New axion bounds derived from the 100-parsec Gaia DR3 white dwarf luminosity function

Using the high-quality 100-parsec Gaia DR3 white dwarf sample and advanced population synthesis models, this study derives a stringent upper limit on the axion-electron coupling ($g_{ae} < 1.68 \times 10^{-13}$ at 95% C.L.), effectively ruling out the previously suggested range of couplings that implied axion-induced cooling.

The Gist

Imagine the universe as a giant, bustling city. In this city, stars are like buildings. Most buildings (stars) are under construction or fully occupied, but White Dwarfs are the old, retired buildings that have stopped producing new energy. They are simply cooling down, like a cup of coffee left on a table, slowly losing heat until they become cold, dark embers.

Scientists have been watching these "cooling coffee cups" for decades to see if they cool down at the expected speed. If they cool down faster than physics predicts, it means something invisible is stealing their heat.

This paper is a detective story about finding that invisible heat thief: a hypothetical particle called the Axion.

The Mystery: The "Ghost" Particle

The Axion is a ghost-like particle that physicists have been hunting for. It's a leading candidate for "Dark Matter," the invisible stuff that holds galaxies together. If axions exist, they might be born inside the hot cores of these white dwarf stars and fly right out, carrying energy away with them. This would make the stars cool down faster than they should.

The Old Detective Work vs. The New High-Tech Lens

In the past, astronomers tried to count these cooling stars using old telescopes and methods that were a bit like trying to count people in a foggy city by only looking at the brightest streetlights. They had to guess how many dimmer stars were hidden in the fog. This led to some confusion: some studies suggested axions were stealing heat (making stars cool faster), while others said they weren't.

The New Game Changer: Gaia
Enter the Gaia mission, a space telescope that acts like a super-precise 3D map of our cosmic neighborhood. The authors of this paper used a specific "100-parsec sample" from Gaia. Think of this as a perfectly curated list of roughly 7,000 white dwarfs within 100 parsecs of Earth. Because they are so close, we know exactly how far away they are and how bright they truly are. There is no fog, no guessing.

The Experiment: Simulating the City

The researchers didn't just look at the data; they built a virtual city inside a computer.

1. The Simulation: They used a "Monte Carlo" code (a fancy way of saying they ran thousands of random simulations) to create a population of fake white dwarfs.
2. The Variables: They programmed these fake stars to cool down in two ways:
   • Standard Mode: Cooling only by normal physics (like a normal cup of coffee).
   • Axion Mode: Cooling by normal physics plus axions stealing heat.
3. The Comparison: They compared their virtual city's "cooling curve" against the real data from the Gaia telescope.

The Big Discovery: The Heat Thief is Gone

Here is the twist: The real stars are cooling exactly as fast as they should without any help from axions.

When the researchers tried to add axions to their simulation, the virtual stars cooled down too fast. They didn't match the real data. In fact, the data was so precise that it ruled out the "mildly improved" fits that previous studies claimed to find.

The Analogy: Imagine you are timing how long it takes a runner to finish a race.

• Old studies said, "The runner is finishing slightly faster than expected! Maybe there's a tailwind (axions) helping them?"
• This new study says, "We have a perfect stopwatch and a clear track. The runner is finishing exactly on time. There is no tailwind. If there were a tailwind, the runner would be way ahead of schedule, but they aren't."

The Verdict: A Tighter Net

Because the data is so clean, the authors could draw a much tighter line around the possible size of the axion.

• They found that if axions exist, they must be extremely weak at interacting with electrons (specifically, the coupling constant $g_{ae}$ must be less than $1.68 \times 10^{-13}$).
• This is one of the strongest limits we have ever set. It effectively says, "If axions are stealing heat from these stars, they are doing a terrible job."

Why This Matters

This paper is a victory for precision. By using the high-quality "100-parsec" sample from Gaia and running sophisticated computer simulations that account for every possible error, the authors cleared up a decade of confusion.

They proved that the "ghost" axion isn't hiding in the cooling of our nearest white dwarfs. While this doesn't mean axions don't exist at all, it tells us exactly where not to look and forces physicists to refine their theories. It's like telling a detective, "The thief isn't in the kitchen; go check the basement."

