Gas phase ODEs
The rate of change of the concentration of a gas phase species due to a single two body reaction is
where is the rate of that reaction in units of . Since we work in fractional abundances rather than concentrations, we can remove factors of since
For reactions between involving only a single body such as ionization by a cosmic ray, we have
The total rate of change of the fractional abundance of a species due to gas phase reactions is then just the sum of these terms for all reactions where it is a product minus the sum of all reactions where it is a reactant.
As a rule, any part of a reaction ODE which does not depend on abundance (eg the rate itself) is calculated between timesteps by the subroutine
rates.f90. The abundances are include in the ODE calculation itself so they can be updated between steps by the solver.
Gas phase chemistry in UCLCHEM uses the UMIST12 database. This is a database listing reactants and products with up to three rate constants which we label for thousands of gas phase reactions. We briefly list here the way in which the rates in the above equations are calculated for each reaction type and more information can be found in McElroy et al. 2013
Two Body Reactions use the Kooji-Arrhenius equation.
Cosmic Ray Protons
Cosmic Ray induced photons
where is the cosmic ray ionization rate in units of 1.3 10 s, E is the efficiency with which cosmic rays cause ionization, is the dust grain albedo, is the attenuated UV field.