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: Solubility Equilibrium Isotope Effects of Noble Gases in Water: Theory and Observations
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
Award ID(s):
1923915
PAR ID:
10497332
Author(s) / Creator(s):
; ;
Publisher / Repository:
Journal of Physical Chemistry B
Date Published:
Journal Name:
The Journal of Physical Chemistry B
Volume:
127
Issue:
45
ISSN:
1520-6106
Page Range / eLocation ID:
9802 to 9812
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Physical Chemistry Chemical Physics)
    The potential energy curves (PECs) for the homonuclear He–He, Ar–Ar, Cu–Cu, and Si–Si dimers, as well as heteronuclear Cu–He, Cu–Ar, Cu–Xe, Si–He, Si–Ar, and Si–Xe dimers, are obtained in quantum Monte Carlo (QMC) calculations. It is shown that the QMC method provides the PECs with an accuracy comparable with that of the state-of-the-art coupled cluster singles and doubles with perturbative triples corrections [CCSD(T)] calculations. The QMC data are approximated by the Morse long range (MLR) and (12-6) Lennard-Jones (LJ) potentials. The MLR and LJ potentials are used to calculate the deflection angles in binary collisions of corresponding atom pairs and transport coefficients of Cu and Si vapors and their mixtures with He, Ar, and Xe gases in the range of temperature from 100 K to 10 000 K. It is shown that the use of the LJ potentials introduces significant errors in the transport coefficients of high-temperature vapors and gas mixtures. The mixtures with heavy noble gases demonstrate anomalous behavior when the viscosity and thermal conductivity can be larger than that of the corresponding pure substances. In the mixtures with helium, the thermal diffusion factor is found to be unusually large. The calculated viscosity and diffusivity are used to determine parameters of the variable hard sphere and variable soft sphere molecular models as well as parameters of the power-law approximations for the transport coefficients. The results obtained in the present work include all information required for kinetic or continuum simulations of dilute Cu and Si vapors and their mixtures with He, Ar, and Xe gases. 
    more » « less
  2. ; ; ; ; ; ; (, Nature Communications)
    Abstract Constraining the magnitude of past hydrological change may improve understanding and predictions of future shifts in water availability. Here we demonstrate that water-table depth, a sensitive indicator of hydroclimate, can be quantitatively reconstructed using Kr and Xe isotopes in groundwater. We present the first-ever measurements of these dissolved noble gas isotopes in groundwater at high precision (≤0.005‰ amu−1; 1σ), which reveal depth-proportional signals set by gravitational settling in soil air at the time of recharge. Analyses of California groundwater successfully reproduce modern groundwater levels and indicate a 17.9 ± 1.3 m (±1 SE) decline in water-table depth in Southern California during the last deglaciation. This hydroclimatic transition from the wetter glacial period to more arid Holocene accompanies a surface warming of 6.2 ± 0.6 °C (±1 SE). This new hydroclimate proxy builds upon an existing paleo-temperature application of noble gases and may identify regions prone to future hydrological change. 
    more » « less
  3. ; ; ; (, Atmospheric Measurement Techniques)
    null (Ed.)
    Abstract. The atmospheric He/N2 ratio is expected to increase due to the emission of He associated with fossil fuels and isexpected to also vary in both space and time due to gravitational separationin the stratosphere. These signals may be useful indicators of fossil fuelexploitation and variability in stratospheric circulation, but directmeasurements of He/N2 ratio are lacking on all timescales. Here wepresent a high-precision custom inlet system for mass spectrometers thatcontinuously stabilizes the flow of gas during sample–standard comparisonand removes all non-noble gases from the gas stream. This enablesunprecedented accuracy in measurement of relative changes in the helium molefraction, which can be directly related to the 4He/N2 ratio usingsupplementary measurements of O2/N2, Ar/N2 and CO2.Repeat measurements of the same combination of high-pressure tanks using ourinlet system achieves a He/N2 reproducibility of∼ 10 per meg (i.e., 0.001 %) in 6–8 h analyses. This compares to interannual changesof gravitational enrichment at ∼ 35 km in the midlatitudestratosphere of order 300–400 per meg and an annual tropospheric increasefrom human fossil fuel activity of less than ∼ 30 per meg yr−1 (bounded by previous work on helium isotopes). The gettering andflow-stabilizing inlet may also be used for the analysis of other noble-gasisotopes and could resolve previously unobserved seasonal cycles inKr/N2 and Xe/N2. 
    more » « less
  4. ; ; ; ; ; (, Global Biogeochemical Cycles)
    Abstract The noble gas signature of incoming Pacific Bottom Water (PBW), when compared to North Atlantic Deep Water, indicates the addition of 450 ± 70 GT a−1glacial melt water to form AABW and subsequently PBW. The downstream evolution of this signature between the southern (20°S to equator) and northern (25°–45°N) bottom waters indicates a decrease in sea level pressure around Antarctica over the past two millennia. Vertical profiles of noble gases in the deep Pacific show exponential relationships with depth with scale heights identical to temperature and salinity. Unlike the other noble gases, helium isotopes show evidence of mid‐depth injection of non‐atmospheric helium. Using observed deviations from exponential behavior, we quantify its magnitude and isotope ratio. There is a clear latitude trend in the isotope ratio of this added helium that decreases from a high exceeding 9 RA(atmospheric3He/4He ratio) in the south to around 8 RAnear the equator. North of 30–40°N, it systematically decreases northward to a low of ∼2 RAnorth of 50°N. This decline results from a combination of northward decline in seafloor spreading, release of radiogenic helium from increased sediment thickness, and the possible emission of radiogenic helium through cold seeps along the Alaskan and North American margins. Finally, we derive an improved method of computing the excess helium isotope concentrations and that the distributions of bottom water3HeXS/4HeXSare consistent with what is known about bottom water flow patterns and the input of low3He/4He sedimentary helium. 
    more » « less
  5. ; ; ; ; ; ; (, The Journal of Chemical Physics)
    We investigate attachment of slow electrons (0–10 eV) to naphthalene (Np) clusters in a crossed beam experiment. Supersonic expansions under different conditions using different buffer gases generate the clusters: in He, Ne, and low pressure Ar, neat (Np)N clusters are formed, while we also observe mixed clusters of naphthalene with rare-gas atoms in co-expansion with Ar above 0.5 bar and with Kr. Negatively charged (Np)n− and Rgm(Np)n− (Rg = Ar, Kr) clusters are analyzed by mass spectrometry, and electron energy dependent ion yields are measured. We show that the smallest stable naphthalene complex with an excess electron, the dimer (Np)2− anion, cannot be formed in a binary collision of a free electron with (Np)2 dimer, nor with (Np)3 trimer. Evaporation of a weakly bound Ar atom(s) from a mixed ArM(Np)2 cluster following electron attachment leads to the dimer (Np)2− anion. Larger (Np)n−, n > 3, transient cluster anions decay via evaporation of an Np unit on a timescale of tens of microseconds. The self-scavenging process opens around 6 eV, where a naphthalene unit is electronically excited by the incoming electron, which is slowed down and trapped. However, the transient negative ion is efficiently stabilized only in the mixed clusters, from which Ar atom(s) can be evaporated. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email