MEASURE

4. measure.collision - Collision models and parameters

Contains classes that represent the collision models available in MEASURE, and methods for calculating various collision parameters. Each collision model provides a method generateCollisionMatrix() that generates the collision matrix \matrix{M}_\mathrm{coll} / \omega = \matrix{P} - \matrix{I} corresponding to the collisional energy transfer probability function P(E, E^\prime) for that model. The available collision models are:

4.1. Collision Models

class measure.collision.CollisionModel

A base class for collision models. To create a custom collision model, derive from this class and implement the generateCollisionMatrix() method, which returns the collision matrix for the collision model you are implementing.

Note

As with all collision models, you can only specify either the deactivating direction or the activating direction of the collisional transfer probabilities function P(E, E^\prime), as the other is constrained by detailed balance.

4.1.1. Single Exponential Down Model

class measure.collision.SingleExponentialDownModel

Refactoring of collision and reaction modules to full Cython syntax. A single exponential down collision model, based around the collisional energy transfer probability function

P(E, E^\prime) = C(E^\prime) \exp \left( - \frac{E^\prime - E}{\alpha} \right) \hspace{40pt} E < E^\prime

where the parameter \alpha = \left< \Delta E_\mathrm{d} \right> represents the average energy transferred in a deactivating collision. This is the most commonly-used collision model, simply because it only has one parameter to determine. The parameter \alpha is specified using the equation

\alpha = \alpha_0 \left( \frac{T}{T_0} \right)^n

where \alpha_0 is the value of \alpha at temperature T_0 in K. Set the exponent n to zero to obtain a temperature-independent value for \alpha.

Attribute Type Description
alpha double The average energy transferred in a deactivating collision in J/mol
generateCollisionMatrix

Generate and return the collision matrix \matrix{M}_\mathrm{coll} / \omega = \matrix{P} - \matrix{I} corresponding to this collision model for a given set of energies Elist in J/mol, temperature T in K, and isomer density of states densStates.

getAlpha

Return the value of the \alpha parameter at temperature T in K. The \alpha parameter represents the average energy transferred in a deactivating collision \left< \Delta E_\mathrm{d} \right>, and has units of J/mol.

4.2. Collision Parameters

measure.collision.calculateCollisionFrequency()

Calculate the Lennard-Jones collision frequency for a given species with a dictionary of bath gases and their mole fractions bathGas at a given temperature T in K and pressure P in Pa. The Lennard-Jones model is generally a slight underestimate, but reasonable enough. If the bath gas is a mixture, arithmetic means are used to compute its effective Lennard-Jones \sigma parameter and molecular weight, while a geometric mean is used to calculate its effective Lennard-Jones \epsilon parameter.

measure.collision.calculateCollisionEfficiency()

Calculate an efficiency factor for collisions, particularly useful for the modified strong collision method. The collisions involve the given species with density of states densStates in mol/J corresponding to energies Elist in J/mol, ground-state energy E0 in J/mol, and first reactive energy Ereac in J/mol. The collisions occur at temperature T in K and are described by the collision model collisionModel, which currently must be a SingleExponentialDownModel object. The algorithm here is implemented as described by Chang, Bozzelli, and Dean [Chang2000].

[Chang2000]A. Y. Chang, J. W. Bozzelli, and A. M. Dean. Z. Phys. Chem. 214, p. 1533-1568 (2000). doi: 10.1524/zpch.2000.214.11.1533