删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

惯性约束聚变装置中靶面光场特性的统计表征方法

本站小编 Free考研考试/2021-12-29

摘要:在激光驱动的惯性约束聚变装置中, 常采用多种束匀滑手段对焦斑的时空特性进行调控. 光传输链路中涉及的光学元件众多、传输变换复杂, 往往导致光传输模型复杂, 且在运用衍射光学方法分析焦斑形态和特征时面临大量的数据处理和计算, 致使出现计算量大、计算效率低等问题, 亟需寻求快速而简便的新方法来描述焦斑的统计特征. 本文利用光场特性的统计表征方法对靶面光场进行表征, 采用圆型复数高斯随机变量直接描述靶面光场的统计特征, 并基于典型焦斑评价参数对衍射光学方法和统计表征方法得到的远场焦斑进行了对比和分析. 结果表明, 采用衍射光学方法和统计表征方法获得的焦斑的瞬时特征基本一致, 其时间积分的远场焦斑有所不同, 但仍可进一步采用相关系数来表征其远场焦斑的时间变化特征.
关键词: 统计光学/
惯性约束聚变/
束匀滑/
焦斑

English Abstract


--> --> -->
激光间接驱动的惯性约束聚变(inertial confinement fusion, ICF)装置中, 靶面辐照均匀性是影响靶丸有效压缩的关键因素[1]. 目前已发展多种束匀滑技术, 如一维光谱角色散匀滑技术[2] (1D smoothing by spectral dispersion, 1D-SSD)、位相板[3] (phase plate, PP)和偏振匀滑[4, 5] (pol- arization smoothing, PS)等, 对靶面光斑的时域、空域和偏振特性进行控制. 然而, 由于在对靶面光斑特征进行计算分析时通常需要对大量的光场分布数据进行计算处理, 进而采用各种评价指标对其进行表征, 因而大规模的并行计算必不可少. 例如, 对于不含时间变量, 空间采样数N = 1024 的单束光双精度处理, 需要的储存量约为16 Gbits; 对于含时间变量, 时间采样数Nt = 64的10束光双精度处理, 需要的储存量为10 T. 然而幸运的是, 目前已发展了多种方法对靶面光强的统计特性进行表征, 其中最具代表性的是在兆焦耳激光(laser mega Joule, LMJ)装置的研究中报道了多光束下焦斑的时空统计特征, 并对不同束匀滑方式下焦斑的尺度、脉宽和速度的统计分布进行了分析[6-9]. 此外, 在束匀滑技术中, 由于连续位相板(continuous phase plate, CPP)的位相分布由随机数多次迭代获得[10], 尽管不同CPP得到的焦斑细节不尽相同, 但其均具有类似的统计特征. 因此, 激光束经过CPP的传输可视为激光束经过随机表面的传输过程. 尽管靶面散斑在不同位置的光强变化具有很强的随机性, 而且不同束匀滑方式得到的靶面散斑随机分布不相同, 但其整体分布却满足一定的统计规律. 因此, 焦斑不能用常规的函数进行描述与表征, 而记录每个点的光强十分繁琐且占用空间大, 需要寻求有效的表征方法. 光场特性的统计表征方法可以描述这种随机过程的统计特征, 通过合理的假设和公式推导, 根据近场光场的统计特性, 避开从近场到远场的数值计算, 进而采用满足一定统计规律的解析表达式直接对靶面光场进行描述与表征, 且无需对靶面光场的每个点进行数据处理.
在以美国国家点火装置为代表的ICF装置中, 往往采用光谱角色散匀滑技术(smoothing by spectral dispersion, SSD), CPP和偏振控制联用的束匀滑方案[11, 12]. SSD的基本原理[13, 14]是利用光栅对时间相位调制后的光束进行色散, 使激光束在远场的散斑产生扫动, 从而在一定积分时间内抹平焦斑的强度调制. 为了实现对焦斑的超快速匀滑, 还提出了基于光克尔效应的径向匀滑(radial smoothing, RS)方案, 即利用光克尔效应实现焦斑尺寸在ps量级的超快变化, 从而抹平焦斑强度调制[15, 16]. 然而, 其焦斑的统计特性尚不清楚. 本文采用光场特性的统计表征方法对靶面光场进行表征, 利用圆型复数高斯随机变量对靶面光强的统计特征进行描述, 并对上述两种束匀滑方案下的瞬时和时间平均下焦斑的统计特征进行了分析.
2
2.1.衍射积分模型
-->以NIF中单束激光为例, 激光束经过光谱角色散、位相板, 最后传输到焦平面[17], 如图1(a)所示. 图1(b)则给出了激光束在径向匀滑方案中的传输示意图.
图 1 激光束传输和聚焦光路 (a) SSD+CPP; (b) RS+CPP
Figure1. Transmission and focusing light path of laser beam: (a) SSD+CPP; (b) RS+CPP.

