正式场与物质相互作用的量子理论激光物理研究生课件

1处理激光问题三个层次的理论1．速率方程理论2．半经典理论3．全量子理论2全量子力学方程半经典方程速率方程对泵浦和弛豫过程取平均忽略掉所有的相位关系用来研究激光线宽、强度的起伏、相干性、光子统计等用来研究阈值条件、输出功率等(连续运转、脉冲运转、调Q激光器)用来研究频率牵引和推斥、粒子数的脉动、相位锁定、超短脉冲、相干光学瞬态过程等三个层次理论之间的关系10.1量子力学的三种图象对同一个物理内容，可以存在多种不同的数学描述方式，这些不同的描述方式是完全等价的。量子力学对微观系统状态及其运动规律存在三种等价的描述方式，称之为图像(picture)，或表象，或绘景，它们是：1．Schrödinger图像2．Heisenberg图像3．相互作用(Interaction)图像4在量子力学中，可观测量不是力学量算符和态矢本身，而是力学量的平均值及其概率分布，它们是随时间演化的。如果把力学量平均值和概率分布随时间的演化，全都归之为态矢随时间的演化，而力学量算符不随时间演化，这种描述方式就是Schrödinger图像；反之，全都归之为力学量算符随时间的演化而态矢保持不变，得到Heisenberg图像；部分归之为态矢变化，部分归之为算符变化，则是相互作用图像。510.1.1Schrödinger图像在该图像中，体系的状态矢量|φ(t)>是随时间t演化的，其演化的方式遵守Schrödinger方程而力学量算符不随时间演化：。力学量平均值随时间的演化由态矢来承载：6令其中算符把t0时刻的态|φ(t0)>变换成t1时刻的态|φ(t1)>，称为时间演化算符，它代表一个连续变换（t0和t1任意），把态矢随时间变化而变化用一个变换算符的作用来体现。7由于概率守恒<φ(t1)|φ(t1)>=<φ(t0)|φ(t0)>，且

故由于Hamiltonian算符是厄米算符，由后面的(1.8)式，可以进一步给出8时间演化算符还满足以下性质满足(1.3)式的算符成为幺正算符，它所代表的的变换称为幺正变换（正交变换可看作是一种特殊的幺正变换）9下面令t1=t，t0=0，且采用简写把由有10代入Schrödinger方程，有由于|φ(0)>是任意的，故代表能量算符的Hamiltonian算符不显含t，(1.7)式有以下形式解1110.1.2Heisenberg图像在下面，Schrödinger图像和Heisenberg图像下的态矢和力学量算符分别带有上标S和H。在Heisenberg图像中，力学量平均值随时间的演化，完全归之于力学量算符随时间的演化，而态矢保持不变。对于力学量平均值，有12其中分别是Heisenberg图像下的态矢和力学量算符。不显含时间的Hamiltonian算符，在两种图像下是相等的，这是因为Hamiltonian算符与时间演化算符是对易的。13因此，对Schrödinger图像下的态矢和力学量算符，利用演化算符进行幺正变换，可以得到Heisenberg图像下的态矢和力学量算符。这种幺正变换不改变态矢内积和力学量平均值，不改变算符之间的对易关系，因此不改变物理内容，两种图像等价。假设力学量算符不显含t，即利用(1.9)式，有14我们有利用(1.7)式，即15总之，在Heisenberg图像中，态矢不随时间演化，而力学量算符是随时间演化的，其演化的方式遵守Heisenberg方程。于是，我们得到在Heisenberg图像下，力学量算符随时间演化的Heisenberg方程。16当一个量子系统的Hamiltonian算符可以分解成两部分：10.1.3相互作用图像其主要部分不含时间（通常是自由部分），而微扰部分只对系统产生较小的影响（通常是相互作用部分），这时就可以采用相互作用图像。相互作用图像下的态矢和算符（带上标I），可由Schrödinger图像下的态矢和算符作如下幺正变换得到：17其中的幺正变换算符是由Hamiltonian算符的主要部分来定义的时间演化算符，它同样满足前面给出的演化算符的一切性质。由(1.13)式中的第一式有18如果Hamiltonian算符的主要部分和微扰部分对易，即有则相互作用图像下的态矢又可以表达为19利用(1.13)－(1.15)式以及Schrödinger方程不难验证，在相互作用图像下，态矢和算符分别满足以下方程（算符不显含时间）：由定义(1.13)，相互作用图像下Hamiltonian算符的微扰项与自由项分别为20因此，在相互作用图像下，态矢和算符都随时间演化，其中态矢的演化遵从Schrödinger方程，且由Hamiltonian算符中的相互作用项（微扰项）推动；算符的演化遵从Heisenberg方程，且由Hamiltonian算符中的自由项推动21以上三种图像是对同一物理内容的不同描述方式，在物理本质上是相互等价的。例如，在三种图像中，算符之间的对易关系不会变，算符的平均值不会变，态矢之间的内积不会变，测不准关系不会变，等等。如果对未微扰系统（）已经有充分了解，加上微扰之后，取相互作用图像是合适的，此时算符的运动方程由未微扰系统的Heisenberg方程来描述，它的解是熟悉的、已知的，而态矢量的运动方程只含一个影响较小的微扰算符，便于近似求解。2210.2辐射场与原子的相互作用10.2.1Schrödinger图像下的Hamiltonian算符EbEanħω|b>|a>(n+1)ħω单模辐射场与二能级原子构成的系统考虑由光场和原子共同组成的系统，其中原子是二能级的，上下能级本征态分别是|a>和|b>，分别对应能量本征值Ea=ћωa和Eb=ћωb。23上下能级的本征态矢量|a>和|b>满足正交归一和完备性关系，例如<a|a>=<b|b>=1,<a|b>=<b|a>=0。在以|a>和|b>作为基矢量的表象下，

