Solar radio emissions and ultralight dark matter

Ultralight axions and dark photons are well-motivated dark matter candidates. Inside the plasma, once the mass of ultralight dark matter candidates equals the plasma frequency, they can resonantly convert into electromagnetic waves, due to the coupling between the ultralight dark matter particles and the standard model photons. The converted electromagnetic waves are monochromatic. In this article, we review the development of using radio detectors to search for ultralight dark matter conversions in the solar corona and solar wind plasma.


Introduction
There is overwhelming evidence from astronomical and cosmological observations supporting that about one-fourth of the energy density in today's universe is composed of dark matter. However, the particle nature of dark matter is still mysterious. Because of the null results of searching for WIMPs [1][2][3], more and more attention has been attracted to ultralight dark matter candidates such as the QCD axion [4], axion-like particles [5], and dark photons [6][7][8][9]. Originally, the theory of QCD axion was put forward to resolve the strong CP problem via the known Peccei-Quinn mechanism [10][11][12][13]. Soon after that, it was shown that, as a bonus, QCD axions are suitable dark matter candidates [4]. In addition, many models such as the string theory will give rise to axion-like particles which have similar features as QCD axions in coupling with the standard model of particle physics. Similarly, axion-like particles can also serve as dark matter candidates. (In the rest of the paper, we will call both the QCD axion and axion-like particles the axions). There are multiple ways in the early Universe of generating axions: both the misalignment mechanism [14][15][16] and the decay of axionic topological defects such as domain walls and strings [17,18] can produce an appropriate amount of axions as dark matter. On the other hand, dark photon is usually a U(1) gauge boson which kinetically mixes with the standard model photons with a tiny coupling strength [19][20][21][22][23][24]. Dark photons are also good dark matter candidates [6][7][8][9] and can be generated in multiple ways such as parametric resonances [25][26][27][28][29][30], inflationary fluctuations [9,31,32,32,33,[33][34][35][36][37][38][39][40], misalignment mechanism [7,8,[41][42][43], and cosmic strings decay [44]. Current measurements show that the local dark matter density in the solar system, ρ DM , is about 0.3 GeV/cm 3 [45,46]. In this paper, we assume that the dark matter is saturated either by axions or dark photons.
The interaction between axions and electromagnetic fields makes it possible for the axion dark matter to convert to electromagnetic waves in a magnetic field. If plasma exists, the interaction between the electromagnetic field and the charged particles in the plasma can induce an effective mass for transverse photon, which equals the plasma frequency ω p . Thus, in the region that ω p = m DM , the axion DM will resonantly convert to electromagnetic waves. For dark photon dark matter, the plasma will also regenerate the mixing of the dark photon arXiv:2304.01056v1 [hep-ph] 3 Apr 2023 field and the photon field. Thus, the once diagonalized fields can convert to each other again inside the plasma. Therefore the dark photon dark matter can also convert into electromagnetic waves in the region ω p = m DM . The energy of the converted photon equals the energy of the dark matter particle, which can be written as m DM + E kDM , where E kDM is the kinetic energy of each dark matter particle. In the galactic halo, we have E kDM ∼ 10 −6 m DM . Thus, we expect a monochromatic electromagnetic wave signal produced by the ultralight dark matter conversions. To take advantage of the resonant conversion to search for ultralight dark matter, we need first to find a large plasma. Secondly, the region between the conversion region and the detector must be transparent for the converted electromagnetic wave to propagate to the detector. A perfect plasma that we can use is the Sun's corona. The distribution of the temperature and plasma density in the Sun's corona is shown in Fig. 1a, where one can calculate that ω p ranges from about 10 −6 eV to 10 −7 eV. Thus, if the mass of axion dark matter or dark photon dark matter falls in this region, it can produce monochromatic electromagnetic wave signals. The frequency in this region is about 10 MHz − 1 GHz, which is in the radio wave band. There is already a lot of radio observation apparatus toward the Sun, such as the Low Frequency Array (LOFAR) [47], having collected plenty of solar radio data, in which we can search for ultralight dark matter conversion signals. There will also be future projects, such as the Square Kilometer Array (SKA) [48] and the upgraded Five Hundred Meter Aperture Spherical Radio Telescope (FAST), which may provide more precise solar data in the future.
The earth's ionosphere also provides the plasma for the ultralight dark matter conversion. However, the frequency range is about 4-6 MHz, which is polluted by human activities, and therefore cannot be used to search ultralight dark matter. Also, because of the screening effect of the earth's ionosphere, there are no terrestrial radio telescopes with frequency ranges below about 10 MHz. Thus, if the frequency of the ultralight dark matter (m DM /(2π)) is smaller than 10 MHz, we need to go outer space to search for the conversion signal. The reason the conversion can still happen is that solar wind provides charged particles in the inter-region between the Sun and the Earth. In Fig. 1b, one can see that ω p in the region between the Sun and the earth is all the way down to 10 −10 eV. Solar probe satellites with long-wave radio spectrometers, such as the STEREO satellites [49] and the Parker Solar Probe [50], have been sent out to observe solar activities, and have collected a lot of data. Thus, one can search for the conversion signal in these observation data.
Similar ideas have been proposed to search for radio waves from axion resonant conversion in the mega magnetic field around the neutron stars [51][52][53] and white dwarfs [54][55][56][57]. On the other hand, solar radio emissions are also related to other dark matter models. For instance, the impulse radio events observed by the Murchison Widefield Array (MWA) [58] may be explained by the annihilation of axion-quark-nugget dark matter [59][60][61][62][63] in the solar corona, the energy released from which may account for the solar corona heating mystery [64][65][66].
In this article, we review the recent developments in searching for ultralight dark matter in solar radio observation. In Sec. 2, we review particle physics models for axions and dark photons, focusing on their interactions with the electromagnetic wave. In Sec. 3, we present the calculation of the conversion probabilities for axions and dark photons. In Sec. 4, we review the calculation of the propagation of the converted electromagnetic waves. We discuss the sensitivities of the radio telescopes for such resonant signals in Sec. 5. We summarize the results in Sec. 6.

