Unbibium

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; 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.

DC Tunneling Ionization

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

$w=4{\omega }_a\frac{E_a}{\left|E\right|}\exp [-\frac{2}{3}\frac{E_a}{\left|E\right|}]$

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

$E_a=\frac{m^2{e}^5}{(4\pi {ϵ}_0{)}^3{\hslash }^4}$,

${\omega }_a=\frac{m{e}^4}{(4\pi {ϵ}_0{)}^2{\hslash }^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 an 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 fields oscillation.                                                 Multiphoton and tunnelionization  of an atom or a moleculedescribe 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 findingthe time evolution of the state of the atom which is described by the Schrödinger equation.

File:Tunnel ionization 3.png

Combined potential of an atom and a uniform laser field. At distances r<r0, the potential of the laser can be neglected, while at distances with r>r0 the Coulomb potential is negligible compared to the potential of the laser field. The electron emerges from under the barrier at r=Rc. Ei is the ionization potential of the atom.

When the intensity of 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. He modeled the MPIprocess 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 it is observed from 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. included the Coulomb interaction at larger internuclear distances. Their model (which we call 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. 

On the experimental side, many experiments have been carried out on the MPI of rare gas atoms using strong laser pulses, both through measuring the total ion yield and the kinetic energy of the electrons. Here, we only consider the experiments designed to measure the total ion yield. Among these experiments,we highlight those by  Chin et al.,   Augst et al. and Auguste et al. Chin et al.  used  a10.6 mm 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 near IR or visible region of frequencies. These findings  weakened the suspicion on the applicability of models basically founded on the assumption of a structurless atom. 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 fit 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 findingthe time evolution of the stateof the atomwhich 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

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

where $\mathbf{E}(t)$ is the electric field of the laser and 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 $\sqrt{2E_i}.\delta (\mathbf{r})$ ( $E_i$ the magnitude of the ionization potential of the atom), the probability current $\mathbf{J}(\mathbf{r},t)$ is calculated. Then, the total MPI rate from short range potential for linear polarization, $W(\mathbf{E},\omega )$, is found from

$W(\mathbf{E},\omega )=li{m}_{x->\mathrm{\infty }}\underset{0}{\overset{\frac{2\pi }{\omega }}{\int }}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\mathbf{J}(\mathbf{r},t)dzdydt$

where $\omega$ is the frequency of the laser, which is assumed to be polarized in the direction of the x axis. The effect of the ionic potential, which behaves like $\frac{Z}{r}$ (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

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

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

$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{)}^{n^{*2}-\left|m\right|-3/2}(1+\gamma {)}^2{)}^{|m/2|+3/4}{A}_m(\omega ,\gamma )e^{-(2(2E_i{)}^{\frac{3}{2}}/F)g(\gamma )}$

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

$f_{lm}=\frac{(2l+1)(l+|m|{)}^{!}}{2^m|m{|}^{!}(l-\left|m\right|{)}^{!}}$

$g(\gamma )=\frac{3}{2\gamma }(1+\frac{1}{2{\gamma }^2}sin{h}^{-1}(\gamma )-\frac{\sqrt{1+{\gamma }^2}}{2\gamma })$

$|C_{n^{*}{l}^{*}}{|}^2=\frac{2^{2n^{*}}}{n^{*}\mathrm{\Gamma }(n^{*}+l^{*}+1)\mathrm{\Gamma }(n^{*}{l}^{*})}$

The coefficient $A_m(\omega ,\gamma )$ is given by

$A_m(\omega ,\gamma )=\frac{4}{3\pi }\frac{1}{|m{|}^{!}}\frac{{\gamma }^2}{1+{\gamma }^2}\sum _{n>v}^{\mathrm{\infty }}{w}_m(\sqrt{\frac{2\gamma }{\sqrt{1+{\gamma }^2}}(n-v)}{e}^{-(n-v)\alpha (\gamma )})$

where

$w_m(x)=e^{-x^2}\underset{0}{\overset{x}{\int }}(x^2-y^2{)}^m{e}^{y^2}dy$

$\alpha (\gamma )=2(sin{h}^{-1}(\gamma )-\frac{\gamma }{\sqrt{1+{\gamma }^2}})$

$v=\frac{E_i}{\omega }(1+\frac{2}{{\gamma }^2})$

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

$W_{ADK}=|C_{n^{*}{l}^{*}}{|}^2\sqrt{\frac{6}{\pi }}{f}_{lm}{E}_i(2(2E_i{)}^{\frac{3}{2}}/F{)}^{n^{*2}-\left|m\right|-3/2}{e}^{-(2(2E_i{)}^{\frac{3}{2}}/3F)}$

In practice, the limit for the QST regime is $\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 $\gamma$, since the potential is bent down during half a cycle of the field oscillation and the ratio can be expressed as

$\gamma =\frac{{\tau }_T}{\frac{1}{2}{\tau }_L}$

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

The 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 10.6 mm 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. were able to quantitatively fit the ionization yield of  diatomic molecules interacting with a Ti:Sappire laser pulse. The conclusion of our work was that the MPI rate of a diatomic molecule can be predicted from the PPTmodel  by assuming that the electrontunnels through a barrier given by $\frac{Z_{eff}}{r}$ instead of  barrier $\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, $Z_{eff}$. Using the semi-empirical model for the MPI rate of unsaturated hydrocarbons is feasible. Of course, 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.

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.