Search

We use series expansion methods to study ground- and excited-state properties in the helically ordered phase of spin-1/2 frustrated antiferromagnets on an anisotropic triangular lattice. We calculate the ground state energy, ordering wavevector, sublattice magnetization and one-magnon excitation spectrum for parameters relevant to Cs_2CuCl_4 and Cs_2CuBr_4. Both materials are modeled in terms of a Heisenberg model with spatially anisotropic exchange constants; for Cs_2CuCl_4 we also take into account the additional Dzyaloshinskii-Moriya (DM) interaction. We compare our results for Cs_2CuCl_4 with unpolarized neutron scattering experiments and find good agreement. In particular, the large quantum renormalizations of the one-magnon dispersion are well accounted for in our analysis, and inclusion of the DM interaction brings the theoretical predictions for the ordering wavevector and the magnon dispersion closer to the experimental results., Comment: 10 pages, 8 figures

We use series expansion methods to calculate the dispersion relation of the one-magnon excitations for the spin-1/2 triangular-lattice nearest-neighbor Heisenberg antiferromagnet above a three-sublattice ordered ground state. Several striking features are observed compared to the classical (large-S) spin-wave spectra. Whereas at low energies the dispersion is only weakly renormalized by quantum fluctuations, significant anomalies are observed at high energies. In particular, we find roton-like minima at special wave-vectors and strong downward renormalization in large parts of the Brillouin zone, leading to very flat or dispersionless modes. We present detailed comparison of our calculated excitation energies in the Brillouin zone with the spin-wave dispersion to order 1/S calculated recently by Starykh, Chubukov, and Abanov [cond-mat/0608002]. We find many common features but also some quantitative and qualitative differences. We show that at temperatures as low as 0.1J the thermally excited rotons make a significant contribution to the entropy. Consequently, unlike for the square lattice model, a non-linear sigma model description of the finite-temperature properties is only applicable at extremely low temperatures., Comment: 14 pages, 11 figures

We investigate the critical behaviour of the spectral weight of a single quasiparticle, one of the key observables in experiment, for the particular case of the transverse Ising model.Series expansions are calculated for the linear chain and the square and simple cubic lattices. For the chain model, a conjectured exact result is discovered. For the square and simple cubic lattices, series analyses are used to estimate the critical exponents. The results agree with the general predictions of Sachdev., Comment: 4 pages, 3 figures

We calculate the excitation spectrum and spectral weights of the alternating antiferromagnetic-ferromagnetic spin-half Heisenberg chain with exchange couplings $J$ and $-|\lambda|J$ as a power series in $\lambda$. For small $|\lambda|$, the gapped one-particle spectrum has a maximum at $k=0$ and there is a rich structure of bound (and anti-bound) states below (and above) the 2-particle continuum. As $|\lambda|$ is increased past unity the spectrum crosses over to the Haldane regime, where the peak shifts away from $k=0$, the one particle states merge with the bottom of the continuum near $k=0$, and the spectral weights associated with the one-particle states become very small. Extrapolation of the spectrum to large $|\lambda|$ confirms that the ground state energy and excitation gap map onto those of the spin-one chain., Comment: 4 pages, 8 figures

Motivated by the possible charge ordering in the recently discovered superconductor Na_xCoO_2 yH_2O, which at filling x=1/3 would correspond to a half-filled honeycomb lattice, we investigate the single-hole dynamics of the t-J model on this lattice. Using exact diagonalization, series expansion and the self-consistent Born approximation, we calculate the quasi-particle dispersion, bandwidth and residues and compare our findings with the well established results for the square lattice. Given the similarities between both lattices, it is not surprising that we find qualitatively similar features for the honeycomb lattice, namely an almost flat band along the edge of the magnetic Brillouin zone and well defined hole pockets around the corners. However, we report a considerable disagreement between the three methods used concerning the bandwidth and discuss possible origins of this discrepancy., Comment: 11 pages, 10 figures

