激光与光电子学进展, 2019, 56 (4): 042701, 网络出版: 2019-07-31   

基于腔结构的可控量子纠缠 下载: 1084次

Controllable Quantum Entanglement Based on Cavity Structure
陆繁 *
作者单位
江南大学理学院, 江苏 无锡 214122
摘要
提出了一种基于V型三能级原子和腔量子电动力学系统的量子纠缠态制备方案。在不同耦合参量条件下,数值分析了任意两腔之间的纠缠演化规律。研究结果表明,当原子以恒定速度通过空腔时,通过改变耦合腔间的耦合参量,可实现对腔间纠缠的控制。利用此方案可同时实现多分量量子纠缠态的产生和控制。
Abstract
A scheme for the preparation of quantum entangled states is proposed based on a V-type three-level atom and cavity quantum electrodynamics system. The dynamics of cavity-cavity entanglement is numerically analyzed under different coupling parameters. The research results show that the cavity-cavity entanglement can be mediated by varying the coupling parameters between cavities when the atom successively passes through the empty cavities at a constant speed. Moreover, this scheme can be utilized for the simultaneous generation and control of multi-component quantum entanglement states.

1 引言

量子纠缠是量子信息和量子计算的重要资源。腔量子电动力学系统是一种分布信息的重要工具,这类系统可用于产生纠缠态。在过去的十年里,科研工作者们已经提出了多种基于腔结构产生纠缠态的方案 [ 1 - 4 ] 。Su等 [ 5 ] 提出一种利用原子的自发辐射和腔的衰减产生最大纠缠态的方案。Ziane等 [ 6 ] 提出了通过控制原子与腔相互作用的时间生成多体纠缠的方案。Lin等 [ 7 ] 提出了利用原子与强失谐腔相互作用产生两体三维最大纠缠态的方案。但是,由于实验条件的限制和来自环境的干扰,纠缠量子位的消相干成为了实现量子技术的主要障碍。消相干影响基于腔量子电动力学系统的量子计算机模型的实现 [ 8 ] 。量子纠缠极易受环境的扰动,在噪音环境中会快速衰减,从而导致纠缠猝死 [ 9 ] 。与纠缠猝死相反, Ficek等 [ 10 ] 发现了纠缠的突然产生。同时,实验上也证明纠缠猝死后能够复苏 [ 11 ] 。此外,Francica等 [ 12 ] 提出用量子芝诺效应抑制纠缠退化,Man等 [ 13 ] 通过调节失谐或腔与附加腔间的耦合强度保护两量子位系统间的纠缠。Wiseman等 [ 14 ] 利用弱测量和反转测量的方法保护多粒子间的纠缠。徐玉虎等 [ 15 ] 发现不同跃迁频率下利用对称失谐可极大地保护纠缠。

在纠缠态的应用方面,量子隐形传态需要用到长时间处于纠缠态的共享量子比特对,而作为许多量子算法基本特征的量子并行性来源于量子状态的相干叠加。消相干不利于量子信息处理, 若在消相干出现时对其进行调控,则可保持量子系统原有的信息。因此,了解量子系统在各种环境下的动力学性质,通过优化耦合参量对系统进行调控,寻找调控和保护纠缠的方法对于量子信息的处理至关重要。闫丽 [ 16 ] 研究了两个二能级原子与共同热库发生相互作用的系统纠缠动力学性质。邢贵超等 [ 17 ] 研究了与热库耦合的光腔内三原子间的纠缠动力学。马蓉等 [ 18 ] 研究了运动原子与场相互作用模型中的量子关联。

为了调控纠缠,本文设计了一个既能生成纠缠态,又能通过参数调节控制纠缠的装置。在这项工作中,受文献[ 19 ]的启发,使用两对相互耦合的无损耗单模腔组成四腔结构模型,令一个 V型三能级原子以恒定速度相继通过耦合腔系统,V型三能级原子的两个偶极跃迁分别与其通过的腔模发生共振耦合。应用该模型研究腔-腔之间的耦合参量对腔模之间纠缠的影响。

2 理论模型

