Ising Variables Research Articles

Correlation equalities and some upper bounds for the coupling constant implying area decay of Wilson loop for Z 3 lattice gauge theories

Correlation identities are obtained for a Z 3 lattice gauge theory where the bonds of the plaquettes are decorated by generalized three-state Ising variables. Making use of correlation inequalities, we obtain rigorous bounds to the area decay of the Wilson loop observable in a range of the coupling parameter larger than those obtained from mean field theory considerations.

To the Theory of Decaying Turbulence

We have found an infinite dimensional manifold of exact solutions of the Navier-Stokes loop equation for the Wilson loop in decaying Turbulence in arbitrary dimension d>2. This solution family is equivalent to a fractal curve in complex space Cd with random steps parametrized by N Ising variables σi=±1, in addition to a rational number pq and an integer winding number r, related by ∑σi=qr. This equivalence provides a dual theory describing a strong turbulent phase of the Navier-Stokes flow in Rd space as a random geometry in a different space, like ADS/CFT correspondence in gauge theory. From a mathematical point of view, this theory implements a stochastic solution of the unforced Navier-Stokes equations. For a theoretical physicist, this is a quantum statistical system with integer-valued parameters, satisfying some number theory constraints. Its long-range interaction leads to critical phenomena when its size N→∞ or its chemical potential μ→0. The system with fixed N has different asymptotics at odd and even N→∞, but the limit μ→0 is well defined. The energy dissipation rate is analytically calculated as a function of μ using methods of number theory. It grows as ν/μ2 in the continuum limit μ→0, leading to anomalous dissipation at μ∝ν→0. The same method is used to compute all the local vorticity distribution, which has no continuum limit but is renormalizable in the sense that infinities can be absorbed into the redefinition of the parameters. The small perturbation of the fixed manifold satisfies the linear equation we solved in a general form. This perturbation decays as t−λ, with a continuous spectrum of indexes λ in the local limit μ→0. The spectrum is determined by a resolvent, which is represented as an infinite product of 3⊗3 matrices depending of the element of the Euler ensemble.

Optimal Thresholds for Fracton Codes and Random Spin Models with Subsystem Symmetry.

Fracton models provide examples of novel gapped quantum phases of matter that host intrinsically immobile excitations and therefore lie beyond the conventional notion of topological order. Here, we calculate optimal error thresholds for quantum error correcting codes based on fracton models. By mapping the error-correction process for bit-flip and phase-flip noises into novel statistical models with Ising variables and random multibody couplings, we obtain models that exhibit an unconventional subsystem symmetry instead of a more usual global symmetry. We perform large-scale parallel tempering MonteCarlo simulations to obtain disorder-temperature phase diagrams, which are then used to predict optimal error thresholds for the corresponding fracton code. Remarkably, we found that the X-cube fracton code displays a minimum error threshold (7.5%) that is much higher than 3D topological codes such as the toric code (3.3%), or the color code (1.9%). This result, together with the predicted absence of glass order at the Nishimori line, shows great potential for fracton phases to be used as quantum memory platforms.

Machine learning for percolation utilizing auxiliary Ising variables.

Machine learning for phase transition has received intensive research interest in recent years. However, its application in percolation still remains challenging. We propose an auxiliary Ising mapping method for the machine learning study of the standard percolation as well as a variety of statistical mechanical systems in correlated percolation representation. We demonstrate that unsupervised machine learning is able to accurately locate the percolation threshold, independent of the spatial dimension of system or the type of phase transition, which can be first-order or continuous. Moreover, we show that, by neural network machine learning, auxiliary Ising configurations for different universalities can be classified with a high confidence level. Our results indicate that the auxiliary Ising mapping method, despite its simplicity, can advance the application of machine learning in statistical and condensed-matter physics.

Dynamic phase transition and universality in a quasi 2D system: Bilayer Ising/Blume-Capel ferromagnet on a honeycomb lattice

Monte Carlo simulation technique is applied on a magnetic bilayer system composed of spin-1/2 Ising and spin-1 Blume-Capel spin variables distributed on a honeycomb lattice. Assuming an indirect exchange interaction between adjacent monolayers, we investigate the influence of single ion anisotropy and number of nonmagnetic layers between the magnetic sheets on the critical behavior of the system. Equilibrium phase diagram in a reduced pseudo critical temperature (Tc/J) versus reduced single ion anisotropy (DB/J) plane has been determined. In the presence of alternating magnetic field, dynamic phase transition properties have been elucidated. The form of the oscillating magnetic field is assumed to be in a square wave form. By performing a detailed finite-size scaling analysis, we have calculated the critical half period t1/2c at which the system undergoes a dynamic phase transition between dynamically ordered and disordered states. Extracted critical exponent ratios β/ν and γ/ν suggest that observed dynamic phase transition for the present multilayer system falls into the same universality class as the equilibrium two-dimensional Ising model.