|a>=1|a>+0|b>，|b>=0|a>+1|b>。将一个矢量用基矢量展开时，展开系数即是该矢量的坐标，由坐标构成的列矩阵，就是矢量的矩阵表示。

因此|a>和|b>在其自身表象下的矩阵表示为24上升算符和下降算符定义上升算符和下降算符如下：25即上升算符把下能级本征态|b>变为上能级本征态|a>，下降算符则把上能级本征态|a>变为下能级本征态|b>。显然上升算符和下降算符互为复共轭转置，即互为厄米共轭。利用和态矢的正交归一性，易证26在原子的量子力学状态|ψ>=Ca|a>+Cb|b>下，上升算符和下降算符的平均值与原子的密度矩阵元对应(书上pp.106-107)：2.系统的Hamiltonian算符由单模光场和二能级原子组成的系统，其总的Hamiltonian算符包括自由光场的贡献(带下标f)、纯原子的贡献(带下标a)，以及光场与原子之间的相互作用的贡献(带下标af，相互作用采用电偶极矩近似)，即有假设单模光场的频率为Ω，二能级原子上下能级本征态分别是|a>和|b>，分别对应能量本征值Ea=ћωa和Eb=ћωb，则有28其中原子的Hamiltonian算符在它自身的表象下，还可以表达为29显然有以下本征方程考虑到原子的固有电偶极矩(平均)为零，即故光场与原子之间相互作用项对应的矩阵为30即其中定义光场与原子之间的耦合系数g为31于是32于是相互作用项又可以表达成考虑到上升算符和下降算符满足即33在上面相互作用包含的四项中

表示原子从上能级跃迁到下能级，同时吸收光场的一个光子；×表示原子从下能级跃迁到上能级，同时吸收光场的一个光子；√表示原子从上能级跃迁到下能级，同时光场增加一个光子；√表示原子从下能级跃迁到上能级，同时光场增加一个光子。×34我们考虑的是由光场和原子构成的孤立系统，无外界作用，为了满足能量守恒，相互作用项只能取为：综上所述，在全量子化理论下，由光场和原子构成的系统的总能量算符为3510.2.2相互作用图像下的相互作用能在上面得到的由光场和原子构成的系统的总哈密顿算符中，前两项对应系统的自由项（定态哈密顿算符），第三项对应相互作用项（微扰项），即有36因此，采用相互作用图象时，变换算符为前面已经讲过，相互作用图象下的态矢和算符可由(1.13)式得到37我们考虑的系统其状态由光场状态和原子状态共同决定（光场状态用单模光子数态|n>描述，原子状态用上下能级本征态|a>和|b>描述）。态矢对应概率振幅，根据同时发生事件的概率相乘原理，可以用|a,n>≡|a>|n>表示原子处于上能级而光场光子数为n的本征态，而用|b,n+1>≡|b>|n+1>表示原子处于下能级而光场光子数为(n+1)的本征态，显然38其中展开系数Ca,n

