Search

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

Nonequilibrium phase transition properties of a mixed Ising ferrimagnetic model consisting of spin-1/2 and spin-3/2 on a square lattice under the existence of a time dependent oscillating magnetic field have been investigated by making use of Monte Carlo simulations with single-spin flip Metropolis algorithm. A complete picture of dynamic phase boundary and magnetization profiles have been illustrated and the conditions of a dynamic compensation behavior have been discussed in detail. According to our simulation results, the considered system does not point out a dynamic compensation behavior, when it only includes the nearest-neighbor interaction, single-ion anisotropy and an oscillating magnetic field source. As the next-nearest-neighbor interaction between the spins-1/2 takes into account and exceeds a characteristic value which sensitively depends upon values of single-ion anisotropy and only of amplitude of external magnetic field, a dynamic compensation behavior occurs in the system. Finally, it is reported that it has not been found any evidence of dynamically first-order phase transition between dynamically ordered and disordered phases, which conflicts with the recently published molecular field investigation, for a wide range of selected system parameters., Comment: 10 pages, 7 figures

Nonequilibrium phase transition properties of the $\pm J$ Ising model under a time dependent oscillating perturbation are investigated within the framework of effective field theory for a two-dimensional square lattice. After a detailed analysis, it is found that the studied system exhibits unusual and interesting behaviors such as reentrant phenomena, and the competition between ferromagnetic and antiferromagnetic exchange interactions gives rise to destruct the dynamic first order phase transitions as well as dynamic tricritical points. Furthermore, according to N\'{e}el nomenclature, the magnetization profiles have been found to obey Q-type, L-type and P-type classification schemes under certain conditions. Finally, it is observed that the treatment of critical percolation with applied field amplitude strongly depends upon the frequency of time varying external field., Comment: 8 pages, 6 figures

Monte Carlo simulation based on Metropolis algorithm has been used with a great success to analyze the dynamic phase transition properties of a single spherical core-shell nanoparticle system with a spin-3/2 core surrounded by a spin-1 shell layer with antiferromagnetic interface coupling under the influence of a time dependent oscillating magnetic filed. It has been found that the dynamic phase boundaries strongly depend on the Hamiltonian parameters such as for the high amplitude and period values of the external field, the phase transition temperature sharply changes whereas it tends to slowly alter as the reduced magnitude of interlayer parameter $\mathrm{J_{int}/J_{sh}}$ increases. Moreover, it is observed that the magnetization curves of the particle have been found to obey P-type, N-type and Q-type classification schemes under certain conditions. Much effort has also been paid to the influence of the particle size on the thermal and magnetic properties of the particle. Finally, a comparison of our observations with those of recently published study including dynamic treatments of a nanocubic core-shell system is represented and the findings indicate that there exists a qualitatively good agreement with some relatively distinct differences., Comment: This paper has been withdrawn by the author because it has been extended and updated

The relaxation and complex magnetic susceptibility treatments of a spin-3/2 Blume-Capel model with quenched random crystal field on a two dimensional square lattice are investigated by a method combining the statistical equilibrium theory and the thermodynamics of linear irreversible processes. Generalized force and flux are defined in irreversible thermodynamics limit. The kinetic equation for the magnetization is obtained by using linear response theory. Temperature and also crystal field dependencies of the relaxation time are obtained in the vicinity of phase transition points. We found that the relaxation time exhibits divergent treatment near the order-disorder phase transition point as well as near the isolated critical point whereas it displays cusp behavior near the first order phase transition point. In addition, much effort has been devoted to investigation of complex magnetic susceptibility response of the system to changing applied field frequencies and it is observed that the considered disordered magnetic system exhibits unusual and interesting behaviors. Furthermore, dynamical mean field critical exponents for the relaxation time and complex magnetic susceptibility are calculated in order to formulate the critical behavior of the system. Finally, a comparison of our observations with those of recently published studies is represented and it is shown that there exists a qualitatively good agreement., Comment: 13 pages, 8 figures

