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

自发参量下转换制备纠缠光子对的特性 下载: 2061次

Property of Entangled Photon Pairs Generated via Spontaneous Parametric Down Conversion
作者单位
1 中国科学技术大学自动化系, 安徽 合肥 230027
2 北京卫星信息工程研究所天地一体化信息技术国家重点实验室, 北京 100086
摘要
通过II型自发参量下转换制备纠缠光子对,在Mathematica环境下分析频率纠缠、频率关联、量子干涉特性及脉冲和晶体参数对特性的影响。结果表明,连续光抽运可以得到最大的频率纠缠度和干涉可见度;脉冲光频宽一定时,随着非线性晶体厚度的增大,双光子的联合光谱变窄,频率纠缠度增大,不可区分性减小,干涉可见度减小;非线性晶体厚度一定时,随着脉冲光频宽的减小,双光子的联合光谱变窄,频率纠缠度增大,不可区分性增大,干涉可见度增大;选取不同的联合光谱函数参数,可以得到具有频率反关联、不关联和正关联特性的双光子。
Abstract
The type-II spontaneous parametric down conversion is used to prepare the entangled photon pairs. In the Mathematica environment, the properties of frequency entanglement, frequency correlation, quantum interference and the influences of parameters for pulse and crystal on the properties are analyzed. The results show that the continuous optical pumping can be used to obtain the maximum frequency entanglement and interference visibility. When the spectral width of pulsed laser is constant, as the nonlinear crystal thickness increases, the combined two-photon spectra become narrow, the frequency entanglement increases, the indistinguishability decreases, and the interference visibility decreases. In contrast, when the nonlinear crystal thickness is constant, as the spectral width of pulsed laser decreases, the combined two-photon spectra become narrow, the frequency entanglement increases, the indistinguishability increases, and the interference visibility increases. For different parameters of the combined spectral function, a photon pair with frequency anti-correlation, non-correlation and positive correlation can be obtained.

1 引言

纠缠光源被广泛应用于量子信息[1]、量子度量衡[2]和量子成像[3]等领域。目前,产生纠缠光子源的有效方法是自发参量下转换(SPDC)[4],其利用高频强光作用非线性晶体,同时产生两个低频纠缠光子对,分别为信号光和闲置光。利用该方法可产生偏振[5]、时间-能量[6]和频率[7]等各种光子自由度形式的纠缠源,且双光子产生的光场具有宽带光谱分布。自发参量下转换的非线性晶体一般为一块负单轴晶体,采用一束非寻常光作为抽运光,如果产生一对偏振相同的信号-闲置光子对[8],则称为I类相位匹配;如果产生的信号-闲置光子对有正交偏振[9],则称为II类相位匹配。在I类和II类相位匹配过程中,如果下转换产生的信号-闲置光具有相同的中心频率(波长),则该过程为简并的;反之则为非简并的。1995年,Kwiat等[10]首次提出通过抽运II型BBO(β相偏硼酸钡晶体,β-BaB2O4)晶体,完成极化纠缠光子对的制备,1999年又提出通过抽运两块I型BBO晶体的方法,实现了高亮度的偏振纠缠光子源[11]。2008年,Baek等[12]研究了第I类SPDC产生的频率简并和非简并纠缠光子对的光谱特性。2009年,卢宗贵等[13]研究了第I类SPDC产生的纠缠双子光谱分布特性。通过SPDC产生的纠缠光源对量子通信、高维量子态、光与物质量子界面和量子度量衡等量子应用研究有重要意义。

纠缠双光子之间的联合光谱可以看作是双光子频率的概率分布,能够反映出双光子的特点。通过研究双光子联合光谱函数在不同参数下的光谱形状,了解双光子特性及不同参数对光谱性能的影响。在自发参量下转换过程中,由于连续光和脉冲光的频宽不同,采用连续光和脉冲光作为抽运源各有利弊。连续光的相干性比脉冲光好,但产生的信号-闲置光子对在抽运光较长的相干时间范围内随机产生,具有较大的时间不确定性[14],导致很难将连续光应用到一些量子过程中,如多光子纠缠态的产生等[15]。脉冲光作抽运源的自发参量下转换过程则可以克服这些缺点。尤其当取频宽很小的抽运脉冲时,纠缠光子对也在抽运脉冲极短的持续时间内产生,此时抽运脉冲可精确控制纠缠光子对的产生,从而有效地降低探测光信号过程中环境噪声对光信号的影响,提高信噪比。因此,在进行光子干涉测量的相关实验中,通常选择脉冲光作为抽运光源。

纠缠光子对逐渐被应用于量子通信和量子计算等领域。笔者等在进行量子定位导航的应用中,需要在地面站接收一对带有延时的纠缠光子对。为了更好地获取和拟合合理信息,本文研究了产生纠缠光子对的信号特性。基于II型自发参量下转换技术制备纠缠光子对,通过仿真实验,研究纠缠双光子中的信号光和闲置光的单光子光谱特性,并根据双光子的联合光谱函数,研究以连续光抽运和在不同的脉冲光频宽和晶体厚度情况下以脉冲光抽运产生双光子的频率纠缠、频率关联和量子干涉特性,以及不同参数对其性能的影响。

2 II型自发参量下转换光子对的产生

典型的II型自发参量下转换制备纠缠光子对的实验装置如图1所示,其中,图1(a)是采用共线匹配的方式产生偏振纠缠双光子,激光器产生一束抽运光入射到单片II型非线性晶体中,抽运光与II型非线性晶体相互作用,经过偏振分束器(PBS)后,输出两束偏振相互正交的信号光和闲置光,此时完成纠缠光源的制备。单光子探测器(D1和D2)把从偏振分束器输出的光信号转为电信号,从中可获得制备出的包括光子数量在内的纠缠光子特性。将D1和D2输出的电信号送入符合计数器(CC)中,CC对信号光和闲置光进行符合判决,输出符合脉冲,并统计单位时间内符合脉冲数,得到符合计数。

