Search

We present a comprehensive numerical study of dynamic phase transitions in the two-dimensional kinetic Ising model under a non-antisymmetric time-dependent magnetic field including a sinusoidal term and a second harmonic component. We demonstrate that the expected antisymmetric property and the scaling behavior of the order parameter are maintained using the recently proposed generalized conjugate field approach. Via a detailed finite-size scaling analysis we compute, for zero-bias field, the set of critical exponents suggesting that the Ising universality class is conserved, even in the absence of half-wave antisymmetry in the time-dependent magnetic field. Our results verify up-to-date experimental observations and provide a deeper understanding of non-equilibrium phase transitions, establishing a broader framework for exploring symmetry-breaking phenomena in driven magnetic systems., Comment: 9 pages, 7 figures, submitted for publication

We investigate the dynamical critical behavior of the two-dimensional three-state Potts model with single spin-flip dynamics in equilibrium. We focus on the mean-squared deviation of the magnetization $M$ (MSD$_{M}$) as a function of time, as well as on the autocorrelation function of $M$. Our simulations reveal the existence of two crossover behaviors at times $\tau_1 \sim L^{z_1}$ and $\tau_2 \sim L^{z_2}$, separating three dynamical regimes. MSD$_{M}$ appears to shift from ordinary diffusion in the first regime, to anomalous diffusion in the second, and finally to be constant in the third regime. The magnetization autocorrelation function on the other hand is found to fluctuate between exponential decay, stretched-exponential decay, and then again exponential decay along these three regimes. This behavior is in agreement with the one reported recently for the two-dimensional Ising ferromagnet [Phys. Rev. E {\bf 108}, 034118 (2023)], indicating that the existence of two dynamic critical exponents is not a peculiarity of the Ising model itself. A comparison of both MSD$_{M}$ and the magnetization's autocorrelation function suggests that within our numerical accuracy the exponents $z_1$ and $z_2$ are shared between the Ising and three-state Potts models at least for the particular case of single spin-flip dynamics studied here, even though their equilibrium universality classes are clearly distinct. Continuity of MSD$_{M}$ requires that $\alpha (z_2 - z_1) = \gamma/\nu - z_1$, in which $\alpha$ is the anomalous exponent in the intermediate regime. Since the ratio $\gamma/\nu$ is not shared between the two models, it follows that $\alpha$ is not shared either, an aspect well verified in our simulations. Finally, we also discuss the relevance of our main findings using another useful observable, namely the line magnetization $M_{l}$., Comment: 7 pages, 5 figures, to be published in Phys. Rev. E

In this paper, we study the thermodynamic properties of spin-$1/2$ antiferromagnetic Heisenberg ladders by means of the stochastic series expansion quantum Monte Carlo technique. This includes the thermal properties of the specific heat, uniform and staggered susceptibilities, spin gap, and structure factor. Our numerical simulations are probed over a large ensemble of random realizations in a wide range of disorder strengths $r$, from the clean ($r=0$) case up to the diluted ($r \rightarrow 1$) limit, and for selected choices of number of legs $L_y$ per site. Our results show some interesting phenomena, like the presence of crossing points in the temperature plane for both the specific heat and uniform susceptibility curves which appear to be universal in $r$, as well as a variable dependence of the spin gap in the amount of disorder upon increasing $L_y$., Comment: 9 pages, 6 figures, final version to be published in Eur. Phys. J. B

We investigate the dynamical critical behavior of the two- and three-dimensional Ising model with Glauber dynamics in equilibrium. In contrast to the usual standing, we focus on the mean-squared deviation of the magnetization $M$, MSD$_M$, as a function of time, as well as on the autocorrelation function of $M$. These two functions are distinct but closely related. We find that MSD$_M$ features a first crossover at time $\tau_1 \sim L^{z_{1}}$, from ordinary diffusion with MSD$_M$ $\sim t$, to anomalous diffusion with MSD$_M$ $\sim t^\alpha$. Purely on numerical grounds, we obtain the values $z_1=0.45(5)$ and $\alpha=0.752(5)$ for the two-dimensional Ising ferromagnet. Related to this, the magnetization autocorrelation function crosses over from an exponential decay to a stretched-exponential decay. At later times, we find a second crossover at time $\tau_2 \sim L^{z_{2}}$. Here, MSD$_M$ saturates to its late-time value $\sim L^{2+\gamma/\nu}$, while the autocorrelation function crosses over from stretched-exponential decay to simple exponential one. We also confirm numerically the value $z_{2}=2.1665(12)$, earlier reported as the single dynamic exponent. Continuity of MSD$_M$ requires that $\alpha(z_{2}-z_{1})=\gamma/\nu-z_1$. We speculate that $z_{1} = 1/2$ and $\alpha = 3/4$, values that indeed lead to the expected $z_{2} = 13/6$ result. A complementary analysis for the three-dimensional Ising model provides the estimates $z_{1} = 1.35(2)$, $\alpha=0.90(2)$, and $z_{2} = 2.032(3)$. While $z_{2}$ has attracted significant attention in the literature, we argue that for all practical purposes $z_{1}$ is more important, as it determines the number of statistically independent measurements during a long simulation., Comment: 6 pages, 6 figures, version published in Physical Review E

