Научная статья на тему 'Двухполевые космологические модели с повторным ускоренным расширением Вселенной'

Двухполевые космологические модели с повторным ускоренным расширением Вселенной Текст научной статьи по специальности «Физика»

CC BY
72
18
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИНФЛЯЦИЯ / СКАЛЯРНОЕ ПОЛЕ / КИРАЛЬНЫЕ КОСМОЛОГИЧЕСКИЕ МОДЕЛИ / INFLATION / SCALAR FIELD / CHIRAL COSMOLOGICAL MODELS

Аннотация научной статьи по физике, автор научной работы — Фомин И.В.

Целью настоящей работы является рассмотрение космологических моделей, содержащих два скалярных поля. Одно поле является инфлатоном, второе представляет собой источник наблюдаемого ускоренного расширения Вселенной в настоящее время. Рассматриваемой в работе модели ставится в соответствие киральная космологическая модель, для которой определена метрика внутреннего пространства полей.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Two-Field Cosmological Models with a Second Accelerated Expansion of the Universe

The purpose of this paper is to consider cosmological models that contain two scalar fields. One field is inflaton, the other is the source of the observed accelerated expansion of universe at the present time. For the model that is considered in this work the corresponding chiral cosmological model was found and the metric of the space of the internal field is defined.

Текст научной работы на тему «Двухполевые космологические модели с повторным ускоренным расширением Вселенной»

Двухполевые космологические модели с повторным ускоренным расширением

Вселенной

И. В. Фомин"

Московский государственный технический университет имени Н. Э. Баумана, кафедра физики. Россия, 105005, Москва, 2-я Бауманская ул., д. 5.

Статья поступила 16.03.2018, принята к публикации 11.04.2018.

Целью настоящей работы является рассмотрение космологических моделей, содержащих два скалярных поля. Одно поле является инфлатоном, второе представляет собой источник наблюдаемого ускоренного расширения Вселенной в настоящее время. Рассматриваемой в работе модели ставится в соответствие киральная космологическая модель, для которой определена метрика внутреннего пространства полей.

Ключевые слова: инфляция, скалярное поле, киральные космологические модели. УДК: 524.8. РАСЗ: 98.80.-k.

ВВЕДЕНИЕ

Стандартная теория космологической инфляции, основанная на теории гравитации Эйнштейна, успешно описывает состояние ускоренного расширения Вселенной на ранних стадиях эволюции. Также стандартная инфляционная космология объясняет происхождение первичных неоднородно стей и предсказывает их спектр [1-3].

Механизм начального инфляционного сценария, подразумевающего стадию ускоренного расширения ранней Вселенной, может быть описан в предположении о существовании скалярного поля (инфлатона), также ускоренное расширение современной Вселенной [4] связывают с существованием темной энергии, которая может быть объяснена новой фундаментальной постоянной (космологической константой), модификацией гравитации на сверхбольших масштабах, наличием нового легкого поля (например, газа Чаплыгина), также квинтэссенцией (скалярным полем с положительной кинетической энергией), фантомной энергией (скалярным полем с отрицательной кинетической энергией) и к-эссенцией (скалярным полем с нетривиальной кинетической энергией) [5, 6]. Отметим, что вышеперечисленные скалярные поля, моделирующие темную энергию, также рассматриваются в качестве инфлатона в различных моделях космологической инфляции [7].

Модели к-эссенции могут быть получены из действия Борна—Инфельда в теории струн [8] и предлагались для описания стадии космологической инфляции и темной энергии [9, 10]. Альтернативным механизмом для получения инфляционного и повторного ускорения Вселенной также являются модели квинтэссенции [11], в которых основной вклад в динамику дает потенциал скалярного поля.

Различие между моделями квинтэссенции и к-эс-сенции заключается в том, что последний класс моделей содержит неканонические кинетические члены в лагранжиане. В этом смысле модели квинтэссенции можно рассматривать как частный случай к-эссенции. Также обсуждаются модели к-эссенции, в которой лагранжиан содержит только кинетическую часть, то есть функцию производных скалярного поля, и не зависит явно от самого поля [12].

