Critical phenomena in the Random Ising Model.

Albert Brown

Published Jan 1, 1986• 0 citations

Abstract

The modern techniques of field theory applied to critical phenomena, are briefly discussed, with particular emphasis on the use of the e-expansion to extract the asymptotic behaviour in the critical region. The dimensional regularisation and minimal subtraction method of renormalisation, Introduced by 't Hooft and Veltman, is outlined. These methods, along with a Monte Carlo simulation on a highly parallel computer, are utilised in extracting information pertaining to the behaviour near criticality of the random, site-diluted Ising model. The universal ratio of the susceptibility amplitude above and below the critical temperature is determined to one higher order than previously in an ɛ-expansion, within the framework of the replica formalism. The calculation is performed using a Taylor expansion, in the number of replicas, of the transverse propagator appearing In the Hamiltonian, about the longitudinal. This higher order correction goes in the right direction for agreement with experiment. An extension of the formalism to embrace the presence of two sets of replica spins is used to derive the o(ɛ¹/²) term in an expansion of the additional static - correlation function, C⁽ˢ⁾(q), arising from the non-interchangeability of the thermal and configurational averages In the random model. The appearance of C(sl(q) in the structure factor Implies that results of neutron scattering experiments to measure the susceptibility amplitude ratio must be re-Interpreted. Finally, analysis of the properties of block spin variables In a Monte Carlo simulation Is used to determine the phase diagram and critical exponents In two-dimensions. Possible Interpretations of the observed evolution of the exponent ratio, ß/υ, under Increasing dilution are discussed.

Related papers

Monte-Carlo simulation study of problems of quantum field theory and critical phenomena.

Jae-Kwon. Kim

Jan 1, 19920 citations

Critical behavior of the random-field Ising model

Misha Gofman, Joan Adler, Amnon Aharony A +2 more

Physical review. B, Condensed matterMar 1, 199659 citations

Three-dimensional randomly dilute Ising model: Monte Carlo results

Pasquale Calabrese, Víctor Martín-Mayor, Andrea Pelissetto C +1 more

arXiv:cond-mat/0306272Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topicsSep 30, 200353 citations

Application of Renormalization-Group Techniques to Random Magnetic Systems.

Ciriyam Jayaprakash

Feb 1, 19790 citations

Criticality of the random-site Ising model: Metropolis, Swendsen-Wang and Wolff Monte Carlo algorithms

Dmytro Ivaneyko, Jaroslav Ilnytskyi, Bertrand Berche +1 more

arXiv:cond-mat/0501291Condensed Matter PhysicsJan 1, 200518 citations

Critical behavior of the three-dimensional random-field Ising model: Two-exponent scaling and discontinuous transition

Heiko Rieger

arXiv:cond-mat/9503041Physical review. B, Condensed matterSep 1, 1995138 citations

Critical phenomena in the Random Ising Model.

Albert Brown

Published Jan 1, 1986• 0 citations

OpenAlex

Abstract

Related papers

Monte-Carlo simulation study of problems of quantum field theory and critical phenomena.

Jae-Kwon. Kim

In chapter one, we explain briefly the continuum limit, scaling, and high temperature expansion of critical phenomena, Monte Carlo algorithms and fitting. In chapter two, different continuum limits of the Ising model in dimensions (D) 2, 3 and 4 are investigated numerically. The data indicate that triviality occurs for D = 4 and fails for D < 4 in each limit. In chapter three, a relation between the critical exponents of the leading and confluent scaling terms is derived using the finite size scaling argument. We also determine the new scaling variable of the 4D Ising model based on a new Monte Carlo simulation data. In chapter four, a Monte Carlo study of two dimensional diluted Ising systems is reported. It is shown that regular dilution does not affect critical exponents, while a random one does, with critical exponents varying continuously with impurity concentration. The importance of fluctuations in producing such effects is emphasized. In chapter five, a different point of view regarding the critical exponent of the specific heat of the 3D Ising model is presented. Based on the analysis of high temperature expansion, finite size scaling and Monte Carlo data in the symmetric phase of the 3D Ising model, it is shown that logarithmic scaling behavior of specific heat is more consistent than power scaling behavior.In chapter one, we explain briefly the continuum limit, scaling, and high temperature expansion of critical phenomena, Monte Carlo algorithms and fitting. In chapter two, different continuum limits of the Ising model in dimensions (D) 2, 3 and 4 are investigated numerically. The data indicate that triviality occurs for D = 4 and fails for D < 4 in each limit. In chapter three, a relation between the critical exponents of the leading and confluent scaling terms is derived using the finite size scaling argument. We also determine the new scaling variable of the 4D Ising model based on a new Monte Carlo simulation data. In chapter four, a Monte Carlo study of two dimensional diluted Ising systems is reported. It is shown that regular dilution does not affect critical exponents, while a random one does, with critical exponents varying continuously with impurity concentration. The importance of fluctuations in producing such effects is emphasized. In chapter five, a different point of view regarding the critical exponent of the specific heat of the 3D Ising model is presented. Based on the analysis of high temperature expansion, finite size scaling and Monte Carlo data in the symmetric phase of the 3D Ising model, it is shown that logarithmic scaling behavior of specific heat is more consistent than power scaling behavior.