研究一个V型三能级原子依次通过耦合光腔(a)和耦合光腔(b)时,耦合参量对光腔之间纠缠的影响。物理模型示意图如 图1 所示,耦合光腔(a)由光腔c 1 和光腔c 2 耦合构成,耦合强度为 J1 ;耦合光腔(b)由光腔c 3 和光腔c 4 耦合构成,耦合强度为 J2 。V型三能级原子的两个激发态能级|e>与能级| i >之间偶极禁戒。能级|g>为基态,能级|e>和能级| i >到能级|g>的跃迁频率分别为 ω1 , ω2 。当原子通过腔c 1 时,腔c 1 的腔模与原子能级|e>↔|g>间跃迁共振,腔c 1 与原子的相互作用强度为 k1 ;当原子通过腔c 3 时,腔c 3 的腔模与原子能级| i >↔|g>间跃迁共振,腔c 3 与原子的相互作用强度为 k2 ;取原子能级|g>为零能量参考点。假设在原子进入光腔时,腔c 1 ,c 2 ,c 3 ,c 4 均处于真空态,原子处在相干叠加态,满足

|φatom0>=α|e>+β|i>,(1)

式中: α , β 满足归一化条件 α2 + β2 =1。原子经过耦合腔系统(a),初态可表示为

|φ0(0)>=α|e>+β|i>|0c10c2>(2)

在旋波近似下,原子与耦合腔系统(a)所构成的子系统的哈密顿量为( h- =1)

H1=ω1|e><e|+ω2|i><i|+ω1a1a1+ω1a2a2+k1a1|g><e|+a1|e><g|+J1a1a2+a1a2,(3)

式中: a1 ( a1 )为腔c 1 的湮灭(产生)算符; a2 ( a2 )为腔c 2 的湮灭(产生)算符。忽略子系统的损耗,即原子与耦合腔系统(a)构成一个封闭系统,系统的总能量守恒,考虑该子系统中最多只有一个激子。为书写方便,在表示中略去腔场的上标,态矢量从左向右依次代表原子、腔c 1 和腔c 2 。系统的态矢量集合{|e00>,|g00>,| i 00>,|g01>,|g10>}构成了该子系统的一个完备子空间。在此表象中,哈密顿量 H1 的本征值谱为

{0,ω1,ω1-J12+k12,ω1+J12+k12,ω2}(4)

原子经过腔c 1 后,系统的状态可表示为

|φ1(t)>=A1|e00>+A2|g00>+A3|i00>+A4|g01>+A5|g10>(5)

将(2)、(3)、(5)式代入薛定谔方程 jh- (∂| φ ( t )>)/∂ t = H | φ ( t )>,其中j 2 =-1, H 为系统的哈密顿量,| φ ( t )>为 t 时刻系统的状态。

求解得到

A1=αJ12J12+k12+k12J12+k12cosJ12+k12t,(6)A2=0,(7)A3=β,(8)A4=αJ1k1J12+k12cosJ12+k12t-1,(9)A5=-jαk1J12+k12sinJ12+k12t(10)

原子经过耦合腔(b),与腔c 3 发生相互作用,系统的初态可以表示为

φ'0(0)>=|φ1(t)>|0c30c4>(11)

在旋波近似下,原子与耦合腔(b)组成系统的哈密顿量为

H2=ω1|e><e|+ω2|i><i|+ω2a3a3+ω2a4a4+k2a3|g><i|+a3|i><g|+J2a3a4+a3a4(12)

式中: a3 ( a3 )为腔c 3 的湮灭(产生)算符; a4a4 为腔c 4 的湮灭(产生)算符。

原子经过腔c 3 后,同样,整个系统处于单激发态,这时系统的状态可以表示为

|φ2t>=B1|e0000>+B2|g0100>+B3|g1000>+B4|i0000>+B5|g0001>+B6|g0010>(13)

类似地,基矢中标记顺序从左到右依次为原子、腔c 1 、腔c 2 、腔c 3 、腔c 4 。将(5)~(13)式代入薛定谔方程可得

