Magnetic correlations in three-dimensional ising spin glasses

By a recursive method numerically exact free energies are calculated forL×L×M Ising lattices with random bonds andL=4, 4M10, applying free boundaries in the direction where the lattice is less small and otherwise periodic boundary conditions. Both for the±J model and the gaussian model the specific heat is in fair agreement with Monte Carlo results obtained for much larger lattices. However, the correlation function [S ₀ S _{R T} ² ]_av is found to decay exponentially with distanceR [for 1R9] at temperatures far below the apparent freezing temperatures of the Monte Carlo simulations, implying that there is no nonzero Edwards-Anderson order parameter in equilibrium. This behavior is qualitatively different from Mattis spin glasses (or Ising ferromagnets) where even smaller lattices show pronounced magnetic order at low temperatures. As the Monte Carlo results give evidence for a nonzero Edwards-Anderson order parameter (for not too long observation times), which is fairly independent of lattice size down to sizes of 4³, we suggest that Edwards-Anderson ordering is a nonequilibrium phenomenon visible only in studying dynamic properties.     

Monte Carlo calculations on the gauge Potts model

Theq-state gauge Potts modelP _q ind-dimensions has been studied using Monte Carlo techniques. Ford=2 no phase transitions were detected. TheP ₂ model ford=3 shows a second order phase transition. On the other hand, all thed=3 (q2) andd=4 cases studied show first order phase transitions. In these cases, it was possible to estimate transition coupling parameters as well as latent heat. For selected cases, a study of the behavior of the Wilson loop factor was done.Partially supported by CONICET, CIC Pcia. de Buenos Aires and SUBCYT, Argentina     

Comparative monte Carlo study of ising spin glasses in two to five dimensions

Monte Carlo simulations of the Edwards-Anderson-Ising spin glass with Gaussian distribution of nearest-neighbor exchange forces in four and five dimensions are performed to check the speculation thatd=4 is the lower critical dimensionality. In contrast to this expectation we find no qualitative difference at all to the results in two and three dimensions. We still find that on not too long time-scales there is an apparently rather well defined freezing temperatureT _f , where the susceptibility has a cusp, and belowT _f nonzero order parametersq, can be found as ford=2, 3. But even ford=5 the decay of the Edwards-Anderson order parameter belowT _f is found to be consistent with a logarithmic variation over several decades of observation time. The possible interpretations of this result are discussed. Our data thus suggest that either there is no equilibrium phase transition in all these dimensions, or more likely that a phase transition exists for 2d5 but the properties of the ordered phase may be rather peculiar.Sonderforschungsbereich 125 Aachen-Jülich-Köln     

Correlation inequalities for vector spin models

Correlation inequalities forn-vector spin models (n 2) are reviewed. A relatively simple and unified derivation of the inequalities is achieved, using duplicate variable methods, for spin dimensionalitiesn=2 (plane rotator model),n=3 (classical Heisenberg model), andn=4. Although correlation inequalities are lacking forn > 4, new proofs are presented for the comparison inequalities relating correlations for systems with arbitrary spin dimensionality to corresponding correlations for systems with low spin dimensionality (n = 1 or 2).Research supported by National Science Foundation under Grant DMR 76-23071.     

Annealedn-vectorp-spin model

A disorderedn-vector model withp spin interactions is introduced and studied in mean field theory for the annealed case. We present complete solutions for the casesn=2 andn=3, and have obtained explicit order parameter equations for all the stable solutions for arbitraryn. For alln andp we find one stable high-temperature phase and one stable low-temperature phase. The phase transition is of first order. Forn=2, it is continuous in the order parameters for p4 and has a jump discontinuity in the order parameters ifp>4. Forn=3, it has a jump discontinuity in the order parameters for allp.     

Some recent developments in spin glasses

I give some experimental and theoretical background to spin glasses, and then discuss the nature of the phase transition in spin glasses withvector spins. Results of Monte Carlo simulations of the Heisenberg spin glass model in three dimensions are presented. A finite-size scaling analysis of the correlation length of the spins and chiralities shows that there is a single, finite-temperature transition at which both spins and chiralities order.     

Dynamical interaction of antiferromagnetic subsystems: a neutron scattering study of the spinwave spectrum of the garnet Fe2Ca3(GeO4)3

