Научная статья на тему 'Двумерная локализация атомных населенностей в четырехуровневых квантовых системах'

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

CC BY
111
28
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПРОСТРАНСТВЕННАЯ ЛОКАЛИЗАЦИЯ / "TRIPOD"-СХЕМА / АТОМНЫЕ НАСЕЛЕННОСТИ / SPATIAL LOCALIZATION / TRIPOD SCHEME / ATOMIC POPULATIONS

Аннотация научной статьи по физике, автор научной работы — Ефремова Екатерина Александровна, Гордеев Максим Юрьевич, Рождественский Юрий Владимирович

Исследуется один из аспектов фундаментальной задачи взаимодействия лазерного излучения с веществом пространственное перераспределение атомных населенностей под действием полей нескольких бегущих волн. Впервые показана возможность двумерной пространственной локализации атомных населенностей под действием поля только бегущих волн в «tripods-конфигурации квантовых состояний. Три бегущие волны, которые распространяются в одной плоскости под углами 120° друг к другу, формируют в этой плоскости систему стоячих волн, под действием которой и происходит локализация атомных населенностей. При этом степень локализации населенностей, в принципе, может составлять сотые доли длины волны падающего оптического излучения. Показано, что необходимым условием зависимости населенностей от пространственных координат в плоскости XY является возбуждение центрального перехода "tripod"-системы полем разнонаправленных линейно поляризованных бегущих волн, а возникающие при этом двумерные зависимости могут иметь сложные структуры, такие как, например, «двойные кратеры».

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

Похожие темы научных работ по физике , автор научной работы — Ефремова Екатерина Александровна, Гордеев Максим Юрьевич, Рождественский Юрий Владимирович

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

TWO-DIMENSIONAL LOCALIZATION OF ATOMIC POPULATIONS IN FOUR-LEVEL QUANTUM SYSTEMS

Abstr3ct. The paper deals with investigation of one aspect of fundamental problem of laser radiation interaction with the matter. This problem is spatial localization of atomic populations due to fields impact of few running waves. We are the first to propose in our work two-dimensional spatial localization of atomic populations in medium with tripod-like configuration of levels under the field influence of running waves only. Three running waves, propagating along one plane 120o angle-wise to each other, form the system of standing waves in this plane. Atomic populations can be localized in the field of these standing waves. Moreover, the degree of such localization can make up hundredth parts of the wavelength of the incident optical radiation. It is shown that an excitation of the central transition of the tripod-like system using a field of multidirectional linearly polarized running waves is the necessary condition of the population dependence from spatial coordinates in the XY plane. The two-dimensional shapes that appear in this system can have very complicated structure such as “double craters”.

Текст научной работы на тему «Двумерная локализация атомных населенностей в четырехуровневых квантовых системах»

1

ФОТОНИКА И ОПТОИНФОРМАТИКА PHOTONICS AND OPTOINFORMATICS

УДК 535:621.373.826]:539

ДВУМЕРНАЯ ЛОКАЛИЗАЦИЯ АТОМНЫХ НАСЕЛЕННОСТЕЙ В ЧЕТЫРЕХУРОВНЕВЫХ КВАНТОВЫХ СИСТЕМАХ Е.А. Ефремова"1, М.Ю. Гордеевa, Ю.В. Рождественскийa

a

а Университет ИТМО, Санкт-Петербург, Россия, mxmgordeev@gmail.com

Аннотация. Исследуется один из аспектов фундаментальной задачи взаимодействия лазерного излучения с веществом - пространственное перераспределение атомных населенностей под действием полей нескольких бегущих волн. Впервые показана возможность двумерной пространственной локализации атомных населенностей под действием поля только бегущих волн в «tripod»-конфигурации квантовых состояний. Три бегущие волны, которые распространяются в одной плоскости под углами 120° друг к другу, формируют в этой плоскости систему стоячих волн, под действием которой и происходит локализация атомных населенностей. При этом степень локализации населенно-стей, в принципе, может составлять сотые доли длины волны падающего оптического излучения. Показано, что необходимым условием зависимости населенностей от пространственных координат в плоскости XI является возбуждение центрального перехода «tripod»-системы полем разнонаправленных линейно поляризованных бегущих волн, а возникающие при этом двумерные зависимости могут иметь сложные структуры, такие как, например, «двойные кратеры».

Ключевые слова: пространственная локализация, «tripod»-схема, атомные населенности.

a ITMO University, Saint Petersburg, Russia, mxmgordeev@gmail.com

