-
Zhang[25]在真空QED基础上发展出了激光场中的QED。真空QED是描述电磁相互作用的精确理论,但只能处理少体问题,例如一个电子与一个光子的Compton散射。激光中的QED则能够处理电子与大量光子构成的激光束的散射问题。在自然单位制(
$ \hbar = c = 1 $ )下,圆偏振的激光场能够用经典矢势描述:$$ {\boldsymbol{A}}({\boldsymbol{x}},\,{t}) = A\big\{{\boldsymbol{x}}_0\cos[k(z-t)]+{\boldsymbol{y}}_0\sin[k(z-t)]\big\}\text{。} $$ (1) 一个沿
$ -z $ 轴前进的电子与单位体积的上述沿$ z $ 轴传播的激光正碰过程中自发辐射一个光子$({\boldsymbol{k}}',\,{\boldsymbol{e}}_i')$ 至立体角$ {\rm{d}}\varOmega_{k'} $ 内的微分散射截面为[25]$$\begin{array}{l} \dfrac{{\rm{d}}\varSigma}{{\rm{d}}\varOmega_{k'}} ({\boldsymbol{p}} , \,\sigma , \,{\boldsymbol{k}};{\boldsymbol{p}}' , \,\sigma' , \,{\boldsymbol{k}}' , \,{\boldsymbol{e}}_i') = \\ \ \ \ \ \ \dfrac{\alpha k'^2}{8{\text{π}} m{\cal{N}} k|p_z|(E-p_z)(E+m)(E'+m)} \times \\ \ \ \ \ \ \displaystyle \sum\limits_{{\cal{N}} = -\infty}^{\infty}\Bigg|\displaystyle \sum\limits_{\nu = 0 , \,\pm1}{{J}}_{{\cal{N}}-\nu}(p_{\perp}'R')\times \\ \ \ \ \ \ \Big[\delta_{\sigma , \,\sigma'}F_i^{(\nu)}{\rm{e}}^{-{\rm{i}}{\cal{N}}\varphi_{k'}}+ \delta_{\sigma , \,-\sigma'}G_i^{(\nu)}{\rm{e}}^{{\rm{i}}(\sigma-{\cal{N}})\varphi_{k'}}\Big]\Bigg|^2, \end{array}$$ (2) 其中:
$\sigma,\,\sigma'$ 分别为表征初末态电子自旋的量子数,取值均为$ \pm 1 $ ;$E,\,{\boldsymbol p} = p_z{\boldsymbol z}_0$ 分别为初始电子的能量、动量;$ E^2 = p_z^2+m^2 $ ;$E',\,{\boldsymbol p}'$ 分别是末态电子的能量、动量;$ E'^2 = p_z'^2+p_{\perp}'^2+m^2 $ ;${\boldsymbol{k}},\, {\boldsymbol{k}}'$ 分别为初末态光子的波矢;$ {\boldsymbol{e}}_i' $ 是垂直于波矢$ {\boldsymbol{k}}' $ 的极化矢量;$ R' = \frac{eA}{k(E'-p_z')} $ ;量子数$ {\cal{N}} $ 可以理解为激光场中参与Compton背散射的光子数。系数$F_i^{(\nu)}, \,G_i^{(\nu)}$ 详见文献[25]。2015年KEK在激光电子Compton背散射实验中,244个电子束团与激光束相互作用的时间内,探测到约1 756.8个能量为29 keV的光子[12]。利用式(2)中的微分散射截面可以计算该实验中的散射光子数目
$ N_{\rm{ph}} $ :$$ N_{\rm{ph}} = \int{\rm{d}}\varOmega_{k'}\boldsymbol\cdot \overline{\cfrac{{\rm{d}}\varSigma}{{\rm{d}}\varOmega_{k'}}}\boldsymbol\cdot t \boldsymbol\cdot N_{\rm{e}} \approx 1.55\times 10^9, $$ (3) 其中
$ \overline{\frac{{\rm{d}}\varSigma}{{\rm{d}}\varOmega_{k'}}} $ 为对任意极化方向的平均微分散射截面,即需要将式(2)中的微分散射截面对电子和光子的极化方向末态求和、初态求平均;$ {\rm{d}}\varOmega_{k'} $ 为探测器覆盖的立体角;$ t $ 为244个电子束团与激光束的相互作用时间;$ N_{\rm{e}} = $ 每个电子束团的电荷量$ \times 244/e $ ,为244个电子束团中的电子数目。同样地,利用实验室系下的KN公式[25]
$$ \begin{array}{l} \cfrac{{\rm{d}}\sigma}{{\rm{d}}\varOmega_{k'}}({\boldsymbol{p}} , \,{\boldsymbol{k}} , \, {\boldsymbol{e}}_i;\, {\boldsymbol{p}}' , \, {\boldsymbol{k}}' , \,{\boldsymbol{e}}_i') = \cfrac{\alpha {k'}^2}{4k^2|p_z|(E-p_z)}\ \,\boldsymbol\cdot \\ \left[\cfrac{(E-p_z)k}{(E-p_z\cos\theta)k'}+\cfrac{(E-p_z\cos\theta)k'}{(E-p_z)k}-2+4\big|{\boldsymbol{e}}\boldsymbol\cdot {\boldsymbol{e}}_i'\big|^2\right] , \, \end{array} $$ (4) 也可计算该实验中的辐射光子数目
$ N_{\rm{ph}} $ :$$ N_{\rm{ph}} = \int {\rm{d}}\varOmega_{k'} \boldsymbol\cdot\overline{\cfrac{{\rm{d}}\sigma}{{\rm{d}}\varOmega_{k'}}}\boldsymbol\cdot \rho_{\nu0}\boldsymbol\cdot t \boldsymbol\cdot N_{\rm{e}} \approx 1.79\times 10^9\text{,} $$ (5) 其中:
$ \theta $ 是光子的散射角度;$ \rho_{\nu 0} $ 为单位体积的激光中的光子数目。能够看到,由式(2)计算出的散射光子数目
$ 1.55\times10^9 $ 与由KN公式算出的$ 1.79\times10^9 $ 接近,但两者均比实验中实际探测到的光子数目1 756.8大很多。这是由于激光电子散射理论[25]假定的入射激光为完全相干的零温激光,而实验中的激光并非如此理想的情形,无法达到理想的激光强度。同时,在计算激光中光子数密度$ \rho_{\nu 0} $ 时同样也将激光束理想化,而实际的激光具有一定的展宽,因此实际的入射光子数密度比理论值低。实验激光束在垂直于传播方向的平面上强度具有空间分布的效应,该效应也会导致这一差异。另外X光子、$ \gamma $ 光子穿透性很强,数目不易精确测量致使探测器存在一定的探测效率。这些都会造成辐射光子数目的理论值高于实际测量值。 -
激光场中的QED[25]提供了一个精确描述激光电子相互作用的理论框架。但是文献[25]中只计算了理想情形,即一束理想的电子束和理想的平面波激光正碰。然而实际的激光束有一定展宽、电子束有能散,激光器中也存在着不同光子模式的竞争。如何研究更接近实际情形的激光电子散射过程?
激光电子Compton散射模型中的微分散射截面式(2)给出了从确定系统初态
$ (p, \,\sigma, \,k) $ 到确定末态$ (p', \,\sigma',\, k') $ 的跃迁概率;实际情形的激光电子散射若要形成激光,则需要计算不同初态跃迁到某一确定光子末态的跃迁概率,也就是计算不同频率的激光和不同能量的电子,适当组合散射得到单一电磁模式光子的跃迁概率。这就需要引入激光电子初态态密度,将初态电子和初态光子的能量分布考虑进来。激光电子正碰过程中的能动量守恒分别为[25]
$$ E'-E+\frac{e^2 A^2}{2}\left(\frac{1}{E'-p_z'}-\frac{1}{E-p_z}\right)+k'-{\cal{N}}k = 0, $$ (6) $$ {\boldsymbol{p}}'-{\boldsymbol{p}}+\frac{e^2 A^2}{2k}\left(\frac{1}{E'-p_z'}-\frac{1}{E-p_z}\right){\boldsymbol{k}}+{\boldsymbol{k}}'-{\cal{N}}{\boldsymbol{k}} = 0\text{。} $$ (7) 从散射的角度看,初态的电子光子发生Compton散射得到末态的电子光子;从能级跃迁的角度,可以看作自由电子在不同能级间跃迁,进而辐射出光子。如果只有单个的自发跃迁过程发生便是自发辐射。自发辐射产生的末态光子的存在会诱发受激辐射,受激辐射将使得末态同一状态的光子数加倍。因此单位时间内受激辐射过程带来的同一末态光子数目增量将正比于通过自发辐射产生的同一末态光子数。
在激光电子Compton散射过程中,从某一确定系统初态自发跃迁到任意末态的跃迁速率
$ \varGamma_{\rm{c}} $ 可写为$$ \varGamma_{\rm{c}} = (2{\text{π}}/\hbar)|\langle f|H'|i\rangle|^2\rho_f = \overline{\cfrac{{\rm{d}}\varSigma}{{\rm{d}}\varOmega_{k'}}}\,{\rm{d}}\varOmega_{k'}\text{;} $$ 而由任意系统初态自发跃迁得到同一光子末态的跃迁速率
$ \varGamma_{\rm{d}} $ 为[26]$$ \begin{split} \varGamma_{\rm{d}} &= (2{\text{π}}/\hbar)|\langle f|H'|i\rangle|^2\rho_i \\& = \widetilde{\cfrac{{\rm{d}}\varSigma}{{\rm{d}}\varOmega_{k'}}}\,{\rm{d}}\varOmega_{k'}\rho_i/\rho_f , \end{split}$$ (8) 其中
$ \widetilde{\frac{{\rm{d}}\varSigma}{{\rm{d}}\varOmega_{k'}}} $ 为将式(2)中的微分散射截面对激光束电子束的极化方向初态求和、末态求平均;$ \rho_f $ 是末态态密度;初态态密度$ \rho_i $ 为激光电子散射系统的守恒能量附近的一个小的能量区间$$ E_s = E-E'+\frac{e^2 A^2}{2}\left(\frac{1}{E-p_z}-\frac{1}{E'-p_z'}\right)-k'+{\cal{N}}k $$ 中包含的能够导致所需末态并且满足能动量守恒条件的所有可能的不同态的组合数目[26]:
$$ \rho_i = \int {\rm{d}}E\ \rho_e(E)\boldsymbol\cdot\frac{{\rm{d}}k}{{\rm{d}}E_s}\boldsymbol\cdot\rho_{\nu}(k)\text{。} $$ (9) 为了得到
$ {\rm{d}}k/{\rm{d}}E_s $ ,需要利用动量守恒条件(7)将$ E_s $ 写成初始电子光子的能量动量、末态光子能量的函数。即$$ E_s = -\dfrac{e^2A^2}{2(p_z+{\cal{N}}k-k_z'-p_z')+\frac{e^2A^2}{E-p_z}}-k'+E-p_z+k_z' , $$ (10) 其中
$$ \begin{split} & p_z' = \dfrac{{ \mathbb{K}}(2k_{\perp}'^2+2m^2+e^2A^2)+ e^2A^2 \sqrt{\mathbb{K}^2+k_{\perp}'^2+m^2+e^2A^2}}{2(k_{\perp}'^2+m^2+e^2A^2)},\\& \mathit{ \mathbb{K}} = \frac{e^2A^2}{2(E-p_z)}+p_z-k_z'+{\cal{N}}k, \\[-10pt] \end{split}$$ (11) 那么便能够得到
$$ \frac{{\rm{d}}k}{{\rm{d}}E_s} = \frac{1}{{\cal{N}}}\boldsymbol\cdot \frac{\sqrt{\mathbb{K}^2+k_{\perp}'^2+m^2+e^2A^2}}{\sqrt{\mathbb{K}^2+k_{\perp}'^2+m^2+e^2A^2}+\mathbb{K}}\text{。} $$ (12) 注意到
$ {\mathbb{A}}^2 = k_{\perp}'^2+m^2+e^2A^2 $ ,$\dfrac{{\rm{d}}k}{{\rm{d}}E_s}$ 可以写成$$ \frac{{\rm{d}}k}{{\rm{d}}E_s} = \frac{1}{{\cal{N}}}\boldsymbol\cdot \frac{1}{\mathbb{A}^2}(\mathbb{K}^2+\mathbb{A}^2-\mathbb{K}\sqrt{\mathbb{K}^2+\mathbb{A}^2})\text{。} $$ 当背景激光场的强度低于
$ 10^{19}\; {\rm{W/m^2}} $ 时,初始电子能量不低于几个MeV时,有$ \mathbb{K}^2\gg \mathbb{A}^2 $ 成立,那么$$ \frac{{\rm{d}}k}{{\rm{d}}E_s} \approx \frac{1}{2{\cal{N}}}\text{。} $$ (13) 假设初始电子、光子能量均服从Gauss分布
$$ \rho_e(E) = (1/{\text{π}}^{1/2})(1\big/\varDelta E)\mathrm{exp}\big[-(E-E_0)^2\big/\varDelta E^2\big] , $$ (14) $$ \rho_\nu(k) = (1/{\text{π}}^{1/2})(1\big/\varDelta k)\mathrm{exp}\big[-(k-k_0)^2\big/\varDelta k^2\big] , $$ (15) 激光电子系统的初态态密度可写为
$$ \begin{split} \rho_i & = \dfrac{1}{2{\cal{N}}\boldsymbol\cdot {\text{π}} \boldsymbol\cdot\varDelta E \varDelta k}\times \\& \ \ \ \displaystyle \int {\rm{d}}E\ \mathrm{exp}\left[-\dfrac{(E-E_0)^2}{\varDelta E^2}\right]\mathrm{exp}\left[-\dfrac{(k-k_0)^2}{\varDelta k^2}\right]\text{。} \end{split}$$ (16) 需要注意,这里初态态密度
$ \rho_i $ 中的初始电子和初始光子的能量已不再是两个相互独立的变量,两者将通过能动量守恒条件(6),(7)及对末态光子能量的选取互相制约。因此,我们需将式(16)中初始光子能量$ k $ 写作初始电子能量$ E $ 和末态光子能量$ k' $ 的函数:$$ k = \cfrac{1}{{\cal{N}}}\boldsymbol\cdot \dfrac{\left(E-p_z\cos\theta \right)k'+\frac{e^2A^2}{2(E-p_z)}\left(1-\cos\theta\right)k'}{E-p_z-(1-\cos\theta)k'}\text{。} $$ (17) 为简化计算,我们只考虑每个散射过程中激光束中只有一个光子被散射(或说被吸收)的过程[27],即
$ {\cal{N}} = 1 $ 。另外,由前灯效应可知实验室系下散射的光子集中在电子前进方向上的一个窄圆锥内。若散射光子方向沿着入射电子方向与其相伴而行则更有利于受激辐射的发生,因此取$ \theta = {\text{π}} \ (即p_\perp' = k_\perp' = 0) $ ,那么有$$ k = \dfrac{(m^2+e^2A^2)k'}{(E-p_z)^2-2(E-p_z)k'}, $$ 其中确定末态光子的能量
$ k' $ 将取决于式(14)与式(15)中初始电子光子的中心能量$E_0, k_0$ 。那么单位能量区间内的初始态密度为$$ \begin{array}{l} \rho_i = \dfrac{1}{2 \text{π} \boldsymbol\cdot\varDelta E \varDelta k}\times \\ \displaystyle \int {\rm{d}}E\ \mathrm{exp}\left \{-\dfrac{(E-E_0)^2}{\varDelta E^2}-\dfrac{\left[ \frac{(m^2+e^2A^2)k'}{(E-p_z)^2-2(E-p_z)k'}-k_0\right]^2}{\varDelta k^2}\right\}\text{。} \end{array} $$ (18) 利用两体Feynman公式[28]计算激光电子系统单位能量区间单位体积的末态态密度:
$$ \begin{split} \rho_f &= (2 \text{π})^{-3}\boldsymbol\cdot k'E'\cfrac{|{\boldsymbol{k}}'|^3 {\rm{d}}\varOmega_{k'}}{(E+k){\boldsymbol{k}}'\boldsymbol\cdot {\boldsymbol{k}}'-k'({\boldsymbol{k}}'+{\boldsymbol{p}}')\boldsymbol\cdot{\boldsymbol{k}}'} \\ & = (2\text{π})^{-3}\boldsymbol\cdot \cfrac{E'k'^2 {\rm{d}}\varOmega_{k'}}{E+k-k'-p'}, \end{split} $$ (19) 那么到确定光子末态的跃迁速率(8)将有如下的形式:
$$ \begin{split} \varGamma_{\rm{d}}= & \dfrac{\alpha {\text{π}}^2}{m|p_z|E'(E+m)(E'+m)}\dfrac{E+k-k'-p'}{k(E-p_z)}\,\boldsymbol\cdot \\& \ \ \ \ \left[\dfrac{eA(p_z-E-m)(p_z'-E'-m)}{2}\right]^2 \,\boldsymbol\cdot \\& \ \ \ \ \left[\left(\dfrac{1}{E-p_z}\right)^2+\left(\cfrac{1}{E'-p_z'}\right)^2\right]\boldsymbol\cdot \rho_i \,, \end{split} $$ (20) 其中初态态密度
$ \rho_i $ 将通过数值积分得到。 -
一定长度上的增益定义为平均每个受激光子穿过该长度的过程中系统中因受激辐射额外散射出的光子数目[26]。受激辐射会使末态光子数目加倍,因此该过程的增益可通过自发辐射产生的光子数目计算。
假设在某一位置
$ x $ 处单位体积内的受激光子数目为$ n $ ,一个位于$ x $ 处的受激光子在一段时间后来到位置$ x+{\rm{d}}x $ 处,这时单位体积内的受激光子数为$ n+{\rm{d}}n $ 。那么在这段时间$ {\rm{d}}x $ 内平均1个受激光子带来的额外光子数目为$$ \frac{{\rm{d}}n}{n} = \varGamma_{\rm{d}}\boldsymbol\cdot{\rm{d}}x\boldsymbol\cdot \rho^{}_{e0}, $$ (21) 其中
$\rho^{}_{e0}$ 是单位体积的电子数,那么$$ n = n^{}_0\exp(\varGamma_{\rm{d}}\rho^{}_{e0}x)\text{。} $$ (22) 激光波荡器FEL中平均1个受激光子穿过相互作用长度
$ L $ 上的增益为$$ G = 10\lg(n/n^{}_0) = 10\lg {\rm{e}}\boldsymbol\cdot \varGamma_{\rm{d}}\rho^{}_{e0}\boldsymbol\cdot L , $$ (23) 这里的增益
$ G $ 为最大可能增益,未考虑受激吸收的影响[26]。 -
2015年KEK在激光电子Compton背散射实验中测到了X光子[12]。2019年Hawaii大学实验组利用FEL作为入射激光光源在Compton背散射实验中同样测到了X光子[13]。2012年KEK利用电子与激光束的散射得到了
$ \gamma $ 光子[14-15]。但是这些实验中出射光都没能形成相干光。如果按照这些实验中的参数搭建激光波荡器自由电子激光器,我们能够计算出其中的最大可能增益$ G $ ,列于表1的前三行。为了比较,我们同时计算了世界上第一台磁波荡器X射线自由电子激光器[1](XFEL)中的最大可能增益,见表1中的第四行。其中第四行第四列等效的入射光波长为表 1 各种类型波荡器FEL中的增益
波荡器类型 $E_0$/MeV $\varDelta E/E_0$ $\lambda_0$/nm $\varDelta k/k_0$ $k'$/keV $\rho^{}_{e0}/{\rm{cm^{-3} } }$ $ G $ $\widetilde{G}$ Ref. 激光 $ 40.0 $ $10^{-3}$ 1 047 $0.01$ $29.0$ $2.65\times10^{12}$ $1.78\times10^{-13}$ $2.61\times 10^{-15}$ [12] 激光 $42.2$ $10^{-3}$ 3 100 $0.01$ $10.9$ $2.47\times10^{14}$ $8.33\times10^{-9}$ $3.62\times 10^{-10}$ [13] 激光 $1.28\times10^3$ $10^{-3}$ 1 032 $0.01$ $2.95\times10^4$ $9.65\times10^{14}$ $5.26\times 10^{-17}$ $7.64\times 10^{-19}$ [14] 磁 $13.6\times10^3$ $10^{-4}$ ${4.28\times 10^8}^{\ast}$ $0.10^{\dagger}$ $8.22$ $7.95\times10^{17}$ $1.01\times10^{4}$ $6.04\times 10^{7}$ [1] 射频 $4.5$ $10^{-3}$ $3\times10^6$ $0.08$ $ 0.13\times 10^{-3} $ $8.00\times10^{17}$ $1.74\times10^6$ $ 7.32\times10^{7}$ − *磁波荡器中等效入射光波长由磁波荡器的参数利用式(24)计算得到;†磁波荡器中入射光的展宽为估计值。 $$ \lambda_0 = (1+\beta)\lambda_{\rm{u}}(1+K^2/2), $$ (24) 其中:
$ \beta $ 是电子的速度;$ K = 3.5 $ 是无量纲的波荡器参数;$ \lambda_{\rm{u}} = $ 3 cm是波荡器周期[1]。表1中前三行计算的是由式(1)描述的一束沿
$ z $ 轴传播的电磁波与沿$ -z $ 轴运动的电子束正碰产生的最大可能增益,这里的电磁波是一束行波。永久磁铁构造的磁波荡器也可以看成电磁场,但一动一静。由于目前实验上还没有通过行波电磁场与电子散射成功得到激光,因此我们把已经成功得到激光的磁波荡器XFEL中的最大可能增益作为标准,探讨用行波电磁场和电子束正碰得到相干光的可能性。表1中的
$ G $ 为相互作用长度$ L = 1 $ cm时的最大可能增益。实际上,第一台XFEL中经过60 m的长度之后才产生了相干X光[1],60 m相当于约140个等效入射光波长。所以我们也计算了经过$ L = 140 $ 个入射光波长的最大可能增益(表1中的$ \widetilde{G} $ ),即对于电子分别经过0.15, 0.44, 0.14 mm、60 m长度上的增益。表1中计算前三行的激光波荡器FEL的增益时用到的实验参数列于表2中,计算第四行磁波荡器XFEL中的增益所用参数见表3。根据表2中所列的入射激光的功率与焦斑尺寸可计算激光强度$ I $ ,其与式(1)描述激光束的相干振幅$ A $ 之间有如下关系[29]:$$ \frac{eA}{mc} = \sqrt{\frac{\alpha \lambda_c \lambda^2 I}{\pi mc^3}}, $$ (25) 其中
$ \lambda $ 为激光波长,$ \lambda_{\rm c} $ 为电子的Compton波长。表 2 计算中所用的激光电子散射实验参数
实验 入射电子能量 每电子束团电荷量 入射激光波长 入射激光功率 激光焦斑尺寸 电子束尺寸 碰撞持续时间 2015 KEK[12] 40.00 MeV 61.5 pC 1 047 nm $4.1$ kW 直径$10$ μm 直径$10$ μm 1.5 μs 2019 Hawaii[13] $42.20$ MeV $50.0$ pC $3.1 $ μm $0.3 $ MW 直径$30$ μm $42~\mu {\rm m}\times 64$ μm 2.0 ps 2012 KEK[14] 1.28 GeV $0.5\times 10^{10}$e 1 032 nm 35.0 kW $26~\mu {\rm m}\times 38$ μm $110~\mu {\rm m}\times10$ μm 0.2 ms 表 3 计算中所用第一台XFEL的实验参数[1]
磁波荡器周期 磁波荡器长度 磁感应强度峰值 波荡器参数 入射电子能量 电流峰值 电子束尺寸 3 cm 132 m 1.25 T 3.5 13.6 GeV 3 kA 直径$10 $ μm 由表1可以看到,经过140个入射光波长长度后激光波荡器FEL中的增益远低于第一台XFEL中的增益。可以说,以上述的激光电子散射的实验参数设置难以获得较大数量的单一模式受激光子。这是否意味着利用激光电子正碰无法得到相干光呢?我们扩大参数空间寻找可能的高增益参数组合,找到了能得到较高增益的电磁(射频)波荡器FEL参数如表1最后一行所列。这时140个入射电磁波波长长度上的增益不低于第一台XFEL中的增益。也就是说利用能量为4.5 MeV,电流为3 kA,束流直径10 μm的电子束与波长为3 mm,功率为300 kW,束流直径为3.5 mm的毫米波正碰,那么经过42 cm相互作用长度之后有可能得到波长为9.73 μm的相干光。
我们注意到,这组参数中的入射光波长已经是射频波段。一方面,不论是磁波荡器、电磁波荡器(射频波、激光用作波荡器),物理上都可以转化为Compton或Thomson散射过程处理计算[9, 11]。相干射频波作为电磁波,也能够用式(1)描述,其与电子的相互作用,也是可以利用本文中的理论框架计算的。另一方面,从我们计算的激光电子散射中的最大可能增益看来,我们仍无法断言光波波段激光与电子散射一定能够作为相干光源;但从找到的高增益参数看,利用射频波电子散射有可能产生相干光。
近几年来,美国SLAC实验室一直在研究利用振荡电磁场产生短周期波荡器替换磁波荡器搭建小型FEL[11, 30]。2014年该实验室利用等效周期为13.9 mm的毫米波波荡器实现了种子型相干谐波产生[11],但未能得到相干光;2019年该实验室设计了周期为1.75 mm,最小波束孔径为2.375 mm的射频波荡器,波荡器参数
$ K = 0.1 $ 时所需功率为1.4 MW,对于铜材质的波荡器腔,该功率下可运行250 ns[30]。该射频波荡器周期非常接近我们找到的高增益光波荡器的周期。我们的计算表明,能量为4 MeV,电流强度为3 kA,束流直径10 μm的电子束,若与SLAC实验室的波长为1.75 mm,波束孔径2.375 mm,功率1.4 MW的射频波正碰, 当相互作用长度达24.5 cm后,则可能形成波长为7.20 μm的相干光。
-
摘要: 借鉴自由电子激光(FEL)发展之初Madey对自由电子激光器中受激辐射引起的增益的讨论,通过在激光场中的量子电动力学(QED)的模型中引入激光电子系统初态态密度以及由不确定的系统初态到确定光子末态的跃迁速率,推导了激光电子正碰过程中受激辐射至单一电磁模式产生的最大可能增益。采用成功得到X射线或
$\gamma $ 射线光子的三个激光电子Compton背散射实验的实验参数计算了激光电子散射过程中的最大可能增益,与第一台X射线自由电子激光(XFEL)中的最大可能增益作比较,进而对激光电子散射作为激光光源的可行性进行评估。计算结果表明,现有的能够得到X射线光子或$\gamma $ 射线光子的激光电子散射实验中的最大可能增益远低于第一台XFEL中的。本工作未能找到合适的激光电子参数以获得比第一台XFEL中更高的最大可能增益,但是在入射电磁波位于射频波段范围内找到了能够实现较高增益的参数组合。Abstract: Following Madey’s calculations on the possible gain induced by stimulated radiation in Free Electron Laser(FEL), we derive the possible maximum gain in laser electron scattering basing on Quantum Electrodynamic(QED) in laser field via introducing the initial state density of laser electron system and the transition rate into a definite final photon state in laser electron scattering. The parameters of laser electron scattering experiments where X-ray or$ \gamma $ -ray photons were succesfully obtained are used to calculate the possible maximum gain in laser electron scattering, and the results are compared with that in the first lasing X-ray FEL(XFEL). The calculations show that if laser undulator FEL were built according to the existing experimental conditions in those laser electrons scattering experiments, the possible maximum gain would be much lower than that in the first lasing XFEL. While we have found appropriate parameters to achieve a relatively high gain when the wavelength of incident light is in radio frequency range.-
Key words:
- gain /
- laser electron scattering /
- stimulated emission /
- laser undulator /
- free electron laser
-
表 1 各种类型波荡器FEL中的增益
波荡器类型 $E_0$/MeV $\varDelta E/E_0$ $\lambda_0$/nm $\varDelta k/k_0$ $k'$/keV $\rho^{}_{e0}/{\rm{cm^{-3} } }$ $ G $ $\widetilde{G}$ Ref. 激光 $ 40.0 $ $10^{-3}$ 1 047 $0.01$ $29.0$ $2.65\times10^{12}$ $1.78\times10^{-13}$ $2.61\times 10^{-15}$ [12] 激光 $42.2$ $10^{-3}$ 3 100 $0.01$ $10.9$ $2.47\times10^{14}$ $8.33\times10^{-9}$ $3.62\times 10^{-10}$ [13] 激光 $1.28\times10^3$ $10^{-3}$ 1 032 $0.01$ $2.95\times10^4$ $9.65\times10^{14}$ $5.26\times 10^{-17}$ $7.64\times 10^{-19}$ [14] 磁 $13.6\times10^3$ $10^{-4}$ ${4.28\times 10^8}^{\ast}$ $0.10^{\dagger}$ $8.22$ $7.95\times10^{17}$ $1.01\times10^{4}$ $6.04\times 10^{7}$ [1] 射频 $4.5$ $10^{-3}$ $3\times10^6$ $0.08$ $ 0.13\times 10^{-3} $ $8.00\times10^{17}$ $1.74\times10^6$ $ 7.32\times10^{7}$ − *磁波荡器中等效入射光波长由磁波荡器的参数利用式(24)计算得到;†磁波荡器中入射光的展宽为估计值。 表 2 计算中所用的激光电子散射实验参数
实验 入射电子能量 每电子束团电荷量 入射激光波长 入射激光功率 激光焦斑尺寸 电子束尺寸 碰撞持续时间 2015 KEK[12] 40.00 MeV 61.5 pC 1 047 nm $4.1$ kW 直径$10$ μm 直径$10$ μm 1.5 μs 2019 Hawaii[13] $42.20$ MeV $50.0$ pC $3.1 $ μm $0.3 $ MW 直径$30$ μm $42~\mu {\rm m}\times 64$ μm 2.0 ps 2012 KEK[14] 1.28 GeV $0.5\times 10^{10}$e 1 032 nm 35.0 kW $26~\mu {\rm m}\times 38$ μm $110~\mu {\rm m}\times10$ μm 0.2 ms 表 3 计算中所用第一台XFEL的实验参数[1]
磁波荡器周期 磁波荡器长度 磁感应强度峰值 波荡器参数 入射电子能量 电流峰值 电子束尺寸 3 cm 132 m 1.25 T 3.5 13.6 GeV 3 kA 直径$10 $ μm -
[1] EMMA P, AKRE R, ARTHUR J, et al. Nature Photonics, 2010, 4(9): 641. doi: 10.1038/nphoton.2010.176 [2] ISHIKAWA T, AOYAGI H, ASAKA T, et al. Nature Photonics, 2012, 6(8): 540. doi: 10.1038/nphoton.2012.141 [3] KANG H S, MIN C K, HEO H, et al. Nature Photonics, 2017, 11(11): 708. doi: 10.1038/s41566-017-0029-8 [4] GINZBURG V L. Rev Mod Phys, 2004, 76: 981. doi: 10.1103/RevModPhys.76.981 [5] VALI V, VALI W. Proceedings of the IEEE, 1963, 51(1): 182. doi: 10.1109/PROC.1963.1677 [6] BALDWIN G C, SOLEM J C. Rev Mod Phys, 1997, 69: 1085. doi: 10.1103/RevModPhys.69.1085 [7] BERTOLOTTI M, SIBILIA C. Applied Physics, 1979, 19(1): 127. doi: 10.1007/BF00900548 [8] AVETISSIAN H K, AVETISSIAN A K, MKRTCHIAN G F. Phys Rev Lett, 2014, 113: 023904. doi: 10.1103/PhysRevLett.113.023904 [9] DI PALMA E, CECCUZZI S, RAVERA G L, et al. Applied Sciences, 2021, 11(20): 9499. [10] SHINTAKE T. Experimental Results of Microwave Undulator[C]//LINDAU I E, TATCHYN R O. Insertion Devices for Synchrotron Sources: Volume 0582. SPIE, 1986: 336. [11] TANTAWI S, SHUMAIL M, NEILSON J, et al. Phys Rev Lett, 2014, 112: 164802. doi: 10.1103/PhysRevLett.112.164802 [12] SHIMIZU H, AKEMOTO M, ARAI Y, et al. Nucl Instr and Meth A, 2015, 772: 26. [13] NIKNEJADI P, KOWALCZYK J M D, HADMACK M R, et al. Phys Rev Accel Beams, 2019, 22: 040704. doi: 10.1103/PhysRevAccelBeams.22.040704 [14] BONIS J, CHICHE R, CIZERON R, et al. Journal of Instrumentation, 2012, 7(01): P01017. doi: 10.1088/1748-0221/7/01/p01017 [15] CHAIKOVSKA I, CASSOU K, CHICHE R, et al. Scientific Reports, 2016, 6(1): 36569. doi: 10.1038/srep36569 [16] MILBURN R H. Phys Rev Lett, 1963, 10: 75. doi: 10.1103/PhysRevLett.10.75 [17] PANTELL R, SONCINI G, PUTHOFF H. IEEE Journal of Quantum Electronics, 1968, 4(11): 905. doi: 10.1109/JQE.1968.1074989 [18] AVETISSIAN H K, MKRTCHIAN G F. Phys Rev E, 2002, 65: 046505. doi: 10.1103/PhysRevE.65.046505 [19] AVETISSIAN H K, MKRTCHIAN G F. Phys Rev ST Accel Beams, 2007, 10: 030703. doi: 10.1103/PhysRevSTAB.10.030703 [20] BONIFACIO R. Nucl Instr and Meth A, 2005, 546(3): 634. [21] BONIFACIO R, PIOVELLA N, COLA M, et al. Nucl Instr and Meth A, 2007, 577(3): 745. [22] ELIASSON B, SHUKLA P K. Phys Rev E, 2011, 83: 046407. doi: 10.1103/PhysRevE.83.046407 [23] ELIASSON B, LIU C S. Journal of Plasma Physics, 2013, 79(06): 995. doi: 10.1017/S0022377813000779 [24] SON S, JOON MOON S. Physics of Plasmas, 2012, 19(6): 063102. doi: 10.1063/1.4714798 [25] ZHANG Q R. Chinese Physics B, 2014, 23(1): 010306. doi: 10.1088/1674-1056/23/1/010306 [26] MADEY J M J. Journal of Applied Physics, 1971, 42(5): 1906. doi: 10.1063/1.1660466 [27] DI PIAZZA A, MÜLLER C, HATSAGORTSYAN K Z, et al. Rev Mod Phys, 2012, 84: 1177. doi: 10.1103/RevModPhys.84.1177 [28] FEYNMAN R P. Quantum Electrodynamics: Volume 3[M]. Boulder: Westview Press, 1998: 94. [29] ZHANG Q R. Chinese Physics B, 2015, 24(5): 054208. doi: 10.1088/1674-1056/24/5/054208 [30] TOUFEXIS F, TANTAWI S G. Phys Rev Accel Beams, 2019, 22: 120701. doi: 10.1103/PhysRevAccelBeams.22.120701