Strong universality and algebraic scaling in two-dimensional Ising spin glasses

At zero temperature, two-dimensional Ising spin glasses are known to fall into several universality classes. Here we consider the scaling at low but nonzero temperatures and provide numerical evidence that eta approximately equal 0 and nu approximately equal 3.5 in all cases, suggesting a unique universality class. This algebraic (as opposed to exponential) scaling holds, in particular, for the +/- J model, with or without dilutions, and for the plaquette diluted model. Such a picture, associated with an exceptional behavior at T = 0, is consistent with a real space renormalization group approach. We also explain how the scaling of the specific heat is compatible with the hyperscaling prediction.     

Nonequilibrium relaxation and critical aging for driven Ising lattice gases

We employ Monte?Carlo simulations to study the nonequilibrium relaxation of driven Ising lattice gases in two dimensions. Whereas the temporal scaling of the density autocorrelation function in the nonequilibrium steady state does not allow a precise measurement of the critical exponents, these can be accurately determined from the aging scaling of the two-time autocorrelations and the order parameter evolution following a quench to the critical point. We obtain excellent agreement with renormalization group predictions based on the standard Langevin representation of driven Ising lattice gases.     

Absence of an Almeida-Thouless line in three-dimensional spin glasses

We present results of Monte Carlo simulations of the three-dimensional Edwards-Anderson Ising spin glass in the presence of a (random) field. A finite-size scaling analysis of the correlation length shows no indication of a transition, in contrast with the zero-field case. This suggests that there is no Almeida-Thouless line for short-range Ising spin glasses.     

Performance limitations of flat-histogram methods

We determine the optimal scaling of local-update flat-histogram methods with system size by using a perfect flat-histogram scheme based upon the exact density of states of 2D Ising models. The typical tunneling time needed to sample the entire bandwidth does not scale with the number of spins N as the minimal N2 of an unbiased random walk in energy space. While the scaling is power law for the ferromagnetic and fully frustrated Ising model, for the +/-J nearest-neighbor spin glass the distribution of tunneling times is governed by a fat-tailed Fréchet extremal value distribution that obeys exponential scaling. Furthermore, the shape parameters of these distributions indicate that statistical sample means become ill defined already for moderate system sizes within these complex energy landscapes.     

Random-bond Ising chain in a transverse magnetic field: A finite-size scaling analysis

We investigate the zero-temperature quantum phase transition of the randombond Ising chain in a transverse magnetic field. Its critical properties are identical to those of the McCoy-Wu model, which is a classical Ising model in two dimensions with layered disorder. The latter is studied via Monte Carlo simulations and transfer matrix calculations and the critical exponents are determined with a finite-size scaling analysis. The magnetization and susceptibility obey conventional rather than activated scaling. We observe that the order parameter and correlation function probability distribution show a nontrivial scaling near the critical point, which implies a hierarchy of critical exponents associated with the critical behavior of the generalized correlation lengths.     

Improved phenomenological renormalization schemes

An analysis is made of various methods of phenomenological renormalization based on finite-dimensional scaling equations for inverse correlation lengths, the singular part of the free energy density, and their derivatives. The analysis is made using two-dimensional Ising and Potts lattices and the three-dimensional Ising model. Variants of equations for the phenomenological renormalization group are obtained which ensure more rapid convergence than the conventionally used Nightingale phenomenological renormalization scheme. An estimate is obtained for the critical finite-dimensional scaling amplitude of the internal energy in the three-dimensional Ising model. It is shown that the two-dimensional Ising and Potts models contain no finite-dimensional corrections to the internal energy so that the positions of the critical points for these models can be determined exactly from solutions for strips of finite width. It is also found that for the two-dimensional Ising model the scaling finite-dimensional equation for the derivative of the inverse correlation length with respect to temperature gives the exact value of the thermal critical index.     

Interface Motion in Disordered Ferromagnets

We consider numerically the depinning transition in the random-field Ising model. Our analysis reveals that the three and four dimensional model displays a simple scaling behavior whereas the five dimensional scaling behavior is affected by logarithmic corrections. This suggests that d = 5 is the upper critical dimension of the depinning transition in the random-field Ising model. Furthermore, we investigate the so-called creep regime (small driving fields and temperatures) where the interface velocity is given by an Arrhenius law.     

Long-term fluctuations in globally coupled phase oscillators with general coupling: finite size effects

We investigate the diffusion coefficient of the time integral of the Kuramoto order parameter in globally coupled nonidentical phase oscillators. This coefficient represents the deviation of the time integral of the order parameter from its mean value on the sample average. In other words, this coefficient characterizes long-term fluctuations of the order parameter. For a system of N coupled oscillators, we introduce a statistical quantity D, which denotes the product of N and the diffusion coefficient. We study the scaling law of D with respect to the system size N. In other well-known models such as the Ising model, the scaling property of D is D～O(1) for both coherent and incoherent regimes except for the transition point. In contrast, in the globally coupled phase oscillators, the scaling law of D is different for the coherent and incoherent regimes: D～O(1/N(a)) with a certain constant a>0 in the coherent regime and D～O(1) in the incoherent regime. We demonstrate that these scaling laws hold for several representative coupling schemes.     

Wang-Landau study of the 3D Ising model with bond disorder

We implement a two-stage approach of the Wang-Landau algorithm to investigate the critical properties of the 3D Ising model with quenched bond randomness. In particular, we consider the case where disorder couples to the nearest-neighbor ferromagnetic interaction, in terms of a bimodal distribution of strong versus weak bonds. Our simulations are carried out for large ensembles of disorder realizations and lattices with linear sizes L in the range L=8-64L=8{-}64. We apply well-established finite-size scaling techniques and concepts from the scaling theory of disordered systems to describe the nature of the phase transition of the disordered model, departing gradually from the fixed point of the pure system. Our analysis (based on the determination of the critical exponents) shows that the 3D random-bond Ising model belongs to the same universality class with the site- and bond-dilution models, providing a single universality class for the 3D Ising model with these three types of quenched uncorrelated disorder.     

