Динамика двухуровневых атомов, взаимодействующих с тепловым полем, для перепутанных начальных атомных состояний

УДК 130.145

Вестник СамГУ. 2015. № 6

ФИЗИКА

Е.К. Башкиров, А.С. Соловьева, М.С. Мастюгин1

ДИНАМИКА ДВУХУРОВНЕВЫХ АТОМОВ, ВЗАИМОДЕЙСТВУЮЩИХ С ТЕПЛОВЫМ ПОЛЕМ, ДЛЯ ПЕРЕПУТАННЫХ НАЧАЛЬНЫХ АТОМНЫХ

СОСТОЯНИЙ2

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

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

Введение

Возможность перепутывания квантовых состояний представляет собой одну из наиболее ярких черт квантовой механики. Перепутанные состояния играют ключевую роль в физике квантовых вычислений [1]. В контексте квантовой электродинамики резонаторов были предложены многочисленные схемы генерация перепутывания состояний двух кубитов за счет взаимодействия с общим полем резонатора [2; 3]. Тепловые поля обычно рассматриваются в качестве состояний, которые содержат минимум информации. Взаимодействие квантовых систем с такими полями обычно приводит к декогеренции. Однако Ким и соавторы [4] недавно показали, что такие некоррелированные поля могут приводить к перепутыва-

-1© Башкиров Е.К., Соловьева А.О., Мастюгин М.С., 2015

Башкиров Евгений Константинович (bash@samsu.ru), Соловьева Анастасия Сергеевна (solovievaas@mail.ru), Мастюгин Михаил Сергеевич (mastyugin.mikhail@mail.ru), кафедра общей и теоретической физики, Самарский государственный университет, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.

2Работа выполнена при финансовой поддержке Министерства образования и науки РФ (государственное задание № 1394).

нию состояний атомов в резонаторе. Они исследовали динамику двух идентичных двухуровневых атомов, резонансно взаимодействующих с модой теплового поля в идеальном резонаторе, и показали возможность генерации атомных перепутанных состояний как для чистых, так и для смешанных состояний атомов. При этом было показано, что степень перепутывания атомов существенно зависит от их начального состояния. В работах [5-8] было выявлено, что тепловое поле может также индуцировать высокую степень перепутывания как для случая неидентичных атомов, так и для многофотонных переходов в атомах. Агуиар с соавторами [9] обобщили результаты работы [4], приняв во внимание прямое диполь-дипольное взаимодействие атомов. При этом было показано, что дипольное взаимодействие между атомами может привести к существенному возрастанию максимальной степени атомного перепутывания. Недавно в ряде работ [10-12] было показано, что степень атомного перепутывания двухкубитных систем, индуцированного тепловым полем, может существенно возрастать за счет наличия начальной атомной когерентности. При этом авторы уделили основное внимание исследованию перепу-тывания кубитов для начальных неперепутанных атомных состояний. Представляет интерес рассмотреть динамику атомного перепутывания, возникающего за счет взаимодействия с тепловым шумом, в случае перепутанных начальных состояний атомов и с учетом их прямого диполь-дипольного взаимодействия. Для исследования атомного перепутывания в указанной модели мы найдем точную динамику, используя подход, отличный от развитого в работах [9; 10] и базирующийся на применении полного набора "одетых состояний" исходного гамильтониана.

1. Модель и решение уравнения эволюции

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

2

н =

¿=1

