Angelos Michaelides

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.

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.

On the accuracy of density-functional theory exchange-correlation functionals for H bonds in small water clusters. II. The water hexamer and van der Waals interactions

Second order Møller-Plesset perturbation theory at the complete basis set limit and diffusion quantum Monte Carlo are used to examine several low energy isomers of the water hexamer. Both approaches predict the so-called prism to be the lowest energy isomer, followed by cage, book, and cyclic isomers. The energies of the four isomers are very similar, all being within 10-15 meV/H(2)O. These reference data are then used to evaluate the performance of several density-functional theory exchange-correlation (xc) functionals. A subset of the xc functionals tested for smaller water clusters [I. Santra et al., J. Chem. Phys. 127, 184104 (2007)] has been considered. While certain functionals do a reasonable job at predicting the absolute dissociation energies of the various isomers (coming within 10-20 meV/H(2)O), none predict the correct energetic ordering of the four isomers nor does any predict the correct low total energy isomer. All xc functionals tested either predict the book or cyclic isomers to have the largest dissociation energies. A many-body decomposition of the total interaction energies within the hexamers leads to the conclusion that the failure lies in the poor description of van der Waals (dispersion) forces in the xc functionals considered. It is shown that the addition of an empirical pairwise (attractive) C(6)R(-6) correction to certain functionals allows for an improved energetic ordering of the hexamers. The relevance of these results to density-functional simulations of liquid water is also briefly discussed.

On the accuracy of density-functional theory exchange-correlation functionals for H bonds in small water clusters: Benchmarks approaching the complete basis set limit

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

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 molecules in the gas phase.

Perspective: How good is DFT for water?

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

A density functional theory study of hydroxyl and the intermediate in the water formation reaction on Pt

Density functional theory has been used to study the adsorption of hydroxyl at low and high coverages and also to investigate the nature of the intermediate in the H2O formation reaction on Pt(111). At low coverages [1/9 of a monolayer (ML) to 1/3 ML] OH binds preferentially at bridge and top sites with a chemisorption energy of ∼2.25 eV. At high coverages (1/2 ML to 1 ML) H bonding between adjacent hydroxyls causes: (i) an enhancement in OH chemisorption energy by about 15%; (ii) a strong preference for OH adsorption at top sites; and (iii) the formation of OH networks. The activation energy for the diffusion of isolated OH groups along close packed rows of Pt atoms is 0.1 eV. This low barrier coupled with H bonding between neighboring OH groups indicates that hydroxyls are susceptible to island formation at low coverages. Pure OH as well as coadsorbed OH and H can be ruled out as the observed low temperature intermediate in the water formation reaction. Instead we suggest that the intermediate consists ...

A one-dimensional ice structure built from pentagons

Heterogeneous ice nucleation has a key role in fields as diverse as atmospheric chemistry and biology. Ice nucleation on metal surfaces affords an opportunity to watch this process unfold at the molecular scale on a well-defined, planar interface. A common feature of structural models for such films is that they are built from hexagonal arrangements of molecules. Here we show, through a combination of scanning tunnelling microscopy, infrared spectroscopy and density-functional theory, that about 1-nm-wide ice chains that nucleate on Cu(110) are not built from hexagons, but instead are built from a face-sharing arrangement of water pentagons. The pentagon structure is favoured over others because it maximizes the water-metal bonding while maintaining a strong hydrogen-bonding network. It reveals an unanticipated structural adaptability of water-ice films, demonstrating that the presence of the substrate can be sufficient to favour non-hexagonal structural units.

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 nonlocal 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. Because of the greater polarizability of the substrate metal atoms, nonlocal 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 icelike wetting layers are thermodynamically stable.

Improved description of soft layered materials with van der Waals density functional theory

The accurate description of van der Waals forces within density functional theory is currently one of the most active areas of research in computational physics and chemistry. Here we report results on the structural and energetic properties of graphite and hexagonal boron nitride, two layered materials where interlayer binding is dominated by van der Waals forces. Results from several density functionals are reported, including the optimized Becke88 van der Waals (optB88-vdW) and the optimized PBE van der Waals (optPBE-vdW) (Klimeš et al 2010 J. Phys.: Condens. Matter 22 022201) functionals. Where comparison to experiment and higher-level theory is possible, the results obtained from the two new van der Waals density functionals are in good agreement. An analysis of the physical nature of the interlayer binding in both graphite and hexagonal boron nitride is also reported.