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

Ultralow-loading platinum-cobalt fuel cell catalysts derived from imidazolate frameworks

Achieving high catalytic performance with the lowest possible amount of platinum is critical for fuel cell cost reduction. Here we describe a method of preparing highly active yet stable electrocatalysts containing ultralow-loading platinum content by using cobalt or bimetallic cobalt and zinc zeolitic imidazolate frameworks as precursors. Synergistic catalysis between strained platinum-cobalt core-shell nanoparticles over a platinum-group metal (PGM)-free catalytic substrate led to excellent fuel cell performance under 1 atmosphere of O₂ or air at both high-voltage and high-current domains. Two catalysts achieved oxygen reduction reaction (ORR) mass activities of 1.08 amperes per milligram of platinum (A mg_Pt ^-1) and 1.77 A mg_Pt ^-1 and retained 64% and 15% of initial values after 30,000 voltage cycles in a fuel cell. Computational modeling reveals that the interaction between platinum-cobalt nanoparticles and PGM-free sites improves ORR activity and durability.

Identifying Few-Molecule Water Clusters with High Precision on Au(111) Surface

Revealing the nature of a hydrogen-bond network in water structures is one of the imperative objectives of science. With the use of a low-temperature scanning tunneling microscope, water clusters on a Au(111) surface were directly imaged with molecular resolution by a functionalized tip. The internal structures of the water clusters as well as the geometry variations with the increase of size were identified. In contrast to a buckled water hexamer predicted by previous theoretical calculations, our results present deterministic evidence for a flat configuration of water hexamers on Au(111), corroborated by density functional theory calculations with properly implemented van der Waals corrections. The consistency between the experimental observations and improved theoretical calculations not only renders the internal structures of absorbed water clusters unambiguously, but also directly manifests the crucial role of van der Waals interactions in constructing water-solid interfaces.

Decisive role of nuclear quantum effects on surface mediated water dissociation at finite temperature

Water molecules adsorbed on inorganic substrates play an important role in several technological applications. In the presence of light atoms in adsorbates, nuclear quantum effects (NQEs) influence the structural stability and the dynamical properties of these systems. In this work, we explore the impact of NQEs on the dissociation of water wires on stepped Pt(221) surfaces. By performing ab initio molecular dynamics simulations with van der Waals corrected density functional theory, we note that several competing minima for both intact and dissociated structures are accessible at finite temperatures, making it important to assess whether harmonic estimates of the quantum free energy are sufficient to determine the relative stability of the different states. We thus perform ab initio path integral molecular dynamics (PIMD) in order to calculate these contributions taking into account the conformational entropy and anharmonicities at finite temperatures. We propose that when adsorption is weak and NQEs on the substrate are negligible, PIMD simulations can be performed through a simple partition of the system, resulting in considerable computational savings. We then calculate the full contribution of NQEs to the free energies, including also anharmonic terms. We find that they result in an increase of up to 20% of the quantum contribution to the dissociation free energy compared with the harmonic estimates. We also find that the dissociation process has a negligible contribution from tunneling but is dominated by zero point energies, which can enhance the rate of dissociation by three orders of magnitude. Finally we highlight how both temperature and NQEs indirectly impact dipoles and the redistribution of electron density, causing work function changes of up to 0.4 eV with respect to static estimates. This quantitative determination of the change in the work function provides a possible approach to determine experimentally the most stable configurations of water oligomers on the stepped surfaces.

Regium-π bonds: An Unexplored Link between Noble Metal Nanoparticles and Aromatic Surfaces

The ability of metal clusters involving elements from group 11 (Ag, Cu, Au) to favorably interact with π systems of different size and electronic nature was evaluated at the PBE0-D3/def2-TZVPP//PBE0-D3/def2-TZVP level of theory. The M₉ clusters (M=Cu, Ag, Au) were used as σ-hole and σ-lump donors, and benzene, trifluorobenzene, and hexafluorobenzene as aromatic rings. In addition, the study was expanded to the analysis of extended π systems by using naphthalene and anthracene as well as their corresponding perfluorinated derivatives. Furthermore, Bader's theory of Atoms in Molecules as well as natural bonding orbital and spin-density calculations were used to further investigate and characterize the regium-π and σ-lump complexes described herein. Apparently, regium-π bonds have not previously been described in the literature and may be of great importance in the understanding of organocatalytic processes involving aromatic substrates as well as in the design of new materials based on this novel subclass of σ-hole bonding.

