激光与光电子学进展, 2023, 60 (1): 0102002, 网络出版: 2022-12-05   

通过奇偶谐波谱重构不对称平面分子的构型 下载: 705次

Probing the Structure of Asymmetric Planar Molecules Using Odd-Even High Harmonics
作者单位
延安大学物理与电子信息学院,陕西 延安 716000
摘要
通过求解含时薛定谔方程的数值方法,研究了不同分子构型的平面分子H32+的高次谐波产生,结果表明:奇偶谐波同时释放,且奇偶谐波的相对产量对取向角和分子构型都很灵敏。基于上述结果,提出了利用奇偶高次谐波谱探测原子核位置的简单方法,从而可以通过奇偶高次谐波谱重构取向不对称平面分子的键长和键角,同时利用高次谐波谱为阿秒尺度探测多中心分子结构提供帮助。
Abstract
The generation of high-order harmonics from planar molecules H32+ with different molecular configurations is studied using the numerical solution of the time-dependent Schr?dinger equation. The results show that planar molecules H32+ emit both odd and even harmonics, and the relative yields of odd and even harmonics are sensitive to orientation angles and molecular configuration. A simple method is proposed, based on the above results, to detect the position of atomic nucleus via odd and even high-order harmonic spectra. Bond lengths and angles from oriented asymmetric planar molecules can be reconstructed using odd-even harmonics at different angles. Meanwhile, probing the polycentric molecular structure in attoseconds using a high-order harmonic spectrum is helpful.

1 引言

在过去的20年,原子和分子高次谐波产生1-2在强激光与物质相互作用的实验和理论研究中一直是一个热门课题。研究高次谐波的兴趣部分源于其在阿秒科学中的重要应用3-4,包括电子结构和动力学的阿秒探测5-6。可通过三步模型来理解原子分子的高次谐波产生过程7:1)电子通过隧穿电离;2)电子在激光场中传播;3)自由电子返回束缚态并伴随高能光子的释放。

近年来,两中心对称分子的高次谐波产生得到广泛研究,如H2+8-9N210-11。结果表明,一个显著的最小值出现在H2+的高次谐波谱上8,这个最小值得到了广泛的理论和实验研究12-14,并在分子轨道层析成像中显示了其重要意义10。对于HeH2+15-17CO18-19NO20等两中心不对称分子,其对称性破缺导致奇偶谐波产生。研究表明,奇偶谐波谱具有不同谱的性质21-22,并携带不同靶的信息23-24。在高次谐波谱中,奇偶谐波具有重要的应用前景:利用奇偶高次谐波谱与分子结构的关系,重构不对称分子轨道19;在不对称库仑势的作用下,利用奇偶高次谐波谱,探究电子的超快动力学25-26。之前主要研究方向为取向双原子分子的奇偶谐波谱,最近学者们越来越关注激光与多原子分子的相互作用,如平面分子H3+27H32+28D3+29。对于平面三角形分子H32+,奇偶谐波都会释放,由于分子轨道不对称性,平面分子H32+的情况更加复杂,并可能出现一些新效应。本课题组之前已通过奇偶谐波谱探测过多中心线性分子的结构30,为了将奇偶高次谐波谱应用于不同构型的多中心不对称分子,需对此问题进行详细的理论研究。

通过求解含时薛定谔方程的数值方法,详细研究了不同取向角及不同分子构型下三维取向的平面分子H32+的奇偶高次谐波产生。结果表明:奇偶谐波谱对平面取向角、空间取向角及核间距比较灵敏,因此奇偶谐波的相对产量及交点随着平面取向角、空间取向角和分子构型的变化而变化。在实验中,平面分子H32+是通过平面分子H3+释放一个电子得到的,但是其不稳定。在平面分子H32+不稳定的情况下,化学反应中的三角键长在近飞秒时间尺度上是变化的28,模拟表明不同分子构型下的瞬时键长可通过不同取向角下奇偶谐波来重构。

