Using the standard Quantum Monte Carlo technique for the Hubbard model, I present here a numerical investigation of the hole propagation in a Quantum Antiferromagnet. The calculation is very well stabilized, using selected sized systems and special use of the trial wavefunction that satisfy the “close shell condition” in presence of an arbitrarily weak Zeeman magnetic field, vanishing in the thermodynamic limit. It will be shown in a forthcoming publication1 that the presence of this magnetic field does not affect thermodynamic properties for the half filled system. Then I have used the same selected sizes for the one hole ground state. I have investigated the question of vanishing or nonvanishing quasiparticle weight, in order to clarify whether the Mott insulator should behave just as conventional insulator with an upper and lower Hubbard band. By comparing the present finite size scaling with several techniques predicting a finite quasiparticle weight (see Fig.1) the data seem more consistent with a vanishing quasiparticle weight, i.e. , as recently suggested by P.W. Anderson2 the Hubbard-Mott insulator should be characterized by non-trivial excitations which cannot be interpreted in a simple quasi-particle picture. However it cannot be excluded , based only on numerical grounds, that a very small but non vanishing quasiparticle weight should survive in the thermodynamic limit.
Using the standard Quantum Monte Carlo technique for the Hubbard model, I present here a numerical investigation of the hole propagation in a Quantum Antiferromagnet. The calculation is very well stabilized, using selected sized systems and special use of the trial wavefunction that satisfy the “close shell condition” in presence of an arbitrarily weak Zeeman magnetic field, vanishing in the thermodynamic limit. It will be shown in a forthcoming publication1 that the presence of this magnetic field does not affect thermodynamic properties for the half filled system. Then I have used the same selected sizes for the one hole ground state. I have investigated the question of vanishing or nonvanishing quasiparticle weight, in order to clarify whether the Mott insulator should behave just as conventional insulator with an upper and lower Hubbard band. By comparing the present finite size scaling with several techniques predicting a finite quasiparticle weight (see Fig.1) the data seem more consistent with a vanishing quasiparticle weight, i.e. , as recently suggested by P.W. Anderson2 the Hubbard-Mott insulator should be characterized by non-trivial excitations which cannot be interpreted in a simple quasi-particle picture. However it cannot be excluded , based only on numerical grounds, that a very small but non vanishing quasiparticle weight should survive in the thermodynamic limit.
Using the standard quantum Monte Carlo technique for the Hubbard model, I present here a numerical investigation of the hole propagation in a quantum antiferromagnet. The calculation is very well stabilized, using selected-size systems and a special choice of the trial wave function that satisfies the ``closed-shell condition'' in the presence of an arbitrarily weak Zeeman magnetic field. First it is shown that the presence of this magnetic field does not affect the thermodynamic properties of the half-filled system. Then I use the same selected sizes for the one-hole ground state. I investigate the question of vanishing or nonvanishing quasiparticle weight, in order to clarify whether the Mott insulator behaves as a conventional insulator with an upper and lower Hubbard band. By comparing the present finite-size scaling with several techniques that predict a finite quasiparticle weight, the data seem more consistent with a vanishing quasiparticle weight, i.e., the Mott-Hubbard insulator should be characterized by nontrivial excitations that cannot be interpreted in a simple picture. However, one cannot exclude, based only on numerical grounds, a very small but nonvanishing quasiparticle weight. In the last part I give some theoretical arguments to explain the results of the Monte Carlo simulation.
We study a ground-state ansatz for the single-hole-doped $t\text{\ensuremath{-}}J$ model in two dimensions via a variational Monte Carlo method. Such a single-hole wave function possesses finite angular momenta generated by hidden spin currents, which give rise to a novel ground-state degeneracy in agreement with recent exact diagonalization (ED) and density matrix renormalization group (DMRG) results. We further show that the wave function can be decomposed into a quasiparticle component and an incoherent momentum distribution in excellent agreement with the DMRG results up to an $8\ifmmode\times\else\texttimes\fi{}8$ lattice. Such a two-component structure indicates the breakdown of Landau's one-to-one correspondence principle, and in particular, the quasiparticle spectral weight vanishes by a power law in the large sample size limit. By contrast, turning off the phase string induced by the hole hopping in the so-called $\ensuremath{\sigma}\ifmmode\cdot\else\textperiodcentered\fi{}t\text{\ensuremath{-}}J$ model, a conventional Bloch-wave wave function with a finite quasiparticle spectral weight can be recovered, also in agreement with the ED and DMRG results. The present study shows that a singular effect already takes place in the single-hole-doped Mott insulator, by which the bare hole is turned into a non-Landau quasiparticle with translational-symmetry breaking. Generalizations to pairing and finite doping are briefly discussed.
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/t from 2 to 8, for a range of densities including half-filling and n=0.3,0.5,0.6,0.75, and 0.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.
An analytic approach is presented of electronic one-particle spectra of the one-band Hubbard model at half filling in the antiferromagnetic phase. Starting from the strong-coupling regime U\ensuremath{\gg}:t, a projection technique is used to set up self-consistent coupled equations for the electron Green function, which are valid down to values U\ensuremath{\approx}t. The self-consistent equation for the hole propagator is a direct generalization of the one found from the t-J model. This gives further support to the ``string'' picture, where propagation of holes creates strings of overturned spins with which the holes interact. In the present work hopping of holes (or electrons) with up spin on the down sublattice is also taken into acount, as well as transitions between the lower and upper Hubbard bands. These are shown to change significantly the incoherent part of the t-J model spectra, by smearing out the shake-off peaks, reminiscent of higher bound string states due to multispin scattering. Coherent (quasiparticle) peaks exist at the band edges, on both sides of the insulating gap. We show that with decreasing U the quasiparticle concept loses its meaning for wave vectors at the center of the magnetic Brillouin zone (MBZ). For large values of U the dispersion of the quasiparticle is found to scale with its band width, which is of order J. Extrema are always found at k=(\ensuremath{\pi}/2,\ensuremath{\pi}/2). The weight of the quasiparticle at this k value decreases logarithmically with increasing U. In the strong-coupling limit the spectrum tends to be symmetric, i.e., to become an even function of the frequency around the chemical potential, for any wave vector. For small values of U the dispersion at the edge of the MBZ flattens away, as expected when approaching the noninteracting limit. The spectral function in this regime, for wave vectors away from the edge of the MBZ, is concentrated mainly on one side of the chemical potential.
Solid state physics, a very broad field in natural science, aims to explain physical properties of materials modeled as a many-body ensemble of the order of ∼ 10^23 interacting particles. For interacting metallic systems, one of the most powerful tools for unraveling the consequences of the collective interaction is the Fermi liquid theory. While still being an approximation, it excellently models many materials especially in the low-temperature and low-energy regime by treating themany-body problem as being equivalent to effectively non-interacting ”particles”, which is done by introducing the concept of quasiparticles. During the last few decades, several exotic phases of matter have been discovered, which originate from strong Coulomb interactions between the electrons, whose description is beyond the Fermi liquid theory. This thesis addresses the question of the validity of Fermi liquid theory once a single hole is introduced into an antiferromagnet with completely localized spins due to the strong interactions. Here, the most-studied case concerns the so-called 2D t − J model, which hosts a well-defined quasiparticle, the spin polaron. Focusing on a more tractable case, the t − Jz model without quantum spin fluctuations, we study the stability of the spin polaron solution once the nextnearest neighbor hopping t′ is introduced into the system. We focus on observations at the Γ-point and the X-point of systems with fixed Jz = 0.4t and varying t′ ∈ [−0.5t, 0.5t]. The model is numerically solved using exact diagonalization on finite-size lattices of up to 32 sites. A study of the quasiparticle weight under finite-size scaling reveals an “anomalous” regime at the Γ-point fort′ ≤ −0.3t, indicated by a vanishing quasiparticle weight. In order to characterize the low-energy dynamics within this “anomalous” regime, we investigated the spin correlation centered around the hole, as well as the magnon number distribution of a single-doped 20-site system. Based on the results, we present a preliminary picture of the antiferromagnetic ground state doped with one hole within that “anomalous” regime (t′ ≤ −0.3t) that deviates from the standard spin polaron as the cloud of magnetic fluctuations, typically localized around the hole, delocalizes over the whole finite-size lattice. Consequently, such systems show a weakly suppressed antiferromagnetic configuration under single-hole doping.
Dynamic properties of the three-dimensional single-band Hubbard model are studied using quantum Monte Carlo combined with the maximum entropy technique. At half-filling, there is a clear gap in the density of states and well-defined quasiparticle peaks at the top (bottom) of the lower (upper) Hubbard band. We find an antiferromagnetically induced weight above the naive Fermi momentum. Upon hole doping, the chemical potential \ensuremath{\mu} moves to the top of the lower band where a robust peak is observed. Results are compared with spin-density-wave mean-field and self-consistent Born approximation results, and also with the infinite-dimensional (D=\ensuremath{\infty}) Hubbard model, and experimental photoemission for three-dimensional transition-metal oxides. \textcopyright{} 1996 The American Physical Society.