Lifshitz analysis of dispersion forces based on quantitative reflection electron energy loss spectroscopy

HYPOTHESIS

The energy loss experienced by an electron while moving through a solid is determined by the optical properties of the surrounding. Hence, quantitative Reflection Electron Energy Loss Spectroscopy (REELS) should allow the determination of optical data required for the calculation of Hamaker coefficients using Lifshitz theory. This approach might improve the accuracy of calculated Hamaker coefficients and should also enable to harness the unique capabilities of REELS to analyse nanostructured surfaces and thin-films with great spatial resolution and surface sensitivity.

EXPERIMENTS

REELS spectra of a survey of insulating polymers and of metal-like Ti_0.23Al_0.32N_0.44 (TiAlN) were measured, the complex dielectric functions determined and the corresponding Hamaker coefficients across vacuum and water calculated. The sensitivity of the quantification procedure towards typical systematic errors was investigated. For polystyrene the results were comparatively analysed using vacuum ultraviolet spectroscopy (VUV).

FINDINGS

The accuracy especially of the non-retarded vacuum Hamaker constants of the polymers was increased when compared to VUV reflectance spectroscopy due to the higher spectral range of REELS. Furthermore, a new correction procedure for the intricate case of unresolved inelastic losses in the REELS spectrum, such as encountered in the case of TiAlN, could be developed using spectroscopic ellipsometry as a complementary mean.

σ-Holes and σ-lumps direct the Lewis basic and acidic interactions of noble metal nanoparticles: introducing regium bonds

Using local DFT-based probes for electrostatic as well as charge transfer/polarization interactions, we are able to characterize Lewis basic and acidic sites on copper, silver and gold nanoparticles.

Ordered hydroxyls on Ca3Ru2O7(001)

As complex ternary perovskite-type oxides are increasingly used in solid oxide fuel cells, electrolysis and catalysis, it is desirable to obtain a better understanding of their surface chemical properties. Here we report a pronounced ordering of hydroxyls on the cleaved (001) surface of the Ruddlesden-Popper perovskite Ca₃Ru₂O₇ upon water adsorption at 105 K and subsequent annealing to room temperature. Density functional theory calculations predict the dissociative adsorption of a single water molecule (E_ads = 1.64 eV), forming an (OH)_ads group adsorbed in a Ca-Ca bridge site, with an H transferred to a neighboring surface oxygen atom, O_surf. Scanning tunneling microscopy images show a pronounced ordering of the hydroxyls with (2 × 1), c(2 × 6), (1 × 3), and (1 × 1) periodicity. The present work demonstrates the importance of octahedral rotation and tilt in perovskites, for influencing surface reactivity, which here induces the ordering of the observed OH overlayers.

As ternary perovskite-type oxides are increasingly used in fuel cells and catalysis, greater understanding of their surface chemical properties is required. Here the authors report a pronounced ordering of hydroxyls on the cleaved (001) surface of Ca₃Ru₂O₇ perovskite induced by O-octahedral rotation and tilt.

Comparing van der Waals DFT methods for water on NaCl(001) and MgO(001)

