J1-J2 Heisenberg Model Research Articles

Indicators of quantum phase transitions in the one-dimensional J1-J2 Heisenberg model: Global quantum discord and multipartite nonlocality analysis

We here investigate the ground state phase transitions of the one-dimensional J1-J2 Heisenberg model. The phase transition is characterized by analyzing the quantum discord, global quantum discord (GQD), bipartite nonlocality, and multipartite nonlocality of the ground state. Our findings reveal that these correlations serve as effective indicators of the phase transition point at J2/J1=0.5. Specifically, quantum discord and bipartite nonlocality display more pronounced critical behavior in small-sized systems. Conversely, for large-scale systems, GQD and multipartite nonlocality offer a more practical and reliable method for detecting many-body driven phase transitions. Additionally, several linear relationships among these correlations are identified.

A simple linear algebra identity to optimize large-scale neural network quantum states

Neural-network architectures have been increasingly used to represent quantum many-body wave functions. These networks require a large number of variational parameters and are challenging to optimize using traditional methods, as gradient descent. Stochastic reconfiguration (SR) has been effective with a limited number of parameters, but becomes impractical beyond a few thousand parameters. Here, we leverage a simple linear algebra identity to show that SR can be employed even in the deep learning scenario. We demonstrate the effectiveness of our method by optimizing a Deep Transformer architecture with 3 × 105 parameters, achieving state-of-the-art ground-state energy in the J1–J2 Heisenberg model at J2/J1 = 0.5 on the 10 × 10 square lattice, a challenging benchmark in highly-frustrated magnetism. This work marks a significant step forward in the scalability and efficiency of SR for neural-network quantum states, making them a promising method to investigate unknown quantum phases of matter, where other methods struggle.

Extending the variational quantum eigensolver to finite temperatures

We present a variational quantum thermalizer (VQT), called quantum-VQT (qVQT), which extends the variational quantum eigensolver to finite temperatures. The qVQT makes use of an intermediate measurement between two variational circuits to encode a density matrix on a quantum device. A classical optimization provides the thermal state and, simultaneously, all associated excited states of a quantum mechanical system. We demonstrate the capabilities of the qVQT for two different spin systems. First, we analyze the performance of qVQT as a function of the circuit depth and the temperature for a one-dimensional Heisenberg chain. Second, we use the excited states to map the complete, temperature dependent phase diagram of a two-dimensional J1–J2 Heisenberg model. Numerical experiments on both quantum simulators and real quantum hardware demonstrate the efficiency of our approach, which can be readily applied to study various quantum many-body systems at finite temperatures on currently available noisy intermediate-scale quantum devices.

The Ising Nematic in the J1–J2 Heisenberg Model on a Square Lattice in a Self-Consistent Spin-Wave Theory

The Ising Nematic in the J1–J2 Heisenberg Model on a Square Lattice in a Self-Consistent Spin-Wave Theory

Coupled Cluster and Perturbation Theories Based on a Cluster Mean-Field Reference Applied to Strongly Correlated Spin Systems.

We introduce perturbation and coupled-cluster theories based on a cluster mean-field reference for describing the ground state of strongly correlated spin systems. In cluster mean-field, the ground state wave function is written as a simple tensor product of optimized cluster states. The cluster language and the mean-field nature of the ansatz allow for a straightforward improvement which uses perturbation theory and coupled-cluster to account for intercluster correlations. We present benchmark calculations on the 1D chain and 2D square J1-J2 Heisenberg model, using cluster mean-field, perturbation theory, and coupled-cluster. We also present an extrapolation scheme that allows us to compute thermodynamic limit energies accurately. Our results indicate that, with sufficiently large clusters, the correlated methods (cPT2, cPT4, and cCCSD) can provide a relatively accurate description of the Heisenberg model in the regimes considered, which suggests that the methods presented can be used for other strongly correlated systems. Some ways to improve upon the methods presented in this work are discussed.

Momentum-resolved time evolution with matrix product states

We introduce a method based on matrix product states (MPS) for computing spectral functions of (quasi) one-dimensional spin chains, working directly in momentum space in the thermodynamic limit. We simulate the time evolution after applying a momentum operator to an MPS ground state by working with the momentum superposition of a window MPS. We show explicitly for the spin-1 Heisenberg chain that the growth of entanglement is smaller in momentum space, even inside a two-particle continuum, such that we can attain very accurate spectral functions with relatively small bond dimension. We apply our method to compute spectral lineshapes of the gapless XXZ chain and the square-lattice J1-J2 Heisenberg model on a six-leg cylinder.

Gapless quantum spin liquid and global phase diagram of the spin-1/2 [formula omitted] square antiferromagnetic Heisenberg model

