摘要:达到高能量密度物理状态后, 光束在介质中的传播行为与经典光学研究范畴相比, 会出现一些新现象. 比如在各向同性介质内可出现光束传播方向改变的现象. 另一方面, 高能量密度物理实验中由高功率激光器产生的束匀滑光束较为常见. 本文分析了空间和时间束匀滑光束在各向同性等离子体传播中出现束偏折现象的机制和条件, 并利用三维激光等离子体相互作用程序LAP3D进行了验证. 模拟表明只有当同时发生成丝不稳定性和存在离子声速量级的横向流时束匀滑激光才会产生显著的束偏折现象.
关键词: 高能量密度物理/
光束传播/
束匀滑光束/
束偏折

English Abstract

全文HTML

--> --> -->

在经典光学范围内, 光束在各向同性介质内的传播方向不会发生改变, 传播方向改变的现象通常发生在不同介质的交界面处. 但当能量驱动源如激光, 使介质能量增加并达到高能量密度状态后, 激光强度足以改变介质折射率, 在这种介质内也会发生光束传播方向改变的现象(束偏折)^[1]. 另一方面在高能量密度物理研究中, 随着近年来高功率激光器技术的发展, 采用激光束辐照产生高能量密度状态成为一种常用的方式. 许多采用大型高功率激光器的实验, 通常输出束匀滑光束^[2-8]. 针对高能量密度状态下, 束匀滑光束发生束偏折的机制和条件, 本文给出了相关研究结果.
束匀滑光束按光束匀滑方式可分为空间束匀滑和时间束匀滑两种. 它们的共同特点是光束包络较为平滑, 但内部结构较为复杂, 内部包含数千个小光斑, 纵向长度为$8{F^2}{\lambda _0}$^[1,2], 横向长度为$2 F{\lambda _0}$^[1,2], 这里F为透镜光学F数, ${\lambda _0}$为激光波长. 分析束匀滑光束发生束偏折的条件, 要从光束内部多个小光束相互作用角度考虑, 涉及到多体相互作用, 仅应用理论分析较为困难. 如果不考虑光束内部结构, 仅从整个光束角度考虑, 其发生束偏折的条件在一定程度上可类比单个光束, 如高斯型光束发生束偏折的机制和条件, 本文从这个角度展开了对束匀滑光束发生偏折现象所需条件的分析和研究.
实验研究^[9,10]和近期利用LAP3D程序对高斯型光束的模拟研究^[11]都表明, 高斯型激光束发生显著偏折需满足两个条件: 发生成丝不稳定性和存在声速量级的等离子体横向流^[11]. 前者在聚焦光斑处形成密度凹坑并改变此区域折射率, 后者使这种密度空间结构能在横向整体移动. 据此认为研究束匀滑光束发生束偏折的条件可进一步细化为研究以下两点: 束匀滑光束发生成丝不稳定性的条件以及声速量级的等离子体横向流对成丝不稳定性的影响.
本文根据以上两个条件对束匀滑光束偏折现象的产生机制和条件进行了研究, 给出了用于模拟三维束匀滑光束偏折现象的数值模拟程序LAP3D的物理建模, 分析了横向流大小对不稳定性发展的影响, 给出了研究空间束匀滑光束和时间束匀滑光束发生束偏折条件的数值模拟结果和分析.

2

-->

2.1.光束传播模型

为描述光束在等离子体中的传输和等离子体响应, 三维激光等离子体相互作用程序LAP3D分别采用了激光包络方程和双温流体方程组^[12,13]的物理模型. 由于束偏折现象发生与否跟光束成丝现象紧密相关, 因此在激光传播模型中只考虑衍射项和折射项, 忽略激光受激散射项, 对应的包络化激光传播方程为

$\left( {\frac{\partial }{{\partial t}} + {v_{\rm{g}}}\frac{\partial }{{\partial z}} - \frac{{{\rm{i}}{c^2}}}{{2{\omega _0}}}D_ \bot ^2 + \nu } \right)E = \frac{{ - {\rm{i}}2{\text{π}}{{\rm{e}}^2}}}{{{m_{\rm e}}{\omega _0}}}\delta {n_{\rm e{\rm{f}}}}E,$

