2.3光与物质相互作用的全量子理论在本节,我们将以量子化辐射场与两能级原子的相互作用为例来阐述光与物质相互作用的全量子理论。
2.3.1原子系统与光波场的总哈密顿在半经典理论中,单电子原子与辐射场的相互作用哈密顿为:E r e H H HF A ⋅-+=(2.47)其中A H 和F H 分别代表无相互作用时的原子和辐射场的能量,r 代表电子的位置矢量,E 代表辐射场的振幅。
当辐射场也被量子化后,我们有:ii ii ii A E i E H σ∑∑==(2.48a)∑+=+kk k k F a a H )2/1(ν(2.48b) ∑∑==ji ij ij ji j j i e r e ,,σμ(2.48c) ∑++=kk k k k a a E E )(ε(2.48d)其中+k a 和k a 分别代表光子的产生和湮灭算符,j i ij =σ代表原子跃迁算符,j e ij =μ代表电偶极矩阵元,2/10)2/(V E k k εν =。
于是,我们得到全量子理论中的哈密顿:∑∑∑∑+++++=ji kk k ij ij k iii i kk k k a a g E a a H ,)(σσν(2.49)其中 /)(k k ij ij k E g εμ⋅-=。
在此,我们已从第一项中略去了零点能。
对于一个两能级原子,考虑到ba ab μμ=,我们可令bak ab k k g g g ==,于是方程(2.49)可进一步简化为:∑∑+++++++=kk k ba ab k bb b aa a kk k k a a g E E a a H ))(()(σσσσν(2.50)若我们令bb aa z σσσ-=,ab σσ=+,ba σσ=-,考虑到ω =-b a E E 和1=+bb aa σσ,并略去常数能量因子2/)(b a E E +,则方程(2.50)变为:∑∑+-++++++=kk k k z k k k k a a g a a H ))((21σσωσν(2.51)在上式中,相互作用能由四项组成。
其中-+σk a 项描述原子由上能级跃迁至下能级同时产生一个k 模式光子的过程,+σk a 项描述与其相反的过程;-σk a 项描述原子由上能级跃迁至下能级同时消灭一个k 模式光子的过程,++σk a 项描述与其相反的过程。
注意,在前两个过程中能量是守恒的,但在后两个过程中能量不守恒,因此需将-σk a 项和++σk a 项略去,这相应于半经典理论中的旋转波近似。
于是,我们有:∑∑-++++++=kk k k z k k k k a a g a a H )(21σσωσν(2.52)下面,我们考虑单模量子场与一个两能级原子的相互作用。
略去耦合因子kg 的下脚标后,我们有:10H H H +=(2.53a) 2/0z a a H ωσν +=+(2.53b) )(1-+++=σσa a g H(2.53c)在相互作用图像下求解原子与场的相互作用更为方便,因此下面我们求出相互作用图像下的哈密顿:)(/1/00t i t i t iH t iH I e a ae g e H e H ∆--+∆+-+==σσ(2.54)其中νω-=∆单模场相对于原子跃迁的失谐。
2.3.2几率振幅法和旋转波近似在相互作用图象下,系统状态函数)(t ψ的运动方程为:ψψI H ti =∂∂(2.55)对于一个由两能级原子和单模场组成的系统,若n a ,(n b ,)表示原子处于上(下)能级a (b ),而光波场有n 个光子的状态,则:∑+=nn b n a n b t c n a t c t ],)(,)([)(,,ψ(2.56)相互作用能(2.54)可引起系统在n a ,和n b ,间的跃迁,因此我们要考虑幅度n a c ,和1,+n b c 的演化规律。
将(2.54)和(2.56)代入方程(2.56)可得:1,,.1+∆+-=n b t i n a c e n ig c (2.57a) n a t i n b c e n ig c ,1,.1∆-++-=(2.57b)方程(2.57)与我们在半经典理论中求得的(2.27)非常相似,只是在此我们需要将光波场的状态也考虑进来。
考虑到系统的初始条件,方程(2.57)的一般解为:2/1,,,2sin )0(122sin 2cos )0()(t i n n b n n n n n a n a e t c n ig t i t c t c ∆+⎭⎬⎫⎩⎨⎧⎪⎭⎫ ⎝⎛ΩΩ+-⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛ΩΩ∆-⎪⎭⎫ ⎝⎛Ω=(2.58a)2/,1,1,2sin )0(122sin 2cos )0()(t i n n a n n n n n b n b et c n ig t i t c t c ∆-++⎭⎬⎫⎩⎨⎧⎪⎭⎫ ⎝⎛ΩΩ+-⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛ΩΩ∆+⎪⎭⎫ ⎝⎛Ω= (2.58b)其中,)1(4222++∆=Ωn g n 。
如果原子最初位于上能级,即:)0()0(,n n a c c =,0)0(1,=+n b c ,则:2/,2sin 2cos )0()(t i n n n n n a e t i t c t c ∆⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛ΩΩ∆-⎪⎭⎫ ⎝⎛Ω=(2.59a)2/1,2sin 12)0()(t i n n n n b et n ig c t c ∆-+⎪⎭⎫ ⎝⎛ΩΩ+-= (2.59b)方程(2.58)或(2.59)是一组完整的解,因为关于量子化的场和原子的所有信息均可由它们获得。
在t 时刻单模场中有n 个光子的几率为2,2,)()()(t c t c n p n b n a +=,即:⎪⎭⎫ ⎝⎛Ω⎪⎪⎭⎫ ⎝⎛Ω+⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫ ⎝⎛Ω⎪⎪⎭⎫ ⎝⎛Ω∆+⎪⎭⎫ ⎝⎛Ω=----2sin 4)0(2sin 2cos )0()(122121,1222t n g t t n p n n n n n n n nn ρρ (2.60)其中,)0(nn ρ表示在0=t 时刻单模光波场中有n 个光子的几率,即:!)0()0(2n e n c nnn nn -==ρ (2.61)另一个我们需要关注的量是能级间的反转:[]∑-=nn b n a t c t c t W 2,2,)()()((2.62)由方程(2-58)和(2-62),易得:()∑∞=⎥⎦⎤⎢⎣⎡ΩΩ++Ω∆=02222cos )1(4)0()(n n n n nn t n g t W ρ (2.63)假如单模场的初始态为真空态(0)0(n nn δρ=),则:[]{}t g g gt W 2/1222222)4(cos 441)(+∆+∆+∆=(2.64)显然,即使0)0(n nn δρ=,能级间的反转仍有拉比振荡发生,因为此时有单模场导致的自发辐射发生。
这明显不同于我们在半经典理论中得出的结论:没有驱动场的时候,位于上能级的原子不能跃迁至下能级。
2.3.3自发辐射的Weisskopf-Weigner 理论我们知道,实际上原子的自发辐射是由真空辐射场的众多不同模式导致的,因此有必要考虑模式连续的辐射场与原子的相互作用。
模式连续的辐射场与原子的相互作用哈密顿为:∑+=-+kt i k k I c h e a r g H k .].)([)(0*νωσ(2.65)其中)exp()(00r k i g r g k k ⋅-=,0r 代表原子的位置坐标。
假定在0=t 时刻,原子处于激发态a ,而辐射场处于真空态0,则在t 时刻原子的态矢为:∑+=kk k b a b t c a t c t 1,)(0,)()(,ψ(2.66)其中1)0(=a c ,0)0(=b c 。
由薛定谔方程我们可求得几率幅)(t c a 和)(t c a 的运动方程:∑--=kk b t i k a t c e r g i t c k )()()(,)(0*.νω(2.67a))()()()(0,.t c e r ig t c a t i k k b k νω---=(2.67b)对方程(2.67b )积分并代入方程(2.67a ),我们有:⎰∑---=ta t t i kk a t c e dt r g t c k 0'))(('20.)()()('νω (2.68)假定辐射场的不同模式在频率上紧密相连,则我们可以用积分来代替对k 的求和,即:⎰⎰⎰∑∞→020203sin )2(2dk k d d Vk ππθθφπ (2.69)其中V 是量子体积。
考虑到θμεν22020cos 2)(ab k k Vr g =(θ是ab μ和k ε之间的夹角),我们可得:⎰⎰--∞-=ta t t i k k ab a tc e dtd c t c k 0'))(('033022.)(6)2(4)('νωννεπμ (2.70)考虑到辐射光强主要集中在原子跃迁频率ω附近,我们可以用3ω代替3k ν并将k ν的积分下限扩展至∞-。
此时,)(2)()(2)(0'''))((0'''t c t c t t dt ed t c dt a ta t t i k ta k ππδννω=-=⎰⎰⎰∞∞---(2.71)于是,我们得到:)(2)(.t c t c a a Γ-=(2.72a) )exp()(2t t c a aa Γ-==ρ(2.72b)其中,32303441cabμωπε-=Γ为弛豫常数。
下面,我们计算在自发辐射过程中辐射场的状态。
我们首先计算系数)(,t c k b 。
将)(t c a 的解代入方程(2.67b ),得:⎥⎦⎤⎢⎣⎡Γ+--=-=Γ---Γ---⎰2/)(1)()()(2/)(002/)('0,.''i e r g edt r ig t c k t t i k tt t i k k b k k ωννωνω(2.73a)k kk t t i rk i k t i e eg ba et k 12/)(10,)(2/)(2/∑⎥⎦⎤⎢⎣⎡Γ+--+=Γ---⋅-Γ-ωνψνω(2.73b)下面,我们引入辐射场的状态函数:k kk rk i k i e g 12/)(0∑Γ+-=⋅-ωνγ(2.74)显然,当1-Γ>>t 时,我们有0γψb →,这是具有不同波矢的单光子态的线性叠加。