ALP production through non-linear Compton scattering in intense fields

We derive production yields for massive pseudo-scalar and scalar axion-like-particles (ALPs), through non-linear Compton scattering of an electron in the background of low- and high-intensity electromagnetic fields. In particular, we focus on electromagnetic fields from Gaussian plane wave laser pulses. A detailed study of the angular distributions and effects of the scalar and pseudo-scalar masses is presented. It is shown that ultra-relativistic seed electrons can be used to produce scalars and pseudo-scalars with masses up to the order of the electron mass. We briefly discuss future applications of this work towards lab-based searches for light beyond-the-Standard-Model particles.


Introduction
The existence of light exotic pseudo-scalar or scalar particles is highly motivated by some of the most widely studied extensions of the Standard Model (SM).For example, the strong-CP problem of QCD can be elegantly solved through the Peccei-Quinn mechanism [1] which predicts the existence of a new spin-0, CP-odd particle called the axion 1 .Many other extensions of the SM also contain spontaneously broken U (1) global symmetries and light spin-0 Goldstone bosons as a result.We broadly refer to these light spin-0 particles as axion-like-particles (ALPs) irrespective of their CP properties.The interactions of these particles are described by the Lagrangian densities where the ± refers to the even/odd CP property of the ALP, φ represents the ALP field, F µν is the electromagnetic field strength tensor, and F µν = 1 2 µνηδ F ηδ is the dual field strength tensor.A vast amount of work has been done studying the experimental signatures 2 of these fields in low-energy lab-based experiments (light-shining-through-wall (LSW) experiments) [5], solar experiments [6], dark matter and stellar evolution experiments [7], beam dumps [8], rare meson decays [9], and in high energy collider experiments [10].In the LSW experiments, and in most other ALP searches, the aim is to produce and detect ALPs through ALP-photon conversion mediated by the g φγγ coupling (for a review on theoretical aspects of LSW experiments see [11]).
Unique opportunities for observing beyond-the-SM (BSM) processes are also found through the study of quantum electrodynamics (QED) in high-intensity fields.The calculation of scattering matrix elements in high intensity fields cannot be performed using a perturbative expansion in the coupling and instead one must use non-perturbative solutions of the Dirac equation to describe the interaction between the electromagnetic field and the electron.The most popular of these solutions, of which very few are known, is the Volkov solution for an electron in a plane wave background [12].For reviews on these methods we refer the reader to [13,14] and we list recent developments in the field in [15].Given the increasing availability of high-intensity lasers from recent and upcoming experiments [16,17] the study of QED in intense fields to observe both SM and BSM phenomena is of crucial importance.So far, using high-intensity lasers to probe BSM phenomena has mainly be studied theoretically, and then through the coupling to ALPs [18] and mini-charged paticles [19] to photons.In this paper we study interactions involving electrons and ALPs in intense electromagnetic fields, focusing on the production of ALPs via Compton scattering from electrons in intense laser pulses.We will make some reasonable assumptions on the parameters in the calculation; the first being that the laser photons have an energy of κ 0 = 1.55eV (corresponding to a wavelength of 800nm), and the second being that the electrons can have momenta of up to a few GeV in optical set-ups.(In colliding 47 GeV electrons with a ps optical laser pulse, the SLAC E144 experiment [20] is an example of combining particle accelerator and laser pulse technology.)We will assume that a bunch of electrons interact incohorently with the external field, and for that reason restrict ourselves to processes involving single electron seeds.(Optical set-ups typically deal with bunches of the order of ∼ 10 8 electrons [17].) We begin the paper with a study of the Compton production of ALPs in a head-on collision between the seed electron and a low-intensity laser pulse.Due to the large number of photons, even in this low-intensity example, the laser pulse can still be treated as a classical background field and we expand the electron wavefunction perturbatively in a small intensity parameter, ξ 1.This dimensionless intensity parameter represents the work done by the external field over the Compton wavelength of an electron, in units of the external field photon energy and so in some way quantifies the number of photons interacting at a time with an electron.ξ will be defined quantitatively at the beginning of the next section.We assume the laser background to have a Gaussian pulse shape, however we also assume that the pulse duration is much larger than the photon wavelength, allowing us to approximate the electron as being in a monochromatic background.After obtaining analytical expressions for the production yield of the ALPs in both a linearly and circularly polarised laser pulse we study the total yield and angular distribution of the emitted ALPs for various ALP masses.We then move on to the study of electron-ALP interactions in high-intensity fields, i.e. ξ 1.Using the Volkov solution for the electron wave-function we take the limit of a constant-crossed background field and calculate the production yield of the ALPs via non-linear Compton scattering.A similar calculation for the emission of a massless pseudoscalar was performed in [21], where bounds on the ALP properties were derived using astrophysical constraints.Employing the Locally Constant Field Approximation (LCFA), see for example [22,23], we use this result to approximate the production yield of ALPs when the background electromagnetic field has a non-trivial profile -such as a Gaussian or that of a focussed laser pulse.Using these solutions, we present a detailed analysis of the energy and angular distribution of the production yield for the ALPs.We study the effects of having a non-zero mass term for the ALPs and perform a comparison between the properties of scalar and pseudo-scalar production.Finally we conclude and discuss how this work can be used in studies of lab-based searches for light ALPs which probe the g φe coupling.

