High-Energy Astrochemistry

• High-energy astrochemistry is defined by the study of molecular transformations driven by ionizing radiation and energetic particles, enabling non-thermal chemistry in both gas and ice phases.
• It investigates processes such as cosmic-ray, X-ray, UV, and shock-induced ionization, where secondary electron cascades trigger extensive chemical networks beyond conventional thermal reactions.
• Integrated models coupling radiative transfer with chemical networks highlight key tracers (e.g., H3+ and HCO+) and validate laboratory experiments on ice radiolysis, informing our understanding of molecular evolution in varied astronomical settings.

High-energy astrochemistry is the study of molecular formation, destruction, and transformation under the influence of ionizing radiation and energetic particles—cosmic rays, X-rays, fast electrons, gamma rays, UV photons, and shocks—through processes fundamentally distinct from thermal or UV-driven chemistry. The unifying criterion is that each irradiating particle or photon deposits sufficient energy to generate a cascade of low-energy secondary electrons (with $N_e \gg 1$), catalyzing extensive non-thermal chemical networks in both gas and ice phases (Gaches et al., 10 Dec 2025). This regime controls critical pathways for molecular ionization, dissociation, heating, physical desorption, and the emergence of chemical complexity in environments ranging from cold prestellar cores to galactic nuclei, and is central to explaining observed abundances of complex organic molecules in regions not accessible to conventional thermal or surface chemistry.

1. Physical Processes and Energy Sources

High-energy astrochemistry is triggered by energetic agents with the ability to ionize and dissociate molecular species far from equilibrium:

• Cosmic rays (CRs): High-energy protons, electrons, and heavier ions accelerated in supernova remnants, colliding-wind binaries, or bow shocks of massive stars. With energies from MeV to PeV, CRs penetrate deep into clouds and ices, driving ionization, heating, and secondary electron cascades. The total ionization rate per H$_2$ molecule is $$\zeta_{\rm H_2} = \int J(E)\,\sigma_{\rm ion,H_2}(E) dE$$, with typical $\zeta_{\rm H_2} \sim 10^{-17} - 10^{-14}$ s$^{-1}$ depending on environment (Becker, 2013, Gaches et al., 10 Dec 2025).
• X-rays: Hard X-rays from AGN or massive YSOs ionize and dissociate atoms and molecules via K-shell photoionization, with secondary Auger electrons producing further ionizations. The primary photoionization rate for species $i$ is $$\zeta_{i,\rm pe} = \int \sigma_{i,\rm pe}(E)\frac{F_X(E)}{E}\,dE$$ (Gaches et al., 10 Dec 2025, Harada, 2017).
• UV photons: Far-UV (6–13.6 eV) fields from O/B stars drive photodissociation and photoionization in photon-dominated regions (PDRs). The unshielded photodissociation rate is $$k_{\rm pd,i} = \int_{v_0}^\infty (4\pi J_\nu/h\nu)\,\sigma_{\rm pd,i}(\nu)\,d\nu$$ (Becker, 2013).
• Shocks: Mechanical energy from supernovae or winds induces adiabatic or radiative shocks, promoting endothermic reactions and sputtering of dust/ices. Gas temperatures $T_2 \sim 10^3$–$10^4$ K enable neutral–neutral reactions with activation barriers (e.g., C + H$_2$ and O + H$_2$) (Becker, 2013, Harada, 2017).
• Radionuclide decay: Radioactive isotopes ($^{26}$Al, $^{60}$Fe) synthesized in massive stars decay, providing local ionization and heating, especially in embedded star/disk systems (Becker, 2013).

In all these regimes, the energy deposition results in a copious production of secondary electrons, which allocate energy via ionizations, excitations, and radical generation, enhancing chemical complexity especially in shielded and cold environments not directly accessible to UV photons (Gaches et al., 10 Dec 2025).

2. Chemistry in the Gas Phase: Ionization, Heating, and Nonthermal Pathways

Ionizing radiation initiates rich ion-molecule and radical chains, bypassing activation barriers that restrict neutral-neutral pathways at low temperatures:

• H$_3^+$ formation: The archetypal process, H$_2$ + CR $\rightarrow$ H$_2^+$ + e$^-$; H$_2^+$ + H$_2$ $\rightarrow$ H$_3^+$ + H. H$_3^+$ serves as a universal ionization and proton-transfer tracer (Becker, 2013, Gaches et al., 10 Dec 2025).
• Helium and carbon networks: CR or X-ray ionization of He produces He$^+$, which destroys CO: He$^+$ + CO $\rightarrow$ C$^+$ + O + He. Subsequent C$^+$ + H$_2$ and H$_3^+$ + C reactions lead to the formation of CH$^+$ and hydrocarbon cations (Gaches et al., 10 Dec 2025).
• Oxygen hydrogenation: CR/X-rays produce O$^+$, leading to chains: O$^+$ $\rightarrow$ OH$^+$ $\rightarrow$ H$_2$O$^+$ $\rightarrow$ H$_3$O$^+$, terminating in neutral water upon dissociative recombination (Gaches et al., 10 Dec 2025, Becker, 2013).
• NH$_3$ synthesis: N$^+$ produced by CR/X-rays is hydrogenated to NH$_3^+$, then to NH$_3$ via NH$_4^+$ dissociative recombination.
• Thermal balance: Each ionization deposits 10–20 eV as heat, leading to $$\Gamma_{\rm CR} \sim n_H \zeta_{\rm CR} \langle E_h \rangle$$ (Harada, 2017). This process maintains elevated gas temperatures, particularly in high-$\zeta$ environments such as starburst nuclei and CR-dominated regions (CRDRs).

