Source author record

D. J. J. Farnell

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

21works
2topics
4close collaborators

Published work

21 publications

preprint2016arXiv

Ground-State Ordering of the J1J_1--J2J_2 Model on the Simple Cubic and Body-Centered Cubic Lattices

The J1J_1--J2J_2 Heisenberg model is a "canonical" model in the field of quantum magnetism in order to study the interplay between frustration and quantum fluctuations as well as quantum phase transitions driven by frustration. Here we apply the Coupled Cluster Method (CCM) to study the spin-half J1J_1--J2J_2 model with antiferromagnetic nearest-neighbor bonds J1>0J_1 >0 and next-nearest-neighbor bonds J2>0J_2 >0 for the simple cubic (SC) and body-centered cubic (BCC) lattices. In particular, we wish to study the ground-state ordering of these systems as a function of the frustration parameter p=z2J2/z1J1p=z_2J_2/z_1J_1, where z1z_1 (z2z_2) is the number of nearest (next-nearest) neighbors. We wish to determine the positions of the phase transitions using the CCM and we aim to resolve the nature of the phase transition points. We consider the ground-state energy, order parameters, spin-spin correlation functions as well as the spin stiffness in order to determine the ground-state phase diagrams of these models. We find a direct first-order phase transition at a value of p=0.528p = 0.528 from a state of nearest-neighbor Néel order to next-nearest-neighbor Néel order for the BCC lattice. For the SC lattice the situation is more subtle. CCM results for the energy, the order parameter, the spin-spin correlation functions and the spin stiffness indicate that there is no direct first-order transition between ground-state phases with magnetic long-range order, rather it is more likely that two phases with antiferromagnetic long-range are separated by a narrow region of a spin-liquid like quantum phase around p=0.55p=0.55. Thus the strong frustration present in the J1J_1--J2J_2 Heisenberg model on the SC lattice may open a window for an unconventional quantum ground state in this three-dimensional spin model.

preprint2015arXiv

Ground-state properties of the triangular-lattice Heisenberg antiferromagnet with arbitrary spin quantum number ss

We apply the coupled cluster method to high orders of approximation and exact diagonalizations to study the ground-state properties of the triangular-lattice spin-ss Heisenberg antiferromagnet. We calculate the fundamental ground-state quantities, namely, the energy e0e_0, the sublattice magnetization MsubM_{\rm sub}, the in-plane spin stiffness ρsρ_s and the in-plane magnetic susceptibility χχ for spin quantum numbers s=1/2,1,…,smaxs=1/2, 1, \ldots, s_{\rm max}, where smax=9/2s_{\rm max}=9/2 for e0e_0 and MsubM_{\rm sub}, smax=4s_{\rm max}=4 for ρsρ_s and smax=3s_{\rm max}=3 for χχ. We use the data for s≥3/2s \ge 3/2 to estimate the leading quantum corrections to the classical values of e0e_0, MsubM_{\rm sub}, ρsρ_s, and χχ. In addition, we study the magnetization process, the width of the 1/3 plateau as well as the sublattice magnetizations in the plateau state as a function of the spin quantum number ss.

preprint2014arXiv

The spin-1/2 square-lattice J_1-J_2 model: The spin-gap issue

We use the coupled cluster method to high orders of approximation in order to calculate the ground-state energy, the ground-state magnetic order parameter, and the spin gap of the spin-1/2 J_1-J_2 model on the square lattice. We obtain values for the transition points to the magnetically disordered quantum paramagnetic phase of J_2^{c1}=0.454J_1 and J_2^{c2}= 0.588 J_1. The spin gap is zero in the entire parameter region accessible by our approach, i.e. for J_2 \le 0.49J_1 and J_2 > 0.58J_1. This finding is in favor of a gapless spin-liquid ground state in this parameter regime.

preprint2013arXiv

Quantum s = 1/2 Antiferromagnets on Archimedean Lattices: The Route from Semiclassical Magnetic Order to Nonmagnetic Quantum States

