CAMP: Special CMAQ Reaction for H2O2

For the reactions:

\[\ce{ HO2 + HO2 ->[k] H2O2 }\]

\[\ce{ HO2 + HO2 + H2O ->[k] H2O2 }\]

CMAQ rate constants are calculated as:

\[ k=k_1+k_2[\mbox{M}] \]

where \(k_1\) and \(k_2\) are Arrhenius rate constants with \(D=300\) and \(E=0\), and \([\mbox{M}]\) is the concentration of air [Gipson].

Input data for CMAQ \(\ce{H2O2}\) reactions have the following format:

{
  "type" : "CMAQ_H2O2",
  "k1_A" : 5.6E-12,
  "k1_B" : -1.8,
  "k1_C" : 180.0,
  "k2_A" : 3.4E-12,
  "k2_B" : -1.6,
  "k2_C" : 104.1,
  "time unit" : "MIN",
  "reactants" : {
    "spec1" : {},
    "spec2" : { "qty" : 2 },
    ...
  },
  "products" : {
    "spec3" : {},
    "spec4" : { "yield" : 0.65 },
    ...
  }
}

The key-value pairs reactants, and products are required. Reactants without a qty value are assumed to appear once in the reaction equation. Products without a specified yield are assumed to have a yield of 1.0.

The two sets of parameters beginning with k1_ and k2_ are the Arrhenius parameters for the \(k_1\) and \(k_2\) rate constants, respectively. When not present, _A parameters are assumed to be 1.0, _B to be 0.0, and _C to be 0.0.

The unit for time is assumed to be s, but inclusion of the optional key-value pair time unit = MIN can be used to indicate a rate with min as the time unit.