Hydrogen bonds and van der waals forces in ice at ambient and high pressures

Complex-Valued K-Means Clustering of Interpolative Separable Density Fitting Algorithm for Large-Scale Hybrid Functional Enabled Ab Initio Molecular Dynamics Simulations within Plane Waves

K-means clustering, as a classic unsupervised machine learning algorithm, is the key step to select the interpolation sampling points in interpolative separable density fitting (ISDF) decomposition for hybrid functional electronic structure calculations. Real-valued K-means clustering for accelerating the ISDF decomposition has been demonstrated for large-scale hybrid functional enabled ab initio molecular dynamics (hybrid AIMD) simulations within plane-wave basis sets where the Kohn-Sham orbitals are real-valued. However, it is unclear whether such K-means clustering works for complex-valued Kohn-Sham orbitals. Here, we propose an improved weight function defined as the sum of the square modulus of complex-valued Kohn-Sham orbitals in K-means clustering for hybrid AIMD simulations. Numerical results demonstrate that the K-means algorithm with a new weight function yields smoother and more delocalized interpolation sampling points, resulting in smoother energy potential, smaller energy drift, and longer time steps for hybrid AIMD simulations compared to the previous weight function used in the real-valued K-means algorithm. In particular, we find that this improved algorithm can obtain more accurate oxygen-oxygen radial distribution functions in liquid water molecules and a more accurate power spectrum in crystal silicon dioxide compared to the previous K-means algorithm. Finally, we describe a massively parallel implementation of this ISDF decomposition to accelerate large-scale complex-valued hybrid AIMD simulations containing thousands of atoms (2,744 atoms), which can scale up to 5,504 CPU cores on modern supercomputers.

Diffusion Monte Carlo method for barrier heights of multiple proton exchanges and complexation energies in small water, ammonia, and hydrogen fluoride clusters

Proton exchange reactions are of key importance in many processes in water. However, it is nontrivial to achieve reliable barrier heights for multiple proton exchanges and complexation energies in hydrogen-bonded systems theoretically. Performance of the fixed-node diffusion quantum Monte Carlo (FN-DMC) with the single-Slater-Jastrow trial wavefunction on total energies, barrier heights of multiple proton exchanges, and complexation energies of small water, ammonia, and hydrogen fluoride clusters is investigated in this study. Effects of basis sets and those of locality approximation (LA), T-move approximation (T-move), and determinant localization approximation (DLA) schemes in dealing with the nonlocal part of pseudopotentials on FN-DMC results are evaluated. According to our results, diffuse basis functions are important in achieving reliable barrier heights and complexation energies with FN-DMC, although the cardinal number of the basis set is more important than diffuse basis functions on total energies of these systems. Our results also show that the time step bias with DLA and LA is smaller than T-move; however, the time step bias of DMC energies with respect to time steps using the T-move is roughly linear up to 0.06 a.u., while this is not the case with LA and DLA. Barrier heights and complexation energies with FN-DMC using these three schemes are always within chemical accuracy. Taking into account the fact that T-move and DLA are typically more stable than LA, FN-DMC calculations with the T-move or DLA scheme and basis sets containing diffuse basis functions are suggested for barrier heights of multiple proton exchanges and complexation energies of hydrogen-bonded clusters.

Waste to Medicine: Evidence from Computational Studies on the Modulatory Role of Corn Silk on the Therapeutic Targets Implicated in Type 2 Diabetes Mellitus

Type 2 diabetes mellitus (T2DM) is characterized by insulin resistance and/or defective insulin production in the human body. Although the antidiabetic action of corn silk (CS) is well-established, the understanding of the mechanism of action (MoA) behind this potential is lacking. Hence, this study aimed to elucidate the MoA in different samples (raw and three extracts: aqueous, hydro-ethanolic, and ethanolic) as a therapeutic agent for the management of T2DM using metabolomic profiling and computational techniques. Ultra-performance liquid chromatography-mass spectrometry (UP-LCMS), in silico techniques, and density functional theory were used for compound identification and to predict the MoA. A total of 110 out of the 128 identified secondary metabolites passed the Lipinski's rule of five. The Kyoto Encyclopaedia of Genes and Genomes pathway enrichment analysis revealed the cAMP pathway as the hub signaling pathway, in which ADORA1, HCAR2, and GABBR1 were identified as the key target genes implicated in the pathway. Since gallicynoic acid (-48.74 kcal/mol), dodecanedioc acid (-34.53 kcal/mol), and tetradecanedioc acid (-36.80 kcal/mol) interacted well with ADORA1, HCAR2, and GABBR1, respectively, and are thermodynamically stable in their formed compatible complexes, according to the post-molecular dynamics simulation results, they are suggested as potential drug candidates for T2DM therapy via the maintenance of normal glucose homeostasis and pancreatic β-cell function.

