The dynamic phase transitions of periodically driven ferromagnets present a complex interplay between internal relaxation dynamics and external driving forces. In this study, a highly optimized Metropolis-Hastings Monte Carlo simulation is developed to model a two-dimensional Ising ferromagnet subjected to an oscillating magnetic field. Following exact thermodynamic validation of the static critical exponent of 1.75, the system is pushed out of equilibrium to investigate two competing mechanisms of phase disruption: structural spatial disorder and temporal frequency variation. Computational results demonstrate that quenched non-magnetic impurities act as localized pinning centers for expanding domain walls, fundamentally "softening" the material and drastically shrinking the dynamically ordered parameter space. Conversely, scaling the driving frequency outpaces localized spin relaxation, effectively bypassing the chaotic limit cycle and expanding the dynamically ordered phase. This dichotomy offers a comprehensive mathematical profile of non-equilibrium phase boundaries in disordered magnetic lattices.
The dynamic phase transitions of periodically driven ferromagnets present a complex interplay between internal relaxation dynamics and external driving forces. In this study, a highly optimized Metropolis-Hastings Monte Carlo simulation is developed to model a two-dimensional Ising ferromagnet subjected to an oscillating magnetic field. Following exact thermodynamic validation of the static critical exponent of 1.75, the system is pushed out of equilibrium to investigate two competing mechanisms of phase disruption: structural spatial disorder and temporal frequency variation. Computational results demonstrate that quenched non-magnetic impurities act as localized pinning centers for expanding domain walls, fundamentally "softening" the material and drastically shrinking the dynamically ordered parameter space. Conversely, scaling the driving frequency outpaces localized spin relaxation, effectively bypassing the chaotic limit cycle and expanding the dynamically ordered phase. This dichotomy offers a comprehensive mathematical profile of non-equilibrium phase boundaries in disordered magnetic lattices.
This study presents a comprehensive Monte Carlo simulation of the two-dimensional (2-D) Ising model using the Metropolis algorithm to investigate critical phenomena and thermodynamic behavior in spin-lattice systems. The model, implemented in MATLAB with periodic boundary conditions, explores equilibrium properties such as magnetization, internal energy, specific heat, and susceptibility across a range of temperatures. By employing various initial spin configurations—ordered and random—the simulations demonstrate the system's ergodicity and convergence to thermal equilibrium. Key results include a sharp decline in magnetization and a pronounced peak in specific heat near the critical temperature, consistent with second-order phase transition behavior. The simulation captures microscopic domain evolution, highlighting the transition from ferromagnetic to paramagnetic phases as thermal fluctuations increase. The study further evaluates algorithmic efficiency, discusses the impact of lattice size on statistical accuracy, and proposes improvements using advanced cluster algorithms and parallel computing frameworks. The findings validate theoretical predictions from Onsager's solution and underscore the versatility of Monte Carlo techniques in modeling collective behavior in magnetic systems. The simulation framework offers a robust foundation for analyzing critical dynamics and extends its relevance to broader applications in material science, computational physics, and complex systems modeling.
Among the many non-equilibrium processes encountered in nature we deal with two different but related aspects. One is the non-equilibrium relaxation process that is at the origin of \'aging phenomena••, and the other one is a non-equilibrium phase transition, called ••dynamic phase transition••. One of the main purposes of our research is to explore more realistic situations than studied previously. Indeed, in the study of aging phenomena certain kinds of disorder effects are considered, and we introduce the ••surface•• as a spatial boundary to the system undergoing the dynamic phase transition. In order to observe these processes as clearly as possible, we study in both cases simple spin systems. Using Monte Carlo simulations we first investigate aging in three-dimensional Ising spin glasses as well as in two-dimensional Ising models with disorder quenched to low temperatures. The time-dependent dynamical correlation length L(t) is determined numerically and the scaling behavior of various two-time quantities as a function of L(t)/L(s) is discussed where t and s are two different times. For disordered Ising models deviations of L(t) from algebraic growth law show up. The generalized scaling forms as a function of L(t)/L(s) reveal a generic simple aging scenario for Ising spin glasses as well as for disordered Ising ferromagnets. We also study the local critical phenomena at a dynamic phase transition by means of numerical simulations of kinetic Ising models with surfaces subjected to a periodic oscillating field. We examine layer-dependent quantities, such as the period-averaged magnetization per layer Q(z) and the layer susceptibility ¥ö(z), and determine local critical exponents through finite size scaling. Both for two and three dimensions, we find that the values of the surface exponents differ from those of the equilibrium critical surface. It is revealed that the surface phase diagram of the non-equilibrium system is not identical to that of the equilibrium system in three dimensions.
To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.
We study the critical dynamics arising from a time-dependent periodic homogenous source coupled to the order-parameter field, which drives a classical ferromagnetic system across a continuous transition. For this purpose, we consider the paradigmatic two-dimensional (2D) Ising model in the presence of a periodic magnetic field h(t)=−Acos(2πt/P), evolving under a purely relaxational dynamics at the critical temperature. We show that the periodic driving gives rise to a peculiar dynamic scaling behavior in the thermodynamic limit, arising from a nontrivial interplay among the time t, the amplitude A and period P of h(t). The relevant scaling variables are τ=t/P and σ=APκ, with κ=yh/z, where yh=(d+2−η)/2 is the critical dimension of the magnetic field, and z is dynamic exponent for the critical relaxational dynamics (κ≈0.865 for the 2D Ising model). The dynamic scaling behaviors of the magnetization and bond-energy density show an oscillatory behavior around a smooth curve which approaches a large-τ stationary behavior. We also briefly discuss the dynamic behavior of an Ising system driven across the critical point by a periodic time-varying temperature at zero magnetic field.
The two-dimensional ferromagnetic Ising model in the presence of a propagating magnetic field wave (with a well-defined frequency and wavelength) was studied by Monte Carlo simulation. This study differs from similar earlier studies, because in earlier studies the oscillating magnetic field was considered to be uniform in space. The time averaged magnetization over a full cycle (the time period) of the propagating magnetic field acted as the dynamic order parameter. The dynamical phase transition was observed. The temperature variation of the dynamic order parameter, the mean square deviation of the dynamic order parameter, the dynamic specific heat and the derivative of the dynamic order parameter were studied. The mean square deviation of the dynamic order parameter and the dynamic specific heat show sharp maxima near the transition point. The derivative of the dynamic order parameter shows a sharp minimum near the transition point. The transition temperature is found to depend also on the speed of propagation of the magnetic field wave.
Many problems in statistical physics resist exact solution once the number of interacting particles becomes large, and the two-dimensional Ising model is a classic case where numerical simulation fills that gap. This paper reports a Monte Carlo study of the two-dimensional Ising model using the Metropolis algorithm, undertaken to examine how magnetization, internal energy, specific heat, and magnetic susceptibility change with temperature. Spin configurations were generated on a square lattice under periodic boundary conditions, with each proposed spin flip accepted or rejected according to the Boltzmann weight of the resulting energy change. The simulation reproduces the expected transition from an ordered ferromagnetic state at low temperature to a disordered paramagnetic state at high temperature, with magnetization falling sharply and both specific heat and susceptibility peaking near a critical temperature of Tc ≈ 2.3–2.4 (in units of J/kB), consistent with the exact value obtained analytically for this system. These results confirm that the Metropolis algorithm, despite its simplicity, remains a dependable tool for locating second-order phase transitions and extracting thermodynamic quantities in lattice spin systems, and the same computational framework can be adapted to more elaborate models used in condensed matter and materials research.