где а+ = \+)¿¿(—\^а- = —^¿(+\ - операторы переходов для возбужденного и основного состояний г-го двухуровневого атома (¿=1,2), а+ и а - опера-

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

Будем полагать, что в начальный момент времени атомы находятся в чистом состоянии, а поле - в одномодовом тепловом состоянии. В настоящей работе мы ограничимся рассмотрением атом-атомных перепутанных состояний. Для описания динамики перепутывания атомов необходимо вычислить редуцированную атомную матрицу плотности рагот(Ь). А. Перес и Хородецкие доказали [13; 14], что необходимым и достаточным условием неперепутанности состояний двух ку-битов является то, что некая дополнительная матрица рТ1от, полученная путем частичной перестановки индексов (транспонирование по переменным одного ку-бита) матрицы плотности рагот(Ь), имеет только неотрицательные собственные

значения. Тогда для определения степени атом-атомного перепутывания можно воспользоваться параметром Переса — Хородецких, который определим как

е = -2^2

где ¡ - отрицательные собственные значения транспонированной по переменным одного кубита атомной матрицы плотности pT1om. Отсутствие у частично транспонированной по переменным одного кубита атомной матрицы плотности (2) отрицательных собственных значений означает, как уже было сказано выше, что атомы в любой момент времени остаются неперепутанными, в этом случае полагают е = 0. Если параметр перепутывания е > 0, то атомы находятся в перепутанном состоянии. Максимальная степень перепутывания атомов соответствует значению параметра перепутывания е = 1.

Рассмотрим теперь динамику параметра перепутывания (2) для различных чистых перепутанных состояний атомной подсистемы белловского типа

или

где cos © и sin © - амплитуды и ф - относительная фаза состояний, входящих в правые части формул (3) и (4), а также одномодового теплового состояния резо-наторного поля, описываемого матрицей плотности

PF

n

Здесь Ш - одномодовое фоковское состояние поля и вероятности pn задаются для теплового поля формулой pn = nn/ (nn + 1)n+ , где n - среднее число фотонов в моде. Для этой цели предварительно найдем решение временного уравнения Шредингера для модели с гамильтонианом (1) в случае чистого фоковского состояния поля и различных начальных чистых состояний кубитов. Для начальной волновой функции полной системы |Ф(0)) указанное решение можно представить в виде

i

где |Ф(0)) = Y1 СД0)|Ф^ и |Фi) и Ei - собственные функции и собственные зна-

i

чения гамильтониана (1) соответственно.

В двухатомном базисе —, — ,n+2), |+, — ,n+1), —, +,n+1), |+, +, n) собственные функции (так называемые "одетые состояния") гамильтониана (1) есть

|*in) = tUXi^—, —,n + 2) + Xi2n|+, —,n + 1) + +Xi3n—, +,n + 1) + Xnn^, +,n))

где

- константы нормировки для собственных функций и

Xiin = 0, Xi2n = —1, X13 = 1, X14 = 0,

X21n = — :, X22n = 0, X23n = 0, X24n = 1,

a/2 + n

2 + n а + Q

Л31и — ы , X32n

A/rrwr+n' 4V/1m'

^ _1(1+ п)(а + 0„) + (а + Q„)(2 + n + а(а + 0„)/2

X33n — ö-и !r.-;-ö--), X34 — 1

2 л/1 + n (6 + 4n + а2 + aQn

„ _ 2 + n v _v _ а

X41n — ,, , X42n — X43n — -... ,, TZ.---4--,

л/1 + W2 + n V1 + n

X44 — 1,

где Qn — V24 + 16n + а2.

Соответствующие собственные значения имеют вид

Ein — 0, E2n — а, E3n — —(1/2)(а + Qn), E4n — -(1/2)(а - Qn), где а — J/g.

В двухатомном базисе \ —, —, 1), |+, —, 0), \ —, +, 0} собственные функции и собственные значения гамильтониана (1) можно представить в виде

|*1> — ^[\+, —,0} — \—, +,0}], Ei/h — —а;

4

\^2> —

\^2>

8 + (ш + а)2 4

2} —

— + а) — , —, 1} + \+, —, 0} + \ —, +, 0} —1(ш — а)\—, —, 1} + \+, —, 0} + \ —, +, 0}

E2/h — — — а); E3/h — !(ш + а);

8+ (ш — а)2 где ш — а/8 + а2.

Наконец, гамильтониан (1) имеет еще одну собственную функцию вида

\Фо} — \—, —, 0},

соответствующую энергии Eo — 0.

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

\—, — ,n + 2}, \+, — ,n + 1}, \—, +, n + 1}, \+, +,n}.

Предположим, что рассматриваемая система в начальный момент времени находится в чистом состоянии \+, —,n + 1}. В этом случае с использованием соотношений (6)-(9) временную волновую функцию можно записать в виде