In short: We looked at the universe's oldest cooling stars with our sharpest eyes, and they told us, "We are cooling just fine on our own. No invisible heat thieves here."

Technical Summary

Here is a detailed technical summary of the paper "New axion bounds derived from the 100-parsec Gaia DR3 white dwarf luminosity function."

1. Problem Statement

The axion is a hypothetical particle proposed to solve the strong-CP problem in QCD and appears in various extensions of the Standard Model (e.g., string theory). If axions couple to electrons, they can be produced efficiently in the hot, degenerate cores of white dwarfs (WDs) via electron-ion Bremsstrahlung. This process acts as an additional cooling channel, accelerating the cooling rate of WDs.

The White Dwarf Luminosity Function (WDLF)—the distribution of WDs as a function of their luminosity—is a sensitive probe of this cooling rate. Previous studies using magnitude-limited samples (e.g., from SDSS and SuperCOSMOS) and semi-analytical modeling suggested that axion-electron couplings ($g_{ae}$) in the range $0.7 \times 10^{-13} < g_{ae} < 2.1 \times 10^{-13}$provided a slightly better fit to the observed WDLF than the standard model ($g_{ae}=0$). However, these earlier constraints suffered from:

• Observational Biases: Magnitude-limited samples require complex corrections (e.g., $1/V_{max}$) that introduce uncertainties.
• Modeling Simplifications: Previous analyses often relied on semi-analytical approaches that did not fully propagate observational errors or account for the specific selection functions of modern surveys.
• Star Formation Rate (SFR) Uncertainty: The shape of the WDLF is heavily dependent on the history of star formation in the Galactic disk, which is poorly constrained.

This work aims to derive robust, updated constraints on $g_{ae}$ using the high-quality, volume-limited Gaia DR3 100-parsec white dwarf sample and a rigorous Monte Carlo population synthesis framework.

2. Methodology

A. Stellar Models and Physics

• Code: The authors used LPCODE, a state-of-the-art stellar evolution code.
• Physics Inputs: The models include detailed non-gray model atmospheres, conductive and radiative opacities (OPAL), neutrino cooling, and time-dependent element diffusion (H, He, C, N, O, Ne).
• Axion Emission: The dominant cooling mechanism considered is axion Bremsstrahlung in a highly degenerate, strongly correlated plasma ($\Gamma \gg 1$). The emission rate ($\epsilon_{BD}$) was calculated using the latest theoretical prescriptions from [32], which improve upon earlier recipes by accounting for ion correlations.
• Model Grid: The authors generated a grid of 360 white dwarf evolutionary tracks. These cover masses from $0.48$to$1.05 M_\odot$and various metallicities ($Z = 0.0001$to$0.02$). Crucially, they computed tracks for axion-electron coupling constants ($g_{ae}$) ranging from$0$to$8.4 \times 10^{-13}$.

B. Population Synthesis

• Monte Carlo Simulation: A custom code generated synthetic populations of the Galactic thin disk.
  • Initial Mass Function (IMF): Salpeter distribution ($N(m) \propto m^{-2.35}$).
  • Star Formation Rate (SFR): A constant SFR over a 10.6 Gyr history was adopted. The authors demonstrated that for the brightest WDs ($M_G \lesssim 11.5$), the WDLF shape is largely independent of the SFR, making this region robust against SFR uncertainties.
  • IFMR: An Initial-to-Final Mass Relation from [70] was used to map progenitor masses to WD masses.
• Observational Forward-Modeling: To compare theory with data, the synthetic populations were processed through a "forward-modeling" pipeline:
  • Conversion of physical parameters to Gaia passband magnitudes ($G, G_{BP}, G_{RP}$).
  • Incorporation of realistic Gaia DR3 astrometric and photometric errors and selection functions.
  • Application of the same classification criteria (DA vs. non-DA) used in the observational sample.

C. Data Selection