(t)和Cb,n+1(t)假定为时间的缓变函数。相互作用图像下的相互作用能为在相互作用图像下，利用系统定态哈密顿算符的两个本征态|a,n>和|b,n+1>，可把系统的一般态矢表示为39把(2.16)式代入(2.19)式，可得书上p.188的(10.2.29)式(作业：由(10.2.22)式推出证明(10.2.29)式），即下面的(2.20)式其中，ω0为二能级原子的共振跃迁频率：为便于推导，可对(2.19)的第二式取它的原始表达式，再去掉不满足能量守恒的项。4010.3原子发射和吸收的跃迁几率EbEanħΩ|b>|a>(n+1)ħΩ单模辐射场与二能级原子构成的系统ω0=(Ea-Eb)/ħ=ωa-ωb|n>|n+1>41对于由单模辐射场与一个二能级原子构成的相互作用系统，上节里，我们已经在相互作用图像下，给出了系统态矢的表达式(2.18)和系统Hamiltonian算符中的相互作用项表达式(2.20)，即有（假定光场的初始光子数为n或者(n+1)）10.3.1系统状态随时间演化的方程42在相互作用图像下，系统态矢随时间的演化，是由系统总能量算符中的相互作用项推动的，即有以下方程将(2.18)、(2.20)和(2.21)式代入上式，再利用定态能量算符本征态的正交归一关系，可以求得（书上p.189方程(10.3.5)和(10.3.6)）。43上面的方程，就是系统处于态|a,n>和|b,n+1>的几率振幅随时间变化的方程。10.3.2状态演化方程的几种解使用迭代法求出Ca,n

(t)和Cb,n+1(t)的各阶近似解，也可以求出它们的强信号解，从而确定系统处于态|a,n>和|b,n+1>上的几率。441.受激吸收几率假定系统在初始t=0时刻处在|b,n+1>的状态即初始时刻系统中的原子处于下能级而光场有(n+1)个光子。在t＞0时刻，由于光场与原子之间的相互作用，原子有一定的几率跃迁到上能级而光场相应减少一个光子，即代表该过程的几率振幅Ca,n

(t)≠0for

t＞0.45为了具体地求出t＞0时的几率振幅Ca,n

(t)，对(3.2)式的第一式积分，取一阶微扰近似，有系统从|b,n+1>态跃迁到|a,n>态的几率为46系统受激吸收的几率正比于光场的初始光子数(n+1)，因此若初始光子数为零，原子就不会从下能级跃迁到上能级共振吸收时，Ω=ω0，几率最大。2.受激发射几率和自发发射几率假定系统在初始t=0时刻处在|a,n>的状态，即原子处于上能级而光场有n个光子的状态同理，对(3.2)式的第二式积分，且取一阶微扰近似，有47系统从|a,n>态跃迁到|b,n+1>态的几率为由于初始光场的光子数为n，故上式右端可分为跟初始光场中n个光子相关的部分和跟n个光子无关的部分，即48(3.5)中第一式与初始光场的光子数n成正比，是n个光子诱发的受激发射(stimulatedemission)几率，n=0时该几率为0；第二式跟初始光场的光子数无关，即使初始光场处于光子数为0的真空态，该项仍然存在，故属于自发发射(spontaneousemission)几率。49初始光场的光子数相同的情况下，同一系统的受激发射几率与受激吸收几率相同，都是正比于初始光场的光子数自发发射几率等于一个光子引发的受激发射几率。自发发射在全量子理论中完全是自然出现的，这是因为量子场存在真空起伏，从而存在零点能，零点能对处于上能级的原子产生干扰，使其产生向下能级的跃迁，这个过程中能量守恒（反之不存在“自发吸收”过程，因为违背能量守恒）。而经典电磁场不存在零点能，不能自然地解释自发发射共振发射时，Ω=ω0，所有跃迁过程几率最大。503.自发发射(3.5)式的第二式表达的是一个模式上的自发发射几率。所有模式（频率）上总的自发发射几率，等于对各个模式上的自发发射几率求和（频率连续分布时，对应连续求和，即积分）。体积V内从频率Ω～(Ω+dΩ)之间的光频范围内的模式数为（已考虑到两种偏振状态）51所有模式上总的自发发射几率为52令x=(Ω-ω0)t/2，考虑到被积函数中(sinx/x)2是一个在x=0

(即Ω=ω0)附近的尖峰函数，且被积函数在Ω>>ω0的区域上积分贡献近似为零，因此对被积函数中的Ω2项取近共振近似Ω=ω0并提出积分号之外。考虑到上式第三个等号右端积分的下限通常满足-ω0t/2<<-1，为方便令-ω0t/2=-∞

，且利用公式53最后求得所有模式上总的自发发射几率为把代入(3.6)式，对空间求平均（相当于用1/2代替sin2kz），并且取近共振近似Ω=ω0，得到54（个人体会：在初始时刻t≈0，-ω0t/2≈0，则得到的结果只有(3.6)-(3.8)式结果的一半）单位时间内的总自发发射几率（即自发发射速率）为554.拉比强信号解前面为了计算受激发射和受激吸收几率，在微扰近似下求解(3.2)式，这种近似只适用于相互作用很微弱的情形。当相互作用项很强时，就只能用非微扰方法求解。下面就用非微扰方法求拉比强信号解56对(3.2)式再次对时间求导，并且以原来的(3.2)式代入，可以得到(3.9)式是两个常系数齐次微分方程，其特征方程是：s2+ps+q=0，对于第一个方程，p=i(Ω-ω0)，对于第二个方程，p=-i(Ω-ω0)；对于两个方程都有q=g2(n+1)。先考虑(3.9)式的第一个方程，其特征方程的根是：57于是(3.9)式的第一个方程的通解为58同理，对于(3.9)式的第二个方程，只需把上面特征方程中的p=i(Ω-ω0)换成p=-i(Ω-ω0)，其他求解过程完全一样，总之(3.9)式的两个方程的通解为：(3.2)式把(3.10)式中四个系数A,B,C,D联系起来。例如，把(3.10)式代入(3.2)式的第一个方程，即59可得这样，我们只需要系数A和B来表达解Ca,n

(t)和Cb,n+1(t)。上述解是(3.2)式的严格解析解，适用于任何情形，而前面的微扰解只适用于弱信号情形，这里的解可用于分析拉比强信号解，其中的常数A,B由系统的初始条件确定。当辐射场与原子共振时，Ω=ω0，上述解化简为60于是有61假设系统初始时刻t=0处于下能级、光场有(n+1)个光子，即于是在任意时刻t时，有622.假设系统初始时刻t=0处于上能级、光场有n+1个光子，同理有结论：在强信号作用下，当辐射场与原子发生共振Ω=ω0时，原子总是在上下两个能级之间不断地反复跃迁，跃迁的频率正比于光场数的平方根，即辐射场越强，状态的反转越快。6310.5激光器的库理论本节将讲述激光器的全量子理论。为此需要对激光器建立一个简化的等效模型6465辐射场原子系统损耗机制泵浦过程原子与场相互作用激光系统环境66假定场为单模辐射场假定原子是二能级的且是静止的假定场的振荡频率等于原子的中心频率把激光器中的辐射场作为研究对象，而把其他的一切看作是对场产生影响的背景条件（称作“库”），把库对场产生的影响归结为两类：增益和损耗。因而库可以等效地由位于上能级的激活原子和位于下能级的损耗原子组成，前者模拟激光振荡中的增益机制，其中包含了外界泵浦的作用；后者模拟激光振荡中的损耗机制，把一切损耗都等效为损耗原子对激光的吸收作用。

单模激光场增益饱和非线性，近似到四级微扰激活原子损耗原子|b>|a>增益损耗线性损耗，近似到二级微扰激光器的全量子理论模型68假设原子与辐射场在t=t0时刻开始发生耦合，由原子和辐射场构成的系统的密度算符，满足Schrödinger图像下的运动方程：10.5.1场与物质密度算符的运动方程其中在Schrödinger图像下，原子与辐射场构成的系统其总Hamiltonian算符满足(2.12)式，即69将Schrödinger图像下的密度算符方程变换到相互作用图象下的密度算符方程，有其中相互作用图像下的相互作用能量算符由(2.16)给出，即为其中Ω是光场频率，ω0=ωa-ωb是原子的中心频率，ћωa和ћωb对应原子的上下能级能量。70对(5.1)式两边从初始时刻t=t0开始进行关于时间的积分，便得到以下递推关系例如，把(5.2)式的第二式代入第一式右边，得71这个过程可以重复进行下去，得到这称作“微扰级数展开”，因为只有相互作用能足够小，相当于在系统上附加的一个“微扰”时，展开的无穷级数才收敛。此种求解方法只适用于微扰情形，故被称为“微扰方法”。72对于以上微扰展开，令（为方便去掉上标I）：其中称为第n级微扰（n=0,1,2…)，且有73利用密度算符，可以把力学量算符的平均值表达如下（其中tr代表求迹运算）10.5.2辐射场的约化密度矩阵考虑原子－辐射场总系统的密度算符，设场的能量本征态为|fn>，原子的能量本征态为|am>，则体系总能量本征态为这两个态矢的并矢式：|am,fn>≡|am>|fn>。以{|am,fn>}为基矢组，则一个只与场相关的算符平均值为74其中注意到与场相关的算符只作用于场的本征态，因而可以与原子的本征态互换位置。定义总系统的约化密度算符为75因此，当只计算与场相关的算符平均值时，可以先把总密度算符对原子变量求迹，得到仅与场相关的约化密度算符，使得计算简化（形式上跟只有场存在时的计算是一样的）。若光场和原子两个系统独立，则总系统的密度算符，等于两个系统各自的密度算符之积：此时上述约化密度算符，就等于光场的密度算符。注：如果让总密度算符对应经典概论统计中的联合分布，则约化密度算符对应边缘分布。下面只考虑共振场情形（即Ω=ω0=ωa-ωb）。前面已经表明，在二能级原子的能量本征态|a>、|b>下，利用上升和下降算符，可把相互作用能表达为（共振时它在相互作用图像下和在Schrödinger图像下的表达式是一样的，下面符号表达上已经去掉上标I）10.5.3总密度算符的求解77原子的一般状态|φ>用它的能量本征态|a>和|b>展开为原子的密度算符为78增益原子初始时刻处于上能级，即有|φ(t0)>=|a>，因此初始时刻的密度算符为于是初始时刻场与原子整体系统的密度算符为79于是把(5.7)和(5.9)式代入(5.4)式，可以求出各级近似。首先一级近似为80再把求得的上式(5.10)和(5.9)式继续代入(5.4)式，可以求出二级近似同理，可以给出三级和四级近似，见书上p.196，方程(10.5.15)和(10.5.16)。可看出偶次近似的非对角元为零，奇次近似的对角元为零。以上都是在原子能量本征态下的矩阵表示。81实例：对于初始时刻t0泵浦到上能级|a>的激活原子,假定它的寿命为τ，即平均来说大约经过时间τ之后，原子跃迁到下能级|b>，并发射一个光子到辐射场中，使得原来的辐射场获得了增益。由于辐射场的密度算符（对应原子－辐射场总系统的约化密度算符）描述了光子数占有概率，因此在时间τ内该密度算符的变化反映了辐射场的增益贡献。从t=t0到t=t0+τ时刻，增益使得辐射场密度算符的增量82其中对原子－辐射场系统的总密度算符，计算近似到四级（考虑增益饱和非线性），即取于是把上式代入(5.12)式，得83最后一式中的（t=t0+τ）由书上p.196方程(10.5.16)给出。把书上p.196方程(10.5.13)－(10.5.16)代入上式，可得从t=t0到t=t0+τ时刻，增益使得辐射场密度算符的增量为84反之，对于库中的损耗原子，它们在初始时刻t0处于下能级|b>，即初始时刻原子－辐射场系统的总密度算符为考虑损耗原子很多时的非饱和吸收，计算中只取二级线性近似，同理，从t=t0到t=t0+τ时刻，损耗原子使得辐射场的密度算符（总系统的约化密度算符）变化为85(5.13)和(5.15)两式分别代表了一个增益原子和一个损耗原子所带来的对辐射场的密度算符（即原子－辐射场总系统的约化密度算符）的贡献8610.5.4相互作用图像下约化密度算符的运动方程(5.13)和(5.15)式分别代表了一个增益原子和一个损耗原子所带来的对辐射场的密度算符的贡献，现在研究存在大量增益原子和损耗原子的情形下该密度算符的运动方程。假设λa和λb分别代表单位时间内注入到相互作用区域中的增益原子和损耗原子数量。由于大量原子的平均作用，可以认为在t0~t0+τ这段时间间隔内，场的变化是很小的，从而取近似87令τf表示辐射场开始出现明显变化的时间，则在时间间隔τ<Δt<<τf内，原子－辐射场系统的约化密度算符（纯辐射场的密度算符）的平均变化率，等于增益原子与损耗原子的贡献之和将(5.13)和(5.15)式代入上式右边，并且利用(5.16)式的近似，可得88讨论：以上计算是基于(5.16)式的“粗粒化近似”(coarse-grainedapproximation)，即假定在时间间隔Δt<τ内辐射场基本不变化，而辐射场随时间变化的微分，对应它在时间间隔τ<Δt<<τf内的平均变化，就好像一维时间是以τ为最小单元构成的一样，把场在这个时间单元以内的变化细节通过平均给抹掉。这种近似带来方便，但是不能描述场的某些噪音。8910.5.5光子数表象下辐射场密度矩阵元运动方程用|φ>表示单模辐射场的态矢，|n>(n=0,1,2…)表示辐射场的粒子数本征态，辐射场的密度算符（原子－辐射场系统的约化密度算符）及其在粒子数表象下的矩阵表示为（为方便计，把下标f去掉）显然ρnn代表光场有n个光子的概率90密度矩阵元随时间的变化为把(5.17)式代入上式右边，并且利用(5.13)和(5.15)，以及91得到场的密度矩阵元的运动方程92对角元ρnn代表辐射场中有n个光子的概率，对上面(5.20)式令m=n，得到对角元ρnn随时间的变化为以上给出的是微扰近似解。严格的解析解，即非微扰的强信号解，前人也已经给出，即93当nB/A<<1时，对(5.22)式取近似则(5.22)式变成(5.21)式。在下一节中我们将知道，条件nB/A<<1相当于光子数远小于激光器震荡阀值，因而是弱信号近似。94(5.21)式描述场处于态|n>上的几率变化，它可以重新表达为或95由于ρnn代表场处于态|n>上的几率，(5.23)式右边的正数项，让场处于态|n>的几率增加，表示光场由其他态向态|n>跃迁的过程；负数项让场处于态|n>的几率减小，表示光场由态|n>向其他态跃迁的过程