The nature of the zero-temperature phase diagram of the spin-1/2J1-J2 Heisenberg model on a square lattice has been debated in the past three decades, and it remains one of the fundamental problems unsettled in the study of quantum many-body theory. By using the state-of-the-art tensor network method, specifically, the finite projected entangled pair state (PEPS) algorithm, to simulate the global phase diagram of the J1-J2 Heisenberg model up to 24×24 sites, we provide very solid evidences to show that the nature of the intermediate nonmagnetic phase is a gapless quantum spin liquid (QSL), whose spin-spin and dimer-dimer correlations both decay with a power law behavior. There also exists a valence-bond solid (VBS) phase in a very narrow region 0.56≲J2/J1≤0.61 before the system enters the well known collinear antiferromagnetic phase. We stress that we make the first detailed comparison between the results of PEPS and the well-established density matrix renormalization group (DMRG) method through one-to-one direct benchmark for small system sizes, and thus give rise to a very solid PEPS calculation beyond DMRG. Our numerical evidences explicitly demonstrate the huge power of PEPS for highly frustrated spin systems. Finally, an effective field theory is also proposed to understand the physical nature of the discovered gapless QSL and its relation to deconfined quantum critical point (DQCP).

Topologically different spin disorder phases of the J1-J2 Heisenberg model on the honeycomb lattice

Searching for spin liquids on the honeycomb J1-J2 Heisenberg model has been attracting great attention in the past decade. In this Paper we investigate the topological properties of the J1-J2 Heisenberg model by introducing nearest-neighbour and next-nearest-neighbour bond parameters. We find that there exist two topologically different phases in the spin disordered regime 0.2<J2/J1<0.5: for J2/J1<0.32, the system is a zero-flux spin liquid which is topological trivial and gapless; for J2/J1>0.32, it is a pi-flux chiral spin liquid, which is topological nontrivial and gapped. These results suggest that there exist two topologically different spin disorder phases in honeycomb J1-J2 Heisenberg model.

Two-dimensional frustrated J1−J2 model studied with neural network quantum states

The use of artificial neural networks to represent quantum wave-functions has recently attracted interest as a way to solve complex many-body problems. The potential of these variational parameterizations has been supported by analytical and numerical evidence in controlled benchmarks. While approaching the end of the early research phase in this field, it becomes increasingly important to show how neural-network states perform for models and physical problems that constitute a clear open challenge for other many-body computational methods. In this paper we start addressing this aspect, concentrating on a presently unsolved model describing two-dimensional frustrated magnets. Using a fully convolutional neural network model as a variational ansatz, we study the frustrated spin-1/2 J1-J2 Heisenberg model on the square lattice. We demonstrate that the resulting predictions for both ground-state energies and properties are competitive with, and often improve upon, existing state-of-the-art methods. In a relatively small region in the parameter space, corresponding to the maximally frustrated regime, our ansatz exhibits comparatively good but not best performance. The gap between the complexity of the models adopted here and those routinely adopted in deep learning applications is, however, still substantial, such that further improvements in future generations of neural-network quantum states are likely to be expected.

Thermal entanglement in a five-qubit XXZ Heisenberg spin chain with the next nearest neighboring interaction

In the study of thermal entanglement of the Heisenberg spin chain model, one usually considers only the spin interaction between the nearest neighboring qubits. Actually, a generalized Heisenberg model, so-called J1-J2 Heisenberg model, which is constructed by considering the fact that not only the nearest neighboring but also the next nearest neighboring spin interaction also plays an important role. In J1-J2 Heisenberg model, due to the next nearest neighboring spin interaction, the frustration effect can occur and has an important influence on the magnetic properties of the model. In this paper we investigate the thermal entanglement of a five-qubit XXZ Heisenberg spin chain with the next nearest neighboring interaction in a magnetic field. Using the numerical method, we calculate the pairwise concurrences of the nearest neighbouring qubits and the next nearest neighboring qubits, abbreviated as C12 and C13 respectively. The numerical results show that the frustration parameter α has an important effect on the pairwise thermal entanglement. Moreover, C12 and C13 have different variations with the change of the frustration parameter α. Meanwhile, it is found that the temperature, magnetic field, Dzyaloshinkii-Moriya (DM) interaction and anisotropic parameter also have great effects on the thermal entanglement. The increasing of temperature can reduce the thermal entanglement. The magnetic field can enhance the thermal entanglement between both two nearest and next nearest neighboring qubits, but when the magnetic field becomes strong enough, only the thermal entanglement between the two nearest neighboring qubits is suppressed. A certain extent of DM interaction can enhance the thermal entanglement between the two nearest neighboring qubits. But for the next nearest neighboring qubits, without the magnetic field, the increasing of DM interaction mainly enlarge the entanglement vanishing area of frustration parameter α. When the system changes from anisotropic to isotropic state, the entanglement vanishing area also changes obviously for C12 and C13. Thus, we can choose appropriate magnetic field strength, temperature, frustration parameter, DM interaction parameter and anisotropic parameter to effectively control and enhance the thermal entanglement of the system.

