A microcanonical Monte Carlo algorithm has been developed to calculate the density of states of a two-dimensional Blume-Capel model in zero field. The full density of states is calculated for a 16*16 lattice as a function of the two spin summations which appear in the Hamiltonian. This permits the partition function and related thermodynamic functions to be evaluated for any temperature and any single site anisotropy from a single data set. Results are presented for the pseudo-critical transition temperature and the maximum in the specific heat as a function of the single site anisotropy parameter. The technique will allow the region near the tricritical point to be explored in detail and will permit the future determination of the zeros of the partition function. It should be noted that the microcanonical sampling method of Lee (1990) was found to be unusably slow for the S=1 problem in regions were the density of states had a significant slope and the microcanonical sampling described in this work is a significant improvement on his sampling method.
A microcanonical Monte Carlo algorithm has been developed to calculate the density of states of a two-dimensional Blume-Capel model in zero field. The full density of states is calculated for a 16*16 lattice as a function of the two spin summations which appear in the Hamiltonian. This permits the partition function and related thermodynamic functions to be evaluated for any temperature and any single site anisotropy from a single data set. Results are presented for the pseudo-critical transition temperature and the maximum in the specific heat as a function of the single site anisotropy parameter. The technique will allow the region near the tricritical point to be explored in detail and will permit the future determination of the zeros of the partition function. It should be noted that the microcanonical sampling method of Lee (1990) was found to be unusably slow for the S=1 problem in regions were the density of states had a significant slope and the microcanonical sampling described in this work is a significant improvement on his sampling method.
The microcanonical partition function ${\ensuremath{\Omega}}_{\ensuremath{\alpha},N}(E)$ of a three-dimensional ferromagnetic simple cubic Blume-Capel model is calculated for several coupling ratios $\ensuremath{\alpha}$ close to and on the first-order side of the tricritical point. To this end a single Monte Carlo simulation of a suitably extended partition function for systems of $L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L$ spins with $L\ensuremath{\in}{8,10,12,14,18}$ was performed. A finite system analog of the latent heat is introduced in order to define a tricritical point in finite systems as well. An empirical scaling of its coordinates to the thermodynamic limit yields ${\ensuremath{\alpha}}^{(t)}=2.84479\ifmmode\pm\else\textpm\fi{}0.00030$ for the tricritical coupling ratio, and ${k}_{B}{T}^{(t)}/J=1.4182\ifmmode\pm\else\textpm\fi{}0.0055$ for the tricritical transition temperature. The results are compared with values obtained on fcc lattices.
We present a modification of the Rose-Machta algorithm [N. Rose and J. Machta, Phys. Rev. E 100, 063304 (2019)2470-004510.1103/PhysRevE.100.063304] and estimate the density of states for a two-dimensional Blume-Capel model, simulating 10^{5} replicas in parallel for each set of parameters. We perform a finite-size analysis of the specific heat and Binder cumulant, determine the critical temperature along the critical line, and evaluate the critical exponents. The obtained results are in good agreement with those previously obtained using various methods-Markov chain Monte Carlo simulation, Wang-Landau simulation, transfer matrix, and series expansion. The simulation results clearly illustrate the typical behavior of specific heat along the critical lines and through the tricritical point.
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.
We investigate the tricritical scaling behavior of the two-dimensional spin-1 Blume-Capel model by using the Wang-Landau method of measuring the joint density of states for lattice sizes up to 48×48 sites. We find that the specific heat deep in the first-order area of the phase diagram exhibits a double-peak structure of the Schottky-like anomaly appearing with the transition peak. The first-order transition curve is systematically determined by employing the method of field mixing in conjunction with finite-size scaling, showing a significant deviation from the previous data points. At the tricritical point, we characterize the tricritical exponents through finite-size-scaling analysis including the phenomenological finite-size scaling with thermodynamic variables. Our estimation of the tricritical eigenvalue exponents, yt=1.804(5), yg=0.80(1), and yh=1.925(3), provides the first Wang-Landau verification of the conjectured exact values, demonstrating the effectiveness of the density-of-states-based approach in finite-size scaling study of multicritical phenomena.
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 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.