Ultralight dark matter
Axion is a pseudoscalar particle whose interaction with the electromagnetic part can be written as The n e profile is based on the in situ measurements of the Parker Solar Probe [69] and then fitted by the relation n e (r) ∼ 3.3 × 10 5 (r/R ) −2 + 4.1 × 10 6 (r/R ) −4 + 8.0 × 10 7 (r/R ) −6 cm −3 given in [70]. The T profile is the relation T(r) ∼ 418 eV × (r/R ) −0.74 [69] based on the in situ measurements. a is the axion field and m a is its mass. F µν = ∂ µ A ν − ∂ ν A µ is the field strength of the standard model photon andF µν = ε µναβ /2 · F αβ is the dual field strength. g aγγ is the coupling between axion and standard model photon. The interaction term in (1) can be further simplified to −g aγγ aE · B in terms of the electric and magnetic fields.
Dark photon is a U(1) gauge boson that, after electroweak symmetry breaking, interacts with the standard model photon via the following Lagrangian: is the field strength of dark photon and m A is its mass. is the kinetic mixing strength between dark photon and photon. To study more conveniently the conversion between dark photon and photon, we can diagonalize (2) in the interaction basis with the following field redefinition: and (2) becomes In solar plasma, photons will acquire an effective mass equal to the plasma frequency, m A = ω p . The plasma is charge neutral consisting of electrons and ions. Ignoring the contributions of ions due to their large masses, we only care about the plasma frequency of electrons, which in the non-relativistic case is α is the fine structure constant, m e is the electron mass, and n e is the electron density. Therefore, for photons in plasma, we should add an effective mass term L ⊃ 1/2 · ω 2 p A 2 in the Lagrangians (1) and (4).

