Biswajit Santra

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.

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.

Advanced capabilities for materials modelling with Quantum ESPRESSO

Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudo-potential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement theirs ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.Quantum EXPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum EXPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water

In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local van der Waals/dispersion (vdW) interactions, via a fully self-consistent density-dependent dispersion correction, and an approximate treatment of nuclear quantum effects, via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H2O)128 significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, SOO(Q), and corresponding oxygen-oxygen radial distribution function, gOO(r), that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of a hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water also yields higher-order correlation functions, such as the oxygen-oxygen-oxygen triplet angular distribution, POOO(θ), and therefore the degree of local tetrahedrality, as well as electrostatic properties, such as the effective molecular dipole moment, that are in much better agreement with experiment.

Structure and energetics of benzene adsorbed on transition-metal surfaces: density-functional theory with van der Waals interactions including collective substrate response

The adsorption of benzene on metal surfaces is an important bench- mark system for hybrid inorganic/organic interfaces. The reliable determination of the interface geometry and binding energy presents a significant challenge for both theory and experiment. Using the Perdew-Burke-Ernzerhof (PBE), PBE+vdW (van der Waals) and the recently developed PBE+vdW surf (density- functional theory with vdW interactions that include the collective electronic response of the substrate) methods, we calculated the structures and energet- ics for benzene on transition-metal surfaces: Cu, Ag, Au, Pd, Pt, Rh and Ir. Our calculations demonstrate that vdW interactions increase the binding en- ergy by more than 0.70eV for physisorbed systems (Cu, Ag and Au) and by an even larger amount for strongly bound systems (Pd, Pt, Rh and Ir). The collective response of the substrate electrons captured via the vdW surf method plays a significant role for most substrates, shortening the equilibrium dis- tance by 0.25A for Cu and decreasing the binding energy by 0.27eV for Rh. The reliability of our results is assessed by comparison with calculations using the random-phase approximation including renormalized single excitations,

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 [B. Santra, J. Klimeš, D. Alfè, A. Tkatchenko, B. Slater, A. Michaelides, R. Car, and M. Scheffler, Phys. Rev. Lett. 107, 185701 (2011)] which showed that vdW forces become increasingly important at high pressures, we find here that all vdW inclusive methods considered improve the relative energies and transition pressures of the high-pressure ice phases compared to those obtained with semi-local or hybrid xc functionals. However, we also find that significant discrepancies between experiment and the vdW inclusive approaches remain in the cohesive properties of the various phases, causing certain phases to be absent from the phase diagram. Therefore, room for improvement in the description of water at ambient and high pressures remains and we suggest that because of the stern test the high pressure ice phases pose they should be used in future benchmark studies of simulation methods for water.

Coupled cluster benchmarks of water monomers and dimers extracted from density-functional theory liquid water: The importance of monomer deformations

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

Local structure analysis in ab initio liquid water

Within the framework of density functional theory, the inclusion of exact exchange and non-local van der Waals/dispersion (vdW) interactions is crucial for predicting a microscopic structure of ambient liquid water that quantitatively agrees with experiment. In this work, we have used the local structure index (LSI) order parameter to analyse the local structure in such highly accurate ab initio liquid water. At ambient conditions, the LSI probability distribution, P(I ), was unimodal with most water molecules characterised by more disordered high-density-like local environments. With thermal excitations removed, the resultant bimodal P(I ) in the inherent potential energy surface (IPES) exhibited a 3:1 ratio between high-density- and low-density-like molecules, with the latter forming small connected clusters amid the predominant population. By considering the spatial correlations and hydrogen bond network topologies among water molecules with the same LSI identities, we demonstrate that the signatures of the experimentally observed low- and high-density amorphous phases of ice are present in the IPES of ambient liquid water. Analysis of the LSI autocorrelation function uncovered a persistence time of ∼ 4 ps – a finding consistent with the fact that natural thermal fluctuations are responsible for transitions between these distinct yet transient local aqueous environments in ambient liquid water.