We report on large-scale Wang-Landau Monte Carlo simulations of the critical behavior of two spin models in two- (2d) and three-dimensions (3d), namely the 2d random-bond Ising model and the pure 3d Blume-Capel model at zero crystal-field coupling. T
he numerical data we obtain and the relevant finite-size scaling analysis provide clear answers regarding the universality aspects of both models. In particular, for the random-bond case of the 2d Ising model the theoretically predicted strong universalitys hypothesis is verified, whereas for the second-order regime of the Blume-Capel model, the expected $d=3$ Ising universality is verified. Our study is facilitated by the combined use of the Wang-Landau algorithm and the critical energy subspace scheme, indicating that the proposed scheme is able to provide accurate results on the critical behavior of complex spin systems.
Using high-precision Monte Carlo simulations and finite-size scaling we study the effect of quenched disorder in the exchange couplings on the Blume-Capel model on the square lattice. The first-order transition for large crystal-field coupling is sof
tened to become continuous, with a divergent correlation length. An analysis of the scaling of the correlation length as well as the susceptibility and specific heat reveals that it belongs to the universality class of the Ising model with additional logarithmic corrections observed for the Ising model itself if coupled to weak disorder. While the leading scaling behavior in the disordered system is therefore identical between the second-order and first-order segments of the phase diagram of the pure model, the finite-size scaling in the ex-first-order regime is affected by strong transient effects with a crossover length scale $L^{ast} approx 32$ for the chosen parameters.
We report a high-precision numerical estimation of the critical exponent $alpha$ of the specific heat of the random-field Ising model in four dimensions. Our result $alpha = 0.12(1)$ indicates a diverging specific-heat behavior and is consistent with
the estimation coming from the modified hyperscaling relation using our estimate of $theta$ via the anomalous dimensions $eta$ and $bar{eta}$. Our analysis benefited form a high-statistics zero-temperature numerical simulation of the model for two distributions of the random fields, namely a Gaussian and Poissonian distribution, as well as recent advances in finite-size scaling and reweighting methods for disordered systems. An original estimate of the critical slowing down exponent $z$ of the maximum-flow algorithm used is also provided.
We study the effect of interfacial phenomena in two-dimensional perfect and random (or disordered) $q$-state Potts models with continuous phase transitions, using, mainly, Monte Carlo techniques. In particular, for the total interfacial adsorption, t
he critical behavior, including corrections to scaling, are analyzed. The role of randomness is scrutinized. Results are discussed applying scaling arguments and invoking findings for bulk critical properties. In all studied cases, i.e., $q = 3$, $4$, and $q = 8$, the spread of the interfacial adsorption profiles is observed to increase linearly with the lattice size at the bulk transition point.
We enlighten some critical aspects of the three-dimensional ($d=3$) random-field Ising model from simulations performed at zero temperature. We consider two different, in terms of the field distributio
We study the $pm J$ three-dimensional Ising model with a longitudinal anisotropic bond randomness on the simple cubic lattice. The random exchange interaction is applied only in the $z$ direction, whereas in the other two directions, $xy$ - planes, w
e consider ferromagnetic exchange. By implementing an effective parallel tempering scheme, we outline the phase diagram of the model and compare it to the corresponding isotropic one, as well as to a previously studied anisotropic (transverse) case. We present a detailed finite-size scaling analysis of the ferromagnetic - paramagnetic and spin glass - paramagnetic transition lines, and we also discuss the ferromagnetic - spin glass transition regime. We conclude that the present model shares the same universality classes with the isotropic model, but at the symmetric point has a considerably higher transition temperature from the spin-glass state to the paramagnetic phase. Our data for the ferromagnetic - spin glass transition line are supporting a forward behavior in contrast to the reentrant behavior of the isotropic model.
We investigate the dependence of the critical Binder cumulant of the magnetization and the largest Fortuin-Kasteleyn cluster on the boundary conditions and aspect ratio of the underlying square Ising lattices. By means of the Swendsen-Wang algorithm,
we generate numerical data for large system sizes and we perform a detailed finite-size scaling analysis for several values of the aspect ratio $r$, for both periodic and free boundary conditions. We estimate the universal probability density functions of the largest Fortuin-Kasteleyn cluster and we compare it to those of the magnetization at criticality. It is shown that these probability density functions follow similar scaling laws, and it is found that the values of the critical Binder cumulant of the largest Fortuin-Kasteleyn cluster are upper bounds to the values of the respective order-parameters cumulant, with a splitting behavior for large values of the aspect ratio. We also investigate the dependence of the amplitudes of the magnetization and the largest Fortuin-Kasteleyn cluster on the aspect ratio and boundary conditions. We find that the associated exponents, describing the aspect ratio dependencies, are different for the magnetization and the largest Fortuin-Kasteleyn cluster, but in each case are independent of boundary conditions.
We study the Blume-Capel model on the square lattice. To allow for wetting and interfacial adsorption, the spins on opposite boundaries are fixed in two different states, +1 and -1, with reduced couplings at one of the boundaries. Using mainly Monte
Carlo techniques, of Metropolis and Wang-Landau type, phase diagrams showing bulk and wetting transitions are determined. The role of the non-boundary state, 0, adsorbed preferably at the interface between -1 and +1 rich regions, is elucidated.
We implement a two-stage approach of the Wang-Landau algorithm to investigate the critical properties of the 3D Ising model with quenched bond randomness. In particular, we consider the case where disorder couples to the nearest-neighbor ferromagneti
c interaction, in terms of a bimodal distribution of strong versus weak bonds. Our simulations are carried out for large ensembles of disorder realizations and lattices with linear sizes $L$ in the range $L=8-64$. We apply well-established finite-size scaling techniques and concepts from the scaling theory of disordered systems to describe the nature of the phase transition of the disordered model, departing gradually from the fixed point of the pure system. Our analysis (based on the determination of the critical exponents) shows that the 3D random-bond Ising model belongs to the same universality class with the site- and bond-dilution models, providing a single universality class for the 3D Ising model with these three types of quenched uncorrelated disorder.
We investigate, by means of extensive Monte Carlo simulations, the magnetic critical behavior of the three-dimensional bimodal random-field Ising model at the strong disorder regime. We present results in favor of the two-exponent scaling scenario, $
bar{eta}=2eta$, where $eta$ and $bar{eta}$ are the critical exponents describing the power-law decay of the connected and disconnected correlation functions and we illustrate, using various finite-size measures and properly defined noise to signal ratios, the strong violation of self-averaging of the model in the ordered phase.
