First-order phase transition and tricritical scaling behavior of the Blume-Capel model: A Wang-Landau sampling approach

Numerical simulation of a two-dimensional Blume-Capel ferromagnet in an oscillating magnetic field with a constant bias

We perform a numerical study of the kinetic Blume-Capel (BC) model to find if it exhibits the metamagnetic anomalies previously observed in the kinetic Ising model for supercritical periods [P. Riego et al., Phys. Rev. Lett. 118, 117202 (2017)0031-900710.1103/PhysRevLett.118.117202; G. M. Buendía et al., Phys. Rev. B 96, 134306 (2017)2469-995010.1103/PhysRevB.96.134306]. We employ a heat-bath Monte Carlo (MC) algorithm on a square lattice in which spins can take values of ±1,0, with a nonzero crystal field, subjected to a sinusoidal oscillating field in conjunction with a constant bias. In the ordered region, we find an equivalent hysteretic response of the order parameters with its respective conjugate fields between the kinetic and the equilibrium model. In the disordered region (supercritical periods), we observed two peaks, symmetrical with respect to zero bias, in the susceptibility and scaled variance curves, consistent with the numerical and experimental findings on the kinetic Ising model. This behavior does not have a counterpart in the equilibrium model. Furthermore, we find that the peaks occur at higher values of the bias field and become progressively smaller as the density of zeros, or the amplitude of the oscillating field, increases. Using nucleation theory, we demonstrate that these fluctuations, as in the Ising model, are not a critical phenomenon, but that they are associated with a crossover between a single-droplet (SD) and a multidroplet (MD) magnetization switching mechanism. For strong (weak) bias, the SD (MD) mechanism dominates. We also found that the zeros concentrate on the droplets' surfaces, which may cause a reduced interface tension in comparison with the Ising model [M. Schick et al., Phys. Rev. B 34, 1797 (1986)0163-182910.1103/PhysRevB.34.1797]. Our results suggest that metamagnetic anomalies are not particular to the kinetic Ising model, but rather are a general characteristic of spin kinetic models, and provide further evidence that the equivalence between dynamical phase transitions and equilibrium ones is only valid near the critical point.

Blume-Capel model analysis with a microcanonical population annealing method

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.

Two-dimensional dilute Baxter-Wu model: Transition order and universality

We investigate the critical behavior of the two-dimensional spin-1 Baxter-Wu model in the presence of a crystal-field coupling Δ with the goal of determining the universality class of transitions along the second-order part of the transition line as one approaches the putative location of the multicritical point. We employ extensive Monte Carlo simulations using two different methodologies: (i) a study of the zeros of the energy probability distribution, closely related to the Fisher zeros of the partition function, and (ii) the well-established multicanonical approach employed to study the probability distribution of the crystal-field energy. A detailed finite-size scaling analysis in the regime of second-order phase transitions in the (Δ,T) phase diagram supports previous claims that the transition belongs to the universality class of the four-state Potts model. For positive values of Δ, we observe the presence of strong finite-size effects, indicative of crossover effects due to the proximity of the first-order part of the transition line. Finally, we demonstrate how a combination of cluster and heat-bath updates allows one to equilibrate larger systems, and we demonstrate the potential of this approach for resolving the ambiguities observed in the regime of Δ≳0.

Universality in the two-dimensional dilute Baxter-Wu model

We study the question of universality in the two-dimensional spin-1 Baxter-Wu model in the presence of a crystal field Δ. We employ extensive numerical simulations of two types, providing us with complementary results: Wang-Landau sampling at fixed values of Δ and a parallelized variant of the multicanonical approach performed at constant temperature T. A detailed finite-size scaling analysis in the regime of second-order phase transitions in the (Δ,T) phase diagram indicates that the transition belongs to the universality class of the four-state Potts model. Previous controversies with respect to the nature of the transition are discussed and attributed to the presence of strong finite-size effects, especially as one approaches the pentacritical point of the model.

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.

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.

Tricritical point in the mixed-spin Blume-Capel model on three-dimensional lattices: Metropolis and Wang-Landau sampling approaches

We investigate the mixed-spin Blume-Capel model with spin-1/2 and spin-S (S=1, 2, and 3) on the simple cubic and body-centered cubic lattices with single-ion-splitting crystal field (Δ) by using the Metropolis and the Wang-Landau Monte Carlo methods. We show that the two methods are complementary: The Wang-Landau algorithm is efficient to construct phase diagrams and the Metropolis algorithm allows access to large-sized lattices. By numerical simulations, we prove that the tricritical point is independent of S for both lattices. The positions of the tricritical point in the phase diagram are determined as [Δ_{t}/J=2.978(1); k_{B}T_{t}/J=0.439(1)] and [Δ_{t}/J=3.949(1); k_{B}T_{t}/J=0.854(1)] for the simple cubic and the body-centered cubic lattices, respectively. A very strong supercritical slowing down and hysteresis were observed in the Metropolis update close to first-order transitions for Δ>Δ_{t} in the body-centered cubic lattice. In addition, for both lattices we found a line of compensation points, where the two sublattice magnetizations have the same magnitude. We show that the compensation lines are also S independent.