Jan 1, 19920 citations

Critical behavior of the random-field Ising model

Misha Gofman, Joan Adler, Amnon Aharony A +2 more

We study the critical properties of the random field Ising model in general dimension d using high-temperature expansions for the susceptibility, \ensuremath{\chi}=${\ensuremath{\sum}}_{\mathit{j}}$ [〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$ ${\mathrm{〉}}_{\mathit{T}}$ -〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$ ${\mathrm{〉}}_{\mathit{T}}$ 〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$ ${\mathrm{〉}}_{\mathit{T}}$ ${]}_{\mathit{h}}$ and the structure factor, G=${\ensuremath{\sum}}_{\mathit{j}}$[〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$ ${\mathrm{〉}}_{\mathit{T}}$ ${]}_{\mathit{h}}$ , where 〈 ${\mathrm{〉}}_{\mathit{T}}$ indicates a canonical average at temperature T for an arbitrary configuration of random fields and [ ${]}_{\mathit{h}}$ indicates an average over random fields. We treated two distributions of random fields, the bimodal in which each ${\mathit{h}}_{\mathit{i}}$ =\ifmmode\pm\else\textpm\fi{} ${\mathit{h}}_{0}$ and a Gaussian distribution in which each ${\mathit{h}}_{\mathit{i}}$ has variance ${\mathit{h}}_{0}^{2}$ . We obtained series for \ensuremath{\chi} and G in the form ${\ensuremath{\sum}}_{\mathit{n}=1,15}$ ${\mathit{a}}_{\mathit{n}}$ (g,d)(J/T ${)}^{\mathit{n}}$ , where J is the exchange constant and the coefficients ${\mathit{a}}_{\mathit{n}}$ (g,d) are polynomials in g\ensuremath{\equiv} ${\mathit{h}}_{0}^{2}$ / ${\mathit{J}}^{2}$ and in d. We assume that as T approaches its critical value, ${\mathit{T}}_{\mathit{c}}$ , one has \ensuremath{\chi}\ensuremath{\sim}(T- ${\mathit{T}}_{\mathit{c}}$ ${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}}$ and G\ensuremath{\sim}(T- ${\mathit{T}}_{\mathit{c}}$ ${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}$. For dimensions above d=2 we find a range of values of g for which the critical exponents obtained from our series seem not to depend on g.For large values of g our results show a g dependence which is attributable to either a tricritical point or a first-order transition. All our results for critical exponents suggest that \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}, in agreement with the two-exponent scaling picture. In addition we have also constructed series for the amplitude ratio, A=(G/${\mathrm{\ensuremath{\chi}}}^{2}$)( ${\mathit{T}}^{2}$ )/( ${\mathit{gJ}}^{2}$ ). We find that A approaches a constant value as T\ensuremath{\rightarrow} ${\mathit{T}}_{\mathit{c}}$ (consistent with \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}) with A\ensuremath{\approxeq}1. It appears that A is somewhat larger for the bimodal than for the Gaussian model, in agreement with a recent analysis at high d. \textcopyright{} 1996 The American Physical Society.We study the critical properties of the random field Ising model in general dimension d using high-temperature expansions for the susceptibility, \ensuremath{\chi}=${\ensuremath{\sum}}_{\mathit{j}}$ [〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$${\mathrm{〉}}_{\mathit{T}}$-〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{〉}}_{\mathit{T}}$〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$${\mathrm{〉}}_{\mathit{T}}$${]}_{\mathit{h}}$ and the structure factor, G=${\ensuremath{\sum}}_{\mathit{j}}$[〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$${\mathrm{〉}}_{\mathit{T}}$${]}_{\mathit{h}}$, where 〈${\mathrm{〉}}_{\mathit{T}}$ indicates a canonical average at temperature T for an arbitrary configuration of random fields and [${]}_{\mathit{h}}$ indicates an average over random fields. We treated two distributions of random fields, the bimodal in which each ${\mathit{h}}_{\mathit{i}}$=\ifmmode\pm\else\textpm\fi{}${\mathit{h}}_{0}$ and a Gaussian distribution in which each ${\mathit{h}}_{\mathit{i}}$ has variance ${\mathit{h}}_{0}^{2}$. We obtained series for \ensuremath{\chi} and G in the form ${\ensuremath{\sum}}_{\mathit{n}=1,15}$${\mathit{a}}_{\mathit{n}}$(g,d)(J/T${)}^{\mathit{n}}$, where J is the exchange constant and the coefficients ${\mathit{a}}_{\mathit{n}}$(g,d) are polynomials in g\ensuremath{\equiv}${\mathit{h}}_{0}^{2}$/${\mathit{J}}^{2}$ and in d. We assume that as T approaches its critical value, ${\mathit{T}}_{\mathit{c}}$, one has \ensuremath{\chi}\ensuremath{\sim}(T-${\mathit{T}}_{\mathit{c}}$${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}}$ and G\ensuremath{\sim}(T-${\mathit{T}}_{\mathit{c}}$${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}$. For dimensions above d=2 we find a range of values of g for which the critical exponents obtained from our series seem not to depend on g.For large values of g our results show a g dependence which is attributable to either a tricritical point or a first-order transition. All our results for critical exponents suggest that \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}, in agreement with the two-exponent scaling picture. In addition we have also constructed series for the amplitude ratio, A=(G/${\mathrm{\ensuremath{\chi}}}^{2}$)(${\mathit{T}}^{2}$)/(${\mathit{gJ}}^{2}$). We find that A approaches a constant value as T\ensuremath{\rightarrow}${\mathit{T}}_{\mathit{c}}$ (consistent with \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}) with A\ensuremath{\approxeq}1. It appears that A is somewhat larger for the bimodal than for the Gaussian model, in agreement with a recent analysis at high d. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matterMar 1, 199659 citations