We use high order linked cluster series to investigate the hard core boson model on the triangular lattice, at zero temperature. Our expansions, in powers of the hopping parameter $t$, probe the spatially ordered `solid' phase and the transition to a uniform superfluid phase. At the commensurate fillings $n=1/3, 2/3$ we locate a quantum phase transition point at $(t/V)_c\simeq 0.208(1)$, in good agreement with recent Monte Carlo studies. At half-filling ($n=1/2$) we find evidence for a solid phase, which persists to $t/V\simeq 0.06$., Comment: 13 pages, 7 figures

We use series expansions to study the excitation spectra of spin-1/2 antiferromagnets on anisotropic triangular lattices. For the isotropic triangular lattice model (TLM) the high-energy spectra show several anomalous features that differ strongly from linear spin-wave theory (LSWT). Even in the Neel phase, the deviations from LSWT increase sharply with frustration, leading to roton-like minima at special wavevectors. We argue that these results can be interpreted naturally in a spinon language, and provide an explanation for the previously observed anomalous finite-temperature properties of the TLM. In the coupled-chains limit, quantum renormalizations strongly enhance the one-dimensionality of the spectra, in agreement with experiments on Cs_2CuCl_4., Comment: 4 pages, 5 figures. New Fig. 3 with higher-order series data, paper shortened, references updated, one added (Ref. 28), minor changes otherwise. Published version

Spin and hole excitation spectra and spectral weights are calculated for the half-filled Hubbard model, as a function of $t/U$. We find that the high energy spin spectra are sensitive to charge fluctuations. The energy difference $\Delta(\pi,0)- \Delta(\pi/2,\pi/2)$, which is negative for the Heisenberg model, changes sign at a fairly small $t/U\approx 0.053(5)$. The hole bandwidth is proportional to $J$, and considerably larger than in the t-J models. It has a minimum at ($\pi/2,\pi/2$) and a very weak dispersion along the antiferromagnetic zone boundary. A good fit to the measured spin spectra in La$_2$CuO$_4$ at $T=10K$ is obtained with the parameter values $U=3.1{\rm eV}$, $t=0.35{\rm eV}$., Comment: 4 pages, 5 figures

We have extended our previous series studies of quantum antiferromagnets at zero temperature by computing the one-magnon dispersion curves and various structure factors for the linear chain, square and simple cubic lattices. Many of these results are new; others are a substantial extension of previous work. These results are directly comparable with neutron scattering experiments and we make such comparisons where possible., Comment: 15 pages, 12 figures, revised version

We present the temperature dependence of the uniform susceptibility of spin-half quantum antiferromagnets on spatially anisotropic triangular-lattices, using high temperature series expansions. We consider a model with two exchange constants, $J_1$ and $J_2$ on a lattice that interpolates between the limits of a square-lattice ($J_1=0$), a triangular-lattice ($J_2=J_1$), and decoupled linear chains ($J_2=0$). In all cases, the susceptibility which has a Curie-Weiss behavior at high temperatures, rolls over and begins to decrease below a peak temperature, $T_p$. Scaling the exchange constants to get the same peak temperature, shows that the susceptibilities for the square-lattice and linear chain limits have similar magnitudes near the peak. Maximum deviation arises near the triangular-lattice limit, where frustration leads to much smaller susceptibility and with a flatter temperature dependence. We compare our results to the inorganic materials Cs$_2$CuCl$_4$ and Cs$_2$CuBr$_4$ and to a number of organic molecular crystals. We find that the former (Cs$_2$CuCl$_4$ and Cs$_2$CuBr$_4$) are weakly frustrated and their exchange parameters determined through the temperature dependence of the susceptibility are in agreement with neutron-scattering measurements. In contrast, the organic materials are strongly frustrated with exchange parameters near the isotropic triangular-lattice limit., Comment: 10 pages, 9 figures and 1 table, revised version

We estimate, using high-temperature series expansions, the transition temperatures of the spin 1/2, 1 and 3/2 Heisenberg ferromagnet and antiferromagnet in 3-dimensions. The manner in which the difference between Curie and Neel temperatures vanishes with increasing spin quantum number is investigated., Comment: 10 pages

We use linked-cluster series expansions, both at T=0 and high temperature, to analyse the phase structure of the spin-$\half$ Heisenberg antiferromagnet with competing first and second-neighbor interactions on the 3-dimensional body-centred-cubic lattice. At zero temperature we find a first-order quantum phase transition at $J_2/J_1 \simeq 0.705 \pm 0.005$ between AF$_1$ (Ne\'el) and AF$_2$ ordered phases. The high temperature series yield quite accurate estimates of the bounding critical line for the AF$_1$ phase, and an apparent critical line for the AF$_2$ phase, with a bicritical point at $J_1/J_2\simeq 0.71$, $kT/J_1\simeq 0.34$. The possibility that this latter transition is first-order cannot be excluded., Comment: 10 pages, 4 figures

Motivated by the geometry of spins in the material CaCu$_2$O$_3$, we study a two-layer, spin-half Heisenberg model, with nearest-neighbor exchange couplings J and \alpha*J along the two axes in the plane and a coupling J_\perp perpendicular to the planes. We study these class of models using the Stochastic Series Expansion (SSE) Quantum Monte Carlo simulations at finite temperatures and series expansion methods at T=0. The critical value of the interlayer coupling, J_\perp^c, separating the N{\'e}el ordered and disordered ground states, is found to follow very closely a square root dependence on $\alpha$. Both T=0 and finite-temperature properties of the model are presented., Comment: 9 pages, 11 figs., 1 table

We develop a linked cluster method to calculate the spectral weights of many-particle excitations at zero temperature. The dynamical structure factor is expressed as a sum of exclusive structure factors, each representing contributions from a given set of excited states. A linked cluster technique to obtain high order series expansions for these quantities is discussed. We apply these methods to the alternating Heisenberg chain around the dimerized limit ($\lambda=0$), where complete wavevector and frequency dependent spectral weights for one and two-particle excitations (continuum and bound-states) are obtained. For small to moderate values of the inter-dimer coupling parameter $\lambda$, these lead to extremely accurate calculations of the dynamical structure factors. We also examine the variation of the relative spectral weights of one and two-particle states with bond alternation all the way up to the limit of the uniform chain ($\lambda=1$). In agreement with Schmidt and Uhrig, we find that the spectral weight is dominated by 2-triplet states even at $\lambda=1$, which implies that a description in terms of triplet-pair excitations remains a good quantitative description of the system even for the uniform chain., Comment: 26 pages, 17 figures

Linked cluster series expansions about the Ising limit are used to study ground state preperties, viz. ground state energy, magnetization and excitation spectra, for mixed spin S=(1/2,1) quantum ferrimagnets on simple bipartite lattices in 1, 2, and 3-dimensions. Results are compared to second-order spin wave theory and, in general, excellent agreement is obtained., Comment: 5 figures, 12 pages, revtex

We use perturbative series expansions about a staggered dimerized ground state to compute the ground state energy, triplet excitation spectra and spectral weight for a one-dimensional model in which each site has an $S=\case 1/2$ spin ${\bf S}_i$ and a pseudospin ${\bf T}_i$, representing a doubly degenerate orbital. An explicit dimerization is introduced to allow study of the confinement of spinon excitations. The elementary triplet represents a bound state of two spinons, and is stable over much of the Brillouine zone. A special line is found in the gapped spin-liquid phase, on which the triplet excitation is dispersionless. The formation of triplet bound states is also investigated., Comment: 9 pages, 9 figures

We develop strong-coupling series expansion methods to study two-particle spectra of quantum lattice models. At the heart of the method lies the calculation of an effective Hamiltonian in the two-particle subspace. We explicitly consider an orthogonality transformation to generate this block diagonalization, and find that maintaining orthogonality is crucial for systems where the ground state and the two-particle subspace are characterized by identical quantum numbers. We discuss the solution of the two-particle Schr\"odinger equation by using a finite lattice approach in coordinate space or by an integral equation in momentum space. These methods allow us to precisely determine the low-lying excitation spectra of the models at hand, including all two-particle bound/antibound states. Further, we discuss how to generate series expansions for the dispersions of the bound/antibound states. These allow us to employ series extrapolation techniques, whereby binding energies can be determined even when the expansion parameters are not small. We apply the method to the (1+1)D transverse Ising model and the two-leg spin-$\case 1/2$ Heisenberg ladder. For the latter model, we also calculate the coherence lengths and determine the critical properties where bound states merge with the two-particle continuum., Comment: 10 figures, revtex

We use recently developed strong-coupling expansion methods to study the two-particle spectra for the frustrated alternating Heisenberg model, consisting of an alternating nearest neighbor antiferromagnetic exchange and a uniform second neighbor antiferromagnetic exchange. Starting from the limit of weakly coupled dimers, we develop high order series expansions for the effective Hamiltonian in the two-particle subspace. In the limit of a strong applied dimerization, we calculate accurately various properties of singlet and triplet bound states and quintet antibound states. We also develop series expansions for bound state energies in various sectors, which can be extrapolated using standard methods to cases where the external bond-alternation goes to zero. We study the properties of singlet and triplet bound states in the latter limit and suggest a crucial role for the bound states in the unbinding of triplets and deconfinement of spin-half excitations., Comment: 17 figures, revtex