Orbitals from local RDMFT: Are they Kohn-Sham or natural orbitals?

Variational minimization scheme for the one-particle reduced density matrix functional theory in the ensemble N-representability domain

The one-particle reduced density-matrix (1-RDM) functional theory is a promising alternative to density-functional theory (DFT) that uses the 1-RDM rather than the electronic density as a basic variable. However, long-standing challenges such as the lack of the Kohn-Sham scheme and the complexity of the pure N-representability conditions are still impeding its wild utilization. Fortunately, ensemble N-representability conditions derived in the natural orbital basis are known and trivial such that almost every functional of the 1-RDM is actually natural orbital functional, which does not perform well for all the correlation regimes. In this work, we propose a variational minimization scheme in the ensemble N-representable domain that is not restricted to the natural orbital representation of the 1-RDM. We show that splitting the minimization into the diagonal and off-diagonal parts of the 1-RDM can open the way toward the development of functionals of the orbital occupations, which remains a challenge for the generalization of site-occupation functional theory in chemistry. Our approach is tested on the uniform Hubbard model using the Müller and the Töws-Pastor functionals, as well as on the dihydrogen molecule using the Müller functional.

Effective local potentials for density and density-matrix functional approximations with non-negative screening density

A way to improve the accuracy of the spectral properties in density functional theory (DFT) is to impose constraints on the effective, Kohn-Sham (KS), local potential [J. Chem. Phys. 136, 224109 (2012)]. As illustrated, a convenient variational quantity in that approach is the "screening" or "electron repulsion" density, ρrep, corresponding to the local, KS Hartree, exchange and correlation potential through Poisson's equation. Two constraints, applied to this minimization, largely remove self-interaction errors from the effective potential: (i) ρrep integrates to N - 1, where N is the number of electrons, and (ii) ρrep ≥ 0 everywhere. In this work, we introduce an effective "screening" amplitude, f, as the variational quantity, with the screening density being ρrep = f2. In this way, the positivity condition for ρrep is automatically satisfied, and the minimization problem becomes more efficient and robust. We apply this technique to molecular calculations, employing several approximations in DFT and in reduced density matrix functional theory. We find that the proposed development is an accurate, yet robust, variant of the constrained effective potential method.

Density inversion method for local basis sets without potential auxiliary functions: inverting densities from RDMFT

A density inversion method is presented, to obtain the constrained, optimal, local potential that has a prescribed asymptotic behaviour and reproduces optimally any given ground-state electronic density. This work builds...

Electronic Properties of Carbon Nanobelts Predicted by Thermally-Assisted-Occupation DFT

Accurate prediction of properties of large-scale multi-reference (MR) electronic systems remains difficult for traditional computational methods (e.g., the Hartree-Fock theory and Kohn-Sham density functional theory (DFT)). Recently, thermally-assisted-occupation (TAO)-DFT has been demonstrated to offer reliable description of electronic properties of various large-scale MR electronic systems. Consequently, in this work, TAO-DFT is used to unlock the electronic properties associated with C-Belt[n] (i.e., the carbon nanobelts containing n fused 12-membered carbon rings). Our calculations show that for all the system sizes reported (n = 4-24), C-Belt[n] have singlet ground states. In general, the larger the size of C-Belt[n], the more pronounced the MR character of ground-state C-Belt[n], as evident from the symmetrized von Neumann entropy and the occupation numbers of active TAO-orbitals. Furthermore, the active TAO-orbitals are delocalized along the circumference of C-Belt[n], as evident from the visualization of active TAO-orbitals.

Quantum Information and Algorithms for Correlated Quantum Matter

Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.

Kinetic-Energy Density-Functional Theory on a Lattice

We present a kinetic-energy density-functional theory and the corresponding kinetic-energy Kohn–Sham (keKS) scheme on a lattice and show that, by including more observables explicitly in a density-functional approach, already simple approximation strategies lead to very accurate results. Here, we promote the kinetic-energy density to a fundamental variable alongside the density and show for specific cases (analytically and numerically) that there is a one-to-one correspondence between the external pair of on-site potential and site-dependent hopping and the internal pair of density and kinetic-energy density. On the basis of this mapping, we establish two unknown effective fields, the mean-field exchange-correlation potential and the mean-field exchange-correlation hopping, which force the keKS system to generate the same kinetic-energy density and density as the fully interacting one. We show, by a decomposition based on the equations of motions for the density and the kinetic-energy density, that we can construct simple orbital-dependent functionals that outperform the corresponding exact-exchange Kohn–Sham (KS) approximation of standard density-functional theory. We do so by considering the exact KS and keKS systems and comparing the unknown correlation contributions as well as by comparing self-consistent calculations based on the mean-field exchange (for the effective potential) and a uniform (for the effective hopping) approximation for the keKS and the exact-exchange approximation for the KS system, respectively.

Natural occupation numbers in two-electron quantum rings

Natural orbitals (NOs) are central constituents for evaluating correlation energies through efficient approximations. Here, we report the closed-form expression of the NOs of two-electron quantum rings, which are prototypical finite-extension systems and new starting points for the development of exchange-correlation functionals in density functional theory. We also show that the natural occupation numbers for these two-electron paradigms are in general non-vanishing and follow the same power law decay as atomic and molecular two-electron systems.