Electronic properties of the dimerized one-dimensional Hubbard model using lattice density-functional theory

TL;DR

This paper applies lattice density-functional theory to study the electronic properties of the dimerized one-dimensional Hubbard model, introducing a new approximation for the interaction-energy functional and analyzing its effectiveness across different correlation regimes.

Contribution

The paper develops a simple explicit approximation for the interaction-energy functional in lattice DFT and applies it to analyze ground-state and excitation properties of dimerized chains.

Findings

The approximation accurately reproduces known solutions and numerical results.

LDFT effectively captures the crossover from weak to strong correlations.

Charge gaps and ground-state properties are reliably predicted across parameter ranges.

Abstract

The dimerized one-dimensional Hubbard model is studied in the framework of lattice density-functional theory (LDFT). The single-particle density matrix gamma_{ij} with respect to the lattice sites is considered as basic variable. The corresponding interaction-energy functional W[gamma_{ij}] is defined by Levy's constrained search. Exact numerical results are obtained for W(gamma_{12},gamma_{23}) where gamma_{12} = gamma_{i,i+1} for i odd and gamma_{23} = gamma_{i,i+1} for i even are the nearest-neighbor density-matrix elements along the chain. The domain of representability of gamma_{ij} and the functional dependence of W(gamma_{12},gamma_{23}) are analyzed. A simple, explicit approximation to W(gamma_{12},gamma_{23}) is proposed, which is derived from scaling properties of W, exact dimer results, and known limits. Using this approximation, LDFT is applied to determine ground-state properties and charge-excitation gaps of finite and infinite dimerized chains as a function of the Coulomb-repulsion strength U/t and of the alternation delta t of the hopping integrals t_{ij} (t_{ij} = t + or - delta t). The accuracy of the method is demonstrated by comparison with available exact solutions and accurate numerical calculations. Goals and limitations of the present approach are discussed particularly concerning its ability to describe the crossover from weak to strong electron correlations.