Nonsequential double ionization of atoms driven by spatially inhomogeneous laser fields

Li Ying-Bin; Zhang Ke; Chen Hong-Mei; Kang Shuai-Jie; Li Zheng-Fa; Cheng Jian-Guo; Wu Yin-Meng; Zhai Chun-Yang; Tang Qing-Bin; Xu Jing-Kun; Yu Ben-Hai

doi:10.7498/aps.72.20230548

Abstract

Authors and contacts

FullText HTML

1. 引　言
2. 浅海全波动混响理论
2.1 浅海海底粗糙界面散射
2.2 浅海海底粗糙界面混响
3. 海底反射系数的三参数模型
4. 海底粗糙界面反向散射模型
4.1 角度特性
4.2 强度特性
5. 结　论

References

Cited By

Abstract

Using a three-dimensional classical ensemble method, we investigate the nonsequential double ionization (NSDI) of xenon atoms from the near infrared wavelength to the mid-infrared wavelength in spatially inhomogeneous laser fields, and compare the results with those from spatially homogeneous laser fields. The results show that the NSDI probability curves from spatially inhomogeneous laser field and spatially homogeneous laser field at short wavelength are similar to each other. With the laser wavelength increasing, NSDI at the high intensities is more and more suppressed for spatially inhomogeneous laser field.Compared with the result from the spatially homogeneous laser field, the final emission angle of two electrons from the NSDI exhibits a very strongly correlated characteristic in the spatially inhomogeneous field, especially at a longer laser wavelength, the final emission angles of two electrons are almost both concentrated around

${0^\circ }$ , meaning that the two electrons are always emitted into the same direction parallelly. Moreover, effective recollision of the NSDI is always dominated by first return of the first electron from the near infrared to the mid-infrared inhomogeneous laser fields, however, the transition from the first return dominance to the second return dominance occurs in the spatially homogeneous laser fields. Further, we reveal the more details of the ultrafast dynamics of the correlated electrons in the spatially inhomogeneous laser field by back-tracing the classical trajectories of NSDI.

Keywords:

PACS:

03.67.Dd	(Quantum cryptography and communication security)
03.65.Ud	(Entanglement and quantum nonlocality)
03.67.Hk	(Quantum communication)

Authors and contacts

College of Physics and Electronic Engineering, Xinyang Normal University, Xinyang 464000, China

Corresponding author: Xu Jing-Kun, 2777394244@qq.com ; Yu Ben-Hai, hnyubenhai@163. com

Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 12004323, 12074329, 12247208, 12104389) and the Nanhu Scholars Program for Young Scholars of Xinyang Normal University, China.

FullText HTML

1. 引　言

海底作为浅海混响的主要散射源, 浅海海底界面混响是很多学者探究的课题. 浅海混响根据其形成过程可以分为传播和散射两个过程. 入射声传播和散射声传播过程同属信道声传播问题, 研究较为完善. 而散射过程则比较复杂, 散射强度的测量也较难实现, 尤其是小掠射角的反向散射的测量, 所以针对散射的研究相对比较受限. 在混响研究伊始, 各国学者均采用经验的散射模型描述海底对声场的散射作用. 由于其形式较为简单, 运算速度较快, 在主动声纳预报系统中有很大的优势^[1-4], 至今仍被延续使用. 但是该类型的散射模型也存在着很大的局限性, 很难分析海底散射机理. 最初的经验散射模型是借鉴光学理论中提出的半无限自由空间中的Lambert散射定律给出的海底散射模型, 是海底散射强度与平面波掠射角之间的关系. 而在波导环境中, 由于频散效应的存在, 平面波散射理论不再适用于声场处理, 遂采用简正波理论, 避免了频散产生的多途现象. 然而简正模态的掠射角不同于平面波的掠射角, 因此对波导环境中的声散射问题不能够直接采用经验的Lambert散射模型. 近年来, 越来越多的学者开始从海底散射的物理机理着手建立散射模型. 国际上关于物理散射模型的建模方法包括有限元方法^[5]、Kirchhoff近似方法^[6]、微扰近似方法^[7-10]. 国内关于物理散射模型的研究最具代表性的是中国科学院声学研究所的高天赋和尚尔昌. Shang等^[11]在Gao^[12]和Tang^[13]的研究基础上提出了全波动混响理论, 并对该模型进行了一系列的发展. 2001年, Gao等^[14]提出了依据混响数据反演海底反向散射矩阵的方法. Wu等^[15,16]将海底反射模型引入到浅海混响模型中简化了混响衰减特性与海底沉积层参数之间的关系, 为海底地声参数的反演提供了一种新方法.

本文在尚尔昌和吴金荣的研究基础上, 对反向散射模型进行进一步分析. 通过引入海底反射系数的相移参数, 描述海底对声场的散射作用, 分析海底反向散射模型的角度特性和强度特性. 在实际应用允许的误差范围内, 进一步分析海底粗糙界面反向散射模型的角度特性和强度特性, 并与经验散射模型进行对比, 说明两者之间的异同.

2. 浅海全波动混响理论

在浅海波导环境中, 海底粗糙界面是海底混响的主要散射源. 在图1所示的浅海波导环境中, 平坦海底的平均海深为 ${H_0}$ , 任意水平位置 $r$ 的海底粗糙界面相对于平均海深的起伏高度为 $\eta (r)$ , 其量值远小于 ${H_0}$ , 均值为0, 方差为 $\sigma _\eta ^{\rm{2}}$ . 海底粗糙界面对声场的散射强度远小于入射声场的强度, 即满足弱散射的条件. 海底为半无限均匀介质, 其声速、密度和声吸收系数分别为 ${c_b},{\rho _b}\,{\text{和}}\,\alpha _b^{(w)}$ (声吸收系数的单位是 ${{{\rm{dB}}}/{\rm{\lambda }}}$ ). 海面为平坦的自由边界, 不考虑其对声场的散射作用.

$图 1 浅海波导环境\r\nFig. 1. Shallow water waveguide environment.$

图 1 浅海波导环境

Fig. 1. Shallow water waveguide environment.

下载: 全尺寸图片幻灯片

2.1 浅海海底粗糙界面散射

全波动混响理论的基础思想是将海底粗糙界面对声场的散射作用看作是“二次声源”向外辐射声能量, 从而将海底粗糙界面的散射作为“声源”向波导环境中辐射声能量^[11]. 该“声源”的声源强度与入射声强度成正比关系. 在水平均匀的波导环境中, 单位简谐点源入射声场满足波动方程为

${\nabla ^2}{u_{{\rm{w,b}}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}} + k_{{\rm{w,b}}}^{\rm{2}}{u_{{\rm{w,b}}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}} = 4{\rm{{\text π} \delta (}}{{{R}}_1} - {{{R}}_0}{\rm{)}},$

(1)

海底粗糙界面的“二次声源”满足

$\begin{split} & {\left. {\frac{{\partial {u_{\rm{w}}}({ R_1},{ R_0})}}{{\partial z}} - \frac{{\partial {u_{\rm{b}}}({ R_1},{ R_0})}}{{\partial z}}} \right|_{z = {H_0}}}\\ = & V({R_1}){G^{\rm{i}}}({ R_1},{ R_0}), \end{split}$

(2)

$\begin{split} & {\left. {{\rho _{\rm{w}}}{u_{\rm{w}}}({ R_1},{ R_0}) - {\rho _{\rm{b}}}{u_{\rm{b}}}({ R_1},{ R_0})} \right|_{z = {H_0}}}\\ = & p({ R_1}){G^{\rm{i}}}({ R_1},{ R_0}), \end{split}$

(3)

${G^{\rm{i}}}({{{R}}_1},{{{R}}_0}) = \sqrt {\frac{{2{\rm{\pi }}i}}{{{k_{\rm{w}}}r}}} \sum\limits_{m = 1}^M {{\phi _m}({z_0}){\phi _m}{\rm{(}}{H_0}{\rm{)}}{{\rm{e}}^{{\rm{i}}{k_m}r - {\beta _m}r}}} ,$

(4)

$V({{{R}}_1}) = \left(k_{\rm{w}}^{\rm{2}} - \frac{{k_{\rm{b}}^{\rm{2}}}}{\kappa }\right)\eta {\rm{(}}{r_1}{\rm{)}} + \left(1 - \frac{1}{\kappa }\right){\nabla _ \bot } \cdot {\rm{(}}\eta ({r_1}{\rm{)}}{\nabla _ \bot }{\rm{)}},$

(5)

