X-ray Searches for Axions from Super Star Clusters Christopher Dessert,1, 2, 3 Joshua W. Foster,1, 2, 3 and Benjamin R. Safdi1, 2, 3 1Leinweber Center for Theoretical Physics, Department of Physics, University of Michigan, Ann Arbor, MI 48109 U.S.A. 2Berkeley Center for Theoretical Physics, University of California, Berkeley, CA 94720, U.S.A. 3Theoretical Physics Group, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, U.S.A. (Dated: March 3, 2021) Axions may be produced in abundance inside stellar cores and then convert into observable X-rays in the Galactic magnetic fields. We focus on the Quintuplet and Westerlund 1 super star clusters, which host large numbers of hot, young stars including Wolf-Rayet stars; these stars produce axions efficiently through the axion-photon coupling. We use Galactic magnetic field models to calculate the expected X-ray flux locally from axions emitted from these clusters. We then combine the axion model predictions with archival Nuclear Spectroscopic Telescope Array (NuSTAR) data from 10 - 80 keV to search for evidence of axions. We find no significant evidence for axions and constrain the −12 −1 −11 axion-photon coupling gaγγ . 3:6 × 10 GeV for masses ma . 5 × 10 eV at 95% confidence. Ultralight axion-like particles that couple weakly to ordinary matter are natural extensions to the Standard 30 Model. For example, string compactifications often pre- 2 dict large numbers of such pseudo-scalar particles that in- 20 teract with the Standard Model predominantly through dimension-five operators [1,2]. If an axion couples to 1 10 quantum chromodynamics (QCD) then it may also solve the strong CP problem [3{6]; in this work we refer to [arcmin] both the QCD axion and axion-like particles as axions. 0 0 0 Axions may interact electromagnetically through the operator = g aF F~µν =4, where a is the axion 10 aγγ µν Dec - Dec − L − 1 field, F is the electromagnetic field-strength tensor, with − ~ 20 10 - 80 keV cts (back. sub.) F its Hodge dual, and gaγγ is the dimensionful coupling − constant of axions to photons. This operator allows both 2 the production of axions in stellar plasmas through the − 30 − Primakoff Process [7,8] and the conversion of axions 2 1 0 1 2 to photons in the presence of static external magnetic − − RA- RA0 [arcmin] fields. Strong constraints on gaγγ for low-mass axions come from the CERN Axion Solar Telescope (CAST) ex- Figure 1. The stacked and pixelated background-subtracted periment [9], which searches for axions produced in the count data (10 - 80 keV) from the NuSTAR observations of Solar plasma that free stream to Earth and then convert the Quintuplet SSC. The locations of the stars are indicated in to X-rays in the magnetic field of the CAST detector. black, while the 90% energy containment region for emission −11 CAST has excluded axion couplings gaγγ & 6:6 10 associated with the SSC is indicated by the black circle, ac- −1 × GeV for axion masses ma . 0:02 eV at 95% con- counting for the NuSTAR point spread function (PSF). RA0 fidence [9]. Primakoff axion production also opens a and DEC0 denote the locations of the cluster center. We find arXiv:2008.03305v2 [hep-ph] 1 Mar 2021 new pathway by which stars may cool, and strong lim- no evidence for axion-induced emission from this SSC, which −11 −1 would follow the spatial counts template illustrated in the its (gaγγ 6:6 10 GeV at 95% confidence for . × inset panel. ma . keV) are derived from observations of the horizon- tal branch (HB) star lifetime, which would be modified in the presence of axion cooling [10]. In this work, we produce some of the strongest con- convert into photons in the Galactic magnetic fields, lead- −9 straints to-date on gaγγ for ma . 10 eV through X-ray ing to signatures observable with space-based X-ray tele- observations with the Nuclear Spectroscopic Telescope scopes such as NuSTAR. We analyze archival NuSTAR Array (NuSTAR) telescope [11] of super star clusters data from the Quintuplet SSC near the Galactic Center (SSCs). The SSCs contain large numbers of hot, young, (GC) along with the nearby Westerlund 1 (Wd1) cluster −12 −1 and massive stars, such as Wolf-Rayet (WR) stars. We and constrain gaγγ . 3:6 10 GeV at 95% con- −11× show that these stars, and as a result the SSCs, are highly fidence for ma . 5 10 eV. In Fig.1 we show the efficient at producing axions with energies 10{100 keV locations of the stars× within the Quintuplet cluster that through the Primakoff process. These axions∼ may then are considered in this work on top of the background- 2 subtracted NuSTAR counts, from 10 - 80 keV, with the lar metallicity is Z 0:02) [31, 32]. Note that higher point-spread function (PSF) of NuSTAR also indicated. metallicities generally≈ lead to the stars entering the WR In the Supplementary Material (SM) we show that ob- classifications sooner, when their cores are cooler. Rota- servations of the Arches SSC yield similar but slightly tion may also cause certain massive stars to be classified weaker limits. as WR stars at younger ages. We model the initial ro- Our work builds upon significant previous efforts to tation distribution as a Gaussian distribution with mean use stars as laboratories to search for axions. Some of µrot and standard deviation σrot for non-negative rota- the strongest constraints on the axion-matter couplings, tion speeds [33, 34]. Refs. [33, 34] found µrot 100 km/s ≈ for example, come from examining how HB, white dwarf and σrot 140 km/s, but to assess systematic uncertain- ≈ (WD), red giant, and neutron star (NS) cooling would ties we vary µrot between 50 and 150 km/s [33]. be affected by an axion [10, 12{20]. When the stars We draw initial stellar velocities from the velocity dis- have large magnetic fields, as is the case for WDs and tribution described above (from 0 to 500 km/s) and ini- NSs, the axions can be converted to X-rays in the stellar tial stellar masses from the Kroupa initial mass func- magnetospheres [21{24]. Intriguingly, in [23, 24] observa- tion [35] (from 15 to 200 M ). We use MESA to evolve tions of the Magnificent Seven nearby isolated NSs found the stars from pre-main-sequence (pre-MS){before core evidence for a hard X-ray excess consistent with the hydrogen ignition{to near-supernova. At each time step expected axion spectrum from nucleon bremsstrahlung. we assign each stellar model a spectroscopic classifica- This work extends these efforts by allowing the axions to tion using the definitions in [36, 37]. We then construct convert to X-rays not just in the stellar magnetic fields an ensemble of models for each spectroscopic classifica- but also in the Galactic magnetic fields [25{27]. tion by joining together the results of the different sim- Axion production in SSCs.|During helium burning, par- ulations that result in the same classification for stellar ticularly massive stars may undergo considerable mass ages within the age range for star formation in the clus- loss, especially through either rotation or binary inter- ter; for Quintuplet, this age range is between 3.0 and action, which can begin to peel away the hydrogen en- 3.6 Myr [38]. Note that each simulation generally pro- velope, revealing the hot layers underneath and revers- vides multiple representative models, taken at different ing the cooling trend. Stars undergoing this process are time steps. In total we compute 105 models per stellar known as WR stars, and these stars are the most im- classification. portant in our analyses. If the star has a small (<40% Quintuplet hosts 71 stars of masses & 50M , with a abundance) remaining hydrogen envelope, it is classified substantial WR cohort [38]. In particular it has 14 WC + as a WNh star; at <5% hydrogen abundance it is clas- WN stars, and we find that these stars dominate the pre- −12 −1 sified as a WN star; otherwise, it is classified as WC or dicted axion flux. For example, at gaγγ = 10 GeV WO, which indicates the presence of >2% carbon, and we compute that the total axion luminosity from the SSC +0:7 35 oxygen, respectively, in the atmosphere. (with Z = 0:035 and µrot = 150 km/s) is 2:1 10 −0:4 × Axions are produced through the photon coupling gaγγ erg/s, with WC + WN stars contributing 70% of that in the high-mass stars in SSCs through the Primakoff flux. Note that the uncertainties arise from∼ performing process γ+(e−;Z) a+(e−;Z). This process converts a multiple (500) draws of the stars from our ensembles of stellar photon to an! axion in the screened electromagnetic representative models. In the 10 - 80 keV energy range field of the nucleons and electrons. The massive stars relevant for NuSTAR the total luminosity is 1:7+0:4 1035 −0:3× are high-temperature and low-density and therefore form erg/s. We take Z = 0:035 and µrot = 150 km/s because nonrelativistic nondegenerate plasmas. The Primakoff these choices lead to the most conservative limits. For emission rate was calculated in [8, 28] as a function of example, taking the metallicity at the lower-end of our temperature, density, and composition, and is described range (Z = 0:018) along with µrot = 100 km/s the pre- in detail in the SM. dicted 10 - 80 keV flux increases by 60%. At fixed ∼ To compute the axion luminosity in a given star, we Z = 0:035 changing µrot from 150 km/s to 100 km/s in- use the stellar evolution code Modules for Experiments creases the total luminosity (over all energies) by 10%, ∼ in Stellar Astrophysics (MESA) [29, 30] to find, at any though the luminosity in the 10 - 80 keV range is virtually particular time in the stellar evolution, radial profiles of unaffected.