Abstract. The paper deals with investigation of one aspect of fundamental problem of laser radiation interaction with the matter. This problem is spatial localization of atomic populations due to fields impact of few running waves. We are the first to propose in our work two-dimensional spatial localization of atomic populations in medium with tripod-like configuration of levels under the field influence of running waves only. Three running waves, propagating along one plane 120o angle-wise to each other, form the system of standing waves in this plane. Atomic populations can be localized in the field of these standing waves. Moreover, the degree of such localization can make up hundredth parts of the wavelength of the incident optical radiation. It is shown that an excitation of the central transition of the tripod-like system using a field of multidirectional linearly polarized running waves is the necessary condition of the population dependence from spatial coordinates in the XY - plane. The two-dimensional shapes that appear in this system can have very complicated structure such as "double - craters".

Keywords: spatial localization, tripod - scheme, atomic populations.

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

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

Следует особо подчеркнуть, что в англоязычной научной литературе вместо термина «локализация атомных населенностей» («localization of atom population») используется термин «локализация атомов» («atom localization»), что ведет к определенному несоответствию между физическим явлением и его названием. Дело в том, что в данном случае речь идет именно о локализации населенностей, т.е. в результате взаимодействия атома с полем оптического излучения происходит пространственное перераспределе-

TWO-DIMENSIONAL LOCALIZATION OF ATOMIC POPULATIONS IN FOUR-LEVEL QUANTUM SYSTEMS E.A. Efremovaa, M.Yu. Gordeeva, Yu.V. Rozhdestvenskya

1

Введение

ние населенностей атомных состояний. Таким образом, начальное пространственное распределение атомов w(x, г = 0) с учетом всех населенностей внутренних состояний р,. (х, г = 0) до w¡n (х) и после взаимодействия wf (х) сохранит свой вид, в то время как пространственные зависимости населенностей после взаимодействия р,. (х, г) изменятся:

win (х) = w( х, г = 0) = £ р,. (х, г = 0) = £ р,. (х, г) = w( х, г) = wf (х).

/ /

Выше предполагалось, что под действием оптического излучения атом не меняет своего трансляционного состояния. Другими словами, атом считается настолько тяжелым, что можно пренебречь энергией отдачи Ел = Й2к2 /2т, (к = 1/Х) (где И - постоянная Планка, к - волновое число излучаемого атомом фотона, т - масса излучающего атома, X - длина волны излучаемого атомом фотона) по сравнению с кинетической энергией атома ЕК = р2 /2т (р = ту).

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

, -1 жизни г ~ у :

5х > у-1 = - Е1 ~3-10-31, п Йу

где считалось, что ЕК ~10-2 Йу для сильных электродипольных оптических переходов в атомах, у - скорость спонтанной релаксации атома, 1 - длина волны падающего излучения. Влияние скорости отдачи на конечную пространственную ширину локализованных населенностей также существенно при попытке выделить из общего ансамбля только атомы, которые находятся в определенном внутреннем состоянии, поскольку такое выделение также может быть реализовано с помощью оптических методов.

Таким образом, пространственная локализация населенностей внутренних состояний атома действительно является результатом неоднородной на длине волны оптического излучения оптической накачки. Физически ясно, что основной проблемой здесь является выбор схемы взаимодействия атомов с полем излучения. При этом выбор схемы взаимодействия весьма ограничен, поскольку для эффективной оптической накачки на нижние состояния системы желательно иметь только один возбужденный верхний уровень. В результате все рассмотренные схемы атомных состояний, использованные для исследований локализации населенностей, в той или иной степени основываются на трехуровневой Л -системе. Так, в работах [1-3] предложена одномерная локализация населенностей для атомов в Л-конфигурации. При этом обобщение на случай двух пространственных измерений нетривиально вследствие того, что в уравнениях для элементов матрицы плотности следует явно учитывать поляризации световых волн, что выражается в учете правил отбора для действующих световых волн. По правилам отбора на переходе с разностью магнитных моментов ДМ=±1 может действовать только поле с круговой поляризацией, в случае же, когда разность магнитных моментов между уровнями ДМ=0, действовать на переходе может лишь поле с линейной поляризацией. В связи с этим для того, чтобы использовать поля как с линейной, так и с круговой поляризацией для локализации населенностей в двух измерениях, в [4] было предложено использовать «Мро&ьконфигурацию атомных состояний.

В настоящей работе впервые представлено получение пространственной двумерной локализации населенностей под действием оптических полей бегущих волн также в четырехуровневой «Мро&ь системе. При этом три бегущие волны, которые распространяются в одной плоскости под углами 120° друг к другу, формируют стоячие волны, которые, взаимодействуя со средним переходом «Мро&ь системы (рис. 1, а), обеспечивают пространственную локализацию атомных населенностей.

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