Для реализации инфляционного сценария и наблюдаемого ускоренного расширения Вселенной также рассматриваются модели с двумя скалярными полями, одно из которых является инфлатоном, второе определяет темную энергию. В рамках этого подхода возможны различные виды кинетического и потенциального взаимодействия между полем инфлатона и скалярным полем, моделирующим темную энергию на инфляционной стадии, а также взаимодействия темной материи и темной энергии на последующих стадиях эволюции Вселенной [13-16].

В настоящей работе рассматривается модель космологической инфляции с двумя взаимодействующими скалярными полями. Первое убывающее поле определяет космологическую динамику на ранней стадии эволюции Вселенной. Темная энергия определяется полем к-эссенции с линейным лагранжианом кинетической части, что позволяет сопоставить ему некоторое эффективное скалярное поле. Данному подходу ставится в соответствие киральная космологическая модель [17], объединяющая взаимодействия кинетического и потенциального типа между скалярными полями.

1. МОДЕЛИ КОСМОЛОГИЧЕСКОЙ ИНФЛЯЦИИ С ОДНИМ СКАЛЯРНЫМ ПОЛЕМ

Запишем действие для модели инфляции со скалярным полем х в виде

S = j dxAJ-~g(R + 1 + Y(х)),

(1)

где У (х) — потенциал скалярного поля х.

Космологическая динамика определяется уравнениями Эйнштейна и полевым уравнением, которые для действия (1) в пространственно плоской Вселенной Фридмана—Робертсона—Уокера (ФРУ) записываются следующим образом [18]:

3 H2 = 1X2 + Y (х),

н = -1X2, 2

•• , ои- , dY(X) п

X + 3HX+--д— = 0.

дх

(2)

(3)

(4)

E-mail: ingvor@inbox.ru

Классификация методов точных решений этой системы уравнений, также метод генерирования новых

а

космологических решений из известных и сопоставления с приближенными решениями рассматривались в работе [18].

Действие для моделей к-эссенции с полем 9 записывается как [19]

5 = | + р(9,Х)),

(5)

3Н2 =р-2Хр х,

н = Хр,х,

1 д

Х = 2Хр,х.

р(9,Х ) = К (9)Х + К (9).

5 =у ^ж4^^^1^+2ф

где кинетическое взаимодействие между полями ф и 9 можно рассматривать как взаимодействие между полями ф и х, определяемое потенциалом и[ф, 9(х)].

Уравнения динамики для действия (15) во Вселенной ФРУ записываются следующим образом:

3Н2 = 1 ф2 + V(ф) + р - 2Хр,х + и[ф, 9(х)], (16)

где р(9, X) — функция скалярного поля 9 и кинетического члена Х = — 9,м9'м, кинетическая энергия поля, в таком случае определяется как —Х/2.

Уравнения космологической динамики для таких моделей

д

ф + 3Нф + — (V(ф) + и[ф, ^(х)]) = 0,

[р,х + 2Хр,хх ] 9 + 3Нр,х 9+ 1д

(17)

(6) (7)

[р,х + 2Хр,хх] 9 + 3Нр,х9 +2 д^ [р — 2Хр,х] = 0.

(8)

Соответствие между этими моделями определяется следующим условием [19]:

3Н2 + Н = У (х) = р(9, X) — Х^^Х1. (9)

В таком случае скалярное поле х будет функцией поля 9 и кинетического члена Х. Из уравнений (3) и (7) запишем

(10)

+ 2 д9 (р — 2Хр,х + и[ф,9(х)])=0. (18)

С линейным лагранжианом поля к-эссенции р(<9>, Х) = К(9)Х + К(9) уравнения принимают вид

3Н2 = 1 ф2 + V(ф) — К(9)992 + К(9) + и[ф, 9(х)],

(19)

д

ф + 3Нф + — (V (ф) + и [ф,9(х)]) = 0, (20) 1 + 92 д

9 + 3^9+^9) д9 (К(9) + и[ф,9(х)]) = 0. (21)

Из соответствия между лагранжианом к-эссенции р(9, Х) = К(у>)Х + К(9) и лагранжианом скалярного поля £ = 2х2 — У (х) с обычным кинетическим членом получим

После подстановки х в уравнение (4), с учетом условия (9), получаем

[р,х + Хр,хх] 9 + 3Нр,х9 +2 д9 [р — 2Хр,х] = 0.