假设输入光场为具有振幅调制和位相畸变的超高斯光束, 其光场的时空分布表达式为
$\begin{split}E\left( {x,y,t} \right) = &\; {E_0}\left( t \right)\left( {{\rm{1 + }}{\sigma _{{\rm{random}}}}} \right)\\ & \times{\rm{exp}}\left[ { - \left( {\frac{{{x^{2N}}}}{{{w^{2N}}}} + \frac{{{y^{2N}}}}{{{w^{2N}}}}} \right)} \right]{\rm{exp}}\left( {{\rm{i}}{\varphi _{{\rm{initial}}}}} \right),\end{split}$
式中E0(t)为电场强度, ${\sigma _{{\rm{random}}}}$为随机振幅调制, w为激光束束腰, N为空间超高斯阶数, ${\varphi _{{\rm{initial}}}}$为初始位相畸变, xy为近场坐标.
在常规SSD匀滑方案中, 临界色散对应的色循环数为1, 因此取色循环数Nc = 1[18], 在经光谱角色散匀滑、三倍频和连续相位板后, 近场光场E1near可表示为
$\begin{split}{E_{{\rm{1near}}}}\left( {x,y,t} \right) =\;&E\left( {x,y,t} \right){\rm{exp}}\left[ {{\rm{i}}3{\omega _0}t} \right.\\& \left. { +\; {\rm{i}}3\delta {\rm{sin}}\left( {2{\text{π}}{v_{\rm{m}}}t + \alpha X} \right)} \right]{\rm{exp}}\left( {{\rm{i}}{\varphi _{{\rm{CPP}}}}} \right),\end{split}$
式中E为超高斯光场分布, ${\omega _0}$为入射光中心角频率, δ为调制深度, vm为调制频率, α为色散系数, ${\varphi _{{\rm{CPP}}}}$为CPP对光场的附加相位.
在RS方案中, 经过径向匀滑装置、三倍频和连续相位板后, 激光束的近场光场可表示为
${E_{{\rm{2near}}}}\left( {x,y,t} \right) = E\left( {x,y,t} \right){\rm{exp}}\left[ {{\rm{i}}\left( {3{\varphi _{{\rm{RS}}}} + {\varphi _{{\rm{CPP}}}}} \right)} \right],$
式中${\varphi _{{\rm{RS}}}}$为光克尔介质引入的球面位相调制. 可利用含时变光场柯林斯公式[19]计算激光束在靶面的瞬时光场分布, 进而计算积分时间内的靶面光强分布. 于是, 远场瞬时光场分布为
$\begin{split} {E_{\rm{f}}}\left( {{x_{\rm{f}}},{y_{\rm{f}}},t} \right) =\; & \frac{{\exp \left( {{\rm{i}}{k_{3\omega }}L} \right)}}{{{\rm i}{\lambda _{3\omega }}B}} \int\limits_{ - \infty }^\infty \int\limits_{ - \infty }^\infty {E_{j{\rm{near}}}}\left( {x,y,t} \right)\\ &\times\exp \left[ { - \frac{{{\rm{i}}{k_{3\omega }}}}{B}\left( {x{x_{\rm{f}}} + y{y_{\rm{f}}}} \right)} \right] {\rm{d}}x{\rm{d}}y, \\ &\left( {j = 1,2} \right),\end{split}$
式中${k_{3\omega }} = 2{\text{π}}/{\lambda _{3\omega }}$为三倍频后的激光束的波数; L = 2f, f为聚焦透镜的焦距; B = f/β0, β0为扩束系统的扩束比; ${x_{\rm{f}}}$, ${y_{\rm{f}}}$为远场坐标.
在积分时间Δt内的靶面平均光强可表示为
${I_{\rm{f}}}\left( {{x_{\rm{f}}},{y_{\rm{f}}}} \right) = \frac{1}{{\Delta t}}\int_0^{\Delta t} {{{\left| {{E_{\rm{f}}}\left( {{x_{\rm{f}}},{y_{\rm{f}}},t} \right)} \right|}^2}} {\rm{d}}t,$
式中Δt为积分时间.
2
2.2.光场特性的统计表征方法
-->在ICF装置中, 常采用CPP对激光束的焦斑进行空间整形[20]. CPP的随机性主要体现在随机种子数上, 而确定性主要体现在相位滤波函数上. 相位滤波过程中通过改变滤波截止频率可以获得不同最小空间周期的CPP[21]. CPP的位相分布由随机数多次迭代获得, 保留了一定的随机特征, 同时其位相梯度也具有一定的确定性. 采用最小空间周期一定而随机数种子不同的CPP, 用相同的统计方法提取不同CPP的位相统计特征, 其位相的统计分布如图2所示.
图 2 不同随机数种子得到的CPP的位相统计分布
Figure2. Statistical distribution of the phase of CPP obtained from different random number seeds.

