Assessment of DFT for endohedral complexes' dipole moment: PNO-LCCSD-F12 as a reference method

State-of-the-art local correlation methods enable affordable gold standard quantum chemistry for up to hundreds of atoms

In this feature, we review the current capabilities of local electron correlation methods up to the coupled cluster model with single, double, and perturbative triple excitations [CCSD(T)], which is a gold standard in quantum chemistry. The main computational aspects of the local method types are assessed from the perspective of applications, but the focus is kept on how to achieve chemical accuracy (i.e., <1 kcal mol-1 uncertainty), as well as on the broad scope of chemical problems made accessible. The performance of state-of-the-art methods is also compared, including the most employed DLPNO and, in particular, our local natural orbital (LNO) CCSD(T) approach. The high accuracy and efficiency of the LNO method makes chemically accurate CCSD(T) computations accessible for molecules of hundreds of atoms with resources affordable to a broad computational community (days on a single CPU and 10-100 GB of memory). Recent developments in LNO-CCSD(T) enable systematic convergence and robust error estimates even for systems of complicated electronic structure or larger size (up to 1000 atoms). The predictive power of current local CCSD(T) methods, usually at about 1-2 order of magnitude higher cost than hybrid density functional theory (DFT), has become outstanding on the palette of computational chemistry applicable for molecules of practical interest. We also review more than 50 LNO-based and other advanced local-CCSD(T) applications for realistic, large systems across molecular interactions as well as main group, transition metal, bio-, and surface chemistry. The examples show that properly executed local-CCSD(T) can contribute to binding, reaction equilibrium, rate constants, etc. which are able to match measurements within the error estimates. These applications demonstrate that modern, open-access, and broadly affordable local methods, such as LNO-CCSD(T), already enable predictive computations and atomistic insight for complicated, real-life molecular processes in realistic environments.

DFT and spatial confinement: a benchmark study on the structural and electrical properties of hydrogen bonded complexes

An extended set of 37 exchange correlation functionals, representing different DFT approximations, has been evaluated on a difficult playground represented by the dipole moment (μ_z), polarizability (α_zz), first hyperpolarizability (β_zzz), and the corresponding interaction-induced electrical properties (Δμ_z, Δα_zz, Δβ_zzz) of spatially confined hydrogen bonded (HB) dimers. A two-dimensional harmonic oscillator potential was used to exert the effect of spatial restriction. The performance of DFT methods in predicting hydrogen bond lengths in the studied molecular complexes upon confinement has also been examined. The data determined using a high-level CCSD(T) method serve as a reference. The conducted analyses allow us to conclude that methods rooted in DFT constitute a precise tool for the calculation of μ_z and α_zz as well as Δμ_z and Δα_zz, as most of the tested functionals provide results affected by rather small relative errors. On the other hand, an accurate description of the nonlinear optical response of the studied HB systems remains a great challenge for most of the analyzed DFT functionals, both in vacuum and in the presence of an analytical confining potential. Some of the tested DFT methods are found to be prone to catastrophic failure in the prediction of β_zzz as well as Δβ_zzz. The obtained results indicate that there is no great chasm in performance between functionals belonging to different DFT approximations or functionals including different amount of Hartree-Fock exchange when the values of dipole moment and first hyperpolarizability as well as the corresponding interaction-induced electrical properties are considered. However, a higher fraction of Hartree-Fock exchange improves the quality of predictions of α_zz and Δα_zz. Additionally, it has been shown that only three functionals from the examined set, namely B2PLYP, B3LYP and ωB97X-D, provide highly accurate structural parameters for the investigated systems. Of significant importance is the conclusion that the ωB97X-D functional, representing a modern and highly parametrized range-separated hybrid, demonstrates the most coherent behavior, showing rather small deviations from the reference data in the case of μ_z, α_zz, Δμ_z and Δα_zz as well as the structural parameters of the studied HB dimers. Moreover, our results indicate that the presence of spatial confinement has a rather small effect on the performance of DFT methods.