(11)

Уравнения (8) и (11) записываются одинаково при условии р,хх = 0.

Таким образом, лагранжиан р(у, Х) определяется кинетическим членом Х и функцией скалярного поля

К = К (9):

3Н2 = 2 ф2 + 2 X2 + w (ф, х),

ф+3Нф+ д^ (ф, х) 0 ф + 3Нф +--—-= 0,

.. + 3Н. + дШ (ф, х) 0

х + 3Нх9 +----= 0,

дх

(22)

(23)

(24)

где Ш(ф,х) = V(ф) + У(х) + и(ф,х), и[ф,9(х)] =

= и (ф, х).

Эта система может быть записана в виде [13]

(12)

Н 2(ф, х) = 1Ш (ф,х)+

Для рассматриваемого лагранжиана р(9, Х) из уравнений (9) и (10) получим соотношения:

У (х) = К (9), (13)

X = ±1 \/ 2К (9^9. (14)

Исходя из этих соотношений, будем рассматривать уравнения динамики поля к-эссенции как уравнения динамики скалярного поля х.

2. УРАВНЕНИЯ КОСМОЛОГИЧЕСКОЙ ДИНАМИКИ

Рассмотрим двухкомпонентную космологическую модель, определяемую действием вида

(я 1

АдН(ф,хП 2 / дН(ф, х) V дф ^ + ^ дх

ф 2 дН(ф, х) . 2 дН(ф, х)

ф = —2--, X = —2-

2

+ 3

дф

дх

, (25)

(26)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Теперь воспользуемся следующим представлением параметра Хаббла Н(ф,х) = Нх(ф) + Н2(х), которое часто используется при анализе космологических моделей, содержащих более одного скалярного поля (см., например, [13]).

В таком случае система уравнений (22)-(24) приводится к виду

[Н(ф)+Н2(х)]2 = 3Ш(ф, х)+

2 ' 2Т

— V(ф) — р(9,Х) — и[ф,9(х)] ), (15)

2

+ 3

дН!(ф^2 + ^дН2(х) и

дф

V дх

ф 2 дНх(ф) . 2 дН2(х) ф = —2——, X = —2-

дф

д

(27)

(28)

2

Далее из уравнения (27) получим выражения для потенциалов самодействия скалярных полей V(ф) и У (х) и потенциала их взаимодействия и(ф, х) через параметр Хаббла

3H2

3Я22

У1(ф) + 1 ф2 = V (ф) + 2

V2(x) + 1 XX2 = Y (х) + 2

6Я1Я2 = U (ф,х).

dH1

(29)

(30)

(31)

Решения этой системы уравнений будут определять космологическую динамику для конкретных космологических моделей.

3. МОДЕЛЬ КОСМОЛОГИЧЕСКОЙ ИНФЛЯЦИИ

Рассмотрим потенциалы скалярного поля V(ф) = = V) ехр ^ —^ и поля к-эссенции К (у) = .

Из уравнений (13)-(14) получаем потенциал У(х) = = В2 ехр (—-X) эффективного скалярного поля х > 0, положительное значение поля выбирается для соответствия рассматриваемой модели однополевой космологической модели с квинтэссенцией [20] в частном случае. Таким образом, потенциал взаимодействия полей будет и(ф, х) = Л ехр ^ — — . Потенциал Ш(ф, х) определяется как

Ш(ф, х) = V) ехр ( А) + В2 ехр (—+

+ЛЧ А—!)• (32)

Для случая А = В, ф = х и V) = В2 получим модель инфляции с одним скалярным полем и потенциалом

вида Ш(ф) = V) еов^—^, который рассматривался

ранее для описания динамики Вселенной в работе [20].

В общем случае точные решения системы уравнений записываются следующим образом:

