We present a method to invert a given density and find the Kohn-Sham (KS) potential in Density Functional Theory (DFT) which shares that density. Our method employs the concept of screening density, which is naturally constrained by the inversion pro
cedure and thus ensures the density being inverted leads to a smooth KS potential with correct asymptotic behaviour. We demonstrate the applicability of our method by inverting both local (LDA) and non-local (Hartree-Fock and Coupled Cluster) densities; we also show how the method can be used to mitigate the effects of self-interactions in common DFT potentials with appropriate constraints on the screening density.
In the constrained minimization method of Gidopoulos and Lathiotakis (J. Chem. Phys. 136, 224109), the Hartree exchange and correlation Kohn-Sham potential of a finite $N$-electron system is replaced by the electrostatic potential of an effective cha
rge density that is everywhere positive and integrates to a charge of $N-1$ electrons. The optimal effective charge density (electron repulsion density, $rho_{rm rep}$) and the corresponding optimal effective potential (electron repulsion potential $v_{rm rep}$) are obtained by minimizing the electronic total energy in any density functional approximation. The two constraints are sufficient to remove the self-interaction errors from $v_{rm rep}$, correcting its asymptotic behavior at large distances from the system. In the present work, we describe, in complete detail, the constrained minimization method, including recent refinements. We also assess its performance in removing the self-interaction errors for three popular density functional approximations, namely LDA, PBE, and B3LYP, by comparing the obtained ionization energies to their experimental values for an extended set of molecules. We show that the results of the constrained minimizations are almost independent of the specific approximation with average percentage errors 15%, 14%, 13% for the above DFAs respectively. These errors are substantially smaller than the corresponding errors of the plain (unconstrained) Kohn-Sham calculations at 38%, 39%, and 27% respectively. Finally, we showed that this method correctly predicts negative values for the HOMO energies of several anions.
Some of the most spectacular failures of density-functional and Hartree-Fock theories are related to an incorrect description of the so-called static electron correlation. Motivated by recent progress on the N-representability problem of the one-body
density matrix for pure states, we propose a way to quantify the static contribution to the electronic correlation. By studying several molecular systems we show that our proposal correlates well with our intuition of static and dynamic electron correlation. Our results bring out the paramount importance of the occupancy of the highest occupied natural spin-orbital in such quantification.
In the present work, a method for the study of the structural deformations of two dimensional planar structures under uniaxial strain is presented. The method is based on molecular mechanics using the original stick and spiral model and a modified on
e which includes second nearest neighbor interactions for bond stretching. As we show, the method allows an accurate prediction of the structural deformations of any two dimensional planar structure as a function of strain, along any strain direction in the elastic regime, if structural deformations are known along specific strain directions, which are used to calculate the stick and spiral model parameters. Our method can be generalized including other strain conditions and not only uniaxial strain. We apply this method to graphene and we test its validity, using results obtained from {it ab initio} Density Functional Theory calculations. What we find is that the original stick and spiral model is not appropriate to describe accurately the structural deformations of graphene in the elastic regime. However, the introduction of second nearest neighbor interactions provides a very accurate description.
The concept of correlation is central to all approaches that attempt the description of many-body effects in electronic systems. Multipartite correlation is a quantum information theoretical property that is attributed to quantum states independent o
f the underlying physics. In quantum chemistry, however, the correlation energy (the energy not seized by the Hartree-Fock ansatz) plays a more prominent role. We show that these two different viewpoints on electron correlation are closely related. The key ingredient turns out to be the energy gap within the symmetry-adapted subspace. We then use a few-site Hubbard model and the stretched H$_2$ to illustrate this connection and to show how the corresponding measures of correlation compare.
We consider necessary conditions for the one-body-reduced density matrix (1RDM) to correspond to a triplet wave-function of a two electron system. The conditions concern the occupation numbers and are different for the high spin projections, $S_z=pm
1$, and the $S_z=0$ projection. Hence, they can be used to test if an approximate 1RDM functional yields the same energies for both projections. We employ these conditions in reduced density matrix functional theory calculations for the triplet excitations of two-electron systems. In addition, we propose that these conditions can be used in the calculation of triplet states of systems with more than two electrons by restricting the active space. We assess this procedure in calculations for a few atomic and molecular systems. We show that the quality of the optimal 1RDMs improves by applying the conditions in all the cases we studied.
Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the
optimal orbitals satisfy a single electron Schrodinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.
Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Colemans ensemble $N$-representability conditions. Recently, the topic of pure-state $N$-representability conditions, also known as generalized Pauli constraints
, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble $N$-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.
