The Ising model with nearest-neighbor interactions on a two-dimensional (2D) square lattice is one of the simplest models for studying ferro-magnetic to para-magnetic transitions. Extensive results are available in the literature for this model, which has become a paradigm for the study of magnetic phase transitions in materials, both theoretically and numerically. After a brief review of the main results obtained with a classical computer, we show how to implement on the D- Wave quantum annealer a more complex Ising model with the addition of competing antiferromagnetic interactions between the diagonal next-to-nearest neighbors with two coupling constants J1 and J2. The dynamics of this system, owing to frustration, are richer than those of the simple Ising model and exhibit a third striped (or antiferromagnetic) phase in addition to the ferro- and para-magnetic phases. In this work, we observed all three phases on the D-Wave hardware, studied the behavior of the solution with different annealing parameters, such as the chain strength and annealing time, and showed how to identify the phase transition by varying the ratio between the ferromagnetic and antiferromagnetic couplings. The same system is studied on a classical computer, with the possibility of taking into account the temperature (fixed on D-Wave) as a free parameter and to explore the full phase diagram: some comparative conclusions with D-Wave are drawn.
We study the frustrated Ising model on the two-dimensional
L×L square lattice with ferromagnetic (FM) nearest-neighbor and antiferromagnetic diagonal-neighbor interactions using the D-wave quantum annealing machine (D-QAM) with 5000+ qubits composed on structure of the Pegasus graph. As the former Monte Carlo and mean field results, we find the FM to stripe order phase transition, through observations of the magnetization
M, energy, magnetic susceptibility and structure factor. We also analyze probability which occurs any
M at a given interaction for many quantum annealing shots to estimate the shape of objective function
f. The only one value of
M with specific phase is observed in the regions far from phase transition for many quantum annealing shots, while several values of
M with different possibilities are appeared in the regimes of phase transition. We guess that
f in the regimes of phase transition retains the multi-modal structure with several local minimums, due to the strong degeneracies caused by frustrations. Finally, we discuss fail of the quantum annealing simulations, through analysis of the number of the chains, defined as the same variable with
N-qubits, as a function of
L.
In this work the ordering processes in classical frustrated Ising models and the ground-state behavior of frustrated quantum spin models are investigated. In the first part the phase transition from a high-temperature paramagnetic phase to a antiferromagnetic ground state is analyzed for the square-lattice Ising model with competing antiferromagnetic interactions between nearest- and next-nearest-neighbor spins. In the case of spatially isotropic interactions the phase transition to a collinear state is of special interest. Using Monte-Carlo simulations on finite lattices as well as field-theoretical methods for the continuous model it is shown that depending on the frustration three different sectors can be distinguished that show (i) Ising-universal behavior (ii) Ashkin-Teller-like non-universal behavior (varying critical exponents) and (iii) non-critical first-order behavior. For spatially anisotropic interactions the emergence of an intermediate incommensurate phase is verified by means of Monte-Carlo simulations which has been predicted by earlier mean-field calculations. The nature of this phase is studied analyzing spin-spin correlation functions and its ordering wave-vector reveals a strong temperature and size dependence. In the second part the influence of quantum fluctuations on classical antiferromagnetic ground states are investigated in anisotropic quantum spin-1/2 Heisenberg models with competing interactions. Using quantum Monte-Carlo simulations, exact diagonalization and perturbation theory the ground-state phase diagrams are calculated for a square lattice with nearest- and next-nearest-neighbor interactions and a honeycomb lattice with up to third-nearest-neighbor interactions whereby the fluctuations are chosen to be non-frustrating. For both cases the stability regions of classical magnetic states are identified and an intermediate phase is detected that shows no sign of long-range magnetic order. Further calculations also rule out valence-bond-crystal or topological order for this state, thus a quantum disordered ground state is stabilized in these models.
One of the surprising, and potentially very useful, capabilities of analog quantum computers, such as D-Wave quantum annealers, is sampling from the Boltzmann, or Gibbs, distribution defined by a classical Hamiltonian. In this study, we thoroughly examine the ability of D-Wave quantum annealers to sample from the Boltzmann distribution defined by a canonical type of competing magnetic frustration
J1-
J2 model; the 1-dimensional ANNNI (axial next-nearest-neighbor Ising) model. Boltzmann sampling error rate is quantified for standard linear-ramp anneals ranging from
5 nanosecond annealing times up to
2000 microseconds on two different D-Wave quantum annealing processors. Interestingly, we find some analog hardware parameters which result in a very high accuracy (down to a TVD of
0.0003) and low temperature sampling (down to
β=32.2) in a frustrated region of the ANNNI model magnetic phase diagram. This bolsters the viability of current analog quantum computers for thermodynamic sampling applications of highly frustrated magnetic spin systems.
We consider the thermal phase transition from a paramagnetic to a stripe-antiferromagnetic phase in the frustrated two-dimensional square-lattice Ising model with competing interactions
J1<0 (nearest neighbor, ferromagnetic) and
J2>0 (second neighbor, antiferromagnetic). The striped phase breaks a
Z4 symmetry and is stabilized at low temperatures for
g=J2/∣J1∣>1/2. Despite the simplicity of the model, it has proved difficult to precisely determine the order and the universality class of the phase transitions. This was done convincingly only recently by Jin et al. [Phys. Rev. Lett. 108, 045702 (2012)]. Here, we further elucidate the nature of these transitions and their anomalies by employing a combination of cluster mean-field theory, Monte Carlo simulations, and transfer-matrix calculations. The
J1-
J2 model has a line of very weak first-order phase transitions in the whole region
1/2<g<g∗, where ${g}^{*}=0.67\ifmmode\pm\else\textpm\fi{}0.01$. Thereafter, the transitions from
g=g∗ to $g\ensuremath{\rightarrow}\ensuremath{\infty}$ are continuous and can be fully mapped, using universality arguments, to the critical line of the well-known Ashkin-Teller model from its four-state Potts point to the decoupled Ising limit. We also comment on the pseudo-first-order behavior at the Potts point and its neighborhood in the Ashkin-Teller model on finite lattices, which in turn leads to the appearance of similar effects in the vicinity of the multicritical point
g∗ in the
J1-
J2 model. The continuous transitions near
g∗ can therefore be mistaken for first-order transitions, and this realization was the key to understanding the paramagnetic-striped transition for the full range of
g>1/2. Most of our conclusions are based on Monte Carlo calculations, while the cluster mean-field and transfer-matrix results provide useful methodological benchmarks for weakly first-order behaviors and Ashkin-Teller criticality.
The phase transitions in two-dimensional antiferromagnetic Ising model are studied on a square lattice by taking the interactions of next-nearest neighbors into account. The model is investigated on basis of the replica Monte Carlo algorithm and the histogrammic analysis of data. The diagram of the critical temperature dependence on an interaction value of next-nearest neighbors is plotted. The studied model reveals the phase transition of second order and the frustration points depending on relative interaction values and magnetic field. A effect, heat capacity splitting near the frustration points, is found.
In the present work, we investigate the phase transitions of the Ising model with interactions between first (J1), second (J2) and third (J3) neighbours in the square lattice with transverse magnetic field. In this study, we adopt antiferromagnetic interactions between first and second neighbours and consider third-neighbours interactions to be both ferromagnetic and antiferromagnetic. The description of the classical and quantum phase transitions of the model is carried out by adopting the cluster mean-field approximation with four sites. As a result, we identified that strong enough third-neighbor interactions lead to the disappearance of tricriticality at the boundary between the superantiferromagnetic (SAF) and paramagnetic (PM) phases. In particular, tricriticality is more sensitive to the presence of ferromagnetic interactions, disappearing for |J3| ă 0.3|J1| in the absence of transverse field. In the presence of antiferromagnetic J3, we find first-order transitions between the degenerate staggered dimer phase and the PM phase. Furthermore, we find a change in the nature of the SAF-PM phase transitions introduced by quantum fluctuations. This phenomenon, called quantum annealed criticality (QAC), consists of classical first-order phase transitions that become second-order phase transitions in the presence of a strong enough transverse field. Our results allow establishing a range of parameters J2 and J3 for which it is possible to find QAC. The analysis of the model criticality allows concluding that a strong enough J3 interaction eliminates the QAC phenomenon. Therefore, our results suggest that attempts to realize QAC in a square lattice system should avoid strong interactions between third neighbors, mainly if these interactions are ferromagnetic.