Spiral order from orientationally correlated random bonds in classical XY models

We discuss the stability of ferromagnetic long-range order in three-dimensional classical XY ferromagnets upon substitution of a small subset of equally oriented bonds by impurity bonds, on which the ferromagnetic exchange J_perp > 0 is replaced by a strong antiferromagnetic coupling J_imp < 0. In the presence of a single impurity bond, once the absolute value of the frustrating coupling J_imp < 0 exceeds a threshold J_c > 0, the ground state becomes two-fold degenerate, corresponding to either clockwise or anticlockwise canting of the spins in the vicinity of the impurity bond. In the presence of a small concentration of impurity bonds, the effective low-energy Hamiltonian is that of Ising variables encoding the sense of rotation of the local canting around the impurities. Those degrees of freedom interact through a dipolar interaction mediated by spin waves. A ferromagnetic Ising ground state indicates the instability of the XY ferromagnet towards a spiral state with a wave vector proportional to the concentration of impurity bonds. To analyze under which circumstances such a ground state arises, we study first impurities forming superlattices. For a subclass of those, we can rigorously establish the existence of spiral order. For another class of superlattices, the Ising variables order ferromagnetically in planes perpendicular to the orientation of impurity bonds, but antiferromagnetically parallel to it, which results in a fan-like XY ground state. Second, we consider the case when the impurity bonds are randomly distributed on the three-dimensional host lattice according to a Poisson process. We show the phenomenon of spiral order by disorder with an ordering wave vector proportional to the impurity concentration. The analytical predictions are confirmed by Monte Carlo simulations and are relevant for magnetic materials such as YBaCuFeO_5.

A Short Path Quantum Algorithm for Exact Optimization

We give a quantum algorithm to exactly solve certain problems in combinatorial optimization, including weighted MAX-2-SAT as well as problems where the objective function is a weighted sum of products of Ising variables, all terms of the same degree D; this problem is called weighted MAX-ED-LIN2. We require that the optimal solution be unique for odd D and doubly degenerate for even D; however, we expect that the algorithm still works without this condition and we show how to reduce to the case without this assumption at the cost of an additional overhead. While the time required is still exponential, the algorithm provably outperforms Grover's algorithm assuming a mild condition on the number of low energy states of the target Hamiltonian. The detailed analysis of the runtime dependence on a tradeoff between the number of such states and algorithm speed: fewer such states allows a greater speedup. This leads to a natural hybrid algorithm that finds either an exact or approximate solution.

Semiconductor of spinons: from Ising band insulator to orthogonal band insulator

We use the ionic Hubbard model to study the effects of strong correlations on a two-dimensional semiconductor. The spectral gap in the limit where on-site interactions are zero is set by the staggered ionic potential, while in the strong interaction limit it is set by the Hubbard U. Combining mean field solutions of the slave spin and slave rotor methods, we propose two interesting gapped phases in between: (i) the insulating phase before the Mott phase can be viewed as gapping a non-Fermi liquid state of spinons by the staggered ionic potential. The quasi-particles of underlying spinons are orthogonal to physical electrons, giving rise to the ‘ARPES-dark’ state where the ARPES gap will be larger than the optical and thermal gap. (ii) The Ising insulator corresponding to ordered phase of the Ising variable is characterized by single-particle excitations whose dispersion is controlled by Ising-like temperature and field dependences. The temperature can be conveniently employed to drive a phase transition between these two insulating phases where Ising exponents become measurable by ARPES and cyclotron resonance. The rare earth monochalcogenide semiconductors where the magneto-resistance is anomalously large can be a candidate system for the Ising band insulator. We argue that the Ising and orthogonal insulating phases require strong enough ionic potential to survive the downward renormalization of the ionic potential caused by Hubbard U.

Unbinding slave spins in the Anderson impurity model

