Search

Search for a new quantum state of matter emerging in a crystal is one of recent trends in condensed matter physics. For magnetic materials, geometrical frustration and high magnetic field are two key ingredients to realize it: a conventional magnetic order is possibly destroyed by competing interactions (frustration) and is replaced by an exotic state that is characterized in terms of quasiparticles, that are magnons, and the magnetic field can control the density and chemical potential of the magnons. Here we show that a synthetic copper mineral, Cd-kapellasite, comprising a kagome lattice made of corner-sharing triangles of Cu2+ ions carrying spin-1/2 exhibits an unprecedented series of fractional magnetization plateaux in ultrahigh magnetic fields up to 160 T, which may be interpreted as crystallizations of emergent magnons localized on the hexagon of the kagome lattice. Our observation reveals a novel type of particle physics realized in a highly frustrated magnet., Comment: 9 pages, 10 figures

We present results of our Monte Carlo simulation of the Ising-O(3) model on the two-dimensional (2D) and quasi-2D lattices. This model is an effective classical model for the stacked square-lattice $J_1$-$J_2$ Heisenberg model where the nearest neighbor ($J_1$) and next-nearest neighbor ($J_2$) couplings are frustrated and we assume that $J_2$ is dominant. We find an Ising ordered phase where the O(3) spins remain disordered in a moderate quasi-2D region. There is a single first order transition for a sufficiently large 3D coupling. Our results can provide a qualitative explanation on the experiments on ferropnictides, namely observed sequence and orders of the structural and magnetic transitions, in terms of the ratio between the inter-layer and intra-layer coupling., Comment: 12 pages, 11 figures; Corrected minor misprints; published in Phys. Rev. B

We use quantum Monte Carlo (QMC) simulations to study the combined effects of harmonic confinement and temperature for bosons in a two dimensional optical lattice. The scale invariant, finite temperature, state diagram is presented for the Bose-Hubbard model in terms of experimental parameters -- the particle number, confining potential and interaction strength. To distinguish the nature of the spatially separated superfluid, Mott Insulator and normal Bose liquid phases, we examine the local density, compressibility, superfluid density and Green's function. In the annular superfluid rings, as the width of the ring decreases, the long range superfluid correlations start to deviate from an equivalent homogeneous 2D system. At zero temperature, the correlation decay is intermediate between 1D and 2D, while at finite temperature, the decay is similar to that in 1D at a much lower temperature. The calculations reveal shortcomings of the local density approximation (LDA) in describing superfluid properties of trapped bosons. We also present the finite temperature phase diagram for the homogeneous two dimensional Bose-Hubbard model. We compare our state diagram with the results of a recent experiment at NIST on a harmonically trapped 2D lattice [Phys. Rev. Lett. 105, 110401 (2010)], and identify a finite temperature effect in the experiment., Comment: 13 pages, 9 figures

We study the finite-temperature transition to the $m=1/2$ magnetization plateau in a model of interacting $S=1/2$ spins with longer range interactions and strong exchange anisotropy on the geometrically frustrated Shastry-Sutherland lattice. This model was shown to capture the qualitative features of the field-induced magnetization plateaus in the rare-earth tetraboride, ${\rm TmB_4}$. Our results show that the transition to the plateau state occurs via two successive transitions with the two-dimensional Ising universality class, when the quantum exchange interactions are finite, whereas a single phase transition takes place in the purely Ising limit. To better understand these behaviors, we perform Monte Carlo simulations of the classical generalized four-state chiral clock model and compare the phase diagrams of the two models. Finally, we estimate a parameter set that can explain the magnetization curves observed in ${\rm TmB_4}$. The magnetic properties and critical behavior of the finite-temperature transition to the $m=1/2$ plateau state are also discussed., Comment: 11 pages, 15 figures

We study two-dimensional Heisenberg antiferromagnets with additional multi-spin interactions which can drive the system into a valence-bond solid state. For standard SU(2) spins, we consider both four- and six-spin interactions. We find continuous quantum phase transitions with the same critical exponents. Extending the symmetry to SU(N), we also find continuous transitions for N=3 and 4. In addition, we also study quantitatively the cross-over of the order-parameter symmetry from Z4 deep inside the valence-bond-solid phase to U(1) as the phase transition is approached., Comment: 4 pages, 5 figures