In the present work, we investigate the effects of disorder on the reversal time ($\tau$) of classical anisotropic Heisenberg ferromagnets in three dimensions by means of Monte Carlo simulations. Starting from the pure system, our analysis suggests that $\tau$ increases with increasing anisotropy strength. On the other hand, for the case of randomly distributed anisotropy, generated from various statistical distributions, a set of results is obtained: (i) For both bimodal and uniform distributions the variation of $\tau$ with the strength of anisotropy strongly depends on temperature. (ii) At lower temperatures, the decrement in $\tau$ with increasing width of the distribution is more prominent. (iii) For the case of normally distributed anisotropy, the variation of $\tau$ with the width of the distribution is non-monotonic, featuring a minimum value that decays exponentially with the temperature. Finally, we elaborate on the joint effect of longitudinal ($h_z$) and transverse ($h_x$) fields on $\tau$, which appear to obey a scaling behavior of the form $\tau h_z^{n} \sim f(h_x)$., Comment: 18 pages preprint style, 10 figures, including Appendix, version to be published in Phys Rev E

We study the question of universality in the two-dimensional spin-$1$ Baxter-Wu model in the presence of a crystal field $\Delta$. We employ extensive numerical simulations of two types, providing us with complementary results: Wang-Landau sampling at fixed values of $\Delta$ 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 $(\Delta, T)$ phase diagram indicates that the transition belongs to the universality class of the $4$-state Potts model. Previous controversies with respect to the nature of the transition are discussed and possibly attributed to the presence of strong finite-size effects, especially as one approaches the pentacritical point of the model., Comment: 10 pages, 12 figures, 1 table, version to be published in Phys. Rev. E

We study the transient behavior of damage propagation in the two-dimensional spin-$1$ Blume-Capel model using Monte Carlo simulations with Metropolis dynamics. We find that, for a particular region in the second-order transition regime of the crystal field--temperature phase diagram of the model, the average Hamming distance decreases exponentially with time in the weakly damaged system. Additionally, its rate of decay appears to depend linearly on a number of Hamiltonian parameters, namely the crystal field, temperature, applied magnetic field, but also on the amount of damage. Finally, a comparative study using Metropolis and Glauber dynamics indicates a slower decay rate of the average Hamming distance for the Glauber protocol., Comment: 17 pages Latex including 8 captioned pdf figures, J. Stat. Phys. (in press) 2022

We investigate aspects of universality in the two-dimensional (2D) spin-$1$ Baxter-Wu model in a crystal field $\Delta$ using a parallel version of the multicanonical algorithm employed at constant temperature $T$. A detailed finite-size scaling analysis in the continuous regime of the $\Delta-T$ phase diagram of the model indicates that the transition belongs to the universality class of the $4$-state Potts model. The presence of first-order-like finite-size effects that become more pronounced as one approaches the pentacritical point of the model is highlighted and discussed., Comment: 6 pages, 6 figures, XXXII IUPAP Conference on Computational Physics

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., Comment: 25 pages (APS preprint style), 13 figures, 1 table

