An official website of the United States government
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
The ab initio potential energy curves of atom pairs and transport properties of high-temperature vapors of Cu and Si and their mixtures with He, Ar, and Xe gases
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
- Award ID(s):
- 2148653
- PAR ID:
- 10448698
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 25
- Issue:
- 6
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 4872 to 4898
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Polyelectrolyte complexes, PECs, are glassy and brittle when dry but may be plasticized with water. Though hydrated PECs contain a high proportion of water, many still exhibit a glass transition in the 0 to 100 oC range. The apparently unique effectiveness of water as a plasticizer of PECs has been an obstacle to further developments in applications and in fundamental studies of PEC properties. In this work it is shown that formamide is an excellent and even superior solvent for plasticizing PECs, substantially decreasing glass transition temperatures relative to those of hydrated PECs when formamide is used as a solvent instead. The affinities of PECs for water and formamide, indicated by the (exothermic) enthalpies of solvent swelling of dry PECs, are comparable. Ion transport dynamics revealed similar lifetimes, about 1 ns, of charge pairs within a PEC solvated with water compared to formamide, despite the differences in their dielectric constants. Ion transport dynamics, which depend on the mobility of pendant groups, have lower cooperativity than those of the polymer backbone. The use of formamide is a significant experimental variable for reducing the glass transition temperature/viscosity of complexed polyelectrolytes and can turn a solidlike hydrated complex into a fluidlike coacervate.more » « less
-
The results of a combined experimental and computational investigation of the structural evolution of Au 81 Si 19 , Pd 82 Si 18 , and Pd 77 Cu 6 Si 17 metallic glass forming liquids are presented. Electrostatically levitated metallic liquids are prepared, and synchrotron x-ray scattering studies are combined with embedded atom method molecular dynamics simulations to probe the distribution of relevant structural units. Metal–metalloid based metallic glass forming systems are an extremely important class of materials with varied glass forming ability and mechanical processibility. High quality experimental x-ray scattering data are in poor agreement with the data from the molecular dynamics simulations, demonstrating the need for improved interatomic potentials. The first peak in the x-ray static structure factor in Pd 77 Cu 6 Si 17 displays evidence for a Curie–Weiss type behavior but also a peak in the effective Curie temperature. A proposed order parameter distinguishing glass forming ability, [Formula: see text], shows a peak in the effective Curie temperature near a crossover temperature established by the behavior of the viscosity, T A .more » « less
-
The linear viscoelastic response, LVR, of a hydrated polyelectrolyte complex coacervate, PEC, was evaluated over a range of frequencies, temperatures, and salt concentrations. The PEC was a nearly-stoichiometric blend of a quaternary ammonium poly([3-(methacrylamido)propyl]trimethylammonium chloride), PMAPTAC, and poly(2-acrylamido-2-methyl-1-propanesulfonic acid sodium salt), PAMPS, an aliphatic sulfonate, selected because they remain fully charged over the conditions of use. Narrow molecular weight distribution polyelectrolytes were prepared using fractionation techniques. A partially deuterated version of PMAPTAC was incorporated to determine the coil radius of gyration, Rg, within PECs using small angle neutron scattering. Chain dimensions were determined to be Gaussian with a Kuhn length of 2.37 nm, which remained constant from 25 to 65 0C. The LVR for a series of matched molecular weight PECs, mostly above the entanglement threshold, exhibited crossovers of modulus versus frequency classically attributed to the reptation time, relaxation between entanglements, and the relaxation of a Kuhn length of units (the “monomer” time). The scaling for zero shear viscosity, η0, versus chain length N, was η0 ~ N3.1, in agreement with “sticky reptation” theory. The lifetime and activation energy, Ep, of a pair between polyanion and polycation repeat units, Pol+Pol-, were determined from diffusion coefficients of salt ions within the PEC. The activation energy for LVR of salt-free PECs was 2Ep, showing that the key mechanism limiting the dynamics of undoped PECs is pair exchange. An FTIR technique was used to distinguish whether SCN- acts as a counterion or a co-ion within PECs. Doping of PECs with NaSCN breaks Pol+Pol- pairing efficiently, which decreases effective crosslinking and decreases viscosity. An equation was derived that quantitatively predicts this effect.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/d2cp04981c
