I-Chun Lin

Structure and Dynamics of Liquid Water from ab Initio Molecular Dynamics—Comparison of BLYP, PBE, and revPBE Density Functionals with and without van der Waals Corrections

We investigate the accuracy provided by different treatments of the exchange and correlation effects, in particular the London dispersion forces, on the properties of liquid water using ab initio molecular dynamics simulations with density functional theory. The lack of London dispersion forces in generalized gradient approximations (GGAs) is remedied by means of dispersion-corrected atom-centered potentials (DCACPs) or damped atom-pairwise dispersion corrections of the C6R(-6) form. We compare results from simulations using GGA density functionals (BLYP, PBE, and revPBE) with data from their van der Waals (vdW) corrected counterparts. As pointed out previously, all vdW-corrected BLYP simulations give rise to highly mobile water whose softened structure is closer to experimental data than the one predicted by the bare BLYP functional. Including vdW interactions in the PBE functional, on the other hand, has little influence on both structural and dynamical properties of water. Augmenting the revPBE functional with either damped atom-pairwise dispersion corrections or DCACP evokes opposite behaviors. The former further softens the already under-structured revPBE water, whereas the latter makes it more glassy. These results demonstrate the delicacy needed in describing weak interactions in molecular liquids.

Importance of van der Waals interactions in liquid water.

We present ab initio molecular dynamics studies on liquid water using density functional theory in conjunction with either dispersion-corrected atom-centered potentials or empirical van der Waals corrections. Our results show that improving the description of van der Waals interactions in DFT-GGA leads to a softening of liquid waters structure with higher mobility. The results obtained with dispersion-corrected atom-centered potentials are especially encouraging. In particular, the radial distribution functions are in better agreement with experiment, and the self-diffusion coefficient increases by more than three-fold compared with the one predicted by the BLYP functional. This work demonstrates that van der Waals interactions are essential in fine-tuning both structural and dynamical properties of liquid water.

Weakly bonded complexes of aliphatic and aromatic carbon compounds described with dispersion corrected density functional theory

Interaction energies and structural properties of van der Waals complexes of aliphatic hydrocarbons molecules and crystals of aromatic hydrocarbon compounds are studied using density functional theory augmented with dispersion corrected atom centered potentials (DCACPs). We compare the performance of two sets of DCACPs, (a) DCACP-MP2, a correction for carbon only, generated using MP2 reference data and a penalty functional that includes only equilibrium properties and (b) DCACP-CCSD(T), a set that has been calibrated against CCSD(T) reference data using a more elaborate penalty functional that explicitly takes into account some long-range properties and uses DCACP corrections for hydrogen and carbon atoms. The agreement between our results and high level ab initio or experimental data illustrates the transferability of the DCACP scheme for the gas and condensed phase as well as for different hybridization states of carbon. The typical error of binding energies for gas-phase dimers amounts to 0.3 kcal/mol. This work demonstrates that only one DCACP per element is sufficient to correct for weak interactions in a large variety of systems, irrespective of the hybridization state.

Describing weak interactions of biomolecules with dispersion-corrected density functional theory

Interaction energies of the biomolecules in the JSCH-2005 database are calculated with density functional theory using the exchange-correlation functional BLYP augmented with dispersion-corrected atom-centered potentials (DCACPs). The results are in excellent agreement with extrapolated CCSD(T) complete basis set limit references with unsigned mean errors of less than 1.6 kcal mol(-1). Geometry optimisations all reach stable configurations that are close to the MP2-optimised reference geometries.

Hydrogen Bonding Described Using Dispersion-Corrected Density Functional Theory

In recent works, dispersion-corrected atom-centered potentials (DCACPs) were developed as a method to account for long-range dispersion forces between molecules in density functional theory calculations within the generalized gradient approximation (GGA). Here, we test the ability of DCACPs to improve the GGA treatment of hydrogen-bonded systems. We assessed both BLYP and dispersion-corrected BLYP with respect to benchmark calculations for the hydrogen bond lengths and binding energies of 20 complexes containing the elements C, H, N, O, and S. Benchmark data included geometries calculated using MP2 and CCSD(T) and binding energies using W2, W1, CBS-QB3, and other CCSD(T) extrapolation schemes. With respect to benchmark methods, dispersion-corrected BLYP exhibited a mean signed error of 0.010 A in the hydrogen bond length and a mean relative error of 5.1% in the hydrogen bond binding energy. By comparison, uncorrected BLYP exhibited error statistics of 0.036 A and 15.9%, respectively. We conclude that DCACPs robustly improve the BLYP description of hydrogen-bonded systems at small additional computational cost. New benchmark geometries (MP2/aug-cc-pVTZ) and new benchmark binding energies (W1) are presented for seven complexes, and the remaining benchmark data were taken from previous literature.

Accurate DFT Descriptions for Weak Interactions of Molecules Containing Sulfur

Dispersion corrected atom centered potentials (DCACPs) have been shown to significantly improve the density functional theory (DFT) description of weak interactions. In this work, we have calibrated a DCACP for sulfur in combination with the widely used Generalized Gradient Approximation (GGA) BLYP, thereby augmenting the existing library of DCACPs for the first- and second-row elements H, C, N, O, and rare gases. Three weakly bound complexes as well as elemental (orthorhombic) sulfur are used as test cases to evaluate the transferability of the DCACP to different chemical environments. It is found that the sulfur DCACP systematically improves the agreement of DFT-calculated weak interactions with respect to MP2 and CCSD(T) level results.

Dispersion Corrected Atom-Centered Potentials for Phosphorus.

Dispersion-corrected atom-centered potentials (DCACPs) for the element phosphorus were generated and tested for the BLYP, BP, and PBE generalized gradient approximations of the exchange-correlation functional. The accuracy and transferability of the DCACPs were tested by evaluating the interaction energy of different weakly bound molecular systems (P2, PH3, and PN dimers). These results were compared to reference CCSD(T) calculations and standard density functional theory (DFT). The DCACP were also tested in the case of condensed phase systems. Specifically, the density and cohesive energies of β-white and black phosphorus were estimated and compared to available experimental data. Our results show an overall strong improvement both at the qualitative and quantitative level, with respect to uncorrected generalized gradient approximation DFT results for all three functionals. In particular, BLYP-corrected results show the maximal transferability, reporting for all systems a deviation from CCSD(T) results of less than 1% in the predicted binding energies.