Materials in Extreme Environments

We are studying carbon and water bearing solids and fluids under pressure, to improve our understanding of the physical and chemical behavior of carbon and water at extreme conditions, as found in the deep interiors of Earth and other planets. We are also investigating how tweaking and modifying chemical bonds with temperature and pressure may lead to the discovery of new materials and novel properties. In particular, applying pressure to materials may lead to the synthesis of new structures, recoverable once pressure is released.

High Pressure Talk summarizing the work of our group in the field of materials under extreme conditions, November 21st, 2022: [video here]

Salty water under pressure

The investigation of salts in water at extreme conditions is crucial to understanding the properties of aqueous fluids in the Earth. We conducted first principles (FP) and classical molecular dynamics simulations of NaCl in the dilute limit, at temperatures and pressures relevant to the Earth’s upper mantle. Our first principles simulations were performed by coupling SSAGES and Qbox. Similar to ambient conditions, we observe two metastable states of the salt: the contact (CIP) and the solvent-shared ion-pair (SIP), which are entropically and enthalpically favored, respectively. We find that the free energy barrier between the CIP and SIP minima increases at extreme conditions, and that the stability of the CIP is enhanced in FP simulations, consistent with the decrease of the dielectric constant of water. The minimum free energy path between the CIP and SIP becomes smoother at high pressure, and the relative stability of the two configurations is affected by water self-dissociation, which can only be described properly by FP simulations. We also computed photoelectron spectra and thermal conductivity.

CO2 in water at extreme conditions

We used first-principles molecular dynamics simulations to study carbonates and carbon dioxide dissolved in water at pressures (P) and temperatures (T) approximating the conditions of the Earth’s upper mantle. Contrary to popular geochemical models assuming that molecular CO2(aq) is the major carbon species present in water under deep earth conditions, we found that at 11 GPa and 1000 K carbon exists almost entirely in the forms of solvated carbonate (CO3)2- and bicarbonate (HCO3)- ions, and that even carbonic acid (H2CO3(aq)) is more abundant than CO2(aq). Furthermore, our simulations revealed that ion pairing between Na+ and (HCO3)-/(CO3)2- is greatly affected by P-T conditions, decreasing with increasing pressure at 800~1000 K. Our results suggest that in the Earth’s upper mantle, water-rich geo-fluids transport a majority of carbon in the form of rapidly interconverting (CO3)2- and (HCO3)- ions, not solvated CO2(aq) molecules. We are also developing a strategy, based on first principles simulations, to determine ratios of Raman scattering cross-sections of aqueous species under extreme conditions, thus providing a key quantity that can be used, in conjunction with Raman measurements, to predict chemical speciation in aqueous fluids, including carbonates.

Vibrational spectroscopy and ionic conductivity of water under pressure

We carried out first-principles simulations to study the diffusivity, vibrational properties, and conductivity of water in a controversial region of its phase diagram (1,000 K, 10-20 GPa). Our results provide insight into water's dissociation mechanism, the origin of its large ionic conductivity, and the spectroscopic signatures of ionic species. We predicted Raman and infrared spectra at several conditions, which will serve as a guide for future experiments.

Molecular anvils for mechanochemistry under pressure

By combining experiments (carried out by a team at Stanford and Slac) and computations, we designed molecular anvils and demonstrated hydrostatic-pressure-driven redox reactions in metal-organic chalcogenides incorporating molecular elements with heterogeneous compressibility: we showed that bending of bond angles or shearing of adjacent chains activates the metal-chalcogen bonds leading to formation of elemental metal. Our results reveal an unexplored mechanism and enable new possibilities for high-specificity mechanosynthesis.

Tetrahedrally coordinated carbonates in Earth’s lower mantle

Only a small fraction of our planet’s total carbon budget is found at the surface. In fact, Earth’s mantle is thought to be the largest carbon reservoir. Carbonates, and in particular ferromagnesite ((Mg,Fe)CO3), are likely candidates for deep-Earth carbon storage and therefore play a key role in the deep carbon cycle. We recently reported unequivocal evidence of tetrahedrally coordinated carbon in high pressure carbonates, obtained by combined experimental and theoretical studies. In particular, we identified a unique vibrational signature present only in the high-pressure phase of magnesite, and thus a new carbon-oxygen bond formed under pressure.[DCO link]

The refractive index and electronic gap of water and ice increase with increasing pressure

Experimentally, it is not yet possible to measure absorption processes taking place in water and ice in diamond anvil cells. The band gap of diamond is smaller than that of water and ice, at least up to 30 GPa, and at high temperatures water becomes corrosive. We used ab initio molecular dynamics simulations and electronic structure calculations to show that both the refractive index and the electronic gap of water and ice increase with increasing pressure (up to 30 GPa), contrary to previous assumptions and contrary to the results of simple, widely used models. Subtle electronic effects, related to the nature of inter-band transitions and band edge localization under pressure, are responsible for this apparently anomalous behavior. [DCO link]

Water and carbonates under pressure

Water is a major component of fluids in the Earth’s mantle, where its properties are substantially different from those at ambient conditions. At the pressures and temperatures of the mantle, experiments on aqueous fluids are challenging, and several fundamental properties of water are poorly known; e.g., its dielectric constant has not been measured. Using ab initio molecular dynamics, we computed the dielectric constant of water under the conditions of the Earth’s upper mantle, and we predicted the solubility products of carbonate minerals.

Hydrocarbons under pressure

The thermodynamic and kinetic properties of hydrocarbons at high pressures and temperature are important for understanding carbon reservoirs and fluxes in the Earth. This is one of the of foci of the Deep Carbon Observatory project, funded by Sloan. Our combined ab-initio MD and free energy calculations reveal how higher hydrocarbons may be formed from methane at deep Earth pressures and temperature. This work provides a basis for understanding experiments that demonstrated polymerization of methane to form high hydrocarbons, and earlier methane forming reactions under pressure (in the news).

Melting of ice

We combined first-principles molecular dynamics simulations with a two phase approach to determine the melting temperature of the ice-VII phase in the range of 10–50 GPa. Our computed melting temperatures are consistent with existing diamond anvil cell experiments. We find that for pressures between 10 and 40 GPa, ice melts as a molecular solid. For pressures above 45 GPa, there is a sharp increase in the slope of the melting curve due to the presence of molecular dissociation and proton diffusion in the solid before melting. The onset of significant proton diffusion in ice-VII as a function of increasing temperature is found to be gradual and bears many similarities to that of a type-II superionic solid.

Phase transitions in nitrogen

Two independent experiments have reported the presence of a maximum in the nitrogen melting curve, below 90 GPa, however the position and the interpretation of the origin of such maximum differ. By means of ab initio MD simulations and thermodynamic integration techniques, we have determined the phase diagram of nitrogen in the range between 20 and 100 GPa. We find a maximum in the melting line, related to a transformation in the liquid, from molecular N2 to polymeric nitrogen accompanied by an insulator-to-metal transition.