We show that a generic single-orbital Anderson impurity model, lacking for instance any kind of particle-hole symmetry, can be exactly mapped without any constraint onto a resonant level model coupled to two Ising variables, which reduce to one if the hybridisation is particle-hole symmetric. The mean-field solution of this model is found to be stable to unphysical spontaneous magnetisation of the impurity, unlike the saddle-point solution in the standard slave-boson representation. Remarkably, the mean-field estimate of the Wilson ratio approaches the exact value $R_\text{W}=2$ in the Kondo regime.

Classical Spin Liquid Instability Driven By Off-Diagonal Exchange in Strong Spin-Orbit Magnets.

We show that the off-diagonal exchange anisotropy drives Mott insulators with strong spin-orbit coupling to a classical spin liquid regime, characterized by an infinite number of ground states and Ising variables living on closed or open strings. Depending on the sign of the anisotropy, quantum fluctuations either fail to lift the degeneracy down to very low temperatures, or select noncoplanar magnetic states with unconventional spin correlations. The results apply to all 2D and 3D tricoordinated materials with bond-directional anisotropy and provide a consistent interpretation of the suppression of the x-ray magnetic circular dichroism signal reported recently for β-Li_{2}IrO_{3} under pressure.

Spanning forests and OSP(N|2M) -invariant σ-models**Dedicated to our friend Tony Guttmann on the occasion of his 70th birthday

The present paper is part of our ongoing work on supersymmetric σ-models, their relation with the Potts model at q = 0 and spanning forests, and the rigorous analytic continuation of the partition function as an entire function of N − 2M, a feature first predicted by Parisi and Sourlas in the 1970s. Here we accomplish two main steps. First, we analyze in detail the role of the Ising variables that arise when the constraint in the model is solved, and we point out two situations in which the Ising and forest variables decouple. Second, we establish the analytic continuation for the model in some special cases: when the underlying graph is a forest, and for the Nienhuis action on a cubic graph. We also make progress in understanding the series-parallel graphs.

Anisotropic four-state clock model in the presence of random fields.

A four-state clock ferromagnetic model is studied in the presence of different configurations of anisotropies and random fields. The model is considered in the limit of infinite-range interactions, for which the mean-field approach becomes exact. Both representations of Cartesian spin components and two Ising variables are used, in terms of which the physical properties and phase diagrams are discussed. The random fields follow bimodal probability distributions and the richest criticality is found when the fields, applied in the two Ising systems, are not correlated. The phase diagrams present new interesting topologies, with a wide variety of critical points, which are expected to be useful in describing different complex phenomena.

Pseudolikelihood Decimation Algorithm Improving the Inference of the Interaction Network in a General Class of Ising Models

In this Letter we propose a new method to infer the topology of the interaction network in pairwise models with Ising variables. By using the pseudolikelihood method (PLM) at high temperature, it is generally possible to distinguish between zero and nonzero couplings because a clear gap separate the two groups. However at lower temperatures the PLM is much less effective and the result depends on subjective choices, such as the value of the ℓ1 regularizer and that of the threshold to separate nonzero couplings from null ones. We introduce a decimation procedure based on the PLM that recursively sets to zero the less significant couplings, until the variation of the pseudolikelihood signals that relevant couplings are being removed. The new method is fully automated and does not require any subjective choice by the user. Numerical tests have been performed on a wide class of Ising models, having different topologies (from random graphs to finite dimensional lattices) and different couplings (both diluted ferromagnets in a field and spin glasses). These numerical results show that the new algorithm performs better than standard PLM.

Engineering Ising-XY spin-models in a triangular lattice using tunable artificial gauge fields

Magnetism plays a key role in modern science and technology, but still many open questions arise from the interplay of magnetic many-body interactions. Deeper insight into complex magnetic behaviour and the nature of magnetic phase transitions can be obtained from, for example, model systems of coupled XY and Ising spins. Here, we report on the experimental realization of such a coupled system with ultracold atoms in triangular optical lattices. This is accomplished by imposing an artificial gauge field on the neutral atoms, which acts on them as a magnetic field does on charged particles. As a result, the atoms show persistent circular currents, the direction of which provides an Ising variable. On this, the tunable staggered gauge field, generated by a periodic driving of the lattice, acts as a longitudinal field. Further, the superfluid ground state presents strong analogies with the paradigm example of the fully frustrated XY model on a triangular lattice. A quantum gas trapped in an optical lattice of triangular symmetry can now be driven from a paramagnetic to an antiferromagnetic state by a tunable artificial magnetic field.

Digitizing the Neveu-Schwarz model on the lightcone worldsheet

