🔭 astrophysics

Theoretical study of the synthesis of hydrogen, helium and lithium in the outermost region of the accretion disks of the Vela X-1 binary system

This theoretical study utilizes nuclear lattice simulations to analyze the transmutation of free neutrons into hydrogen, helium, and lithium isotopes within the outermost regions of the Vela X-1 accretion disk across a temperature range of 10⁶ to 10⁷ K, estimating neutron lifetimes and resulting mass fraction abundances.

Original authors: Loidel Puentes-Milián

Published 2026-06-30✓ Author reviewed ⓘ

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Loidel Puentes-Milián

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a cosmic kitchen where a neutron star and a giant blue star are locked in a dance. Between them swirls a massive, spinning disk of gas and dust called an accretion disk. This is the setting for the study of Vela X-1, a famous binary system in our universe.

While we often think of stars as the only places where new elements are cooked up (like a giant cosmic oven), this paper asks a different question: What happens in the cooler, outer edges of this swirling disk?

Here is a simple breakdown of what the author, Loidel Puentes-Milián, discovered by running thousands of computer simulations.

1. The Cosmic "Lego" Set

Think of the disk as a box of loose Lego bricks. The main bricks are protons (hydrogen) and helium. But there are also a few loose, unstable pieces called free neutrons.

In this "kitchen," the temperature isn't uniform. The center is scorching hot, but the outer edges (where this study focused) are cooler, ranging from 1 million to 10 million degrees. Even though that sounds hot to us, it's actually the "chill zone" of the disk.

The author built a digital model—a "nuclear lattice"—to watch how these Lego bricks interact. They wanted to see how the free neutrons would crash into other bricks to build new things like deuterium (heavy hydrogen), tritium (super-heavy hydrogen), and lithium.

2. The "Neutron Lifespan" Mystery

Free neutrons are like unstable jellybeans; they don't last long on their own. They usually decay quickly. The study asked: How long do these neutrons survive in the outer disk before they disappear or get used up?

The Surprise: The author found that the neutrons live for about the same amount of time (roughly 1,100 to 1,200 seconds) regardless of whether the spot is slightly warmer or slightly cooler.

The Analogy: Imagine a crowded dance floor. Usually, if the music gets faster (hotter), people move around more and bump into each other more often. You'd expect the neutrons to get "used up" faster in the hot spots. But here, the presence of lithium acts like a giant magnet. Lithium grabs the neutrons so efficiently that it doesn't matter how fast the music is playing; the neutrons get snatched up at the same steady pace. This makes the "lifespan" of the neutrons very consistent across the whole outer region.

3. The "Hot vs. Cold" Chemistry

The study tracked what happened to the different types of atoms after the neutrons were all used up.

The Heat Sensitive Crew (Deuterium, Tritium, Helium-3): These are like delicate glass figurines. In the hotter parts of the disk, they get smashed to pieces easily. As the temperature goes up, their numbers drop significantly.
The Tough Cookies (Lithium): Lithium is like a rock. It doesn't care much about the heat. The study found that the amount of Lithium-6 and Lithium-7 remained almost the same whether the disk was cool or hot. Because Lithium-6 grabs neutrons so fast, it turns into Lithium-7 almost immediately, keeping their ratios steady.

4. The "Electron Balance"

The disk is a plasma, which means it's a soup of charged particles. For the soup to stay stable, the number of positive charges (protons) must equal the number of negative charges (electrons).

The study confirmed that the ratio of protons to electrons stays constant, just as physics demands. However, because protons are much heavier than electrons (like a bowling ball vs. a ping-pong ball), the mass ratio isn't 1-to-1. The simulation showed that for every unit of electron mass, there are about 1,000 units of proton mass, which makes perfect sense.

5. The Final Recipe

After running 5,000 different scenarios for each temperature, the author concluded what the outer edge of the Vela X-1 disk is actually made of once the chaos settles:

The Main Ingredients: Mostly Protons (Hydrogen) and Helium-4.
The Traces: Tiny amounts of Lithium.
The Missing Pieces: Very little Deuterium, Tritium, or Helium-3, especially in the warmer spots, because the heat destroys them.

Summary

In simple terms, this paper is a theoretical recipe book for the outer edge of a cosmic disk. It tells us that while the heat of the disk can destroy fragile atoms like tritium and deuterium, the presence of lithium acts as a stabilizer for neutrons. The result is a region dominated by simple hydrogen and helium, with a steady, unchanging amount of lithium, regardless of the temperature fluctuations. It's a snapshot of how the universe builds (and breaks) the lightest elements in the wild environment around a neutron star.

Technical Summary: Theoretical Study of Hydrogen, Helium, and Lithium Synthesis in the Outer Accretion Disks of Vela X-1

Problem Statement
While nucleosynthesis is traditionally associated with stellar interiors and explosive events, recent theoretical frameworks suggest that accretion regions around compact objects may harbor thermodynamic and irradiation conditions sufficient to activate nuclear reaction channels distinct from those in main-sequence stars. Specifically, the outer regions of accretion disks (temperatures $10^6$ – $10^7$ K) may act as laboratories for light particle synthesis and destruction. However, the specific impact of these environments on light isotopes (hydrogen, helium, and lithium) and the dynamics of charge carriers (free electrons) remain relatively unexplored. This study addresses the gap in understanding the transmutation of free neutrons into isotopes of hydrogen, helium, and lithium within the outermost region of the Vela X-1 binary system's accretion disk.