In this work, a range of van der Waals type density functionals are applied to the H₂O/NaCl(001) and H₂O/MgO(001) interface systems to explore the effect of an explicit dispersion treatment. The functionals we use are the self-consistent vdW functionals vdW-DF, vdW-DF2, optPBE-vdW, optB88-vdW, optB86b-vdW, and vdW-DF-cx, as well as the dispersion-corrected PBE-TS and PBE-D2 methods; they are all compared with the standard PBE functional. For both NaCl(001) and MgO(001), we find that the dispersion-flavoured functionals stabilize the water-surface interface by approximately 20%-40% compared to the PBE results. For NaCl(001), where the water molecules remain intact for all overlayers, the dominant contribution to the adsorption energy from "density functional theory dispersion" stems from the water-surface interactions rather than the water-water interactions. The optPBE-vdW and vdW-DF-cx functionals yield adsorption energies in good agreement with available experimental values for both NaCl and MgO. To probe the strengths of the perturbations of the adsorbed water molecules, we also calculated water dipole moments and found an increase up to 85% for water at the MgO(001) surface and 70% at the NaCl(001) surface, compared to the gas-phase dipole moment.

Molecular and dissociative adsorption of water and hydrogen sulfide at perfect and defective Cu(110) surfaces

The adsorption of H₂O and H₂S onto Cu(110) surfaces lead to considerably different geometries and surface reconstruction and relaxation mechanisms.

First-principles analysis of the stability of water on oxidised and reduced CuO(111) surfaces

A first-principles density functional theory investigation, with the inclusion of the Hubbard + U correction (PBE + U) on Cu-3d states, on the interaction of water with a CuO(111) surface.

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.

Capturing Solvation Effects at a Liquid/Nanoparticle Interface by Ab Initio Molecular Dynamics: Pt201 Immersed in Water

Solvation can substantially modify the adsorption properties of heterogeneous catalysts. Although essential for achieving realistic theoretical models, assessing such solvent effects over nanoparticles is challenging from a computational standpoint due to the complexity of those liquid/metal interfaces. This effect is investigated by ab initio molecular dynamics simulations at 350 K of a large platinum nanoparticle immersed in liquid water. The first solvation layer contains twice as much physisorbed water molecules above the terraces, than chemisorbed ones located only at edges and corners. The solvent stabilizes the binding energy of chemisorbates: 66% of the total gain comes from interactions with physisorbed molecules and 34% from the influence of bulk liquid.

Water at Interfaces

The interfaces of neat water and aqueous solutions play a prominent role in many technological processes and in the environment. Examples of aqueous interfaces are ultrathin water films that cover most hydrophilic surfaces under ambient relative humidities, the liquid/solid interface which drives many electrochemical reactions, and the liquid/vapor interface, which governs the uptake and release of trace gases by the oceans and cloud droplets. In this article we review some of the recent experimental and theoretical advances in our knowledge of the properties of aqueous interfaces and discuss open questions and gaps in our understanding.

Understanding and Designing the Gold-Bio Interface: Insights from Simulations

Gold nanoparticles (AuNPs) are an integral part of many exciting and novel biomedical applications, sparking the urgent need for a thorough understanding of the physicochemical interactions occurring between these inorganic materials, their functional layers, and the biological species they interact with. Computational approaches are instrumental in providing the necessary molecular insight into the structural and dynamic behavior of the Au-bio interface with spatial and temporal resolutions not yet achievable in the laboratory, and are able to facilitate a rational approach to AuNP design for specific applications. A perspective of the current successes and challenges associated with the multiscale computational treatment of Au-bio interfacial systems, from electronic structure calculations to force field methods, is provided to illustrate the links between different approaches and their relationship to experiment and applications.

Shape Evolution of Metal Nanoparticles in Water Vapor Environment

The structures of the metal nanoparticles are crucial for their catalytic activities. How to understand and even control the shape evolution of nanoparticles under reaction condition is a big challenge in heterogeneous catalysis. It has been proved that many reactive gases hold the capability of changing the structures and properties of metal nanoparticles. One interesting question is whether water vapor, such a ubiquitous environment, could induce the shape evolution of metal nanoparticles. So far this question has not received enough attention yet. In this work, we developed a model based on the density functional theory, the Wulff construction, and the Langmuir adsorption isotherm to explore the shape of metal nanoparticle at given temperature and water vapor pressure. By this model, we show clearly that water vapor could notably increase the fraction of (110) facets and decrease that of (111) facets for 3-8 nm Cu nanoparticles, which is perfectly consistent with the experimental observations. Further investigations indicate the water vapor has different effects on the different metal species (Cu, Au, Pt, and Pd). This work not only helps to understand the water vapor effect on the structures of metal nanoparticles but also proposes a simple but effective model to predict the shape of nanoparticles in certain environment.

