Manipulating Twisted Electrons in Strong-Field Ionization

We investigate the discrete orbital angular momentum (OAM) of photoelectrons freed in strongfield ionization. We use these `twisted' electrons to provide an alternative interpretation on existing experimental work of vortex interferences caused by strong field ionization mediated by two counterrotating circularly polarized pulses separated by a delay. Using the strong field approximation, we derive an interference condition for the vortices. In computations for a neon target we find very good agreement of the vortex condition with photoelectron momentum distributions computed with the strong field approximation, as well as with the time-dependent methods Qprop and R-Matrix. For each of these approaches we examine the OAM of the photoelectrons, finding a small number of vortex states localized in separate energy regions. We demonstrate that the vortices arise from the interference of pairs of twisted electron states. The OAM of each twisted electron state can be directly related to the number of arms of the spiral in that region. We gain further understanding by recreating the vortices with pairs of twisted electrons and use this to determine a semiclassical relation for the OAM. A discussion is included on measuring the OAM in strong field ionization directly or by employing specific laser pulse schemes as well as utilizing the OAM in time-resolved imaging of photo-induced dynamics.

We investigate the discrete orbital angular momentum (OAM) of photoelectrons freed in strongfield ionization. We use these 'twisted' electrons to provide an alternative interpretation on existing experimental work of vortex interferences caused by strong field ionization mediated by two counterrotating circularly polarized pulses separated by a delay. Using the strong field approximation, we derive an interference condition for the vortices. In computations for a neon target we find very good agreement of the vortex condition with photoelectron momentum distributions computed with the strong field approximation, as well as with the time-dependent methods Qprop and R-Matrix. For each of these approaches we examine the OAM of the photoelectrons, finding a small number of vortex states localized in separate energy regions. We demonstrate that the vortices arise from the interference of pairs of twisted electron states. The OAM of each twisted electron state can be directly related to the number of arms of the spiral in that region. We gain further understanding by recreating the vortices with pairs of twisted electrons and use this to determine a semiclassical relation for the OAM. A discussion is included on measuring the OAM in strong field ionization directly or by employing specific laser pulse schemes as well as utilizing the OAM in time-resolved imaging of photo-induced dynamics.

I. INTRODUCTION
Phase vortices in light and matter are the product of the orbital angular momentum (OAM) carried by free particles resulting in a rotating wavefront. The OAM of a free particle is a quantum number that can be measured and manipulated, with many of the same properties as spin, yet it has only been observed in experiments fairly recently in photons (3 decades ago) [1] and even later for electrons (1 decade ago) [2] (see [3,4] for reviews). It has been found that these twisted photons or electrons have huge potential in imaging, and recent experiments with electron vortex beams have revealed their role in chiral energy-loss spectroscopy and magnetic dichroism [3].
The study of twisted light and electrons in strong field processes is a very new but rapidly developing topic. In high-harmonic generation (HHG) it has been shown that a twisted driving field can be used to produce twisted UV light [5]. Later this was made time-varying [6] and further extended to the generation of UV light with torus knot topology [7]. The effect of twisted electrons has been investigated in HHG [8]. However, the effect of OAM on photoelectrons in strong-field ionization has not been studied in detail.
Recently, the OAM content of photoelectrons was computed using the strong field approximation (SFA). This work focused on reaching exceedingly high OAM values for quasi-relativistic field intensities [9] and terahertz fields [10]. Twisted electrons have also been studied in the context of rescattering with a well defined OAM [11]. However, currently there is very little work focused on the generation, measurement and manipulation of photoelectrons carrying OAM for strong field ionization. That said, there is a still a lot to learn from related strong field work. Through conservation of momentum, photoelectrons ionized by a linearly polarized laser field will adopt the same OAM as the magnetic quantum number of the bound state of origin [11]. The linear field does not alter the magnetic quantum number or the photoelectron's OAM, and there is a direct mapping between the two. Furthermore, it has been demonstrated that circularly polarized laser fields preferentially tunnel ionize p electrons, whose magnetic quantum number m dictates that they are 'counter-rotating' with respect to the laser field [12,13]. Depending on the helicity of the laser field, either m = 1 or m = −1 electrons dominate the signal, while ionization from orbitals with m = 0 is strongly suppressed. Thus, given the interaction with the magnetic quantum number, tunnel ionization via circularly polarized light may be a route to generate, control, and measure photoelectrons with non-trivial OAM.
A more explicit example of twisted electrons in multiphoton processes is the formation of interference vortices [14][15][16][17][18][19][20][21][22][23]. Computations of the momentum distribution of photoelectrons ionized via two time-delayed, counter-rotating circularly-polarized laser fields reveal a Fermat spiral interference pattern. Initially treated theoretically for singlephoton ionization [14] and later two-photon ionization [15] of helium, such spirals have since been observed in both single and double ionization of molecules [16,17]. Experimental verification came in Refs. [18,19], where the spiral interference pattern was demonstrated for three-photon ionization of potassium. Calculations using the multielectron ab initio R-Matrix with time dependence (RMT) method fully corroborated this effect [23], and revealed the threedimensional characteristics of the interference vortices. Experimental progress in this domain continues to gather pace, and vortex measurements at mid-infrared wavelengths are being carried out [24]. Interference vortices are indicative of photoelectrons carrying OAM, as the interference of two states with differing OAM produces such spirals. The above-mentioned multiphoton studies have provided explanations in terms of the magnetic quantum number in the bound states and scattering states to which the photoelectron is promoted, but do not go as far as discussing the OAM of the free electron or its measurement and control. In this work we examine the interference vortices in the non-perturbative regime, where a description of promotion to scattering states will not hold. However, analysis of coand counter-rotating electrons [12,13] is valid and provides additional insight.
Twisted electrons and the OAM distributions are the focus of this article, and with these tools we present an alternative understanding of the vortices. Using the SFA [25], as well as the time-dependent Schrödinger equation (TDSE)-based solvers, Qprop [26,27] and RMT [28][29][30], we compute both the photoelectron's momentum distribution and its constituent OAM components. The resulting momentum-space vortices may be understood by deriving an interference condition using the SFA, and by recreating the same pattern using the interference between states of different values of OAM. We also identify a semi-classical relationship that demonstrates a link between the OAM and above threshold ionization peaks of the photoelectron. Understanding the OAM in photoelectrons, which undergo strong-field ionization via circularly-polarized light, imbues the seemingly plain momentum distributions with previously unseen structure revealed only by such interferometric schemes. This opens the question: is there a generalised way to measure, directly or indirectly, the OAM of photoelectrons in strong-field experiments? Given the high correlation between the photoelectron OAM and the quantum numbers of the initial, bound state, consideration of the OAM could be extremely useful in relating physical observables to the bound state and disentangling interferences in time-resolved measurements. Furthermore, the potential for the outgoing twisted photoelectron states to be chiral could aid the detection and spectrographic measurement of chiral molecular targets [31].
The article is structured as follows. In section Sec. II we outline the theoretical background for the three methods employed in this article, the SFA (Sec. II A), the single active electron (SAE) 3D-TDSE solver Qprop (Sec. II B), and the multielectron TDSE solver RMT (Sec. II C). In Sec. II D we present the methodology for computing the OAM distributions for each model. Next we focus on the interference vortices themselves. In Sec. III, we present the derivation of the vortex interference condition using the SFA (Sec. III A) and compare the results of all methods and that condition (Sec. III B). While in Sec. IV we turn our attention to the orbital angular momentum (OAM), in Sec. IV A we compare the photoelectron OAM distribution from all methods and in Sec. IV B we use this to construct the interference vortex entirely from twisted electron states. In Sec. V we make our concluding statements and discuss potential direct and indirect measurement schemes. We employ atomic units throughout (denoted a.u.), where the elementary charge, electron mass and the Planck constant are set to one, e = m =h = 1.

Laser field
We consider a bicircular field comprised of two time-delayed, 400-nm pulses, described by the vector potential where ± denotes the rotation of the circular field. Each field is given by where A 0 = 2 U p is the peak vector potential strength, which is related to the peak electric field amplitude by A 0 = E 0 /ω. The laser frequency is ω = 0.114 a.u. (corresponding to a 400 nm wavelength), N c is the total number of laser cycles, and δ is the time delay between the pulses, which we will typically take to be δ = 2πNc ω so that the two field envelopes do not overlap. Each pulse reaches an associated laser peak intensity of 2.5 × 10 13 W/cm 2 and the field retains this profile for time t ∈ [0, 2πN c /ω], and is zero otherwise.

A. Strong Field Approximation
Our starting point in the SFA [25] is the transition amplitude for direct ATI from the initial bound state |ψ 0 to a final Volkov state with drift momentum p given by [32][33][34] where d(p, t ) = p + A(t )|V |Ψ 0 , |Ψ 0 is the ground state of the target and the action is given by Here, I p is the ionization potential of our target. We employ the saddle point approximation, seeking stationary action for the integration variable t , Now the probability distribution can be computed from Eq. (3) as where the prefactor c(p, t s , t), derived from application of the saddle point approximation and also includes the t independent phase from Eq. (3), is given by In order to capture the essential physics of the interference vortices, we make some additional assumptions. The minimal requirement is to have interference between two photoelectron orbits deriving from each pulse. Thus, in the SFA model we take only two ionization events from near the peak of each laser pulse, A + (t) and A − (t − δ), but otherwise neglect the pulse envelopes to allow for an analytic description. Thus, in the SFA model, the vector potentials are defined by where the ± denotes the direction of rotation of the field. In this way we are able to separate the two contributions of ionization at the times t + and t − + δ, given by In Eq. (9) we have written the solutions in spherical momentum coordinates p, θ and φ as they are the most natural for this system. The transition amplitude can then be written as a sum of the two separate solutions The actions S + (p), S − (p) and prefactors c + (p, t), c − (p, t), d + (p) and d − (p) have the times t + or t − + δ substituted into them such that and Following the prescription of [12] we split the prefactor d(p, t) into two parts where k = p + A(t) and k(p, t), θ k (p, t) and φ k (p, t) are the spherical coordinates of k. The function f nlm (k, θ) is dependent on the specific bound state, which in our case is the ground state of neon, while φ k is the so-called the tunneling angle [12], which leads to the enhanced ionization of 'counter-rotating' electrons. It will continue to play an important role in the description of the interference vortices.

B. Qprop
The 3D-TDSE calculations were performed using the latest version of Qprop [27], where a very accurate method for the calculation of photoelectron spectra (PES), dubbed i-SURFV, is implemented. In short, Qprop is a velocity gauge 3D-TDSE solver that allows studies within both the single active electron (SAE) approximation and manyelectron systems via the solution of the time-dependent Kohn-Sham equations. From a computational viewpoint, the velocity gauge formulation of the strong laser-matter problem appears to be more efficient, providing results of accuracy equivalent to the length gauge with much less computational effort. For our SAE model of the Ne atom, we have employed the model potential of Ref. [35] and started the time-dependent simulations from the 2p 0 , 2p −1 and 2p +1 bound states. These states were computed via imaginary time propagation. The pulse parameters are those given in Sec. II 1 and, specifically, the effective potential has the form where Z = 1. For neon the coefficients a i are a 1 = 8.069, a 2 = 2.148, a 3 = −3.570, a 4 = 1.986, a 5 = 0.931, a 6 = 0.602 [35], which gives the correct ionization potential of I p = 0.79 a.u. In this computation we considered angular momenta up to L = 14 and all magnetic sublevels.

C. R-Matrix with time dependence
The R-matrix with time dependence method (RMT) is an ab initio approach that solves the time-dependent Schrödinger equation for single ionization of multielectron atoms, ions, and molecules in arbitrarily-polarized, strong laser fields [28][29][30]. The method divides position space into two distinct regions. The inner region, which usually extends to around 20 Bohr radii (20 a.u.), confines the target nucleus. In this region, a many-body wave function is treated, accounting for electron exchange and correlation. The outer region, which extends to larger distances, contains a single ejected electron that interacts with both the laser field and the singly-ionized residual target.
Our description of the neon target is provided in a previous work [36]. In short, we couple a single electron to the Ne + residual ion, expanding the wave function in LM L Sπ symmetries up to L = 99, retaining both the 1s 2 2s 2 2p 5 and 1s 2 2s2p 6 Ne + residual ion states. The photoelectron wave function is calculated on a finite-difference radial grid in the outer region, which extends to 3400 a.u.. The laser pulse configuration is that given in Sec. II 1. The wave function is propagated in time for a total of 35 fs, which gives sufficient time for the ionizing wavepacket to extend to large distances within the outer region. The photoelectron momentum distribution is obtained by performing a Fourier transform of the ejected-electron radial wavefunction.

D. Computing OAM Distributions
In order to compute the orbital angular momentum distributions we will take two approaches. Firstly, given a transition amplitude such as that in the SFA we can transform it into the 'OAM basis' by computing its Fourier series coefficients [9] M lv (p || , Plotting the interference vortex condition over the polarization pxpy-plane of the laser field. The three spirals are plotted using 3 different sets of n from Eq. (24); the first 7 valid n, the next 8 and then the next 9 are used, from left to right. The minimum valid value of n is 37 for our field, target parameters and φ = π. Using a delay δ = 4 × 2π/ω, an angular frequency of ω = 0.114 a.u., Up = 0.027 a.u., corresponding to a peak laser intensity of I0 = 2.5 × 10 13 W/cm 2 . The ionization potential is Ip = 0.79 a.u. for a neon target.
However, given that this is a transformation into cylindrical co-ordinates, and many TDSE solvers employ a basis of spherical harmonics, some computational power can be saved by converting a distribution expanded in such a basis, so that the transformation of Eq. (15) reduces to where tan θ = p ⊥ /p || , and L is the maximum value of the azimuthal quantum number l required for adequate convergence.

III. VORTEX INTERFERENCE
A. Vortex interference condition from the standard SFA implementation In this section we derive a condition for the interference vortices from the SFA description. Such conditions can help to give some insight into the physical dynamics of the process. The starting point to derive the condition is Eq. (10). We can demonstrate that the prefactors c + (p, t) and c − (p, t) are the same and given by Thus, c ± (p, t) can be factored out and will not play a role in the interference. For the prefactor d ± (p), co-ordinates k and θ k of f nlm (k, θ k ) can be shown to be the same when evaluated at either t + or t − + δ. From the saddle point Eq. (5) we can determine Thus, we can deduce Hence, f nlm (k, θ k ) can be factored out and we are left with the e imφ k term, which (as presented in [12]) will act to weight the contribution from each laser pulse depending on the value of m. For non-zero m this will lead to a loss of contrast of the interference fringes as electrons 'co-rotating' with the field will be preferentially ionized. Now we may write the transition amplitude as where C contains all the prefactors. The term e im∆φ k (p) acts to switch off or blur the interference in the way described.
The interference condition can be derived by maximizing the term e i∆S(p) , where Setting ∆S(p) = 2πn, where n is an integer we can find its maxima and thus a condition for the interference vortices This can be re-arranged to give the equation of a spiral From Eq. (24) for any particular angle φ there is a minimum value of n, below which there will be no real solutions for the radius p. This is reminiscent of similar interference conditions that may be derived in the SFA e.g. to describe above-threshold ionization (ATI) peaks [37]. The condition Eq. (24) is capable of describing a varying number of spiral arms. In FIG. 1 we plot the spiral for different values of n using both the positive and negative branches.
In panels (a), (b) and (c) we use the first 7, next 8 and the following 9 integer values of n, respectively and both branches of Eq. (24). Using more and higher values of n leads the spiral arms filling a larger region, with some arms becoming crowded in the center, more evenly spaced for higher p and finally developing gaps for the highest p. Thus, this suggests that in specific regions there will be a particular number of evenly spaced spiral arms (similar to Fermat spirals), which increases with p. In short, we expect the interference vortices to have an increasing number of spiral arms with p.

B. Computation of interference vortices
In this section we present numerical calculations of strong field ionization of neon via the two counter-rotating fields described above. In FIG. 2  The interference vortex condition is superimposed on all plots and matches very well, clearly capturing the core behaviour. It is notable that such a simplified SFA model is able to capture the features of the vortex. There are however differences, such as the hole in the center of the SFA results, which is less visible in the RMT and Qprop outcomes. This is likely due to the Coulomb force, which will decelerate the outgoing electron wavepacket and thus reduce its final momentum closing the hole. This can be roughly accounted for in the SFA if the radial final momentum is taken to be lower than the initial, which will reduce the size of the hole. One clear aberration is the discontinuity along φ = 0 in the SFA results. In the SFA φ is fixed between ±π yielding a discontinuity. In the current construction the interference vortices can only be reproduced properly if φ can vary over a range larger than 2π in some regions, as it does in the interference condition.
In the RMT and Qprop momentum distribution calculation the spherical harmonic expansion was limited to to angular momenta l > 3. The rationale for this is that the spiral interference arises from an ionization process involving at least 7 photons, which typically leads to population of states close to l = 7. Population in states of low angular momenta are likely due to processes involving intermediate excited states, which are not accounted for by the SFA treatment. At the specific field parameters used here (peak intensity of 2. We find that the intermediate state effects are sensitive to the laser pulse parameters, and their contribution may be reduced if alternative parameters are used. We have verified that for higher intensities (10 14 W/cm 2 ) or long pulses (10-cycle sin 2 pulses), the low values of angular momenta do not play a significant role, and the interference vortices arising from at least 7-photon ionization from the ground state are the dominant effect. In fact very similar spirals can be produced with QProp and RMT using all angular momenta for an intensity of 10 14 W/cm 2 . This is aided by the fact the spiral fringes do not vary considerably with intensity in this regime. For longer pulses, the interference vortices are visible along with the characteristic ATI rings/ peaks, with the number of arms of each spiral increasing with the order of the ATI processes. A similar effect was observed in lower-order ATI processes in recent experiments [18].
Although it is difficult to resolve by eye, the number of spiral arms does vary, in FIG. 2, increasing with momentum. We have verified this by performing Fourier analysis on the momentum distributions across rings of fixed cylindrical radial momentum p ⊥ so that we can count the number of spiral arms intersecting each ring. This is also corroborated by the very good match with the interference condition. Thus, very good agreement can be achieved with a relatively simple SFA model, where many features can be identified. However, if we now examine the OAM of the outgoing photoelectron we can achieve additional insight and understanding. while positive OAM is dotted), and it is the interference between these that leads to the interference vortices. The RMT and Qprop OAM distributions are very similar but they differ somewhat from the SFA. In particular, there is a slight shift in the peak positions and relative heights. The highest peak in the SFA is for l v = ±9, while it is l v = ±8 for RMT and Qprop. This is related to the larger 'hole' in the SFA results (FIG. 2) and most likely can be traced to Coulomb distortions. There is also a shift of around half a photon's energy between the SFA peaks and those from RMT and Qprop. On each panel in FIG. 3  final OAM. For m = −1 electrons (not shown), the opposite tendency appears, and states of positive OAM dominate.

B. Building interference vortices from twisted electrons
Using the distributions presented in FIG. 3 the interference vortices can be reconstructed. This provides particular insight into the meaning of the OAM in a strong field context. To recreate this interference, for the case of m = 0, we will consider a distribution of outgoing photoelectrons with two opposite values of OAM ±l v . For an initial state of general m one should consider two states with OAM m ± l v . This can be written as Here |ψ m±lv,p (t) is electron vortex state with OAM m ± l v and momentum p = (p || , p ⊥ ) and w(p) provides a weighting over the cylindrical momentum coordinates. Now, if we project this state onto a 3D plane wave momentum state we can examine the interference in the p x p y -plane, The m cancels in the interference fringes, as they depend on the difference of the values of OAM. We have used the following momentum representation of Bessel electron vortex states [4] p|ψ Thus, we identify the interference between the two vortex states as leading to the condition This result generalizes the expressions derived in Refs. [14,15] for one-photon and two-photon ionization to the strong field regime, and describes a Fermat spiral with 2l v arms. In FIG. 4 [(a-b) and (d-e)], we plot examples l v = 1 and l v = 2, for different values of δ, which inversely sets the growth of the radius with respect to the angle φ. These spirals are reminiscent of those presented in Refs. [14,15], and observe a similar delay dependence. In panels (c) and (f) we plot a reconstruction of the full interference vortices by combining pairs of OAM in different momentum regions as indicated by the distributions given in FIG. 3. The SFA and RMT distributions for m = 0 are recreated in panels (c) and (f), respectively. This interference construction will also hold for m = 1, the only difference being positive OAM are suppressed. This explains why the original spiral condition works so well for both cases. For each pair of OAM values included the weighting function w(p) is set as a Gaussian over p ⊥ to reflect the position, heights and widths of the respective OAM distributions. This clearly demonstrates the link between interference vortices and the OAM of the photoelectrons. However, we have two separate spiral conditions. From the interference of photoelectrons of opposite OAM, Eq. (28) leads to a condition for Fermat spirals with a fixed number of arms determined by the OAM. But due to the localization of OAM of the photoelectron, this condition is only valid in specific momentum regions and the number of arms of the interference vortex does not remain fixed. Alternatively, from the SFA, Eq. (23), we derived a condition capable of describing a varying number of arms. It is, however, possible to reconcile these two, seemingly disparate descriptions. By employing the semiclassical relation in Eq. (23) we recover the same form as Eq. (28), where δ = 2πN c /ω is the number of laser cycles in the delay. Thus, there is a p-dependent semiclassical expression for the OAM, l v (p), which maps directly onto the ATI equation and can be interpreted as each additional photon contributing ±1 to the OAM. The integer n is shifted by N c × l(0), where l(0) corresponds to lowest possible OAM given l(0) photons are required to overcome the barrier I p + U p . The expression given by Eq. (30) maps the SFA condition onto vortex Fermat spirals. It demonstrates, via this semiclassical relation, in specific regions where l(p) is an integer, that Eq. (30) will behave like the condition given by Eq. (28).

V. DISCUSSION AND CONCLUSIONS
In the previous sections the main dynamics of the interference vortices in the strong field regime were captured and new light was shed on their formation by analysing the orbital angular momenta (OAM). In this section we will place the interference vortices and photoelectron OAM in a wider context, making parallels with other systems as well as discussing the potential to measure the OAM and exploit it in time-resolved imaging.
Vortices and spirals in a physical system are indicative of rotational dynamics and symmetries. When we talk about vortices in quantum matter the physics of liquids is the first example that comes to mind [38]. In the present context, however, the analogy to Bose Einstein condensates (BEC) of dilute atomic gases is more appropriate [39]. Vortices in BEC can be created by rotations (like in a bucket of liquid Helium) [40], that is "stirring with an optical spoon" [41] or a "phase imprinting" [42]. While the vortices of charges higher than ±1 are not stable in these systems, rotating BECs form beautiful Abrikosov lattices of charge one vortices. Still, interference of matter waves was, from the very beginning, proposed as a way to observe topological defects [43][44][45], stimulated by non-linear optics [46][47][48][49]. Interference of a vortex state with a plain BEC leads to the famous fork patterns [38,39]. The interference of two BECs with vortices leads to spiral patterns similar to those discussed here [38,39], except that i) the patterns in BEC are formed in the density profile in real space and ii) they are formed by a macroscopic quantum state, and are thus robust with respect to decoherence.
The question of how to measure the OAM of photoelectrons in strong field ionization may also be answered via interference. The main difference here is that interference structures are formed in the electron momentum space. The interference vortices are in principle measurable in experiments with angular resolution of the electrons created in ionization. Indeed, in the multi-photon regime the experimental measurements clearly showed the Fermat spirals [18,19]. Still, the problem is how sensitive the interference patterns may be to the parameters of the system. Will incoherent effects in the strong field regime blur the interference? Preliminary simulations of focal averaging in our regime of strong and non-perturbative laser fields, indicate that the feasibility of observing the interference pattern is robust. The interference fringes do not vary considerably with intensity until U p > I p , which for neon at λ = 400 nm occurs at intensities greater than 7 × 10 14 W/cm 2 . Going to higher intensities or longer wavelength would perhaps require the use of pre-designed laser pulses with a flat top, for instance [50,51], using more sophisticated methods. It is also important that the strong field ionization leading to the interference vortices is the dominant process. For lower intensities and multielectron targets, coupling to excited states can occur. Among other effects this could cause superposed interference vortices from both the ground and excited states and the combination can be difficult to disentangle. For larger intensities and longer pulses these effects can be minimized. For more complex targets such as molecules contributions from multiple states or molecular centres may also present difficulties as again the superposition of many interference vortices, associated with different ionization potentials and magnetic quantum numbers, could make them difficult to resolve.
For the case of circular fields we have presented a clear path to determine the OAM of the photoelectrons, via a second field of opposite helicity. What about elliptical, linear or more complex fields? In the case of linear fields, as they do not interact with the OAM, the same scheme as discussed in this work would not be possible. It would, however, for any field be possible to produce a well-characterized reference for measurement of OAM via interference using a circular pulse. Interference of pairs of twisted electron states will always yield a spiral in the p x p y -plane with the number of arms being equal to the difference between the two OAM values. Thus, a careful choices of circular pulses following any target pulse (after a delay) could reveal the OAM distribution. Exploiting the localisation of OAM in momentum space for a circular field would lead to different interference vortices in different regions, in effect 'scanning' the OAM of the photoelectrons ionized from the target pulse. However, this method would start to break down if multiple OAMs populate the same region of momentum space for a particular target pulse. Multiple vortices may be difficult to resolve and lead to interferences too complex to analyse.
Another approach is to measure the OAM of the outgoing photoelectrons directly. This would require significant alteration to standard strong field experimental setups, such as velocity map imaging (VMI) [52,53] or a reaction microscope (ReMi) detection system [54][55][56], so needs some justification. However, it would not only be a much more generalized and robust way to determine the OAM. Direct measurement could also allow for measurement of incompatible observables, which is not obviously possible using interferometric schemes via tailored laser fields. Incompatible observables are exploited in various aspects of quantum information and metrology including quantum key distribution, Bell inequalities and quantum sensing [57,58]. Additionally, the measurement of OAM could enhance existing strong field procedures. For example, it would allow direct measurement of the quantum magnetic number m of the initial state. Furthermore, it has already be shown that recolliding twisted photoelectrons are sensitive to chiral molecular targets [11], which may be exploited for imaging. This may also have implications for photoelectron holography [59], where recollision/ reinteraction with the target plays an important role and chiral phases could be revealed. The inclusion of fields with high ellipticity and recollision probability, such as bi-circular fields, would enable manipulation of the OAM while retaining strong interaction with the target. Thus, OAM measurements could open up a range of possibilities for strong field research.
So what is the state of the art in OAM measurement and how can this be applied to strong field systems? In [3,4] a review of state of measurement in OAM is given. A range of methodologies are available, an example of a typical approach is using diffraction via holograms or particularly shaped apertures, converting the phase information into spatial information. The first/most well-known method uses a fork phase mask [60,61] but there are difficulties such as the the limited transparency of the mask. In [62] the OAM sorter developed transforms the OAM into spatial position, reminiscent of a Stern-Gerlach style measurement of spin. Such a scheme could be envisaged in a strong field set-up. This would forgo the usual momentum information from the p x p y -plane in favour of a p ⊥ and l measurement. In theory this could be achieved in existing detection systems with the addition the hologram/phase mask and corrector as described in the OAM sorter [62]. This would map one of the dimensions in the p x p y -plane to the OAM, leading to different spots along the detector for different values of OAM. However, there are a great many practical considerations including the kinetic energy of the photoelectrons, propagation distance, detector and mask efficiencies to name a small few. In trying to combine these systems together it may be that they are incompatible without major alterations. However, the additional control and information that may be gleaned from strong field systems as a result certainly makes it worth consideration.
Because the OAM gives a direct measurement of the quantum magnetic number as well as providing sensitivity to chiral molecules [11,31], it is a useful tool for photo-induced time-resolved measurements. The first thought may be in exploiting OAM for such measurements through the incident light such as in twisted attosecond pulses. This is a nontrivial task but progress has been made in high-harmonic generation table top sources [63,64] as well in free-electron lasers [65]. However, the focus of this work is on the OAM of the outgoing photoelectron, which may be exploited via detection. As such there would be many possibilities for pump-probe schemes. For example the two counter-rotating circular fields employed in this work could each act as a probe, while two pumps separated by the same delay could be i) XUV attosecond pulses or ii) few-cycle IR fields. In the case of i) we have a situation reminiscent of the RABBITT (reconstruction of attosecond harmonic beating by interference of two-photon transitions) and attosecond streaking techniques [66,67]. In both cases, i) and ii), the idea would be to use VMI to exploit the interference, counting the number of spirals to determine the OAM. In order to apply such methods the interference vortices should be further studied in more complex targets to address the above-mentioned reservations for molecules. Direct measurement of the OAM could allow for any pump probe configuration as well as potentially removing the difficulties of using more complex targets. The ability to employ a variety pump-probe scheme with any ellipticity could allow for non-trivial coupling of the laser pulses to the OAM and novel time-resolved measurements of photo-induced dynamics.
In conclusion, we have presented a description of interference vortices in the strong field regime. Using the strong field approximation we have captured the main physical mechanism, deriving an interference condition, which closely matches calculations performed with the time-dependent Schrödinger equation-based solvers Qprop and R-Matrix with time dependence. We can explain the blurring of interference in the m = 1 case by the asymmetric yield from circular fields of opposite helicity, which is an exclusively strong field effect. By examining the system using orbital angular momentum (OAM) of twisted photoelectron states, we find a new interpretation of the electron vortices in terms of the interference of pairs of OAM states. Not only do we find good agreement between all models here but we also uncover a semiclassical relationship for the OAM by considering the interference between pairs of vortex states and linking it to the condition derived using the SFA. This has consequences, opening up the possibility to measure the OAM directly, or indirectly using interferometric schemes. The OAM of photoelectrons can open the possibility to a range of novel measurements and control in strong field systems.