# Classical simulations including electron correlations for sequential double ionization

## Abstract

With a classical ensemble model that including electron correlations during the whole ionization process, we investigated strong-field sequential double ionization of Ar by elliptically polarized pulses at the quantitative level. The experimentally observed intensity-dependent three-band or four-band structures in the ion momentum distributions are well reproduced with this classical model. More importantly, the experimentally measured ionization time of the second electrons [A. N. Pfeiffer et al., Nature Phys. 7, 428 (2011)], which can not be predicted by the standard independent-electron model, is quantitatively reproduced by this fully classical correlated model. The success of our work encourages classical description and interpretation of the complex multi-electron effects in strong field ionization where nonperturbative quantum approaches are currently not feasible.

###### pacs:

32.80.Rm, 31.90.+s, 32.80.FbAmong various intense laser-induced phenomena, strong field double ionization (DI) is one of the most important and fundamental processes. During the past decades, a great number of experimental as well as theoretical studies have been performed on this area. It has been known that DI proceeds either sequentially or nonsequentially. In nonsequential double ionization (NSDI), the second electron is ionized by the recollision of the first tunneled electron (1). Because of this recollision, the two electrons from NSDI exhibit a highly correlated behavior (2); (3); (4); (5); (6). In sequential double ionization (SDI), it is usually assumed that no correlation exists between the two electrons, and thus the ionization of the electrons can be treated as two independent tunneling-ionization steps. However, this assumption has been called into doubt by recent experiments (7); (8). In Ref. (7), it has been shown that there is a clear angular correlation between the two electrons from SDI, which implies that the successive ionization steps are not independent in SDI. In Ref. (8), the authors found that the ionization time of the second electron from SDI is much earlier than the prediction of the independent-electron model. These observations declare that the electron correlations in SDI should be reexamined carefully.

Theoretically, an accurate description of the electron correlations in DI needs full quantum theory. However, due to the enormous computational demand, the solution of the time-dependent schrödinger equation of the two-electron system in the strong field, especially in the case of elliptically polarized laser fields, is not feasible currently. Instead, a fully classical treatment of the two- and multi-electron systems proposed by Eberly et al. has been well established (9); (10); (11); (12). During the past decade, this model has been successful in exploring the strong-field ionization processes at the qualitative level. However, it fails when made a quantitative comparison with experiments. For example, in ref. (13) it has been shown that the saturation intensity for SDI of Ar in the classical simulation is much higher than the experimental data (8), and the amplitude of oscillation in the ratio of the parallel to antiparallel emitted SDI counts as a function of intensity is much larger than the experimental observation (14). It remains a question: can a classical treatment describe the strong field processes at the quantitative level (15)? In this Letter, we gave an affirmative answer. With a delicate modification to Eberly’s model, we performed a first quantitative simulation on strong field SDI of Ar in the close-to-circular laser fields. Our numerical results exhibit intensity-dependent three-band or four-band structures in the ion momentum distributions, which are consistent with the recent experimental results (8). Especially, the experimentally measured ionization time of the second electron (8), which can not be predicted by the standard independent-electron model, is reproduced surprisingly well with this purely classical model which fully takes into account the electron correlations during the entire SDI process. The quantitative agreement between our numerical results and experimental data indicates that a classical treatment is a good approximation in describing the strong-field SDI and has the potential to shed light on the subtle multi-electron effect in strong-field double and multiple ionizations.

In the classical picture of a two-electron system, one electron often drops deeply into the nuclear potential well, leading to the autoionization of the other electron. In Eberly’s classical model (9); (10); (16), a soft-core potential is introduced for the ion-electron interaction to avoid autoionization. However, in this soft-core potential classical model (SPCM), the first and second ionization potentials can not be matched with those of the investigated target. In fact, the first electron often ionizes leaving the second electron with an energy much lower than the second ionization potential of the target (17); (18). In tunneling and over-the-barrier ionization, the ionization rate is very sensitive to the ionization energy. Thus, the SPCM may be deficient in describing strong field double ionization that occurs through tunneling or over-the-barrier escape. However, this deficiency can be avoided with the Heisenberg-core potential, which not only prevents autoionization but also gives the ground-configuration energies of the multi-electron atoms (19). In the past decades, this potential has been extensively employed in the classical investigations of atomic and molecular collisions and laser-matter interactions (20); (21). Here, we employ the Heisenberg-core potential instead of the soft-core potential in the classical two-electron model to study SDI of Ar by the elliptical laser pulses.