2 数值方法

所研究的三核平面分子H32+的哈密顿函数为H(t)=p22+V(r)+rE(t)pr为电子的动量和位矢。采用的是原子单位=e=me=1。其中平面分子H32+被假设位于xoy平面内,激光场在取向Ω=θ,φ方向是线极化的。在三维模拟中,平面分子H32+坐标系统缩略示意如图1所示。分子系统的质心与坐标系的原点o重合,θ为激光极化与z轴的夹角,而φ为激光极化的方位角,R1R2R3为核间距。采用的软心库仑势为V(r)=-Zξ+r-Ra2-Zξ+r-Rb2-Zξ+r-Rc2ξ=0.5为用来避免库仑奇点的软心参数,Z为有效电荷,通过调整Z获得期望的基态能量,RaRbRc为原点与3个核之间的矢。E(t)=ef(t)Esin ω0t为外电场,f(t)为包络函数,eeθ,φ=sin θcos φ,sin θsin φ,cos φ为沿着激光极化方向的单位矢,Eω0分别为外电场的振幅及频率。

图 1. 三维模拟中平面分子H32+的坐标系统示意图

Fig. 1. Schematic diagram of coordinate system of planar molecules H32+ in 3D simulation

下载图片 查看所有图片

在计算中采用10个周期的激光脉冲,其中3个光学周期是线性上升的,另外7个光学周期保持不变。所采用的激光强度和波长分别为I=8×1014 W/cm2λ=800  nm。含时薛定谔方程iψ˙t=Htψtψ˙(t)为波函数ψ(t)的导数)通过谱方法31数值求解,每个周期的时间格点是2048,对于精确的模拟,采用的空间格点为Lx×Ly×Lz=410×51×51 a.u.,空间步长为Δx=Δy=Δz=0.4 a.u.。每个时间步长后,将Ht)的含时薛定谔方程的ψt乘以一个窗函数,用以吸收边界处的连续态波包。沿着x轴的窗函数有下列形式:当xx0时,F(x)=1;当xx0时,F(x)=cos18πx-x0Lx-2x0,采用x0=Lx8作为沿着x轴方向吸收方案的边界,y轴、z轴相似处理。得到含时薛定谔方程的波ψt,高次谐波谱的干涉项F||ω可表示为

F||(ω)=ψ(t)|e||V|ψ(t)exp(iωt)dt

式中:e||为激光极化方向,ω为释放光子频率;V为对库仑势求梯度。为了便于比较,此处不同分子构型的平面分子H32+具有相同的电离能Ip=1.285 a.u.

3 解析方法

基于强场近似原理32,由于势的对称性,原子和对称分子仅释放奇次谐波,奇次谐波的产生路径可简单表示为〈p|er|0〉〈0|er|pr,此公式意味着电子从相同的初始态0电离出来并重新与之复合,其中初始态0具有确定的宇称,在高次谐波产生过程中,偶极矩阵元〈p|exr|0〉i与电子沿着激光极化方向ex电离有关,而偶极矩阵元〈0|exr|pr与电子沿着e方向复合有关。

