Search

Biphasic interfaces are complex but fascinating regimes that display a number of properties distinct from those of the bulk. The CO2–H2O 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 CO2–H2O interface, particularly concerning the interfacial tension and phase behavior of CO2 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 CO2–H2O interface. Interfacial tensions are predicted from 1 to 500 bars and found to be in close agreement with experiment at pressures for which experimental data are available. Structural analyses indicate the buildup of an adsorbed, saturated CO2 film forming at a low pressure (20 bars) with properties similar to those of the bulk liquid, but preferential perpendicular alignment with respect to the interface. The CO2 monolayer buildup 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. [ABSTRACT FROM AUTHOR]

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 20 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 characterize 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, revPBE-D3 and RSCAN reproduce reference absolute lattice energies with the smallest error, while 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. [ABSTRACT FROM AUTHOR]

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 Ih. 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. [ABSTRACT FROM AUTHOR]

The O vacancy (Ov) formation energy, EOv, 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 EOv 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. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rock salt MgO, producing the first accurate and well-converged determinations of EOv with this method. These reference values are used to benchmark exchange–correlation functionals in DFT, and we find that all the studied functionals underestimate EOv, 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. [ABSTRACT FROM AUTHOR]

Jiang and co-workers use high resolution STM to investigate the reaction and self-assembly of (3,6-dibromo-9,10-phenanthrenequinone, or DBPQ) molecules on Ag (100) and Ag (110) surfaces in order to understand the mechanism of bottom-up assembly on surfaces.[31] They show that, through the inclusion of multiple functional groups within a precursor molecule, it becomes possible to fabricate new low-dimensional materials with unique chemical, physical, and electronic properties. Herbst and Fransson consider the core-valence separation approximation that is often used in the calculation of core-level spectra.[5] They show how to quantify the errors in this approximation, thereby opening the door to error-quantified predictions relevant to x-ray spectroscopy. 153(16), 164108 (2020).10.1063/5.0019557 5 M. F. Herbst and T. Fransson, "Quantifying the error of the core-valence separation approximation", J. Chem. Phys. Zhu and co-workers tackle this problem for a model system containing a 2D semiconductor heterojunction and show convincingly the efficient hot electron transfer from photoexcited MoTe SB 2 sb to WS SB 2 sb .[30] This finding provides important insight into the competition between hot electron cooling and transfer at 2D semiconductor interfaces and suggests an intriguing possibility for the exploration of hot electron devices. [Extracted from the article]