Н1(ф) = Д2 ехр (А) , Н2(х) = В2 ехр (— ,

D2

A2Vo

3A2 - 4' + R2 •

H(ф' х) = D2 exp (A) + B2 exp (-1) ,

B2 =

Л

(33)

(34)

ф(г) = A ln

/ A2

V2D2(t + C)

^(t) = [2(t + C)]-B/2' A2 B2

6D2'

' x(t) = B ln (2(t + C))' (35)

H (t)

+

2D2(t + C) 2(t + C)'

a2 1 л

a(t) = a1a2 = a0(t + C) 2d2 +12d2

(36)

где A, B, C, D, Л — постоянные свободные параметры модели.

Теперь рассмотрим случай, когда ф ^ х, и масштабный фактор на стадии инфляции запишем следующим образом:

j(t) ~ a0(t + C) 2D2 .

(37)

Далее определим параметры медленного скатывания

Ö = е -

H

"H2

е

aa 1 —^

a2 Ha

2Не

2HH

2D2 A2 2D2 A2

const'

= const.

(38)

(39)

Рассматривая квантовые флуктуации скалярного поля ф, исходя из условия ¿ф ^ ¿х, в линейном порядке теории космологических возмущений, в соответствии с результатами, полученными в работах [21, 22] (отметим, что в работе [21] расчет параметров космологических возмущений производился на основе точных решений уравнений космологической динамики, без использования приближения медленного скатывания), определим основные параметры космологических возмущений на пересечении радиуса Хаббла (к = аН):

PR(k) = -

Н 4

h2

A6

8п2Н 8п2е 64D6n2(t + C)2 '

h2

a4

64D4 n2(t + C )2:

ns(k) - 1

4H - ^д

H + H2

J - 2е 1 - е

2

4D2

ng (k)

2H

A2 - 2D2 4D2

H + H2

1-е

Pg . 8D

r = —— = 4е =-

Г Pr 4е A2

2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

A2 - 2D2 4(ns - 1) ns - 3 '

(40)

(41)

(42)

(43)

(44)

где Тп, Тд, п^ и по — спектральные индексы (наклоны спектров) скалярных и тензорных возмущений, г — тензорно-скалярное отношение.

Далее определим число е-фолдов как функцию времени

N (t)=/

Hdt = 2D2 ln(t + C)'

(45)

причем различие значений числа e-фолдов на завершении стадии инфляции и в ее начале оценивается как N(te) - N(ti) « 50 - 60 [23].

Также запишем измеренные параметры космологических возмущений и ограничения на их значения по данным наблюдений анизотропии реликтового излучения спутником PLANCK [23]

109PS = 2.142 ± 0.049' ns = 0.9667 ± 0.0040'

Pt = rps'

r < 0.112.

(46)

Верификация рассматриваемой модели на основе наблюдательных ограничений может быть проведена посредством выбора свободных параметров А и Д, также время завершения стадии космологической инфляции ¿е можно оценить из выражения (45), рассматривая N =0.

2

2

r=

0.25 -

0.20

0.15

0.10

0.05

0.00

т--1-1-1-1-

Planck TT+lowP+lensing (ANeff = 0.39)

\ \ ACDM Planck TT+lowP

+lensing+ext

0.25 -

0.20

0.15

0.10

0.05

0.00

T-1-1-

Planck TT+lowP Planck TT+lowP+BKP +lensing+ext

1.00

1.00

ns, наклон спектра скалярных возмущении

Рисунок. Зависимость зависимость r = r(ns) для модели с потенциалом (32) на завершении стадии инфляции, наложенная на данные наблюдений спутника PLANCK с доверительными вероятностями 68% и 95% (использован рис. 21 из [23], в которой также представлено подробное описание анализа данных наблюдений). Диаграммы красного цвета соответствуют комбинации наблюдений анизотропии реликтового излучения для больших l ^ 30 и малых значений мультипольных коэффициентов вращательно-инвариантной величины Cl, в терминах которой обычно производится построение спектров мощности космологических возмущений. Синий цвет соответствуют комбинированному анализу предыдущих и дополнительных наблюдений по барионным акустическим осцилляциям, линзированию, измерению параметра Хаббла Ho. Серый цвет соответствует изменению диаграммы с учетом дополнительных релятивистских степеней свободы с ANeff = 0.39. На диаграмме справа дополнительно учитываются измерения B-моды поляризации реликтового излучения BICEP2/Keck Array и PLANCK (выделено зеленым цветом) [23])

Тем не менее эффективным методом проверки соответствия космологических моделей наблюдениям, независимо от выбора свободных параметров, является построение диаграммы г = г(пз) и ее наложение на ограничения, обусловленные комбинацией наблюдений анизотропии и поляризации реликтового излучения, барионных акустических осцилляций, данных по гравитационному линзированию и измерениям значения параметра Хаббла Н0 в современную эпоху [23].

На рисунке приведена зависимость г = г(пз), из которой видно, что рассматриваемая модель соответствует наблюдаемым ограничениям на параметры космологических возмущений лучше, чем модель с квадратичным потенциалом V гс ф2, и близка к модели с линейным потенциалом V гс ф.

Постинфляционная космологическая динамика включает стадию преобладания излучения с масштабным фактором а(£) = а0(£ + С)1/2 и стадию преобладания вещества а(£) = а0(£ + С)2/3.

В современную эпоху начинается стадия повторного ускоренного расширения Вселенной с масштабным фактором

a(t) ~ a0(t + C)3

2 + -

3 +

(47)

для которой поле к-эссенции 9 (или эффективное скалярное поле ) моделирует темную энергию, определяющую динамику Вселенной на этой стадии.

4. КИРАЛЬНАЯ КОСМОЛОГИЧЕСКАЯ МОДЕЛЬ

Нелинейные сигма-модели с потенциалом взаимодействия (киральные космологические модели) часто используют для описания различных эпох эволюции Вселенной (см. обзор [17]). В таких моделях присутствует внутреннее пространство целей (полей), метрика которого определяет взаимодействие между скалярными полями.

Запишем действие для нелинейной сигма-модели [17]

5 = У ¿ж4^ (Я + 2дмФ3д"" — Пф, х)) ,

