Single hole dynamics in the one-dimensional - model

We present a new finite-temperature quantum Monte Carlo algorithm to compute imaginary-time Green functions for a single hole in the t-J model on non-frustrated lattices. Spectral functions are obtained with the Maximum Entropy method. Simulations of the one-dimensional case show that a simple charge-spin separation Ansatz is able to describe the overall features of the spectral function such as the bandwidth and the compact support of the spectral function, over the whole energy range for values of J / t from 1/3 to 4. This is contrasted with the two-dimensional case. The quasiparticle weight Z_k is computed on lattices up to L =128 sites in one dimension, and scales as . Received 15 February 2000     

Spin susceptibility of interacting electrons in one dimension: Luttinger liquid and lattice effects

The temperature-dependent uniform magnetic susceptibility of interacting electrons in one dimension is calculated using several methods. At low temperature, the renormalization group reveals that the Luttinger liquid spin susceptibility approaches zero temperature with an infinite slope in striking contrast with the Fermi liquid result and with the behavior of the compressibility in the absence of umklapp scattering. This effect comes from the leading marginally irrelevant operator, in analogy with the Heisenberg spin 1/2 antiferromagnetic chain. Comparisons with Monte Carlo simulations at higher temperature reveal that non-logarithmic terms are important in that regime. These contributions are evaluated from an effective interaction that includes the same set of diagrams as those that give the leading logarithmic terms in the renormalization group approach. Comments on the third law of thermodynamics as well as reasons for the failure of approaches that work in higher dimensions are given. Received 2 March 1999     

Charge density waves and bond order waves in a quarter filled extended Hubbard ladder

We investigate the phase diagram of a quarter filled Hubbard ladder with nearest-neighbor Coulomb repulsion using bosonization and renormalization group approach. Focusing on the strong-repulsion regime, we discuss the effect of an interchain exchange interaction J and interchain repulsion V on the possible ground states of the system and charge order configurations. Since the spin excitations always possess a gap, we find competing bond-order wave and charge density wave phases as possible ground states of the ladder model. We discuss the elementary excitations in these various phases and point an analogy between the excitations on some of these phases and those of a Kondo-Heisenberg insulator. We also study the order of the quantum phase transitions between the different ground states of the system. We obtain second order transitions in the Ising or SU(2)₂ universality class or first order transitions. We map the complete phase diagram in the J –V plane by integrating perturbative renormalization-group equations. Finally, we discuss the effect of doping away from half-filling and the effect of an applied magnetic field.     

Luttinger liquid with asymmetric dispersion

We present an extension of the Tomonaga-Luttinger model in which left and right-moving particles have different Fermi velocities. We derive expressions for one-particle Green's functions, momentum-distributions, density of states, charge compressibility and conductivity as functions of both the velocity difference ε and the strength of the interaction β. This allows us to identify a novel restricted region in the parameter space in which the system keeps the main features of a Luttinger liquid but with an unusual behavior of the density of states and the static charge compressibility κ. In particular κ diverges on the boundary of the restricted region, indicating the occurrence of a phase transition. Received 20 May 2002 / Received in final form 23 August 2002 Published online 19 November 2002     

Competition of disorder and interchain hopping in a two-chain Hubbard model

We study the interplay of Anderson localization and interaction in a two chain Hubbard ladder allowing for arbitrary ratio of disorder strength to interchain coupling. We obtain three different types of spin gapped localized phases depending on the strength of disorder: a pinned 4k _F Charge Density Wave (CDW) for weak disorder, a pinned 2k _F CDW^π for intermediate disorder and two independently pinned single chain 2k _F CDW for strong disorder. Confinement of electrons can be obtained as a result of strong disorder or strong attraction. We give the full phase diagram as a function of disorder, interaction strength and interchain hopping. We also study the influence of interchain hopping on localization length and show that localization is enhanced by a small interchain hopping but suppressed by a large interchain hopping. Received 6 April 2001     

Charge localization in organic conductors (TM)X: The influence of anion ordering