图2可知, 采用不同随机数种子设计得到CPP的统计位相分布大致吻合, 满足正态分布的统计特性. 位相梯度是影响焦斑分布的关键物理量, 能够较好地反映光学元件引入畸变波前的低频特性. 为了进一步分析不同CPP位相分布之间的误差, 采用均方根梯度(gradient root-mean square, GRMS)表征CPP面形对远场光场的影响[22]. 根据梯度的概念, CPP沿x, y方向上的位相梯度函数gx(x, y), gy(x, y)分别是CPP沿x, y方向的一阶偏微分:
${g_x}\left( {x,y} \right) = \frac{{\partial {\varphi _{{\rm{CPP}}}}\left( {x,y} \right)}}{{\partial x}},$
${g_y}\left( {x,y} \right) = \frac{{\partial {\varphi _{{\rm{CPP}}}}\left( {x,y} \right)}}{{\partial y}}.$
则CPP的总面形梯度为
$g\left( {x,y} \right) = \sqrt {{g_x}{{\left( {x,y} \right)}^2} + {g_y}{{\left( {x,y} \right)}^2}} ,$
g(x, y)的均方根值为CPP的均方根梯度GRMSCPP[22]:
${\rm{GRM}}{{\rm{S}}_{{\rm{CPP}}}}\!=\!\sqrt {\frac{{\left\{ {\displaystyle\sum\limits_{i = 1}^n {\displaystyle\sum\limits_{j = 1}^m {{{\left[ {{g_{i,j}}\left( {x,y} \right) \!- \!\bar g\left( {x,y} \right)} \right]}^2}} } } \right\}}}{{\left( {n - 1} \right)\left( {m - 1} \right)}}},$
式中$\bar g(x,y)$为位相梯度的平均值; ${g_{i,j}}\left( {x,y} \right)$为离散点对应的位相梯度值; n, m分别为离散点对应的行数和列数.
计算最小空间周期相同而随机数种子不同的CPP的GRMS值, 得到3个CPP的GRMS值分别为0.4540, 0.4480和0.4479, 即采用最小空间周期相同而随机数种子不同的CPP得到的GRMS值基本一致. 由此可见, 不同随机数种子设计得到不同CPP的位相分布在一定误差范围内满足相同的统计规律. 这说明了CPP位相分布的随机性与确定性并存. 对于这种不完全随机的统计分布特征, 更有利于我们从其中提取出规律性, 可进一步根据CPP面形的统计特征推导出靶面光场的统计特征.
CPP的面形分布是连续且随机的, 因而可以将其看作一个表面高度为随机函数的衍射光学元件[23]. 因此, 激光束经过CPP汇聚至靶面可视为光源照明粗糙表面产生散射光的过程, 可以采用统计光学的理论模型对焦斑特性进行统计分析.
图3所示, 将CPP看成许多小单元构成的位相元件, 且对激光束的附加位相满足某一类统计分布. 当对位相板划分的单元数足够多时, 可用圆型复数高斯随机变量对瞬时焦平面进行描述:
图 3 经CPP调制后激光束汇聚至靶面的过程
Figure3. The process of the laser beam converged to the target plane after the modulation of CPP.

