We consider a Hubbard-like model of strongly interacting spinless fermions and hardcore bosons on a square lattice, such that nearest neighbor occupation is forbidden. Stripes (lines of holes across the lattice forming antiphase walls between ordered domains) are a favorable way to dope this system below half filling. The problem of a single stripe can be mapped to fermions in one dimension, which allows understanding of its elementary excitations and calculation of the stripe's effective mass for transverse vibrations. Using Lanczos exact diagonalization, we investigate the excitation gap and dispersion of a hole on a stripe, and the interaction of two holes. We also study the interaction of a few stripes, finding that they repel and the interaction energy decays with stripe separation as if they are hardcore particles moving in one (transverse) direction. To determine the stability of an array of stripes against phase separation into a particle-rich phase and a hole-rich liquid, we evaluate the liquid's equation of state, finding that the stripe-array is not stable for bosons but is possibly stable for fermions.
Spinless fermions on a lattice with nearest-neighbor repulsion serve as a toy-version Hubbard model and have a symmetry-broken even/odd superlattice at half filling. At infinite repulsion, doped holes form charged stripes that are antiphase walls ( as noted by Mila [Phys. Rev. B
49, 14 047 (1994)]). Exact-diagonalization data for systems up to 36 sites around 1/4 filling, and also for one or two holes added to a stripe of length up to 12, indicate stability of the stripe-array state against phase separation. In the boson version of the model, the same behavior can be stabilized by addition of a term with four creation/annihilation operators.
We investigate the phase diagram of spinless fermions on a square lattice with nearest-neighbor interaction, using the recently developed projective truncation approximation in Green's function equation of motion. For attractive interaction, the ground state is in a homogeneous
p+ip superconducting (SC) phase at high or low fermion densities. Near half filling is a phase separation (PS) between the SC phases. Allowing inhomogeneous solution, we obtain
p-wave SC domains with positive interface energy. As temperature increases, the SC phases transit into normal phases above
Tsc, generating a homogeneous normal phase (far away from
n=1/2), or a PS between normal phases with different densities (close to
n=1/2). Further increasing temperature to
Tps, the PS disappears and the particle-hole symmetry of the Hamiltonian is recovered. For repulsive interaction, depending on the filling, the ground state is in charge-ordered phase (half filling), charge-disordered phase (large hole/particle doping), or PS between them (weak doping). At finite temperature, the regime of charge order phase moves to finite
V and extends to finite doping regime.
A theoretical investigation of boson versions of the $t\ensuremath{-}J$ and $t\ensuremath{-}{J}_{z}$ models on the square lattice is carried out, by means of Green's function Monte Carlo simulations. Accurate ground-state energy estimates as a function of hole doping are obtained, allowing one to investigate the stability of the uniform phase against separation of the system into hole-rich and hole-free phases. In the boson $t\ensuremath{-}{J}_{z}$ model, such a separation is found to occur for arbitrarily small values of
Jz, at sufficiently low hole doping. Phase separation is suppressed in the boson $t\ensuremath{-}J$ model, which features a uniform ground state at any doping, for $J/t\ensuremath{\lesssim}1.5.$ Relevance of this study to the corresponding fermion models is discussed. Fermi statistics enhances the tendency toward phase separation; in particular, phase separation at low doping is predicted in the fermion $t\ensuremath{-}{J}_{z}$ model at any
Jz>0. The possible formation of stripes of holes is investigated for systems featuring both periodic and cylindrical boundary conditions. No evidence of a striped ground state is found in either the $t\ensuremath{-}J$ or $t\ensuremath{-}{J}_{z}$ boson models.
We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near
1/8 doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about
1/5 doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.
We consider an extended Hubbard model of interacting fermions on a lattice. The fermion kinetic energy corresponds to a tight-binding Hamiltonian with nearest-neighbor (-t) and next-nearest-neighbor (
t∣IH) hopping matrix elements. In addition to the on-site Hubbard interaction (U) we also consider a nearest-neighbor repulsion (V). We obtain the zero-temperature phase diagram of our model within the Hartree-Fock approximation. We consider ground states having charge- and spin-density-wave ordering as well as states with orbital antiferromagnetism or spin nematic order. The latter two states correspond to particle-hole binding with ${\mathrm{d}}_{{\mathrm{x}}^{2}\mathrm{\ensuremath{-}}{\mathrm{y}}^{2}}$ symmetry in the charge and spin channels, respectively. For
t∣IH=0, only the charge-density-wave and spin-density-wave states are energetically stable. For nonzero
t∣IH, we find that orbital antiferromagnetism (or spin nematic) order is stable over a finite portion of the phase diagram at weak coupling. This region of stability is seen to grow with increasing values of
t∣IH.