First-principles study of the structure of water layers on flat and stepped Pb electrodes

On the basis of perodic density functional theory (DFT) calculations, we have addressed the geometric structures and electronic properties of water layers on flat and stepped Pb surfaces. In contrast to late d-band metals, on Pb(111) the energy minimum structure does not correspond to an ice-like hexagonal arrangement at a coverage of 2/3, but rather to a distorted structure at a coverage of 1 due to the larger lattice constant of Pb. At stepped Pb surfaces, the water layers are pinned at the step edge and form a complex network consisting of rectangles, pentagons and hexagons. The thermal stability of the water layers has been studied by using ab initio molecular dynamics simulations (AIMD) at a temperature of 140 K. Whereas the water layer on Pb(111) is already unstable at this temperature, the water layers on Pb(100), Pb(311), Pb(511) and Pb(711) exhibit a higher stability because of stronger water-water interactions. The vibrational spectra of the water layers at the stepped surfaces show a characteristic splitting into three modes in the O-H stretch region.

Identical Binding Energies and Work Functions for Distinct Adsorption Structures: Olympicenes on the Cu(111) Surface

Reliability is one of the major concerns and challenges in designing organic/inorganic interfaces for (opto)electronic applications. Even small structural differences for molecules on substrates can result in a significant variation in the interface functionality, due to the strong correlation between geometry, stability, and electronic structure. Here, we employed state-of-the-art first-principles calculations with van der Waals interactions, in combination with atomic force microscopy experiments, to explore the interaction mechanism for three structurally related olympicene molecules adsorbed on the Cu(111) surface. The substitution of a single atom in the olympicene molecule switches the nature of adsorption from predominantly physisorptive character [olympicene on Cu(111)], to an intermediate state [olympicene-derived ketone on Cu(111)], then to chemisorptive character [olympicene radical on Cu(111)]. Despite the remarkable difference in adsorption structures (by up to 0.9 Å in adsorption height) and different nature of bonding, the olympicene, its ketone, and its radical derivatives have essentially identical binding energies and work functions upon interaction with the metal substrate. Our findings suggest that the stability and work functions of molecular adsorbates could be rendered insensitive to their adsorption structures, which could be a useful property for (opto)electronic applications.

Neural network molecular dynamics simulations of solid–liquid interfaces: water at low-index copper surfaces

Molecular dynamics simulation of the water–copper interface have been carried out using high-dimensional neural network potential based on density functional theory.

Adsorption of water and ethanol on noble and transition-metal substrates: a density functional investigation within van der Waals corrections

We report an extensive density functional theory investigation of water and ethanol adsorption on several Cu-, Pt-, and Au-based substrates including substrates with low-coordinated sites due to intermixing of Pt–Cu and Pt–Au in the topmost surface.

Non-covalent adsorption of amino acid analogues on noble-metal nanoparticles: influence of edges and vertices

First-principles calculations on nanoscale-sized noble metal nanoparticles demonstrate that planes, edges and vertices show different noncovalent adsorption propensities depending on the adsorbate functional group.

Initial stages of water solvation of stepped platinum surfaces

Steps act as anchoring points for water adsorption and dominate water structures on stepped platinum surfaces.

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

How Does Water Wet a Surface?

