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, as the length of the barrier that electrons have to pass decreases, the electrons can escape from the atom's potential more easily. Tunneling Ionization is a quantum mechanical phenomenon, since in the classical picture an electron does not have sufficient energy to overcome the potential barrier of the atom.

When the atom is in a DC external field, the Coulomb potential barrier is lowered and the electron has an increased, non-zero probability of tunnelling through the potential 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] using parabolic coordinates. This provides a simplified physical system that given it proper exponential dependence of the ionization rate on the applied external field. When ${\displaystyle E<<E_{a}}$, 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 atomic 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.

A more physically meaningful form for the ionization rate above can be obtained by noting that the Bohr radius and Hydrogen atom ionization energy are given by

${\displaystyle a_{0}={\frac {4\pi \epsilon _{0}\hbar ^{2}}{me^{2}}}}$,

${\displaystyle E_{ion}=R_{H}={\frac {me^{4}}{8\epsilon _{0}^{2}h^{2}}}}$,

where ${\displaystyle R_{H}\approx \mathrm {13.6\,eV} }$ is the Rydberg energy. Then, the parameters ${\displaystyle E_{a}}$ and ${\displaystyle \omega _{a}}$ can be written as

${\displaystyle E_{a}={\frac {2R_{H}}{ea_{0}}}}$, ${\displaystyle \omega _{a}={\frac {2R_{H}}{\hbar }}}$.

so that the total ionization rate can be rewritten

${\displaystyle w=8{\frac {R_{H}}{\hbar }}{\frac {2R_{H}/a_{0}}{\left|eE\right|}}\exp \left[-{\frac {4}{3}}{\frac {R_{H}/a_{0}}{\left|eE\right|}}\right]}$.

This form for the ionization rate ${\displaystyle w}$ emphasizes that the characteristic electric field needed for ionization ${\displaystyle E_{a}={\frac {2E_{ion}}{ea_{0}}}}$ is proportional to the ratio of the ionization energy ${\displaystyle E_{ion}}$ to the characteristic size of the electron's orbital ${\displaystyle a_{0}}$. Thus, atoms with low ionization energy (such as alkali metals) with electrons occupying orbitals with high principal quantum number ${\displaystyle n}$ (i.e. far down the periodic table) ionize most easily under a DC field. Furthermore, for a Hydrogenic atom, the scaling of this characteristic ionization field goes as ${\displaystyle Z^{3}}$, where ${\displaystyle Z}$ is the nuclear charge. This scaling arises because the ionization energy scales as ${\displaystyle \propto Z^{2}}$ and the orbital radius as ${\displaystyle \propto Z^{-1}}$. More accurate and general formulas for the tunneling from Hydrogen orbitals can also be obtained.[3]

As an empirical point of reference, the characteristic electric field ${\displaystyle E_{a}}$ for the ordinary Hydrogen atom is about ${\displaystyle 51\mathrm {\frac {Volt}{Angstrom}} }$ (or ${\displaystyle 5.1\cdot 10^{3}\,\mathrm {MV/cm} }$) and the characteristic frequency ${\displaystyle \omega _{a}}$ is ${\displaystyle 4.1\cdot 10^{4}\,\mathrm {THz} }$.

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<r_{0}}$, 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.[4] 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[5]). 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.[6][7] 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.[8]

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.,[9] Augst et al.[10] and Auguste et al.[11] 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 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.[12] 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.[13] 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.[14] 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.[15][16] 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.[17] 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.[18]

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[19] times are imaginary (corresponding to the decay of the wavefunction under the barrier). In a recent publication[20] 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.[21][22]

References

