The dimerized one-dimensional Hubbard model is studied in the framework of lattice density-functional theory (LDFT). The single-particle density matrix ${\gamma }_{\mathrm{ij}}$ with respect to the lattice sites is considered as the basic variable. The corresponding interaction-energy functional $W\left[{\gamma }_{\mathrm{ij}}\right]$ 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\mathrm{}}$ for odd i and ${\gamma }_{23}={\gamma }_{i,i+1\mathrm{}}$ for even i are the nearest-neighbor density-matrix elements along the chain. The domain of representability of ${\gamma }_{\mathrm{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_{\mathrm{ij}}$ ${(t}_{\mathrm{ij}}=t\pm \delta t).\mathrm{}$ 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.

• Received 1 August 2002

DOI:https://doi.org/10.1103/PhysRevB.67.035115