One atom or molecule binds to another through various types of bond, the strengths of which range from several meV to several eV. Although some computational methods can provide accurate descriptions of all bond types, those methods are not efficient enough for many studies (for example, large systems, ab initio molecular dynamics and high-throughput searches for functional materials). Here, we show that the recently developed non-empirical strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) within the density functional theory framework predicts accurate geometries and energies of diversely bonded molecules and materials (including covalent, metallic, ionic, hydrogen and van der Waals bonds). This represents a significant improvement at comparable efficiency over its predecessors, the GGAs that currently dominate materials computation. Often, SCAN matches or improves on the accuracy of a computationally expensive hybrid functional, at almost-GGA cost. SCAN is therefore expected to have a broad impact on chemistry and materials science.

The methods of density-functional perturbation theory may be used to calculate various physical response properties of insulating crystals including elastic, dielectric, Born charge, and piezoelectric tensors. These and other important tensors may be defined as second derivatives of the total energy with respect to atomic-displacement, electric-field, or strain perturbations, or as mixed derivatives with respect to two of these perturbations. The resulting tensor quantities tend to be coupled in complex ways in polar crystals, giving rise to a variety of variant definitions. For example, it is generally necessary to distinguish between elastic tensors defined under different electrostatic boundary conditions, and between dielectric tensors defined under different elastic boundary conditions. Here, we describe an approach for computing all of these various response tensors in a unified and systematic fashion. Applications are presented for two materials, wurtzite ZnO and rhombohedral BaTiO3, at zero temperature. 62.20.Dc, 77.22.Ch, 77.65.Bn, 71.15.Mb

Water is of the utmost importance for life and technology. However, a genuinely predictive ab initio model of water has eluded scientists. We demonstrate that a fully ab initio approach, relying on the strongly constrained and appropriately normed (SCAN) density functional, provides such a description of water. SCAN accurately describes the balance among covalent bonds, hydrogen bonds, and van der Waals interactions that dictates the structure and dynamics of liquid water. Notably, SCAN captures the density difference between water and ice Ih at ambient conditions, as well as many important structural, electronic, and dynamic properties of liquid water. These successful predictions of the versatile SCAN functional open the gates to study complex processes in aqueous phase chemistry and the interactions of water with other materials in an efficient, accurate, and predictive, ab initio manner.water | ab initio theory | hydrogen bonding | density functional theory | molecular dynamics W ater is arguably the most important molecule for life and is involved in almost all biological processes. Without water, life, as we know it, would not exist, earning water the pseudonym "matrix of life," among others (1). Despite the apparent simplicity of an H2O molecule, water in the condensed phase displays a variety of anomalous properties that originate from its complex structure. In an ideal arrangement, water molecules form a tetrahedral network of hydrogen (H) bonds with each vertex being occupied by a water molecule. This tetrahedral network is realized in the solid phase ice Ih, but thermal fluctuations disrupt the H-bond network in the liquid state, with the network fluctuating on picosecond to nanosecond timescales. Due to the complexity of the H-bond network and its competition with thermal fluctuations, a precise molecular-level understanding of the structure of liquid water remains elusive. Major challenges lie in unambiguously capturing the atomic-scale fluctuations in water experimentally. Current approaches such as time-resolved spectroscopy (2, 3) and diffraction measurements (4, 5) may be able to resolve changes on picosecond timescales but rely on interpretation through models, which often cannot describe all of the details of liquid water with quantitative accuracy. Not surprisingly, the nature of the H-bond network in liquid water continues to be at the center of scientific debate, and advances in both experiment and theory are needed, especially with regard to quantitative modeling of aqueous phase chemistry.Ab initio molecular dynamics (AIMD) simulation (6) is an ideal approach for modeling the condensed phases of water across the phase diagram and aqueous phase chemistry using quantum mechanical principles (7-11), although for some applications, such as the study of liquid vapor phase equilibria (12), Monte Carlo methods are better suited. In particular, Kohn-Sham density functional theory (DFT) (13)-used to model the system in its electronic ground state-provides an efficient framework that enables the si...

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 vdW/dispersion interactions, via a fully self-consistent density-dependent dispersion correction, and approximate nuclear quantum effects (aNQE), 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, vdW, and aNQE as resulting from a large-scale AIMD simulation of (H2O)128 at the PBE0+vdW level of theory, 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 as 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 the PBE0 hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of nonlocal vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water obtained at the PBE0+vdW+aNQE level of theory 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 also in much better agreement with experiment.

The spinel cobalt oxide Co 3 O 4 is a magnetic semiconductor containing cobalt ions in Co 2+ and Co 3+ oxidation states. We have studied the electronic, magnetic and bonding properties of Co 3 O 4 using density functional theory (DFT) at the Generalized Gradient Approximation (GGA), GGA+U, and PBE0 hybrid functional levels. The GGA correctly predicts Co 3 O 4 to be a semiconductor, but severely underestimates the band gap. The GGA+U band gap (1.96 eV) agrees well with the available experimental value (~ 1.6 eV), whereas the band gap obtained using the PBE0 hybrid functional (3.42 eV) is strongly overestimated. All the employed exchange-correlation functionals predict 3 unpaired d electrons on the Co 2+ ions, in agreement with crystal field theory, but the values of the magnetic moments given by GGA+U and PBE0 are in closer agreement with the experiment than the GGA value, indicating a better description of the cobalt localized d states. Bonding properties are studied by means of Maximally Localized Wannier Functions (MLWFs). We find d-type MLWFs on the cobalt ions, as well as Wannier functions with the character of sp 3 d bonds between cobalt and oxygen ions.Such hybridized bonding states indicate the presence of a small covalent component in the primarily ionic bonding mechanism of this compound.