The Hamiltonian of the two-electron atom in the Heisenberg-core potential classical model (HPCM) is (atomic units are used throughout this Letter until stated otherwise):

(1) |

where and are the position and canonical momentum of the th electron, respectively. is the Heisenberg-core potential, which is expressed as (19):

(2) |

The parameter indicates the rigidity of the Heisenberg core and is chosen to be 2. For a given , the parameter is chosen to match the second ionization potential of the target, i.e., it is set to make the minimum of the one-electron Hamiltonian [] equal to the second ionization potential of Ar (-1.01 a.u.). For =2 we obtain =1.225. For the two-electron atom, in refs. (19); (20); (21), the ground-state energy of the system is obtained by minimizing the Hamiltonian . In that configuration, the two electrons locate at opposite sides of the nucleus, and they are stationary to each other (19). Differently, in our calculation, we employed the approach used by Eberly et al. (16) to determine the ground-state energy of the atom, i.e., by inputting an energy that equals to the sum of the first and the second ionization potentials of the target (-1.59 a.u. for Ar). The initial distributions of the ground-state atom in the phase space are obtained with the approach in SPCM (10); (11). Figures. 1(a) and 1(b) display the position distributions of the ground-state atom. We state that in our calculation, the second ionization potential is related to the minimum of while the ground-state energy of the atom is determined by the “input” energy, which is higher than the minimal value of . Thus, the two electrons are not fixed at two certain points but distributed in a finite region of the phase space [see figs. 1(a) and 1(b)]. This relaxation of the position of the electron pairs in the phase space allows sufficient electron correlations in the initial state and during the ionization of the first electron.

The Hamiltonian of the two-electrons in the presence of the laser field is

(3) |

where is the electric field of the laser pulses. The evolution of the system in the laser field is determined by the following equations:

(4) |

The electric field is given as , where f(t)= is the field envelope. , and are the laser frequency, the ellipticity and carrier-envelope phase (CEP), respectively. denotes the pulse duration (FWHM).

First of all, we perform a calculation to test how well the first and second ionization potentials of Ar have been produced by this HPCM. The two-electron atoms are exposed to a 33-fs, 788-nm laser pulse with the intensity I=4.0 PW/cm and ellipticity . We trace the energy evolutions of the two electrons and record the energy of the second electron when the first electron is ionized. This energy is assumed to be the second ionization potential of the model atom, which is shown in fig. 1(c). The first ionization potential of the model atom is also shown in fig. 1(c), obtained by subtracting the second ionization potential from the initial energy of the two-electron atom (-1.59 a.u.). It is clearly shown that the first and second ionization potentials of the model atom are about -0.57 a.u. and -1.02 a.u. (almost at the bottom of Hamiltonian ), respectively, very close to the realistic first and second ionization potentials of Ar.

A recent experiment has shown that in SDI of Ar by strong elliptical laser pulses the ratio of the parallel and antiparallel electron emissions along the minor elliptical axis exhibits an oscillating behavior as a function of laser intensity (14). This behavior has been predicted by the SPCM and explained as a multi-electron effect beyond independent-electron assumption (13). We display our calculations with the HPCM in fig. 2, where the oscillating behavior is also clear. It implies that the electron correlations in SDI are included and represented in this HPCM. One can clearly see from fig. 2 that at the high laser intensities the oscillating curve is a bit below 1, meaning the two electrons prefer to emit into the opposite hemispheres. This behavior well agrees with the experimental data (see fig. 2(a) of ref. (14)).