不对称分子高次谐波的产生是一个多通道过程,例如本研究中的模型平面分子H32+,基态0没有确定的宇称,结果奇偶谐波谱都会释放,这种情况下基态可以分解为偶宇称和奇宇称,即|=0〉=|0g〉+|0u〉。因此不对称分子奇次谐波的释放由两条路径产生贡献,可以表示为Lggo=〈pu|er|0gi〈0g|er|purLuuo=〈pg|er|0ui〈0u|er|pgr,第一条路径Lggo对应电子从初始态的偶宇称|0g〉电离并与之复合,第二条路径Luuo对应电子从初始态的奇宇称|0u〉电离并与之复合;类似地,当电子从基态的不同宇称电离并返回到基态的不同宇称时,偶次谐波就会释放,通过Luuo=〈pu|er|0gi〈0u|er|pgr Luuo=〈pg|er|0ui〈0g|er|pur表示这两条偶次路径。数值模拟表明,基态的偶宇称|0g〉更容易通过隧穿电离,因此,这条路径LggoLgue在奇(偶)次谐波的释放中起主导作用。此外,在电离过程中,激光场使库仑势发生了严重扭曲,导致电离过程对分子结构不敏感。对于高次谐波的产生,只有返回过程与分子结构紧密相关。对于奇次谐波产生路径LggoDodd=〈0g|er|pur;对于偶次谐波产生路径LgueDeven=〈0u|er|pgr。现在探讨三核平面分子H32+的奇偶高次谐波谱与奇偶束缚态-连续态跃迁偶极矩之间的密切关系。假设平面分子的基态波函数为ϕ0r=r0=Nfexp-κra+exp-κrb+exp-κrcra=r-R1rb=r-R2rc=r-R3κ=2IpNf为归一化系数,Ip为基态0的电离能,故在沿着激光偏振方向,与奇次谐波释放有关的偶极矩项可以写为

Doddω,θ,φGoddω,θ,φdrϕ0rersinpkr

式中:Goddω,θ,φ=cospkR1cos θ+cospkR2cos θ+cospkR3cos θ为干涉项,在沿着激光偏振方向,与偶次谐波释放有关的偶极矩项可以写为

Devenω,θ,φGevenω,θ,φdrϕ0rersinpkr

式中:Gevenω,θ,φ=sinpkR1cos θ+sinpkR2cos θ+sinpkR3cos θ为干涉项。连续态由平面波函数|exp(ipkr)〉近似得到,动量大小pk=2ω,这里考虑了库仑效应10

基于强场近似,沿激光偏振方向相应偶极矩Dodd(even)和高次谐波谱Sodd(even)之间的定性关系1019

Soddω,θ,φ=Nθ,φaωDoddω,θ,φSevenω,θ,φ=Nθ,φaωDevenω,θ,φ

式中:aωθφ的不敏感连续电子的光谱振幅;Nθ,φ为电离率。

4 数值结果和讨论

平面分子H32+在不同的平面和空间取向角下沿着激光偏振方向发射的奇偶高次谐波谱如图2所示,图2中所采用的核间距为R1=1.8 a.u.R2=1.6 a.u.R3=2 a.u.,且每列数据具有相同的空间取向角θ,结果显示在每列的顶部。

图 2. 三维取向平面分子H32+的奇偶谐波谱和相应偶极矩的比较。(a)θ=90oφ=20o的奇偶谐波谱;(b)θ=90oφ=50o的奇偶谐波谱;(c)θ=90oφ=20oφ=50o的奇偶偶极矩;(d)θ=45oφ=20o的奇偶谐波谱;(e)θ=45oφ=50o的奇偶谐波谱;(f)θ=45oφ=20oφ=50o的奇偶偶极矩

Fig. 2. Comparison of odd and even harmonic spectra and corresponding dipoles of 3D orientation planar molecules H32+. (a) θ=90o,φ=20o odd and even harmonic spectra; (b) θ=90o,φ=50o odd and even harmonic spectra; (c) θ=90o, φ=20o and φ=50o odd and even dipoles; (d) θ=45o, φ=20o odd and even harmonic spectra; (e) θ=45o, φ=50o odd and even harmonic spectra; (f) θ=45o, φ=20o and φ=50o odd and even dipoles

下载图片 查看所有图片