Such non-thermal pathways dominate in environments where thermal activation is suppressed, explaining the presence of iCOMs (interstellar complex organic molecules) in cold ($T \sim 10$ K) gas (Gaches et al., 10 Dec 2025). The electron abundance is governed by $x_e \sim \sqrt{\zeta/\alpha n_H}$, with $\alpha$ the recombination rate.

3. Ice Radiolysis and Non-Thermal Surface Chemistry

Energetic radiation penetrates molecular ices, initiating a sequence of radiolytic transformations:

• Mechanisms: Ionization plus dissociative recombination, dissociative excitation, excitation to stable states, and desorption into gas phase (Gaches et al., 10 Dec 2025). Each interaction triggers a cascade of secondary electrons, leading to bond cleavage and radical mobilization.
• Radiolysis rates: The rate for a product $i$ induced by radiation $k$ is $$k_{i,k} = G_i\,\left(S_{e,k}/100~{\rm eV}\right)\,\phi_k$$, with $G_i$ the yield per 100 eV, $S_{e,k}$ the stopping power, and $\phi_k$ the flux (Gaches et al., 10 Dec 2025). Sputtering rates at high $S_e$ scale as $Y_{\rm sput} \sim Y_e^0 S_e^2$.
• Laboratory constraints: Experiments using keV–MeV electrons and ions demonstrate robust formation of simple and complex products such as CO$_2$, H$_2$CO, CH$_4$, HNCO, and carbon-chain oxides from H$_2$:CO: N$_2$ ices, with conversion yields of 5–10% for electrons, and yields 1–2 orders of magnitude lower for Ly-α photons due to inefficient absorption (Martin-Domenech et al., 2020). Methanol radiolysis yields CH$_4$, H$_2$CO, CO, CO$_2$, and radicals (CH$_2$OH, CH$_3$O, HCO) (Herczku et al., 2021, Mifsud et al., 2021).

Laboratory Parameter	Value or Range	Reference
Electron flux	$\sim 2.4 \times 10^{11}$ e$^-$ cm$^{-2}$ s$^{-1}$	(Martin-Domenech et al., 2020)
Product yields	5–10% (e$^-$) for CO$_2$, H$_2$CO, CH$_4$, etc.	(Martin-Domenech et al., 2020)
CH$_3$OH cross-section	$\sigma_{\rm CH_3OH} \sim 3 \times 10^{-18}$ cm$^2$	(Mifsud et al., 2021)

These experiments validate the inclusion of bulk ice radiolysis in astrochemical models, extending the surface hydrogenation paradigm to non-thermal contexts and yielding COMs in cold cores (Gaches et al., 10 Dec 2025).

4. Integrated Modeling Frameworks and Numerical Treatment

Modern modeling of high-energy astrochemistry entails coupling radiative transfer, physical dynamics, and chemical networks:

• General formulation: Each species $i$ evolves as

$$\frac{dn_i}{dt} = \sum_{j,k}k_{j,k} n_j n_k + \sum_j \Gamma_j n_j - n_i \sum_l k_{i,l} n_l - \alpha_i n_i n_e$$

with $\Gamma_j$ representing CR and X-ray induced rates (Gaches et al., 10 Dec 2025, Becker, 2013).

• Chemical networks: Reaction databases (UMIST, KIDA) are updated to include electron-impact, CR- and X-ray-driven reactions, with cross-sections from sources such as ALeCS (Gaches et al., 10 Dec 2025, Bisbas et al., 2023).
• Surface and ice processes: Codes such as Nautilus, UCLCHEM, 3D-PDR support radiolytic chemistry and cosmic-ray-induced desorption (Gaches et al., 10 Dec 2025, Becker, 2013).
• Physical-radiative coupling: Users solve the steady-state or time-dependent coupled equations for density, temperature, ionizing flux, and molecular abundances, including heating and cooling balances ($\Gamma_{\rm CR}$, $\Gamma_X$, $\Lambda_{\rm line}$, $\Lambda_{dg}$) (Harada, 2017).
• Treatment of gradients: Accurate propagation of cosmic-ray and X-ray fluxes in heterogeneous media requires consideration of attenuation, secondary electron distributions, and Monte Carlo approaches for electron cascade handling (Gaches et al., 10 Dec 2025).

These frameworks are essential for interpreting observables and linking high-energy feedback to molecular abundances.

