Dissociative recombination and ion-pair formation in $\mathrm{HeH^+}$ isotopologues: A time-dependent wave-packet study including rotational coupling

This study employs time-dependent wave-packet propagation to demonstrate that including a large manifold of resonant states and rotational couplings significantly enhances the dissociative recombination and resonant ion-pair formation cross sections for $\mathrm{HeH^+}$ isotopologues, thereby highlighting the critical role of multi-state nonadiabatic effects in accurately modeling electron-molecule collisions in astrophysical plasmas.

The Gist

The Big Picture: A Cosmic First Step

Imagine the early universe as a giant, empty construction site. Before stars and galaxies could form, the very first "building block" had to be created. Scientists believe that block was a molecule made of one helium atom and one hydrogen atom stuck together, called HeH+. It's like the "first brick" of the universe.

However, this first brick is fragile. It's constantly being hit by tiny, fast-moving particles called electrons. When an electron hits the HeH+ molecule, two things can happen:

1. Dissociative Recombination (DR): The electron sticks to the molecule, causing it to instantly shatter into a helium atom and a hydrogen atom.
2. Resonant Ion-Pair Formation (RIP): The electron hits the molecule, causing it to split into two charged pieces: a positive helium ion and a negative hydrogen ion.

This paper is a detailed computer simulation of exactly how these collisions happen.

The New Approach: A Bigger Net and More Spin

Previous scientists tried to simulate these crashes, but they were looking at the problem through a narrow keyhole. They only watched a few specific "paths" the molecule could take and ignored how the molecule spins.

The authors of this paper built a much more sophisticated simulation. Think of it like upgrading from a simple fishing rod to a massive, high-tech net.

• The Bigger Net (More States): Instead of watching just a few paths, they tracked 23 different electronic states (different ways the electrons inside the molecule can arrange themselves). This is like checking 23 different escape routes instead of just one.
• The Spin (Rotational Coupling): They also included how the molecule spins as it flies. Imagine a spinning top; if it spins fast, it might wobble and change direction. The authors realized that this "wobble" (rotational coupling) helps the molecule find new ways to break apart that previous models missed.

What They Found: The Breakup is Faster Than We Thought

When they ran their new, more complex simulation, they found something surprising: The molecule breaks apart much more easily than we previously thought.

• The "Shatter" Rate: The probability of the molecule breaking (the cross-section) is significantly higher in their new model. It's like realizing that a glass vase is actually made of a much more brittle material than we thought; it shatters with a much lighter tap.
• The Spin Matters: They found that the spinning motion of the molecule acts like a bridge, helping the electrons jump between different energy levels and making the breakup more likely.
• The "Heavy" vs. "Light" Effect: They tested different versions of the molecule (using heavier or lighter isotopes, like swapping regular hydrogen for "heavy" hydrogen). They found a clear rule: The lighter the molecule, the faster it breaks.
  • Analogy: Imagine two runners on a track. The lighter runner (lighter isotope) runs so fast that they sprint past the "danger zone" before they can trip. The heavier runner (heavier isotope) moves slower, giving them more time to trip and fall (break apart). Wait, actually, the paper says the opposite for the result: The lighter molecules break apart more often because they move so fast through the critical zone that they successfully escape before the electron can bounce back off. It's a race against time where the faster runner wins the "breakup" more often.

Two Ways to Look at the Same Thing

The authors ran the simulation in two different mathematical "languages" (Adiabatic and Diabatic).

• Adiabatic: This is like watching a movie where the scenery changes smoothly as the characters move.
• Diabatic: This is like watching the same movie but focusing on the characters' internal states changing instantly.
  They found that while both languages tell the same story, they highlight different details. In one language, certain types of spins (called $2\Sigma$) are the main heroes causing the breakup. In the other, different spins ($2\Pi$) play a bigger role at lower speeds.

Why This Matters for the Universe

The paper concludes that because the molecule breaks apart more easily than we thought, it might not survive as long in the early universe as some old models predicted.

• The Cosmic Balance: If HeH+ breaks apart too quickly, there might be less of it floating around in space than we thought.
• The "First Brick" Status: Since HeH+ is considered the first molecule in the universe, knowing exactly how fast it gets destroyed helps astronomers understand the chemistry of the early cosmos, the gas clouds between stars, and the glowing shells around dying stars (planetary nebulae).

Summary

In short, this paper says: "We built a better, more detailed computer model of how the universe's first molecule gets destroyed by electrons. We found that it breaks apart much more easily than we thought, especially when it's spinning and when it's made of lighter ingredients. This means we need to update our maps of the early universe to account for this faster destruction."

Technical Summary

Technical Summary: Dissociative Recombination and Ion-Pair Formation in HeH+ Isotopologues

Problem Statement
Dissociative recombination (DR) and resonant ion-pair (RIP) formation are critical electron-driven processes governing the evolution of molecular plasmas, particularly in primordial and astrophysical environments. The helium hydride ion (HeH+), recognized as the first molecular ion formed in the Universe, serves as a primary benchmark for understanding these mechanisms. While previous theoretical studies have utilized multichannel quantum defect theory (MQDT), R-matrix methods, and time-dependent wave-packet approaches, earlier dynamical models often suffered from limitations: they were restricted to a small subset of resonant states and frequently omitted rotational couplings between different electronic symmetries. Consequently, discrepancies remained between theoretical predictions and experimental storage-ring measurements, which indicated significantly higher recombination rates than simple direct-capture models could explain. This study addresses the need for a comprehensive theoretical investigation that incorporates a large manifold of coupled electronic states and explicit rotational interactions to accurately describe electron–molecule collision processes in HeH+ and its isotopologues.