H2O2(s) and H2O2·2H2O(s) crystals compared with ices: DFT functional assessment and D3 analysis

The H2O and H2O2 molecules resemble each other in a multitude of ways as has been noted in the literature. Here, we present density functional theory (DFT) calculations for the H2O2(s) and H2O2·2H2O(s) crystals and make selected comparisons with ice polymorphs. The performance of a number of dispersion-corrected density functionals-both self-consistent and a posteriori ones-are assessed, and we give special attention to the D3 correction and its effects. The D3 correction to the lattice energies is large: for H2O2(s) the D3 correction constitutes about 25% of the lattice energy using PBE, much more for RPBE, much less for SCAN, and it primarily arises from non-H-bonded interactions out to about 5 Å.The large D3 corrections to the lattice energies are likely a consequence of several effects: correction for missing dispersion interaction, the ability of D3 to capture and correct various other kinds of limitations built into the underlying DFT functionals, and finally some degree of cell-contraction-induced polarization enhancement. We find that the overall best-performing functionals of the twelve examined are optPBEvdW and RPBE-D3. Comparisons with DFT assessments for ices in the literature show that where the same methods have been used, the assessments largely agree.

Impact of crystalline structure on the digestibility of amylopectin-based starch-lipid complexes

In this study, chain-elongated waxy corn starch (mWCS) was complexed with lauric acid (LA) to produce starch-lipid complexes (mWCS@LA) with a mixture of B- and V-type crystalline structures. Results from in vitro digestion showed that mWCS@LA had higher digestibility than mWCS, and the logarithm of slope plots of mWCS@LA revealed a two-stage digestion pattern, with digestion rate of the first stage (k₁ = 0.038 min^-1) being much higher than that of the following stage (k₂ = 0.0116 min^-1). The complexation between the long branch chains of mWCS and LA formed amylopectin-based V-type crystallites that were rapidly hydrolyzed during the first stage. The digesta isolated from the second stage of digestion had a B-type crystallinity of 52.6 %, and starch chains with degree of polymerization of 24-28 mainly contributed to the formation of the B-type crystalline structure. The results from the present study reveal that the B-type crystallites were more resistant to amylolytic hydrolysis than the amylopectin-based V-type crystallites.

Extending density functional theory with near chemical accuracy beyond pure water

Density functional simulations of condensed phase water are typically inaccurate, due to the inaccuracies of approximate functionals. A recent breakthrough showed that the SCAN approximation can yield chemical accuracy for pure water in all its phases, but only when its density is corrected. This is a crucial step toward first-principles biosimulations. However, weak dispersion forces are ubiquitous and play a key role in noncovalent interactions among biomolecules, but are not included in the new approach. Moreover, naïve inclusion of dispersion in HF-SCAN ruins its high accuracy for pure water. Here we show that systematic application of the principles of density-corrected DFT yields a functional (HF-r²SCAN-DC4) which recovers and not only improves over HF-SCAN for pure water, but also captures vital noncovalent interactions in biomolecules, making it suitable for simulations of solutions.

Capturing Individual Hydrogen Bond Strengths in Ices via Periodic Local Vibrational Mode Theory: Beyond the Lattice Energy Picture

