CC BY
25
УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Том 152, кн. 2
Физико-математические пауки
2010
УДК 530.145
ОСОБЕННОСТИ ДИССИПАТИВНОЙ ДИНАМИКИ Л-АТОМА ПОД ДЕЙСТВИЕМ КВАНТОВАННОГО ПОЛЯ В СЖАТОМ СОСТОЯНИИ
АЛ. Гельман,, В.А. Миронов
Аннотация
Числеппо квантовым методом Мопте-Карло и аналитически методом Гейзепберга Ланжевена исследована днссипативная динамика трехуровневого Л-атома под действием квантованного электромагнитного поля в режиме электромагпитпо-гшдуцироваппой прозрачности (ЭИП). Подробно рассмотрено влияние уровня сжатия квантованного света па параметры атомной системы, отвечающие за процессы релаксации и шум в полосе прозрачности ЭИП. Показано, что квантовые свойства света, такие, как сжатие, могут передаваться атому в дашгом режиме. Исследована задача распространения квантованного поля в сжатом когерентном состоянии в среде и найдена зависимость уровня сжатия поля па выходе из среды от параметров па входе и свойств среды. Обсуждается применение метода Мопте-Карло к численному моделированию задачи распространения в свете недавних экспериментов по прохождению сжатого вакуума в рассмотренной среде.
Ключевые слова: сжатый свет, квантовые скачки, квантовый шум. квантовая память.
Введение
Недавние успехи в физике квантовой информации [1] привели к развитию таких новых приложений, как квантовые вычисления, квантовая криптография, телепор-тация. квантовые компьютеры. Одна из основных проблем процессы диссипации в атомных системах, предлагаемых для реализации квантовых протоколов. Диссипация вызывает флуктуации и квантовый шум [2] и может значительно повлиять на динамику системы. Существует два основных подхода к анализу этой проблемы: метод оператора плотности и метод Гейзепберга Ланжевена [2. 3]. Оба метода предполагают описание ансамбля систем и дают усредненное описание динамики системы. Развитый за последние годы метод Монте-Карло расчета волновой функции системы [4] позволяет аналитически описывать единичные реализации процессов в квантовых системах (например, один ион в магнитной ловушке) и по-новому взглянуть на физику процессов в системе. Разработанный изначально для аналитического описания диссипативной динамики одиночных квантовых систем, метод стал мощным способом численного расчета квантового кинетического уравнения как в форме Линдблада [2] (в наиболее общем виде уравнения для матрицы плотности в марковском приближении), так и для иемарковских систем.
Взаимодействие квантованного электромагнитного (ЭМ) поля со средой трехуровневых Л-атомов достаточно хорошо изучено за последние несколько лет [5, 8]. Эффект электромагнитно-индуцированной прозрачности (ЭИП) (см. [9. 10]) может быть реализован в такой системе при малой, вплоть до нескольких фотонов, интенсивности пробной волны, что позволяет реализовать квантовую память в системе и сделать вывод о перспективности при реализации квантовых вычислений. В отличие от работ, где рассматривается процесс прохождения поля через такую
\Ъ)
И
Л
сроду, в настоящей работе подробно исследуются особенности отклика среды на квантованное поле в наиболее общем, сжатом когерентном состоянии. Рассмотрено влияние квантово-статистичоских свойств света на среду и показана возможность снижения шумов параметров атомной системы, влияющих па процессы хранения и передачи квантовой информации. Задача решается аналитически методом Гейзенберга Ланжовона и численно квантовым методом Монте-Карло. Результаты численного моделирования сопоставляются с аналитическим расчетом методом: между ними наблюдается хорошее соответствие. Первый эксперимент со светом в сжатом вакуумном состоянии был проведен в 2004 г. [11]. Недавние экспериментальные работы подтверждают возрастающий интерес в этом направлении [12. 13].
1. Аналитическое рассмотрение методом Гейзенберга — Ланжевена
Рассмотрим трехуровневый атом (рис. 1), на переход |а) — |е) которого действует сильное классическое лазерное поле Ep с частот ой шр. На переход |а) — |b) действует слабое квантованное поле (сигнальное) Es = e(iMs)V(t)e-l^st + э.с., где а оператор уничтожения, подчиняющийся бозонным коммутационным соотношениям, е(ш ) = л/2ттПш/V, V - объем квантования, для простоты выбранный как объем системы, Д = uab — ws - отстройка от резонанса. Гамильтониан взаимодействия V в резонансном приближении имеет вид
V = —h(gaabci(t)e-iu-t + B.c.) — h(aacO(t)e-i"pt + э.с.), (1)
где О - частота Раби поля накачки, g - константа связи, а^ - операторы атомных переходов. Видно, что в гамильтониан взаимодействия явно входит оператор поля
Es
операторов определяется уравнением Гейзенберга Ланжевена [2, 3]
да i
-|jr = -lijVij + [Н, + Fij, (2)
где Yij ~ скорости релаксации, Fj - дельта-коррелированные операторы шума Ланжевена с корреляционными функциями (Fj (t)F+ (t')) = DjS(t — t'). Первый член в правой части (2) отражает возможность спонтанного излучения в атоме. Наличие шумового члена является проявлением флуктуационно-диссипационной теоремы.
В первом порядке теории возмущений по малому параметру е = gV/О получаем из (2) отделившуюся систему уравнений
а ba = (—Y + «Д)аьа + igV + iüabc + Fba,
(3)
а bc = (—Г + iA)abc + iO*aba + Fbc-
Учтем квантовый характер сигнального поля. Если на вход подается пробное поле в когерентном состоянии, (V)COh = а, где а пропорциональна амплитуде поля,
то. переходя к усредненным величинам в (3). получаем:
= (-7 + ¿Д)(стЬа) + гда + Ш(стЬс)
(4)
(а ьс) = (-Г + гД)(стьс) + гП*' '
Стационарное решение для поляризации на переходе |а) — |Ь) запишется в виде
гда{Т + гА)
(7 + гД)(Г + гД) + |ПРГ И
Этот результат, полученный нами из полностью квантового рассмотрения, соответствует классическому анализу эффекта ЭИП (см. [9. 10]). Для неклассических состояний поля, например фоковского с определенным числом частиц |п) или сжатого вакуума [2], (У) = 0 и среднее значение поляризации па переходе (5) (аЬа) = 0. При этом для описания поведения системы нам недостаточно знать средние значения. необходимо также исследовать флуктуации отклонения рассматриваемых параметров от среднего. Для этого введем операторы ба^ = а^ — (ац), бУ = У— (а). Используя преобразование Фурье получаем из (3):
йаьаМ =-(7 + гД)(г + гд) + |прр-' (6)
где Д = Д—и, для итератора У выполняется нормировка: [У(£), У(£')] = Ь/еб^Ь—Ь'), Ь — длина квантования. Таким образом, мы можем определить, как изменяется не только среднее значение поляризации па переходе |а) — |Ь), но и ее флуктуация. Существует два основных вклада в |а) — |Ь): ЭМ-иоля, причем зависимость от параметров системы имеет вид, аналогичный (5), и шума добавляются члены, определяемые лаижевеиовскими операторами, связанными с диссипативными процессами в системе. При данном подходе мы но учитывали изменение статистики сигнального поля при прохождении через среду.
Для описания процессов релаксации удобно использовать операторы, аналогичные квадратурам X = (У+У+)/2, у = (У — а+)/2г для ЭМ-поля, а\ = (аЬа + а+)/2, а2 = (аЬа — а+а)/2^, и операторы аз, а4 для аЬс. Их дисперсию можно получить в общем виде для любого квантованного поля из (4) (6), пользуясь флуктуациоино-диссипациоппой теоремой [2, 3], например, для а2 имеем:
(5а1) = ^^((<5Х2)(|П2|Г + Д27 + 7Г2)2 + ((5У2)Д2(Г2 + Д2 - \П2\)2 +
+ (бХбУ + бУбХ )(|П2|Г + Д27 + 7Г2)Д(Г2 + Д2 — |П2|)) +
|П2|Г + Д27 + 7Г2 + 2|с1е1;|2 ' (П
где |det|2 = (|02| — Д2 + Г7)2 + Д2(Г + 7)2 . При выводе мы учитывали, что флуктуации операторов шума и поля не коррелированны.
Далее вычисления проводим для света в сжатом когерентном состоянии, для которого имеем [2]: (6Х2)8Ч2 = е-2г/4, (бУ2)8Ч2 = е2г/4, (бХбУ + бУбХ)Щ2 = = 0, считая параметр сжатия £ = твгв действительным (в = 0 для определенности). Дисперсия оператора X в сжатом когерентном состоянии может быть много меньше, чем в когерентном состоянии в зависимости от степени сжатия т: (бХ2)щ2 ^ (бХ2)СоЬ = 1/4- Оказывается, этим свойством света можно воспользоваться и уменьшить дисперсию параметров атомной системы (2), например,
Рис. 2. Дисперсия оператора о\ в зависимости от нормированной расстройки Д/7. Пробное поле в когерентном состоянии г = 0, (пунктирная линия), в сжатом состоянии г = = 0.5 (штриховая линия) иг = 1.5 (сплошная линия). Крестиком отмечены результаты численного моделирования (подробнее см. разд. 2). Параметры Г = 10-37, П = 2/З7, д2Ь/а1 = 1
для {Sal) имеем с учетом нормировки:
(6*1) = ^^ (е2г(|П2|Г + 7(А2 + Г2))2 + е~2г Д2(Г2 + Д2 - |П2|)2) +
|^2|Г + 7(Д2 + Г2) + 2jdetp ' (8)
Рассмотрение данного вопроса стимулировано результатами анализа аналогичных систем в резонаторе [5. 6]. где было показано, что существует возможность полного переноса сжатия на флуктуирующие параметры среды эффект спинового сжатия (spin squeezing). В таких процессах принципиальную роль играют параметры (7) и (8). полученные нами аналитически, несмотря на то что напрямую они не измеримы. Очевидно при этом, что изменение состояния света при прохождении через среду связано именно с взаимодействием со средой и обратным воздействием среды на поле. В разд. 3 этот вопрос будет рассмотрен подробно.
2. Численное моделирование квантовым методом Монте-Карло
С помощью унитарного преобразования U = e(iUs°ab+iUp°ac)t избавимся от явной зависимости от времени в гамильтониане: (1)
H = HAacc + hAaaa - Ъд(ааъа + э.с.) - h(aacS! + э.с.). (9)
Уравнение для оператора плотности в форме Линдблада [2] примет вид
Р= -i/h[H,p]+ LabP + LacP + LcbP- (10)
Операторы Lj действуют следующим образом:
LabP = Yab/2(2abapaab - a ааР - pa аа),
LacP = Yac/2(2acapaac - a ааР - pa аа), LcbP = r(accpabb + abbpacc)-
Для дальнейшего решения методом Монте-Карло уравнение (10) должно быть переписано [4]: р = -i/h[Heñ-p — pH+ ] JjР-, где введены Иев - эффективный неэрмитов гамильтониан, Jj - оператор скачка, которые имеют вид
Heff = H — ih/2(lab + Yac — Г/2 )<7аа,
Jj Р = Yj Cj pj, j = 1 2, 3, ci = aba, C2 = ffca, C3 = acc — ^bb,
Yi = Yab, Y2 = Yac, Y3 = Г/2.
Предполагается, что над исследуемой системой проводятся измерения, позволяющие учесть процессы спонтанного излучения (например, с помощью фотодетектора ). Изменение наших знаний о поле излучения при регистрации фотона либо отсутствии фотоотсчета непосредственно приводит к изменению волновой функции системы. Тогда динамика квантовой системы представляется как непрерывное изменение векторов состояния в отсутствие фотоотсчетов под действием неэрмн-тового гамильтониана Heg:
±т)) = -i/mcS т.)), (и)
которое прерывается скачками, переводящими волновую функцию в определенное состояние при регистрации фотоотсчета
| ф(г + Ai)) = л Cj | Ф(г)), (12)
где Л - нормировочная константа, At - шаг по времени, через который реги-
cj
ci c2
ма с уровня |а), переводя волновую функцию системы в состояния |b) и |c) при скачке соответственно. Оператор c3 отвечает за дефазировку на низкочастотном переходе |c) — |b). Вероятности скачков Pj иод действием операторов Cj зависят от времени и могут быть вычислены [4]:
Pj (t)= Yj Atm)lc+cj | ф). (13)
Одномодовое сжатое когерентное состояние может быть получено действием оператора смещения и оператора сжатия па вакуумное состояние |0) [2, 3]:
|в,С) = exp(ea+ + ¡3*a) exp[(£*a2 — С«+2)/2]|0), (14)
где С характеризует сжатие, в ~ когерентное состояние. Состояние |в, С) может быть разложено по фоковским состояниям с определенным числом фотонов |n)
у. (e^thr)"^ —(I/312 —e~¿e/32t.h r)/2 jj ( íi^'2 \ , n) (15)
где Hn - полиномы Эрмита.
При достаточно сильном сжатии распределение числа фотонов p(n) является узким и быстро спадающим [3], поэтому оказывается достаточным использовать конечное число членов разложения (15) для получения требуемой точности. Волновую функцию системы «атом + иоле» разложим по базису |ф(t)) = i n c¿,n(t)|i, n), где i = a, b, c - атомный уровень, n = 0,1, 2,... - число фотонов в моде поля. Мы получим систему уравнений для амплитуд вероятности ci;j (t), описывающую динамику атома при отсутствии скачков. Вероятности квантовых скачков определяются
в соответствии с (13). Мы исследуем среднее значение и дисперсию атомных операторов <jjj , характеризующих параметры атома в данной задаче, в стационарном состоянии, вычисленные в предыдущем разделе методом Гейзснберга Лаижевеиа аналитически. На основе метода Монте-Карло они могут быть просто вычислены и, например, для ai, имеют вид
(ох) = (VKOki = +К.С, (16)
n
{o-D = №)\*i№)) = + K»l2> <^2> = - b>2- (17)
n
Для проверки теоретических результатов мы рассмотрели поле в сжатом когерентном состоянии со средним числом фотонов (и) = 16 и степенью сжатия r = 1.5 (рис. 2). Такие степени сжатия легко получают в эксперименте [12, 13]. Как показали наши численные исследования, в данном случае достаточно ограничиться числом фотонов и = 100 в разложении (15), размерность системы есть N = 3 х 100 = 300 . Для соблюдения точности расчетов, дающей статистическую ошибку расчета среднего и среднеквадратического отклонения (öai) < 0.01, число реализаций метода выбрано равным 100. t = 10y-1 есть время, при котором система переходит в стационарное состояние, шаг по времени равен 0.05y-1 . Сравнение результатов, полученных методом Монте-Карло и Гейзснберга Ланжевена представлено на рис. 2 для перехода |а) — |b) и степеней сжатия r = 0.5 и r = 1.5. Видно, что результаты численного расчета и аналитического решения находятся в хорошем соответствии. Преимуществом численного решения является возможность учета статистики фотонов квантованного света при взаимодействии со средой, исследования поля произвольной статистики, что практически невозможно сделать методом Гейзснберга Ланжевена аналитически и требует большого времени при численном расчете элементов матрицы плотности. Важно отметить, что использование метода Монте-Карло позволило нам учесть квантовый шум, возникающий при диссипативном взаимодействии системы с резервуаром, что принципиально для правильного решения задачи. Кроме того, этим методом удобно было бы исследовать изменение статистики фотонов при распространении поля в среде либо в стационарном режиме при взаимодействии со средой в резонаторе.
3. Распространение пробной волны в режиме ЭИП
Рассмотрим процесс прохождение сигнального поля через среду в режиме ЭИП. Пусть L - длина среды. От операторов для отдельного атома перейдем к операторам, непрерывно зависящим от переменной z путем усреднения
z zkeNz z zkeNz
по интервалу Az , содержащему число частиц Nz ^ 1 , так что на Az медленно меняющуюся амплитуду сигнального поля можно считать постоянной. Здесь aj -оператор атомного перехода для отдельного атома. Тогда гамильтониан взаимо-
N
V = ~J ^N(z)(h(g<7abEa(z,t)+H.c.) + h((7acn(z,t)e^t-^k''z+3.c.), (18)
где от суммирования по числу частиц мы перешли к интегрированию. В качестве сигнальной волны будем рассматривать широкополосный сжатый вакуум [2].
Для описания свойств поля используем квадратурный оператор
X(г, ш) = (а(г, ш8 + ш) + а+(г, ш8 - ш))/2,
аналогичный рассмотренному в разд. 1. спектр флуктуаций с учетом нормировки имеет вид
Б(ш) = е/Ь ! <ш' {5Х(ш)5Х(ш')).
В представлении Гейзенберга укороченное уравнение для медленно меняющейся амплитуды сигнального поля в среде запишется в виде
д д \
— +с—\а(г,$={дМаьа(з,$. (19)
Применим преобразование Фурье и перейдем к уравнению для флуктуаций параметров системы в (19):
д г-/ ч ъш ¡дМ
— оа(г,и>) = —оа(г,и>)-\--осгъа[2,ш). (20)
дг е е
Подставим сюда ЗаЬа(г,ш) из (6), тогда получим:
д
— 6а(г, ш) = к(ш)6а(г, ш) + ш), дг
где
ш д2М (Г-Н А)
к(ш) =--
Р (г,ш) =
с с (7 + «Д)(Г + гД) + |П|2 ' гдЫ (Г + (К)Ръа{г, и) + из)
с (7 + гД)(Г + гД) +|П|2
Ограничимся случаем точного резонанса, когда юа,ъ — и>е=0. Это оправдано, так как в Д расстройка и сдвиг частоты, связанный с немонохроматичностыо. дают одинаковый вклад. Для спектра флуктуаций получим:
=2Яек(си)8(г,си) + ^О(си), (21)
где шумовой член имеет вид
д2М2 |П2|Г + ш27 + 7 Г2
В(ш)
2е2 (7Г - ш2 + |П2|)2 + ш2(Г + 7)2'
Выражение для В(ш) получено с помощью флуктуационно-диссипационной теоремы в приближении идеального режима ЭИП и слабой пробной волны, когда основная часть атомов находится в состоянии |Ь). Решение (21) в данном приближении можно легко найти
Б(Ь,и;)=3(0,и;)Т-^(Т-1), (22)
где Т = е2ЬПек(ш) - коэффициент прохождения.
Х
выходе из среды от параметров поля на входе и свойств среды. Из него видно, что принципиальным является учет флуктуаций в системе (второй член в (21)).
Рис. 3. Коэффициент прохождения Т(ш) (слева) и спектр сжатия £(ш) в децибелах (справа), в зависимости от частоты ш (МГц) для различных скоростей релаксации в среде; Г = 10- 7 - сплошная линия, Г = 10- 7 - пунктирная линия, Г = 10- 7 - штриховая линия. Штрих-пунктирная линия - уровень сжатия на входе; 7 = 6 МГц, О = 27, д2ЬЫ/с7 = 102 (см. [11])
Их учет позволяет не только сохранять коммутационные соотношения между операторами, но и получить правильное решение (22). Видно, что при ш = 0, Г = 0, то есть при идеальных условиях ЭИП, Ие к(ш) = 0 и спектр па выходе имеет вид Я(Ь, ш) = Я(0, ш), то есть нет зависимости от флуктуации. Это согласуется с представлением о том, что при идеальном режиме верхний уровень совсем не заселен, {&аа) = 0. Тогда, так как в системе отсутствует диссипация, пет и флуктуаций, и спектр сигнала не модифицируется. Если на вход подается вакуум Я(0, ш) = 1/4, то па выходе — спектр вакуумных флуктуаций: Я(Ь,ш) = 1/4. Уровень сжатия на выходе из среды как отношение уровня флуктуаций сжатого света к вакууму п(ш) = Я(Ь,ш)/Ятас(ш) представлен на рис. 3. Прохождение сжатого вакуума через среду в режиме ЭИП было недавно исследован экспериментально [12, 13]. Сжатый вакуум в качестве пробной волны подавался на вход среды в присутствии волны накачки с помощью балансного гомодинного детектора и спектроаналпзато-ра измерялся уровень сжатия на выходе. Данные, полученные по аналитическим зависимостям (21), (22), показаны на рис. 3 для параметров, аналогичных параметрам работы [11], и качественно совпадают с результатами эксперимента [11]. Таким образом, аналитически и численно найденные в разд. 1, 2 параметры среды могут быть проверены экспериментально.
Заключение
В работе рассмотрено взаимодействие квантованного ЭМ-поля в сжатом состоянии со средой Л-атомов в режиме ЭИП. Нами исследованы воздействие поля па среду и возможность передачи квантовых параметров поля, таких, как сжатие, атому. Поскольку квантовые параметры атома, отвечающие за взаимодействие с квантованным светом, измерить напрямую невозможно, мы использовали метод Монте-Карло прямого численного моделирования, позволяющий определить необходимые характеристики, учесть квантовый шум произвольной статистики и перспективный для дальнейших исследований. Дополнительно нами использовался аналитический метод Гейзенберга Лаижевена. Мы рассматривали свет в наиболее общем, сжатом когерентном состоянии. Использование данного состояния позволило получить решение аналитически и подтвердить результаты численного моделирования. Нами найдены аналитические зависимости дисперсии операторов атомных переходов, характеризующих параметры атома, от статистических
параметров света при взаимодействии. Показано, что квантовые свойства света, такие, как сжатие, могут передаваться атому в таком режиме. Рассмотрен способ экспериментальной проверки предложенных зависимостей.
Авторы выражают благодарность АЛ!. Сатанину и М.В. Федорову за полезные обсуждения.
Работа выполнена при поддержке РФФИ (проект Л*1' 09-02-97059).
Summary
A.J. Gelman, V.A. Mironov. Peculiarities of Dissipative Dynamics of Л-atom Driven by a Quantized Field in Squeezed State.
The dissipative dynamics of the three-level Л-atom driven by a quantized electromagnetic field under the electromagnet.ically induced transparency (EIT) conditions is investigated both analytically by the Heisenberg Langevin method and numerically using the quantum Monte Carlo wave function method (QMC method). The impact of the squeezing level of quantized light 011 the atomic system parameters, which determine the relaxation processes and noise in EIT transparency window, is considered in detail. It is shown that under these conditions quantum properties of light (squeezing) can be transferred to an atom. The work also considers the problem of the propagation of quantum field in squeezed coherent state in the medium of atoms and analyzes the dependence of squeezing level of light at the medium output 011 the parameters at the input and properties of the medium. The application of the QMC method to numerical simulation of the propagation problem is discussed in the light of the recent experiments 011 the propagation of squeezed vacuum in the medium under examination.
Key words: squeezed light, quantum jumps, quantum noise, quantum memory.
Литература
1. Нильсен M., Чанг И. Квантовые вычисления и квантовая информация. М.: Мир, 2006. 211 с.
2. Gardiner С. W., Zoiler P. Quantum noise. Berlin: Springer, 2000. 438 р.
3. Скалли М.О., Зубайри M.G. Квантовая оптика. М.: Физматлит, 2003. 512 с.
4. Pleniu М.В., Knight P.L. The quantum-jump approach to dissipative dynamics in quantum optics // Rev. Mod. Pliys. 1998. V. 70, No 1. P. 101 143.
5. Danian A., Pinard M. Quantum-state transfer between fields and atoms in elect.romag-netically induced transparency // Pliys. Rev. A. 2004. V. 69. P. 043810.
6. Vernac L., Pinard M., Giacubinu E. Spin squeezing in two-level systems // Pliys. Rev. A. 2000. V. 62, No 6. P. 063812-1 063812-9.
7. Mazets I.E. Adiabatic propagation of quantized light pulses in an atomic medium with the tripod level configuration // ЖЭТФ. 2006. T. 130, № 3. C. 421 425.
8. Гельман А.И., Миронов В.А. Подавление шума в атомной системе под действием поля в сжатом когерентном состоянии // ЖЭТФ. 2010. Т. 137, Л' 4. С. 627 636.
9. Lukin M.D. // Colloquium: Trapping and manipulating photon states in atomic ensembles // Rev. Mod. Pliys. 2003. V. 75. P. 457 472.
10. Harris S.E. Electromagnetically Induced Transparency // Pliys. Today. 1997. V. 50, No 7. P. 36 42.
11. Akamatsu D., Akiba К., Kozuma М. Electromagnetically Induced Transparency with Squeezed Vacuum // Pliys. Rev. Lett. 2004. V. 92, No 20. P. 203602-1 203602-4.
12. Honda K., Akamatsu D., Arikawa M. et al. Storage and Retrieval of a Squeezed Vacuum // Phys. Rev. Lett. 2008. V. 100, No 9. P. 093601-1 093601-4.
13. Appel J., Figueroa E., Kurystuv D. et al. Quantum Memory for Squeezed Light // Phys. Rev. Lett. 2008. V. 100, No 9. P. 093602-1 093602-4.
Поступила в редакцию 09.12.09
Гельман Александр Иосифович младший научный сотрудник Института прикладной физики РАН, г. Нижний Новгород. Е-шаП: деЛтап варр1. иы-ппоь. ги
Миронов Вячеслав Александрович кандидат физико-математических паук, старший научный сотрудник Института прикладной физики РАН, г. Нижний Новгород. Е-шаП: fravinanQappl.sci-nnov.ru