Revisiting the field-driven edge transition of the tricritical two-dimensional Blume-Capel model

We reconsider the tricritical Blume-Capel model on the square lattice with a magnetic field acting on the open boundaries in one direction. Periodic boundary conditions are applied in the other direction. We apply three types of Monte Carlo algorithms, local Metropolis updates, and cluster algorithms of the Wolff and geometric type, adapted to the symmetry properties of the model. Statistical analyses of the bulk magnetization, the bulk Binder ratio, the edge magnetization, and the connected product of the edge and bulk magnetizations lead to new results confirming the presence of a singular edge transition at H_{sc}≈0.68, as we reported earlier [Phys. Rev. E 71, 026109 (2005)PLEEE81539-375510.1103/PhysRevE.71.026109]. We provide a plausible answer concerning a discrepancy between the behavior of the edge Binder ratio reported in that work and our new results.

Monte Carlo study of the interfacial adsorption of the Blume-Capel model

We investigate the scaling of the interfacial adsorption of the two-dimensional Blume-Capel model using Monte Carlo simulations. In particular, we study the finite-size scaling behavior of the interfacial adsorption of the pure model at both its first- and second-order transition regimes, as well as at the vicinity of the tricritical point. Our analysis benefits from the currently existing quite accurate estimates of the relevant (tri)critical-point locations. In all studied cases, the numerical results verify to a level of high accuracy the expected scenarios derived from analytic free-energy scaling arguments. We also investigate the size dependence of the interfacial adsorption under the presence of quenched bond randomness at the originally first-order transition regime (disorder-induced continuous transition) and the relevant self-averaging properties of the system. For this ex-first-order regime, where strong transient effects are shown to be present, our findings support the scenario of a non-divergent scaling, similar to that found in the original second-order transition regime of the pure model.

Smallest neural network to learn the Ising criticality

Learning with an artificial neural network encodes the system behavior in a feed-forward function with a number of parameters optimized by data-driven training. An open question is whether one can minimize the network complexity without loss of performance to reveal how and why it works. Here we investigate the learning of the phase transition in the Ising model and find that having two hidden neurons can be enough for an accurate prediction of critical temperature. We show that the networks learn the scaling dimension of the order parameter while being trained as a phase classifier, demonstrating how the machine learning exploits the Ising universality to work for different lattices of the same criticality within a single set of trainings in one lattice geometry.

Phase diagram of the J_{1}-J_{2} frustrated anisotropic antiferromagnet with spin S=1 on the quadratic lattice

In the paper the phase diagram of J_{1}-J_{2} frustrated antiferromagnet with spin S=1 and single-ion anisotropy is studied on the planar quadratic lattice in the cluster approximation. The Bogolyubov inequality is adopted for the Gibbs energy calculation for the case of 2×2 and 4×4 clusters. On this basis, the ranges of existence of the antiferromagnetic, superantiferromagnetic, and paramagnetic phases are investigated for the antiferromagnetic nearest-neighbor (J_{1}<0) and next-nearest-neighbor (J_{2}<0) interactions. In particular, the occurrence of tricritical and triple points is discussed and a comparison between the results for 2×2 and 4×4 clusters is made. The results are also compared with the classical MFA method, adopted here for the model in question, as well as with selected literature results for particular choices of interaction parameters.

Dynamic phase transitions in the presence of quenched randomness

We present an extensive study of the effects of quenched disorder on the dynamic phase transitions of kinetic spin models in two dimensions. We undertake a numerical experiment performing Monte Carlo simulations of the square-lattice random-bond Ising and Blume-Capel models under a periodically oscillating magnetic field. For the case of the Blume-Capel model we analyze the universality principles of the dynamic disordered-induced continuous transition at the low-temperature regime of the phase diagram. A detailed finite-size scaling analysis indicates that both nonequilibrium phase transitions belong to the universality class of the corresponding equilibrium random Ising model.

Dynamic phase transition of the Blume-Capel model in an oscillating magnetic field