Local stretching force constants derived from periodic local vibrational modes at the vdW-DF2 density functional level have been employed to quantify the intrinsic hydrogen bond strength of 16 ice polymorphs, ices Ih, II, III, IV, V, VI, VII, VIII, IX, XI, XII, XIII, XIV, XV, XVII, and XIX, that are stable under ambient to elevated pressures. Based on this characterization on 1820 hydrogen bonds, relationships between local stretching force constants and structural parameters such as hydrogen bond length and angle were identified. Moreover, different bond strength distributions, from uniform to inhomogeneous, were observed for the 16 ices and could be explained in relation to different local structural elements within ices, that is, rings, that consist of different hydrogen bond types. In addition, criteria for the classification of hydrogen bonds as strong, intermediate, and weak were introduced. The latter was used to explore a different dimension of the water-ice phase diagram. These findings will provide important guidelines for assessing the credibility of new ice structures.

Superionicity, disorder, and bandgap closure in dense hydrogen chloride

Hydrogen bond networks play a crucial role in biomolecules and molecular materials such as ices. How these networks react to pressure directs their properties at extreme conditions. We have studied one of the simplest hydrogen bond formers, hydrogen chloride, from crystallization to metallization, covering a pressure range of more than 2.5 million atmospheres. Following hydrogen bond symmetrization, we identify a previously unknown phase by the appearance of new Raman modes and changes to x-ray diffraction patterns that contradict previous predictions. On further compression, a broad Raman band supersedes the well-defined excitations of phase V, despite retaining a crystalline chlorine substructure. We propose that this mode has its origin in proton (H⁺) mobility and disorder. Above 100 GPa, the optical bandgap closes linearly with extrapolated metallization at 240(10) GPa. Our findings suggest that proton dynamics can drive changes in these networks even at very high densities.

Energies, structures, and harmonic frequencies of small water clusters from the direct random phase approximation

The binding energies, structures, and vibrational frequencies of water clusters up to 20 molecules are computed at the direct random phase approximation (RPA) level of theory and compared to theoretical benchmarks. Binding energies of the WATER27 set, which includes neutral and positively and negatively charged clusters, are predicted to be too low in the complete basis set limit by an average of 7 kcal/mol (9%) and are worse than the results from the best density functional theory methods or from the Møller-Plesset theory. The RPA shows significant basis set size dependence for binding energies. The order of the relative energies of the water hexamer and dodecamer isomers is predicted correctly by the RPA. The mean absolute deviation for angles and distances for neutral clusters up to the water hexamer are 0.2° and 0.6 pm, respectively, using quintuple-ζ basis sets. The relative energetic order of the hexamer isomers is preserved upon optimization. Vibrational frequencies for these systems are underestimated by several tens of wavenumbers for large basis sets, and deviations increase with the basis set size. Overall, the direct RPA method yields accurate structural parameters but systematically underestimates binding energies and shows strong basis set size dependence.

Interactions between large molecules pose a puzzle for reference quantum mechanical methods

Quantum-mechanical methods are used for understanding molecular interactions throughout the natural sciences. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are state-of-the-art trusted wavefunction methods that have been shown to yield accurate interaction energies for small organic molecules. These methods provide valuable reference information for widely-used semi-empirical and machine learning potentials, especially where experimental information is scarce. However, agreement for systems beyond small molecules is a crucial remaining milestone for cementing the benchmark accuracy of these methods. We show that CCSD(T) and DMC interaction energies are not consistent for a set of polarizable supramolecules. Whilst there is agreement for some of the complexes, in a few key systems disagreements of up to 8 kcal mol^-1 remain. These findings thus indicate that more caution is required when aiming at reproducible non-covalent interactions between extended molecules.

Ammonium fluoride's analogy to ice: Possibilities and limitations

Ammonium fluoride, NH₄F, is often seen as an analog to ice, with several of its solid phases closely resembling known ice phases. While its ionic and hydrogen-ordered nature puts topological constraints on the ice-like network structures it can form, it is not clear what consequences these constraints have for NH₄F compound formation and evolution. Here, we explore computationally the reach and eventual limits of the ice analogy for ammonium fluoride. By combining data mining of known and hypothetical ice networks with crystal structure prediction and density functional calculations, we explore the high-pressure phase diagram of NH₄F and host-guest compounds of its hydrides. Pure NH₄F departs from ice-like behavior above 80 GPa with the emergence of close-packed ionic structures. The predicted stability of NH₄F hydrides shows that NH₄F can act as a host to small guest species, albeit in a topologically severely constraint configuration space. Finally, we explore the binary NH₃-HF chemical space, where we find candidate structures for several unsolved polyfluoride phases; among them is the chemical analog to H₂O₂ dihydrate.

