We study the ground state properties of the t–J model on two and three leg ladders. The phase diagram of the anisotropic couplings along the rungs and legs is investigated to find the opening of a charge transfer gap which appears as a charge density plateau. In a perturbative approach we start from isolated rungs at zero leg couplings. Then the stability of the charge density plateau is examined by first order perturbation theory in terms of an effective Hamiltonian using the numerical Lanczos method on finite size ladders. The results of this approach are compared with the exact diagonalization of finite systems. We have found that an improvement of our results depends on the initial clusters for zeroth order perturbation which cannot be recovered by higher order perturbation in rung clusters. For the two leg ladder case the application of 4-site (2 × 2) plaquettes helped us to overcome earlier flaws when the charge density is greater than one-half. We have also addressed the generalization of our method to trace phase separation in this model, where a linear dependence of ground state energy versus charge density gives evidence of a phase separated state.
We explore the ground-state phase diagram of the
S=1/2 two-leg ladder with different leg interactions. The
xy and
z components of the leg interactions between nearest-neighbor spins in the
a (
b) leg are respectively denoted by
Jl,a and
ΔlJl,a (
Jl,b and
ΔlJl,b). On the other hand, the
xy and
z components of the uniform rung interactions are respectively denoted by
ΓrJr and
Jr. In the above,
Δl and
Γr are the
XXZ-type anisotropy parameters for the leg and rung interactions, respectively. This system has a frustration when
Jl,aJl,b<0 irrespective of the sign of
Jr. The phase diagram on the
Δl (
∣Δl∣≤1.0) versus
Jl,b (
−2.0≤Jl,b≤3.0) plane in the case where
Jl,a=0.2,
Jr=−1.0, and
Γr=0.5 is determined numerically. We employ the physical consideration, and the level spectroscopy and phenomenological renormalization-group analyses of the numerical date obtained by the exact diagonalization method. The resultant phase diagram contains the ferromagnetic, Haldane, N{é}el, nematic Tomonaga-Luttinger liquid (TLL), partial ferrimagnetic, and
XY1 phases. Interestingly enough, the nematic TLL phase appears in the strong-rung unfrustrated region as well as in the strong-rung frustrated one. We perform the first-order perturbational calculations from the strong rung coupling limit to elucidate the characteristic features of the phase diagram. Furthermore, we make the density-matrix renormalization-group calculations for some physical quantities such as the energy gaps, the local magnetization, and the spin correlation functions to supplement the reliability of the phase diagram.
We have numerically investigated the doped t-J ladder using exact diagonalization. We have studied both the limit of strong interchain coupling and isotropic coupling. The ladder scales to the Luther-Emery liquid regime in the strong interchain coupling limit. In this strong coupling limit there is a simple picture of the excitation spectrum that can be continued to explain the behavior at isotropic coupling. At J=0 we have indications of a ferromagnetic ground state. At a large J/t the ladder is phase separated into holes and a Heisenberg ladder. At intermediate coupling the ground state shows hole pairing with a modified d-wave symmetry. The excitation spectrum separates into a limited number of quasiparticles which carry charge +\ensuremath{\Vert}e\ensuremath{\Vert} and spin 1/2 and a triplet magnon mode. At half filling the former vanish but the latter evolves continuously into the magnon band of the spin liquid. At low doping the quasiparticles form a dilute Fermi gas with a strong attraction but simultaneously the Fermi wave vector, as would be measured in photoemission, is large. The dynamical structure factors are calculated and are found to be very similar to calculations on two-dimensional clusters. \textcopyright{} 1996 The American Physical Society.
Using numerical results from density matrix renormalization group (DMRG) calculations for the t-J model, on systems as large as 10\ifmmode\times\else\texttimes\fi{}7, we examine the structure of the one and two hole ground states in ladder systems and in two dimensional clusters. A simple theoretical framework is used to explain why holes bind in pairs in two-dimensional antiferromagnets. For the case J/t=0.5, which we have studied, the hole pairs reside predominantly on a 2\ifmmode\times\else\texttimes\fi{}2 core plaquette with the probability that the holes are on diagonal sites greater than nearest-neighbor sites. There is a strong singlet bond connecting the spins on the two remaining sites of the plaquette. We find that a general characteristic of dynamic holes in an antiferromagnet is the presence of frustrating antiferromagnetic bonds connecting next-nearest-neighbor sites across the holes. Pairs of holes bind in order to share the frustrating bonds. At low doping, in addition to hole pairs, there are two additional low-energy structures which spontaneously form on certain finite systems. The first is an undoped L\ifmmode\times\else\texttimes\fi{}2 spin-liquid region, or ladder. The second is a hole moving along a one dimensional chain of sites. At higher doping we expect that hole pairing is always favored.
The ground-state phase diagram of a two-leg fermionic dipolar ladder with inter-site interactions is studied using density matrix renormalization group (DMRG) techniques. We use a state-of-the-art implementation of the DMRG algorithm and finite size scaling to simulate large system sizes with high accuracy. We also consider two different model systems and explore stable phases in half and quarter filling factors. We find that in the half filling, the charge and spin gaps emerge in a finite value of the dipole-dipole and on-site interactions. In the quarter filling case, s-wave superconducting state, charge density wave, homogenous insulating and phase separation phases occur depend on the interaction values. Moreover, in the dipole-dipole interaction, the D-Mott phase emerges when the hopping terms along the chain and rung are the same, whereas, this phase has been only proposed for the anisotropic Hubbard model. In the half filling case, on the other hand, there is either charge-density wave or charged Mott order phase depends on the orientation of the dipole moments of the particles with respect to the ladder geometry.