。λa和λb分别代表单位时间内注入的增益原子和损耗原子数量，由(5.17)式知A,B∝λa,C∝λb，故含系数A、B的项跟增益原子数成正比，含系数C的项与损耗原子数成正比。有了这些准备，我们来分析(5.23)或(5.23)’式右边各项的物理意义。9697系数A跟单位时间内注入的增益原子数λa成正比，(5.23)式第1项和第4项表示在场的作用下，增益原子对场的线性增益作用，其中第1项表示处于态|n-1>的光场，由于增益原子受激发射而跃迁到态|n>，使得场处于态|n>的几率ρnn增加(故第1项为正)，此贡献与场处于初态|n-1>的几率ρn-1,n-1成正比，这是因为增益原子受激发射几率与外场光子数成正比；98第4项表示处于态|n>的光场，由于增益原子受激发射而跃迁到态|n+1>，使得场处于初态|n>的几率ρnn减小(故第4项为负)，此贡献与场处于初态|n>的几率ρnn成正比，这是因为原子受激发射几率与外场光子数成正比。n=0时，第4项不等于零，表示处于真空态|0>的场，由于增益原子自发发射而跃迁到态|1>，因此系数A是与自发发射有关的量。99系数B跟A一样，与增益原子数λa成正比，但是跟含A的线性增益项相反，(5.23)式中含系数B的项是非线性损耗项（跟粒子数平方成正比），这是由原子增益饱和造成的非线性效应。接近增益饱和时，原来处于上能级的增益原子，纷纷跃迁为处于下能级的损耗原子，所以非线性损耗与初始的增益原子数成正比。故系数B反映原子的增益饱和性质。100系数C跟单位时间内注入的损耗原子数λb成正比，(5.23)式第3项和第6项表示在场的作用下，损耗原子对场的线性损耗作用，其中第3项表示处于态|n+1>的光场，由于损耗原子受激吸收而跃迁到态|n>，使得场处于态|n>的几率ρnn增加(故第3项为正)，此贡献与场处于初态|n+1>的几率ρn+1,n+1成正比，这是因为损耗原子受激吸收几率与外场光子数成正比；101第6项表示处于态|n>的光场，由于损耗原子受激吸收而跃迁到态|n-1>，使得场处于态|n>的几率ρnn减小(故第6项为负)，此贡献与场处于初态|n>的几率ρnn成正比，这是因为损耗原子受激吸收几率与外场光子数成正比。n=0时第6项为0，即处于真空态|0>的光场，其零点能不足以引起损耗原子的“自发吸收”。102辐射场的平均光子数及其随时间的变化为把(5.23)式代入(5.25)式右端，并利用(5.24)式，可得103在激光器中，当泵浦速率超过阀值时，腔内光子数急剧增加，使得<n2>远大于<n>，于是上式可简化为由于光强In与平均光子数<n>成正比，上式可与半经典理论中得到的单模场运动方程相比较104即有可见，(A-C)与线性净时间增益系数αn对应，B与自饱和系数βn对应，同时平均光子数运动方程比单模场运动方程多出一项正常数项A，使得<n>=0时，腔内光子数仍然可以随时间增加而增加，即腔内辐射场可以从光子数为零的场建立起来，这是自发发射导致的，故常数A与自发发射相关。因此全量子理论能够揭示激光器的自激振荡性质。10510.6激光的光子统计本节将讨论激光场在阀值以下和阀值以上时的光子统计规律，揭示激光振荡时辐射场从热辐射场演变为相干辐射场的过程与性质型106在(5.23)’式中，前三项是场在态|n-1>和|n>之间跃迁，后三项是场在态|n+1>和|n>之间跃迁，它们共同决定了场处于态|n>上的几率ρnn随时间的变化率。在稳态震荡的条件下，ρnn保持不变：dρnn/dt=0(光场的平均光子数保持不变)，此时场在态|n-1>和|n>之间以及在态|n+1>和|n>之间的跃迁，分别处于动态平衡状态，从而有107对(6.1)式中的第一个方程作变换(n→n+1)则得到第二个方程，因此(6.1)式中的两个方程其实对应同一个方程，它们都是反映光场在两个相邻能态之间的几率流量为0。从中得到以下递推关系从而有108利用归一化条件109备课时，在讲稿本上写上p.200本人给的笔注，以便给学生板书回顾一下有激光光场在稳定状态下，处于各个光子数态|n>的几率不再随时间变化，(6.3)式即为描述此时的激光光场的光子分布的。但要注意，这里只是近似到四级微扰下的解110激光振荡在阀值以下时，单位时间内注入的增益原子少于损耗原子，即λa<λb，此时远未达到饱和状态，故系数B=0，且有10.6.1激光振荡在阀值以下的情况此时(6.3)式变为由于A/C<1，n=0时的几率ρ00最大，而n↑时，ρnn↓，ρnn

随n的变化曲线由图(6.1)中的实线表示111图(6.1)ρnn随n的变化曲线图实线：阀值以下；点划线：阀值；虚线：阀值以上nρnn112此时(6.4)式变为于是(6.6)式变为上式还可以重新写成113由(6.8)式可求得平均光子数为令x=A/C<1，利用114有把(6.10)代入(6.9)式，可得0<A/C<1，令k是Boltzmann常数，T是绝对温度，于是115下面证明，阀值以下的激光场的光子分布是热辐射场的光子分布。热辐射场是混合系综，假设场处于态|n>的几率为pn，则热辐射场的密度算符为其中pn由Boltzmann公式给出116由于场处于态|n>的几率pn又可以写成上式代入(6.13)式，得到密度算符为117在以态|n>为基矢的光子数表象下，热辐射场的平均光子数为118即热辐射场的平均光子数为用上式表示光场占有|n>态的几率，有119对比(6.17)和(6.11)，以及(6.16)和(6.12)，可知工作在阀值以下的激光光场，其光子分布与热辐射场的分布相同，它是若干模式的统计系综分布。10.6.2激光振荡在阀值处的情况在阀