B1=αJ12J12+k12+k12J12+k12cosJ12+k12t,(14)B2=αJ1k1J12+k12cosJ12+k12t-1,(15)B3=-jαk1J12+k12sinJ12+k12t,(16)B4=βJ22J22+k22+k22J22+k22cosJ22+k22t,(17)B5=βJ2k2J22+k22cosJ22+k22t-1,(18)B6=-jβk2J22+k22sinJ22+k22t(19)

为研究光腔之间的纠缠,这里仅考虑原子状态为| g >时腔场的状态,对原子求迹得腔场的约化态,可得

|φ(t)>=B2|0100>+B3|1000>+B5|0001>+B6|0010>(20)

将(20)式归一化得到

|φ(t)>=B2N|0100>+B3N|1000>+B5N|0001>+B6N|0010>,(21)

式中: N = B22+B32+B52+B62

图 1. 物理模型示意图

Fig. 1. Schematic of physical model

下载图片 查看所有图片

3 数值模拟与结果分析

通过数值模拟的方法研究系统参量 k1 , k2 , J1 , J2 对光腔之间纠缠度的影响。采用Wootters [ 20 ] 提出的并发度度量各腔之间的纠缠,并发度定义为

c(ρ)=max{0,λ1-λ2-λ3-λ4},(22)

式中: λ i ( i =1,2,3,4)为矩阵 ρ '= ρ ( σ y 􀱋 σ y ) ρ* ×( σ y 􀱋 σ y )的本征值且 λ1 > λ2 > λ3 > λ4 ,这里 ρ 为系统的密度矩阵, ρ*ρ 的复共轭矩阵, σ y = 0-jj0 为泡利矩阵。并发度0≤ c ( ρ )≤1,当 c ( ρ )=1时,系统处于最大纠缠态,而 c ( ρ )=0时系统处于可分离态。

对于类贝尔态系统,态密度矩阵可表示为

ρ=a00z0by00y*c0z*00d,(23)

式中: a , b , c , d 为非负实数, a + b + c + d =1; y , z 为复数。此类态密度矩阵的矩阵元呈X型,因此称为X态,对于X态,其并发度公式简化为

c(ρ)=2max{0,z-bc,y-ad}(24)

通过求迹的方法可以得到腔c i 与腔c j 构成的子系统的约化密度矩阵。对于(21)式所表示的量子系统,腔c 1 和腔c 2 子系统的约化密度矩阵可通过对腔c 3 和腔c 4 子系统求迹的方法得到

<00|φ(t)>=B2N|01>+B3N|10><11|φ(t)>=0<01|φ(t)>=B5N|00><10|φ(t)>=B6N|00>(25)

容易得到X型密度矩阵为

ρc1-c2=1N2B52+B620000B22B2B3*00B2*B3B3200000(26)

将(26)式与(23)式进行比较,并代入(24)式,计算得到腔c 1 与腔c 2 间并发度的表达式为

cρc1-c2=21N2B2B3*(27)

同理,可得其他两腔间并发度表达式为

cρc2-c3=21N2B2B6*,(28)cρc1-c4=21N2B3B5*,(29)cρc3-c4=21N2B5B6*,(30)cρc1-c3=21N2B3B6*,(31)cρc2-c4=21N2B2B5*(32)

由(1)式可知原子初始时处于两个高能态的叠加态,以下模拟中选取参数 α = β = 2/ 2,即原子初始处于两高能态的最大纠缠态。

首先,考虑原子仅与腔c 1 相互作用,即 k1 ≠0, k2 =0时,腔c 1 -腔c 2 之间纠缠随原子与腔c 1 相互作用时间的变化。当 k1 =0 . 2, k2 =0时,不同耦合常数 J1 下,腔c 1 -腔c 2 之间的纠缠随相互作用时间的变化如 图2 所示。

图 2. 当 k1 =0.2, k2 =0时,不同耦合常数 J1 下腔c 1 与腔c 2 间的纠缠随相互作用时间的变化。(a) J1 较大时;(b) J1 较小时

Fig. 2. Entanglement between cavities c 1 and c 2 versus interaction time for different coupling constants J1 when k1 =0.2, k2 =0. (a) Larger J1 ; (b) smaller J1

下载图片 查看所有图片

