Quantum Monte Carlo Loop Algorithm for the t-J Model

We propose a generalization of the Quantum Monte Carlo loop algorithm to the t-J model by a mapping to three coupled six-vertex models. The autocorrelation times are reduced by orders of magnitude compared to the conventional local algorithms. The method is completely ergodic and can be formulated directly in continuous time. We introduce improved estimators for simulations with a local sign problem. Some first results of finite temperature simulations are presented for a t-J chain, a frustrated Heisenberg chain, and t-J ladder models.Comment: 22 pages, including 12 figures. RevTex v3.0, uses psf.te

Numerical renormalization-group study of spin correlations in one-dimensional random spin chains

We calculate the ground-state two-spin correlation functions of spin-1/2 quantum Heisenberg chains with random exchange couplings using the real-space renormalization group scheme. We extend the conventional scheme to take account of the contribution of local higher multiplet excitations in each decimation step. This extended scheme can provide highly accurate numerical data for large systems. The random average of staggered spin correlations of the chains with random antiferromagnetic (AF) couplings shows algebraic decay like $1/r^2$, which verifies the Fisher's analytic results. For chains with random ferromagnetic (FM) and AF couplings, the random average of generalized staggered correlations is found to decay more slowly than a power-law, in the form close to $1/\ln(r)$. The difference between the distribution functions of the spin correlations of the random AF chains and of the random FM-AF chains is also discussed.Comment: 14 pages including 8 figures, REVTeX, submitted to Physical Review

Lightly Doped t-J Three-Leg Ladders - an Analog for the Underdoped Cuprates

The three-leg ladder has one odd-parity and two even-parity channels. At low doping these behave quite differently. Numerical calculations for a t-J model show that the initial phase upon hole doping has two components - a conducting Luttinger liquid in the odd-parity channel, coexisting with an insulating (i.e. undoped) spin liquid phase in the even-parity channels. This phase has a partially truncated Fermi surface and violates the Luttinger theorem. This coexistence of conducting fermionic and insulating paired bosonic degrees of freedom is similar to the recent proposal of Geshkenbein, Ioffe, and Larkin for the underdoped spin-gap normal phase of the cuprates. A mean field approximation is derived which has many similarities to the numerical results. One difference however is an induced hole pairing in the odd-parity channel at arbitrary small dopings, similar to that proposed by Geshkenbein, Ioffe, and Larkin for the two-dimensional case. At higher dopings, we propose that a quantum phase transition will occur as holes enter the even-parity channels, resulting in a Luther-Emery liquid with hole pairing with essentially d-wave character. In the mean field approximation a crossover occurs which we interpret as a reflection of this quantum phase transition deduced from the numerical results.Comment: RevTex, 36 pages with 16 figure

Phase diagram and hidden order for generalized spin ladders

We investigate the phase diagram of antiferromagnetic spin ladders with additional exchange interactions on diagonal bonds by variational and numerical methods. These generalized spin ladders interpolate smoothly between the $S=1/2$ chain with competing nn and nnn interactions, the $S=1/2$ chain with alternating exchange and the antiferromagnetic $S=1$ chain. The Majumdar-Ghosh ground states are formulated as matrix product states and are shown to exhibit the same type of hidden order as the af $S=1$ chain. Generalized matrix product states are used for a variational calculation of the ground state energy and the spin and string correlation functions. Numerical (Lanczos) calculations of the energies of the ground state and of the low-lying excited states are performed, and compare reasonably with the variational approach. Our results support the hypothesis that the dimer and Majumdar-Ghosh points are in the same phase as the af $S=1$ chain.Comment: 23 pages, REVTEX, 7 figure

Monte Carlo Study of the Separation of Energy Scales in Quantum Spin 1/2 Chains with Bond Disorder

