Nonequilibrium kinetics of excess defect generation and dynamic scaling in the Ising spin chain under slow cooling

We investigate the relaxation dynamics of a one-dimensional Ising chain via Glauber kinetic Monte Carlo simulations, when the system is cooled slowly from infinite temperature to zero temperature with different cooling protocols. The main quantity of interest is the excess defect density that represents the total defect density minus the equilibrium defect density at varying temperatures. We find that, for three classes of cooling protocols, the time dependence of the excess defect density for various cooling speed shows a dynamic scaling behavior that largely encompasses the Kibble-Zurek mechanism as well as Krapivsky's calculation of the final defect density at zero temperature. We also find distinct features in the behavior of the dynamic scaling when the temperature approaches in a power-law fashion to zero temperature and the excess defect density reaches a peak at finite temperature, where the scaling of the excess defect density at its peak and that at zero temperature exhibits different asymptotic behavior due to different logarithmic corrections. This characteristic behavior can be attributed to the exponential divergence of the relaxation time near zero temperature. While the qualitative theoretical predictions by Krapivsky on the asymptotic exponents are confirmed, the asymptotic predictions for amplitudes overestimate by up to 40% the simulation results.