The adsorption and reactions of water on surfaces has attracted great interest, as water is involved in many physical and chemical processes at interfaces. On metal surfaces, the adsorption energy of water is comparable to the hydrogen bond strength in water. Therefore, the delicate balance between the water-water and the water-metal interaction strength determines the stability of water structures. In such systems, kinetic effects play an important role and many metastable states can form with long lifetimes, such that the most stable state may not reached. This has led to difficulties in the theoretical prediction of water structures as well as to some controversial results. The direct imaging using scanning tunneling microscopy (STM) in ultrahigh vacuum at low temperatures offers a reliable means of understanding the local structure and reaction of water molecules, in particular when interpreted in conjunction with density functional theory calculations. In this Account, a selection of recent STM results on the water adsorption and dissociation on close-packed metal surfaces is reviewed, with a particular focus on Ru(0001). The Ru(0001) surface is one where water adsorbs intact in a metastable state at low temperatures and where partially dissociated layers are formed at temperatures above ∼150 K. First, we will describe the structure of intact water clusters starting with the monomer up to the monolayer. We show that icelike wetting layers do not occur on close-packed metal surfaces but instead hydrogen bonded layers in the form of a mixture of pentagonal, hexagonal, and heptagonal molecular rings are observed. Second, we will discuss the dissociation mechanism of water on Ru(0001). We demonstrate that water adsorption changes from dissociative to molecular as a function of the oxygen preadsorbed on Ru. Finally, we briefly review recent STM experiments on bulk ice (Ih and Ic) and water adsorption on insulating thin films. We conclude with an outlook illustrating the manipulation capabilities of STM in respect to probe the proton and hydrogen dynamics in water clusters.

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.

The gold/ampicillin interface at the atomic scale

In the fight against antibiotic resistance, gold nanoparticles (AuNP) with antibiotics grafted on their surfaces have been found to be potent agents. Ampicillin-conjugated AuNPs have been thus reported to overcome highly ampicillin-resistant bacteria. However, the structure at the atomic scale of these hybrid systems remains misunderstood. In this paper, the structure of the interface between an ampicillin molecule AMP and three flat gold facets Au(111), Au(110) and Au(100) has been investigated with numerical simulations (dispersion-corrected DFT). Adsorption energies, bond distances and electron densities indicate that the adsorption of AMP on these facets goes through multiple partially covalent bonding. The stability of the AuNP/AMP nanoconjugates is explained by large adsorption energies and their potential antibacterial activity is discussed on the basis of the constrained spatial orientation of the grafted antibiotic.

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.

Dispersion corrected RPBE studies of liquid water

The structure of liquid water has been addressed by ab initio molecular dynamics simulations based on density functional theory. Exchange-correlation effects have been described by the popular PBE and RPBE functionals within the generalized gradient approximation as these functionals also yield satisfactory results for metals which is important to model electrochemical interfaces from first principles. In addition, dispersive interactions are included by using dispersion-corrected schemes. It turns out that the dispersion-corrected RPBE functional reproduces liquid water properties quite well in contrast to the PBE functional. This is caused by the replacement of the over-estimated directional hydrogen-bonding in the PBE functional by non-directional dispersive interactions.

Molecular adsorption at Pt(111). How accurate are DFT functionals?

Recent DFT functionals including van der Waals contributions are benchmarked against micro calorimetric measurement for chemisorption on Pt(111). OptPBE-vdW and PBE-dDsC are found to be the most accurate, especially for unsaturated molecules where large errors can be committed with other functionals as PBE.

The accurate calculation of the band gap of liquid water by means of GW corrections applied to plane-wave density functional theory molecular dynamics simulations

Knowledge about the intrinsic electronic properties of water is imperative for understanding the behaviour of aqueous solutions that are used throughout biology, chemistry, physics, and industry.

Modeling adsorption and reactions of organic molecules at metal surfaces