Three-dimensional randomly dilute Ising model: Monte Carlo results

Pasquale Calabrese, Víctor Martín-Mayor, Andrea Pelissetto C +1 more

We perform a high-statistics simulation of the three-dimensional randomly dilute Ising model on cubic lattices L3 with L< or =256. We choose a particular value of the density, x=0.8, for which the leading scaling corrections are suppressed. We determine the critical exponents, obtaining nu=0.683(3), eta=0.035(2), beta=0.3535(17), and alpha=-0.049(9), in agreement with previous numerical simulations. We also estimate numerically the fixed-point values of the four-point zero-momentum couplings that are used in field-theoretical fixed-dimension studies. Although these results somewhat differ from those obtained using perturbative field theory, the field-theoretical estimates of the critical exponents do not change significantly if the Monte Carlo result for the fixed point is used. Finally, we determine the six-point zero-momentum couplings, relevant for the small-magnetization expansion of the equation of state, and the invariant amplitude ratio R(+)(xi) that expresses the universality of the free-energy density per correlation volume. We find R(+)(xi)=0.2885(15).We perform a high-statistics simulation of the three-dimensional randomly dilute Ising model on cubic lattices L3 with L< or =256. We choose a particular value of the density, x=0.8, for which the leading scaling corrections are suppressed. We determine the critical exponents, obtaining nu=0.683(3), eta=0.035(2), beta=0.3535(17), and alpha=-0.049(9), in agreement with previous numerical simulations. We also estimate numerically the fixed-point values of the four-point zero-momentum couplings that are used in field-theoretical fixed-dimension studies. Although these results somewhat differ from those obtained using perturbative field theory, the field-theoretical estimates of the critical exponents do not change significantly if the Monte Carlo result for the fixed point is used. Finally, we determine the six-point zero-momentum couplings, relevant for the small-magnetization expansion of the equation of state, and the invariant amplitude ratio R(+)(xi) that expresses the universality of the free-energy density per correlation volume. We find R(+)(xi)=0.2885(15).

arXiv:cond-mat/0306272Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topicsSep 30, 200353 citations

Application of Renormalization-Group Techniques to Random Magnetic Systems.