One-dimensional Heisenberg spin 1/2 chains with random ferro- and antiferromagnetic bonds are realized in systems such as $Sr_3 CuPt_{1-x} Ir_x O_6$. We have investigated numerically the thermodynamic properties of a generic random bond model and of a realistic model of $Sr_3 CuPt_{1-x} Ir_x O_6$ by the quantum Monte Carlo loop algorithm. For the first time we demonstrate the separation into three different temperature regimes for the original Hamiltonian based on an exact treatment, especially we show that the intermediate temperature regime is well-defined and observable in both the specific heat and the magnetic susceptibility. The crossover between the regimes is indicated by peaks in the specific heat. The uniform magnetic susceptibility shows Curie-like behavior in the high-, intermediate- and low-temperature regime, with different values of the Curie constant in each regime. We show that these regimes are overlapping in the realistic model and give numerical data for the analysis of experimental tests.Comment: 7 pages, 5 eps-figures included, typeset using JPSJ.sty, accepted for publication in J. Phys. Soc. Jpn. 68, Vol. 3. (1999

Spin Waves in Random Spin Chains

We study quantum spin-1/2 Heisenberg ferromagnetic chains with dilute, random antiferromagnetic impurity bonds with modified spin-wave theory. By describing thermal excitations in the language of spin waves, we successfully observe a low-temperature Curie susceptibility due to formation of large spin clusters first predicted by the real-space renormalization-group approach, as well as a crossover to a pure ferromagnetic spin chain behavior at intermediate and high temperatures. We compare our results of the modified spin-wave theory to quantum Monte Carlo simulations.Comment: 3 pages, 3 eps figures, submitted to the 47th Conference on Magnetism and Magnetic Material

Thermodynamics of Random Ferromagnetic Antiferromagnetic Spin-1/2 Chains

Using the quantum Monte Carlo Loop algorithm, we calculate the temperature dependence of the uniform susceptibility, the specific heat, the correlation length, the generalized staggered susceptibility and magnetization of a spin-1/2 chain with random antiferromagnetic and ferromagnetic couplings, down to very low temperatures. Our data show a consistent scaling behavior in all the quantities and support strongly the conjecture drawn from the approximate real-space renormalization group treatment.A statistical analysis scheme is developed which will be useful for the search of scaling behavior in numerical and experimental data of random spin chains.Comment: 13 pages, 13 figures, RevTe

Quantum magnetism in the stripe phase: bond- versus site order

It is argued that the spin dynamics in the charge-ordered stripe phase might be revealing with regards to the nature of the anomalous spin dynamics in cuprate superconductors. Specifically, if the stripes are bond ordered much of the spin fluctuation will originate in the spin sector itself, while site ordered stripes require the charge sector as the driving force for the strong quantum spin fluctuations.Comment: 4 pages, 3 figures, LaTe

Entropy Driven Dimerization in a One-Dimensional Spin-Orbital Model

We study a new version of the one-dimensional spin-orbital model with spins S=1 relevant to cubic vanadates. At small Hund's coupling J_H we discover dimerization in a pure electronic system solely due to a dynamical spin-orbital coupling. Above a critical value J_H, a uniform ferromagnetic state is stabilized at zero temperature. More surprisingly, we observe a temperature driven dimerization of the ferrochain, which occurs due to a large entropy released by dimer states. This dynamical dimerization seems to be the mechanism driving the peculiar intermediate phase of YVO_3.Comment: 5 pages, 4 figure

SU(4) Spin-Orbital Two-Leg Ladder, Square and Triangle Lattices

Based on the generalized valence bond picture, a Schwinger boson mean field theory is applied to the symmetric SU(4) spin-orbital systems. For a two-leg SU(4) ladder, the ground state is a spin-orbital liquid with a finite energy gap, in good agreement with recent numerical calculations. In two-dimensional square and triangle lattices, the SU(4) Schwinger bosons condense at (\pi/2,\pi/2) and (\pi/3,\pi/3), respectively. Spin, orbital, and coupled spin-orbital static susceptibilities become singular at the wave vectors, twice of which the bose condensation arises at. It is also demonstrated that there are spin, orbital, and coupled spin-orbital long-range orderings in the ground state.Comment: 5 page