No Arabic abstract
The cluster-in-molecule (CIM) local correlation approach with an accurate distant pair correlation energy correction is presented. For large systems, the inclusion of distant pair correlation energies is essential for the accurate predictions of absolute correlation energies and relative energies. Here we propose a simple and efficient scheme for evaluating the distant pair correlation energy correction. The corrections can be readily extracted from electron correlation calculations of clusters with almost no additional effort. Benchmark calculations show that the improved CIM approach can recover more than 99.97% of the conventional correlation energy. By combining the CIM approach with the domain based local pair natural orbital (DLPNO) local correlation approach, we have provided accurate binding energies at the CIM-DLPNO-CCSD(T) level for a test set consisting of eight weakly bound complexes ranging in size from 200 to 1027 atoms. With these results as the reference data, the accuracy and applicability of other electron correlation methods and a few density functional methods for large systems have been assessed.
A recently proposed local self-interaction correction (LSIC) method [Zope textit{et al.} J. Chem. Phys., 2019,{bf 151}, 214108] when applied to the simplest local density approximation provides significant improvement over standard Perdew-Zunger SIC (PZSIC) for both equilibrium properties such as total or atomization energies as well as properties involving stretched bond such as barrier heights. The method uses an iso-orbital indicator to identify the single-electron regions. To demonstrate the LSIC method, Zope textit{et al.} used the ratio $z_sigma$ of von Weizsacker $tau_sigma^W$ and total kinetic energy densities $tau_sigma$, ($z_sigma = tau_sigma^W/tau_sigma$) as a scaling factor to scale the self-interaction correction. The present work further explores the LSIC method using a simpler scaling factor as a ratio of orbital and spin densities in place of the ratio of kinetic energy densities. We compute a wide array of both, equilibrium and non-equilibrium properties using the LSIC and orbital scaling methods using this simple scaling factor and compare them with previously reported results. Our study shows that the present results with simple scaling factor are comparable to those obtained by LSIC($z_sigma$) for most properties but have slightly larger errors. We furthermore study the binding energies of small water clusters using both the scaling factors. Our results show that LSIC with $z_{sigma}$ has limitation in predicting the binding energies of weakly bonded system due to the inability of $z_{sigma}$ to distinguish weakly bonded region from slowly varying density region. LSIC when used with density ratio as a scaling factor, on the other hand, provides good description of water cluster binding energies, thus highlighting the appropriate choice of iso-orbital indicator.
A full coupled-cluster expansion suitable for sparse algebraic operations is developed by expanding the commutators of the Baker-Campbell-Hausdorff series explicitly for cluster operators in binary representations. A full coupled-cluster reduction that is capable of providing very accurate solutions of the many-body Schrodinger equation is then initiated employing screenings to the projection manifold and commutator operations. The projection manifold is iteratively updated through the single commutators $leftlangle kappa right| [hat H,hat T]left| 0 rightrangle$ comprised of the primary clusters $hat T_{lambda}$ with substantial contribution to the connectivity. The operation of the commutators is further reduced by introducing a correction, taking into account the so-called exclusion principle violating terms, that provides fast and near-variational convergence in many cases.
Potential energy surfaces of the hydrogen molecular ion H$_2^+$ in the Born-Oppenheimer approximation are computed by means of the Riccati-Pade method (RPM). The convergence properties of the method are analyzed for different states. The equilibrium internuclear distance, as well as the corresponding electronic plus nuclear energy, and the associated separation constants, are computed to 40 digits of accuracy for several bound states. For the ground state the same parameters are computed with more than 100 digits of accuracy. Additional benchmark values of the electronic energy at different internuclear distances are given for several additional states. The software implementation of the RPM is given under a free software license. The results obtained in the present work are the most accurate available so far, and further additional benchmarks are made available through the software provided.
We present a near-linear scaling formulation of the explicitly-correlated coupled-cluster singles and doubles with perturbative triples method (CCSD(T)$_{overline{text{F12}}}$) for high-spin states of open-shell species. The approach is based on the conventional open-shell CCSD formalism [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] utilizing the domain local pair-natural orbitals (DLPNO) framework. The use of spin-independent set of pair-natural orbitals ensures exact agreement with the closed-shell formalism reported previously, with only marginally impact on the cost (e.g. the open-shell formalism is only 1.5 times slower than the closed-shell counterpart for the $text{C}_text{160}text{H}_{text{322}}$ n-alkane, with the measured size complexity of $approx1.2$). Evaluation of coupled-cluster energies near the complete-basis-set (CBS) limit for open-shell systems with more than 550 atoms and 5000 basis functions is feasible on a single multi-core computer in less than 3 days. The aug-cc-pVTZ DLPNO-CCSD(T)$_{overline{text{F12}}}$ contribution to the heat of formation for the 50 largest molecules among the 348 core combustion species benchmark set [J. Klippenstein et al., J. Phys. Chem. A 121, 6580 (2017)] had root-mean-square deviation (RMSD) from the extrapolated CBS CCSD(T) reference values of 0.3 kcal/mol. For a more challenging set of 50 reactions involving small closed- and open-shell molecules [G. Knizia et al., J. Chem. Phys. 130, 054104 (2009)] the aug-cc-pVQ(+d)Z DLPNO-CCSD(T)$_{overline{text{F12}}}$ yielded a RMSD of $sim$0.4 kcal/mol with respect to the CBS CCSD(T) estimate.
Aiming to remedy the incorrect asymptotic behavior of conventional semilocal exchange-correlation (XC) density functionals for finite systems, we propose an asymptotic correction scheme, wherein an exchange density functional whose functional derivative has the correct (-1/r) asymptote can be directly added to any semilocal density functional. In contrast to semilocal approximations, our resulting exchange kernel in reciprocal space exhibits the desirable singularity of the type O(-1/q^2) as q -> 0, which is a necessary feature for describing the excitonic effects in non-metallic solids. By applying this scheme to a popular semilocal density functional, PBE [J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)], the predictions of the properties that are sensitive to the asymptote are significantly improved, while the predictions of the properties that are insensitive to the asymptote remain essentially the same as PBE. Relative to the popular model XC potential scheme, our scheme is significantly superior for ground-state energies and related properties. In addition, without loss of accuracy, two closely related schemes are developed for the efficient treatment of large systems.