Search

Water confined in nanoscale cavities plays a crucial role in everyday phenomena in geology and biology, as well as technological applications at the water-energy nexus. However, even understanding the basic properties of nano-confined water is extremely challenging for theory, simulations, and experiments. In particular, determining the melting temperature of quasi-one-dimensional ice polymorphs confined in carbon nanotubes has proven to be an exceptionally difficult task, with previous experimental and classical simulations approaches report values ranging from $\sim 180 \text{ K}$ up to $\sim 450 \text{ K}$ at ambient pressure. In this work, we use a machine learning potential that delivers first principles accuracy to study the phase diagram of water for confinement diameters $ 9.5 < d < 12.5 \text{ \AA}$. We find that several distinct ice polymorphs melt in a surprisingly narrow range between $\sim 280 \text{ K}$ and $\sim 310 \text{ K}$, with a melting mechanism that depends on the nanotube diameter. These results shed new light on the melting of ice in one-dimension and have implications for the operating conditions of carbon-based filtration and desalination devices.

Biphasic interfaces are complex but fascinating regimes that display a number of properties distinct from those of the bulk. The CO$_2$-H$_2$O interface, in particular, has been the subject of a number of studies on account of its importance for the carbon life cycle as well as carbon capture and sequestration schemes. Despite this attention, there remain a number of open questions on the nature of the CO$_2$-H$_2$O interface, particularly concerning the interfacial tension and phase behavior of CO$_2$ at the interface. In this paper, we seek to address these ambiguities using ab initio-quality simulations. Harnessing the benefits of machine-learned potentials and enhanced statistical sampling methods, we present an ab initio-level description of the CO$_2$-H$_2$O interface. Interfacial tensions are predicted from 1-500 bar and found to be in close agreement with experiment at the pressures for which experimental data is available. Structural analyses indicate the build-up of an adsorbed, saturated CO$_2$ film forming at low pressure (20 bar) with properties similar to those of the bulk liquid, but preferential perpendicular alignment with respect to the interface. CO$_2$ monolayer build-up coincides with a reduced structuring of water molecules close to the interface. This study highlights the predictive nature of machine-learned potentials for complex macroscopic properties of biphasic interfaces, and the mechanistic insight obtained into carbon dioxide aggregation at the water interface is of high relevance for geoscience, climate research, and materials science., Comment: 14 pages, 10 figures

Many of graphene's remarkable properties are intrinsically linked to its inherent ripples. Defects, whether naturally present or artificially introduced, are known to have a strong impact on the rippling of graphene. However, how defects alter ripple dynamics in two-dimensional (2D) materials in general, and graphene in particular, remains largely unexplored. Here, using machine learning-driven molecular dynamics simulations, we reveal a fundamental connection between defect concentration and ripple dynamics in freestanding graphene sheets. Specifically, we find that at a critical concentration of approximately 0.1%, dynamic rippling undergoes a transition from freely propagating to static ripples. This is in quantitative alignment with the experimentally observed turning point in the non-monotonic scaling of the Young's modulus and emphasises the critical interplay between defects and material dynamics. Our work not only unveils the significant impact of defects on rippling dynamics in graphene but also paves the way to design two-dimensional devices with tailored properties., Comment: 8 pages, 20 pages supplementary information

CO$_2$ capture using carbon-based materials, particularly graphene and graphene-like materials, is a promising strategy to deal with CO$_2$ emissions. However, significant gaps remain in our understanding of the molecular-level interaction between CO$_2$ molecules and graphene, particularly, in terms of chemical bonding and electron transfer. In this work, we employ random structure search and density functional theory to understand the adsorption of CO$_2$ molecules on Ca, Sr, Na, K, and Ti decorated graphene surfaces. Compared to the pristine material, we observe enhanced CO$_2$ adsorption on the decorated graphene surfaces. Particularly on group 2 metals and titanium decorated graphene, CO$_2$ can be strongly chemisorbed as a bent CO$_2$ anion or as an oxalate, depending on the number of CO$_2$ molecules. Electronic structure analysis reveals the adsorption mechanism to involve an ionic charge transfer from the metal adatom to the adsorbed CO$_2$. Overall, this study suggests that reducing CO$_2$ to oxalate on group 2 metals and titanium metal-decorated graphene surfaces is a potential strategy for CO$_2$ storage., Comment: Main manuscript and Supplementary information provided

Calculating sublimation enthalpies of molecular crystal polymorphs is relevant to a wide range of technological applications. However, predicting these quantities at first-principles accuracy -- even with the aid of machine learning potentials -- is a challenge that requires sub-kJ/mol accuracy in the potential energy surface and finite-temperature sampling. We present an accurate and data-efficient protocol based on fine-tuning of the foundational MACE-MP-0 model and showcase its capabilities on sublimation enthalpies and physical properties of ice polymorphs. Our approach requires only a few tens of training structures to achieve sub-kJ/mol accuracy in the sublimation enthalpies and sub 1 % error in densities for polymorphs at finite temperature and pressure. Exploiting this data efficiency, we explore simulations of hexagonal ice at the random phase approximation level of theory at experimental temperatures and pressures, calculating its physical properties, like pair correlation function and density, with good agreement with experiments. Our approach provides a way forward for predicting the stability of molecular crystals at finite thermodynamic conditions with the accuracy of correlated electronic structure theory., Comment: 9 pages; 4 figures

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation.

Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferability from one chemical system to the next. Here, using the state-of-the-art MACE architecture we introduce a single general-purpose ML model, trained on a public database of 150k inorganic crystals, that is capable of running stable molecular dynamics on molecules and materials. We demonstrate the power of the MACE-MP-0 model - and its qualitative and at times quantitative accuracy - on a diverse set problems in the physical sciences, including the properties of solids, liquids, gases, chemical reactions, interfaces and even the dynamics of a small protein. The model can be applied out of the box and as a starting or "foundation model" for any atomistic system of interest and is thus a step towards democratising the revolution of ML force fields by lowering the barriers to entry., Comment: 119 pages, 63 figures, 37MB PDF