In figs. 3(a) and 3(d) we display the ion momentum distributions in the polarization plane for SDI of Ar by 33-fs laser pulses with laser intensities 1.0 PW/cm and 4.0 PW/cm, respectively. The ellipticity and wavelength are respectively 0.77 and 788 nm, same as those in the recent experiment (8). In the direction of the minor elliptical axis ( axis), the distribution exhibits a three-band structure at the relatively low laser intensity while a four-band structure at the relatively high laser intensity. In the direction of the major elliptical axis ( axis), as shown in figs. 3(c) and 3(f), the spectra show a gaussian shape for both intensities. These results agree well with the experiment (8). The origin of the four-band structure is interpreted as due to the different values of the electric field at which the two electrons are released (22); (23). Because of the ellipticity, the electrons preferentially ionize along the major axis, leading to the electrons with final momenta along the minor axis. The two outer bands in fig. 3(d) correspond to the events where the two electrons emit into the parallel directions whereas the two inner bands result from the events where the two electrons release into the antiparallel directions. At the relatively low laser intensity, the momentum amplitudes of the two electrons are almost the same because both electrons are ionized around the pulse centre (8). Thus the antiparallel electron emissions result in the nearly zero momentum of ion, leading to the three-band structure in ion momentum distribution [fig. 3(a)]. The above picture for the intensity-dependent ion momentum spectra has been confirmed by back analyzing the SDI trajectories of our calculations (not shown here).

Figure. 4(a) displays the correlated radial momentum distribution of the electron pairs from SDI of Ar. The laser parameters are the same as those in fig. 3(d). Note that in this calculation the focal volume effect assuming a Gaussian beam profile has been considered. The momentum in fig. 4(a) is defined as , which is an injective function of time under the condition that the electron ionized before the peak of the pulse (8). The repulsion behavior along the diagonal means that the two electrons achieve different final momenta at the end of the pulses, implying the different release times of the two electrons. In ref. (8), the ionization times of the two electrons in SDI are read from the electron’s final momentum . It is found that the ionization time for the first electron agrees well with the prediction of the independent-electron model. However, the ionization of the second electrons occurs much earlier than the prediction of the independent-electron model. With the same procedure as that in ref. (8), we extract the release times of two electrons in our classical simulations where the electron correlations are included during the entire process of SDI. The results are shown in fig. 4(b), where the experimental data from ref. (8) are also displayed. Note that for our numerical results, the laser intensities in fig. 4(b) are scaled with a constant factor of 0.82, which is well in the experimental uncertainty range. Surprisingly, the release times from our classical simulations agree excellently well with the experimental data, both for the first and the second electrons. For comparison, we have repeated the calculations above with the SPCM. The ionization times of the first and the second electrons from SPCM [as shown in fig. 4(b)] deviate seriously from the experimental data. Here, we only displayed the results at intensities above 4.0 PW/cm because in the SPCM the DI yield is very low at relatively low intensities (the DI probability does not exceed 1% until the laser intensity reaches 4.0 PW/cm). During the past decades, classical methods have been widely employed to describe the strong field phenomena. However, they are standing at the qualitative point of view. As far as we know, our calculation is the first successful example of quantitative investigation of the strong field DI with a classical treatment.

For the theoretical treatment of the independent-electron model in ref. (8), the Coulomb interaction between the ion and the escaping electron is neglected. Here we estimate how strong the Coulomb interaction affects the release-time reading. It has been demonstrated that the angular distributions of the electrons can serve as a signal to estimate the importance of the Coulomb interaction on the electrons (24); (25). In fig. 5 we display the angular spectra of the first and the second electrons. Note that the angular spectra peak at 180 (and 0) if there is no interaction between the ion and the escaping electron. Fig. 5 shows that the spectrum of the first electron (the red curves) peaks at an angle slightly deviating from 180, indicating a weak Coulomb correction on the electron trajectory (24); (25). For the second electron (the green curves), however, the spectrum peaks almost at 180, implying negligible Coulomb interaction between the ion and escaping electron. These results indicate that the influence of the Coulomb attraction on the second electron is even weaker than that on the first electron. Thus, it confirms that the Coulomb attraction has negligible contribution to the deviation between the experimental measurement and independent-electron prediction for the ionization time of the second electron reported in Ref. (8).

It has been proposed that inelastic tunneling (26); (27), where the first electron escapes leaving the core in an excited state, is important in double and multiple ionization at high laser intensities. However, the most typical picture of the ionization process revealed by our classical model is that the first electron ionizes leaving the second electron at the ground sate of Ar, as shown in fig. 1(d) [see the inset of fig. 1(d)]. Note that the second electron stays at a state with an energy higher than -1.01 a.u. after the first ionization for a small part of the SDI trajectories shown in figure 1(c). This “excitation” is responsible for the oscillation behavior of the parallel-antiparallel ratio (see fig. 2) (13). However, it is not responsible for the early release time of the second electron because the second ionization time shown in fig. 4(b) does not change when these trajectories are excluded. Another reason for the deviation of the independent-electron prediction of the second ionization from the measured data (8) possibly ascribes to the empirical formula for ionization rate (28), which may be inaccurate in the experimental condition where nonadiabatical effect is important to the ionization process (29). Then, it is not clear to what extent the deviation of the second ionization time comes from the inaccuracy of the ionization rate formula and to what extent the electron correlations influence the ionization of the second electron. These questions call for more detailed investigations.