图 1. II型自发参量下转换制备纠缠光子对的实验装置示意图。(a)共线匹配;(b)非共线匹配

Fig. 1. Schematic of experimental device for generation of entangled photon pairs via type-II spontaneous parametric down conversion. (a) Collinear matching; (b) non-collinear matching

下载图片 查看所有图片

图1(b)利用II型非线性晶体,采用非共线匹配方式直接产生双光子纠缠态,抽运光与晶体相互作用所输出的两束光通过半波片HWP0出射后,再分别经过双折射晶体C1和C2,沿着C1出射的光在经过半波片HWP1和四分之一波片QWP1后进入偏振器P1和单光子探测器D1中,沿着C2出射的光直接由偏振器P2和单光子探测器D2接收,其中偏振器P1和P2均由两个堆叠的偏振分束器组成,HWP0、C1和C2可补偿两束光因双折射造成的离散效应,HWP1和QWP1可制备其他Bell态[10]

对于II型自发参量下转换,其过程可以表示为e→o + e,其中输入的抽运光为e光,产生的输出信号光为o光,另一个输出闲置光为e光,这两个光子的偏振方向互相垂直。其中o光被称为寻常光,其在晶体传播中的折射率固定不变;e光被称为非寻常光,是偏振方向垂直于o光振动的光,在传播中的折射率随入射光方向发生变化。当非共线匹配时,II型SPDC过程在频率简并和频率非简并情况下光场的空间结构如图2所示,其中,抽运光的分布为两个相交的圆锥面,其上圆为e光,下圆为o光;当共线匹配时,这两个圆锥面是相切的。当相位匹配切割角逐渐增大时,两个圆锥面都向抽运光方向靠拢,发生交叠,此时两圆交叉的两点可能是e光也可能是o光,如果其中一个为e光,则另一个一定为o光,这样在两方向上的一对光子便形成偏振纠缠的双光子态。调整合适的抽运光方向与晶体光轴的夹角,使两个锥面相交,两锥体交线的两个方向各在抽运光的一侧,则沿这两个方向出射的下转换光子即处在偏振纠缠态,其表达式为:|ψ>=1/2[|H1V2>+exp(iϕ)|V1H2>],式中:H1H2为水平偏振态,V1V2为垂直偏振态,角度ϕ可以通过附加一个双折射相移器或小幅度旋转下转换晶体本身来改变大小,其在0到π间变化,互换不同偏振可产生4种Bell态,即|ψ+>=1/ 2[|H1V2>+|V1H2>]、|ψ->=1/ 2[|H1V2>-|V1H2>]、|ψ+>=1/ 2[|H1H2>+|V1V2>]和|ψ->=1/ 2[|H1H2>-|V1V2>]。

图 2. II型SPDC光场示意图

Fig. 2. Illustration of type-II SPDC light field

下载图片 查看所有图片

3 双光子联合光谱以及单光子光谱与抽运频宽和晶体厚度的关系

根据量子理论,自发参量下转换所产生的纠缠光子对波函数|ψ>可表示为[16]

|ψ>=CdωsdωiΦ(ωs,ωi)εp(ωs+ωi)as(ωs)ai(ωi)|0>,(1)

式中:C为常系数;ωsωi分别为信号光和闲置光的频率, as(ωs)和 ai(ωi)分别为信号光和闲置光的产生算符;|0>为真空态;Φ(ωs,ωi)为相位匹配函数,其表达式为

Φ(ωs,ωi)=sinc(ΔL/2),(2)

式中:sinc(ΔL/2)= 2sin(ΔL/2)ΔL;L为非线性晶体厚度;Δ为径向相位失配量,其表达式为

Δ=ks(ωs)+ki(ωi)-kp(ωp),(3)

其中ks(ωs)、ki(ωi)和kp(ωp)分别为信号光、闲置光和抽运光的光波数,且有kr(ωr)=nrωr/c(r=p,s,i),nr为折射率,c为真空中的光速,ωr为光频率。当Δ=0即kp(ωp)=ks(ωs)+ki(ωi)时,双光子满足动量守恒即相位匹配条件,此时自发参量下转换过程的效率最大。

(1)式中,εp(ωs+ωi)为抽运线性函数,其表达式为

εpωs+ωi=exp[-(ωs+ωi-Ωp)2/σp2],(4)

式中:σpΩp分别为抽运光的频宽和中心频率。

抽运线性函数是抽运脉冲函数时域分布的傅里叶变换,满足能量守恒,即下转换过程仅允许产生光子对的频率之和等于抽运频率,即ωp=ωs+ωi的情况。自发参量下转换产生的双光子的联合光谱函数S(ωs,ωi)是相位匹配函数和抽运线性函数乘积的平方,即

S(ωs,ωi)=Φ(ωs,ωi)εp(ωs+ωi)2(5)

S(ωs,ωi)可以看作光子频率的概率分布,体现了联合光谱的强度。

直接测量双光子联合光谱,并将其与双光子的频率纠缠特性联系起来,可实现纠缠量化。Fedorov等[17]提出在脉冲抽运时信号光子与闲置光子间的频率纠缠度:将Schmidt数值R近似为单光子光谱宽度与双光子联合光谱宽度的比值,其表达式为

R=ΔνsΔνc,(6)

式中:Δνc为双光子联合光谱的频率宽度,即联合光谱图形在νi=νs轴上的频率宽度;Δνs为信号光s的单光子光谱的半峰全宽(FWHM),即峰值高度一半处的峰宽度;R在0到∞范围内波动。