The purpose of this article is to extend the light-cone worldsheet lattice description of string theory to include the Neveu-Schwarz model. We model each component of the fermionic worldsheet field by a critical Ising model. We show that a simple choice of boundary conditions for the Ising variables leads to the half-integer modes required by the model. We identify the $G$-parity operation within the Ising model and formulate the procedure for projecting onto the even $G$-parity sector. We construct the lattice version of the three open string vertex, with the necessary operator insertion at the interaction point. We sketch a formalism for summing planar open string multiloop amplitudes, and we discuss prospects for numerically summing them. If successful, the methods described here could provide an alternative to lattice gauge theory for computations in large $N$ QCD.

A schematic model for molecular affinity and binding with Ising variables

After discussing the relevance of statistical physics in molecular recognition processes, we present a schematic model for ligand-receptor association based on an Ising chain. We discuss the possible behaviors of the affinity when the stiffness of the ligand increases. We also consider the case of flexible receptors. A variety of interesting behaviors is obtained, including some affinity modulation upon bond hardening or softening. The affinity of a ligand for its receptor is shown to depend on the details of its rigidity profile, and we question the possibility of encoding information in the rigidities as well as in the shape. An exhaustive study of the selectivity of patterns with length n < 8 is carried out. Connection with other spin models, in particular spin glasses is mentioned in the conclusion.

Topological superfluid 3He-B in magnetic field and ising variable

The topological superfluid 3He-B provides many examples of the interplay of symmetry and topology. Here we consider the effect of magnetic field on topological properties of 3He-B. Magnetic field violates the time reversal symmetry. As a result, the topological invariant supported by this symmetry ceases to exist; and thus the gapless fermions on the surface of 3He-B are not protected any more by topology: they become fully gapped. Nevertheless, if perturbation of symmetry is small, the surface fermions remain relativistic with mass proportional to symmetry violating perturbation—magnetic field. The 3He-B symmetry gives rise to the Ising variable I = ±1, which emerges in magnetic field and which characterizes the states of the surface of 3He-B. This variable also determines the sign of the mass term of surface fermions and the topological invariant describing their effective Hamiltonian. The line on the surface, which separates the surface domains with different I, contains 1 + 1 gapless fermions, which are protected by combined action of symmetry and topology.

Spin-lattice model of magnetoelectric transitions inRbCoBr3

Extensive Monte-Carlo simulations are performed to analyze a recent neutron-diffraction experiment on a distorted triangular lattice compound ${\text{RbCoBr}}_{3}$. We consider a spin-lattice model where both spin and lattice are Ising variables. This model explains well successive magnetic and dielectric transitions observed in the experiment. The exchange interaction parameters and the spin-lattice coupling are estimated. It is found that the spin-lattice coupling is important to explain the slow growth of a ferrimagnetic order. The present simulations were made possible by developing a modified Monte-Carlo algorithm, which accelerates slow Monte-Carlo dynamics of quasi-one-dimensional frustrated systems.

Spin and chiral stiffness of the XY spin glass in two dimensions

We analyse the zero-temperature behaviour of theXY Edwards–Anderson spin glass model on a square lattice. A newly developed algorithmcombining exact ground-state computations for Ising variables embedded into the planarspins with a specially tailored evolutionary method, resulting in the genetic embeddedmatching (GEM) approach, allows for the computation of numerically exact ground statesfor relatively large systems. This enables a thorough reinvestigation of the long-standingquestions of (i) extensive degeneracy of the ground state and (ii) a possible decoupling ofspin and chiral degrees of freedom in such systems. The new algorithm together withappropriate choices for the considered sets of boundary conditions and finite-size scalingtechniques allows for a consistent determination of the spin and chiral stiffness scalingexponents.

Effective Hamiltonian for the pyrochlore antiferromagnet: Semiclassical derivation and degeneracy

In the classical pyrochlore lattice Heisenberg antiferromagnet, there is a macroscopic continuous ground-state degeneracy. We study the semiclassical limit of large spin length $S$, keeping only the lowest-order (in $1∕S$) correction to the classical Hamiltonian. We perform a detailed analysis of the spin-wave modes, and using a real-space loop expansion, we produce an effective Hamiltonian, in which the degrees of freedom are Ising variables representing fluxes through loops in the lattice. We find a family of degenerate collinear ground states, related by gaugelike ${Z}_{2}$ transformations, and provide bounds for the order of the degeneracy. We further show that the theory can readily be applied to determine the ground states of the Heisenberg Hamiltonian on related lattices and to field-induced collinear magnetization plateau states.