Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical review. B, Condensed matter

Response of Ising systems to oscillating and pulsed fields: Hysteresis, ac, and pulse susceptibility

Muktish Acharyya, Bikas K. Chakrabarti

Published Sep 1, 1995• 232 citations
Publisher

Abstract

We have studied, using Monte Carlo (MC) simulation for ferromagnetic Ising systems in one to four dimensions and solving numerically the mean-field (MF) equation of motion, the nature of the response magnetization m(t) of an Ising system in the presence of a periodically varying external field [h(t)=${\mathit{h}}_{0}$cos(\ensuremath{\omega}t)]. From these studies, we determine the m-h loop or hysteresis loop area A(=\ensuremath{\oint}mdh) and the dynamic order parameter Q(=\ensuremath{\oint}mdt) and investigate their variations with the frequency (\ensuremath{\omega}) and amplitude (${\mathit{h}}_{0}$) of the applied external magnetic field and the temperature (T) of the system. The variations in A are fitted to a scaling form, assumed to be valid over a wide range of parameter (\ensuremath{\omega},${\mathit{h}}_{0}$,T) values, and the best-fit exponents are obtained in all three dimensions (D=2,3,4). The scaling function is Lorentzian in the MF case and is exponentially decaying, with an initial power law, for the MC cases. The dynamic phase boundary (in the ${\mathit{h}}_{0}$-T plane) is found to be frequency dependent and the transition (from Q\ensuremath{\ne}0 for low T and ${\mathit{h}}_{0}$ to Q=0 for high T and ${\mathit{h}}_{0}$) across the boundary crosses over from a discontinuous to a continuous one at a tricritical point. These boundaries are determined in various cases.We find that the response can be generally expressed as m(t)=P(\ensuremath{\omega}(t-${\mathrm{\ensuremath{\tau}}}_{\mathrm{eff}}$)) where P denotes a periodic function with the same frequency \ensuremath{\omega} of the perturbing field and ${\mathrm{\ensuremath{\tau}}}_{\mathit{e}\mathit{f}\mathit{f}}$(${\mathit{h}}_{0}$,\ensuremath{\omega},T) denotes the effective delay. We established that this effective delay ${\mathrm{\ensuremath{\tau}}}_{\mathit{e}\mathit{f}\mathit{f}}$ of the response is the crucial term and it practically determines all the above observations for A, Q, etc. Investigating the nature of the in-phase (\ensuremath{\chi}\ensuremath{'}) and the out-of-phase (\ensuremath{\chi}\ensuremath{'}\ensuremath{'}) susceptibility, defined as \ensuremath{\chi}\ensuremath{'}=(${\mathit{m}}_{0}$/${\mathit{h}}_{0}$)cos(\ensuremath{\varphi}) and \ensuremath{\chi}\ensuremath{'}\ensuremath{'}=(${\mathit{m}}_{0}$/${\mathit{h}}_{0}$)sin(\ensuremath{\varphi}); \ensuremath{\varphi}=\ensuremath{\omega}${\mathrm{\ensuremath{\tau}}}_{\mathit{e}\mathit{f}\mathit{f}}$ [and ${\mathit{m}}_{0}$ is the amplitude of m(t)], we find that the loop area A is directly given by \ensuremath{\chi}\ensuremath{'}\ensuremath{'} and also the temperature variation of \ensuremath{\chi}\ensuremath{'}\ensuremath{'}, at fixed \ensuremath{\omega} and ${\mathit{h}}_{0}$, gives a prominent peak at the dynamic transition point. We have also studied the behavior of the response magnetization by the application of a short-duration (compared with the relaxation time) pulsed magnetic field.Here, we observed (both in the MC and MF cases) that the width ratio (of the half-width and the width of the response magnetization and of the pulsed field, respectively) and the susceptibility (the ratio of excess magnetization over its equilibrium value and the height of the pulsed field) both show sharp peaks at the order-disorder (ferromagnetic-paramagnetic) transition point. We have also studied similar response behavior of Ising systems in the presence of time-varying longitudinal and transverse fields, solving numerically the mean-field equation of motion. We have again studied here the nature of the dynamic phase transition and the behavior of the ac susceptibility (both longitudinal and transverse) across the dynamic phase boundary. For a short-duration pulse of the transverse field, the width ratio and the pulse susceptibility are again seen to diverge at the order-disorder transition point.

DOI
10.1103/physrevb.52.6550

Related papers

Dynamic response of an Ising system to a pulsed field

Muktish Acharyya, Jayanta K. Bhattacharjee, Bikas K. Chakrabarti

The dynamical response to a pulsed magnetic field has been studied here both using Monte Carlo simulation and by solving numerically the mean-field dynamical equation of motion for the Ising model. The ratio ${\mathrm{R}}_{\mathrm{p}}$ of the response magnetization half-width to the width of the external field pulse has been observed to diverge and pulse susceptibility ${\mathrm{\ensuremath{\chi}}}_{\mathrm{p}}$ (ratio of the response magnetization peak height and the pulse height) gives a peak near the order-disorder transition temperature ${\mathrm{T}}_{\mathrm{c}}$ (for the unperturbed system). The Monte Carlo results for the Ising system on a square lattice show that ${\mathrm{R}}_{\mathrm{p}}$ diverges at ${\mathrm{T}}_{\mathrm{c}}$, with the exponent \ensuremath{\nu}z\ensuremath{\cong}2.0, while ${\mathrm{\ensuremath{\chi}}}_{\mathrm{p}}$ shows a peak at ${\mathrm{T}}_{\mathrm{c}}^{\mathrm{e}}$, which is a function of the field pulse width \ensuremath{\delta}t. A finite-size (in time) scaling analysis shows that ${\mathrm{T}}_{\mathrm{c}}^{\mathrm{e}}$=${\mathrm{T}}_{\mathrm{c}}$+C(\ensuremath{\delta}t${)}^{\mathrm{\ensuremath{-}}1\mathrm{/}\mathrm{x}}$, with x=\ensuremath{\nu}z\ensuremath{\cong}2.0. The mean-field results show that both the divergence of R and the peak in ${\mathrm{\ensuremath{\chi}}}_{\mathrm{p}}$ occur at the mean-field transition temperature, while the peak height in ${\mathrm{\ensuremath{\chi}}}_{\mathrm{p}}$\ensuremath{\sim}(\ensuremath{\delta}t${)}^{\mathrm{y}}$, y\ensuremath{\cong}1 for small values of \ensuremath{\delta}t. These results also compare well with an approximate analytical solution of the mean-field equation of motion.

arXiv:cond-mat/9611182•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Mar 1, 1997•18 citations

Dynamic magnetization-reversal transition in the Ising model

Arkajyoti Misra, Bikas K. Chakrabarti

We report the results of mean field and the Monte Carlo study of the dynamic magnetization-reversal transition in the Ising model, brought about by the application of an external field pulse applied in opposition to the existing order before the application of the pulse. The parameters controlling the transition are the strength ${h}_{p}$ and the duration $\ensuremath{\Delta}t$ of the pulse. In the mean field case, an approximate analytical expression is obtained for the phase boundary which agrees well with that obtained numerically in the small $\ensuremath{\Delta}t$ and large $T$ limit. The order parameter of the transition has been identified, and is observed to vary continuously near the transition. The order parameter exponent $\ensuremath{\beta}$ was estimated both for the mean field $(\ensuremath{\beta}=1)$ and the Monte Carlo $(\ensuremath{\beta}=0.90\ifmmode\pm\else\textpm\fi{}0.02$ in two dimensions) cases. The transition shows a ``critical slowing down'' type behavior near the phase boundary with diverging relaxation time. The divergence was found to be logarithmic in the mean field case and exponential in the Monte Carlo case. The finite size scaling technique was employed to estimate the correlation length exponent $\ensuremath{\nu}$ $(=1.5\ifmmode\pm\else\textpm\fi{}0.3$ in two dimensions) in the Monte Carlo case.

arXiv:cond-mat/9805289•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Oct 1, 1998•18 citations

High-frequency magnetic field on crystal field diluted S = 1 Ising system: magnetic relaxation near continuous phase transition points

Gül Gülpınar, Rıza Erdem

Magnetization relaxation and the steady state response of the S = 1 Ising model with random crystal field to a time varying magnetic field with a frequency ω is modelled and studied here by a method that combines the statistical equilibrium theory with the theory of irreversible thermodynamics. The method offers information on the relaxation time (τ) of the system as well as the temperature (θ) and ω dependencies of the complex (AC or dynamical) susceptibility (i.e., χ(ω) = χ′(ω) − iχ″(ω)). The so-called low- and high-frequency regions are separated by τ because τ−1 → 0 as θ approaches the critical temperatures (θc). One can choose to keep the frequency ω fixed and observe the low-frequency behaviors followed by the high-frequency behaviors when θ → θc. It is shown that χ(ω) exhibits different behaviors in low- and high-frequency regimes that are separated by the quantity ωτ: χ′(ω) converges to static susceptibility and χ″(ω) → 0 for ωτ ≪ 1. However, in the high-frequency region where ωτ ≫ 1, χ′(ω) vanishes and χ″(ω) displays a peak at the critical temperature (θc). Besides the above, the logarithm of the susceptibility components versus log(ω) is also plotted. From these plots, one plateau (a step-like) region and a shifted peak with rising temperature is observed for the real and imaginary parts, respectively.

Canadian Journal of Physics•Apr 11, 2018•4 citations

Competing field pulse induced dynamic transition in Ising Models

Arnab Chatterjee, Bikas K. Chakrabarti

We discuss dynamic magnetization-reversal phenomena in the Ising model under a finite-duration external magnetic field, which competes with the existing order for . The nature of the phase boundary is estimated from the mean-field (MF) equation of motion. The susceptibility and relaxation time diverge at the MF phase boundary. A Monte Carlo study also shows divergence of the relaxation time and order-parameter fluctuations at the phase boundary. The fourth-order cumulant shows two distinct behaviors. For low temperatures and pulse durations, the value of the cumulant at the crossing point for different system sizes is far less than that for the static transition in the same dimension. This suggests a new universality class for the dynamic transition. For higher temperatures and pulse durations, the transition falls in a MF-like weak-singularity universality class.

Phase Transitions•May 1, 2004•14 citations

Hysteresis and self-organized criticality in the O(N) model in the limit N to infinity

Deepak Dhar, P B Thomas

The authors consider the response of the ferromagnetic N-vector model to a sinusoidally varying external magnetic field in the large-N limit. In all dimensions d>2, they show that at low frequencies omega , and small amplitudes H 0 of the field, the area of the hysteresis loop scales as (H 0 omega ) 1/2 with logarithmic corrections. At very high frequencies, the area varies as H 0 2 / omega . They find that for any H 0 there is a dynamical phase transition separating these two frequency regimes. They determine numerically the critical frequency as a function of the field strength. In the high-frequency phase the magnetization is predominantly transverse to the external magnetic field.

Journal of Physics A Mathematical and General•Oct 7, 1992•91 citations

Ising model in an oscillating magnetic field: Mean-field theory

Michael F. Zimmer

I consider the dynamics of a soft-spin Ising magnet, subject to the time-varying external magnetic field h(t)=h cos\ensuremath{\Omega}t. The system is modeled with a time-dependent Ginzburg-Landau equation and is studied at a mean-field level. The time-averaged magnetization (scrM) acts as the order parameter and divides the temperature-h plane into two phases. In contrast to a previous result that predicted a dynamical tricritical point separating a line of continuous and discontinuous transitions in scrM, I find that the transition is always continuous. The previous work utilized Glauber dynamics which, it is here argued, becomes inapplicable near the phase boundary line due to a critical slowing down.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jun 1, 1993•62 citations
We have studied, using Monte Carlo (MC) simulation for ferromagnetic Ising systems in one to four dimensions and solving numerically the mean-field (MF) equation of motion, the nature of the response magnetization m(t) of an Ising system in the presence of a periodically varying external field [h(t)=h0{\mathit{h}}_{0}h0​cos(\ensuremath{\omega}t)]. From these studies, we determine the m-h loop or hysteresis loop area A(=\ensuremath{\oint}mdh) and the dynamic order parameter Q(=\ensuremath{\oint}mdt) and investigate their variations with the frequency (\ensuremath{\omega}) and amplitude (h0{\mathit{h}}_{0}h0​) of the applied external magnetic field and the temperature (T) of the system. The variations in A are fitted to a scaling form, assumed to be valid over a wide range of parameter (\ensuremath{\omega},h0{\mathit{h}}_{0}h0​,T) values, and the best-fit exponents are obtained in all three dimensions (D=2,3,4). The scaling function is Lorentzian in the MF case and is exponentially decaying, with an initial power law, for the MC cases. The dynamic phase boundary (in the h0{\mathit{h}}_{0}h0​-T plane) is found to be frequency dependent and the transition (from Q\ensuremath{\ne}0 for low T and h0{\mathit{h}}_{0}h0​ to Q=0 for high T and h0{\mathit{h}}_{0}h0​) across the boundary crosses over from a discontinuous to a continuous one at a tricritical point. These boundaries are determined in various cases.We find that the response can be generally expressed as m(t)=P(\ensuremath{\omega}(t-\ensuremathτeff{\mathrm{\ensuremath{\tau}}}_{\mathrm{eff}}\ensuremathτeff​)) where P denotes a periodic function with the same frequency \ensuremath{\omega} of the perturbing field and \ensuremathτeff{\mathrm{\ensuremath{\tau}}}_{\mathit{e}\mathit{f}\mathit{f}}\ensuremathτeff​(h0{\mathit{h}}_{0}h0​,\ensuremath{\omega},T) denotes the effective delay. We established that this effective delay \ensuremathτeff{\mathrm{\ensuremath{\tau}}}_{\mathit{e}\mathit{f}\mathit{f}}\ensuremathτeff​ of the response is the crucial term and it practically determines all the above observations for A, Q, etc. Investigating the nature of the in-phase (\ensuremath{\chi}\ensuremath{'}) and the out-of-phase (\ensuremath{\chi}\ensuremath{'}\ensuremath{'}) susceptibility, defined as \ensuremath{\chi}\ensuremath{'}=(m0{\mathit{m}}_{0}m0​/h0{\mathit{h}}_{0}h0​)cos(\ensuremath{\varphi}) and \ensuremath{\chi}\ensuremath{'}\ensuremath{'}=(m0{\mathit{m}}_{0}m0​/h0{\mathit{h}}_{0}h0​)sin(\ensuremath{\varphi}); \ensuremath{\varphi}=\ensuremath{\omega}\ensuremathτeff{\mathrm{\ensuremath{\tau}}}_{\mathit{e}\mathit{f}\mathit{f}}\ensuremathτeff​ [and m0{\mathit{m}}_{0}m0​ is the amplitude of m(t)], we find that the loop area A is directly given by \ensuremath{\chi}\ensuremath{'}\ensuremath{'} and also the temperature variation of \ensuremath{\chi}\ensuremath{'}\ensuremath{'}, at fixed \ensuremath{\omega} and h0{\mathit{h}}_{0}h0​, gives a prominent peak at the dynamic transition point. We have also studied the behavior of the response magnetization by the application of a short-duration (compared with the relaxation time) pulsed magnetic field.Here, we observed (both in the MC and MF cases) that the width ratio (of the half-width and the width of the response magnetization and of the pulsed field, respectively) and the susceptibility (the ratio of excess magnetization over its equilibrium value and the height of the pulsed field) both show sharp peaks at the order-disorder (ferromagnetic-paramagnetic) transition point. We have also studied similar response behavior of Ising systems in the presence of time-varying longitudinal and transverse fields, solving numerically the mean-field equation of motion. We have again studied here the nature of the dynamic phase transition and the behavior of the ac susceptibility (both longitudinal and transverse) across the dynamic phase boundary. For a short-duration pulse of the transverse field, the width ratio and the pulse susceptibility are again seen to diverge at the order-disorder transition point.
The dynamical response to a pulsed magnetic field has been studied here both using Monte Carlo simulation and by solving numerically the mean-field dynamical equation of motion for the Ising model. The ratio Rp{\mathrm{R}}_{\mathrm{p}}Rp​ of the response magnetization half-width to the width of the external field pulse has been observed to diverge and pulse susceptibility ${\mathrm{\ensuremath{\chi}}}_{\mathrm{p}}$ (ratio of the response magnetization peak height and the pulse height) gives a peak near the order-disorder transition temperature Tc{\mathrm{T}}_{\mathrm{c}}Tc​ (for the unperturbed system). The Monte Carlo results for the Ising system on a square lattice show that Rp{\mathrm{R}}_{\mathrm{p}}Rp​ diverges at Tc{\mathrm{T}}_{\mathrm{c}}Tc​, with the exponent \ensuremath{\nu}z\ensuremath{\cong}2.0, while ${\mathrm{\ensuremath{\chi}}}_{\mathrm{p}}$ shows a peak at Tce{\mathrm{T}}_{\mathrm{c}}^{\mathrm{e}}Tce​, which is a function of the field pulse width \ensuremath{\delta}t. A finite-size (in time) scaling analysis shows that Tce{\mathrm{T}}_{\mathrm{c}}^{\mathrm{e}}Tce​=Tc{\mathrm{T}}_{\mathrm{c}}Tc​+C(\ensuremath{\delta}t${)}^{\mathrm{\ensuremath{-}}1\mathrm{/}\mathrm{x}}$, with x=\ensuremath{\nu}z\ensuremath{\cong}2.0. The mean-field results show that both the divergence of R and the peak in ${\mathrm{\ensuremath{\chi}}}_{\mathrm{p}}$ occur at the mean-field transition temperature, while the peak height in ${\mathrm{\ensuremath{\chi}}}_{\mathrm{p}}$\ensuremath{\sim}(\ensuremath{\delta}t)y{)}^{\mathrm{y}})y, y\ensuremath{\cong}1 for small values of \ensuremath{\delta}t. These results also compare well with an approximate analytical solution of the mean-field equation of motion.
We report the results of mean field and the Monte Carlo study of the dynamic magnetization-reversal transition in the Ising model, brought about by the application of an external field pulse applied in opposition to the existing order before the application of the pulse. The parameters controlling the transition are the strength hp{h}_{p}hp​ and the duration $\ensuremath{\Delta}t$ of the pulse. In the mean field case, an approximate analytical expression is obtained for the phase boundary which agrees well with that obtained numerically in the small $\ensuremath{\Delta}t$ and large TTT limit. The order parameter of the transition has been identified, and is observed to vary continuously near the transition. The order parameter exponent $\ensuremath{\beta}$ was estimated both for the mean field $(\ensuremath{\beta}=1)$ and the Monte Carlo $(\ensuremath{\beta}=0.90\ifmmode\pm\else\textpm\fi{}0.02$ in two dimensions) cases. The transition shows a ``critical slowing down'' type behavior near the phase boundary with diverging relaxation time. The divergence was found to be logarithmic in the mean field case and exponential in the Monte Carlo case. The finite size scaling technique was employed to estimate the correlation length exponent $\ensuremath{\nu}$ $(=1.5\ifmmode\pm\else\textpm\fi{}0.3$ in two dimensions) in the Monte Carlo case.
Magnetization relaxation and the steady state response of the S = 1 Ising model with random crystal field to a time varying magnetic field with a frequency ω is modelled and studied here by a method that combines the statistical equilibrium theory with the theory of irreversible thermodynamics. The method offers information on the relaxation time (τ) of the system as well as the temperature (θ) and ω dependencies of the complex (AC or dynamical) susceptibility (i.e., χ(ω) = χ′(ω) − iχ″(ω)). The so-called low- and high-frequency regions are separated by τ because τ−1 → 0 as θ approaches the critical temperatures (θc). One can choose to keep the frequency ω fixed and observe the low-frequency behaviors followed by the high-frequency behaviors when θ → θc. It is shown that χ(ω) exhibits different behaviors in low- and high-frequency regimes that are separated by the quantity ωτ: χ′(ω) converges to static susceptibility and χ″(ω) → 0 for ωτ ≪ 1. However, in the high-frequency region where ωτ ≫ 1, χ′(ω) vanishes and χ″(ω) displays a peak at the critical temperature (θc). Besides the above, the logarithm of the susceptibility components versus log(ω) is also plotted. From these plots, one plateau (a step-like) region and a shifted peak with rising temperature is observed for the real and imaginary parts, respectively.
We discuss dynamic magnetization-reversal phenomena in the Ising model under a finite-duration external magnetic field, which competes with the existing order for . The nature of the phase boundary is estimated from the mean-field (MF) equation of motion. The susceptibility and relaxation time diverge at the MF phase boundary. A Monte Carlo study also shows divergence of the relaxation time and order-parameter fluctuations at the phase boundary. The fourth-order cumulant shows two distinct behaviors. For low temperatures and pulse durations, the value of the cumulant at the crossing point for different system sizes is far less than that for the static transition in the same dimension. This suggests a new universality class for the dynamic transition. For higher temperatures and pulse durations, the transition falls in a MF-like weak-singularity universality class.
The authors consider the response of the ferromagnetic N-vector model to a sinusoidally varying external magnetic field in the large-N limit. In all dimensions d>2, they show that at low frequencies omega , and small amplitudes H 0 of the field, the area of the hysteresis loop scales as (H 0 omega ) 1/2 with logarithmic corrections. At very high frequencies, the area varies as H 0 2 / omega . They find that for any H 0 there is a dynamical phase transition separating these two frequency regimes. They determine numerically the critical frequency as a function of the field strength. In the high-frequency phase the magnetization is predominantly transverse to the external magnetic field.
I consider the dynamics of a soft-spin Ising magnet, subject to the time-varying external magnetic field h(t)=h cos\ensuremath{\Omega}t. The system is modeled with a time-dependent Ginzburg-Landau equation and is studied at a mean-field level. The time-averaged magnetization (scrM) acts as the order parameter and divides the temperature-h plane into two phases. In contrast to a previous result that predicted a dynamical tricritical point separating a line of continuous and discontinuous transitions in scrM, I find that the transition is always continuous. The previous work utilized Glauber dynamics which, it is here argued, becomes inapplicable near the phase boundary line due to a critical slowing down.