ALP production in a low-intensity laser pulse
In this section we study the Compton production of a scalar or pseudo-scalar from an electron in a low intensity external electromagnetic field.The external electromagnetic field is parametrised by where ξ = eA 0 /m e is the dimensionless intensity parameter, and f (x) is a pulse shape which describes the spatial dependence of the vector potential (the parameter ξ can be defined in a gauge-and Lorentzinvariant manner using the stress-energy tensor, see [24]).The low-intensity regime is then defined by the intensity parameter satisfying ξ 1.We choose a linearly polarised external field in the µ = (0, 1, 0, 0) µ direction, and label the perpendicular polarisation as ˜ µ = (0, 0, 1, 0) µ .The photon momenta are described by κ µ = κ 0 (1, 0, 0, 1) µ , and we define a dimensionless phase ϕ = κ • x which we use to parametrise the position of the electron wavefunction with respect to the external field.In the low intensity regime, the effects of this field on the electron wavefunction can be treated perturbatively, i.e. we expand the wavefunction to first order in a µ where the zeroth order part contains the free electron wavefunction and the first order part contains the interaction.Thus we have ψ p,r =ψ (0)  p,r + ψ (1)   p,r where the G(x−y) = 0|T ψ(x)ψ(y)|0 = e +iε e −ip•(x−y) is the fermionic propagator.This method approximates that only one photon is absorbed by the electron prior to the emission of the pseudo-scalar or scalar particle.The matrix element for the process e − → φ + e − can be written as where p, p , and k are the momenta of the incoming and outgoing electrons and the outgoing scalar or pseudo-scalar, respectively.The r and r label the spinor indices of the incoming and outgoing electrons.Then p,r , and (2.4) The matrix element for scalar production is obtained by replacing γ 5 by the spinor identity matrix as indicated by the structure of the interaction in Eq. 1.1.The outgoing wavefunction of the scalar or pseudo-scalar field is The probability can then be written in the form where T ± contains the traces over the spinor indices for either the scalar or the pseudo-scalar interaction.Note here that we have used the lightfront coordinates for the particle momenta, a description of which can be found in Appendix A. The spatial dependence of the external field now enters through the Fourier transform of the pulse shape, f , and its argument follows from momentum conservation imposed on the momentum of the recoiling electron, where the last expression is simply the on-shell condition.

The pulse shape and the monochromatic limit
To obtain the expression in Eq. 2.6 we began by Fourier transforming the profile f (x) as where κ µ = κ 0 (1, 0, 0, 1) µ describes a plane wave trajectory for the photon with κ 0 being the photon energy.We suppose that the pulse shape for f (x) is Gaussian with respect to the phase ϕ, i.e.
2 −iϕ (2.9) where Φ = κ • τ is a pulse duration with τ µ = τ 0 (1, 0, 0, 0) µ in the lab frame, and the terms linear in ϕ in the exponent describe the oscillations of the plane wave with frequency κ 0 .From this we can calculate f (2.10) and insert it into Eq.2.6.This results in a complicated expression which can be simplified by assuming that Φ 1.Using this we arrive at the monochromatic (or long-pulse) limit in Eq. 2.6, the long pulse limit implies that the incoming electron absorbs or emits one photon of fixed energy to or from the external field.After these manipulations the probability can be written as 2.13) where one of the integrals over d 2 k ⊥ dk − will be used to enforce the delta function condition.