其中, E为激光电场强度; ${v_{\rm{g}}}$ 为激光波包传播群速度; c为真空中激光光速; ${\omega _0}$为激光频率; $\nu $为逆轫致吸收率; ${m_{\rm e}}$为电子质量; $\delta {n_{\rm e{\rm{f}}}}$为电子扰动数密度; $D_ \bot ^2$为广义衍射算子^[12], 在下文推导中为简化用$\nabla _ \bot ^2$来代替.
2 -->

2.2.流速对束偏折(成丝)发展影响的分析

从前期对单个高斯光束偏折的模拟研究^[11]发现, 产生束偏折的空间位置与离子声速大小有以下关系: 当横向流为离子声速时, 发生束偏折时的光束传播距离最短; 当横向流速逐渐小于离子声速时, 发生束偏折的位置也逐渐“靠后”, 即光束传播距离逐渐增加. 当横向流速远大于离子声速后, 束偏折现象消失. 通过下面分析可知, 这实际上是横向流速大小对成丝不稳定性发展影响的一种体现.
设沿激光传播方向为z方向, 存在横向(沿x方向)流速为${u_ \bot }$的等离子体流场, 考虑流体方程组中的质量连续方程:

${\partial \rho }/\partial t + \nabla \cdot ({ u}) = 0, $

和动量方程:

$\dfrac{\partial { S}}{\partial t} + \nabla \cdot ({ {S u}}) + \nabla Q = - \nabla (p + p_{\rm e}) - \rho \nabla \psi , $

其中, ${{S}}$为电子动量, ${{u}}$为流速, ${{Q}}$为人为粘性应力张量, $\psi $为有质动力势, p和${p_e}$为离子压力和电子压力.
对上面两个方程线性化, 设$E = {E_0} + \delta E$, $n = {n_0} + \delta {{{n}}_{\rm ef}}$, 则质量连续方程和动量方程经化简后, 可得$\delta {n_{e{\rm{f}}}}$满足的密度涨落方程为

$\begin{split}&\left( {\frac{{{\partial ^2}}}{{\partial {t^2}}} + u_ \bot ^2\frac{\partial }{{\partial {x^2}}} - c_{\rm{s}}^2{\nabla ^2}} \right)\delta {n_{\rm e{\rm{f}}}}\\ =& \;\frac{{Z{n_{\rm e}}{e^2}}}{{mM{\omega _0}}}{\nabla ^2}\left( {{E_0} \cdot \delta E} \right),\end{split}$

其中, ${c_{\rm{s}}}$为离子声速, M为离子质量, 对(1)式和(4)式做傅里叶分析可得色散关系为

$\begin{split} &{\omega ^2} - k_ \bot ^2\left( {c_{\rm{s}}^2 - u_ \bot ^2} \right) = \frac{{k_ \bot ^2v_{{\rm{os}}}^2\omega _{{\rm{pi}}}^2}}{4}\\\; & \times\left[ {\frac{1}{{D\left( {\omega \!-\! {\omega _0},k \!-\! {k_0}} \right)}} \!+\! \frac{1}{{D\left( {\omega \!+\! {\omega _0},k \!+\! {k_0}} \right)}}} \right].\end{split}$

考虑$\omega = {\rm{i}}\gamma \ll {\omega _0}$并且${k} \cdot { k} _0 = 0$, 化简得到光束传播过程中存在横向小扰动及有横向流场发生成丝不稳定性对应的空间增长率:

$\frac{{{\gamma ^2}}}{{k_ \bot ^2}} = {{\frac{1}{8}\frac{{{n_{\rm e}}}}{{{n_{{\rm{cr}}}}}}{{\left( {\frac{{{v_{{\rm{os}}}}}}{{{v_{\rm e}}}}} \right)}^2}k_ \bot ^2}\bigg/{\left( {1 - \frac{{u_ \bot ^2}}{{c_{\rm{s}}^2}}} \right)}} - \frac{{k_ \bot ^4}}{{4k_0^2}}, $