1. Corkum, P. B. (1993-09-27). "Plasma perspective on strong field multiphoton ionization". Physical Review Letters. American Physical Society (APS). 71 (13): 1994–1997. doi:10.1103/physrevlett.71.1994. ISSN 0031-9007. PMID 10054556.
2. L.D. Landau and E.M. Lifshitz, Quantum Mechanics (Pergamon, New York, 1965), 2nd ed., pg 276.
3. Yamabe, Tokio; Tachibana, Akitomo; Silverstone, Harris J. (1977-09-01). "Theory of the ionization of the hydrogen atom by an external electrostatic field". Physical Review A. 16 (3): 877–890. doi:10.1103/PhysRevA.16.877.
4. Keldysh L V 1965 Soviet Phys. JETP 2354
5. Wolkow, D. M. (1935). "Über eine Klasse von Lösungen der Diracschen Gleichung". Zeitschrift für Physik (in German). Springer Science and Business Media LLC. 94 (3–4): 250–260. doi:10.1007/bf01331022. ISSN 1434-6001. S2CID 123046147.
6. Perelemov A M, Popov V S and Terent'ev M V 1966 SovietPhys. JETP, 23 924
7. Perelemov A M and Popov V S 1967 Soviet Phys.JETP, 25 336
8. Ammosov M V, Delone N B and Krainov V P 1986 SovietPhys. JETP, 64 1191
9. Chin, S L; Yergeau, F; Lavigne, P (1985-04-28). "Tunnel ionisation of Xe in an ultra-intense CO₂ laser field (10¹⁴ W cm⁻²) with multiple charge creation". Journal of Physics B: Atomic and Molecular Physics. IOP Publishing. 18 (8): L213–L215. doi:10.1088/0022-3700/18/8/001. ISSN 0022-3700.
10. Augst, S.; Meyerhofer, D. D.; Strickland, D.; Chint, S. L. (1991-04-01). "Laser ionization of noble gases by Coulomb-barrier suppression". Journal of the Optical Society of America B. The Optical Society. 8 (4): 858. doi:10.1364/josab.8.000858. ISSN 0740-3224.
11. Auguste, T; Monot, P; Lompre, L A; Mainfray, G; Manus, C (1992-10-28). "Multiply charged ions produced in noble gases by a 1 ps laser pulse at lambda =1053 nm". Journal of Physics B: Atomic, Molecular and Optical Physics. IOP Publishing. 25 (20): 4181–4194. doi:10.1088/0953-4075/25/20/015. ISSN 0953-4075.
12. Larochelle, S; Talebpour, A; Chin, S L (1998-03-28). "Non-sequential multiple ionization of rare gas atoms in a Ti:Sapphire laser field". Journal of Physics B: Atomic, Molecular and Optical Physics. IOP Publishing. 31 (6): 1201–1214. doi:10.1088/0953-4075/31/6/008. ISSN 0953-4075.
13. CHIN, S. L. (2004). "From multiphoton to tunnel ionization". Advances in Multi-Photon Processes and Spectroscopy. 16. WORLD SCIENTIFIC. pp. 249–271. doi:10.1142/9789812796585_0003. ISBN 978-981-256-031-5. ISSN 0218-0227.
14. Walsh, T D G; Decker, J E; Chin, S L (1993-02-28). "Tunnel ionization of simple molecules by an intense CO2laser". Journal of Physics B: Atomic, Molecular and Optical Physics. IOP Publishing. 26 (4): L85–L90. doi:10.1088/0953-4075/26/4/002. ISSN 0953-4075.
15. Talebpour, A; Larochelle, S; Chin, S L (1998-01-28). "Suppressed tunnelling ionization of the molecule in an intense Ti:sapphire laser pulse". Journal of Physics B: Atomic, Molecular and Optical Physics. IOP Publishing. 31 (2): L49–L58. doi:10.1088/0953-4075/31/2/003. ISSN 0953-4075.
16. Talebpour, A.; Yang, J.; Chin, S.L. (1999). "Semi-empirical model for the rate of tunnel ionization of N₂ and O₂ molecule in an intense Ti:sapphire laser pulse". Optics Communications. Elsevier BV. 163 (1–3): 29–32. doi:10.1016/s0030-4018(99)00113-3. ISSN 0030-4018.
17. Talebpour, A; Larochelle, S; Chin, S L (1998-06-28). "Multiphoton ionization of unsaturated hydrocarbons". Journal of Physics B: Atomic, Molecular and Optical Physics. IOP Publishing. 31 (12): 2769–2776. doi:10.1088/0953-4075/31/12/012. ISSN 0953-4075.
18. Jaron-Becker, A. (2012). "Molecular Dynamics in Strong Laser Fields". IEEE Journal of Selected Topics in Quantum Electronics. Institute of Electrical and Electronics Engineers (IEEE). 18 (1): 105–112. doi:10.1109/jstqe.2011.2108271. ISSN 1077-260X. S2CID 16703524.
19. Büttiker, M.; Landauer, R. (1982-12-06). "Traversal Time for Tunneling". Physical Review Letters. American Physical Society (APS). 49 (23): 1739–1742. doi:10.1103/physrevlett.49.1739. ISSN 0031-9007.
20. Landsman, Alexandra S.; Weger, Matthias; Maurer, Jochen; Boge, Robert; Ludwig, André; et al. (2014-11-14). "Ultrafast resolution of tunneling delay time". Optica. The Optical Society. 1 (5): 343. arXiv:1301.2766. doi:10.1364/optica.1.000343. ISSN 2334-2536.
21. Fertig, H. A. (1990-11-05). "Traversal-Time Distribution and the Uncertainty Principle in Quantum Tunneling". Physical Review Letters. American Physical Society (APS). 65 (19): 2321–2324. doi:10.1103/physrevlett.65.2321. ISSN 0031-9007. PMID 10042518.
22. Yamada, Norifumi (2004-10-18). "Unified Derivation of Tunneling Times from Decoherence Functionals". Physical Review Letters. American Physical Society (APS). 93 (17): 170401. doi:10.1103/physrevlett.93.170401. ISSN 0031-9007. PMID 15525052.