Liquid water contains the building blocks of diverse ice phases

Water molecules can arrange into a liquid with complex hydrogen-bond networks and at least 17 experimentally confirmed ice phases with enormous structural diversity. It remains a puzzle how or whether this multitude of arrangements in different phases of water are related. Here we investigate the structural similarities between liquid water and a comprehensive set of 54 ice phases in simulations, by directly comparing their local environments using general atomic descriptors, and also by demonstrating that a machine-learning potential trained on liquid water alone can predict the densities, lattice energies, and vibrational properties of the ices. The finding that the local environments characterising the different ice phases are found in water sheds light on the phase behavior of water, and rationalizes the transferability of water models between different phases.

Molecular understanding of water is challenging due to the structural complexity of liquid water and the large number of ice phases. Here the authors use a machine-learning potential trained on liquid water to demonstrate the structural similarity of liquid water and that of 54 real and hypothetical ice phases.

Hydration of NH_{4}^{+} in Water: Bifurcated Hydrogen Bonding Structures and Fast Rotational Dynamics

Understanding the hydration and diffusion of ions in water at the molecular level is a topic of widespread importance. The ammonium ion (NH_{4}^{+}) is an exemplar system that has received attention for decades because of its complex hydration structure and relevance in industry. Here we report a study of the hydration and the rotational diffusion of NH_{4}^{+} in water using ab initio molecular dynamics simulations and quantum Monte Carlo calculations. We find that the hydration structure of NH_{4}^{+} features bifurcated hydrogen bonds, which leads to a rotational mechanism involving the simultaneous switching of a pair of bifurcated hydrogen bonds. The proposed hydration structure and rotational mechanism are supported by existing experimental measurements, and they also help to rationalize the measured fast rotation of NH_{4}^{+} in water. This study highlights how subtle changes in the electronic structure of hydrogen bonds impacts the hydration structure, which consequently affects the dynamics of ions and molecules in hydrogen bonded systems.

Squeezing Oil into Water under Pressure: Inverting the Hydrophobic Effect

The molecular structure of dense homogeneous fluid water-methane mixtures has been determined for the first time using high-pressure neutron-scattering techniques at 1.7 and 2.2 GPa. A mixed state with a fully H-bonded water network is revealed. The hydration shell of the methane molecules is, however, revealed to be pressure-dependent with an increase in the water coordination between 1.7 and 2.2 GPa. In parallel, ab initio molecular dynamics simulations have been performed to provide insight into the microscopic mechanisms associated with the phenomenon of mixing. These calculations reproduce the observed phase change from phase separation to mixing with increasing pressure. The calculations also reproduce the experimentally observed structural properties. Unexpectedly, the simulations show mixing is accompanied by a subtle enhancement of the polarization of methane. Our results highlight the key role played by fine electronic effects on miscibility and the need to readjust our fundamental understanding of hydrophobicity to account for these.

Self-interaction error overbinds water clusters but cancels in structural energy differences

We gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi-Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction-corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/[Formula: see text] from about 42 meV/[Formula: see text] for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction-corrected SCAN. These results indicate that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

Experimental and Simulation Insights into Local Structure and Luminescence Evolution in Eu3+-Doped Nanocrystals under High Pressure

Tremendous effort has been devoted to tailoring structure-correlated properties, especially for the luminescence of lanthanide nanocrystals (NCs). High pressure has been demonstrated as a decent way to tune the performance of lanthanide NCs; however, little attention has been paid to the local structure evolution accompanied by extreme compression and its effect on luminescence. Here, we tailor the local structure around lanthanide ions with pressure in β-NaGdF₄ NCs, in which Eu³⁺ ions were doped as optical probes for local structure for the sensitive electric dipole transition. As the pressure increases, the intensity ratio of the ⁵D₀ → ⁷F₂ to ⁵D₀ → ⁷F₁ transition decreases monotonically from 2.04 to 0.81, implying a higher local symmetry around Eu³⁺ ions from compression. In situ X-ray diffraction demonstrates that the sample maintains the hexagonal structure up to 33.5 GPa, and density functional theory calculations reveal the tendency of the local structure to vary under high pressure.

