Proton behaviour, structure and elasticity of serpentine at high-pressure

Abstract Serpentines occur in the subduction zone settings, both along the slab and within the mantle wedge, they are candidates for transporting water in to the deep earth. Their presence is manifested by serpentine mud volcanoes, high electrical conductivities, magnetic and seismic anomalies. Using theoretical methods, we predict a pressure induced structural transformations in serpentine. The transformations are related to the behavior of the silicate framework and misfit between octahedral and tetrahedral layers. As the structure is compressed, the octahedral layer and tetrahedral layers are compressed at different rates. At 7 GPa, the misfit between the layers vanishes. This causes non-linear pressure dependence of tetrahedral rotational angle. This is also manifested by the onset of anomalous pressure dependence of the elastic constants c 11 , c 33 , c 12 , c 13 . Beyond 7 GPa, the misfit between the layers grows again reaching extremum at 22 GPa. This is also manifested by discontinuity in average Si–O bond length, volume of tetrahedron and re-orientation of hydroxyl vector. The symmetry of the crystal-structure however, remains unaffected. Evidence of pressure-induced hydrogen bonding is absent in serpentine, as evident from reduction of O–H bond length upon compression. Results of compression for the low-pressure regime ( P K 0  = 79 GPa, K 0 ′ = 12 and K 0 ″ = − 2, where K is the bulk modulus, prime indicates pressure derivatives, and O refers to zero pressure. Our best estimates of K 0 , K 0 ′ and the Gruneisen parameter, γ at 300 K and zero pressure based on our results are: 61 GPa, 17, and 0.77, respectively. At low pressures, serpentine structure is anisotropic with c 11  ~ 2.4 ×  c 33 . The pressure derivative of elastic constants ( ∂c ij / ∂P ) are such, that around 22 GPa c 11 ~ c 33 . An elastic instability ( c 66  50 GPa) is also noted. The elastic constant tensor reveals large acoustic anisotropy (41% in V P ) and seismic wave velocities that are significantly higher than those inferred from experiments on serpentinites.

[1]  R. Hemley,et al.  High-pressure infrared spectra of talc and lawsonite , 2006 .

[2]  K. Katahara Clay Mineral Elastic Properties , 1996 .

[3]  D. Miller,et al.  Mantle wedge water contents estimated from seismic velocities in partially serpentinized peridotites , 2003 .

[4]  R. Jeanloz,et al.  Deep-Focus Earthquakes and Recycling of Water into the Earth's Mantle , 1991, Science.

[5]  M. Mottl,et al.  Mariana blueschist mud volcanism: Implications for conditions within the subduction zone , 1999 .

[6]  A. Yamagishi,et al.  First-principles studies on the elastic constants of a 1:1 layered kaolinite mineral , 2005 .

[7]  L. A. Reinen Seismic and aseismic slip indicators in serpentinite gouge , 2000 .

[8]  D. Chadi,et al.  Special points for Brillouin-zone integrations , 1977 .

[9]  N. Christensen,et al.  Serpentinites, Peridotites, and Seismology , 2004 .

[10]  H. Mao,et al.  Sound Velocities of Hot Dense Iron: Birch's Law Revisited , 2005, Science.

[11]  R. Wentzcovitch,et al.  Invariant molecular-dynamics approach to structural phase transitions. , 1991, Physical review. B, Condensed matter.

[12]  T. Duffy,et al.  Single-crystal elasticity of wadsleyites, β -Mg 2 SiO 4 , containing 0.37-1.66 wt.% H 2 O (DOI:10.1016/j.epsl.2007.10.045) , 2008 .

[13]  L. Stixrude,et al.  High-pressure proton disorder in brucite , 2004 .

[14]  G. Kresse,et al.  Ab initio molecular dynamics for liquid metals. , 1993 .

[15]  S. Fukura,et al.  Pressure-induced change in the compressional behavior of the O-H bond in chrysotile: A Raman high-pressure study up to 4.5 GPa , 2007 .

[16]  J. Parise,et al.  Elastic properties of hydrous ringwoodite (γ-phase) in Mg2SiO4 , 1998 .

[17]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[18]  R. Hyndman,et al.  An inverted continental Moho and serpentinization of the forearc mantle , 2002, Nature.

[19]  N. Christensen Elasticity of ultrabasic rocks , 1966 .

[20]  T. Duffy,et al.  Elasticity of enstatite and its relationship to crystal structure , 1986 .

[21]  B. Reynard,et al.  Equation of state of antigorite, stability field of serpentines, and seismicity in subduction zones , 2006 .

[22]  L. Benco,et al.  HARTREE-FOCK STUDY OF PRESSURE-INDUCED STRENGTHENING OF HYDROGEN BONDING IN LIZARDITE-1T , 1998 .

[23]  T. Duffy,et al.  Single-crystal elasticity of brucite, Mg(OH)2, to 15 GPa by Brillouin scattering , 2006 .

[24]  G. Kresse,et al.  Optimized norm-conserving pseudopotentials , 1992 .

[25]  H. Mao,et al.  Pressure-induced amorphization of crystalline silica , 1988, Nature.

[26]  D. Wallace,et al.  Thermodynamics of Crystals , 1972 .

[27]  J. Tsuchiya,et al.  First-principles study of hydrogen bond symmetrization of phase D under high pressure , 2005 .

[28]  G. Kresse,et al.  Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set , 1996 .

[29]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[30]  R. Pellenq,et al.  Atomistic calculations of structural and elastic properties of serpentine minerals: the case of lizardite , 2006 .

[31]  Kazushige Obara,et al.  Nonvolcanic Deep Tremor Associated with Subduction in Southwest Japan , 2002, Science.

[32]  F. Birch The velocity of compressional waves in rocks to 10 kilobars: 1. , 1960 .

[33]  T. Yagi,et al.  Amorphization of Serpentine at High Pressure and High Temperature , 1996, Science.

[34]  G. V. Gibbs,et al.  Quadratic Elongation: A Quantitative Measure of Distortion in Coordination Polyhedra , 1971, Science.

[35]  M. Roderick,et al.  The cause of decreased pan evaporation over the past 50 years. , 2002, Science.

[36]  Volker Heine,et al.  The Pseudopotential Concept , 1970 .

[37]  R. Hemley,et al.  Theory of minerals at high pressure , 1998 .

[38]  J. Rodgers,et al.  Ab initio elastic properties of talc from 0 to 12 GPa: Interpretation of seismic velocities at mantle pressures and prediction of auxetic behaviour at low pressure , 2008 .

[39]  D. Lauretta,et al.  A Nebular Origin for Chondritic Fine-Grained Phyllosilicates , 2003, Science.

[40]  M. O'keeffe,et al.  On Si—O—Si configurations in silicates , 1978 .

[41]  L. Stixrude,et al.  Hydrogen incorporation in stishovite at high pressure and symmetric hydrogen bonding in δ-AlOOH , 2004 .

[42]  J. Zussman,et al.  The crystal structure of the serpentine mineral, lizardite Mg3Si2O5(OH)4 , 1965 .

[43]  F. Mauri,et al.  First-principles calculation of the infrared spectrum of lizardite , 2002 .

[44]  L. Stixrude Talc under tension and compression: Spinodal instability, elasticity, and structure , 2002 .

[45]  Price,et al.  Ab initio molecular dynamics with variable cell shape: Application to MgSiO3. , 1993, Physical review letters.

[46]  L. Dubrovinsky,et al.  The high-pressure behaviour of the 10 Å phase: A spectroscopic and diffractometric study up to 42 GPa , 2006 .

[47]  P. Meredith,et al.  Simulation of Subduction Zone Seismicity by Dehydration of Serpentine , 2002, Science.

[48]  Renata M. Wentzcovitch,et al.  High‐pressure elastic properties of major materials of Earth's mantle from first principles , 2001 .

[49]  Lars Stixrude,et al.  Thermodynamics of mantle minerals – I. Physical properties , 2005 .

[50]  P. Comodi,et al.  High-pressure structural study of muscovite , 1995 .

[51]  N. Christensen Ophiolites, seismic velocities and oceanic crustal structure , 1978 .

[52]  J. Holloway,et al.  Experimental determination of elastic properties of talc to 800°C, 0.5 GPa; calculations of the effect on hydrated peridotite, and implications for cold subduction zones , 2000 .

[53]  C. Viti,et al.  Hydrogen positions and thermal expansion in lizardite-1T from Elba: A low-temperature study using Rietveld refinement of neutron diffraction data , 1996 .

[54]  B. Grobéty Polytypes and higher-order structures of antigorite: A TEM study , 2003 .

[55]  E. Knittle,et al.  The equation of state, amorphization, and high-pressure phase diagram of muscovite , 1994 .

[56]  T. Ahrens,et al.  Shock wave equation of state of serpentine to 150 GPa: Implications for the occurrence of water in the Earth's lower mantle , 1991 .

[57]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[58]  A. Kleppe,et al.  The effect of pressure upon hydrogen bonding in chlorite: A Raman spectroscopic study of clinochlore to 26.5 GPa , 2003 .

[59]  D. Vanderbilt,et al.  Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. , 1990, Physical review. B, Condensed matter.

[60]  Holtz,et al.  Effect of pressure on the Raman vibrational modes of layered aluminosilicate compounds. , 1993, Physical review. B, Condensed matter.

[61]  M. Lazzeri,et al.  Elasticity of serpentines and extensive serpentinization in subduction zones , 2007 .

[62]  R. A. Robie,et al.  Thermodynamic properties of minerals and related substances at 298.15 K and 1 bar (10[5] pascals) pressure and at higher temperatures , 1995 .

[63]  P. Dera,et al.  Novel pressure‐induced phase transformations in hydrous layered materials , 2002 .

[64]  F. Birch,et al.  Finite strain isotherm and velocities for single‐crystal and polycrystalline NaCl at high pressures and 300°K , 1978 .

[65]  M. Drábek,et al.  ‘‘Cs-tetra-ferri-annite:’’ High-pressure and high-temperature behavior of a potential nuclear waste disposal phase , 1999 .

[66]  Robert M. Hazen,et al.  Bulk modulus—volume relationship for cation-anion polyhedra , 1979 .

[67]  J. Nye Physical Properties of Crystals: Their Representation by Tensors and Matrices , 1957 .

[68]  Simon M. Peacock,et al.  Serpentinization of the forearc mantle , 2003 .

[69]  B. Reynard,et al.  P–V Equations of State and the relative stabilities of serpentine varieties , 2006 .

[70]  L. Stixrude,et al.  The 10 Å phase at high pressure by first principles calculations and implications for the petrology of subduction zones , 2007 .

[71]  Robert M. Hazen,et al.  The crystal structures and compressibilities of layer minerals at high pressure; II, Phlogopite and chlorite , 1978 .

[72]  Lars Stixrude,et al.  First-principles study of illite–smectite and implications for clay mineral systems , 2002, Nature.

[73]  J. Rodgers,et al.  Predicted elastic properties of the hydrous D phase at mantle pressures: Implications for the anisotropy of subducted slabs near 670-km discontinuity and in the lower mantle , 2007 .

[74]  B. Reynard,et al.  High-Pressure Creep of Serpentine, Interseismic Deformation, and Initiation of Subduction , 2007, Science.

[75]  S. Clark,et al.  Equation of state measurements of chlorite, pyrophyllite, and talc , 2002 .

[76]  R. Jeanloz,et al.  Static compression of Ca(OH)2 at room temperature: Observations of amorphization and equation of sta , 1990 .

[77]  M. Mellini,et al.  Effects of pressure on the structure of lizardite-1T , 1989 .

[78]  N. H. March,et al.  Theory of the inhomogeneous electron gas , 1983 .

[79]  Y. Fei,et al.  Static compression of the hydrous magnesium silicate phase D to 30 GPa at room temperature , 1999 .

[80]  B. Reynard,et al.  High-pressure behaviour of serpentine minerals: a Raman spectroscopic study , 2004 .