We report a method for determining the spatial dependence of the magnetic exchange coupling, $dJ/dr$, from magnetostriction measurements of a quantum magnet. The organic Ni $S = 1$ system NiCl$_2$-4SC(NH$_2$)$_2$ exhibits lattice distortions in response to field-induced canted antiferromagnetism between $H_{c1} = 2.1$ T and $H_{c2} = 12.6$ T. We are able to model the magnetostriction in terms of uniaxial stress on the sample created by magnetic interactions between neighboring Ni atoms along the c-axis. The uniaxial strain is equal to $(1/E)dJ_c/dx_c < S_{\bf r} \cdot S_{{\bf r}+ {\bf e}_c} >$, where $E$, $J_c$, $x_c$ and ${\bf e}_c$ are the Young's modulus, the nearest neighbor (NN) exchange coupling, the variable lattice parameter, and the relative vector between NN sites along the c-axis. We present magnetostriction data taken at 25 mK together with Quantum Monte Carlo calculations of the NN spin-spin correlation function that are in excellent agreement with each other. We have also measured Young's modulus using resonant ultrasound, and we can thus extract $dJ_c/dx_c = 2.5$ K/$\AA$, yielding a total change in $J_c$ between $H_{c1}$ and $H_{c2}$ of 5.5 mK or 0.25% in response to an 0.022% change in length of the sample.

NiCl$_2$-4SC(NH$_2$)$_2$ (DTN) is a quantum S=1 chain system with strong easy-pane anisotropy and a new candidate for the Bose-Einstein condensation of the spin degrees of freedom. ESR studies of magnetic excitations in DTN in fields up to 25 T are presented. Based on analysis of the single-magnon excitation mode in the high-field spin-polarized phase and previous experimental results [Phys. Rev. Lett. 96, 077204 (2006)], a revised set of spin-Hamiltonian parameters is obtained. Our results yield $D=8.9$ K, $J_c=2.2$ K, and $J_{a,b}=0.18$ K for the anisotropy, intrachain, and interchain exchange interactions, respectively. These values are used to calculate the antiferromagnetic phase boundary, magnetization and the frequency-field dependence of two-magnon bound-state excitations predicted by theory and observed in DTN for the first time. Excellent quantitative agreement with experimental data is obtained.

We show that the spatial dimensionality of the quantum critical point associated with Bose--Einstein condensation at T=0 is reduced when the underlying lattice comprises a set of layers coupled by a frustrating interaction. Our theoretical predictions for the critical temperature as a function of the chemical potential correspond very well with recent measurements in BaCuSi$_{2}$O$_{6}$ [S. E. Sebastian \textit{et al}, Nature \textbf{411}, 617 (2006)]., Comment: 5 pages, 2 figures

Competition between electronic ground states near a quantum critical point (QCP) - the location of a zero-temperature phase transition driven solely by quantum-mechanical fluctuations - is expected to lead to unconventional behaviour in low-dimensional systems. New electronic phases of matter have been predicted to occur in the vicinity of a QCP by two-dimensional theories, and explanations based on these ideas have been proposed for significant unsolved problems in condensed-matter physics, such as non-Fermi-liquid behaviour and high-temperature superconductivity. But the real materials to which these ideas have been applied are usually rendered three-dimensional by a finite electronic coupling between their component layers; a two-dimensional QCP has not been experimentally observed in any bulk three-dimensional system, and mechanisms for dimensional reduction have remained the subject of theoretical conjecture. Here we show evidence that the Bose-Einstein condensate of spin triplets in the three-dimensional Mott insulator BaCuSi2O6 provides an experimentally verifiable example of dimensional reduction at a QCP. The interplay of correlations on a geometrically frustrated lattice causes the individual two-dimensional layers of spin-1/2 Cu2+ pairs (spin dimers) to become decoupled at the QCP, giving rise to a two-dimensional QCP characterized by power law scaling distinctly different from that of its three-dimensional counterpart. Thus the very notion of dimensionality can be said to acquire an 'emergent' nature: although the individual particles move on a three-dimensional lattice, their collective behaviour occurs in lower-dimensional space., Comment: 14 pages, 4 figures

We report an experimental determination of the phase boundary between a quantum paramagnetic state and the proposed spin Bose-Einstein condensate of triplons in the spin gap compound BaCuSi2O6. The ordering temperature is related to the proximity to a quantum critical point at the lower critical magnetic field H_c1 = 23.52 +/- 0.03T by a power law parameterized by critical exponent nu. We obtain an experimental estimate of nu = 0.63 +/- 0.03 which is in good agreement with the mean field prediction of nu = 2/3 for the 3D XY model, used to describe the Bose condensation of a 3D dilute interacting Bose gas., Comment: 4 pages, 3 figures; added reference

Besides being an ancient pigment, BaCuSi$_2$O$_6$ is a quasi-2D magnetic insulator with a gapped spin dimer ground state. The application of strong magnetic fields closes this gap creating a gas of bosonic spin triplet excitations called triplons. The topology of the spin lattice makes BaCuSi$_2$O$_6$ an ideal candidate for studying the Bose-Einstein condensation of triplons as a function of the external magnetic field, which acts as a chemical potential. In agreement with quantum Monte Carlo numerical simulations, we observe a distinct lambda-anomaly in the specific heat together with a maximum in the magnetic susceptibility upon cooling down to liquid Helium temperatures., Comment: published on August 20, 2004