Основные уравнения

Рассмотрим теперь подробно энергетическую схему уровней атома в «Мро&ьконфигурации (рис. 1, а). В этом случае атомная система состоит из трех нижних состояний, оптические переходы между которыми запрещены, и одного верхнего уровня. На оптически разрешенных переходах | т >-1 4 >, т = 1 - 3, действуют поля с частотами Раби g1, g2 и g3 и отстройками от резонансных значений частот переходов Д1, Д2, Д3 соответственно. На рис. 1, б, показана ориентация полей. Видно, что три бегущие волны с одинаковыми частотами Раби g2, действующие на переходе |2 > - |4 >, распространяются в плоскостиХУ под углом 120° друг к другу, в то время как поля на переходах |1 >-| 4 >, |3 >-| 4 > с часто-

тами Раби g1 и g3 распространяются в отрицательном и положительном направлении оси Z соответственно. В результате полное поле, которое действует на «tripods-систему (рис. 1), можно записать в виде

E = e2E2 (cos (га2/ - k!,r) + cos (ra2i - k2r) + cos (ra2i - k2r)) + e1 E1 cos (ra1i + к1 z) + e3E3 cos (ra3i - к3z), (1) где волна с ю1, к1 распространяется в отрицательном, волна с ю3, к3 - в положительном направлении оси Z, а в плоскости XY действуют три волны с одинаковыми частотами, но разными направлениями волновых векторов k'2 л k2 = 120o (i Ф j; i, j = 1,2,3).

Рис. 1. (а): Схема квантовых уровней атома в «ШроС»-конфигурации. Переходы | т >-| 4 >, (т = 1 -3) разрешены в дипольном приближении, а переходы между нижними состояниями запрещены. Частоты Раби световых волн и отстройки возбуждающих полей относительно частот переходов ют4, (т = 1 -3) обозначены gm и Дт (т = 1 -3) соответственно. (б): Взаимная ориентация оптических полей. В плоскости

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

При этом, согласно правилам отбора, вектора е13 задают круговую поляризацию в положительном и отрицательном направлении соответственно, а е2 определяет линейную поляризацию.

Далее для описания динамики системы на рис. 1 используем уравнения для элементов матрицы плотности р.. (х, у, г) в виде

щ. = [ н, р]. + /г. р. (2)

с гамильтонианом в представлении взаимодействия в виде Н = Н0 +¥ , где Н0 задает внутреннее со-

1 3

стояние системы без возмущения, а V = —V (ё/Ч • е..)Е1 определяет взаимодействие с полем оптического излучения (1) для переходов | т >-|4 >, т = 1 -3, с матричным элементом оператора дипольного взаимодействия d/4. В выражении (2) матрица Г.. задает скорости релаксации элементов р. (х, у, г). При

этом скорость релаксации диагональных матричных элементов (т.е. населенностей) определяется естественной шириной 2у верхнего возбужденного состояния системы 2у = у1 + у2 + у3, а скорости релаксации недиагональных матричных элементов Г. (/ Ф .) могут, наряду со скоростью естественного распада, учитывать и другие возможные механизмы уширения (столкновения, конечную ширину спектра возбуждающих полей и пр.). В результате уравнения для элементов матрицы плотности ¡5.. (х, у, г) системы

уровней в «tripod»-конфигурации имеют следующий вид:

Ф11 = ^(Р14 - Р41 ) + "У1Р44, ф22 = g2(Р24 - Р42 ) + /У2Р44 , Ф33 = gз(Pз4 - Р43) + "УзР44 , ф44 = gl(P41 - Р14 ) + g2(Р42 - Р24 )+ gз(P43 - Р34) -/(Т1 + У2 + У3 )P44,

."Р14 = gl(Pп - Р 44 ) + g2Pl2+ gзPl3 + (Д1 /Г14 )Р14 , (3)

ф24 = g2 (Р22 - Р44 ) + glP21+ gзP23 + (Д 2 -/Г24)P24, Фз4 = gз(p33 - Р44) + &Рз1 + g 2Р32 + (Д3 /Г34 )Рз4 , /Р12 = g2Pl4 - glP42 + (Д1 -Д2 )Р12 , "Р13 = gзPl4 - g 1Р43 + (Д1 -Д3)Pl3, "Р23 = gзP24 - g2P43 + (Д2 - Д3 )P23,

где ру- = р* (, Ф У). При выводе системы уравнений (3) мы пренебрегли членами, содержащими временные осцилляции на удвоенной оптической частоте (резонансное приближение) и использовали так называемое приближение вращающейся волны, которое состоит в замене недиагональных матричных элементов: оптических рт4 = рт4е'Л"' (т = 1 - 3) и низкочастотных когерентностей