Methodology
The research employs a theoretical nuclear lattice model to simulate the evolution of nuclear abundances under non-equilibrium conditions.

System Parameters: The simulations focus on the outer region of the Vela X-1 accretion disk, utilizing a density of $\rho = 10^{-15}$ g/cm $^3$ and a temperature range of $T = 10^6$ to $10^7$ K.
Nuclear Network: A network of 27 nuclear reactions was constructed, encompassing neutron and proton capture, electron capture, and beta decays involving neutrons ( $n$ ), protons ( $p$ ), deuterium ( $d$ ), tritium ( $t$ ), helium-3 ( $^3$ He), helium-4 ( $^4$ He), lithium-6 ( $^6$ Li), lithium-7 ( $^7$ Li), and electrons ( $e^-$ ). Endothermic reactions with MeV thresholds were excluded as energetically unfeasible for the target region.
Initial Conditions: Initial mass fractions were derived from the composition of the companion star (HD 77581), dominated by protons ( $X_p \approx 0.7$ ) and helium-4 ( $X_{^4He} \approx 0.28$ ), with trace amounts of lithium and negligible tritium. To account for variability, 5,000 sets of initial conditions were randomly generated for each of the 10 temperature steps, normalized to satisfy mass conservation ( $\sum X_i = 1$ ).
Numerical Integration: The temporal evolution of molar abundances was solved using a system of 9 coupled differential equations. Integration was performed using a 4th-order Adams-Moulton predictor-corrector algorithm, initialized by a 4th-order Runge-Kutta method. Simulations ran until free neutrons were depleted ( $X_n < 10^{-8}$ ) or a maximum time of $1.5 \times 10^4$ s was reached.

Key Contributions and Results
The study provides quantitative estimates of free neutron lifetimes and the final mass fraction abundances of light isotopes across the specified temperature range.

Free Neutron Lifetimes:
- The average lifetime of free neutrons was estimated using exponential decay fits (single and double exponential).
- Finding: The lifetimes exhibit quasi-homogeneity across the temperature range, showing almost complete independence from temperature.
- Mechanism: This stability is attributed to the presence of lithium-6. The reaction rates for neutron capture by lithium-6 (reactions 24 and 25) are temperature-independent and significantly higher than other neutron-consuming reactions, thereby dominating the depletion timescale.
Relative Isotopic Distributions:
- Hydrogen Isotopes: The ratio of protons to deuterium ( $X_p/X_d$ ) increases with temperature, indicating that deuterium abundance decreases as temperature rises. The ratio of protons to tritium ( $X_p/X_t$ ) remains extremely high ( $10^{10}$ – $10^{11}$ ), confirming the scarcity of tritium production in the outer disk compared to inner regions.
- Helium Isotopes: The ratio of helium-4 to helium-3 ( $X_{^4He}/X_{^3He}$ ) increases with temperature, suggesting that higher temperatures favor the destruction of helium-3.
- Lithium Isotopes: The ratio of lithium-7 to lithium-6 ( $X_{^7Li}/X_6Li$ ) remains approximately unity (1:1) across most configurations and temperatures.
- Electron-Proton Ratio: The mass fraction ratio $X_p/X_{e^-}$ remains quasi-independent of temperature, clustering around $10^3$ . This reflects the mass difference between protons and electrons while maintaining plasma charge neutrality.
Absolute Abundance Distributions:
- Temperature Dependence: Deuterium, tritium, and helium-3 abundances decrease as temperature increases. Tritium shows a temperature-independent minimum abundance floor.
- Lithium Homogeneity: The mass fraction abundances of lithium-6 and lithium-7 remain nearly constant across all temperatures. This is a direct consequence of the temperature-independent reaction rates governing lithium-6 consumption and production.
- Stability of Major Species: Protons and helium-4 show no perceptible changes in abundance during the neutron depletion phase, confirming their dominance and stability in the outer disk environment.

Significance and Claims
The paper claims to extend previous theoretical frameworks by explicitly incorporating lithium synthesis and the dynamic role of free electrons in accretion disk nucleosynthesis. The primary significance of the work lies in establishing a quantitative link between nuclear microphysics and electronic interactions in non-equilibrium accretion environments.

The authors conclude that the outermost layer of the Vela X-1 accretion disk is chemically stable regarding major species (protons and helium-4) but acts as a selective filter for light isotopes:

High temperatures efficiently destroy deuterium, tritium, and helium-3.
Lithium isotopes are largely unaffected by the thermal conditions of the disk plasma due to specific reaction kinetics.
Free neutrons persist for a consistent duration regardless of the specific zone in the outer disk, driven by lithium interactions.

Ultimately, the study posits that the outer accretion disk of Vela X-1 is composed principally of protons and helium-4, with traces of lithium, and serves as a natural laboratory where specific nuclear reaction channels (particularly those involving lithium) dictate the chemical evolution of light elements under extreme, yet non-explosive, conditions.