根据频率一致纠缠源的判据[18],通过比较纠缠光子对的信号光子与闲置光子频谱之和的频宽ΔΛ+和信号光子与闲置光子频谱之差的频宽ΔΛ-的大小可判断频率的关联特性。在双光子联合光谱上,ΔΛ+为联合光谱图形在正对角线νi=νs上的频率宽度即ΔνcΛ-为联合光谱图形在负对角线νi=-νs上的频率宽度。若ΔΛ+Λ-,则双光子为频率反关联;若ΔΛ+Λ-,则双光子为频率不关联;若ΔΛ+Λ-,则双光子为频率正关联。

还可通过对双光子联合光谱的积分运算来体现量子干涉可见度,对信号光子和闲置光子可区分性的度量直接体现在干涉可见度上。双光子的干涉可见度V可以表示为[16]

V=dωsdωiA(ωs,ωi)A(ωi,ωs)dωsdωiA(ωs,ωi)2,(7)

式中:A(ωs,ωi)=Φ(ωs,ωi)εp(ωs+ωi)为双光子联合光谱的幅度函数,当且仅当 A(ωs,ωi)= A(ωi,ωs),即联合光谱的幅度函数关于正对角线νi=νs对称时,干涉可见度V达到最大值1,两个单光子的光谱完全重合,不可区分[19];两个单光子光谱的差异越大,不可区分性越小,双光子的干涉可见度也越小。

推导双光子的联合光谱以及单光子光谱与抽运频宽和晶体厚度的关系,设抽运光、信号光和闲置光的频率ωr分别为

ωr=Ωr+νr,(8)

式中:Ωr为对应光场的中心频率;νr为对应光场频率的偏移量,即差分频率。根据能量守恒,差分频率满足νs+νi=νp

将光波数kr(νr)=nr(Ωr+νr)/c按级数展开取一阶项,忽略二阶以上项,可得[20]

kr(νr)=Kr+νrμr(Ωr),(9)

式中:Kr为展开式的常数项,满足零阶分量的动量守恒,即Ks+Ki=Kp;μr(Ωr)为沿晶体传播对应光场的光束在中心频率Ωr处测得的群速度。

将(9)式代入(3)式,可得到径向相位失配量Δ与差分频率νr的关系,即

Δ=Ks+Ki-Kp+νsμs(Ωs)+νiμi(Ωi)-νpμp(Ωp),(10)

将其化简为

Δ=-(νpD++0.5ν-D),(11)

式中:νp=ωp-Ωp=ωs+ωi-Ωp;ν-=νi-νs=ωi-ωs+Ω;Ω=Ωs-Ωi为信号光和闲置光中心频率的偏移量;D为闲置光与信号光在晶体中的逆群速度之差,其与μr(Ωr)之间的关系为D=1i(Ωi)-1/μs(Ωs);D+为信号光和闲置光在晶体中传播的平均逆群速度与抽运光传播的逆群速度差,其表达式为D+=0.5[1/μs(Ωs)+1/μi(Ωi)]-1/μp(Ωp);DD+的大小、正负均取决于晶体的自身性质,即对特定的晶体如BBO随中心频率(波长)的变化而各不相同。

将(11)式代入(2)式,可得

Φ(ωs,ωi)=sinc{0.5L[-(ωs+ωi-Ωp)D+-0.5D(ωi-ωs+Ω)]}(12)

Ω≠0时,即信号光和闲置光的中心频率不同,则两个输出信号频率是非简并的,将(12)式中ωsωi交换后,可得

Φ(ωi,ωs)=sinc{0.5L[-(ωi+ωs-Ωp)D+-0.5D(ωs-ωi-Ω)]},(13)

此时|Φ(ωs,ωi)|≠|Φ(ωi,ωs)|,因此相位匹配函数不对称。因为抽运线性函数εp(ωs+ωi)能够决定可用于下转换泵能量的范围,相位匹配函数Φ(ωs,ωi)能够决定分配给两个下变频光子的能量,所以相位匹配函数不对称时会导致能量分配不均,进而体现为信号光和闲置光的单光子光谱各不相同,可区分。

Ω=0时,即信号光和闲置光的中心频率相等,则频率简并,设Ωs=Ωi= Ω-,又因为Ωp=Ωs+Ωi,所以Ωp=2 Ω-,此时(12)式转化为

Φ(ωs,ωi)=sinc{0.5L[-(ωs+ωi-2Ω-)D+-0.5D(ωi-ωs)]}(14)

同时当D+=0时,将ωsωi交换后,相位匹配函数关于ωsωi对称,即|Φ(ωs,ωi)|=|Φ(ωi,ωs)|,下转换所产生的两个单光子的光谱完全相同,此时的双光子是完全不可区分的。

将(8)式关于信号光频率ωs和闲置光频率ωi的表达式代入(14)式,可得相位匹配函数Φ(ωs,ωi)与信号光差分频率νs和闲置光差分频率νi的关系为

Φ(Ω-+νs,Ω-+νi)=sinc{0.5L[-(νs+νi)D+-0.5D(νi-νs)]}(15)

以同样的方式,可将(4)式的抽运线性函数εp(ωs+ωi)转化为

εp2Ω-+νs+νi)=exp[-(νs+νi)2/σp2],(16)

将(15)和(16)式代入(5)式,可得中心频率差值Ω=0频率简并情况下,自发参量下转换产生的双光子的联合光谱函数为

S(νs,νi)=exp[-(νs+νi)2/σp2]sinc{0.5L[-(νs+νi)D+-0.5D(νi-νs)]}2(17)

由(17)式可以看出,影响下转换产生的双光子联合光谱性能的参数有抽运光的频宽σp与非线性晶体厚度L

