Ising Universality Class Research Articles

Critical Casimir forces between defects in the 2D Ising model

An exact statistical mechanical derivation is given of the critical Casimir interactions between two defects in a planar lattice-gas Ising model. Each defect is a finite group of nearest-neighbor spins with modified coupling constants. Such a system can be regarded as a model of a binary liquid mixture with the molecules confined to a membrane and the defects mimicking protein inclusions embedded into the membrane. As suggested by recent experiments, certain cellular membranes appear to be tuned to the proximity of a critical demixing point belonging to the two-dimensional Ising universality class. Therefore one can expect the emergence of critical Casimir forces between membrane inclusions. These forces are governed by universal scaling functions, which we derive for simple defects. We prove that the scaling law appearing at criticality is the same for all types of defects considered here.

Charge order in an interacting monolayer under transverse bias

A monolayer of molecules or quantum dots sandwiched between electrodes can be driven out of equilibrium by the application of a bias voltage between the electrodes. We study charge ordering, i.e., the spontaneous formation of a charge density wave, and the perpendicular current in such a system within a master-equation approach augmented by mean-field and classical Monte Carlo methods. Our approach is suitable for weak tunneling between the monolayer and the electrodes. For a square lattice with nearest-neighbor Coulomb repulsion, we present a comprehensive study of the zero-temperature phases controlled by the on-site energy, the bias voltage, and the degeneracy of the occupied single-site state. One of the most interesting results is the prediction of a conducting charge-density-wave phase that only occurs at a finite bias voltage. We also study the universality classes of the phase transitions towards charge-ordered states at zero and nonzero temperatures. While all transitions at $T>0$ and some at $T=0$ belong to the two-dimensional Ising universality class, we also find an absorbing-to-active phase transition in the $\mathbb{Z}_2$ symmetric directed percolation (DP2) class at $T=0$.

Structure and dynamics of binary liquid mixtures near their continuous demixing transitions