$p{\rm{(}}{{{R}}_1}{\rm{)}} = {\rm{(}}{\rho _{\rm{w}}} - {\rho _{\rm{b}}}{\rm{)}}\eta {\rm{(}}{r_1}{\rm{)}}\frac{\partial }{{\partial z}},$

(6)

${k_{\rm{w}}} = \frac{{2{{\text π}}f}}{{{c_{\rm{H}}}}},\;\;\;\;\;\;{k_{\rm{b}}} = \frac{{2{{\text π}}f}}{{{c_{\rm{b}}}}},\;\;\;\;\;\;\kappa = \frac{{{\rho _{\rm{b}}}}}{{{\rho _{\rm{w}}}}},$

(7)

其中, ${{{R}}_0} = {\rm{(}}{r_{\rm{s}}},{z_0}{\rm{)}}$ , ${{{R}}_1} = ({r_1},{H_0})$ 分别是声源和散射点的空间位置, 两者之间的水平距离差为 $r = |{r_{\rm{s}}} - {r_1}|$ ; $u_{\rm{w}}, u_{\rm{b}}$ 分别是水介质和沉积层介质中“二次源”的辐射声场(下脚标w和b分别代表水介质和沉积层介质); ${k_{{\rm{w,b}}}}$ 是水介质和沉积层介质中的波数, 与声源频率 $f$ 和介质声速有关; ${c_{\rm{H}}} =$ ${c_{\rm{w}}}{\rm{(}}{H_0}{\rm{)}}$ 是海深处的声速, 沉积层介质和水介质的密度比用 $\kappa$ 表示; ${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ 是波导环境中初级声场的格林函数. 根据简正波理论的思想, 在浅海波导环境中, ${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ 可以表示为简正模态 ${\phi _m}({H_0})$ (实际上是 ${\phi _m}(z)$ , $z$ 是接收深度, 在这里的接收深度为海深 ${H_0}$ )叠加的形式, 如(4)式所示, k_m和 ${\beta _m}$ 分别是简正模态本征值的实部和虚部. $V({{{R}}_1})$ 和 $p({{{R}}_1})$ 则分别是质点振速算子和压力算子, 作用到初级声场的格林函数 ${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ , 得到海底粗糙界面的“二次速度源”的声源强度 $V{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ 和“二次压力源”的声源强度 $p{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)}}$ . 在瑞利参数较小的情况下, Bass给出了(5)式和(6)式的近似表达式. ${\nabla _ \bot }$ 为水平梯度算子,

${\nabla _ \bot } = \frac{\partial }{{\partial r}}.$

(8)

根据格林定理, 水平均匀波导环境中, “二次速度源”在接收点 ${{R}} = {\rm{(}}{r_{\rm{r}}},z{\rm{)}}$ 处的辐射声场可以通过格林函数 $G_{{\rm{w,b}}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$ 给出形式解. 辐射声场的格林函数满足

$\begin{aligned} & {\nabla ^2}G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} + k_{\rm{w}}^{\rm{2}}G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} = 0\;\;\;z < {H_0}\\ & {\nabla ^2}G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} + k_{\rm{b}}^2G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} = 0\;\;\;z > {H_0} \end{aligned},$

(9)

同时满足边界条件

$\begin{gathered} {\left. {\frac{{\partial G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}}}{{\partial z}} - \frac{{\partial G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}}}{{\partial z}}} \right|_{z = {H_0}}} = {\rm{\delta (}}{{R}} - {{{R}}_1}{\rm{)}} \\ {\left. {\rho G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} - {\rho _{\rm{b}}}G_{\rm{b}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}} \right|_{z = {H_0}}} = 0 \\ \end{gathered} .$

(10)

在接收点 “二次速度源”的辐射声场满足形式解,

$u_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} = \int {G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}V{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)d}}{{{R}}_1}} .$

(11)

“二次压力源”的辐射声场同样有(11)式的形式解,

$u_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} = \int {G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}p{\rm{(}}{{{R}}_1}{\rm{)}}{G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{{R}}_0}{\rm{)d}}{{{R}}_1}} ,$

(12)

其中, $G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$ 是“二次压力源”的辐射声场的格林函数. 根据边界处法向振速与压力之间的关系, 不难得到 $G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$ 与 $G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}$ 满足

$G_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}} = - \frac{1}{{{\rho _{\rm{w}}}}}\frac{\partial }{{\partial z}}G_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1}{\rm{)}}.$

(13)

叠加“二次速度源”和“二次压力源”的辐射声场, 得到海底粗糙界面向波导环境中辐射的总声场

${u_{\rm{w}}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} = u_{\rm{w}}^{\rm{V}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}} + u_{\rm{w}}^{\rm{p}}{\rm{(}}{{R}},{{{R}}_1},{{{R}}_0}{\rm{)}}.$

(14)

将(11)—(13)式代入(14)式, 得到海底粗糙界面散射声场的一般形式解

$\begin{split} & {u_{\rm{w}}}( R,{ R_1},{ R_0})\\ = & \int [G_{\rm{w}}^{\rm{V}}( R,{ R_1})V({ R_1})\\ & - \frac{1}{{{\rho _{\rm{w}}}}}\frac{\partial }{{\partial z}}G_{\rm{w}}^{\rm{V}}( R,{ R_1})p({ R_1})]{G^{\rm{i}}}({ R_1},{ R_0}){\rm{d}}{ R_1}. \end{split}$

(15)

实际上, 入射声场的格林函数、散射场的格林函数以及初级声场的格林函数有一致的形式解, 忽略他们之间的差异,

${G^{\rm{i}}}{\rm{(}}{{{R}}_1},{{R}}{\rm{)}} = G_{\rm{w}}^{\rm{V}}{\rm{(}}{{{R}}_1},{{R}}{\rm{)}} = G{\rm{(}}{{{R}}_1},{{R}}{)},$

(16)

所以对(15)式进行简化得到

$\begin{split} & {u_{\rm{w}}}( R,{ R_1},{ R_0})\\ = & \int_S {[G( R,{ R_1})V({ R_1})} \\ & - \frac{1}{{{\rho _{\rm{w}}}}}\frac{{\partial G( R,{ R_1})}}{{\partial z}}p({ R_1})]G({ R_1},{ R_0}){\rm{d}}{ R_1}. \end{split}$

(17)

对于本地混响, 声源和接收均在同一水平位置. 为了方便计算, 通常令其位于过源点的垂直轴线上, 即 ${r_{\rm{s}}} = {r_{\rm{r}}} = 0$ . 将(4)—(6)式代入(17)式, 并通过分部积分法得到简谐点源的海底粗糙界面的反向散射声场, 将其描述为简正模态的叠加形式^[11],

$\begin{split} & {u_{\rm{w}}}(r,z)\\ = & \displaystyle\frac{{2\pi i}}{{{k_{\rm{w}}}r}}\sum\limits_m^M {\sum\limits_n^M {{\phi _m}({z_0}){\phi _m}({H_0})} } \cdot C_{mn}^\eta K_{mn}^\eta ({k_m},{k_n})\\ & \times {\phi _n}({H_0}){\phi _n}(z){{\rm{e}}^{ - ({\beta _m} + {\beta _n})r}} \end{split}$

(18)

其中

$C_{mn}^\eta = k_{\rm{w}}^{\rm{2}} - \frac{{k_{\rm{b}}^2}}{\kappa } + {\rm{(}}1 - \frac{1}{\kappa }{\rm{)}}{k_m}{k_n} + \frac{{1 - \kappa }}{{{\kappa ^2}}}{\gamma _m}{\gamma _n},$

(19)

$K_{mn}^\eta {\rm{(}}{k_m},k{}_n{\rm{)}} = \int {\eta {\rm{(}}r'{\rm{)}}{{\rm{e}}^{{\rm{i(}}{k_m} + {k_n}{\rm{)}}r'}}{\rm{d}}r'} .$

(20)

2.2 浅海海底粗糙界面混响

实际上(18)式是平坦波导环境中海底粗糙界面散射场的稳定解形式. 对于声源脉冲信号s(t), 其功率谱 $F\left( \omega \right)$ 满足Fourie变换

$F{\rm{(}}\omega {\rm{)}} = \int_{ - \infty }^{ + \infty } {s{\rm{(}}t{\rm{)}}{{\rm{e}}^{ - {\rm{i}}\omega t}}{\rm{d}}t},\tag{21a}$

(21a)

$s{\rm{(}}t{\rm{)}} = \frac{1}{{2{\rm{\pi }}}}\int_{ - \infty }^{ + \infty } {F{\rm{(}}\omega {\rm{)}}{{\rm{e}}^{{\rm{i}}\omega t}}{\rm{d}}\omega }.\tag{21b}$

