# This file is part of BurnMan - a thermoelastic and thermodynamic toolkit for the Earth and Planetary Sciences
# Copyright (C) 2012 - 2017 by the BurnMan team, released under the GNU
# GPL v2 or later.
"""
Murakami_etal_2012
^^^^^^^^^^^^^^^^^^
Minerals from Murakami et al. (2012) supplementary table 5 and references therein, V_0 from
Stixrude & Lithgow-Bertolloni 2005. Some information from personal communication with Murakami.
"""
from __future__ import absolute_import
from ..classes import mineral_helpers as helpers
from ..classes.mineral import Mineral
[docs]
class mg_perovskite(Mineral):
def __init__(self):
self.params = {
"equation_of_state": "slb2",
"V_0": 24.45e-6, # S & L-B 2005
"K_0": 281e9,
"Kprime_0": 4.1,
"G_0": 173e9,
"Gprime_0": 1.56,
"molar_mass": 0.100,
"n": 5,
"Debye_0": 1070.0,
"grueneisen_0": 1.48,
"q_0": 1.4,
"eta_s_0": 2.4,
}
Mineral.__init__(self)
[docs]
class mg_perovskite_3rdorder(Mineral):
def __init__(self):
self.params = {
"equation_of_state": "slb3",
"V_0": 24.45e-6, # S & L-B 2005
"K_0": 281e9,
"Kprime_0": 4.1,
"G_0": 171.42e9,
"Gprime_0": 1.83,
"molar_mass": 0.100,
"n": 5,
"Debye_0": 1070.0,
"grueneisen_0": 1.48,
"q_0": 1.4,
"eta_s_0": 2.4,
}
Mineral.__init__(self)
[docs]
class fe_perovskite(Mineral):
def __init__(self):
self.params = {
"equation_of_state": "slb2",
"V_0": 24.607e-6,
"K_0": 251.9e9,
"Kprime_0": 4.01,
"G_0": 164.7e9,
"Gprime_0": 1.58,
"molar_mass": 0.102,
"n": 5,
"Debye_0": 1054.0,
"grueneisen_0": 1.48,
"q_0": 1.4,
"eta_s_0": 2.4,
}
Mineral.__init__(self)
[docs]
class mg_periclase(Mineral):
def __init__(self):
self.params = {
"equation_of_state": "slb2",
"V_0": 11.24e-6, # S & L-B 2005
"K_0": 161e9,
"Kprime_0": 3.9,
"G_0": 131e9,
"Gprime_0": 1.92,
"molar_mass": 0.0403,
"n": 2,
"Debye_0": 773.0, # S& L-B 2005
"grueneisen_0": 1.5,
"q_0": 1.5, # S&L-B 2005
"eta_s_0": 3.0,
}
Mineral.__init__(self)
[docs]
class fe_periclase(helpers.HelperSpinTransition):
def __init__(self):
helpers.HelperSpinTransition.__init__(
self, 63.0e9, fe_periclase_LS(), fe_periclase_HS()
)
[docs]
class fe_periclase_3rd(helpers.HelperSpinTransition):
def __init__(self):
helpers.HelperSpinTransition.__init__(
self, 63.0e9, fe_periclase_LS(), fe_periclase_HS()
)
[docs]
class fe_periclase_HS(
Mineral
): # From Murakami's emails, see Cayman for details, represents Mg# = .79
def __init__(self):
self.params = {
"equation_of_state": "slb2",
"V_0": 11.412e-6,
"K_0": 159.1e9,
"Kprime_0": 4.11,
"G_0": 105.43e9,
"Gprime_0": 1.773,
"molar_mass": 0.047,
"n": 2,
"Debye_0": 706.0,
"grueneisen_0": 1.45,
"q_0": 1.5,
"eta_s_0": 2.54,
}
Mineral.__init__(self)
[docs]
class fe_periclase_LS(
Mineral
): # From Murakami's emails, see Cayman for details, represents Mg# = .79
def __init__(self):
self.params = {
"equation_of_state": "slb2",
"V_0": 11.171e-6,
"K_0": 170.0e9,
"Kprime_0": 4.00,
"G_0": 116.34e9,
"Gprime_0": 1.668,
"molar_mass": 0.047,
"n": 2,
"Debye_0": 706.0,
"grueneisen_0": 1.45,
"q_0": 1.5,
"eta_s_0": 2.54,
}
Mineral.__init__(self)
[docs]
class fe_periclase_HS_3rd(Mineral):
def __init__(self):
self.params = {
"equation_of_state": "slb3",
"V_0": 11.412e-6,
"K_0": 159.1e9,
"Kprime_0": 4.11,
"G_0": 129.35e9,
"Gprime_0": 1.993,
"molar_mass": 0.0469,
"n": 2,
"Debye_0": 706.0,
"grueneisen_0": 1.45,
"q_0": 1.5,
"eta_s_0": 2.54,
}
Mineral.__init__(self)
[docs]
class fe_periclase_LS_3rd(Mineral):
def __init__(self):
self.params = {
"equation_of_state": "slb3",
"V_0": 11.171e-6,
"K_0": 170.0e9,
"Kprime_0": 4.0,
"G_0": 151.67e9,
"Gprime_0": 1.754,
"molar_mass": 0.0469,
"n": 2,
"Debye_0": 706.0,
"grueneisen_0": 1.45,
"q_0": 1.5,
"eta_s_0": 2.54,
}
Mineral.__init__(self)
mg_bridgmanite = mg_perovskite
fe_bridgmanite = fe_perovskite
mg_bridgmanite_3rdorder = mg_perovskite_3rdorder