Oxygen vacancies are ubiquitous in TiO2 and play key roles in catalysis and magnetism applications. Despite being extensively investigated, the electronic structure of oxygen vacancies in TiO2 remains controversial both experimentally and theoretically. Here, we report a study of a neutral oxygen vacancy in TiO2 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 phases of TiO2. Specifically, embedded coupled cluster with singles, doubles, and perturbative triples 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 into the electronic structure of the oxygen vacancy in TiO2, clarifying earlier controversies and potentially inspiring future studies of defects with correlated wave function theories. [ABSTRACT FROM AUTHOR]

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 optimization of the many-body smooth overlap of atomic positions 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 [V. L. Deringer and G. Csányi, Phys. Rev. B 95, 094203 (2017)] with the high numerical accuracy necessary for crystalline graphene [Rowe et al., 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. [ABSTRACT FROM AUTHOR]

Due to their current and future technological applications, including realization of water filters and desalination membranes, water adsorption on graphitic sp2-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 noncovalent 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 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 semilocal exchange-correlation functionals, nonlocal (Fock) exchange contributions, and long-range van der Waals (vdW) treatments, is tested by computing the deviations from the reference interaction energies. The calculated mean absolute deviations on the WaC18 set depend crucially on the DFT approach, ranging from 135 meV for local density approximation (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. [ABSTRACT FROM AUTHOR]

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 (PPs) to substantially improve efficiency. In order to deal with nonlocal terms of PPs, 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 automation of FN-DMC and its much easier application in large systems. [ABSTRACT FROM AUTHOR]

Ice nucleation plays a significant role in a large number of natural and technological processes, but it is challenging to investigate experimentally because of the small time scales (ns) and short length scales (nm) involved. On the other hand, conventional molecular simulations struggle to cope with the relatively long time scale required for critical ice nuclei to form. One way to tackle this issue is to take advantage of free energy or path sampling techniques. Unfortunately, these are computationally costly. Seeded molecular dynamics is a much less demanding alternative that has been successfully applied already to study the homogeneous freezing of water. However, in the case of heterogeneous ice nucleation, nature’s favourite route to form ice, an array of suitable interfaces between the ice seeds and the substrate of interest has to be built, and this is no trivial task. In this paper, we present a Heterogeneous SEEDing (HSEED) approach which harnesses a random structure search framework to tackle the ice-substrate challenge, thus enabling seeded molecular dynamics simulations of heterogeneous ice nucleation on crystalline surfaces. We validate the HSEED framework by investigating the nucleation of ice on (i) model crystalline surfaces, using the coarse-grained mW model, and (ii) cholesterol crystals, employing the fully atomistic TIP4P/ice water model. We show that the HSEED technique yields results in excellent agreement with both metadynamics and forward flux sampling simulations. Because of its computational efficiency, the HSEED method allows one to rapidly assess the ice nucleation ability of whole libraries of crystalline substrates—a long-awaited computational development in, e.g., atmospheric science. [ABSTRACT FROM AUTHOR]

The delicate balance between hydrogen bonding and van der Waals interactions determines 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. In this quest, 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 longer chained alcohols prefer to form zigzag chains held together by extended hydrogen 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 hydrogen-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 hydrogen bonded aggregates and paves the way for rationalizing 2D chiral supramolecular assembly. [ABSTRACT FROM AUTHOR]

The design and production of novel 2-dimensional materials have seen great progress in the last decade, prompting further exploration of the chemistry of such materials. Doping and hydrogenating graphene are an experimentally realised method of changing its surface chemistry, but there is still a great deal to be understood on how doping impacts on the adsorption of molecules. Developing this understanding is key to unlocking the potential applications of these materials. High throughput screening methods can provide particularly effective ways to explore vast chemical compositions of materials. Here, alchemical derivatives are used as a method to screen the dissociative adsorption energy of water molecules on variousBNdoped topologies of hydrogenated graphene. The predictions from alchemical derivatives are assessed by comparison to density functional theory. This screening method is found to predict dissociative adsorption energies that span a range of more than 2 eV, with a mean absolute error <0.1 eV. In addition, we show that the quality of such predictions can be readily assessed by examination of the Kohn-Sham highest occupied molecular orbital in the initial states. In this way, the root mean square error in the dissociative adsorption energies of water is reduced by almost an order of magnitude (down to ~0.02 eV) after filtering out poor predictions. The findings point the way towards a reliable use of first order alchemical derivatives for efficient screening procedures. [ABSTRACT FROM AUTHOR]

Non-bonded potentials are included in most force fields and therefore widely used in classical molecular dynamics simulations of materials and interfacial phenomena. It is commonplace to truncate these potentials for computational efficiency based on the assumption that errors are negligible for reasonable cutoffs or compensated for by adjusting other interaction parameters. Arising from a metadynamics study of the wetting transition of water on a solid substrate, we find that the influence of the cutoff is unexpectedly strong and can change the character of the wetting transition from continuous to first order by creating artificial metastable wetting states. Common cutoff corrections such as the use of a force switching function, a shifted potential, or a shifted force do not avoid this. Such a qualitative difference urges caution and suggests that using truncated non-bonded potentials can induce unphysical behavior that cannot be fully accounted for by adjusting other interaction parameters. [ABSTRACT FROM AUTHOR]

Molecular adsorption on surfaces plays an important part in catalysis, corrosion, desalination, and various other processes that are relevant to industry and in nature. As a complement to experiments, accurate adsorption energies can be obtained using various sophisticated electronic structure methods that can now be applied to periodic systems. The adsorption energy of water on boron nitride substrates, going from zero to 2-dimensional periodicity, is particularly interesting as it calls for an accurate treatment of polarizable electrostatics and dispersion interactions, as well as posing a practical challenge to experiments and electronic structure methods. Here, we present reference adsorption energies, static polarizabilities, and dynamic polarizabilities, for water on BN substrates of varying size and dimension. Adsorption energies are computed with coupled cluster theory, fixed-node quantum Monte Carlo (FNQMC), the random phase approximation, and second order Møller-Plesset theory. These wavefunction based correlated methods are found to agree in molecular as well as periodic systems. The best estimate of the water/h-BN adsorption energy is -107±7 meV from FNQMC. In addition, the water adsorption energy on the BN substrates could be expected to grow monotonically with the size of the substrate due to increased dispersion interactions, but interestingly, this is not the case here. This peculiar finding is explained using the static polarizabilities and molecular dispersion coefficients of the systems, as computed from time-dependent density functional theory (DFT). Dynamic as well as static polarizabilities are found to be highly anisotropic in these systems. In addition, the many-body dispersion method in DFT emerges as a particularly useful estimation of finite size effects for other expensive, many-body wavefunction based methods. [ABSTRACT FROM AUTHOR]

Prompted by this observation, the reference dataset used to fit the GAP-20 model was recomputed, using DFT calculations with moreconsistent reciprocal-space sampling than was used in the original training data. (b) bcc atomization energy for updated GAP-20 model (GAP-20U); the better regularization of the model reduces the magnitude of the fluctuations, and the bcc atomization energy remains high. (a) bcc atomization energy for the original GAP-20 model, showing a bcc carbon minimum below the atomization energy of graphite. [Extracted from the article]

The formation of ice affects many aspects of our everyday life as well as important technologies such as cryotherapy and cryopreservation. Foreign substances almost always aid water freezing through heterogeneous ice nucleation, but the molecular details of this process remain largely unknown. In fact, insight into the microscopic mechanism of ice formation on different substrates is difficult to obtain even if state-of-the-art experimental techniques are used. At the same time, atomistic simulations of heterogeneous ice nucleation frequently face extraordinary challenges due to the complexity of the water-substrate interaction and the long time scales that characterize nucleation events. Here, we have investigated several aspects of molecular dynamics simulations of heterogeneous ice nucleation considering as a prototypical ice nucleating material the clay mineral kaolinite, which is of relevance in atmospheric science. We show via seeded molecular dynamics simulations that ice nucleation on the hydroxylated (001) face of kaolinite proceeds exclusively via the formation of the hexagonal ice polytype. The critical nucleus size is two times smaller than that obtained for homogeneous nucleation at the same supercooling. Previous findings suggested that the flexibility of the kaolinite surface can alter the time scale for ice nucleation within molecular dynamics simulations. However, we here demonstrate that equally flexible (or non flexible) kaolinite surfaces can lead to very different outcomes in terms of ice formation, according to whether or not the surface relaxation of the clay is taken into account. We show that very small structural changes upon relaxation dramatically alter the ability of kaolinite to provide a template for the formation of a hexagonal overlayer ofwater molecules at the water-kaolinite interface, and that this relaxation therefore determines the nucleation ability of this mineral. [ABSTRACT FROM AUTHOR]

Novel uses for 2-dimensional materials like graphene and hexagonal boron nitride (h-BN) are being frequently discovered especially for membrane and catalysis applications. Still however, a great deal remains to be understood about the interaction of environmentally and industrially relevant molecules such as water with these materials. Taking inspiration from advances in hybridising graphene and h-BN, we explore using density functional theory, the dissociation of water, hydrogen, methane, and methanol on graphene, h-BN, and their isoelectronic doped counterparts: BN doped graphene and C doped h-BN. We find that doped surfaces are considerably more reactive than their pristine counterparts and by comparing the reactivity of several small molecules, we develop a general framework for dissociative adsorption. From this a particularly attractive consequence of isoelectronic doping emerges: substrates can be doped to enhance their reactivity specifically towards either polar or non-polar adsorbates. As such, these substrates are potentially viable candidates for selective catalysts and membranes, with the implication that a range of tuneable materials can be designed. [ABSTRACT FROM AUTHOR]

Kohn-Sham density functional theory (DFT) has become established as an indispensable tool for investigating aqueous systems of all kinds, including those important in chemistry, surface science, biology, and the earth sciences. Nevertheless, many widely used approximations for the exchangecorrelation (XC) functional describe the properties of pure water systems with an accuracy that is not fully satisfactory. The explicit inclusion of dispersion interactions generally improves the description, but there remain large disagreements between the predictions of different dispersioninclusive methods.We present here a review of DFT work on water clusters, ice structures, and liquid water, with the aim of elucidating how the strengths and weaknesses of different XC approximations manifest themselves across this variety of water systems. Our review highlights the crucial role of dispersion in describing the delicate balance between compact and extended structures of many different water systems, including the liquid. By referring to a wide range of published work, we argue that the correct description of exchange-overlap interactions is also extremely important, so that the choice of semi-local or hybrid functional employed in dispersion-inclusive methods is crucial. The origins and consequences of beyond-2-body errors of approximate XC functionals are noted, and we also discuss the substantial differences between different representations of dispersion. We propose a simple numerical scoring system that rates the performance of different XC functionals in describing water systems, and we suggest possible future developments [ABSTRACT FROM AUTHOR]

The assembly of complex structures in nature is driven by an interplay between several intermolecular interactions, from strong covalent bonds to weaker dispersion forces. Understanding and ultimately controlling the self-assembly of materials requires extensive study of how these forces drive local nanoscale interactions and how larger structures evolve. Surface-based self-assembly is particularly amenable to modeling and measuring these interactions in well-defined systems. This study focuses on 2-butanol, the simplest aliphatic chiral alcohol. 2-butanol has recently been shown to have interesting properties as a chiral modifier of surface chemistry; however, its mode of action is not fully understood and a microscopic understanding of the role non-covalent interactions play in its adsorption and assembly on surfaces is lacking. In order to probe its surface properties, we employed high-resolution scanning tunneling microscopy and density functional theory (DFT) simulations. We found a surprisingly rich degree of enantiospecific adsorption, association, chiral cluster growth and ultimately long range, highly ordered chiral templating. Firstly, the chiral molecules acquire a second chiral center when adsorbed to the surface via dative bonding of one of the oxygen atom lone pairs. This interaction is controlled via the molecule's intrinsic chiral center leading to monomers of like chirality, at both chiral centers, adsorbed on the surface. The monomers then associate into tetramers via a cyclical network of hydrogen bonds with an opposite chirality at the oxygen atom. The evolution of these square units is surprising given that the underlying surface has a hexagonal symmetry. Our DFT calculations, however, reveal that the tetramers are stable entities that are able to associate with each other by weaker van der Waals interactions and tessellate in an extended square network. This network of homochiral square pores grows to cover the whole Au(111) surface. Our data reveal that the chirality of a simple alcohol can be transferred to its surface binding geometry, drive the directionality of hydrogen-bonded networks and ultimately extended structure. Furthermore, this study provides the first microscopic insight into the surface properties of this important chiral modifier and provides a well-defined system for studying the network's enantioselective interaction with other molecules. [ABSTRACT FROM AUTHOR]

Ice formation is one of the most common and important processes on earth and almost always occurs at the surface of a material. A basic understanding of how the physicochemical properties of a material's surface affect its ability to form ice has remained elusive. Here, we use molecular dynamics simulations to directly probe heterogeneous ice nucleation at a hexagonal surface of a nanoparticle of varying hydrophilicity. Surprisingly, we find that structurally identical surfaces can both inhibit and promote ice formation and analogous to a chemical catalyst, it is found that an optimal interaction between the surface and the water exists for promoting ice nucleation. We use our microscopic understanding of the mechanism to design a modified surface in silico with enhanced ice nucleating ability. C 2015 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License. [ABSTRACT FROM AUTHOR]

Coarse grained molecular dynamics simulations are presented in which the sensitivity of the ice nucleation rate to the hydrophilicity of a graphene nanoflake is investigated. We find that an optimal interaction strength for promoting ice nucleation exists, which coincides with that found previously for a face centered cubic (111) surface. We further investigate the role that the layering of interfacial water plays in heterogeneous ice nucleation and demonstrate that the extent of layering is not a good indicator of ice nucleating ability for all surfaces. Our results suggest that to be an efficient ice nucleating agent, a surface should not bind water too strongly if it is able to accommodate high coverages of water. [ABSTRACT FROM AUTHOR]

Despite a recent flurry of experimental and simulation studies, an accurate estimate of the interaction strength of water molecules with hexagonal boron nitride is lacking. Here, we report quantum Monte Carlo results for the adsorption of a water monomer on a periodic hexagonal boron nitride sheet, which yield a water monomer interaction energy of -84 ± 5 meV. We use the results to evaluate the performance of several widely used density functional theory (DFT) exchange correlation functionals and find that they all deviate substantially. Differences in interaction energies between different adsorption sites are however better reproduced by DFT. [ABSTRACT FROM AUTHOR]

Density functional theory (DFT) studies of weakly interacting complexes have recently focused on the importance of van der Waals dispersion forces, whereas the role of exchange has received far less attention. Here, by exploiting the subtle binding between water and a boron and nitrogen doped benzene derivative (1,2-azaborine) we show how exact exchange can alter the binding conformation within a complex. Benchmark values have been calculated for three orientations of the water monomer on 1,2-azaborine from explicitly correlated quantum chemical methods, and we have also used diffusion quantum Monte Carlo. For a host of popular DFT exchange-correlation functionals we show that the lack of exact exchange leads to the wrong lowest energy orientation of water on 1,2-azaborine. As such, we suggest that a high proportion of exact exchange and the associated improvement in the electronic structure could be needed for the accurate prediction of physisorption sites on doped surfaces and in complex organic molecules. Meanwhile to predict correct absolute interaction energies an accurate description of exchange needs to be augmented by dispersion inclusive functionals, and certain non-local van der Waals functionals (optB88- and optB86b-vdW) perform very well for absolute interaction energies. Through a comparison with water on benzene and borazine (B3N3H6) we show that these results could have implications for the interaction of water with doped graphene surfaces, and suggest a possible way of tuning the interaction energy. [ABSTRACT FROM AUTHOR]

Methanol is a versatile chemical feedstock, fuel source, and energy storage material. Many reactions involving methanol are catalyzed by transition metal surfaces, on which hydrogen-bonded methanol overlayers form. As with water, the structure of these overlayers is expected to depend on a delicate balance of hydrogen bonding and adsorbate-substrate bonding. In contrast to water, however, relatively little is known about the structures methanol overlayers form and how these vary from one substrate to another. To address this issue, herein we analyze the hydrogen bonded networks that methanol forms as a function of coverage on three catalytically important surfaces, Au(111), Cu(111), and Pt(111), using a combination of scanning tunneling microscopy and density functional theory. We investigate the effect of intermolecular interactions, surface coverage, and adsorption energies on molecular assembly and compare the results to more widely studied water networks on the same surfaces. Two main factors are shown to direct the structure of methanol on the surfaces studied: the surface coverage and the competition between the methanol-methanol and methanol-surface interactions. Additionally, we report a new chiral form of buckled hexamer formed by surface bound methanol that maximizes the interactions between methanol monomers by sacrificing interactions with the surface. These results serve as a direct comparison of interaction strength, assembly, and chirality of methanol networks on Au(111), Cu(111), and Pt(111) which are catalytically relevant for methanol oxidation, steam reforming, and direct methanol fuel cells. [ABSTRACT FROM AUTHOR]

High quality reference data from diffusion Monte Carlo calculations are presented for bulk sI methane hydrate, a complex crystal exhibiting both hydrogen-bond and dispersion dominated interactions. The performance of some commonly used exchange-correlation functionals and all-atom point charge force fields is evaluated. Our results show that none of the exchange-correlation functionals tested are sufficient to describe both the energetics and the structure of methane hydrate accurately, while the point charge force fields perform badly in their description of the cohesive energy but fair well for the dissociation energetics. By comparing to ice Ih, we show that a good prediction of the volume and cohesive energies for the hydrate relies primarily on an accurate description of the hydrogen bonded water framework, but that to correctly predict stability of the hydrate with respect to dissociation to ice Ih and methane gas, accuracy in the water-methane interaction is also required. Our results highlight the difficulty that density functional theory faces in describing both the hydrogen bonded water framework and the dispersion bound methane. [ABSTRACT FROM AUTHOR]

Exploring the role of van der Waals (vdW) forces on the adsorption of molecules on extended metal surfaces has become possible in recent years thanks to exciting developments in density functional theory (DFT). Among these newly developed vdW-inclusive methods, interatomic vdW approaches that account for the nonlocal screening within the bulk [V. G. Ruiz, W. Liu, E. Zojer, M. Scheffler, and A. Tkatchenko, Phys. Rev. Lett. 108, 146103 (2012)] and improved nonlocal functionals [J. Klimeš, D. R. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 022201 (2010)] have emerged as promising candidates to account efficiently and accurately for the lack of long-range vdW forces in most popular DFT exchange-correlation functionals. Here we have used these two approaches to compute benzene adsorption on a range of close-packed (111) surfaces upon which it either physisorbs (Cu, Ag, and Au) or chemisorbs (Rh, Pd, Ir, and Pt).We have thoroughly compared the performance between the two classes of vdW-inclusive methods and when available compared the results obtained with experimental data. By examining the computed adsorption energies, equilibrium distances, and binding curves we conclude that both methods allow for an accurate treatment of adsorption at equilibrium adsorbate-substrate distances. To this end, explicit inclusion of electrodynamic screening in the interatomic vdW scheme and optimized exchange functionals in the case of nonlocal vdW density functionals is mandatory. Nevertheless, some discrepancies are found between these two classes of methods at large adsorbate-substrate separations. [ABSTRACT FROM AUTHOR]

The dissolution of NaCl in water is one of the most common everyday processes, yet it remains poorly understood at the molecular level. Here we report the results of an extensive density functional theory study in which the initial stages of NaCl dissolution have been examined at low water coverages. Our specific approach is to study how the energetic cost of moving an ion or a pair of ions to a less coordinated site at the surface of various NaCl crystals varies with the number of water molecules adsorbed on the surface. This "microsolvation" approach allows us to study the dependence of the defect energies on the number of water molecules in the cluster and thus to establish when and where dissolution becomes favorable. Moreover, this approach allows us to understand the roles of the individual ions and water molecules in the dissolution process. Consistent with previous work we identify a clear preference for dissolution of Cl ions over Na ions. However, the detailed information obtained here leads to the conclusion that the process is governed by the higher affinity of the water molecules to Na ions than to Cl ions. The Cl ions are released first as this exposes more Na ions at the surface creating favorable adsorption sites for water. We discuss how this mechanism is likely to be effective for other alkali halides. [ABSTRACT FROM AUTHOR]

Density-functional theory (DFT) has been widely used to study water and ice for at least 20 years. However, the reliability of different DFT exchange-correlation (xc) functionals for water remains a matter of considerable debate. This is particularly true in light of the recent development of DFT based methods that account for van der Waals (vdW) dispersion forces. Here, we report a detailed study with several xc functionals (semi-local, hybrid, and vdW inclusive approaches) on ice Ih and six proton ordered phases of ice. Consistent with our previous study [B. Santra, J. Klimesˇ, D. Alfè, A. Tkatchenko, B. Slater, A. Michaelides, R. Car, and M. Scheffler, Phys. Rev. Lett. 107, 185701 (2011)] which showed that vdW forces become increasingly important at high pressures, we find here that all vdW inclusive methods considered improve the relative energies and transition pressures of the high-pressure ice phases compared to those obtained with semi-local or hybrid xc functionals. However, we also find that significant discrepancies between experiment and the vdW inclusive approaches remain in the cohesive properties of the various phases, causing certain phases to be absent from the phase diagram. Therefore, room for improvement in the description of water at ambient and high pressures remains and we suggest that because of the stern test the high pressure ice phases pose they should be used in future benchmark studies of simulation methods for water. [ABSTRACT FROM AUTHOR]

Quantum Monte Carlo calculations with the diffusion Monte Carlo (DMC) method have been used to compute the binding energy curves of hydrogen on benzene, coronene, and graphene. The DMC results on benzene agree with both Mo\ller-Plessett second order perturbation theory (MP2) and coupled cluster with singles, doubles, and perturbative triples [CCSD(T)] calculations, giving an adsorption energy of ∼25 meV. For coronene, DMC agrees well with MP2, giving an adsorption energy of ∼40 meV. For physisorbed hydrogen on graphene, DMC predicts a very small adsorption energy of only 5 ± 5 meV. Density functional theory (DFT) calculations with various exchange-correlation functionals, including van der Waals corrected functionals, predict a wide range of binding energies on all three systems. The present DMC results are a step toward filling the gap in accurate benchmark data on weakly bound systems. These results can help us to understand the performance of current DFT based methods, and may aid in the development of improved approaches. [ABSTRACT FROM AUTHOR]

The impact of quantum nuclear effects on hydrogen (H-) bond strength has been inferred in earlier work from bond lengths obtained from path integral molecular dynamics (PIMD) simulations. To obtain a direct quantitative assessment of such effects, we use constrained-centroid PIMD simulations to calculate the free energy changes upon breaking the H-bonds in dimers of HF and water. Comparing ab initio simulations performed using PIMD and classical nucleus molecular dynamics (MD), we find smaller dissociation free energies with the PIMD method. Specifically, at 50 K, the H-bond in (HF)2 is about 30% weaker when quantum nuclear effects are included, while that in (H2O)2 is about 15% weaker. In a complementary set of simulations, we compare unconstrained PIMD and classical nucleus MD simulations to assess the influence of quantum nuclei on the structures of these systems. We find increased heavy atom distances, indicating weakening of the H-bond consistent with that observed by direct calculation of the free energies of dissociation. [ABSTRACT FROM AUTHOR]

To understand the performance of popular density-functional theory exchange-correlation (xc) functionals in simulations of liquid water, water monomers and dimers were extracted from a PBE simulation of liquid water and examined with coupled cluster with single and double excitations plus a perturbative correction for connected triples [CCSD(T)]. CCSD(T) reveals that most of the dimers are unbound compared to two gas phase equilibrium water monomers, largely because monomers within the liquid have distorted geometries. Of the three xc functionals tested, PBE and BLYP tend to predict too large dissociation energies between monomers within the dimers. We show that this is because the cost to distort the monomers to the geometries they adopt in the liquid is systematically underestimated with these functionals. PBE0 reproduces the CCSD(T) monomer deformation energies very well and consequently the dimer dissociation energies much more accurately than PBE and BLYP. Although this study is limited to water monomers and dimers, the results reported here may provide an explanation for the overstructured radial distribution functions routinely observed in BLYP and PBE simulations of liquid water and are of relevance to water in other phases and to other associated molecular liquids. [ABSTRACT FROM AUTHOR]

Density functional theory (DFT) molecular dynamics simulations of a thin (∼15 Å) water film on NaCl(001) have been performed, with the aim of understanding the structural and dynamic properties of this important interfacial water system. The interaction of the water film with the surface orders the water molecules in the immediate vicinity of the interface. This is reflected by oscillations in the planar-averaged water density distribution along the surface normal that extend to about 8 Å from the surface. The interaction with the substrate leaves many of the water molecules in the immediate vicinity with broken hydrogen bonds and as a consequence considerably reduced dipole moments. Indeed a clear correlation between the number of hydrogen bonds which a water molecule is involved in and its dipole moment for both water on NaCl and bulk water is observed. How the DFT results obtained here compare to those obtained with various empirical potentials is briefly discussed. [ABSTRACT FROM AUTHOR]

We address the nature of the bond between water molecules and metal surfaces through a systematic density-functional theory (DFT) study of H2O monomer adsorption on a series of close-packed transition metal surfaces: Ru(0001), Rh(111), Pd(111), and Ag(111). Aiming to understand the origin behind energetic and structural trends along the 4d series we employ a range of analysis tools such as the electron reactivity function, decomposition of densities of states, electron density differences, and inspection of individual Kohn–Sham orbitals. The results obtained from our DFT calculations allow us to rationalize the bonding between water and transition metal surfaces as a balance of covalent and electrostatic interactions. A frontier orbital scheme based on so-called two-center four-electron interactions between the molecular orbitals of H2O—mainly the 1b1— and d-band states of the surface proves incisive in understanding these systems. [ABSTRACT FROM AUTHOR]

Weak noncovalent interactions such as van der Waals and hydrogen bonding are ubiquitous in nature, yet their accurate description with electronic structure theories is challenging. Here we assess the ability of a variety of theories to describe a water-benzene binding energy curve. Specifically, we test Hartree–Fock, second-order Mo\ller–Plesset perturbation theory, coupled cluster, density functional theory with several exchange-correlation functionals with and without empirical vdW corrections, and quantum Monte Carlo (QMC). Given the relative paucity of QMC reports for noncovalent interactions, it is interesting to see that QMC and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are in very good agreement for most of the binding energy curve, although at short distances there are small deviations on the order of 20 meV. [ABSTRACT FROM AUTHOR]

Within the framework of the GGA+U implementation of density functional theory, we investigate atomistic and electronic structures of Au adsorbed on the stoichiometric and the defective CeO2{111} surfaces, in the latter of which either O or Ce vacancies are presented. We show that on the stoichiometric surface, the most stable adsorption site of Au is not on the top of the outermost O atoms, as previously suggested, but on a bridgelike site in which the Au directly binds to two O atoms. We suggest that on both sites, the original empty Ce 4f states near the Fermi level facilitate the oxidation of the Au; the preference of the Au for being on the bridgelike site is due to the larger O 2p-d(Au) mixing, accompanied by more significant electron redistributions. On the reduced surface with O vacancies, the most stable adsorption site of Au is near the vacancy position. Unlike that on the stoichiometric surface, strong ionic bonding character exists between Au and Ce, as the former becomes Auδ- due to the occupation of the 6s(Au) orbitals. Upon substitution for one of the Ce atoms in the lattice, the Au possesses a much stronger positive charge than that in other cases. We find that although Au is strongly bonded when it is at the Ce vacancy site, the overall binding (i.e., with the Ce vacancy formation energy being taken into account) is weaker than that for Au adsorbed at the stoichiometric surface. [ABSTRACT FROM AUTHOR]

The ability of several density-functional theory (DFT) exchange-correlation functionals to describe hydrogen bonds in small water clusters (dimer to pentamer) in their global minimum energy structures is evaluated with reference to second order Mo\ller-Plesset perturbation theory (MP2). Errors from basis set incompleteness have been minimized in both the MP2 reference data and the DFT calculations, thus enabling a consistent systematic evaluation of the true performance of the tested functionals. Among all the functionals considered, the hybrid X3LYP and PBE0 functionals offer the best performance and among the nonhybrid generalized gradient approximation functionals, mPWLYP and PBE1W perform best. The popular BLYP and B3LYP functionals consistently underbind and PBE and PW91 display rather variable performance with cluster size. [ABSTRACT FROM AUTHOR]

Biphasic interfaces are complex but fascinating regimes that display a number of properties distinct from those of the bulk. The CO2–H2O 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 CO2–H2O interface, particularly concerning the interfacial tension and phase behavior of CO2 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 CO2–H2O interface. Interfacial tensions are predicted from 1 to 500 bars and found to be in close agreement with experiment at pressures for which experimental data are available. Structural analyses indicate the buildup of an adsorbed, saturated CO2 film forming at a low pressure (20 bars) with properties similar to those of the bulk liquid, but preferential perpendicular alignment with respect to the interface. The CO2 monolayer buildup 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. [ABSTRACT FROM AUTHOR]

This Special Topic section on Advanced Electronic Structure Methods for Solids and Surfaces contains a collection of research papers that showcase recent advances in the high accuracy prediction of materials and surface properties. It provides a timely snapshot of a growing field that is of broad importance to chemistry, physics, and materials science. [ABSTRACT FROM AUTHOR]

The interaction of water molecules with metal surfaces is typically weak and as a result van der Waals (vdW) forces can be expected to be of importance. Here we account for the systematic poor treatment of vdW forces in most popular density functional theory exchange-correlation functionals by applying accurate non-local vdW density functionals. We have computed the adsorption of a variety of exemplar systems including water monomer adsorption on Al(111), Cu(111), Cu(110), Ru(0001), Rh(111), Pd(111), Ag(111), Pt(111), and unreconstructed Au(111), and small clusters (up to 6 waters) on Cu(110). We show that non-local correlations contribute substantially to the water-metal bond in all systems, whilst water-water bonding is much less affected by non-local correlations. Interestingly non-local correlations contribute more to the adsorption of water on the reactive transition metal substrates than they do on the noble metals. The relative stability, adsorption sites, and adsorption geometries of competing water adstructures rarely differ when comparing results obtained with semi-local functionals and the non-local vdW density functionals, which explains the previous success of semi-local functionals in characterizing adsorbed water structures on a number of metal surfaces. [ABSTRACT FROM AUTHOR]

Electron dispersion forces play a crucial role in determining the structure and properties of biomolecules, molecular crystals, and many other systems. However, an accurate description of dispersion is highly challenging, with the most widely used electronic structure technique, density functional theory (DFT), failing to describe them with standard approximations. Therefore, applications of DFT to systems where dispersion is important have traditionally been of questionable accuracy. However, the last decade has seen a surge of enthusiasm in the DFT community to tackle this problem and in so-doing to extend the applicability of DFT-based methods. Here we discuss, classify, and evaluate some of the promising schemes to emerge in recent years. A brief perspective on the outstanding issues that remain to be resolved and some directions for future research are also provided. [ABSTRACT FROM AUTHOR]

Oxygen vacancies are ubiquitous in TiO2 and play key roles in catalysis and magnetism applications. Despite being extensively investigated, the electronic structure of oxygen vacancies in TiO2 remains controversial both experimentally and theoretically. Here, we report a study of a neutral oxygen vacancy in TiO2 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 phases of TiO2. Specifically, embedded coupled cluster with singles, doubles, and perturbative triples 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 into the electronic structure of the oxygen vacancy in TiO2, clarifying earlier controversies and potentially inspiring future studies of defects with correlated wave function theories. [ABSTRACT FROM AUTHOR]

Biphasic interfaces are complex but fascinating regimes that display a number of properties distinct from those of the bulk. The CO2-H2O 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 CO2-H2O interface, particularly concerning the interfacial tension and phase behavior of CO2 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 CO2-H2O interface. Interfacial tensions are predicted from 1 to 500 bars and found to be in close agreement with experiment at pressures for which experimental data are available. Structural analyses indicate the buildup of an adsorbed, saturated CO2 film forming at a low pressure (20 bars) with properties similar to those of the bulk liquid, but preferential perpendicular alignment with respect to the interface. The CO2 monolayer buildup 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., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

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 phases of TiO 2 . Specifically, embedded coupled cluster with singles, doubles, and perturbative triples 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 into 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.