(21b)

根据Fourie变换的性质, 得到入射信号为s(t)的海底散射声场为

$u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{(}}r,z{\rm{)}} = \frac{1}{{2{\rm{\pi }}}}\int_{ - \infty }^{ + \infty } {F{\rm{(}}\omega {\rm{)}}{u_{\rm{w}}}{\rm{(}}r,z{\rm{)}}{{\rm{e}}^{{\rm{i}}\omega t}}{\rm{d}}\omega } .$

(22)

将声源信号的功率谱函数 $F\left( \omega \right)$ 和(18)式代入(22)式, 同时取水平波束在中心频率处的两阶泰勒级数展开, 进而得到脉宽为 $\tau$ 声源信号s(t)的反向散射声场,

$\begin{split} & u_{\rm{w}}^{s(t)}(r,z)\\ = & \frac{{2\pi i}}{{{k_{\rm{w}}}r}}s(t - {t_0})\sum\limits_m^M {\sum\limits_n^M {{\phi _m}({z_0}){\phi _m}({H_0}) } } \\ & \times C_{mn}^\eta K_{mn}^\eta \cdot {\phi _n}({H_0}){\phi _n}(z){{\rm{e}}^{ - ({\beta _m} + {\beta _n})r}} \end{split}$

(23)

其中, t₀是信号的传播时间, 与散射环的水平距离满足 ${t_0} \approx {{2r}/{{c_{\rm{H}}}}}$ .

混响是能够同一时刻被接收、来自各方向散射回波的叠加. 对于海底粗糙界面混响而言, 同一时刻接收到的散射回波来自于宽度为 $\Delta r = {{{c_{\rm{H}}}\tau }/2}$ 的散射圆环, 其内径为 ${r_{{\rm{in}}}} = r - \displaystyle\frac{{{c_{\rm{H}}}\tau }}{4}$ , 外径为 ${r_{{\rm{ex}}}} = r +$ $\displaystyle\frac{{{c_{\rm{H}}}\tau }}{4}$ . 对(23)式计算强度, 并在散射环面积上进行积分得到海底粗糙界面混响平均强度

${I_\eta }{\rm{(}}r,z{\rm{)}} = \int_A {\langle u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{[}}u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{]}}^{*}\rangle } {\rm{d}}A,$

(24)

其中, ${\rm{[}}u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{]}}^*$ 是 $u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}$ 的复共轭.

假设散射环面积内的海底界面起伏高度满足各向同性的原则, 对(24)式进一步简化为

${I_\eta }{\rm{(}}r,z{\rm{)}} = 2{\rm{\pi }}r\int_{{r_{{\rm{in}}}}}^{{r_{{\rm{ex}}}}} {\langle u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{[}}u_{\rm{w}}^{s{\rm{(}}t{\rm{)}}}{\rm{]}}^{*}\rangle {\rm{d}}r} .$

(25)

将(23)式代入(25)式, 并假设脉宽远小于混响时间t, 将声源强度近似为

${E_0} = {s^2}\left( {t - {t_0}} \right) \cdot \tau .$

(26)

在这里, 只考虑其非相干特性, 得到短脉冲的非相干混响平均强度

$\begin{split} {I_\eta }(r,z) = & {E_0} \cdot {\left(\frac{{2{\text{π}}}}{{{k_{\rm{w}}}r}}\right)^2}\\ & \times {\text{π}} rc\sum\limits_m^M {\sum\limits_n^M {\phi _m^2({z_0})\phi _m^2({H_0}) \cdot S_{mn}^\eta } } \\ & \times \phi _n^2({H_0})\phi _n^2(z){{\rm{e}}^{ - 2({\beta _m} + {\beta _n})r}}, \end{split}$

(27)

$S_{mn}^\eta = {\left[ {C_{mn}^\eta } \right]^2}{S^\eta }\left( {2{k_0}} \right),$

(28)

$\begin{split} & {S^\eta }(2{k_0})\\ = & \int_{{r_{{\rm{in}}}}}^{{r_{{\rm{ex}}}}} {\int_{{r_{{\rm{in}}}}}^{{r_{{\rm{ex}}}}} {\langle \eta (r')\eta (r'')\rangle {{\rm{e}}^{{\rm{i}}({k_m} + {k_n})r' - {\rm{i}}({k_m} + {k_n})r''}}} } \\ & \times {\rm{d}}r''{\rm{d}}r'\; = \;\sigma _\eta ^{\rm{2}}{P^\eta }(2{k_0}), \end{split}$

(29)

其中, ${S^\eta }(2{k_0})$ 是海底界面的粗糙度谱, 与海底沉积层的地声参数无关, 只是声源频率的函数; $S_{mn}^\eta$ 是海底粗糙界面引起的入射模态和散射模态之间的耦合系数, 反映海底粗糙界面的反向散射能力, 由海底界面粗糙度谱和地声参数决定. 定义

$\Theta _{mn}^\eta = \phi _m^2({H_0})S_{mn}^\eta \phi _m^2({H_0})$

(30)

是区别于经验散射模型的物理散射模型. 它结合了波导环境特性和简正波理论的思想, 建立了受格林函数约束的散射模型, 明确了海底地声参数以及海底粗糙界面与海底反向散射的定量关系. 因此, 地声参数的准确度直接影响到海底反向散射函数的精确度. 通常情况下, 很难直接获取大面积的地声数据, 通过反演得到的地声参数又受地声模型的影响较大, 而且反演的未知参数较多, 导致结果存在很大的不确定性. 从反射系数的角度思考, 海底对声场的影响主要体现在海底反射系数. 海底反射系数的三参数模型与海底地声模型无关, 且参数较少, 利用海底反射系数代替地声参数作为输入参数, 能很大程度上简化海底粗糙界面的反向散射模型.

3. 海底反射系数的三参数模型

尚尔昌在1979年提出了海底反射系数的三参数模型^[17], 随后给出了海底小掠射角的反射系数幅值和相移与掠射角之间的量化关系, 该参数与地声模型无关. 海底反射系数的幅值用Q参数表示

$\ln \left| {V\left( \theta \right)} \right| = - Q\theta ,\;\;\theta \to 0,$

(31)

海底反射系数的相移用P参数表示

$\arg \left[ {V\left( \theta \right)} \right] = -{\text{π}} + P/\theta ,$

(32)

$\theta$ 是海底界面处入射声波的掠射角; Q参数反映海底声吸收效果, 体现声能量的衰减; P参数则描述各阶简正模态在海深处的能量, 即控制各阶简正模态的水平波数. 海底反射系数的相移参数P与地声参数之间的转换关系满足^[18,19]

${P^{{\rm{Low}}}} = \frac{{2\kappa }}{{\sqrt {1 - {\upsilon ^{ - 2}}} }},\;\;{\rm{ }}\theta < \theta ^{\rm{*}},$

(33)

${{{P^{{\rm{Crit}}}} = {\text{π}}}/{\theta ^{*}}},\;\;{\rm{ }}\theta \to \theta ^{\rm{*}},$

(34)

其中, $\upsilon = {{{c_{\rm{b}}}}/{{c_{\rm{H}}}}}$ 表示沉积层介质和水介质的声速比; $\kappa$ 是密度比, 如(7)式所示; $\theta ^{*}$ 是临界掠射角.

(33)式和(34)式分别给出了不同海底掠射角时海底反射系数的相移参数的不同形式. ${P^{{\rm{Low}}}}$ 是甚小掠射角处的海底反射系数相移参数, 是海底沉积层和水介质的声速比和密度比的函数; ${P^{Crit}}$ 是临界角附近海底反射系数的相移参数, 是临界角的函数.

海底反射系数的幅值参数Q与地声参数之间的转换关系满足^[18,19]

$Q = 0.036 \cdot \frac{{{\upsilon ^{ - 1}}\kappa }}{{{{\left[ {1 - {\upsilon ^{ - 2}}} \right]}^{3/2}}}} \cdot \alpha _b^{\left( w \right)},$

(35)

其中, $\alpha _b^{\left( w \right)}$ 是海底沉积层的声吸收系数, 单位是 $dB/\lambda$ .

根据Snell折射定律, $\upsilon$ 与临界掠射角之间满足

${\rm{cos(}}\theta ^{*}{\rm{)}} = {{{c_{\rm{H}}}}/{{c_{\rm{b}}}}} = {\upsilon ^{ - 1}}.$

(36)

