From Microscopic Interactions to Density Functionals

Sandra Kemler

Published Jan 1, 2017• 2 citations

Abstract

Density functional theory provides the basis for uncounted studies of ground-state properties of many-body systems. However, the connection between the energy density functional and the underlying microscopic interactions of a given theory is not fully understood. In this thesis we use renormalization-group techniques in combination with density functional theory to study many-body systems from microscopic interactions. We apply our formalism to different one-dimensional systems. We start with systems of identical fermions interacting via a short-range repulsive and long-range attractive interaction which serves as a simple one-dimensional toy model for nuclei. After that we study systems of spin-1/2 fermions, where we assume interactions only between fermions with different spins and consider both a non-local interaction as in the previous case and a contact interaction. In particular, the contact interaction plays a prominent role for ultracold Fermi gases. We calculate ground-state properties such as the energy, density, intrinsic density and density correlation functions and compare our results to values obtained from other approaches. Moreover, we show how energies of excited states and the absolute square of the ground-state wavefunction can be extracted from the density-density correlation function. The relation between our formalism and conventional density functional theory as well as many-body perturbation theory is discussed which may help to guide the development of ab initio functionals for quantitative studies of nuclei in the future.

DOI: 10.26083/tuprints-00006013

Related papers

Formation of selfbound states in a one-dimensional nuclear model—a renormalization group based density functional study

Sandra Kemler, Martin Pospiech, Jens Braun

Abstract In nuclear physics, density functional theory (DFT) provides the basis for state-of-the art studies of ground-state properties of heavy nuclei. However, the direct relation of the density functional underlying these calculations and the microscopic nuclear forces is not yet fully understood. We present a combination of DFT and renormalization group (RG) techniques which allows to study selfbound many-body systems from microscopic interactions. We discuss its application with the aid of systems of identical fermions interacting via a long-range attractive and short-range repulsive two-body force in one dimension. We compute ground-state energies, intrinsic densities, and density correlation functions of these systems and compare our results to those obtained from other methods. In particular, we show how energies of excited states as well as the absolute square of the ground-state wave function can be extracted from the correlation functions within our approach. The relation between many-body perturbation theory and our DFT-RG approach is discussed and illustrated with the aid of the calculation of the second-order energy correction for a system of N identical fermions interacting via a general two-body interaction. Moreover, we discuss the control of spuriously emerging fermion self-interactions in DFT studies within our framework. In general, our approach may help to guide the development of energy functionals for future quantitative DFT studies of heavy nuclei from microscopic interactions.

arXiv:1606.04388Journal of Physics G Nuclear and Particle PhysicsDec 5, 201625 citations

Density-functional theory in one dimension for contact-interacting fermions

Rudolph J Magyar, Kieron Burke

A density-functional theory is developed for fermions in one dimension, interacting via a $\ensuremath{\delta}$ function. Such systems provide a natural testing ground for questions of principle, as the local-density approximation should be highly accurate since for this interaction type the exchange contribution to the local-density approximation is intrinsically self-interaction-free. The exact-exchange contribution to the total energy is a local functional of the density. A local-density approximation for correlation is obtained using perturbation theory and Bethe ansatz results for the one-dimensional contact-interacting uniform Fermi gas. The ground-state energies are calculated for two finite systems, the analogs of helium and of Hooke's atom. The local-density approximation is shown to be excellent as expected.

arXiv:cond-mat/0402108Physical Review ASep 21, 200461 citations

Towards a renormalization group approach to density functional theory—general formalism and case studies

Sandra Kemler, Jens Braun

We discuss a two-particle point-irreducible (2PPI) approach to many-body physics which relies on a renormalization group (RG) flow equation for the associated effective action. In particular, the general structure and properties of this RG flow equation are analyzed in detail. Moreover, we discuss how our 2PPI-RG approach relates to density functional theory and argue that it can in principle be used to study ground-state properties of non-relativistic many-body systems from microscopic interactions, such as (heavy) nuclei. For illustration purposes, we use our formalism to compute the ground-state properties of two toy models. Communicated by Achim Schwenk

arXiv:1304.1161Journal of Physics G Nuclear and Particle PhysicsJul 10, 201333 citations

Ground-state energy, density profiles, and momentum distribution of attractively interacting 1D Fermi gases with hard-wall boundaries: a Monte Carlo study

Shill,, Drut,, Porter Porter +1 more

Motivated by the realization of hard-wall boundary conditions in experiments with ultracold atoms, we investigate the ground-state properties of spin-1/2 fermions with attractive interactions in a one-dimensional box.We use lattice Monte Carlo methods to determine essential quantities like the energy, which we compute as a function of coupling strength and particle number in the regime from few to many particles.Many-fermion systems bound by hard walls display non-trivial density profiles characterized by so-called Friedel oscillations (which are similar to those observed in harmonic traps).In non-interacting systems, the characteristic length scale of the oscillations is set by (2k F ) -1 , where k F is the Fermi momentum, while repulsive interactions tend to generate Wignercrystal oscillations of period (4k F ) -1 .Based on the non-interacting result, we find a remarkably simple parametrization of the density profiles of the attractively interacting case, which we generalize to the one-body density matrix.While the total momentum is not a conserved quantity in the presence of hard walls, the magnitude of the momentum does provide a good quantum number.We are therefore able to provide a detailed characterization of the (quasi-)momentum distribution, which displays rather robust discontinuity at the Fermi surface.In addition, we determine the spatially varying on-site density-density correlation, which in turn yields Tan's contact density and, upon integration, Tan's contact.As is well known, the latter fully determines the short-range correlations and plays a crucial role in a multitude of equilibrium and non-equilibrium sum rules.

Apr 18, 20200 citations

Development and applications of density matrix functional theory of generalized Hubbard Models

Matthieu Saubanère

In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.

Jan 25, 20120 citations

Superfluid Local Density Approximation: A Density Functional Theory Approach to the Nuclear Pairing Problem

Aurel Bulgac

I describe the foundation of a Density Functional Theory approach to include pairing correlations, which was applied to a variety of systems ranging from dilute fermions, to neutron stars and finite nuclei. Ground state properties as well as properties of excited states and time-dependent phenomena can be achieved in this manner within a formalism based on microscopic input.

arXiv:1204.2207WORLD SCIENTIFIC eBooksMar 1, 20131 citation