5. Observational Diagnostics and Tracers

Distinct chemical signatures in line emission and molecular abundances provide empirical diagnostics of high-energy astrochemistry:

• Ionization tracers: H$_3^+$, HCO$^+$, HOC$^+$, OH$^+$, and H$_3$O$^+$. The HCO$^+$/HOC$^+$ ratio (range: 10–150 in NGC 253) and elevated HOC$^+$ abundance up to $6 \times 10^{-10}$ uniquely trace high CR ionization or intense UV in starburst nuclei (Harada et al., 2021).
• Molecular ions and transitions: Far-infrared lines of CH$^+$, OH$^+$, H$_2$O$^+$ (Herschel); CO SLEDs (J=1–0 to 10–9); [C II] 158 μm and C I emission; molecular ratios such as SiO/CH$_3$OH for shocks, CN/HCN for PDRs, HCN/HCO$^+$ for XDR/CRDR conditions (Harada, 2017, Becker, 2013, Bisbas et al., 2023).
• Abundance patterns: In $\alpha$-enhanced (C/O subsolar) and high-$\zeta_{\rm CR}$ environments, CO and C become underabundant and CO(1–0) remains the most robust H$_2$ tracer, while [C I] can become “dark” ([CI]-dark galaxies) (Bisbas et al., 2023).
• High-energy tracers in external galaxies: Enhanced abundances of H$_3$O$^+$, SO, and HCO$^+$ in starbursts, and CN (PDR), SiO (shock), and HC$_3$N (dense gas) in mergers and nuclear regions (Harada, 2017, Harada et al., 2021).

Observational constraints now routinely reach the spatial and spectral sensitivity needed to distinguish between PDR, XDR, CRDR, and MDR chemo-thermal regimes.

6. Environmental Variation: Starbursts, AGN, and Metallicity Effects

The balance of ionizing fluxes and chemical outcomes is highly dependent on environmental context:

• Starburst galaxies and CRDRs: Elevated $\zeta_{\rm CR}$ (up to $10^{-14}$–$10^{-13}$ s$^{-1}$), high UV ($G_0 \sim 10^4$), and shocks combine. For example, NGC 253 observations require $\zeta \gtrsim 10^{-14}$ s$^{-1}$ or $G_0 \gtrsim 10^3$ to reproduce observed HOC$^+$ enhancement and low [HCO$^+$]/[HOC$^+$] ratios (Harada et al., 2021).
• XDRs and AGN: X-ray luminosities up to $10^{45}$ erg s$^{-1}$ drive deep ionization and heating, uniform to large columns. Observed line ratios (e.g., HCN/HCO$^+$, doubly-ionized atoms such as C$^{++}$) and SLEDs differentiate AGN-induced chemistry from star formation–dominated regions (Harada, 2017).
• Low-metallicity and $\alpha$-enhanced systems: Suppressed carbon abundance ($[\rm C/O] < 0$) and high $\zeta_{\rm CR}$ environments produce strong differential effects, with CO and C phases diminished and CO(1–0) remaining the most stable molecular tracer (Bisbas et al., 2023). The [CII]/FIR deficit in ULIRGs is only partially explained by low [C/O], implicating additional thermal and radiative processes.

These insights are crucial for extragalactic surveys and modeling the ISM/CGM across cosmic time.

7. Frontier Questions and Developments

Several challenges and directions shape the field:

• Modeling of cosmic-ray and X-ray transport: Accurate inclusion of magnetic mirroring, propagation, and attenuation in heterogeneous media is needed. Current models extend beyond simple $\sigma j$ integrals (Gaches et al., 10 Dec 2025).
• Inclusion of secondary electron cascades: Complete treatment of electron-induced processes in both gas and ice phases is advancing with detailed cross-section databases (e.g., ALeCS) and on-the-spot/Monte Carlo methods (Gaches et al., 10 Dec 2025).
• Surface and ice radiolysis: General adoption of bulk radiolysis in chemical codes is required for capturing the buildup of complex molecules in cold environments, especially with mixed ices and at low temperatures.
• Microphysical modeling: ab initio and molecular-dynamics simulations (REAXFF, Geant4-DNA) are beginning to provide mechanistic underpinning for yields, branching ratios, and the influence of ice structure (Gaches et al., 10 Dec 2025).
• Non-thermal excitation and diagnostics: Emission from H$_2$ lines excited by secondary UV (Prasad–Tarafdar effect) is becoming an incisive probe of $\zeta$, but is not yet widely incorporated in chemical models.
• Origins of chirality and spin chemistry: Laboratory and theoretical studies are poised to test whether spin-polarized secondary electrons from grains could produce chiral asymmetry, potentially linking high-energy astrochemistry to the emergence of homochirality (Gaches et al., 10 Dec 2025).

With the synergy of JWST/ALMA observations, next-generation experiments (INFRA-ICE, ICA, VIZSLA), and computational modeling, the full impact of high-energy astrochemistry on the molecular evolution of the universe is becoming accessible.