By considering the quench dynamics of two-dimensional frustrated Ising models through numerical simulations, we investigate the dynamical critical behavior on the multicritical Nishimori point (NP). We calculate several dynamical critical exponents, namely, the relaxation exponent z_{c}, the autocorrelation exponent λ_{c}, and the persistence exponent θ_{c}, after a quench from the high temperature phase to the NP. We confirm their universality with respect to the lattice geometry and bond distribution. For a quench from a power-law correlated initial state to the NP, the aging dynamics are much slower. We also look up the issue of multifractality during the critical dynamics by investigating different moments of the spatial correlation function. We observe a single growth law for all the length scales extracted from different moments, indicating that the equilibrium multifractality at the NP does not affect the dynamics.
By considering the quench dynamics of two-dimensional frustrated Ising models through numerical simulations, we investigate the dynamical critical behavior on the multicritical Nishimori point (NP). We calculate several dynamical critical exponents, namely, the relaxation exponent z_{c}, the autocorrelation exponent λ_{c}, and the persistence exponent θ_{c}, after a quench from the high temperature phase to the NP. We confirm their universality with respect to the lattice geometry and bond distribution. For a quench from a power-law correlated initial state to the NP, the aging dynamics are much slower. We also look up the issue of multifractality during the critical dynamics by investigating different moments of the spatial correlation function. We observe a single growth law for all the length scales extracted from different moments, indicating that the equilibrium multifractality at the NP does not affect the dynamics.
The ±J Ising model is a simple frustrated spin model, where the exchange couplings independently take the discrete value -J with probability p and +J with probability 1-p. It is especially appealing due to its connection to quantum error correcting codes. Here, we investigate the nonequilibrium critical behavior of the two-dimensional ±J Ising model, after a quench from different initial conditions to a critical point T_{c}(p) on the paramagnetic-ferromagnetic (PF) transition line, especially above, below, and at the multicritical Nishimori point (NP). The dynamical critical exponent z_{c} seems to exhibit nonuniversal behavior for quenches above and below the NP, which is identified as a preasymptotic feature due to the repulsive fixed point at the NP, whereas for a quench directly to the NP, the dynamics reaches the asymptotic regime with z_{c}≃6.02(6). We also consider the geometrical spin clusters (of like spin signs) during the critical dynamics. Each universality class on the PF line is uniquely characterized by the stochastic Loewner evolution with corresponding parameter κ. Moreover, for the critical quenches from the paramagnetic phase, the model, irrespective of the frustration, exhibits an emergent critical percolation topology at the large length scales.
The random-bond Ising model on the square lattice has several disordered critical points, depending on the probability distribution of the bonds. There are a finite-temperature multicritical point, called the Nishimori point, and a zero-temperature fixed point, for both a binary distribution where the coupling constants take the values ± J and a Gaussian disorder distribution. Inclusion of dilution in the ± J distribution ( J = 0 for some bonds) gives rise to another zero-temperature fixed point which can be identified with percolation in the non-frustrated case ( J ≥ 0). We study these fixed points using numerical (transfer matrix) methods. We determine the location, critical exponents, and central charge of the different fixed points and study the spin–spin correlation functions. Our main findings are the following. (1) We confirm that the Nishimori point is universal with respect to the type of disorder, i.e. we obtain the same central charge and critical exponents for the ± J and Gaussian distributions of disorder. (2) The Nishimori point, the zero-temperature fixed point for the ± J and Gaussian distributions of disorder, and the percolation point in the diluted case all belong to mutually distinct universality classes. (3) The paramagnetic phase is re-entrant below the Nishimori point, i.e. the zero-temperature fixed points are not located exactly below the Nishimori point, either for the ± J distribution, or for the Gaussian distribution.
We investigate the critical behavior of the random-bond $\ifmmode\pm\else\textpm\fi{}J$ Ising model on a square lattice at the multicritical Nishimori point in the $T\text{\ensuremath{-}}p$ phase diagram, where T is the temperature and p is the disorder parameter (p=1 corresponds to the pure Ising model). We perform a finite-size scaling analysis of high-statistics Monte Carlo simulations along the Nishimori line defined by $2p\ensuremath{-}1=\text{tanh}(1/T)$, along which the multicritical point lies. The multicritical Nishimori point is located at ${p}^{\ensuremath{\ast}}=0.890\text{ }81(7)$, ${T}^{\ensuremath{\ast}}=0.9528(4)$, and the renormalization-group dimensions of the operators that control the multicritical behavior are y1=0.655(15) and y2=0.250(2); they correspond to the thermal exponent $\ensuremath{\nu}\ensuremath{\equiv}1/{y}_{2}=4.00(3)$ and to the crossover exponent $\ensuremath{\phi}\ensuremath{\equiv}{y}_{1}/{y}_{2}=2.62(6)$.
We consider the three-dimensional $\ifmmode\pm\else\textpm\fi{}J$ model defined on a simple cubic lattice and study its behavior close to the multicritical Nishimori point, where the paramagnetic-ferromagnetic, the paramagnetic-glassy, and the ferromagnetic-glassy transition lines meet in the $T\text{\ensuremath{-}}p$ phase diagram (p characterizes the disorder distribution and gives the fraction of ferromagnetic bonds). For this purpose, we perform Monte Carlo simulations on cubic lattices of size $L\ensuremath{\leqslant}32$ and a finite-size-scaling analysis of the numerical results. The magnetic-glassy multicritical point is found at p∗=0.76820(4), along the Nishimori line given by $2p\ensuremath{-}1=\mathrm{tanh}(J∕T)$. We determine the renormalization-group dimensions of the operators that control the renormalization-group flow close to the multicritical point, y1=1.02(5), y2=0.61(2), and the susceptibility exponent $\ensuremath{\eta}=\ensuremath{-}0.114(3)$. The temperature and crossover exponents are $\ensuremath{\nu}=1∕{y}_{2}=1.64(5)$ and $\ensuremath{\phi}={y}_{1}∕{y}_{2}=1.67(10)$, respectively. We also investigate the model-A dynamics, obtaining the dynamic critical exponent z=5.0(5).
We study the dynamical scaling properties of the two-dimensional fully frustrated Ising model by simulated annealing. The model has a highly degenerate, critical ground state with well known equilibrium properties. To study its out-of-equilibrium behavior, we develop an annealing protocol that is convenient for extending the Kibble-Zurek mechanism to the case of an exponentially divergent correlation length in the limit of the temperature T→Tc=0. We find excellent agreement with the scaling forms derived, including a dependence on the initial temperature that is not present in the conventional Kibble-Zurek mechanism. The dynamic exponent z is consistent among all quantities studied; fits of its dependence on the initial temperature give z=1.99(2) (susceptibility), z=2.072(8) (excess energy), and z=2.09(3) (Edwards-Anderson order parameter). These values are close to z=2, which is consistent with previous studies of relaxation after instantaneous quenches and of time dependent correlations in equilibrium, including the estimate z=2.004(9) from equilibrium autocorrelation functions.
We explore the critical properties of the recently discovered finite-time dynamical phase transition in the nonequilibrium relaxation of Ising magnets after a temperature quench. The transition is characterized by a sudden switch in the relaxation dynamics and it occurs at a sharp critical time. While previous works have focused either on mean-field interactions or on investigating the critical time, we analyze the critical fluctuations at the phase transition in the nearest-neighbor Ising model on a square lattice using Monte Carlo simulations. By means of a finite-size scaling analysis, we extract the critical exponents for the transition. In two spatial dimensions, we find that the exponents are consistent with those of the two-dimensional Ising universality class when the system is initially in the vicinity of the critical point. For initial temperatures below the critical one, however, the critical exponents differ from the Ising-exponents, indicating a distinct dynamical critical phenomenon.