Abstract:We present a modification of the Rose-Machta algorithm (Phys. Rev. E 100 (2019) 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 results obtained are in good agreement with those obtained previously 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.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is, in the two - dimensional Blume - Capel model, through the improved Microcanonical Population Annealing (MCPA) method, to study the finite - size behaviors of the Density of States (DoS), Specific Heat ($C$) and Binder Cumulant ($U_4$) of the system, and to determine the critical temperature and critical exponents. Specifically, the author focuses on:
1. **Improved algorithm**: The paper proposes an improved Rose - Machta algorithm. By introducing the Ceiling Energy ($E_c$) and Floor Energy ($E_f$), it estimates the Density of States of the system. This method can cover the entire energy spectrum more accurately, thereby improving the estimation accuracy of the system state.
2. **Analysis of critical behavior**: Through the finite - size analysis of Specific Heat and Binder Cumulant, the critical behavior of the system on the second - order phase transition line, tricritical point and first - order phase transition line is studied. In particular, the author focuses on how the maximum value of Specific Heat changes with the system size and how these changes reflect different phase transition types.
3. **Critical temperature and critical exponents**: Determine the critical temperature $T_c$ and critical exponents, especially the critical exponent $\alpha$ of Specific Heat and the critical exponent $\nu$ of the correlation length. The author estimates the critical temperature by the position of the maximum value of Specific Heat and the position of the minimum value of Binder Cumulant, and extracts the critical exponents by fitting these data.
4. **Comparison with existing results**: Verify the accuracy of the new method. By comparing the results with those obtained by previous different methods (such as Markov Chain Monte Carlo simulation, Wang - Landau simulation, transfer matrix method and series expansion, etc.), ensure the effectiveness and reliability of the new method.
### Main research contents
- **Model introduction**: The Blume - Capel model is a lattice model containing spins $\sigma_i$. The spin can take three values: - 1, 0, 1. The Hamiltonian of the model is:
\[
H=-J\sum_{\langle i, j\rangle}\sigma_i\sigma_j+\Delta\sum_i\sigma_i^2
\]
where $J$ is the ferromagnetic coupling constant and $\Delta$ is the crystal field.
- **Microcanonical population annealing algorithm**: The improved MCPA algorithm generates copies of the initial random configuration and gradually adjusts the energy upper limit and lower limit to cover the entire energy spectrum. The key steps of the algorithm include generating the initial state, performing MCMC steps, adjusting the energy upper limit and lower limit, calculating the Culling Fraction, etc.
- **Finite - size analysis of Specific Heat and Binder Cumulant**: By calculating the finite - size behaviors of Specific Heat ($C$) and Binder Cumulant ($U_4$), the critical behavior of the system on different phase transition lines is studied. Special attention is paid to how the maximum value of Specific Heat $C_{\text{max}}$ changes with the system size $L$ and how these changes reflect different phase transition types.
- **Estimation of critical temperature and critical exponents**: Estimate the critical temperature $T_c$ by the position of the maximum value of Specific Heat and the position of the minimum value of Binder Cumulant, and extract the critical exponents by fitting these data. Especially near the tricritical point, study how the maximum value of Specific Heat changes with the system size to verify whether it conforms to the expected power - law behavior.
### Conclusion
The paper successfully estimates the Density of States, Specific Heat and Binder Cumulant of the two - dimensional Blume - Capel model through the improved MCPA algorithm, and determines the critical temperature and critical exponents. The research results are compared with those obtained by previous different methods.