Photon enhancement in a homogeneous axion dark matter background

We study the propagation of photons in a homogeneous axion dark matter background. When the axion decay into two photons is stimulated, the photon field exhibits a parametric instability in a small bandwidth centered on one half of the axion mass. We estimate analytically the enhancement for both coherent and non-coherent axion fields and we find that this effect could be relevant in the context of miniclusters and galactic halos.


Introduction
Several ideas beyond the standard model of elementary particles have been performed to explain the dark matter of the universe at a fundamental scale. One of the most serious candidates is the QCD axion which was originally proposed in order to explain the conservation of the discrete symmetries P and C P in strong interactions [1][2][3]. In this theoretical framework, the axion is a pseudo-Goldston boson arising from the spontaneous breaking of the PQ (Peccei-Quinn) symmetry defined at some energy scale f a . Below the QCD energy scale, the axion acquires its mass m a by nonperturbative effects. It is related to PQ symmetry breaking scale by m a 6 × 10 −6 eV 10 12 GeV As the axion mass, all the couplings with the standard model are inversely proportional to f a and since the PQ symmetry breaking is supposed to be very high ( f a ∼ 10 12 GeV), the axion is weakly coupled to the standard model particles.
Besides this feeble interaction, one of the most important features of axion physics is that they could be efficiently produced in the early universe by non-thermal mechanisms [4][5][6], such as the misalignment or from topological defect decays [7][8][9]. Under this context, they can be abundant and a e-mail: ariel.arza@ufl.edu cold enough to be an excellent dark matter candidate. On the other hand, pseudo-Nambu-Goldstone bosons emerging from the spontaneous breaking of global symmetries appear in several extensions of the Standard Model, and generically in all string compactifications [10]. They share many of the axion features, such as their small mass (also suppressed by the scale of the spontaneous symmetry breaking), their weak coupling to the standard model particles, and very importantly, they can also be produced in the early universe by the same non-thermal processes as axions so they are cold dark matter candidates too [11,12]. Due to all these similarities they are widely known as axion-like particles (ALPs), even though they do not feature an a priori relationship between the decay constant and their mass. One interesting attribute of axions (from now on we will refer as axions for both QCD axions and ALPs) from the experimental and observational point of view, is its coupling to two photons. The effective lagrangian of the axion-photon interaction is where a is the axion field and g its coupling to two photons. F μν is the electromagnetic strength tensor where A μ is the electromagnetic field.F μν is the dual of the electromagnetic strength tensor and is defined asF Axions have not been experimentally observed yet, so their existence is still in suspense, but because of their cosmological role as dark matter candidates, the experimental efforts to search for them in the corresponding parameter space have been intense. In the laboratory, axions are searched for using the coupling to two photons described by Eq. (2). Especially, static magnetic fields have been strongly implemented [13][14][15]. Some running experiments, such as haloscopes searches [16], already have access to the parameter space where they are cold dark matter candidates, and several new proposals plan to do so in the near future (ALPS-II [17], IAXO [18], HAYSTAC [19,20], MADMAX [21,22], among others).
In this paper we are interested in studying the effects on an incident electromagnetic plane wave traveling within a homogeneous axion dark matter background. We find an small window where the photon field experiences parametric resonance, we estimate the rate of the enhancement in coherent and non-coherent axion background and finally, we apply our results to astrophysical structures. Related topics were discussed recently in References [23][24][25], but they do not consider non-coherent axion fields. Studies about the propagation of photons in an axion dark matter background had already been performed a few years ago in order to account for their effects in optical experiments [26,27], however the instability was not considered.
The article is organized as follows: in Sect. 2 we obtain amplitude equations for states excited by the process a → 2γ and we analyze them numerically. Then, we solve the linear regime analytically finding the instability window and enhancement rate. In Sect. 3 we estimate the results obtained in Sect. 2 for coherent and non-coherent homogeneous axion fields. In Sect. 4 we discuss the possible astrophysical relevance in axion miniclusters and galactic halos.

