Skip to content

Calliope.oxygen fugacity

oxygen_fugacity

OxygenFugacity(model='oneill')

log10 oxygen fugacity as a function of temperature

Source code in src/calliope/oxygen_fugacity.py 
14
15
def __init__(self, model='oneill'):
    self.callmodel = getattr(self, model)

fischer(T)

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

Source code in src/calliope/oxygen_fugacity.py 
26
27
28
def fischer(self, T):
    """Fischer et al. (2011) IW (FeO equation of state, EPSL 304, 496)"""
    return 6.94059 - 28.1808 * 1e3 / T

oneill(T)

O'Neill and Eggins (2002) IW

Source code in src/calliope/oxygen_fugacity.py 
30
31
32
33
34
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)