Influence of multi-photon excitation on the atomic above-threshold ionization*
Tian Yuan-Yea),b), Wang Chun-Chenga),b), Li Su-Yua),b), Guo Fu-Minga),b), Ding Da-Juna),b), Wim-G Roeterdinka),b), Chen Ji-Genc), Zeng Si-Liangd),e), Liu Xue-Shena),b), Yang Yu-Juna),b)†
Institute of Atomic and Molecular Physics, Jilin University, Changchun 130012, China
Jilin Provincial Key Laboratory of Applied Atomic and Molecular Spectroscopy, Jilin University, Changchun 130012, China
Department of Physics and Materials Engineering, Taizhou University, Taizhou 318000, China
Data Center for High Energy Density Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
Science and Technology Computational Physics Laboratory, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China

Corresponding author. E-mail: yangyj@jlu.edu.cn

*Project supported by the National Basic Research Program of China (Grant No. 2013CB922200), the National Natural Science Foundation of China (Grants Nos. 11274141, 11034003, 11304116, 11274001, and 11247024), and the Jilin Provincial Research Foundation for Basic Research, China (Grant No. 20140101168JC).

Abstract

Using the time-dependent pseudo-spectral scheme, we solve the time-dependent Schrödinger equation of a hydrogen-like atom in a strong laser field in momentum space. The intensity-resolved photoelectron energy spectrum in above-threshold ionization is obtained and further analyzed. We find that with the increase of the laser intensity, the above-threshold ionization emission spectrum exhibits periodic resonance structure. By analyzing the population of atomic bound states, we find that it is the multi-photon excitation of bound state that leads to the occurrence of this phenomenon, which is in fairly good agreement with the experimental results.

Keyword: 32.80.Rm; 42.50.Hz; time-dependent pseudo-spectral scheme; above-threshold ionization; resonance structure
1. Introduction

With the rapid development of laser technology, the intensity of a laser field can be obtained to be close to or even larger than the intensity of the Coulomb field that the valence electrons in an atom can feel.[1] As the atoms and molecules interact with this intense laser pulse, the above-threshold ionization (ATI), high-order harmonic generation (HHG), non-sequential ionization, and other strong field phenomena can be observed.[26] ATI is one of the typical unperturbed physical phenomena.

Recently, quite a few new phenomena have been found in research on ATI. With the development of the ultrashort laser technology, the carrier envelope phase (CEP) of the laser pulse plays a very important role in the process of ionization. In experiments, the ATI photoelectron spectrum can be used to measure the CEP of the laser pulse, [7, 8] and detect the internal structures of atoms and molecules.[912] Utilizing the ultrashort intense infrared laser pulse[13, 14] to interact with atoms and molecules, new singular structures have been found in the low-energy part of ATI photoelectron spectrum, [1517] which is induced by the mid-infrared laser pulse and Coulomb long-range interaction.[18, 19] The development of the free electron laser and HHG technology provides a chance to study the interaction between high-frequency laser pulses. In the ATI photoelectron spectrum of the hydrogen atom under the action of high-frequency intense laser pulse it is found that with the increase of the intensity of the incident laser pulse, the peak shifts and a multi-peak structure occurs.[20, 21] This phenomenon can be attributed to the dynamic Stark effect and the interference between two electronic wave-packets.[2224] Li et al. discussed the resonance ionization via intensity-resolved ATI spectrum.[25] They used a quantum-trajectory Monte Carlo model to reproduce the measured photoelectron angular distribution of ATI.[26]

As the laser pulse interacts with the atom, the ground state electron may resonate with higher excited states, and a Freeman resonant peak can be observed in the low-energy region of photoelectron emission spectrum.[27, 28] Previous study has shown that this phenomenon is more sensitive to the CEP of the laser pulse, which can be used to detect the initial CEP of the laser pulse.[29, 30] This paper focuses on the study of the relationship between the resonance phenomenon, the structure of the ATI, and the population of the bound states of the system, which are observed within specific intensities.

To illuminate the above problems, according to the momentum space, we calculate the photoelectron spectra of hydrogen-like atom in different intensities by numerically solving the three-dimensional time-dependent Schrö dinger equation (TDSE), and find that with the increase of the intensity of incident laser pulse, the ATI peak structure changes from a single peak to a multi-peak and the resonant peak structure appears. However, if we further increase the intensity of the incident laser pulse, the ATI peak changes back to a single peak structure, and the resonant peak structure disappears, exhibiting a periodic resonance structure. By analyzing the population of the bound states, it is found that there exists a direct relation between the population of highly excited state and the resonance structure.