Simplifying the expressions for the yields
The function T ± can be written as where ± corresponds to scalar and pseudo-scalar, respectively.Before writing these factors it is useful to note that momentum conservation implies where we recall that r is the variable from the Fourier transformation of f (ϕ) defined in Eq. 2.12.
Using these relations we can write the factors in the trace as We see that when the electrons collide head-on with the laser pulse i.e. ( • p) = (˜ • p) = 0, the expressions simplify greatly and the scalar and pseudo-scalar kinematics differ only in the γ ± terms.
From here onwards, we focus solely on the case in which ( •p) = (˜ •p) = 0.The next step is to perform the k integrals.The first point to note is that the integrations force r = 1, because the negative r = −1 solution is kinematically forbidden.The second point to note is that we can do the integrals either in d 2 k ⊥ dk − or in d 2 k ⊥ dk 3 .And lastly the third point to note is that the integrand is independent of ˜ • k, apart from through r which is a function of ( We define • k = m e ρ cos φ and ˜ • k = m e ρ sin φ with 0 ≤ φ ≤ 2π and ρ ≥ 0, such that the delta function simplifies to .18)We can now write the probability as To simplify the expressions we define k − = vp − , m φ = δm e , and κ • p = η p m 2 e .Performing the ρ and φ integrals we have (2.20) Note that we require g > 0 to obtain real solutions, putting a limit on v which can be written as (2.21)This ensures v < 1 and thus the expressions for T ± never become singular.In the case where we use a circularly polarised external field we find that the result is the exact same.This is only true however in the limit where • p = ˜ • p = 0.As mentioned in the introduction, κ 0 = 1.55eV and |p| O (10)GeV, from which it follows that η p 1. (This can be seen from

Angular distributions and effects of a non-zero ALP mass
In Figure 2 we plot the total yield for pseudo-scalar and scalar production as a function of the ALP mass for various seed electron momenta in the MeV range 3 .We refer to the quantity P as the production yield as it represents the number of ALPs expected to be emitted while the electron is in the external field.When dealing with the external field P can obviously not represent a probability since if Φ is large enough P can be larger than 1.In a realistic experimental set-up the interaction would take place between a laser pulse and a bunch of approximately 10 8 electrons.We can see from Figure 2 that the range of ALP masses probed by this interaction depends entirely on the energy of the initial seed electron, with the photon energy, κ 0 , being fixed to 1.55 eV.We can see that the production yield cuts off quite sharply as δ reaches a critical value that depends on the seed electron energy.With larger | p| values the cut-off increases approximately linearly, and | p| ∼ O(GeV) allows the interaction to probe m φ ∼ 2 × 10 4 κ 0 ∼ O(m e ) masses.It is useful to look at the scattering in lab-frame polar coordinates (see Appendix A.1 for details on the co-ordinate transformation).The exact expressions of the yield in these co-ordinates are lengthy and we refrain from presenting them here, although the reader can easily deduce them from the information already provided.In the monochromatic limit there is a direct relationship between the ALPs' energy and the polar angle at which they are emitted.This is irrespective of the scalar or pseudo-scalar nature of the particle and we have plotted the relationship for various seed electron energies in Figure 3. From this plot we can see that for | k| κ 0 the polar angle of emission is close to zero, which corresponds to emission parallel to the momentum of the incoming laser photon.For larger ALP momentum the polar angle shifts closer to π, which corresponds to emission parallel to the momentum of the incoming electron.As the seed electron energy increases the polar distribution becomes more sharply localised towards the polar angle of the incoming electron4 .Permitting the ALP to have a non-zero mass drastically alters the properties of the emission, as can be seen in Figure 4.In particular, a non-zero mass alters the polar angle at which the ALP is emitted and restricts it to lay closer to the polar angle along which the initial seed electron travels.Note that with κ 0 = 1.55 eV we only reach η p = 1 when | p| ∼ 40 GeV.We have not plotted the | k| distribution of the emitted ALPs for the simple reason that the distributions are approximately constant over the ranges depicted in Figure 3, a feature which persists even for seed electron energies of O(GeV).

ALP production in a constant-crossed field
If the intensity of the external field is such that ξ O(1), the approximation that the electron absorbs only one photon before emitting the pseudo-scalar or scalar particle breaks down.All orders of photon exchange with the electron must be included for a calculation to be consistent.In this section we calculate the production yield for pseudo-scalar and scalar particles via non-linear Compton scattering in a constant-crossed external field.This result is particularly important as in the rest frame of an ultra-relativistic particle all electromagnetic fields resemble a constant-crossed field [13].
We will use the notation and structures of the previous section, and introduce new ideas along the way.The relevant non-linear parameter that we use in the study of high-intensity fields is χ q = ξη q , sometimes referred to as the quantum nonlinearity parameter, which is equal to the work done by the external field over a Compton wavelength in units of the electron rest energy.Assuming κ 0 = 1.55 e and ξ ∼ O(100) we have χ q η q and to generate an O(1) χ p for the seed electron we only require | p| ∼ 100 ξ GeV.In comparison, for non-relativistic electrons with | p| O(MeV) and the same laser parameters, χ ∼ 6 × 10 −3 .
To calculate the scattering matrix when ξ ∼ O(1) or larger we must use the non-perturbative Volkov solution of the Dirac equation for an electron in a plane-wave background, The dynamics of the background field is described by the ¡ κ ¡ a and the S p terms, which are non-linear in the vector potential: and a µ (z) is the external plane-wave EM field defined as with κ µ = κ 0 (1, 0, 0, 1) µ being the photon momentum vector.Despite using a plane-wave solution of the Dirac equation the constant field limit can be taken in integrated expressions through κ 0 → 0. We will see that the result for the total yield in the constant field is independent of κ 0 and the limit is trivial.We will start with the calculation for pseudo-scalar production and then present the result for scalar production, this allows us to describe the calculation in more detail.The matrix element is written in a similar way: The Γ − (ϕ) function arises from the Volkov solution and contains all the spinor and external field dependence Performing similar steps to the previous section and transforming to light-front coordinates we find that the total probability is given by originating from the Fourier transform of Γ † − Γ − .The function T (p, k, p , ϕ, ϕ ) contains the spinor traces that arise after taking the spin sum.The α and β functions from the above expression are given by The trace element of the spin sum is contained within the T (p, k, p , ϕ, ϕ ) function, which can be simplified using the momentum conservation relations in Eq. 2.15 to find The spatial integrals can be computed exactly in the constant field limit, corresponding to g(x) = x.
To perform the ϕ and ϕ integrals a change of variables is useful, and we choose The integrals can be performed exactly using the integral identities in appendix B and the total probability can be written as where we recall δ = m φ /m e .The above result is independent of • k.This is precisely because we have chosen the constant-crossed field background, which we will discuss in more detail in the next section.Using the identities listed in Appendix B the ˜ • k integral can be performed exactly and we find where The calculation for the production of a scalar proceeds analogously, apart from the γ 5 operator is replaced with the spinor identity matrix in the interaction between the φ and fermion fields.The fermion trace for the scalar field production is equal to that in the pseudo-scalar trace plus a factor of 8m 2 e , and the terms in the exponent of the integrand in the spin sum remain unchanged.Integrating the (ϕ + , ϕ − ) variables we find that the total probability can be written as where the form of ∆ and Σ in the χ q notation can be found in Eq. 3.12.Integrating over d(˜ • k) we then find, It is important to notice here that the argument of the Airy function is the same as in the pseudo-scalar case, which follows from the kinematics of the collision.

ALP production in high-intensity backgrounds
In the constant-crossed field calculation there appears a divergent integral in • k.However, this can be reinterpreted in the following way.The integral over the electron's phase co-ordinate ϕ, performed at the amplitude level, leads to the Airy functions at the level of the probability.The contribution from the Airy functions occurs mainly when they have an argument of the order of unity or less.This corresponds to a finite region of the electron's trajectory in ϕ.This so-called "coherence interval" [13] becomes ever smaller as ξ increases.In the limit ξ → ∞, the relevant part of the electron trajectory corresponds to the stationary phase: Therefore, there is a one-to-one mapping between the electron's stationary phase (representing its classical trajectory) and the value of • k at which an ALP is emitted.This means the divergent integral in • k can be reinterpreted as an integration over the electron's phase ϕ * as it propagates through the background.Writing the probabilities for the scalar and pseudo-scalar emissions as a probability per unit phase using d( where the difference in the two lies in the pre-factor of the Ai 1 term.To calculate the probability of emission in a non-trivial external field we then use the Locally Constant Field Approximation (LCFA) and make the replacement where g (ϕ * ) is the profile of the electric field.
In Appendix A.1 we discuss the transformation from d( • k)dk − to polar coordinates and the consequences for light-front invariants.The same transformation can also be used here to obtain information on the angular distribution of the emitted ALPs in a non-trivial external field background.Using Eq. 4.1 we can write .4)where g ( • k) is simply the external electric field profile written in terms of the ALP momenta and the ellipsis corresponds to the integrands in Eq. 4.2 written with the replacement χ p,k → χ p,k (ϕ * ).Note that the Eq.4.4 does not explicitly depend on the non-linearity parameter ξ as this only enters through the χ-parameters, also the ratios χ k χp = k − p − remain independent of both ξ and g(ϕ * ).

Yield distributions in a constant field background
We can study properties of the yield distribution for scalars and pseudo-scalars in a constant field by taking the external field profile to be constant over some finite distance, i.e.
with Θ(x) being the Heaviside step function and L being some finite phase.The yield distribution in χ k remains constant with ϕ * thus we can sample the distribution at one point to examine its behaviour, this is shown in Figure 5 where we have set g φe = 1.In this figure we display the yield distributions in χ k for various seed electron energies, where the probabilities have been re-scaled for purposes of comparison.There is a clear difference between the scalar and pseudo-scalar χ k distributions.
The distribution for pseudo-scalar production is peaked away from zero for all values of χ p,0 , with the maximum of the distribution moving closer to χ k,0 = χ p,0 for larger χ p,0 .For χ p,0 1 the distribution in χ k for scalar production is peaked at zero, whereas for larger values of χ p,0 the distribution becomes peaked away from zero and begins to resemble the distribution for pseudo-scalar production.

Yield distributions in a Gaussian background
We now study the emission of pseudo-scalars and scalars from an electron in an external field with a Gaussian profile described by with Φ being the pulse duration in units of inverse κ 0 .We choose the duration for the high-intensity laser pulse to be 100 fs throughout this chapter (corresponding to Φ ≈ 300).The distribution in χ k,0 now has a non-trivial dependence on the phase ϕ * .To show this we have plotted the pseudo-scalar and scalar yields for various seed electron energies in Figures 6 and 7, respectively, where we have set g φe = 1.In the pseudo-scalar case we see that the distribution is localised at a point which for low seed electron energies is at χ k,0 χ p,0 , but for larger seed electron energies moves towards χ k,0 χ p,0 .This is in contrast to the scalar case in which the distribution is always localised around χ k,0 0, where larger seed electron energies increase the spread of the distribution towards χ k,0 χ p,0 .Similar behaviour can also be seen in the constant field case, Figure 5.Note that these distributions are symmetric around ϕ * = 0, which corresponds to the central peak of the Gaussian profile in g (ϕ * ).In Figures 8 and 9 we show, using Eq.A.8, how this translates to distributions in the polar angle and energy of the emitted ALP in lab frame coordinates.For both the pseudo-scalar and scalar yields the higher energy particles are emitted at smaller polar angles for small seed electron energies.Larger seed electron energies results in the ALPs being emitted at larger angles, becoming parallel with the direction of propagation of incoming electrons for very large seed electron energies.These distributions are similar to the constant-crossed field case in that the scalar yield is peaked at | k| 0 while the pseudo-scalar yield is peaked at a non-zero | k| determined by the seed electron energy.

Total yields and ALP mass dependence
Thus far in this section we have taken the ALP to be massless, i.e. δ = 0.For the mass effects to become significant in the high-intensity regime we require δ 0.1.This contrasts with the results we found in the low intensity external field in Section 2, where the range of masses probed by the interaction was limited by the energy of the photons in the background field.In this section we will study what happens to the differential yield and the total yield when we allow for sizeable ALP masses.
We start by plotting the total integrated yield for a Gaussian background field against the seed electron χ p,0 for various values of δ in Figure 10.The analogous plot for a constant field background has the same features but at a different magnitude, scaled by how long the electron is taken to interact with the constant field (which is of course, formally infinite).In Figure 10 we have allowed for a much larger range in χ p,0 than in previous plots, this is only done because the full picture of the effects related to a non-zero ALP mass on the total yield only become apparent at these larger values of χ p,0 .The sharp peak for the massless scalar can be resolved in the logarithmic plot in Figure 11 below.We can see that there is a steep exponential increase in the production rate as χ p,0 is increased to χ p,0 ∼ 0.3 followed by a gentler exponential decrease.It is also interesting to see the effects of δ = 0 on the χ k distribution of the emitted ALP in a constant field, we have plotted this in Figure 12 where χ p,0 has been fixed to 1.  Lastly we show in Figure 13 the effects of a non-zero δ on the angular distribution of emitted ALPs in a Gaussian background with χ p,0 = 1.In all of these cases we see that for larger values of the ALP mass the properties of scalar and pseudo-scalar emission become very alike.This is hinted at by the expressions in Eq. 4.2 where we see that the difference between the analytic formulas for pseudo-scalar and scalar emission lies in the δ-dependent pre-factor of the Ai 1 function.

ALP production in a high-intensity laser pulse
The external field backgrounds studied in the previous subsections serve as a useful test case for the physics of ALP production in more general external fields.However in a laser-based experimental set-up the strong electromagnetic field will have a carrier wave frequency.Assuming that the photons are linearly polarised, an example of this is to take where we have assumed a Gaussian pulse shape.In using this pulse shape, we expect the results to be similar to that in the previous section where we simply had a Gaussian profile.The main difference in a laser pulse background is that the cosine modulation results in a modulation of the differential yield in ϕ * .We can see an example of this in Figure 14 where we plot the χ k distribution for pseudo-scalar emission with χ p,0 = 1 and m φ = 0.For Figure 14 we choose a shorter pulse duration of 10 fs such that the modulation effects are more apparent, as for longer pulses the wavelength in units of ϕ * becomes very small in comparison the the pulse duration in units of ϕ * .We can see that in the case of the laser pulse background the overall shape of the contours matches that in the case of just a Gaussian background.It is also important to note that the peak values in the differential yield are also the same.This tells us that the total yield in both cases will be similar, and in fact the total yield in the case of a laser pulse only differs from that in the Gaussian case by a numerical factor ∼ 50%.The effects of a non-zero mass and varying seed electron momentum mirrors those in Figure 10 and 11 for the Gaussian case.

Conclusions, analysis, and outlook
In the interest of future lab-based experimental ALP searches a detailed study of ALP production via Compton scattering in low-and high-intensity electromagnetic fields, with a particular focus on laser pulses, was performed.Particular properties of the production yields depend strongly on the CP nature of the ALP, through which its coupling to electrons is determined.The basic set-up that we envisage is an electron colliding almost head-on with a laser pulse.For optical lasers, photons have energies of the order of eV, and in an optical set-up, the electron momentum could be anywhere from keV to several GeV.
In Section 2 the production yields for ALPs were derived in the case of a low-intensity laser pulse interacting with an electron.That is, larger ALP masses suppress the production yield, with the cutoff on the ALP mass imposed by having a non-zero production yield being largely determined by the energy of the photons in the background field.The angular distributions in Section 2.3 show that the ALPs with larger energies are emitted in the direction in which the incoming seed electron propagates, whereas lower energy ALPs are emitted off-axes.We can estimate the number of axions emitted in an interaction between a low-intensity laser pulse and an electron bunch as where W parametrises the kinematical and axion mass dependencies, N e × 10 8 is the number of electrons per bunch, and Φ = Φ × 10 3 .Typical expectations from a lab-based experimental set-up are that N e and Φ would be O(1) numbers.For scalar production we find that W 0.066 for light ALPs regardless of what the electron momentum is, however this number begins to decrease for larger ALP masses.To maintain W 0.066 for larger axion masses, up to ∼ 0.1m e , it is only required to increase the momentum of the electrons interacting with the pulse.However, for m φ > m e , W becomes suppressed.For pseudo-scalars the situation is much different.For light ALP masses W strongly depends on the momenta of the incoming electrons, where for electrons at rest W ∼ 10 −12 .For ultra-relativistic electrons however (| p| ∼ 10 4 m e ) the W factor for the pseudo-scalars becomes similar to that for the scalars.The behaviour of W for the ALPs can be determined from Figure 2.
In Section 3 the production yields for ALPs in a high-intensity constant-crossed electromagnetic field were derived.In Section 4 the LCFA was used to translate this result to an approximation for the production yield in non-trivial field configurations, such as a laser pulse.The sensitivity of the yield to the CP nature of the ALP and the mass of the ALP was studied in detail in  It was found that the CP-even states naturally have a larger production yield than the CP-odd states, however as the axion mass is increased the total production yields become qualitatively similar both in magnitude and in their sensitivity to the seed electron energy.An increase in the ALP mass always leads to a reduction in the production yield, however this is not as drastic as in the low-intensity regime.
Here we find that the production yield remains sizeable even for ALPs with masses greater than that of the electron.The reason for this lies in the fact that the electron can absorb many photons before emitting an ALP, thus increasing the available energy for ALP production.Also, regardless of the ALP mass both production yields become similar in magnitude for larger seed electron energies.This is quite similar to the behaviour seen in the low-intensity case.The most interesting effect comes from the angular and momentum distributions of the emitted ALPs -using a Gaussian pulse for the background field, it was demonstrated that the angular distributions were strongly dependent on the CP nature of the ALP.We found that CP-odd ALPs have a momentum distribution peaked at some non-zero value determined by the seed electron energy, and the CP-even ALPs have a momentum distribution peaked at zero.However, allowing the ALP mass to increase to O(m e ) we find that the angular and momentum distributions for the CP-even and -odd ALPs become virtually indistinguishable, differing only in magnitude.Due to the non-trivial dependence of the production yield on ξ and Φ in a highintensity laser background we cannot factorise these as we have in Eq 5.1 for the low intensity laser pulse.However, as can be seen from the previous section the production yields are typically at least an order of magnitude larger than in the case of a low-intensity laser pulse, with a larger range of accessible ALP masses.
These results constitute the first analysis of spin-0 particle production via non-linear Compton scattering in intense laser pulses.In this paper we have not only described how theoretical predictions for such a process are calculated, but we have already obtained information on the characteristics of the production mechanism which will be crucial to understanding how to detect these ALPs in a lab-based experimental set-up.

Figure 1 :
Figure 1: A schematic diagram of the interaction taking place between the laser background and electron, leading to a scattered electron and ALP.

Figure 2 :
Figure 2: The total yield is plotted against the mass of the emitted ALP for various seed electron energies.

Figure 3 :
Figure 3: Relationship between the energy of the emitted ALP and the polar angle at which it is emitted in lab frame.These plots are for m φ = 0.

Figure 4 :
Figure 4: Relationship between the energy of the emitted ALP and the polar angle at which it is emitted in lab frame for | p| = me and various ALP masses.

Figure 5 :
Figure 5: The relationship between the differential yield and the χ k parameter is plotted for different seed electron energies with ξ0 = 100 and g φe = 1 5 .

Figure 6 :Figure 7 :
Figure 6: The pseudo-scalar yield, P−, is plotted as a function of χ k /χp and the phase ϕ * in a Gaussian background with a pulse duration of 100 fs and m φ = 0.

Figure 8 :Figure 9 :
Figure 8: The pseudo-scalar yield, P−, is plotted as a function of the angle and momenta of the emitted ALP in a Gaussian background with m φ = 0.

Figure 10 :
Figure 10: The pseudo-scalar and scalar yields, P±, are plotted as a function of χp,0 for various ALP masses, with m φ = δme, in a Gaussian background.

Figure 11 :
Figure 11: The scalar yield, P+ is plotted as a function of χp,0 for δ = 0 in a Gaussian background.

Figure 12 :
Figure 12: The pseudo-scalar (left) and scalar (right) yields, P±, are plotted as a function of the χ k /χp and the phase ϕ * in a Gaussian background for two different ALP masses.Also, χp,0 has been fixed to 1.

Figure 13 :
Figure 13: The pseudo-scalar (left) and scalar (right) yields, P±, are plotted as a function of the polar angle and momenta of the emitted ALP in a Gaussian background for two different ALP masses.Also, χp,0 has been fixed to 1.

Figure 14 :
Figure 14: The pseudo-scalar yield, P−, is plotted as a function of ϕ * and χ k /χp for δ = 0 in (a) a Gaussian pulse background and (b) a laser pulse background, with a pulse duration of 10 fs and κ 0 = 1.55 eV.