A New Method to Calculate the Spin-Glass Order Parameter of the Two-Dimensional +/-J Ising Model

A new method to numerically calculate the $n$th moment of the spin overlap of the two-dimensional $\pm J$ Ising model is developed using the identity derived by one of the authors (HK) several years ago. By using the method, the $n$th moment of the spin overlap can be calculated as a simple average of the $n$th moment of the total spins with a modified bond probability distribution. The values of the Binder parameter etc have been extensively calculated with the linear size, $L$, up to L=23. The accuracy of the calculations in the present method is similar to that in the conventional transfer matrix method with about $10^{5}$ bond samples. The simple scaling plots of the Binder parameter and the spin-glass susceptibility indicate the existence of a finite-temperature spin-glass phase transition. We find, however, that the estimation of $T_{\rm c}$ is strongly affected by the corrections to scaling within the present data ($L\leq 23$). Thus, there still remains the possibility that $T_{\rm c}=0$, contrary to the recent results which suggest the existence of a finite-temperature spin-glass phase transition.Comment: 10 pages,8 figures: final version to appear in J. Phys.

Criticality in the two-dimensional random-bond Ising model

The two-dimensional (2D) random-bond Ising model has a novel multicritical point on the ferromagnetic to paramagnetic phase boundary. This random phase transition is one of the simplest examples of a 2D critical point occurring at both finite temperatures and disorder strength. We study the associated critical properties, by mapping the random 2D Ising model onto a network model. The model closely resembles network models of quantum Hall plateau transitions, but has different symmetries. Numerical transfer matrix calculations enable us to obtain estimates for the critical exponents at the random Ising phase transition. The values are consistent with recent estimates obtained from high-temperature series.Comment: minor changes, 7 pages LaTex, 8 postscript figures included using epsf; to be published Phys. Rev. B 55 (1997

Ground states with cluster structures in a frustrated Heisenberg chain

We examine the ground state of a Heisenberg model with arbitrary spin S on a one-dimensional lattice composed of diamond-shaped units. A unit includes two types of antiferromagnetic exchange interactions which frustrate each other. The system undergoes phase changes when the ratio $\lambda$ between the exchange parameters varies. In some phases, strong frustration leads to larger local structures or clusters of spins than a dimer. We prove for arbitrary S that there exists a phase with four-spin cluster states, which was previously found numerically for a special value of $\lambda$ in the S=1/2 case. For S=1/2 we show that there are three ground state phases and determine their boundaries.Comment: 4 pages, uses revtex.sty, 2 figures available on request from [email protected], to be published in J. Phys.: Cond. Mat

Shape Effects of Finite-Size Scaling Functions for Anisotropic Three-Dimensional Ising Models

The finite-size scaling functions for anisotropic three-dimensional Ising models of size $L_1 \times L_1 \times aL_1$ ($a$: anisotropy parameter) are studied by Monte Carlo simulations. We study the $a$ dependence of finite-size scaling functions of the Binder parameter $g$ and the magnetization distribution function $p(m)$. We have shown that the finite-size scaling functions for $p(m)$ at the critical temperature change from a two-peak structure to a single-peak one by increasing or decreasing $a$ from 1. We also study the finite-size scaling near the critical temperature of the layered square-lattice Ising model, when the systems have a large two-dimensional anisotropy. We have found the three-dimensional and two-dimensional finite-size scaling behavior depending on the parameter which is fixed; a unified view of 3D and 2D finite-size scaling behavior has been obtained for the anisotropic 3D Ising models.Comment: 6 pages including 11 eps figures, RevTeX, to appear in J. Phys.

High Temperature Expansion Study of the Nishimori multicritical Point in Two and Four Dimensions

We study the two and four dimensional Nishimori multicritical point via high temperature expansions for the $\pm J$ distribution, random-bond, Ising model. In $2d$ we estimate the the critical exponents along the Nishimori line to be $\gamma=2.37\pm 0.05$, $\nu=1.32\pm 0.08$. These, and earlier $3d$ estimates $\gamma =1.80\pm 0.15$, $\nu=0.85\pm 0.08$ are remarkably close to the critical exponents for percolation, which are known to be $\gamma=43/18$, $\nu=4/3$ in $d=2$ and $\gamma=1.805\pm0.02$ and $\nu=0.875\pm 0.008$ in $d=3$. However, the estimated $4d$ Nishimori exponents $\gamma=1.80\pm 0.15$, $\nu=1.0\pm 0.1$, are quite distinct from the $4d$ percolation results $\gamma=1.435\pm 0.015$, $\nu=0.678\pm 0.05$.Comment: 5 pages, RevTex, 3 postscript files; To appear in Physical Review

Ground-State Phase Diagram of the Two-Dimensional Quantum Heisenberg Mattis Model

The two-dimensional $S=1/2$ asymmetric Heisenberg Mattis model is investigated with the exact diagonalization of finite clusters. The N\'eel order parameter and the spin glass order parameter can be smoothly extrapolated to the thermodynamic limit in the antiferromagnetic region, as in the pure Heisenberg antiferromagnet. The critical concentration of the N\'eel phase is consistent with that of the two-dimensional Ising Mattis model, and the spin glass order parameter increases monotonously as the ferro-bond concentration increases. These facts suggest that quantum fluctuation does not play an essential role in two-dimensional non-frustrated random spin systems. KEYWORDS: quantum spin system, ground state, randomness, Mattis model, N\'eel order, spin glass orderComment: 10 pages, LaTeX, 6 compressed/uuencoded postscript figures, J. Phys. Soc. Jpn. 65 (1996) No. 2 in pres

Aging Relation for Ising Spin Glasses

We derive a rigorous dynamical relation on aging phenomena -- the aging relation -- for Ising spin glasses using the method of gauge transformation. The waiting-time dependence of the auto-correlation function in the zero-field-cooling process is equivalent with that in the field-quenching process. There is no aging on the Nishimori line; this reveals arguments for dynamical properties of the Griffiths phase and the mixed phase. The present method can be applied to other gauge-symmetric models such as the XY gauge glass.Comment: 9 pages, RevTeX, 2 postscript figure

Vicinal Surface with Langmuir Adsorption: A Decorated Restricted Solid-on-solid Model

We study the vicinal surface of the restricted solid-on-solid model coupled with the Langmuir adsorbates which we regard as two-dimensional lattice gas without lateral interaction. The effect of the vapor pressure of the adsorbates in the environmental phase is taken into consideration through the chemical potential. We calculate the surface free energy $f$, the adsorption coverage $\Theta$, the step tension $\gamma$, and the step stiffness $\tilde{\gamma}$ by the transfer matrix method combined with the density-matrix algorithm. Detailed step-density-dependence of $f$ and $\Theta$ is obtained. We draw the roughening transition curve in the plane of the temperature and the chemical potential of adsorbates. We find the multi-reentrant roughening transition accompanying the inverse roughening phenomena. We also find quasi-reentrant behavior in the step tension.Comment: 7 pages, 12 figures (png format), RevTeX 3.1, submitted to Phys. Rev.