Coulomb-Corrected Strong Field Approximation without Singularities and Branch Points
The domain of validity of the Coulomb-Corrected Strong Field Approximation (CCSFA) is going to be analyzed in relation to the semi-classical dynamics of electrons during ionization of hydrogen-like targets. Our analysis is limited to ionization driven by Ti-sapphire laser pulses with intensities up to roughly Wcm. For such parameters, the effects related to radiation pressure are small and the laser field can be described in the dipole approximation. By applying the Magnus expansion for the exact retarded electron propagator we obtain an effective action which is free of Coulomb singularities and branch points when the complex-time trajectories are used. Furthermore, we show that the classical action is exactly recovered as the asymptotic limit of its effective counterpart. The applicability of such limit is also discussed.
The great English mathematician and philosopher Bertrand Russell (1872-1970) has once written that ’all exact science is dominated by the idea of approximation.’ One could continue his line of thought by saying that the most beautiful aspects of all exact science are the approximations. It has to be remembered, however, that any approximation has its domain of validity, which is inherently related to it and equally important in its applications. The strong-field physics (SFP) is, in fact, one of the best illustrations of the above statements. Recent developments in this field have led to the creation of new branch of modern science, the attophysics BK (); IK (). One of the most prominent tools in the SFP, the Strong-Field Approximation (SFA) in photoionization, was originally introduced by Keldysh Keldysh () in the length gauge and further developed by Faisal Faisal () and Reiss Reiss () for other forms of the Schrödinger equation (i.e., in the Kramers-Henneberger frame kh1 (); kh2 () or in the velocity gauge, respectively). The common feature of those approaches is the approximation of the exact scattering state of the photoelectron by the Volkov solution Volkov (); EJK1998 (). Subsequently, such approximation has been extended to treat the non-relativistic laser-assisted scattering processes Fedorov1 () and its relativistic counterpart Fedorov2 (), where the relativistic form of the Volkov solution EKK2009 (); DiPiazza2012 (); DiPiazza2018 () is fully exploited.
Since the SFA does not offer a proper explanation of various experimental results, it requires further developments. For instance, an attempt to incorporate the interaction of photoelectrons with the parent ion by means of the Coulomb-Volkov state was undertaken in Tsoar (); Kornev (); Duchateau (). This procedure was successfully applied to the analysis of ionization driven by elliptically polarized laser fields Ferrante (). (Note that similar investigations were recently presented for bi-circular laser fields Milo2018 ().) Further generalizations, which fully account for the low-frequency approximation, were considered in Refs. K1 () and K2 () for the scattering and ionization processes, respectively. This approach, called the Coulomb-Volkov Strong-Field Approximation (CVSFA), was recently generalized in PU2018 () by incorporating the density functional theory for the initial bound state of many-electron atoms, or in ZZ2018 () by applying the parabolic quasi-Sturmian-Floquet approach. Also, the scattering states in the high-frequency approximation were explored Gav1 (); Gav2 (). An alternative approach to the CVSFA, which we call the Born-Series Strong-Field Approximation (BSSFA), is to apply the Born expansion to the final electron scattering state in both the binding potential and the laser field. Here, the first two terms of the Born series were successfully used in studies of the re-scattering process in ionization Milo2018 (); Resc1 ().
Another path of theoretical explorations is related to the eikonal approximation eik1 (); eik2 (); eik3 () and its generalizations eik4 (); KK2015 (); VKK2016 (). In our further discussion, we shall refer to those approaches as the Coulomb-Corrected Strong-Field Approximation (CCSFA) ccsfa1 (); ccsfa2 (); ccsfa3 (); ccsfa4 (). In connection to these theoretical investigations is the purely classical analysis of laser-assisted atomic or molecular processes classdyn () such as ionization classion1 (); classion2 (), plasma dynamics in atomic clusters classion2a (), and re-collisions classresc (). In fact, those works were stimulated by the successful application of the classical dynamics to ionization of Rydberg states by microwave fields classrf1 (); classrf2 ().
The two groups of approaches mentioned above (the CVSFA and BSSFA vs. the CCSFA and the classical analysis) are contradictory when the behavior of the electron wave-packet close to the Coulomb singularity is analyzed. For the first group of approximations (and actually for the full quantum treatment; see, e.g., Nurhuda (); num1 ()) this singularity does not present any conceptual or numerical problem. In contrast, the classical analysis shows severe numerical problems when the electron trajectories approach the singularity (as discussed in classion2 ()). Furthermore, in the CCSFA, when the method of complex-time electron trajectories is applied, the Coulomb singularity leads to singular branch points. Such branch points do not exist in the BSSFA or in the exact numerical solution of the Schrödinger equation. The aim of our investigations is to determine the origin of these problems and to derive the effective CCSFA (related to the generalized eikonal approximation and eikonal perturbation theory eik4 (); KK2015 (); VKK2016 ()) which is free of Coulomb singularities and branch points.
Note that the CCSFA and the classical dynamics follow from the quantum theory in the limit when the Planck constant goes to 0 or the time evolution is very short. However, it must be taken into account that and time have physical dimensions. Thus, in order to provide a sensible physical meaning of this limit, one has to construct a dimensionless parameter out of the Plank constant, the time of evolution, and other relevant physical quantities, such that the CCSFA or the classical dynamics are recovered in such limit. This is going to be done below for the ionization of a hydrogen-like ion driven by laser fields of moderate intensity. Such condition guarantees that the dipole approximation is valid throughout our analysis. For instance, if the Ti-sapphire laser field (wavelength 800 nm and frequency such that ) is considered, the laser pulse should not exceed intensities of the order of . For larger intensities, relativistic effects related to the radiation pressure press0 (); KKpress () (see, also press1 (); press2 (); press3 (); press4 (); press5 (); press6 ()) and relativistic mass corrections CKKvortex () become already visible.
Our analysis below is based on the Magnus expansion, as opposed to the more general eikonal perturbation theory developed in eik4 (); KK2015 (); VKK2016 (). The Magnus expansion magnus1 (); magnus2 (); magnus2a (); magnus3 () allows one to construct an approximate exponential representation of the propagator of the system. Every order of such expansion corresponds to a partial re-summation of infinite terms of an ordinary Born series. Next, we apply this approximation up to the linear term with respect to the binding potential and show under which conditions the classical dynamics is restored. This determines the domain of validity of the ordinary CCSFA ccsfa1 (); ccsfa2 (); ccsfa3 (); ccsfa4 () and it allows us to remove the artificial Coulomb singularities and branch points. We show that the latter arise from the application of the relevant asymptotic expansion beyond its domain of validity. We demonstrate that the effective CCSFA, in the appropriate limits, leads to the exact Born expansion or to the classical (or complex-time) dynamics. Moreover, it appears that the effective action is complex and accounts for the electron wave packet spreading in the course of its time evolution. Due to this fact the electron trajectory is also complex even for real times and real initial/boundary conditions, but such that in the classical limit its imaginary part vanishes. Note that the method presented here and in eik4 (); KK2015 (); VKK2016 () allows to incorporate further terms, which are nonlinear with respect to the static binding potential, into the effective action.
Ii General theory and Magnus expansion
We start with general statements about the non-relativistic quantum-mechanical evolution of an electron in arbitrary electromagnetic fields. Such evolution, from an initial time up to a final time , is determined by the retarded propagator . If the Hamiltonian of the system, in terms of the ‘primed’ coordinates, is , the propagator satisfies the differential equation (see, e.g., Ref. KK2015 ())
As many studies of ionization relate to the momentum or energy distribution of photoelectrons, it is most convenient to analyze the propagator in momentum space. For this reason, we introduce its Fourier transform, , calculated over the non-primed coordinate,
and, as it can be checked, fulfills the retardation and initial conditions,
As mentioned in magnus2 (), the Magnus expansion finds applications in many subfields of physics including atomic, molecular, and particle physics, quantum electrodynamics, etc. Here, we shall use it to solve the homogeneous equation,
where is an a priori unknown function. Such function can be expanded in the power series
where is a small and real parameter (see below). The challenge, however, is to determine each term in the sum (7). Below, we illustrate this for the case when can also be represented as a Born series.
First, we separate the total Hamiltonian of the system into two parts, , where describes a relatively weak interaction compared to the main contributing Hamiltonian . This is stressed by a small parameter . Next, we also assume that the transformed propagator can be represented as a Born series with respect to . Specifically, up to the first order in the perturbation, we write it down as
Therefore, up to the first order in , we obtain
By comparing Eqs. (8) and (10), we conclude that and . Hence, by performing the Born and Magnus expansions of the transformed propagator simultaneously, the terms can be determined. However, it has to be noted that the complexity of the procedure increases with , as for higher orders in a larger number of factors, arising from the series expansion in (9), contribute to the particular term.
Up to now our analysis has been very general. In the following, the Hamiltonian will describe an electron interacting with an oscillating laser field, which in the dipole approximation is defined by the vector potential , and with a general time-dependent binding potential, , in the Kramers-Henneberger frame.
Iii Magnus and Born expansions in the Kramers-Henneberger frame
Our aim is to derive an approximate expression for the propagator in the Kramers-Henneberger (KH) frame. Such frame is particularly suitable to perform our calculations, as the interaction with the scalar and vector potentials is reduced to a single term in the Hamiltonian.
where and is the free particle Hamiltonian. is the so-called displacement vector, and it relates to the oscillating electric field such that . This, in turn, suggests that determines the classical trajectory of a free electron in the light field, and corresponds to its acceleration kh2 (). In this paper, it is assumed that the laser pulse acts over a finite period of time (from up to ), which means that for and . We also choose the vector potential such that , and assume that the displacement vanishes at times .
We proceed with expanding the propagator in a Born series. This is done by using repeatedly the Lippmann-Schwinger equation,
where is the free-particle propagator. In the first-order Born approximation, we obtain that
Note that such general iterative procedure leads to a series expansion of the total propagator in powers of the binding potential. This is the essence of the BSSFA.
Since we are interested in deriving an expression for the propagator in momentum space, we perform the Fourier transform of in (13) with respect to . In addition, since our main focus is the photoionization of atoms/ions, the binding potential is static. This means that and we can use
As a result,
where one can derive that . Moreover, in Eq. (15), we have introduced a classical trajectory of an electron in the laser field, . Note that the second term in the square bracket in Eq. (15) can be related to an effective potential experienced by the electron,
Having this in mind, we rewrite (15) in a more compact form,
which represents the Magnus expansion of the propagator in momentum space (10) with
It allows us to formally write down that, up to the linear term in the potential,
This expression is fundamental for our further analysis. As we will demonstrate shortly, the fact that the propagator depends on an effective potential rather than on the classical one, has very important consequences in the photoionization dynamics.
iii.1 Effective CCSFA
In order to define our effective CCSFA, let us modify the classical trajectory of the electron in the laser field by a correction that is linear in . Namely, we introduce an effective trajectory which satisfies the same boundary conditions as the classical trajecotry; namely, and . By making a transformation from the KH to the laboratory frame, we obtain that, up to linear terms in , the propagator (19) is given by
where the subscript stands for the velocity gauge and where , for fixed and , is a function of and . The functional in (21) is given by
which can be recognized as the classical action of an electron interacting with the laser field and the scalar effective potential . This, in turn, implies that the dynamics of the system is governed by the effective (and complex) Newton equation
, being a complex function, differs from the classical potential as it contains important quantum corrections (see, Sec. IV). Furthermore, the effective potential (16) depends explicitly on the initial time , meaning that the Newton equation (23) contains information about past events. In other words, contrary to its classical counterpart, the effective quantum-mechanical evolution is non-Markovian.
Now, we present for the spherically-symmetric Coulomb potential . Here, is the fine-structure constant, is the atomic number, and is the speed of light. From the definition (16), it follows that (see, Refs. KK2015 (); VKK2016 ())
where . The mathematical and physical consequences of this expression are analyzed below.
Iv Classical limit of the effective CCSFA and its range of validity
The propagator in the effective CCSFA [Eq. (19)] depends on the effective potential (24), which differs from the pure Coulomb one. As a result, the effective interaction already accounts for the spreading of the electron wave packet and quantum diffusion effects KK2015 (); VKK2016 (). Such effects, while negligible at infinitesimally-short time intervals, become relevant at the characteristic time-scale at which photoionization takes place (see below). In addition, is an entire function of , i.e., it does not have singularities or branch points when the electron trajectories come back to the parent ion. This is in contrast to the pure Coulomb potential, as it is singular at the origin of coordinates. To see this better, consider the power-series expansion of the error function (see, Ref. olver2010 ()) in Eq. (24),
Thus, the effective potential can be expressed as an absolutely convergent power series in the whole space,
and, in particular, , provided that . On the other hand, for large arguments, the effective potential can be analyzed by making use of the asymptotic expansion of the complementary error function, (see, Ref. olver2010 ()). In doing so, we find that for large (or, equivalently, for large ), , i.e., the Coulomb interaction is recovered when the argument of the effective potential is large enough.
In the upper panel of Fig. 1 we show the real (blue line) and imaginary (green line) parts of together with the curve (red line). While the first function determines the behavior of the effective potential [Eq. (24)], the second one defines the classical Coulomb interaction. It can be seen that the imaginary part of goes fast to zero for large values of . In contrast, behaves in a very similar way as , when . This is in agreement with our previous analysis which established that, in the limit of large arguments, the effective potential and the classical one coincide. On the other hand, for small values of the situation changes drastically. First, is nonzero and remains finite. Second, the general behavior of and differ considerably: while the former remains bounded and reaches values near to one, the latter increases rapidly and presents a singularity at . This agrees with our earlier conclusions that the effective potential is an entire function of and does not contain singularities or branch points. The dramatic differences in the behavior of the two potentials near the origin of coordinates () can be attributed to the inherent quantum-mechanical nature of photoionization; diffusion and spreading of the wave packet play an important role in the process and they are neglected in a classical treatment or CCSFA.
As we have shown, only the asymptotic form of (which is valid provided that is large) has singularities and branch points at small values of . A very similar situation is met in the theory of special functions. As an example, consider the Bessel function , which is an entire function of (cf., Fig. 2). Its power series expansion,
is absolutely convergent and well-behaved for all complex arguments . However, its asymptotic behavior, for , might indicate that exhibits the singularity and the branch point at . This is not the case because the asymptotic form is not applicable there. We illustrate this in Fig. 2, where (blue line) and its asymptotic expansion (red line) are presented.
Now, an important question arises: What are the arguments of for which the effective potential can be well approximated by its asymptotic (classical) form? This can happen for real such that . In the lower panel of Fig. 1 we present the real and imaginary parts of this function (blue and green lines, respectively). Hence, one can estimate that the condition is quite well fulfilled for , i.e., the classical form of the effective potential can be applied provided that
This condition is satisfied if either:
For fixed and , the particle mass is very large. In fact, such statement is the essence of the Born-Oppenheimer approximation, i.e., it is assumed that the dynamics of massive particles can be approximated by the classical treatment.
For fixed and , the trajectories do not come back close to the origin of coordinates. This limit is justified, for instance, in the ionization of Rydberg states (RS) by microwave fields. If the distance between the electron and the nucleus is approximately equal to the Bohr radius of the RS (with principal quantum number ), and corresponds to its orbital period, the classical limit is applicable provided that [see, Eq. (28)]. Thus, for the electron dynamics is determined by the classical evolution classrf1 (); classrf2 ().
In scattering processes the condition (28) is equivalent to , where is the electron momentum and is the impact parameter. This is, in fact, the condition of applicability of the traditionally-used eikonal approximation. Hence, such approximation is valid provided that and/or are sufficiently large and it excludes the case of backward scattering.
It follows from here that the asymptotic expansion of the effective potential can be applied, for instance, to ionization driven by elliptically polarized laser fields, as the electron trajectories do not return to the parent ion. The situation is different when a linearly polarized pulse is considered. In such case, the electron can come back to the origin of coordinates at times close to half of the laser cycle, i.e., , where is the carrier frequency of the laser field. This means that the asymptotic limit (i.e., the classical theory or the CCSFA) is only justified provided that , in atomic units. For a typical frequency of eV a.u., we have or, equivalently, a.u. This means that the electron trajectory should be far away from the parent ion. The range of applicability of the classical description of ionization is improved for X-ray laser pulses. For instance, for eV the electron can approach the origin of coordinates up to a distance of the order of , where is the Bohr radius.
Both CCSFA and the classical approach are asymptotic limits of the quantum theory based on the Schrödinger equation. Therefore, they can be used only within their domains of validity. We suggest that, beyond those domains, either the methods based on the Born expansion or the effective CCSFA with the regular complex effective action should be used. Our analysis has been limited to the case when the effective action is a linear functional of the binding potential. However, by considering higher terms in the Born expansion or the eikonal perturbation theory, one can derive further modifications of the effective CCSFA which are consistent with the quantum theory.
This work is supported by the National Science Centre (Poland) under Grant No. 2014/15/B/ST2/02203.
- This work was presented at the 27th Annual International Laser Physics Workshop (Nottingham, July 16-20, 2018).
- Brabec T and Krausz F 2000 Rev. Mod. Phys. 72 545
- Krausz F and Ivanov M 2009 Rev. Mod. Phys. 81 163
- Keldysh L V 1964 Zh. Eksp. Teor. Fiz. 47 1945; 1965 Sov. Phys. JETP 20 1307
- Faisal F H M 1973 J. Phys. B 6 L89
- Reiss H R 1980 Phys. Rev. A 22 1786
- Kramers H A 1956 Collected Scientific Papers (Amsterdam: North-Holland) p 262
- Henneberger W C 1968 Phys. Rev. Lett. 21 838
- Wolkow D V 1935 Z. Phys. 94 250
- Ehlotzky F, Jaroń A and Kamiński J Z 1998 Phys. Rep. 297 63
- Bunkin F V and Fedorov M V 1965 Zh. Eksp. Teor. Fiz. 49 1215; 1966 Sov. Phys. JETP 22 844
- Denisov M M and Fedorov M V 1967 Zh. Eksp. Teor. Fiz. 53 1340; 1968 Sov. Phys. JETP 26 779
- Ehlotzky F, Krajewska K and Kamiński J Z 2009 Rep. Prog. Phys. 72 046401
- Di Piazza A, Müller C, Hatsagortsyan K Z and Keitel C H 2012 Rev. Mod. Phys. 84 1177
- Di Piazza A 2018 Phys. Rev. D 97 056028
- Jain M and Tzoar N 1978 Phys. Rev. A 18 538
- Kornev A S and Zon B A 2002 J. Phys. B: At. Mol. Opt. Phys. 35 2451
- Duchateau G 2018 Eur. Phys. J. Plus 133 186
- Basile S, Trombetta F and Ferrante G 1988 Phys. Rev. Lett. 61 2435
- Gazibegović-Busuladžić A, Becker W and Milošević D B 2018 Opt. Express 26 12684
- Kamiński J Z 1986 Phys. Scr. 34 770
- Kamiński J Z, Jaroń A and Ehlotzky F 1996 Phys. Rev. A 53 1756
- Pyak P E and Udachenko V I 2018 J. Phys. B: At. Mol. Opt. Phys. 51 065001
- Zaytzev A S, Zaytzev S A, Ancarani L U and Kouzakov K A 2018 Phys. Rev. A 97 043417
- Offerhaus M J, Kamiński J Z and Gavrila M 1985 Phys. Lett. A 112 151
- Gavrila M 2002 J. Phys. B: At. Mol. Opt. Phys. 35 R147
- Becker W, Lohr A and Kleber M 1995 Quantum Semiclass. Opt. 7 423
- Choudhury B J and Bakar B S 1974 J. Phys. B 7 L137
- Choudhury B J and Bakar B S 1975 J. Phys. B 8 L228
- Zon B A 1975 J. Phys. B 8 L86
- Kamiński J Z 1984 Acta Phys. Pol. A 66 517
- Cajiao Vélez F, Krajewska K and Kamiński J Z 2015 Phys. Rev. A 91 053417
- Cajiao Vélez F, Krajewska K and Kamiński J Z 2016 J. Phys.: Conf. Ser. 691 012003
- Popruzhenko S V, Mur V D, Popov V S and Bauer D 2008 Phys. Rev. Lett. 101 193003
- Pisanty E and Ivanov M 2016 Phys. Rev. A 93 043408
- Maxwell A S, Al-Jawahiry A, Das T and Figueira de Morisson Faria C 2017 Phys. Rev. A 96 023420
- Shvetsov-Shilovski N I and Lein M 2018 Phys. Rev. A 97 013411
- Grozdanov T, Grujić P and Krstić P 1989 Classical Dynamics in Atomic and Molecular Physics (Singapore: World Scientific)
- Gajda M, Grochmalicki J, Lewenstein M and Rza̧żewski K 1992 Phys. Rev. A 46 1638
- Schmitz H, Boucke K and Kull H J 1998 Phys. Rev. A 57 467
- El-Khawaldeh A and Kull H J 2017 Phys. Rev. A 95 043401
- Mauger F, Chandre C and Uzer T 2010 Phys. Rev. Lett. 105 083002
- Leopold J G and Parsival I C 1978 Phys. Rev. Lett. 41 944
- Leopold J G and Parsival I C 1979 J. Phys. B 12 709
- Nurhuda M 2004 Comp. Phys. Commun. 162 1
- Majorosi S and Czirják A 2016 Comp. Phys. Commun. 208 9
- Bethune-Waddell M and Chau K J 2015 Rep. Prog. Phys. 78 122401
- Krajewska K and Kamiński J Z 2015 Phys. Rev. A 92 043419
- Smeenk C T L, Arissian L, Zhou B, Mysyrowicz A, Villeneuve D M, Staudte A and Corkum P B 2011 Phys. Rev. Lett. 106 193002
- Titi A S and Drake G W F 2012 Phys. Rev. A 85 041404(R)
- Reiss H R 2013 Phys. Rev. A 87 033421
- Chelkowski S, Bandrauk A D and Corkum P B 2014 Phys. Rev. Lett. 113 263005
- Ivanov I A 2015 Phys. Rev. A 91 043410
- He P-L, Lao D and He F 2017 Phys. Rev. Lett. 118 163203
- Cajiao Vélez F, Krajewska K and Kamiński J Z 2018 Phys. Rev. A 97 043421
- Magnus W 1954 Comm. Pure Appl. Math. 7 649
- Mananga E S and Charpentier T 2016 Phys. Rep. 609 1
- Blanes S, Casas F, Oteo J A and Ros J 2009 Phys. Rep. 470 151
- Kuwahara T, Moria T and Saito K 2016 Ann. Phys. 367 96
- Olver F W J, Lozier D W, Boisvert R F and Clark C W 2010 NIST Handbook of Mathematical Functions (New York: Cambridge)
- Dirac P A M 1933 Phys. Z. Sowjetunion 3 64
- Feynman R P 1948 Rev. Mod. Phys. 20 367
- Kleinert H 2006 Path Integrals in Quantum Mechanics, Statistics, Polimer Physics, and Financial Markets (Singapore: World Scientific)