CAMP: Special CMAQ Reaction for OH+HNO3

For the reaction:

\[\ce{ OH + HNO3 ->[k] NO3 + H2O }\]

CMAQ rate constants are calculated as:

\[ k=k_0+(\frac{k_3[\mbox{M}]}{1+k_3[\mbox{M}]/k_2}) \]

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

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

{
  "type" : "CMAQ_OH_HNO3",
  "k0_A" : 5.6E-12,
  "k0_B" : -1.8,
  "k0_C" : 180.0,
  "k2_A" : 3.4E-12,
  "k2_B" : -1.6,
  "k2_C" : 104.1,
  "k3_A" : 3.2E-11,
  "k3_B" : -1.5,
  "k3_C" : 92.0,
  "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 three sets of parameters beginning with k0_, k2_, and k3_, are the Arrhenius parameters for the \(k_0\), \(k_2\) and \(k_3\) 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.