Cosmic rays may be linked to the formation of volatiles necessary for prebiotic chemistry. We explore the effect of cosmic rays in a hydrogen-dominated atmosphere, as a proof-of-concept that ion-neutral chemistry may be important for modelling hydrog
en-dominated atmospheres. In order to accomplish this, we utilize Monte Carlo cosmic ray transport models with particle energies of $10^6$ eV $< E < 10^{12}$ eV in order to investigate the cosmic ray enhancement of free electrons in substellar atmospheres. Ion-neutral chemistry is then applied to a Drift-Phoenix model of a free-floating giant gas planet. Our results suggest that the activation of ion-neutral chemistry in the upper atmosphere significantly enhances formation rates for various species, and we find that C$_2$H$_2$, C$_2$H$_4$, NH$_3$, C$_6$H$_6$ and possibly C$_{10}$H are enhanced in the upper atmospheres because of cosmic rays. Our results suggest a potential connection between cosmic ray chemistry and the hazes observed in the upper atmospheres of various extrasolar planets. Chemi-ionization reactions are briefly discussed, as they may enhance the degree of ionization in the cloud layer.
Context: Observations of small carbon-bearing molecules such as CCH, C4H, c-C3H2, and HCO in the Horsehead Nebula have shown these species to have higher abundances towards the edge of the source than towards the center. Aims: Given the determinati
on of a wide range of values for zeta (s-1), the total ionization rate of hydrogen atoms, and the proposal of a column-dependent zeta(N_H), where N_H is the total column of hydrogen nuclei, we desire to determine if the effects of zeta(N_H) in a single object with spatial variation can be observable. We chose the Horsehead Nebula because of its geometry and high density. Method: We model the Horsehead Nebula as a near edge-on photon dominated region (PDR), using several choices for zeta, both constant and as a function of column. The column-dependent zeta functions are determined by a Monte Carlo model of cosmic ray penetration, using a steep power-law spectrum and accounting for ionization and magnetic field effects. We consider a case with low-metal elemental abundances as well as a sulfur-rich case. Results: We show that use of a column-dependent zeta(N_H) of 5(-15) s-1 at the surface and 7.5(-16) s-1 at Av = 10 on balance improves agreement between measured and theoretical molecular abundances, compared with constant values of zeta.
