Skip to content

Calliope.oxygen fugacity

oxygen_fugacity

OxygenFugacity(model=DEFAULT_FO2_MODEL)

log10 oxygen fugacity as a function of temperature

Source code in src/calliope/oxygen_fugacity.py 
21
22
def __init__(self, model=DEFAULT_FO2_MODEL):
    self.callmodel = getattr(self, model)

fischer(T)

Fischer et al. (2011) IW (FeO equation of state, EPSL 304, 496).

The coefficients are the p = 1 bar isoline of their Figure 6, obtained by integrating their Equation 2.

Source code in src/calliope/oxygen_fugacity.py 
33
34
35
36
37
38
39
def fischer(self, T):
    """Fischer et al. (2011) IW (FeO equation of state, EPSL 304, 496).

    The coefficients are the p = 1 bar isoline of their Figure 6,
    obtained by integrating their Equation 2.
    """
    return 6.94059 - 28.1808 * 1e3 / T

oneill(T)

O'Neill and Eggins (2002) IW

Source code in src/calliope/oxygen_fugacity.py 
41
42
43
44
45
def oneill(self, T):
    """O'Neill and Eggins (2002) IW"""
    # 8.31441 reproduces O'Neill & Eggins (2002) Eq. 11 verbatim;
    # do not replace with constants.R_gas (8.31446...).
    return 2 * (-244118 + 115.559 * T - 8.474 * T * np.log(T)) / (np.log(10) * 8.31441 * T)