Quantum-mechanical calculations based on Kohn–Sham density functional theory (DFT) played a significant role in accurately predicting various aspects of materials behavior over the past decade. The Kohn–Sham approach to DFT reduces the many-body Schrodinger (eigen value) problem of interacting electrons into an equivalent problem of noninteracting electrons in an effective mean field that is governed by electron-density. Despite the reduced computational complexity of Kohn–Sham DFT, large-scale DFT calculations are still computationally very demanding with the resulting computational complexity scaling cubically with number of atoms in a given materials system. Numerical algorithms with reduced scaling behavior which are robust, computationally efficient and scalable on parallel computing architectures are always desirable to enable simulations at larger scales and on more complex systems. Following this line of thought, this study explores the use of tensor structured methods for ab-initio numerical solution of Kohn–Sham equations arising in DFT calculations. Earlier studies on tensor-structured methods have been quite successful in the accurate calculation of Hartree and the nonlocal exchange operators arising in the Hartree–Fock equations. A recent investigation of low-rank Tucker-type decomposition of the electron-density of large aluminum clusters (obtained from the finite-element discretization of orbital free DFT) shows the exponential decay of approximation error with respect to Tucker rank (number of tensor-basis functions in Tucker type representation). The results also indicate a smaller Tucker rank for the accurate representation of the electron density and is only weakly dependent on the system sizes studied. The promising success of tensor-structured techniques in resolving the electronic structure of material systems has enabled us to take a step further. In this study, we propose a systematic way of computing a globally adapted Tucker-type basis for solving the Kohn–Sham DFT problem by using a separable approximation of the Kohn–Sham Hamiltonian. Further, the resulting Kohn–Sham eigenvalue problem is projected into the aforementioned Tucker basis and is solved for ground-state energy using a self-consistent field iteration. The rank of the resulting Tucker representation and the computational complexity of the calculation are examined on representative benchmark examples involving metallic and insulating systems.
In this work we discuss the extraction of mean field single particle Hamiltonians from a many body wave function of a fermionic system. While we primarily discuss lattice models we also look at neon in an augccpvdz basis set as an example. Our approach allows us to discuss the result of a many particle wave function in terms of a non-interacting description. In contrast to density functional theory approaches on the lattice this approach allows the extraction of appropriate kinetic terms. In result the extracted mean field theories are closer to the physics of the problem under investigation. Therefore, the extracted mean field Hamiltonians may provide an improved starting point for perturbative approaches. In addition, the technique can be used to decide whether a density matrix renormalization group calculation for interacting fermions has converged to the true ground state.
By utilizing the knowledge that a Hamiltonian is a unique functional of its ground-state density, the following fundamental connections between densities and Hamiltonians are revealed: Given that ${\ensuremath{\rho}}_{\ensuremath{\alpha}}, {\ensuremath{\rho}}_{\ensuremath{\beta}},\dots{},{\ensuremath{\rho}}_{\ensuremath{\omega}}$ are ground-level densities for interacting or noninteracting Hamiltonians
H1,H2,…,HM (
M arbitrarily large) with local potentials
v1,
v2,
…,
vM, but given that we do not know which $\ensuremath{\rho}$ belongs with which
H, the correct mapping is possible and is obtained by minimizing $\ensuremath{\int}d\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}} [{v}_{1}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}){\ensuremath{\rho}}_{\ensuremath{\alpha}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})+{v}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}){\ensuremath{\rho}}_{\ensuremath{\beta}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})+\ensuremath{\cdots}{v}_{M}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}){\ensuremath{\rho}}_{\ensuremath{\omega}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})]$ with respect to optimum permutations of the $\ensuremath{\rho}$'s among the
v's. A tight rigorous bound connects a density to its interacting ground-state energy via the one-body potential of the interacting system and the Kohn-Sham effective one-body potential of the auxiliary noninteracting system. A modified Kohn-Sham effective potential is defined such that its sum of lowest orbital energies equals the true interacting ground-state energy. Moreover, of all those effective potentials which differ by additive constants and which yield the true interacting ground-state density, this modified effective potential is the most invariant with respect to changes in the one-body potential of the true Hamiltonian. With the exception of the occurrence of certain linear dependencies,
a density will not generally be associated with any ground-state wave function (is not wave function
v representable) if that density can be generated by a special linear combination of three or more densities that arise from a common set of degenerate ground-state wave functions. Applicability of the "constrained search" approach to density-functional theory is emphasized for non-
v-representable as well as for
v-representable densities. In fact, a particular constrained ensemble search is revealed which provides a general sufficient condition for non-
v representability by a wave function. The possible appearance of noninteger occupation numbers is discussed in connection with the existence of non-
v representability for some Kohn-Sham noninteracting systems.