1. Introduction
During the past decades, strong-field ionization measurements in atoms and molecules have led to numerous insights into the electron dynamics on short time scales. In particular, several nonlinear optical processes, such as the above-threshold ionization (ATI, [
1,
2]), tunneling ionization, high-order harmonic generation (HHG, [
3,
4]), or the nonsequential double ionization (NSDI, [
5]) have attracted much interest and can be readily controlled by tailoring the temporal shape and duration of ultrashort laser pulses. In ATI, for example, the energy and momentum distributions of photoelectrons are often recorded for different targets and (short) laser pulses of different frequency
, intensity
I, polarization
, pulse duration (i.e., number of laser cycles,
), or by even steering the carrier–envelope phase
. In contrast to the detailed modeling of the driving laser pulse, however, the target atoms are typically described in rather a simplified manner, and especially the initial state of the
active electron is often just taken as a hydrogenic
state [
6,
7]. Because of this and further simplifications in modeling the target atoms, many observations are still understood only
qualitatively so far.
Figure 1a displays the prototypical geometry and observables of an ATI experiment. Here, atoms are exposed to an intense driving laser pulse with given intensity
I, wavelength
, ellipticity
or, perhaps, even a superposition of such light fields. A detector D records the photoelectrons that are emitted due to the interaction of laser pulse with the target atoms. Routinely, the photoelectron energy distributions are recorded at a fixed detector position within the
polarization plane (
Figure 1b). The observed photoelectron energy spectra then exhibit ATI peaks that are spaced by the photon energy. If the detector position is varied within the polarization plane, azimuthal angular distributions can be recorded for photoelectrons of selected energy (
Figure 1c). These angular distributions strongly depend not only on the shape of the driving laser pulse but also on the outgoing electron wave in the potential of the photoion. In addition, the full photoelectron momentum distributions are often measured within the polarization plane as shown in
Figure 1d.
While any reliable theory of the strong-field (ionization) processes from above must have its roots in the time-dependent Schrödinger equation, a direct (numerical) integration of this equation becomes unfeasible already for the three-dimensional motion of a single (active) electron in a static soft-core potential, not to speak about the many-electron nature of most targets [
8]. Therefore, a number of analytical methods have been developed as well and nowadays provide good insights into the underlying electron dynamics. In particular, the strong-field approximation (SFA) [
9,
10,
11] provides an efficient single-electron treatment and has become, despite several limitations in its original form, a very valuable tool for computing the ATI and HHG spectra for a wide range of laser parameters and targets [
12,
13,
14]. Here, however, a re-formulation of the SFA in terms of partial waves and spherical tensors [
15] is applied and help adopt this method towards modern strong-field measurements. This re-formulation enables one to incorporate all central features of the incident laser pulse as well as the electronic structure of the target atoms.
To support the analysis of different strong-field measurements, this work reports an implementation of the (direct) SFA amplitude in its partial-wave representation within the framework of
Jac, the Jena Atomic Calculator [
16]. This toolbox, which facilitates the (relativistic) computation of atomic structures and processes [
17,
18], has been expanded here in order to model the initial-bound and final-Volkov states in the computation of strong-field amplitudes. Apart from the active-electron waves, however, our implementation below is flexible also in choosing the polarization, shape (envelope) and even the CEP phase of the driving laser pulse. Indeed, all these features have been found to be (very) crucial to further adopt the theoretical modeling of strong-field ionization processes to ongoing experiments.
The paper is structured as follows: After a brief discussion of the SFA and the direct amplitude in terms of partial waves in
Section 2, emphasis is placed on the implementation within the framework of
Jac as well as the role of appropriate data structures for simplifying the communication
with and
within the program.
Section 3 then explains and discusses how the different energy, angular and momentum distributions can be obtained quite readily by just specifying the initial and final levels of the target atom as well as the parameters of the laser pulse. This includes the choice of atomic potential and the Volkov states in the evaluation of the amplitudes. Finally, a short summary and conclusions are given in
Section 4, with emphasis on possible extensions of the code towards rescattering phenomena, or the computation of harmonic spectra.
2. Strong-Field Amplitudes and Probabilities
2.1. Brief Account of the Strong-Field Approximation
The SFA has been known as perhaps the most straight avenue for modeling strong-field ionization processes and for analyzing most of the associated ATI spectra and momentum distributions. In this approximation, a (so-called)
active electron is assumed to undergo a transition from its initial bound state
into the laser-dressed continuum
owing to its interaction with the laser pulse (cf.
Figure 1a), while the motion of all other electrons of the target is typically assumed to be unaffected. Not much need to be said here about the basic SFA theory, which can be found in various texts [
12,
13]. In this approximation, the probability for the strong-field ionization of atoms and for finding a photoelectron with asymptotic momentum
at the detector,
can then be expressed in terms of transition amplitudes as [
6]
and where
refers to the laser–electron interaction,
the time evolution and
to the
atomic potential as seen by the outgoing electron.
Indeed, the two (strong-field) amplitudes and can be readily interpreted in terms of a (re-) scattered photoelectron and are often referred to as the direct and rescattering amplitudes, respectively. In this work, we shall focus especially on the direct amplitude that describes those photoelectrons which are directly released from the target atom by the laser potential, , and then freely propagate within the laser field as Volkov solution , until they reach the detector. Indeed, this amplitude often provides a good approximation for most strong-field ionization processes and, in particular, if the laser field is not linearly polarized. Typically, the following assumptions are made to further simplify the amplitude (2):
The initial-bound state is entirely determined by the atomic potential and is not affected by the laser field.
The photoelectron with asymptotic momentum arrives as plane-wave at the detector, i.e., .
Once the electron is released from the atom, the atomic potential does not affect its (electronic) motion within the continuum.
Often, moreover, a Coulomb potential has been applied, in line with the ionization potential of the target atoms, and the initial state has been taken just as ground state orbital in this (Coulomb) potential. With these assumptions in mind, the momentum distribution of the photoelectrons can then be expressed by a closed (analytical) formula. Obviously, however, these assumptions neglect both a proper representation of the initial state of the atoms as well as the (static) potential of the photoion upon the outgoing electron wave (continuum) and, hence, quite major parts of the electronic structure of the target atoms.
Several, if not most, of these assumptions can be easily released, if the initial and final states are consequently expressed in terms of partial waves as typical for atomic structure theory [
15]. In such a partial wave expansion, the representation of the initial bound and (Volkov) continuum states can be incorporated along with the parameterization of the short an intense laser pulses. It is
this representation of the (direct) transition amplitude
on which we shall focus in the implementation below and which paves the way for extending the strong-field theory towards the study of non-dipole contributions in light-atom interactions as well as towards many-particle correlations in strong-field ionization processes.
2.2. Partial-Wave Representation of Strong-Field Amplitudes
In the derivation of the SFA amplitude ((2) and (3)), indeed, no additional assumptions have to be made about the form of the atomic potential
or about the potentials
and
of the driving laser field, but which—of course—affect the Volkov states
. If, for the coupling of the radiation field, we restrict ourselves to the dipole approximation [
] and the velocity gauge [
], the vector potential of an elliptically-polarized laser pulse can be written in terms of its spherical tensor components as
and where
denotes the (real-valued) amplitude,
the pulse envelope,
the fundamental frequency, and where
refers to the carrier–envelope phase of the laser pulse. In this notation, moreover, the (complex) polarization unit vector
defines the orientation of the polarization ellipse in terms of the ellipticity
. The vector potential (
4) therefore implies already all the properties
and
of the laser field which can be controlled experimentally.
A partial-wave representation of the amplitudes ((2) and (3)) also enables one to adapt both the initial bound state
and the final continuum state to the target potential of interest [
15,
19]. This is readily achieved, for instance, by using (self-consistent) solutions from atomic structure theory. For the outgoing photoelectron
moreover, one only needs to replace the partial wave
in the expansion of a plane-wave Volkov state by the corresponding solutions of either a Coulomb–Volkov or distorted-Volkov state in order to account for a realistic potential of the target, including the associated Coulombic and non-Coulombic phase shifts [
20,
21]. Here, we shall provide only a brief discussion of the theory, just enough to follow our implementation below, while all further details are given in Refs. [
15,
19].
For the direct SFA transition amplitude (2), Equations (23) and (24) of reference [
15] display a rather lengthy formula that is written in a basis of well-defined total angular momenta for the initially bound and the final (photo-) electron. This expression depends explicitly on the Volkov phase and all the parameters of the driving laser pulse, and it also accounts for the spatial dependence of the active electron in terms of the reduced one-particle matrix elements
of the momentum operator, as typical for atomic structure theory. The advantage of such an expression in a spherical basis arises from the—prior and separate—integration over all radial and spherical coordinates. This expression therefore also enables one to readily incorporate and discuss different contributions from the electron–photon interaction and the representation of the active electron without any need to re-derive the transition amplitude(s) for every target potential and/or laser (pulse) configuration separately. Below, we shall focus especially upon realistic (single-electron) initial states and the improved representation of the continuum for the outgoing photoelectron. It is
this partial-wave representation of the direct SFA amplitude which makes the present extension an integral part of the
Jac toolbox and which goes well beyond of what has been used originally in the SFA. However, since the expression in Ref. [
15] is still restricted to the electric-dipole approximation, it neither accounts for the magnetic field nor any spatial dependence of the electric field, though this can likely be done as well (cf.
Section 4).
2.3. Implementation of the (Direct) Strong-Field Amplitude
Like in atomic structure theory, a partial-wave representation of all (strong-field) amplitudes enables one to deal quite independently with atomic potentials, the Volkov state for the outgoing electron or the laser–electron interaction in terms of the given laser parameters. Such a representation also helps identify the building blocks for computing the photoelectron spectra and/or momentum distributions, as they are observed experimentally. In the Jac toolbox, we therefore aim to distinguish between the target and laser parameters as the input of a computations, and the generated observables (spectra) as output. A simple access to all individual input parameters will enable us to compute the strong-field ionization amplitudes in quite different approximations. Although the current implementation is still restricted to the direct SFA amplitude in velocity gauge, the same or very similar building blocks will occur if other amplitudes or gauges are to be considered in the future. The use of partial waves and spherical tensors even ensures that these amplitudes can be readily combined with atomic structure codes to further include electronic correlations and relativistic contributions to the strong-field ionization studies.
Not much need to be said about
Jac itself, the Jena Atomic Calculator [
16] that supports atomic (structure) calculations of different kind and complexity and that has been summarized at various places [
22,
23,
24]. Apart from energies and wave functions for open-shell atoms and ions, this toolbox also helps compute a good number of their excitation and decay processes. With the design and implementation of
Jac [
25], we moreover aim for establishing a
descriptive language that (i) is simple enough for both, a seldom or more frequent use of this toolbox, (ii) emphasizes the underlying atomic physics and, furthermore, (iii) avoids most technical slang as common to many other electronic structure codes. The implementation of
Jac is based on Julia [
26,
27], a recently developed programming language for scientific computing, and supports its use without much prior knowledge about neither the language nor the code itself.
Figure 2 displays a few selected atomic processes that are presently supported by
Jac, and which are shown together with useful features and control parameters for calculating strong-field amplitudes. The set of parameters in the right panel of
Figure 2 hereby indicates how between different pulses we shall distinguish in these computations geometries and/or gauges for the coupling of the radiation field, and as far they have been worked out until the present. In particular, the initial bound and final Volkov states of the (photo-) electron just appear in the reduced matrix elements
and can be taken either as a hydrogenic orbital, scaled upon the ionization potential, or as realistic one-electron wave function. Here, the continuum orbitals are generated in the static potential of the photoion and with energies as measured at the detector [
28].
Special care has to be taken about the envelope of the laser pulses. In a spherical-wave expansion, this envelope enters the direct amplitude in terms of (so-called) pulse-shape integrals
and
, cf. Ref. [
15]. These one-dimensional, (time) integrals are often obtained numerically but can be evaluated also analytically for continuous beams and a few other forms of the envelope. In our implementation below, the envelope (shape) of the laser pulse is accessed by a proper (abstract) data type, quite similar to the frequency, intensity, number of cycles and the polarization of the incident pulses. In typical applications of
Jac, one needs to select these parameters based on the given setup of the experiment and different practical considerations in order to keep the computations feasible.
2.4. Data Types for Modeling Photoelectron Distributions and Above-Threshold Experiments
From a physics viewpoint, we normally wish to trace back the simulation of different spectra and photoelectron distributions to just computing the (direct) SFA amplitude from above, though for specifically selected target atoms, approximation of wave functions and parameters of the laser pulse. Obviously, this requires simple access to all these data as well as special care to bring them together with the internal calls of the program. To facilitate the communication with and the data transfer within the program, the Jac toolbox is built upon a large number of data structures in order to specify useful and frequently recurring objects in such computations, and which also establish their language elements. Two prominent examples for such data structures, that frequently appear in atomic structure theory, are an Orbital for specifying the quantum numbers and radial components of single-electron orbital functions, or a Level for the full representation of an approximate initial or final bound state of the target atoms, and which itself comprises all information about the orbitals, the coupling of the angular momenta and the mixing of the many-electron target states. These target states are typically obtained self-consistently in a Dirac–Fock–Slater potential and hence are based on orbitals in line with the given target. Jac’s explicit set of data structures has been enlarged for the present update of the code by several types and now helps compute, analyze and explore the desired photoelectron spectra for different laser pulses and targets. In total, there are about 250 of these data structures in Jac, though most of them remain hidden to the user.
To model different strong-field ionization measurements, we wish (and need) to characterize especially the incident laser pulse in terms of its frequency, intensity, shape and polarization as well as the observables (spectra) to be simulated. In addition, we wish to control the target potentials and representation of the Volkov states in the strong-field amplitude. All this input is very central to the implementation and must be readily accessible by the given hierarchy of data structures. While we shall not explain these structures in all detail here,
Figure 3 displays a few of them from
Jac’s
Pulse and
StrongField modules. The abstract type
Pulse.AbstractEnvelope (middle panel), for example, just deals with the envelope of the laser pulse and comprises various concrete types for specifying a particular shape, pulse duration or number of cycles. Similarly, the data type
StrongField.AbstractSFAObservable enables one to specify the observable of choice and its resolution. All this information about the observable, target and pulse parameters finally define (an instance of) a
StrongField.Computation (lower panel), and which can be utilized in
Jac analogue to the previously implemented
Atomic.Computation [
16,
29] or
Cascade.Computation [
30,
31].
Finally,
Table 1 displays several other data structures that are relevant as well for the computation and analysis of strong-field photoelectron distributions. They are explained only in brief, while further details can be obtained from
Jac’s User Guide [
25] or by just using Julia’s help facilities [
32]. The definition and hierarchy of these data structures however nicely illustrate how they help implement different strong-field ionization scenarios and, hence, a wide range of potential applications in atomic and atto-second physics. In the next section, we make use of these data types to simulate various energy, angular and momentum distributions for a krypton target.
3. Energy and Momentum Distributions for Atomic Krypton
In the literature, the SFA has been frequently applied for comparing the energy and momentum distributions with experiments and for studying pulses and targets of quite a different sort. In these computations, more often than not, the active electron has initially been assumed to be in a hydrogenic
state, and by just matching the ionization potential to the target of interest. However, such a simple approach provides only little insight into the role that the target atoms play in strong-field ionization. Here, we wish to demonstrate that our partial-wave representation of the SFA amplitude enables us to adopt the initial-bound and final-Volkov states to realistic target potentials. We also show how the ATI spectra and momentum distributions can be obtained for pulses of different intensity, polarization and pulse duration. All these computations are performed by applying the
Jac toolbox [
16], which integrates the electronic structure and a good deal of atomic processes within a single computational framework, and which has now been expanded to facilitate the simulation of strong-field ionization distributions. For the sake of convenience, all simulations below are performed for krypton (
14 eV) and a right-circularly polarized,
cycle driving laser pulse with wavelength
nm, intensity
W/cm
and carrier–envelope phase
. Here, we shall not compare our implementation with experiment or previous computations, which have been done recently for a number of other targets [
19].
3.1. Above-Threshold Energy Spectra
Often, the observed ATI spectra can be qualitatively reproduced by simply using the SFA and plane-wave Volkov continuum states, since the peak structure of these photoelectron spectra itself arises from the interaction of the (quasi-) free electron with the electric field of the ionizing laser pulse. For these reasons, most energy spectra also exhibit distinct peaks, which are just spaced by the photon energy of the incident laser beam. These peaks become easily visible by measuring the photoelectron energy for a fixed azimuthal angle along some line in the polarization plane. Besides the selected laser parameters, these energy spectra depend of course also on the target atoms as well as on how the photoelectrons are described on their way to the detector, including the Volkov continuum and, possibly, even a re-scattering of the photoelectrons.
The ATI energy spectra of the strong-field photoelectrons can be modeled also by the present extension to the
Jac toolbox.
Figure 4 displays the (Julia) input which needs to be prepared by the user and which enables one to calculate such spectra for different targets and pulses. In this input, we have assumed that the
ground level of atomic krypton and the
final levels of the photoion have been computed before by the
Jac toolbox and are just provided by the variables
initialLevel and
finalLevel. Here, we make use of a slightly larger charge of the nucleus in order to adopt the
ionization potential to experiment. To characterize the laser pulse, moreover, we provide the wavelength, intensity and carrier–envelope phase and assume a
envelope as well as a right-circularly polarized plane-wave beam. Some of these given parameters first need to be converted to atomic units in order to be applicable in the computation of the field amplitude. We also specify here the velocity gauge and the electric-dipole (E1) interaction, even if these parameters must not be changed in the present implementation. The choice of a hydrogenic orbital with scaled nuclear charge can be made by a boolean in the
StrongField.Settings().
In the input above, we finally also specify as observable an SfaEnergyDistribution(), and which is to be calculated for and (i.e., along the x-axis), and for 200-electron energies in the interval eV. All this input together determines the (strong-field) computation comp::StrongField.Computation and can be readily adopted to many other experimental scenarios. All that is needed in Jac is to perform(comp, output=true) this computation, and where the optional parameter output=true just tells Jac to return the calculated data (amplitudes) to the user for printing and post-processing.
Figure 5 displays the photoelectron energy spectra, emitted along the
x-axis, for a krypton target and a right-circularly polarized laser pulse. The left panel shows the spectra as obtained for a computed with a hydrogenic
initial wave function with adopted ionization potential and for a plane-wave Volkov continuum (black-solid curve) as well as a Coulomb–Volkov continuum (red dashed curve). On the right panel, in contrast, the spectra are computed for an initial
orbital of neutral krypton and a plane-wave Volkov continuum (black-solid curve), a Coulomb–Volkov continuum (red-dashed curve) as well as a distorted-Volkov continuum (blue-dotted curve). In all these computations, a right-circularly polarized
pulse of wavelength
nm, intensity
W/cm
, carrier–envelope phase
and with just
cycles has been utilized.
Input quite similar to
Figure 4 can be employed also for studying the angle and momentum distributions of photoelectrons for different laser pulses and targets. While no further input data will be shown below, we refer the reader for details to the User Guide and the online documentation of the
Jac program. Moreover, rather moderate changes to this input will be sufficient in the future to expand the
StrongField module to other gauges, amplitudes or many-electron features. While such an expansion of the code indeed appears straightforward, major effort will still be needed for its implementation and testing.
3.2. Photoelectron Angular Distribution for Elliptically-Polarized Laser Pulses
In the electric-dipole (E1) approximation, the angular distribution of the photoelectrons is restricted to the
polarization plane and just reflects at fixed photoelectron energy the ionization probability in Equation (
1) for different azimuthal angles
. If, moreover, the laser field dominates the electron dynamics in the continuum, the observed photoelectron angular distribution (PAD) should also reflect the symmetry of the vector potential of the laser beam. In practice, however, a
Coulomb asymmetry in the PAD was (first) observed by Goreslavski et al. [
33] in the ATI of xenon gas targets and, since then, has been found to be a valuable testbed for improving theory. For lithium, argon and xenon, for example, the SFA theory was shown to reproduce this asymmetry, if a target-specific initial orbital function is chosen along with a distorted-Volkov continuum for the active electron [
34].
Figure 6 displays different photoelectron angular distributions in the polarization plane (
) for a krypton target. Angular distributions are shown for elliptically-polarized laser pulses with
(left panel) and
(right panel) at fixed photoelectron energy
according to the third ATI peak in
Figure 5. Different approximations are compared for these angular distributions: a hydrogenic
initial orbital together with a plane-wave Volkov continuum (black-solid curves); a self-consistent initial
orbital of neutral krypton together with a Coulomb–Volkov continuum (red long-dashed curves); the same initial
orbital but together with a distorted-Volkov continuum (blue-dashed curves). All these distributions are normalized on their maximum, while all other laser parameters are the same as in
Figure 5. Indeed, a self-consistent
orbital of neutral krypton together with a Coulomb–Volkov continuum (red long-dashed curves) leads to a clear rotation of the PAD as mentioned above. Moreover, the PAD no longer exhibits an inversion symmetry with regard to the origin because of the short duration of the laser pulse. If, in addition, the Coulomb–Volkov continuum is replaced by an distorted-Volkov continuum (blue-dashed curves), and which accounts for an outgoing electron in the potential of the Kr
photoion, the rotation angle still changes rather remarkably. In Ref. [
19], it was shown that such a distorted-Volkov continuum (often) leads for different targets to better agreement with experiment.
3.3. Photoelectron Momentum Distribution for Few-Cycle Laser Pulses
Theoretical photoelectron momentum distributions (PMD) have been calculated in the literature by means of quite different methods, and by making use of even a larger number of case-specific modifications to these methods. Generally, the experimentally observed symmetries of the PMD cannot be explained so readily by just applying a plane-wave Volkov continuum [
33], but can be improved further if the Coulomb potential of the residual ion is taken into account. In our implementation of the SFA direct amplitude, this is achieved by replacing the plane-wave Volkov continuum by either Coulomb–Volkov or distorted-Volkov states. For the low-energy photoelectrons with (say)
, the ionization probability is then often enhanced by up to an order of magnitude, if the ionic charge just increases from
to 1. This has been explained by the attractive Coulomb potential of the residual ion that pulls the electron back to the ion and hence reduces its kinetic energy. The low-energy part of the ATI spectra can be further improved by adding a short-range potential to the (long-range) Coulomb potential and by making use of distorted-Volkov states.
Figure 7 shows the photoelectron momentum distributions in the polarization plane (
) for the strong-field ionization of a krypton target. Momentum distributions are shown for circularly-polarized laser beams with three different CEP phases:
(left panel),
(middle panel) and
(right panel) and by applying a self-consistent initial
orbital of neutral krypton together with a plane-wave Volkov continuum. All further laser parameters are the same as in
Figure 5. Obviously, the PMD in this figure exhibits a (very) clear rotation since the photoelectrons are preferably emitted in the polarization plane along the maximum of the vector potential [
6], and which changes with the carrier envelope phase
. It will be interesting to explain with
Jac in future work how the Coulomb asymmetry and the CEP dependence act together upon the angular or momentum distributions, and, especially of the initial-bound and distorted-Volkov continuum states of different atomic targets, are taken properly into account. In these studies, both the Coulomb and short-range interactions can be easily incorporated into the continuum by just replacing the radial wave functions of the active electron.
4. Conclusions and Outlook
Up to the present, the SFA has been found as perhaps the most powerful method for predicting or analyzing the electron dynamics in strong-field ionization. Often, this approximation helps describe features in the observed electron distributions even quantitatively, if the initial-bound and final-continuum states of the photoelectron are well adopted to the target atoms, and if combined with a proper parameterization of the laser field. With the present implementation of the direct SFA amplitude into the Jac toolbox, this method can now be applied to different targets and strong-field scenarios. In particular, the implementation of the SFA in the partial-wave representation enables us to readily control (and replace) the wave functions and various details about the laser–electron interaction. It also enables us to extend this implementation for incorporating further interactions and mechanisms into the modeling.
Detailed calculations are performed for a krypton target as well as for different ATI spectra and PMD. These examples clearly show how the target potential affects the photoelectrons on their way to the detector and, hence, all the observed spectra. In particular, we have demonstrated how the electronic structure of the atomic targets can be taken into account in the representation of the active electron and how the dynamics of the outgoing electron can be readily controlled by applying different approximations for the Volkov continuum. Moreover, the use of partial waves and spherical tensor operators facilitates a simpler comparison of different pulse shapes and how they influence the observed ATI spectra and PMD.
Several extensions to the SFA are still desirable and appear feasible within a framework, which is based on a partial-wave representation of the associated strong-field amplitudes. While further effort will be needed to decompose these amplitudes into a form, suitable for computations, a few useful extensions concern:
Non-dipole interactions: For spatially-structured light fields, non-dipole contributions to the Volkov continuum usually arise from the spatially dependent Volkov phase [
35,
36,
37], and which need first to be expressed into a partial-wave representation in order to become applicable within
Jac. These
non-dipole terms beyond the widely used E1 approximation capture the combined-electric
and magnetic-fields upon the electron dynamics [
38,
39]. Their implementation into the
Jac toolbox will help predict the energy and momentum shifts at long wavelengths of the driving fields.
Coupling of the radiation field: Apart from the (commonly applied) velocity gauge, the direct amplitude can be also implemented in length gauge. This leads to more complicated pulse-shape integrals that also comprise the reduced matrix elements of the momentum operator, since the kinetic momentum then needs to be replaced by the (time-dependent) canonical momentum. While such an implementation requires further work, the direct SFA amplitude in length gauge was shown to provide more accurate results for the ionization of non-spherical
electrons [
40].
Rescattering amplitude: For laser pulses with proper polarization, the electrons are known to be partly
rescattered by the photoion, which then leads to processes, such as high-order ATI, the non-sequential double ionization, or even to high-order harmonics beyond the well-known
cut-off law [
41]. A partial-wave representation of the rescattering amplitude (2) is currently worked out and can be applied to account for realistic rescattering potentials.
High-harmonic generation: Similar to the rescattering ATI amplitude above, a recombination amplitude needs to be computed in order to obtain the dipole moment of emitted high-harmonic radiation. For modeling HHG, again, we expect to benefit from a re-formulation of the dipole amplitude in terms of partial waves and from including realistic initial and continuum orbitals [
42,
43].
Role of bound states: The coupling of the ground and continuum states to other excited (bound) states has been analyzed in the literature for just a (very) few selected HHG spectra [
44]. A partial-wave representation of the SFA amplitudes facilitates the coupling to excited states of the target and may help explain the formation and influence of (autoionizing) resonances in the HHG plateau.
Many-electron effects: A consequent partial-wave decomposition of all strong-field amplitudes enables one to incorporate many-electron contributions beyond a (spherical) short-range potential into the formalism. Apart from the self-consistent field and the mixing of important configurations, this also refers to the treatment of the multipole contributions (higher than E1), if the corresponding many-electron matrix elements are utilized [
45,
46].
Nonsequential double ionization (NSDI): When the photoelectron returns to the photoion, the electron can scatter inelastically under the ionization of a second electron. Theoretically, the NSDI is typically described semi-classically by using excitation and/or ionization cross sections for the second (ionizing) step of the process [
47,
48]. A partial-wave representation of all associated quantum SFA amplitude facilitates a coherent treatment of this nonlinear ionization process.
For all these desirable extensions, the partial-wave representation of the SFA [
15], and its implementation in
Jac provides a straight and conceivably the best way to advance theory and the light–atom interaction in strong fields.