Phase Transition of Ice at High Pressures and Low Temperatures

The behavior of ice under extreme conditions undergoes the change of intermolecular binding patterns and leads to the structural phase transitions, which are needed for modeling the convection and internal structure of the giant planets and moons of the solar system as well as H2O-rich exoplanets. Such extreme conditions limit the structural explorations in laboratory but open a door for the theoretical study. The ice phases IX and XIII are located in the high pressure and low temperature region of the phase diagram. However, to the best of our knowledge, the phase transition boundary between these two phases is still not clear. In this work, based on the second-order Møller–Plesset perturbation (MP2) theory, we theoretically investigate the ice phases IX and XIII and predict their structures, vibrational spectra and Gibbs free energies at various extreme conditions, and for the first time confirm that the phase transition from ice IX to XIII can occur around 0.30 GPa and 154 K. The proposed work, taking into account the many-body electrostatic effect and the dispersion interactions from the first principles, opens up the possibility of completing the ice phase diagram and provides an efficient method to explore new phases of molecular crystals.

Exotic Hydrogen Bonding in Compressed Ammonia Hydrides

Hydrogen-rich compounds attract significant fundamental and practical interest for their ability to accommodate diverse hydrogen bonding patterns and their promise as superior energy storage materials. Here, we report on an intriguing discovery of exotic hydrogen bonding in compressed ammonia hydrides and identify two novel ionic phases in an unusual stoichiometry NH₇. The first is a hexagonal R3̅ m phase containing NH₃-H⁺-NH₃, H^-, and H₂ structural units stabilized above 25 GPa. The exotic NH₃-H⁺-NH₃ unit comprises two NH₃ molecules bound to a proton donated from a H₂ molecule. Above 60 GPa, the structure transforms to a tetragonal P4₁2₁2 phase comprising NH₄⁺, H^-, and H₂ units. At elevated temperatures, fascinating superionic phases of NH₇ with part-solid and part-liquid structural forms are identified. The present findings advance fundamental knowledge about ammonia hydrides at high pressure with broad implications for studying planetary interiors and superior hydrogen storage materials.

Stability of Hydrogen Hydrates from Second-Order Møller-Plesset Perturbation Theory

The formation of gas hydrates and clathrates critically depends on the interaction between the host water network and the guest gas species. Density functional calculations can struggle to quantitatively capture these dispersion-type interactions. Here, we report wave function-based calculations on hydrogen hydrates that combine periodic Hartree-Fock with a localized treatment of electronic correlation. We show that local second-order Møller-Plesset perturbation theory (LMP2) reproduces the stability of the different filled-ice-like hydrates in excellent agreement with experimental data. In contrast to various dispersion-corrected density functional theory implementations, LMP2 correctly identifies the pressures needed to stabilize the C₀, C₁, and C₂ hydrates and does not find a spurious region of stability for an ice-I_h-based dihydrate. Our results suggest that LMP2 or similar approaches can provide quantitative insights into the mechanisms of formation and eventual decomposition of molecular host-guest compounds.

Filling Ices with Helium and the Formation of Helium Clathrate Hydrate

We have formed the long-sought He-clathrate. This was achieved by refilling helium into ice XVI, opening a new synthesis route for exotic forms of clathrate hydrates. The process was followed by neutron diffraction; structures and cage fillings were established. The stabilizing attractive van der Waals interactions are enhanced by multiple cage fillings with theoretically up to four helium atoms per large cage and up to one per small cage; He-clathrate hydrates can be considered as a solid-state equivalent of the clustering of small apolar entities dissolved in the liquid state of water. Unlike most other guests, helium easily enters and leaves the water cages at temperatures well below 100 K, hampering applications as a gas storage material. Despite the weak dispersive interactions, the inclusion of helium has a very significant effect on lattice constants; this is also established for helium inclusion in ice Ih and suggests that lattice parameters are arguably the most sensitive measure to gauge dispersive water-gas interactions.

