Abstract:One of the most widely used and computationally efficient models that accounts for London dispersion interactions within density functional theory (DFT) is the D3 dispersion correction model. In this work, we demonstrate that this model can induce the appearance of unphysical minima on the potential energy surface (PES) when the coordination number of atoms changes. Optimizing to these artifactual structures can lead to significant errors in determining the interaction energy between two molecules and in estimating the thermodynamic properties of the system. In several specific examples, such as Kuratowski-type H2-NiKur and H2-PdKur clusters, these local minima exhibited extremely high PES curvature, resulting in incorrect estimations of harmonic frequencies and significant overestimations of zero-point energy and enthalpy values. Although such erroneous behavior of the D3 model is relatively rare, it can occur across a wide range of chemical species, including molecules like the [Li(C6H6)]+ complex and the dispiro(acridan)-substituted pyracene (DSAP) molecule. Our analysis reveals that the root of the problem lies in the definition of the AB atomic-pair dependent C6AB coefficients in the D3 model. To address this issue, we propose a re-parameterization of the D3 model by introducing a modified C6AB functional form that now depends on the specific pair of considered atoms. This new model, termed D3-Smooth (or D3S for short), is designed to smooth out the PES associated with the dispersion correction. By doing so, we demonstrate that D3S eliminates unphysical local minima while maintaining the quite satisfactory accuracy of the parent D3 method in interaction energy benchmark sets. For example, the RMS difference between using the D3(BJ) correction to B3LYP and the D3S(BJ) correction across the large MGCDB84 dataset of nearly 5000 data points is only 0.12 kJ/mol. Similar results are obtained for every other D3-corrected functional tested. Consistent with this result, no significant improvement could be obtained for the B3LYP-D3S(0) correction by reoptimizing the damping function.

Removing the deadwood from DFT/MRCI wave functions: The p-DFT/MRCI method

A perturbative approximation to DFT/MRCI: DFT/MRCI(2)

A DFT/MRCI Hamiltonian Parameterized Using Only Ab Initio Data: I. Valence Excited States

Calculation of quasi-diabatic states within the DFT/MRCI(2) framework: the QD-DFT/MRCI(2) method

Externally-Contracted Multireference Configuration Interaction Method Using a DMRG Reference Wave Function

Excited State Structure and Minimum Energy Conical Intersection Optimization Using DFT/MRCI

Minimum Energy Conical Intersection Optimization Using DFT/MRCI(2)

Time-Reversal Symmetry in RDMFT and pCCD with Complex-Valued Orbitals

Cluster decomposition of full configuration interaction wave functions: a tool for chemical interpretation of systems with strong correlation

Combining density based dynamical correlation with a reduced density matrix strong correlation description

Reduced density matrix functional theory from an ab initio seniority-zero wave function: Exact and approximate formulations along adiabatic connection paths

A DFT/MRCI Hamiltonian Parameterized Using Only Ab Initio Data: II. Core-Excited States

Enhancing Reduced Density Matrix Functional Theory Calculations by Coupling Orbital and Occupation Optimizations

Density Functional Theory Transformed into a One-electron Reduced Density Matrix Functional Theory for the Capture of Static Correlation

Renormalized-Residue-Based Multireference Configuration Interaction Method for Strongly Correlated Systems

Improved on-the-Fly MCTDH Simulations with Many-Body-Potential Tensor Decomposition and Projection Diabatization

A New View on Density Corrected DFT: Can One Get a Better Answer for a Good Reason?

Diagrammatic Multiplet-Sum Method (MSM) Density-Functional Theory (DFT): Investigation of the Transferability of Integrals in "Simple" DFT-Based Approaches to Multi-Determinantal Problems

RADE: A Reduced Approach to Density-Functional Expansion

Smoother Semiclassical Dispersion for Density Functional Theory via D3S: Understanding and Addressing Unphysical Minima in the D3 Dispersion Correction Model

Density matrices in full configuration interaction quantum Monte Carlo: Excited states, transition dipole moments and parallel distribution