Conversion in solar plasma
The equations of motions corresponding to the Lagrangians (1) and (4) for the axion and dark photon cases can both be written in the form [71] where X denotes either axion a or dark photon A . For the axion case, we havẽ and for the dark photon case, we havẽ Since the external environment in the solar plasma mostly depends on the radial direction r (this is true for the dominant term ∆ AA which is determined by electron density according to (5)), we only keep the spatial derivative in the radial direction in the equation of motion (6) while ignoring trivial dynamics in the other two spatial dimensions. Some differences between the axion case and the dark photon case are as follows. In comparison with the dark photon case, the conversion between axion and photon requires the presence of magnetic field B T transverse to the axion trajectory. The solar magnetic field provides the source of B T . Due to the pseudoscalar feature of axion, the converted photons A in (6) is polarized in the direction of B T [56,71]. For the dark photon case, (6) describes converted photons in two transverse modes perpendicular to the propagation trajectory.
We assume the solutions to the equation of motion (6) are in the wave form: where ω is the energy of wave, k = ω 2 − m 2 X is the momentum, and k r is the momentum component in the radial direction. For resonant conversion to happen, we should have m A = m X such that both ω and k of A and X match before and after the conversion. In the case we consider here where the dark matter velocity is non-relativistic v ∼ 10 −3 c, we have ω m X . Then, near the resonant layer in solar plasma, we have the following relation [71]: The last approximated relation holds because the parameters of solar plasma, e.g., plasma frequency, vary slowly in a length scale much larger than the wavelength 1/k r such that Therefore, using the relation (10), the equation of motion (6) can be simplified to The off-diagonal matrix M 1 is suppressed by the coupling constant, either or g aγγ , and thus can be treated as a perturbative term. If we start with a X particle, (A, X) = (0, 1) initially, then solving the equation of motion (12) to first order of M 1 , the probability that X eventually converts to a photon A, namely, (A, X) = (1, 0), is 6 of 18 We can solve the integration using the saddle-point method where we have used the expression of norm to remove the redundant phase factor. r c is the position when h(r) reaches its local maximum/minimum h (r c ) = 0, that is, the position of resonant X-A conversion when ω p (r c ) = m X . Only the second derivative plays an important role, the expression of which is The length of resonant layer δr c can be estimated as the length that makes h(r) ≈ h(r c ) The resonant layer is very thin, δr c /r c ∼ 10 −6 , for the typical values r c = R , ω p = 2π · 100 MHz, and v r = 10 −3 c.
Applying the saddle-point approximation (15), we solve (14) and obtain Then, for each case we have and Note that we have an extra prefactor 2/3 in the dark photon case. This is because the converted photons polarize in three directions with equal probability but the longitudinal mode cannot propagate out of the plasma. In comparison, the photons converted from axions polarize only in one transverse mode parallel to B T as mentioned above. Furthermore, the conversion probabilities for the two cases are related via the following replacement This relation is useful when we know the conversion rate of one case and directly calculate that of the other case.

Propagation of the converted photons
After the radio photons are converted from axion or dark photon, they will be refracted, absorbed, and scattered by thermal electrons in the plasma. We discuss these effects during the photon propagation in the following.

