啁啾激光場中氫分子離子電離過程的理論研究

張彩萍，苗向陽

（山西師范大學(xué)物理與信息工程學(xué)院，臨汾041004）

1 引 言

近年，基于原子在強激光場中動力學(xué)研究的不斷成熟，及激光技術(shù)的不斷提高，越來越多的學(xué)者將目光投向分子在超快強激光中的動力學(xué)行為［1-3］.其中氫分子離子作為最簡單的研究體系已被廣泛應(yīng)用到分子在強激光場中電離機制的研究中，并取得豐碩的成果.強激光場中的電離理論主要包含電荷共振增強電離（Charge resonance enhanced ionization）［3］，多光子電離（Multiphoton ionization）［4］， 遂 穿 電 離（Tunneling ionization）［4］， 閾 上 電 離（Above threshold ioniza-tion）［5］，躍壘電離（Over barrier ionization）.Bandrauk［6］等人利用激光泵浦技術(shù)在非B-O（Born-Oppenheimer）近似下研究處于疊加態(tài)的氫分子離子的光電離機制.Ben-Itzhak［7］等人在B-O 近似下研究超短激光脈沖作用下氫分子離子的電荷共振增強電離機制并與Zuo［3］等人的理論解釋做了比較.而啁啾場中有關(guān)電離過程的研究較少.本文通過使用超短啁啾脈沖（5fs，790nm）進一步探究氫分子離子在強激光場中的電離動力學(xué).

2 理論方法

本文在計算中通過執(zhí)行LZH-DICP 程序［8］，在非B-O 近似下數(shù)值求解了氫分子離子一維含時薛定諤方程：其中體系的哈密頓量可表示為

其中高斯外包絡(luò)函數(shù)為

程序中初始波函數(shù)通過對角化在分離格點表象（DVR）［8，9］中構(gòu)建的矩陣中得到，波函數(shù)隨時間的演化采用二階劈裂算符［8-11］的方法：

T 為核與電子的動能算符，V 表示體系所有的勢能.為了節(jié)約時間，將核與電子分開傳播得

3 結(jié)果和討論

圖1描述了在不同激光強度（4.0×1014W／cm2-1.0×1015W／cm2）下氫分子離子的電離幾率隨啁啾參數(shù)的變化，圖中每一點為響應(yīng)啁啾參數(shù)下最大的電離幾率C（t）

其中jz為流算符［12，13］

圖1 氫分子離子電離幾率在不同激光強度下隨啁啾參數(shù)β的變化.四條虛線代表的激光強度分別為4.0×1014 W／cm2（三角形），6.0×1014 W／cm2（正方形），8.0×1014 W／cm2（鉆石形）和1.0×1015 W／cm2（圓形）Fig.1 The ionization probabilities of H+2 with the change of chirp parameterβat different laser intensities.And the four dotted lines represent different laser intensities：4.0×1014 W／cm2（triangles），6.0×1014 W／cm2（squares），8.0×1014 W／cm2（diamonds）and 1.0×1015 W／cm2（circles）

圖2 （a）和（b）分別表示β=0和β=6時激光場的變化.（c）和（d）分別表示平均核間距在相應(yīng)啁啾參數(shù)下隨時間的變化.其中實線代表激光強度為1.0×1015 W／cm2，點劃線代表6.0×1014 W／cm2Fig.2 （a）and（b）show the electronic laser fields atβ=0andβ=6，respectively.（c）and（d）present the corresponding nuclear motion R（t）with the laser intensities of 1.0×1015 W／cm2（solid line）and 6.0×1014 W／cm2（dot-dash line），respectively

本文首先從激光作用后核運動的角度分析啁啾脈沖對氫分子離子電離過程的影響.圖2（a）和（b）分別描述β=0和β=6時激光場隨時間的變化，隨著啁啾參數(shù)去值得增大，B 區(qū)域的次峰值個數(shù)不斷增加即隨著啁啾參數(shù)取值的增大激光頻率不斷增加.而圖1中電離幾率從β=0到β=6不斷增加主要歸因于B 區(qū)域激光頻率的不斷增加.根據(jù)公式