The Bernal-Fowler ice rules stipulate that each water molecule in an ice crystal should form four hydrogen bonds. However, in extreme or constrained conditions, the arrangement of water molecules deviates from conventional ice rules, resulting in properties significantly different from bulk water. In this study, we employ machine learning-driven first-principles simulations to identify a new stabilization mechanism in nanoconfined ice phases beyond conventional ice rules. Instead of forming four hydrogen bonds, nanoconfined crystalline ice can form a quasi-one-dimensional hydrogen-bonded structure that exhibits only two hydrogen bonds per water molecule. These structures consist of strongly hydrogen-bonded linear chains of water molecules that zig-zag along one dimension, stabilized by van der Waals interactions that stack these chains along the other dimension. The unusual interplay of hydrogen bonding and van der Waals interactions in nanoconfined ice results in atypical proton behavior such as potential ferroelectric behavior, low dielectric response, and long-range proton dynamics., Comment: main text: 14 pages, 6 figures

The extent of ion pairing in solution is an important phenomenon to rationalise transport and thermodynamic properties of electrolytes. A fundamental measure of this pairing is the potential of mean force (PMF) between solvated ions. The relative stabilities of the paired and solvent shared states in the PMF and the barrier between them are highly sensitive to the underlying potential energy surface. However direct application of accurate electronic structure methods is challenging, since long simulations are required. We develop wavefunction based machine learning potentials with the Random Phase Approximation (RPA) and second order Moller-Plesset (MP2) perturbation theory for the prototypical system of Na and Cl ions in water. We show both methods in agreement, predicting the paired and solvent shared states to have similar energies (within 0.2 kcal/mol). We also provide the same benchmarks for different DFT functionals as well as insight into the PMF based on simple analyses of the interactions in the system.

The structure of oxide-supported metal nanoclusters plays an essential role in their sharply enhanced catalytic activity over bulk metals. Simulations provide the atomic-scale resolution needed to understand these systems. However, the sensitive mix of metal-metal and metal-support interactions which govern their structure puts stringent requirements on the method used, requiring going beyond standard density functional theory (DFT). The method of choice is coupled cluster theory [specifically CCSD(T)], but its computational cost has so far prevented applications. In this work, we showcase two approaches to make CCSD(T) accuracy readily achievable in oxide-supported nanoclusters. First, we leverage the SKZCAM protocol to provide the first benchmarks of oxide-supported nanoclusters, revealing that it is specifically metal-metal interactions that are challenging to capture with DFT. Second, we propose a CCSD(T) correction ($\Delta$CC) to the metal-metal interaction errors in DFT, reaching comparable accuracy to the SKZCAM protocol at significantly lower cost. This forges a path towards studying larger systems at reliable accuracy, which we highlight by identifying a ground state structure in agreement with experiments for Au$_{20}$ on MgO; a challenging system where DFT models have yielded conflicting predictions.

Nucleation in small volumes of water has garnered renewed interest due to the relevance of pore condensation and freezing under conditions of low partial pressures of water, such as in the upper troposphere. Molecular simulations can in principle provide insight on this process at the molecular scale that is challenging to achieve experimentally. However, there are discrepancies in the literature as to whether the rate in confined systems is enhanced or suppressed relative to bulk water at the same temperature and pressure. In this study, we investigate the extent to which the size of the critical nucleus and the rate at which it grows in thin films of water are affected by the thickness of the film. Our results suggest that nucleation remains bulk-like in films that are barely large enough accommodate a critical nucleus. This conclusion seems robust to the presence of physical confining boundaries. We also discuss the difficulties in unambiguously determining homogeneous nucleation rates in nanoscale systems, owing to the challenges in defining the volume. Our results suggest any impact on a film's thickness on the rate is largely inconsequential for present day experiments., Comment: 26 pages, 6 figures in main script, 6 in ESI, submitted to Faraday Discussions for the "Water at interfaces Faraday Discussion" in September 2023

Low-density amorphous ice (LDA) is one of the most common solid materials in the Universe and influences a myriad of cosmological phenomena. In fundamental research, it is the potential key material for understanding the many famous anomalies of liquid water. Despite its significance, the structural nature of LDA remains heavily debated. It is unclear if LDA is a glassy state representing a liquid state or a heavily disordered crystal. Here we show that the experimental structure factor of LDA is best reproduced computationally by small ice grains with amorphous regions at the grain boundaries. Moreover, LDA materials made through different routes crystallise to give ice I with different stacking characteristics. These results suggest LDA is not truly amorphous but instead a partially crystalline material. Going forward, great care is needed in concluding from diffraction data if materials are truly glassy or instead partially crystalline states.