The Fe³⁺ ions in the garnet Ca₃Fe₂Ge₃O₁₂ form two identical antiferromagnetic subsystems. The interaction between the two subsystems is vanishing within molecular field approximation forq=0. A coupling appears due to the spin fluctuations. The dynamics of the system is described by the Hamiltonian for a Heisenberg antiferromagnet. Symmetry requirements impose two exchange parameters between the sublattices (nearest neighbours)J ₁ in the direction of the 3-fold axis andJ' ₁ in the other three space diagonals. The interaction within each sublattice (second nearest neighbours) is described by the exchange parameterJ ₂. The measured spin wave dispersion curves for the three principal symmetry directions are very well reproduced by a model calculation withJ ₁=-0.909(9) K,J' ₁=-0.307(8) K andJ ₂=-0.615(2)K. The observed intensities are in agreement with predictions from the model. Forq0 the model predicts two acoustic branches going towards zero frequency. A calculation beyond linear spin wave theory forq=0 predicts a quantum gap for the lower acoustic branch. This gap has been found at 0.033(4) THz. An anisotropy gap of 0.007 THz has been taken from the literature.     

Excitations in Heisenberg spin glasses

Within the RPA approach forT=0, the excitations of the Heisenberg spin glass system Eu _x Sr_1–x S are studied by numerical methods, using a continued fraction algorithm. Both the density of statesg(E) and also the spectral functionS(q,E) are calculated for systems with (16)³ sites, withx=0.4, 0.5, and 0.6 (spin glass phase), and also forx0.7 (ferromagnetic phase). Forq-vectors within the (1,1,1) plane,S(q,E) shows magnon peaks even in the spin glass phase, over the whole range ofq. However, these peaks are quite broad, and there is considerable intensity at small energies even for largeq, leading to a finite intercept ofg(E) forE0. Over a large temperature range, the specific heat is approximately linear inT forx0.7.     

On the critical dynamics of disordered spin systems

The dynamic critical behaviour of spin systems with quenched impurities, and of amorphous spin systems as characterized by the additional presence of random anisotropy directions, is studied by renormalization group methods to second order in=4–d. For the Halperin-Hohenberg-Ma model with purely relaxational dynamics it is concluded that in three dimensions (d=3) the critical slowing down should be enhanced by impurities for systems with Ising type statics, whereas there is no change forXY- and Heisenberg systems. For amorphous systems, however, the critical dynamics should change also in theXY- and Heisenberg cases. Furthermore, it is concluded that additional conserved, but noncritical modes become always statically decoupled from the order parameter for systems with impurities, but not for amorphous systems. Thus, for the impure system, the energy density mode and the asymmetric models of Halperin, Hohenberg and Siggia are ruled out. But the effects of dynamic coupling remain: Especially, the relationz=d/2 for the dynamic exponent of planar and isotropic antiferromagnets is modified for impure or amorphous systems.     

Monte Carlo studies of the internal energy and specific heat of a classical Heisenberg spin glass

We present data obtained from Monte Carlo studies of the Edwards-Anderson model of a classical Heisenberg spin glass. The internal energy is calculated for a 10 × 10 × 10 array. The specific heat, obtained from the temperature derivative of the internal energy, has a rounded peak at slightly less than one half the mean field transition temperature.     

Thermodynamics of the spin-S antiferromagnetic Heisenberg chain

The numerical solution of the Bethe ansatz equations of an integrableSU (2)-invariant generalization of the spin-S antiferromagnetic Heisenberg chain in zero magnetic field is presented. The thermodynamics is obtained numerically. The temperature dependence of the entropy, specific heat and susceptibility is presented forS5/2. The results are compared to those of then-channel Kondo problem with a spin-S impurity withn=2S.     

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.     

Classical Heisenberg antiferromagnets with nearest and next-nearest neighbor interactions on the face-centered cubic lattice: a model for EuTe?

Magnetic properties of the Heisenberg antiferromagnet with spin quantum numberS on the face-centered cubic lattice are studied as function of temperature and magnetic field, using molecular field approximation and Monte Carlo methods. In order to model Europiumtelluride, we use isotropic exchange interactions between nearest- and nextnearest neighbors; the values of these exchange constants are taken from experiments. In addition, a pseudo-dipolar anisotropy (truncated after the next-nearest neighbor distance) is included; the molecular field calculations also are performed with the full dipolar of real EuTe in two respects: the structure in zero magnetic field involves 8 sublattices in the model rather than only two; the bicritical point, above which in the temperatureT magnetic fieldH plane the spin flop phase appears, occurs atH=0 in the model rather than at nonzero field. Possible additional interactions responsible for these discrepancies are discussed. Applying finite size scaling techniques we give also a preliminary analysis of the critical behavior of the model.     