$\begin{split}A\left( {x,y} \right) = & \sum\limits_{j = 1}^N {\sum\limits_{l = 1}^N {{A_0}\left( {jh,lh} \right)} }\\ &\times\exp \!\left( {{\rm{i}}{\phi _{i,j}}} \right)\exp\!\left[\!{ - \frac{{{\rm{i}}2{\text{π}}}}{{\lambda F}}\left( {jhx + lhy} \right)} \right]\!,\end{split}$
式中将位相板划分成N × N个小单元, 不同小单元之间的振幅分布和位相分布不同; j, l分别是CPP横向和纵向划分的不同单元数, h是每个子单元的口径, F = f/h为系统的F数, f为透镜焦距.
在利用圆型复数高斯随机变量对瞬时焦斑进行描述时, (10)式满足以下统计特征: 位相分布?满足与连续位相板位相分布相似的统计分布规律, 而振幅概率密度函数则满足Rayleigh分布, 即
${P_{{A_0}}}\left( {{A_0}} \right) = \frac{{{A_0}}}{{{\sigma ^2}}}\exp \left( {\frac{{{A_0}}}{{2{\sigma ^2}}}} \right).$
光强概率密度函数满足负指数分布:
${P_I}\left( I \right) = \frac{1}{{\overline I }}\exp \left( { - \frac{I}{{\overline I }}} \right).$
束匀滑方式、打靶构型、单束和集束方式均会影响散斑的统计分布, 因而在对瞬时焦平面光场进行描述时, (10)式中的A0?满足的统计分布会随着不同的束匀滑方式发生变化. 对于不同的束匀滑方式, 如采用KPP和RPP等空域束匀滑手段时, 需要根据相应的相位板面形的统计特性分析靶面光场的光强与位相统计特性. 对于实际打靶过程中由于光学器件性能等原因引起的与理论设计之间的偏差, 则需要针对具体的情况对统计光学表征模型加以完善和修正.
激光束经过位相板整形后再经过透镜聚焦, 不妨定义焦平面处电场的相关函数和相关系数[24]
${\varGamma _A}\left( {{x_1},{y_1};{x_2},{y_2}} \right) = \overline {A\left( {{x_1},{y_1}} \right){A^ * }\left( {{x_2},{y_2}} \right)} ,$
${\mu_A}\left( {\Delta x,\Delta y} \right) = {{{\varGamma _A}\left( {\Delta x,\Delta y} \right)}/{{\varGamma _A}\left( {0,0} \right)}},$
式中(x1, y1)和(x2, y2)是焦平面上的任意两点; Δx = x1x2, Δy = y1y2, 焦平面场的相关函数等于透镜处光强的傅里叶变换, 即满足
$\begin{split}{\varGamma _A}\left( {\Delta x,\Delta y} \right) =\;& \frac{\kappa }{{{\lambda ^2}{f^2}}}\iint {I\left( {\alpha ,\beta } \right)}\\&\times\exp \left[ { - {\rm{j}}\frac{{2{\text{π}}}}{{\lambda f}}\left( {\alpha \Delta x + \beta \Delta y} \right)} \right]{\rm{d}}\alpha {\rm{d}}\beta ,\end{split}$
式中I (α, β)为透镜前的近场光强分布, κ是量纲为长度平方的常量. 焦平面的功率谱密度函数PSDI (νx, νy)代表强度涨落功率在二维频率平面上的分布, 它等于近场光强自相关函数的傅里叶变换, 即靶面光场的PSD可由近场光强的自相关求得. 由自相关定理, PSDI (νx, νy)可表示为

${\rm{PS}}{{\rm{D}}_I}\left( {{\nu _x},{\nu _y}} \right) = {\overline I ^2}\left[ {{\text{δ}} \left( {{\nu _x},{\nu _y}} \right) + {{\left( {\lambda f} \right)}^2}\frac{{\displaystyle\iint {I\left( {\alpha ,\beta } \right)I\left( {\alpha + \lambda f{\nu _x},\beta + \lambda f{\nu _y}} \right){\rm{d}}\alpha {\rm{d}}\beta }}}{{{{\left[ {\displaystyle\iint {I\left( {\alpha ,\beta } \right)}{\rm{d}}\alpha {\rm{d}}\beta } \right]}^2}}}} \right],$