式中, $\gamma $为成丝增长率, ${k_ \bot }$为横向扰动波数, ${k_0}$为入射光束波数, ${v_{{\rm{os}}}}$和${v_{\rm e}}$分别为电子在光场中的抖动速度和热速度.
对于在各向同性等离子体中传播的高斯光束, 可认为在给定参数下(6)式中物理量${v_{\rm e}}$、${v_{{\rm{os}}}}$、${n_{\rm e}}$、${n_{{\rm{cr}}}}$、${k_0}$、${c_{\rm{s}}}$均为定值. 在光束的扰动波数${k_ \bot }$范围内, 以${k_{{\rm{opt}}}}$为例考察(6)式中横向流大小对不稳定性发展的影响. 这里${k_{{\rm{opt}}}}$为没有横向流时, 在此条件下使成丝不稳定性增长率取最大值时对应的扰动波数. 进一步分析(6)式可知, 当横向流速为离子声速时, 对应不稳定增长率为极大值; 当流速逐渐小于离子声速时, 不稳定性增长率逐渐减小, 发生成丝的空间位置也逐渐“靠后”; 当流速大于离子声速时, 成丝发展受抑制.
结合本节第一段所描述的模拟结果^[11], 上面分析表明, 对高斯型光束, 横向流对其束偏折的影响本质上是横向流对成丝不稳定性发展的影响. 考虑到束匀滑光束包络内部包含多个小焦斑, 其光强分布可近似等同于高斯光束, 因此可推断束匀滑光束发生显著的束偏折现象也需满足发生成丝不稳定性和存在声速量级横向流的条件. 本文下面两节将采用数值模拟研究的方法来进一步分析和判断.

激光器产生的入射光束通过位于聚焦透镜前的相位板后, 被分割成多个子光束, 在远场叠加形成空间束匀滑光束. 本文模拟中应用连续相位板(continuous phase plate, CPP), 它使子光束间相位差在0到${\text{π}}$之间连续变化, 对比随机相位板, 不仅在远场可产生更为平滑的包络分布, 也使焦斑中能量提高到入射光能量的95%^[14].
焦平面内经CPP后生成的空间束匀滑光束的电场幅值分布可表示为下式:

$\begin{split} E(x',y') =\; & \frac{1}{{{\lambda _0}f}}\int E(x,y)t(x,y)\\ & \times\exp \left[ { - {\rm{i}}\frac{{2{\text{π}}}}{{\lambda f}}\left( {xx' + yy'} \right)} \right] {\rm{d}}x{\rm{d}}y,\end{split}$

其中, f为透镜焦距, $E(x, y)$为入射光复振幅, $t\left( {x, y} \right)$为CPP的相位信息.
LAP3D程序采用笛卡尔坐标系, 模拟三维大尺度上发生的空间束匀滑光束偏折现象. 模拟中激光从(x, y)面入射, 沿z方向传播. (x, y, z)方向上对应的模拟尺度分别为$\left( {2048{\lambda _0}, 2048{\lambda _0}, 2048{\lambda _0}} \right)$, 单位为激光波长${\lambda _0}$, ${\lambda _0} = $ 0.351 μm. 模拟空间内等离子体均匀分布, 初始等离子体密度为$0.{\rm{1 }}{n_{\rm{c}}}$, 电子温度${T_e}$为$3\;{\rm{ keV}}$, 离子温度${T_{\rm{i}}}$为$1\;{\rm{ keV}}$, 离子电离度$Z = 1$.
模拟中采用光学F数为11, 长轴约为200 μm, 短轴约为100 μm的椭圆焦斑Φ200, 椭圆焦斑内的电场幅值分布如文献[15]中图1所示. 根据空间束匀滑光束入射面平均光强大小, 设计了5个模型. 首先模拟了在没有横向等离子体流的条件下, 空间束匀滑光束在均匀等离子体中的传播行为, 重点考察是否发生了光束成丝现象. 然后引入声速量级的横向流, 再次计算上述5个模型, 重点考察是否发生了光束偏折现象. 等离子体横向流由x方向左边界流入, 从x方向右边界流出, 垂直于激光传播方向. 这些模拟结果如表1所示.
图 1 不同激光强度下空间束匀滑光束对应的束偏折模拟结果　(a) Φ200模型2加横向流; (b) Φ200模型5加横向流. 图中横纵坐标对应模拟空间坐标z和x, 其量纲为激光波长. 横向流速等于离子声速
Figure1. Beam deflection simulation results at different incident intensity: (a) Transverse flow and average intensity lower than filamentation threshold; (b) transverse flow and average intensity higher than filamentation threshold. x and yaxes of two figures corresponding to xand z axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length. The transverse flow speed equals ion sound speed.

