Yukiyasu Ozeki

Phase diagram and critical exponents of the ±J ising model in finite dimensions by Monte Carlo renormalization group

We have performed large-scale numerical simulations for the ± J Ising model on the simple cubic and the square lattices with asymmetric weight of ferromagnetic and antiferromagnetic bonds. The Monte Carlo renormalization group method is used to estimate T c , γ/ν, β/ν and ν at various points along the phase boundary between ferromagnetic and paramagnetic phases including the tricritical point. The obtained values of T c and the resulting phase diagrams of the simple cubic and the square lattices are much more accurate than those estimated so far. In the simple cubic lattice, we find evidence for the possibility that universality of the ferromagnetic critical exponents holds along the critical line except at the tricritical point, while weak universality holds even at the tricritical point. Weak universality seems to hold also in the square lattice.

Phase Diagram of the ± J Ising Model in Two Dimensions

The ± J Ising model with asymmetric probability weight of ± bonds in two dimensions is studied in the p - T plane by the numerical transfer-matrix method. We calculate finite-width susceptibilities and finite-width correlation lengths of long strips of various sizes up to 14×10 5 . Finite-size scaling analysis suggests that the random antiphase state exists adjacent to the ferromagnetic phase. The boundary separating these phases is vertical (parallel to the temperature axis) in the phase diagram. The weak universality seems to hold along the ferromagnetic critical line.

Phase diagram of gauge glasses

The authors introduce a general class of random spin systems which are symmetric under local gauge transformations. Their model is a generalization of the usual Ising spin glass and includes the Zq, XY, and SU (2) gauge glasses. For this general class of systems, the internal energy and an upper bound on the specific heat are calculated explicitly in any dimensions on a special line in the phase diagram. Although the line intersects a phase boundary at a multicritical point, the internal energy and the bound on the specific heat are found to be written in terms of a simple function. They also show that the boundary between the ferromagnetic and nonferromagnetic phases is parallel to the temperature axis in the low-temperature region of the phase diagram. This means the absence of re-entrant transitions. All these properties are derived by simple applications of gauge transformations of spin and randomness degrees of freedom.

Ground State Properties of the ± J Ising Model in Two Dimensions

We investigate the ground state properties of the ± J Ising model with asymmetric probability weight in two dimensions. We propose a numerical transfer matrix method which enables us to calculate directly the ground-state properties such as correlation functions, defect energy and entropy at T =0. We find that the spin glass correlation function decays in a power law like \(\tilde{g}(r)\sim{r}^{-\tilde{\eta}}\) with \(\tilde{\eta}\simeq0.21\) in the range \(0.5\leq{p}\lesssim0.85\), suggesting the weak universality of spin glass-the concentration of + J bond is denoted by p . At the ferromagnetic phase boundary locating around p =0.89, the ferromagnetic correlation function decays as \(g(r)\sim{r}^{-\tilde{\eta}}\) with η≃0.240±0.016, which suggests the weak universality of ferromagnetic critical exponent from the non-random critical point ( p =1, T =2.269···) to the ground sate counterpart. The behavior of correlation functions, and also that of defect energies indicate the existence of an extra-phase, c...

Nonequilibrium Relaxation of Fluctuations of Physical Quantities

The nonequilibrium relaxation (NER) process of fluctuations of physical quantities is studied simulationally. It is shown that the NER method, which is a convenient technique for studying the phase and transition, is useful not only for identifying the phase and locating the transition point, but also for estimating both static and dynamical exponents. As an example, the cubic-lattice ferromagnetic Ising model is analyzed. The transition inverse-temperature is estimated to be K c =0.2216595(15). The exponents of correlation length, magnetization, and specific heat are estimated to be ν=0.635(5), β= 0.325(5) and α= 0.14(2), respectively. The dynamical exponent z is estimated to be 2.055(10).

Multicritical dynamics for the ? J Ising Model

Exact dynamical properties are discussed for the ? J Ising model around the multicritical point (MCP). It is found that the relation of relaxations between the non-equilibrium remanent magnetization m(t) and the equilibrium autocorrelation function q(t) at the MCP is different from that at the pure critical point. The dynamic critical exponent for the ferromagnetic ordering defined by and that for the spin glass ordering defined by become identical at the MCP. Accurate numerical calculations for them are performed in two and three dimensions using the non-equilibrium relaxation analysis. The MCP is located at and the exponent is estimated as for the square lattice. They are estimated as and for the simple cubic lattice.

Interfacial approach tod-dimensional ±J Ising models in the neighborhood of the ferromagnetic phase boundary

Two- and three-dimensional ±J Ising models in the neighborhood of the ferromagnetic phase (FP) boundary in the concentration-temperature (p-T) plane are studied, investigating the size dependence of interfacial free energies calculated by a transfer matrix method. Thep andT dependences of two stiffness exponents relevant to the FP and the nonferromagnetic ordered phase lead to the following results in two dimensions, giving a unified view. It is confirmed that the random antiphase state (RAS) exists in contact with the vertical FP boundary. Spatial fluctuations are dominant near the vertical boundary, which is separated by the Nishimori line from the remaining FP boundary governed by thermal fluctuations. The RAS is a kind of Mattis spin glass such that it changes to the FP smoothly with nonsingular physical connectivity, but with a percolation singularity of its ferromagnetic part. Universal finite-size critical amplitudes are consistent with them. Results in three dimensions give only suggestions which are similar to the two-dimensional results. These results suggest important insight into spin-glass properties in higher dimensions.

Ground-State Long-Range Order in the Two-Dimensional XXZ Model

We prove the existence of Neel-type long-range order in the ground state of the spin-1/2 XXZ model with Δ (exchange anisotropy)>1.67 on the square lattice. We further show the existence of long-range order for Δ >1.10 and 0≦ Δ <0.59 by assuming monotonicity of nearest neighbor correlations as functions of the system size. The assumption of monotonicity is supported by numerical calculations.

Nonequilibrium relaxation method

The nonequilibrium relaxation (NER) method is a numerical technique to analyse equilibrium phase transitions. One can estimate the transition point and critical exponents calculating relaxations of order parameter and fluctuations in NER processes from nonequilibrium initial states to the equilibrium ones. It is applied to the second-order transition as well as the first-order one; the precise estimation of transition temperature for the latter case is achieved with the mixed phase initialization technique. Since the equilibration is not necessary in the simulation, one can analyse systems with large sizes which can be recognized as in the thermodynamic limit up to the observation time. This leads to the analysis being more efficient for slowly relaxing systems with frustrations and randomness. It can be extended to the Kosterlitz?Thouless transition and the spin glass transition cases by the use of the finite-time scaling analysis. It is also extended to quantum systems.

Non-Universal Critical Behavior in the Ferromagnetic Transition of the ± J Ising Model

Ferromagnetic (FM) transition of the ± J Ising model is studied for square and simple-cubic lattices by Monte Carlo simulations. Non-equilibrium relaxation (NER) process of magnetization from the all-aligned spin state is analyzed. The transition temperature is estimated at several values of FM bond concentrations to provide the boundary of the FM phase. The estimated NER exponent of magnetization shows non-universal behavior and it becomes larger smoothly when the transition temperature becomes higher.