Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA well-scaling natural orbital theory
77
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree-Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the one-particle density matrix.
rate research
Read More
We introduce a novel energy functional for ground-state electronic-structure calculations. Its fundamental variables are the natural spin-orbitals of the implied singlet many-body wave function and their joint occupation probabilities. The functional derives from a sequence of controlled approximations to the two-particle density matrix. Algebraic scaling of computational cost with electron number is obtainable in general, and Hartree-Fock scaling in the seniority-zero version of the theory. Results obtained with the latter version for saturated small molecular systems are compared with those of highly-accurate quantum-chemical computations. The numerical results are variational, capturing most of the correlation energy from equilibrium to dissociation. Their accuracy is considerably greater than that obtainable with current density-functional theory approximations and with current functionals of the one-particle density matrix only.
Semi-local approximations to the density functional for the exchange-correlation energy of a many-electron system necessarily fail for lobed one-electron densities, including not only the familiar stretched densities but also the less familiar but closely-related noded ones. The Perdew-Zunger (PZ) self-interaction correction (SIC) to a semi-local approximation makes that approximation exact for all one-electron ground- or excited-state densities and accurate for stretched bonds. When the minimization of the PZ total energy is made over real localized orbitals, the orbital densities can be noded, leading to energy errors in many-electron systems. Minimization over complex localized orbitals yields nodeless orbital densities, which reduce but typically do not eliminate the SIC errors of atomization energies. Other errors of PZ SIC remain, attributable to the loss of the exact constraints and appropriate norms that the semi-local approximations satisfy, and suggesting the need for a generalized SIC. These conclusions are supported by calculations for one-electron densities, and for many-electron molecules. While PZ SIC raises and improves the energy barriers of standard generalized gradient approximations (GGAs) and meta-GGAs, it reduces and often worsens the atomization energies of molecules. Thus PZ SIC raises the energy more as the nodality of the valence localized orbitals increases from atoms to molecules to transition states. PZ SIC is applied here in particular to the SCAN meta-GGA, for which the correlation part is already self-interaction-free. That property makes SCAN a natural first candidate for a generalized SIC.
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.
Machine learning is used to approximate the kinetic energy of one dimensional diatomics as a functional of the electron density. The functional can accurately dissociate a diatomic, and can be systematically improved with training. Highly accurate self-consistent densities and molecular forces are found, indicating the possibility for ab-initio molecular dynamics simulations.
The semilocal meta generalized gradient approximation (MGGA) for the exchange-correlation functional of Kohn-Sham (KS) density functional theory can yield accurate ground-state energies simultaneously for atoms, molecules, surfaces, and solids, due to the inclusion of kinetic energy density as an input. We study for the first time the effect and importance of the dependence of MGGA on the kinetic energy density through the dimensionless inhomogeneity parameter, $alpha$, that characterizes the extent of orbital overlap. This leads to a simple and wholly new MGGA exchange functional, which interpolates between the single-orbital regime, where $alpha=0$, and the slowly varying density regime, where $alpha approx 1$, and then extrapolates to $alpha to infty$. When combined with a variant of the Perdew-Burke-Erzerhof (PBE) GGA correlation, the resulting MGGA performs equally well for atoms, molecules, surfaces, and solids.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