模型	光斑平均强度/${\rm{W} } \cdot {\rm{c} }{ {\rm{m} }^{ {\rm{ - 2} } } }$	成丝现象	束偏折现象
1	4.30 × 1013	否	否
2	3.86 × 1014	否	否
3	1.07 × 1015	是	是
4	2.11 × 1015	是	是
5	3.49 × 1015	是	是

表1空间束匀滑光束Φ200成丝和束偏折现象模拟结果
Table1.Simulation results for filamentaion and beam deflection in the case of CPP smoothed beam Φ200.

对于空间束匀滑激光在均匀等离子体传播过程中能否发生成丝现象, 可根据下面公式^[15]来判断,

$ {I_{{\rm{av}}}} \leqslant {\left( {0.8{\text{π}}\frac{{{n_{\rm e}}}}{{{n_{\rm{c}}}}}\frac{{{F^2}\lambda _0^2}}{{{T_{\rm e}}({\rm{keV}})}}} \right)^{ - 1}}, $

其中, ${I_{{\rm{av}}}}$, ${\lambda _0}$, ${n_{\rm e}}$, ${n_{\rm{c}}}$, ${T_{\rm e}}$和F分别为光斑平均光强、激光波长、电子密度、等离子体临界密度、电子温度和光学F数. 当空间束匀滑入射光束平均光强大于${I_{{\rm{av}}}}$时, 空间束匀滑光束在传播过程中会发生成丝现象; 反之, 小于此值时, 空间束匀滑光束不成丝. 根据上文给出的模拟参数, 可知表1中模型3, 4, 5可发生成丝不稳定性情况, LAP3D模拟结果也证实了这一点.
对于束偏折现象的模拟表明, 如果空间束匀滑光束发生了成丝现象, 则引入声速量级的横向流后, 光束会发生显著的束偏折现象; 反之, 如果光束没有发生成丝现象, 则不会出现束偏折现象. 如图1所示, 模型2中没有发生激光成丝现象, 加横向流时对应的模拟结果中光束传播方向没有变化; 而模型5中发生了激光成丝现象, 加横向流后对应模拟结果显示光束传播方向随传播距离增加明显偏向等离子体流方向, 即发生了束偏折现象, 这些模拟结果与物理预期相同.
此外从空间束匀滑光束入射面和出射面内光斑电场幅值空间分布情况也可判断光束是否发生了明显的束偏折现象. 如图2所示, 对比Φ200模型5加横向流时模拟结果给出的电场幅值在入射面和出射面内的分布, 可以发现在出射面内光斑整体偏向等离子体流的方向, 表明发生了束偏折现象.
图 2 对比空间束匀滑光束发生束偏折时Φ200模型5加横向流的入射面和出射面内光斑电场幅值分布　(a)入射面; (b)出射面. 图中横纵坐标对应模拟空间坐标x和y, 其量纲为激光波长. 横向流速等于离子声速
Figure2. Comparison of spatial distribution of laser electric field between laser entrance and exit planes as beam deflection presents: (a) Laser entrance plane; (b) laser exit plane. x and y axes of two figures corresponding to x and y axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length. The transverse flow speed equals ion sound speed.

