Stars have been recognized as optimal laboratories to probe axion properties. In the last decades there have been significant advances in this field due to a better modeling of stellar systems and accurate observational data. In this work we review the current status of constraints on axions from stellar physics. We focus in particular on the Sun, globular cluster stars, white dwarfs and (proto)-neutron stars.

In the presence of a magnetic field, axions can convert into photons and vice versa. The phenomenology of the conversion is captured by a system of two coupled Klein-Gordon equations, which, assuming that the axion is relativistic, is usually recast into a pair of first-order Schrödinger-like equations. In such a limit, focusing on a constant magnetic field and plasma frequency, the equations admit an exact analytic solution. The relativistic limit significantly simplifies the calculations and, therefore, it is widely used in phenomenological applications. In this work, we discuss how to evaluate the axion-photon system evolution without relying on such relativistic approximation. In particular, we give an exact analytical solution, valid for any axion energy, in the case that both the magnetic field and plasma frequency are constant. Moreover, we devise an analytic perturbative expansion that allows for tracking the conversion probability in a slightly inhomogeneous magnetic field or plasma frequency, whose characteristic scale of variation is much larger than the typical axion-photon oscillation length. Finally, we discuss a case of resonant axion-photon conversion giving useful simplified formulae that might be directly applied to dark matter axions converting in neutron star magnetospheres.

We propose a novel class of dark matter (DM) candidates in the form of a heavy composite axionlike particle (ALP) with highly suppressed electromagnetic interactions populating vast yet unexplored domains in the ALP parameter space. This is achieved for the first time in the simplest dark confining gauge theory yielding a new composite glueball ALP (GALP) DM coupling-mass relation found in terms of two distinct fundamental scales-the large dark fermion mass scale and the dynamical scale of dark confinement. The presence of a heavy fermion portal between the visible (photons) and dark (GALPs) sectors ensures a strong radiative suppression of the GALP-photon coupling naturally without any fine-tuning. The observable features of heavy GALP DM in a minimal realization are controlled by only three physical parameters. Our work paves the road for a novel research field exploring the theory and phenomenology of composite ALPs in multimessenger astrophysics and cosmology.

We investigate the potential of IAXO and its intermediate version, BabyIAXO, to detect axions produced in core-collapse supernovae (SNe). Our study demonstrates that these experiments have realistic chances of identifying SN axions, offering crucial insights into both axion physics and SN dynamics. IAXO's sensitivity to SN axions allows for the exploration of regions of the axion parameter space inaccessible through solar observations. In addition, in the event of a nearby SN, d ∼ 𝒪(100) pc, and sufficiently large axion couplings,  g_aγ ≳ 10^-11 GeV^-1, IAXO could have a chance to significantly advance our understanding of axion production in nuclear matter and provide valuable information about the physics of SNe, such as pion abundance, the equation of state, and other nuclear processes occurring in extreme environments.

Core-collapse supernovae (SNe) provide a unique environment to study feebly interacting particles (FIPs) such as axionlike particles (ALPs), sterile neutrinos, and dark photons (DPs). This paper focuses on heavy FIPs produced in SNe, whose decay produces electrons and positrons, generating observable secondary signals during their propagation and annihilation. We focus on the in-flight annihilation of positrons, which emerge as the most significant contribution to the resulting 𝛾-ray spectrum. Using data from COMPTEL and EGRET, we derive the most stringent bounds on the FIP-electron couplings for heavy ALPs, sterile neutrinos, and DPs. These results strengthen existing bounds by one to two orders of magnitude, depending on the FIP model.

We obtain novel constraints on new scalar fields interacting with Standard Model fermions through Lorentz-violating couplings, bridging searches for scalar particles and Lorentz-symmetry tests. These constraints arise from torsion-balance experiments, magnetometer searches, and an excessive energy loss in Red Giant stars. Torsion-balance experiments impose stringent constraints, benefitting from large macroscopic sources including the Sun and Earth. Magnetometer-based searches, which detect pseudo-magnetic fields through spin precession, offer additional limiting power to low-mass scalar fields. Meanwhile, observations of Red Giant stars place strong limits on additional energy loss mechanisms, extending these constraints to higher scalar mass ranges and a wider range of Lorentz-violating couplings. Combining data from laboratory experiments and astrophysical observations, this approach strengthens constraints on Lorentz-violating interactions and paves the way for future investigations into physics beyond the Standard Model.

The γ-ray emission originating from in-flight annihilation (IA) of positrons is a powerful observable for constraining high-energy positron production from exotic sources. By comparing diffuse γ-ray observations of INTEGRAL, COMPTEL, and EGRET to theoretical predictions, we set the most stringent constraints on electrophilic feebly interacting particles, thereby proving IA as a valuable probe of new physics. In particular, we extensively discuss the case of MeV-scale sterile neutrinos, where IA sets the most stringent constraints, excluding |Uμ4|2≳10-13 and |Uτ4|2≳2×10-13 for sterile neutrinos mixed with μ and τ neutrinos, respectively. These constraints improve existing limits by more than an order of magnitude. We briefly discuss the application of these results to a host of exotic positron sources such as dark photons, axionlike particles, primordial black holes, and sub-GeV dark matter.

We propose a new theory of dark matter based on axion physics and cosmological phase transitions. We show that theories in which a gauge coupling increases through a first-order phase transition naturally result in ‘axion relic pockets’: regions of relic false vacua stabilised by the pressure from a kinematically trapped, hot axion gas. Axion relic pockets provide a viable and highly economical theory of dark matter: the macroscopic properties of the pockets depend only on a single parameter (the phase transition temperature). We describe the formation, evolution and present-day properties of axion relic pockets, and outline how their phenomenology is distinct from existing dark matter paradigms. We briefly discuss how laboratory experiments and astronomical observations can be used to test the theory, and identify gamma-ray observations of magnetised, dark-matter-dense environments as particularly promising.

Sterile neutrinos with masses up to (100)  MeV can be copiously produced in a supernova (SN) core through the mixing with active neutrinos. In this regard, the SN 1987A detection of neutrino events has been used to put constraints on active-sterile neutrino mixing, exploiting the well-known SN cooling argument. We refine the calculation of this limit including neutral current interactions with nucleons, which constitute the dominant channel for sterile neutrino production. We also include, for the first time, the charged current interactions between sterile neutrinos and muons, relevant for the production of sterile neutrinos mixed with muon neutrinos in the SN core. Using the recent modified luminosity criterion, we extend the bounds to the case where sterile states are trapped in the stellar core. Additionally, we study the decays of heavy sterile neutrinos, affecting the SN explosion energy and possibly producing a gamma-ray signal. We also illustrate the complementarity of our new bounds with cosmological bounds and laboratory searches.

High-frequency gravitational waves (f≳1 MHz) may provide a unique signature for the existence of exotic physics. The lack of current and future gravitational-wave experiments sensitive at those frequencies leads to the need of employing different indirect techniques. Notably, one of the most promising ones is constituted by graviton-photon conversions in magnetic fields. In this work, we focus on conversions of a gravitational-wave background into photons inside the Milky Way magnetic field, taking into account the state-of-the-art models for both regular and turbulent components. We discuss how graviton-to-photon conversions may lead to imprints in the cosmic photon background spectrum in the range of frequencies f∼109-1026 Hz, where the observed photon flux is widely explained by astrophysics emission models. Hence, the absence of any significant evidence for a diffuse photon flux induced by graviton-photon conversions allows us to set stringent constraints on the gravitational-wave strain hc, strengthening current astrophysical bounds by ∼1-2 orders of magnitude in the whole range of frequencies considered.