Methodology
The authors employ a time-dependent wave-packet propagation method to solve the nuclear dynamics of the system. The study is characterized by three main methodological advancements:

1. Extended Electronic Manifold: Nuclear dynamics are treated on a manifold of 23 coupled resonant electronic states, encompassing $^2\Sigma$, $^2\Pi$, and $^2\Delta$ symmetries. This represents a significant expansion over previous wave-packet studies that utilized only a few states.
2. Explicit Rotational Coupling: The Hamiltonian explicitly includes rotational coupling terms (L-uncoupling) between states of different symmetries ($\Sigma$-$\Pi$ and $\Pi$-$\Delta$). These couplings are derived using selection rules and the pure precision approximation, allowing for population transfer between symmetry manifolds that was previously neglected.
3. Dual Representations: Calculations are performed in both adiabatic and strictly diabatic representations.
   • In the adiabatic representation, non-adiabatic couplings are treated via second-order derivative terms.
   • In the diabatic representation, the authors apply a strict adiabatic-to-diabatic transformation of the 23 coupled states. This approach preserves the character of the states (covalent vs. ion-pair) across avoided crossings.
   • For RIP formation, two specific diabatization schemes are compared: a two-state transformation (Method 1) and the full 23-state transformation (Method 2).

The initial wave packet is constructed from the vibrational ground state ($v=0$) of the HeH+ ion. Propagation is carried out using a Crank-Nicholson algorithm with a local complex potential (the "boomerang" model) to account for autoionization widths ($\Gamma$). Reaction cross sections are computed via half-Fourier transforms of the wave packet in the asymptotic region. Thermal rate coefficients are subsequently derived by averaging these cross sections over a Maxwell-Boltzmann distribution of electron energies.

Key Contributions and Results
The study yields several significant findings regarding the dynamics of HeH+ isotopologues ($^3$HeH+, $^4$HeH+, $^3$HeD+, $^4$HeD+, $^3$HeT+, $^4$HeT+):

• Enhanced DR Cross Sections: The inclusion of the large 23-state manifold and rotational couplings results in DR cross sections that are systematically larger than those reported in earlier theoretical works (e.g., Larson et al., Sarpal et al.). The new results show better agreement with the upper limits of experimental data from storage-ring experiments.
• Symmetry Contributions:
  • $^2\Sigma$ Dominance: In the diabatic representation, $^2\Sigma$ states dominate the recombination dynamics across most of the energy range (0–50 eV), attributed to strong coupling with the ionic ground state and favorable Franck-Condon factors.
  • Role of $^2\Pi$: At lower collision energies, $^2\Pi$ states contribute significantly, particularly in the adiabatic representation, due to rotationally induced mixing between $\Sigma$ and $\Pi$ manifolds.
  • $^2\Delta$ Contribution: The contribution of $^2\Delta$ states remains relatively small, as they are only indirectly connected to the ground state via rotational interactions.
• Resonant Ion-Pair (RIP) Formation: Both diabatization schemes predict RIP cross sections significantly larger than previous models, especially at low collision energies. Method 2 (strict 23-state diabatization) yields larger cross sections than Method 1 (two-state), indicating that the cumulative effect of multiple avoided crossings is critical for efficient population transfer to the ion-pair channel.
• Isotopic Effects: A clear inverse dependence of cross section magnitude on reduced mass is observed. Lighter isotopologues exhibit larger cross sections and thermal rate coefficients. This is attributed to faster nuclear motion in lighter species, which increases the probability of reaching the dissociation region before autoionization occurs.
• Thermal Rate Coefficients: The calculated thermal rate coefficients decrease monotonically with increasing electron temperature, consistent with experimental trends. The results fall within the envelope of rotationally resolved data (Novotný et al.) for rotational levels $j=0$ to $j \ge 3$, suggesting that the dominant physics is captured by the electronic dynamics and rotational couplings included in the model.

Significance and Claims
The authors assert that their work underscores the critical importance of multi-state coupling and non-adiabatic effects in accurately modeling electron–molecule collisions. By extending the electronic-state manifold and explicitly including rotational couplings, the study provides a more complete description of the available dissociation pathways.

The findings have specific implications for astrophysics and plasma modeling:

• Astrophysical Modeling: Since HeH+ is a key species in primordial gas chemistry and planetary nebulae (e.g., NGC 7027), the larger DR rate coefficients predicted here suggest that HeH+ may be destroyed more efficiently than assumed in earlier astrochemical models. This implies that updated abundance calculations may be necessary to interpret astronomical observations and understand the chemical evolution of the early Universe.
• Methodological Benchmark: The study serves as a stringent test case for modern wave-packet methodologies, demonstrating that accurate modeling of HeH+ requires a rigorous treatment of the interplay between electronic structure, non-adiabatic couplings, and rotational interactions.

The authors acknowledge that while their local complex potential approach captures the dominant mechanisms, it neglects interference effects between overlapping resonances and indirect recombination via high-lying Rydberg states. They suggest that future rigorous approaches based on MQDT or modern R-matrix methodologies could provide further benchmarks, but affirm that the current model successfully captures the primary dynamics governing the fragmentation of HeH+ isotopologues.