р12 = р.е'^ -Д2)' , р1з = рпе'(Д1 -Дз)' , р2з = р2зе'(Д2-Дз)'.

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

1 л/3 1 ^/з

к2г =-к2y, к\г = ^к2у- —к2х к32г = ^к2у + —к2х. (4)

Таким образом, выражения (4) явно демонстрируют, что в плоскости ХУ три бегущие разнонаправленные волны действительно задают систему стоячих волн, которые обеспечивают координатную зависимость населенностей состояний квантовой системы.

Результаты вычислений и обсуждения

Система уравнений для элементов матрицы плотности (3) полностью определяет динамику квантовой системы в «Мро&ьконфигурации для поля (1). Для исследования пространственной зависимости населенностей (т.е. диагональных элементов матрицы плотности) мы используем стационарное решение уравнений (3), которое может быть получено в пренебрежении временными производными в левых частях по сравнению с членами в правых частях, содержащих скорости распада. Физически это означает, что рассматриваются значения элементов матрицы плотности на временах, значительно превышающих t >> у , I у .

На рис. 2 представлены пространственные зависимости населенностей всех состояний системы в «Мро&ьконфигурации. Видно, что в нашем случае действительно имеет место двумерная (в плоскости ХУ) локализация населенностей для всех состояний системы. При выбранных параметрах только населенность р11 демонстрирует пространственную зависимость в виде «холмов», а для остальных населенностей имеют место «кратеры». При этом в определенных точках плоскости ХУ максимальные значения населенностей первых двух состояний могут достигать единицы. Это означает, что в этих точках сосредоточена вся населенность квантовой системы.

С другой стороны, среди пространственных структур на рис. 2 наиболее интересными являются зависимости населенностей уровней |2 >, |3 > и |4 > р22,р33,р44. Так, пространственные зависимости в плоскости ХУ населенностей третьего и четвертого уровней р33, р44 демонстрируют двойную структуру -так называемые «двойные кратеры», в отличие от обычного «кратера» (рис. 2, в). При этом вид сверху наглядно показывает сложную структуру пространственных зависимостей (рис. 2, е, з). Наибольший интерес здесь представляет график населенности второго уровня р22 . Полученный кратер имеет крайне узкие стенки, что означает высокую степень локализации, которая для данных параметров составляет 0,04 длины волны падающего излучения.

Заключение

Исследована возможность двумерной локализации населенности в «Мро&ьсистеме при взаимодействии с полем бегущих световых волн. Найдено, что необходимым условием зависимости населенно-стей от пространственных координат в плоскости ХУ является возбуждение центрального перехода «М-ро&»-системы полем разнонаправленных линейно поляризованных бегущих волн. Возникающие при этом двумерные зависимости могут иметь сложные структуры, такие как, например, «двойные кратеры».

Отметим, что для реализации двумерной локализации не важен способ получения пространственной зависимости линейно поляризованного поля, резонансного центральному переходу «Мро&ьсистемы. Исходя из этого, полученное двумерное распределение населенностей можно рассматривать как своеобразную визуализацию распределения интенсивности светового поля. Характеристики такого «экрана» определяются параметрами условий, при которых можно наблюдать перераспределение населенностей.

В нашем случае инертность, т.е. быстрота отклика на воздействие, характеризуется временем опТ -1 -7

тической накачки t —— = у ~10 с (8 = у), а минимальный размер «зерна», т.е. минимальный размер

8

области, который может быть отражен таким «экраном», составляет 8х > уя у-1 ~3 10 3 X ~2нм . При

этом возможность получения такого разрешения непосредственно связана со значениями скоростей атомов и может быть реализована на атомных ансамблях, охлажденных до температуры энергии отдачи Тя = ЕК / кв = й2к2 / 2ткв , где кв - постоянная Больцмана.

1,0 0,8 0,6

Р 0,4 Pli '

ky

-6 -4 -2 0 2 4 6 kx

б

-6 -4

-2 0 kx

jy

1,000 0,750 0,625 0,375 0,250 0,000

M 1,000 0,750 0,625 2 0,375 0,250 4 0,000

2 4 6

0,18 0,14 0,1 Рзз0,06 0,02

-4 kx

ky

4

д

0,18 0,14 0,10 P44 0,06 0,02

-4 kx

6 4

2 1,000

0 0,750

£ 0,625

-2 0,375 0,250

-4 0,000

I ' г

-6 -4 -2 0 kx

2 4

6

4

2 —1,000

1 0 ky 0,750

0,625

-2 0,375

--24 0,250 0,000

r

-6 -4