（c）和（d）圖分別給出了相應(yīng)啁啾參數(shù)下I=1.0×1015W／cm2（實線）和I=6.0×1014W／cm2（點劃線）平均核間距R（t）隨時間的變化.從（c）圖中可以看出在激光作用前氫分子離子的核間距保持在2.8a.u.，當t=10fs時，氫分子離子開始受到激光的激發(fā)，兩核運動加快核間距不斷增大，而且在激光場的上升區(qū)域核間距的增大較快，當t＞18fs即在激光場的下降區(qū)域核間距的增加相對較慢.如圖所示當激光作用結(jié)束后，核間距在較高激光強度下（I=1.0×1015W／cm2）由2.8a.u.增加到7.3a.u.，較低強度下（I=6.0×1014W／cm2）由2.8a.u.增加到5.7a.u..核間距的越大導(dǎo)致核對電子的束縛越小，因而I=1.0×1015W／cm2下氫分子離子的電離幾率（0.55）要遠大于I=6.0×1014W／cm2的電離幾率（0.2）.（d）圖為引入啁啾脈沖后β=6時，平均核間距隨時間變化，對比（c）與（d）可發(fā)現(xiàn)引入啁啾脈沖后，核波包的擴散更顯著，核間距迅速增大.兩種激光強度下的核間距分別為12.3a.u.（實線）和11.7a.u.（點劃線），因而如圖1所示在相應(yīng)的激光參數(shù)下電離幾率明顯增強.從而進一步從核運動的角度解釋了啁啾脈沖可增強氫分子離子的電離過程.

圖3 不同激光強度不同啁啾參數(shù)下電子波包密度隨時間的變化.第一列和第二列所對應(yīng)的啁啾參數(shù)分別為β=0和β=6.（a）和（c）所對應(yīng)的激光強度為1.0×1015 W／cm2，（b）和（d）相應(yīng)的激光強度為6.0×1014 W／cm2Fig.3 The time evolution of the electronic probability densities at different intensities withβ=0（first column）andβ=6（second column）.The intensities are 1.0×1015 W／cm2 for（a）and（c），and 6.0×1014 W／cm2 for（b）and（d）

為了進一步從電子運動的角度分析氫分子離子的電離過程，根據(jù)公式

圖3描述了不同激光強度下電子波包密度隨時間的變化.第一列與第二列所對應(yīng)的啁啾參數(shù)為β=0與β=6，且（a）和（c）相應(yīng)激光強度為I=1.0×1015W／cm2，（b）和（d）相應(yīng)的激光強度為I=6.0×1014W／cm2.如圖所示，隨著激光場周期性的改變，電子在電場力的作用下同樣做周期性的運動.在無啁啾激光場中，從（a）圖中可以看出在t=18fs，即達到激光的峰值強度1.0×1015W／cm2時，氫分子離子受到強激光的作用，大部分電子被瞬間電離，電子波包密度隨之減少.與（a）圖比較可發(fā)現(xiàn)，當峰值強度為6.0×1014W／cm2時，在18fs到20fs這一時間段內(nèi)仍由相當數(shù)量的電子存在，說明較高的激光強度可有效使電子脫離原子核的束縛，提高電離幾率，與圖1中β=0時的數(shù)值結(jié)果相吻合；同理，當氫分子離子處于啁啾激光場中時，隨著啁啾參數(shù)取值的增大，電子在激光場中的運動也發(fā)生明顯的改變，如圖2（b）中，在B區(qū)域內(nèi)基本不存在電子的分布，即受到啁啾激光的激發(fā)后，氫分子離子吸收高頻光子達到電離態(tài)，電子被瞬間電離.比較（a）（或（c））與（b）（或（d））可以發(fā)現(xiàn)啁啾激光場中的電離過程發(fā)生的更快些，在適當?shù)倪眳?shù)下，相應(yīng)的幾率也會有明顯的增加.

4 結(jié) 論