The dynamic and static critical behavior of a family of binary Lennard-Jones liquid mixtures, close to their continuous demixing points (belonging to the so-called model H' dynamic universality class), are studied computationally by combining semi-grand canonical Monte Carlo simulations and large-scale molecular dynamics (MD) simulations, accelerated by graphic processing units (GPU). The symmetric binary liquid mixtures considered cover a variety of densities, a wide range of compressibilities, and various interactions between the unlike particles. The static quantities studied here encompass the bulk phase diagram (including both the binodal and the λ-line), the correlation length, and the concentration susceptibility, of the finite-sized systems above the bulk critical temperature Tc, the compressibility and the pressure at Tc. Concerning the collective transport properties, we focus on the Onsager coefficient and the shear viscosity. The critical power-law singularities of these quantities are analyzed in the mixed phase (above Tc) and non-universal critical amplitudes are extracted. Two universal amplitude ratios are calculated. The first one involves static amplitudes only and agrees well with the expectations for the three-dimensional Ising universality class. The second ratio includes also dynamic critical amplitudes and is related to the Einstein-Kawasaki relation for the interdiffusion constant. Precise estimates of this amplitude ratio are difficult to obtain from MD simulations, but within the error bars our results are compatible with theoretical predictions and experimental values for model H'. Evidence is reported for an inverse proportionality of the pressure and the isothermal compressibility at the demixing transition, upon varying either the number density or the repulsion strength between unlike particles.

Phase transitions and multicritical behavior in the Ising model with dipolar interactions.

In this work, the phase diagram of the ferromagnetic Ising model with dipole interactions is revisited with the aim of determining the nature of the phase transition between stripe-ordered phases with width n (h_{n}) and tetragonal liquid (TL) phases. Extensive Monte Carlo simulations are performed in order to study the short-time dynamic behavior of the observables for selected values of the ratio between the ferromagnetic exchange and dipolar constants δ. The obtained results indicate that the h_{1}-TL phase transition line is continuous up to δ=1.2585, while for the h_{2}-TL line a weak first-order character is found in the interval 1.2585≤δ≤1.36 and becomes continuous for 1.37≤δ≤1.9. This result suggests the existence of a tricritical point close to δ=1.37. When it is appropriate, the complete set of critical exponents is obtained, and in all the studied cases they depend on δ but do not belong to the Ising universality class. Furthermore, short-time dynamic studies reveal that at the point where the mentioned lines meet the h_{1}-h_{2} line, i.e., at δ=1.2585, the critical phase corresponding to the h_{1}-TL transition coexists with the h_{2} phase.

Truncated Hilbert space approach to the 2d ϕ 4 theory

We apply the massive analogue of the truncated conformal space approach to study the two dimensional $\phi^{4}$ theory in finite volume. We focus on the broken phase and determine the finite size spectrum of the model numerically. We interpret the results in terms of the Bethe-Yang spectrum, from which we extract the infinite volume masses and scattering matrices for various couplings. We compare these results against semiclassical analysis and perturbation theory. We also analyze the critical point of the model and confirm that it is in the Ising universality class.

Critical properties of the eight-vertex model in a field

The general eight-vertex model on a square lattice is studied numerically by using the Corner Transfer Matrix Renormalization Group method. The method is tested on the symmetric (zero-field) version of the model, the obtained dependence of critical exponents on model's parameters is in agreement with Baxter's exact solution and weak universality is verified with a high accuracy. It was suggested long time ago that the symmetric eight-vertex model is a special exceptional case and in the presence of external fields the eight-vertex model falls into the Ising universality class. We confirm numerically this conjecture in a subspace of vertex weights, except for two specific combinations of vertical and horizontal fields for which the system still exhibits weak universality.

Behavior and finite-size effects of the sixth order cumulant in the three-dimensional Ising universality class**Supported by Fund Project of Sichuan Provincial Department of Education (16ZB0339), Fund Project of Chengdu Technological University for Doctor (2016RC004), Major State Basic Research Development Program of China (2014CB845402) and National Natural Science Foundation of China (11405088, 11221504)

The high-order cumulants of conserved charges are suggested to be sensitive observables to search for the critical point of Quantum Chromodynamics (QCD). This has been calculated to the sixth order in experiments. Corresponding theoretical studies on the sixth order cumulant are necessary. Based on the universality of the critical behavior, we study the temperature dependence of the sixth order cumulant of the order parameter using the parametric representation of the three-dimensional Ising model, which is expected to be in the same universality class as QCD. The density plot of the sign of the sixth order cumulant is shown on the temperature and external magnetic field plane. We found that at non-zero external magnetic field, when the critical point is approached from the crossover side, the sixth order cumulant has a negative valley. The width of the negative valley narrows with decreasing external field. Qualitatively, the trend is similar to the result of Monte Carlo simulation on a finite-size system. Quantitatively, the temperature of the sign change is different. Through Monte Carlo simulation of the Ising model, we calculated the sixth order cumulant of different sizes of systems. We discuss the finite-size effects on the temperature at which the cumulant changes sign.

Effective theory and breakdown of conformal symmetry in a long-range quantum chain

We deal with the problem of studying the symmetries and the effective theories of long-range models around their critical points. A prominent issue is to determine whether they possess (or not) conformal symmetry (CS) at criticality and how the presence of CS depends on the range of the interactions. To have a model, both simple to treat and interesting, where to investigate these questions, we focus on the Kitaev chain with long-range pairings decaying with distance as power-law with exponent α. This is a quadratic solvable model, yet displaying non-trivial quantum phase transitions. Two critical lines are found, occurring respectively at a positive and a negative chemical potential. Focusing first on the critical line at positive chemical potential, by means of a renormalization group approach we derive its effective theory close to criticality. Our main result is that the effective action is the sum of two terms: a Dirac action SD, found in the short-range Ising universality class, and an “anomalous” CS breaking term SAN. While SD originates from low-energy excitations in the spectrum, SAN originates from the higher energy modes where singularities develop, due to the long-range nature of the model. At criticality SAN flows to zero for α>2, while for α<2 it dominates and determines the breakdown of the CS. Out of criticality SAN breaks, in the considered approximation, the effective Lorentz invariance (ELI) for every finite α. As α increases such ELI breakdown becomes less and less pronounced and in the short-range limit α→∞ the ELI is restored. In order to test the validity of the determined effective theory, we compared the two-fermion static correlation functions and the von Neumann entropy obtained from them with the ones calculated on the lattice, finding agreement. These results explain two observed features characteristic of long-range models, the hybrid decay of static correlation functions within gapped phases and the area-law violation for the von Neumann entropy. The proposed scenario is expected to hold in other long-range models displaying quasiparticle excitations in ballistic regime. From the effective theory one can also see that new phases emerge for α<1. Finally we show that at every finite α the critical exponents, defined as for the short-range (α→∞) model, are not altered. This also shows that the long-range paired Kitaev chain provides an example of a long-range model in which the value of α where the CS is broken does not coincide with the value at which the critical exponents start to differ from the ones of the corresponding short-range model. At variance, for the second critical line, having negative chemical potential, only SAN (SD) is present for 1<α<2 (for α>2). Close to this line, where the minimum of the spectrum coincides with the momentum where singularities develop, the critical exponents change where CS is broken.

Bak–Tang–Wiesenfeld model on the square site-percolation lattice

The Bak–Tang–Wiesenfeld (BTW) model is considered on the site-diluted square lattice, tuned by the occupancy probability p. Various statistical observables of the avalanches are analyzed in terms of p, e.g. the fractal dimension of their exterior frontiers, gyration radius, loop lengths and Green’s function. The model exhibits critical behavior for all amounts of p, and the exponents of the statistical observables are analyzed. We find a distinct universality class at , which is unstable towards a p = 1 (BTW) fixed point. This universality class displays some common features such as a two-dimensional (2D) Ising universality class, e.g. the fractal dimension of loops in the thermodynamic limit is which is compatible with the fractal dimension of geometrical spin clusters of the 2D critical Ising model (with ).

Conformal perturbation of off-critical correlators in the 3D Ising universality class

Thanks to the impressive progress of conformal bootstrap methods we have now very precise estimates of both scaling dimensions and OPE coefficients for several 3D universality classes. We show how to use this information to obtain similarly precise estimates for off-critical correlators using conformal perturbation. We discuss in particular the $< \sigma (r) \sigma (0) >$, $< \epsilon (r) \epsilon (0) >$ and $< \sigma (r) \epsilon (0) >$ two point functions in the high and low temperature regimes of the 3D Ising model and evaluate the leading and next to leading terms in the $s = t r^{\Delta_{t}}$ expansion, where $t$ is the reduced temperature. Our results for $< \sigma (r) \sigma (0) >$ agree both with Monte Carlo simulations and with a set of experimental estimates of the critical scattering function.

Classification of magnons in rotated ferromagnetic Heisenberg model and their competing responses in transverse fields

Competing orders is a general concept to describe various quantum phases and transitions in various materials. One efficient way to investigate competing orders is to first classify different class of excitations in a given quantum phase, then study their competing responses under various external probes. This strategy may not only lead to deep understanding of the quantum phase itself, but also its deep connections to various other quantum phases nearby. We implement this approach by studying the Rotated Ferromagnetic Heisenberg model (RFHM) in two different transverse fields $h_x$ and $h_z$ which can be intuitively visualized as studying spin-orbit couplings (SOC) effects in 2d Ising or anisotropic XY model in a transverse field. At a special SOC class, it was known that the RFHM at a zero field owns an exact ground state called Y-x state. It supports non only the commensurate C-C$_0$ and C-C$_{\pi} $ magnons, but also the in-commensurate C-IC magnons. These magnons are non-relativistic, not contained in the exact ground state, so need to be thermally excited. Their dramatic response under the longitudinal $ h_y $ field was recently worked out by the authors. Here we find they respond very differently under the two transverse fields. Any $h_x$ ($h_z$) changes the collinear Y-x state to a canted co-planar YX-x (YZ-x) state which suffers quantum fluctuations. The C-C$_0$, C-C$_{\pi} $ and C-IC magnons sneak into the quantum ground state, become relativistic and play leading roles even at $ T=0 $. We map out the boundaries among the C-C$_0$, C-C$_{\pi} $ and C-IC magnons, especially the detailed evolution of the C-IC magnons inside the canted phases. As $h_x$ ($h_z$) increases further, the C-C$_0$ magnons always win the competition and emerge as the seeds to drive a transition from the YX-x (YZ-x) to the X-FM ( Z-FM ) which is shown to be in the 3d Ising universality class.

Molecular dynamics simulation of a binary mixture near the lower critical point.

2,6-lutidine molecules mix with water at high and low temperatures but in a wide intermediate temperature range a 2,6-lutidine/water mixture exhibits a miscibility gap. We constructed and validated an atomistic model for 2,6-lutidine and performed molecular dynamics simulations of 2,6-lutidine/water mixture at different temperatures. We determined the part of demixing curve with the lower critical point. The lower critical point extracted from our data is located close to the experimental one. The estimates for critical exponents obtained from our simulations are in a good agreement with the values corresponding to the 3D Ising universality class.

Universal fidelity near quantum and topological phase transitions in finite one-dimensional systems

We study the quantum fidelity (groundstate overlap) near quantum phase transitions of the Ising universality class in one dimensional (1D) systems of finite size L. Prominent examples occur in magnetic systems (e.g. spin-Peierls, the anisotropic XY model), and in 1D topological insulators of any topologically nontrivial Altland-Zirnbauer-Kitaev universality class. The rescaled fidelity susceptibility is a function of the only dimensionless parameter LM, where 2M is the gap in the fermionic spectrum. We present analytic expressions for the fidelity susceptibility for periodic and open boundaries conditions with zero, one or two edgestates. The latter are shown to have a crucial impact and alter the susceptibility both quantitatively and qualitatively. We support our analytical solutions with numerical data.

Dimerization transitions in spin-1 chains

We study spontaneous dimerization transitions in a Heisenberg spin-1 chain with additional next-nearest neighbor (NNN) and 3-site interactions using extensive numerical simulations and a conformal field theory analysis. We show that the transition can be second order in the WZW SU(2)$_2$ or Ising universality class, or first-order. We argue that these features are generic because of a marginal operator in the WZW SU(2)$_2$ model, and because of two topologically distinct non-dimerized phases with or without edge states. We also provide explicit numerical evidence of conformal towers of singlets inside the spin gap at the Ising transition. Implications for other models are briefly discussed.

Noise-induced absorbing phase transition in a model of opinion formation

In this work we study a 3-state (+1, −1, 0) opinion model in the presence of noise and disorder. We consider pairwise competitive interactions, with a fraction p of those interactions being negative (disorder). Moreover, there is a noise q that represents the probability of an individual spontaneously change his opinion to the neutral state. Our aim is to study how the increase/decrease of the fraction of neutral agents affects the critical behavior of the system and the evolution of opinions. We derive analytical expressions for the order parameter of the model, as well as for the stationary fraction of each opinion, and we show that there are distinct phase transitions. One is the usual ferro–paramagnetic transition, that is in the Ising universality class. In addition, there are para-absorbing and ferro-absorbing transitions, presenting the directed percolation universality class. Our results are complemented by numerical simulations.

Tricritical Ising phase transition in a two-ladder Majorana fermion lattice

We introduce a two-ladder lattice model with interacting Majorana fermions that could be realized on the surfaces of a topological insulator film. We study this model by a combination of analytical and numerical techniques and find a phase diagram that features both gapless and gapped phases as well as interesting phase transitions including a quantum critical point in the tricritical Ising (TCI) universality class. The latter occurs at an intermediate coupling strength at a meeting point of a first-order transition line and an Ising critical line and is known to be described by a superconformal field theory with central charge $c =\frac{7}{10}$. We discuss the experimental feasibility of constructing the model and tuning parameters to the vicinity of the TCI point where signatures of the elusive supersymmetry can be observed.

Majority-vote model on spatially embedded networks: Crossover from mean-field to Ising universality classes

We study through Monte Carlo simulations and finite-size scaling analysis the nonequilibrium phase transitions of the majority-vote model taking place on spatially embedded networks. These structures are built from an underlying regular lattice over which directed long-range connections are randomly added according to the probability P_{ij}∼r^{-α}, where r_{ij} is the Manhattan distance between nodes i and j, and the exponent α is a controlling parameter [J. M. Kleinberg, Nature (London) 406, 845 (2000)NATUAS0028-083610.1038/35022643]. Our results show that the collective behavior of this system exhibits a continuous order-disorder phase transition at a critical parameter, which is a decreasing function of the exponent α. Precisely, considering the scaling functions and the critical exponents calculated, we conclude that the system undergoes a crossover among distinct universality classes. For α≤3 the critical behavior is described by mean-field exponents, while for α≥4 it belongs to the Ising universality class. Finally, in the region where the crossover occurs, 3<α<4, the critical exponents are dependent on α.

A numerical study of planar arrays of correlated spin islands

We present our analysis of a system of interacting islands of XY spins on a triangular lattice that has been introduced a few years ago by Eley et al. to account for the phenomenology in experiments on tunable arrays of proximity coupled long superconductor-normal metal-superconductor junctions. The main features of the model are the separation of a local and a global interaction energy scale and the mesoscopic character of the spin islands. Upon lowering the temperature the model undergoes two crossovers corresponding to an increasing phase coherence on a single island and to the onset of global coherence across the array; the latter is a thermodynamical phase transition in the Ising universality class. The dependence of the second transition on the island edge-to-edge spacing is related to the proximity-effect of the coupling constant.

Roberge-Weiss endpoint at the physical point ofNf=2+1QCD

We study the phase diagram of $N_f = 2+1$ QCD in the $T - \mu_B$ plane and investigate the critical point corresponding to the onset of the Roberge-Weiss transition, which is found for imaginary values of $\mu_B$. We make use of stout improved staggered fermions and of the tree level Symanzik gauge action, and explore four different sets of lattice spacings, corresponding to $N_t = 4,6,8,10$, and different spatial sizes, in order to assess the universality class of the critical point. The continuum extrapolated value of the endpoint temperature is found to be $T_{\rm RW} = 208(5)$ MeV, i.e. $T_{\rm RW}/T_c \sim 1.34(7)$, where $T_c$ is the chiral pseudocritical temperature at zero chemical potential, while our finite size scaling analysis, performed on $N_t = 4$ and $N_t = 6$ lattices, provides evidence for a critical point in the $3d$ Ising universality class.

Nonequilibrium and nonhomogeneous phenomena around a first-order quantum phase transition

We consider non-equilibrium phenomena in a very simple model that displays a zero-temperature first-order phase transition. The quantum Ising model with a four-spin exchange is adopted as a general representative of first-order quantum phase transitions that belong to the Ising universality class, such as for instance the order-disorder ferroelectric transitions, and possibly first-order T = 0 Mott transitions. In particular, we address quantum quenches in the exactly solvable limit of infinite connectivity and show that, within the coexistence region around the transition, the system can remain trapped in a metastable phase, as long as it is spatially homogeneous so that nucleation can be ignored. Motivated by the physics of nucleation, we then study in the same model static but inhomogeneous phenomena that take place at surfaces and interfaces. The first order nature implies that both phases remain locally stable across the transition, and with that the possibility of a metastable wetting layer showing up at the surface of the stable phase, even at T=0. We use mean-field theory plus quantum fluctuations in the harmonic approximation to study quantum surface wetting.