Visualization-based analysis of structural and dynamical properties of simulated hydrous silicate melt
Tóm tắt
We have explored first-principles molecular dynamics simulation data for hydrous MgSiO3 liquid (with 10 wt% water) to gain insight into its structural and dynamical behavior as a function of pressure (0–150 GPa) and temperature (2,000–6,000 K). By visualizing/analyzing a number of parameters associated with short- and mid-range orders, we have shown that the melt structure changes substantially on compression. The speciation of the water component at low pressures is dominated by the isolated structures (with over 90% hydrogen participated) consisting of hydroxyls, water molecules, O–H–O bridging and four-atom (O–H–O–H and H–O–H–O) groups, where every oxygen atom may be a part of polyhedron or free (i.e., bound to only magnesium atom). Hydroxyls favor polyhedral sites over magnesium sites whereas molecular water is almost entirely bound to magnesium sites, and also interpolyhedral bridging (Si–O–H–O–Si) dominates other types of bridging. Water content is shown to enhance and suppress, respectively, the proportions of hydroxyls and molecular water. As compression increases, these isolated structures increasingly combine with each other to form extended structures involving a total of five or more O and H atoms and also containing threefold coordination species, which together consume over 80% hydrogen at the highest compression studied. Our results show that water lowers the mean coordination numbers of different types including all cation–anion environments. The hydrous melt tends to be more tetrahedrally coordinated but with the Si–Si network being more disrupted compared to the anhydrous melt. Protons increase the content of non-bridging oxygen and decrease the contents of bridging oxygen as well as oxygen triclusters (present at pressures above 10 GPa). The calculated self-diffusion coefficients of all atomic species are enhanced in the presence of water compared to those of the anhydrous melt. This is consistent with the prediction that water depolymerizes the melt structure at all pressures. Our analysis also suggests that proton diffusion involves two processes—the transfer of H atoms (requiring the rupture and formation of O–H bonds) and the motion of hydroxyls as hydrogen carriers (requiring the rupture and formation of Si–O and/or Mg–O bonds). Both the processes are operative at low compression whereas only the first process is operative at high compression.
Tài liệu tham khảo
Angell CA, Cheeseman PA, Tamaddon S (1982) Pressure enhancement of ion mobilities in liquid silicates from computer-simulation studies to 800-kilobars. Science 218:885–887
Bercovici D, Karato S (2003) Whole mantle convection and the transition-zone water filter. Nature 425:39–44
Bhattarai D, Karki BB (2009) Atomistic visualization: space–time multiresolution integration of data analysis and rendering. J Mol Graph Model 27:951–968
Ceperley DM, Alder BJ (1980) Ground state of the electron gas by a stochastic method. Phys Rev Lett 45:566–569
Closmann C, Williams Q (1995) In-situ spectroscopic investigation of high-pressure hydrated (Mg, Fe)SiO3 glasses—OH vibrations as a probe of glass structure. Am Mineral 80:201–212
de Koker N, Stixrude L, Karki BB (2008) Thermodynamics, structure, dynamics, and melting of Mg2SiO4 liquid at high pressure. Geochim Cosmochim Acta 72:1427–1441
Farber DL, Williams Q (1996) An in situ Raman spectroscopic study of Na2Si2O5 at high pressures and temperatures: structure of compressed liquids and glasses. Am Mineral 81:273–283
Funamori N, Yamamoto S, Yagi T, Kikegawa T (2004) Exploratory studies of silicate melt structure at high pressures and temperatures by in situ X-ray diffraction. J Geophys Res 109:B03203
Giordano D, Romano DB, Dingwell DB, Poe B, Behrens H (2004) The combined effects of water and fluorine on the viscosity of silicic magmas. Geochim Cosmochim Acta 68:5159
Goldman N, Fried LE, Kuo IFW, Mundy CJ (2005) Bonding in the superionic phase of water. Phys Rev Lett 94:217801
Goncharov AF, Struzhkin VV, Mao HK, Hemley RJ (1999) Raman spectroscopy of dense H2O and the transition to symmetric hydrogen bonds. Phys Rev Lett 83:1998–2001
Inoue T (1994) Effect of water on melting phase relations and melt composition in the system Mg2SiO4–MgSiO3–H2O up to 15 GPa. Phys Earth Planet Inter 85:237–263
Karki BB, Bhattarai D, Stixrude L (2007) First principles simulations of liquid silica: structural and dynamical behavior at high pressure. Phys Rev B 76:104205
Kresse G, Furthmuller J (1996) Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput Mater Sci 6:15–50
Kresse G, Hafner J (1994) Norm-conserving and ultrasoft pseudopotentials for first-row and transition-elements. J Phys Condens Matter 6:8245–8257
Kubicki JD, Lasaga AC (1991) Molecular dynamics simulations of pressure and temperature effects on MgSiO3 and Mg2SiO4 melts and glasses. Phys Chem Miner 17:661–673
Lacks DJ, Rear DB, Van Orman JA (2007) Molecular dynamics investigation of viscosity, chemical diffusivities and partial molar volumes of liquids along the MgO–SiO2 joins as functions of pressure. Geochim Cosmochim Acta 71:1312–1323
Lange RA (1994) The effect of H2O, CO2 and F on the density and viscosity of silicate melts. Rev Mineral 30:331–369
Lange RA, Carmichael ISE (1987) Densities of Na2O–K2O–MgO–FeO–Fe2O3–Al2O3–TiO2–SiO2 liquids—new measurements and derived partial molar properties. Geochim Cosmochim Acta 51:2931–2946
Lee SK, Li JF, Cai YQ, Hiraoka N, Eng PJ, Okuchi T, Mao HK, Meng Y, Hu MY, Chow P, Shu J, Li B, Fukui H, Lee BH, Kim HN, Yoo CS (2008) X-ray Raman scattering study of MgSiO3 glass at high pressure: implication for triclustered MgSiO3 melt in Earth’s mantle. Proc Natl Acad Sci USA 105:7925–7929
Lesher CE, Hervig RL, Tinker D (1996) Self diffusion of network formers (silicon and oxygen) in naturally occurring basaltic liquid. Geochim Cosmochim Acta 60:405–413
Mao WL, Mao HK, Meng Y, Eng PJ, Hu MY, Chow P, Cai YQ, Shu J, Hemley RJ (2006) X-ray-induced dissociation of H2O and formation of O2–H2 alloy at high pressure. Science 314:636–638
Martin GB, Spera FJ, Ghiorso MS, Nevins D (2009) Structure, thermodynamic and transport properties of molten Mg2SiO4: molecular dynamics simulations and model EOS. Am Mineral 94:693–703
Matsukage K, Zhicheng J, Karato S (2005) Density of hydrous silicate melt at the conditions of Earth’s deep mantle. Nature 438:488–491
Mermin ND (1965) Thermal properties of inhomogeneous electron gas. Phys Rev 137:A1441–A1443
Mookherjee M, Stixrude L, Karki BB (2008) Hydrous silicate melt at high pressure. Nature 452:983–986
Ochs FA, Lange RA (1999) The density of hydrous magmatic liquids. Science 283:1314–1317
Pohlmann M, Benoit M, Kob W (2004) First-principles molecular dynamics simulations of a hydrous silica melt: structural properties and hydrogen diffusion mechanism. Phys Rev B 70:184209
Revenaugh J, Sipkin SA (1994) Seismic evidence for silicate melt atop the 410-km mantle discontinuity. Nature 369:474–476
Rustad JR, Yuen DA, Spera FJ (1990) Molecular-dynamics of liquid SiO2 under high-pressure. Phys Rev A 42:2081–2089
Sakamaki T, Suzuki A, Ohtani E (2006) Stability of hydrous melt at the base of the Earth’s upper mantle. Nature 439:192–194
Shen A, Keppler H (1995) Infrared spectroscopy of hydrous silicate melts to 1000°C and 10 kbar: direct observation of H2O speciation in a diamond-anvil cell. Am Mineral 80:1335–1338
Shen A, Keppler H (1997) Direct observation of complete miscibility the albite-H2O system. Nature 385:710–712
Stixrude L, Karki BB (2005) Structure and freezing of MgSiO3 liquid in Earth’s lower mantle. Science 310:297–299
Stixrude L, Lithgow-Bertelloni C (2005) Thermodynamics of mantle minerals. I. Physical properties. Geophys J Int 162:610–632
Stolper EM (1982) The speciation of water in silicate melts. Geochim Cosmochim Acta 46:2609–2620
Tinker D, Lesher CE, Hutcheon (2003) Self-diffusion of Si and O in diopside–anorthite melt at high pressure. Geochim Cosmochim Acta 67:133–142
Tsuneyuki S, Matsui Y (1995) Molecular dynamics study of pressure enhancement of ion mobilities in liquid silica. Phys Rev Lett 74:3197–3200
Tyburczy JA, Duffy TS, Ahrens TJ, Lange MA (1991) Shock wave equation of state of serpentine to 150 GPa: implications for the occurrence of water in the Earth’s mantle. J Geophys Res 96:18011–18027
Wan TK, Duffy TS, Scandolo S, Car R (2007) First principles study of density, viscosity and diffusion coefficients of liquid MgSiO3 at conditions of the Earth’s deep mantle. J Geophys Res 112:03208. doi:10.1029/2005JB004135
Wang Y, Perdew J (1991) Correlation hole of the spin-polarized electron gas with exact small-wave-vector and high-density scaling. Phys Rev B 44:13298–13307
Waseda Y, Toguri JM (1990) Structure of silicate melts determined by X-ray diffraction. In: Marumo F (ed) Dynamics process of material transport and transformation in earth’s interior. Terra Scientific Publishing Company, Tokyo, pp 37–51
Wasserman EA, Yuen DA, Rustad JR (1993) Molecular dynamics study of the transport properties of perovskite melts under high temperature and pressure conditions. Earth Planet Sci Lett 114:373–384
Watson EB (1981) Diffusion in magmas at depth in the Earth: the effects of pressure and dissolved H2O. Earth Planet Sci Lett 52:291–301
Wentzcovitch RM, Martins JL, Allen PB (1992) Energy versus free-energy conservation in first-principles molecular dynamics. Phys Rev B 45:11372–11374
Whittington A, Richet P, Holtz F (2000) Water and the viscosity of depolymerized aluminosilicate melts. Geochim Cosmochim Acta 64:3725
Williams Q, Jeanloz R (1988) Spectroscopic evidence for pressure-induced coordination changes in silicate-glasses and melts. Science 239:902–905
Wolf GH, Durben DJ, McMillan PF (1990) High-pressure Raman spectroscopic study of sodium tetrasilicate (Na2Si4O9) galls. J Chem Phys 93:2280–2288
Xue X, Stebbins JF, Kanzaki M, McMillan PF, Poe B (1991) Pressure induced silicon coordination and tetrahedral structural changes in alkali silicate melts up to 12 GPa: NMR, Raman, and infrared spectroscopy. Am Mineral 76:8–26
Yamada Y, Inoue T, Irifune T (2004) Melting of enstatite from 13 to 18 GPa under hydrous conditions. Phys Earth Planet Inter 147:45–56
Yamada A, Inoue T, Urakawa A, Funakoshi K, Funamori N, Kikegawa T, Ohfuji H, Irifune T (2007) In situ X-ray experiment on the structure of hydrous Mg–Silicate melt under high pressure and high temperature. Geophys Res Lett 34:L10303
Zhang Y (1999) H2O in rhyolitic glasses and melts: measurement, speciation, solubility and diffusion. Rev Geophys 37:493–516
Zhang Y, Stolper EM (1991) Water diffusion in a basaltic melt. Nature 351:306–309
Zhang G, Guo G, Refson K, Zhao Y (2004) Finite-size effects at both high and low temperatures in molecular dynamics calculations of the self-diffusion coefficient and viscosity of liquid silica. J Phys Condens Matter 16:9127–9135
Zhang Y, Xu Z, Zhu M, Wang H (2007) Silicate melt properties and volcanic eruptions. Rev Geophys 45:RG4004. doi:10.1029/2006RG000216