Ciriyam Jayaprakash

"Renorma1ization-group methods have been applied in the study of quenched random magnetic systems in recent years. We begin with a brief review of second-order phase transitions in pure, homogeneous systems and also of the .. renorma1ization group framework. Then we provide an introduction to quenched random magnetic systems. Next, momentum-space methods and position-space techniques as applied to quenched random magnets are outlined and compared. Grinstein and Luther applied the Wilson-Fisher E-expansion to random n-vector models; Khme1'nitsky discovered that the random Ising model (n = 1) possessed a ""random"" fixed point of 0(squareroot(e)). This fixed point was found to have one marginal and one irrelevant operator. We have investigated the stability of this fixed point using Ca11an-Symanzik equations and renorma1ized perturbation theory. We find the fixed point stable in the next order; we have also obtained critical exponents to one higher order. Next, position-space techniques are used to study some simple model systems. In addition to critical exponents, global thermodynamic properties are determined. These calculations are based on the Migda1-Kadanoff approximate recursion relations suitably generalized to the inhomogeneous case. Firstly we study the randomly bond-dilute two-dimensional nearest- neighbor Ising model on a square lattice. Calculations give both thermal and magnetic exponents associated with the percolative fixed point. Differential recursion relations yield a phase diagram which is in quantitative agreement with all known results. Curves for the specific heat, percolation probability, and magnetization are displayed. The critical region of the specific heat becomes unobservably narrow well above the percolation threshold Pc. This provides a possible explanation for the apparent specific-heat rounding in certain experiments. We then study the Edwards-Anderson model of a spin glass. The current theoretical situation, which is far from satisfactory at present, is briefly reviewed. We treat the spin-l/2 Ising model with independently random nearestneighbor interactions in dimensionalities d = 2, 3, and 4. The phase diagram, which is in qualitative agreement with mean-field results, exhibits paramagnetic, ferromagnetic, antiferromagnetic, and spin-glass phases. The spinglass and paramagnetic phases meet along an extended second-order phase boundary, which terminates in two tricritical points. Critical and tricritical exponents are calculated. The spin-glass specific-heat exponent turns out to be large and negative, compatibly with recent experiments which show a rounded specific heat anomaly. Global specific-heat curves are also displayed for d = 2.""Renorma1ization-group methods have been applied in the study of quenched random magnetic systems in recent years. We begin with a brief review of second-order phase transitions in pure, homogeneous systems and also of the .. renorma1ization group framework. Then we provide an introduction to quenched random magnetic systems. Next, momentum-space methods and position-space techniques as applied to quenched random magnets are outlined and compared. Grinstein and Luther applied the Wilson-Fisher E-expansion to random n-vector models; Khme1'nitsky discovered that the random Ising model (n = 1) possessed a ""random"" fixed point of 0(squareroot(e)). This fixed point was found to have one marginal and one irrelevant operator. We have investigated the stability of this fixed point using Ca11an-Symanzik equations and renorma1ized perturbation theory. We find the fixed point stable in the next order; we have also obtained critical exponents to one higher order. Next, position-space techniques are used to study some simple model systems. In addition to critical exponents, global thermodynamic properties are determined. These calculations are based on the Migda1-Kadanoff approximate recursion relations suitably generalized to the inhomogeneous case. Firstly we study the randomly bond-dilute two-dimensional nearest- neighbor Ising model on a square lattice. Calculations give both thermal and magnetic exponents associated with the percolative fixed point. Differential recursion relations yield a phase diagram which is in quantitative agreement with all known results. Curves for the specific heat, percolation probability, and magnetization are displayed. The critical region of the specific heat becomes unobservably narrow well above the percolation threshold Pc. This provides a possible explanation for the apparent specific-heat rounding in certain experiments. We then study the Edwards-Anderson model of a spin glass. The current theoretical situation, which is far from satisfactory at present, is briefly reviewed. We treat the spin-l/2 Ising model with independently random nearestneighbor interactions in dimensionalities d = 2, 3, and 4. The phase diagram, which is in qualitative agreement with mean-field results, exhibits paramagnetic, ferromagnetic, antiferromagnetic, and spin-glass phases. The spinglass and paramagnetic phases meet along an extended second-order phase boundary, which terminates in two tricritical points. Critical and tricritical exponents are calculated. The spin-glass specific-heat exponent turns out to be large and negative, compatibly with recent experiments which show a rounded specific heat anomaly. Global specific-heat curves are also displayed for d = 2."