Frequency evolutions of hysteresis loop area and hysteresis tools such as remanence and coercivity of a kinetic Ising model in the presence of quenched bond dilution are investigated in detail. The kinetic equation describing the time dependence of the magnetization is derived by means of effective-field theory with single-site correlations. It is found that the frequency dispersions of hysteresis loop area, remanence and coercivity strongly depend on the quenched bond randomness, as well as applied field amplitude and oscillation frequency. In addition, the shape of the hysteresis curves for a wide variety of Hamiltonian parameters is studied and some interesting behaviors are found. Finally, a comparison of our observations with those of recently published studies is represented and it is shown that there exists a qualitatively good agreement., Comment: 11 pages, 8 figures

We have presented dynamic phase transition features and stationary-state behavior of a ferrimagnetic small nanoparticle system with a core-shell structure. By means of detailed Monte Carlo simulations, a complete picture of the phase diagrams and magnetization profiles have been presented and the conditions for the occurrence of a compensation point $T_{comp}$ in the system have been investigated. According to N\'{e}el nomenclature, the magnetization curves of the particle have been found to obey P-type, N-type and Q-type classification schemes under certain conditions. Much effort has been devoted to investigation of hysteretic response of the particle and we observed the existence of triple hysteresis loop behavior which originates from the existence of a weak ferromagnetic core coupling $J_{c}/J_{sh}$, as well as a strong antiferromagnetic interface exchange interaction $J_{int}/J_{sh}$. Most of the calculations have been performed for a particle in the presence of oscillating fields of very high frequencies and high amplitudes in comparison with exchange interactions which resembles a magnetic system under the influence of ultrafast switching fields. Particular attention has also been paid on the influence of the particle size on the thermal and magnetic properties, as well as magnetic features such as coercivity, remanence and compensation temperature of the particle. We have found that in the presence of ultrafast switching fields, the particle may exhibit a dynamic phase transition from paramagnetic to a dynamically ordered phase with increasing ferromagnetic shell thickness., Comment: 12 pages, 12 figures

Nonequilibrium behavior and dynamic phase transition properties of a kinetic Ising model under the influence of periodically oscillating random-fields have been analyzed within the framework of effective field theory (EFT) based on a decoupling approximation (DA). Dynamic equation of motion has been solved for a simple cubic lattice ($q=6$) by utilizing a Glauber type stochastic process. Amplitude of the sinusoidally oscillating magnetic field is randomly distributed on the lattice sites according to bimodal and trimodal distribution functions. For a bimodal type of amplitude distribution, it is found in the high frequency regime that the dynamic phase diagrams of the system in temperature versus field amplitude plane resemble the corresponding phase diagrams of pure kinetic Ising model. Our numerical results indicate that for a bimodal distribution, both in the low and high frequency regimes, the dynamic phase diagrams always exhibit a coexistence region in which the stationary state (ferro or para) of the system is completely dependent on the initial conditions whereas for a trimodal distribution, coexistence region disappears depending on the values of system parameters., Comment: 11 pages, 11 figures

In this work, the relaxation process of the spin-3/2 Blume-Capel model with quenched random crystal field on a two dimensional square lattice has been investigated by a method which combines the statistical equilibrium theory and the thermodynamic of linear irreversible processes. Generalized force and flux have been defined in irreversible thermodynamics limit. The kinetic equation for the magnetization has been obtained by using linear response theory. In order to understand effect of the first-order phase transition on the relaxation process, low temperature regions have been analyzed. And also, dependence on the temperature and random crystal field of the relaxation process has been investigated near the second-order and isolated critical phase transition points. Moreover, evidence of the reentrant behavior has been observed on the relaxation process. Finally, critical exponents of the relaxation time and magnetization have been found by benefiting from log-log plot., Comment: This paper has been withdrawn by the author because it has been extended and updated

This investigation employs a square lattice and Glauber dynamics methodology to probe the effects of diluteness in the crystal-field interaction in a Blume-Capel Ising system under an oscillating magnetic-field. Fourteen different phase diagrams have been observed in temperature-magnetic-field space as the concentration of the crystal-field interactions is varied. Besides, a comparison is given with the results of the pure spin-1 Ising systems., Comment: 25 pages, 24 figures