单光子光谱区别于双光子联合光谱,可以单独反映双光子各自的变化,通过研究单光子光谱可以了解单光子的变化情况。将(17)式的联合光谱函数分别沿信号光和闲置光的差分频率方向作投影,得到信号光的单光子光谱Ss(νs)和闲置光的单光子光谱Si(νi)为

Ss(νs)=dνiS(νs,νi)Si(νi)=dνsS(νs,νi)(18)

将(8)式代入(12)式,得到相位匹配函数Φ(ωs,ωi)与信号光差分频率νs和闲置光差分频率νi的关系为

Φ(Ωs+νs,Ωi+νi)=sinc{0.5L[-(νs+νi)D+-0.5D(νi-νs+Ω)D]}(19)

以同样的方式,将(4)式的抽运线性函数εp(ωs+ωi)转化为

εp(Ωs+νs+Ωi+νi)=exp[-(νs+νi)2/σp2],(20)

将(19)式和(20)式代入(5)式,得到中心频率差值Ω≠0频率非简并的情况下自发参量下转换产生的双光子的联合光谱函数为

S(νs,νi)=exp[-(νs+νi)2/σp2]·sinc{0.5L[-(νs+νi)D+-0.5D(νi-νs+Ω)]}2,(21)

将(21)式代入(18)式,得到Ss(νs)和Si(νi)分别为

Ss(νs)=dνiexp[-(νs+νi)2/σp2]·sinc{0.5L[-(νs+νi)D+-0.5D(νi-νs+Ω)]}2Si(νi)=dνsexp[-(νs+νi)2/σp2]·sinc{0.5L[(νs+νi)D+-0.5D(νi-νs+Ω)]}2(22)

下面将通过对自发参量下转换过程的仿真实验来分析双光子的联合光谱和单光子光谱在不同参数下表现出的频率纠缠、频率关联以及量子干涉特性。

4 不同参数下的光谱特性实验及其结果分析

详细分析了信号光和闲置光的单光子光谱特性、连续光抽运时双光子的光谱特性、脉冲光抽运时脉冲的频宽和晶体的厚度对双光子光谱特性的影响,以及具有特殊频率关联的双光子间特性。其中所有的仿真实验均是在Mathematica软件上完成的。采用绘图函数Plot对信号光和闲置光的单光子光谱作图,采用绘制三维图形函数Plot3D对相位匹配函数Φ( Ω-+νs, Ω-+νi)、抽运线性函数εp(2 Ω-+νs+νi)和双光子联合光谱函数S(νs,νi)作三维空间图;采用密度函数Density Plot作二维俯视图,绘图点Plot Points为50,绘图范围PlotRange为Full,颜色函数Color Function设为“Sunset Colors”,在Mathematica中将各个函数的表达式代入Plot/Plot3D/Density Plot[Ss(νs)/Si(νi)/Φ( Ω-+νs, Ω-+νi)/εp(2 Ω-+νs+νi)/S(νs,νi),νsνiσpD+DL]中进行仿真实验绘图。

4.1 单光子光谱特性分析

采用脉冲光作为自发参量下转换的抽运光进行两部分实验:1) 在信号光和闲置光的中心频率差值Ω≠0,即频率非简并的情况下,观察参数Ω的变化对单光子光谱特性的影响;2) 在信号光和闲置光的中心频率差值Ω=0,即频率简并的情况下,观察参数D+的变化对单光子光谱特性的影响。

通过(22)式Ss(νs)和Si(νi)所示的光谱函数可作出单光子光谱图,首先在Mathematica中使用数值积分函数NIntegrate对双光子联合光谱函数中信号光差分频率νs(或闲置光差分频率νi)进行积分,然后使用Plot函数绘制出单光子光谱函数随闲置光差分频率νi(或信号光差分频率νs)变化的图像。在实验中,选择脉冲光的频宽为σp=1 nm,晶体厚度L=1 mm,D+=-1.82×10-13 s/mm,D=1×10-13 s/mm[21]。在第一个实验中,分别取Ω=5 THz和Ω=-5 THz,将所有参数代入(22)式中,Ω>0和Ω<0时信号光子和闲置光子的单光子光谱Ss(νs)和Si(νi)分别如图3(a)和(b)所示。在第二个实验中取Ω=0,分别使D+=-1.82×10-13 s/mm和D+=0[22],其他参数不变,将所有参数代入(22)式中,信号光子和闲置光子的单光子光谱Ss(νs)和Si(νi)如图3(c)和(d)所示。在图3中蓝色实线为信号光(s)的单光子光谱,红色虚线为闲置光(i)的单光子光谱。

图3(a)和(b)可以看出,当Ω≠0时,频率非简并,此时实验不满足相位匹配条件,即Δ≠0。信号光和闲置光的单光子光谱的对称轴存在偏移,且偏移量为Ω,此时信号光和闲置光在频率上是可区分的,并且这种偏移的程度越大,信号光和闲置光的单光子光谱可区分性越大。从图3(c)和(d)可以看出,当Ω=0时,频率简并,即信号光和闲置光的中心频率相同,所以两个单光子光谱对称轴均在差分频率为零的轴线处。另一方面,当D+≠0时,相位匹配函数Φ(ωs,ωi)不对称,从图3(a)、(b)和(c)可以看出,双光子的单光子光谱幅度各异,是可区分的;当D+=0时,信号光和闲置光的单光子光谱完全重合,是不可区分的。

图 3. 不同参数下信号光(s)和闲置光(i)的单光子光谱图。(a) Ω>0;(b) Ω<0;(c) D+≠0;(d) D+=0

Fig. 3. Single photon spectra of signal light (s) and idle light (i) under different parameters. (a) Ω>0; (b) Ω<0; (c) D+≠0; (d) D+=0

下载图片 查看所有图片