Feb 1, 19790 citations

Criticality of the random-site Ising model: Metropolis, Swendsen-Wang and Wolff Monte Carlo algorithms

Dmytro Ivaneyko, Jaroslav Ilnytskyi, Bertrand Berche +1 more

We apply numerical simulations to study the criticality of the 3D Ising model with random site quenched dilution. The emphasis is given to the issues not being discussed in detail before. In particular, we attempt a comparison of different Monte Carlo techniques, discussing regions of their applicability and advantages/disadvantages depending on the aim of a particular simulation set. Moreover, besides evaluation of the critical indices we estimate the universal ratio $Γ^+/Γ^-$ for the magnetic susceptibility critical amplitudes. Our estimate $Γ^+/Γ^- = 1.67 \pm 0.15$ is in a good agreement with the recent MC analysis of the random-bond Ising model giving further support that both random-site and random-bond dilutions lead to the same universality class.We apply numerical simulations to study the criticality of the 3D Ising model with random site quenched dilution. The emphasis is given to the issues not being discussed in detail before. In particular, we attempt a comparison of different Monte Carlo techniques, discussing regions of their applicability and advantages/disadvantages depending on the aim of a particular simulation set. Moreover, besides evaluation of the critical indices we estimate the universal ratio $Γ^+/Γ^-$ for the magnetic susceptibility critical amplitudes. Our estimate $Γ^+/Γ^- = 1.67 \pm 0.15$ is in a good agreement with the recent MC analysis of the random-bond Ising model giving further support that both random-site and random-bond dilutions lead to the same universality class.

arXiv:cond-mat/0501291Condensed Matter PhysicsJan 1, 200518 citations

Critical behavior of the three-dimensional random-field Ising model: Two-exponent scaling and discontinuous transition

Heiko Rieger

In extensive Monte Carlo simulations the phase transition of the random-field Ising model in three dimensions is investigated. The values of the critical exponents are determined via finite-size scaling. For a Gaussian distribution of the random fields it is found that the correlation length \ensuremath{\xi} diverges with an exponent \ensuremath{\nu}=1.1\ifmmode\pm\else\textpm\fi{}0.2 at the critical temperature and that \ensuremath{\chi}\ensuremath{\sim}${\ensuremath{\xi}}^{2\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}}$ with \ensuremath{\eta}=0.50\ifmmode\pm\else\textpm\fi{}0.05 for the connected susceptibility and ${\mathrm{\ensuremath{\chi}}}_{\mathrm{dis}}$\ensuremath{\sim}${\ensuremath{\xi}}^{4\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}$ with \ensuremath{\eta}\ifmmode\bar\else\textasciimacron\fi{}=1.03\ifmmode\pm\else\textpm\fi{}0.05 for the disconnected susceptibility. Together with the amplitude ratio A=${\mathrm{lim}}_{\mathit{T}\ensuremath{\rightarrow}\mathit{T}\mathit{c}}$${\mathrm{\ensuremath{\chi}}}_{\mathrm{dis}}$/${\mathrm{\ensuremath{\chi}}}^{2}$( ${\mathit{h}}_{\mathit{r}}$ /T ${)}^{2}$ being close to one this gives further support for a two-exponent scaling scenario implying \ensuremath{\eta}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\eta}. The magnetization behaves discontinuously at the transition, i.e., \ensuremath{\beta}=0. However, no divergence for the specific heat and in particular no latent heat is found. Also the probability distribution of the magnetization does not show a multipeak structure that would be characteristic for the phase-coexistence at first-order phase-transition points.In extensive Monte Carlo simulations the phase transition of the random-field Ising model in three dimensions is investigated. The values of the critical exponents are determined via finite-size scaling. For a Gaussian distribution of the random fields it is found that the correlation length \ensuremath{\xi} diverges with an exponent \ensuremath{\nu}=1.1\ifmmode\pm\else\textpm\fi{}0.2 at the critical temperature and that \ensuremath{\chi}\ensuremath{\sim}${\ensuremath{\xi}}^{2\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}}$ with \ensuremath{\eta}=0.50\ifmmode\pm\else\textpm\fi{}0.05 for the connected susceptibility and ${\mathrm{\ensuremath{\chi}}}_{\mathrm{dis}}$\ensuremath{\sim}${\ensuremath{\xi}}^{4\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}$ with \ensuremath{\eta}\ifmmode\bar\else\textasciimacron\fi{}=1.03\ifmmode\pm\else\textpm\fi{}0.05 for the disconnected susceptibility. Together with the amplitude ratio A=${\mathrm{lim}}_{\mathit{T}\ensuremath{\rightarrow}\mathit{T}\mathit{c}}$${\mathrm{\ensuremath{\chi}}}_{\mathrm{dis}}$/${\mathrm{\ensuremath{\chi}}}^{2}$(${\mathit{h}}_{\mathit{r}}$/T${)}^{2}$ being close to one this gives further support for a two-exponent scaling scenario implying \ensuremath{\eta}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\eta}. The magnetization behaves discontinuously at the transition, i.e., \ensuremath{\beta}=0. However, no divergence for the specific heat and in particular no latent heat is found. Also the probability distribution of the magnetization does not show a multipeak structure that would be characteristic for the phase-coexistence at first-order phase-transition points.