式中δ函数对应于平均强度$\bar I$所贡献的零频率离散功率, 后一项则表示靶面光强的强度随空间频率的分布.
对于表征靶面光场平整度的评价函数δRMS, 其定义式为
${\delta _{{\rm{RMS}}}} \!=\! {\left\{ {\frac{{\displaystyle\iint_A {{{\left[ {{I_\nu}\left( {x,y} \right) - {I_{{\rm{obj}}}}\left( {x,y} \right)} \right]}^2}{\rm{d}}x{\rm{d}}y}}}{{\displaystyle\iint_A {{{\left[ {{I_{{\rm{obj}}}}\left( {x,y} \right)} \right]}^2}{\rm{d}}x{\rm{d}}y}}}} \right\}^{{1}/{2}}},$
式中A为焦斑的面积, Iν为焦平面的平均光强, Iobj为焦平面上不同点处的光强.
焦斑的光通量对比度的表达式为
$C{\rm{ = }}\frac{1}{{{I_\nu }}}\sqrt {{{\iint_{\left( {x,y} \right) \in A} {\left[ {{I_{{\rm{obj}}}}\left( {x,y} \right) - {I_\nu }} \right]^2}}}{\rm{d}}x{\rm{d}}y/A} .$
从统计光学的角度表述光通量对比度, 则有
$C = \frac{{{\sigma _I}}}{{\overline I }},$
式中${\sigma _I}$为靶面光强的标准差, $\bar I$是靶面光强的均值. 当目标函数为超高斯分布, 且只对中心平顶区域做统计时, 满足柯西施瓦茨不等式中等号成立的条件, 将(17)式中分子和分母同除以焦斑面积A, 可认为δRMS = C.
由Parseval定理, 可得到RMS与PSD的关系:
$\delta _{{\rm{RMS}}}^2 = \iint {{\rm{PS}}{{\rm{D}}_I}\left( {{\nu _x},{\nu _y}} \right)}{\rm{d}}{\nu _x}{\rm{d}}{\nu _y} = {C^2}, $
即光通量对比度的值的平方在一定条件下等于PSD的积分.
分析(20)式可知, PSD曲线的面积等于RMS的平方, 也在一定条件下等于光通量对比度的平方. 通过PSD的值与频率间隔相乘再求和可以得到RMS以及光通量对比度. 在建立了靶面光场的衍射积分模型和统计光学模型后, 我们将利用上述模型对激光束在近场和远场的瞬时、积分特性进行分析.
基于上述理论模型, 我们对靶面光场特性的统计表征方法的可行性进行分析, 并采用衍射积分模型中常规评价函数与统计表征方法中评价函数对不同束匀滑方式的焦斑均匀性进行分析和对比.
2
3.1.光场特性的统计表征方法表征瞬时焦斑的有效性验证
-->为了验证采用圆型复数高斯随机变量描述靶面光场的可行性, 先对靶面光场的统计特性进行分析. 输入光场的参数为: 光束束腰半径w = 186 mm; 超高斯阶数N = 6; 中心角频率ω0 = 1.79 × 1015 Hz; 光束波长λ0 = 1053 nm; CPP的PV值为7.3λ0. 由于CPP面形分布的随机性, 光束通过CPP后在靶面形成散斑. 对该散斑的光强和位相的统计特性进行分析, 典型结果如图4所示.
图 4 激光束经过CPP整形后靶面光强和位相统计特征 (a) CPP整形后的靶面光强分布; (b) CPP整形后的靶面振幅分布; (c) CPP位相与远场位相统计分布
Figure4. The statistical characteristics of the laser beam's intensity and phase of the target plane after CPP's shaping: (a) Intensity distribution of the target plane after CPP's reshaping; (b) amplitude distribution of the target plane after CPP's shaping; (c) statistical distribution of CPP's phase and far field phase.

图4(a)(b)可知, 激光束经过CPP整形之后, 在靶面形成散斑光强的统计分布近似为负指数分布, 振幅的统计分布近似为瑞利分布, 与圆型复数高斯随机变量满足相同的统计特征. 在图4(c)中, 拟合曲线服从正态分布, 可见CPP位相的统计分布大致为正态分布, 且CPP的位相分布的统计特性与靶面散斑的位相分布的统计特性大致一致. 因此, 圆型复数高斯随机变量中的位相分布?应该与CPP位相的统计分布一致, 即满足正态分布.
我们进一步利用衍射积分模型对激光束通过CPP后的靶面光强分布进行数值模拟, 并对由圆型复数高斯随机变量描述的靶面光强分布进行统计分析. 将连续位相板视为由512 × 512位相单元组成, 进而采用常规评价参数对两种方法得到的瞬时焦斑的光强和位相两个方面进行比较, 并运用FOPAI曲线比较两者靶面瞬时光强的不同峰值功率占总功率的份额, 采用位相统计分布规律来比较数值求解的远场位相分布与光场特性的统计表征方法得到的远场位相分布. 典型结果如图5所示.
图 5 数值求解与的瞬时远场光强特性比较 (a)瞬时焦斑光强FOPAI对比; (b)数值求解远场位相与解析求解远场位相统计特性
Figure5. Comparison of characteristics of instantaneous far-field intensity solved by numerical analysis and that Solved by analytical solution: (a) FOPAI's comparison instantaneous focal spot intensity; (b) statistical characteristics of numerical solution far-field phase and analytical solution far-field phase.