We investigate ground states of ss=1/2 Heisenberg antiferromagnets on the eleven two-dimensional (2D) Archimedian lattices by using the coupled cluster method. Magnetic interactions and quantum fluctuations play against each other subtly in 2D quantum magnets and the magnetic ordering is thus sensitive to the features of lattice topology. Archimedean lattices are those lattices that have 2D arrangements of regular polygons and they often build the underlying magnetic lattices of insulating quasi-two-dimensional quantum magnetic materials. Hence they allow a systematic study of the relationship between lattice topology and magnetic ordering. We find that the Archimedian lattices fall into three groups: those with semiclassical magnetic ground-state long-range order, those with a magnetically disordered (cooperative quantum paramagnetic) ground state, and those with a fragile magnetic order. The most relevant parameters affecting the magnetic ordering are the coordination number and the degree of frustration present.

preprint2012arXiv

Ground-state phases of the frustrated spin-1/2 J1J_{1}--J2J_{2}--J3J_{3} Heisenberg ferromagnet (J1<0J_{1}<0) on the honeycomb lattice with J3=J2>0J_{3}=J_{2}>0

We study the ground-state (gs) properties of the frustrated spin-1/2 J1J_{1}--J2J_{2}--J3J_{3} Heisenberg model on a honeycomb lattice with ferromagnetic (FM) nearest-neighbor (J1=−1J_{1}=-1) exchange and frustrating antiferromagnetic (AFM) next-nearest-neighbor (J2>0J_{2}>0) and next-next-nearest-neighbor (J3>0J_{3}>0) exchanges, for the case J3=J2J_{3}=J_{2}. We use the coupled cluster method in high orders of approximation, complemented by the exact diagonalization of a lattice with 32 sites, and calculate the gs energy, magnetic order parameter, and spin-spin correlation functions. We find a quantum phase transition between regions characterized by FM order and a form of AFM ("striped") collinear order at J2c≈0.1095±0.0005J^{c}_{2} \approx 0.1095 \pm 0.0005. We compare results for the FM case (with J1=−1J_{1}=-1) to previous results for the corresponding AFM case (with J1=+1J_{1}=+1). While the magnetic order parameters behave similarly for the FM and the AFM models for large values of the frustration parameter J2J_{2}, there are considerable differences between them for J2/∣J1∣≲0.6J_{2}/|J_{1}| \lesssim 0.6. For example, the quasiclassical collinear magnetic long-range order for the AFM model (with J1=+1J_{1}=+1) breaks down at J2c2≈0.60J^{c_{2}}_{2} \approx 0.60, whereas the "equivalent" point for the FM model (with J1=−1J_{1}=-1) occurs at J2c≈0.11J^{c}_{2} \approx 0.11. Unlike in the AFM model (with J1=+1J_{1}=+1), where a plaquette valence-bond crystal phase intrudes between the two corresponding quasiclassical AFM phases (with Néel and striped order) for J2c1<J2<J2c2J^{c_{1}}_{2} < J_{2} < J^{c_{2}}_{2}, with J2c1≈0.47J^{c_{1}}_{2} \approx 0.47, we find no clear indications in the FM model for an intermediate magnetically disordered phase between the phases exhibiting FM and striped order. Instead, the evidence points strongly to a direct first-order transition between the two ordered phases of the FM model.

preprint2012arXiv

Phase diagram of a frustrated Heisenberg antiferromagnet on the honeycomb lattice: the J1J_{1}--J2J_{2}--J3J_{3} model

