Show table of references in a web browser or get individual references for species. Check self consistency of individual entries in database.
1 2 3 4 5 6 7
character, numeric, or list; bibliographic reference key(s)
logical, keep duplicated references?
dataframe, equations-of-state parameters in the format of
character, physical state of species
character, property of interest (Cp or V)
logical, return the difference between calculated and tabulated values?
dataframe, containing G, H and S, in the format of
character, path to a file
character, name of compound(s) in group additivity calculation
thermo.refs with default arguments uses
browseURL to display the sources of thermodynamic data in
thermo()$refs, with the URLs in that table showing as hyperlinks in the browser.
key is character, the citation information for those reference
keys (including URLs) are returned.
key is numeric, the values refer to the species in those rows of
thermo()$OBIGT, and the citation information for each listed reference (
thermo()$OBIGT$ref2) is returned.
key is a list, it is interpreted as the result of a call to
subcrt, and the citation information for each species involved in the calculation is returned.
Only unique references are returned, unless
keep.duplicates is TRUE.
In that case, a single reference for each species is returned, ignoring anything in
checkEOS compares heat capacity and volume calculated from equation-of-state parameters with reference (tabulated) values at 25 \degC and 1 bar and prints a message and returns the calculated value if tolerance is exceeded.
The Helgeson-Kirkham-Flowers equations of state parameters are in
eos, which is a data frame with columns (and column names) in the same format as
property can be one of Cp or
The code only distinguishes between
states of aq and all others.
The default tolerances, given in
thermo()$opt$V.tol, are 1 cal/K.mol for Cp and 1 cm3/mol for V.
ret.diff is TRUE, the differences are returned irrespective of their values, and no messages are printed.
checkGHS compares G (standard molal Gibbs energy of formation from the elements) calculated from H (standard molal enthalpy of formation) and S (standard molal entropy) with reference (tabulated) values of G at 25 \degC and 1 bar.
A message is printed and the calculated difference is returned if it exceeds the value given in
thermo()$opt$G.tol, which has a default value of 100 cal/mol.
The calculation requires that G, H and S, and the chemical formula of the species all be present.
check.OBIGT is a function to check self-consistency of each entry in the thermodynamic database, using
The output is a table listing only species that exceed at least one of the tolerance limits, giving the species index (rownumber in 'thermo()$OBIGT'), species name and state, and DCp, DV and DG, for the calculated differences (only those above the tolerances are given).
Values of DCp and DG are given in the units present in the data files.
This function is used to generate the file found at
dumpdata returns all of the available data, from both the default and optional data files, or writes it to a file if
file is not NULL.
The format is the same as
thermo$OBIGT, except for a single prepended column named source, giving the source of the data (OBIGT refers to the default database, and DEW, SLOP98, and SUPCRT92 are the optional data files).
RH2OBIGT implements a group additivity algorithm for standard molal thermodynamic properties and equations of state parameters of crystalline and liquid organic molecules from Richard and Helgeson, 1998. The names of the
compounds and their physical
state are searched for in the indicated
file, that also contains chemical formulas and group stoichiometries; the names of the groups are stored in the column names of this file, and must be present in
thermo$OBIGT. The default
extdata/thermo/RH98_Table15.csv) includes data taken from Table 15 of Richard and Helgeson, 1998 for high molecular weight compounds in crystalline and liquid states. An error is produced if any of the
state combinations is not found in the
file, if any of the group names for a given
state combination is not found in
thermo()$OBIGT, or if the chemical formula calculated from group additivity (with the aid of
as.chemical.formula) is not identical to that listed in the
The values returned (
mod.OBIGT are the rownumbers of the affected species.
Richard, L. and Helgeson, H. C. (1998) Calculation of the thermodynamic properties at elevated temperatures and pressures of saturated and aromatic high molecular weight solid and liquid hydrocarbons in kerogen, bitumen, petroleum, and other organic matter of biogeochemical interest. Geochim. Cosmochim. Acta 62, 3591–3636. https://doi.org/10.1016/S0016-7037(97)00345-1
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
# citation information for Helgeson et al., 1998 thermo.refs("HOK+98") # two references for alanine thermo.refs(info("alanine")) # three references for species in the reaction s <- subcrt(c("O2","O2"),c("gas","aq"),c(-1,1)) thermo.refs(s) ## Not run: ## marked dontrun because it opens a browser # show the contents of thermo()$refs thermo.refs() ## End(Not run) ## calculate thermodynamic properties of organic compounds ## using group additivity, after Richard and Helgeson, 1998 RH2OBIGT()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.