-2 0 2 4 6 kx

ж з

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

Рис. 2. Координатные зависимости населенностей состояний «ШроС»-системы: р11(х,у) (а); р22(х,у) (в); р33(х,у) (д); р44(х,у) (ж). Справа (б, г, з, е) представлен вид сверху для каждой населенности. Параметры взаимодействия и релаксационные константы системы были выбраны в виде Y1,2,з = Y, Г14 = Г24 = Г34 = 1 = Y, 92 = 3Y, дз = 2Y, Д1 = Y, = 4Y, Д1 = -2Y

6

а

в

г

6

6

е

Литература

1. Qamar S., Zhu S.-Y, Zubairy M.S. Atom localization via resonance fluorescence // Physical Review A -Atomic, Molecular, and Optical Physics. 2000. V. 61. N 6. P. 1-5.

2. Agarwal G.S., Kapale K.T. Subwavelength atom localization via coherent population trapping // Journal of Physics B: Atomic, Molecular, and Optical Physics. 2006. V. 39. N 17. P. 3437-3446.

3. Xu J., Hu X.-M. Sub-half-wavelength localization of an atom via trichromatic phase control // Journal of Physics B: Atomic, Molecular, and Optical Physics. 2007. V. 40. N 7. P. 1451-1459.

4. Ivanov V., Rozhdestvensky Y. Two-dimensional atom localization in a four-level tripod system in laser fields // Physical Review A - Atomic, Molecular, and Optical Physics. 2010. V. 81. N 3. Art. N 033809.

5. Qamar S., Zhu S.-Y., Zubairy M.S. Precision localization of single atom using Autler-Townes microscopy // Optics Communications. 2000. V. 176. N 4. P. 409-416.

6. Ghafoor F., Qamar S., Zubairy M.S. Atom localization via phase and amplitude control of the driving field // Physical Review A - Atomic, Molecular, and Optical Physics. 2002. V. 65. N 4. P. 0438191-0438198.

7. Paspalakis E., Knight P.L. Localizing an atom via quantum interference // Physical Review A - Atomic, Molecular, and Optical Physics. 2001. V. 63. N 7. P. 065802/1-065802/4.

8. Paspalakis E., Terzis A.F., Knight P.L. Quantum interference induced sub-wavelength atomic localization // Journal of Modern Optics. 2005. V. 52. N 12. P. 1685-1694.

9. Liu C., Gong S., Cheng D., Fan X., Xu Z. Atom localization via interference of dark resonances // Physical Review A - Atomic, Molecular, and Optical Physics. 2006. V. 73. N 2. Art. N 025801.

10. Carreno F., Anton M.A. Gradient echo memory in a tripod-like dense atomic medium // Optics Communications. 2010. V. 283. P. 4787-4795.

11. Gornyi M.B., Matisov B.G., Rozhdestvenskii Yu.V. Coherent population trapping in an optically dense medium // Sov. Phys. JETP. 1989. V. 68. N 4. P. 728-732.

12. Агапьев Б. Д., Горный М.Б., Матисов Б.Г., Рождественский Ю.В. Когерентное пленение населенно-стей в квантовых системах // Успехи физических наук. 1993. Т. 163. № 9. С. 1-36.

13. Скалли М.О., Зубайри М.С. Квантовая оптика : Пер. с англ. М.: ФИЗМАТЛИТ, 2003. 510 с.

14. Стенхольм С. Основы лазерной спектроскопии : Пер. с англ. М.: Мир, 1987. 312 с.

15. Аллен Л., Эберли Дж. Оптический резонанс и двухуровневые атомы : Пер. с англ. М.: Мир, 1978.

223 с.

Ефремова Екатерина

Александровна

Гордеев Максим Юрьевич

Maxim Yu. Gordeev

Рождественский Юрий Владимирович

Ekaterina A Efremova

Yuri V. Rozhdestvensky

кандидат физико-математических наук, доцент, СПбГУ, Санкт-Петербург, Россия, Efremova.kate@gmail.com

аспирант, Университет ИТМО, Санкт-Петербург, Россия, mxmgordeev@gmail. com

доктор физико-математических наук, ведущий научный сотрудник, Университет ИТМО, Санкт-Петербург, Россия, rozd-yu@mail.ru

PhD, Associate professor, Saint Petersburg State University, Saint Petersburg, Russia, Efremova.kate@gmail.com

postgraduate, Saint Petersburg State University, Saint Petersburg, Russia, mxmgordeev@gmail.com

D.Sc., Senior scientific researcher, ITMO University, Saint Petersburg, Russia, rozd-yu@mail.ru

Принято к печати 08.05.14 Accepted 08.05.14

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