Φ200对应的模拟结果表明, 影响空间束匀滑光束明显发生束偏折的物理条件与高斯光束时的主要物理因素相一致, 即空间束匀滑光束满足成丝条件和存在流速在离子声速量级的等离子体横向流.

采用时间束匀滑的典型应用为谱色散匀滑化技术, 它利用光栅、三倍频晶体, 展宽入射光束频率, 然后使光束通过相位板和聚焦透镜, 形成远场的谱色散匀滑(smoothing by spectral dispersion, SSD)光束, 其传播主要特征表现为光束在横向上的“扫描”行为^[16]. SSD光束在焦平面上的光强分布为

$\begin{split}& E(x,y)= \\ \; &{E_0}\exp \left( {{\rm{i}}3\omega t} \right)\sum\limits_n {{J_n}\left( {3\delta } \right)} \frac{{\sin \left( {n\gamma + q} \right)}}{{n\gamma + q}}\frac{{\sin p}}{p}\sum\limits_{KL}\\ &\times{\exp \left[ {{\rm{i}}\left( {n{\Omega _m}t - 2n\gamma L - 2Lq - 2Kp + {\phi _{KL}}} \right)} \right]} ,\end{split} $

式中, $3\delta $为调制振幅; ${\varOmega _m}$为调制频率; $(p, q) = $$(x, y){{k\varDelta }/({2 f}})$, $x, y$为焦平面上的空间坐标, k是三倍频后的波数, $\varDelta $为相邻相位板面元之间的距离; ${\phi _{KL}}$为$\left( {K, L} \right)$通过相位板面元造成的相移; $\gamma = {{\alpha \varDelta }/2}$, 是由不同频率以不同角度入射造成的位移. 修改(9)式可得LAP3D程序中采用的入射面上光强分布的数值计算模型为^[16]

$E(x,y) = \sum\limits_n {{\rm{ETEMP}}\left( {x,y,n} \right)} {J_n}\left( {3\delta } \right)\exp \left( {{\rm{i}}n{\varOmega _m}t} \right), $

其中

$\begin{split}&{\rm{ETEMP}}(x,y,n) \\= \; &{E_0}{\rm{sinc}}[0.5{k_{0y}}(n \cdot {N_{\rm{c}}}F\lambda + y)] \cdot {\rm{sinc}}[0.5{k_{0x}}x]\\ & \times \sum\limits_{KL}\exp[ - {\rm{i}}L{k_{0y}}(n \cdot {N_{\rm{c}}}F\lambda + y) \\ &- {\rm{i}}K{k_{0x}}x + {\rm{i}}{\varphi _{KL}}] ,\end{split}$

这里, ${k_{0 y}} \equiv \dfrac{{2{\text{π}}}}{{F{N_y}\lambda }}$, ${k_{0 x}} \equiv \dfrac{{2{\text{π}}}}{{F{N_x}\lambda }}$; ${N_x}$, ${N_y}$分别为沿x, y方向上相位板的面元个数. 相较于(9)式, LAP3D程序计算每一时间步采用(11)式, 可大幅节约SSD光束传播过程的计算时间并提高程序并行计算能力^[16].
仿照上节研究空间束匀滑光束出现束偏折现象条件的思路, 可认为对于SSD匀滑光束, 发生束偏折的前提是发生成丝不稳定性. 因此这里首先考察SSD光束发生成丝的条件. 美国****Kruer等^[17-21]研究分析了有一定带宽的激光抑制不稳定性发展的条件, 指出当不稳定性增长率为${\gamma _0}$的波引入带宽$\Delta \omega $时, 不稳定性引发的振幅增长可被削弱为没有带宽时的${{{\gamma _0}}/{\Delta \omega }}$倍, 即

$\big\langle {\tilde f} \big\rangle = \beta \exp \left( {\frac{{\gamma _0^2t}}{{\Delta \omega }}} \right),$

