BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

D. J. J. Farnell

D. J. J. Farnell contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elcond-mat.stat-mech

Trust snapshot

Quick read

Trust 21 - Emerging
21works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

21 published item(s)

preprint2016arXiv

Ground-State Ordering of the $J_1$--$J_2$ Model on the Simple Cubic and Body-Centered Cubic Lattices

The $J_1$--$J_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 $J_1$--$J_2$ model with antiferromagnetic nearest-neighbor bonds $J_1 >0$ and next-nearest-neighbor bonds $J_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=z_2J_2/z_1J_1$, where $z_1$ ($z_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.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.55$. Thus the strong frustration present in the $J_1$--$J_2$ Heisenberg model on the SC lattice may open a window for an unconventional quantum ground state in this three-dimensional spin model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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-$s$ Heisenberg antiferromagnet. We calculate the fundamental ground-state quantities, namely, the energy $e_0$, the sublattice magnetization $M_{\rm sub}$, the in-plane spin stiffness $ρ_s$ and the in-plane magnetic susceptibility $χ$ for spin quantum numbers $s=1/2, 1, \ldots, s_{\rm max}$, where $s_{\rm max}=9/2$ for $e_0$ and $M_{\rm sub}$, $s_{\rm max}=4$ for $ρ_s$ and $s_{\rm max}=3$ for $χ$. We use the data for $s \ge 3/2$ to estimate the leading quantum corrections to the classical values of $e_0$, $M_{\rm sub}$, $ρ_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 $s$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

We investigate ground states of $s$=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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ground-state phases of the frustrated spin-1/2 $J_{1}$--$J_{2}$--$J_{3}$ Heisenberg ferromagnet ($J_{1}<0$) on the honeycomb lattice with $J_{3}=J_{2}>0$

We study the ground-state (gs) properties of the frustrated spin-1/2 $J_{1}$--$J_{2}$--$J_{3}$ Heisenberg model on a honeycomb lattice with ferromagnetic (FM) nearest-neighbor ($J_{1}=-1$) exchange and frustrating antiferromagnetic (AFM) next-nearest-neighbor ($J_{2}>0$) and next-next-nearest-neighbor ($J_{3}>0$) exchanges, for the case $J_{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 (&#34;striped&#34;) collinear order at $J^{c}_{2} \approx 0.1095 \pm 0.0005$. We compare results for the FM case (with $J_{1}=-1$) to previous results for the corresponding AFM case (with $J_{1}=+1$). While the magnetic order parameters behave similarly for the FM and the AFM models for large values of the frustration parameter $J_{2}$, there are considerable differences between them for $J_{2}/|J_{1}| \lesssim 0.6$. For example, the quasiclassical collinear magnetic long-range order for the AFM model (with $J_{1}=+1$) breaks down at $J^{c_{2}}_{2} \approx 0.60$, whereas the &#34;equivalent&#34; point for the FM model (with $J_{1}=-1$) occurs at $J^{c}_{2} \approx 0.11$. Unlike in the AFM model (with $J_{1}=+1$), where a plaquette valence-bond crystal phase intrudes between the two corresponding quasiclassical AFM phases (with Néel and striped order) for $J^{c_{1}}_{2} < J_{2} < J^{c_{2}}_{2}$, with $J^{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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase diagram of a frustrated Heisenberg antiferromagnet on the honeycomb lattice: the $J_{1}$--$J_{2}$--$J_{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 $J_{1}$--$J_{2}$--$J_{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 $J_{1} \equiv 1$ we find four magnetically ordered phases in the parameter window $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 $J_{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) $(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) $(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) $(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.

0 citations0 reviewsNo code linkedOpen access
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 $J_{1}>0$ along two of the three lattice directions and $J_{2} \equiv κJ_{1} > 0$ along the third. In the classical limit the ground-state (GS) phase for $κ< 1/2$ has collinear ferrimagnetic (Néel$&#39;$) order where the $J_2$-coupled chain spins are ferromagnetically ordered in one direction with the remaining spins aligned in the opposite direction, while for $κ> 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<κ<κ_{c_1}$ for the Néel$&#39;$ state and for (at least part of) the region $κ>κ_{c_2}$ for the canted phase. However, they are now separated by a paramagnetic phase without either sort of magnetic order in the region $κ_{c_1} < κ< κ_{c_2}$, which includes the isotropic kagome point $κ= 1$ where the stable GS phase is now believed to be a topological ($\mathbb{Z}_2$) spin liquid. Our best numerical estimates are $κ_{c_1} = 0.515 \pm 0.015$ and $κ_{c_2} = 1.82 \pm 0.03$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The frustrated Heisenberg antiferromagnet on the checkerboard lattice: the $J_{1}$--$J_{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) $J_{1}$--$J_{2}$ system on the checkerboard lattice, with nearest-neighbor exchange bonds $J_{1}>0$ and next-nearest-neighbor bonds $J_{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 $κ< κ_{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)$ in a large-$s$ analysis (where $s$ 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$) 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 $κ= κ_{c_1}$, and the fourfold AFM states become infinitely susceptible to PVBC ordering at $κ= κ_{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 $κ= κ_{c_2}$. We thus find a Neel-ordered gs phase for $κ<κ_{c_1}$, a PVBC-ordered phase for $κ_{c_1} < κ< κ_{c_2}$, and a CDVBC-ordered phase for $κ> κ_{c_2}$. Both transitions are probably direct ones, although we cannot exclude very narrow coexistence regions confined to $0.79 \lesssim κ\lesssim 0.81$ and $1.20 \lesssim κ\lesssim 1.22$ respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The frustrated Heisenberg antiferromagnet on the honeycomb lattice: $J_{1}$--$J_{2}$ model

We study the ground-state (gs) phase diagram of the frustrated spin-1/2 $J_{1}$--$J_{2}$ antiferromagnet with $J_{2}=κJ_1>0$ ($J_{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 $κ_{c_1}<κ<κ_{c_2}$, where $κ_{c_1} \approx 0.207 \pm 0.003$ and $κ_{c_2} \approx 0.385 \pm 0.010$. The transition at $κ_{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 \lesssim κ\lesssim 0.24$), while that at $κ_{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$. The spiral phases that are present classically for all values $κ> 1/6$ are absent for all $κ\lesssim 1$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The magnetization process of the spin-one triangular-lattice Heisenberg antiferromagnet

We apply the coupled cluster method and exact diagonalzation to study the uniform susceptibility and the ground-state magnetization curve of the triangular-lattice spin-1 Heisenberg antiferromagnet. Comparing our theoretical data for the magnetization curve with recent measurements on the s=1 triangular lattice antiferromagnet Ba3NiSb2O9 we find a very good agreement.

0 citations0 reviewsNo code linkedOpen access
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 $J_{1}$-$J_{2}$-$J_{3}$ antiferromagnet with $J_{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 $κ< κ_{c_{1}} \approx 0.47$, a collinear striped AFM phase for $κ> κ_{c_{2}} \approx 0.60$, and a paramagnetic PVBC phase for $κ_{c_{1}} \lesssim κ\lesssim κ_{c_{2}}$. The transition at $κ_{c_{2}}$ appears to be of first-order type, while that at $κ_{c_{1}}$ is continuous. Since the Néel and PVBC phases break different symmetries our results favor the deconfinement scenario for the transition at $κ_{c_{1}}$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 $J$-$J&#39;$ model interpolating between both systems by varying $J&#39;/J$ from $J&#39;/J=0$ (bounce limit) to $J&#39;/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 &#34;pure&#34; bounce ($J&#39;/J=0$) and maple-leaf ($J&#39;/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 \le J&#39;/J \lesssim J&#39;_c/J $ and that the magnetic order parameter varies only weakly with $J&#39;/J$. At $J&#39;_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&#39;$ 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&#39;/J \gtrsim 1.07$ and another one at 2/3 of saturation emerging only at large $J&#39;/J \gtrsim 3$.

0 citations0 reviewsNo code linkedOpen access
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 \ge \frac 12$ for both the ground and excited states, and independently to &#34;generalised&#34; 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 \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 $J_1$--$J_2$ at the Majumdar-Ghosh point $J_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 &#34;special solutions&#34; 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 CaV$_4$O$_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 &#34;brute-force&#34; approach; i.e., one in which none of the translational or point-group symmetries of the finite lattice are used.

0 citations0 reviewsNo code linkedOpen access
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=0$), which interpolates between Heisenberg antiferromagnets (HAF&#39;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 $J_{1} \geq 0$ along two of the equivalent directions and $J_{2} \geq 0$ along the third. Sites connected by $J_{2}$ bonds are themselves connected to the missing NN non-kagome sites of the triangular lattice by bonds of strength $J_{1}&#39; \geq 0$. When $J_{1}&#39;=J_{1}$ and $J_{2}=0$ the model reduces to the square-lattice HAF. The magnetic ordering of the system is investigated and its $T=0$ phase diagram discussed. Results for the kagome HAF limit are among the best available.

0 citations0 reviewsNo code linkedOpen access
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 $J_{1} > 0$, but only half the next-nearest-neighbor bonds are present with a strength $J_{2} \equiv κJ_{1} > 0$. The bonds are arranged such that on the $2 \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 $κ_{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 $M$ approaches a nonzero value $M_{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κ$ also appears to be continuous, or very nearly so, at the critical point $κ_{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 $κ_{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.

0 citations0 reviewsNo code linkedOpen access
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|.

0 citations0 reviewsNo code linkedOpen access
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 `LSUB$m$&#39; 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$ for the square-lattice antiferromagnet, which is also in agreement with the results of other approximate methods (e.g., $χ=0.0669$ via QMC). Our estimate for the range of the extent of the ($M/M_s=$)$\frac 13$ magnetisation plateau for the triangular-lattice antiferromagnet is $1.37< λ< 2.15$, which is in good agreement with results of spin-wave theory ($1.248 < λ< 2.145$) and exact diagonalisations ($1.38 < λ< 2.16$). The CCM value for the in-plane magnetic susceptibility per site is $χ=0.065$, which is below the result of the spin-wave theory (evaluated to order 1/S) of $χ_{SWT}=0.0794$.

0 citations0 reviewsNo code linkedOpen access
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 $J_{1}$--$J_{2}&#39;$ model, defined on an anisotropic 2D lattice. With respect to an underlying square-lattice geometry the model contains antiferromagnetic ($J_{1} > 0$) bonds between nearest neighbors and competing ($J_{2}&#39;>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 $J_{2}&#39; \equiv κJ_{1}$ bonds along parallel chains and $J_{1}$ bonds providing an interchain coupling. Hence, the model interpolates between a spin-half HAF on the square lattice at one extreme ($κ= 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$. 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 $κ_{c_{1}} = 0.80 \pm 0.01$, by contrast with the corresponding second-order transition between the equivalent classical states at $κ_{\rm cl} = 0.5$. We also find strong evidence for a second critical point at $κ_{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 $J_{1}$--$J_{2}&#39;$ model that did not, however, evaluate the corresponding critical point.

0 citations0 reviewsNo code linkedOpen access
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, $s$. This new ``general-$s$&#39;&#39; 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$. 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 $Δ\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$.

0 citations0 reviewsNo code linkedOpen access
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 $z$-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 LSUB$m$ scheme in which we retain all $k$-body correlations defined on all possible locales of $m$ adjacent lattice sites ($k \le m$). The ``raw&#39;&#39; CCM LSUB$m$ 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 LSUB$m$ results to the limit $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.

0 citations0 reviewsNo code linkedOpen access