Scaling properties of geometric parallelization

We present a universal scaling law for all geometrically parallelized computer simulation algorithms. For algorithms with local interaction laws we calculate the scaling exponents for zero and infinite lattice size. The scaling is tested on local (cellular automata, Metropolis Ising) as well as cluster (Swendsen-Wang) algorithms. The practical aspects of the scaling properties lead to a simple recipe for finding the optimum number of processors to be used for the parallel simulation of a particular system.     

Pair function for the rectangular Ising ferromagnet

A representation of the pair correlation function for the rectangular Ising model in zero magnetic field is derived using a new spinor technique; this enables the scaling limit to be established, as well as several analytical properties of the scaling functions.Supported by the Fonds National Suisse de la Recherche Scientifique, by the National Science Foundation Grant No. PHY76-17191, and by the National Research Council of Canada Grant No. NRC A9344     

Scaling relationship between effective critical exponents throughout the crossover region in thin Ising films

The Monte Carlo (MC) approach is used to check the validity of the scaling relationship for the effective critical exponents in thin Ising films. We investigate this relationship not just in the critical region but throughout the crossover to the expected two-dimensional behavior. Our results indicate that this scaling relationship is very well-fulfilled throughout the entire crossover temperature region, as predicted by a previous renormalization group analysis. The two-dimensional universality class of Ising films is confirmed by means of data collapsing plots for plates with increasing L, up to L=100. The evolution of the maximum value of the effective critical exponents with film thickness is discussed. Received 22 April 1999     

Nagel scaling, relaxation, and universality in the kinetic ising model on an alternating isotopic chain

The dynamic critical exponent and the frequency and wave-vector dependent susceptibility of the kinetic Ising model on an alternating isotopic chain with Glauber dynamics are examined. The analysis provides a connection between a microscopic model and the Nagel scaling curve originally proposed to describe dielectric susceptibility measurements of several glass-forming liquids. While support is given to the hypothesis relating the Nagel scaling to multiple relaxation processes, it is also found that the scaling function may exhibit plateau regions and does not hold for all temperatures.     

On the validity of scaling theory for the anisotropic ising model

When scaling theory is applied to successive terms of a perturbation expansion for the anisotropic Ising model, a consistent gap index of 1.75 is predicted. We have recently derived a series expansion for this model and present here an analysis of the series which suggests that the scaling theory results are not obeyed.     

Finite-size scaling for the ising model on the Möbius strip and the klein bottle

We study the finite-size scaling properties of the Ising model on the M?bius strip and the Klein bottle. The results are compared with those of the Ising model under different boundary conditions, that is, the free, cylindrical, and toroidal boundary conditions. The difference in the magnetization distribution function p(m) for various boundary conditions is discussed in terms of the number of the percolating clusters and the cluster size. We also find interesting aspect-ratio dependence of the value of the Binder parameter at T = T(c) for various boundary conditions. We discuss the relation to the finite-size correction calculations for the dimer statistics.     

Quantum Ising Criticality in Dimerized XY Spin-1/2 Chain

In the present paper we propose a spin-1/2 chain which provides an exactly solvable example to studythe Ising criticality with the central charge c = 1/2. We performthe diagonalization of this model in the presence ofmagnetic field. From the full energy spectrum, the central charge and the scaling dimensions are given at the criticalpoint. The results show evidently that the quantum Ising criticality exists in such a system.     

Biorthogonal quantum criticality in non-Hermitian many-body systems

We develop the perturbation theory of the fidelity susceptibility in biorthogonal bases for arbitrary interacting non-Hermitian many-body systems with real eigenvalues. The quantum criticality in the non-Hermitian transverse field Ising chain is investigated by the second derivative of the ground-state energy and the ground-state fidelity susceptibility. We show that the system undergoes a second-order phase transition with the Ising universal class by numerically computing the critical points and the critical exponents from the finite-size scaling theory. Interestingly, our results indicate that the biorthogonal quantum phase transitions are described by the biorthogonal fidelity susceptibility instead of the conventional fidelity susceptibility.     

Computer simulation of the critical behavior of 3D disordered ising model

The critical behavior of the disordered ferromagnetic Ising model is studied numerically by the Monte Carlo method in a wide range of variation of concentration of nonmagnetic impurity atoms. The temperature dependences of correlation length and magnetic susceptibility are determined for samples with various spin concentrations and various linear sizes. The finite-size scaling technique is used for obtaining scaling functions for these quantities, which exhibit a universal behavior in the critical region; the critical temperatures and static critical exponents are also determined using scaling corrections. On the basis of variation of the scaling functions and values of critical exponents upon a change in the concentration, the conclusion is drawn concerning the existence of two universal classes of the critical behavior of the diluted Ising model with different characteristics for weakly and strongly disordered systems.     

Finite-size scaling for near-critical continuum fluids at constant pressure

We consider the application of finite-size scaling methods to isothermal-isobaric (constant-NpT) simulations of pure continuum fluids. A finite-size scaling ansatz is made for the dependence of the relevant scaling operators on the particle number. To test the proposed scaling form, constant pressure simulations of the Lennard-Jones fluid at its liquid-vapour critical point are performed. The critical scaling operator distributions are obtained and their scaling with particle number is found to be consistent with the proposed behaviour. The forms of these scaling distributions are shown to be identical to their Ising model counterparts. The relative merits of employing the constant-NpT and grand canonical (constant-μVT) ensembles for simulations of fluid critically are also discussed.