We review recent findings on spin glass models. Both the equilibrium properties and the dynamic properties are covered. We focus on progress in theoretical, in particular numerical, studies, while its relationship to real magnetic materials is also mentioned., Comment: Chapter 6 in ``Frustrated Spin Systems'' edited by H.T.Diep

This is a reply to the comment by A. W. Sandvik (cond-mat/0010433) on our paper Phys. Rev. Lett. 84, 4204 (2000). We show that his data do not conflict with our data nor with our conclusions., Comment: RevTeX, 1 page; Revised version

Critical behavior of the quantum phase transition of a site-diluted Heisenberg antiferromagnet on a square lattice is investigated by means of the quantum Monte Carlo simulation with the continuous-imaginary-time loop algorithm. Although the staggered spin correlation function decays in a power law with the exponent definitely depending on the spin size $S$, the correlation-length exponent is classical, i.e., $\nu=4/3$. This implies that the length scale characterizing the non-universal quantum phase transition is nothing but the mean size of connected spin clusters., Comment: 4 pages, 3 figures

We study the critical behavior of Ising quantum magnets with broadly distributed random couplings (J), such that $P(\ln J) \sim |\ln J|^{-1-\alpha}$, $\alpha>1$, for large $|\ln J|$ (L\'evy flight statistics). For sufficiently broad distributions, $\alpha<\alpha_c$, the critical behavior is controlled by a line of fixed points, where the critical exponents vary with the L\'evy index, $\alpha$. In one dimension, with $\alpha_c=2$, we obtaind several exact results through a mapping to surviving Riemann walks. In two dimensions the varying critical exponents have been calculated by a numerical implementation of the Ma-Dasgupta-Hu renormalization group method leading to $\alpha_c \approx 4.5$. Thus in the region $2<\alpha<\alpha_c$, where the central limit theorem holds for $|\ln J|$ the broadness of the distribution is relevant for the 2d quantum Ising model., Comment: 10pages, 13figures, final form

The quantum critical behavior and the Griffiths-McCoy singularities of random quantum Ising ferromagnets are studied by applying a numerical implementation of the Ma-Dasgupta-Hu renormalization group scheme. We check the procedure for the analytically tractable one-dimensional case and apply our code to the quasi-one-dimensional double chain. For the latter we obtain identical critical exponents as for the simple chain implying the same universality class. Then we apply the method to the two-dimensional case for which we get estimates for the exponents that are compatible with a recent study in the same spirit., Comment: 10 pages LaTeX, eps-figures and PTP-macros included. Proceedings of the ICCP5, Kanazawa (Japan), 1999

Ground-state magnetic properties of the diluted Heisenberg antiferromagnet on a square lattice are investigated by means of the quantum Monte Carlo method with the continuous-time loop algorithm. It is found that the critical concentration of magnetic sites is independent of the spin size S, and equal to the two-dimensional percolation threshold. However, the existence of quantum fluctuations makes the critical exponents deviate from those of the classical percolation transition. Furthermore, we found that the transition is not universal, i.e., the critical exponents significantly depend on S., Comment: RevTeX, 4 pages including 5 EPS figures

We study the quantum phase transition in the two-dimensional random Ising model in a transverse field by Monte Carlo simulations. We find results similar to those known analytically in one-dimension. At the critical point, the dynamical exponent is infinite and the typical correlation function decays with a stretched exponential dependence on distance. Away from the critical point there are Griffiths-McCoy singularities, characterized by a single, continuously varying exponent, z', which diverges at the critical point, as in one-dimension. Consequently, the zero temperature susceptibility diverges for a RANGE of parameters about the transition., Comment: 4 pages RevTeX, 3 eps-figures included

A cluster algorithm formulated in continuous (imaginary) time is presented for Ising models in a transverse field. It works directly with an infinite number of time-slices in the imaginary time direction, avoiding the necessity to take this limit explicitly. The algorithm is tested at the zero-temperature critical point of the pure two-dimensional (2d) transverse Ising model. Then it is applied to the 2d Ising ferromagnet with random bonds and transverse fields, for which the phase diagram is determined. Finite size scaling at the quantum critical point as well as the study of the quantum Griffiths-McCoy phase indicate that the dynamical critical exponent is infinite as in 1d., Comment: 4 pages RevTeX, 3 eps-figures included

With the help of EXACT ground states obtained by a polynomial algorithm we compute the domain wall energy at zero-temperature for the bond-random and the site-random Ising spin glass model in two dimensions. We find that in both models the stability of the ferromagnetic AND the spin glass order ceases to exist at a UNIQUE concentration p_c for the ferromagnetic bonds. In the vicinity of this critical point, the size and concentration dependency of the first AND second moment of the domain wall energy are, for both models, described by a COMMON finite size scaling form. Moreover, below this concentration the stiffness exponent turns out to be slightly negative \theta_S = -0.056(6) indicating the absence of any intermediate spin glass phase at non-zero temperature., Comment: 7 pages Latex, 5 postscript-figures included

We have studied the three-dimensional Ising spin glass with a $\pm J$ distribution by Monte Carlo simulations. Using larger sizes and much better statistics than in earlier work, a finite size scaling analysis shows quite strong evidence for a finite transition temperature, $T_c$, with ordering below $T_c$. Our estimate of the transition temperature is rather lower than in earlier work, and the value of the correlation length exponent, $\nu$, is somewhat higher. Because there may be (unknown) corrections to finite size scaling, we do not completely rule out the possibility that $T_c = 0$ or that $T_c$ is finite but with no order below $T_c$. However, from our data, these possibilities seem less likely., Comment: Postscript file compressed using uufiles. The postscript file is also available by anonymous ftp at ftp://chopin.ucsc.edu/pub/sg3d.ps

We present cluster Monte Carlo algorithms for the $XYZ$ quantum spin models. In the special case of $S=1/2$, the new algorithm can be viewed as a cluster algorithm for the 8-vertex model. As an example, we study the $S=1/2$ $XY$ model in two dimensions with a representation in which the quantization axis lies in the easy plane. We find that the numerical autocorrelation time for the cluster algorithm remains of the order of unity and does not show any significant dependence on the temperature, the system size, or the Trotter number. On the other hand, the autocorrelation time for the conventional algorithm strongly depends on these parameters and can be very large. The use of improved estimators for thermodynamic averages further enhances the efficiency of the new algorithms., Comment: 20 pages including 3 figures and 3 tables, accepted for publication in J. Stat. Phys

We generalize the Fortuin-Kasteleyn (FK) cluster representation of the partition function of the Ising model to represent the partition function of quantum spin models with an arbitrary spin magnitude in arbitrary dimensions. This generalized representation enables us to develop a new cluster algorithm for the simulation of quantum spin systems by the worldline Monte Carlo method. Because the Swendsen-Wang algorithm is based on the FK representation, the new cluster algorithm naturally includes it as a special case. As well as the general description of the new representation, we present an illustration of our new algorithm for some special interesting cases: the Ising model, the antiferromagnetic Heisenberg model with $S=1$, and a general Heisenberg model. The new algorithm is applicable to models with any range of the exchange interaction, any lattice geometry, and any dimensions., Comment: 46 pages, 10 figures, to appear in J.Stat.Phys

We discuss the development of cluster algorithms from the viewpoint of probability theory and not from the usual viewpoint of a particular model. By using the perspective of probability theory, we detail the nature of a cluster algorithm, make explicit the assumptions embodied in all clusters of which we are aware, and define the construction of free cluster algorithms. We also illustrate these procedures by rederiving the Swendsen-Wang algorithm, presenting the details of the loop algorithm for a worldline simulation of a quantum $S=$ 1/2 model, and proposing a free cluster version of the Swendsen-Wang replica method for the random Ising model. How the principle of maximum entropy might be used to aid the construction of cluster algorithms is also discussed., Comment: 25 pages, 4 figures, to appear in Phys.Rev.E

We introduce two novel quantum Monte Carlo methods and employ them to study the superfluid-insulator transition in a two-dimensional system of hard-core bosons. One of the methods is appropriate for zero temperature and is based upon Green's function Monte Carlo; the other is a finite-temperature world-line cluster algorithm. In each case we find that the dynamical exponent is consistent with the theoretical prediction of $z=2$ by Fisher and co-workers., Comment: Revtex, 10 pages, 3 figures (postscript files attached at end, separated by %%%%%% Fig # %%%%%, where # is 1-3). LA-UR-94-2708

Two cluster algorithms, based on constructing and flipping loops, are presented for worldline quantum Monte Carlo simulations of fermions and are tested on the one-dimensional repulsive Hubbard model. We call these algorithms the loop-flip and loop-exchange algorithms. For these two algorithms and the standard worldline algorithm, we calculated the autocorrelation times for various physical quantities and found that the ordinary worldline algorithm, which uses only local moves, suffers from very long correlation times that makes not only the estimate of the error difficult but also the estimate of the average values themselves difficult. These difficulties are especially severe in the low-temperature, large-$U$ regime. In contrast, we find that new algorithms, when used alone or in combinations with themselves and the standard algorithm, can have significantly smaller autocorrelation times, in some cases being smaller by three orders of magnitude. The new algorithms, which use non-local moves, are discussed from the point of view of a general prescription for developing cluster algorithms. The loop-flip algorithm is also shown to be ergodic and to belong to the grand canonical ensemble. Extensions to other models and higher dimensions is briefly discussed., Comment: 36 pages, RevTex ver.3