Computer investigations of the three-dimensional Ising model

We have been studying the three-dimensional Ising model using some finite-size scaling ideas. The simulation is done by a fast microcanonical method. Here we present our results for the critical exponents and.     

Interface motion in a two-dimensional Ising model with a field

We determine by Monte Carlo simulations the width of an interface between the stable phase and the metastable phase in a two-dimensional Ising model with a magnetic field, in the case of nonconversed order parameter (Glauber dynamics). At zero temperature, the width increases ast with–1/3, as predicted by earlier theories. As temperature increases, the value of the effective exponent that we measure decreases toward the value 1/4, which is the value in the absence of magnetic field.     

Anisotropy of the interface tension of the three-dimensional Ising model

We determine the interface tension for the 100, 110 and 111 interface of the simple cubic Ising model with nearest-neighbour interaction using novel simulation methods. To overcome the droplet/strip transition and the droplet nucleation barrier we use a newly developed combination of the multimagnetic algorithm with the parallel tempering method. We investigate a large range of inverse temperatures to study the anisotropy of the interface tension in detail.     

Series and Monte Carlo study of high-dimensional Ising models

Ising models in high dimensions are used to compare high-temperature series expansions with Monte Carlo simulations. Simulations of the magnetization on four-, six-, and seven-dimensional hypercubic lattices give consistent values of the critical temperature from both equilibrium and nonequilibrium data ford=6 and 7. We tabulate 15 terms of series expansions for the susceptibility for generald and giveJ/k _B T _c =0.092295 (3) and 0.077706 (2) ford=6 and 7. In contrast to five dimensions, where earlier series found nonanalytic scaling corrections, for d=6 and 7 the leading scaling correction may be analytic inT-T _c . In most cases these expansions gave more accurate results than these simulations.     

Fixed point behavior of cumulants in the three-dimensional Ising universality class

High-order cumulants and factorial cumulants of conserved charges are suggested for the study of the critical dynamics in heavy-ion collision experiments. In this paper, using the parametric representation of the three-dimensional Ising model which is believed to belong to the same universality class as quantum chromo-dynamics, the temperature dependence of the second- to fourth-order (factorial) cumulants of the order parameter is studied. It is found that the values of the normalized cumulants are independent of the external magnetic field at the critical temperature, which results in a fixed point in the temperature dependence of the normalized cumulants. In finite-size systems simulated using the Monte Carlo method, this fixed point behavior still exists at temperatures near the critical. This fixed point behavior has also appeared in the temperature dependence of normalized factorial cumulants from at least the fourth order. With a mapping from the Ising model to QCD, the fixed point behavior is also found in the energy dependence of the normalized cumulants (or fourth-order factorial cumulants) along different freeze-out curves.     

Phase transitions and universality in nonequilibrium steady states of stochastic Ising models

We present results of direct computer simulations and of Monte Carlo renormalization group (MCRG) studies of the nonequilibrium steady states of a spin system with competing dynamics and of the voter model. The MCRG method, previously used only for equilibrium systems, appears to give useful information also for these nonequilibrium systems. The critical exponents are found to be of Ising type for the competing dynamics model at its second-order phase transitions, and of mean-field type for the voter model (consistent with known results for the latter).     

Surface tension from finite-volume vacuum tunneling in the 3D Ising model

We measure the surface tension in the broken phase of the 3D Ising model at a temperatureT=0.955T _c with two different methods which are taken from quantum field theory in finite volumes. Both methods rely on finite-size effects close to the phase transition. The first one measures from the size dependence of the vacuum tunneling energy, which is determined by the decay of a correlation, giving=0.030. The second one extracts from the size dependence of the rate of flip events and its corresponding correlation time. It leads to=0.027. Both values agree reasonably with other calculations.     

Machine learning phase transitions of the three-dimensional Ising universality class

Exploration of the QCD phase diagram and critical point is one of the main goals in current relativistic heavy-ion collisions. The QCD critical point is expected to belong to a three-dimensional (3D) Ising universality class. Machine learning techniques are found to be powerful in distinguishing different phases of matter and provide a new way to study the phase diagram. We investigate phase transitions in the 3D cubic Ising model using supervised learning methods. It is found that a 3D convolutional neural network can be trained to effectively predict physical quantities in different spin configurations. With a uniform neural network architecture, it can encode phases of matter and identify both second- and first-order phase transitions. The important features that discriminate different phases in the classification processes are investigated. These findings can help study and understand QCD phase transitions in relativistic heavy-ion collisions.     

Critical dynamics of the three-dimensional Ising model: A Monte Carlo study