An investigation of the different contributions leading to charge localization in a 1/2 or 1/4 filled band 1D conductor has been conducted through a study of transport properties in the solid solution [(TMTSF)_1-x (TMTTF) _x]₂ReO₄. The existence of an ordering transition of the anions allows to identify two contributions to the electronic potential with wave vector 4k_F. A dominant on-site 4k_F potential besides the bond contribution is revealed when Umklapp scattering is pertinent via the weakening of the localization arising at the (0, 1/2, 1/2) anion ordering which is stabilized under pressure in the compound [(TMTSF) _0.5 (TMTTF)_0.5]₂ReO₄ at variance with the enhancement of localization observed in the homomolecular (TMTTF)₂ReO ₄ material. Received: 13 May 1998 / Revised: 8 July 1998 / Accepted: 9 July 1998     

Dynamical mean-field theory of a simplified double-exchange model

Simplified double-exchange model including transfer of the itinerant electrons with spin parallel to the localized spin in the same site and the indirect interaction J of kinetic type between localized spins is comprihensively investigated. The model is exactly solved in infinite dimensions. The exact equations describing the main ordered phases (ferromagnetic and antiferromagnetic) are obtained for the Bethe lattice with (z is the coordination number) in analytical form. The exact expression for the generalized paramagnetic susceptibility of the localized-spin subsystem is also obtained in analytical form. It is shown that temperature dependence of the uniform and the staggered susceptibilities has deviation from Curie-Weiss law. Dependence of Curie and Néel temperatures on itinerant-electron concentration is discussed to study instability conditions of the paramagnetic phase. Anomalous temperature behaviour of the chemical potential, the thermopower and the specific heat is investigated near the Curie point. It is found for J=0 that the system is unstable towards temperature phase separation between ferromagnetic and paramagnetic states. A phase separation connected with antiferromagnetic and the paramagnetic phases can occur only at . Zero-temperature phase diagram including the phase separation between ferromagnetic and antiferromagnetic states is given. Received 28 May 1999 and Received in final form 14 July 1999     

Phase separation in the strongly correlated Falicov-Kimball model in infinite dimensions

Phase separation in the strongly correlated Falicov-Kimball model in infinite dimensions is examined. We show that the phase separation can occur for any values of the interaction constant J^* when the site energy of the localized electrons is equal to zero. Electron-poor regions always have homogeneous state and electron-rich regions have chessboard state for , chessboard state or homogeneous state in dependence upon temperature for 0<J ^* <0.03 and homogeneous state for J ^* =0. For J ^* =0 and T=0, phase separation (segregation) occurs at .The obtained results are exact for the Bethe lattice with infinite number of the nearest neighbours. Received 1 December 1998 and Received in final form 12 April 1999     

Thermal evolution of hole dynamics in the large-dimensional model

The dynamics of a single hole in the t - J model is solved exactly for all temperature, T, in the limit of large spatial dimensions, , using the Feenberg renormalized perturbation series. We focus in particular on single-particle spectra, together with optical and static hole conductivities. Explicit results are illustrated for a Bethe lattice, and exemplify the continuous thermal evolution of the underlying string picture from the T =0 string-pinned limit through to the paramagnetic phase. Quenched site-disorder is also readily incorporated, exact results thereby being obtained for the interplay between disorder and thermally-induced hole dynamics. Received 16 September 1998     

Density-matrix renormalization study of the Hubbard model [4]on a Bethe lattice

The half-filled Hubbard model on the Bethe lattice with coordination number z=3 is studied using the density-matrix renormalization group (DMRG) method. Ground-state properties such as the energy per site E, average local magnetization , its fluctuations and various spin correlation functions are determined as a function of the Coulomb interaction strength U/t. The local magnetic moments increase monotonically with increasing Coulomb repulsion U/t showing antiferromagnetic order between nearest neighbors []. At large U/t, is strongly reduced with respect to the saturation value 1/2 due to exchange fluctuations between nearest neighbors (NN) spins [ for ]. shows a maximum for U/t=2.4-2.9 that results from the interplay between the usual increase of with increasing U/t and the formation of important permanent moments at large U/t. While NN sites show antiferromagnetic spin correlations that increase with increasing Coulomb repulsion, the next NN sites are very weakly correlated over the whole range of U/t. The DMRG results are discussed and compared with tight-binding calculations for U=0, independent DMRG studies for the Heisenberg model and simple first-order perturbation estimates. Received 8 February 1999 and Received in final form 14 June 1999     

to -12pt [-12pt]to 16ptRapid Note Wavefunctions for the Luttinger liquid