将(36)式代入(34)式, 临界角附近的海底反射系数的相移参数 ${P^{Crit}}$ 与声速之间满足

$\upsilon = {\rm{co}}{{\rm{s}}^{ - 1}}{\rm{(}}\frac{{\rm{\pi }}}{{{P^{{\rm{Crit}}}}}}{\rm{)}}.$

(37)

联合(37)式和(33)式得到甚小掠射角的海底反射系数的相移参数 ${P^{Low}}$ 与密度之间满足关系式

$\kappa = \frac{{{P^{{\rm{Low}}}}}}{2}\sin \left(\frac{{\text{π}}}{{{P^{{\rm{Crit}}}}}}\right).$

(38)

(37)式和(38)式分别是海底反射系数的参数与等效的均匀半无限海底介质的地声参数之间的转换关系.

4. 海底粗糙界面反向散射模型

常用的经验散射模型如(39)式所示:

${\Theta _E} = {\mu _E}{\sin ^l}\theta {\sin ^k}\varphi ,$

(39)

其中, ${\mu _{\rm{E}}}$ 是经验散射系数, 反映反向散射强度; $\theta$ 是入射掠射角; $\varphi$ 是散射掠射角. 当 $l = k = 1$ , 是经验的Lambert散射模型. 该类型的散射模型是可分离的经验海底反向散射模型, 不可分离的经验散射模型(参见附录)也被用来描述海底反向散射.

(30)式的海底反向散射模型是根据波导环境建立的散射模型, 明确了入射简正模态、散射简正模态以及海底地声参数对反向散射的强度特性和角度特性的影响. 浅海波导环境中简正模态是受格林函数严格约束的稳定解, 与海底沉积层的介质参数和水文环境有关. 在掠射角小于临界角的条件下, 忽略各阶简正模态的水平波数(掠射角)之间的差异, 即 ${k_m} \approx {k_n} \approx {k_{\rm{w}}}$ , 则(19)式可以近似描述为与模态无关的量.

$\begin{split} C_{mn}^\eta = & \left[1 \!- \!\frac{{k_{\rm{b}}^{\rm{2}}}}{{k_{\rm{w}}^2\kappa }} \!+\! \left(1 - \frac{1}{\kappa }\right)\frac{{{k_m}{k_n}}}{{k_{\rm{w}}^2}} + \frac{{1 - \kappa }}{{{\kappa ^2}}}\frac{{{\gamma _m}{\gamma _n}}}{{k_{\rm{w}}^2}}\right]k_{\rm{w}}^2 \\ & \approx \;\left[1 - \frac{{k_{\rm{b}}^{\rm{2}}}}{{k_{\rm{w}}^{\rm{2}}\kappa }} + \left(1 - \frac{1}{\kappa }\right) + \frac{{1 - \kappa }}{{{\kappa ^2}}}\frac{{k_{\rm{w}}^2 - k_{\rm{b}}^{\rm{2}}}}{{k_{\rm{w}}^{\rm{2}}}}\right]k_{\rm{w}}^{\rm{2}} \\ & = \;\left[2 - 2{\kappa ^{ - 1}} + {\kappa ^{ - 2}} - {\left(\kappa \upsilon \right)^{ - 2}}\right]k_{\rm{w}}^2. \end{split}$

(40)

将(37)式和(38)式代入(40)式, 用海底反射系数的相移参数代替地声参数, 得到

$C_{mn}^\eta = k_{\rm{w}}^2 \cdot \varsigma {\rm{(}}P{\rm{)}},$

(41)

其中

$\varsigma {\rm{(}}P{\rm{)}} \approx 2 - \frac{4}{{P{\rm{sin(}}{{\text{π}}/P}{\rm{)}}}} + \frac{4}{{{P^2}}}.$

(42)

在这里忽略 ${P^{Low}}$ 和 ${P^{Crit}}$ 之间的差异, 文献[18]中也说明了这种近似的合理性.

同时, 海底反射系数的相移参数 $P$ 决定了各阶简正模态在海底粗糙界面平均深度处的能量, 即

${\phi _m}({H_0}) = \sqrt {\frac{2}{{{H_0}}}} \sin \left(\frac{{P{\theta _m}}}{2}\right).$

(43)

将(41)式和(43)式代入(30)式得到远距离(小掠射角)条件下的海底反向散射核函数,

${\Theta _P} = \frac{4}{{H_0^2}}S(2{k_0}){[k_w^2 \cdot \varsigma (P)]^2} \cdot {\sin ^2}\left(\frac{{P\theta }}{2}\right){\sin ^2}\left(\frac{{P\varphi }}{2}\right),$

(44)

其中 $\theta = {\theta _m}$ 为入射掠射角, $\varphi = {\theta _n}$ 为散射掠射角.

采用海底反射系数的相移参数 $P$ 描述海底反射系数, 明确了海底对声散射的物理机理, 不同于经验反向散射模型.

4.1 角度特性

在远距离条件下, (44)式表明, 海底反向散射的角度特性受海底反射系数的相移参数 $P$ 的影响, 不同于经验散射模型. 在临界角附近,

$\mathop {\lim }\limits_{\theta \to \theta *} \frac{{P\theta }}{2} = \frac{{\text{π}}}{2}.$

(45)

以 $\displaystyle\frac{{P\theta }}{2}$ 为宗量, 不能对 $\sin \left(\displaystyle\frac{{P\theta }}{2}\right)$ 进行小角度近似, 所以在临界角附近, (44)式所描述的反向散射模型随掠射角并不满足线性变化的关系.

以 ${{P\theta }/2} = 1$ 为分界点, 当 ${{P\theta }/2} < 1$ 时, 为甚小掠射角范围, 小角度近似引起的海底反向散射强度小于3 dB, 在实际应用中可忽略不计; 当 ${{P\theta }/2} > 1$ 时, 为临界角附近, 小角度近似使海底反向散射强度出现大于3 dB的误差, 实际应用中不可忽略. 所以, 在随后的分析中, 对反向散射模型的强度特性和角度特性进行分段处理.

当 ${{P\theta }/2} > 1$ , 即 ${2/P} < \theta < \theta *$ 时, (44)式不能进行小角度近似, 其随角度的变化关系受P参数加权.

$\begin{split}\Theta _P^{Crit} = & \mu _P^{Crit}{\sin ^2}\left(\frac{{P\theta }}{2}\right){\sin ^2}\left(\frac{{P\varphi }}{2}\right) \\ & \propto {\sin ^2}\left(\frac{{P\theta }}{2}\right){\sin ^2}\left(\frac{{P\varphi }}{2}\right),\end{split}$

(46)

其中

$\mu _P^{Crit} = \frac{4}{{H_0^2}}S(2{k_0}){[k_0^2 \cdot \varsigma (P)]^2}.$

(47)

当 ${{P\theta }/2} < 1$ , 即 $\theta < {2/P}$ 时, 对(44)式的角度项进行小角度近似,

$\sin \left(\frac{{P\theta }}{2}\right) \approx \frac{{P\theta }}{2}.$

(48)

将(48)式代入(44)式得到

$\Theta _P^{Low} = \mu _P^{Low} \cdot {\theta ^2}{\varphi ^2} \propto {\theta ^2}{\varphi ^2},$

(49)

其中

$\mu _P^{Low} = \frac{{{P^4}}}{{4H_0^2}}S(2{k_0}){[k_0^2 \cdot \varsigma (P)]^2}.$

(50)

对(39)式的经验散射模型进行小角度近似得到

${\Theta _E} \approx {\mu _E}{\theta ^l}{\varphi ^k} \propto {\theta ^l}{\varphi ^k}.$

(51)

当 $l = k = 2$ 时, 比较(51)式和(49)式具有相同的角度关系. 所以当掠射角 $\theta < {2/P}$ 时, 基于物理散射机理的反向散射模型与 $l = k = 2$ 的可分离经验散射模型具有一致的角度特性.

4.2 强度特性

通过4.1节的分析, 掠射角不同时, 海底反向散射的角度特性不同, 同时也反映了海底反向散射系数的差异.

(47)式指出, 在临界角附近, 海底反向散射系数由海底反射系数的相移参数P、海底界面粗糙度谱, 海深以及声源信号的频率决定. 忽略临界角附近以及甚小掠射角的P参数之间的差异^[18], 则(42)式中等式右边第三项相比前两项为小量, 可以忽略不计. 第二项中采用小角度近似,

${\rm{sin}}\frac{{\text{π}}}{P} \approx \frac{{\text{π}}}{P},$

(52)