We use the coupled cluster method in high orders of approximation to make a comprehensive study of the ground-state (GS) phase diagram of the spin-1/2 J1J_{1}--J2J_{2}--J3J_{3} model on a two-dimensional honeycomb lattice with antiferromagnetic (AFM) interactions up to third-nearest neighbors. Results are presented for the GS energy and the average local on-site magnetization. With the nearest-neighbor coupling strength J1≡1J_{1} \equiv 1 we find four magnetically ordered phases in the parameter window J2,J3∈[0,1]J_{2},J_{3} \in [0,1], namely the Néel (N), striped (S), and anti-Néel (aN) collinear AFM phases, plus a spiral phase. The aN phase appears as a stable GS phase in the classical version of the model only for values J3<0J_{3}<0. Each of these four ordered phases shares a boundary with a disordered quantum paramagnetic (QP) phase, and at several widely separated points on the phase boundaries the QP phase has an infinite susceptibility to plaquette valence-bond crystalline order. We identify all of the phase boundaries with good precision in the parameter window studied, and we find three tricritical quantum critical points therein at: (a) (J2c1,J3c1)=(0.51±0.01,0.69±0.01)(J_{2}^{c_1},J_{3}^{c_1})=(0.51 \pm 0.01,0.69 \pm 0.01) between the N, S, and QP phases; (b) (J2c2,J3c2)=(0.65±0.02,0.55±0.01)(J_{2}^{c_2},J_{3}^{c_2})=(0.65 \pm 0.02,0.55 \pm 0.01) between the S, spiral, and QP phases; and (c) (J2c3,J3c3)=(0.69±0.01,0.12±0.01)(J_{2}^{c_3},J_{3}^{c_3})=(0.69 \pm 0.01,0.12 \pm 0.01) between the spiral, aN, and QP phases.

preprint2012arXiv

Spin-1/2 Heisenberg antiferromagnet on an anisotropic kagome lattice