Recent experiments have shown that the out-of-plane dielectric constant of water confined in nanoslits of graphite and hexagonal boron nitride (hBN) is vanishingly small. Despite extensive effort based mainly on classical force-field molecular dynamics (FFMD) approaches, the origin of this phenomenon is under debate. Here we used ab initio molecular dynamics simulations (AIMD) and AIMD-trained machine learning potentials to explore the structure and electronic properties of water confined inside graphene and hBN slits. We found that the reduced dielectric constant arises mainly from the anti-parallel alignment of the water dipoles in the perpendicular direction to the surface in the first two water layers near the solid interface. Although the water molecules retain liquid-like mobility, the interfacial layers exhibit a net ferroelectric ordering and constrained hydrogen-bonding orientations which lead to much reduced polarization fluctuations in the out-of-plane direction at room temperature. Importantly, we show that this effect is independent of the distance between the two confining surfaces of the slit, and it originates in the spontaneous polarization of interfacial water. Our calculations also show no significant variations in the structure and polarization of water near graphene and hBN, despite their different electronic structures. These results are important as they offer new insight into a property of water that plays a critical role in the long-range interactions between surfaces, the electric double-layer formation, ion solvation and transport, as well as biomolecular functioning.

Life on Earth depends upon the dissolution of ionic salts in water, particularly NaCl. However, an atomistic scale understanding of the process remains elusive. Simulations lend themselves conveniently to studying dissolution since they provide the spatio-temporal resolution that can be difficult to obtain experimentally. Nevertheless, the complexity of various inter- and intra-molecular interactions require careful treatment and long time scale simulations, both of which are typically hindered by computational expense. Here, we use advances in machine learning potential methodology to resolve for the first time at an ab initio level of theory the dissolution mechanism of NaCl in water. The picture that emerges is that of a steady ion-wise unwrapping of the crystal preceding its rapid disintegration, reminiscent of crumbling. The onset of crumbling can be explained by a strong increase in the ratio of the surface to volume of the crystal. Overall, dissolution is comprised of a series of highly dynamical microscopic sub-processes, resulting in an inherently stochastic mechanism. These atomistic level insights now pave the way for a general understanding of dissolution mechanisms in other crystals, and the methodology is primed for more complex systems of recent interest such as water/salt interfaces under flow and salt crystals under confinement.

Friction at water-carbon interfaces remains a major puzzle with theories and simulations unable to explain experimental trends in nanoscale waterflow. A recent theoretical framework -- quantum friction (QF)-- proposes to resolve these experimental observations by considering nonadiabatic coupling between dielectric fluctuations in water and graphitic surfaces. Here, using a classical model that enables fine-tuning of the solid's dielectric spectrum, we provide evidence from simulations in general support of QF. In particular, as features in the solid's dielectric spectrum begin to overlap with water's librational and Debye modes, we find an increase in friction in line with that proposed by QF. At the microscopic level, we find that this contribution to friction manifests more distinctly in the dynamics of the solid's charge density than that of water. Our findings suggest that experimental signatures of QF may be more pronounced in the solid's response rather than liquid water's.

Hydrogen is a key player in global strategies to reduce greenhouse gas emissions. In order to make hydrogen a widely-used fuel, we require more efficient methods of storing it than the current standard of pressurized cylinders. An alternative method is to adsorb H$_2$ in a material and avoid the use of high pressures. Among many potential materials, layered materials such as graphene present a practical advantage as they are lightweight. However, graphene and other 2D materials typically bind H$_2$ too weakly to store it at the typical operating conditions of a hydrogen fuel cell. Modifying the material, for example by decorating graphene with adatoms, can strengthen the adsorption energy of H$_2$ molecules, but the underlying mechanisms are still not well understood. In this work, we systematically screen alkali and alkaline earth metal decorated graphene sheets for the adsorption of hydrogen gas from first principles, and focus on the mechanisms of binding. We show that there are three mechanisms of adsorption on metal decorated graphene and each leads to distinctly different hydrogen adsorption structures. The three mechanisms can be described as weak van der Waals physisorption, metal adatom facilitated polarization, and Kubas adsorption. Among these mechanisms, we find that Kubas adsorption is easily perturbed by an external electric field, providing a way to tune H$_2$ adsorption.

Ice is one of the most important and interesting molecular crystals exhibiting a rich and evolving phase diagram. Recent discoveries mean that there are now twenty distinct polymorphs; a structural diversity that arises from a delicate interplay of hydrogen bonding and van der Waals dispersion forces. This wealth of structures provides a stern test of electronic structure theories, with Density Functional Theory (DFT) often not able to accurately characterise the relative energies of the various ice polymorphs. Thanks to recent advances that enable the accurate and efficient treatment of molecular crystals with Diffusion Monte Carlo (DMC), we present here the DMC-ICE13 dataset; a dataset of lattice energies of 13 ice polymorphs. This dataset encompasses the full structural complexity found in the ambient and high-pressure molecular ice polymorphs and when experimental reference energies are available our DMC results deliver sub-chemical accuracy. Using this dataset we then perform an extensive benchmark of a broad range of DFT functionals. Of the functionals considered, we find revPBE-D3 and RSCAN to reproduce reference absolute lattice energies with the smallest error, whilst optB86b-vdW and SCAN+rVV10 have the best performance on the relative lattice energies. Our results suggest that a single functional achieving reliable performance for all phases is still missing, and that care is needed in the selection of the most appropriate functional for the desired application. The insights obtained here may also be relevant to liquid water and other hydrogen bonded and dispersion bonded molecular crystals.

