We develop an functional framework to compute one-body quasiprobabilities for both fermionic and bosonic systems. Since the key variable is a quasidensity, this theory circumvents the problems of finding the Pauli potential or approximating the kinetic energy that are known to limit the applicability of standard orbital-free density functional theory. We present a set of strategies to (a) compute the one-body Wigner quasiprobability in an orbital-free manner from the knowledge of the universal functional and (b) obtain those functionals from the functionals of the one-body reduced density matrix (1-RDM). We find that the universal functional of optical lattices results from a translation, a contraction, and a rotation of the corresponding functional of the 1-RDM, revealing the strong connection between these two functional theories. Furthermore, we relate the key concepts of and . Published by the American Physical Society 2024
We develop an functional framework to compute one-body quasiprobabilities for both fermionic and bosonic systems. Since the key variable is a quasidensity, this theory circumvents the problems of finding the Pauli potential or approximating the kinetic energy that are known to limit the applicability of standard orbital-free density functional theory. We present a set of strategies to (a) compute the one-body Wigner quasiprobability in an orbital-free manner from the knowledge of the universal functional and (b) obtain those functionals from the functionals of the one-body reduced density matrix (1-RDM). We find that the universal functional of optical lattices results from a translation, a contraction, and a rotation of the corresponding functional of the 1-RDM, revealing the strong connection between these two functional theories. Furthermore, we relate the key concepts of and . Published by the American Physical Society 2024
Quantum Fisher information (QFI) is a central concept in quantum sciences used to quantify the ultimate precision limit of parameter estimation, detect quantum phase transitions, witness genuine multipartite entanglement, or probe nonlocality. Despite this widespread range of applications, computing the QFI value of quantum many-body systems is, in general, a very demanding task. Here we combine ideas from functional theories and quantum information to develop a functional framework for the QFI of fermionic and bosonic ground states. By relying upon the constrained-search approach, we demonstrate that the QFI matrix terms can universally be determined by the one-body reduced density matrix (1-RDM), thus avoiding the use of exponentially large wave functions. Furthermore, we show that QFI functionals can be determined from the universal 1-RDM functional by calculating its derivatives with respect to the coupling strengths, thus becoming the generating functional of the QFI. We showcase our approach with the Bose-Hubbard model and present exact analytical and numerical QFI functionals. Our results provide the first connection between the one-body reduced density matrix functional theory and the quantum Fisher information. Published by the American Physical Society 2024
We showcase the advantages of orbital-free density-potential functional theory (DPFT), a more flexible variant of Hohenberg-Kohn density functional theory. DPFT resolves the usual trouble with the gradient-expanded kinetic energy functional by facilitating systematic semiclassical approximations in terms of an effective potential energy that incorporates all interactions. With the aid of two systematic approximation schemes we demonstrate that DPFT is not only scalable, universally applicable in both position and momentum space, and allows kinetic and interaction energy to be approximated consistently, but can also compete with highly accurate, yet restricted, methods. As two- and three-dimensional geometries are extensively covered elsewhere, our focus here is on one-dimensional settings, with semiclassical observables systematically derived from both the Wigner function formalism and a split-operator approach. The high quality of our results for Fermi gases in Morse potentials invites the use of DPFT for describing more exotic systems, such as trapped large-spin fermion mixtures with contact or dipole-dipole interactions.
Time-dependent orbital-free density functional theory (DFT) is an efficient method for calculating the dynamic properties of large-scale quantum systems due to the low computational cost compared to standard time-dependent DFT. We formalize this method by mapping the real system of interacting fermions onto a fictitious system of noninteracting bosons. The dynamic Pauli potential and associated kernel emerge as key ingredients of time-dependent orbital-free DFT. Using the uniform electron gas as a model system, we derive an approximate frequency-dependent Pauli kernel. Pilot calculations suggest that space nonlocality is a key feature for this kernel. Nonlocal terms arise already in the second-order expansion with respect to unitless frequency and reciprocal space variable ($\frac{\ensuremath{\omega}}{q\phantom{\rule{0.16em}{0ex}}{k}_{F}}$ and 2kFq, respectively). Given the encouraging performance of the proposed kernel, we expect it will lead to more accurate orbital-free DFT simulations of nanoscale systems out of equilibrium. Additionally, the proposed path to formulate nonadiabatic Pauli kernels presents several avenues for further improvements which can be exploited in future work to improve the results.
From drug design to material science, computational chemistry has become a valuable tool for many fields of applications in the past decades. Quantum chemistry in particular, can provide accurate descriptions of molecules and solids by relying on principles of quantum mechanics with the Schrödinger equation at its core. Unfortunately this equation becomes prohibitively expensive to solve for systems made of more than a few electrons and one has to rely on other approaches. Among the most popular is Density Functional Theory, relatively inexpensive and exact in principle. One has, however, to introduce approximations in practice, deteriorating its accuracy for challenging, so-called strongly correlated, systems. On the contrary, other methods like CASSCF remain reliable for these systems, but although cheaper than the solution of the Schrödinger equation, remain too expensive for extended systems. There is therefore room for methods in between, that could handle large correlated systems with reasonable accuracy. One such method is the One Body Reduced Density Matrix (1-RDM) Functional Theory (RDMFT), based on the fact that an electronic system can, in principle exactly, be described by its 1-RDM only. Although its theory dates back to the 1970s, computational hindrances have prevented this method from shining so far. Indeed, computations in RDMFT rely on a constrained optimization of the eigenvalues and eigenvectors of the 1-RDM, that can be reformulated as an unconstrained optimization by parameterizing the latter but in both cases suffers from convergence problems. In this thesis, the role of this parametrization is analyzed and motivates a now second order optimization algorithm, which improves the convergence speed by two orders of magnitude w.r.t. the state of the art algorithm. This understanding of the problem also inspires an alternative approach with the 1-RDM as the optimization variable. This however requires the development of new approximations and a first investigation in that direction, based on Machine Learning is performed. Lastly, a generalization of the method to excited states is attempted, highlighting some key computational challenges.
The one-body reduced density matrix γ plays a fundamental role in describing and predicting quantum features of bosonic systems, such as Bose-Einstein condensation. The recently proposed reduced density matrix functional theory for bosonic ground states establishes the existence of a universal functional F[γ] that recovers quantum correlations exactly. Based on a novel decomposition of γ, we have developed a method to design reliable approximations for such universal functionals: our results suggest that for translational invariant systems the constrained search approach of functional theories can be transformed into an unconstrained problem through a parametrization of an Euclidian space. This simplification of the search approach allows us to use standard machine-learning methods to perform a quite efficient computation of both F[γ] and its functional derivative. For the Bose-Hubbard model, we present a comparison between our approach and Quantum Monte Carlo.
Density functional theory constitutes the workhorse of modern electronic structure calculations due to its favourable computational cost despite the fact that it usually fails to describe strongly correlated systems. A particularly promising approach to overcome those difficulties is reduced density matrix functional theory (RDMFT): It abandons the complexity of the N-particle wave function and at the same time explicitly allows for fractional occupation numbers. It is the goal of this thesis to initiate and establish a bosonic RDMFT for both ground state and excited state energy calculations. Motivated by the Onsager and Penrose criterion which identifies RDMFT as a particularly suitable approach to describe Bose-Einstein condensates (BECs), we derive the universal functional for a homogeneous BEC in the Bogoliubov regime. Remarkably, the gradient of the universal functional is found to diverge repulsively in the regime of complete condensation. This introduces the new concept of a BEC force, which provides a universal explanation for quantum depletion since it is merely based on the geometry of quantum states. In the second part of the thesis, we propose and work out an ensemble RDMFT targeting excitations in bosonic quantum systems. This endeavour further highlights the potential of convex analysis for the development of functional theories in the future. Indeed by resorting to several concepts from convex analysis, we succeeded to provide a comprehensive foundation of w-ensemble RDMFT for bosons which is further based on a generalization of the Ritz variational principle and a constrained search formalism. In particular, we solve the emerging N-representability problem leading to a generalization of Pauli's famous exclusion principle to bosonic mixed states.