We use the coupled cluster method to study the zero-temperature properties of an extended two-dimensional Heisenberg antiferromagnet formed from spin-1/2 moments on an infinite spatially anisotropic kagome lattice of corner-sharing isosceles triangles, with nearest-neighbor bonds only. The bonds have exchange constants J1>0J_{1}>0 along two of the three lattice directions and J2≡κJ1>0J_{2} \equiv κJ_{1} > 0 along the third. In the classical limit the ground-state (GS) phase for κ<1/2κ< 1/2 has collinear ferrimagnetic (Néel′') order where the J2J_2-coupled chain spins are ferromagnetically ordered in one direction with the remaining spins aligned in the opposite direction, while for κ>1/2κ> 1/2 there exists an infinite GS family of canted ferrimagnetic spin states, which are energetically degenerate. For the spin-1/2 case we find that quantum analogs of both these classical states continue to exist as stable GS phases in some regions of the anisotropy parameter κκ, namely for 0<κ<κc10<κ<κ_{c_1} for the Néel′' state and for (at least part of) the region κ>κc2κ>κ_{c_2} for the canted phase. However, they are now separated by a paramagnetic phase without either sort of magnetic order in the region κc1<κ<κc2κ_{c_1} < κ< κ_{c_2}, which includes the isotropic kagome point κ=1κ= 1 where the stable GS phase is now believed to be a topological (Z2\mathbb{Z}_2) spin liquid. Our best numerical estimates are κc1=0.515±0.015κ_{c_1} = 0.515 \pm 0.015 and κc2=1.82±0.03κ_{c_2} = 1.82 \pm 0.03.

preprint2012arXiv

The frustrated Heisenberg antiferromagnet on the checkerboard lattice: the J1J_{1}--J2J_{2} model

We study the ground-state (gs) phases of the spin-half anisotropic planar pyrochlore (or crossed chain) model using the coupled cluster method (CCM). The model is a frustrated antiferromagnetic (AFM) J1J_{1}--J2J_{2} system on the checkerboard lattice, with nearest-neighbor exchange bonds J1>0J_{1}>0 and next-nearest-neighbor bonds J2≡κJ1>0J_{2} \equiv κJ_{1} > 0. Using various AFM classical ground states as CCM model states we present results for their gs energy, average on-site magnetization, and susceptibilities to plaquette valence-bond crystal (PVBC) and crossed-dimer valence-bond crystal (CDVBC) ordering. We show that the state with Neel ordering is the gs phase for κ<κc1≈0.80±0.01κ< κ_{c_1} \approx 0.80 \pm 0.01, but that none of the fourfold set of AFM states selected by quantum fluctuations at O(1/s)O(1/s) in a large-ss analysis (where ss is the spin quantum number) from the infinitely degenerate set of AFM states that form the gs phase for the classical version of the model (for κ>1κ>1) survives the quantum fluctuations to form a stable magnetically-ordered gs phase for the spin-half case. The Neel state becomes susceptible to PVBC ordering at or very near to κ=κc1κ= κ_{c_1}, and the fourfold AFM states become infinitely susceptible to PVBC ordering at κ=κc2≈1.22±0.02κ= κ_{c_2} \approx 1.22 \pm 0.02. In turn, we find that these states become infinitely susceptible to CDVBC ordering for all values of κκ above a certain critical value at or very near to κ=κc2κ= κ_{c_2}. We thus find a Neel-ordered gs phase for κ<κc1κ<κ_{c_1}, a PVBC-ordered phase for κc1<κ<κc2κ_{c_1} < κ< κ_{c_2}, and a CDVBC-ordered phase for κ>κc2κ> κ_{c_2}. Both transitions are probably direct ones, although we cannot exclude very narrow coexistence regions confined to 0.79≲κ≲0.810.79 \lesssim κ\lesssim 0.81 and 1.20≲κ≲1.221.20 \lesssim κ\lesssim 1.22 respectively.

preprint2012arXiv

The frustrated Heisenberg antiferromagnet on the honeycomb lattice: J1J_{1}--J2J_{2} model

We study the ground-state (gs) phase diagram of the frustrated spin-1/2 J1J_{1}--J2J_{2} antiferromagnet with J2=κJ1>0J_{2}=κJ_1>0 (J1>0J_{1}>0) on the honeycomb lattice, using the coupled-cluster method. We present results for the ground-state energy, magnetic order parameter and plaquette valence-bond crystal (PVBC) susceptibility. We find a paramagnetic PVBC phase for κc1<κ<κc2κ_{c_1}<κ<κ_{c_2}, where κc1≈0.207±0.003κ_{c_1} \approx 0.207 \pm 0.003 and κc2≈0.385±0.010κ_{c_2} \approx 0.385 \pm 0.010. The transition at κc1κ_{c_1} to the Néel phase seems to be a continuous deconfined transition (although we cannot exclude a very narrow intermediate phase in the range 0.21≲κ≲0.240.21 \lesssim κ\lesssim 0.24), while that at κc2κ_{c_2} is of first-order type to another quasiclassical antiferromagnetic phase that occurs in the classical version of the model only at the isolated and highly degenerate critical point κ=1/2κ= 1/2. The spiral phases that are present classically for all values κ>1/6κ> 1/6 are absent for all κ≲1κ\lesssim 1.

preprint2011arXiv

The frustrated Heisenberg antiferromagnet on the honeycomb lattice: A candidate for deconfined quantum criticality

We study the ground-state (gs) phase diagram of the frustrated spin-1/2 J1J_{1}-J2J_{2}-J3J_{3} antiferromagnet with J2=J3=κJ1J_{2} = J_{3} =κJ_1 on the honeycomb lattice, using coupled-cluster theory and exact diagonalization methods. We present results for the gs energy, magnetic order parameter, spin-spin correlation function, and plaquette valence-bond crystal (PVBC) susceptibility. We find a Néel antiferromagnetic (AFM) phase for κ<κc1≈0.47κ< κ_{c_{1}} \approx 0.47, a collinear striped AFM phase for κ>κc2≈0.60κ> κ_{c_{2}} \approx 0.60, and a paramagnetic PVBC phase for κc1≲κ≲κc2κ_{c_{1}} \lesssim κ\lesssim κ_{c_{2}}. The transition at κc2κ_{c_{2}} appears to be of first-order type, while that at κc1κ_{c_{1}} is continuous. Since the Néel and PVBC phases break different symmetries our results favor the deconfinement scenario for the transition at κc1κ_{c_{1}}.

preprint2011arXiv

The Heisenberg antiferromagnet on the kagome lattice with arbitrary spin: A high-order coupled cluster treatment

Starting with the sqrt{3} x sqrt{3} and the q=0 states as reference states we use the coupled cluster method to high orders of approximation to investigate the ground state of the Heisenberg antiferromagnet on the kagome lattice for spin quantum numbers s=1/2,1,3/2,2,5/2, and 3. Our data for the ground-state energy for s=1/2 are in good agreement with recent large-scale density-matrix renormalization group and exact diagonalization data. We find that the ground-state selection depends on the spin quantum number s. While for the extreme quantum case, s=1/2, the q=0 state is energetically favored by quantum fluctuations, for any s>1/2 the sqrt{3} x sqrt{3} state is selected. For both the sqrt{3} x sqrt{3} and the q=0 states the magnetic order is strongly suppressed by quantum fluctuations. Within our coupled cluster method we get vanishing values for the order parameter (sublattice magnetization) M for s=1/2 and s=1, but (small) nonzero values for M for s>1. Using the data for the ground-state energy and the order parameter for s=3/2,2,5/2, and 3 we also estimate the leading quantum corrections to the classical values.

preprint2011arXiv

The spin-half Heisenberg antiferromagnet on two Archimedian lattices: From the bounce lattice to the maple-leaf lattice and beyond

We investigate the ground state of the two-dimensional Heisenberg antiferromagnet on two Archimedean lattices, namely, the maple-leaf and bounce lattices as well as a generalized JJ-J′J' model interpolating between both systems by varying J′/JJ'/J from J′/J=0J'/J=0 (bounce limit) to J′/J=1J'/J=1 (maple-leaf limit) and beyond. We use the coupled cluster method to high orders of approximation and also exact diagonalization of finite-sized lattices to discuss the ground-state magnetic long-range order based on data for the ground-state energy, the magnetic order parameter, the spin-spin correlation functions as well as the pitch angle between neighboring spins. Our results indicate that the "pure" bounce (J′/J=0J'/J=0) and maple-leaf (J′/J=1J'/J=1) Heisenberg antiferromagnets are magnetically ordered, however, with a sublattice magnetization drastically reduced by frustration and quantum fluctuations. We found that magnetic long-range order is present in a wide parameter range 0≤J′/J≲Jc′/J0 \le J'/J \lesssim J'_c/J and that the magnetic order parameter varies only weakly with J′/JJ'/J. At Jc′≈1.45JJ'_c \approx 1.45 J a direct first-order transition to a quantum orthogonal-dimer singlet ground state without magnetic long-range order takes place. The orthogonal-dimer state is the exact ground state in this large-J′J' regime, and so our model has similarities to the Shastry-Sutherland model. Finally, we use the exact diagonalization to investigate the magnetization curve. We a find a 1/3 magnetization plateau for J′/J≳1.07J'/J \gtrsim 1.07 and another one at 2/3 of saturation emerging only at large J′/J≳3J'/J \gtrsim 3.

preprint2010arXiv

High-Order Coupled Cluster Method (CCM) Formalism 3: Finite-Size CCM

Recent developments of high-order CCM have been to extend existing formalism and codes to s≥12s \ge \frac 12 for both the ground and excited states, and independently to "generalised" expectation values for a wide range of one- and two-body spin operators. An advantage of the CCM is that the Goldstone linked-cluster theorem is obeyed at all levels of approximation and so it provides results in the infinite lattice limit N→∞N \to \infty from the outset. However, recent results have also shown that the CCM can provide exact (symmetry-breaking) results for the spin-half linear-chain J1J_1--J2J_2 at the Majumdar-Ghosh point J2/J1=0.5J_2/J_1=0.5 by identifying special solutions of the CCM equations for the usual Néel model state. Interestingly, the CCM provides exact (non-symmetry-breaking) results for systems in which small magnetic clusters become de-coupled from each other when the bonds connecting them tend to zero. These exact results involve the identification of "special solutions" of the CCM equations for the Néel state. An example of this is given by a spin-half system with nearest-neighbour bonds for an underlying lattice corresponding to the magnetic material CaV4_4O9_9 (CAVO) in which one of the two different types of bonds on the lattice tend to zero. Larger finite-sized systems may be considered by appropriate choice of the unit cell and the bonds on it. We show here that exact diagonalisation results for ground-state energy and excitation energy gap for the spin-half and spin-one linear Heisenberg model on chains of length up to N=12 sites for s=1/2 and N=6 sites for s=1 with periodic boundary conditions are reproduced exactly using high-order CCM via this "brute-force" approach; i.e., one in which none of the translational or point-group symmetries of the finite lattice are used.

preprint2010arXiv

Magnetic order in a spin-1/2 interpolating kagome-square Heisenberg antiferromagnet

The coupled cluster method is applied to a spin-half model at zero temperature (T=0T=0), which interpolates between Heisenberg antiferromagnets (HAF's) on a kagome and a square lattice. With respect to an underlying triangular lattice the strengths of the Heisenberg bonds joining the nearest-neighbor (NN) kagome sites are J1≥0J_{1} \geq 0 along two of the equivalent directions and J2≥0J_{2} \geq 0 along the third. Sites connected by J2J_{2} bonds are themselves connected to the missing NN non-kagome sites of the triangular lattice by bonds of strength J1′≥0J_{1}' \geq 0. When J1′=J1J_{1}'=J_{1} and J2=0J_{2}=0 the model reduces to the square-lattice HAF. The magnetic ordering of the system is investigated and its T=0T=0 phase diagram discussed. Results for the kagome HAF limit are among the best available.

preprint2010arXiv

Magnetic order on a frustrated spin-1/2 Heisenberg antiferromagnet on the Union Jack lattice

We use the coupled cluster method (CCM) to study the zero-temperature phase diagram of a 2D frustrated spin-half antiferromagnet, the so-called Union Jack model. It is defined on a square lattice such that all nearest-neighbor bonds are present with a strength J1>0J_{1} > 0, but only half the next-nearest-neighbor bonds are present with a strength J2≡κJ1>0J_{2} \equiv κJ_{1} > 0. The bonds are arranged such that on the 2×22 \times 2 unit cell they form the pattern of the Union Jack flag. We find strong evidence for a first phase transition between a Néel phase and a canted ferrimagnetic phase at a critical coupling κc1=0.66±0.02κ_{c_{1}} = 0.66 \pm 0.02. At the transition the energy and its first derivative seem continuous, thus providing a typical scenario of a second-order transition, although a weakly first-order transition cannot be excluded. By contrast, the average on-site magnetization MM approaches a nonzero value Mc1=0.195±0.005M_{c_{1}}=0.195 \pm 0.005 on both sides of the transition, which is more typical of a first-order transition. The slope dM/dκdM/dκ also appears to be continuous, or very nearly so, at the critical point κc1κ_{c_{1}}. We find strong evidence that the canted phase becomes unstable at large values of κκ, and hence we have also used the CCM with a model collinear semi-stripe-ordered ferrimagnetic state in which alternating rows (and columns) are ferromagnetically and antiferromagnetically ordered. We find tentative evidence, based on the relative energies of the two states, for a second (first-order) phase transition between the canted and semi-stripe-ordered states at a large value of the coupling parameter around κc2≈125±5κ_{c_{2}} \approx 125 \pm 5. This prediction, however, is based on an extrapolation of the CCM results for the canted state into regimes where the CCM equations at any level of approximation beyond the lowest have no solutions.

preprint2010arXiv

The frustrated spin-1/2 J1-J2 Heisenberg ferromagnet on the square lattice: Exact diagonalization and Coupled-Cluster study

We investigate the ground-state magnetic order of the spin-1/2 J1-J2 Heisenberg model on the square lattice with ferromagnetic nearest-neighbor exchange J1<0 and frustrating antiferromagnetic next-nearest neighbor exchange J2>0. We use the coupled-cluster method to high orders of approximation and Lanczos exact diagonalization of finite lattices of up to N=40 sites in order to calculate the ground-state energy, the spin-spin correlation functions, and the magnetic order parameter. We find that the transition point at which the ferromagnetic ground state disappears is given by J2^{c1}=0.393|J1| (exact diagonalization) and J2^{c1}=0.394|J1| (coupled-cluster method). We compare our results for ferromagnetic J1 with established results for the spin-1/2 J1-J2 Heisenberg model with antiferromagnetic J1. We find that both models (i.e., ferro- and antiferromagnetic J1) behave similarly for large J2, although significant differences between them are observed for J2/|J1| \lesssim 0.6. Although the semiclassical collinear magnetic long-range order breaks down at J2^{c2} \approx 0.6J1 for antiferromagnetic J1, we do not find a similar breakdown of this kind of long-range order until J2 \sim 0.4|J1| for the model with ferromagnetic J1. Unlike the case for antiferromagnetic J1, if an intermediate disordered phase does occur between the phases exhibiting semiclassical collinear stripe order and ferromagnetic order for ferromagnetic J1 then it is likely to be over a very small range below J2 \sim 0.4|J1|.

preprint2009arXiv

High-Order Coupled Cluster Method Study of Frustrated and Unfrustrated Quantum Magnets in External Magnetic Fields

We apply the coupled cluster method (CCM) in order to study the ground-state properties of the (unfrustrated) square-lattice and (frustrated) triangular-lattice spin-half Heisenberg antiferromagnets in the presence of external magnetic fields. Here we determine and solve the basic CCM equations by using the localised approximation scheme commonly referred to as the `LSUBmm' approximation scheme and we carry out high-order calculations by using intensive computational methods. We calculate the ground-state energy, the uniform susceptibility, the total (lattice) magnetisation and the local (sublattice) magnetisations as a function of the magnetic field strength. Our results for the lattice magnetisation of the square-lattice case compare well to those results of QMC for all values of the applied external magnetic field. We find a value for magnetic susceptibility of χ=0.070χ=0.070 for the square-lattice antiferromagnet, which is also in agreement with the results of other approximate methods (e.g., χ=0.0669χ=0.0669 via QMC). Our estimate for the range of the extent of the (M/Ms=M/M_s=)13\frac 13 magnetisation plateau for the triangular-lattice antiferromagnet is 1.37<λ<2.151.37< λ< 2.15, which is in good agreement with results of spin-wave theory (1.248<λ<2.1451.248 < λ< 2.145) and exact diagonalisations (1.38<λ<2.161.38 < λ< 2.16). The CCM value for the in-plane magnetic susceptibility per site is χ=0.065χ=0.065, which is below the result of the spin-wave theory (evaluated to order 1/S) of χSWT=0.0794χ_{SWT}=0.0794.

preprint2009arXiv

Magnetic order in a spin-half interpolating square-triangle Heisenberg antiferromagnet

Using the coupled cluster method we study the zero-temperature phase diagram of a spin-half Heisenberg antiferromagnet (HAF), the so-called J1J_{1}--J2′J_{2}' model, defined on an anisotropic 2D lattice. With respect to an underlying square-lattice geometry the model contains antiferromagnetic (J1>0J_{1} > 0) bonds between nearest neighbors and competing (J2′>0J_{2}'>0) bonds between next-nearest neighbors across only one of the diagonals of each square plaquette, the same diagonal in every square. Considered on an equivalent triangular-lattice geometry the model may be regarded as having two sorts of nearest-neighbor bonds, with J2′≡κJ1J_{2}' \equiv κJ_{1} bonds along parallel chains and J1J_{1} bonds providing an interchain coupling. Hence, the model interpolates between a spin-half HAF on the square lattice at one extreme (κ=0κ= 0) and a set of decoupled spin-half chains at the other (κ→∞κ\to \infty), with the spin-half HAF on the triangular lattice in between at κ=1κ= 1. We find strong evidence that quantum fluctuations favor a first-order transition from quasiclassical Néel order to a quantum helical state at a first critical point at κc1=0.80±0.01κ_{c_{1}} = 0.80 \pm 0.01, by contrast with the corresponding second-order transition between the equivalent classical states at κcl=0.5κ_{\rm cl} = 0.5. We also find strong evidence for a second critical point at κc2=1.8±0.4κ_{c_{2}} = 1.8 \pm 0.4 where another first-order transition occurs, this time from the quantum helical phase to a collinear stripe-ordered phase. This latter result provides quantitative verification of a recent qualitative prediction of Starykh and Balents [Phys.\ Rev. Lett. {\bf 98}, 077205 (2007)] for the J1J_{1}--J2′J_{2}' model that did not, however, evaluate the corresponding critical point.

preprint2003arXiv

Coupled Cluster Method Calculations Of Quantum Magnets With Spins Of General Spin Quantum Number

We present a new high-order coupled cluster method (CCM) formalism for the ground states of lattice quantum spin systems for general spin quantum number, ss. This new ``general-ss'' formalism is found to be highly suitable for a computational implementation, and the technical details of this implementation are given. To illustrate our new formalism we perform high-order CCM calculations for the one-dimensional spin-half and spin-one antiferromagnetic {\it XXZ} models and for the one-dimensional spin-half/spin-one ferrimagnetic {\it XXZ} model. The results for the ground-state properties of the isotropic points of these systems are seen to be in excellent quantitative agreement with exact results for the special case of the spin-half antiferromagnet and results of density matrix renormalisation group (DMRG) calculations for the other systems. Extrapolated CCM results for the sublattice magnetisation of the spin-half antiferromagnet closely follow the exact Bethe Ansatz solution, which contains an infinite-order phase transition at Δ=1Δ=1. By contrast, extrapolated CCM results for the sublattice magnetisation of the spin-one antiferromagnet using this same scheme are seen to go to zero at Δ≈1.2Δ\approx 1.2, which is in excellent agreement with the value for the onset of the Haldane phase for this model. Results for sublattice magnetisations of the ferrimagnet for both the spin-half and spin-one spins are non-zero and finite across a wide range of ΔΔ, up to and including the Heisenberg point at Δ=1Δ=1.

preprint2000arXiv

High-Order Coupled Cluster Method Calculations for the Ground- and Excited-State Properties of the Spin-Half XXZ Model

In this article, we present new results of high-order coupled cluster method (CCM) calculations, based on a Néel model state with spins aligned in the zz-direction, for both the ground- and excited-state properties of the spin-half {\it XXZ} model on the linear chain, the square lattice, and the simple cubic lattice. In particular, the high-order CCM formalism is extended to treat the excited states of lattice quantum spin systems for the first time. Completely new results for the excitation energy gap of the spin-half {\it XXZ} model for these lattices are thus determined. These high-order calculations are based on a localised approximation scheme called the LSUBmm scheme in which we retain all kk-body correlations defined on all possible locales of mm adjacent lattice sites (k≤mk \le m). The ``raw'' CCM LSUBmm results are seen to provide very good results for the ground-state energy, sublattice magnetisation, and the value of the lowest-lying excitation energy for each of these systems. However, in order to obtain even better results, two types of extrapolation scheme of the LSUBmm results to the limit m→∞m \to \infty (i.e., the exact solution in the thermodynamic limit) are presented. The extrapolated results provide extremely accurate results for the ground- and excited-state properties of these systems across a wide range of values of the anisotropy parameter.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.