实验结果表明:1) 当频率非简并即Ω≠0时,自发参量下转换过程仍能发生,但不满足相位匹配条件,存在相位失配量Δ,因此自发参量下转换的效率下降,产生的双光子的频率纠缠度不高;2) 当Ω=0且D+=0时,若相位匹配函数对称,信号光和闲置光的单光子光谱完全重合,不可区分;若相位匹配函数不对称,两个单光子的光谱存在差异,是可区分的。因此,在实际实验中,应该尽量减少相位匹配函数的不对称性,以提高双光子的不可区分性和干涉可见度。

4.2 连续抽运光下的双光子的特性分析

连续光是指抽运光频宽σp→0的情况。当连续光作为自发参量下转换的光源时,进行相位匹配函数、抽运线性函数、双光子联合光谱以及单光子光谱的仿真实验。实验中所选取参数分别为:L=1 mm,D+=-1.82×10-13 s/mm,D=1×10-13 s/mm。令信号光和闲置光的差分频率νsνi在[-8 THz,8 THz]范围内变化。将各参数代入(15)式和(16)式,得到以信号光、闲置光差分频率为变量的相位匹配函数及抽运线性函数,如图4(a1)和(b1)所示,其二维俯视图分别如图4(a2)和(b2)所示;将各参数代入(17)式,得到如图4(c)所示的双光子联合光谱函数图,其二维俯视图如图4(c3)所示;将各参数代入(22)式,令νsνi在[-20 THz,20 THz]范围内变化,得到如图4(d)所示的单光子光谱图。

图4(a1)中可以看出,相位匹配函数整体呈“山”字状,结合图4(a2)可以看出,相位匹配函数关于νi=D+νs轴对称,且在该对称轴上函数达到最大值,将(15)式Φ( Ω-+νs, Ω-+νi)中的信号光差分频率νs和闲置光差分频率νi进行交换可以得到Φ( Ω-s, Ω-+νi),因为D+≠0,所以Φ( Ω-+νs, Ω-+νi)≠Φ( Ω-+νi, Ω-+νs),相位匹配函数Φ( Ω-+νs, Ω-+νi)关于νi=νs轴不对称,即函数不具有对称性。从图4(b1)可见,抽运线性函数的空间图形呈倒“T”状,从图4(b2)可见,εp(2 Ω-+νs+νi)仅在负对角线νi=-νs轴上取得函数最大值,因为连续光的频宽σp趋近于0,所以其函数图像的宽度很窄,表现为负对角线上的一条细线。

图 4. 抽运源为连续光时的各函数图形。(a1)(a2)相位匹配函数的三维图和二维图;(b1)(b2)抽运线性函数的三维图和二维图;(c1)(c2)联合光谱函数的三维图和二维图;(d)单光子光谱图

Fig. 4. Source function graphs when continuous light used as pump. (a1) (a2) 3D and 2D plots of phase matching function; (b1) (b2) 3D and 2D plots of pumping linear function; (c1) (c2) 3D and 2D plots of combined spectral function; (d) single photon spectra

下载图片 查看所有图片

由于双光子的联合光谱是由抽运线性函数与相位匹配函数乘积的平方得到的,所以从图4(c1)可以看到,联合光谱函数的三维图形呈峰状,从图4(c2)可以看到,联合光谱函数是一条沿着抽运线性函数方向、且与相位匹配函数的νi=D+νs轴坐标范围一致的细线。从图4(c2)中可以明显观察到,横坐标信号光差分频率νs越大,纵坐标闲置光的差分频率νi越小,联合光谱图形在νi=νs上的频率宽度为ΔΛ+,在νi=-νs上的频率宽度为ΔΛ-,测量得到ΔΛ+→0,ΔΛ-≈12.43 THz,明显ΔΛ+Λ-,说明以连续光作为自发参量下转换抽运光时产生的信号光子和闲置光子满足频率反关联特性。以连续光作为自发参量下转换过程的抽运源时,如图4(c2)所示的联合光谱图像关于νi=νs轴对称,所以其幅度函数A(ωs,ωi)满足 A(ωs,ωi)= A(ωi,ωs),由(7)式可知,干涉可见度V=1,达到最大值。

图4(d)中可以看出,信号光(s)和闲置光(i)的单光子光谱是重合的,这是因为虽然相位匹配函数Φ(ωs,ωi)不对称,但σp接近零,则抽运线性函数图像可以看作是宽度接近零的细线,且关于νi=νs轴对称,所以其被相位匹配函数截断的细线即双光子联合光谱函数S(νs,νi)的图像也关于νi=νs轴对称,由(18)式可知,Ss(νs)=Si(νi),故两者的单光子光谱完全重叠不可区分。从图4(d)可见,信号光(s)的单光子光谱的半峰全宽为Δνs≈5.17 THz,因为连续光的频宽σp→0,图4(c2)中联合光谱宽度Δνc→0,因此R→∞,达到最大的频率纠缠度。

4.3 脉冲光的频宽和晶体厚度对双光子的特性影响

当采用脉冲光作为抽运光源进行自发参量下转换时,脉冲光的频宽σp和非线性晶体厚度L分别为抽运线性函数εp(2 Ω-+νs+νi)和相位匹配函数Φ( Ω-+νs, Ω-+νi)的重要参数,通过分别固定晶体厚度时改变脉冲光频宽和固定脉冲光频宽时改变晶体厚度,确定二者对双光子的联合光谱、单光子光谱以及双光子频率纠缠度和干涉可见度的影响。

1) 固定晶体的厚度L,改变脉冲光的频宽σp

