Chemical sensitivity to the ratio of the cosmic-ray ionization rates of He and H2 in dense clouds

Valentine Wakelam,Eric Herbst,Franck Selsis,Gérard Massacrier

DOI: https://doi.org/10.1051/0004-6361%3A20065472

2006-08-25

Abstract:Aim: To determine whether or not gas-phase chemical models with homogeneous and time-independent physical conditions explain the many observed molecular abundances in astrophysical sources, it is crucial to estimate the uncertainties in the calculated abundances and compare them with the observed abundances and their uncertainties. Non linear amplification of the error and bifurcation may limit the applicability of chemical models. Here we study such effects on dense cloud chemistry. Method: Using a previously studied approach to uncertainties based on the representation of rate coefficient errors as log normal distributions, we attempted to apply our approach using as input a variety of different elemental abundances from those studied previously. In this approach, all rate coefficients are varied randomly within their log normal (Gaussian) distribution, and the time-dependent chemistry calculated anew many times so as to obtain good statistics for the uncertainties in the calculated abundances. Results: Starting with so-called ``high-metal'' elemental abundances, we found bimodal rather than Gaussian like distributions for the abundances of many species and traced these strange distributions to an extreme sensitivity of the system to changes in the ratio of the cosmic ray ionization rate zeta\_He for He and that for molecular hydrogen zeta\_H2. The sensitivity can be so extreme as to cause a region of bistability, which was subsequently found to be more extensive for another choice of elemental abundances. To the best of our knowledge, the bistable solutions found in this way are the same as found previously by other authors, but it is best to think of the ratio zeta\_He/zeta\_H2 as a control parameter perpendicular to the ''standard'' control parameter zeta/n\_H.

Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **To determine whether the gas - phase chemical model can explain the multiple molecular abundances observed in astrophysical sources under uniform and time - invariant physical conditions, and to study the extreme sensitivity of these models to the ratio of cosmic - ray ionization rates of helium (He) and molecular hydrogen (H₂) (\(\frac{\zeta_{He}}{\zeta_{H2}}\)) and its impact on the chemical abundance distribution**.

### Specific problem description

1. **Model applicability assessment**:
- Researchers need to assess whether the gas - phase chemical model can explain the observed molecular abundances under uniform and time - invariant physical conditions.
- By comparing the calculated molecular abundances with the observational results and their uncertainties, researchers hope to understand the applicability and limitations of the model.

2. **Non - linear amplification effect and bifurcation phenomenon**:
- Non - linear amplification effect and bifurcation phenomenon may limit the application of chemical models. Researchers are concerned about how these effects affect the chemistry of dense clouds.
- In particular, researchers have found that the abundance distributions of some species show a bimodal rather than a Gaussian distribution, indicating that the system is very sensitive to changes in certain parameters.

3. **Impact of the ratio of cosmic - ray ionization rates of helium and molecular hydrogen**:
- Researchers are particularly concerned about the impact of the ratio of cosmic - ray ionization rates of helium (He) and molecular hydrogen (H₂) \(\frac{\zeta_{He}}{\zeta_{H2}}\) on chemical abundances.
- It has been found that changes in this ratio can lead to the bistability phenomenon of the system, that is, there are two stable chemical abundance solutions under the same set of parameters.

### Research methods

- **Monte Carlo method**:
- Use the Monte Carlo method to calculate the theoretical error caused by the uncertainties of rate coefficients and physical conditions.
- The uncertainty of the rate coefficient is represented by a log - normal distribution, and all rate coefficients randomly change within their respective log - normal distributions.

- **Parameter variation analysis**:
- Study the impact of different elemental abundances (such as high - metal, low - metal and intermediate - metal abundances) on sensitivity.
- Analyze the impact of changes in the \(\frac{\zeta_{He}}{\zeta_{H2}}\) ratio on the abundances of various species, especially oxygen (O₂), carbon (C), water (H₂O), electrons (e⁻) and H₃⁺ ions, etc.

### Results and conclusions

- **Bimodal distribution and bistability phenomenon**:
- It has been found that the abundance distributions of many species show a bimodal rather than a Gaussian distribution, indicating that the system is extremely sensitive to changes in \(\frac{\zeta_{He}}{\zeta_{H2}}\).
- This sensitivity may lead to the bistability phenomenon, that is, there are two stable chemical abundance solutions under the same set of parameters.

- **Importance of control parameters**:
- \(\frac{\zeta_{He}}{\zeta_{H2}}\) is considered an important control parameter, similar to the standard control parameter \(\frac{\zeta}{n_H}\).
- Different choices of elemental abundances will affect the manifestation and degree of sensitivity, but \(\frac{\zeta_{He}}{\zeta_{H2}}\) is always a key factor.

Through these studies, the author aims to reveal the hidden complexity and non - linear behavior in the dense - cloud chemical model and provide new perspectives and tools for future astrochemical research.