使得(42)式中的第二项近似与P参数无关, 所以 $\varsigma (P) \approx 0.73$ . 图2也表明了 $\varsigma (P)$ 随P参数的变化不是很明显. 所以, 在临界角附近的海底反向散射系数近似与海底沉积层介质无关, 只是海底界面粗糙度谱、海深以及声源频率的函数.

$\mu _P^{Crit} \approx \frac{{2.11}}{{H_0^2}}S(2{k_0})k_0^4.$

(53)

$图 2 $\zeta $随P参数的变化\r\nFig. 2. $\zeta $ varied with P parameters near critical angle.$

图 2

$\zeta$

随P参数的变化

Fig. 2.

$\zeta$

varied with P parameters near critical angle.

下载: 全尺寸图片幻灯片

对于甚小掠射角的海底反向散射系数满足(50)式. 结合 $\varsigma (P)$ 近似为常数的结果, (50)式可以近似为与 ${P^4}$ 成正比, 即

$\mu _P^{Low} \approx {0.13}\frac{{{P^4}}}{{H_0^2}}S(2{k_0})k_0^4 \propto {P^4}.$

(54)

图3以第一类海底为例说明了这种近似结果的合理性. 图中给出了甚小掠射角的海底反向散射系数随 $P$ 参数的变化曲线(仿真参数: 海深为50 m, 声源频率为600 Hz, 海底粗糙界面的标准差为 $0.1\;{\rm{m}}$ , 相关尺度为 ${\rm{10}}\;{\rm{m}}$ , 计算得到Goff-Jordan谱^[20]为－32.8 dB)以及对其进行 $P$ 参数的四次方拟合结果, 两者的变化趋势基本一致. 所以, 在甚小掠射角的条件下, 海底反向散射系数与海底反射系数的相移参数 $P$ 密切相关, 即与海底沉积层的声速比和密度比有关, 与声吸收系数无关.

$图 3 海底甚小掠射角反向散射系数对P参数的依赖\r\nFig. 3. Relationship between P and bottom backscattering coefficient at very low grazing angle.$

图 3 海底甚小掠射角反向散射系数对P参数的依赖

Fig. 3. Relationship between P and bottom backscattering coefficient at very low grazing angle.

下载: 全尺寸图片幻灯片

另外, 从混响平均强度的角度分析, 同样说明了这种近似的合理性. 海底界面粗糙度谱函数采用图3的仿真参数, 海底是第一类均匀半无限介质, 声速为 $18{\rm{36}}\;{\rm{m/s}}$ , 密度为 ${\rm{2}}{\rm{.03}}\;{\rm{g/c}}{{\rm{m}}^{\rm{3}}}$ , 声吸收系数为 $0.88\;{\rm{dB/}}\lambda$ , 海深为50 m, 水介质为1500 m/s的等声速剖面的均匀水体. 声源频率为600 Hz, 发射深度为10 m, 接收深度为30 m. 计算水平散射距离从5 km到50 km的混响平均强度衰减曲线. 图4中实线是根据(30)式仿真的非近似混响平均强度衰减曲线; “○”是根据(46)式和(53)式仿真的临界角的近似结果, 与非近似的混响平均强度衰减曲线在13 km以前完全重合; 随着水平距离的增加, 两者相差逐渐增大; “ $\bullet$ ”是根据(49)式和(54)式仿真的甚小掠射角的近似结果, 在20 km以前, 与非近似结果之间的差异随水平距离的增大而减小, 在20 km以后两者基本吻合. 在散射距离较近时, 掠射角相对较大, 临界角附近近似的海底反向散射模型与真实结果更接近; 在散射距离较远时, 掠射角相对较小, 甚小掠射角的近似结果与真实结果更接近.

$图 4 海底反向散射模型的比较\r\nFig. 4. Compare with bottom backscattering model with different approximate.$

图 4 海底反向散射模型的比较

Fig. 4. Compare with bottom backscattering model with different approximate.

下载: 全尺寸图片幻灯片

5. 结　论

在浅海混响平均强度模型中, 经验散射模型在混响特性分析中存在明显的局限性, 而现有的物理散射模型受地声模型的影响较大. 本文以全波动混响理论的物理散射模型为基础, 结合海底反射系数的三参数模型, 将海底反射系数的相移参数等效代替地声参数, 描述海底对声场的散射作用, 简化了海底反向散射模型. 通过理论分析, 明确了海底掠射角以 ${2/P}$ 为分界点, 海底反向散射的角度特性和强度特性对海底反射系数的相移参数存在不同的依赖关系. 在掠射角满足 $\theta < {2/P}$ 时, 海底反向散射的角度特性近似描述为与可分离经验散射模型的角度特性一致, 与 $P$ 参数无关; 而其散射系数则近似描述为与P⁴线性增强. 当海底掠射角满足 ${2/P} < \theta < \theta ^{*}$ 时, 海底反向散射的角度特性是受 $P$ 参数加权的, 即受海底地声参数的影响; 而其散射系数则近似为与 $P$ 参数无关. 所以不同掠射角范围, 海底对反向散射声场的强度特性和角度特性的贡献不同. 掠射角较大时海底对反向散射声场的影响主要体现在其角度特性; 掠射角非常小时, 海底的影响主要体现在其强度特性.

附录A　经验反向散射模型

经验反向散射模型有多种形式, 主要体现在其角度特性之间的差异. 最常用的海底散射模型是可分离的散射模型,

$\Theta = {\mu _{\rm{E}}}{\rm{si}}{{\rm{n}}^l}\theta {\rm{si}}{{\rm{n}}^k}\varphi, \tag{A1}$

(A1)

其中, ${\mu _{\rm{E}}}$ 是海底反向散射系数; $l,k$ 取不同的值对应不同的散射模型. 常见的Lambert散射模型则是(A1)式中 $l = k = 1$ 时的散射模型. 该模型最初是Mackenzie^[1]在处理深海海底混响模型时由光学的散射原理引入, 给出经验的海底反向散射模型, 随后根据实验数据与理论结果对比得到 ${\mu _{\rm{E}}}$ 的分贝值大约在－27 dB左右, 也就是说 ${\mu _{\rm{E}}} \approx {10^{ - 2.7}}$ . 当 $l = 1,k = 0$ 表示散射源强正比于 $\mu {\rm{sin}}\theta$ 的均匀散射; 当 $l = 0,k = 0$ 表示与角度无关的均匀散射; 同时也存在 $l = k = 2$ 的反向散射模型.

在海底界面大尺度不均匀波导环境中, 小掠射角和基尔霍夫近似的条件下的海底散射模型可以表示为

$\Theta = {\mu _{\rm{E}}}{\rm{si}}{{\rm{n}}^l}\left(\frac{{\theta + \varphi }}{2}\right) \;\;l = 0,\frac{1}{2},1,\frac{3}{2},....\tag{A2}$

(A2)

在海底界面小尺度不均匀波导环境中, 小掠射角和基尔霍夫近似的条件下的海底散射模型可以表示为

$\Theta = {\mu _{\rm{E}}}{\rm{si}}{{\rm{n}}^l}\left[{\rm cos}^{- 1}\left(\frac{{{\rm{cos}}\theta + {\rm{cos}}\varphi }}{2}\right)\right]\;\;l = 0,\frac{1}{2},1,\frac{3}{2},....\tag{A3}$

(A3)

另外一种不可分离的散射模型, 同样被用来描述海底反向散射过程

$\Theta = {\mu _{\rm{E}}}\frac{{{\rm{sin}}\theta {\rm{sin}}\varphi }}{{{\rm{sin}}\theta + {\rm{sin}}\varphi }}.\tag{A4}$

(A4)

通常无论是可分离的经验反向散射强度还是不可分离的经验反向散射强度, 它们的共同点在于海底反向散射强度特性与角度特性相互分离, 相互独立.

References