• Sample: The 100-parsec volume-limited sample from Gaia DR3 [28].
• Constraints: The analysis focused on DA white dwarfs (hydrogen atmosphere) with masses between 0.5 and 0.7 $M_\odot$ and absolute magnitudes $M_G < 11.5$.
  • This mass range excludes He-core WDs and massive ONe-core WDs (minimizing binary evolution contamination).
  • The magnitude cut ($M_G < 11.5$) ensures the sample is dominated by recent star formation, rendering the WDLF shape insensitive to the SFR history.

D. Statistical Analysis

• Chi-Squared Test: The authors compared the observed WDLF counts with the synthetic models using Pearson's $\chi^2$ statistic.
• Binning: Multiple binning schemes (uniform in magnitude and uniform in star count) were tested to ensure robustness.
• Confidence Levels: Constraints were derived at 68%, 95%, and 99.7% confidence levels (C.L.) based on the resulting $p$-values.

3. Key Contributions

1. First Application of Gaia DR3 100-pc Sample to Axion Bounds: This is the first study to utilize the high-purity, volume-limited Gaia DR3 sample for axion constraints, eliminating the biases inherent in previous magnitude-limited catalogs.
2. Rigorous Forward-Modeling: Unlike previous semi-analytical studies, this work uses a full Monte Carlo population synthesis code that explicitly folds in observational errors and selection effects, providing a more faithful comparison between theory and observation.
3. Updated Bremsstrahlung Rates: The implementation of the latest axion emission formulas (accounting for strong ion correlations) ensures the theoretical cooling rates are physically accurate.
4. Decoupling SFR Uncertainty: By isolating the bright end of the WDLF ($M_G < 11.5$), the authors demonstrated that the constraints on axion cooling are independent of the uncertain Galactic star formation history, a major source of error in previous works.

4. Results

• Disfavored Axion Coupling: The analysis strongly disfavors the "improved fit" suggested by earlier studies. The data is best described by the standard model ($g_{ae} = 0$).
• New Upper Limits: The study establishes the following upper limits on the axion-electron coupling constant ($g_{ae}$):
  • 68% C.L.: $g_{ae} < 1.12 \times 10^{-13}$
  • 95% C.L.: $g_{ae} < 1.68 \times 10^{-13}$
  • 99.7% C.L.: $g_{ae} < 2.52 \times 10^{-13}$
• Mass Constraints: Translating these limits to the QCD axion mass ($m_a$) using the DFSZ model relation ($g_{ae} \propto m_a \cos^2 \beta$), the constraints correspond to:
  • $m_a \cos^2 \beta < 4$ meV (68% C.L.)
  • $m_a \cos^2 \beta < 6$ meV (95% C.L.)
  • $m_a \cos^2 \beta < 9$ meV (99.7% C.L.)
• Fit Degradation: The quality of the fit deteriorates abruptly for $g_{ae} > 0.7 \times 10^{-13}$, indicating that even modest axion couplings are inconsistent with the observed bright WD population.

5. Significance

• Tighter Constraints: The derived limit of $g_{ae} < 1.68 \times 10^{-13}$ (95% C.L.) is among the strongest available astrophysical bounds, significantly improving upon previous constraints derived from older, less precise datasets.
• Resolution of Discrepancy: The paper resolves the tension with earlier works that found a preference for non-zero axion couplings. The authors attribute the previous "improvement" to simplifying assumptions in modeling and the lower quality of magnitude-limited samples, which were susceptible to SFR uncertainties and selection biases.
• Methodological Benchmark: The study sets a new standard for using stellar populations to constrain particle physics, demonstrating that volume-limited samples combined with forward-modeling Monte Carlo simulations are essential for robust results.
• Complementarity: While the globular cluster 47 Tuc provides even tighter constraints (due to a well-known SFR and distance), the Galactic disk analysis offers a complementary probe that is less susceptible to cluster-specific systematics (e.g., crowding, mass segregation).

In conclusion, this work provides a robust, statistically rigorous exclusion of significant axion-induced cooling in white dwarfs, pushing the upper bound on the axion-electron coupling to $1.68 \times 10^{-13}$ at 95% confidence.