The three-dimensional XY model is investigated in the presence of a uniform magnetic field applied in the $X$-direction. The nearest neighbour intraplanar interaction is considered ferromagnetic, and the interplanar nearest neighbour interaction is chosen to be antiferromagnetic. Starting from a high-temperature initial random spin configuration, the equilibrium phase of the system at any finite temperature was achieved by cooling the system using the Monte Carlo single spin-flip Metropolis algorithm with a random updating rule. The components of total magnetisation and the sublattice magnetisations were calculated. The variance of the antiferromagnetic order parameter and the susceptibility have been calculated. In a specific range of relative strengths of interactions (antiferromagnetic/ferromagnetic) and the applied magnetic field, the system shows the equilibrium phase transitions at different temperatures. The phase diagrams (in the field-temperature plane) were obtained for different values of the relative interaction strengths. The ordered region bounded by the phase boundary was found to increase as the ratio of relative interaction strength increased. Furthermore, the maximum value of the susceptibility ($\chi^m_{ay}$) was found to increase with the system size ($L$). For $\chi^m_{ay} \sim L^{{{\gamma} \over {\nu}}}$, the exponent $\gamma/\nu$ has been estimated to be 2.10$\pm0.11$ ., Comment: 10 pages Latex and 8 captioned eps figures

We present an extensive Monte Carlo investigation of the metastable lifetime through the reversal of the magnetization of spin-$s$ Ising and Blume-Capel models, where $s=\{1/2, 1, 3/2, 2, 5/2, 3, 7/2\}$. The mean metastable lifetime (or reversal time) is studied as a function of the applied magnetic field and for both models is found to obey the Becker-Doring theory, as was initially developed for the case of $s=1/2$ Ising ferromagnet within the classical nucleation theory. Moreover, the decay of metastable volume fraction nicely follows the Avrami's law for all values of $s$ and for both models considered., Comment: 21 preprint-style pages, 7 figures, 3 tables; version to be published in Phys. Rev. E

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 the 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 the 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 smoothening of the transition to second-order with the presence of strong scaling corrections., Comment: 11 pages, 13 figures, minor correction to references appearing in Fig. 1, to be published in Phys. Rev. E

A two-leg quenched random bond disordered antiferromagnetic spin$-1/2$ Heisenberg ladder system is investigated by means of stochastic series expansion (SSE) quantum Monte Carlo (QMC) method. Thermal properties of the uniform and staggered susceptibilities, the structure factor, the specific heat and the spin gap are calculated over a large number of random realizations in a wide range of disorder strength. According to our QMC simulation results, the considered system has a special temperature point at which the specific heat take the same value regardless of the strength of the disorder. Moreover, the uniform susceptibility is shown to display the same character except for a small difference in the location of the special point. Finally, the spin gap values are found to decrease with increasing disorder parameter and the smallest gap value found in this study is well above the weak coupling limit of the clean case., Comment: 15 pages, 9 figures

The nonequilibrium behaviour of a core-shell nanoparticle has been studied by Monte- Carlo simulation. The core consists of Ising spins of $\sigma=1/2$ and the shell contains Ising spins of $S=1$. The interactions within the core and in the shell are considered ferromagnetic but the interfacial interaction between core and shell is antiferromagnetic. The nanoparticle system is kept in open boundary conditions and is driven by randomly varying (in time but uniform over the space) magnetic field. Depending on the width of the randomly varying field and the temperature of the system, the core, shell and total magnetization varies in such a manner that the time averages vanish for higher magnitude of the width of random field, exhibiting a dynamical symmetry breaking transitions. The susceptibilities get peaked at two different temperatures indicating nonequilibrium multiple transitions. The phase boundaries of the nonequilibrium multiple transitions are drawn in the plane formed by the axes of temperature and the width of the randomly varying field. Furthermore, the effects of the core and shell thicknesses on the multiple transitions have been discussed., Comment: 14 pages Latex including 8 captioned figures

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., Comment: 28 pages, 12 figures, version to be published in Physical Review E (minor title correction)

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 non-equilibrium phase transition belongs to the universality class of the equilibrium Ising model and allow us to construct a dynamic phase diagram, in analogy to 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., Comment: 27 pages, 1 table, 15 figures, minor corrections (updated figures 3 and 4)

Magnetic response of the spin-$1/2$ cylindrical nanowire to the propagating magnetic field wave has been investigated by means of Monte Carlo simulation method based on Metropolis algorithm. The obtained microscopic spin configurations suggest that the studied system exhibits two types of dynamical phases depending on the considered values of system parameters: Coherent propagation of spin bands and spin-frozen or pinned phases, as in the case of the conventional bulk systems under the influence of a propagating magnetic field. By benefiting from the temperature dependencies of variances of dynamic order parameter, internal energy and the derivative of dynamic order parameter of the system, dynamic phase diagrams are also obtained in related planes for varying values of the wavelength of the propagating magnetic field. Our simulation results demonstrate that as the strength of the field amplitude is increased, the phase transition points tend to shift to the relatively lower temperature regions. Moreover, it has been observed that dynamic phase boundary line shrinks inward when the value of wavelength of the external field decreases., Comment: 7 pages, 4 figures