In conclusion, we have investigated SDI of Ar by the elliptical laser pulses with the HPCM in which the electron correlation is included during the entire process. The experimental observed ion momentum spectra and oscillating behavior of the ratio of the antiparallel to parallel electron emissions are well reproduced by this classical model. Especially, the measured ionization time of the second electron, which strongly deviates from prediction of the standard independent-electron model, is excellently reproduced by the HPCM. The quantitative agreement between our classical calculations and experimental results provides strong supports to the classical treatment of the multi-electron processes induced by strong laser fields, which is currently indispensable because the nonperturbative quantum treatments of the complex effect are not feasible.

We acknowledge helpful discussions with Dr. A. N. Pfeiffer, Dr. X. Wang and Prof. J. H. Eberly. We thanks Dr. Pfeiffer for providing us with their experimental data. This work was supported by the National Science Fund (No.60925021 and No. 11004070) and the 973 Program of China (No. 2011CB808103).

### Footnotes

- Corresponding author: lupeixiang@mail.hust.edu.cn

### References

- P. B. Corkum, Phys. Rev. Lett. 71, 1994 (1993).
- Th. Weber et al., Nature 405, 658 (2000).
- A. Rudenko et al., Phys. Rev. Lett. 93, 253001 (2004).
- A. Becker, R. Dörner, and R. Moshammer, J. Phys. B 38, S753 (2005).
- Y. Zhou, Q. Liao, and P. Lu, Opt. Express 18, 16025 (2010).
- Y. Zhou et al., Opt. Express 19, 2301 (2011).
- A. Fleischer et al., Phys. Rev. Lett. 107, 113003 (2011).
- A. N. Pfeiffer et al., Nature Phys. 7, 428 (2011).
- Phay J. Ho et al., Phys. Rev. Lett. 94, 093002 (2005).
- S. L. Haan et al., Phys. Rev. Lett. 97, 103008 (2006).
- Y. Zhou, Q. Liao, and P. Lu, Phys. Rev. A 80, 023412 (2009).
- Y. Zhou, Q. Liao, and P. Lu, Phys. Rev. A 82, 053402 (2010).
- X. Wang and J. H. Eberly, arXiv:1102.0221v1 (2011)
- A. N. Pfeiffer et al., New. J. Phys. 13, 093008 (2011).
- K. Ueda and K. L. Ishikawa, Nature Phys. 7, 371 (2011).
- R. Panfili, S. L. Haan, and J. H. Eberly, Phys. Rev. Lett. 89, 113001 (2002).
- D. Bauer, Phys. Rev. A 56, 3028 (1997).
- Y. Zhou et al., Chin. Phys. Lett. 25, 3950 (2008).
- C. L. Kirschbaum and L. Wilets, Phys. Rev. A 21, 834 (1980).
- D. Zajfman and D. Maor, Phys. Rev. Lett. 56, 320 (1986).
- P. B. Lerner, K. J. LaGattuta, and J. S. Cohen, Phys. Rev. A 49, R12 (1994).
- C. M. Maharjan et al., Phys. Rev. A 72, 041403R(2005).
- X. Wang and J. H. Eberly, Phys. Rev. Lett. 103, 103007 (2009).
- S. V. Popruzhenko, G. G. Paulus, and D. Bauer, Phys. Rev. A 77, 053409 (2008).
- A. N. Pfeiffer et al., Nature Phys. 8, 76 (2012).
- A. S. Kornev, E. B. Tulenko, and B. A. Zon, Phys. Rev. A 68, 043414 (2003).
- W. A. Bryan et al., Nature Physics 2, 379 (2006).
- X. M. Tong and C. D. Lin, J. Phys. B 38, 2593 (2005).
- I. Barth and O. Smirnova, Phys. Rev. A 84, 063415 (2011).