\Ф(*)} — Zi2n\ — , —,n + 2} + Z22,n\ + , —,n +1} + Z32,n\ — , +,n +1} + Z42,n\ + , +,n}, (10) где

Zi2,n — e-lElnt/% ein Y21n Xiin + e-lE2nt/% Y22nX21n +

+e-1E3nt/h ы Y23n X31n + e-lE4nt/R e4n Y24n Х4Ш, Z22,n — e-lElnt/n £1n Y21n X12n + e-lE2nt/K^ bn Y22n X22n +

+e-lEsnt/h £3n Y23n X32n + e-lE4nt/R ^4n Y24n X42n, Z32,n — e-lElnt/h £1n Y21n X13n + e-lE2nt/n Y22n X23n + +e-lEsnt/h £3n Y23n X33n + e-lE4nt/R ^4n Y24n X43n,

242,п — е-гВ1п*1Л а„ У21п ХХАп + е-1^1п У22п Х24п+ +е-гЕ3п1/П ы У23п Хз4п + в-гЕ4п1/П ^п У*4п Х44П

и У1]п — £,jnXjin

Если же начальное состояние полной системы выбрать в виде |+, —, 0), то временная функция примет вид

где

£12 — — (2/ ш) е-га1/2 sin(wt/2), 222 — (1/2) \егаг + с-гаг/2 [сов(ш^2) — (%а/ш)$т(шг/2)]} ,

г32 — (1/2) {—е1а* + в-1а*/2 [cos(шt/2) — (га/ш)вт(шг/2)]} .

Если мы выберем в начальный момент времени волновую функцию в виде I —, +,п + 1), то ее временная эволюция будет описываться вектором состояния

Выражения для величин ^ можно получить из формул для , если в

заменить величины У2п на У3т, где г — 1, 2, 3,4.

Для начального состояния —, +, 0) временная функция имеет вид

|Ф(*)) — 21з|—, —, 1) + 223+, —, 0) + 2зз—, +, 0),

где 213 — , 223 — 2з2, 2зз — 222

Для начального состояния |+, +, п) временная волновая функция есть

Выражения для величин 2ц ,п можно получить из формул для ,п , если в (10) заменить величины У2т на Уцп, где г — 1, 2, 3, 4.

Для начального состояния —, —,п + 2) в последующие моменты времени волновая функция может быть представлена как

№

Выражения для величин 2.14,п можно получить из формул для , если в

заменить величины У2п на У4п где г — 1, 2, 3,4.

Для начального состояния |—, —, 1) временная функция имеет вид

где

214 — (1/ш) е-га1/2 [ш сов(шг/2) + га вт(шг/2)], 224 — 234 — — (2г/ш)е-га1/2 sin(шt/2).

Наконец, начальное состояние —, —, 0) не эволюционирует с течением времени. Теперь вернемся к решению поставленной задачи. Пусть атомы приготовлены в одном из перепутанных состояний вида (3) или (4), а поле - в одномодовом состоянии (5). Тогда используя соотношения (6)-(15), мы можем найти явный вид решения квантового уравнения Лиувилля для временной матрицы плотности полной системы, включающей в себя два атома и моду поля. Проводя ее усреднение по переменным поля, мы можем получить редуцированную атомную матрицу плотности. В двухатомном базисе |+, +), |+, —), —, +), —, —) она имеет вид

рл(г) =