Multipartite quantum correlations in the extended J1–J2 Heisenberg model

Multipartite entanglement and the maximum violation of Bell inequalities are studied in finite clusters of spins in an extended J1–J2 Heisenberg model at zero temperature. The ensuing highly frustrated states will unveil a rich structure for different values of the corresponding spin–spin interaction strengths. The interplay between nearest-neighbors, next-nearest neighbors and further couplings will be explored using multipartite correlations. The model is relevant to certain quantum annealing computation architectures where an all-to-all connectivity is considered.

Spin–spin correlation length in a two-dimensional frustrated magnet and its relation to doping

In a spherically symmetric self-consistent approach (SSSA), the spin-1/2 J1–J2 Heisenberg model on a two-dimensional square lattice is considered for two-time retarded spin–spin Green’s functions. The spin excitation spectrum, ω(q), and spin gaps at symmetric points are obtained for the entire J1–J2 diagram, i.e., for any ϕ, J1 = cosϕ, and J2 = sinϕ. The structure factor cq and the correlation length ξ at finite temperature are calculated in the entire range of parameters. A radical difference in the behavior of the system in the upper, frustrated (0 ⩽ ϕ ⩽ π), and the lower, nonfrustrated (π ⩽ ϕ ⩽ 2π), regions of the diagram is demonstrated. In the latter region, there is a first-order phase transition that is unique on the phase diagram. For a weakly frustrated antiferromagnet (J1 > J2 > 0), the results obtained are compared with the experimental dependence of ξ on temperature and doping level. A correspondence rule is proposed between frustration in a spin model and the doping of an antiferromagnet with holes.

Ground state of the two-dimensional spin-1/2 J1-J2 Heisenberg model within the Jordan–Wigner fermionization

Ground state of the two-dimensional spin-1/2 <i>J</i>₁-<i>J</i>₂ Heisenberg model within the Jordan–Wigner fermionization

High-temperature Expansion for Frustrated Magnets: Application to the J1-J2 Model on the BCC Lattice

We present the high-temperature expansion up to 11th order for the specific heat C and the uniform susceptibility χ0 and up to 9th order for the structure factor SQ of the frustrated spin-half J1-J2 Heisenberg model on the BCC lattice. We consider ferromagnetic as well as antiferromagnetic nearest-neighbor exchange J1 and frustrating antiferromagnetic next-nearest- neighbor exchange J2. We discuss the influence of frustration on the temperature dependence of these quantities. Furthermore, we use the HTE series to determine the critical temperature Tc as a function of the frustration parameter J2.

A frustrated spin-1 J1-J2 Heisenberg antiferromagnet: An anisotropic planar pyrochlore model

The zero-temperature ground-state (GS) properties and phase diagram of a frustrated spin-1 J1-J2 Heisenberg model on the checkerboard square lattice are studied, using the coupled cluster method. We consider the case where the nearest-neighbour exchange bonds have strength J1 > 0 and the next-nearest-neighbour exchange bonds present (viz., in the checkerboard pattern of the planar pyrochlore) have strength J2 = κJ1 > 0. We find significant differences from both the spin-1/2 and classical versions of the model. We find that the spin-1 model has a first phase transition at κC1 ≈ 1.00 ± 0.01 (as does the classical model at κcl = 1) between two antiferromagnetic phases, viz., a quasiclassical Néel phase (for κ < κC1) and one of the infinitely degenerate family of quasiclassical phases (for k > κC1) that exists in the classical model for κ > κc1, which is now chosen by the order by disorder mechanism as (probably) the "doubled Néel" (or Néel*) state. By contrast, none of this family survives quantum fluctuations to form a stable GS phase in the spin-1/2 case. We also find evidence for a second quantum critical point at κC2 ≈ 2.0 ± 0.5 in the spin-1 model, such that for κ > κC2 the quasiclassical (Néel*) ordering melts and a nonclassical phase appears, which, on the basis of preliminary evidence, appears unlikely to have crossed-dimer valence-bond crystalline (CDVBC) ordering, as in the spin-1/2 case. Unlike in the spin-1/2 case, where the Néel and CDVBC phases are separated by a phase with plaquette valence-bond crystalline (PVBC) ordering, we find very preliminary evidence for such a PVBC state in the spin-1 model for all κ > κC2.

Entanglement spectroscopy of SU(2)-broken phases in two dimensions

