In the two-dimensional Ising model weak random surface field is predicted to be a marginally irrelevant perturbation at the critical point. We study this question by extensive Monte Carlo simulations for various strength of disorder. The calculated effective (temperature or size dependent) critical exponents fit with the field-theoretical results and can be interpreted in terms of the predicted logarithmic corrections to the pure systems critical behaviour.
The statistics of critical spin-spin correlation functions in Ising systems with non-frustrated disorder are investigated on a strip geometry, via numerical transfer-matrix techniques. Conformal invariance concepts are used, in order to test for logarithmic corrections to pure power-law decay against distance. Fits of our data to conformal-invariance expressions, specific to logarithmic corrections to correlations on strips, give results with the correct sign, for the moments of order $n=0-4$ of the correlation-function distribution. We find an interval of disorder strength along which corrections to pure-system behavior can be decomposed into the product of a known $n$-dependent factor and an approximately $n$-independent one, in accordance with predictions. A phenomenological fitting procedure is proposed, which takes partial account of subdominant terms of correlation-function decay on strips. In the low-disorder limit, it gives results in fairly good agreement with theoretical predictions, provided that an additional assumption is made.
Using high-precision Monte-Carlo simulations based on a parallel version of the Wang-Landau algorithm and finite-size scaling techniques we study the effect of quenched disorder in the crystal-field coupling of the Blume-Capel model on the square lattice. We mainly focus on the part of the phase diagram where the pure model undergoes a continuous transition, known to fall into the universality class of the pure Ising ferromagnet. A dedicated scaling analysis reveals concrete evidence in favor of the strong universality hypothesis with the presence of additional logarithmic corrections in the scaling of the specific heat. Our results are in agreement with an early real-space renormalization-group study of the model as well as a very recent numerical work where quenched randomness was introduced in the energy exchange coupling. Finally, by properly fine tuning the control parameters of the randomness distribution we also qualitatively investigate the part of the phase diagram where the pure model undergoes a first-order phase transition. For this region, preliminary evidence indicate a smoothening of the transition to second-order with the presence of strong scaling corrections.
The three-dimensional bimodal random-field Ising model is investigated using the N-fold version of the Wang-Landau algorithm. The essential energy subspaces are determined by the recently developed critical minimum energy subspace technique, and two implementations of this scheme are utilized. The random fields are obtained from a bimodal discrete $(pmDelta)$ distribution, and we study the model for various values of the disorder strength $Delta$, $Delta=0.5, 1, 1.5$ and 2, on cubic lattices with linear sizes $L=4-24$. We extract information for the probability distributions of the specific heat peaks over samples of random fields. This permits us to obtain the phase diagram and present the finite-size behavior of the specific heat. The question of saturation of the specific heat is re-examined and it is shown that the open problem of universality for the random-field Ising model is strongly influenced by the lack of self-averaging of the model. This property appears to be substantially depended on the disorder strength.
We provide a non-trivial test of supersymmetry in the random-field Ising model at five spatial dimensions, by means of extensive zero-temperature numerical simulations. Indeed, supersymmetry relates correlation functions in a D-dimensional disordered system with some other correlation functions in a D-2 clean system. We first show how to check these relationships in a finite-size scaling calculation, and then perform a high-accuracy test. While the supersymmetric predictions are satisfied even to our high-accuracy at D=5, they fail to describe our results at D=4.
The one-parametric Wang-Landau (WL) method is implemented together with an extrapolation scheme to yield approximations of the two-dimensional (exchange-energy, field-energy) density of states (DOS) of the 3D bimodal random-field Ising model (RFIM). The present approach generalizes our earlier WL implementations, by handling the final stage of the WL process as an entropic sampling scheme, appropriate for the recording of the required two-parametric histograms. We test the accuracy of the proposed extrapolation scheme and then apply it to study the size-shift behavior of the phase diagram of the 3D bimodal RFIM. We present a finite-size converging approach and a well-behaved sequence of estimates for the critical disorder strength. Their asymptotic shift-behavior yields the critical disorder strength and the associated correlation lengths exponent, in agreement with previous estimates from ground-state studies of the model.