第30卷 第3期 哈尔滨师范大学自然科学学报 NATURAL SCIENCES JOURNAL OF HARBIN NORMAL UNIVERSITY 耦合双原子与两个单模腔场的 光子数分布和自发辐射特性 杨金宝,吕树臣,隋 欣 (哈尔滨师范大学) 【摘要】研究由两个全同的二能级原子分别同时与两个耦合微腔构成的系 统,腔场处于弱相干状态,忽略了两个原子之间的相互作用.采用理论推导和数值 计算相结合的方法,研究系统的平均光子数和腔场自发辐射光谱的性质.讨论了腔 场的衰减率与平均光子数的关系,以及腔场的耦合强度对弱相干场耦合腔系统自 发辐射光谱的影响.研究结果表明:当腔场的衰减率增加时,在两个耦合腔场的非 相干泵浦率和衰减率共同作用下,腔场的平均光子数减少,并随衰减率的增加而产 生了一个真空拉比震荡.当C 腔场的本征频率 变大时,发射光谱S( )双波峰 整体向高频段飘移.另一方面,当C 腔场的耦合强度g 增加时,高频段和低频段 的发射峰分别向高频段和低频段移动. 【关键词】量子主方程;原子一腔耦合;平均光子数;自发辐射 量子理论的J—C模型,采用量子主方程¨ 的方 0 引言 近年来,在J—c模型 建立后,腔量子电 动力学(cQED)作为探索微观物理世界的重要工 法,研究了在旋波近似下,由两个全同的二能级 原子分别同时与两个耦合微腔构成的系统中,平 均光子数和腔场自发辐射光谱的相关特性,讨论 了腔场的衰减率与平均光子数的关系,以及腔场 具,在量子计算 j、量子通信 和量子点激 光器 等领域得到了广泛的应用,原子与腔场之 间的相互作用受到了人们的普遍关注,随着研究 的耦合强度对弱相干场耦合腔系统自发辐射光 谱的影响. 的深入,人们对二能级原子与相干场、单模场,多 模场等腔场相互作用的问题进行了大量的实验 和研究,发现了许多与腔场有关的非经典效应, 1 物理模型和理论计算 该文采用文献 的物理模型,研究两个独 立的全同二能级Rydberg原子( 和A )同时和 例如光子的反群聚效应和亚泊松分布 J、真空 拉比l9 、以及原子的“坍缩和回复”现 象 等.另一方面,不少学者对双原子与腔场相 互作用产生的发射光谱做了进一步的研究,获得 了显著的成果.例如宋同强、冯健等¨ 研究了非 等同的双原子与单模场相互作用时,两个原子发 射光谱的相关变化.他们的研究结果表明:在真 两个纠缠的高Q耦合微腔(C.和c )构成的系 统,如图1所示. 在旋波近似下,双原子一微腔耦合系统的 主方程 为 ddo = [日,p]+ (p) (1) 空场作用情况下,拉比震荡会产生4个“坍缩和 系统的哈密顿量可以写成 2 1 2 1 回复”的发射峰,在弱场作用下,发射峰变为5 个;在强场作用下,发射峰变为3个.该文应用全 收稿日期:2014—03—25 日=∑— y 0 。 +∑告击 + 第3期 耦合双原子与两个单模腔场的光子数分布和自发辐射特性 119 图1 两个二能级原子与两个单模腔场的耦合模型 2 1 ∑如i(alo' ̄ +口 )+g3(0 口:+a1口 ) I:l ‘_ (2) 式(2)中 为第i个单模腔场的本征频率, 为 第i个原子的跃迁频率,g 为第i个原子与单模腔 场之间的耦合强度,g 是两个腔场之间的耦合强 度,a 和a是腔模的产生、湮灭算符, +和 一是 原子的泡利赝算符.式(1)中超算符L(p)¨ 用 林德布拉德形式表示为 L(p)=k1(2口 p01一aI口 p—pa10 )+ k2(2a;pa2一a20 p-pa2024-)+k3(2a1p口 一ala1P —papa1)+k4(2a2p0 一口 02p—papa2) (3) 式(3)中k 、k:分别为两个腔场的非相干泵浦 率 ,k,、k 分别为两个腔场的衰减率.由方程 (3),可以得到以下微分方程组: ( ) Ig, fg】0 1 ) ) 『 l + 一 ● I ) l 留。 +七0 喀 2一 『g0:I I 0 J ( ) ( ’) l 0 ( )’) O (4) ll( )l)]py,+毛一 『gl fg3( 口I)J l=l 0 。七 I 1『2~七4J l×l((( ) q)]《口I)jI (5) 由方程组(4)可以研究系统的平均光子数与腔 场的非相干泵浦率和腔场的衰减率的变化关系. 依据Wiener—Khintchine理论‘】 ,在稳态时发 射光谱L20 的表达式可以成 r∞ .s( )=2Re J e“ ”(0 (t)口) dt (6) J 0 由方程(5)和(6)可以研究系统的自发辐射光 谱.该文重点研究第一个耦合腔场.在单光子关 联近似( 口>=一(0)和费米统计近似[ 一, +]+=1下,由式(4)和(5),可以来求解耦合系 统稳态时,平均光子数为 <口 a1> s l + 可]I,~ala ̄ 一 3 一 (8) 其中m= 1一 2,n=k1+k2一k3一 4. 联合式(5)和(6),可以求解双原子一腔耦 合系统中,由双原子辐射产生的发射光谱,由耦 合微腔场产生的发射光谱,以及整个系统产生的 发射光谱. 2 数值计算与结果分析 利用(4)式和(5)式,可以计算腔场处于稳 态时,耦合系统的平均光子数,结果如图2~图 .1 量≤g qd蠢 ∞ ∞ ∞ ∞ 鲫 ∞ ∞ ∞ 0 ∞ 帕 5所示. ÷_r1. 10 一 .厂 . ● ●_ _ -I ~ 。~~“ 。~~ 。~’ 图2 稳态时,耦合腔场的平均光子数(a;a > 随 腔场非相干泵浦率k。的变化曲线 取定C 腔场的衰减率为k =0.03,C 腔场 的非相干泵浦率和衰减率分别为 =0.02和 =0.04,两个腔场之间的耦合强度为g =10k . 图2中两个耦合腔场的本征频率之间的关系为 y 一y:=l0,由公式(7)、(8)可以得到在稳态 时,耦合腔场c 的平均光子数(0 0 ) 随腔场非 相干泵浦率k 的变化曲线.如图2所示,在两个 耦合腔场的非相干泵浦率和衰减率共同作用下, 耦合腔场c 的平均光子数(a;a ) 随腔场非相 干泵浦率k 的增加而产生了一个真空拉比震 荡 .当k <k 时,为二能级原子A 发射光子 的过程,腔场平均光子数(。 a ) 随非相干泵浦 率k 的增加而增加;当 :k。时,腔场的平均光 子数(a;a ) 出现了一个真空拉比震荡,为二能 级原子A 从下能级激发到上能级,并快速吸收 120 哈尔滨师范大学自然科学学报 2014年第3O卷 光子的过程,腔场内的平均光子数急速减少;当 后 > 。时,平均光子数(a?a ) 随腔场非相干泵 浦率 的增加而增加,为二能级原子A 发射光 子,腔场的平均光子数回复的过程. b 藿 薹 塞 图3 稳态时,在两个腔场的本征频率分别取不同 差值时,耦合腔场C。的平均光子数(。 0 > 随腔场 非相干泵浦率 的变化曲线,图为 一Y:=10、20、 30的情况 取定C,腔场的衰减率为 =0.03,C 腔场 的非相干泵浦率和衰减率分别为 :=0.02和 =0.04,两个腔场之间的耦合强度为g,=10k . 由公式(7)、(8)可以得到在稳态时,两个耦合腔 场的本征频率分别取不同差值,耦合腔场c 的 平均光子数( a ) 随腔场非相干泵浦率 的 变化曲线.如图3所示,耦合腔场c 的平均光子 数(。 0 ) 随腔场非相干泵浦率 的增加而产 生了一个真空拉比震荡,当两个耦合腔场C 和 c2之间的本征频率之差增加时,耦合腔场c,的 平均光子数(n a ) 峰值增加,真空拉比震荡加 I. 矗 D 宣 至 皇 呈 磊 0 图4 稳态时,耦合腔场C 的平均光子数<n a ) 随腔场衰减率b的变化曲线 取定C 腔场的非相干泵浦率 =0.01,C 腔场的非相干泵浦率和衰减率分别为 =0.02 和 =0.04,两个腔场之间的耦合强度为g,: 10k .图4中两个腔场的本征频率之间的关系为 一 ::10,由公式(7)、(8)可以得到在稳态 时,耦合腔场c 的平均光子数(0 0 ) 随腔场衰 减率 的变化曲线.如图4所示,在两个耦合腔 场的非相干泵浦率和衰减率共同作用下,耦合腔 场C 的平均光子数(口 。。) 随腔场衰减率 的 增加而产生了一个真空拉比震荡,震荡趋势与图 2中震荡趋势相反.当 < 时,为二能级原子 A 吸收光子的过程,腔场平均光子数(口 a ) 随 衰减率 。的增加而减少;当 ,= ,时,腔场的平 均光子数(。 0 ) 出现了一个真空拉比震荡,为 二能级原子A 从上能级跃迁至下能级,并快速 发射出光子的过程,腔场的平均光子数急速增 加;当 > 时,平均光子数(。 0 ) 继续随腔 场衰减率 的增加而减少,最后趋于稳定,所得 结论与图2相符. 1∞4 图5 在C,腔场中,不同本征频率Y,条件下,系统 发射光谱S(Ca))随二能级原子A 跃迁频率 的变化 曲线,图为y1为998,999、1(300、1001、1002时的情况 取定C 腔场的非相干泵浦率和衰减率分别 为 =0.01、 =0.03,C 腔场的非相干泵浦率 和衰减率分别为 =0.02、 =0.04,C 腔场的 本征频率为y:=1000 Hz,在C,腔场中二能级原 子A。和单模腔场的耦合强度为g。=0.1,两个腔 场之间的耦合强度为g =0.2.图5中c 腔场具 有不同的本征频率 时,由公式(5)和(6)可以 得到,系统发射光谱S( )随二能级原子A。跃迁 频率 的变化曲线. 如图5所示,在真空场作用情况下,拉比震 荡会产生2个“坍缩和回复”的发射峰.当c 腔 场本征频率 变大时,发射光谱S( )双波峰整 体向高频段飘移. 取定C 腔场的非相干泵浦率和衰减率分别 为 =0.01、 =0.03,C腔场的非相干泵浦率 和衰减率分别为 =0.02、 =0.04,两个腔场 的本征频率分别为 =y:=1000 Hz,腔场之间 的耦合强度为g =0.2.图6中C 腔场具有不同 的二能级原子A 与腔场的耦合强度g 时,由公 式(5)和(6)可以得到,系统发射光谱s(∞)随 第3期 耦合双原子与两个单模腔场的光子数分布和自发辐射特性 l21 图6 在C 腔场中,不同耦合强度g,条件下,系统 发射光谱S( )随二能级原子A 跃迁频率to 的变 换曲线,图为g.为0.1、0.2、0.3、0.4、0.5时的情况 二能级原子A 跃迁频率to 的变换曲线. 如图6所示,当c 腔场的耦合强度g 变大 时,系统发射光谱S( )低频段所对应的波峰向 低频段移动,高频段对应的波峰向高频段移动. 3 结束语 研究了两个二能级原子与两个单模腔场耦 合,具有不同的衰减率和耦合强度时,系统的平 均光子数和腔场自发辐射光谱的相关特性.通过 量子主方程的推导,得到一组关于平均光子数的 封闭的一阶非线性微分方程组,利用计算机软件 得到方程组的解析解,并绘制了关于系统平均光 子数和腔场自发辐射光谱的变换曲线.研究结果 表明:在旋波近似下,当腔场的衰减率增加时,在 两个耦合腔场的非相干泵浦率和衰减率共同作 用下,腔场的平均光子数减少,并随衰减率的增 加而产生了一个真空拉比震荡.腔场本征频率变 大,发射光谱S( )双波峰整体向高频段飘移.另 一方面,当腔场的耦合强度增加时,高频段和低 频段的发射峰分别向高频段和低频段移动. 参考文献 [1]Wang Z C,Shen FH,eta1.Sudden deathand sudden birth of atomic entanglement in doulde Jaynes—Cummings model composed of two coupled cavities.Chinese Journal of Quan- rum Electronics(量子电子学报),2011,28(1):65—72. [2]Wu Chunwang,Yang Han,Zhong Xiaojun,et a1.One—way quantum computation with c ircuit quantum electmdynamics [J].Phys Rev A,2010,81(3):034301. [3]Zhou Y L,Li C z.Robust quantum gates via photon tirggering electromagnetieally induced transparency[J].Phys.Rev. A,2011,84(4):044304. [4]Tang Zhaohong,Li Gaoxiang.Entanglement of a two—mode ifeld in a collective three—level atomic system.Phys Rev A. 2011,84:063801. [5] Wei Xiong,Liu Ye.Tunable optimal quantum cloning ma. chines with trapped atoms.Journal of the Optical Society of A- mefica B,2012,29:901—905. [6] Che Chi,Liu Qingfeng,Ma Jing,et a1.Displacement damage effects on the characteristics of quantum dot lasers.Acta Phys Sin,2013,62:094219. [7] Wang Yueyuan,Wang Jicheug,Liu Shutian.The Wigner function and phase properties of superposition of two coherent states with the vacuum state.Chin Phys B,2010,19(7): 074206. [8] Kang Dongpeng,Liao Qinghong,Muhammad Ashfaq Ahamd, et a1.Entropy squeezing of an atom with a k—photon in the Jaynes Cummings mode1.Chin Phys B,2010,19(1): 014206. [9] Peter E,Senellart P,Martrou D,et a1.Exciton—photon strong—coupling regime for a single quantum dot embedded in a microcavity[J].Phys Rev Lett,2005,95(6):067401. [10] Tan Lei,Zhang Yuqing,Wu Mingliu.Quantum phase trnasi— tions for two coupled sites with dipole—coupled effective Jaynes—Cummings mode1.Phys Rev A,2011,84:063816. Wang Fei,Xiao Ming.Output sideband quantum correlations with nonadiabatie elimination.Acta Optica Sinica,2012,32 (12):1227001. [12] Narozhny N B.Sanchez—Mondragon J J.Eberly J H.Coher- ence Versus Incoherence:Collapse and Re—vival in a Simple Quantum Mode1.Phys Rev A,1981,23(1):236—247. [13] Song Tongqiang,Feng Jian,Wang Wenzheng et a1.DynamicM properties of the atoms and the field in the Raman interaction of two atoms with a single—mode cavity field[J].Acta Physica Sinica,1995,44(7):1056—1063(in Chinese). [14] Christina Ginzel,Hans—Jurgen Briegel,Ullrich Martini, Quantum optical master equations:The one—atom laser. Phys Rev A,1993,48:732. [15] 王雪迎.两原子与两耗散腔场量子体系中线性熵的研究. 东北师范大学硕士学位论文,2010. [16] Laussy F P,Valle E D,Tejedor C.Strong coupling of quan- tum dots in microcavities[J].Phys Rev Lett,2008,101(8): 083601. [17] Haroche S,Raimond J M.Exploring the Quantum:Atoms, Cavities and Photons.Oxford University Press,2006. [18] Cao Gui,Yu Zhongyuan,Liu Yumin,et a1.Influences of pure dephasing and incoherent pumping in a coupled quantum dot cavity system and its application.Phys B,201 1,406:3805 —3809. [19] Zhang Wenfu,Chang Heng M A.Some issues on AR models ofr wind field simulation.Chinese Journal of Computational Mechanics,2009,26(1):124—130. [2O] Tan Rong,Li Gaoxinag,Zbigniew Ficek.Squeezed single—at- om laser in a photonic crysta1.Phys Rev A,2008,78: 023833. [21] Yao Pe ̄un,Pathak P K,Illes E,Hughcs S.Nonlinear photolu— minescence spectra from a quantum—dot—cavity system:In— terplay of pump——induced stimulated emission and anharmon— ic cavity QED.Phys Rev B,2010,81:033309. 122 哈尔滨师范大学自然科学学报 2014年第3O卷 The Photon Number Distribution and the Spontaneous Radiation Coupling with Two Atom——cavity Systems Yang Jinbao,Lv Shuchen,Sui Xin (Harbin Normal University) Abstract:In this paper,the system that consists of two identical two—level atoms simultaneously with two coupled single—mode cavity was researched.The cavity field in weak coherent state was considered and the interaction between the two atoms was ignored.Combined the theoretical analysis and numerical computation,the properties of the mean photon number and the spontaneous radiation spectrum are studied. The relationship between the decay rate and the mean photon number of cavity field were discussed,and the influence of the intensity of coupled cavity field to the spontaneous radiation spectrum of the weak coherent coupled cavity field were discussed.The results show that the mean photon number of the cavity was decreased and a vacuum rabi oscillation was observed as the cavity decay rate was increased under the incoherent pumping rate and the decay rate working together on two coupled cavities.When the cavity field intinsic ̄equency became large,the emission spectra of double peak integrallry drifted towards high  ̄equency.On the other hand,while the coupling strength of the cavity field increased,the emission peak of high ̄equency moved towards to the high ̄equency side and low ̄equency moved towards to the low  ̄equency side. Keywords:Quantum master equation;Atom—cavi ̄coupled;The mean photon number;Spontaneous emission spectrum (责任编辑:季春阳)