arXiv:cond-mat/9503041Physical review. B, Condensed matterSep 1, 1995138 citations

{\mathrm{〉}}_{\mathit{T}}

-〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$

{\mathrm{〉}}_{\mathit{T}}

〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$

{\mathrm{〉}}_{\mathit{T}}

{]}_{\mathit{h}}

and the structure factor, G=${\ensuremath{\sum}}_{\mathit{j}}$[〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$

{\mathrm{〉}}_{\mathit{T}}

{]}_{\mathit{h}}

, where 〈

{\mathrm{〉}}_{\mathit{T}}

indicates a canonical average at temperature T for an arbitrary configuration of random fields and [

{]}_{\mathit{h}}

indicates an average over random fields. We treated two distributions of random fields, the bimodal in which each

{\mathit{h}}_{\mathit{i}}

=\ifmmode\pm\else\textpm\fi{}

{\mathit{h}}_{0}

and a Gaussian distribution in which each

{\mathit{h}}_{\mathit{i}}

has variance

{\mathit{h}}_{0}^{2}

. We obtained series for \ensuremath{\chi} and G in the form ${\ensuremath{\sum}}_{\mathit{n}=1,15}$

{\mathit{a}}_{\mathit{n}}

(g,d)(J/T

{)}^{\mathit{n}}

, where J is the exchange constant and the coefficients

{\mathit{a}}_{\mathit{n}}

(g,d) are polynomials in g\ensuremath{\equiv}

{\mathit{h}}_{0}^{2}

{\mathit{J}}^{2}

and in d. We assume that as T approaches its critical value,

{\mathit{T}}_{\mathit{c}}

, one has \ensuremath{\chi}\ensuremath{\sim}(T-

{\mathit{T}}_{\mathit{c}}

${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}}$ and G\ensuremath{\sim}(T-

{\mathit{T}}_{\mathit{c}}

${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}$. For dimensions above d=2 we find a range of values of g for which the critical exponents obtained from our series seem not to depend on g.For large values of g our results show a g dependence which is attributable to either a tricritical point or a first-order transition. All our results for critical exponents suggest that \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}, in agreement with the two-exponent scaling picture. In addition we have also constructed series for the amplitude ratio, A=(G/${\mathrm{\ensuremath{\chi}}}^{2}$)(

{\mathit{T}}^{2}

)/(

{\mathit{gJ}}^{2}

). We find that A approaches a constant value as T\ensuremath{\rightarrow}

{\mathit{T}}_{\mathit{c}}

(consistent with \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}) with A\ensuremath{\approxeq}1. It appears that A is somewhat larger for the bimodal than for the Gaussian model, in agreement with a recent analysis at high d. \textcopyright{} 1996 The American Physical Society.

Γ^+/Γ^-

for the magnetic susceptibility critical amplitudes. Our estimate

Γ^+/Γ^- = 1.67 \pm 0.15

is in a good agreement with the recent MC analysis of the random-bond Ising model giving further support that both random-site and random-bond dilutions lead to the same universality class.

{\mathit{h}}_{\mathit{r}}

{)}^{2}

being close to one this gives further support for a two-exponent scaling scenario implying \ensuremath{\eta}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\eta}. The magnetization behaves discontinuously at the transition, i.e., \ensuremath{\beta}=0. However, no divergence for the specific heat and in particular no latent heat is found. Also the probability distribution of the magnetization does not show a multipeak structure that would be characteristic for the phase-coexistence at first-order phase-transition points.