图2中:φ为平面取向角;odd为奇次谐波谱;even为偶次谐波谱;20° oddφ=20°的奇次固有偶极矩,20° evenφ=20°的偶次固有偶极矩,50o oddφ=50°的奇次固有偶极矩,50° evenφ=50°的偶次固有偶极矩。对于大角度和中间角θ及所有角φ,通过比较发现不同构型的平面分子H32+的偶次谐波的发射较强(除了由于分子的几何形状φ=60°的等边三角形情况)。总体上偶次谐波随空间取向角θ的减小而减小,当θ=0°时,偶次谐波产量消失,且φ=90°的偶次谐波产量比φ=0°时偶次谐波产量高,奇偶谐波谱的相对产量显示出明显的角依赖。例如图2(a)中,对于低阶谐波谱而言,偶次谐波的产量低于奇次谐波的产量;对于中间和高阶次谐波谱则相反:因此奇偶高次谐波谱有一个明显的交点。随着空间取向角减小及平面取向角增大,交点的位置向更高阶的谐波方向偏移,且随着核间距和激光波长的变化,奇偶谐波的相对产量也会发生变化,导致奇偶谐波谱的交点会向其他谐波阶方向偏移。

谐波谱和相应偶极矩之间进一步的比较表明:在不同角度下,图2中奇偶谐波的相对产量与相关偶极矩的预测一致,见式(4)。例如,在θ=90°的情况下,图2(a)中φ=20°图2(b)中φ=50°的奇偶谐波谱交点与图2(c)中相应偶极矩曲线的交点一致。这些现象也适用于其他空间取向角θ和平面取向角φ及不同分子构型的平面分子。通过以上分析,无论偶次谐波的产量比奇次谐波高,或与奇次谐波产量差不多,奇偶高次谐波谱在相对产量和交点方面都与相应的偶极矩很好地吻合。

可从奇偶高次谐波谱上知道奇偶谐波谱交点的能量ω,根据高次谐波谱Soddeven与相应偶极矩Doddeven之间的定性关系,再通过奇偶偶极矩Doddeven的交点来重构多中心不对称平面分子H32+的键长和键角,即Dodd=Deven。又因为积分项drexp-κre||·rsinpk·r为原子的偶极矩,而奇偶偶极矩的积分项是相等的,因此可通过干涉项Goddeven来重构不对称平面分子H32+的键长和键角,即Godd=Geven。假定奇偶谐波谱在θ1φ1θ2φ2θ3φ3三组角上的交点分别位于3个不同的谐波阶,能量分别为ω1ω2ω3,因此有

Goddω1,θ1,φ1=Gevenω1,θ1,φ1Goddω2,θ2,φ2=Gevenω2,θ2,φ2Goddω3,θ3,φ3=Gevenω3,θ3,φ3

所有干涉项GoddGeven都为R1R2R3的函数,式(5)的解给出了R1R2R3的预测。为检测上述方案的适用性,后续讨论中关注3种平面三角形构型:图3中的等边三角形R1=R2=R3=2 a.u.图4中的等腰三角形R1=R2=2 a.u.R3=1.6 a.u.图5中的任意三角形R1=1.8 a.u.R2=1.6 a.u.R3=2  a.u.ω1ω2ω3分别对应三组角θ1,φ1θ2,φ2θ3,φ3奇偶谐波交点的谐波能量。在模拟中,把空间取向角从θ=90°改变到θ=45°,平面取向角从φ=0°改变到φ=90°,并在不同交点及不同平面和空间角度下,对通过式(5)重构的键长进行平均,相关结果如图3~图5所示。

图 3. 平面等边三角形构型的结果比较

Fig. 3. Comparison of results for planar equilateral triangle configuration

下载图片 查看所有图片

图 4. 平面等腰三角形构型的结果比较

Fig. 4. Comparison of results for planar isosceles triangle configuration

下载图片 查看所有图片

图 5. 平面任意三角形构型的结果比较

Fig. 5. Comparison of results for planar scalene triangle configuration

下载图片 查看所有图片

图3~图5中:accurate result为精确结果;reconstructed result为基于式(5)的重构结果。