这里, $\tilde f$对应不稳定性增长的幅值, $\beta $由初始条件确定. 此式表明当激光带宽远大于不稳定性增长率${\gamma _0}$时, 可有效抑制不稳定性增长. 由此可认为当SSD光束带宽大于成丝不稳定性增长率时, 成丝不稳定性的发展会受到抑制.
为考察SSD带宽对其在等离子体传播过程中发生束偏折的影响, 用LAP3D模拟了三种不同调制频率的SSD光束在均匀等离子体中的传播, 这里取调制频率分别为${10^{ - 3}}{\omega _0}$, ${10^{ - 4}}{\omega _0}$和${10^{ - 5}}{\omega _0}$. 模拟中取入射SSD光束光斑平均光强为$2.11 \times $$ {10^{15}}\;{{\rm{W}}/{{\rm{c}}{{\rm{m}}^2}}}$, 其余模拟参数与上节中算例相同.
根据模拟参数, 可分析哪种调制频率会出现束偏折现象. 成丝不稳定性最大增长率${\gamma _{\max }}$为

${\gamma _{{\rm{max}}}} = \frac{1}{8}{\left( {\frac{{{v_{{\rm{os}}}}}}{{{v_{\rm e}}}}} \right)^2}\frac{{{\omega _{{\rm{p}}\rm e}}}}{{{k_0}{c^2}}},$

其中, ${v_{{\rm{os}}}}$, ${v_{\rm e}}$, ${\omega _{{\rm{p}}e}}$, ${k_0}$, c分别为电子在光场中的抖动速度、电子热速度、等离子体频率、光束波矢和真空光速. 根据上文模拟参数可知成丝不稳定性最大增长率${\gamma _{\max }} \approx 2.56 \times {10^{{\rm{ - }}4}}{\rm{ }}{\omega _0}$. SSD光束带宽$\Delta \omega $为

$\Delta \omega = 2 \cdot \delta \cdot {f_{{\rm{SSD}}}},$

其中, $\delta $为调制深度, 模拟中取$\delta = 12$, ${f_{{\rm{SSD}}}}$为调制频率. 由本文算例所取的SSD光束调制频率可知, 当SSD光束调制频率${f_{{\rm{SSD}}}}$为${10^{ - 3}}{\rm{ }}{\omega _0}$和${10^{ - 4}}{\rm{ }}{\omega _0}$时, SSD光束带宽均远大于算例对应的成丝最大增长率${\gamma _{\max }} \approx 2.56 \times {10^{ - 4}}{\rm{ }}{\omega _0}$. 但当SSD光束调制频率${f_{{\rm{SSD}}}}$为${10^{ - 5}}{\rm{ }}{\omega _0}$时, SSD光束带宽小于成丝最大增长率.
根据上面的理论分析, 可以预期在取${10^{ - 3}}{\rm{ }}{\omega _0}$和${10^{ - 4}}{\rm{ }}{\omega _0}$这两种调制频率时, 时间束匀滑光束可以抑制成丝不稳定性; 相应束偏折现象也不会出现. 当SSD光束调制频率为${10^{ - 5}}{\rm{ }}{\omega _0}$时, 可发生成丝不稳定性; 当同时存在声速量级横向流时会出现显著的束偏折现象. 下面给出的模拟结果也证实了这些判断.
图3给出了当SSD光束调制频率为${10^{ - 3}}{\rm{ }}{\omega _0}$且没有横向流时的传播过程的模拟结果. 模拟发现, 在给定的传播距离内(约2000个激光波长), 与空间束匀滑光束成直线传播相比, SSD光束在空间传播中表现出明显的“扫描”行为, 这和已有模拟结果^[16]和理论预期相符.
图 3 调制频率为${10^{ - 3}}{\omega _0}$的时间束匀滑光束传播行为　(a)对应11000激光周期; (b)对应13750激光周期. 图中横纵坐标对应模拟空间坐标z和y, 其量纲为激光波长
Figure3. Propagation of SSD beam at modulation frequency of 10^–3ω₀: (a) Corresponding simulation result at 11000 th laser periods; (b) corresponding simulation result at 13750 th laser periods. x and y axes of two figures corresponding to z and y axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length.