选取参数分别为:L=1 mm,D+=-1.82×10-13 s/mm,D=1×10-13 s/mm,脉冲光频宽σp分别取0.5,2,5 nm,令信号光和闲置光的差分频率νsνi在[-8 THz,8 THz]范围内变化,将各参数代入(16)式,得到以νsνi为变量的抽运线性函数三维图,依次如图5(a1)、(a2)和(a3)所示;将各参数代入(17)式,得到双光子联合光谱三维图,依次如图5(b1)、(b2)和(b3)所示,其对应的俯视图如图5(c1)、(c2)和(c3)所示;将各参数代入(22)式,令νsνi在[-20 THz,20 THz]范围内变化,得到双光子的单光子光谱依次如图5(d1)、(d2)和(d3)所示。

图 5. 不同脉冲光频宽下的抽运线性函数三维图、联合光谱函数三维图及二维图和单光子光谱图。(a1) (b1) (c1) (d1) σp=0.5 nm;(a2) (b2) (c2) (d2) σp=2 nm;(a3) (b3) (c3) (d3) σp=5 nm

Fig. 5. 3D plots of pumping linear functions, 3D plots and 2D plots of combined spectral functions, single photon spectra under different spectral widths of pulsed laser. (a1) (b1) (c1) (d1) σp=0.5 nm; (a2) (b2) (c2) (d2) σp=2 nm; (a3) (b3) (c3) (d3) σp=5 nm

下载图片 查看所有图片

采用不同的脉冲光频宽进行仿真实验时,晶体厚度L固定的情况下相位匹配函数Φ( Ω-+νs, Ω-+νi)的幅度不变,由于DD+L的取值与连续光作为抽运光时所选参数相同,此时相位匹配函数的三维图同图4(a1)一样,从图5(a1)、(a2)和(a3)可以看出,脉冲光频宽σp取值越大,抽运线性函数εp(2 Ω-+νs+νi)在负对角线νi=-νs附近的频率宽度越大,所以εp(2 Ω-+νs+νi)和Φ( Ω-+νs, Ω-+νi)在空间相交重叠的范围也越大。因此, S(νs,νi)的图形越宽越长,即Δνc越大,但联合光谱图形的频率范围受相位匹配函数中νi=D+νs轴坐标范围的限制,在图5(c1)、(c2)和(c3)中表现为联合光谱图形越来越靠近白色细线,即νi=D+νs轴,直到在这条轴附近的一定范围后,图形的频率宽度和范围将不再随σp的增大而增大。图5(d1)、(d2)和(d3)中所有信号光(s)和闲置光(i)的单光子光谱均是可区分的,这是因为D+≠0,相位匹配函数Φ(ωs,ωi)不对称,所以分配给信号光和闲置光的能量不均,双光子的单光子光谱不一致,是可区分的,且随着σp的变大,参与频率下转换的抽运光子变多,用于下转换的泵能量频率范围变大,信号光和闲置光分配能量差异变大,因此单光子光谱图形越来越宽即Δνs越大,双光子之间的单光子光谱差异越大。

不同脉冲光频宽下的频率纠缠度与干涉可见度如表1所示,其中,双光子光谱的频率宽度Δνc可通过图5(c1)、(c2)和(c3)测得,信号光s的单光子光谱的半峰全宽Δνs可以通过图5(d1)、(d2)和(d3)测得,根据(6)式求出不同脉冲光频宽σp下的频率纠缠度,根据(7)式可以求出不同脉冲光频宽σp下的干涉可见度V,使用Mathematica中数值积分函数NIntegrate对(7)式的分子和分母进行二重积分,可得σp=0.5 nm时,V≈4.748/4.922≈0.96;σp=2 nm时,V≈12.436/16.098≈0.77;σp=5 nm时,V≈15.892/22.769≈0.69。从表1可以看出,随着σp的变大,Δνc和Δνs均变大,而频率纠缠度R变小,双光子的干涉可见度变小。

表 1. 不同脉冲光频宽下的频率纠缠度与干涉可见度

Table 1. Frequency entanglement and interference visibility under different spectral widths of pulsed laser

σp /nmΔνc /THzΔνs /THzRV
0.50.675.868.750.96
21.136.946.140.77
52.2612.585.570.69

查看所有表

2) 固定脉冲光的频宽σp,改变晶体的厚度L

选取各参数分别为:脉冲光的频宽为σp=1 nm,D+=-1.82×10-13 s/mm,D=1×10-13 s/mm,晶体厚度L分别取0.5,1,3 mm,令信号光和闲置光的差分频率νsνi在[-8 THz,8 THz]范围内变化,将σp=1 nm代入(16)式,得到以νsνi为变量的抽运线性函数三维图,如图6(a)所示;将参数D+D和不同的L代入(15)式,得到以νsνi为变量的相位匹配函数三维图,依次如图6(b1)、(b2)和(b3)所示;将各参数代入(17)式,依次得到如图6(c1)、(c2)和(c3)所示的双光子联合光谱函数三维图,其对应的俯视图如图6(d1)、(d2)和(d3)所示;将各参数代入(22)式,令νsνi在[-20 THz,20 THz]范围内变化,依次得到如图6(e1)、(e2)和(e3)所示的单光子光谱图。

采用不同厚度的非线性晶体进行仿真实验时,脉冲光频宽σp固定为1 nm,抽运线性函数εp(2 Ω-+νs+νi)的三维图如图6(a)所示,从图6(b1)、(b2)和(b3)可以看出,L值越大,相位匹配函数Φ( Ω-+νs, Ω-+νi)在νi=D+νs轴附近的频率宽度越小,εp(2 Ω-+νs+νi)和Φ( Ω-+νs, Ω-+νi)在空间相交重叠的范围也越小,因此,S(νs,νi)图形越窄越短,即Δνc越小。因为由抽运线性函数产生的下转换光场被限制在相位匹配函数的νi=D+νs轴附近,所以相位匹配函数Φ( Ω-+νs, Ω-+νi)频率宽度越小,下转换产生的双光子联合光谱图形越靠近νi=D+νs轴,该轴与联合光谱图形相交的频率范围不变。图6(e1)、(e2)和(e3)中所有信号光(s)和闲置光(i)的单光子光谱均是可区分的,D+≠0使相位匹配函数Φ(ωs,ωi)不对称,因此双光子在下转换时分配的能量不均,双光子的单光子光谱不一致,是可区分的;L值变大,能量分配差异变大,双光子之间的单光子光谱差异随之变大,且Δνs变小。

