# 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': .100,
'n': 5,
'Debye_0': 1070.,
'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': .100,
'n': 5,
'Debye_0': 1070.,
'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': .102,
'n': 5,
'Debye_0': 1054.,
'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': .0403,
'n': 2,
'Debye_0': 773., # 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': .047,
'n': 2,
'Debye_0': 706.,
'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': .047,
'n': 2,
'Debye_0': 706.,
'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': .0469,
'n': 2,
'Debye_0': 706.,
'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': .0469,
'n': 2,
'Debye_0': 706.,
'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