Relativistic Tunneling Ionization in Crossed Petawatt Laser Beams

A scheme is proposed for a two-beam experiment on multi-petawatt laser systems aimed at searching for relativistic effects in electron tunneling from deep levels of multiply charged ions. It is shown that the most efficient relativistic ionization will occur in the field of a standing wave that arises in the case of counter-propagating beams. Estimates are given showing that the observation of relativistic effects during the ionization of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s$$\end{document}-states of hydrogen-like ions with charge \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z = 80{-} 90$$\end{document} is potentially possible at laser intensity in the beam overlap domain \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \simeq {\kern 1pt} $$\end{document}1025 W/cm2.

Multiphoton and tunneling ionization of atoms and molecules in intense laser fields has been observed experimentally since the mid-1960s [1]. Shortly after the discovery of multiphoton absorption of highpower laser radiation, an essentially nonlinear interaction regime was observed, in which the number of photons involved in each ionization event exceeds the minimum required for the transition to the continuum [2]. The theory of nonlinear ionization in the field of a plane electromagnetic wave was formulated in [3][4][5][6], where analytical formulas for the ionization rate and momentum distributions of photoelectrons were derived within the strong field approximation, when the Coulomb interaction of the emitted electron with the atomic residual is not accounted for. The Coulomb factor in the ionization rate was found in [7] in the limit of a slowly varying field, and then in [8] for an arbitrary laser frequency. The results obtained in the pioneering works [3][4][5][6][7] were significantly advanced in the subsequent decades. At present, the theory of nonlinear ionization of atoms and molecules is widely used to describe the effects of the interaction of highpower laser radiation with matter. The current state of the art in the theory of multiphoton and tunneling ionization of atomic systems in strong fields is presented in reviews [9][10][11][12][13][14].
As the maximal intensity of laser radiation increases, the study of the ionization of deep levels of heavy atoms and multiply charged ions becomes topical. Such processes are interesting in particular in view of their possible use for direct measurement of the laser intensity at the focus (see [15][16][17] and references therein). Besides, the possibility of observing relativistic effects in ionization is of fundamental interest. The photoelectron motion in a laser field becomes relativistic at: (1) Here, and ω are the amplitude and frequency of the laser wave, m and e are the electron mass and the absolute value of its charge, respectively; c is the speed of light. In the following, atomic units are used with . For infrared laser fields with wavelengths μm and , this relativistic threshold is reached at W/cm 2 . Kinematic effects in the spectra of photoelectrons and of radiation emitted by these electrons are well studied (see, for example, [18][19][20] and references therein). At intensities W/cm 2 currently available at many laser facilities, photoelectrons quickly become relativistic after ionization, although the process of their detachment from an atom remains determined by the laws of non-relativistic quantum mechanics, and ionization rates can be calculated with high accuracy using formulas from [3,[5][6][7]. This extension of the applicability of the nonrelativistic theory into the domain of ultrarelativistic intensities is mainly due to the fact that the binding energies of the levels ionized in such fields remain small compared to the rest energy of the electron. As a result, relativistic corrections to the sub-barrier electron trajectory and the respective corrections to the ionization probability appear minor even at intensities W/cm 2 .
At intensities W/cm 2 , currently available or expected to achieve in the coming years with new multi-petawatt laser systems [21][22][23][24], the probability of the electron detachment from an atom is determined not only by the electric but also by the magnetic field of the laser wave, and relativistic corrections to the electron trajectory become significant during the process of transition from the bound state to the continuum. Qualitatively, the regime of this transition is determined by the value of the Keldysh parameter [3]: At , ionization can be considered as multiphoton absorption, while in the opposite limit it is described as tunneling through a slowly varying potential barrier. As the laser intensity increases, the effective value of the Keldysh parameter , upon reaching which the level is rapidly ionized, decreases as . For lasers with a wavelength of 1 μm and atoms in the ground state, the tunneling regime is reached starting from about 10 14 W/cm 2 . Thus, at record intensities >10 23 W/cm 2 considered in this work, ionization proceeds in a deep tunneling regime, making it possible to use the constant field approximation to calculate the ionization probability with a very high accuracy. Presently, for the first time, conditions arise for the experimental observation of the relativistic tunneling effect in atoms.
The theory of ionization of atoms in fields of relativistic intensity determined by the condition (1) was first formulated in [4]. Relativistic tunneling effect was theoretically studied in [25,26], as well as in a number of other later works (see reviews [13,27] and references therein). These papers show that the probability of ionization of a deep-lying atomic level, whose energy is determined taking into account relativistic corrections, can differ significantly, by several orders of magnitude, from that calculated by non-relativistic formulas [5][6][7]. The ratio of ionization rates calculated by non-relativistic and relativistic formulas can be both greater or less than unity [9]. In a recent paper [28], this point was investigated in more detail, and it was shown that the effect of a magnetic field on the probability of tunneling in heavy ions, associated with the Zeeman splitting of levels, can be detected at intensities W/cm 2 . To find effects associated with relativistic corrections to the sub-barrier motion of an electron, even higher intensities exceeding 10 24 W/cm 2 are needed. Achieving such intensities will become possible, in particular, by superimposing several multi-petawatt laser beams in the focal region. Such a scheme for obtaining extreme light fields is proposed at the XCELS [21] facility, which is currently under construction. The maximal number of beams is 12 at a power of 15 PW in each beam. In the electromagnetic field that arises when the laser beams 10 10 × 23 > 2 10 J cross, the electric E and magnetic H fields are no longer equal in amplitude and can also be non-orthogonal. Expressions for the probability of the relativistic tunneling effect in a static electromagnetic field of various configurations were obtained in [25,26] and can be used, after necessary generalizations, to analyze possible multibeam experiments. In this paper, we present such an analysis for the case of ionization of the s-states of multiply charged ions in the field of two crossed laser beams, i.e., the simplest multibeam configuration that can be implemented at the XCELS facility. It is shown that for a fixed power of each of the beams, the most significant relativistic effect can be achieved using the counter-propagating geometry.
Let us consider the field resulting from a coherent overlap of two laser beams linearly polarized in the same direction. Describing each of the beams in the plane wave approximation with the same amplitude values of electric and magnetic fields , , where are unit vectors specifying the directions of wave propagation, and are the wave vectors, we obtain that the maximum value of the electric field is . In this case, the magnetic field amplitude is . Here θ is the half angle between the directions of the wave vectors. The electric and magnetic fields remain orthogonal. Taking into account the highly nonlinear dependence of the ionization rate on the electric field strength, we examine this quantity in those parts of the focus that will make a dominant contribution to the formation of multiply charged ions. These are parts of space where the amplitude of the electric field is maximum.
Consider tunneling from an s-state of a multiply charged ion with an ionization potential . In this case, there are no linear Zeeman shifts in the magnetic field, and the spin effects remain insignificant even in extremely strong fields [29]. In the semiclassical approximation [9,25,26], the ionization rate can be presented in the form (3) where the classical action is calculated along the trajectory corresponding to the zero electron velocity at the barrier exit at : Here, ϕ and A are the scalar and vector potentials specifying the electromagnetic field of the laser wave, and the complex time instant of the motion start is determined from the condition that the initial energy of the photoelectron is equal to that of the bound state.
When the condition is satisfied, the time of subbarrier motion is a small fraction of the optical period, , and the field can be treated as constant; in this case , . The electron trajectory satisfies the equations of motion in complex time with the Lorentz force, and is determined from the equation: Omitting the details of calculations, part of them are contained in [25,26], we present a result applicable at : where is the nonrelativistic ionization rate [5,6,9]. For s-states it is given by (7) The reduced electric field and the effective quantum number are equal to (8) Z is the atomic residual charge. The asymptotic coefficient of the single-electron wave function at large distances from the atomic core can be determined numerically or using a semiclassical asymptotic [9]. For , time in the formulas (6), (7) can be considered as a parameter. Expressions (6), (7) for the ionization rate are obtained in the semiclassical approximation. Their applicability conditions require that the barrier width remains large compared to the electron de Broglie wavelength . This means that the value of the reduced field (8) is . In this case, taking into account the fact that the ionization of deep albeit weakly relativistic levels, , is considered, the field strength upon reaching which the ionization rate (6) becomes significant, grows as . Thus, for the state of hydrogen-like xenon with I p = 41300 eV ≈ 1518 a.u. the value of the reduced field, which provides almost complete ionization in a single period of the laser wave, , is reached at 10 3 a.u. corresponding to the intensity 10 23 W/cm 2 . Note that for s-levels, where , the above-barrier ionization regime [9] is typically not reached [30], and the above given semiclassical formulas (6), (7) remain quantitatively accurate. In the case of ionization of p-states with it may be necessary to take into ν ≈ 2 account corrections arising due to the suppression of the potential barrier (see [30] and references therein).
We will be interested to analyze the difference between the relativistic and nonrelativistic ionization rates determined by the exponential factor in (6). In particular, the dependence of this exponential factor on the angle between the beams will be examined. In [9,25,26] a semiclassical expression was obtained for the probability of relativistic tunneling for the case of perpendicular and unequal fields under the condition , which corresponds to Lorentz ionization [31]. Formulas for the opposite case we are interested in can be obtained from [26] by the method of analytic continuation. A detailed derivation and analysis of a more general case on nonorthogonal fields will be presented elsewhere. The polynomial in the exponent of (6) gives the first term of the expansion of the exact expression for the imaginary part of the action in a small parameter for the case of orthogonal fields under the condition . This formula is sufficient for making estimates that allow one to establish the conditions necessary for experimental observation of relativistic effects in laser fields of extreme intensity.
In order to formulate these conditions, note the following.
1. The probability of tunneling ionization is a sharp function of the reduced field F, therefore, for a given electron level with an ionization potential , tunneling effectively occurs in a narrow range of fields near [15]. The value of this interval depends logarithmically on the laser pulse duration, which makes it possible to provide a fairly universal estimate of the field threshold for observing tunneling ionization. At , the complete ionization of levels with any values of the azimuthal and magnetic quantum numbers will occur within one optical period (or even faster) of a laser pulse with a wavelength of 1 μm [15]. For laser pulses with duration of fs, which corresponds to the planned parameters of the XCELS laser facility, a gas target density of 10 15 cm -3 and a focus volume of , several tens of ions (a number sufficient for registration of multiply charged ions by modern detectors) will be produced at . When going to , the ionization rate will drop by 5 orders of magnitude so that the probability of generating even one ion per pulse becomes negligible. Thus, the value can be taken as the lower threshold for the observation of tunneling ionization.
2. The restriction on the reduced field F specifies the maximum value of the binding energy of the level, whose ionization can, in principle, be observed at a certain maximum value of the laser intensity at the focus. Intensities of ~10 23 W/cm 2 [24] have been achieved so far. With optimistic expectations, intensi- = 0.02 F = 0.015 F = 0.02 F ties of ~10 25 W/cm 2 can be achieved at 100-petawatt class facilities XCELS in Russia and Station of Extreme Light in China. According to the above estimates, in the field with such an intensity it is possible to observe tunneling from levels with ionization potentials of about 120 keV, which corresponds to the ground states of hydrogen-like radon or thorium ions. In this case so that the expansion of the action in this parameter quickly converges, and the approximate expression (6) can be used for estimates.
In the case of two parallel beams, the coefficient in front of in (6) is , for (beams propagating orthogonally) , and for counterpropagating beams -. As a result, in the transition from ionization in two parallel or nearly collinear beams to the case of counter-propagating beams, the ionization rate and the yield of fully ionized atoms with binding energies of states of about 100 keV should increase by approximately a factor of 2. The increase in the ionization probability is mainly due to relativistic corrections to the sub-barrier electron trajectory. Thus, an experimental measurement of the increase in the yield of multiply charged ions in the situation described above could be interpreted as an observation of the influence of relativistic effects on tunneling in a strong external field. It should be noted that the laser field intensity remains also ultrarelativistic for , that is, for parallel beam propagation or in the case of a single beam. The relativistic correction to the tunneling probability in this case appears smaller, because the magnetic field deforms the subbarrier electron trajectory, making it two-dimensional, which leads to an increase in the imaginary part of the action (4) and to the suppression of the ionization probability [26]. Relativistic corrections to the action due to the electric field, on the contrary, increase the probability. In the case of a crossed field, these corrections are almost completely compensated, and the net effect appears strongly suppressed.
As one can see from the above estimates, the relativistic effects in tunneling turn out to be relatively weak even at high intensities [15,28]. A relativistic enhancement of the ionization probability in the field of a standing wave by an order of magnitude compared to the ionization probability in parallel beams can be expected at intensities W/cm 2 , which is hardly possible to achieve in the coming years.
Let us discuss a possible experimental realization. The ion yield must be measured in a rarefied gas jet at concentrations <10 15 cm -3 , which allow one to avoid the influence of plasma fields, the appearance of bremsstrahlung, and collisional ionization. The requirements to an experimental setup aimed to observe tunneling ionization of deep atomic levels in a strong laser field are described in detail in [32,33]. The measured value is the number of ions of the max- imum charge state depending on the mutual orientation of the beams at fixed energies of the latter. A significant problem for extracting the contribution of relativistic ionization to the yield of ions of the maximum charge state will be the change in the structure of the focus with variation in the angle between the beams. As the angle θ increases, the part of space where the electric field amplitude takes values close to the maximum will decrease. This effect may happen to be considerably more significant than the increase in the number of ions due to the suppression of the magnetic field in the counter-propagating beams. The influence of a change in the effective focal volume can be eliminated by measuring the yield of ions of smaller charge at lower intensities, for which relativistic effects are clearly insignificant, and the ionization probability is determined by the amplitude of the electric field only. Thus, an experiment in which the angle between the beams changes from a small value to (orthogonal beams) and then to (counter-propagating beams) must be carried out at two values of laser power that ensure the achievement of an intensity of ≈10 25 W/cm 2 , at which hydrogen-like ions Rn and Th are ionized and, for example, at W/cm 2 when the shells of the same atoms are ionized. In the second case, the measured ratio of the number of ions of charge for different angles between the beams, , etc., will then make it possible to calibrate the measurement results for and to extract the contribution of relativistic effects.
The proposed experiment can also be performed in a scheme with the number of beams >2. Considering that the maximum ionization rate corresponds to the case of a zero magnetic field, such a scheme in its most preferable form should use an even number of mutually counter-propagating beams. Note that the proposed experimental scheme is feasible only at a high degree of coherence of the overlapping beams. It is not yet clear what degree of coherence can be achieved on multi-beam laser systems with multi-petawatt power. N format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.