An official website of the United States government
Information resulting from a comprehensive investigation into the intrinsic strengths of hydrated divalent magnesium clusters is useful for elucidating the role of aqueous solvents on the Mg2+ ion, which can be related to those in bulk aqueous solution. However, the intrinsic Mg–O and intermolecular hydrogen bond interactions of hydrated magnesium ion clusters have yet to be quantitatively measured. In this work, we investigated a set of 17 hydrated divalent magnesium clusters by means of local vibrational mode force constants calculated at the ωB97X-D/6-311++G(d,p) level of theory, where the nature of the ion–solvent and solvent–solvent interactions were interpreted from topological electron density analysis and natural population analysis. We found the intrinsic strength of inner shell Mg–O interactions for [Mg(H2O)n]2+ (n = 1–6) clusters to relate to the electron density at the bond critical point in Mg–O bonds. From the application of a secondary hydration shell to [Mg(H2O)n]2+ (n = 5–6) clusters, stronger Mg–O interactions were observed to correspond to larger instances of charge transfer between the lp(O) orbitals of the inner hydration shell and the unfilled valence shell of Mg. As the charge transfer between water molecules of the first and second solvent shell increased, so did the strength of their intermolecular hydrogen bonds (HBs). Cumulative local vibrational mode force constants of explicitly solvated Mg2+, having an outer hydration shell, reveal a CN of 5, rather than a CN of 6, to yield slightly more stable configurations in some instances. However, the cumulative local mode stretching force constants of implicitly solvated Mg2+ show the six-coordinated cluster to be the most stable. These results show that such intrinsic bond strength measures for Mg–O and HBs offer an effective way for determining the coordination number of hydrated magnesium ion clusters.
more »
« less
Structural and Thermodynamic Properties of Magnesium-Rich Liquids at Ultrahigh Pressure
We explore the structural properties of Mg, MgO, and MgSiO3 liquids from ab initio computer simulations at conditions that are relevant for the interiors of giant planets, stars, shock compression measurements, and inertial confinement fusion experiments. Using path-integral Monte Carlo and density functional theory molecular dynamics, we derive the equation of state of magnesium-rich liquids in the regime of condensed and warm dense matter, with densities ranging from 0.32 to 86.11 g cm−3 and temperatures from 20,000 K to 5 × 108 K. We study the electronic structure of magnesium as a function of density and temperature and the correlations of the atomic motion, finding an unexpected local maximum in the pair correlation functions that emerges at high densities which decreases the coordination number of elemental magnesium and reveals a higher packing. This phenomenon is not observed in other magnesium liquids, which maintain a rather constant coordination number.
more »
« less
- Award ID(s):
- 2020249
- PAR ID:
- 10431055
- Date Published:
- Journal Name:
- Minerals
- Volume:
- 13
- Issue:
- 7
- ISSN:
- 2075-163X
- Page Range / eLocation ID:
- 885
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Thermophysical properties of silicate liquids under extreme conditions are critical for understanding the accretion and evolution of super‐Earths and sub‐Neptunes. The thermal equation of state and viscosity of silicate liquids determine the adiabatic profiles and dynamics of magma oceans. However, these properties are challenging to constrain at elevated pressures in experiments. Here, we perform ab initio molecular dynamics simulations of MgSiO3liquid across a wide range of pressures (0–1,200 GPa) and temperatures (2200–14000 K) and analyze its structure, the Grüneisen parameter, and viscosity. Our results reveal the clear temperature and pressure dependence of the Grüneisen parameter, which varies synchronously with the O‐O coordination number. The Grüneisen parameter shifts from positive to negative temperature dependence between ∼20 and 70 GPa, corresponding to a peak in the O‐O coordination number and SiO5abundance. Initially, the Grüneisen parameter increases with pressure and then decreases, showing limited temperature dependence above ∼300 GPa, where its behavior resembles that of solids. Furthermore, we determine the adiabat and viscosity profiles of magma oceans in super‐Earths and sub‐Neptunes. The results suggest that the mantles of super‐Earths and sub‐Neptunes may solidify either from the bottom up or at pressures of ∼120–150 GPa, depending on the curvature of the mantle melting line. The low viscosity of magma oceans likely enhances convective currents and facilitate efficient differentiation. These thermophysical properties, now quantified up to terapascal pressures, enable updates to the mass‐radius relation of magma ocean exoplanets, showing notable differences compared to their solid counterparts.more » « less
-
While the excess chemical potential is the key quantity in determining phase diagrams, its direct computation for high-density liquids of long polymer chains has posed a significant challenge. Computationally, the excess chemical potential is calculated using the Widom insertion method, which involves monitoring the change in internal energy as one incrementally introduces individual molecules into the liquid. However, when dealing with dense polymer liquids, inserting long chains requires generating trial configurations with a bias that favors those at low energy on a unit-by-unit basis: a procedure that becomes more challenging as the number of units increases. Thus, calculating the excess chemical potential of dense polymer liquids using this method becomes computationally intractable as the chain length exceeds N ≥ 30. Here, we adopt a coarsegrained model derived from the integral equation theory for which inserting long polymer chains becomes feasible. The integral equation theory of coarse graining (IECG) represents a polymer as a sphere or a collection of blobs interacting through a soft potential. We employ the IECG approach to compute the excess chemical potential using Widom’s method for polymer chains of increasing lengths, extending up to N = 720 monomers, and at densities reaching up to ρ = 0.767 g/cm3. From a fundamental perspective, we demonstrate that the excess chemical potentials remain nearly constant across various levels of coarse graining, offering valuable insights into the consistency of this type of procedure. Ultimately, we argue that current Monte Carlo algorithms, originally designed for atomistic simulations, such as configurational bias Monte Carlo (CBMC) methods, can significantly benefit from the integration of the IECG approach, thereby enhancing their performance in the study of phase diagrams of polymer liquids.more » « less
-
Mg 2 GeO 4 is important as an analog for the ultrahigh-pressure behavior of Mg 2 SiO 4 , a major component of planetary interiors. In this study, we have investigated magnesium germanate to 275 GPa and over 2,000 K using a laser-heated diamond anvil cell combined with in situ synchrotron X-ray diffraction and density functional theory (DFT) computations. The experimental results are consistent with the formation of a phase with disordered Mg and Ge, in which germanium adopts eightfold coordination with oxygen: the cubic, Th 3 P 4 -type structure. DFT computations suggest partial Mg-Ge order, resulting in a tetragonal I 4 ¯ 2 d structure indistinguishable from I 4 ¯ 3 d Th 3 P 4 in our experiments. If applicable to silicates, the formation of this highly coordinated and intrinsically disordered phase may have important implications for the interior mineralogy of large, rocky extrasolar planets.more » « less
-
The utilization of metallic anodes holds promise for unlocking high gravimetric and volumetric energy densities and is pivotal to the adoption of ‘beyond Li’ battery chemistries. Much of the promise of magnesium batteries stems from claims regarding their lower predilection for dendrite growth. Whilst considerable effort has been invested in the design of novel electrolytes and cathodes, detailed studies of Mg plating are scarce. Using galvanostatic electrodeposition of metallic Mg from Grignard reagents in symmetric Mg–Mg cells, we establish a phase map characterized by disparate morphologies spanning the range from fractal aggregates of 2D nanoplatelets to highly anisotropic dendrites with singular growth fronts and nanowires entangled in the form of mats. The effects of electrolyte concentration, applied current density, and coordinating ligands have been explored. The study demonstrates a complex range of electrodeposited morphologies including canonical dendrites with shear moduli conducive to penetration through typical polymeric separators. We further demonstrate a strategy for mitigating Mg dendrite formation based on the addition of molecular Lewis bases that promote nanowire growth through selective surface coordination.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/min13070885