We employ numerical simulations and finite-size scaling techniques to investigate the properties of the dynamic phase transition that is encountered in the Blume-Capel model subjected to a periodically oscillating magnetic field. We mainly focus on the study of the two-dimensional system for various values of the crystal-field coupling in the second-order transition regime. Our results indicate that the present nonequilibrium phase transition belongs to the universality class of the equilibrium Ising model and allow us to construct a dynamic phase diagram, in analogy with the equilibrium case, at least for the range of parameters considered. Finally, we present some complementary results for the three-dimensional model, where again the obtained estimates for the critical exponents fall into the universality class of the corresponding three-dimensional equilibrium Ising ferromagnet.

Universality from disorder in the random-bond Blume-Capel model

Using high-precision Monte Carlo simulations and finite-size scaling we study the effect of quenched disorder in the exchange couplings on the Blume-Capel model on the square lattice. The first-order transition for large crystal-field coupling is softened to become continuous, with a divergent correlation length. An analysis of the scaling of the correlation length as well as the susceptibility and specific heat reveals that it belongs to the universality class of the Ising model with additional logarithmic corrections which is also observed for the Ising model itself if coupled to weak disorder. While the leading scaling behavior of the disordered system is therefore identical between the second-order and first-order segments of the phase diagram of the pure model, the finite-size scaling in the ex-first-order regime is affected by strong transient effects with a crossover length scale L^{*}≈32 for the chosen parameters.

Computation of a Theoretical Membrane Phase Diagram and the Role of Phase in Lipid-Raft-Mediated Protein Organization

Lipid phase heterogeneity in the plasma membrane is thought to be crucial for many aspects of cell signaling, but the physical basis of participating membrane domains such as "lipid rafts" remains controversial. Here we consider a lattice model yielding a phase diagram that includes several states proposed to be relevant for the cell membrane, including microemulsion-which can be related to membrane curvature-and Ising critical behavior. Using a neural-network-based machine learning approach, we compute the full phase diagram of this lattice model. We analyze selected regions of this phase diagram in the context of a signaling initiation event in mast cells: recruitment of the membrane-anchored tyrosine kinase Lyn to a cluster of transmembrane IgE-FcεRI receptors. We find that model membrane systems in microemulsion and Ising critical states can mediate roughly equal levels of kinase recruitment (binding energy ∼ -0.6 k_B T), whereas a membrane near a tricritical point can mediate a much stronger kinase recruitment (-1.7 k_B T). By comparing several models for lipid heterogeneity within a single theoretical framework, this work points to testable differences between existing models. We also suggest the tricritical point as a new possibility for the basis of membrane domains that facilitate preferential partitioning of signaling components.

Partition function zeros of the p-state clock model in the complex temperature plane

We investigate the partition function zeros of the two-dimensional p-state clock model in the complex temperature plane by using the Wang-Landau method. For p=5, 6, 8, and 10, we propose a modified energy representation to enumerate exact irregular energy levels for the density of states without any binning artifacts. Comparing the leading zeros between different p's, we provide strong evidence that the upper transition at p=6 is indeed of the Berezinskii-Kosterlitz-Thouless (BKT) type in contrast to the claim of the previous Fisher zero study [Phys. Rev. E 80, 042103 (2009)10.1103/PhysRevE.80.042103]. We find that the leading zeros of p=6 at the upper transition collapse onto the zero trajectories of the larger p's including the XY limit while the finite-size behavior of p=5 differs from the converged behavior of p≥6 within the system sizes examined. In addition, we argue that the nondivergent specific heat in the BKT transition is responsible for the small partition function magnitude that decreases exponentially with increasing system size near the leading zero, fundamentally limiting access to large systems in search for zeros with an estimator under finite statistical fluctuations.

Discovering phases, phase transitions, and crossovers through unsupervised machine learning: A critical examination

We apply unsupervised machine learning techniques, mainly principal component analysis (PCA), to compare and contrast the phase behavior and phase transitions in several classical spin models-the square- and triangular-lattice Ising models, the Blume-Capel model, a highly degenerate biquadratic-exchange spin-1 Ising (BSI) model, and the two-dimensional XY model-and we examine critically what machine learning is teaching us. We find that quantified principal components from PCA not only allow the exploration of different phases and symmetry-breaking, but they can distinguish phase-transition types and locate critical points. We show that the corresponding weight vectors have a clear physical interpretation, which is particularly interesting in the frustrated models such as the triangular antiferromagnet, where they can point to incipient orders. Unlike the other well-studied models, the properties of the BSI model are less well known. Using both PCA and conventional Monte Carlo analysis, we demonstrate that the BSI model shows an absence of phase transition and macroscopic ground-state degeneracy. The failure to capture the "charge" correlations (vorticity) in the BSI model (XY model) from raw spin configurations points to some of the limitations of PCA. Finally, we employ a nonlinear unsupervised machine learning procedure, the "autoencoder method," and we demonstrate that it too can be trained to capture phase transitions and critical points.