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.
Attention:The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 7:00 AM ET to 7:30 AM ET on Friday, April 24 due to maintenance. We apologize for the inconvenience.

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. ; ; ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract Viscosity of silicate melts governs magma transport and influences mantle dynamics, yet effects of pressure and water on melt viscosity remain poorly understood. Here, we report in situ falling‐sphere viscosity measurements on diopside (Di) melts with 0–3 wt.% H2O along the liquidus up to 7 GPa and 2103 K using synchrotron X‐ray radiography. By incorporating our hydrous melt data into a previously validated model for the dry system, the effects of pressure, temperature, and H2O contents on Di melt viscosity can be satisfactorily captured by the function: whereT*is the homologous temperature,xH2Ois the molar % H2O,η0 = 8.90 (1.50) × 10−8 Pa s,b0 = 3.02 (0.10), andH*(P) = 15.72 (0.03)−0.35 (0.01)·P + 1.07 (0.07) × 10−2·P2−1.19 (0.14) × 10−4 P3, ×10−3 GPa−1. Adding 3 wt.% H2O systematically reduces viscosity by ∼0.7 log units. For both dry and hydrous melts, viscosity along the liquidus decreases monotonically with increasing pressure, suggesting that moderate hydration may not significantly alter the compressional behavior of Di melts. Combining the Di viscosity model with models for feldspar and olivine, we simulated the viscosity of analog basaltic magmas under mantle conditions. Increasing H2O content from 0 to 3wt.% raises mobility of basaltic magma increases by >1 order of magnitude. In hot plume settings, the mobility further increases by a factor of 30 relative to typical ambient mantle. Assuming a simple percolation model, the increased mobility corresponds to faster melt ascent in mantle plumes that could, in part, explain the voluminous magmatism of large igneous provinces. 
    more » « less
  3. ; ; ; ; (, 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
  4. ; ; (, 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
  5. ; ; ; ; (, 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
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email