We study the pure and random-bond versions of the square lattice ferromagnetic Blume-Capel model, in both the first-order and second-order phase transition regimes of the pure model. Phase transition temperatures, thermal and magnetic critical exponents are determined for lattice sizes in the range L = 20-100 via a sophisticated two-stage numerical strategy of entropic sampling in dominant energy subspaces, using mainly the Wang-Landau algorithm. The second-order phase transition, emerging under random bonds from the secondorder regime of the pure model, has the same values of critical exponents as the two-dimensional Ising universality class, with the effect of the bond disorder on the specific heat being well described by doublelogarithmic corrections, our findings thus supporting the marginal irrelevance of quenched bond randomness. On the other hand, the second-order transition, emerging under bond randomness from the first-order regime of the pure model, has a distinctive universality class with = 1.30͑6͒ and ␤ / = 0.128͑5͒. These results amount to a strong violation of universality principle of critical phenomena, since these two second-order transitions, with different sets of critical exponents, are between the same ferromagnetic and paramagnetic phases. Furthermore, the latter of these two sets of results supports an extensive but weak universality, since it has the same magnetic critical exponent ͑but a different thermal critical exponent͒ as a wide variety of two-dimensional systems with and without quenched disorder. In the conversion by bond randomness of the first-order transition of the pure system to second order, we detect, by introducing and evaluating connectivity spin densities, a microsegregation that also explains the increase we find in the phase transition temperature under bond randomness.

Dominant energy subspaces of statistical systems are defined with the help of restrictive conditions on various characteristics of the energy distribution, such as the probability density and the fourth order Binder's cumulant. Our analysis generalizes the ideas of the critical minimum energy subspace (CRMES) technique, applied previously to study the specific heat's finite-size scaling. Here, we illustrate alternatives that are useful for the analysis of further finite-size anomalies and the behavior of the corresponding dominant subspaces is presented for the two-dimensional (2D) Baxter-Wu and the 2D and 3D Ising models. In order to show that a CRMES technique is adequate for the study of magnetic anomalies, we study and test simple methods which provide the means for an accurate determination of the energy-order-parameter (E,M) histograms via Wang-Landau random walks. The 2D Ising model is used as a test case and it is shown that high-level Wang-Landau sampling schemes yield excellent estimates for all magnetic properties. Our estimates compare very well with those of the traditional Metropolis method. The relevant dominant energy subspaces and dominant magnetization subspaces scale as expected with exponents alpha/nu and gamma/nu, respectively. Using the Metropolis method we examine the time evolution of the corresponding dominant magnetization subspaces and we uncover the reasons behind the inadequacy of the Metropolis method to produce a reliable estimation scheme for the tail regime of the order-parameter distribution.

We investigate and contrast, via entropic sampling based on the Wang-Landau algorithm, the effects of quenched bond randomness on the critical behavior of two Ising spin models in 2D. The random bond version of the superantiferromagnetic (SAF) square model with nearest-and next-nearest-neighbor competing interactions and the corresponding version of the simple Ising model are studied and their general universality aspects are inspected by a detailed finite-size scaling (FSS) analysis. We find that, the random bond SAF model obeys weak universality, hyperscaling, and exhibits a strong saturating behavior of the specific heat due to the competing nature of interactions. On the other hand, for the random Ising model we encounter some difficulties for a definite discrimination between the two well-known scenarios of the logarithmic corrections versus the weak universality. Yet, a careful FSS analysis of our data favors the field-theoretically predicted logarithmic corrections.

The effects of bond randomness on the phase diagram and critical behavior of the square lattice ferromagnetic Blume-Capel model are discussed. The system is studied in both the pure and disordered versions by the same efficient two-stage Wang-Landau method for many values of the crystal field, restricted here in the second-order phase transition regime of the pure model. For the random-bond version several disorder strengths are considered. We present phase diagram points of both pure and random versions and for a particular disorder strength we locate the emergence of the enhancement of ferromagnetic order observed in an earlier study in the ex-first-order regime. The critical properties of the pure model are contrasted and compared to those of the random model. Accepting, for the weak random version, the assumption of the double logarithmic scenario for the specific heat we attempt to estimate the range of universality between the pure and random-bond models. The behavior of the strong disorder regime is also discussed and a rather complex and yet not fully understood behavior is observed. It is pointed out that this complexity is related to the ground-state structure of the random-bond version.

We have applied Sullivan's model to study the interfacial properties of a spherical drop embedded in a one-component fluid using the mean field approximation. We have examined the effect of drop size (expressed by the equimolar radius Re) on the density profile, pressure tensor, pressure difference Delta p across the interface, surface tension, Tolman's length, and density, pressure, and normal component of pressure tensor at the centre of the drop. Delta p is found to be an ill defined quantity in the sense that it can be defined in various ways, whose results coincide for large drops. The surface tension is a non-monotonic function of Re; it increases slowly from its flat interface value as Re decreases until a maximum value is attained, then it decreases rapidly. For small supersaturations, the drops are under tension and compression while at large supersaturations they are only under tension. The results of this theory are compared qualitatively with previous molecular dynamics simulations and theoretical calculations for drops.

No abstract

No abstract

We present a variation of the N-fold way algorithm, which improves efficiency when one combines the Wang–Landau method with the N-fold way. It is shown that the new version of the N-fold way algorithm has good performance when used for importance sampling and compared with the usual N-fold version or the Metropolis algorithm. The new N-fold algorithm combined with the Wang–Landau method is applied to the square Ising model using a multi-range approach. A comparative study is presented for all these algorithms, Wang–Landau and the two combined versions with the N-fold way. The role of boundary effects is discussed.