Extensive Monte Carlo simulations have been performed to analyze the dynamical behavior of the three-dimensional Ising model with local dynamics. We have studied the equilibrium correlation functions and the power spectral densities of odd and even observables. The exponential relaxation times have been calculated in the asymptotic one-exponential time region. We find that the critical exponentz=2.09 ±0.02 characterizes the algebraic divergence with lattice size for all observables. The influence of scaling corrections has been analyzed. We have determined integrated relaxation times as well. Their dynamical exponentz _int agrees withz for correlations of the magnetization and its absolute value, but it is different for energy correlations. We have applied a scaling method to analyze the behavior of the correlation functions. This method verifies excellent scaling behavior and yields a dynamical exponentz _scal which perfectly agrees withz.     

Critical wetting in the square Ising model with a boundary field

The Ising square lattice with nearest-neighbor exchangeJ>0 and a free surface at which a boundary magnetic fieldH ₁ acts has a second-order wetting transition. We study the surface excess magnetization and the susceptibility ofL×M lattices by Monte Carlo simulation and probe the critical behavior of this wetting transition, applying finite-size scaling methods. For the cases studied, the results are not consistent with the presumably exactly known values of the critical exponents, because the asymptotic critical region has not yet been reached. Implication of our results for critical wetting in three dimensions and for the application of the present model to adsorbed wetting layers at surface steps are briefly discussed.Alexander von Humboldt-Fellow     

The roughening transition of the three-dimensional Ising interface: A Monte Carlo study

We study the roughening transition of an interface in an Ising system on a 3D simple cubic lattice using a finite-size scaling method. The particular method has recently been proposed and successfully tested for various solid-on-solid models. The basic idea is the matching of the renormalization-groupflow of the interface with that of the exactly solvable body-centered cubic solid-on-solid model. We unambiguously confirm the Kosterlitz-Thouless nature of the roughening transition of the Ising interface. Our result for the inverse transition temperatureK _r=0.40754(5) is almost two orders of magnitude more accurate than the estimate of Mon, Landau, and Stauffer.     

Fast vectorized algorithm for the Monte Carlo simulation of the random field Ising model

An algorithm for the simulation of the 3-dimensional random field Ising model with a binary distribution of the random fields is presented. It uses multi-spin coding and simulates 64 physically different systems simultaneously. On one processor of a Cray YMP it reaches a speed of 184 million spin updates per second. For smaller field strength we present a version of the algorithm that can perform 242 million spin updates per second on the same machine.     

Finite-size effects in surface tension. I. Fluctuating interfaces

We consider a two-dimensional Ising cylinder of circumferenceM and heightN, with a floating interface introduced by the appropriate boundary conditions. An exact analysis of the finite-size effects in surface tension is given and the scaling function for all temperatures is calculated. The results are compared with the Monte Carlo data of Mon and Jasnow.On leave from: Department of Theoretical Chemistry, Oxford University, Oxford, OX1 3UB, England.     

Complete wetting in the three-dimensional transverse Ising model

We consider a three-dimensional Ising model in a transverse magnetic fieldh and a bulk fieldH. An interface is introduced by an appropriate choice of boundary conditions. At the point (H=0,h=0) spin configurations corresponding to different positions of the interface are degenerate. By studying the phase diagram near this multiphase point using quantum mechanical perturbation theory, we show that the quantum fluctuations, controlled byh, split the multiphase degeneracy giving rise to an infinite sequence of layering transitions.     

二维伊辛模型的畴壁动力学

表面可以改变纳米磁性薄膜的结构和相变温度,畴壁动力学由此成为研究的重点.本文采用动力学蒙特卡罗模拟方法,对二维Ising模型磁畴界面的非平衡动力学展开数值研究.系统初态设为半正半负,即由完全有序但自旋取向完全相反的两部分组成,其间的磁畴壁随时间生长.通过对磁化标度形式的分析,发现畴壁内外的动力学标度形式虽然相同,但临界...     

二维伊辛模型的畴壁动力学

表面可以改变纳米磁性薄膜的结构和相变温度,畴壁动力学由此成为研究的重点。本文采用动力学蒙特卡罗模拟方法,对二维Ising模型磁畴界面的非平衡动力学展开数值研究。系统初态设为半正半负,即由完全有序但自旋取向完全相反的两部分组成,其间的磁畴壁随时间生长。通过对磁化标度形式的分析,发现畴壁内外的动力学标度形式虽然相同,但临界指数在数值上却存在很大差异,相差一个 =1,这是由初始条件导致的。     

Surface Critical Behavior of Two-Dimensional Dilute Ising Models

Ising models with nearest neighbor ferromagnetic random couplings on a square lattice with a (1, 1) surface are studied, using Monte Carlo techniques and a star-triangle transformation method. In particular, the critical exponent of the surface magnetization is found to be close to that of the perfect model, Β₁ = 1/2. The crossover from surface to bulk critical properties is discussed.