Electromagnetic plane wave in a homogeneous axion background
The equations of motion for photons in an homogeneous axion background can be written in the Coulomb gauge as Let's consider an incident linearly polarized and plane electromagnetic wave with frequency ω as shown in Fig. 1. Inside the region filled with dark matter axions, the potential vector can be written as where the amplitudes A T (t) and A R (t) vary in time much slower than e −iωt . On the other hand, the axion field has the form a(t) = α(t)e −im a t + c.c., where the time dependent amplitude α(t) is introduced due to the feedback from the potential vector. The most interesting scenario is where the incident wave stimulates the process a → 2γ . In that case, the conservation of momentum and energy in the axion rest frame requires the condition ω = m a /2 to be satisfied.
Equations (3) and (4) become where we have neglected high frequency oscillation terms.
The complete behavior of this system is shown in Fig. 2, where we have chosen the initial conditions A T (0) = α(0)/2 and A R (0) = 0 (the initial wave travels forward). We can see that the system is unstable increasing the electromagnetic field considerably. In order to compute analytically this instability, we focus in the linear behavior (see dashed black line in Fig. 2) which corresponds to the limit of a constant axion energy density. In such a case, the axion amplitude α(t) is considered as a slow oscillating function, which can be written as where ρ is the axion energy density and corresponds to small deviations of the axion frequency. Defining A T = a T e − i 2 t and A R = a R e − i 2 t , Eqs. (6) and (7) leads to We can see that for there is a parametric resonance where a T and a R grow exponentially in time being their complete solutions where σ = g 2 Inside the instability window (11), the potential vector grows exponentially in time according to Eqs. (12) and (13), and the electromagnetic signals (e.g. the power) will be enhanced at the rate e γ t , where The signal is a spectral line centered at m a /2 with a bandwidth of ω = g √ ρ/2 [because of the parametric resonance window (11)] and amplitude where d P 0 /dω is the power spectrum of the incident beam.
To finish this section, we want to make some comments about the validity of our results in the linear approximation. To estimate the moment in which the exponential behavior is no longer valid, we apply ∂ t to Eq. (6). Using Eqs. (7) and (8) we find We can see that the instability works until a time in which After that, the right-hand side of (17) becomes negative and A T starts to oscillate. In the linear regime, after a sufficient time t, we can approximate |A R | 2 ≈ 1 4 |A T (0)| 2 e γ t and |α| 2 ≈ |α(0)| 2 . Condition (18) translates to Since the maximum enhancement is closely related with the validity of the linear regime, it must be proportional to A = |α(0)| 2 /|A T (0)| 2 . In Fig. 2 we use A = 4, however it could be much bigger considering that |α(0)| 2 is given by the dark matter energy density.

Enhancement in a coherent and non-coherent axion background
Consider a dispersive axion background with frequency distribution F(ω a ) and total energy density ρ T . This distribution is normalized as and characterized by a bandwidth ω a given by where δv is the velocity dispersion. On the other hand, we have seen in the previous section that Eq. (11) defines the bandwidth δω a = g 2ρ (22) in the axion frequency where axions produce parametric resonance. When δω a ∼ ω a or greater, the axion field is coherent enough to matching ρ = ρ T . In this case, Eq. (15) becomes However, when δω a ω a , the axion field is considered as non-coherent for our purposes, and only a small axion population will be engaged with the instability. The effective energy density of this population is given by Combining Eqs. (22) and (24), we find that the effective density is Replacing (25) into (14), and assuming that F(ω a ) ∼ ω −1 a , we have

Possible astrophysical relevance
In order to figure out the parameter space (m a , g) where this effect could be important, we focus in two different astrophysical scenarios characterized by a large dark matter energy density; axion miniclusters and galactic halos. As a criterion, we establish that all the dark matter is composed by axions (or ALPs) and that parametric resonance plays a significative role when where d corresponds to the distance traveled by the incident photon within the axion dark matter background, which is  (m a , g). The gray region is the space of parameters excluded by experiments (see Ref. [30]). Condition (27) is satisfied above the red and purple line for axion miniclusters and galactic halos assuming that the caustic ring model is correct. These lines were plotted using Eqs. (32) and (33) actually the maximum time during which the parametric resonance can occur. Using some typical physical parameters, we plot in Fig. 3 the regions in the parameter space where condition (27) is satisfied for both axion miniclusters and galactic halos. More details of these results will be explained in the following.

Axion miniclusters
Axion miniclusters have been well motivated if inflation occurs before PQ phase transition [28,29]. In this case, inhomogeneities of the axion fluid form bound objects by gravitational instabilities. The typical total energy densities of the axion miniclusters is about ρ mc ∼ 10 −18 g/cm 3 [9]. For axion miniclusters we assume that the axion background is coherent enough, for instance it could be considered as a Bose Einstein condensate [31,32]. In this case the exponential rate is just γ ∼ g √ ρ mc /2 [see Eq. (23)]. The typical size is where R(t) is the scale factor and t eq ∼ 2 × 10 12 s is the cosmological time of matter-radiation equality. t 1 corresponds to the cosmological time when the axion field starts to oscillate and is defined by m a (t 1 )t 1 ∼ 1. In the most general scenario (see [12]), the axion energy density today can be expressed as ρ a,0 0.58ρ CDM m a 10 −5 eV where ρ CDM is the today dark matter energy density. From Eq. (29) we can find t 1 assuming that all the dark matter are axions, we have Condition (27) gives us

Galactic halos
The caustic ring halo model [33] could be very interesting in this context because predicts regions in galaxies with large energy densities and very small velocity dispersion. For instance, it is estimated that the 5th caustic ring in the Milky Way, which is located very close to the sun, features an energy density of 1.5 ×10 −24 g/cm 3 and a velocity dispersion about 53 m/s or less in a extension of about 130 pc [34,35]. In addition, there is no reason to think that other caustic rings in galactic halos can not have larger densities and widths.
Considering reasonable parameters such as an energy density ρ h ∼ 1 GeV/cm 3 , a velocity dispersion δv h ∼ 10 −8 c (where c is the speed of light) and a width d h ∼ 1 kpc, condition (27) This expression was found using Eq. (26) which is more appropriate in this context where the axion field is noncoherent. As a final remark, It is very important to make clear that in this manuscript we are only considering homogeneous axion backgrounds. Inhomogeneities could change partially or drastically our results.