Warning
This documentation is for an old version of Cantera. You can find docs for newer versions here.
Cantera supports a number of different types of reactions, including several types of homogeneous reactions, surface reactions, and electrochemical reactions. For each, there is a corresponding entry type. The simplest entry type is reaction, which can be used for any homogeneous reaction that has a rate expression that obeys the law of mass action, with a rate coefficient that depends only on temperature.
All of the entry types that define reactions share some common features. These are described first, followed by descriptions of the individual reaction types in the following sections.
The reaction equation determines the reactant and product stoichiometry. A relatively simple parsing strategy is currently used, which assumes that all coefficient and species symbols on either side of the equation are delimited by spaces:
2 CH2 <=> CH + CH3 # OK
2 CH2<=>CH + CH3 # OK
2CH2 <=> CH + CH3 # error
CH2 + CH2 <=> CH + CH3 # OK
2 CH2 <=> CH+CH3 # error
The incorrect versions here would generate “undeclared species” errors and would halt processing of the input file. In the first case, the error would be that the species 2CH2 is undeclared, and in the second case it would be species CH+CH3.
Whether the reaction is reversible or not is determined by the form of the equality sign in the reaction equation. If either <=> or = is found, then the reaction is regarded as reversible, and the reverse rate will be computed from detailed balance. If, on the other hand, => is found, the reaction will be treated as irreversible.
The rate coefficient is specified with an embedded entry corresponding to the rate coefficient type. At present, the only implemented type is the modified Arrhenius function
which is defined with an Arrhenius entry:
rate_coeff = Arrhenius(A=1.0e13, n=0, E=(7.3, 'kcal/mol'))
rate_coeff = Arrhenius(1.0e13, 0, (7.3, 'kcal/mol'))
As a shorthand, if the rate_coeff field is assigned a sequence of three numbers, these are assumed to be \((A, n, E)\) in the modified Arrhenius function:
rate_coeff = [1.0e13, 0, (7.3, 'kcal/mol')] # equivalent to above
The units of the preexponential factor A can be specified explicitly if desired. If not specified, they will be constructed using the quantity, length, and time units specified in the units directive. Since the units of A depend on the reaction order, the units of each reactant concentration (different for bulk species in solution, surface species, and pure condensedphase species), and the units of the rate of progress (different for homogeneous and heterogeneous reactions), it is usually best not to specify units for A, in which case they will be computed taking all of these factors into account.
Note: if \(n \ne 0\), then the term \(T^n\) should have units of \(K^n\), which would change the units of A. This is not done, however, so the units associated with A are really the units for \(k_f\) . One way to formally express this is to replace \(T^n\) by the nondimensional quantity \([T/(1 K)]^n\).
An optional identifying string can be entered in the ID field, which can then be used in the reactions field of a phase or interface entry to identify this reaction. If omitted, the reactions are assigned ID strings as they are read in, beginning with '0001', '0002', etc.
Note that the ID string is only used when selectively importing reactions. If all reactions in the local file or in an external one are imported into a phase or interface, then the reaction ID field is not used.
Certain conditions are normally flagged as errors by Cantera. In some cases, they may not be errors, and the options field can be used to specify how they should be handled.
Normally, when a reaction is imported into a phase, it is checked to see that it is not a duplicate of another reaction already present in the phase, and an error results if a duplicate is found. But in some cases, it may be appropriate to include duplicate reactions, for example if a reaction can proceed through two distinctly different pathways, each with its own rate expression. Another case where duplicate reactions can be used is if it is desired to implement a reaction rate coefficient of the form:
While Cantera does not provide such a form for reaction rates, it can be implemented by defining N duplicate reactions, and assigning one rate coefficient in the sum to each reaction. If the 'duplicate' option is specified, then the reaction not only may have a duplicate, it must. Any reaction that specifies that it is a duplicate, but cannot be paired with another reaction in the phase that qualifies as its duplicate generates an error.
The reaction entry is used to represent homogeneous reactions with pressureindependent rate coefficients and mass action kinetics. Examples of reaction entries that implement some reactions in the GRIMech 3.0 natural gas combustion mechanism [3] are shown below:
units(length = 'cm', quantity = 'mol', act_energy = 'cal/mol')
...
reaction( "O + H2 <=> H + OH", [3.87000E+04, 2.7, 6260])
reaction( "O + HO2 <=> OH + O2", [2.00000E+13, 0.0, 0])
reaction( "O + H2O2 <=> OH + HO2", [9.63000E+06, 2.0, 4000])
reaction( "O + HCCO <=> H + 2 CO", [1.00000E+14, 0.0, 0])
reaction( "H + O2 + AR <=> HO2 + AR", [7.00000E+17, 0.8, 0])
reaction( "HO2 + C3H7 <=> O2 + C3H8", [2.55000E+10, 0.255, 943])
reaction( "HO2 + C3H7 => OH + C2H5 + CH2O", [2.41000E+13, 0.0, 0])
A threebody reaction is a gasphase reaction of the form:
Here M is an unspecified collision partner that carries away excess energy to stabilize the AB molecule (forward direction) or supplies energy to break the AB bond (reverse direction).
Different species may be more or less effective in acting as the collision partner. A species that is much lighter than A and B may not be able to transfer much of its kinetic energy, and so would be inefficient as a collision partner. On the other hand, a species with a transition from its ground state that is nearly resonant with one in the AB* activated complex may be much more effective at exchanging energy than would otherwise be expected.
These effects can be accounted for by defining a collision efficiency \(\epsilon\) for each species, defined such that the forward reaction rate is
where
where \(C_k\) is the concentration of species k. Since any constant collision efficiency can be absorbed into the rate coefficient \(k_f(T)\), the default collision efficiency is 1.0.
A threebody reaction may be defined using the three_body_reaction entry. The equation string for a threebody reaction must contain an 'M' or 'm' on both the reactant and product sides of the equation.
Some examples from GRIMech 3.0 are shown below:
three_body_reaction( "2 O + M <=> O2 + M", [1.20000E+17, 1, 0],
" AR:0.83 C2H6:3 CH4:2 CO:1.75 CO2:3.6 H2:2.4 H2O:15.4 ")
three_body_reaction( "O + H + M <=> OH + M", [5.00000E+17, 1, 0],
efficiencies = " AR:0.7 C2H6:3 CH4:2 CO:1.5 CO2:2 H2:2 H2O:6 ")
three_body_reaction(
equation = "H + OH + M <=> H2O + M",
rate_coeff = [2.20000E+22, 2, 0],
efficiencies = " AR:0.38 C2H6:3 CH4:2 H2:0.73 H2O:3.65 "
)
As always, the field names are optional if the field values are entered in the declaration order.
A falloff reaction is one that has a rate that is firstorder in [M] at low pressure, like a threebody reaction, but becomes zeroorder in [M] as [M] increases. Dissociation / association reactions of polyatomic molecules often exhibit this behavior.
The simplest expression for the rate coefficient for a falloff reaction is the Lindemann form [2]:
In the lowpressure limit, this approaches \(k0{\rm [M]}\), and in the highpressure limit it approaches \(k_\infty\).
Defining the nondimensional reduced pressure:
The rate constant may be written as
More accurate models for unimolecular processes lead to other, more complex, forms for the dependence on reduced pressure. These can be accounted for by multiplying the Lindemann expression by a function \(F(T, P_r)\):
This expression is used to compute the rate coefficient for falloff reactions. The function \(F(T, P_r)\) is the falloff function, and is specified by assigning an embedded entry to the falloff field.
A widelyused falloff function is the one proposed by Gilbert et al. [1]:
The Troe directive requires specifying the first three parameters \((A, T_3, T_1)\). The fourth parameter, \(T_2\), is optional, defaulting to 0.0.
This falloff function is based on the one originally due to Stewart et al. [4], which required three parameters \((a, b, c)\). Kee et al. [5] generalized this function slightly by adding two more parameters \((d, e)\). (The original form corresponds to \(d = 1, e = 0\).) Cantera supports the extended 5parameter form, given by:
In keeping with the nomenclature of Kee et al. [5], we will refer to this as the “SRI” falloff function. It is implemented by the SRI directive.
For these reactions, the rate falls off as the pressure increases, due to collisional stabilization of a reaction intermediate. Example:
which competes with:
Like falloff reactions, chemicallyactivated reactions are described by blending between a “low pressure” and a “high pressure” rate expression. The difference is that the forward rate constant is written as being proportional to the low pressure rate constant:
and the optional blending function F may described by any of the parameterizations allowed for falloff reactions. Chemicallyactivated reactions can be defined using the chemically_activated_reaction directive.
An example of a reaction specified with this parameterization:
chemically_activated_reaction('CH3 + OH (+ M) <=> CH2O + H2 (+ M)',
kLow=[2.823201e+02, 1.46878, (3270.56495, 'cal/mol')],
kHigh=[5.880000e14, 6.721, (3022.227, 'cal/mol')],
falloff=Troe(A=1.671, T3=434.782, T1=2934.21, T2=3919.0))
In this example, the units of \(k_0\) (kLow) are m^3/kmol/s and the units of \(k_\infty\) (kHigh) are 1/s.
The pdep_arrhenius class represents pressuredependent reaction rates by logarithmically interpolating between Arrhenius rate expressions at various pressures. Given two rate expressions at two specific pressures:
The rate at an intermediate pressure \(P_1 < P < P_2\) is computed as
Multiple rate expressions may be given at the same pressure, in which case the rate used in the interpolation formula is the sum of all the rates given at that pressure. For pressures outside the given range, the rate expression at the nearest pressure is used.
An example of a reaction specified in this format:
pdep_arrhenius('R1 + R2 <=> P1 + P2',
[(0.001315789, 'atm'), 2.440000e+10, 1.04, 3980.0],
[(0.039473684, 'atm'), 3.890000e+10, 0.989, 4114.0],
[(1.0, 'atm'), 3.460000e+12, 0.442, 5463.0],
[(10.0, 'atm'), 1.720000e+14, 0.01, 7134.0],
[(100.0, 'atm'), 7.410000e+30, 5.54, 12108.0],
[(100.0, 'atm'), 1.900000e+15, 0.29, 8306.0])
The first argument is the reaction equation. Each subsequent argument is a sequence of four elements specifying a pressure and the Arrhenius parameters at that pressure.
Class chebyshev_reaction represents a phenomenological rate coefficient \(k(T,P)\) in terms of a bivariate Chebyshev polynomial. The rate constant can be written as:
where \(\alpha_{tp}\) are the constants defining the rate, \(\phi_n(x)\) is the Chebyshev polynomial of the first kind of degree \(n\) evaluated at \(x\), and
are reduced temperature and reduced pressures which map the ranges \((T_\mathrm{min}, T_\mathrm{max})\) and \((P_\mathrm{min}, P_\mathrm{max})\) to \((1, 1)\).
A Chebyshev rate expression is specified in terms of the coefficient matrix \(\alpha\) and the temperature and pressure ranges. An example of a Chebyshev rate expression where \(N_T = 6\) and \(N_P = 4\) is:
chebyshev_reaction('R1 + R2 <=> P1 + P2',
Tmin=290.0, Tmax=3000.0,
Pmin=(0.001, 'atm'), Pmax=(100.0, 'atm'),
coeffs=[[1.44280e+01, 2.59970e01, 2.24320e02, 2.78700e03],
[ 2.20630e+01, 4.88090e01, 3.96430e02, 5.48110e03],
[2.32940e01, 4.01900e01, 2.60730e02, 5.04860e03],
[2.93660e01, 2.85680e01, 9.33730e03, 4.01020e03],
[2.26210e01, 1.69190e01, 4.85810e03, 2.38030e03],
[1.43220e01, 7.71110e02, 1.27080e02, 6.41540e04]])
Note that the Chebyshev polynomials are not defined outside the interval \((1,1)\), and therefore extrapolation of rates outside the range of temperatures and pressure for which they are defined is strongly discouraged.
References
[1]  R. G. Gilbert, K. Luther, and J. Troe. Ber. Bunsenges. Phys. Chem., 87:169, 1983. 
[2] 

[3]  Gregory P. Smith, David M. Golden, Michael Frenklach, Nigel W. Moriarty, Boris Eiteneer, Mikhail Goldenberg, C. Thomas Bowman, Ronald K. Hanson, Soonho Song, William C. Gardiner, Jr., Vitali V. Lissianski, , and Zhiwei Qin. GRIMech version 3.0, 1997. see http://www.me.berkeley.edu/gri_mech. 
[4]  P. H. Stewart, C. W. Larson, and D. Golden. Combustion and Flame, 75:25, 1989. 
[5]  (1, 2) R. J. Kee, F. M. Rupley, and J. A. Miller. ChemkinII: A Fortran chemical kinetics package for the analysis of gasphase chemical kinetics. Technical Report SAND898009, Sandia National Laboratories, 1989. 