Calculate the chemical affinities of formation reactions of species.
numeric, zero or more named arguments, used to identify the variables of interest in the calculations
character, denoting the property to be calculated. Default is A, for chemical affinity of formation reactions of species of interest
list, output from
character. This argument is used to identify a conserved basis species (or PBB) in a chemical activity buffer. Default is PBB.
numeric, indices of proteins in
numeric, logarithms of activities of proteins identified in
affinity calculates the chemical affinities of reactions to form the species of interest from the basis species.
The equation used to calculate chemical affinity A is A=RT*ln(K/Q), where K denotes the equilibrium constant of the reaction and Q stands for the activity product of the species in the reaction.
The calculation of chemical affinities relies on the current definitions of the
basis species and
species of interest.
Calculations are possible at single values of temperature, pressure, ionic strength and chemical activities of the basis species, or as a function of one or more of these variables.
Zero, one, or more leading arguments to the function identify which of the chemical activities of basis species, temperature, pressure and/or ionic strength to vary.
The names of each of these arguments may be the formula of any of the basis species of the system, or T, P, pe, pH, Eh, or IS (but names may not be repeated).
The names of charged basis species such as K+ and SO4-2 should be quoted when used as arguments.
The value of each argument is of the form
c(min, max) or
c(min, max, res) where
max refer to the minimimum and maximum values of variable identified by the name of the argument, and
res denotes the resolution, or number of points along which to do the calculations;
red is assigned a default value of 128 if it is missing.
For any arguments that refer to basis species, the numerical values are the logarithms of activity (or fugacity for gases) of that basis species.
If T, P, and/or IS are not among the
vars, their constant values can be supplied in
IS (in mol kg^-1).
The units of T and P are those set by
P.units (on program start-up these are °C and bar, respectively).
sout, if provided, replaces the call to
subcrt, which can greatly speed up the calculations if this intermediate step is stored by other functions (e.g.,
exceed.Ttr is passed to
subcrt so that the properties of mineral phases beyond their transition temperatures can optionally be calculated.
If one or more buffers are assigned to the definition of
basis species, the logarithms of activities of these basis species are taken from the buffer (see
loga.protein arguments can be used to compute the chemical affinities of formation reactions of proteins that are not in the current
iprotein contains the indices (rownumbers) of desired proteins in
This uses some optimizations to calculate the properties of many proteins in a fraction of the time it would take to calculate them individually.
When the length(s) of the variables is(are) greater than 3, the function enters the transect mode of operation. In this mode of operation, instead of performing the calculations on an n-dimensional grid, the affinities are calculated on a transect of changing T, P, and/or chemical activity of basis species.
The function can also be used to calculate other thermodynamic properties of formation reactions.
Valid properties are A or NULL for chemical affinity, logK or logQ for logarithm of equilibrium constant and reaction activity product, or any of the properties available in
subcrt except for rho.
The properties returned are those of the formation reactions of the species of interest from the basis species.
It is also possible to calculate the properties of the species of interest themselves (not their formation reactions) by setting the
property to G.species, Cp.species, etc.
Except for A, the properties of proteins or their reactions calculated in this manner are restricted to nonionized proteins.
A list, elements of which are
sout output from
property name of the calculated property (A for chemical affinity),
species definition of basis species and species of interest in effect at runtime,
P temperature and pressure, in the system units of Kelvin and bar, set to
numeric() (length=0) if either one is a variable,
vars the names of the variables,
vals the values of the variables (a list, one element for each variable),
values the result of the calculation (a list, one element for each species, with names taken from the species index in
The elements of the lists in
values are arrays of n dimensions, where n is the number of variables.
The values of chemical affinity of formation reactions of the species are returned in dimensionless units (for use with decimal logarithms, i.e., A/2.303RT).
Names other than T or P in
vars generally refer to basis species, and the corresponding
vals are the logarithms of activity or fugacity.
However, if one or more of
pH is among the variables of interest,
vals holds the values of the those variables as indicated.
Helgeson, H. C., Richard, L, McKenzie, W. F., Norton, D. L. and Schmitt, A. (2009) A chemical and thermodynamic model of oil generation in hydrocarbon source rocks. Geochim. Cosmochim. Acta 73, 594–695. http://dx.doi.org/10.1016/j.gca.2008.03.004
ionize.aa, activated if proteins are among the species of interest, H+ is in the basis, and
thermo$opt$ionize is TRUE.
equilibrate for using the results of
affinity to calculate equilibrium activities of species, and
diagram to plot the results.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## set up a system and calculate ## chemical affinities of formation reactions basis(c("SiO2", "MgO", "H2O", "O2"), c(-5, -5, 0, 999)) species(c("quartz","enstatite","forsterite")) # chemical affinities (A/2.303RT) at 25 deg C and 1 bar affinity() # at higher temperature and pressure affinity(T=500, P=2000) # at 25 temperatures and pressures affinity(T=c(500, 1000, 5), P=c(1000, 5000, 5)) # equilibrium constants of formation reactions affinity(property="logK") # standard molal Gibbs energies of species, # user units (default: cal/mol) affinity(property="G.species") # standard molal Gibbs energies of reactions affinity(property="G") # a T,P-transect # (fluid pressure from Helgeson et al., 2009 Fig. 7) affinity(T=c(25, 110, 115, 215), P=c(11, 335, 500, 1450))