模拟发现,原子与耦合光腔(a)中腔c 1 的局域耦合诱导产生了腔c 1 -腔c 2 之间的纠缠,在一定的耦合时间范围内,耦合参量 J1 越大,腔c 1 -腔c 2 之间纠缠量越大。长时相互作用下,如 图2 (a)所示,耦合参量 J1 越大,腔c 1 -腔c 2 之间纠缠变化周期越短,振荡越剧烈;而如 图2 (b)所示,耦合参量 J1 越小,腔c 1 -腔c 2 之间的纠缠变化周期越长,在一个周期内出现两个峰值的时间间隔越短,纠缠死亡时间越长,纠缠曲线呈现M型。当 J1 =0时,腔c 1 与腔c 2 相互独立,即腔c 1 -腔c 2 间纠缠为零。其中的物理机制为原子与光腔之间的局域耦合诱导产生了原子与光腔之间的纠缠,光腔内光子存在的概率一定,而腔-腔耦合作用越强,腔-腔之间光子交换速度越快,从而表现为腔-腔之间的相干性呈周期性变化。因而调节腔-腔耦合参量 J1 ,可以调控腔c 1 -腔c 2 之间的纠缠量。

进一步地,考虑原子先后以相同的速度经过腔c 1 和腔c 3 ,即原子与腔c 1 作用的时间 t1 等于其与腔c 3 作用的时间 t2 ,研究 J2 =0时,不同耦合强度 J1 对腔c 1 -腔c 2 、腔c 1 -腔c 3 、腔c 2 -腔c 3 之间纠缠的影响。 图3 为不同耦合常数 J1 下,腔间纠缠随相互作用时间的变化。

图 3. 不同耦合常数 J1 下腔间纠缠随相互作用时间的变化。(a) k1 = k2 ≠0,腔c 1 、腔c 3 间的纠缠;(b) k1 = k2 =0.1,腔c 1 、腔c 3 间的纠缠;(c) k1 = k2 =0.1,腔c 2 、腔c 3

Fig. 3. Entanglement between cavities versus interaction time for different coupling constants J1 . (a) Entanglement between cavities c 1 and c 3 when k1 = k2 ≠0; (b) entanglement between cavities c 1 and c 3 when k1 = k2 =0.1; (c) entang

下载图片 查看所有图片

原子以相同速度先后经过腔c 1 和腔c 3 ,选取参数 J1 = J2 =0,经过这一过程,如 图3 (a)所示,腔c 1 与腔c 3 处于最大纠缠态,将(15)、(16)、(18)、(19)式代入(31)式同样可得 cρc1-c3 ≡1,即处于相干叠加态的原子将其相干性完全转移成腔c 1 -腔c 3 之间的纠缠。如 图3 (b)所示,在一定的耦合作用时间范围内,耦合参量 J1 越大,腔c 1 -腔c 3 之间的纠缠越小,而如 图3 (c)、(d)所示,腔c 2 -腔c 3 之间的纠缠和腔c 1 -腔c 2 之间的纠缠随着 J1 的增大而增大。由这一现象可以看出,引进附加腔,调节腔间耦合参量可调控纠缠量。有趣的是,在 图3 (b)和 图3 (c)上可以找到一个与耦合常数 J1 无关的时间点,使得腔c 1 -腔c 3 之间和腔c 2 -腔c 3 之间的纠缠量为零。因为从(19)式可以看出, J2 =0时, B6 仅与 k2 有关,与 J1 无关,所以选取合适时间点可使 B6 ≡0,从而使得(28)式和(31)式为0,即存在与 J1 无关的耦合时间点使腔c 1 -腔c 3 之间和腔c 2 -腔c 3 之间的纠缠为0,而此时腔c 1 -腔c 2 之间的纠缠量不为0,改变 J1 可得到不同的腔c 1 -腔c 2 之间的纠缠量。

最后,考虑原子以相同速度先后经过两对对称的耦合腔,即 k1 = k2 ≠0, J1 = J2 ≠0。 图4 为不同耦合常数 J1 下,腔c i -腔c j (i≠j,i,j=1,2,3,4)之间的纠缠随相互作用时间的变化。