[1]	Lein M, Hay N, Velotta R, Marangos J P, Knight P L 2002 Phys. Rev. Lett. 88 183903 Google Scholar
[2]	Ghimire S, Reis D A 2019 Nat. Phys. 15 10
[3]	Yang B, Schafer K J, Walker B, Kulander K C, Agostini P, DiMauro L F 1993 Phys. Rev. Lett. 71 3770 Google Scholar
[4]	Milošević D B, Paulus G G, Becker W 2003 Opt. Express 11 1418 Google Scholar
[5]	Busuladi M, Gazibegovi-Busuladi A, Miloevi D B, Becker W 2008 Phys. Rev. Lett. 100 203003 Google Scholar
[6]	He M R, Li Y, Zhou Y M, Li M, Cao W, Lu P X 2018 Phys. Rev. Lett. 120 133204 Google Scholar
[7]	Zhou Y, Tolstikhin O I, Morishita T 2016 Phys. Rev. Lett. 116 173001 Google Scholar
[8]	Rudenko A, Jesus V, Ergler T, Zrost K, Ullrich J 2007 Phys. Rev. Lett. 99 263003 Google Scholar
[9]	Quan W, Hao X L, Wang Y L, Chen Y J, Yu S G, Xu S P, Xiao Z L, Sun R P, Lai X Y, Hu S L, Liu M Q, Shu Z, Wang X D, Li W D, Becker W, Liu X J, Chen J 2017 Phys. Rev. A 96 032511 Google Scholar
[10]	Huang C, Zhou Y M, Zhang Q B, Lu P X 2013 Opt. Express 21 11382 Google Scholar
[11]	Tong A H, Liao Q, Zhou Y M, Lu P X 2010 Opt. Express 18 9064 Google Scholar
[12]	Shaaran T, Augstein B B, Figueira D 2011 Phys. Rev. A 84 013429 Google Scholar
[13]	Wang Y L, Xu S P, Chen Y J, Kang H P, Lai X Y, Quan W, Liu X J, Hao X L, Li W D, Hu S L, Chen J, Becker W, Chu W, Yao J P, Zeng B, Cheng Y, Xu Z Z 2017 Phys. Rev. A 95 063415 Google Scholar
[14]	Corkum P B 1993 Phys. Rev. Lett. 71 1994 Google Scholar
[15]	Feuerstein B, Moshammer R, Fischer D, Dorn A, Schroter C D, Deipenwisch J, Crespo Lopez-Urrutia J R, Hohr C, Neumayer P, Ullrich J, Rottke H, Trump C, Wittmann M, Korn G, Sandner W 2001 Phys. Rev. Lett. 87 043003 Google Scholar
[16]	Wang Y L, Xu S P, Quan W, Gong C, Lai X Y, Hu S L, Liu M Q, Chen J, Liu X J 2016 Phys. Rev. A 94 053412 Google Scholar
[17]	Huang C, Zhong M M, Wu Z M 2016 Opt. Express 24 28361 Google Scholar
[18]	苏杰, 刘子超, 廖健颖, 李盈傧, 黄诚 2022 物理学报 71 193201 Google Scholar Su J, Liu Z C, Liao J Y, Li Y B, Huang C 2022 Acta Phys. Sin. 71 193201 Google Scholar
[19]	Herink G, Solli D R, Guide M, Ropers C 2012 Nature 483 190 Google Scholar
[20]	Ortmann L, Landsman A S 2018 Phys. Rev. A 97 023420 Google Scholar
[21]	Husakou A, Im S J, Herrmann J 2011 Phys. Rev. A 83 043839 Google Scholar
[22]	Yavuz I, Bleda E A, Altun Z, Topcu T 2012 Phys. Rev. A 85 013416 Google Scholar
[23]	Ciappina M F, Biegert J, Quidant R, Lewenstein M 2012 Phys. Rev. A 85 033828 Google Scholar
[24]	Ortmann L, Pérez-Hernández J A, Ciappina M F, Schötz J, Chacón A, Zeraouli G, Kling M F, Roso L, Lewenstein M, Landsman A S 2017 Phys. Rev. Lett. 119 053204 Google Scholar
[25]	Ciappina M F, Pérez-Hernández J, Shaaran T, Roso L, Lewenstein M 2014 Phys. Rev. A 89 013409 Google Scholar
[26]	Chacón A, Ortmann L, Cucchietti F, Suárez N, Pérez-Hernández J A, Ciappina M F, Landsman A S, Lewenstein M 2017 Appl. Phys. B 123 116 Google Scholar
[27]	Xu J K, Li Y B, Zhou Y B, Chen Y M, Li M, Yu B H, Lu P X 2022 Opt. Express 30 15951 Google Scholar
[28]	Liang J T, Zhou Y M, Liao Y J, Jiang W C, Li M, Lu P X 2022 Ultrafast Sci. 2022 9842716 Google Scholar
[29]	Ho P J, Eberly J H 2005 Phys. Rev. Lett. 95 193002 Google Scholar
[30]	Haan S L, Breen L, Karim A, Eberly J H 2006 Phys. Rev. Lett. 97 103008 Google Scholar
[31]	Li Y B, Xu J K, Yu B H, Wang X 2020 Opt. Express 28 7341 Google Scholar
[32]	Chen Z H, Su J, Zeng X, Huang X F, Li Y B, Huang C 2021 Opt. Express 29 29576 Google Scholar
[33]	Li Y B, Yu B H, Tang Q B, Wang X, Hua D Y, Tong A H, Jiang C H, Ge G X, Li Y C, Wan J G 2016 Opt. Express 24 6469 Google Scholar
[34]	Li Y B, Qin L L, Liu A H, Zhang K, Tang Q B, Zhai C Y, Xu J K, Chen S, Yu B H, Chen J 2022 Chin. Phys. Lett. 39 093201 Google Scholar
[35]	李盈傧, 秦玲玲, 陈红梅, 李怡涵, 何锦锦, 史璐珂, 翟春洋, 汤清彬, 刘爱华, 余本海 2022 物理学报 71 043201 Google Scholar Li Y B, Qin L L, Chen H M, Li Y H, He J J, Shi L K, Zhai C, Y Tang Q B, Liu A H, Yu B H 2022 Acta Phys. Sin. 71 043201 Google Scholar
[36]	Ciappina M F, Pérez-Hernández J, Shaaran T, Biegert J, Quidant R, Lewenstein M 2012 Phys. Rev. A 86 023413 Google Scholar
[37]	Wang Z, Lan P F, Luo J H, He L X, Lu P X 2013 Phys. Rev. A 88 063838 Google Scholar
[38]	Chacón A, Ciappina M F, Lewenstein M 2016 Phys. Rev. A 94 043407 Google Scholar
[39]	Gao X Z, Landsman A S, Wang H, S Huang P, Zhang Y P, Wang B, Wang Y S, Cao H B, Fu Y X, Pi L W 2021 New J. Phys. 23 113017 Google Scholar
[40]	Yavuz I, Ciappina M F, Chacón A, Altun Z, Kling M F, Lewenstein M 2016 Phys. Rev. A 93 033404 Google Scholar
[41]	Chen Y K, Zhou Y M, Tan J, Li M, Cao W, Lu P X 2021 Phys. Rev. A 104 043107 Google Scholar
[42]	Li Y, Wang X, Yu B, Tang Q, Wang G, Wan J 2016 Sci. Rep. 6 37413 Google Scholar

Cited By

期刊类型引用(2)

1.	赵宁，江英华，周贤韬. 基于单光子的高效量子安全直接通信方案. 物理学报. 2022(15): 36-41 . 百度学术
2.	张晨涛，石小涛，朱文新，朱金龙，郝向英，金锐博. 利用域排列算法设计铌酸锂晶体实现3μm中红外波段频域纯态单光子源. 物理学报. 2022(20): 81-87 . 百度学术

其他类型引用(1)

图 1 对于不同波长, 空间非均匀激光场 (蓝色方块) 和空间均匀激光场 (红色圆圈) 下DI产率随激光强度的变化　(a) 800 nm; (b) 1600 nm; (c) 3200 nm; (d) 4800 nm

Figure 1. Probabilities of DI as a function of laser intensity for the spatially inhomogeneous (blue squares) and spatially homogeneous (red circles) laser fields at different wavelengths: (a) 800 nm; (b) 1600 nm; (c) 3200 nm; (d) 4800 nm.

DownLoad: Full-Size Img PowerPoint

图 2 空间非均匀激光场和空间均匀激光场中返回电子在碰撞前 0.03T 的动能与再碰撞时间的关系 (激光强度均为4×10¹³ W/cm²)　(a)—(d) ε = 0.003; (e)—(h) ε = 0

Figure 2. Kinetic energy of the returning electron before the recollision 0.03T vs. the recollision time for the spatially inhomogeneous and spatially homogeneous laser fields (Laser intensities are both 4×10¹³ W/cm²): (a)–(d) ε = 0.003; (e)–(h) ε = 0.

DownLoad: Full-Size Img PowerPoint

图 3 关联电子对在激光脉冲结束时的发射角分布　(a)—(d) 空间非均匀激光场; (e)—(h) 空间均匀激光场