不同晶体厚度下的频率纠缠度与干涉可见度如表2所示,其中,双光子光谱的频率宽度Δνc可以通过图6(d1)、(d2)和(d3)测得,信号光s的单光子光谱的半峰全宽Δνs可以通过图6(e1)、(e2)和(e3)测得,根据(6)式求出不同晶体厚度L下的频率纠缠度,根据(7)式可以求出不同晶体厚度L下的干涉可见度V,使用Mathematica中数值积分函数NIntegrate对(7)式的分子和分母进行二重积分,可以得到L=0.5 mm时,V≈13.903/14.025≈0.99;L=1 mm时,V≈5.587/6.651≈0.84;L=3 mm时,V≈2.071/3.514≈0.58。从表2可以看出,随着L值变大,Δνc和Δνs均变小,频率纠缠度R变大,同时双光子的干涉可见度变小。

本小节可以得出以下结论:1) 为了得到尽量大的双光子频率纠缠度和干涉可见度,脉冲光频宽σp应尽量小,σp趋近0时,脉冲光抽运变为连续光抽运,达到最大纠缠度和干涉可见度;2) 确定σp后,为提高频率纠缠度还可以增大晶体厚度L,但L过大时会降低干涉可见度,所以应该选取合适的L

图 6. 不同晶体厚度下的抽运线性函数三维图、相位匹配函数三维图、联合光谱函数三维图及二维图、单光子光谱图。(a) σp=1 nm;(b1) (c1) (d1)(e1) L=0.5 nm;(b2) (c2) (d2)(e2) L=1 mm;(b3) (c3) (d3)(e3) L=3 mm

Fig. 6. 3D plots of pumping linear functions, 3D plots of phase matching functions, 3D plots and 2D plots of combined spectral functions, single photon spectra under different crystal thicknesses. (a) σp=1 nm; (b1) (c1) (d1) (e1) L=0.5 nm; (b2) (c2) (d2) (e2) L=1 mm; (b3) (c3) (d3) (e3) L=3 mm

下载图片 查看所有图片

表 2. 不同晶体厚度下的频率纠缠度与干涉可见度

Table 2. Frequency entanglement and interference visibility under different crystal thicknesses

L /mmΔνc /THzΔνs /THzRV
0.51.7411.206.440.99
10.916.427.050.84
30.292.528.690.58

查看所有表

4.4 具有特殊频率关联双光子的特性分析

采用自发参量下转换技术不仅可以制备以上所述的频率反关联特性的双光子,也可以通过改变双光子联合光谱函数的参数,使脉冲抽运的自发参量下转换过程产生具有频率不关联和正关联特性的双光子[23]。由(15)式可知,当D+=0时,相位匹配函数Φ( Ω-+νs, Ω-+νi)具有对称性,即函数关于νi=νs轴对称,又因为脉冲光的频宽σp和非线性晶体的厚度L不同时,双光子联合光谱图像也不同,因此一定参数下自发参量下转换过程可以制备出具有频率不关联或正关联特性的双光子。

当选取脉冲光的频宽为σp=3 nm,非线性晶体的厚度为L=2 mm,D+=0,D=1×10-13 s/mm时,令信号光和闲置光的差分频率νsνi在[-8 THz,8 THz]范围内变化,将各参数代入(17)式,得到双光子的联合光谱图如图7(a)所示。保持D+=0,D=1×10-13 s/mm,使脉冲光频宽σp=10 nm,晶体厚度L=5 mm[24],令信号光和闲置光的差分频率νsνi在[-8 THz,8 THz]范围内变化,将各参数代入(17)式,此时双光子的联合光谱图如图7(b)所示。

图 7. 不同频率关联的双光子联合光谱。(a)不关联;(b)正关联

Fig. 7. Combined two-photon spectra under different frequency correlation. (a) Non-correlation; (b) positive correlation

下载图片 查看所有图片

图7(a)可以看出,双光子联合光谱图形关于νi=νsνi=-νs都是对称的,测量得到联合光谱图形在正对角线νi=νs上的频率宽度ΔΛ+≈4.52 THz,在负对角线νi=νs上的频率宽度ΔΛ-≈4.52 THz,因为ΔΛ+Λ-,所以信号光子和闲置光子是不关联的;从图7(b)可见,横坐标信号光差分频率νs越大,纵坐标闲置光的差分频率νi也越大,测量得到ΔΛ+≈19.23 THz,ΔΛ-≈2.08 THz,明显ΔΛ+Λ-,所以信号光子和闲置光子满足频率正关联特性。

仿真实验表明,改变双光子联合光谱函数的参数可以使纠缠光子对之间的频率关联特性从反关联向不关联和正关联演化,这对量子密码学多光子纠缠态的产生和多光子干涉实验等具有重要意义。

5 结论

分别对以连续光和脉冲光作为自发参量下转换的抽运光时,双光子的联合光谱和单光子光谱进行了仿真实验,对比研究了不同参数下对纠缠光子对的频率纠缠、频率关联和量子干涉特性的影响,可得到以下结论。1) 连续光作为抽运光时产生的双光子干涉可见度最大,脉冲光作为抽运光时,频宽越小或晶体厚度越小,干涉可见度越大;2) 连续光作为抽运光时产生的双光子频率纠缠度最大,脉冲光作为抽运光时,其频宽越小或晶体厚度越大,频率纠缠度越大;3) 选取合适的脉冲光频宽和晶体厚度可以产生频率负关联、不关联和正关联的双光子。