Atomic units are used throughout this paper unless otherwise indicated.

2. Theory

In order to investigate the photoelectron emission spectrum under the action of intense laser pulses through numerical simulation, we should first calculate the time-dependent wave-function of the electron in an intense laser field. In the coordinate space, the ionized electron will move to a zone which is far from the ions; therefore, in this case, larger computational space and more mesh grids are needed so as to make the results convergent. In this paper, based on the momentum space, the time-dependent pseudo-spectral scheme is used to solve the TDSE of the system, [31, 32] which can reduce the requirements for computational source on the promise of accurate calculation. In the velocity gauge, the coordinate space TDSE in the dipole approximation for an atom interacting with a laser field is given by

where c is the speed of light in vacuum, A(t) denotes the vector potential of the laser field, and U(r) represents the Coulomb potential of the atom. The momentum space wave function Φ (k, t) can be calculated from the Fourier transform of the coordinate space wave function Ψ (r, t) in the whole coordinate space as follows:

By applying Eq.  (2) to Eq.  (1), the TDSE in the momentum space can be obtained, which is an integro– differential equation given by

where V(k, k' ) is the momentum space Coulomb potential defined by

and for hydrogen-like system whose number of nuclear charge is Z, the above expression can be reduced to[33]

By utilizing the time-dependent pseudo-spectral scheme, we can solve Eq.  (3) and obtain the time-dependent wave-function at any moment. Using the wave-function of the atom at the end of the laser pulse, and projecting it onto the eigen-state in the field-free case, we can calculate the excitation and ionization of the atom under the action of an intense laser pulse. After obtaining the final-state wave function of the system, we calculate the information about the ATI by utilizing the scheme proposed by Dionissopoulou et al.[34] At the end of the laser pulse, we can obtain the population of each state by projecting the eigen-state wave-function in the field-free case onto the final-state wave function. By projecting the scattered continuum state

onto the wave function at the end of the laser pulse, if cl (ε , t) is the population amplitude of the normalized continuum state, the angular distribution of photoelectron can be obtained as

For the case of a linearly polarized laser pulse, it can be reduced to

and the corresponding single differential cross section is

3. Results and discussion

The laser pulse used in this paper is linearly polarized, whose incident direction is along the z axis, and whose electric field is

where A0 = E0/ω , E0 and ω are the peak amplitude and frequency of the laser electric field, respectively, and f(t) = sin2(ω t/2τ R) is the envelope of the laser’ s electric field, with τ R being the full width at half maximum (FWHM) of the laser pulse, which is selected as 10 optical cycles, and φ is the CEP, which is set to be 0. The wavelength of the incident laser pulse is 780  nm, and the ionization potential is 11.3  eV.

The intensity of the incident laser pulse changes from 2  TW/cm2 to 100  TW/cm2, and the variations of the total ionization and excitation with the intensity of the incident laser are investigated, and the results are shown in Fig.  1. It can be seen from the figure that with the increase of the incident intensity, the total ionization and excitation of atom increase on the whole; however, the peaked structures appear in some specific intensities. More importantly, for ionization and excitation, there exists a difference between their behaviors of peak structure: as for excitation, its peak structure is more obvious, and the intensities which the peak positions of the ionization and excitation correspond to are different, as shown by the arrow in the figure.

Fig.  1. Variations of the total ionization and excitation probability with the intensity of the incident laser pulse.

To further understand the change law of the ATI peak, we calculate the variation of the population of the bound state with the incident intensity, as shown in Fig.  2. In the figure, we present the variation of the probability of ionized electron with the incident intensity. It can be seen from the figure that the position where the peak of n photons disappears, calculated according to the formula, is in good agreement with that obtained from quantum calculation. More importantly, the disappearance of the ATI peaks continues the previous trend and extends to the negative energy area, and there is high population at the higher excited states as shown in the figure. We can also see from the figure that when population of the bound states is larger, the resonant multi-peak structure can be observed in the corresponding ATI peak, which is shown by the elliptic region in the figure. For the field strength at which the bound state of the atoms is resonance-enhanced, the population of the higher excited state is larger, and the corresponding total excitation of the atom is also enhanced, as shown in Fig.  1. For example, in Fig.  2, as the intensity changes between 60  TW/cm2 and 70  TW/cm2, the apparent multi-photon resonant excitation appears and the peak also appears in the corresponding product of total excitation, as shown by the solid black arrow in Fig.  1. It can be seen from Fig.  1 that the probability peak of the total excitation corresponds to the minimum of the total ionization. Therefore, we can deduce the cause which results in the difference between the intensities that the ionization and excitation peaks correspond to, i.e., as the higher excited states interact with the intense laser pulse, the energy level will move. Consequently, as the incident intensity meets the condition of multi-photon excitation, the excitation probability of an atom increases while the corresponding excitation probability is not that high. Therefore, as the incident intensity increases, it meets the condition of multi-photon ionization and the peaks will appear on the total ionization.

Fig.  2. Variations of population of continuous and bound states with the intensity of the incident laser.

In order to make a comparison with the experimental results, the volume effect of laser pulse is taken into consideration in our calculation.[35] In Fig.  3(a), we present the ATI spectra as the intensities of the laser are 16  TW/cm2, 24  TW/cm2, and 40  TW/cm2, with the volume effect considered, which corresponds to the solid black curve, dash red curve, and dash– dot blue curve in Fig.  3(a), respectively. We can observe the change of ATI peak structure from the figure. When the intensity is 16  TW/cm2, the ATI peak structure has a single-peak. When the intensity is 24  TW/cm2, the ATI peak structure does not have a single-peak but changes into a multi-peak structure. If we further increase the incident intensity, the ATI peak structure changes back to a single-peak and the resonant structure disappears, indicating that this kind of resonant structure strongly depends on the intensity. In the experiments, we observe the ATI spectra of formic acid molecule whose ionization energy is 11.3  eV, [36] as shown in Fig.  3(b). The multi-peak structure is shown by the solid green arrow in Fig.  3(b). It can be seen from the figure that with the increase of the incident intensity, the ATI spectrum becomes identical to the theoretical result. By analyzing the population of the atomic bound state, we find that it is the multi-photon excitation of the bound state, leading to the occurrence of this multi-peak structure.

In order to understand the variations of the excitation and ionization probability with the intensity of the incident laser, we calculate the variation of the photoelectron spectrum with the intensity of the incident laser, as shown in Fig.  4. It can be seen from the figure that with the increase of the incident intensity, the cut-off position increases gradually. Since the influence of AC-STARK effect of laser field on the continuous state is large while that on the ground state is small, more energy is needed to make the atom ionized. The peak position of the ATI photoelectron emission spectrum is EATI = n × ω IpUp, where Ip is the ionization potential of the atom, Up is the pondermotive energy, and n is the number of photons absorbed by the ground state electron. The Up increases linearly with the increase of the incident intensity. As a result, the ATI peak moves towards the low energy side. We present the position where the energy of the first ATI peak is zero as the number of photons absorbed by the ground state electron is n = 8– 10. With the increase of the intensity, the ATI peak exhibits the repetitive structures of single-peak and multi-peak, as shown in the panel on the right-hand side of Fig.  4.

Fig.  3. ATI spectra in different incident intensities, obtained from (a) the theoretical calculations and (b) the experimental measurements.[36]

Fig.  4. Change of photoelectron spectrum with the intensity of the incident laser.

In addition, we investigate the photoelectron spectra of the atom for different wavelengths of the incident laser, as shown in Fig.  5. The wavelengths of the incident laser are selected to be 760  nm (Fig.  5(a)), 780  nm (Fig.  5(b)), and 800  nm (Fig.  5(c)), respectively. It can be seen from the figure that, within the same intensity and photoelectron energy range, different wavelengths show similar results: with the increase of the incident intensity, the ATI peaks move towards the low-energy side, the ATI peak that corresponds to intensity in which the first peak reaches zero energy starts to change from a single-peak structure to a multi-peak structure, as shown by the elliptic region in the figure, indicating that it gradually plays a dominant role in the process of multi-photon excitation of the high Rydberg state. The difference between the change laws of ATI peaks at different wavelengths lies in the fact that the slopes of the changes of ATI peaks with incident intensity are different, reflecting that, in the same intensity, different frequencies correspond to different values of Up.

Fig.  5. Evolutions of the photoelectron emission spectrum with incident intensity when the wavelengths of the incident laser are selected to be (a) 760  nm, (b) 780  nm, and (c) 800  nm, respectively.

4. Conclusions

In this study, based on the momentum space, we numerically solve the TDSE accurately, and systematically investigate the ATI emission spectra at different incident intensities. With the increase of the incident intensity, the ATI emission spectra exhibit periodic resonance structures of single-peak structure and multi-peak structure, respectively. By analyzing the change of the probability of the eigen-state at the end of the laser pulse with the incident intensity, the occurrence of the multi-peak structure can be attributed to the multi-photon excitation. With the volume effect considered: with the change of the incident intensity, the calculated ATI spectrum is identical to the experimental results.

Acknowledgment

We thank the High Performance Computing Center of Jilin University for supercomputer time.

Reference
1 Protopapas M, Keitel C H and Knight P L 1997 Rep. Prog. Phys. 60 389 DOI:10.1088/0034-4885/60/4/001 [Cited within:1] [JCR: 13.232]
2 Agostini P, Fabre F, Mainfray G, Petite G and Rahman N K 1979 Phys. Rev. Lett. 42 1127 DOI:10.1103/PhysRevLett.42.1127 [Cited within:1] [JCR: 7.943]
3 McPherson A, Gibson G, Jara H, Johann U, Luk T S, McIntyre I A, Boyer K and Rhodes C K 1987 J. Opt. Soc. Am. B 4 595 DOI:10.1364/JOSAB.4.000595 [Cited within:1]
4 Walker B, Sheehy B, DiMauro L F, Agostini P, Schafer K J and Kuland er K C 1994 Phys. Rev. Lett. 73 1227 DOI:10.1103/PhysRevLett.73.1227 [Cited within:1] [JCR: 7.943]
5 Chen J, Zeng B, Liu X, Cheng Y and Xu Z Z 2009 New J. Phys. 11 113021 DOI:10.1088/1367-2630/11/11/113021 [Cited within:1] [JCR: 4.063]
6 Tian Y Y, Guo F M and Yang Y J 2013 Acta Phys. Sin. 62 073202 DOI:10.7498/aps.62.073202(in Chinese) [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
7 Tong X M, Hino K and Toshima N 2006 Phys. Rev. A 74 031405 DOI:10.1103/PhysRevA.74.031405 [Cited within:1] [JCR: 3.042]
8 Kling M F, Rauschenberger J, Verhoef A J, Hasovic E, Uphues T, MiloŠević D B, Muller H G and Vrakking M J J 2008 New J. Phys. 10 025024 DOI:10.1088/1367-2630/10/2/025024 [Cited within:1] [JCR: 4.063]
9 Kamta G L and Band rauk A D 2006 Phys. Rev. A 74 033415 DOI:10.1103/PhysRevA.74.033415 [Cited within:1] [JCR: 3.042]
10 Meckel M, Comtois D, Zeidler D, Staudte A, Pavičić D, Band ulet H C, Pépin H, Kieffer J C, Dörner R, Villeneuve D M and Corkum P B 2008 Science 320 1478 DOI:10.1126/science.1157980 [Cited within:1]
11 van der Zwan E V and Lein M 2012 Phys. Rev. Lett. 108 043004 DOI:10.1103/PhysRevLett.108.043004 [Cited within:1] [JCR: 7.943]
12 Blaga C I, Xu J L, Dichiara A D, Sistrunk E, Zhang K, Agostini P, Miller T A, DiMauro L F and Lin C D 2012 Nature 483 194 DOI:10.1038/nature10820 [Cited within:1] [JCR: 38.597]
13 Andriukaitis G, Balči&ūnas T, Ališauskas S, Pugžlys A, Baltuška A, Popmintchev T, Chen M C, Murnane M M and Kapteyn H C 2011 Opt. Lett. 36 2755 DOI:10.1364/OL.36.002755 [Cited within:1] [JCR: 3.385]
14 Agostini P and DiMauro L F 2008 Contemp. Phys. 49 179 DOI:10.1080/00107510802221630 [Cited within:1] [JCR: 2.614]
15 Reiss H R 1980 Phys. Rev. A 22 1786 DOI:10.1103/PhysRevA.22.1786 [Cited within:1]
16 Reiss H R 1987 J. Phys. B: At. Mol. Opt. Phys. 20 L79 DOI:10.1088/0022-3700/20/3/003 [Cited within:1] [JCR: 2.031]
17 Blaga C I, Catoire F, Colosimo P, Paulus G G, Muller H G, Agostini P and Dimauro L F 2009 Nat. Phys. 5 335 DOI:10.1038/nphys1228 [Cited within:1] [JCR: 19.352]
18 Quan W, Lin Z, Wu M, Kang H, Liu H, Liu X, Chen J, Liu J, He X T, Chen S G, Xiong H, Guo L, Xu H, Fu Y, Cheng Y and Xu Z Z 2009 Phys. Rev. Lett. 103 093001 DOI:10.1103/PhysRevLett.103.093001 [Cited within:1] [JCR: 7.943]
19 Liu C P, Nakajima T, Sakka T and Ohgaki H 2008 Phys. Rev. A 77 043411 DOI:10.1103/PhysRevA.77.043411 [Cited within:1] [JCR: 3.042]
20 Demekhin P V and Cederbaum L S 2012 Phys. Rev. Lett. 108 253001 DOI:10.1103/PhysRevLett.108.253001 [Cited within:1] [JCR: 7.943]
21 Demekhin P V, Hochstuhl D and Cederbaum L S 2013 Phys. Rev. A 88 023422 DOI:10.1103/PhysRevA.88.023422 [Cited within:1] [JCR: 3.042]
22 Demekhin P V and Cederbaum L S 2013 Phys. Rev. A 88 043414 DOI:10.1103/PhysRevA.88.043414 [Cited within:1] [JCR: 3.042]
23 Yu C, Fu N, Zhang G Z and Yao J Q 2013 Phys. Rev. A 87 043405 DOI:10.1103/PhysRevA.87.043405 [Cited within:1] [JCR: 3.042]
24 Yu C, Fu N, Hu T, Zhang G Z and Yao J Q 2013 Phys. Rev. A 88 043408 DOI:10.1103/PhysRevA.88.043408 [Cited within:1]
25 Li Q G, Tong X M, Morishita T, Wei H and Lin C D 2014 Phys. Rev. A 89 023421 DOI:10.1103/PhysRevA.89.023421 [Cited within:1] [JCR: 3.042]
26 Li M, Geng J W, Liu H, Deng Y K, Wu C Y, Peng L Y, Gong Q H and Liu Y Q 2014 Phys. Rev. Lett. 112 113002 DOI:10.1103/PhysRevLett.112.113002 [Cited within:1] [JCR: 7.943]
27 Freeman R R, Bucksbaum P H, Milchberg H, Darack S, Schumacher D and Geusic M E 1987 Phys. Rev. Lett. 59 1092 DOI:10.1103/PhysRevLett.59.1092 [Cited within:1] [JCR: 7.943]
28 Wang Y, Zhang J T, Xu Z Z, Wu Y S, Wang J T and Guo D S 2009 Phys. Rev. A 80 053417 DOI:10.1103/PhysRevA.80.053417 [Cited within:1] [JCR: 3.042]
29 Milošević D B, Paulus G G and Becke W 2003 Opt. Express 11 1418 DOI:10.1364/OE.11.001418 [Cited within:1] [JCR: 3.546]
30 Tian Y Y, Wei S S, Guo F M, Li S Y and Yang Y J 2013 Acta Phys. Sin. 62 153202 (in Chinese) DOI:10.7498/aps.62.153202 [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
31 Zhou Z Y and Chu S I 2011 Phys. Rev. A 83 013405 DOI:10.1103/PhysRevA.83.013405 [Cited within:1] [JCR: 3.042]
32 Tian Y Y, Li S Y, Wei S S, Guo F M, Zeng S L, Chen J G and Yang Y J 2014 Chin. Phys. B 23 053202 DOI:10.1088/1674-1056/23/5/053202 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
33 Land au R H 1983 Phys. Rev. C 27 2191 DOI:10.1103/PhysRevC.27.2191 [Cited within:1] [JCR: 3.715]
34 Dionissopoulou S, Mercouris T, Lyras A and Nicolaides C A 1997 Phys. Rev. A 55 4397 DOI:10.1103/PhysRevA.55.4397 [Cited within:1] [JCR: 3.042]
35 Morishita T, Chen Z G, Watanabe S and Lin C D 2007 Phys. Rev. A 75 023407 DOI:10.1103/PhysRevA.75.023407 [Cited within:1] [JCR: 3.042]
36 Wang C C, Tian Y Y, Luo S Z, Roeterdink W G, Yang Y J, Ding D J, Okunishi M, Prumper G, Shimada K, Ueda K and Zhu R 2014 Phys. Rev. A 90 023405 DOI:10.1103/PhysRevA.90.023405 [Cited within:1] [JCR: 3.042]