[Paper Review] Reduction of chemical networks. II. Analysis of the fractional ionisation in protoplanetary discs

D. S. Semenov, D. S. Wiebe+1 more|ArXiv.org|Mar 23, 2004

Astrophysics and Star Formation Studies35 references94 citations

TL;DR

This paper develops reduced chemical networks to accurately model fractional ionisation in protoplanetary discs over 10⁶ years, revealing three distinct ionisation layers: a low-ionisation midplane dominated by cosmic rays and radionuclides (10 species, ~10 reactions), a complex intermediate layer driven by X-rays and surface hydrogenation of carbon chains (100+ species, hundreds of reactions), and a surface layer governed by UV photoionisation (20 species, ~20 reactions). The equilibrium approximation (xₑ = √(ζ/βn_H)) fails in dynamic regimes, especially in the intermediate layer, where non-equilibrium chemistry dominates.

ABSTRACT

(abridged) We analyse the evolution of the fractional ionisation in a steady-state protoplanetary disc with a vertical temperature gradient and with gas-grain chemistry including surface reactions. The ionisation due to stellar X-rays, stellar and interstellar UV radiation, cosmic rays and radionuclide decay is taken into account. Using our reduction schemes as a tool for the analysis, we isolate small sets of chemical reactions that reproduce the evolution of the ionisation degree at representative disc locations with an accuracy of 50%-100%. Column densities of key molecules are calculated and compared to the results of other recent studies and observational data.

Motivation & Objective

• To analyze the time-dependent fractional ionisation in a steady-state protoplanetary disc over 10⁶ years, accounting for multi-source ionisation and gas-grain chemistry.
• To identify and isolate minimal chemical networks that reproduce ionisation evolution with 50%–100% accuracy in different disc regions.
• To challenge the common assumption of ionisation equilibrium (xₑ = √(ζ/βn_H)) by demonstrating its inaccuracy in dynamic, non-equilibrium conditions, especially in the intermediate layer.
• To assess the relevance of surface reactions, particularly hydrogenation of long carbon chains, for ionisation control in the X-ray-irradiated intermediate layer.
• To provide reduced networks for use in MHD modelling, acknowledging limitations in dynamical contexts.

Proposed method

• Application of species-based network reduction techniques (from Paper I) to isolate key species and reactions that control ionisation evolution in representative disc locations.
• Use of the full UMIST 95 network with surface chemistry as the reference, including ionisation from stellar X-rays, UV, cosmic rays, and radionuclide decay.
• Construction of reduced networks with 10–120 species and corresponding reactions, selected to reproduce fractional ionisation within a factor of 2 over 10⁶ years.
• Incorporation of surface reactions on dust grains, particularly for hydrogenation of carbon chains, to model ionisation in the intermediate layer.
• Comparison of reduced network results with equilibrium ionisation (xₑ = √(ζ/βn_H)) and observational data on column densities.
• Evaluation of network performance across three disc layers: midplane (shielded), intermediate (X-ray-irradiated), and surface (UV-irradiated).

Experimental results

Research questions

• RQ1 How do different ionisation sources—cosmic rays, radionuclides, X-rays, and UV—contribute to fractional ionisation in distinct regions of a protoplanetary disc?
• RQ2 To what extent does the standard equilibrium ionisation formula xₑ = √(ζ/βn_H) accurately represent time-dependent ionisation in dynamic disc conditions?
• RQ3 What role do surface reactions, particularly hydrogenation of long carbon chains, play in sustaining ionisation in the X-ray-irradiated intermediate layer?
• RQ4 Why do reduced networks in the intermediate layer require over 100 species and hundreds of reactions, while midplane and surface layers require far fewer?
• RQ5 Can reduced networks be reliably used in MHD simulations of disc evolution, and what are the limitations in dynamically evolving systems?

Key findings

• The midplane, shielded from high-energy radiation, maintains a very low fractional ionisation (≤10⁻¹²) sustained only by cosmic rays and radionuclides, requiring only ~10 species and ~10 reactions for accurate modelling.
• In the intermediate layer, X-rays dominate ionisation, but the chemistry is highly complex due to the formation and reactivity of long carbon chains (up to C6), necessitating reduced networks with over 100 species and hundreds of reactions.
• Surface hydrogenation of carbon chains is identified as a crucial process for ionisation control in the intermediate layer, a factor previously overlooked in simplified models.
• The equilibrium ionisation formula xₑ = √(ζ/βn_H) yields values within a factor of 2 of the true ionisation in the midplane and surface layer, but fails in the intermediate layer where non-equilibrium dynamics dominate.
• In the surface layer, ionisation is primarily governed by C⁺ photoionisation and recombination, requiring a recombination coefficient of 1.4×10⁻¹³T⁻⁰.⁶¹, not the standard β value, and equilibrium is reached in less than 100 years.
• Reduced networks are not directly transferable to dynamical MHD models due to their steady-state derivation, but merging them into a single 120-species network could enable future use in time-dependent simulations.