参考文献

[1] BouwmeesterD, EkertA, ZeilingerA. The physics of quantum information[M]. Berlin, Heidelberg: Springer, 2000: 11- 28.

[2] Giovannetti V, Lloyd S. MacCone L. Advances in quantum metrology[J]. Nature Photonics, 2011, 5(4): 222-229.

[3] D'Angelo M. Kim Y H, Kulik S P, et al. Identifying entanglement using quantum ghost interference and imaging[J]. Physical Review Letters, 2004, 92(23): 233601.

[4] Burnham D C, Weinberg D L. Observation of simultaneity in parametric production of optical photon pairs[J]. Physical Review Letters, 1970, 25(2): 84-87.

[5] Aspect A, Dalibard J, Roger G. Experimental test of Bell's inequalities using time- varying analyzers[J]. Physical Review Letters, 1982, 49(25): 1804-1807.

[6] Kwiat P G, Steinberg A M, Chiao R Y. High-visibility interference in a Bell-inequality experiment for energy and time[J]. Physical Review A, 1993, 47(4): R2472-R2475.

[7] Tapster P R, Rarity J G. Owens P C M. Violation of Bell's Inequality over 4 km of optical fiber[J]. Physical Review Letters, 1994, 73(14): 1923-1926.

[8] 尹娟娟, 俞侃, 包佳祺. 飞秒激光泵浦Ⅰ类BBO晶体中自发参量下转换的研究[J]. 光子学报, 2011, 40(9): 1376-1380.

    Yin J J, Yu K, Bao J Q. Spontaneous parametric down conversion in type-iBBO crystal pumped by femtosecond pulses[J]. Acta Photonica Sinica, 2011, 40(9): 1376-1380.

[9] Kim Y H, Grice W P. Generation of pulsed polarization-entangled two-photon state via temporal and spectral engineering[J]. Journal of Modern Optics, 2002, 49(14/15): 2309-2323.

[10] Kwiat P G, Mattle K, Weinfurter H, et al. New high-intensity source of polarization-entangled photon pairs[J]. Physical Review Letters, 1995, 75(24): 4337-4341.

[11] Kwiat P G, Waks E, White A G, et al. Ultrabright source of polarization-entangled photons[J]. Physical Review A, 1999, 60(2): R773-R776.

[12] Baek S Y, Kim Y H. Spectral properties of entangled photon pairs generated via frequency-degenerate type-I spontaneous parametric down-conversion[J]. Physical Review A, 2008, 77(4): 043807.

[13] 卢宗贵, 刘红军, 景峰, 等. 基于自发参量下转换产生参量荧光的光谱分布特性理论分析[J]. 物理学报, 2009, 58(7): 4689-4696.

    Lu Z G, Liu H J, Jing F, et al. Theoretical analysis of spectral properties of parametric fluorescence via spontaneous parametric down-conversion[J]. Acta Physica Sinica, 2009, 58(7): 4689-4696.

[14] Shi B S, Tomita A. Generation of a pulsed polarization entangled-photon pair via a two-crystal geometry[J]. Physical Review A, 2003, 67(4): 043804.

[15] Bouwmeester D. PanJ W, Daniell M, et al. Observation of three-photon Greenberger-Horne-Zeilinger entanglement[J]. Physical Review Letters, 1999, 82(7): 1345-1349.

[16] Grice W P, Walmsley I A. Spectral information and distinguishability in type-II down-conversion with a broadband pump[J]. Physical Review A, 1997, 56(2): 1627-1634.

[17] Fedorov V, Efremov A, Volkov A, et al. Short-pulse or strong-field breakup processes: a route to study entangled wave packets[J]. Journal of Physics B: Atomic, Molecular and Optical Physics, 2006, 39(13): S467-S483.

[18] Valencia A, Ceré A, Shi X J, et al. Shaping the waveform of entangled photons[J]. Physical Review Letters, 2007, 99(24): 243601.

[19] Ou Z Y. Temporal distinguishability of an N-photon state and its characterization by quantum interference[J]. Physical Review A, 2006, 74(6): 063808.

[20] Keller T E, Rubin M H. Theory of two-photon entanglement for spontaneous parametric down-conversion driven by a narrow pump pulse[J]. Physical Review A, 1997, 56(2): 1534-1541.

[21] Pe ina J, Sergienko A V, Jost B M, et al. . Dispersion in femtosecond entangled two-photon interference[J]. Physical Review A, 1999, 59(3): 2359-2368.

[22] Nasr M B, Di Giuseppe G. Saleh B E A, et al. Generation of high-flux ultra-broadband light by bandwidth amplification in spontaneous parametric down conversion[J]. Optics Communications, 2005, 246(4/5/6): 521-528.

[23] Giovannetti V. MacCone L, Shapiro J H, et al. Extended phase-matching conditions for improved entanglement generation[J]. Physical Review A, 2002, 66(4): 043813.

[24] LutzT, KolenderskiP, JenneweinT. A down-conversion source of positively spectrally correlated and decorrelated telecom photon pairs[C]∥Conference on Lasers and Electro-Optics, June 9-14, 2013, San Jose, California United States. Washington: Optical Society of America.2013: QTh1E. 4.

宋媛媛, 陈鼎, 丛爽. 自发参量下转换制备纠缠光子对的特性[J]. 激光与光电子学进展, 2019, 56(4): 043002. Yuanyuan Song, Ding Chen, Shuang Cong. Property of Entangled Photon Pairs Generated via Spontaneous Parametric Down Conversion[J]. Laser & Optoelectronics Progress, 2019, 56(4): 043002.

引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

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