Critical dynamics of cluster algorithms in the random-bond Ising model.

In the present work, we present an extensive Monte Carlo simulation study on the dynamical properties of the two-dimensional random-bond Ising model. The correlation time τ of the Swendsen-Wang and Wolff cluster algorithms is calculated at the critical point. The dynamic critical exponent z of both algorithms is also measured by using the numerical data for several lattice sizes up to L=512. It is found for both algorithms that the autocorrelation time decreases considerably and the critical slowing-down effect reduces upon the introduction of bond disorder. Additionally, simulations with the Metropolis algorithm are performed, and the critical slowing-down effect is observed to be more pronounced in the presence of disorder, confirming the previous findings in the literature. Moreover, the existence of the non-self-averaging property of the model is demonstrated by calculating the scaled form of the standard deviation of autocorrelation times. Finally, the critical exponent ratio of the magnetic susceptibility is estimated by using the average cluster size of the Wolff algorithm.

Monte Carlo study of the two-dimensional kinetic Blume-Capel model in a quenched random crystal field.

We investigate by means of Monte Carlo simulations the dynamic phase transition of the two-dimensional kinetic Blume-Capel model under a periodically oscillating magnetic field in the presence of a quenched random crystal-field coupling. We analyze the universality principles of this dynamic transition for various values of the crystal-field coupling at the originally second-order regime of the corresponding equilibrium phase diagram of the model. A detailed finite-size scaling analysis indicates that the observed nonequilibrium phase transition belongs to the universality class of the equilibrium Ising ferromagnet with additional logarithmic corrections in the scaling behavior of the heat capacity. Our results are in agreement with earlier works on kinetic Ising models.

Ferromagnetic TM2BC (TM = Cr, Mn) monolayers for spintronic devices with high Curie temperature.

Transition metal boro-carbide (TM2BC) structures crystallize in the layered orthorhombic structure in their bulk phases. In this study, however, we find that TM2BC (TM = Cr, Mn) prefer a tetragonal (t) crystal structure in their monolayer phases due to the occurrence of strong sp2 bonds between the metal and B/C atoms instead of sp3 + sp2 bonds which exist in the orthorhombic phase. The calculated energy difference between the orthorhombic and the tetragonal structures based on density functional theory (DFT) is more than 1 eV per unit cell. In addition, t-Cr2BC and t-Mn2BC monolayers are dynamically and thermally stable with their magnetic metal electronic structures. For further investigations, we combine our DFT calculations with the Monte Carlo simulations and find that both t-TM2BC monolayers show ferromagnetic properties. The calculated Curie temperatures are 846 K and 128 K for t-Cr2BC and t-MnBC sheets, respectively. In addition, we examine the magnetic anisotropy energies (MAE) of t-TM2BC monolayers and find that both structures prefer out-of-plane as the easy axis magnetization direction and the applied electric field can easily modulate the MAE of the monolayers. Our theoretical calculations reveal that t-TM2BC (TM = Cr, Mn) sheets have great potential for the future design of controllable spintronic devices with their tunable MAE properties.

Ising universality in the two-dimensional Blume-Capel model with quenched random crystal field.

Using high-precision Monte Carlo simulations based on a parallel version of the Wang-Landau algorithm and finite-size scaling techniques, we study the effect of quenched disorder in the crystal-field coupling of the Blume-Capel model on a square lattice. We mainly focus on the part of the phase diagram where the pure model undergoes a continuous transition, known to fall into the universality class of a pure Ising ferromagnet. A dedicated scaling analysis reveals concrete evidence in favor of the strong universality hypothesis with the presence of additional logarithmic corrections in the scaling of the specific heat. Our results are in agreement with an early real-space renormalization-group study of the model as well as a very recent numerical work where quenched randomness was introduced in the energy exchange coupling. Finally, by properly fine tuning the control parameters of the randomness distribution we also qualitatively investigate the part of the phase diagram where the pure model undergoes a first-order phase transition. For this region, preliminary evidence indicate a smoothing of the transition to second-order with the presence of strong scaling corrections.