由上一节给出的模拟结果可知在给定的等离子体条件下, 本节算例中的SSD光束强度远大于成丝阈值, 但图3中显示SSD光束并没有出现束发散现象. 这表明当SSD光束取调制频率为${10^{ - 3}} {\omega _0}$时, 抑制了成丝不稳定性. 由此可以预期, 当加入等离子体横向流时, 也不会明显发生束偏折现象, 模拟结果证实了这一点(图4).
图 4 有横向离子声速量级等离子体流时调制频率为${10^{ - 3}}{\omega _0}$的时间束匀滑光束的传播行为　(a)对应11000激光周期; (b)对应13750激光周期. 图中横纵坐标对应模拟空间坐标z和y, 其量纲为激光波长. 横向流速等于离子声速
Figure4. Propagation of SSD beam with transverse flow at modulation frequency of 10^–3 ω₀: (a) Corresponding simulation result at 11000 th laser periods; (b) corresponding simulation result at 13750th laser periods. x and y axes of two figures corresponding to z and y axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length. The transverse flow speed equals ion sound speed.

除调制频率为${10^{ - 3}} {\omega _0}$的情况, 还模拟了SSD光束调制频率为${10^{ - 4}}{\omega _0}$时, 无横向流和有横向流时的激光传播情况. 模拟结果表明, 无横向流时算例没有出现激光成丝现象, 相应模拟流速在离子声速量级的横向流时的算例如图5所示, 也没有出现束偏折现象. 当调制频率为${10^{ - 5}} {\omega _0}$时, 模拟结果表明: 无横向流时, SSD光束会出现成丝现象; 有横向流时, 光束会出现束偏折现象. 由于其调制频率较低, 所以在本文所选的模拟空间尺度内, SSD光束传播行为类似于前一节的空间束匀滑激光, 如束发散时光场空间分布的模拟结果与图1(b)类似, 故这里略去.
图 5 对比时间束匀滑光束在调制频率为${10^{ - 4}}{\rm{ }}{\omega _0}$时的光束传播行为　(a)等离子体横向流速为零; (b)等离子体横向流速等于离子声速. 图中横纵坐标对应模拟空间坐标y和z, 其量纲为激光波长
Figure5. Propagation of SSD beam at modulation frequency of 10^–4ω₀: (a) No transverse flow; (b) the transverse flow speed equals ion sound speed. x and y axes of two figures corresponding to y and z axes of simulation coordinates, respectively. The spatial scale is in unit of laser wave length.

综上理论分析和模拟结果都表明, 当SSD光束带宽远大于成丝稳定性最大增长率时, 通过抑制成丝不稳定性发展, 可抑制束偏折现象出现.

通过对高斯型光束发生束偏折所需条件进行引申, 提出了束匀滑光束发生束偏折的条件, 并应用数值模拟程序进行了验证和研究. 模拟证实, 对于空间束匀滑激光和时间束匀滑激光来说, 只有同时满足发生成丝不稳定性和存在离子声速量级的横向流时, 才会产生显著的束偏折现象. 这表明束匀滑光束发生束偏折的物理成因与高斯型光束类似, 即成丝不稳定性产生的密度空间结构变化——“凹坑”, 以及在等离子体横向流作用下产生横向移动. 结合初步理论分析表明, 横向流速为离子声速时束偏折现象最显著.
研究表明, 在均匀等离子体及存在离子声速量级横向流时, 如果空间束匀滑光束入射光斑平均光强大于成丝阈值, 可引发成丝不稳定性从而发生束偏折现象; 对于时间束匀滑光束, 当其带宽远大于成丝不稳定性最大增长率时, 可以抑制成丝不稳定性从而避免束偏折现象产生, 反之则会发生束偏折现象.