We have elucidated the dynamic phase transition features and finite-size scaling analysis of the triangular lattice system under the presence of a square-wave magnetic field. It has been found that as the value of half-period of the external field reaches its critical value, whose location is estimated by means of Binder cumulant, the system presents a dynamic phase transition between dynamically ordered and disordered phases. Moreover, at the dynamic phase transition point, finite-size scaling of the Monte Carlo results for the dynamic order parameter and susceptibility give the critical exponents $\beta/\nu=0.143\pm0.004$ and $\gamma/\nu=1.766\pm0.036$, respectively. The obtained critical exponents show that present magnetic system belongs to same universality class with the two-dimensional equilibrium Ising model., Comment: 5 pages, 7 figures

We elucidate the finite temperature entanglement properties of $N=9$ qubits Heisenberg $XX$ and $XXZ$ models under the presence of a polarized magnetic field in $xz$ plane by means of concurrence concept. We perform a systematic analysis for a wide range of the system parameters. Our results suggest that the global phase regions which separate the entangled and non-entangled regions sensitively depend upon the spin-spin interaction term of the $z-$ component of two neighboring spins $J_{z}/J_{x}$, temperature as well as polarized magnetic field components. Thereby, we think that polarized magnetic field can be used a control parameter to determine the amount of thermal entanglement between pair of qubits for different temperatures and spin-spin interaction terms. Moreover, it has been found that the nearest-neighbor pair of qubits does not point out a re-entrant type entanglement character when one only deals with the nearest-neighbor pair of qubits. However, as one considers next-nearest neighbor pair of qubits, it is possible to see the evidences of re-entrant type entanglement behaviors., Comment: 6 pages, 4 figures

Dynamic phase transition properties of ferromagnetic thin film system under the influence both bias and time dependent magnetic fields have been elucidated by means of kinetic Monte Carlo simulation with local spin update Metropolis algorithm. The obtained results after a detailed analysis suggest that bias field is the conjugate field to dynamic order parameter, and it also appears to define a phase line between two antiparallel dynamic ordered states depending on the considered system parameters. Moreover, the data presented in this study well qualitatively reproduce the recently published experimental findings where time dependent magnetic behavior of a uniaxial cobalt films is studied in the neighborhood of dynamic phase transition point., Comment: 15 pages, 5 Figures

By means of Monte Carlo simulation with local spin update Metropolis algorithm, we have elucidated non-equilibrium phase transition properties and stationary-state treatment of a disordered binary ferromagnetic alloy of the type $A_{p}B_{1-p}$ on a square lattice. After a detailed analysis, we have found that the system shows many interesting and unusual thermal and magnetic behaviors, for instance, the locations of dynamic phase transition points change significantly depending upon amplitude and period of the external magnetic field as well as upon the active concentration of $A-$ type components. Much effort has also been dedicated to clarify the hysteresis tools, such as coercivity, dynamic loop area as well as dynamic correlations between time dependent magnetizations and external time dependent applied field as a functions of period and amplitude of field as well as active concentration of of $A-$ type components, and outstanding physical findings have been reported in order to better understand the dynamic process underlying present system., Comment: 11 pages, 8 figures. Accepted in Physics Letters A

By means of detailed Monte Carlo (MC) simulations, we have presented dynamic phase transition (DPT) properties of ferromagnetic thin-films. Thermal variations of surface, bulk and total dynamical order parameters (DOP) for a film and total order parameter for the films with different thicknesses have been examined. Opposite regimes of the critical value of reduced exchange interaction (surface to bulk ratio) $R_{c}$ at which the critical temperature becomes independent of film thickness $L$ has been also taken into consideration. The average magnetizations of each layer is reversed in these regimes. Based on the results, we have confirmed that the system represents a crossover behavior in between ordinary to extraordinary transition in the presence of surface exchange enhancement., Comment: 6 pages, 5 figures. arXiv admin note: text overlap with arXiv:1306.3562 by other authors