In magnetically ordered systems the breaking of SU(2) symmetry in the thermodynamic limit is associated with the appearance of a special type of low-lying excitations in finite size energy spectra, the so called tower of states (TOS). In the present work we numerically demonstrate that there is a correspondence between the SU(2) tower of states and the lower part of the ground state entanglement spectrum (ES). Using state-of-the-art DMRG calculations, we examine the ES of the 2D antiferromagnetic J1-J2 Heisenberg model on both the triangular and kagom\'e lattice. At large ferromagnetic J2 the model exhibits a magnetically ordered ground state. Correspondingly, its ES contains a family of low-lying levels that are reminiscent of the energy tower of states. Their behavior (level counting, finite size scaling in the thermodynamic limit) sharply reflects tower of states features, and is characterized in terms of an effective entanglement Hamiltonian that we provide. At large system sizes TOS levels are divided from the rest by an entanglement gap. Our analysis suggests that (TOS) entanglement spectroscopy provides an alternative tool for detecting and characterizing SU(2)-broken phases using DMRG.

Phase transitions in the 2D J 1-J 2 Heisenberg model with arbitrary signs of exchange interactions

The ground state of the J 1-J 2 Heisenberg model with arbitrary signs of exchange is studied for spin S = 1/2 in the case of the two-dimensional (2D) square lattice. The states with different types of spin long-range order (antiferromagnetic checkerboard, stripe, collinear ferromagnetic) as well as the disordered spin liquid states are described in the framework of one analytical approach. In particular, it is shown that the phase transition between the ferromagnetic spin liquid and the ferromagnet with long-range order is of the second order. In the vicinity of such transition, we have found the ferromagnetic state with a rapidly varying condensate function.

The quasi-one-dimensional character of spin waves in K2Fe7Se8

The ground state phase diagram, the physical properties of spin waves and the dynamic structure factor of the new family of 122 iron-vacancy ordered compound K2Fe7Se8 are theoretically investigated based on the J1–J2 Heisenberg model. Four antiferromagnetic phases are found to be stable in the reasonable physical parameter region. The dispersion relations of the spin waves for the newly discovered AFM2 and AFM3 phases are studied using the linearized spin wave theory. Two acoustic spin wave branches and 12 optical spin wave branches are present in both the AFM2 and AFM3 phases, and each of them are twofold degenerate. One more interesting finding is that the spin waves of K2Fe7Se8 have the quasi-one-dimensional character at the high frequency range. To offer theoretical guidance for the further neutron scattering experiments, the inelastic neutron scattering patterns for the AFM2 and AFM3 phases are also investigated.

THEORETICAL MODELS FOR MAGNETIC PROPERTIES OF IRON PNICTIDES

We attempt to describe the magnetic properties of parent pnictide compounds by using both the J1–J2 Heisenberg model and its three-dimensional generalization, the J1–J2–Jc model. We also include spin anisotropy in the XY plane. In order to obtain the average magnetization and spin wave dispersion, we use the Green's functions method for spin operators in the random phase approximation. We obtain estimates for the model parameters by considering the low temperature experimental dispersion for the compounds CaFe 2 As 2 and BaFe 2 As 2 and conclude that theoretical dispersion can fit the experimental one if spatially anisotropic Hamiltonian is used. A good agreement between theory and experiment indicates that the Heisenberg model is applicable to parent pnictides at low temperatures. The applicability of the model for higher temperatures is checked by calculating the Néel temperature for both compounds. It turns out that the model overestimates the measured critical temperature. The Heisenberg model is not applicable to parent pnictides, for temperature comparable to Néel temperature. Our results thus confirm that all the magnetic properties of parent pnictides cannot be described with purely localized degrees of freedom, and that the itinerant magnetism should have an important role in these compounds. All results given in Sec. 3 are general and could be used in description of classes of compounds with spin stripe structure.

Field-induced staggered moment stabilization in frustrated quantum magnets

For low-dimensional frustrated quantum magnets, the dependence of the staggered moment on a magnetic field is nonmonotonic: For small and intermediate fields, quantum fluctuations are gradually suppressed, leading to an increase of the staggered moment as a function of the field strength. For large applied magnetic fields, the classically expected field dependence is recovered, namely a monotonous decrease with increasing field strength. The staggered moment is eventually suppressed when reaching the fully polarized state at the saturation field. The quantitative analysis of this behavior is an excellent tool to determine the frustration parameter of a magnetic compound. We have developed a general finite-size scaling scheme for numerical exact-diagonalization data of low-dimensional frustrated magnets, which we apply to the recently measured field dependence of the magnetic neutron scattering intensity of Cu(pz)2(ClO4)2 in the framework of the S = 1/2 two-dimensional (2D) J 1–J 2 Heisenberg model. We also apply linear spin-wave theory to complement our numerical findings. Our results show that Cu(pz)2(ClO4)2 is a quasi-2D antiferromagnet with intermediate frustration J 2/J 1 = 0.2.