We present a theoretical study on the response properties to an external electric field of strongly correlated one-dimensional metals. Our investigation is based on the recently developed Bethe ansatz local-density approximation (BALDA) to the density functional theory formulation of the Hubbard model. This is capable of describing both Luttinger liquid and Mott-insulator correlations. The BALDA calculated values for the static linear polarizability are compared with those obtained by numerically accurate methods, such as exact (Lanczos) diagonalization and the density-matrix renormalization group, over a broad range of parameters. In general BALDA linear polarizabilities are in good agreement with the exact results. The response of the exact exchange and correlation potential is found to point in the same direction of the perturbing potential. This is well reproduced by the BALDA approach, although the fine details depend on the specific parameterization for the local approximation. Finally we provide a numerical proof for the nonlocality of the exact exchange and correlation functional.
We present a theoretical study on the response properties to an external electric field of strongly correlated one-dimensional metals. Our investigation is based on the recently developed Bethe ansatz local-density approximation (BALDA) to the density functional theory formulation of the Hubbard model. This is capable of describing both Luttinger liquid and Mott-insulator correlations. The BALDA calculated values for the static linear polarizability are compared with those obtained by numerically accurate methods, such as exact (Lanczos) diagonalization and the density-matrix renormalization group, over a broad range of parameters. In general BALDA linear polarizabilities are in good agreement with the exact results. The response of the exact exchange and correlation potential is found to point in the same direction of the perturbing potential. This is well reproduced by the BALDA approach, although the fine details depend on the specific parameterization for the local approximation. Finally we provide a numerical proof for the nonlocality of the exact exchange and correlation functional.
Since the advent of high-Tc cuprate superconductors in 1986, strongly correlated electron systems have attracted much attention. Since the cuprates are essentially two-dimensional, low-dimensional systems have moved into the focus of condensed-matter theory. From a theoretical point of view, one-dimensional systems are of particular interest because there are exact numerical and analytical methods which permit detailed studies and deep insights into the many-body problem. In the first part of this thesis, using the numerically exact methods Exact Diagonalization and the Density-Matrix Renormalization Group (DMRG), we investigate the properties of the Tomonaga–Luttinger liquid which is the generic metallic state of matter in one dimension. In particular, we concentrate on the investigation of the so-called Tomonaga–Luttinger liquid parameter which determines the critical exponent α for the density of states near the Fermi energy. Experimental results for some quasi one-dimensional materials report α≳1, which would imply Kρ≲0.17, a value which cannot be reconciled with the bare Hubbard model where KρH≥0.5, i.e., αH≲1/8. We develop new accurate numerical methods to obtain Kρ and investigate how to obtain such small values for Kρ for slightly doped charge-density-wave insulators. In the second part of this thesis, using the Thermodynamic Bethe Ansatz (TBA) as exact analytical method, we investigate the one-dimensional Hubbard model in the spindisordered regime, which is characterized by the temperature being much larger than the magnetic energy scale but small compared to the Mott–Hubbard gap. Our study is motivated by the controversy about the Mott–Hubbard insulator in infinite dimensions whose ground state is also spin-disordered. In this system the determination of the precise value of the critical interaction strength Uc where the Mott–Hubbard gap closes is still unsolved. Therefore, we provide an example of a Hubbard-type model with a disordered spin background which can be solved exactly. The thermodynamics of our model can be understood in terms of gapped charged excitations with an effective dispersion which are decoupled form the spin degrees of freedom; the latter contribute only entropically. An interpretation of this regime in terms of a putative interacting-electron system at zero temperature leads to a metal-insulator transition at a finite interaction strength above which the gap opens linearly. Our exact results indicate that the strong-coupling expansion of the ground-state energy cannot be used to locate Uc. However, the strong-coupling expansion of the gap permits a reliable extrapolation of the critical interaction strength.
The response of a one-dimensional fermion system is investigated using density functional theory (DFT) within the local-density approximation (LDA) and compared with exact results. It is shown that DFT-LDA reproduces surprisingly well some of the characteristic features of the Luttinger liquid, for example, the dispersion of the collective charge excitations. On the other hand, the approximation fails even qualitatively for details of the dynamic response and for quantities for which backscattering is important, i.e., those quantities which are crucial for an accurate description of transport. In particular, the Drude weight in the presence of a single impurity is discussed.
One-dimensional conductors are a long-standing topic of research with direct applications to organic conductors and mesoscopic rings. The discovery of the ceramic high-temperature superconductors has revitalized the interest in low-dimensional charge and spin fluctuations of highly correlated electron systems. Several mechanisms proposed to explain the high-T c superconductors invoke properties of the two-dimensional Hubbard model, but probably also some one-dimensional aspects are relevant. Numerous one-dimensional models for correlated electrons have been studied with various approximate, asymptotically exact and exact methods. These results lead to the concept of Luttinger liquid for interacting electron gases without excitation gaps (metallic systems). Characteristic of Luttinger liquids are the charge and spin separation, marginal Fermi liquid properties, e.g. the absence of quasiparticles in the vicinity of the Fermi surface, nonuniversal power-law singularities in the one-particle spectral function and the related absence of a discontinuity in the momentum distribution at the Fermi level, the power-law decay of correlation functions for long times and large distances, persistent currents in finite rings, etc. Due to the peculiarities of the phase space in one dimension some of the models have sufficient conserved currents to be completely integrable. We review exact results derived within the framework of Bethe's ansatz for integrable one-dimensional models of correlated electrons. The Bethe-ansatz method is presented by explicitly showing the steps leading to the solution of the N-component electron gas interacting via a δ-function potential (repulsive and attractive interaction), which is probably the simplest model of correlated electrons. Emphasis is given to the procedure to extract the groundstate properties, the classification of states, the excitation spectrum, the thermodynamics and finite size effects, such as critical exponents of correlation functions and persistent currents. The method is then applied to numerous other models, e.g. (i) a two-band model involving attractive and repulsive potentials and crystalline fields splitting the bands, (ii) the traditional Hubbard chain with attractive and repulsive U, (iii) the degenerate Hubbard model with repulsive U, which displays a metal–insulator transition at a finite U, (iv) a two-band Hubbard model with repulsive U, (v) the traditional supersymmetric t–J model (vi) a two-band supersymmetric t–J model with band-splitting and (vii) the N-component supersymmetric t–J model. Finally, results for models with long-range interactions, in particular r-2 and sinh -2(r) potentials, are briefly reviewed.
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.
A density-functional theory (DFT) of lattice fermion models is presented, which uses the single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ as basic variable. A simple, explicit approximation to the interaction-energy functional $W[\ensuremath{\gamma}]$ of the Hubbard model is derived from exact dimer results, scaling properties of $W[\ensuremath{\gamma}]$ and known limits. Systematic tests on the one-dimensional chain show a remarkable agreement with the Bethe-ansatz exact solution for all interaction regimes and band fillings. New results are obtained for the ground-state energy and charge-excitation gap in two dimensions. A successful description of strong electron correlations within DFT is achieved.
In the first part of our theoretical study of correlated atomic wires on substrates, we introduced lattice models for a one-dimensional quantum wire on a three-dimensional substrate and their approximation by quasi-one-dimensional effective ladder models [Abdelwahab et al., preceding paper, Phys. Rev. B 96, 035445 (2017)]. In this second part, we apply this approach to the case of a correlated wire with a Hubbard-type electron-electron repulsion deposited on an insulating substrate. The ground-state and spectral properties are investigated numerically using the density-matrix renormalization group method and quantum Monte Carlo simulations. As a function of the model parameters, we observe various phases with quasi-one-dimensional low-energy excitations localized in the wire, namely, paramagnetic Mott insulators, Luttinger liquids, and spin-1/2 Heisenberg chains. The validity of the effective ladder models is assessed for selected parameters by studying the dependence of results on the number of legs and comparing to the full three-dimensional model. We find that narrow ladder models accurately reproduce the quasi-one-dimensional excitations of the full three-dimensional model but predict only qualitatively whether excitations are localized around the wire or delocalized in the three-dimensional substrate.