Refraction effect-
The behavior of converted photons propagating in the solar plasma obeys the rule of refraction: n(r c ) n(r) = sin θ(r) sin θ(r c ) . (22) θ(r) is the angle of incidence or refraction with respect to the radial direction. n(r) is the refraction index as a function of distance, which is equal to the inverse of phase velocity: At the position of resonant conversion, n(r c ) v ∼ 10 −3 . Once the converted photons leave the resonant layer, the refraction index n(r) quickly increases to values close to 1 as the number density of electrons n e (and thus the plasma frequency, according to (5)) abruptly decreases with the distance. For example, as we can see from Fig. 1a, n e drops by about two orders of magnitude from the height 10 4 km to the height 10 6 km. Therefore, the converted photons will quickly propagate nearly radially θ π/2 after traveling a distance based on the refraction law (22).
Absorption effect-In addition to the refraction effect discussed above, the converted photons could be absorbed by interacting with the plasma along the propagation. The absorption effect is mainly due to the inverse bremsstrahlung process, the rate of which is estimated as [67,72] Γ invB 8πα 3 n e n I 3m 2 e ω 3 2πm e T ln 2T 2 n I is the number density of ions in the plasma. This equation is derived under the condition ω T m e , which is satisfied in both solar corona and solar wind plasma, as can be seen in Fig. 1. Taking into consideration the long-range Coulomb force cut off by the plasma screening effect (Debye screening length λ d ∼ T/(m e ω 2 p )), we have the ln factor in (24). At first glance, it seems that the inverse bremsstrahlung rate ∼ ω −3 will diverge for low-energy photons, but the numerator factor n e n I n 2 e ∝ ω 4 p ∝ ω 4 at the resonant layer actually makes Γ invB decease quickly as ω → 0.
Scattering effect-The converted photons also interact with plasma via Compton scattering, the rate of which is Such interaction can shift the photon energy and propagation direction. We will elaborate more on this effect in the following. To estimate how the converted photons are affected by the effects of absorption and scattering, one can define an attenuation factor by simply adding the two effects together, Γ att = Γ invB + Γ Com [67]. Then, the survival probability of photons after conversion can be estimated as [67] It denotes the probability that a converted photon is neither absorbed nor scattered during propagation. Compared with the inverse bremsstrahlung process (24), the Compton scattering plays a minor role. This can be seen by taking the ratio of (25) to (24) which gives ∼ (m e /T) −3/2 1.
In addition, the propagating photons can scatter with random electron density fluctuations (see e.g., Ref. [73]), which is labeled as 'irregular refractions' in Ref. [74]. Ref. [73] treated the density fluctuations effectively static, and therefore the scatterings do not change the photon energy but only the photon propagating direction. Therefore, (26) represents a simplified and conservative scenario. More realistically, the converted photons will be scattered into different directions and get an angular distribution. As pointed out in Ref. [73], the directions of solar radio emissions will be quickly randomized by scattering effect (mainly due to the electron density fluctuations) near the resonant region where the plasma frequency is equal to the emission frequency, while the large-scale refraction generally tends to shift the propagation path towards the radial direction. Solving exactly the evolution dominated by repeated scatterings and inverse bremsstrahlung processes is difficult, and this problem is usually transformed into the framework of Fokker-Planck equation (see e.g., Refs. [73,[75][76][77]). Essentially, the Fokker-Planck equation can describe the evolution of the phase-space distribution N(r, k, t) of photons, where the refraction effect represented by the gradually-varying n e (r), the absorption effect represented by a collision absorption coefficient, and the scattering effect represented by a diffusion tensor, are included [73,75]. Based on the Fokker-Planck equation, Ref. [73] has developed a numerical Monte Carlo ray-tracing method to study the 3-dimension stochastic propagation process of radio emissions where the anisotropic n e fluctuation has been considered.
Applying the ray-tracing numerical code [73], Ref. [78] has studied the propagation of converted photons in the plasma. The survival probability P sur defined in the numerical study is only related to absorption regardless of whether scattered or not, which is different from that in (26). The scattering effect against the refraction effect is represented by the angular distribution function g(θ) where θ is the angle between the final propagation direction and the radial direction on the last scattering surface (∼ 6R , beyond which the interactions of photons with the plasma are negligible). The numerical results show that P sur 0.5 − 0.07 and g(θ) distributes mostly within the range θ ∼ 0.15 rad, in the radio frequency range 30 − 80 MHz, and both the two effects are stronger for higher frequencies.

Detection
For dark matter particles with typical initial velocity v 0 220 km/s, the solar gravitational effect cannot be neglected when they travel close to the Sun since GM /(R v 2 0 ) 3.45 is an O(1) number. The trajectory of a dark matter particle from far away is hyperbolic in the solar gravitational potential. In the polar coordinates (r, θ) centered at the Sun, the dark matter trajectory can be expressed as where b is the impact factor and e is the eccentricity of trajectory. The tangential velocity v θ and radial velocity v r of dark matter are respectively The dark matter velocity v 0 follows the Maxwellian distribution in the Galactic frame [79][80][81]: and we have ∞ 0 dv 0 f G (v 0 ) = 1. v p is the most probable speed, which is taken as the speed of Local Standard of Rest, i.e., the circular speed about the Galactic center at the solar position, v p v 220 km/s. The distribution (29) should be truncated at v esp which is the Galactic escape speed at the solar position, v esp 544 km/s. However, this modification is exponentially tiny as ∞ v esp dv 0 f G (v 0 ) ≈ 0.66%, which can be safely ignored. Next, we make a Galilean boost to get the velocity distribution in the rest frame of the Sun (see e.g., Ref. [80]): Dark photons or axions in the dark matter wind can reach the resonant-conversion layer r c if their impact parameters b are within Then, the total power of the converted photons per unit of solid angle flowing out of the sphere r c is [67] dP 0 dΩ The prefactor 2 accounts for the fact that the dark matter particles entering and leaving the resonant sphere can both resonantly convert into photons with equal probability. The converted photons propagating inward will be totally reflected as the plasma frequency is larger in deeper region. Then, with the Maxwellian velocity distribution (30) taken into consideration, the emitted power (32) is corrected as [78] dP The converted photons will generate a spike-like feature superposed on the normal solar radio spectrum. Such signals can be detected by the radio telescopes such as SKA and LOFAR when they reach the Earth. For Earth-based radio telescopes, the spectral flux density (in unit of W/m 2 /Hz) received is thus [67,78] The factor β represents the suppression on the final detected signal due to the scattering effect. It depends on the angular distribution g(θ) that requires a detailed ray-tracing numerical simulation [78]. If there is no scattering effect, then β = 1. B is the frequency bandwidth over which the signal distributes. We average the signal over the bandwidth and get the signal per unit of frequency as shown in (34). B should be chosen as the larger one between the telescope resolution B res and the natural bandwidth B DM . Due to the velocity dispersion of dark matter, the converted photons have a natural bandwidth of frequency [67,82]: On the other hand, we need to know the telescope's ability to detect weak signals. The minimal detectable spectral flux density of a telescope is [67,83] S min = SEFD η s n pos Bt obs . (36) η s is the system efficiency of a telescope, which is 0.9 and 1 for SKA and LOFAR respectively [83,84]; n pol = 2 is the number of polarizations of photons; t obs is the observation time. SEFD is the system equivalent (spectral) flux density: T sys is the system temperature of the radio telescope antenna; T nos is the noise temperature induced on the antenna when pointing towards the Sun, which can be estimated as the blackbody temperature of the quiet Sun [67]; A eff is the effective area of the radio telescope array such as SKA and LOFAR. Then, taking the signal (34) equal the telescope detection ability (36), one can find the upper limit of the coupling g aγγ for the axion case or for the dark photon case that can be detected by the telescope. Eq. (38) reflects how large the parameter space of dark photon/axion dark matter can be explored via observing solar radio emissions.
The effective area A eff for SKA1 array is about 10 4 − 10 5 m 2 , while it is about 10 3 m 2 for LOFAR. The antenna system temperatures T sys are also given in Ref. [67] based on the data of Ref. [47,83], ranging from ∼ 10 to 10 4 K for the five frequency bands. Lastly, the telescope resolution B res for SKA1-Low, SKA1-Mid and LOFAR are respectively 1, 3.5 and 197 kHz. All of them are larger than the natural bandwidth of dark matter (35) for the radio frequency bands we are interested in. Therefore, one should choose B = B res when computing the received signal (34).
With the above properties of the radio telescopes known, the sensitivity reaches for the dark photon case based on (38) is plotted in Fig. 2. The observation time is assumed to be 1 hour and 100 hours respectively. As shown in the plot, the sensitivity reaches of the kinetic mixing parameter based on SKA and LOFAR radio telescopes could be several orders of magnitude better than the existing constraints from the cosmic microwave background (CMB) [8,86] and the cavity experiment WISPDMX [87].
The sensitivity reaches obtained in Fig. 2 can be translated to the axion case via the relation (21), and the result is plotted in Fig. 3. As shown in Ref. [88], the magnetic field strength is about 1 − 4 Gauss in the solar corona at 1.05 − 1.35R from the Sun center. To be conservative, we take the transverse magnetic field |B T | = 1 Gauss in (21). As we can see from Fig. 3, our projected sensitivities are better than the Light-Shining-through-a-Wall experiments (LSW, including CROWS, ALPS, and OSQAR [89][90][91]). Also, they are better than the constraint from the CAST experiment in part of the SKA1 frequency range. However, in most of the relevant frequency range, they cannot exceed the constraints from multiple The projected sensitivity reaches of the parameter space of the dark photon dark matter based on the radio telescopes SKA phase 1 and LOFAR, with the observation time assumed to be 1 hour and 100 hours respectively. The plot also shows the existing constraints from CMB [8,86], the cavity experiment WISPDMX [87], and the multiple haloscope-type experiments [8,[100][101][102][103][104]. Figure adapted from Ref. [67].
astrophysical searches in various environments [57,[92][93][94][95][96][97][98] (which is shown as a combined constraint in Fig. 3). The performance of sensitivity reaches in the axion case is not as good as that in the dark photon case. This is because the solar magnetic field is relatively low which suppresses the axion-photon conversion rate as we can see from (21). Figs. 2 and 3 are based on the simplified scenario of the propagation process of the converted photons that took P sur as (26) and assumed β = 1 in (34). Instead, Ref. [78] has numerically simulated the radio photon propagation based on the Monte Carlo ray-tracing method. Also, the Maxwellian velocity distribution has been considered in computing the emitted power (33), although this effect does not make much difference. Furthermore, Ref. [78] has analyzed real data of LOFAR (three samples of 17 min observations at different dates). Finally, Ref. [78] gets a constraint up to ∼ 10 −13 which is about 1 to 2 orders of magnitude better than the CMB constraint [8,86] in the frequency range 30 − 80 MHz. This constraint is weaker than the projections in Fig. 2, mainly because the angular distribution in (34) of converted photons based on the numerical simulation has been considered. A similar constraint on the axion-photon coupling strength g aγγ ∼ 10 −9 GeV −1 has been obtained in the same frequency range [78], which is about 1 to 2 orders of magnitude better than the existing constraints from the LSW experiments [89][90][91], but does not exceed the constraints from the CAST experiment [99] and multiple astrophysical bounds [57,[92][93][94][95][96][97][98].
In addition to the terrestrial facilities such as SKA and LOFAR discussed above, space radio probes are also powerful tools for searching for ultralight dark matter signals from resonant conversion. The STEREO observatory [49,127] consists of two nearly identical satellites, observing the Sun in the 1 AU orbit. It carries a low-frequency radio receiver operating in the frequency range 2.6 kHz − 150 kHz and a higher-frequency receiver operating in the frequency range 125 kHz − 16 MHz. Furthermore, the Parker Solar Probe [50,69] provides new opportunities for observing the conversion signal with the perihelion reaching ∼ 10R from  (21) with |B T | taken as 1 Gauss [88]. Together we show the existing constraints, which can be classified into four categories: the astrophysical searches for axions in various environments, including white dwarfs, neutron stars and pulsars, quasars and blazars, the radio galaxy NGC 1275, and Globular Clusters, etc. [57,[92][93][94][95][96][97][98]; various haloscope experiments aiming to detect axions in the Galactic dark matter halo [100,102,; the helioscope, CAST, aiming to detect axions emitted by the Sun; and Light-Shining-through-a-Wall (LSW) experiments including CROWS, ALPS, and OSQAR [89][90][91]. The existing constraints are also nicely summarized in Ref. [126]. the Sun center. In fact, the Parker satellite is the first spacecraft that can fly into the solar corona. The best advantage for us is that Parker satellite can measure in situ the resonant-conversion events in the solar-wind plasma as shown in Fig. 1b, thanks to its motion between the Sun and 1 AU. The Parker satellite carries a low-frequency radio receiver operating in the frequency range 10 kHz − 1.7 MHz and a higher-frequency radio receiver operating in the frequency range 1.3 MHz − 19.2 MHz. We expect that the Parker satellite has a strong performance in searching for resonant-conversion signals in these frequency ranges. The analysis of space probes such as Paker Solar Probe and STEREO to search for ultralight dark matter is currently in progress by the authors of this Review.

Summary
Solar plasma is an ideal place for detecting ultralight axions or dark photons. When axion dark matter or dark photon dark matter particles cross the solar corona or solar wind, they can resonantly convert into standard model photons when their mass equals the plasma frequency. The plasma frequency of solar corona and solar wind spans from 1 GHz to 0.1 MHz for the distance extending up to 1 AU. Therefore, searching for the radio signals converted from ultralight axions or dark photons in the plasma provides a complementary method to other haloscope studies, to detect these dark matter particles in the mass range ∼ 10 −6 − 10 −10 eV and set stringent constraints on the parameter space of these dark matter models. Based on this mechanism, radio telescopes such as LOFAR and SKA perform well in searching for ultralight dark matter, especially for the dark photon case. Future projected SKA and upgraded FAST are expected to obtain solar radio data with much higher precision. In addition to the terrestrial telescopes, we can use space solar probes such as STEREO and Parker to search for the ultralight dark matter. We conclude that searching for the resonant-conversion radio signals from solar plasma is an efficient method, which is competitive and complementary compared to other direct detection experiments [8,87,[89][90][91][99][100][101][102][103][104] aiming for the ultralight dark matter.