We study the frustrated ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model with ferromagnetic nearest-neighbor coupling ${J}_{1}<0$ and antiferromagnetic next-nearest-neighbor coupling ${J}_{2}>0$ at and close to the $z=4$ quantum critical point (QCP) at ${J}_{1}/{J}_{2}=\ensuremath{-}4$. The ${J}_{1}\ensuremath{-}{J}_{2}$ model plays an important role for recently synthesized chain cuprates as well as in supersymmetric Yang-Mills theories. We study the thermodynamic properties using field theory, a modified spin-wave theory, as well as numerical density-matrix renormalization group calculations. Furthermore, we compare with results for the classical model obtained by analytical methods and Monte Carlo simulations. As one of our main results, we present numerical evidence that the susceptibility at the QCP seems to diverge with temperature $T$ as $\ensuremath{\chi}\ensuremath{\sim}{T}^{\ensuremath{-}1.2}$ in the quantum case, in contrast to the classical model where $\ensuremath{\chi}\ensuremath{\sim}{T}^{\ensuremath{-}4/3}$.
We study the frustrated J1\ensuremath−J2 Heisenberg model with ferromagnetic nearest-neighbor coupling J1<0 and antiferromagnetic next-nearest-neighbor coupling J2>0 at and close to the z=4 quantum critical point (QCP) at J1/J2=\ensuremath−4. The J1\ensuremath−J2 model plays an important role for recently synthesized chain cuprates as well as in supersymmetric Yang-Mills theories. We study the thermodynamic properties using field theory, a modified spin-wave theory, as well as numerical density-matrix renormalization group calculations. Furthermore, we compare with results for the classical model obtained by analytical methods and Monte Carlo simulations. As one of our main results, we present numerical evidence that the susceptibility at the QCP seems to diverge with temperature T as \ensuremathχ\ensuremath∼T\ensuremath−1.2 in the quantum case, in contrast to the classical model where \ensuremathχ\ensuremath∼T\ensuremath−4/3.
Abstract. We studied the thermodynamics of the one-dimensional J1-J2 spin-1/2 Heisenberg chain for ferromagnetic nearest-neighbor bonds J1 < 0 and frustrating antiferromagnetic nextnearest-neighbor bonds J2> 0 using full diagonalization of finite rings and a second-order Green-function formalism. Thereby we focus on J2 < |J1|/4 where the ground state is still ferromagnetic, but the frustration influences the thermodynamic properties. We found that their critical indices are not changed by J2. The analysis of the low-temperature behavior of the susceptibility χ leads to the conclusion that this behavior changes from χ ∝ T −2 at J2 < |J1|/4 to χ ∝ T −3/2 at the quantum-critical point J2 = |J1|/4. Another effect of the frustration is the appearance of an extra low-T maximum in the specific heat Cv(T) for J2 � |J1|/8, indicating its strong influence on the low-energy spectrum. Introduction: In low-dimensional frustrated quantum magnets thermal and quantum fluctuations strongly influence the low-temperature physics [1,2]. Special attention has been paid to one-dimensional (1D) J1-J2 quantum Heisenberg magnets, see Ref. [3] and references therein. Recent experimental studies have shown that edge-shared chain cuprates, such as LiVCuO4, Li(Na)Cu2O2, Li2ZrCuO4, and Li2CuO2 [4–13], represent a family of quantum magnets for
We study the ground-state (GS) properties of the frustrated spin-21J1--J2--J3 Heisenberg model on the two-dimensional honeycomb lattice with ferromagnetic nearest-neighbor (${J}_{1}=\ensuremath{-}1$) exchange and frustrating antiferromagnetic next-nearest-neighbor (J2>0) and next-next-nearest-neighbor (J3>0) exchanges, for the case J3=J2. We use the coupled-cluster method implemented to high orders of approximation, complemented by the Lanczos exact diagonalization of a large finite lattice with 32 sites, in order to calculate the GS energy, magnetic order parameter, and spin-spin correlation functions. In one scenario we find a quantum phase transition point between regions characterized by ferromagnetic order and a form of antiferromagnetic (``striped'') collinear order at ${J}_{2}^{c}\ensuremath{\approx}0.1095\ifmmode\pm\else\textpm\fi{}0.0005$, which is below the corresponding hypothetical transition point at J2cl=71 ($\ensuremath{\approx}0.143$) for the classical version of the model, in which we momentarily ignore the intervening noncollinear spiral phase in the region 101<J2<51. Hence we see that quantum fluctuations appear to stabilize somewhat the collinear antiferromagnetic order in preference to the ferromagnetic order in this model. We compare results for the present ferromagnetic case (with ${J}_{1}=\ensuremath{-}1$) to previous results for the corresponding antiferromagnetic case (with J1=+1). The magnetic order parameter is found to behave similarly for the ferromagnetic and the antiferromagnetic models for large values of the frustration parameter J2. However, there are considerable differences in the behavior of the order parameters for the two models for ${J}_{2}/|{J}_{1}|\ensuremath{\lesssim}0.6$. For example, the quasiclassical collinear magnetic long-range order for the antiferromagnetic model (with J1=+1) breaks down at ${J}_{2}^{{c}_{2}}\ensuremath{\approx}0.60$, whereas the ``equivalent'' point for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) occurs at ${J}_{2}^{c}\ensuremath{\approx}0.11$. Unlike in the antiferromagnetic model (with J1=+1), where a plaquette valence-bond crystal phase intrudes between the two corresponding quasiclassical antiferromagnetic phases (with N\'eel and striped order) for J2c1<J2<J2c2, with ${J}_{2}^{{c}_{1}}\ensuremath{\approx}0.47$, we find no clear indications at all in the ferromagnetic model for an intermediate magnetically disordered phase between the corresponding phases exhibiting ferromagnetic and striped order. Instead the evidence for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) points to one of two scenarios: either there is a direct first-order transition between the two magnetically ordered phases, as mentioned above; or there exists an intervening phase between them in the very narrow range $0.10\ensuremath{\lesssim}{J}_{2}\ensuremath{\lesssim}0.12$, which is probably a remnant of the spiral phase that exists in the classical counterpart of the model over the larger range 101<J2<51.
We study Heisenberg spin-1/2 and spin-1 chains with alternating ferromagnetic (J1F) and antiferromagnetic (J1A) nearest-neighbor interactions and ferromagnetic next-nearest-neighbor interaction (J2F). In this model frustration is induced by the non-zero J2F and the effect of dimerization is realized by changing the strength of J1A for a fixed value of J1F. The model with site spin s behaves like a Haldane spin chain with site spin 2s in the limit of vanishing J2F and large J1F/J1A. We show that the exact ground state of the model can be found along a line in the parameter space. The phase diagram in the space of J1A−J2F is determined using numerical techniques complemented by some analytical calculations. A number of quantities, including the ground state spin, energy gap, structure factor, entanglement entropy, and zero temperature magnetization, are studied to understand the complete phase diagram. An interesting and application-wise important feature of this model is that it exhibits a macroscopic magnetization jump in the presence of a magnetic field; we study this using an effective Hamiltonian.
The thermodynamic quantities (spin-spin correlation functions $\ensuremath{\langle}{\mathbf{S}}_{0}{\mathbf{S}}_{n}\ensuremath{\rangle}$, correlation length $\ensuremath{\xi}$, spin susceptibility $\ensuremath{\chi}$, and specific heat CV) of the frustrated one-dimensional J1-J2 Heisenberg ferromagnet with arbitrary spin quantum number S below the quantum critical point, i.e., for J2<∣J1∣/4, are calculated using a rotation-invariant Green-function formalism and full diagonalization as well as a finite-temperature Lanczos technique for finite chains of up to N=20 sites. The low-temperature behavior of the susceptibility $\ensuremath{\chi}$ and the correlation length $\ensuremath{\xi}$ are well described by $\ensuremath{\chi}=\frac{2}{3}{S}^{4}\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right){T}^{\ensuremath{-}2}+A{S}^{5/2}{\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right)}^{1/2}{T}^{\ensuremath{-}3/2}$ and $\ensuremath{\xi}={S}^{2}\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right){T}^{\ensuremath{-}1}+B{S}^{1/2}{\left(|{J}_{1}|\ensuremath{-}4{J}_{2}\right)}^{1/2}{T}^{\ensuremath{-}1/2}$ with $A\ensuremath{\approx}1.1\ensuremath{\cdots}1.2$ and $B\ensuremath{\approx}0.84\ensuremath{\cdots}0.89$. The vanishing of the factors in front of the temperature at J2=∣J1∣/4 indicates a change of the critical behavior of $\ensuremath{\chi}$ and $\ensuremath{\xi}$ at $T\ensuremath{\rightarrow}0$. The specific heat may exhibit an additional frustration-induced low-temperature maximum when approaching the quantum critical point. Such a maximum appears for S=1/2 and S=1, but was not found for S>1.
The thermodynamics of coupled frustrated ferromagnetic chains is studied within a spin-rotation-invariant Green's-function approach. We consider an isotropic Heisenberg spin-half system with a ferromagnetic in-chain coupling J1<0 between nearest neighbors and a frustrating antiferromagnetic next-nearest-neighbor in-chain coupling J2>0. We focus on the moderate strength of frustration J2<∣J1∣/4 such that the in-chain spin-spin correlations are predominantly ferromagnetic. We consider two interchain couplings (ICs) ${J}_{\ensuremath{\perp},y}$ and ${J}_{\ensuremath{\perp},z}$, corresponding to the two axes perpendicular to the chain, where ferromagnetic as well as antiferromagnetic ICs are taken into account. We discuss the influence of frustration on the ground-state properties for antiferromagnetic ICs, where the ground state is of a quantum nature. The major part of our study is devoted to the finite-temperature properties. We calculate the critical temperature Tc as a function of the competing exchange couplings ${J}_{2},{J}_{\ensuremath{\perp},y},{J}_{\ensuremath{\perp},z}$. We find that for fixed ICs, Tc decreases monotonically with increasing frustration J2, where as ${J}_{2}\ensuremath{\rightarrow}|{J}_{1}|/4$ the Tc(J2) curve drops down rapidly. To characterize the magnetic ordering below and above Tc, we calculate the spin-spin correlation functions $\ensuremath{\langle}{\mathbf{S}}_{0}{\mathbf{S}}_{\mathbf{R}}\ensuremath{\rangle}$, the magnetic order parameter M, the uniform static susceptibility ${\ensuremath{\chi}}_{0}$, as well as the correlation length $\ensuremath{\xi}$. Moreover, we discuss the specific heat CV and the temperature dependence of the excitation spectrum ${\ensuremath{\omega}}_{\mathbf{q}}$. As ${J}_{2}\ensuremath{\rightarrow}|{J}_{1}|/4$, some unusual frustration-induced features were found, such as an increase of the in-chain spin stiffness (in the case of ferromagnetic ICs) or of the in-chain spin-wave velocity (in the case of antiferromagnetic ICs) with growing temperature.
The antiferromagnetic ${J}_{1}\ensuremath{-}{J}_{2}$ model is a spin-1/2 chain with isotropic exchange J1>0 between first neighbors and ${J}_{2}=\ensuremath{\alpha}{J}_{1}$ between second neighbors. The model supports both gapless quantum phases with nondegenerate ground states and gapped phases with $\mathrm{\ensuremath{\Delta}}(\ensuremath{\alpha})>0$ and doubly degenerate ground states. Exact thermodynamics is limited to $\ensuremath{\alpha}=0$, the linear Heisenberg antiferromagnet (HAF). Exact diagonalization of small systems at frustration $\ensuremath{\alpha}$ followed by density matrix renormalization group calculations returns the entropy density $S(T,\ensuremath{\alpha},N)$ and magnetic susceptibility $\ensuremath{\chi}(T,\ensuremath{\alpha},N)$ of progressively larger systems up to N=96 or 152 spins. Convergence to the thermodynamic limit, $S(T,\ensuremath{\alpha})$ or $\ensuremath{\chi}(T,\ensuremath{\alpha})$, is demonstrated down to $T/J\ensuremath{\sim}0.01$ in the sectors $\ensuremath{\alpha}<1$ and $\ensuremath{\alpha}>1$. $S(T,\ensuremath{\alpha})$ yields the critical points between gapless phases with ${S}^{\ensuremath{'}}(0,\ensuremath{\alpha})>0$ and gapped phases with ${S}^{\ensuremath{'}}(0,\ensuremath{\alpha})=0$. The ${S}^{\ensuremath{'}}(T,\ensuremath{\alpha})$ maximum at ${T}^{*}(\ensuremath{\alpha})$ is obtained directly in chains with large $\mathrm{\ensuremath{\Delta}}(\ensuremath{\alpha})$ and by extrapolation for small gaps. A phenomenological approximation for $S(T,\ensuremath{\alpha})$ down to T=0 indicates power-law deviations ${T}^{\ensuremath{-}\ensuremath{\gamma}(\ensuremath{\alpha})}$ from $exp[\ensuremath{-}\mathrm{\ensuremath{\Delta}}(\ensuremath{\alpha})/T]$ with exponent $\ensuremath{\gamma}(\ensuremath{\alpha})$ that increases with $\ensuremath{\alpha}$. The $\ensuremath{\chi}(T,\ensuremath{\alpha})$ analysis also yields power-law deviations, but with exponent $\ensuremath{\eta}(\ensuremath{\alpha})$ that decreases with $\ensuremath{\alpha}$. Spin correlation functions account for $S(T,\ensuremath{\alpha})$ differences between frustration $\ensuremath{\alpha}<1$ within a chain and $\ensuremath{\alpha}>1$ between HAFs on sublattices. $S(T,\ensuremath{\alpha})$ and the spin density $\ensuremath{\rho}(T,\ensuremath{\alpha})=4T\ensuremath{\chi}(T,\ensuremath{\alpha})$ probe the thermal and magnetic fluctuations, respectively, of strongly correlated spin states. Gapless chains have constant $S(T,\ensuremath{\alpha})/\ensuremath{\rho}(T,\ensuremath{\alpha})$ for T<0.10. Remarkably, the ratio decreases (increases) with T in chains with large (small) $\mathrm{\ensuremath{\Delta}}(\ensuremath{\alpha})$.