# Tunnel ionization

Tunnel ionization is a process in which electrons in an atom (or a molecule) pass through the potential barrier and escape from the atom (or molecule). In an intense electric field, the potential barrier of an atom (molecule) is distorted drastically. Therefore, the length of the barrier that electrons have to pass decreases and electrons can escape from the atom (molecule) easily. Tunneling Ionization is a QM phenomenon; as a non-zero probability event for observing a particle escaping from the deformed Coulomb potential barrier, obviously this phenomenon is forbidden by classical laws, as in the classical picture an electron does not have sufficient energy to escape.

When the atom is in an external DC field, the Coulomb potential barrier is lowered and the electron can tunnel through the barrier. In the case of an alternating electric field, the direction of the electric field reverses after the half period of the field. The ionized electron may come back to its parent ion. The electron may recombine with the nucleus (nuclei) and its kinetic energy is released as light (high harmonic generation). If the recombination does not occur, further ionization may proceed by collision between high-energy electrons and a parent atom (molecule). This process is known as non-sequential ionization.[1]

## DC tunneling ionization

Tunneling ionization from the ground state of a Hydrogen atom in an electrostatic (DC) field was solved schematically by Landau.[2] This provides a simplified physical system that given it proper exponential dependence of the ionization rate on the applied external field. When ${\displaystyle E<, the ionization rate for this system is given by:

${\displaystyle w=4\omega _{a}{\frac {E_{a}}{\left|E\right|}}\exp \left[-{\frac {2}{3}}{\frac {E_{a}}{\left|E\right|}}\right]}$

Landau expressed this in units where ${\displaystyle m=e=\hbar =1}$. In SI units the previous parameters can be expressed as:

${\displaystyle E_{a}={\frac {m^{2}e^{5}}{(4\pi \epsilon _{0})^{3}\hbar ^{4}}}}$,
${\displaystyle \omega _{a}={\frac {me^{4}}{(4\pi \epsilon _{0})^{2}\hbar ^{3}}}}$.

The ionization rate is the total probability current through the outer classical turning point. This is found using the WKB approximation to match the ground state hydrogen wavefunction though the suppressed coulomb potential barrier.

## AC electric field

The ionization rate of a hydrogen atom in an alternating electric field, like that of a laser, can be treated, in the appropriate limit, as the DC ionization rate averaged over a single period of the electric field's oscillation. Multiphoton and tunnel ionization of an atom or a molecule describe the same process by which a bounded electron, through the absorption of more than one photon from the laser field, is ionized. The difference between them is a matter of definition under different conditions. They can henceforth be called MPI (multiphoton ionization) whenever the distinction is not necessary. The dynamics of the MPI can be described by finding the time evolution of the state of the atom which is described by the Schrödinger equation.

Combined potential of an atom and a uniform laser field. At distances ${\displaystyle r, the potential of the laser can be neglected, while at distances with ${\displaystyle r>r_{0}}$, the Coulomb potential is negligible compared to the potential of the laser field. The electron emerges from under the barrier at ${\displaystyle r=R_{c}}$. ${\displaystyle E_{i}}$ is the ionization potential of the atom.

When the intensity of the laser is strong, the lowest-order perturbation theory is not sufficient to describe the MPI process. In this case, the laser field on larger distances from the nucleus is more important than the Coulomb potential and the dynamic of the electron in the field should be properly taken into account. The first work in this category was published by Keldysh.[3] He modeled the MPI process as a transition of the electron from the ground state of the atom to the Volkov states (the state of a free electron in the electromagnetic field[4]). In this model, the perturbation of the ground state by the laser field is neglected and the details of atomic structure in determining the ionization probability are not taken into account. The major difficulty with Keldysh's model was its neglect of the effects of Coulomb interaction on the final state of the electron. As is observed from the figure, the Coulomb field is not very small in magnitude compared to the potential of the laser at larger distances from the nucleus. This is in contrast to the approximation made by neglecting the potential of the laser at regions near the nucleus. Perelomov et al.[5][6] included the Coulomb interaction at larger internuclear distances. Their model (which is called the PPT model) was derived for short-range potential and includes the effect of the long-range Coulomb interaction through the first-order correction in the quasi-classical action. In the quasi-static limit, the PPT model approaches the ADK model.[7]

Many experiments have been carried out on the MPI of rare gas atoms using strong laser pulses, through measuring both the total ion yield and the kinetic energy of the electrons. Here, one only considers the experiments designed to measure the total ion yield. Among these experiments are those by Chin et al.,[8] Augst et al.[9] and Auguste et al.[10] Chin et al. used a 10.6 μm CO2 laser in their experiment. Due to the very small frequency of the laser, the tunneling is strictly quasi-static, a characteristic which is not easily attainable using pulses in the near infrared or visible region of frequencies. These findings weakened the suspicion on the applicability of models basically founded on the assumption of a structureless atom. Larochelle et al.[11] have compared the theoretically predicted ion versus intensity curves of rare gas atoms interacting with a Ti:sapphire laser with experimental measurement. They have shown that the total ionization rate predicted by the PPT model fits very well the experimental ion yields for all rare gases in the intermediate regime of Keldysh parameter.

## Analytical formula for the rate of MPI

( be careful, there are a lot of typos in the following section ) The dynamics of the MPI can be described by finding the time evolution of the state of the atom which is described by the Schrödinger equation. The form of this equation in the electric field gauge, assuming the single active electron (SAE) approximation and using dipole approximation, is the following

${\displaystyle i{\frac {\partial }{\partial t}}\Psi (\mathbf {r} ,\,t)=-{\frac {1}{2m}}\nabla ^{2}\Psi (\mathbf {r} ,\,t)+(\mathbf {E} (t)\cdot \mathbf {r} +V(\mathbf {r} ))\Psi (\mathbf {r} ,\,t)}$

where ${\displaystyle \mathbf {E} (t)}$ is the electric field of the laser and ${\displaystyle V(r)}$ is the static Coulomb potential of the atomic core at the position of the active electron. By finding the exact solution of equation (1) for a potential ${\displaystyle {\sqrt {2E_{i}}}.\delta (\mathbf {r} )}$ (${\displaystyle E_{i}}$ the magnitude of the ionization potential of the atom), the probability current ${\displaystyle \mathbf {J} (\mathbf {r} ,t)}$ is calculated. Then, the total MPI rate from short range potential for linear polarization, ${\displaystyle W(\mathbf {E} ,\omega )}$, is found from

${\displaystyle W(\mathbf {E} ,\omega )=\lim _{x\to \infty }\int _{0}^{\frac {2\pi }{\omega }}\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }\mathbf {J} (\mathbf {r} ,t)\,dz\,dy\,dt}$

where ${\displaystyle \omega }$ is the frequency of the laser, which is assumed to be polarized in the direction of the ${\displaystyle x}$ axis. The effect of the ionic potential, which behaves like ${\displaystyle {\frac {Z}{r}}}$ (${\displaystyle Z}$ is the charge of atomic or ionic core) at a long distance from the nucleus, is calculated through first order correction on the semiclassical action. The result is that the effect of ionic potential is to increase the rate of MPI by a factor of

${\displaystyle I_{PPT}=(2(2E_{i})^{\frac {3}{2}}/F)^{n^{*}}}$

Where ${\displaystyle n^{*}=Z/{\sqrt {2E_{i}}}}$ and ${\displaystyle F}$ is the peak electric field of laser. Thus, the total rate of MPI from a state with quantum numbers ${\displaystyle l}$ and ${\displaystyle m}$ in a laser field for linear polarization is calculated to be

${\displaystyle W_{PPT}=I_{PPT}W(\mathbf {E} ,\omega )=|C_{n^{*}l^{*}}|^{2}{\sqrt {\frac {6}{\pi }}}f_{lm}E_{i}(2(2E_{i})^{\frac {3}{2}}/F)^{2n^{*}-|m|-3/2}(1+\gamma ^{2})^{|m/2|+3/4}A_{m}(\omega ,\gamma )e^{-{\frac {2}{3}}g(\gamma )(2E_{i})^{\frac {3}{2}}/F}}$

where ${\displaystyle \gamma ={\frac {\omega {\sqrt {2E_{i}}}}{F}}}$ is the Keldysh`s adiabaticity parameter and ${\displaystyle l^{*}=n^{*}-1}$ .The coefficients ${\displaystyle f_{lm}}$, ${\displaystyle g(\gamma )}$ and ${\displaystyle C_{n^{*}l^{*}}}$ are given by

${\displaystyle f_{lm}={\frac {(2l+1)(l+|m|){!}}{2^{|m|}|m|{!}(l-|m|){!}}}}$
${\displaystyle g(\gamma )={\frac {3}{2\gamma }}((1+{\frac {1}{2\gamma ^{2}}})\sinh ^{-1}(\gamma )-{\frac {\sqrt {1+\gamma ^{2}}}{2\gamma }})}$
${\displaystyle |C_{n^{*}l^{*}}|^{2}={\frac {2^{2n^{*}}}{n^{*}\Gamma (n^{*}+l^{*}+1)\Gamma (n^{*}-l^{*})}}}$

The coefficient ${\displaystyle A_{m}(\omega ,\gamma )}$ is given by

${\displaystyle A_{m}(\omega ,\gamma )={\frac {4}{\sqrt {3\pi }}}{\frac {1}{|m|!}}{\frac {\gamma ^{2}}{1+\gamma ^{2}}}\sum _{n>v}^{\infty }e^{-(n-v)\alpha (\gamma )}w_{m}\left({\sqrt {{\frac {2\gamma }{\sqrt {1+\gamma ^{2}}}}(n-v)}}\right)}$,

where

${\displaystyle w_{m}(x)=e^{-x^{2}}\int _{0}^{x}(x^{2}-y^{2})^{m}e^{y^{2}}\,dy}$
${\displaystyle \alpha (\gamma )=2(\sinh ^{-1}(\gamma )-{\frac {\gamma }{\sqrt {1+\gamma ^{2}}}})}$
${\displaystyle v={\frac {E_{i}}{\omega }}(1+{\frac {1}{2\gamma ^{2}}})}$

The ADK model is the limit of the PPT model when ${\displaystyle \gamma }$ approaches zero (quasi-static limit). In this case, which is known as quasi-static tunnelling (QST), the ionization rate is given by

${\displaystyle W_{ADK}=|C_{n^{*}l^{*}}|^{2}{\sqrt {\frac {6}{\pi }}}f_{lm}E_{i}(2(2E_{i})^{\frac {3}{2}}/F)^{2n^{*}-|m|-3/2}e^{-(2(2E_{i})^{\frac {3}{2}}/3F)}}$.

In practice, the limit for the QST regime is ${\displaystyle \gamma <1/2}$. This is justified by the following consideration.[12] Referring to the Figure, the ease or difficulty of tunneling can be expressed as the ratio between the equivalent classical time it takes for the electron to tunnel out the potential barrier while the potential is bent down. This ratio is indeed ${\displaystyle \gamma }$, since the potential is bent down during half a cycle of the field oscillation and the ratio can be expressed as

${\displaystyle \gamma ={\frac {\tau _{T}}{{\frac {1}{2}}\tau _{L}}}}$,

where ${\displaystyle \tau _{T}}$ is the tunneling time (classical time of flight of an electron through a potential barrier, and ${\displaystyle \tau _{L}}$ is the period of laser field oscillation.

## MPI of molecules

Contrary to the abundance of theoretical and experimental work on the MPI of rare gas atoms, the amount of research on the prediction of the rate of MPI of neutral molecules was scarce until recently. Walsh et al.[13] have measured the MPI rate of some diatomic molecules interacting with a 10.6 μm CO2 laser. They found that these molecules are tunnel-ionized as if they were structureless atoms with an ionization potential equivalent to that of the molecular ground state. Talebpour et al.[14][15] were able to quantitatively fit the ionization yield of diatomic molecules interacting with a Ti:sapphire laser pulse. The conclusion of the work was that the MPI rate of a diatomic molecule can be predicted from the PPT model by assuming that the electron tunnels through a barrier given by ${\displaystyle {\frac {Z_{eff}}{r}}}$ instead of barrier ${\displaystyle {\frac {1}{r}}}$ which is used in the calculation of the MPI rate of atoms. The importance of this finding is in its practicality; the only parameter needed for predicting the MPI rate of a diatomic molecule is a single parameter, ${\displaystyle Z_{eff}}$. Using the semi-empirical model for the MPI rate of unsaturated hydrocarbons is feasible.[16] This simplistic view ignores the ionization dependence on orientation of molecular axis with respect to polarization of the electric field of the laser, which is determined by the symmetries of the molecular orbitals. This dependence can be used to follow molecular dynamics using strong field multiphoton ionization.[17]

## Tunneling time

The question of how long a tunneling particle spends inside the barrier region has remained unresolved since the early days of quantum mechanics. It is sometimes suggested that the tunneling time is instantaneous because both the Keldysh and the closely related Buttiker-Landauer[18] times are imaginary (corresponding to the decay of the wavefunction under the barrier). In a recent publication[19] the main competing theories of tunneling time are compared against experimental measurements using the attoclock in strong laser field ionization of helium atoms. Refined attoclock measurements reveal a real and not instantaneous tunneling delay time over a large intensity regime. It is found that the experimental results are compatible with the probability distribution of tunneling times constructed using a Feynman Path Integral (FPI) formulation.[20][21]

## References

1. ^ Corkum, Paul B. "Plasma perspective on strong field multiphoton ionization."Physical Review Letters 71.13 (1993): 1994.
2. ^ L.D. Landau and E.M. Lifshitz, Quantum Mechanics (Pergamon, New York, 1965), 2nd ed., pg 276.
3. ^ Keldysh L V 1965 Soviet Phys. JETP 2354
4. ^ Volcov D M 1934 Z. Phys. 94 250
5. ^ Perelemov A M, Popov V S and Terent'ev M V 1966 SovietPhys. JETP, 23 924
6. ^ Perelemov A M and Popov V S 1967 Soviet Phys.JETP, 25 336
7. ^ Ammosov M V, Delone N B and Krainov V P 1986 SovietPhys. JETP, 64 1191
8. ^ Chin S L, Yergeau F and Lavigne P 1985 J. Phys.B 18 L213
9. ^ Augst S, Meyerhofer D D, Strickland D and Chin SL 1991 JOSA B 8 858
10. ^ AugusteT, Monot P, Lommpre L A, Mainfray G and Manus C 1992 J. Phys. B 254181
11. ^ Larochelle S, Talebpour A and Chin S L 1998 J.Phys. B 31 1201
12. ^ Chin, S. L. "From multiphoton to tunnel ionization." Advances in Multiphoton Processes and Spectroscopy 16 (2004): 249.
13. ^ Walsh T D G, Decker J E and Chin S L 1993 J. Phys. B 26 L85
14. ^ Talebpour A, Larochelle A and Chin S L 1998 J.Phys. B 31 L49
15. ^ A. Talebpour, J. Yang, S.L. Chin, Semi-empirical model for the rate of tunnel ionization of N2 and O2 molecule in an intense Ti:sapphire laser pulse, Optics Communications, 1999, 163:29–32
16. ^ Talebpour A, Larochelle A and Chin S L, Multiphoton ionization of unsaturated hydrocarbons, 1998 J. Phys. B: At. Mol. Opt. Phys. 31 2769
17. ^ A. Jaron-Becker, Molecular Dynamics in Strong Laser Fields, IEEE JOURNAL OF SELECTED TOPICS IN QUANTUM ELECTRONICS, 2012, 18:105-112
18. ^ Buttiker, M. and R. Landauer, Traversal Time for Tunneling. Physical Review Letters, 1982. 49(23): p. 1739–1742.
19. ^ Landsman, Alexandra, et al. "Tunneling Time in Ultrafast Science is Real and Probabilistic." arXiv preprint arXiv:1301.2766 (2013).
20. ^ Fertig, H.A., Traversal-Time Distribution and the Uncertainty Principle in Quantum Tunneling. Physical Review Letters, 1990. 65(19): p. 2321-2324.
21. ^ Yamada, N., Unified derivation of tunneling times from decoherence functionals.Physical Review Letters, 2004. 93(17).