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, 2012•0 citations