图5可以看出, 在满足统计光学假设的前提下, 数值求解与统计分析得到的瞬时远场光强分布的FOPAI曲线和位相分布的统计特征均能够较好地符合, 由此初步验证了采用圆型复数高斯随机变量对靶面光强分布进行描述和表征是可行的.
2
3.2.不同束匀滑方案下焦斑的统计特征
-->下面利用光场特性的统计表征方法, 对1D-SSD+CPP和RS+CPP两种不同束匀滑方案得到的焦斑进行统计特性分析. 1D-SSD+CPP方案计算时所采用的参数[13]为: 光束束腰半径w = 186 mm; 超高斯阶数N = 6; 中心角频率ω0 = 1.79 × 1015 Hz; 调制深度δ = 2.33; 调制频率ωm = 17 GHz; 光栅色散系数${\rm{d}}\theta /{\rm{d}}\lambda = 2156.8\;{\text{μ}}{\rm{rad}}/{\rm{nm}}$; 光束波长λ0 = 1053 nm; 三倍频后的波长λ = 351 nm; 激光调制带宽Δλ = 0.3 nm; SSD积分时间为20 ps; 透镜焦距f = 7.7 m. RS方案中计算所采用的参数[25]为: 抽运光的束腰宽度w = 148.8 mm; 抽运光的峰值强度Ip = 16 GW/cm2, 子脉冲之间的延迟时间Td = 10 ps; 子脉冲脉宽Tw = 4.5 ps; 光克尔介质选用硝基苯, 其光克尔系数n2 = 2 × 10–18 m2/W, 中心厚度d = 100 mm; 主激光束的计算参数与1D-SSD+CPP方案中的参数相同; 积分时间为20 ps (两个抽运光周期); 光通量对比度的积分区域为焦斑86.5%环围能量. 1D-SSD+CPP和RS+CPP两种束匀滑方案下的瞬时和积分时间后的焦斑PSD曲线如图6所示.
图 6 不同束匀滑方案下瞬时与积分焦斑的统计特性 (a) 1D-SSD+CPP瞬时、积分焦斑及其PSD; (b) RS+CPP瞬时、积分焦斑及其PSD
Figure6. Statistical characteristics of instantaneous and integral focal spots obtained by different beam smoothing schemes: (a) Instantaneous, integral focal spots and their PSD of 1D-SSD+CPP; (b) instantaneous, integral focal spots and their PSD of RS+CPP.

图6可知, 经过CPP整形后得到的两种瞬时焦斑统计特性十分相似但又不尽相同. 经过时间积分后, 1D-SSD+CPP 方案得到的焦斑内部出现了沿光栅色散方向的条纹状强度调制, 而RS+CPP方案因快速变焦使焦斑内部散斑在径向方向上产生扫动, 其焦斑在径向方向更为均匀, 因而时间积分后两种方案得到的焦斑在细节上不再具有相似性. 由此可见, 采用圆型复数高斯随机变量只能对瞬时的靶面光场进行描述. 此外, 由于CPP在设计过程中抑制了部分高频成分, 因而1D-SSD+CPP和RS+CPP两种不同束匀滑方案的瞬时与时间积分后的PSD曲线的中低频部分均能够基本重合, 但高频部分却存在差异. 根据(20)式, 分别对瞬时焦斑与积分时间后的焦斑PSD在频率域进行积分, 进而再对其进行开方. 对于1D-SSD+CPP方案, 瞬时焦斑范围内PSD积分开方值为1.079, 积分时间后焦斑范围内PSD积分开方值为1.067. 对于RS+CPP方案, 瞬时焦斑范围内PSD积分开方值为1.077, 积分时间后焦斑范围内PSD积分开方值为1.056. 由此可见, 上述两种方案的瞬时与积分时间后的焦斑PSD积分再开方值均近似相等.
下面进一步对1D-SSD+CPP和RS+CPP两种方案下的瞬时和积分焦斑统计特性进行了分析(表1).
StatisticsPSD integral square value
of instantaneous
Instantaneous luminous
flux contrast
PSD integral square value
of time integral
Integrated luminous
flux contrast
CPP1.0691.094
Statistical optical0.9790.987
SSD+CPP1.0791.0931.0670.514
RS+CPP1.0771.0901.0560.478


表1瞬时与积分焦斑的PSD积分与光通量对比度的统计关系
Table1.Statistical relationship between PSD integral and luminous flux contrast of instantaneous and integral focal spots.

