BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Angelos Michaelides

Angelos Michaelides appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.chem-phphysics.comp-phcond-mat.softcond-mat.othercond-mat.mes-hallcond-mat.stat-mechcond-mat.dis-nncond-mat.str-elquant-ph

Catalog footprint

What is connected

48works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

48 published item(s)

preprint2023arXiv

Classical quantum friction at water-carbon interfaces

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Can molecular simulations reliably compare homogeneous and heterogeneous ice nucleation?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

How do interfaces alter the dynamics of supercooled water?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Long range ionic and short range hydration effects govern strongly anisotropic clay nanoparticle interactions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Tracking single atoms in a liquid environment

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Water flow in single-wall nanotubes: Oxygen makes it slip, hydrogen makes it stick

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Rapid water diffusion at cryogenic temperatures through an inchworm-like mechanism

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Accurate and Transferable Machine Learning Potential for Carbon

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The color center singlet state of oxygen vacancies in TiO$_2$

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Small polarons and the Janus nature of $\text{TiO}_\text{2}(110)$

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Boosting the accuracy and speed of quantum Monte Carlo: size-consistency and time-step

Diffusion Monte Carlo (DMC) simulations for fermions are becoming the standard to provide high quality reference data in systems that are too large to be investigated via quantum chemical approaches. DMC with the fixed-node approximation relies on modifications of the Green function to avoid singularities near the nodal surface of the trial wavefunction. We show that these modifications affect the DMC energies in a way that is not size-consistent, resulting in large time-step errors. Building on the modifications of Umrigar {\em et al.} and of DePasquale {\em et al.} we propose a simple Green function modification that restores size-consistency to large values of time-step; substantially reducing the time-step errors. The new algorithm also yields remarkable speedups of up to two orders of magnitude in the calculation of molecule-molecule binding energies and crystal cohesive energies, thus extending the horizons of what is possible with DMC.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insight into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that in the last few decades have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state of the art computational methods, by reviewing simulations of e.g. ice nucleation or crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insight into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that in doing so the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that by improving (i.) existing interatomic potentials; and (ii.) currently available enhanced sampling methods, the community can move towards accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evidence for Stable Square Ice from Quantum Monte Carlo