图 4. 不同耦合常数 J1 下腔c i -腔c j (i,j=1,2,3;i≠j)之间的纠缠随相互作用时间的变化。(a) k1 = k2 =0 . 1,腔c 1 、腔c 3 间的纠缠;(b) k1 = k2 =0 . 1,腔c 2 、腔c 4

Fig. 4. Entanglement between cavities c i and c j ( i , j =1, 2, 3; i ≠ j ) versus interaction time for different coupling constants J1 . (a) Entanglement between cavities c 1 and c 3 when k1 =

下载图片 查看所有图片

图4 (a)和(b)所示,腔c 1 -腔c 3 ,腔c 2 -腔c 4 间纠缠随腔间耦合参量的增强而振荡加剧,且腔c 1 -腔c 3 间纠缠特性同腔c 2 -腔c 4 间的纠缠特性严格相反,这是因为原子偶极与腔模的局域耦合诱导产生了腔c 1 -腔c 3 之间的纠缠,而腔c 1 -腔c 2 之间、腔c 3 -腔c 4 之间的耦合又诱导产生了腔c 2 -腔c 4 之间的纠缠,由于耦合参量的对称性,所以腔c 2 -腔c 4 之间的纠缠正好补偿了腔c 1 -腔c 3 之间的纠缠损失。由(31)、(32)式计算可得 c ( ρc1-c3 )+ c ( ρc2-c4 )=1,数值模拟结果如 图4 (d)所示,腔c 1 -腔c 3 ,腔c 2 -腔c 4 之间纠缠量之和为常数1,所以它们的纠缠变化规律严格相反。随着 J1 的增强,腔c i 与腔c j ( i-j ≠0,2)之间的纠缠随时间的变化周期缩短,振荡加剧,而 J1 越小,纠缠在一个周期内出现两个峰值的时间间隔越短,纠缠死亡时间越长,纠缠曲线呈现M型,如 图4 (c)所示。由于耦合参量的对称性, 图4 (c)展示了与 图2 (a)相似的变化规律,其中的机理仍然为引入附加光腔,腔-腔之间的耦合参量越大,腔间光子交换的速度越快,腔间纠缠变化周期越短。当 α = β , k1 = k2 , J1 = J2 时,由(15)~(19)式可得 B2 = B5 , B3 = B6 ,进而由(27)~(32)式得出 c ( ρc1-c2 ), c ( ρc1-c4 ), c ( ρc3-c2 ), c ( ρc3-c4 )相等,所以腔c i 与腔c j ( i-j ≠0,2)纠缠规律完全相同。

4 结论

提出了一种利用V型三能级原子通过腔量子电动力学系统生成纠缠态的方案。研究了耦合腔间的不同耦合参量对腔单元之间纠缠的影响。当原子以恒定速度通过两对耦合腔时,改变耦合腔间的耦合参量,可控制腔与腔之间纠缠的变化规律。研究结果表明,耦合腔间的耦合强度增强,腔间纠缠振荡周期缩短,振荡加剧。特别地,利用对称的耦合腔结构,原子与腔间的局域耦合诱导原子的相干叠加态发生转移,腔-腔之间的耦合诱导腔间的纠缠发生转移。通过引进附加腔,可生成、转移纠缠态,调控纠缠量,实现调节耦合参量产生和控制多分量量子纠缠态的目的。

参考文献

[1] Zheng S B, Guo G C. Efficient scheme for two-atom entanglement and quantum information processing in cavity QED[J]. Physical Review Letters, 2000, 85(11): 2392-2395.

[2] Osnaghi S, Bertet P, Auffeves A, et al. Coherent control of an atomic collision in a cavity[J]. Physical Review Letters, 2001, 87(3): 037902.

[3] Wang G Y, Ai Q, Ren B C, et al. Error-detected generation and complete analysis of hyperentangled Bell states for photons assisted by quantum-dot spins in double-sided optical microcavities[J]. Optics Express, 2016, 24(25): 28444-28458.

[4] Chen Y H, Xia Y, Chen Q Q, et al. Fast and noise-resistant implementation of quantum phase gates and creation of quantum entangled states[J]. Physical Review A, 2015, 91(1): 012325.