Understanding the molecular and electronic structure of electrolytes at interfaces requires an analysis of the interactions between the electrode surface, the ions, and the solvent environment on equal footing. Here, we tackle this challenge by exploring the initial stages of Cs+ hydration on a Cu(111) surface by combining experiment and theory. Remarkably, we observe "inside out" solvation of Cs ions, i.e, their preferential location at the perimeter of the water clusters on the metal surface. In addition, water-Cs complexes containing multiple Cs+ ions are observed to form at these surfaces. Established models based on maximum ion-water coordination and the double layer notion cannot account for this situation and the complex interplay of microscopic interactions is the key to a fundamental understanding., Comment: 5 pages, 3 figures

The chemical behaviour of single metal atoms largely depends on the local coordination environment, including interactions with the substrate and with the surrounding gas or liquid. However, the key instrumentation for studying such systems at the atomic scale generally requires high vacuum conditions, limiting the degree to which the aforementioned environmental parameters can be investigated. Here we develop a new platform for transmission electron microscopy investigation of single metal atoms in liquids and study the dynamic behaviour of individual platinum atoms on the surface of a single layer MoS2 crystal in water. To achieve the record single atom resolution, we introduce a double liquid cell based on a 2D material heterostructure, which allows us to submerge an atomically thin membrane with liquid on both sides while maintaining the total specimen thickness of only ~ 70 nm. By comparison with an identical specimen imaged under high vacuum conditions, we reveal drastic differences in the single atom resting sites and atomic hopping behaviour, demonstrating that in situ imaging conditions are essential to gain complete understanding of the chemical activity of individual atoms. These findings pave the way for in situ liquid imaging of chemical processes with single atom precision., Comment: Fixed Captions for figures 1 and 2 (did not appear in previous submission)

The structure of liquid water in the proximity of an interface can deviate significantly from that of bulk water, with surface-induced structural perturbations typically converging to bulk values at about ~1 nm from the interface. While these structural changes are well established it is, in contrast, less clear how an interface perturbs the dynamics of water molecules within the liquid. Here, through an extensive set of molecular dynamics simulations of supercooled bulk and interfacial water films and nano-droplets, we observe the formation of persistent, spatially extended dynamical domains in which the average mobility varies as a function of the distance from the interface. This is in stark contrast with the dynamical heterogeneity observed in bulk water, where these domains average out spatially over time. We also find that the dynamical response of water to an interface depends critically on the nature of the interface and on the choice of interface definition. Overall these results reveal a richness in the dynamics of interfacial water that opens up the prospect of tuning the dynamical response of water through specific modifications of the interface structure or confining material.

Experimental measurements have reported ultra-fast and radius-dependent water transport in carbon nanotubes which are absent in boron nitride nanotubes. Despite considerable effort, the origin of this contrasting (and fascinating) behaviour is not understood. Here, with the aid of machine learning-based molecular dynamics simulations that deliver first-principles accuracy, we investigate water transport in single-wall carbon and boron nitride nanotubes. Our simulations reveal a large, radius-dependent hydrodynamic slippage on both materials with water experiencing indeed a $\approx 5$ times lower friction on carbon surfaces compared to boron nitride. Analysis of the diffusion mechanisms across the two materials reveals that the fast water transport on carbon is governed by facile oxygen motion, whereas the higher friction on boron nitride arises from specific hydrogen-nitrogen interactions. This work not only delivers a clear reference of unprecedented accuracy for water flow in single-wall nanotubes, but also provides detailed mechanistic insight into its radius and material dependence for future technological application.

The O vacancy (Ov) formation energy, $E_\textrm{Ov}$, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g. the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute $E_\textrm{Ov}$ accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods such as local coupled cluster with single, double and perturbative triple excitations [CCSD(T)]. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rocksalt MgO, producing the first accurate and well-converged determinations of $E_\textrm{Ov}$ with this method. These reference values are used to benchmark exchange-correlation functionals in DFT and we find that all studied functionals underestimate $E_\textrm{Ov}$, with the average error decreasing along the rungs of Jacob's ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates., Comment: Accepted for publication in Journal of Chemical Physics on 2022-03-10

The aggregation of clay particles in aqueous solution is a ubiquitous everyday process of broad environmental and technological importance. However, it is poorly understood at the all-important atomistic level since it depends on a complex and dynamic interplay of solvent-mediated electrostatic, hydrogen-bonding, and dispersion interactions. With this in mind we have performed an extensive set of classical molecular dynamics simulations (included enhanced sampling simulations) on the interactions between model kaolinite nanoparticles in pure and salty water. Our simulations reveal highly anisotropic behaviour in which the interaction between the nanoparticles varies from attractive to repulsive depending on the relative orientation of the nanoparticles. Detailed analysis reveals that at large separation (>1.5 nm) this interaction is dominated by electrostatic effects whereas at smaller separations the nature of the water hydration structure becomes critical. This study highlights an incredible richness in how clay nanoparticles interact, which should be accounted for in e.g. coarse grained models of clay nanoparticle aggregation.

