skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: First-principles computation of diffusional Mg isotope fractionation in silicate melts
Diffusional isotope fractionation occurs in geochemical processes (such as magma mixing, bubble growth, and crystal growth), even at magmatic temperatures. Isotopic mass dependence of diffusion is commonly expressed as Di Dj ¼ mj mi   b , where Di and Dj are diffusion coefficients of two isotopes whose masses are mi and mj. How the dimensionless empirical parameter b depends on temperature, pressure, and composition remains poorly constrained. Here, we conducted a series of first-principles molecular dynamics simulations to evaluate the b factor of Mg isotopes in MgSiO3 and Mg2SiO4 melts using pseudo-isotope method. In particular, we considered interactions between Mg isotopes by simultaneously putting pseudo-mass and normalmass Mg atoms in a simulation supercell. The calculated b for Mg isotopes decreases linearly with decreasing temperature at zero pressure, from 0:158  0:004 at 4000 K to 0:121  0:017 at 2200 K for MgSiO3 melt and from 0:150  0:004 at 4000 K to 0:101  0:012 at 2200 K for Mg2SiO4 melt. Moreover, our simulations of compressed Mg2SiO4 melt along the 3000 K isotherm show that the b value decreases linearly from 0:130  0:006 at 0 GPa to 0:060  0:011 at 17 GPa. Based on our diffusivity results, the empirically established positive correlation between b and solvent-normalized diffusivity (Di/DSi) seems to be applicable only at constant temperatures or in narrow temperature ranges. Analysis of atomistic mechanisms suggests that the calculated b values are inversely correlated with force constants of Mg at a given temperature or pressure. Good agreement between our first principles results with available experimental data suggests that interactions between isotopes of major elements must be considered in calculating b for major elements in silicate melts. Also, we discuss diffusion-controlled crystal growth by considering our calculated b values.  more » « less
Award ID(s):
1764140
PAR ID:
10296624
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geochimica et cosmochimica acta
Volume:
290
ISSN:
0016-7037
Page Range / eLocation ID:
27-40
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Geochimica et cosmochimica acta)
    null (Ed.)
    Diffusional isotope fractionation has been widely used to explain lithium (Li) isotope variations in minerals and rocks. Isotopic mass dependence of Li diffusion can be empirically expressed as , where is the diffusivity of a Li isotope. The knowledge about temperature and compositional dependence of the factor which is essential for understanding diffusion profiles and mechanisms remains unclear. Based on the potential energy and interatomic forces generated by deep neural networks trained with ab initio data, we performed deep potential molecular dynamics (DPMD) simulations of several Li pseudo-isotopes (with mass = 2, 7, 21, 42 g/mol) in albite, hydrous albite, and model basalt melts to evaluate the factor. Our calculated diffusivities for 7Li in albite and model basalt melts at 1800 K compare well with experimental results. We found that in albite melt decreases from at 4000 K to at 1800 K. The presence of water appears to slightly weaken the temperature dependence of , with decreasing from to in hydrous albite melt. The calculated in model basalt melt takes much smaller values, decreasing from at 4000 K to at 1800 K. Our prediction of in albite and hydrous albite melts is in good agreement with experimental data. More importantly, our results suggest that Li isotope diffusion in silicate melts is strongly dependent on melt composition. The temperature and compositional effects on can be qualitatively explained in terms of ionic porosity and the coupled relationship between Li diffusion and the mobility of the silicate melt network. Two types of diffusion experiments are suggested to test our predicted temperature and compositional dependence of . This study shows that DPMD is a promising tool to simulate the diffusion of elements and isotopes in silicate melts. 
    more » « less
  2. ; ; ; ; (, ACS earth and space chemistry)
    In this study, we use f irst-principles molecular dynamics simulations to explore the behavior of anhydrous aluminosilicate melt with a stoichiometry of NaAlSi2O6 up to pressures of ∼30 GPa and temperatures between 2500 and 4000 K. We also examine the effect of water (∼4 wt % H2O) on the equation of state and transport properties of the aluminosilicate melt and relate them to atomistic scale changes in the melt structure. Our results show that water reduces the density and bulk modulus of the anhydrous melt. However, the pressure derivative of the bulk modulus of the hydrous melt is larger than that of the anhydrous melt. The pressure dependence of the transport property exhibits an anomalous behavior. At a pressure of ∼12 GPa, anhydrous aluminosilicate melts exhibit maxima in diffusion and minima in viscosity. Dissolved water in melts also affects both diffusion and viscosity. In hydrous aluminosilicate melts, the maxima in diffusion and the minima in viscosity occur at ∼14 GPa. The anomalous behavior of transport properties is related to the pressure-induced changes in the melt structure. At shallower depths, i.e., up to 100 km, relevant for subduction zone settings, the lower density compounded by the lower viscosity of hydrous aluminosilicate melts is likely to provide buoyancy for upward migration. At greater depths of ∼180−200 km, greater compressibility of the hydrous aluminosilicate melts together with the minimum viscosity could hinder magma migration and may explain the presence of a partial melt layer at the lithosphere−asthenosphere boundary. 
    more » « less
  3. ; ; (, Geochemistry, Geophysics, Geosystems)
    Abstract The continental crust is produced by the solidification of aluminosilicate‐rich magmas which are sourced from deep below the surface. Migration of the magma depends on the density (ρ) contrast to source rocks and the melt viscosity (η). At the surface, these silica‐rich melts are typically sluggish due to highη > 1,000 Pa s. Yet at their source regions, the melt properties are complexly influenced by pressure (P), temperature (T), and water contents (). In this study, we examined the combinedP‐T‐ effects on the behavior of melts with an albite stoichiometry (NaAlSi3O8). We usedfirst‐principlesmolecular dynamics simulations to examine anhydrous (0 wt % H2O) and hydrous (5 wt % H2O) melts. To constrain thePandTeffects, we exploredP ≤ 25 GPa across several isotherms between 2500 and 4000 K. The melts show anomalousP‐ρrelationships at lowP ∼ 0 GPa and highT ≥ 2500 K, consistent with vaporization. At lithospheric conditions, meltρincreases with compression and is well described by a finite‐strain formalism. Water lowers the melt density (ρhydrous < ρanhydrous) but increases the compressibility, that is, 1/Khydrous>1/KanhydrousorKhydrous < Kanhydrous. We also find that the meltηdecreases with pressure and then increases with further compression. Water decreases the viscosity (ηhydrous < ηanhydrous) by depolymerizing the melt structure. The ionic self‐diffusivities are increased by the presence of water. The decreasedρandηby H2O increase the mobility of magma at crustal conditions, which could explain the rapid eruption and migration timescales for rhyolitic magmas as observed in the Chaitén volcano in Chile. 
    more » « less
  4. ; ; ; ; (, Minerals)
    null (Ed.)
    The bridgmanite–akimotoite–majorite (Bm–Ak–Mj or BAM) triple point in MgSiO3 has been measured in large-volume press (LVP; COMPRES 8/3 assembly) and laser-heated diamond anvil cell (LHDAC). For the LVP data, we calculated pressures from the calibration provided for the assembly. For the LHDAC data, we conducted in situ determination of pressure at high temperature using the Pt scale at synchrotron. The measured temperatures of the triple point are in good agreement between LVP and LHDAC at 1990–2000 K. However, the pressure for the triple point determined from the LVP is 3.9 ± 0.6 GPa lower than that from the LHDAC dataset. The BAM triple point determined through these experiments will provide an important reference point in the pressure–temperature space for future high-pressure experiments and will allow mineral physicists to compare the pressure–temperature conditions measured in these two different experimental methods. 
    more » « less
  5. ; ; (, Scientific reports)
    null (Ed.)
    Water (H2O) as one of the most abundant fluids present in Earth plays crucial role in the generation and transport of magmas in the interior. Though hydrous silicate melts have been studied extensively, the experimental data are confined to relatively low pressures and the computational results are still rare. Moreover, these studies imply large differences in the way water influences the physical properties of silicate magmas, such as density and electrical conductivity. Here, we investigate the equation of state, speciation, and transport properties of water dissolved in Mg1-xFexSiO3 and Mg2(1-x)Fe2xSiO4 melts (for x = 0 and 0.25) as well as in its bulk (pure) fluid state over the entire mantle pressure regime at 2000 to 4000 K using first-principles molecular dynamics. The simulation results allow us to constrain the partial molar volume of the water component in melts along with the molar volume of pure water. The predicted volume of silicate melt+water solution is negative at low pressures and becomes zero above 15 GPa. Consequently, the hydrous component tends to lower the melt density to similar extent over much of the mantle pressure regime irrespective of composition. Our results also show that hydrogen diffuses fast in silicate melts and enhances the melt electrical conductivity in a way that differs from electrical conduction in the bulk water. The speciation of the water component varies considerably from the bulk water structure as well. Water is dissolved in melts mostly as hydroxyls at low pressure and as -O-H-O-, -O-H-O-H- and other extended species with increasing pressure. On the other hand, the pure water behaves as a molecular fluid below 15 GPa, gradually becoming a dissociated fluid with further compression. On the basis of modeled density and conductivity results, we suggest that partial melts containing a few percent of water may be gravitationally trapped both above and below the upper mantle-transition region. Moreover, such hydrous melts can give rise to detectable electrical conductance by means of electromagnetic sounding observations. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email