(48)

где У(ф, х) = V(ф) + У(х), — метрический тензор плоского пространства Фридмана—Робертсона— Уокера, Л 3 — метрический тензор пространства полей,

ф1 = ф ф2 = х.

Метрика пространства полей определяется следующим образом

ds2 = hij .

(49)

В работе [17] было показано, что система уравнений Эйнштейна и полевых уравнений в плоском пространстве Фридмана—Робертсона—Уокера для двухкомпо-нентной сигма-модели приводятся к виду

1 Лиф2 + Л^ф* + 2^22 (¿)х2 = —Н, (50)

(51)

y(t) = 3H2 + H.

Сопоставляя с уравнениями (27)-(28), для случая hii = 1, h22 = 1,hi2 = hi2(0, x) получим

hi2 (Ф, X) = 3 H1(0)H2 (X'

2 ( 9И1(ф)\ (сЭИ2 (x) ^ ' 1 d^^ dx J

(52)

В результате компоненты метрики пространства полей

3

hii = 1, h22 = 1, hi2 = h2i = -2 AB, (53)

1

h- = i - - 2AB

hij = 1 -3 AB

1

л

2

Таким образом, рассматриваемой модели космологической инфляции со скалярным полем ф и полем к-эссенции у соответствует киральная космологическая модель с потенциалами самодействия V (ф) = V) ехр ^—^, К (у) = и компонентами (53) метрики внутреннего пространства полей, определяющего взаимодействие между ними.

Поскольку тензор Н^ является постоянным, его возмущения ¿Н^- =0 и взаимодействие между полями ф и у не учитываются при расчете параметров космологических возмущений.

ЗАКЛЮЧЕНИЕ

В данной работе рассматривалась космологическая модель с двумя взаимодействующими каноническим скалярным полем ф и полем к-эссенции у, которому соответствует эффективное скалярное поле х, поскольку в моделях к-эссенции с линейным лагранжианом существует возможность определить соответствие между полем к-эссенции и эффективным скалярным полем.

Введение в модель второго поля к-эссенции обусловлено наблюдаемым ускоренным расширением Вселенной в настоящее время, причем на стадии космологической инфляции выполняется условие ф ^ х. В рамках рассматриваемой модели поле к-эссенции у уменьшается с течением времени, но эффективно соответствующее ему скалярное поле х растет и является источником ускоренного расширения на больших временах.

Также рассматриваемой космологической модели поставлена в соответствие киральная космологическая модель, для которой был найден вид метрического тензора пространства полей Н^-, определяющий потенциальное и кинетическое взаимодействие между полями. В данном случае тензор пространства полей является постоянным, что исключает влияние взаимодействия между полями на спектральные параметры космологических возмущений.

В рамках данного подхода были рассчитаны параметры космологических возмущений на завершении первой (инфляционной) стадии ускоренного расширения, что позволило сопоставить полученные результаты с данными наблюдений. В линейном порядке теории космологических возмущений, при нивелирова-

нии вклада возмущений поля х и тензора hj, на стадии инфляции рассматриваются только адиабатические возмущения скалярного поля ф, то есть спектр возмущений является гауссовым.

Отметим, что аналогичным образом можно построить космологические модели с каноническим и фантомным полями или так называемые квинтомные модели [17], также модели со скалярным (каноническим или фантомным) и тахионным полями, которые можно эффективно описывать посредством некоторого скалярного поля [24], что является тематикой дальнейших исследований автора в данном направлении.

Работа выполнена при финансовой поддержке РФФИ (гранты № 16-02-00488 A и № 16-08-00618 A.)

СПИСОК ЛИТЕРАТУРЫ

1. Starobinsky A.A. // Phys. Lett. B. 1980. 91. P. 99.

2. Guth A.H. // Phys. Rev. D. 1981. 23, N 2. P. 347.

3. LindeA.D. // Phys. Lett. B. 1982. 108. P. 389.

4. Perlmutter S. et al. // Astrophys. J. 1999. 517. P. 565.

5. Сажин М.В., Сажина О. С. // АЖ 2016. 93, № 4. С. 394.

6. Mukhanov V.F., Vikman A. // JCAP. 2006. 0602. P. 004.

7. Baumann D, Peiris H. // Adv. Sci. Lett. 2009. 2. P. 105.

8. Chimento L.P. // Phys. Rev. D. 2004. 69. 123517.

9. Armendariz-Picon C. et al. // Phys. Rev. Lett. 2000. 85. P. 4438.

10. Chiba Т., Okabe Т., Yamaguchi M. // Phys. Rev. D. 2000. 62. 023511.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

11. Peebles P.J.E., Vilenkin A. // Phys. Rev. D. 1999. 59. 063505.

12. de Putter R., Linder E. V. // Astropart. Phys. 2007. 28. P. 263.

13. Byrnes C. T., Tasinato G. // JCAP. 2009. 0908. P. 016.

14. Boehmer C. // AIP Conf. Proc. 2009. 1122. P. 197.

15. Costa F.E.M., Alcaniz J. S., Jain D. // Phys. Rev. D. 2012. 85. 107302.

16. Fomin I. V. // J. Phys. Conf. Ser. 2017. 918, N 1. 012009.

17. Chervon S. V. // Quantum Matter. 2012. 2. P. 71.

18. Фомин И. В., Червон С. В. //Изв. вузов. Физика. 2017. 60, № 3. С. 46.

19. Aguirregabiria J. M., Chimento L. P., Lazkoz R. // Phys. Lett. B. 2005. 631. P. 93.

20. Sahni V., WangL.M. // Phys. Rev. D. 2000. 62. 103517.

21. Chervon S. V., Fomin I. V. // Grav. Cosmol. 2008. 14. P. 163.

22. ЛукашВ.Н. //УФН. 2006. 176. С. 113.

23. Ade P.A.R. et al. [Planck Collaboration] // Astron. Astrophys. 2016. 594. P. A13.

24. Avelino P.P. et al. // Phys. Rev. D. 2010. 82. 063534.

Two-Field Cosmological Models with a Second Accelerated Expansion of the Universe I. V. Fomin

Bauman Moscow State Technical University, Laboratory of Electrodynamics of Moving Media. Moscow 105005, Russia. E-mail: ingvor@inbox.ru.

The purpose of this paper is to consider cosmological models that contain two scalar fields. One field is inflaton, the other is the source of the observed accelerated expansion of universe at the present time. For the model that is considered in this work the corresponding chiral cosmological model was found and the metric of the space of the internal field is defined.

Keywords: inflation, scalar field, chiral cosmological models. PACS: 98.80.-k. Received 16 March 2018.

English version: Moscow University Physics Bulletin. 2018. 73, No. 6. Pp. 696-701.

Сведения об авторе

Фомин Игорь Владимирович — канд. физ.-мат. наук, доцент; e-mail: ingvor@inbox.ru.

i Надоели баннеры? Вы всегда можете отключить рекламу.