In principle, the answer to the posed titular question is undoubtedly 'yes.' But in practice, requisite reference data for homogeneous systems have been obtained with a treatment of intermolecular interactions that is different from that typically employed for heterogeneous systems. In this article, we assess the impact of the choice of truncation scheme when comparing water in homogeneous and inhomogeneous environments. Specifically, we use explicit free energy calculations and a simple mean field analysis to demonstrate that using the 'cut-and-shift' version of the Lennard-Jones potential (common to most simple point charge models of water) results in a systematic increase in the melting temperature of ice I$_{\rm h}$. In addition, by drawing an analogy between a change in cutoff and a change in pressure, we use existing literature data for homogeneous ice nucleation at negative pressures to suggest that enhancements due to heterogeneous nucleation may have been overestimated by several orders of magnitude., Comment: Main article: 9 pages, 7 figures. Supplementary Material: 10 pages, 4 figures

Water diffusion across the surfaces of materials is of importance to disparate processes such as water purification, ice formation, and more. Despite reports of rapid water diffusion on surfaces the molecular-level details of such processes remain unclear. Here, with scanning tunneling microscopy, we observe structural rearrangements and diffusion of water trimers at unexpectedly low temperatures (< 10 K) on a copper surface; temperatures at which water monomers or other clusters do not diffuse. Density functional theory calculations reveal a facile trimer diffusion process involving transformations between elongated and almost cyclic conformers in an inchworm-like manner. These subtle intermolecular reorientations maintain an optimal balance of hydrogen-bonding and water-surface interactions throughout the process. This work shows that the diffusion of hydrogen-bonded clusters can occur at exceedingly low temperatures without the need for hydrogen bond breakage or exchange; findings that will influence Ostwald ripening of ice nanoclusters and hydrogen bonded clusters in general.

