Abstract
Speciation and transport properties of supercritical fluids is critical for understanding their behavior in the Earth’s interior. Here, we report a systematic first principles molecular dynamics simulation study of the structure, speciation, self-difusivity (D), and viscosity (η) of SiO2 melt, NaAlSi3O8 melt, SiO2-H2O and NaAlSi3O8-H2O fluids at 2000–3500 K with 0–70 wt% H2O. Our calculations show that as the water content increases, the proportion of Q0 species (Qn species, where n is the number of bridging oxygens in an individual Si/Al-O polyhedra) increases while Q4 decreases. The proportions of Q1, Q2, and Q3 species first increase and then decrease with increasing water content. The diffusivity sequence for the supercritical SiO2-H2O fluids is DH >DO >DSi, and for the supercritical NaAlSi3O8-H2O fluids, on the whole, is DNa ≈ DH >DO >DAl ≈ DSi. The viscosities of the two systems decrease drastically at the beginning of the increase in water content, and then decrease slowly. We demonstrate that the exponential decrease in the viscosity of polymerized silicate melt with increasing water content is due to a sharp decrease in the proportion of Q4 species and increase in Si-O-H. The typical structural feature of supercritical fluid is that it contains a large amount of easy-to-flow partially polymerized or depolymerized protonated silicate units, which leads to a low viscosity while being enriched in silicate. This feature provides supercritical fluids the potential to transport elements that are hard to migrate in aqueous fluids or hydrous silicate melts, such as high field strength elements.
References cited
Allen, M.P. and Tildesley, D.J. (1989) Computer Simulation of Liquids, 408 p. Oxford University Press.Search in Google Scholar
Anderson, G.M. and Burnham, C.W. (1965) The solubility of quartz in supercritical water. American Journal of Science, 263(6), 494–511.Search in Google Scholar
Audétat, A. and Keppler, H. (2004) Viscosity of fluids in subduction zones. Science, 303, 513–516, https://doi.org/10.1126/science.1092282Search in Google Scholar
Audétat, A. and Keppler, H. (2005) Solubility of rutile in subduction zone fluids, as determined by experiments in the hydrothermal diamond anvil cell. Earth and Planetary Science Letters, 232, 393–402, https://doi.org/10.1016/j.epsl.2005.01.028.Search in Google Scholar
Bajgain, S.K. and Mookherjee, M. (2020) Structure and properties of albite melt at high pressures. ACS Earth & Space Chemistry, 4, 1–13, https://doi.org/10.1021/acsearthspacechem.9b00187.Search in Google Scholar
Bajgain, S.K., Peng, Y., Mookherjee, M., Jing, Z.C., and Solomon, M. (2019) Properties of hydrous aluminosilicate melts at high pressures. ACS Earth & Space Chemistry, 3, 390–402, https://doi.org/10.1021/acsearthspacechem.8b00157Search in Google Scholar
Bockris, J.O., Mackenzie, J.D., and Kitchener, J.A. (1955) Viscous flow in silica and binary liquid silicates. Transactions of the Faraday Society, 51, 1734–1748, https://doi.org/10.1039/tf9555101734.Search in Google Scholar
Bureau, H. and Keppler, H. (1999) Complete miscibility between silicate melts and hydrous fluids in the upper mantle: experimental evidence and geochemical implications. Earth and Planetary Science Letters, 165(2), 187–196.Search in Google Scholar
Caldeweyher, E., Ehlert, S., Hansen, A., Neugebauer, H., Spicher, S., Bannwarth, C., and Grimme, S. (2019) A generally applicable atomic-charge dependent London dispersion correction. The Journal of Chemical Physics, 150, 154122, https://doi.org/10.1063/1.5090222Search in Google Scholar
Chen, W., Xiong, X.L., Wang, J.T., Xue, S., Li, L., Liu, X.C., Ding, X., and Song, M.S. (2018) TiO2 solubility and Nb and Ta partitioning in rutile-silica-rich supercritical fluid systems: Implications for subduction zone processes. Journal of Geophysical Research. Solid Earth, 123, 4765–4782, https://doi.org/10.1029/2018JB015808.Search in Google Scholar
Chen, W., Xiong, X.L., and Takahashi, E. (2021) Zircon solubility in solute-rich supercritical fluids and Zr transfer from slab to wedge in the deep subduction process. Journal of Geophysical Research: Solid Earth, 126(9), e2021JB021970.Search in Google Scholar
de Koker, N.P., Stixrude, L., and Karki, B.B. (2008) Thermodynamics, structure, dynamics, and freezing of Mg2SiO4 liquid at high pressure. Geochimica et Cos-mochimica Acta, 72, 1427–1441, https://doi.org/10.1016/j.gca.2007.12.019Search in Google Scholar
de Koker, N., Karki, B.B., and Stixrude, L. (2013) Thermodynamics of the MgO-SiO2 liquid system in Earth’s lowermost mantle from first principles. Earth and Planetary Science Letters, 361, 58–63, https://doi.org/10.1016/j.epsl.2012.11.026.Search in Google Scholar
Dingwell, D.B., Romano, C., and Hess, K.U. (1996) The effect of water on the viscosity of a haplogranitic melt under P-T-X conditions relevant to silicic volcanism. Contributions to Mineralogy and Petrology, 124, 19–28, https://doi.org/10.1007/s004100050170.Search in Google Scholar
Dingwell, D.B., Hess, K.U., and Romano, C. (1998) Viscosity data for hydrous peraluminous granitic melts: Comparison with a metaluminous model. American Mineralogist, 83, 236–239, https://doi.org/10.2138/am-1998-3-406Search in Google Scholar
Dolejs, D. and Manning, C.E. (2010) Thermodynamic model for mineral solubility in aqueous fluids: Theory, calibration and application to model fluid-flow systems. Geofluids, 10, 20–40.Search in Google Scholar
Dufils, T., Sator, N., and Guillot, B. (2020) A comprehensive molecular dynamics simulation study of hydrous magmatic liquids. Chemical Geology, 533, 119300, https://doi.org/10.1016/j.chemgeo.2019.119300Search in Google Scholar
Eyring, H. (1982) Statistical Mechanics and Dynamics. Wiley.Search in Google Scholar
Flyvbjerg, H. and Petersen, H.G. (1989) Error-estimates on averages of correlated data. The Journal of Chemical Physics, 91, 461–466, https://doi.org/10.1063/1.457480.Search in Google Scholar
Ghosh, D.B. and Karki, B.B. (2011) Diffusion and viscosity of Mg2SiO4 liquid at high pressure from first-principles simulations. Geochimica et Cosmochimica Acta, 75, 4591–4600, https://doi.org/10.1016/j.gca.2011.05.030Search in Google Scholar
Gillan, M. J., Alfè, D., and Michaelides, A. (2016) Perspective: How good is DFT for water? The Journal of Chemical Physics, 144, 130901, https://doi.org/10.1063/1.4944633.Search in Google Scholar
Grimme, S., Antony, J., Ehrlich, S., and Krieg, H. (2010) A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. The Journal of Chemical Physics, 132, 154104, https://doi.org/10.1063/1.3382344Search in Google Scholar
Grimme, S., Ehrlich, S., and Goerigk, L. (2011) Effect of the damping function in dispersion corrected density functional theory. Journal of Computational Chemistry, 32, 1456–1465, https://doi.org/10.1002/jcc.21759Search in Google Scholar
Guissani, Y., and Guillot, B. (1996) A numerical investigation of the liquid-vapor coexistence curve of silica. Journal of Chemical Physics, 104(19), 7633–7644.Search in Google Scholar
Hack, A.C., Thompson, A.B., and Aerts, M. (2007) Phase relations involving hydrous silicate melts, aqueous fluids, and minerals. Reviews in Mineralogy and Geochemistry, 65(1), 129–185.Search in Google Scholar
Harris, K.R. and Woolf, L.A. (2004) Temperature and volume dependence of the viscosity of water and heavy water at low temperatures. Journal of Chemical & Engineering Data, 49, 1064–1069, https://doi.org/10.1021/je049918mSearch in Google Scholar
Hayden, L.A., and Manning, C.E. (2011) Rutile solubility in supercritical NaAlSi3O8- H2O fluids. Chemical Geology, 284(1-2), 74–81.Search in Google Scholar
Hermann, J., Zheng, Y.-F., and Rubatto, D. (2013) Deep fluids in subducted continental crust. Elements (Quebec), 9, 281–287, https://doi.org/10.2113/gselements.9.4.281.Search in Google Scholar
Hunt, J.D. and Manning, C.E. (2012) A thermodynamic model for the system SiO2- H2O near the upper critical end point based on quartz solubility experiments at 500–1100 °C and 5–20 kbar. Geochimica et Cosmochimica Acta, 86, 196–213, https://doi.org/10.1016/j.gca.2012.03.006Search in Google Scholar
Karki, B.B. and Stixrude, L. (2010) First-principles study of enhancement of transport properties of silica melt by water. Physical Review Letters, 104, 215901, https://doi.org/10.1103/PhysRevLett.104.215901Search in Google Scholar
Karki, B.B., Bohara, B., and Stixrude, L. (2011) First-principles study of diffusion and viscosity of anorthite (CaAl2Si2O8) liquid at high pressure. American Mineralogist, 96, 744–751, https://doi.org/10.2138/am.2011.3646Search in Google Scholar
Karki, B.B., Ghosh, D.B., and Bajgain, S.K. (2018) Simulation of silicate melts under pressure. In Y. Kono and C. Sanloup, Eds., Magmas Under Pressure: Advances in High-Pressure Experiments on Structure and Properties of Melts, p. 419–453. Elsevier.Search in Google Scholar
Kennedy, G.C., Wasserburg, G.J., Heard, H.C., and Newton, R.C. (1962) The upper three-phase region in the system SiO2-H2O. American Journal of Science, 260, 501–521, https://doi.org/10.2475/ajs.260.7.501Search in Google Scholar
Kobsch, A., and Caracas, R. (2020) The critical point and the supercritical state of alkali feldspars: Implications for the behavior of the crust during impacts. Journal of Geophysical Research: Planets, 125(9), e2020JE006412.Search in Google Scholar
Kresse, G. and Furthmüller, J. (1996a) Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science, 6, 15–50, https://doi.org/10.1016/0927-0256(96)00008-0Search in Google Scholar
Kresse, G. and Furthmüller, J. (1996b) Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Physical Review B: Condensed Matter, 54, 11169–11186, https://doi.org/10.1103/PhysRevB.54.11169Search in Google Scholar
Kresse, G. and Hafner, J. (1993) Ab initio molecular dynamics for liquid metals. Physical Review B: Condensed Matter, 47, 558–561, https://doi.org/10.1103/PhysRevB.47.558.Search in Google Scholar
Kresse, G. and Hafner, J. (1994) Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. Physical Review B: Condensed Matter, 49, 14251–14269, https://doi.org/10.1103/PhysRevB.49.14251Search in Google Scholar
Li, Y., Pan, H., Liu, Q., Ming, X., and Li, Z. (2022) Ab initio mechanism revealing for tricalcium silicate dissolution. Nature Communications, 13, 1253, https://doi.org/10.1038/s41467-022-28932-2.Search in Google Scholar
Makhluf, A.R., Newton, R.C., and Manning, C.E. (2020) Experimental investigation of phase relations in the system NaAlSi3O8-H2O at high temperatures and pressures: Liquidus relations, liquid-vapor mixing, and critical phenomena at deep crust-upper mantle conditions. Contributions to Mineralogy and Petrology, 175, 76, https://doi.org/10.1007/s00410-020-01711-2Search in Google Scholar
Manning, C.E. (1994) The solubility of quartz in H2O in the lower crust and upper-mantle. Geochimica et Cosmochimica Acta, 58(22), 4831–4839.Search in Google Scholar
Manning, C.E. (2004) The chemistry of subduction-zone fluids. Earth and Planetary Science Letters, 223, 1–16, https://doi.org/10.1016/j.epsl.2004.04.030Search in Google Scholar
Manning, C.E. (2018) Fluids of the lower crust: Deep is different. Annual Review of Earth and Planetary Sciences, 46, 67–97, https://doi.org/10.1146/annurevearth-060614-105224.Search in Google Scholar
Mibe, K., Kanzaki, M., Kawamoto, T., Matsukage, K.N., Fei, Y.W., and Ono, S. (2004) Determination of the second critical end point in silicate-H2O systems using high-pressure and high-temperature X-ray radiography. Geo-chimica et Cosmochimica Acta, 68, 5189–5195, https://doi.org/10.1016/j.gca.2004.07.015.Search in Google Scholar
Mibe, K., Kanzaki, M., Kawamoto, T., Matsukage, K.N., Fei, Y.W., and Ono, S. (2007) Second critical endpoint in the peridotite-H2O system. Journal of Geophysical Research. Solid Earth, 112 (B3), 8.Search in Google Scholar
Mibe, K., Chou, I.M., and Bassett, W.A. (2008) In situ Raman spectroscopic investigation of the structure of subduction-zone fluids. Journal of Geophysical Research. Solid Earth, 113, (B4).Search in Google Scholar
Mibe, K., Kawamoto, T., Matsukage, K.N., Fei, Y., and Ono, S. (2011) Slab melting versus slab dehydration in subduction-zone magmatism. Proceedings of the National Academy of Sciences, 108, 8177–8182, https://doi.org/10.1073/pnas. 1010968108.Search in Google Scholar
Momma, K. and Izumi, F. (2011) VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. Journal of Applied Crystallography, 44, 1272–1276, https://doi.org/10.1107/S0021889811038970Search in Google Scholar
Mysen, B.O. (2007) The solution behavior of H2O in peralkaline aluminosilicate melts at high pressure with implications for properties of hydrous melts. Geochimica et Cosmochimica Acta, 71, 1820–1834, https://doi.org/10.1016/j.gca.2007.01.007.Search in Google Scholar
Mysen, B.O. (2010) Structure of H2O-saturated peralkaline aluminosilicate melt and coexisting aluminosilicate-saturated aqueous fluid determined in-situ to 800 °C and 800 MPa. Geochimica et Cosmochimica Acta, 74, 4123–4139, https://doi.org/10.1016/j.gca.2010.04.024.Search in Google Scholar
Mysen, B.O., Mibe, K., Chou, I.M., and Bassett, W.A. (2013) Structure and equilibria among silicate species in aqueous fluids in the upper mantle: Experimental SiO2-H2O and MgO-SiO2-H2O data recorded in situ to 900 °C and 5.4 GPa. Journal of Geophysical Research. Solid Earth, 118, 6076–6085, https://doi.org/10.1002/2013JB010537.Search in Google Scholar
Nakamura, Y. (1974) The system SiO2-H2O-H2 at 15 kbar. Carnegie Institution of Washington Yearbook, 73, 259–263.Search in Google Scholar
Newton, R.C. and Manning, C.E. (2008) Thermodynamics of SiO2-H2O fluid near the upper critical end point from quartz solubility measurements at 10 kbar. Earth and Planetary Science Letters, 274, 241–249, https://doi.org/10.1016/j.epsl.2008.07.028.Search in Google Scholar
Ni, H.W., Zhang, L., Xiong, X.L., Mao, Z., and Wang, J.Y. (2017) Supercritical fluids at subduction zones: Evidence, formation condition, and physicochemi-cal properties. Earth-Science Reviews, 167, 62–71, https://doi.org/10.1016/j.earscirev.2017.02.006.Search in Google Scholar
Perdew, J.P., Burke, K., and Ernzerhof, M. (1996) Generalized gradient approximation made simple. Physical Review Letters, 77, 3865–3868, https://doi.org/10.1103/PhysRevLett.77.3865.Search in Google Scholar
Romano, C., Poe, B., Mincione, V., Hess, K.U., and Dingwell, D.B. (2001) The viscosities of dry and hydrous XAlSi3O8 (X = Li, Na, K, Ca0.5, Mg0.5) melts. Chemical Geology, 174, 115–132, https://doi.org/10.1016/S0009-2541(00)00311-9.Search in Google Scholar
Shen, A.H. and Keppler, H. (1997) Direct observation of complete miscibility in the albite-H2O system. Nature, 385, 710–712, https://doi.org/10.1038/385710a0.Search in Google Scholar
Spandler, C., Mavrogenes, J., and Hermann, J. (2007) Experimental constraints on element mobility from subducted sediments using high-P synthetic fluid/melt inclusions. Chemical Geology, 239, 228–249, https://doi.org/10.1016/j.chemgeo.2006.10.005.Search in Google Scholar
Spiekermann, G., Wilke, M., and Jahn, S. (2016) Structural and dynamical properties of supercritical H2O-SiO2 fluids studied by ab initio molecular dynamics. Chemical Geology, 426, 85–94, https://doi.org/10.1016/j.chemgeo.2016.01.010Search in Google Scholar
Stalder, R., Ulmer, P., Thompson, A.B., and Gunther, D. (2000) Experimental approach to constrain second critical end points in fluid/silicate systems: Near-solidus fluids and melts in the system albite-H2O. American Mineralogist, 85, 68–77, https://doi.org/10.2138/am-2000-0108Search in Google Scholar
Steele-MacInnis, M. and Schmidt, C. (2014) Silicate speciation in H2O-Na2O-SiO2 fluids from 3 to 40 mol% SiO2, to 600 °C and 2 GPa. Geochimica et Cosmochimica Acta, 136, 126–141, https://doi.org/10.1016/j.gca.2014.04.009Search in Google Scholar
Sun, Y., Zhou, H., Yin, K., Zhao, M., Xu, S., and Lu, X. (2018) Transport properties of Fe2SiO4 melt at high pressure from classical molecular dynamics: Implications for the lifetime of the magma ocean. Journal of Geophysical Research. Solid Earth, 123, 3667–3679, https://doi.org/10.1029/2018JB015452Search in Google Scholar
Sun, Y., Zhou, H., Liu, X., Yin, K., and Lu, X. (2020) Physical state of an early magma ocean constrained by the thermodynamics and viscosity of iron silicate liquid. Earth and Planetary Science Letters, 551, 116556, https://doi.org/10.1016/j.epsl.2020.116556.Search in Google Scholar
Thomas, R., Davidson, P., and Appel, K. (2019) The enhanced element enrichment in the supercritical states of granite-pegmatite systems. Acta Geochimica, 38, 335–349, https://doi.org/10.1007/s11631-019-00319-zSearch in Google Scholar
Wang, Y., Sakamaki, T., Skinner, L.B., Jing, Z., Yu, T., Kono, Y., Park, C., Shen, G., Rivers, M.L., and Sutton, S.R. (2014) Atomistic insight into viscosity and density of silicate melts under pressure. Nature Communications, 5, 3241, https://doi.org/10.1038/ncomms4241Search in Google Scholar
Wang, Q.X., Zhou, D.Y., Li, W.C., and Ni, H.W. (2021) Spinodal decomposition of supercritical fluid forms melt network in a silicate-H2O system. Geochemical Perspectives Letters, 18, 22–26, https://doi.org/10.7185/geochemlet.2119Search in Google Scholar
Xiong, X.L., Keppler, H., Audetat, A., Gudfinnsson, G., Sun, W.D., Song, M.S., Xiao, W.S., and Yuan, L. (2009) Experimental constraints on rutile saturation during partial melting of metabasalt at the amphibolite to eclogite transition, with applications to TTG genesis. American Mineralogist, 94, 1175–1186, https://doi.org/10.2138/am.2009.3158Search in Google Scholar
Yang, F., Hess, K.U., Unruh, T., Mamontov, E., Dingwell, D.B., and Meyer, A. (2017) Intrinsic proton dynamics in hydrous silicate melts as seen by quasielas-tic neutron scattering at elevated temperature and pressure. Chemical Geology, 461, 152–159, https://doi.org/10.1016/j.chemgeo.2017.01.009Search in Google Scholar
Yuan, L., Steinle-Neumann, G., and Suzuki, A. (2020) Structure and density of H2O-rich Mg2SiO4 melts at high pressure from ab initio simulations. Journal of Geophysical Research: Solid Earth, 125(10), e2020JB020365.Search in Google Scholar
Zotov, N. and Keppler, H. (2002) Silica speciation in aqueous fluids at high pressures and high temperatures. Chemical Geology, 184, 71–82, https://doi.org/10.1016/S0009-2541(01)00353-9.Search in Google Scholar
© 2023 by Mineralogical Society of America
Articles in the same Issue
- Heavy halogen compositions of lamprophyres derived from metasomatized lithospheric mantle beneath eastern North China Craton
- Compositional trends in Ba-, Ti-, and Cl-rich micas from metasomatized mantle rocks of the Gföhl Unit, Bohemian Massif, Austria
- Experimental determination of quartz solubility in H2O-CaCl2 solutions at 600–900 °C and 0.6–1.4 GPa
- The use of boron nitride to impose reduced redox conditions in experimental petrology
- Structures and transport properties of supercritical SiO2-H2O and NaAlSi3O8-H2O fluids
- Hydrologic regulation of clay-mineral transformations in a redoximorphic soil of subtropical monsoonal China
- Witness to strain: Subdomain boundary length and the apparent subdomain boundary density in large strained olivine grains
- Libyan Desert Glass: New evidence for an extremely high-pressure-temperature impact event from nanostructural study
- Crystal vs. melt compositional effects on the partitioning of the first-row transition and high field strength elements between clinopyroxene and silicic, alkaline, aluminous melts
- Microbially induced clay weathering: Smectite-to-kaolinite transformation
- Hydrous wadsleyite crystal structure up to 32 GPa
- Multiple fluid sources in skarn systems: Oxygen isotopic evidence from the Haobugao Zn-Fe-Sn deposit in the southern Great Xing’an Range, NE China
- Crocobelonite, CaFe23+(PO4)2O, a new oxyphosphate mineral, the product of pyrolytic oxidation of natural phosphides
- Tetrahedrite-(Ni), Cu6(Cu4Ni2)Sb4S13, the first nickel member of tetrahedrite group mineral from Luobusa chromite deposits, Tibet, China
- New Mineral Names: Heavy metal and minerals from China
- Book Review
Articles in the same Issue
- Heavy halogen compositions of lamprophyres derived from metasomatized lithospheric mantle beneath eastern North China Craton
- Compositional trends in Ba-, Ti-, and Cl-rich micas from metasomatized mantle rocks of the Gföhl Unit, Bohemian Massif, Austria
- Experimental determination of quartz solubility in H2O-CaCl2 solutions at 600–900 °C and 0.6–1.4 GPa
- The use of boron nitride to impose reduced redox conditions in experimental petrology
- Structures and transport properties of supercritical SiO2-H2O and NaAlSi3O8-H2O fluids
- Hydrologic regulation of clay-mineral transformations in a redoximorphic soil of subtropical monsoonal China
- Witness to strain: Subdomain boundary length and the apparent subdomain boundary density in large strained olivine grains
- Libyan Desert Glass: New evidence for an extremely high-pressure-temperature impact event from nanostructural study
- Crystal vs. melt compositional effects on the partitioning of the first-row transition and high field strength elements between clinopyroxene and silicic, alkaline, aluminous melts
- Microbially induced clay weathering: Smectite-to-kaolinite transformation
- Hydrous wadsleyite crystal structure up to 32 GPa
- Multiple fluid sources in skarn systems: Oxygen isotopic evidence from the Haobugao Zn-Fe-Sn deposit in the southern Great Xing’an Range, NE China
- Crocobelonite, CaFe23+(PO4)2O, a new oxyphosphate mineral, the product of pyrolytic oxidation of natural phosphides
- Tetrahedrite-(Ni), Cu6(Cu4Ni2)Sb4S13, the first nickel member of tetrahedrite group mineral from Luobusa chromite deposits, Tibet, China
- New Mineral Names: Heavy metal and minerals from China
- Book Review