Figure 3. Distribution of the emitting angle between the correlated electron pairs at the end of the laser pulses: (a)–(d) Spatially inhomogeneous laser fields; (e)–(h) spatially homogeneous laser fields.

DownLoad: Full-Size Img PowerPoint

图 4 空间非均匀激光场和空间均匀激光场下末态电子的横向动量和纵向动量分布 (激光波长均为4800 nm)　(a) ε = 0.003; (b) ε = 0

Figure 4. Transverse and longitudinal momentum distributions of final state electrons for the spatially inhomogeneous and spatially homogeneous laser fields (Laser wavelengths are both 4800 nm): (a) ε = 0.003; (b) ε = 0.

DownLoad: Full-Size Img PowerPoint

图 5 空间非均匀激光场和空间均匀激光场下飞行时间 (a), (c)及时间延迟 (b), (d)的统计分布　(a), (b) ε = 0.003; (c), (d) ε = 0 (图(b), (d)中竖直实线对应0.25个光周期)

Figure 5. Distributions of travel time (a), (c) and time delay (b), (d) for the spatially inhomogeneous and spatially homogeneous laser fields: (a), (b) ε = 0.003; (c), (d) ε = 0 (Vertical solid line in panels (b) and (d) corresponds to 0.25 cycles).

DownLoad: Full-Size Img PowerPoint

图 6 电子沿z轴正方向 (红线) 和负方向 (蓝线) 电离的演化轨迹. 1st, 2nd和3rd分别代表电子第1次、第2次和第3次返回母离子核, 黑色虚线代表原点处振荡的激光电场　(a), (b) 对应的激光波长分别为 800 nm和3200 nm, 激光强度均为4×10¹³ W/cm²

Figure 6. Evolution trajectories of electron ionization along positive z direction (red lines) and the negative z direction (blue lines). 1st, 2nd and 3rd represent the first, second and third return of the electrons to the parent ion core. The black dashed lines show the oscillating laser electric field at the origin: (a), (b) Correspond to wavelengths of 800 nm and 3200 nm, respectively and the laser intensities are both 4×10¹³ W/cm².

DownLoad: Full-Size Img PowerPoint

图 7 (a), (b)两个电子与母离子核的距离随时间的演化; (c), (d)两个电子沿z轴的速度随时间的演化; (a), (c) 沿z轴负方向电离并发生一次返回的电子运动轨迹; (b), (d) 沿z轴负方向电离发生二次返回的电子运动轨迹

Figure 7. (a), (b) Distance of two electrons from the parent ion core as a function of time; (c), (d) velocity of two electrons along the z direction as a function of time; (a), (c) trajectory of electron ionization along negative z direction with one return; (b), (d) trajectory of electron ionization along negative z direction with secondary return.

DownLoad: Full-Size Img PowerPoint

[1]	Lein M, Hay N, Velotta R, Marangos J P, Knight P L 2002 Phys. Rev. Lett. 88 183903 Google Scholar
[2]	Ghimire S, Reis D A 2019 Nat. Phys. 15 10
[3]	Yang B, Schafer K J, Walker B, Kulander K C, Agostini P, DiMauro L F 1993 Phys. Rev. Lett. 71 3770 Google Scholar
[4]	Milošević D B, Paulus G G, Becker W 2003 Opt. Express 11 1418 Google Scholar
[5]	Busuladi M, Gazibegovi-Busuladi A, Miloevi D B, Becker W 2008 Phys. Rev. Lett. 100 203003 Google Scholar
[6]	He M R, Li Y, Zhou Y M, Li M, Cao W, Lu P X 2018 Phys. Rev. Lett. 120 133204 Google Scholar
[7]	Zhou Y, Tolstikhin O I, Morishita T 2016 Phys. Rev. Lett. 116 173001 Google Scholar
[8]	Rudenko A, Jesus V, Ergler T, Zrost K, Ullrich J 2007 Phys. Rev. Lett. 99 263003 Google Scholar
[9]	Quan W, Hao X L, Wang Y L, Chen Y J, Yu S G, Xu S P, Xiao Z L, Sun R P, Lai X Y, Hu S L, Liu M Q, Shu Z, Wang X D, Li W D, Becker W, Liu X J, Chen J 2017 Phys. Rev. A 96 032511 Google Scholar
[10]	Huang C, Zhou Y M, Zhang Q B, Lu P X 2013 Opt. Express 21 11382 Google Scholar
[11]	Tong A H, Liao Q, Zhou Y M, Lu P X 2010 Opt. Express 18 9064 Google Scholar
[12]	Shaaran T, Augstein B B, Figueira D 2011 Phys. Rev. A 84 013429 Google Scholar
[13]	Wang Y L, Xu S P, Chen Y J, Kang H P, Lai X Y, Quan W, Liu X J, Hao X L, Li W D, Hu S L, Chen J, Becker W, Chu W, Yao J P, Zeng B, Cheng Y, Xu Z Z 2017 Phys. Rev. A 95 063415 Google Scholar
[14]	Corkum P B 1993 Phys. Rev. Lett. 71 1994 Google Scholar
[15]	Feuerstein B, Moshammer R, Fischer D, Dorn A, Schroter C D, Deipenwisch J, Crespo Lopez-Urrutia J R, Hohr C, Neumayer P, Ullrich J, Rottke H, Trump C, Wittmann M, Korn G, Sandner W 2001 Phys. Rev. Lett. 87 043003 Google Scholar
[16]	Wang Y L, Xu S P, Quan W, Gong C, Lai X Y, Hu S L, Liu M Q, Chen J, Liu X J 2016 Phys. Rev. A 94 053412 Google Scholar
[17]	Huang C, Zhong M M, Wu Z M 2016 Opt. Express 24 28361 Google Scholar
[18]	苏杰, 刘子超, 廖健颖, 李盈傧, 黄诚 2022 物理学报 71 193201 Google Scholar Su J, Liu Z C, Liao J Y, Li Y B, Huang C 2022 Acta Phys. Sin. 71 193201 Google Scholar
[19]	Herink G, Solli D R, Guide M, Ropers C 2012 Nature 483 190 Google Scholar
[20]	Ortmann L, Landsman A S 2018 Phys. Rev. A 97 023420 Google Scholar
[21]	Husakou A, Im S J, Herrmann J 2011 Phys. Rev. A 83 043839 Google Scholar
[22]	Yavuz I, Bleda E A, Altun Z, Topcu T 2012 Phys. Rev. A 85 013416 Google Scholar
[23]	Ciappina M F, Biegert J, Quidant R, Lewenstein M 2012 Phys. Rev. A 85 033828 Google Scholar
[24]	Ortmann L, Pérez-Hernández J A, Ciappina M F, Schötz J, Chacón A, Zeraouli G, Kling M F, Roso L, Lewenstein M, Landsman A S 2017 Phys. Rev. Lett. 119 053204 Google Scholar
[25]	Ciappina M F, Pérez-Hernández J, Shaaran T, Roso L, Lewenstein M 2014 Phys. Rev. A 89 013409 Google Scholar
[26]	Chacón A, Ortmann L, Cucchietti F, Suárez N, Pérez-Hernández J A, Ciappina M F, Landsman A S, Lewenstein M 2017 Appl. Phys. B 123 116 Google Scholar
[27]	Xu J K, Li Y B, Zhou Y B, Chen Y M, Li M, Yu B H, Lu P X 2022 Opt. Express 30 15951 Google Scholar
[28]	Liang J T, Zhou Y M, Liao Y J, Jiang W C, Li M, Lu P X 2022 Ultrafast Sci. 2022 9842716 Google Scholar
[29]	Ho P J, Eberly J H 2005 Phys. Rev. Lett. 95 193002 Google Scholar
[30]	Haan S L, Breen L, Karim A, Eberly J H 2006 Phys. Rev. Lett. 97 103008 Google Scholar
[31]	Li Y B, Xu J K, Yu B H, Wang X 2020 Opt. Express 28 7341 Google Scholar
[32]	Chen Z H, Su J, Zeng X, Huang X F, Li Y B, Huang C 2021 Opt. Express 29 29576 Google Scholar
[33]	Li Y B, Yu B H, Tang Q B, Wang X, Hua D Y, Tong A H, Jiang C H, Ge G X, Li Y C, Wan J G 2016 Opt. Express 24 6469 Google Scholar
[34]	Li Y B, Qin L L, Liu A H, Zhang K, Tang Q B, Zhai C Y, Xu J K, Chen S, Yu B H, Chen J 2022 Chin. Phys. Lett. 39 093201 Google Scholar
[35]	李盈傧, 秦玲玲, 陈红梅, 李怡涵, 何锦锦, 史璐珂, 翟春洋, 汤清彬, 刘爱华, 余本海 2022 物理学报 71 043201 Google Scholar Li Y B, Qin L L, Chen H M, Li Y H, He J J, Shi L K, Zhai C, Y Tang Q B, Liu A H, Yu B H 2022 Acta Phys. Sin. 71 043201 Google Scholar
[36]	Ciappina M F, Pérez-Hernández J, Shaaran T, Biegert J, Quidant R, Lewenstein M 2012 Phys. Rev. A 86 023413 Google Scholar
[37]	Wang Z, Lan P F, Luo J H, He L X, Lu P X 2013 Phys. Rev. A 88 063838 Google Scholar
[38]	Chacón A, Ciappina M F, Lewenstein M 2016 Phys. Rev. A 94 043407 Google Scholar
[39]	Gao X Z, Landsman A S, Wang H, S Huang P, Zhang Y P, Wang B, Wang Y S, Cao H B, Fu Y X, Pi L W 2021 New J. Phys. 23 113017 Google Scholar
[40]	Yavuz I, Ciappina M F, Chacón A, Altun Z, Kling M F, Lewenstein M 2016 Phys. Rev. A 93 033404 Google Scholar
[41]	Chen Y K, Zhou Y M, Tan J, Li M, Cao W, Lu P X 2021 Phys. Rev. A 104 043107 Google Scholar
[42]	Li Y, Wang X, Yu B, Tang Q, Wang G, Wan J 2016 Sci. Rep. 6 37413 Google Scholar