CONSPECTUS: The understanding of adsorption and reactions of (large) organic molecules at metal surfaces plays an increasingly important role in modern surface science and technology. Such hybrid inorganic/organic systems (HIOS) are relevant for many applications in catalysis, light-emitting diodes, single-molecule junctions, molecular sensors and switches, and photovoltaics. Obviously, the predictive modeling and understanding of the structure and stability of such hybrid systems is an essential prerequisite for tuning their electronic properties and functions. At present, density-functional theory (DFT) is the most promising approach to study the structure, stability, and electronic properties of complex systems, because it can be applied to both molecules and solids comprising thousands of atoms. However, state-of-the-art approximations to DFT do not provide a consistent and reliable description for HIOS, which is largely due to two issues: (i) the self-interaction of the electrons with themselves arising from the Hartree term of the total energy that is not fully compensated in approximate exchange-correlation functionals, and (ii) the lack of long-range part of the ubiquitous van der Waals (vdW) interactions. The self-interaction errors sometimes lead to incorrect description of charge transfer and electronic level alignment in HIOS, although for molecules adsorbed on metals these effects will often cancel out in total energy differences. Regarding vdW interactions, several promising vdW-inclusive DFT-based methods have been recently demonstrated to yield remarkable accuracy for intermolecular interactions in the gas phase. However, the majority of these approaches neglect the nonlocal collective electron response in the vdW energy tail, an effect that is particularly strong in condensed phases and at interfaces between different materials. Here we show that the recently developed DFT+vdW(surf) method that accurately accounts for the collective electronic response effects enables reliable modeling of structure and stability for a broad class of organic molecules adsorbed on metal surfaces. This method was demonstrated to achieve quantitative accuracy for aromatic hydrocarbons (benzene, naphthalene, anthracene, and diindenoperylene), C60, and sulfur/oxygen-containing molecules (thiophene, NTCDA, and PTCDA) on close-packed and stepped metal surfaces, leading to an overall accuracy of 0.1 Å in adsorption heights and 0.1 eV in binding energies with respect to state-of-the-art experiments. An unexpected finding is that vdW interactions contribute more to the binding of strongly bound molecules on transition-metal surfaces than for molecules physisorbed on coinage metals. The accurate inclusion of vdW interactions also significantly improves tilting angles and adsorption heights for all the studied molecules, and can qualitatively change the potential-energy surface for adsorbed molecules with flexible functional groups. Activation barriers for molecular switches and reaction precursors are modified as well.

Molecular switches from benzene derivatives adsorbed on metal surfaces

Transient precursor states are often experimentally observed for molecules adsorbing on surfaces. However, such precursor states are typically rather short-lived, quickly yielding to more stable adsorption configurations. Here we employ first-principles calculations to systematically explore the interaction mechanism for benzene derivatives on metal surfaces, enabling us to selectively tune the stability and the barrier between two metastable adsorption states. In particular, in the case of the tetrachloropyrazine molecule, two equally stable adsorption states are identified with a moderate and conceivably reversible barrier between them. We address the feasibility of experimentally detecting the predicted bistable behaviour and discuss its potential usefulness in a molecular switch.

Short-lived precursors typically occur before molecules chemisorb on surfaces. Liu et al. predict that for benzene derivatives on metal surfaces, the precursors can be long-lived and the transition to chemisorption states can be reversible, which may be useful in molecular switch applications.

Efficient simulations of the aqueous bio-interface of graphitic nanostructures with a polarisable model

To fully harness the enormous potential offered by interfaces between graphitic nanostructures and biomolecules, detailed connections between adsorbed conformations and adsorption behaviour are needed. To elucidate these links, a key approach, in partnership with experimental techniques, is molecular simulation. For this, a force-field (FF) that can appropriately capture the relevant physics and chemistry of these complex bio-interfaces, while allowing extensive conformational sampling, and also supporting inter-operability with known biological FFs, is a pivotal requirement. Here, we present and apply such a force-field, GRAPPA, designed to work with the CHARMM FF. GRAPPA is an efficiently implemented polarisable force-field, informed by extensive plane-wave DFT calculations using the revPBE-vdW-DF functional. GRAPPA adequately recovers the spatial and orientational structuring of the aqueous interface of graphene and carbon nanotubes, compared with more sophisticated approaches. We apply GRAPPA to determine the free energy of adsorption for a range of amino acids, identifying Trp, Tyr and Arg to have the strongest binding affinity and Asp to be a weak binder. The GRAPPA FF can be readily incorporated into mainstream simulation packages, and will enable large-scale polarisable biointerfacial simulations at graphitic interfaces, that will aid the development of biomolecule-mediated, solution-based graphene processing and self-assembly strategies.