Disordered system withn orbitals per site: 1/n expansion

Averaged Green's functions for a disordered electronic system withn orbitals per site are expanded in powers of 1/n. These expansions should be valid in the region of extended states. The expansion coefficients for the d.c. conductivity are finite for dimensionalityd>2 and diverge asd approaches 2. Similarities of two types of two-particle Green's functions with the transverse and longitudinal susceptibilities of a ferromagnet with broken continuous symmetry are pointed out. Arguments for two being the lower critical dimensionality for the hydrodynamics and the mobility edge are given. Provided our series can be exponentiated we find that no metallic conductivity exists for finiten andd=2 in one of our models. Critical exponents ford infinitesimal above two are given. In this limitv diverges like 1/(d–2) and the conductivity vanishes linearly at the mobility edge.The diagrams of the Green's functions are given in terms of vertices of short-range order and of the two-particle propagators of then= limit. Diagrams withs loops contribute in ordern ^–s . The diagrams can be rearranged so that a number of vertices vanishes like the square of the wavevector. This feature prevents infrared divergencies for the d.c. conductivity ford>2.Work supported in part by a DFG fellowship (R.O.), by the Material Research Laboratory of the National Science Foundation at the University of Chicago (F.W.), and by the Deutsche Forschungsgemeinschaft through the Sonderforschungsbereich 123 (Stochastic Mathematical Models) at the Universität Heidelberg     

Dynamic weighting in simulations of spin systems

We establish the invariance property of a dynamically weighted Monte Carlo process and apply the method to the simulation of spin glasses and Ising models. For the two-dimensional Edwards-Anderson model, we obtain an ergodicity time of O(L^2.44), where L is the linear dimension of the lattice. The results suggest that dynamic weighting is a highly promising new tool for Monte Carlo simulation.     

Single-cluster Monte Carlo dynamics for the Ising model

We present an extensive study of a new Monte Carlo acceleration algorithm introduced by Wolff for the Ising model. It differs from the Swendsen-Wang algorithm by growing and flipping single clusters at a random seed. In general, it is more efficient than Swendsen-Wang dynamics ford>2, giving zero critical slowing down in the upper critical dimension. Monte Carlo simulations give dynamical critical exponentsz _w=0.33±0.05 and 0.44+0.10 ind=2 and 3, respectively, and numbers consistent withz _w=0 ind=4 and mean-field theory. We present scaling arguments which indicate that the Wolff mechanism for decorrelation differs substantially from Swendsen-Wang despite the apparent similarities of the two methods.     

Quenchedn-vectorp-spin model

A disorderedn-vector model withp spin interactions previously introduced is studied for the quenched case by means of the replica method and a generalized Parisi theory. We present formal solutions for generaln andp and then study the casep . The high-temperature solution is stable at all temperatures and there is only one phase transition at a temperatureT _g. Only longitudinal lowtemperature solutions are possible. There is one spin-glass solution, and it is stable for allT _g. The phase transition atT _g is of first order and displays a jump discontinuity in the order parametersq _j ^(L) andd. The spin-glass free energy is temperature dependent forn > 1 while it is constant whenn = 1.     

Quantum statistical monte carlo methods and applications to spin systems

A short review is given concerning the quantum statistical Monte Carlo method based on the equivalence theorem⁽¹⁾ thatd-dimensional quantum systems are mapped onto (d+1)-dimensional classical systems. The convergence property of this approximate tansformation is discussed in detail. Some applications of this geneal appoach to quantum spin systems are reviewed. A new Monte Carlo method, “thermo field Monte Carlo method,” is presented, which is an extension of the projection Monte Carlo method at zero temperature to that at finite temperatures. Invited talk presented at “Frontiers of Quantum Monte Carlo,” Los Alamos National Laboratory, September 3–6, 1985.     

Critical dynamics and statics of uniaxial dipolar magnets

The critical behaviour of uniaxial ferromagnets is studied by an anisotropic renormalization transformation. For dimensionalityd smaller than the critical dimensionalityd _c (d _u ), which depends on the dimensionality of the uniaxial directiond _u , the critical exponents are computed by the-expansion. The critical dynamics is based on a timedependent Ginzburg Landau model, for a non-conserved and for a conserved order parameter. For the experimentally relevant cased=3 andd _u =1 the logarithmic corrections in the frequency dependent susceptibility are computed.Work supported by the Fonds zur Förderung der wissenschaftlichen Forschung