[5] Su S L, Shao X Q, Wang H F, et al. Scheme for entanglement generation in an atom-cavity system via dissipation[J]. Physical Review A, 2014, 90(5): 054302.

[6] Ziane M, El Guerbouz R, Siyouri F Z, et al. Generation of different classes of multipartite entanglement using cavity-QED[J]. Communications in Theoretical Physics, 2018, 69(2): 131-136.

[7] Lin X M, Zhou Z W, Wu Y C, et al. Preparation of two-qutrit entangled state in cavity QED[J]. Chinese Physics Letters, 2005, 22(6): 1318-1320.

[8] Pellizzari T, Gardiner S A, Cirac J I, et al. Decoherence, continuous observation, and quantum computing: a cavity QED model[J]. Physical Review Letters, 1995, 75(21): 3788-3791.

[9] Pandit M, Das S, Roy S S, et al. Effects of cavity-cavity interaction on the entanglement dynamics of a generalized double Jaynes-Cummings model[J]. Journal of Physics B: Atomic, Molecular and Optical Physics, 2018, 51(4): 045501.

[10] Ficek Z. Tana R. Delayed sudden birth of entanglement[J]. Physical Review A, 2008, 77(5): 054301.

[11] Xu J S, Li C F, Gong M, et al. Experimental demonstration of photonic entanglement collapse and revival[J]. Physical Review Letters, 2010, 104(10): 100502.

[12] Francica F, Maniscalco S, Plastina F. Off-resonant quantum Zeno and anti-Zeno effects on entanglement[J]. Physica Scripta, 2010, T140: 014044.

[13] Man Z X, Xia Y J, Franco R L. Cavity-based architecture to preserve quantum coherence and entanglement[J]. Scientific Reports, 2015, 5: 13843.

[14] Wiseman HM , Milburn GJ .Quantum measurement and control [M]. Cambridge : Cambridge University Press , 2009 : 460 .

[15] 徐玉虎, 任学藻, 刘雪莹. 两任意量子比特Rabi模型的纠缠演化特性[J]. 光学学报, 2018, 38(1): 0127001.

    Xu Y H, Ren X Z, Liu X Y. Entanglement evolution characteristics of quantum Rabi models with two arbitrary qubits[J]. Acta Optica Sinica, 2018, 38(1): 0127001.

[16] 闫丽. 两子系统间纠缠演化特性[J]. 激光与光电子学进展, 2017, 54(3): 032701.

    Yan L. Evolution property of entanglement between two subsystems[J]. Laser & Optoelectronics Progress, 2017, 54(3): 032701.

[17] 邢贵超, 夏云杰. 与热库耦合的光学腔内三原子间的纠缠动力学[J]. 物理学报, 2018, 67(7): 070301.

    Xing G C, Xia Y J. Entanglement dynamics of three atoms in optical cavity coupled to reservior[J]. Acta Physica Sinica, 2018, 67(7): 070301.

[18] 马蓉. 艾合买提·阿不力孜, 艾尔肯江·艾木都拉, 等. 运动原子和场相互作用模型中的量子关联[J]. 激光与光电子学进展, 2018, 55(5): 052701.

    Ma R, Ahmad A, Erkinjan H, et al. Quantum correlations in moving atom-field interaction model[J]. Laser & Optoelectronics Progress, 2018, 55(5): 052701.

[19] Yönaç M, Yu T, Eberly J H. Pairwise concurrence dynamics: a four-qubit model[J]. Journal of Physics B: Atomic, Molecular and Optical Physics, 2007, 40(9): S45-S59.

[20] Wootters W K. Entanglement of formation of an arbitrary state of two qubits[J]. Physical Review Letters, 1998, 80(10): 2245-2248.

陆繁. 基于腔结构的可控量子纠缠[J]. 激光与光电子学进展, 2019, 56(4): 042701. Fan Lu. Controllable Quantum Entanglement Based on Cavity Structure[J]. Laser & Optoelectronics Progress, 2019, 56(4): 042701.

本文已被 8 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!