We used first-principle approaches to calculate the equilibrium isotopic fractionation factors of potassium (K) and rubidium (Rb) in a variety of minerals of geological relevance (orthoclase, albite, muscovite, illite, sylvite, and phlogopite). We also used molecular dynamics simulations to calculate the equilibrium isotopic fractionation factors of K in water. Our results indicate that K and Rb form bonds of similar strengths and that the ratio between the equilibrium fractionations of K and Rb is approximately 3–4. Under low-temperature conditions relevant to weathering of continents or alteration of seafloor basalts (∼25 °C), the K isotopic fractionation between solvated K+ and illite (a proxy for K-bearing clays) is +0.24‰, exceeding the current analytical precision, so equilibrium isotopic fractionation can induce measurable isotopic fractionations for this system at low temperature. These findings, however, cannot easily explain why the δ41K value of seawater is shifted by +0.6‰ relative to igneous rocks. Our results indicate that part of the observed fractionation is most likely due to kinetic effects. The narrow range of mean force constants for K and Rb in silicate minerals suggests that phase equilibrium is unlikely to create large K and Rb isotopic fractionations at magmatic temperatures (at least in silicate systems). Kinetic effects associated with diffusion can, however, produce large K and Rb isotopic fractionations in igneous rocks.

The phase diagram of water at extreme conditions plays a critical role in Earth and planetary science, yet remains poorly understood. Here we report a first-principles investigation of the liquid at high temperature, between 11 GPa and 20 GPa-a region where numerous controversial results have been reported over the past three decades. Our results are consistent with the recent estimates of the water melting line below 1,000 K and show that on the 1,000-K isotherm the liquid is rapidly dissociating and recombining through a bimolecular mechanism. We found that short-lived ionic species act as charge carriers, giving rise to an ionic conductivity that at 11 GPa and 20 GPa is six and seven orders of magnitude larger, respectively, than at ambient conditions. Conductivity calculations were performed entirely from first principles, with no a priori assumptions on the nature of charge carriers. Despite frequent dissociative events, we observed that hydrogen bonding persists at high pressure, up to at least 20 GPa. Our computed Raman spectra, which are in excellent agreement with experiment, show no distinctive signatures of the hydronium and hydroxide ions present in our simulations. Instead, we found that infrared spectra are sensitive probes of molecular dissociation, exhibiting a broad band below the OH stretching mode ascribable to vibrations of complex ions.

Perturbations to water, both by ions and confining media, have been the focus of numerous experimental and theoretical studies. Yet, several open questions remain, including the extent to which such perturbations modify the structural and dielectric properties of the liquid. Here, we present a first-principles molecular dynamics study of alkali cations in water (Li+, Na+, and K+) as well as of water and LiCl and KCl solutions under confinement within carbon nanotubes (CNTs) of small diameter (1.1–1.5 nm). Our simulations support the view that the water structure is only modified locally in the presence of cations. We found that molecular polarizabilities are fingerprints of hydrogen bonding modifications, which occur at most up to the second solvation shell for all cations in bulk water. Under confinement, we found that the overall value of the molecular polarizability of water molecules near the surface is determined by the balance of two effects, which are quantitatively different in CNTs of different radii: the presence of broken hydrogen bonds at the surface leads to a decrease in the polarizabilities of water molecules, while the interaction with the CNT enhances polarizabilities. Interestingly, the reduction in dipole moments of interfacial water molecules under confinement is instead driven only by changes in the water structure and not by interfacial interactions. As expected, confinement effects on water molecular polarizabilities and dipole moments are more pronounced in the case of the 1.1 nm CNT.

The molecular interactions and structural behavior of a previously unexplored clathrate system, hydrogen-loaded β-hydroquinone (β-HQ+H2), were investigated under high pressure with synchrotron X-ray diffraction and Raman/infrared spectroscopies. The β-HQ+H2 system exhibits coupling of two independently rare phenomena: multiple occupancy and negative compressibility. The number of H2 molecules per cavity increases from one to three, causing unit cell volume increase by way of unique crystallographic interstitial guest positioning. We anticipate these occupancy-derived trends may be general to a range of inclusion compounds and may aid the chemical and crystallographic design of both high-occupancy hydrogen storage clathrates and novel, variable-composition materials with tunable mechanical properties.