Gaining fundamental understanding of crystal nucleation processes in metal alloys is crucial for the development and design of high-performance materials with targeted properties. Yet, crystallization is a complex non-equilibrium process and, despite having been studied for decades, the microscopic aspects that govern the crystallization mechanism of a material remain to date elusive. Recent evidence shows that spatial heterogeneity in the supercooled liquid, characterised by extended regions with distinctive mobility and order, may be a key microscopic factor that determines the mechanism of crystal nucleation. These findings have advanced our view of the fundamental nature of crystallization, as most research has assumed that crystal clusters nucleate from random fluctuations in a `homogeneous' liquid. Here, by analysing transition path sampling trajectories, we show that dynamical heterogeneity plays a key role in the mechanism of crystal nucleation in an elemental metal, nickel. Our results demonstrate that crystallization occurs preferentially in regions of low mobility in the supercooled liquid, evidencing the collective dynamical nature of crystal nucleation in Ni. In addition, our results show that low mobility regions form before and spatially overlap with pre-ordered domains that act as precursors to the crystal phase that subsequently emerges. Our results show a clear link between dynamical and structural heterogeneity in the supercooled liquid and its impact on the nucleation mechanism, revealing microscopic descriptors that could pave a novel way to control crystallization processes in metals.

Water in nanoscale cavities is ubiquitous and of central importance to everyday phenomena in geology and biology. However, the properties of nanoscale water can be remarkably different from bulk, as shown e.g., by the anomalously low dielectric constant of water in nanochannels [1], near frictionless water flow [2], or the possible existence of a square ice phase [3]. Such properties suggest that nanoconfined water could be engineered for technological applications in nanouidics [4], electrolyte materials [5], and water desalination [6]. Unfortunately, challenges in experimentally characterising water on the nanoscale and the high cost of first-principles simulations have prevented the molecular level understanding required to control the behavior of water. Here we combine a range of computational approaches to enable a first-principles level investigation of a single layer of water within a graphene-like channel. We find that monolayer water exhibits surprisingly rich and diverse phase behavior that is highly sensitive to temperature and the van der Waals pressure acting within the nanochannel. In addition to multiple molecular phases with melting temperatures varying non-monotonically by over 400 degrees with pressure, we predict a hexatic phase, which is an intermediate between a solid and a liquid, and a superionic phase with a high electrical conductivity exceeding that of battery materials. Notably, this suggests that nanoconfinement could be a promising route towards superionic behavior at easily accessible conditions., Comment: 33 pages, 8 figures

Single-Atom Alloys (SAAs) are a special class of alloy surface catalysts that offer well defined, isolated active sites in a more inert metal host. The dopant sites are generally assumed to have little or no influence on the properties of the host metal, and transport of chemical reactants and products to and from the dopant sites is generally assumed to be facile. Here, by performing density functional theory calculations and surface science experiments, we identify a new physical effect on SAA surfaces, whereby adsorption is destabilised by up to 300 meV on host sites within the perimeter of the reactive dopant site. We identify periodic trends for this behaviour, and demonstrate a zone of exclusion around the reactive sites for a range of adsorbates and combinations of host and dopant metals. Experiments confirm an increased barrier for CO diffusion towards the dopant on a RhCu SAA. This effect offers new possibilities for understanding and designing active sites with tunable energetic landscapes surrounding them., Comment: 25 pages, 5 figures; Accepted in The Journal of Physical Chemistry Letters

Graphene's intrinsically corrugated and wrinkled topology fundamentally influences its electronic, mechanical, and chemical properties. Experimental techniques allow the manipulation of pristine graphene and the controlled production of defects which allows to control the atomic out-of-plane fluctuations and, thus, tune graphene's properties. Here, we perform large scale machine learning-driven molecular dynamics simulations to understand the impact of defects on the structure of graphene. We find that defects cause significantly higher corrugation leading to a strongly wrinkled surface. The magnitude of this structural transformation strongly depends on the defect concentration and specific type of defect. Analysing the atomic neighborhood of the defects reveals that the extent of these morphological changes depends on the preferred geometrical orientation and the interactions between defects. While our work highlights that defects can strongly affect graphene's morphology, it also emphasises the differences between distinct types by linking the global structure to the local environment of the defects.

Simulation techniques based on accurate and efficient representations of potential energy surfaces are urgently needed for the understanding of complex aqueous systems such as solid-liquid interfaces. Here, we present a machine learning framework that enables the efficient development and validation of models for complex aqueous systems. Instead of trying to deliver a globally-optimal machine learning potential, we propose to develop models applicable to specific thermodynamic state points in a simple and user-friendly process. After an initial ab initio simulation, a machine learning potential is constructed with minimum human effort through a data-driven active learning protocol. Such models can afterwards be applied in exhaustive simulations to provide reliable answers for the scientific question at hand. We showcase this methodology on a diverse set of aqueous systems with increasing degrees of complexity. The systems chosen here comprise bulk water with different ions in solution, water on a titanium dioxide surface, as well as water confined in nanotubes and between molybdenum disulfide sheets. Highlighting the accuracy of our approach with respect to the underlying ab initio reference, the resulting models are evaluated in detail with an automated validation protocol that includes structural and dynamical properties and the precision of the force prediction of the models. Finally, we demonstrate the capabilities of our approach for the description of water on the rutile titanium dioxide (110) surface to analyze the structure and mobility of water on this surface. Such machine learning models provide a straightforward and uncomplicated but accurate extension of simulation time and length scales for complex systems.

The fundamental understanding of crystallization, in terms of microscopic kinetic and thermodynamic details, remains a key challenge in the physical sciences. Here, by using in situ graphene liquid cell transmission electron microscopy, we reveal the atomistic mechanism of NaCl crystallization from solutions confined within graphene cells. We find that rock salt NaCl forms with a peculiar hexagonal morphology. We also see the emergence of a transitory graphite-like phase, which may act as an intermediate in a two-step pathway. With the aid of density functional theory calculations, we propose that these observations result from a delicate balance between the substrate-solute interaction and thermodynamics under confinement. Our results highlight the impact of confinement on both the kinetics and thermodynamics of crystallization, offering new insights into heterogeneous crystallization theory and a potential avenue for materials design., Comment: Main: 7 pages, 4 figures; Supplemental Material: 12 pages, 7 figures

Oxygen vacancies are ubiquitous in TiO$_2$ and play key roles in catalysis and magnetism applications.Despite being extensively investigated, the electronic structure of oxygen vacancies in TiO$_2$ remains controversial both experimentally and theoretically.Here we report a study of a neutral oxygen vacancy in TiO$_2$ using state-of-the-art quantum chemical electronic structure methods.We find that the ground state is a color center singlet state in both the rutile and the anatase phase of TiO$_2$. Specifically, embedded CCSD(T) calculations find, for an oxygen vacancy in rutile, that the lowest triplet state energy is 0.6 eV above the singlet state, and in anatase the triplet state energy is higher by 1.4 eV. Our study provides fresh insights on the electronic structure of the oxygen vacancy in TiO$_2$, clarifying earlier controversies and potentially inspiring future studies of defects with correlated wave function theories.

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.

We introduce an interatomic potential for hexagonal boron nitride (hBN) based on the Gaussian approximation potential (GAP) machine learning methodology. The potential is based on a training set of configurations collected from density functional theory (DFT) simulations and is capable of treating bulk and multilayer hBN as well as nanotubes of arbitrary chirality. The developed force field faithfully reproduces the potential energy surface predicted by DFT while improving the efficiency by several orders of magnitude. We test our potential by comparing formation energies, geometrical properties, phonon dispersion spectra and mechanical properties with respect to benchmark DFT calculations and experiments. In addition, we use our model and a recently developed graphene-GAP to analyse and compare thermally and mechanically induced rippling in large scale two-dimensional (2D) hBN and graphene. Both materials show almost identical scaling behaviour with an exponent of $\eta \approx 0.85$ for the height fluctuations agreeing well with the theory of flexible membranes. Based on its lower resistance to bending, however, hBN experiences slightly larger out-of-plane deviations both at zero and finite applied external strain. Upon compression a phase transition from incoherent ripple motion to soliton-ripples is observed for both materials. Our potential is freely available online at [http://www.libatoms.org]., Comment: This article has been accepted for publication in the Journal of Physical Chemistry C

We present an accurate machine learning (ML) model for atomistic simulations of carbon, constructed using the Gaussian approximation potential (GAP) methodology. The potential, named GAP-20, describes the properties of the bulk crystalline and amorphous phases, crystal surfaces and defect structures with an accuracy approaching that of direct ab initio simulation, but at a significantly reduced cost. We combine structural databases for amorphous carbon and graphene, which we extend substantially by adding suitable configurations, for example, for defects in graphene and other nanostructures. The final potential is fitted to reference data computed using the optB88-vdW density functional theory (DFT) functional. Dispersion interactions, which are crucial to describe multilayer carbonaceous materials, are therefore implicitly included. We additionally account for long-range dispersion interactions using a semianalytical two-body term and show that an improved model can be obtained through an optimisation of the many-body smooth overlap of atomic positions (SOAP) descriptor. We rigorously test the potential on lattice parameters, bond lengths, formation energies and phonon dispersions of numerous carbon allotropes. We compare the formation energies of an extensive set of defect structures, surfaces and surface reconstructions to DFT reference calculations. The present work demonstrates the ability to combine, in the same ML model, the previously attained flexibility required for amorphous carbon [Phys. Rev. B, 95, 094203, (2017)] with the high numerical accuracy necessary for crystalline graphene [Phys. Rev. B, 97, 054303, (2018)], thereby providing an interatomic potential that will be applicable to a wide range of applications concerning diverse forms of bulk and nanostructured carbon., Comment: The following article has been accepted by The Journal of Chemical Physics. After it is published, it will be found at https://publishing.aip.org/resources/librarians/products/journals/

Polarons are ubiquitous in many semiconductors and have been linked with conductivity and optical response of materials for photovoltaics and heterogeneous catalysis, yet how surface polarons influence adsorption remains unclear. Here, by modelling the surface of rutile titania using density functional theory, we reveal the effect of small surface polarons on water adsorption, dissociation, and hydrogen bonding. On the one hand the presence of such polarons significantly suppresses dissociation of water molecules that are bonded directly to polaronic sites. On the other hand, polarons facilitate water dissociation at certain non-polaronic sites. Furthermore, polarons strengthen hydrogen bonds, which in turn affects water dissociation in hydrogen bonded overlayer structures. This study reveals that polarons at the rutile surface have complex, multi-faceted, effects on water adsorption, dissociation and hydrogen bonding, highlighting the importance of polarons on water structure and dynamics on such surfaces. We expect that many of the physical properties of surface polarons identified here will apply more generally to surfaces and interfaces that can host small polarons, beyond titania.

Due to their current and future technological applications, including realisation of water filters and desalination membranes, water adsorption on graphitic sp$^{2}$-bonded carbon is of overwhelming interest. However, these systems are notoriously challenging to model, even for electronic structure methods such as density functional theory (DFT), because of the crucial role played by London dispersion forces and non-covalent interactions in general. Recent efforts have established reference quality interactions of several carbon nanostructures interacting with water. Here, we compile a new benchmark set (dubbed \textbf{WaC18}), which includes a single water molecule interacting with a broad range of carbon structures, and various bulk (3D) and two dimensional (2D) ice polymorphs. The performance of 28 approaches, including semi-local exchange-correlation functionals, non-local (Fock) exchange contributions, and long-range van der Waals (vdW) treatments, are tested by computing the deviations from the reference interaction energies. The calculated mean absolute deviations on the WaC18 set depends crucially on the DFT approach, ranging from 135 meV for LDA to 12 meV for PBE0-D4. We find that modern vdW corrections to DFT significantly improve over their precursors. Within the 28 tested approaches, we identify the best performing within the functional classes of: generalized gradient approximated (GGA), meta-GGA, vdW-DF, and hybrid DF, which are BLYP-D4, TPSS-D4, rev-vdW-DF2, and PBE0-D4, respectively.

Fixed node diffusion quantum Monte Carlo (FN-DMC) is an increasingly used computational approach for investigating the electronic structure of molecules, solids, and surfaces with controllable accuracy. It stands out among equally accurate electronic structure approaches for its favorable cubic scaling with system size, which often makes FN-DMC the only computationally affordable high-quality method in large condensed phase systems with more than 100 atoms. In such systems FN-DMC deploys pseudopotentials to substantially improve efficiency. In order to deal with non-local terms of pseudopotentials, the FN-DMC algorithm must use an additional approximation, leading to the so-called localization error. However, the two available approximations, the locality approximation (LA) and the T-move approximation (TM), have certain disadvantages and can make DMC calculations difficult to reproduce. Here we introduce a third approach, called the determinant localization approximation (DLA). DLA eliminates reproducibility issues and systematically provides good quality results and stable simulations that are slightly more efficient than LA and TM. When calculating energy differences -- such as interaction and ionization energies -- DLA is also more accurate than the LA and TM approaches. We believe that DLA paves the way to the automization of FN-DMC and its much easier application in large systems.

Hydrogen is the most abundant element in the universe. It is also the lightest and as such the most quantum of the elements, in the sense that quantum tunnelling, quantum delocalisation, and zero-point motion can be important. For practical reasons most computer simulations of materials have not taken such effects into account, rather they have treated atomic nuclei as classical point-like particles. However, it is an exciting time for the theory and simulation of materials and thanks to significant methodological developments over the last few decades, nuclear quantum effects can now be accurately treated in complex materials. In this brief review we discuss our recent studies on the role nuclear quantum effects play in hydrogen containing systems. We give examples of how the quantum nature of the nuclei has a significant impact on the location of the boundaries between phases in high pressure condensed hydrogen, including a dramatic lowering of the solid to liquid melting line. We show how nuclear quantum effects facilitate the dissociative adsorption of molecular hydrogen on solid surfaces and the diffusion of atomic hydrogen across surfaces; effects that are of relevance to the catalytic performance of the surfaces. Finally, we discuss how nuclear quantum effects alter the strength and structure of hydrogen bonds, including the hydrogen bonds in DNA. Overall these studies demonstrate that nuclear quantum effects can manifest in many different, interesting, and at times non-intuitive ways. Whilst historically it has been difficult to know in advance what influence nuclear quantum effects will have, some of the important conceptual foundations have now started to emerge.

When an ice crystal is born from liquid water two key changes occur: (i) the molecules order; and (ii) the mobility of the molecules drops as they adopt their lattice positions. Most research on ice nucleation (and crystallization in general) has focused on understanding the former with less attention paid to the latter. However, supercooled water exhibits fascinating and complex dynamical behavior, most notably dynamical heterogeneity (DH), a phenomenon where spatially separated domains of relatively mobile and immobile particles coexist. Strikingly, the microscopic connection between the DH of water and the nucleation of ice has yet to be unraveled directly at the molecular level. Here we tackle this issue via computer simulations which reveal that: (i) ice nucleation occurs in low-mobility regions of the liquid; (ii) there is a dynamical incubation period in which the mobility of the molecules drops prior to any ice-like ordering; and (iii) ice-like clusters cause arrested dynamics in surrounding water molecules. With this we establish a clear connection between dynamics and nucleation. We anticipate that our findings will pave the way for the examination of the role of dynamical heterogeneities in heterogeneous and solution-based nucleation.

Molecular adsorption on surfaces plays a central role in catalysis, corrosion, desalination, and many other processes of relevance to industry and the natural world. Few adsorption systems are more ubiquitous or of more widespread importance than those involving water and carbon, and for a molecular level understanding of such interfaces water monomer adsorption on graphene is a fundamental and representative system. This system is particularly interesting as it calls for an accurate treatment of electron correlation effects, as well as posing a practical challenge to experiments. Here, we employ many-body electronic structure methodologies that can be rigorously converged and thus provide faithful references for the molecule-surface interaction. In particular, we use diffusion Monte-Carlo (DMC), coupled cluster (CCSD(T)), as well as the random phase approximation (RPA) to calculate the strength of the interaction between water and an extended graphene surface. We establish excellent, sub-chemical, agreement between the complementary high-level methodologies, and an adsorption energy estimate in the most stable configuration of approximately -100\,meV is obtained. We also find that the adsorption energy is rather insensitive to the orientation of the water molecule on the surface, despite different binding motifs involving qualitatively different interfacial charge reorganisation. In producing the first demonstrably accurate adsorption energies for water on graphene this work also resolves discrepancies amongst previously reported values for this widely studied system. It also paves the way for more accurate and reliable studies of liquid water at carbon interfaces with cheaper computational methods, such as density functional theory and classical potentials.

Organic molecules such as steroids or amino acids form crystals that can facilitate the formation of ice -- arguably the most important phase transition on earth. However, the origin of the ice nucleating ability of organic crystals is still largely unknown. Here, we combine experiments and simulations to unravel the microscopic details of ice formation on cholesterol, a prototypical organic crystal widely used in cryopreservation. We find that cholesterol -- which is also a substantial component of cell membranes -- is an ice nucleating agent more potent than many inorganic substrates, including the mineral feldspar (one of the most active ice nucleating material in the atmosphere). Scanning electron microscopy measurements reveal a variety of morphological features on the surfaces of cholesterol crystals: this suggests that the topography of the surface is key to the broad range of ice nucleating activity observed (from -4 to -20 C\degree). In addition, we show via molecular simulations that cholesterol crystals aid the formation of ice nuclei in a unconventional fashion. Rather than providing a template for a flat ice-like contact layer (as found in the case of many inorganic substrates), the flexibility of the cholesterol surface and its low density of hydrophilic functional groups leads to the formation of molecular cages involving both water molecules and terminal hydroxyl groups of the cholesterol surface. We argue that the phenomenal ice nucleating activity of steroids such as cholesterol (and potentially of many other organic crystals) is due to (i) the ability of flexible hydrophilic surfaces to form unconventional ice-templating structures and (ii) the different nucleation sites offered by the diverse topography of the crystalline surfaces., Comment: Chem. Sci., 2018, accepted on 27 Aug 2018

The delicate balance between H-bonding and van der Waals interactions determine the stability,structure and chirality of many molecular and supramolecular aggregates weakly adsorbed on solid surfaces.Yet the inherent complexity of these systems makes their experimental study at the molecular level very challenging.Small alcohols adsorbed on metal surfaces have become a useful model system to gain fundamental insight into the interplay of such molecule-surface and molecule-molecule interactions.Here, through a combination of scanning tunneling microscopy and density functional theory,we compare and contrast the adsorption and self-assembly of a range of small alcohols from methanol to butanol on Au(111).We find that that longer chained alcohols prefer to form zigzag chains held together by extended H-bonded networks between adjacent molecules.When alcohols bind to a metal surface datively via one of the two lone electron pairs of the oxygen atom they become chiral.Therefore,the chain structures are formed by a H-bonded network between adjacent molecules with alternating adsorbed chirality.These chain structures accommodate longer alkyl tails through larger unit cells, while the position of the hydroxyl group within the alcohol molecule can produce denser unit cells that maximize intermolecular interactions.Interestingly,when intrinsic chirality is introduced into the molecule as in the case of 2-butanol the assembly changes completely and square packing structures with chiral pockets are observed. This is rationalized by the fact that the intrinsic chirality of the molecule directs the chirality of the adsorbed hydroxyl group meaning that heterochiral chain structures cannot form.Overall this study provides a general framework for understanding the effect of simple alcohol molecular adstructures on H-bonded aggregates and paves the way for rationalizing 2D chiral supramolecular assembly.