[1]	Ge Zhen-Jie, Su Xu, Bai Li-Hua. Nonsequential double ionization of Ar atoms in counter-rotating two-color elliptically polarized laser fields. Acta Physica Sinica, 2024, 73(9): 093201. doi: 10.7498/aps.73.20231583
[2]	He Tong-Tong, Liu Zi-Chao, Li Ying-Bin, Huang Cheng. Manipulating nonsequential double ionization of atoms by parallel polarized three-color laser fields. Acta Physica Sinica, 2024, 73(16): 163201. doi: 10.7498/aps.73.20240737
[3]	Liu Yi-Jun, Chen Yi-Wei, Zhu Yu-Jian, Huang Yan, An Dong-Dong, Li Qing-Xin, Gan Qi-Kang, Zhu Wang, Song Jun-Wei, Wang Kai-Yuan, Wei Ling-Nan, Zong Qi-Jun, Liu Shuo-Han, Li Shi-Wei, Liu Zhi, Zhang Qi, Xu Ying-Hai, Cao Xin-Yu, Yang Ao, Wang Hao-Lin, Yang Bing, Andy Shen, Yu Ge-Liang, Wang Lei. Isospin polarized Chern insulator state of C = 4 in twisted double bilayer graphene. Acta Physica Sinica, 2023, 72(14): 147303. doi: 10.7498/aps.72.20230497
[4]	Zhong Guo-Hua, Lin Hai-Qing. Aromatic superconductors: Electron-phonon coupling and electronic correlations. Acta Physica Sinica, 2023, 72(23): 237403. doi: 10.7498/aps.72.20231751
[5]	Su Jie, Liu Zi-Chao, Liao Jian-Ying, Li Ying-Bin, Huang Cheng. Intensity-dependent electron correlation in nonsequential double ionization of Ar atoms in counter-rotating two-color elliptically polarized laser fields. Acta Physica Sinica, 2022, 71(19): 193201. doi: 10.7498/aps.71.20221044
[6]	Zeng Xue, Su Jie, Huang Xue-Fei, Pang Hui-Ling, Huang Cheng. Frequency-ratio-dependent ultrafast dynamics in nonsequential double ionization by co-rotating two-color circularly polarized laser fields. Acta Physica Sinica, 2021, 70(24): 243201. doi: 10.7498/aps.70.20211112
[7]	Huang Cheng, Zhong Ming-Min, Wu Zheng-Mao. Intensity-dependent recollision dynamics in strong-field nonsequential double ionization. Acta Physica Sinica, 2019, 68(3): 033201. doi: 10.7498/aps.68.20181811
[8]	Lin Tong, Hu Die, Shi Li-Yu, Zhang Si-Jie, Liu Yan-Qi, Lv Jia-Lin, Dong Tao, Zhao Jun, Wang Nan-Lin. Infrared spectroscopy study of ironbased superconductor Li0.8Fe0.2 ODFeSe. Acta Physica Sinica, 2018, 67(20): 207102. doi: 10.7498/aps.67.20181401
[9]	Zhang Bin, Zhao Jian, Zhao Zeng-Xiu. Multiconfiguration time-dependent Hartree-Fock treatment of electron correlation in strong-field ionization of H2 molecules. Acta Physica Sinica, 2018, 67(10): 103301. doi: 10.7498/aps.67.20172701
[10]	Huang Cheng, Zhong Ming-Min, Wu Zheng-Mao. Nonsequential double ionization of aligned molecules by few-cycle laser pulses at low intensity. Acta Physica Sinica, 2016, 65(8): 083301. doi: 10.7498/aps.65.083301
[11]	Wu Shao-Quan, Fang Dong-Kai, Zhao Guo-Ping. Effect of electronic correlations on magnetotransport through a parallel double quantum dot. Acta Physica Sinica, 2015, 64(10): 107201. doi: 10.7498/aps.64.107201
[12]	Xin Guo-Guo, Zhao Qing, Liu Jie. Maximum correlation at the transition to the saturation regime of nonsequential double ionization. Acta Physica Sinica, 2012, 61(13): 133201. doi: 10.7498/aps.61.133201
[13]	Tong Ai-Hong, Feng Guo-Qiang, Deng Yong-Ju. Dependence of nonsequential double ionization of He on intensity ratio of orthogonal two-color field. Acta Physica Sinica, 2012, 61(9): 093303. doi: 10.7498/aps.61.093303
[14]	Yu Ben-Hai, Li Ying-Bin. Laser intensity dependence of nonsequential double ionization of argon atoms by elliptically polarized laser pulses. Acta Physica Sinica, 2012, 61(23): 233202. doi: 10.7498/aps.61.233202
[15]	Yu Ben-Hai, Li Ying-Bin, Tang Qing-Bin. The nonsequential double ionization of argon atoms with elliptically polarized laser pulse. Acta Physica Sinica, 2012, 61(20): 203201. doi: 10.7498/aps.61.203201
[16]	Tong Ai-Hong, Liao Qing, Zhou Yue-Ming, Lu Pei-Xiang. Internuclear-distance dependence of nonsequential double ionization of H2 in different alignments. Acta Physica Sinica, 2011, 60(4): 043301. doi: 10.7498/aps.60.043301
[17]	Zhang Dong-Ling, Tang Qing-Bin, Yu Ben-Hai, Chen Dong. Nonsequential double ionization of argon atom below the recollision threshold. Acta Physica Sinica, 2011, 60(5): 053205. doi: 10.7498/aps.60.053205
[18]	Tang Qing-Bin, Zhang Dong-Ling, Yu Ben-Hai, Chen Dong. Three-dimensional classical micro-canonical simulation of nonsequential double ionization with a few-cycle laser pulse. Acta Physica Sinica, 2010, 59(11): 7775-7781. doi: 10.7498/aps.59.7775
[19]	Wang Wei, Sun Jia-Fa, Liu Mei, Liu Su. First-principles calculations on the electronic band structure of β-Pyrochlore superconductors AOs2O6 (A=K,Rb,Cs). Acta Physica Sinica, 2009, 58(8): 5632-5639. doi: 10.7498/aps.58.5632
[20]	Li Hong-Yun, Wang Bing-Bing, Jiang Hong-Bing, Chen Jing, Li Xiao-Feng, Liu Jie, Gong Qi-Huang, Fu Pan-Ming. The effect of static electric field on nonsequential double ionization. Acta Physica Sinica, 2008, 57(1): 124-131. doi: 10.7498/aps.57.124

期刊类型引用(2)

1.	赵宁，江英华，周贤韬. 基于单光子的高效量子安全直接通信方案. 物理学报. 2022(15): 36-41 . 百度学术
2.	张晨涛，石小涛，朱文新，朱金龙，郝向英，金锐博. 利用域排列算法设计铌酸锂晶体实现3μm中红外波段频域纯态单光子源. 物理学报. 2022(20): 81-87 . 百度学术

姓名
邮箱
手机号码
标题
留言内容
验证码

Search

Article

留言板