A Monte Carlo method has been used to study a simple cubic Ising ferromagnet in a random quenched magnetic field. The Hamiltonian for this model is ℋ==JΣ(ij)σiσj−ΣiHiσi, where σi,σj=±1, J is the nearest-neighbor interaction constant, and the field Hi=tH is fixed at each site with ti=±1 at random and Σti=0. L×L×L lattices with periodic boundary conditions have been studied for a range of H and T. As expected we find a ferromagnetic ordered state which for small H undergoes a second order phase change to the paramagnetic state with increasing temperature. A finite size scaling analysis of the preliminary data suggests that the critical exponent β is substantially smaller (β∼0.2) than the usual 3-dim Ising value of 0.31. Results obtained for small lattices indicate that below kT/J-2 the transition becomes first order suggesting that a tricritical point appears on the critcal field curve.
A Monte Carlo method has been used to study a simple cubic Ising ferromagnet in a random quenched magnetic field. The Hamiltonian for this model is ℋ==JΣ(ij)σiσj−ΣiHiσi, where σi,σj=±1, J is the nearest-neighbor interaction constant, and the field Hi=tH is fixed at each site with ti=±1 at random and Σti=0. L×L×L lattices with periodic boundary conditions have been studied for a range of H and T. As expected we find a ferromagnetic ordered state which for small H undergoes a second order phase change to the paramagnetic state with increasing temperature. A finite size scaling analysis of the preliminary data suggests that the critical exponent β is substantially smaller (β∼0.2) than the usual 3-dim Ising value of 0.31. Results obtained for small lattices indicate that below kT/J-2 the transition becomes first order suggesting that a tricritical point appears on the critcal field curve.
The effects of locally random magnetic fields are considered in a nonequilibrium Ising model defined on a square lattice with nearest-neighbor interactions. In order to generate the random magnetic fields, we have considered random variables { h } that change randomly with time according to a double-Gaussian probability distribution, which consists of two single-Gaussian distributions, centered at + h o and − h o , with the same width σ. This distribution is very general and can recover in appropriate limits the bimodal distribution ( ) and the single-Gaussian one ( h o = 0). We performed Monte Carlo simulations in lattices with linear sizes in the range L = 32–512. The system exhibits ferromagnetic and paramagnetic steady states. Our results suggest the occurrence of first-order phase transitions between the above-mentioned phases at low temperatures and large random-field intensities h o for some small values of the width σ. By means of finite size scaling, we estimate the critical exponents in the low-field region, where we have continuous phase transitions. In addition, we show a sketch of the phase diagram of the model for some values of σ.
The ferromagnetic Ising spin-(1/2 model in a finite simple-cubic lattice \ensuremath{\Lambda} is studied by Monte Carlo methods when two subsets of the lattice sites in \ensuremath{\Lambda}, say ${\ensuremath{\Omega}}^{+}$ and ${\ensuremath{\Omega}}^{\mathrm{\ensuremath{-}}}$, contain (the same number of) spins fixed at \ifmmode\pm\else\textpm\fi{}1, respectively, the global defect concentration being x\ensuremath{\le}0.25. We study the thermodynamic properties of the model for different choices of \ensuremath{\Omega}=${\ensuremath{\Omega}}^{+}$\ensuremath{\cup}${\ensuremath{\Omega}}^{\mathrm{\ensuremath{-}}}$. A finite-size-scaling analysis reveals that the transition remains second order with pure critical exponents for regularly spaced defects, the critical temperature varying with the symmetry of \ensuremath{\Omega}. Any small randomness in \ensuremath{\Omega}, however, makes the transition weakly first order; the transition becomes more abrupt for defects located fully at random, and the long-range order is suppressed when the numbers of defects in ${\ensuremath{\Omega}}^{+}$ and ${\ensuremath{\Omega}}^{\mathrm{\ensuremath{-}}}$ differ from each other. We also discuss our findings in relation to the random-field and frustration problems.
A Monte Carlo method is applied to simulate the static critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder. Numerical results are presented for the spin concentrations of p = 1.0, 0.95, 0.9, 0.8, 0.6 on L × L × L lattices with L = 20–60 under periodic boundary conditions. The critical temperature is determined by the Binder cumulant method. A finite-size scaling technique is used to calculate the static critical exponents α, β, γ, and ν (for specific heat, susceptibility, magnetization, and correlation length, respectively) in the range of p under study. Universality classes of critical behavior are discussed for three-dimensional diluted systems.
Field-induced phase transitions in the bcc Ising antiferromagnet with nearest-neighbor interactions have been studied using a Monte Carlo method. The properties of $N\ifmmode\times\else\texttimes\fi{}N\ifmmode\times\else\texttimes\fi{}N$ lattices with periodic boundary conditions have been determined for $N\ensuremath{\le}16$ and the size effects analyzed using finite-size scaling theory. We find no field dependence of the critical exponents and only a slow variation of critical amplitudes in small fields in agreement with the "smoothness postulate." At low temperatures the critical-field curve shows a weak, but clear maximum which is about 2% greater than the zero-temperature value Hc. Crossover to noncritical exponential behavior due to thermally induced spin flips is observed as we move away from the low-temperature phase boundary.
In this work, we simulated a quenched random mixed-spin Ising model on the square lattice. The model system consists of two different particles with spins σ = 1/2 (states ±1/2) and S = 1 (states ±1, 0). These particles are randomly distributed on the lattice, and we considered only nearest-neighbor interactions. This model can represent a random magnetic binary alloy AxB1-x, obtained from the high-temperature quenching of a liquid mixture. We performed Monte Carlo simulations for several lattice sizes and temperatures, and we found its critical temperature through the reduced fourth-order cumulant. We also determined the magnetization, the susceptibility, and the specific heat as a function of temperature. We used finite-size scaling arguments to estimate the critical exponents β, γ, and ν of the model. We showed that the quenched model is in the same universality class of the two-dimensional pure Ising model. We also investigated the sample to sample fluctuations that occur in the values of the thermodynamic quantities in order to describe the self-averaging properties at the critical point.
The effects of locally random magnetic fields are considered in a nonequilibrium Ising model defined on a square lattice with nearest-neighbors interactions. In order to generate the random magnetic fields, we have considered random variables {h} that change randomly with time according to a double-gaussian probability distribution, which consists of two single gaussian distributions, centered at +ho and −ho, with the same width σ. This distribution is very general, and can recover in appropriate limits the bimodal distribution (σ → 0) and the single gaussian one (ho = 0). We performed Monte Carlo simulations in lattices with linear sizes in the range L = 32 − 512. The system exhibits ferromagnetic and paramagnetic steady states. Our results suggest the occurence of first-order phase transitions between the above-mentioned phases at low temperatures and large random-field intensities ho, for some small values of the width σ. By means of finite size scaling, we estimate the critical exponents in the low-field region, where we have continuous phase transitions. In addition, we show a sketch of the phase diagram of the model for some values of σ. PACS numbers: 05.10.Ln, 05.50.+q, 64.60.De, 75.10Hk, 75.40.MgFirst-order phase transition in a 2D random-field Ising model with conflicting dynamics2 1.