Tunnel ionization

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 describes 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 multiphoton ionization (MPI) 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.

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 Leonid Keldysh. 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). 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. A. M. Perelomov, V. S. Popov and M. V. Terent'ev included the Coulomb interaction at larger internuclear distances. Their model (which is called the PPT model after their initials) 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 by M. V. Ammosov, N. B. Delone, and V. P. Krainov.

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 S. L. Chin et al., S. Augst et al. and T. Auguste et al. Chin et al. used a 10.6 μm CO₂ laser in their experiment. Due to the very small frequency of the laser, the tunneling is strictly quasi-static, a characteristic that 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. S. Larochelle et al. 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

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{\mathrm{\partial }}{\mathrm{\partial }t}\mathrm{\Psi }(\mathbf{r},t)=-\frac{1}{2m}{\mathrm{\nabla }}^2\mathrm{\Psi }(\mathbf{r},t)+(\mathbf{E}(t)\cdot \mathbf{r}+V(\mathbf{r}))\mathrm{\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_{\text{i}}}.\delta (\mathbf{r})}$ (⁠${\displaystyle E_{\text{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 )=\underset{x\to \mathrm{\infty }}{lim}{\int }_0^{\frac{2\pi }{\omega }}{\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}{\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\mathbf{J}(\mathbf{r},t)dzdydt}$

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 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 semi-classical action. The result is that the effect of ionic potential is to increase the rate of MPI by a factor of

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

Where ${\displaystyle n^{\ast }=Z/\sqrt{2E_{\text{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_{\text{PPT}}=I_{\text{PPT}}W(\mathbf{E},\omega )=|C_{n^{\ast }{l}^{\ast }}{|}^2\sqrt{\frac{6}{\pi }}{f}_{lm}{E}_i(2(2E_{\text{i}}{)}^{\frac{3}{2}}/F{)}^{2n^{\ast }-|m|-3/2}(1+{\gamma }^2{)}^{|m/2|+3/4}{A}_m(\omega ,\gamma )e^{-\frac{2}{3}g(\gamma )(2E_{\text{i}}{)}^{\frac{3}{2}}/F}}$

where ${\displaystyle \gamma =\frac{\omega \sqrt{2E_{\text{i}}}}{F}}$ is the Keldysh's adiabaticity parameter and ⁠${\displaystyle l^{\ast }=n^{\ast }-1}$⁠. The coefficients ${\displaystyle f_{lm}}$, ${\displaystyle g(\gamma )}$ and ${\displaystyle C_{n^{\ast }{l}^{\ast }}}$ 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^{\ast }{l}^{\ast }}{|}^2=\frac{2^{2n^{\ast }}}{n^{\ast }\mathrm{\Gamma }(n^{\ast }+l^{\ast }+1)\mathrm{\Gamma }(n^{\ast }-l^{\ast })}}$

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}^{\mathrm{\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_{\text{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_{\text{ADK}}=|C_{n^{\ast }{l}^{\ast }}{|}^2\sqrt{\frac{6}{\pi }}{f}_{lm}{E}_i(2(2E_{\text{i}}{)}^{\frac{3}{2}}/F{)}^{2n^{\ast }-|m|-3/2}{e}^{-(2(2E_{\text{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. 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 }_{\text{T}}}{\frac{1}{2}{\tau }_{\text{L}}}}$,

where ${\displaystyle {\tau }_{\text{T}}}$ is the tunneling time (classical time of flight of an electron through a potential barrier, and ${\displaystyle {\tau }_{\text{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. have measured the MPI rate of some diatomic molecules interacting with a 10.6 μm CO₂ 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. A. Talebpour et al. 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 Z_{\text{eff}}/r}$ instead of barrier ${\displaystyle 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_{\text{eff}}}$⁠. Using the semi-empirical model for the MPI rate of unsaturated hydrocarbons is feasible. 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.

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 times are imaginary (corresponding to the decay of the wavefunction under the barrier). In a recent publication 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. However, later work in atomic hydrogen has demonstrated that most of the tunneling time measured in the experiment is purely from the long-range Coulomb force exerted by the ion core on the outgoing electron.