图3图4图5中可知,对于所有的平面三角形构型情况,与精确的键长和键角比,基于式(5)的重构方案基本上重构了键长和键角,如:图3中的等边三角形,基于式(5)的重构方案重构了左侧底角A'59.7°和顶角B'58.2°,非常接近精确的键角A=60°B=60°,基于式(5)重构的键长为R1'=2.033 a.u.R2'=2.08 a.u.R3'=2.033 a.u.,与精确的键长R1=R2=R3=2 a.u.也非常接近,这些结果说明该方案是适用的。在其他激光和分子参数下,也验证了通过奇偶高次谐波谱重构原子核位置这个简单方法的正确性。

5 结论

综上,利用奇偶谐波谱的交点以及谐波谱与相应偶极矩的关系,在三维空间取向下重构了平面不对称分子H32+的分子构型。重构的键长和键角与精确的键长和键角符合良好,这些现象也适用于考虑取向效应的情况。由于平面不对称分子在许多化学反应过程中都是活跃的,该结果对探索超快过程中平面不对称分子的结构和动力学具有重要意义。

参考文献

[1] McPherson A, Gibson G, Jara H, et al. Studies of multiphoton production of vacuum-ultraviolet radiation in the rare gases[J]. Journal of the Optical Society of America B, 1987, 4(4): 595-601.

[2] L'Huillier A, Schafer K J, Kulander K C. Theoretical aspects of intense field harmonic generation[J]. Journal of Physics B: Atomic, Molecular and Optical Physics, 1991, 24(15): 3315-3341.

[3] Corkum P B, Krausz F. Attosecond science[J]. Nature Physics, 2007, 3(6): 381-387.

[4] Krausz F, Ivanov M. Attosecond physics[J]. Reviews of Modern Physics, 2009, 81(1): 163-234.

[5] Lépine F, Ivanov M Y, Vrakking M J J. Attosecond molecular dynamics: fact or fiction?[J]. Nature Photonics, 2014, 8(3): 195-204.

[6] Chen J G, Yang Y J, Chen J, et al. Probing dynamic information and spatial structure of Rydberg wave packets by harmonic spectra in a few-cycle laser pulse[J]. Physical Review A, 2015, 91(4): 043403.

[7] Corkum P B. Plasma perspective on strong field multiphoton ionization[J]. Physical Review Letters, 1993, 71(13): 1994-1997.

[8] Lein M, Hay N, Velotta R, et al. Role of the intramolecular phase in high-harmonic generation[J]. Physical Review Letters, 2002, 88(18): 183903.

[9] Li W Y, Dong F L, Yu S J, et al. Ellipticity of near-threshold harmonics from stretched molecules[J]. Optics Express, 2015, 23(24): 31010-31025.

[10] Itatani J, Levesque J, Zeidler D, et al. Tomographic imaging of molecular orbitals[J]. Nature, 2004, 432(7019): 867-871.

[11] Zhou X X, Tong X M, Zhao Z X, et al. Alignment dependence of high-order harmonic generation from N2 and O2 molecules in intense laser fields[J]. Physical Review A, 2005, 72(3): 033412.

[12] Chen Y J, Liu J, Hu B. Reading molecular messages from high-order harmonic spectra at different orientation angles[J]. The Journal of Chemical Physics, 2009, 130(4): 044311.

[13] Chen Y J, Hu B. Role of ionization in orientation dependence of molecular high-order harmonic generation[J]. The Journal of Chemical Physics, 2009, 131(24): 244109.

[14] Li W Y, Wang S, Shi Y Z, et al. Probing the structure of stretched molecular ions with high-harmonic spectroscopy[J]. Journal of Physics B: Atomic, Molecular and Optical Physics, 2017, 50(8): 085003.

[15] Chen Y J, Zhang B. Tracing the structure of asymmetric molecules from high-order harmonic generation[J]. Physical Review A, 2011, 84(5): 053402.

[16] Zhang B, Chen Y J, Jiang X Q, et al. Identifying the interference effect in different harmonic-emission channels from oriented asymmetric molecules[J]. Physical Review A, 2013, 88(5): 053428.