( Р11 Р12 Р13 Р14 ^

* Р12 Р22 Р23 Р24

* Р13 * Р23 Р33 Р34

* ^ Р14 * Р24 * Р34 Р44 у

Элементы матрицы (16) получены нами в явном виде, однако не приведены в работе ввиду чрезмерно громоздкого вида. Соответствующая (16) частично транспонированная по переменным одного кубита атомная матрица плотности рТ есть

( Р11 Р12 р!з р2З \

Р12 Р22 Р14 Р24

Р13 Р14 Р33 Р34

V Р23 Р24 Р34 Р44 )

Вычислив ее собственные значения, мы можем определить параметр Переса — Хородецих (2), позволяющий оценить степень перепутывания состояний атомов.

рТ V

2. Основные результаты

Результаты численного моделирования параметра атом-атомного перепутыва-ния (2) для различных значений параметров модели и начальных атомных состояний вида (3) и (4) приведены на рис. 1-3.

Рис. 1. Зависимость параметра перепутывания (2) от приведенного времени дЬ для

начального атомного состояния вида

Из рисунков хорошо видно, что при определенном выборе параметров начального атомного перепутанного состояния вида (3) или (4) возможно сохранение высокой степени перепутывания атомов в процессе их дальнейшей эволюции даже в случае высоких интенсивностей теплового поля в резонаторе. При этом зависимость степени перепутывания от интенсивности диполь-дипольного взаимодействия является весьма слабой, что представляется существенным преимуществом рассматриваемой схемы по сравнению со схемой, использующей неперепутанные начальные состояния атомов [11]. Это связано с тем, что для того чтобы достичь

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

fit) 1.0*

0.8 0.6 0.4

0.2

0.0.

fn(t) 1.0t

10

15

20

0.8 0.6 0.4 0.2 qt 0.0,

15

20

5 10

а б

