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.


Title: Quantum kinetic energy and isotope fractionation in aqueous ionic solutions
At room temperature, the quantum contribution to the kinetic energy of a water molecule exceeds the classical contribution by an order of magnitude. The quantum kinetic energy (QKE) of a water molecule is modulated by its local chemical environment and leads to uneven partitioning of isotopes between different phases in thermal equilibrium, which would not occur if the nuclei behaved classically. In this work, we use ab initio path integral simulations to show that QKEs of the water molecules and the equilibrium isotope fractionation ratios of the oxygen and hydrogen isotopes are sensitive probes of the hydrogen bonding structures in aqueous ionic solutions. In particular, we demonstrate how the QKE of water molecules in path integral simulations can be decomposed into translational, rotational and vibrational degrees of freedom, and use them to determine the impact of solvation on different molecular motions. By analyzing the QKEs and isotope fractionation ratios, we show how the addition of the Na + , Cl − and HPO 4 2− ions perturbs the competition between quantum effects in liquid water and impacts their local solvation structures.  more » « less
Award ID(s):
1904800 1652960
PAR ID:
10147976
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
ISSN:
1463-9076
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Stalagmites are an important archive of terrestrial climate information. However, there remains questions about the ability of stalagmites to form in oxygen isotopic equilibrium and thus record, in a simple manner, the oxygen isotopic composition and temperature of formational fluids. Recent studies have suggested that the combined application of 48 and 47 carbonate clumped isotope measurements can quantify the extent of kinetic isotope fractionation in stalagmites and thus used to correct for these kinetic isotope effects and solve for the original formation temperatures. Here we measure the 47 and 48 values from 16 different samples of the same stalagmite from central California that spans the deglaciation (11 to 20 kya). Each sample is replicated three to five times. We find that based on these measurements the extent of kinetic fractionation present in the carbonate from this stalagmite is minimal. The temperature calculated from 47 in this stalagmite ranges from 11.6 to 19.9 °C, in agreement with regional reconstructions of temperatures from 47 values of lake carbonates. In contrast, previously published the 18O and 2H values of inclusion fluids (Wortham et al., 2022) from this stalagmite suggest periods of increasing kinetic fractionation of the water isotopes at 13 and 15 ka. These periods have been previously interpreted to be times of a reduction in effective moisture regionally. We suggest by this comparison that the use of both water isotopes and the dual clumped isotope system in stalagmites can aide the interpretation of where kinetic fractionation occurs in the hydrologic and carbonate precipitation system in caves. We will discuss the work’s implications for paleoclimate records from stalagmites and other terrestrial systems in seasonally dry and Mediterranean regions. 
    more » « less
  2. Copper shows limited isotopic variation in equilibrated mantle-derived silicate rocks, but large isotopic fractionation during kinetic processes. For example, lunar and terrestrial samples that have experienced evaporation were found to have an isotopic fractionation of up to 12.5‰ in their 65Cu/63Cu ratios, while komatiites, lherzolites, mid-ocean ridge and ocean island basalts show negligible Cu isotope fractionation as a result of equilibrium partial melting and crystal fractionation. The contrast between the observed magnitudes of equilibrium and kinetic isotope fractionation for Cu calls for a better understanding of kinetic Cu isotope fractionation. One of the mechanisms for creating large kinetic isotopic fractionation even at magmatic temperatures is diffusion. In this study, we performed Cu isotopic measurements on Cu diffusion couple experiments to constrain the beta factor for Cu isotopic fractionation by diffusion. We demonstrate a Monte Carlo approach for the regression and error estimation of the measured isotope profiles, which yielded beta values of 0.16 ± 0.03 and 0.18 ± 0.03 for the two experimental charges measured. Our results are subsequently applied to a quantitative model for the evaporation of a molten sphere to discuss the role of diffusion in affecting the bulk Cu isotopic fractionation between liquid and vapor during evaporation. We apply the model to Cu evaporation experiments and tektite data to show that convection primarily governs mass transport for evaporation during tektite formation. In addition, we show that Cu isotopes can be used as a tool to test the role of kinetics during various magmatic processes such as magmatic sulfide ore deposit formation, porphyry-type ore deposit formation, and fluid-rock interactions. 
    more » « less
  3. This dataset compiles the results of our study of nuclear quantum effects (NQEs) and equilibrium isotope effects (EIEs) of 92 chemically diverse, organic molecular liquids. It contains the average macroscopic properties (density, molar volume, thermal expansion coefficient, isothermal compressibility, dielectric constant, heat of vaporization) and their associated standard errors computed with four independent classical and path-integral molecular dynamics (PIMD) simulations of each system and their deuterated counterparts. All simulations use the Topology Automated Force-Field Interactions (TAFFI) framework to describe the potential energy surface. In addition, the computed nuclear quantum effects and equilibrium isotope effects resulting from comparison of these simulations are included. 
    more » « less
  4. Titanium (Ti) typically exhibits low mobility in geologic fluids due to the low aqueous solubility of common (Fe-)Ti oxide minerals. Consequently, Ti isotope variations (δ49/47Ti, given as δ49Ti) in geologic systems are primarily attributed to magmatic differentiation. Thus, the potential for fluid-mineral fractionation has received less attention. However, ligand-rich fluids are capable of mobilizing Ti as observed in natural systems and laboratory studies. As hydrothermal ore mineralization is commonly associated with ligand-rich brines capable of transporting significant quantities of metals, Ti isotopes may aid in understanding mineralization and alteration in complex hydrothermal systems. Here we present data from computational modeling of various Ti coordination complexes theorized to exist in geologic systems and/or under relevant experimental conditions as well as computed fractionation factors for various Ti-bearing crystalline phases to investigate the basic mechanics of equilibrium fluid-mineral Ti isotope fractionation. These results indicate that equilibrium fluid-mineral Ti isotope exchange between our modeled Ti complexes and phases with 6-coordinated Ti is predicted to generally lead to enrichment of heavy Ti isotopes in the fluid. Because minerals with 6-coordinated Ti (such as magnetite and ilmenite) are the most important reservoirs of Ti in the solid Earth, Ti isotope equilibration between terrestrial rocks and fluids can be generalized to enrich the fluid in heavy Ti isotopes. We also performed magnetite-ülvospinel leaching experiments to investigate fluid-mineral Ti isotope fractionation in this phase. Mineral leaching experiments varying acid strength, leaching temperature, and reaction time with HCl and HF qualitatively support the prediction that the fluid phase will become enriched in heavy Ti isotopes during fluid-mineral interactions that approach equilibrium with Ti-rich magnetite. Additionally, the leaching data also suggest that the fluid becomes slightly enriched in lighter Ti isotopes when Ti exchange is limited—potentially due to kinetic effects. Therefore, magnetite from natural systems may be depleted in heavy Ti isotopes during regenerative mineral replacement involving equilibration with fluids or may possibly become depleted in light Ti isotopes under a kinetic fractionation regime—leading to mineral δ49Ti values that are insufficiently explained by magmatic differentiation or inter-mineral fractionation. These results are a first look at fluid-mineral interactions that may affect Ti isotope fractionation in hydrothermal mineral systems, and Ti isotopes should be further studied as a potential method of understanding aqueous metal transport and tracing alteration in mineral deposits. 
    more » « less
  5. The abundance and isotopic composition of noble gases dissolved in water have many applications in the geosciences. In recent years, new analytical techniques have opened the door to the use of high-precision measurements of noble gas isotopes as tracers for groundwater hydrology, oceanography, mantle geochemistry, and paleoclimatology. These analytical advances have brought about new measurements of solubility equilibrium isotope effects (SEIEs) in water (i.e., the relative solubilities of noble gas isotopes) and their sensitivities to the temperature and salinity. Here, we carry out a suite of classical molecular dynamics (MD) simulations and employ the theoretical method of quantum correction to estimate SEIEs for comparison with experimental observations. We find that classical MD simulations can accurately predict SEIEs for the isotopes of Ar, Kr, and Xe to order 0.01‰, on the scale of analytical uncertainty. However, MD simulations consistently overpredict the SEIEs of Ne and He by up to 40% of observed values. We carry out sensitivity tests at different temperatures, salinities, and pressures and employ different sets of interatomic potential parameters and water models. For all noble gas isotopes, the TIP4P water model is found to reproduce observed SEIEs more accurately than the SPC/E and TIP4P/ice models. Classical MD simulations also accurately capture the sign and approximate magnitude of temperature and salinity sensitivities of SEIEs for heavy noble gases. We find that experimental and modeled SEIEs generally follow an inverse-square mass dependence, which implies that the mean-square force experienced by a noble gas atom within a solvation shell is similar for all noble gases. This inverse-square mass proportionality is nearly exact for Ar, Kr, and Xe isotopes, but He and Ne exhibit a slightly weaker mass dependence. We hypothesize that the apparent dichotomy between He–Ne and Ar–Kr–Xe SEIEs may result from atomic size differences, whereby the smaller noble gases are more likely to spontaneously fit within cavities of water without breaking water–water H-bonds, thereby experiencing softer collisions during translation within a solvation shell. We further speculate that the overprediction of simulated He and Ne SEIEs may result from the neglection of higher-order quantum corrections or the overly stiff representation of van der Waals repulsion by the widely used Lennard-Jones 6–12 potential model. We suggest that new measurements of SEIEs of heavy and light noble gases may represent a novel set of constraints with which to refine hydrophobic solvation theories and optimize the set of interatomic potential models used in MD simulations of water and noble gases. 
    more » « less