Standard bosonization techniques lead to phonon-like excitations in a Luttinger liquid (LL), reflecting the absence of Landau quasiparticles in these systems. Yet in addition to the above excitations some LL are known to possess solitonic states carrying fractional quantum numbers (e.g. the spin 1/2 Heisenberg chain). We have reconsidered the zero modes in the low-energy spectrum of the Gaussian boson LL Hamiltonian both for fermionic and bosonic LL: in the spinless case we find that two elementary excitations carrying fractional quantum numbers allow to generate all the charge and current excited states of the LL. We explicitly compute the wavefunctions of these two objects and show that one of them can be identified with the 1D version of the Laughlin quasiparticle introduced in the context of the Fractional Quantum Hall effect. For bosons, the other quasiparticle corresponds to a spinon excitation. The eigenfunctions of Wen's chiral LL Hamiltonian are also derived: they are quite simply the one dimensional restrictions of the 2D bulk Laughlin wavefunctions. Received 26 January 1999 and Received in final form 21 April 1999     

Spin 3/2 fermions with attractive interactions in a one-dimensional optical lattice: phase diagrams,entanglement entropy,and the effect of the trap

We study spin 3/2 fermionic cold atoms with attractive interactions confined in a one-dimensional optical lattice. Using numerical techniques, we determine the phase diagram for a generic density. For the chosen parameters, one-particle excitations are gapped and the phase diagram is separated into two regions: one where the two-particle excitation gap is zero, and one where it is finite. In the first region, the two-body pairing fluctuations (BCS) compete with the density ones. In the other one, a molecular superfluid (MS) phase, in which bound-states of four particles form, competes with the density fluctuations. The properties of the transition line between these two regions is studied through the behavior of the entanglement entropy. The physical features of the various phases, comprising leading correlations, Friedel oscillations, and excitation spectra, are presented. To make the connection with experiments, the effect of a harmonic trap is taken into account. In particular, we emphasize the conditions under which the appealing MS phase can be realized, and how the phases could be probed by using the density profiles and the associated structure factor. Lastly, the consequences on the flux quantization of the different nature of the pairing in the BCS and MS phases are studied in a situation where the condensate is in a ring geometry.     

Spectral broadening due to long-range Coulomb interactions in the molecular metal TTF-TCNQ

We employ density-functional theory to calculate realistic parameters for an extended Hubbard model of the molecular metal TTF-TCNQ. Considering both intra- and intermolecular screening in the crystal, we confirm the importance of the suspected longer-range Coulomb interactions along the molecular stacks, as well as inter-stack coupling. Contrary to past belief, these terms do not lead to the formation of a Wigner lattice, but simply broaden the spectral function.     

Luttinger liquids with boundaries: Power-laws and energy scales

We present a study of the one-particle spectral properties for a variety of models of Luttinger liquids with open boundaries. We first consider the Tomonaga-Luttinger model using bosonization. For weak interactions the boundary exponent of the power-law suppression of the spectral weight close to the chemical potential is dominated by a term linear in the interaction. This motivates us to study the spectral properties also within the Hartree-Fock approximation. It already gives power-law behavior and qualitative agreement with the exact spectral function. For the lattice model of spinless fermions and the Hubbard model we present numerically exact results obtained using the density-matrix renormalization-group algorithm. We show that many aspects of the behavior of the spectral function close to the boundary can again be understood within the Hartree-Fock approximation. For the repulsive Hubbard model with interaction U the spectral weight is enhanced in a large energy range around the chemical potential. At smaller energies a power-law suppression, as predicted by bosonization, sets in. We present an analytical discussion of the crossover and show that for small U it occurs at energies exponentially (in -1/U) close to the chemical potential, i.e. that bosonization only holds on exponentially small energy scales. We show that such a crossover can also be found in other models. Received 8 February 2000 and Received in final form 25 April 2000     

Strong-coupling perturbation theory of the Hubbard model

The strong-coupling perturbation theory of the Hubbard model is presented and carried out to order (t/U)⁵ for the one-particle Green function in arbitrary dimension. The spectral weight is expressed as a Jacobi continued fraction and compared with new Monte-Carlo data of the one-dimensional, half-filled Hubbard model. Different regimes (insulator, conductor and short-range antiferromagnet) are identified in the temperature-hopping integral (T,t) plane. This work completes a first paper on the subject (Phys. Rev. Lett. 80, 5389 (1998)) by providing details on diagrammatic rules and higher-order results. In addition, the non half-filled case, infinite resummations of diagrams and the double occupancy are discussed. Various tests of the method are also presented. Received 25 October 1999     

