Bond percolation with parallelism restriction on square lattices

As a simple approximation for the ±J spin glass we studied bond percolation on square lattices. However, two neighboring chains of ferromagnetic bonds are required for spins to be regarded as connected. We determine the percolation thresholdp _c =0.8282±0.0002 and the critical exponent =0.75 _–0.05 ^+0.02 for this specific percolation by means of Monte-Carlo simulation on square lattices (up to 150×150).     

Critical temperature for two-dimensional ising ferromagnets with quenched bond disorder

The configuration-averaged free energy of a quenched, random bond Ising model on a square lattice which contains an equal mixture of two types of ferromagnetic bonds J₁ and J₂ is shown to obey the same duality relation as the ordered rectangular model with the same two bond strengths. If the random.system has a single, sharp critical point, the critical temperature T_c must be identical to that of the ordered system, i.e., sinh(2J ₁/kT _c) sinh(2J ₂/kT _c) = 1. Since _c ^(B) = 1/2, we can takeJ ₂ 0 and use Bergstresser-type inequalities to obtain(/dp) exp(–2J ₁/kT_c¦_p=p_c + = 1, in agreement with Bergstresser's rigorous result for the diluted ferromagnet near the percolation threshold.Work supported in part by National Science Foundation Grant No. DMR 76-21703, Office of Naval Research Grant No. N00014-76-C-0106, and National Science Foundation MRL program Grant No. DMR 76-00678.Paper presented at the 37th Yeshiva University Statistical Mechanics Meeting, May 10, 1977.     

On the ground-state threshold in random two-dimensional Ising ±J models

For square, triangular, and for hexagonal lattices there is numerical and theoretical support that the ground-state thresholdp _c between ferro- and paramagnetism in random 2D Ising ±J models, withp as the concentration of antiferromagnetic bonds, is identical top ^*which is characterized by minimal matching properties of frustrated plaquettes. From square lattices of size 100×100 we have got p_c,sq<0.117 by simulations which produced average groundstate magnetizations per spin by means of exact minimal matchings. Moreover, from the squareL×L-lattices treated (L = 10, 20, 50, 100) we obtained the estimatep _c,sq 0.1 which is in agreement with the Grinstein estimatep _c,sq 0.099 andp _c,sq 0.105 by Freund and Grassberger.     

Statistics of self-avoiding walks on random lattices

The statistics of directed self-avoiding walks (SAWs) on randomly bond diluted square lattices have been solved exactly and a computer simulation study of the statistics of ordinary SAWs on diluted square lattices has also been performed. The configurational averaging has been performed here over the logarithms of the distribution functions. We find that the critical behaviour remains unchanged below a certain dilution concentrationp ^*, dependent on the length of the chains considered (p ^*=0 forN), and a crossover to a higher order critical behaviour occurs beyond that point.     

Exact ground states of two-dimensional ±J Ising spin glasses

In this paper we study the problem of finding an exact ground state of a two-dimensional ±J Ising spin glass on a square lattice with nearest neighbor interactions and periodic boundary conditions when there is a concentrationp of negative bonds, withp ranging between 0.1 and 0.9. With our exact algorithm we can determine ground states of grids of sizes up to 50×50 in a moderate amount of computation time (up to 1 hr each) for several values ofp. For the ground-state energy of an infinite spin-glass system withp=0.5 we estimateE _0.5 =–1.4015±0.0008. We report on extensive computational tests based on more than 22,000 experiments.     

Magnetization at corners in two-dimensional Ising models

We investigate the corner spin magnetization of two-dimensional ferromagnetic Ising models in various wedge geometries. Results are obtained for triangular and square lattices by numerical studies on finite wedges using the star-triangle transformation, as well as analytic calculations using corner transfer matrices and a fermionic representation of the row-to-row transfer matrix. The corner magnetizations vanish at the bulk critical temperature T_c with an exponent _c differing from the bulk exponent. For isotropic systems with free edges we find that _c is given simply by _c =/2, where is the angle at the corner. For apex magnetizations of conical lattices we obtain the strikingly similar result _a=/4. These formulas apply equally well to anisotropic lattices if the angle is interpreted as an effective angle, ^eff, determined by the relative strengths of the interactions.     

Proof of the first order phase transition in the Slater KDP model

The Slater KDP model defined on d-dimensional tetrahedral lattices is proved to have a phase transition for which the entropy and energy are discontinuous functions at a transition temperaturekT _c =/ln2, independent of dimensionality.     

Monte Carlo study of phase diagram for correlated site-bond percolation in two dimensions

At the critical point of the square Ising model, the percolation threshold for randomly active bonds between up spins is close top _Bc =0.60 and seems compatible with the predictionp _Bc =1-exp(–2J/k _B T _c )=0.586 of Coniglio and Klein. Longer simulations on larger lattices are necessary for a more precise clarification.     

Spin-wave analysis of easy-axis quantum antiferromagnets on the triangular lattice

We perform the standard spin wave analysis of the triangular Heisenberg quantum antiferromagnet with nearest neighbour coupling. The exchange interaction is taken to be J _{z=J x=J y} (0<1). We give a simple explanation of the non-trivial classical degeneracy pointed out by Miyashita and Kawamura and show that it is removed by quantum fluctuations, but that the degeneracy manifests itself through the appearance of a second gapless spin-wave branch. The existence of a second gapless mode has a drastic influence on the quasiclassical behaviour near the Ising limit: the energy gain with respect to the Ising state energy is found to be linear in , and the reduction of the sublattice magnetization on two of the three sublattices remains finite as 0. These findings essentially invalidate the original qualitative arguments [14] in favour of a spin-liquid ground state of the anisotropic triangular antiferromagnet.     

Extremity of the disordered phase in the Ising model on the Bethe lattice

We prove that the disordered Gibbs distribution in the ferromagnetic Ising model on the Bethe lattice is extreme forTT _c ^SG , whereT _c ^SG is the critical temperature of the spin glass model on the Bethe lattice, and it is not extreme forT _c ^SG .     

First-order phase transitions in finite-size strips with interfaces

Finite-size rounding of the magnetization discontinuity at the magnetic phase transition atH=0 (T<T _c ) in 2d Ising-type strips of sizeL ×L , with ± boundary conditions alongL inducing an interface of lengthL , is studied by phenomenological considerations and transfer matrix techniques. Scaling expressions are derived forL =O(L ) and also in the infinite strip limitL . Most of the results can be extended to the 3d case.     

A simulation of the mixed spin 3–spin 3/2 ferrimagnetic Ising model

The mixed spin 3–spin 3/2 ferrimagnetic Ising model was simulated using cooling algorithm on cellular automaton (CA). The simulations were carried out in the intervals ?4 ≤ D_A/J ≤ 8 and ?4 ≤ D_B/J ≤ 8 for the square lattices with periodic boundary conditions. The ground-state phase diagram of the model has different types of ferrimagnetic phases. Although only the antiferromagnetic nearest-neighbor interaction was contained in the Hamiltonian, the compensation points emerged through D_A/J = 2 at kT/J = 0. The values of the critical exponents (ν, α , β and γ) were estimated within the framework of the finite-size scaling theory and power-law relations for the selected D_A/J values (?2, 0, 1, 2, and 4). The estimated critical exponent values were in good agreement with the universal values of the two-dimensional Ising model (ν = 1, α = α′ = 0, β = 0.125, β′ = 0.875 and γ = γ′ = 1.75).     

Computer simulation of the nonequilibrium critical behavior of disordered two-dimensional Ising systems

We report the results of a computer simulation of the critical relaxation of the magnetization in the two-dimensional Ising model with nonmagnetic impurity atoms frozen at the lattice sites. We assume a square lattice of dimension 400² with spin concentrationsp=1.0, 0.95, 0.9, 0.85, 0.8, 0.75, 0.7. The Monte Carlo and dynamic renormalization group methods are used to determine the dynamical critical indexz as a function ofp: z(p): z(1)=2.24±0.07,z(0.95)=2.24±0.06,z(0.85)=2.38±0.05,z(0.8)=2.51±0.06,z(0.75)=2.66±0.07,z(0.7)=2.88±0.06. It is shown thatz(p) obeys a singular scaling law of the formz=A | ln (p–p _c) |+B withA=0.56±0.07,B=1.62±0.07.Omsk State University. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 8, pp. 83–88, August, 1994.     

Molecular field effects in Mössbauer spectra of FeSb2O4

Mössbauer effect measurements have been made using⁵⁷Fe in FeSb₂O₄. At liquid helium temperature a combined electric quadrupole and magnetic hyperfine interaction is observed withH _eff=185±2 kOe, 1/2e ² qQ=2.94±0.09 mm/sec and =0.37±0.09. The direction ofH _eff is perpendicular to thec axis of the crystal and at 33° to the <110> direction. Thec axis is determined to be the direction of the intermediate principal EFG tensor axis. Calculations are made using a molecular field term in the Hamiltonian for the Fe²⁺ orbitals. The results of these calculations are used to explain the observed values of 1/2e ² qQ and and permit a determination of the ordring of the T_2g orbitals among the T_2g energy levels.     

Update of electroweak parameters fromZ decays

Based on 520 000 fermion pairs accumulated during the first three years of data collection by the ALEPH detector at LEP, updated values of the resonance parameters of theZ are determined to beM _Z =(91.187±0.009) GeV, _Z =(2.501±0.012) GeV, _had ⁰ =(41.60±0.27) nb, andR =20.78±0.13. The corresponding number of light neutrino species isN =2.97±0.05. The forward-backward asymmetry in lepton-pair decays is used to determine the ratio of vector to axial-vector couplings of leptons:g _V ² (M _Z ² )/g _A ² (M _Z ² )=0.0052±0.0016. Combining this with ALEPH measurements of theb andc quark asymmetries and polarization gives sin² _W ^eff =0.2326±0.0013. Assuming the minimal Standard Model, and including measurements ofM _W /M _Z fromp colliders and neutrino-nucleon scattering, the mass of the top quark is GeV.Deceased     

Critical behavior of short range Potts glasses

We study by means of Monte Carlo simulations and the numerical transfer matrix technique the critical behavior of the short rangep=3 state Potts glass model in dimensionsd=2,3,4 with both Gaussian and bimodal (±J) nearest neighbor interactions on hypercubic lattices employing finite size scaling ideas. Ind=2 in addition the degeneracy of the glass ground state is computed as a function of the number of Potts states forp=3, 4, 5 and compared to that of the antiferromagnetic ground state. Our data indicate a transition into a glass phase atT=0 ind=2 with an algebraic singularity, aT=0 transition ind=3 with an essential singularity of the form exp(const.T ^–2), and an algebraic singularity atT0.25 ind=4. We conclude that the lower critical dimension of the present model isd _c =3 or very close to it. Some of the critical exponents are estimated and their respective values discussed.     

Numerical studies on the anderson localization problem

Our study of Anderson's tight binding model for strongly disordered electronic systems is extended to a numerical treatment of thed c-conductivity atT=0. For 100 × 100 square lattices, 129 × 129 triangular lattices, and for diamond lattices with 27,000 sites, the behaviour of is studied as a function of the Fermi energy and the disorder. The calculations are based on the exact eigenfunction representation of the Kubo formula, which is evaluated by the systematic application of recursion algorithms. Our results are in favour of Mott's original suggestion of a minimum metallic conductivity _min, both in two and three dimensions. In two dimensions we find the universal value of _min=(0.11 ±0.02)e ²/. Based on the thesis of J. Stein, Regensburg 1979     

Isotropic majority-vote model on a square lattice

The stationary critical properties of the isotropic majority vote model on a square lattice are calculated by Monte Carlo simulations and finite size analysis. The critical exponents, , and are found to be the same as those of the Ising model and the critical noise parameter is found to beq _c =0.075±0.001.     

Corrections to cluster-radius scaling for branched polymers and percolation

We report analyses of series enumerations for the mean radius of gyration for isotropic and directed lattice animals and for percolation clusters, in two and three dimensions. We allow for the leading correction to the scaling behaviour and obtain estimates of the leading correction-to-scaling exponent . We find -0.640±0.004 and =0.87±0.07 for isotropic animals in 2d, and =0.64±0.06 in 3d. For directed lattice animals we argue that the leading correction has= or= ; we also estimate =0.82±0.01 and 0.69 ±0.01 ind=2, 3 respectively. For percolation clusters at and abovep _c, we find (p _c) =0.58±0.06 and (p>p _c)=0.84±0.09 in 2d, and (p _c)=0.42±0.11 and (p>p _c)=0.41 ±0.09 in 3d.     

Monte Carlo simulation of very large kinetic Ising models

We study the relaxation of Ising models in three and four dimensions aboveT _c , using multi-spin coding for lattices up to 360³ and 40⁴. The nonlinear relaxation time diverges as (T–T _c )^{–1.05±0.04} in three dimensions, where corrections to scaling are taken into account. In four dimensions the effective exponent is about 0.72; logarithmic correction factors make the analysis difficult here. The linear relaxation time for the asymptotic exponential decay is found to be larger, with effective exponents 1.31 (d=2) and 0.97 (d=4). The difference in the linear and nonlinear relaxation exponents is compatible in three dimensions with the orderparameter exponent , as predicted by Racz.Work supported by SFB 125 Aachen-Jülich-KölnWork started at Department de Physique des Systemes Desordonnes, Universite de Provence, Centre St-Jerome, F-13397 Marseille Cedex 13, France