Рис. 2. Зависимость параметра перепутывания (2) от приведенного времени дЬ для начального атомного состояния вида (3). Относительная фаза ф = 9/10п. Параметр диполь-дипольного взаимодействия а =1. Среднее число фотонов в моде п = 1 (а) и п = 10 (б). Параметр в равен: 11/24п (точечная линия), 5/12п (штриховая линия'

и п/4

e(t) 0.010

20

40 60

б

80

10(?

Рис. 3. Зависимость параметра перепутывания

a

Выводы

Таким образом, в настоящей статье нами показана возможность сохранения высокой степени перепутывания кубитов для модели двух атомов, взаимодействующих с тепловым шумом, даже при высоких интенсивностях теплового поля. При этом имеется возможность управления и контроля за степенью перепутывания

за счет изменения параметров модели. При описании динамики перепутывания кубитов важен также учет диссипативных процессов, возникающих за счет взаимодействия кубитов с окружением. Такие взаимодействия приводят к исчезновению квантовой перепутанности кубитов. Учет влияния окружения на динамику перепутывания кубитов, оценка времени декогеренции и рассмотрение других механизмов контроля перепутывания в рамках рассмотренной в настоящей статье модели на основе управляющих уравнений, выведенных ранее в работах [15-19], будет предметом нашего следующего исследования.

Литература

[1] Schumacker D., Westmoreland M. D. Quantum Processes, Systems and Information. New York: Cambridge University Press, 2010. 469 p.

[2] Buluta I., Ashhab F., Nori F. Natural and artificial atoms for quantum computation // Reports on Progress in Physics. 2011. V. 74. P. 104401.

[3] Guhne O., Toth G. Entanglement detection // Physics Reports. 2014. V. 474. P. 1-75.

[4] Kim M.S., Lee J., Ahn D., Knight P.L. Entanglement induced by a single-mode heat environment // Physical Review. 2002. V. A65. P. 040101.

[5] Башкиров Е.К. Перепутанные состояния в системе двух неидентичных атомов, взаимодействующих с тепловым полем // Вестник Самарского государственного университета. Естественнонаучная серия. 2006. № 3. С. 21-29.

[6] Zhou L., Song H.S. Entanglement induced by a single-mode thermal field and criteria for entanglement // Journal of Optics. 2002. V. B4. P. 425-429.

[7] Башкиров Е.К., Ступацкая М.П. Перепутывание двух атомов, взаимодействующих с тепловым электромагнитным полем // Компьютерная оптика. 2011. Т. 35. № 2. P. 243-249.

[8] Башкиров Е.К., Мастюгин М.С. Перепутывание двух сверхпроводящих кубитов, взаимодействующих с двухмодовым тепловым полем // Компьютерная оптика. 2013. Т. 37. № 3. C. 278-285.

[9] Bashkirov E.K., Mastyugin M.S. The dynamics of entanglement in two-atom Tavis-Cummings model with non-degenerate two-photon transitions for four-qubits initial atom-field entangled states // Optics Communications. 2014. V. 313. P. 170-174.

[10] The entanglement of two dipole-dipole coupled atoms in a cavity interacting with a thermal field / L.S. Aguiar [et al.] // Journal of Optics. 2005. V. B7. P. S769-S771.

[11] Hu Y.H., Fang M.F. Coherence-Enhanced Entanglement Induced by a Two-Mode Thermal Field // Communications in Theoretical Physics. 2010. V. 54. P. 421-426.

[12] Башкиров Е.К., Мастюгин М.С. Влияние диполь-дипольного взаимодействия и атомной когерентности на перепутывание двух атомов с вырожденными двухфо-тонными переходами // Оптика и спектроскопия. 2014. Т. 114. № 4. С. 678-683.

[13] Bashkirov E.K., Mastyugin M.S. The influence of atomic coherence and dipole-dipole interaction on entanglement of two qubits with nondegenerate two-photon transitions // Pramana. 2015. V. 84. P. 127-135.

[14] Peres A. Separability Criterion for Density Matrices // Physical Review Letters. 1996. V. 77. P. 1413-1415.

[15] Horodecki M., Horodecki P., Horodecki R. Separability of Mixed States: Necessary and Sufficient Conditions // Physical Letters. 1996. V. A223. P. 1-8.

[16] Bashkirov E.K., Le Kien F., Shumovsky A.S. Collective spontaneous radiation of two atoms in the finite-Q cavity // Physica A: Statistical Mechanics and its Applications. 1990. V. 167. № 3. P. 935-944.

[17] Башкиров Е.К., Мангулова Е.Г. Динамика двух двухуровневых атомов в неидеальном резонаторе // Известия Российской академии наук. Сер.: Физическая. 2000. Т. 64. № 10. С. 2075-2079.

[18] Башкиров Е.К. Спонтанное излучение двух трехуровневых атомов в неидеальном резонаторе // Известия Российской академии наук. Сер.: Физическая. 2004. Т. 68. № 9. С. 1292-1295.

[19] Bashkirov E.K. Dynamics of phonon mode in superradiance regime of laser cooling of crystals // Physics Letters A. 2005. V. 341. № 1-4. P. 345-351.

References

[1] Schumacker D., Westmoreland M.D. Quantum Processes, Systems, and Information. New York, Cambridge University Press, 2010, 469 p. [in English].

[2] Buluta I., Ashhab F., Nori F. Natural and artificial atoms for quantum computation. Reports on Progress in Physics, 2011, Vol. 74, p. 104401 [in English].

[3] Guhne O., Toth G. Entanglement detection. Physics Reports, 2014, Vol. 474, pp. 1-75 [in English].

[4] Kim M.S., Lee J., Ahn D., Knight P.L. Entanglement induced by a single-mode heat environment. Physical Review, 2002, Vol. A65, p. 040101 [in English].

[5] Bashkirov E.K. Entangled states in the system of two nonidentiacl atoms interacting with thermal field. Vestnik Samarskogo gosudarstvennogo universiteta. Estestvennonauchnaia seriia [Vestnik of Samara State University. Natutal Science Series], 2006, no. 3, pp. 21-29 [in Russian].

[6] Zhou L., Song H.S. Entanglement induced by a single-mode thermal field and criteria for entanglement. Journal of Optics, 2002, Vol. B4, pp. 425-429 [in English].

[7] Bashkirov E.K., Stupatskaya M.P. Entanglement of two atoms interacting with thermal field. Komp'iuternaia optika [Computer Optics], 2011, Vol. 35, no. 2, pp. 243-249 [in Russian].

[8] Bashkirov E.K., Mastyugin M.S. Entanglement of two superconducting qubits interacting with two-mode thermal field. Komp'iuternaia optika [Computer Optics], 2013, Vol. 37, no. 3, pp. 278-285 [in Russian].

[9] Bashkirov E.K., Mastyugin M.S. The dynamics of entanglement in two-atom Tavis-Cummings model with non-degenerate two-photon transitions for four-qubits initial atom-field entangled states. Optics Communications, 2014, Vol. 313, pp. 170-174 [in English].

[10] Aguiar L.S., Munhoa P.P., Vidiella-Barranco A., Roversi J.A. The entanglement of two dipole-dipole coupled atoms in a cavity interacting with a thermal field. Journal of Optics, 2005, Vol. B7, pp. S769-S771 [in English].

[11] Hu Y.H., Fang M.F. Coherence-Enhanced Entanglement Induced by a Two-Mode Thermal Field. Communications in Theoretical Physics, 2010, Vol. 54, pp. 421-426 [in English].

[12] Bashkirov E.K., Mastyugin M.S. Influence of dipole-dipole interaction and atomic coherence on entanglement of two atoms with degenerate two-photon transitions. Optika i spektroskopiia [Optics and Spectroscopy], 2014, Vol. 114, no. 4, pp. 678-683 [in Russian].

[13] Bashkirov E.K., Mastyugin M.S. The influence of atomic coherence and dipole-dipole interaction on entanglement of two qubits with nondegenerate two-photon transitions. Pramana, 2015, Vol. 84, pp. 127-135 [in English].

[14] Peres A. Separability Criterion for Density Matrices. Physical Review Letters, 1996, Vol. 77, pp. 1413-1415.

[15] Horodecki M., Horodecki P., Horodecki R. Separability of Mixed States: Necessary and Sufficient Conditions. Physical Letters, 1996, Vol. A223, pp. 1-8 [in English].

[16] Bashkirov E.K., Le Kien F., Shumovsky A.S. Collective spontaneous radiation of two atoms in the finite-Q cavity. Physica A: Statistical Mechanics and its Applications, 1990, Vol. 167, no 3, pp. 935-944 [in English].

[17] Bashkirov E.K., Mangulova E.G. Dynamics of two two-level atoms in finite-Q cavity. Izvestiia Rossiiskoi akademii nauk. Seriia fizicheskaia [Bulletin of the Russian Academy of Sciences: Physics], 2000, Vol. 64, no. 10, pp. 2075-2079 [in Russian].

[18] Bashkirov E.K. Spontaneous radiation of two three-level atoms in lossless cavity. Izvestiia Rossiiskoi akademii nauk. Seriia fizicheskaia [Bulletin of the Russian Academy of Sciences: Physics], 2004, Vol. 68, no. 9, pp. 1292-1295 [in Russian].

[19] Bashkirov E.K. Dynamics of phonon mode in superradiance regime of laser cooling of crystals. Physics Letters A., 2005, Vol. 341, no. 1-4, pp. 345-351 [in English].

E.K. Bashkirov, A.S. Solovieva, M.S. Mastyugin3

DYNAMICS OF TWO-LEVEL ATOMS INTERACTING WITH THERMAL FIELD FOR ENTANGLED INITIAL

ATOMIC STATES4

In this paper the entanglement dynamics for system of two two-level atoms interacting with a mode of thermal electromagnetic field in lossless cavity has been investigated. Initial atomic states are the entangled Bell type states. Using the full set of eigenvectors of the Hamiltonian of the considered model we have derived the exact solution for the density matrix for the whole system. On its basis the reduced atomic density matrix and Peres — Horodetcki parameter have been calculated. Calculating of entanglement parameter shows the possibility of high degree of entanglement even for large intensity of thermal field. Thus there is a possibility of maintenance and control over the degree of entanglement.

Key words: two-level atoms, thermal field, dipole-dipole interaction, lossless cavity, one-photon transitions, atomic coherence, entanglement, atomic density matrix.

Статья поступила в редакцию 9//V/2015. The article received 9//V/2015.

3Bashkirov Eugene Konstantinovich (bash@samsu.ru), Solovieva Anastasya Sergeevna (anssolovieva@mail.ru), Mastyugin Mikhail Sergeevich (mastyugin.mikhail@mail.ru), Department of General and Theoretical Physics, Samara State University, 1, Academician Pavlov Street, Samara, 443011, Russian Federation.

4The work is carried out with the financial support from the Ministry of Education and Science of the Russian Federation (state assignment № 1394).