Flow to strong coupling in the two-dimensional Hubbard model

We extend the analysis of the renormalization group flow in the two-dimensional Hubbard model close to half-filling using the recently developed temperature flow formalism. We investigate the interplay of d-density wave and Fermi surface deformation tendencies with those towards d-wave pairing and antiferromagnetism. For a ratio of next nearest to nearest neighbor hoppings, t'/t = - 0.25, and band fillings where the Fermi surface is inside the Umklapp surface, only the d-pairing susceptibility diverges at low temperatures. When the Fermi surface intersects the Umklapp surface close to the saddle points, d-wave pairing, d-density wave, antiferromagnetic and, to a weaker extent, d-wave Fermi surface deformation susceptibilities grow together when the interactions flow to strong coupling. We interpret these findings as indications for a non-trivial strongly coupled phase with short-ranged superconducting and antiferromagnetic correlations, in close analogy with the spin liquid ground state in the well-understood two-leg Hubbard ladder. Received 23 January 2002     

Tails of the dynamical structure factor of 1D spinless fermions beyond the Tomonaga approximation

We consider one-dimensional (1D) interacting spinless fermions with a non-linear spectrum in a clean quantum wire (non-linear bosonization). We compute diagrammatically the 1D dynamical structure factor, S(ω,q), beyond the Tomonaga approximation focusing on it's tails, |ω| ≫vq, i.e. the 2-pair excitation continuum due to forward scattering. Our methodology reveals three classes of diagrams: two “chiral” classes which bring divergent contributions in the limits ω→±vq, i.e. near the single-pair excitation continuum, and a “mixed” class (so-called Aslamasov-Larkin or Altshuler-Shklovskii type diagrams) which is crucial for the f-sum rule to be satisfied. We relate our approach to the T=0 ones present in the literature. We also consider the case and show that the 2-pair excitation continuum dominates the single-pair one in the range: |q|T/k_F ≪ω±vq ≪T (substantial for q ≪k_F). As applications we first derive the small-momentum optical conductivity due to forward scattering: σ∼1/ω for T ≪ω and σ∼T/ω² for T ≫ω. Next, within the 2-pair excitation continuum, we show that the attenuation rate of a coherent mode of dispersion Ω_q crosses over from , e.g. γ_q ∼|q|³ for an acoustic mode, to , independent of Ω_q, as temperature increases. Finally, we show that the 2-pair excitation continuum yields subleading curvature corrections to the electron-electron scattering rate: , where V is the dimensionless strength of the interaction.     

Electron-doping versus hole-doping in the 2D t-t' Hubbard model

We compare the one-loop renormalization group flow to strong coupling of the electronic interactions in the two-dimensional t-t'-Hubbard model with t' = - 0.3t for band fillings smaller and larger than half-filling. Using a numerical N-patch scheme ( N = 32, ..., 96) we show that in the electron-doped case with decreasing electron density there is a rapid transition from a d _{x2 - y2}-wave superconducting regime with small characteristic energy scale to an approximate nesting regime with strong antiferromagnetic tendencies and higher energy scales. This contrasts with the hole-doped side discussed recently which exhibits a broad parameter region where the renormalization group flow suggests a truncation of the Fermi surface at the saddle points. We compare the quasiparticle scattering rates obtained from the renormalization group calculation which further emphasize the differences between the two cases. Received 19 December 2000 and Received in final form 28 February 2001     

An investigation of the 2D attractive Hubbard model

We present an investigation of the 2D attractive Hubbard model, considered as an effective model relevant to superconductivity in strongly interacting electron systems. We use both hybrid Monte-Carlo simulations and existing hopping parameter expansions to explore the low temperature domain. The increase of the static S-wave pair correlation with decreasing temperature, which depends weakly on the band filling in the explored temperature range, is analyzed in terms of an expected Kosterlitz-Thouless superconducting transition. Using both our data and previously published results, we show that the evidence for this transition is weak: If it exists, its temperature is very low. The number of unpaired electrons remains nearly constant with temperature at fixed attractive potential strength. In contrast, the static magnetic susceptibility decreases fast with temperature, and cannot be related only to pair formation. We introduce a method by which the Padé approximants of the existing series for the susceptibility give sensible results down to rather low temperature region, as shown by comparison with our numerical data. Received: 30 October 1996 / Revised: 23 October 1997 / Accepted: 29 January 1998