Fast and accurate quantum Monte Carlo for molecular crystals

Computational approaches based on the fundamental laws of quantum mechanics are now integral to almost all materials design initiatives in academia and industry. If computational materials science is genuinely going to deliver on its promises, then an electronic structure method with consistently high accuracy is urgently needed. We show that, thanks to recent algorithmic advances and the strategy developed in our manuscript, quantum Monte Carlo yields extremely accurate predictions for the lattice energies of materials at a surprisingly modest computational cost. It is thus no longer a technique that requires a world-leading computational facility to obtain meaningful results. While we focus on molecular crystals, the significance of our findings extends to all classes of materials.

Computer simulation plays a central role in modern-day materials science. The utility of a given computational approach depends largely on the balance it provides between accuracy and computational cost. Molecular crystals are a class of materials of great technological importance which are challenging for even the most sophisticated ab initio electronic structure theories to accurately describe. This is partly because they are held together by a balance of weak intermolecular forces but also because the primitive cells of molecular crystals are often substantially larger than those of atomic solids. Here, we demonstrate that diffusion quantum Monte Carlo (DMC) delivers subchemical accuracy for a diverse set of molecular crystals at a surprisingly moderate computational cost. As such, we anticipate that DMC can play an important role in understanding and predicting the properties of a large number of molecular crystals, including those built from relatively large molecules which are far beyond reach of other high-accuracy methods.

Kinetically Controlled Two-Step Amorphization and Amorphous-Amorphous Transition in Ice

We report the results of in situ structural characterization of the amorphization of crystalline ice Ih under compression and the relaxation of high-density amorphous (HDA) ice under decompression at temperatures between 96 and 160 K by synchrotron x-ray diffraction. The results show that ice Ih transforms to an intermediate crystalline phase at 100 K prior to complete amorphization, which is supported by molecular dynamics calculations. The phase transition pathways show clear temperature dependence: direct amorphization without an intermediate phase is observed at 133 K, while at 145 K a direct Ih-to-IX transformation is observed; decompression of HDA shows a transition to low-density amorphous ice at 96 K and ∼1  Pa, to ice Ic at 135 K and to ice IX at 145 K. These observations show that the amorphization of compressed ice Ih and the recrystallization of decompressed HDA are strongly dependent on temperature and controlled by kinetic barriers. Pressure-induced amorphous ice is an intermediate state in the phase transition from the connected H-bond water network in low pressure ices to the independent and interpenetrating H-bond network of high-pressure ices.

Performance of a nonempirical density functional on molecules and hydrogen-bonded complexes

Recently, Tao and Mo derived a meta-generalized gradient approximation functional based on a model exchange-correlation hole. In this work, the performance of this functional is assessed on standard test sets, using the 6-311++G(3df,3pd) basis set. These test sets include 223 G3/99 enthalpies of formation, 99 atomization energies, 76 barrier heights, 58 electron affinities, 8 proton affinities, 96 bond lengths, 82 harmonic vibrational frequencies, 10 hydrogen-bonded molecular complexes, and 22 atomic excitation energies. Our calculations show that the Tao-Mo functional can achieve high accuracy for most properties considered, relative to the local spin-density approximation, Perdew-Burke-Ernzerhof, and Tao-Perdew-Staroverov-Scuseria functionals. In particular, it yields the best accuracy for proton affinities, harmonic vibrational frequencies, hydrogen-bond dissociation energies and bond lengths, and atomic excitation energies.

Performance of a nonempirical exchange functional from density matrix expansion: comparative study with different correlations

Recently, Tao and Mo proposed a meta-generalized gradient approximation for the exchange–correlation energy with remarkable accuracy for molecules, solids, and surfaces.

The role of van der Waals and exchange interactions in high-pressure solid hydrogen

Our study of the van der Waals interactions in solid molecular hydrogen structures indicates two candidates for phase III.

Modelling pH and potential in dynamic structures of the water/Pt(111) interface on the atomic scale

Modelling liquid structures averages of water in the interface with Pt(111) as grand canonical averages, that are functions of pH and electrode potential, using work functions as the absolute potential scale.

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.