Abstract We show that accurate insights into the critical properties of the Blume–Capel model at two dimensions can be deduced from Monte Carlo simulations, even for small system sizes, when one analyses the behaviour of the zeros of the partition function. The phase diagram of the model displays a line of second-order phase transitions ending at a tricritical point, then a line of first-order transitions. We concentrate on critical and tricritical properties and compare the accuracy achieved via standard finite-size scaling of thermodynamic quantities with that from the zeros analysis. This latter analysis showcases spectacular precision, even for systems as small as 64 spins. We also show that the zeros are very sensitive to subtle crossover effects.
Abstract We show that accurate insights into the critical properties of the Blume–Capel model at two dimensions can be deduced from Monte Carlo simulations, even for small system sizes, when one analyses the behaviour of the zeros of the partition function. The phase diagram of the model displays a line of second-order phase transitions ending at a tricritical point, then a line of first-order transitions. We concentrate on critical and tricritical properties and compare the accuracy achieved via standard finite-size scaling of thermodynamic quantities with that from the zeros analysis. This latter analysis showcases spectacular precision, even for systems as small as 64 spins. We also show that the zeros are very sensitive to subtle crossover effects.
We investigate the location of the critical and tricritical points of the three-dimensional Blume-Capel model by analyzing the behavior of the first Lee-Yang zero, the density of partition function zeros, and higher-order cumulants of the magnetization. Our analysis is conducted through Monte-Carlo simulations, intentionally using only small system sizes. We demonstrate that this approach yields excellent results for studying the critical behavior of the model. Our findings indicate that, at the tricritical point, where logarithmic corrections are anticipated, the numerical results align closely with the theoretical exponents describing these corrections. These expected values are then employed to accurately determine the coordinates of the tricritical point. At the model's critical point, the corrections correspond to those of the three-dimensional Ising model criticality, which we also use to precisely ascertain the critical temperature at zero crystal field. Additionally, we utilize more traditional thermodynamic quantities to validate the self-consistency of our analysis.
Abstract Since the landmark work of Lee and Yang, locating the zeros of the partition function in the complex magnetic-field plane has become a powerful method for studying phase transitions. Fisher later extended this approach to complex temperatures, and subsequent generalizations introduced other control parameters, such as the crystal field. In previous works (Moueddene et al 2024 J. Stat. Mech. 023206; 2024 Phys. Rev. E 110 064144) we applied this framework to the two- and three-dimensional Blume–Capel model–a system with a rich phase structure where a second-order critical line meets a first-order line at a tricritical point. We showed that the scaling of Lee-Yang, Fisher, and crystal-field zeros yields accurate critical exponents even for modest lattice sizes. In the present study, we extend this analysis and demonstrate that simulations need not be performed exactly at the nominal transition point to obtain reliable exponent estimates. Strikingly, small system sizes are sufficient, which not only improves methodological efficiency but also advances the broader goal of reducing the carbon footprint of large-scale computational studies.
We report on numerical simulations of the two-dimensional spin-1 Blume-Capel ferromagnet embedded in a triangular lattice. Utilizing a range of Monte Carlo and finite-size scaling techniques, we explore several critical aspects along the crystal field–temperature ( Δ , T ) transition line. Wang-Landau simulations measuring the joint density of states in combination with the method of field mixing allow us to probe the phase coexistence curve in high resolution, determining the tricritical point ( Δ t , T t ) with improved accuracy and verifying the tricritical exponents. Extensive multicanonical simulations identifying transition points across the phase diagram characterize the Ising universality class for Δ < Δ t with precise determination of thermal and magnetic critical exponents expected in the second-order regime. On the other hand, for Δ > Δ t , a finite-size scaling analysis is dedicated to revealing the first-order signature in the surface tension that linearly increases upon lowering the temperature deeper into the first-order transition regime. Finally, a comprehensive picture of the phase diagram for the model is presented, collecting transition points obtained from the combined numerical approach in this study and previous estimates in the literature.
The Blume-Capel Model as a ϕ6 Field Theory: A Monte Carlo and Flow Model Study
We report Monte Carlo simulations of a continuous spin variant of the Blume-Capel model. We compute the phase diagrams in the crystal-field-temperature plane, as a function of the softness of the spins. We show that normalizing flows, a class of generative models capable of direct sampling without autocorrelation, are effective samplers across the phase diagram, including in the first-order transition region and in the vicinity of the tricritical point. Finite-size scaling analyses of flow samples locate the first- and second-order transitions as well as the tricritical point with accuracy comparable to canonical and extended-ensemble Monte Carlo.
Abstract We perform Monte Carlo simulations, combining both the Wang–Landau and the Metropolis algorithms, to investigate the phase diagrams of the Blume–Capel model on different types of nonregular lattices (Lieb lattice (LL), decorated triangular lattice (DTL), and decorated simple cubic lattice (DSC)). The nonregular character of the lattices induces a double transition (reentrant behavior) in the region of the phase diagram at which the nature of the phase transition changes from first-order to second-order. A physical mechanism underlying this reentrance is proposed. The large-scale Monte Carlo simulations are performed with the finite-size scaling analysis to compute the critical exponents and the critical Binder cumulant for three different values of the anisotropy Δ / J ∈ 0 , 1 , 1.34 ( for LL ) , 1.51 ( for DTL and DSC ) , showing thus no deviation from the standard Ising universality class in two and three dimensions. We report also the location of the tricritical point to considerable precision: (Δ t / J = 1.3457(1); k B T t / J = 0.309(2)), (Δ t / J = 1.5766(1); k B T t / J = 0.481(2)), and (Δ t / J = 1.5933(1); k B T t / J = 0.569(4)) for LL, DTL, and DSC, respectively.
In this paper, we study finite-size effects in the Blume–Capel model through the analysis of the zeros of the partition function. We consider a complete graph and use the behavior of the partition function zeros to elucidate the crossover from effective to asymptotic properties. While in the thermodynamic limit the exact solution yields the asymptotic mean-field behavior, for finite size systems an effective critical behavior is observed. We show that even for large systems the criticality is not asymptotic. We also present insights into how partition function zeros in different complex fields (temperature, magnetic field, crystal field) give different precision and provide us with different parts of the larger picture. This includes the differences between criticality and tricriticality as seen through the lens of Fisher, Lee–Yang, and crystal field zeros.