On the room-temperature phase diagram of high pressure hydrogen: an ab initio molecular dynamics perspective and a diffusion Monte Carlo study

High-pressure structures of solid hydrogen: Insights from ab initio molecular dynamics simulations

Understanding the structural behavior of solid hydrogen under high pressures is crucial for uncovering its unique properties and potential applications. In this study, starting from the phase I of solid hydrogen-free-rotator hcp structure, we conduct extensive ab initio molecular dynamics calculations to simulate the cooling, heating, and equilibrium processes within a pressure range of 80-260 GPa. Without relying on any structure previously predicted, we identify the high-pressure phase structures of solid hydrogen as P21/c for phase II, P6522 for phase III, and BG1BG2BG3 six-layer structure for phase IV, which are different from those proposed previously using the structure-search method. The reasonability of these structures are validated by Raman spectra and x-ray diffraction patterns by comparison with the experimental results. Our results actually show pronounced changes in the c/a ratio between phases I, III, and IV, which hold no brief for the experimental interpretation of an isostructural hcp transformations for phases I-III-IV.

Toward Accurate Adsorption Energetics on Clay Surfaces

Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Hydrogen segregation and its roles in structural stability and metallization: silane under pressure

We present results from first-principles calculations on silane (SiH4) under pressure. We find that a three dimensional P-3 structure becomes the most stable phase above 241 GPa. A prominent structural feature, which separates the P-3 structure from previously observed/predicted SiH4 structures, is that a fraction of hydrogen leaves the Si-H bonding environment and forms segregated H2 units. The H2 units are sparsely populated in the system and intercalated with a polymeric Si-H framework. Calculations of enthalpy of formation suggest that the P-3 structure is against the decomposition into Si-H binaries and/or the elemental crystals. Structural stability of the P-3 structure is attributed to the electron-deficient multicenter Si-H-Si interactions when neighboring silicon atoms are linked together through a common hydrogen atom. Within the multicenter bonds, electrons are delocalized and this leads to a metallic state, possibly also a superconducting state, for SiH4. An interesting outcome of the present study is that the enthalpy sum of SiH4 (P-3 structure) and Si (fcc structure) appears to be lower than the enthalpy of disilane (Si2H6) between 200 and 300 GPa (for all previously predicted crystalline forms of Si2H6), which calls for a revisit of the stability of Si2H6 under high pressure.

Predicting crystal structures and properties of matter under extreme conditions via quantum mechanics: the pressure is on

The role of quantum mechanical calculations in understanding and predicting the behavior of matter at extreme pressures is discussed in this feature contribution.