[17] Lu R F, Yu C, Wang Y H, et al. Control of electron localization to isolate and enhance molecular harmonic plateau in asymmetric HeH2+ system[J]. Physics Letters A, 2014, 378(1/2): 90-94.

[18] Pan Y, Zhao S F, Zhou X X. Generation of isolated sub-40-as pulses from gas-phase CO molecules using an intense few-cycle chirped laser and a unipolar pulse[J]. Physical Review A, 2013, 87(3): 035805.

[19] Chen Y J, Fu L B, Liu J. Asymmetric molecular imaging through decoding odd-even high-order harmonics[J]. Physical Review Letters, 2013, 111(7): 073902.

[20] Wu J, Zeng H, Guo C L. Triple-ionization-induced dissociation of NO in strong laser fields[J]. Physical Review A, 2006, 74(3): 031404.

[21] Bian X B, Bandrauk A D. Multichannel molecular high-order harmonic generation from asymmetric diatomic molecules[J]. Physical Review Letters, 2010, 105(9): 093903.

[22] Heslar J, Telnov D, Chu S I. High-order-harmonic generation in homonuclear and heteronuclear diatomic molecules: exploration of multiple orbital contributions[J]. Physical Review A, 2011, 83(4): 043414.

[23] Frumker E, Hebeisen C T, Kajumba N, et al. Oriented rotational wave-packet dynamics studies via high harmonic generation[J]. Physical Review Letters, 2012, 109(11): 113901.

[24] Kraus P M, Rupenyan A, Wörner H J. High-harmonic spectroscopy of oriented OCS molecules: emission of even and odd harmonics[J]. Physical Review Letters, 2012, 109(23): 233903.

[25] Frumker E, Kajumba N, Bertrand J B, et al. Probing polar molecules with high harmonic spectroscopy[J]. Physical Review Letters, 2012, 109(23): 233904.

[26] Yu S J, Zhang B, Li Y P, et al. Ellipticity of odd-even harmonics from oriented asymmetric molecules in strong linearly polarized laser fields[J]. Physical Review A, 2014, 90(5): 053844.

[27] Lötstedt E, Kato T, Yamanouchi K. D3+and H3+in intense laser fields studied with a quasiclassical model[J]. Physical Review A, 2012, 85(5): 053410.

[28] Lefebvre C, Lu H Z, Chelkowski S, et al. Electron-nuclear dynamics of the one-electron nonlinear polyatomic molecule H32+ in ultrashort intense laser pulses[J]. Physical Review A, 2014, 89(2): 023403.

[29] Lötstedt E, Kato T, Yamanouchi K. Classical dynamics of laser-driven D₃⁺[J]. Physical Review Letters, 2011, 106(20): 203001.

[30] Su N, Yu S J, Li W Y, et al. Probing the structure of multi-center molecules with odd-even high harmonics[J]. Chinese Physics B, 2018, 27(5): 054213.

[31] Feit M D, Jr Fleck J A, Steiger A. Solution of the Schrödinger equation by a spectral method[J]. Journal of Computational Physics, 1982, 47(3): 412-433.

[32] Lewenstein M, PhBalcou, CorkumP. B., et al. Theory of high-harmonic generation by low-frequency laser fields[J]. Physical Review A, 1994, 49(3): 2117-2132.

于术娟, 刘竹琴, 刘艳峰, 李雁鹏. 通过奇偶谐波谱重构不对称平面分子的构型[J]. 激光与光电子学进展, 2023, 60(1): 0102002. Shujuan Yu, Zhuqin Liu, Yanfeng Liu, Yanpeng Li. Probing the Structure of Asymmetric Planar Molecules Using Odd-Even High Harmonics[J]. Laser & Optoelectronics Progress, 2023, 60(1): 0102002.

本文已被 1 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!