本文在非B-O 近似下，通過數(shù)值求解氫分子離子體系在強激光條件下一維含時薛定諤方程，在超短啁啾激光場中進一步研究氫分子離子的電離過程.研究中有效利用流算符計算得到不同激光強度（4.0×1014W／cm2～1.0×1015W／cm2）下的電離幾率，結(jié)果表明啁啾脈沖的使用可有效提高氫分子離子的電離幾率，且隨著激光強度的增大，電離過程對啁啾參數(shù)的依賴性增強.通過計算平均核間距及電子波包密度，從核與電子在激光場中運動的角度，進一步說明啁啾激光強度越大，核運動越顯著核間距越大，對電子的束縛越小，電子的擴散越明顯，從而有效增強氫分子離子的電離過程.通過本文的理論研究，以期人們能夠進一步認識強激光場中小分子體系的電離動力學(xué).

致 謝 感謝韓克利研究員提供的LZH-DICP計算程序以及寶貴的指導(dǎo)意見.

［1］ Chelkowski S，Bandrauk A D，Staudte A，et al.Dynamic nuclear interference structures in the coulomb explosion spectra of a hydrogen molecule in intense laser fields：reexamination of molecular enhanced ionization［J］.Phys.Rev.A，2007，76：013405.

［2］ Uhlmann M，Kunert T，Schmidt R.Molecular alignment of fragmenting H+2and H2in strong laser fields［J］.Phys.Rev.A，2005，72：045402.

［3］ Zuo T，Bandrauk A D.Charge-resonance-enhanced ionization of diatomic molecular ions by intense lasers［J］.Phys.Rev.A，1995，52：R2511.

［4］ Posthumus J H.The dynamics of small molecules in intense laser fields［J］.Rep.Prog.Phys.，2004，67：623.

［5］ Paulus G G，Lindner F，Walther H，et al.Measurement of the phase of few-cycle laser pulses ［J］.Phys.Rev.Lett.，2003，91：253004.

［6］ Bandrauk A D，Chelkowski S，Corkum P B，et al.Attosecond photoionization of a coherent superposition of bound and dissociative molecular states：effect of nuclear motion ［J］.J.Phys.B，2009，42：134001.

［7］ Ben-Itzhak I，Wang P Q，Sayler A M，et al.Elusive enhanced ionization structure for H+2in intense ultrashort laser pulses［J］.Phys.Rev.A，2008，78：063419.

［8］ Lu R F，Zhang P Y，Han K L.Attosecond-resolution quantum dynamics calculations for atoms and molecules in strong laser fields［J］.Phys.Rev.E，2008，77：066701.

［9］ Zhang P Y，Lu R F，Zhang A J，et al.Full six-dimensional nonadiabatic quantum dynamics calculation for the energy pooling reaction O2（aΔ1Δ）+O2（aΔ1Δ）→O2（bΔ1∑）+O2（XΔ3∑）［J］.J.Chem.Phys.，2008，128：091103.

［10］ Han J，JI A P，Miao X Y.An isolated 40attosecond pulse generated by using the combination of a two-color laser pulse and a low-frequency field［J］.J.A.Mol.Phys.，2012，29（4）：650（in Chinese）［韓靜，吉安平，苗向陽.通過雙色激光場和低頻電場合成獲得孤立40阿秒脈沖［J］.原子與分子物理學(xué)報，2012，29（4）：650］

［11］ Liu S S，Miao X Y.Enhancement of high-order harmonic emission by using a coherent superposition in a two-color laser field［J］.J.A.Mol.Phys.，2012，29（5）：881（in Chinese）［劉莎莎，苗向陽.在雙色場中使用相干疊加態(tài)來提高高次諧波的發(fā)射［J］.原子與分子物理學(xué)報，2012，29（5）：881］

［12］ Feuerstein B，Thumm U.On the computation of momentum distributions within wavepacket propagation calculations ［J］.J.Phys.B，2003，36：707.

［13］ Feuerstein B，Thumm U.Fragmentation of H+2in strong 800nm laser pulses：initial-vibrational-state dependence［J］.Phys.Rev.A，2003，67：043405.

［14］ Kamta G L，Bandrauk A D.Phase dependence of enhanced ionization in asymmetric molecules［J］.Phys.Rev.Lett.，2005，94：203003.