In this article we present new formalism for high-order coupled cluster method (CCM) calculations for "generalized" ground-state expectation values and the excited states of quantum magnetic systems with spin quantum number s ≥ 1/2. We use high-order CCM to demonstrate spontaneous symmetry breaking in the spin-half J1-J2 model for the linear chain using the coupled cluster method (CCM). We show that we are able to reproduce exactly the dimerized ground (ket) state at the Majumdar-Ghosh point (J2/J1= 1/2) using a Néel model state. We show that the onset of dimerized phase is indicated by a bifurcation of the nearest-neighbour ket- and bra-state correlation coefficients for the nearest-neighbour Néel model state. We show that ground-state energies are in good agreement with the results of exact diagonalizations of finite-length chains across this entire regime (i. e., J1>0 and J2≤1/2). The effects of the bifurcation point are also observed for the sublattice magnetization for the nearest-neighbour model state. Finally, we use the new formalism for the excited state in order to obtain the excitation energy as a function of J2/J1 for the nearest-neighbour model state by solving up to the LSUB14 level of approximation. We obtain an extrapolated value for the excited-state energy gap of -0.0036 at J2/J1=0.0 and of 0.2310 at J2/J1=0.5. We show that an excitation energy gap opens up at J2/J1 ≈ 0.24, although the gap only becomes large at J2/J1≈ 0.4.
In this article we present new formalism for high-order coupled cluster method (CCM) calculations for "generalized" ground-state expectation values and the excited states of quantum magnetic systems with spin quantum number s ≥ 1/2. We use high-order CCM to demonstrate spontaneous symmetry breaking in the spin-half J1-J2 model for the linear chain using the coupled cluster method (CCM). We show that we are able to reproduce exactly the dimerized ground (ket) state at the Majumdar-Ghosh point (J2/J1= 1/2) using a Néel model state. We show that the onset of dimerized phase is indicated by a bifurcation of the nearest-neighbour ket- and bra-state correlation coefficients for the nearest-neighbour Néel model state. We show that ground-state energies are in good agreement with the results of exact diagonalizations of finite-length chains across this entire regime (i. e., J1>0 and J2≤1/2). The effects of the bifurcation point are also observed for the sublattice magnetization for the nearest-neighbour model state. Finally, we use the new formalism for the excited state in order to obtain the excitation energy as a function of J2/J1 for the nearest-neighbour model state by solving up to the LSUB14 level of approximation. We obtain an extrapolated value for the excited-state energy gap of -0.0036 at J2/J1=0.0 and of 0.2310 at J2/J1=0.5. We show that an excitation energy gap opens up at J2/J1 ≈ 0.24, although the gap only becomes large at J2/J1≈ 0.4.
Recent developments of high-order CCM have been to extend existing formalism and codes to s≥21 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→∞ from the outset. However, recent results have also shown that the CCM can provide exact (symmetry-breaking) results for the spin-half linear-chain J1--J2 at the Majumdar-Ghosh point J2/J1=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 CaV4O9 (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.
The coupled cluster method (CCM) is a well-known method of quantum many-body theory, and in this article we present an application of the CCM to the spin-half ${J}_{1}{\ensuremath{-}J}_{2}$ quantum spin model with nearest- and next-nearest-neighbor interactions on the linear chain and the square lattice. We present results for ground-state expectation values of such quantities as the energy and the sublattice magnetization. The presence of critical points in the solution of the CCM equations, which are associated with phase transitions in the real system, is investigated. Completely distinct from the investigation of the critical points, we also make a link between the expansion coefficients of the ground-state wave function in terms of an Ising basis and the CCM ket-state correlation coefficients. We are thus able to present evidence of the breakdown, at a given value of J2/J1, of the Marshall-Peierls sign rule which is known to be satisfied at the pure Heisenberg point (J2=0) on any bipartite lattice. For the square lattice, our best estimates of the points at which the sign rule breaks down and at which the phase transition from the antiferromagnetic phase to the frustrated phase occurs are, respectively, given by ${J}_{2}{/J}_{1}\ensuremath{\approx}0.26$ and ${J}_{2}{/J}_{1}\ensuremath{\approx}0.61.$
CCM Calculations For The Ground-State Properties Of The One-Dimensional Spin-Half J1--J2 Model: Possible Evidence of Collinear Ordering for J2/J1>21?
In this article we investigate the linear-chain spin-half J1--J2 model by using high-order coupled cluster method (CCM) calculations. We employ three model states, namely, a nearest-neighbour (n.n.) Néel model state in which neighbouring spins are antiparallel, a next-nearest-neighbour (n.n.n.) Néel model state in which next-neighbouring spins are antiparallel, and finally a type of "double spiral" model state with {\it two} pitch angles. For J2/J1≤21, we find that the n.n. Néel model state produces the lowest energies. For J2/J1>21, we find that the stable states for the quantum system are those for the "traditional" spiral state in which the two pitch angles are identical and the collinear n.n.n. Néel model state. We show that ground-state energies for the collinear n.n.n. model state are lower than those of the spiral state in a finite region for J2/J1>21. We produce results for the dimer order parameter by using the n.n.n. Néel model state and these results are in reasonable agreement with previous results of DMRG. Again for the n.n.n. Néel model state, it is shown that the correlation length increases with increasing J2/J1 above J2/J1=21 using this model state and that shift in the peak of the structure function occurs from q=π to q=π/2 as one increases J2/J1. Although these results are not conclusive, they are intriguing because they suggest that a ground state exhibiting collinear order might occur for some region for J2/J1≥21.
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 LSUBm scheme in which we retain all k-body correlations defined on all possible locales of m adjacent lattice sites (k≤m). The ``raw'' CCM LSUBm 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 LSUBm results to the limit m→∞ (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.
The ${J}_{1}\text{\ensuremath{-}}{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}\text{\ensuremath{-}}{J}_{2}$ model with antiferromagnetic nearest-neighbor bonds J1>0 and next-nearest-neighbor bonds J2>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/z1J1, where z1(z2) 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\'eel order to next-nearest-neighbor N\'eel 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}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the sc lattice may open a window for an unconventional quantum ground state in this three-dimensional spin model.
High-Order Coupled Cluster Method (CCM) Formalism 1: Ground- and Excited-State Properties of Lattice Quantum Spin Systems with s>=1/2
The coupled cluster method (CCM) is a powerful and widely applied technique of modern-day quantum many-body theory. It has been used with great success in order to understand the properties of quantum magnets at zero temperature. This is due largely to the application of computational techniques that allow the method to be applied to high orders of approximation using localised approximation schemes, e.g., such as the LSUBm scheme. In this article, the high-order CCM formalism for the ground and excited states of quantum magnetic systems are extended to those with spin quantum number s≥21. Solution strategies for the ket- and bra-state equations are also considered. Aspects of extrapolation of CCM expectation values are discussed and future topics regarding extrapolations are presented.