Recent experiments on ice formed by water under nanoconfinement provide evidence for a two-dimensional (2D) `square ice' phase. However, the interpretation of the experiments has been questioned and the stability of square ice has become a matter of debate. Partially this is because the simulation approaches employed so far (force fields and density functional theory) struggle to accurately describe the very small energy differences between the relevant phases. Here we report a study of 2D ice using an accurate wave-function based electronic structure approach, namely Diffusion Monte Carlo (DMC). We find that at relatively high pressure square ice is indeed the lowest enthalpy phase examined, supporting the initial experimental claim. Moreover, at lower pressures a `pentagonal ice' phase (not yet observed experimentally) has the lowest enthalpy, and at ambient pressure the `pentagonal ice' phase is degenerate with a `hexagonal ice' phase. Our DMC results also allow us to evaluate the accuracy of various density functional theory exchange correlation functionals and force field models, and in doing so we extend the understanding of how such methodologies perform to challenging 2D structures presenting dangling hydrogen bonds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fast diffusion of water nanodroplets on graphene

Diffusion across surfaces generally involves motion on a vibrating but otherwise stationary substrate. Here, using molecular dynamics, we show that a layered material such as graphene opens up a new mechanism for surface diffusion whereby adsorbates are carried by propagating ripples via a motion similar to surfing. For water nanodroplets, we demonstrate that the mechanism leads to exceedingly fast diffusion that is 2-3 orders of magnitude faster than the self-diffusion of water molecules in liquid water. We also reveal the underlying principles that regulate this new mechanism for diffusion and show how it also applies to adsorbates other than water, thus opening up the prospect of achieving fast and controllable motion of adsorbates across material surfaces more generally.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ice Formation on Kaolinite: Insights from Molecular Dynamics Simulations

The formation of ice affects many aspects of our everyday life as well as 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 via state-of-the-art experimental techniques. 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 timescales 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 of water molecules at the water-kaolinite interface, and that this relaxation therefore determines the nucleation ability of this mineral.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs

Despite the inherently quantum mechanical nature of hydrogen bonding, it is unclear how nuclear quantum effects (NQEs) alter the strengths of hydrogen bonds. With this in mind, we use ab initio path integral molecular dynamics to determine the absolute contribution of NQEs to the binding in DNA base pair complexes, arguably the most important hydrogen-bonded systems of all. We find that depending on the temperature, NQEs can either strengthen or weaken the binding within the hydrogen-bonded complexes. As a somewhat counterintuitive consequence, NQEs can have a smaller impact on hydrogen bond strengths at cryogenic temperatures than at room temperature. We rationalize this in terms of a competition of NQEs between low-frequency and high-frequency vibrational modes. Extending this idea, we also propose a simple model to predict the temperature dependence of NQEs on hydrogen bond strengths in general.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Microscopic Mechanism and Kinetics of Ice Formation at Complex Interfaces: Zooming in on Kaolinite

Most ice in nature forms thanks to impurities which boost the exceedingly low nucleation rate of pure supercooled water. However, the microscopic details of ice nucleation on these substances remain largely unknown. Here, we have unraveled the molecular mechanism and the kinetics of ice formation on kaolinite, a clay mineral playing a key role in climate science. We find that the formation of ice at strong supercooling in the presence of this clay is twenty orders of magnitude faster than homogeneous freezing. The critical nucleus is substantially smaller than that found for homogeneous nucleation and, in contrast to the predictions of classical nucleation theory (CNT), it has a strong 2D character. Nonetheless, we show that CNT describes correctly the formation of ice at this complex interface. Kaolinite also promotes the exclusive nucleation of hexagonal ice, as opposed to homogeneous freezing where a mixture of cubic and hexagonal polytypes is observed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stability of Complex Biomolecular Structures: Vander Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects

Biomolecules are complex systems stabilized by a delicate balance of weak interactions, making it important to assess all energetic contributions in an accurate manner. However, it is a priori unclear which contributions make more of an impact. Here, we examine stacked polyglutamine (polyQ) strands, a peptide repeat often found in amyloid aggregates. We investigate the role of hydrogen bond (HB) cooperativity, van der Waals (vdW) dispersion interactions, and quantum contributions to free energies, including anharmonicities through density functional theory and ab initio path integral simulations. Of these various factors, we find that the largest impact on structural stabilization comes from vdW interactions. HB cooperativity is the second largest contribution as the size of the stacked chain grows. Competing nuclear quantum effects make the net quantum contribution small but very sensitive to anharmonicities, vdW, and the number of HBs. Our results suggest that a reliable treatment of these systems can only be attained by considering all of these components.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Interplay of Covalency, Hydrogen Bonding and Dispersion Leads to a Long Range Chiral Network: The Example of 2-Butanol

The assembly of complex structures is driven by an interplay between several intermolecular interactions, from strong covalent bonds to weaker dispersion forces. 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 shows interesting properties as a chiral modifier of surface chemistry, however, its mode of action is not fully understood. In order to probe its surface properties we employed high-resolution scanning tunneling microscopy and 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 reveal that the tetramers are able to associate with each other by weaker van der Waals interactions and tessellate in an extended square network. Our data reveals that the chirality of a simple alcohol can be transferred to its surface binding geometry, drive the directionality of hydrogen bonds 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Many Faces of Heterogeneous Ice Nucleation: Interplay Between Surface Morphology and Hydrophobicity

What makes a material a good ice nucleating agent? Despite the importance of heterogeneous ice nucleation to a variety of fields, from cloud science to microbiology, major gaps in our understanding of this ubiquitous process still prevent us from answering this question. In this work, we have examined the ability of generic crystalline substrates to promote ice nucleation as a function of the hydrophobicity and the morphology of the surface. Nucleation rates have been obtained by brute-force molecular dynamics simulations of coarse-grained water on top of different surfaces of a model fcc crystal, varying the water-surface interaction and the surface lattice parameter. It turns out that the lattice mismatch of the surface with respect to ice, customarily regarded as the most important requirement for a good ice nucleating agent, is at most desirable but not a requirement. On the other hand, the balance between the morphology of the surface and its hydrophobicity can significantly alter the ice nucleation rate and can also lead to the formation of up to three different faces of ice on the same substrate. We have pinpointed three circumstances where heterogeneous ice nucleation can be promoted by the crystalline surface: (i) the formation of a water overlayer that acts as an in-plane template; (ii) the emergence of a contact layer buckled in an ice-like manner; and (iii) nucleation on compact surfaces with very high interaction strength. We hope that this extensive systematic study will foster future experimental work aimed at testing the physiochemical understanding presented herein.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Toward Accurate Adsorption Energetics on Clay Surfaces

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tuning dissociation using isoelectronically doped graphene and hexagonal boron nitride: water and other small molecules

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 elevant 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Water on BN doped benzene: A hard test for exchange-correlation functionals and the impact of exact exchange on weak binding

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 (B$_3$N$_3$H$_6$) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

2D ice from first principles: structures and phase transitions

Despite relevance to disparate areas such as cloud microphysics and tribology, major gaps in the understanding of the structures and phase transitions of low-dimensional water ice remain. Here we report a first principles study of confined 2D ice as a function of pressure. We find that at ambient pressure hexagonal and pentagonal monolayer structures are the two lowest enthalpy phases identified. Upon mild compression the pentagonal structure becomes the most stable and persists up to \textit{ca.} 2 GPa at which point square and rhombic phases are stable. The square phase agrees with recent experimental observations of square ice confined within graphene sheets. We also find a double layer AA stacked square ice phase, which clarifies the difference between experimental observations and earlier force field simulations. This work provides a fresh perspective on 2D confined ice, highlighting the sensitivity of the structures observed to both the confining pressure and width.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Atomistic details of oxide surfaces and surface oxidation: the example of copper and its oxides

The oxidation and corrosion of metals are fundamental problems in materials science and technology that have been studied using a large variety of experimental and computational techniques. Here we review some of the recent studies that have led to significant advances in our atomic-level understanding of copper oxide, one of the most studied and best understood metal oxides. We show that a good atomistic understanding of the physical characteristics of cuprous (Cu$_2$O) and cupric (CuO) oxide and of some key processes of their formation has been obtained. Indeed, the growth of the oxide has been shown to be epitaxial with the surface and to proceed, in most cases, through the formation of oxide nano-islands which, with continuous oxygen exposure, grow and eventually coalesce. We also show how electronic structure calculations have become increasingly useful in helping to characterise the structures and energetics of various Cu oxide surfaces. However a number of challenges remain. For example, it is not clear under which conditions the oxidation of copper in air at room temperature (known as native oxidation) leads to the formation of a cuprous oxide film only, or also of a cupric overlayer. Moreover, the atomistic details of the nucleation of the oxide islands are still unknown. We close our review with a brief perspective on future work and discuss how recent advances in experimental techniques, bringing greater temporal and spatial resolution, along with improvements in the accuracy, realism and timescales achievable with computational approaches make it possible for these questions to be answered in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Friction of Water on Graphene and Hexagonal Boron Nitride from ab initio Methods: Very Different Slippage Despite Very Similar Interface Structures

Friction is one of the main sources of dissipation at liquid water/solid interfaces. Despite recent progress, a detailed understanding of water/solid friction in connection with the structure and energetics of the solid surface is lacking. Here we show for the first time that \textit{ab initio} molecular dynamics can be used to unravel the connection between the structure of nanoscale water and friction for liquid water in contact with graphene and with hexagonal boron nitride. We find that whilst the interface presents a very similar structure between the two sheets, the friction coefficient on boron nitride is $\approx 3$ times larger than that on graphene. This comes about because of the greater corrugation of the energy landscape on boron nitride arising from specific electronic structure effects. We discuss how a subtle dependence of the friction on the atomistic details of a surface, that is not related to its wetting properties, may have a significant impact on the transport of water at the nanoscale, with implications for the development of membranes for desalination and for osmotic power harvesting.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Molecular simulations of heterogeneous ice nucleation. I. Controlling ice nucleation through surface hydrophilicity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Molecular simulations of heterogeneous ice nucleation. II. Peeling back the layers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Understanding corrosion inhibition with van der Waals DFT methods: the case of benzotriazole

The corrosion of materials is an undesirable and costly process affecting many areas of technology and everyday life. As such, considerable effort has gone into understanding and preventing it. Organic molecule based coatings can in certain circumstances act as effective corrosion inhibitors. Although they have been used to great effect for more than sixty years, how they function at the atomic-level is still a matter of debate. In this work, computer simulation approaches based on density functional theory are used to investigate benzotriazole (BTAH), one of the most widely used and studied corrosion inhibitors for copper. In particular, the structures formed by protonated and deprotonated BTAH molecules on Cu(111) have been determined and linked to their inhibiting properties. It is found that hydrogen bonding, van der Waals interactions and steric repulsions all contribute in shaping how BTAH molecules adsorb, with flat-lying structures preferred at low coverage and upright configurations preferred at high coverage. The interaction of the dehydrogenated benzotriazole molecule (BTA) with the copper surface is instead dominated by strong chemisorption via the azole moiety with the aid of copper adatoms. Structures of dimers or chains are found to be the most stable structures at all coverages, in good agreement with scanning tunnelling microscopy results. Benzotriazole thus shows a complex phase behaviour in which van der Waals forces play an important role and which depends on coverage and on its protonation state and all of these factors feasibly contribute to its effectiveness as a corrosion inhibitor.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Water on hexagonal boron nitride from diffusion Monte Carlo

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Benchmarking the performance of Density Functional Theory and Point Charge Force Fields in their Description of sI Methane Hydrate against Diffusion Monte Carlo

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, whilst 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Insight into the description of van der Waals forces for benzene adsorption on transition metal (111) surfaces

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. Klimes, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nature of proton transport in a water-filled carbon nanotube and in liquid water

Proton transport (PT) in bulk liquid water and within a thin water-filled carbon nanotube has been examined with ab initio pathintegral molecular dynamics (PIMD). Barrierless proton transfer is observed in each case when quantum nuclear effects (QNEs) are accounted for. The key difference between the two systems is that in the nanotube facile PT is facilitated by a favorable prealignment of water molecules, whereas in bulk liquid water solvent reorganization is required prior to PT. Configurations where the quantum excess proton is delocalized over several adjacent water molecules along with continuous interconversion between different hydration states reveals that, as in liquid water, the hydrated proton under confinement is best described as a fluxional defect, rather than any individual idealized hydration state such as Zundel, Eigen, or the so-called linear H7O3+ complex along the water chain. These findings highlight the importance of QNEs in intermediate strength hydrogen bonds (HBs) and explain why H+ diffusion through nanochannels is impeded much less than other cations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Accuracy of van der Waals Inclusive Density-Functional Theory Exchange-Correlation Functionals for Ice at Ambient and High Pressures

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 [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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum effects in the diffusion of hydrogen on Ru(0001)

An understanding of hydrogen diffusion on metal surfaces is important, not just for its role in heterogeneous catalysis and hydrogen fuel cell technology, but also because it provides model systems where tunneling can be studied under well-defined conditions. Here we report helium spin-echo measurements of the atomic-scale motion of hydrogen on the Ru(0001) surface between 75 and 250 K. Quantum effects are evident at temperatures as high as 200 K, while below 120 K we observe a tunneling-dominated temperature independent jump rate of 1.9$\times$10$^9$ s$^{-1}$, many orders of magnitude faster than previously seen. Quantum transition state theory calculations based on ab initio path-integral simulations reproduce the temperature dependence of the rate at higher temperatures and predict a crossover to tunneling-dominated diffusion at low temperatures, although the tunneling rate is under-estimated, highlighting the need for future experimental and theoretical studies of hydrogen diffusion on well-defined surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Solvent Induced Proton Hopping at a Water-Oxide Interface

Despite widespread interest, a detailed understanding of the dynamics of proton transfer at interfaces is lacking. Here we use ab initio molecular dynamics to unravel the connection between interfacial water structure and proton transfer for the widely studied and experimentally well-characterized water-ZnO$(10\bar{1}0)$ interface. We find that upon going from a single layer of adsorbed water to a liquid multilayer changes in the structure are accompanied by a dramatic increase in the proton transfer rate at the surface. We show how hydrogen bonding and rather specific hydrogen bond fluctuations at the interface are responsible for the change in the structure and proton transfer dynamics. The implications of this for the chemical reactivity and for the modelling of complex wet oxide interfaces in general are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Microscopic Features of Heterogeneous Ice Nucleation May Affect the Macroscopic Morphology of Atmospheric Ice Crystals

It is surprisingly difficult to freeze water. Almost all ice that forms under "mild" conditions (temperatures > -40 degrees Celsius) requires the presence of a nucleating agent - a solid particle that facilitates the freezing process - such as clay mineral dust, soot or bacteria. In a computer simulation, the presence of such ice nucleating agents does not necessarily alleviate the difficulties associated with forming ice on accessible timescales. Nevertheless, in this work we present results from molecular dynamics simulations in which we systematically compare homogeneous and heterogeneous ice nucleation, using the atmospherically important clay mineral kaolinite as our model ice nucleating agent. From our simulations, we do indeed find that kaolinite is an excellent ice nucleating agent but that contrary to conventional thought, non-basal faces of ice can nucleate at the basal face of kaolinite. We see that in the liquid phase, the kaolinite surface has a drastic effect on the density profile of water, with water forming a dense, tightly bound first contact layer. Monitoring the time evolution of the water density reveals that changes away from the interface may play an important role in the nucleation mechanism. The findings from this work suggest that heterogeneous ice nucleating agents may not only enhance the ice nucleation rate, but also alter the macroscopic structure of the ice crystals that form.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A molecular perspective of water at metal interfaces

Water/solid interfaces are relevant to a broad range of physicochemical phenomena and technological processes such as corrosion, lubrication, heterogeneous catalysis and electrochemistry. Although many fields have contributed to rapid progress in the fundamental knowledge of water at interfaces, detailed molecular-level understanding of water/solid interfaces comes mainly from studies on flat metal substrates. These studies have recently shown that a remarkably rich variety of structures form at the interface between water and even seemingly simple flat surfaces. In this Review we discuss the most exciting work in this area, in particular the emerging physical insight and general concepts about how water binds to metal surfaces. We also provide a perspective on outstanding problems, challenges and open questions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Classical and quantum ordering of protons in cold solid hydrogen under megabar pressures

A combination of state-of-the-art theoretical methods has been used to obtain an atomic-level picture of classical and quantum ordering of protons in cold high-pressure solid hydro-gen. We focus mostly on phases II and III of hydrogen, exploring the effects of quantum nuclear motion on certain features of these phases (through a number of ab initio path-integral molecular dynamics (PIMD) simulations at particular points on the phase diagram). We also examine the importance of van der Waals forces in this system by performing calculations using the optB88-vdW density functional, which accounts for non-local correlations. Our calculations reveal that the transition between phases I and II is strongly quantum in nature, resulting from a competition between anisotropic inter-molecular interactions that restrict molecular rotation and thermal plus quantum fluctuations of the nuclear positions that facilitate it. The transition from phase II to III is more classical because quantum nuclear motion plays only a secondary role and the transition is determined primarily by the underlying potential energy surface. A structure of P21c symmetry with 24 atoms in the primitive unit cell is found to be stable when anharmonic quantum nuclear vibrational motion is included at finite temperatures using the PIMD method. This structure gives a good account of the infra-red (IR) and Raman vibron frequencies of phase II. We find ad-ditional support for C2/c structure as a strong candidate for phase III, since it remains transparent up to 300 GPa, even when quantum nuclear effects are included. Finally, we find that accounting for van der Waals forces improves the agreement between experiment and theory for the parts of the phase diagram considered, when compared to previous work which employed the widely-used Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Perspective: Advances and challenges in treating van der Waals dispersion forces in density functional theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The role of van der Waals forces in water adsorption on metals

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 (DFT) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Understanding the role of ions and water molecules in the NaCl dissolution process

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Benzene Adsorbed on Metals: Concerted Effect of Covalency and van der Waals Bonding

The adsorption of aromatic molecules on metal surfaces plays a key role in condensed matter physics and functional materials. Depending on the strength of the interaction between the molecule and the surface, the binding is typically classified as either physisorption or chemisorption. Van der Waals (vdW) interactions contribute significantly to the binding in physisorbed systems, but the role of the vdW energy in chemisorbed systems remains unclear. Here we study the interaction of benzene with the (111) surface of transition metals, ranging from weak adsorption (Ag and Au) to strong adsorption (Pt, Pd, Ir, and Rh). When vdW interactions are accurately accounted for, the barrier to adsorption predicted by standard density functional theory (DFT) calculations essentially vanishes, producing a metastable precursor state on Pt and Ir surfaces. Notably, vdW forces contribute more to the binding of covalently bonded benzene than they do when benzene is physisorbed. Comparison to experimental data demonstrates that some of the recently developed methods for including vdW interactions in DFT allow quantitative treatment of both weakly and strongly adsorbed aromatic molecules on metal surfaces, extending the already excellent performance found for gas-phase molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantum simulation of low-temperature metallic liquid hydrogen

Experiments and computer simulations have shown that the melt-ing temperature of solid hydrogen drops with pressure above about 65 GPa, suggesting that a liquid state might exist at low temperatures. It has also been suggested that this low temperature liquid state might be non-molecular and metallic, although evidence for such behaviour is lacking. Here, we report results for hydrogen at high pressures using ab initio path-integral molecular dynamics methods, which include a description of the quantum motion of the protons at finite temperatures. We have determined the melting temperature as a function of pressure by direct simulation of the coexistence of the solid and liquid phases and have found an atomic solid phase from 500 to 800 GPa which melts at <200 K. Beyond this and up to pressures of 1200 GPa a metallic atomic liquid is stable at temperatures as low as 50 K. The quantum motion of the protons is critical to the low melting temperature in this system as ab initio simulations with classical nuclei lead to a considerably higher melting temperature of \sim300 K across the entire pressure range considered.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

The first principles approaches, density functional theory (DFT) and quantum Monte Carlo, have been used to examine the balance between van der Waals (vdW) forces and hydrogen (H) bonding in ambient and high pressure phases of ice. At higher pressure, the contribution to the lattice energy from vdW increases and that from H bonding decreases, leading vdW to have a substantial effect on the transition pressures between the crystalline ice phases. An important consequence, likely to be of relevance to molecular crystals in general, is that transition pressures obtained from DFT functionals which neglect vdW forces are greatly overestimated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Van der Waals density functionals applied to solids

The van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)] is a promising approach for including dispersion in approximate density functional theory exchange-correlation functionals. Indeed, an improved description of systems held by dispersion forces has been demonstrated in the literature. However, despite many applications, standard general tests on a broad range of materials are lacking. Here we calculate the lattice constants, bulk moduli, and atomization energies for a range of solids using the original vdW-DF and several of its offspring. We find that the original vdW-DF overestimates lattice constants in a similar manner to how it overestimates binding distances for gas phase dimers. However, some of the modified vdW functionals lead to average errors which are similar to those of PBE or better. Likewise, atomization energies that are slightly better than from PBE are obtained from the modified vdW-DFs. Although the tests reported here are for "hard" solids, not normally materials for which dispersion forces are thought to be important, we find a systematic improvement in cohesive properties for the alkali metals and alkali halides when non-local correlations are accounted for.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

To wet or not to wet? Dispersion forces tip the balance for water-ice on metals

Despite widespread discussion, the role of van der Waals dispersion forces in wetting remains unclear. Here we show that non-local correlations contribute substantially to the water-metal bond and that this is an important factor in governing the relative stabilities of wetting layers and 3D bulk ice. Due to the greater polarizability of the substrate metal atoms, non-local correlations between water and the metal exceed those between water molecules within ice. This sheds light on a long-standing problem, wherein common density functional theory exchange-correlation functionals incorrectly predict that none of the low temperature experimentally characterized ice-like wetting layers are thermodynamically stable.

0 citations0 reviewsNo code linkedOpen access