表1可知, 对于四种不同方案, 瞬时焦斑范围内的PSD积分开方值与瞬时光通量对比度均近似相等, 和理论分析一致. 对于1D-SSD+CPP和RS+CPP两种束匀滑方式时间积分后, 焦斑范围内的PSD积分下降相对较缓, 而光通量对比度下降较快, 两者不再近似相等, 即积分时间后PSD积分与光通量对比度之间不再满足统计等价关系.
2
3.3.近场时间相关性与远场均匀性
-->为了进一步阐明瞬时和积分时间后PSD积分与光通量对比度之间的统计关系, 我们采用近场、远场光强分布的时间相关系数表征SSD+CPP和RS+CPP两种束匀滑方案在不同时刻焦斑光强的关联程度(图7).
图 7 不同束匀滑方案的近场、远场的时间相关特性 (a) 1D-SSD+CPP时间相关特性; (b) RS+CPP时间相关特性
Figure7. The near-field, far-field temporal and spatial correlation characteristics of different beam smoothing schemes: (a) Temporal correlation of 1D-SSD+CPP; (b) temporal correlation of RS+CPP.

图7可以看出, 1D-SSD+CPP与RS+CPP的近场光强时间相关性强, 几乎不随时间变化. 不同时刻的近场光强相关系数均近似为1, 即不同时刻的近场光强相似性很大, 可视为近场光强在不同时刻的分布基本不变. PSD等于近场光强的自相关, 且PSD曲线的面积随积分时间基本不变, 与近场光强时间相关性强相互验证.
在1D-SSD+CPP 方案中, 远场时间相关性随时间增加而迅速减小, 随后在小幅度范围内呈周期性变化, 即不同时刻远场光场的分布不尽相同, 因此, 低相关性的远场光强叠加可使远场光场在时间积分后分布更均匀, 其光通量对比度会随积分时间增大而减小. 而在RS+CPP方案中, 远场光强时间相关性随时间呈周期性变化, 且其周期与抽运光周期一致, 与子脉冲之间的延迟时间${T_{\rm{d}}}$相等, 积分时间后远场光场分布更均匀, 光通量对比度降低. 这一结果与积分时间后PSD积分和光通量对比度之间不再满足统计等价关系的结论相互验证.
为了寻求快速而简便的新方法来描述焦斑的特征, 提出了用于描述和表征靶面焦斑光场特性的统计表征方法. 采用圆型复数高斯随机变量描述靶面光强的统计分布, 且与衍射积分模型得到焦斑的FOPAI曲线和位相分布进行对比, 其统计特征均能较好符合, 说明了光场特性的统计表征方法能够用于描述靶面光强的统计特征. 通过合理的假设和公式推导, 得到了部分评价指标的解析表达式, 如焦斑PSD和焦斑RMS. 在此基础上, 分析了时间积分后的不同束匀滑方案下的远场特性, 得出了焦斑的功率谱密度和光通量对比度之间的统计关系, 说明了光场特性的统计表征方法能够很好地反映焦斑特征. 与衍射积分的方法相比, 光场特性的统计表征方法在数值计算过程中仅需根据CPP面形的统计特性直接获得靶面光场统计分布的解析表达式, 避开了从近场到远场的数值计算过程, 且无需对靶面光场的每个点进行数据处理, 简洁有效且无需大规模数据储存及处理.
相关话题/统计 方案 计算 光学 空间

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • Ce-La-Th合金高压相变的第一性原理计算
    摘要:采用第一性原理计算对Ce0.8La0.1Th0.1在高压下fcc-bct的结构相变、弹性性质及热力学性质进行了研究讨论.通过对计算结果的分析,发现了合金在压力下的相变规律,压强升高到31.6GPa附近时fcc相开始向bct相转变,到34.9GPa时bct相趋于稳定.对弹性模量的计算结果从另一角 ...
    本站小编 Free考研考试 2021-12-29
  • 内在缺陷与Cu掺杂共存对ZnO电磁光学性质影响的第一性原理研究
    摘要:采用基于自旋密度泛函理论的平面波超软赝势方法,研究了Cu掺杂ZnO(简称CuZn)与内在缺陷共存对ZnO电磁光性质的影响.结果表明,Cu是以替位受主的形式掺入的;制备条件对CuZn及内在缺陷的形成起至关重要的作用,富氧条件下Cu掺杂有利于内在缺陷的形成,且CuZn-Oi最易形成;相反在缺氧条件 ...
    本站小编 Free考研考试 2021-12-29
  • GaAs纳米线晶体结构及光学特性
    摘要:采用分子束外延技术在N-型Si(111)衬底上利用自催化生长机制外延砷化镓(GaAs)纳米线,对生长的纳米线进行扫描电子显微镜测试,纳米线垂直度高,长度直径均匀度好.对纳米线进行光致发光(photoluminescence,PL)光谱测试,发现低温10K下两个发光峰P1和P2分别位于1.493 ...
    本站小编 Free考研考试 2021-12-29
  • 螺旋锥束计算机断层成像倾斜扇束反投影滤波局部重建算法
    摘要:螺旋锥束计算机断层成像(CT)作为常用的临床诊断工具,如何尽可能地减少其辐射剂量是热点研究领域之一.局部成像利用准直器减小射线直照区域,能够有效降低CT辐射剂量.然而,局部成像会造成投影数据横向截断,产生局部重建问题.现有螺旋反投影滤波(BPF)算法只能实现局部曲面重建,难以实现局部体区域重建 ...
    本站小编 Free考研考试 2021-12-29
  • 基于瞬态光栅频率分辨光学开关法测量飞秒脉冲的研究
    摘要:超宽光谱的飞秒脉冲测量一直是超快激光领域的重要研究方向之一.常规的飞秒脉冲自相关方法是通过测量自相关倍频信号来获得,而倍频信号具有波长选择性,不同中心波长的飞秒脉冲测量需要更换不同的倍频晶体,十分不方便.因此,提出了一种改进型的瞬态光栅频率分辨光学开关(TG-FROG)方法用于测量飞秒脉冲.该 ...
    本站小编 Free考研考试 2021-12-29
  • 光声光谱检测装置中光声池的数值计算及优化
    摘要:利用光声光谱技术进行痕量气体的检测具有独特的优势,光声池是系统装置中最为重要的核心部件,它决定着整机性能的优劣.以一圆柱形共振型光声池为研究对象,基于声学与吸收光谱学的基本理论,建立了光声池声场激发的数学模型;利用数值模拟方法对光声池空腔结构进行了声学模态仿真,获得了前8阶声学模态值以及声压可 ...
    本站小编 Free考研考试 2021-12-29
  • 高温高压下立方氮化硼和六方氮化硼的结构、力学、热力学、电学以及光学性质的第一性原理研究
    摘要:本文采用基于密度泛函理论的第一性原理平面波赝势和局域密度近似方法,优化了立方和六方氮化硼的几何结构,系统地研究了零温高压下立方和六方氮化硼的几何结构、力学、电学以及光学性质.结构与力学性质研究表明:立方氮化硼的结构更加稳定,两种结构的氮化硼均表现出一定的脆性,而六方氮化硼的热稳定性则相对较差; ...
    本站小编 Free考研考试 2021-12-29
  • W-In体系溶质晶界偏聚行为的第一性原理计算
    摘要:基于第一性原理构建了钨基合金体系的溶质偏聚模型,以W-In体系为例研究了不同浓度下溶质的晶界偏聚行为和成键特征,从电子结构层面揭示了W-In体系的键合作用,预测了W-In体系界面稳定性随溶质浓度的变化规律.结合键布居、电荷密度、差分电荷密度和态密度等电子结构分析,发现了W-In体系中溶质原子在 ...
    本站小编 Free考研考试 2021-12-29
  • 含磷酸胍基间作用的磷酸双乙酸胍晶体电子结构与光学性质研究
    摘要:基于磷酸胍基间作用在L-精氨酸磷酸盐晶体特异性与生物化学功能中的重要角色,已制备了含有磷酸胍基的新晶体磷酸双乙酸胍.本文采用第一性原理计算了磷酸双乙酸胍晶体的电子结构与三个晶向上的光学性质,探讨了其中基团间作用与光学性质的关联.结果表明,磷酸双乙酸胍晶体能隙为4.77eV,远小于磷酸二氢钾晶体 ...
    本站小编 Free考研考试 2021-12-29
  • 100 MeV质子双环双散射体扩束方案设计
    摘要:为满足质子单粒子效应实验对大面积、均匀化质子束流的需求,针对中国原子能科学研究院100MeV质子回旋加速器提供的100MeV质子进行了双环双散射体扩束方案设计.Geant4模拟表明该方案可在2.4m位置产生一个均匀性为±1.89%、半径为8cm的照射野,在5m位置产生一个均匀性为±5.32%、 ...
    本站小编 Free考研考试 2021-12-29