Magnetic Switching of a Stoner-Wohlfarth Particle Subjected to a Perpendicular Bias Field

Characterized by uniaxial magnetic anisotropy, the Stoner-Wohlfarth particle experiences a change in magnetization leading to a switch in behavior when tuned by an externally applied field, which relates to the perpendicular bias component (hperp) that remains substantially small in comparison with the constant switching field (h0). The dynamics of the magnetic moment that governs the magnetic switching is studied numerically by solving the Landau-Lifshitz-Gilbert (LLG) equation using the Mathematica code without any physical approximations; the results are compared with the switching time obtained from the analytic method that intricately treats the non-trivial bias field as a perturbation. A good agreement regarding the magnetic switching time (ts) between the numerical calculation and the analytic results is found over a wide initial angle range (0.01 < θ0 < 0.3), as h0 and hperp are 1.5 × K and 0.02 × K, where K represents the anisotropy constant. However, the quality of the analytic approximation starts to deteriorate slightly in contrast to the numerical approach when computing ts in terms of the field that satisfies hperp > 0.15 × K and h0 = 1.5 × K. Additionally, existence of a comparably small perpendicular bias field (hperp << h0) causes ts to decrease in a roughly exponential manner when hperp increases.


Introduction
Growing research attention has been paid to controlling the magnetization reversal in ferromagnetic multilayer configurations in the past two decades [1][2][3][4][5][6]. To primitively describe the magnetic dynamics of single-domain ferromagnets [7], it appears necessary to revisit the Stoner-Wohlfarth (SW) model that was originally proposed by E. C. Stoner and E. P. Wohlfarth in 1948 [8]. Based on their early work, many additional extensions [9][10][11] of the simplified SW model have undergone intensive investigation motivated by the current progress in computational technology applied to solid-state physics. The concept of "magnetic switching" plays a key role in the reading and writing of information and monitoring the magnetization in magnetic material. In our study, the time-dependent magnetization vector of a particle satisfying the SW model requirement is examined and associated with the applied torque (due to externally applied magnetic fields) that causes magnetic moment reversal to occur. For the crucial importance of investigating the dynamical behavior of magnetization in a SW particle modulated by a significantly small bias field, one notices that discerning how to efficiently manipulate the magnetic state of a primitive single-domain ferromagnet is key to the development of ferromagnetic spintronics. Specifically, there exists two major aspects that motivate our current study. First, the local modulation of magnetization being done by tuning the magnitude of the bias magnetic field should be completed in a more efficient manner than the traditional method which uses electric fields to store information in the computer memory named "dynamic random-access memory" (DRAM), because information storage in such a conventional method is transient and easily lost when turning off the computer. Therefore, newer way of memorizing information such as "magnetic random-access memory" (MRAM), which uses long-lasting ferromagnetism and ferroelectricity, is strongly desired. Second, increasing the magnetic switching rates is required by the disk speed and recorded density with respect to the physical equipment in computer memory part. Accordingly, ferromagnetic materials have become increasingly essential and demonstrate many promising practical applications including nano-scale memory devices and magnetic field probes [12][13][14].
At present, our work discusses the switching behavior of a Stoner-Wohlfarth particle (sometimes called single-domain magnetic grain) for the instance when the external fields that involve the driving field (h 0 ) along with the small bias field (h perp ) are simultaneously present. Among the results, the time it takes for the magnetization vector to experience an angle change, which is a bit smaller than π/2, and how the bias field strength (or initial angle parameter) effects the switching time are extensively explored in this paper.
There exist many computational techniques that can be used to reveal the magnetic dynamics governed by the Landau-Lifshitz-Gilbert (LLG) equation [15], which was developed after the Landau-Lifshitz (LL) equation [16,17]. In general, many codes utilizing Monte Carlo methods, fast Fourier transformations, etc., have been utilized to obtain a numerical simulation of the magnetization occurring in actual magnetic grains, as well as the time evolution of other physical quantities, namely, the vortex profile, dipole-dipole interaction, etc. [18][19][20]. In addition to these programs, which aim to acquire the time evolution of magnetic quantities without any physical approximations, one learns that an analytic solution within the framework of the LLG equation provides another valid and effective method to interpret the actual magnetic behavior of various systems [21][22][23], including our study, which uses appropriate analytic approximations, such as the perturbation theory. Specifically, for our work, we compare the magnetic switching time gained using two different computational methods (numerical and analytical approaches) and examine the apparent agreement (or slight discrepancy) between these approaches. This paper is structured as follows. We briefly introduce the fundamental model (LLG equation) used to explain the time evolution of the magnetic moment in the context of the Bloch equation along with the LL equation, with the goal of establishing the analytic approach (involving appropriate approximations) that applies to the switching time computation. Then, the interplay of the magnetic driving field and the perpendicular bias field occurring in a Stoner-Wohlfarth particle is reviewed and contrasted with our calculated results, which suggest that magnetic switching can be effectively tuned using a significantly small bias field. In the discussion section, the magnetic switching time in the case of axial symmetry (no bias field), as well as in the case of a small perpendicular bias field, is analytically examined, exhibiting fair agreement with the numerically obtained data.

Magnetic Moment Dynamics
For the purpose of characterizing the magnetic dynamics that govern the time evolution of the magnetization vector M, the Bloch equation [24] was initially proposed via relation dM/dt = γ 0 × M × H, where γ 0 described the gyromagnetic ratio while H represented the external field. Then, Landau and Lifshitz extended the Bloch equation to the parameter M subjected to the effective field (H eff ), which arose from the negative functional magnetization derivative with respect to the magnetic energy (−δE m /δM), as well as the dissipation terms, which were related to the loss and relaxation processes in ferromagnetic materials. Consequently, the derived Landau-Lifshitz (LL) equation is as follows: where Landau and Lifshitz assumed that the dissipation was described by a nonlinear term (second term in Equation (1)) of the form restricted by the same effective field, while introducing the coefficient λ and taking into account the gyromagnetic ratio γ 0 . The first term on the right-hand side (RHS) of Equation (1) describes the dissipationless motion of the moment. Notice that the LL equation conserves the magnetic moment |M|= (M·× M) 1/2 and takes the scalar product of M and all terms, except for the nonlinear term appearing in Equation (1), which leads to M·(dM/dt) = 0. Therefore, one finds that |M| = M s = a constant, where M s represents the magnitude of the saturation magnetization under the circumstance of a negligible second term in the RHS of Equation (1). This can be illustrated with the example of a uniaxial anisotropic magnet placed into an external field (H eff ) directed along the same axis. The left panel of the figure below shows that without dissipation, the moment performs a circular motion along the parallels of a unit sphere, so that |M| = a constant is maintained. In the absence of dissipation, the precession of the moment persists for an infinitely long time and is shown in the left panel of Figure 1. conserves the magnetic moment |M|= (M·× M) 1/2 and takes the scalar product of M and all terms, except for the nonlinear term appearing in Equation (1), which leads to M·(dM/dt) = 0. Therefore, one finds that |M| = Ms = a constant, where Ms represents the magnitude of the saturation magnetization under the circumstance of a negligible second term in the RHS of Equation (1). This can be illustrated with the example of a uniaxial anisotropic magnet placed into an external field (Heff) directed along the same axis. The left panel of the figure below shows that without dissipation, the moment performs a circular motion along the parallels of a unit sphere, so that |M| = a constant is maintained. In the absence of dissipation, the precession of the moment persists for an infinitely long time and is shown in the left panel of Figure 1. In real experiments, dissipation losses play a crucial role in regulating the dynamics of the moment. Even a small variation in the dissipation torque can cause the moment to change significantly. As the right panel of Figure 1 shows, the dissipation torque is the only torque that actually pushes M toward the minimum magnetic energy, resulting in a magnetization trajectory that is on the spherical surface due to its constant module |M|, and eventually precesses to the direction which is almost aligned along the −Heff field direction.
The divergence problem arises when significantly large dissipation losses (λ >> 1) are available in the LL equation. To conquer this trouble, Gilbert modified the dissipation term by taking a damping quantity of the form α × Ms −1 × M × (dM/dt) into account. α acts as the Gilbert damping coefficient. Consequently, the motion of a magnetic moment in the context of a redefined damping factor along with effective field was mathematically given by the so-called Landau-Lifshitz-Gilbert (LLG) equation: where γ is another gyromagnetic ratio, and α turns out to be the Gilbert damping parameter. It appears straightforward to verify that we can retrieve the LL equation format from the LLG counterpart. Taking the cross product of the left-hand side (LHS) of Equation (2) with M, and using M norm conservation (M·× (dM/dt) = 0), one can obtain the following: Notice that Ms represents the magnitude of the moment vector M. Inserting M × (dM/dt), which roots from the LHS of Equation (3) into Equation (2), one eventually obtains the following: Here, one may identify γ0 as γ × (1 + α 2 ) −1 and λγ0/|M| as α × γ × Ms −1 × (1+α 2 ) −1 , then both equations (LL and LLG) become identical. In real experiments, dissipation losses play a crucial role in regulating the dynamics of the moment. Even a small variation in the dissipation torque can cause the moment to change significantly. As the right panel of Figure 1 shows, the dissipation torque is the only torque that actually pushes M toward the minimum magnetic energy, resulting in a magnetization trajectory that is on the spherical surface due to its constant module |M|, and eventually precesses to the direction which is almost aligned along the −H eff field direction.
The divergence problem arises when significantly large dissipation losses (λ >> 1) are available in the LL equation. To conquer this trouble, Gilbert modified the dissipation term by taking a damping quantity of the form α × M s −1 × M × (dM/dt) into account. α acts as the Gilbert damping coefficient. Consequently, the motion of a magnetic moment in the context of a redefined damping factor along with effective field was mathematically given by the so-called Landau-Lifshitz-Gilbert (LLG) equation: where γ is another gyromagnetic ratio, and α turns out to be the Gilbert damping parameter. It appears straightforward to verify that we can retrieve the LL equation format from the LLG counterpart. Taking the cross product of the left-hand side (LHS) of Equation (2) with M, and using M norm conservation (M·× (dM/dt) = 0), one can obtain the following: Notice that M s represents the magnitude of the moment vector M. Inserting M × (dM/dt), which roots from the LHS of Equation (3) into Equation (2), one eventually obtains the following: Here, one may identify γ 0 as γ × (1 + α 2 ) −1 and λγ 0 /|M| as α × γ × M s −1 × (1+α 2 ) −1 , then both equations (LL and LLG) become identical.
The next step was to transform Equation (4) into the dimensionless form. Dividing both sides of Equation (2) by γ × M s gave rise to the below expression: with the following definitions: n = M/M s and h eff = γ × H eff, one can rearrange Equation (5) to obtain the dimensionless equation: dn/dt = −n × h eff + α × n × dn/dt. Then, taking the cross product of n and the expression involving dn/dt, one can finally determine the following:

Stoner-Wohlfarth Model
To understand how the presence of magnetic anisotropy leads to the existence of preferred directions for moment M, one realizes that the origin of anisotropy energy E A stems from two factors: The crystalline anisotropy ultimately caused by the spin-orbit interactions in the material [25] and the shape anisotropy caused by the energy of the fields induced by the magnetic moment [26]. The properties of these two types of energy remain complex in multi-domain samples, while the anisotropy energy E A (n) of a single-domain case has the uniaxial form: where K is the anisotropy constant that stays positive, and z denotes a unit vector along the z-axis.
The minimum of the anisotropy is achieved for n = +z or n = −z in the context of a positive K, while the z-axis here is defined as the easy axis of the magnetic material. Considering a single-domain magnet that is placed into an external magnetic field H, one notices that an additional contribution to the energy that is equivalent to (-M·× H) arises. Hence, the total energy is consequently defined by the following: In terms of spherical angles, θ and ϕ are commonly used to parameterize the unit vector n = (sinθcosϕ, sinθsinϕ, cosθ), and the total energy E A is rewritten based on the re-expressed component (n x , n y , n z ) as the following: Such energy, which depends solely on angle θ, can now be plotted as a series of E(θ) graphs in the following (see Figure 2). The qualitative shape of each diagram appears to be different relying on the value of |H| = H z /cosθ. The next step was to transform Equation (4) into the dimensionless form. Dividing both sides of Equation (2) by γ × Ms gave rise to the below expression: With the following definitions: n = M/Ms and heff = γ × Heff, one can rearrange Equation (5) to obtain the dimensionless equation: dn/dt = −n × heff + α × n × dn/dt. Then, taking the cross product of n and the expression involving dn/dt, one can finally determine the following:

Stoner-Wohlfarth Model
To understand how the presence of magnetic anisotropy leads to the existence of preferred directions for moment M, one realizes that the origin of anisotropy energy EA stems from two factors: The crystalline anisotropy ultimately caused by the spin-orbit interactions in the material [25] and the shape anisotropy caused by the energy of the fields induced by the magnetic moment [26]. The properties of these two types of energy remain complex in multi-domain samples, while the anisotropy energy EA(n) of a single-domain case has the uniaxial form: where K is the anisotropy constant that stays positive, and z denotes a unit vector along the z-axis. The minimum of the anisotropy is achieved for n = +z or n = −z in the context of a positive K, while the z-axis here is defined as the easy axis of the magnetic material. Considering a single-domain magnet that is placed into an external magnetic field H, one notices that an additional contribution to the energy that is equivalent to (-M·× H) arises. Hence, the total energy is consequently defined by the following: In terms of spherical angles, θ and φ are commonly used to parameterize the unit vector n = (sinθcosφ, sinθsinφ, cosθ), and the total energy EA is rewritten based on the re-expressed component (nx, ny, nz) as the following: Such energy, which depends solely on angle θ, can now be plotted as a series of E(θ) graphs in the following (see Figure 2). The qualitative shape of each diagram appears to be different relying on the value of |H| = Hz/cosθ. As shown in the figure above, for −K/|M| < |H| < K/|M|, there are two minima of energy in the angle interval (0, π), while for |H| > K/|M| (or |H| < −K/|M|) there exists only one minimum within the same parameter range. If one initially sets the field to a positive value, for instance, |H| > K/|M|, the available energy minimum can only be acquired at θ = 0, which corresponds to the moment As shown in the figure above, for −K/|M| < |H| < K/|M|, there are two minima of energy in the angle interval (0, π), while for |H| > K/|M| (or |H| < −K/|M|) there exists only one minimum within the same parameter range. If one initially sets the field to a positive value, for instance, |H| > K/|M|, the available energy minimum can only be acquired at θ = 0, which corresponds to the moment pointing in the +z direction. By decreasing the field until it eventually approaches its negative counterpart (from |H| = 1.4 × K/|M| to |H| = −1.4 × K/|M| in Figure 2), one can switch the moment from up (cosθ = 1) to down (cosθ = −1). The switching behavior indeed occurs when the energy minimum in the curves at θ = 0 disappears.
It is apparent to deduce from Equation (9) that this critical phenomenon can be numerically captured at the field value |H| = −K/|M|. Continuing to tune the field toward −1.4 × K/|M|, one finds that the moment abruptly jumps into a state with θ = π, (i.e., being directed along −z). In other words, the moment stays in the downstate situation. If the field is increased up from a negative value (e.g., |H| < −K/|M|), the moment remains in the minimum state (θ = π) until the field reaches the value of |H| = +K/|M|, at which point this minimum configuration disappears and the moment eventually jumps back into the "new" stable minimum at θ = 0. Consequently, Figure 3 shows the magnetic energy E as a function of angle θ for the special cases where |H| = ±K/|M|. The |M|(|H|) dependence has a form shown below ( Figure 4). from up (cosθ = 1) to down (cosθ = −1). The switching behavior indeed occurs when the energy minimum in the curves at θ = 0 disappears. It is apparent to deduce from Equation (9) that this critical phenomenon can be numerically captured at the field value |H| = −K/|M|. Continuing to tune the field toward −1.4 × K/|M|, one finds that the moment abruptly jumps into a state with θ = π, (i.e., being directed along -z). In other words, the moment stays in the downstate situation. If the field is increased up from a negative value (e.g., |H| < −K/|M|), the moment remains in the minimum state (θ = π) until the field reaches the value of |H| = +K/|M|, at which point this minimum configuration disappears and the moment eventually jumps back into the "new" stable minimum at θ = 0. Consequently, Figure 3 shows the magnetic energy E as a function of angle θ for the special cases where |H| = ±K/|M|. The |M|(|H|) dependence has a form shown below ( Figure 4).  It is observed that this magnetic hysteresis has a rectangular hysteresis loop, and its width (or coercive field) is K/|M|. A more complicated situation where H has two nonzero components given by H = (Hx, 0, Hz) is of great significance and is analogous to a Stoner-Wohlfarth particle. It is assumed that the longitudinal component Hz stays extremely larger than the perpendicular component Hx which is free to be tuned and is defined as a bias field for later discussion. Currently, the study of the from up (cosθ = 1) to down (cosθ = −1). The switching behavior indeed occurs when the energy minimum in the curves at θ = 0 disappears. It is apparent to deduce from Equation (9) that this critical phenomenon can be numerically captured at the field value |H| = −K/|M|. Continuing to tune the field toward −1.4 × K/|M|, one finds that the moment abruptly jumps into a state with θ = π, (i.e., being directed along -z). In other words, the moment stays in the downstate situation. If the field is increased up from a negative value (e.g., |H| < −K/|M|), the moment remains in the minimum state (θ = π) until the field reaches the value of |H| = +K/|M|, at which point this minimum configuration disappears and the moment eventually jumps back into the "new" stable minimum at θ = 0. Consequently, Figure 3 shows the magnetic energy E as a function of angle θ for the special cases where |H| = ±K/|M|. The |M|(|H|) dependence has a form shown below ( Figure 4).  It is observed that this magnetic hysteresis has a rectangular hysteresis loop, and its width (or coercive field) is K/|M|. A more complicated situation where H has two nonzero components given by H = (Hx, 0, Hz) is of great significance and is analogous to a Stoner-Wohlfarth particle. It is assumed that the longitudinal component Hz stays extremely larger than the perpendicular component Hx which is free to be tuned and is defined as a bias field for later discussion. Currently, the study of the switching behavior of a uniaxial magnet in the presence of a bias field, in addition to a longitudinal field, has been attracting increasing attention [27,28]. For the purpose of seeking the critical value of It is observed that this magnetic hysteresis has a rectangular hysteresis loop, and its width (or coercive field) is K/|M|. A more complicated situation where H has two nonzero components given by H = (H x , 0, H z ) is of great significance and is analogous to a Stoner-Wohlfarth particle. It is assumed that the longitudinal component H z stays extremely larger than the perpendicular component H x which is free to be tuned and is defined as a bias field for later discussion. Currently, the study of the switching behavior of a uniaxial magnet in the presence of a bias field, in addition to a longitudinal field, has been attracting increasing attention [27,28]. For the purpose of seeking the critical value of H x at which the switching occurs supposing that H z is fixed, the functional dependence H z (H x ) describing interplay of these two critical fields is plotted in the following figure and given mathematically as [11,[29][30][31]: Depicting the above relation between the bias field and the fixed field leads to the resulting diagram, which is referred to as the asteroid curve, in which the dashed line in Figure 5 shows how the total field H changes when H z varies itself at a fixed value of H x . describing interplay of these two critical fields is plotted in the following figure and given mathematically as [11,[29][30][31]: Depicting the above relation between the bias field and the fixed field leads to the resulting diagram, which is referred to as the asteroid curve, in which the dashed line in Figure 5 shows how the total field H changes when Hz varies itself at a fixed value of Hx. It is assumed that the longitudinal field Hz is abruptly tuned from the field point A to the one corresponding to the field at point B, as a function of a constant Hx. This is due to the fact that, tangents to the asteroid give rise to magnetization directions with external energy. For a system with a uniaxial anisotropy, the tangents that are closer to the easy axis can result in stable solutions (minimum energy configuration). Specifically, in the initial state A, the magnetic moment is pointing mostly in the +z direction because of a slight deviation from +z due to the presence of a nonzero, yet very small, Hx (Hx << K/|M|). After an abrupt change of Hz to the value confined at state B, the "up equilibrium" no longer exists, and the moment rotates to the "down equilibrium"; however, the moment is directed almost along −z (again, a small deviation from -z is present due to Hx). An important characteristic of this switching process is numerically manifested by the switching time ts, since such parameter plays a crucial role in technological applications, pertinent to which the single-domain particles are used to store digital information. In other words, a smaller switching time means information is written into memory equipment faster. In this paper, the switching time of a uniaxial Stoner-Wohlfarth particle at a significantly small bias field was investigated by virtue of an analytical method along with a numerical calculation.

Switching Behavior under Axial Symmetry
Before performing practical calculations to determine the magnetic switching time of an SW particle subjected to an additional bias field, it was necessary to convert the derived vector Equation (6) within the small damping approximation (α << 1) to a system of two scalar equations on (θ, φ).  It is assumed that the longitudinal field H z is abruptly tuned from the field point A to the one corresponding to the field at point B, as a function of a constant H x . This is due to the fact that, tangents to the asteroid give rise to magnetization directions with external energy. For a system with a uniaxial anisotropy, the tangents that are closer to the easy axis can result in stable solutions (minimum energy configuration). Specifically, in the initial state A, the magnetic moment is pointing mostly in the +z direction because of a slight deviation from +z due to the presence of a nonzero, yet very small, H x (H x << K/|M|). After an abrupt change of H z to the value confined at state B, the "up equilibrium" no longer exists, and the moment rotates to the "down equilibrium"; however, the moment is directed almost along −z (again, a small deviation from −z is present due to H x ). An important characteristic of this switching process is numerically manifested by the switching time t s , since such parameter plays a crucial role in technological applications, pertinent to which the single-domain particles are used to store digital information. In other words, a smaller switching time means information is written into memory equipment faster. In this paper, the switching time of a uniaxial Stoner-Wohlfarth particle at a significantly small bias field was investigated by virtue of an analytical method along with a numerical calculation.

Switching Behavior under Axial Symmetry
Before performing practical calculations to determine the magnetic switching time of an SW particle subjected to an additional bias field, it was necessary to convert the derived vector Equation (6) within the small damping approximation (α << 1) to a system of two scalar equations on (θ, ϕ). Figure 6 introduces the two vectors e θ and e ϕ , both of which are tangent to the unit sphere. By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθsinφ, −sinθ) and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordinate system, including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the form: where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equation (6) is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent vectors appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with respect to the angles, one finally finds the following: We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in the context where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K × cos 2 θ − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle variables are given: Here, conducting an integration of Equation (14) over θ and t subjected to a provided initial angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching time. The switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment positioned close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) sufficient to induce a switch. A representation of the dependence of θ on t (in the unit of α × K) according to the model described mathematically in Equation (14) is shown in Figure 7. As the switching process progressed, different curves reflecting angle θ (given in the unit of radian) versus time, arose due to the distinct values of θ0. By the end of the switching process, the moment had approached the -z direction, and the angle θ approached π. By differentiating the expressions n = (sinθcosϕ, sinθsinϕ, cosθ), e θ = (cosθcosϕ, cosθsinϕ, −sinθ) and e ϕ = (−sinϕ, cosϕ, 0) with respect to time t while making use of the right-hand coordinate system, including a standard unit basis (e.g., e θ × e ϕ = n), one can transform Equation (6) into the form: By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθsinφ, −sinθ) and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordinate system, including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the form: where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equation (6) is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent vectors appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with respect to the angles, one finally finds the following: dϕ/dt = sin −1 θ × (əESW/əθ) -α × sin −2 θ × (əESW/əφ) We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in the context where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K × cos 2 θ − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle variables are given: Here, conducting an integration of Equation (14) over θ and t subjected to a provided initial angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching time. The switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment positioned close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) sufficient to induce a switch. A representation of the dependence of θ on t (in the unit of α × K) according to the model described mathematically in Equation (14) is shown in Figure 7. As the switching process progressed, different curves reflecting angle θ (given in the unit of radian) versus time, arose due to the distinct values of θ0. By the end of the switching process, the moment had approached the -z direction, and the angle θ approached π. By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθsinφ, −sinθ) and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordinate system, including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the form: where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equation (6) is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent vectors appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with respect to the angles, one finally finds the following: dϕ/dt = sin −1 θ × (əESW/əθ) -α × sin −2 θ × (əESW/əφ) We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in the context where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K × cos 2 θ − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle variables are given: Here, conducting an integration of Equation (14) over θ and t subjected to a provided initial angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching time. The switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment positioned close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) sufficient to induce a switch. A representation of the dependence of θ on t (in the unit of α × K) according to the model described mathematically in Equation (14) is shown in Figure 7. As the switching process progressed, different curves reflecting angle θ (given in the unit of radian) versus time, arose due to the distinct values of θ0. By the end of the switching process, the moment had approached the -z direction, and the angle θ approached π. By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cos and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coord including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS o is approximated as dn/dt. Next, by equating the coefficients pertinent to both ta appearing in the LHS of Equation (11) and rearranging the appropriate time derivative to the angles, one finally finds the following: We now investigated the switching behavior in the case of axial symmetry (Hx = 0) where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1 − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle given: Here, conducting an integration of Equation (14) over θ and t subjected to a pr angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switchi switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic mome close to the +z direction and applies a negative magnetic field -H0 (aligned along −z induce a switch. A representation of the dependence of θ on t (in the unit of α × K) acc model described mathematically in Equation (14) is shown in Figure 7. As the swit progressed, different curves reflecting angle θ (given in the unit of radian) versus time the distinct values of θ0. By the end of the switching process, the moment had appr direction, and the angle θ approached π. By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθco and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) where ESW is identified as the energy pertaining to an SW particle of interest, and the is approximated as dn/dt. Next, by equating the coefficients pertinent to b appearing in the LHS of Equation (11) and rearranging the appropriate time deri to the angles, one finally finds the following: We now investigated the switching behavior in the case of axial symmetry (H where the magnetic energy was assumed to depend solely on the θ parameter (E − |M| × |H| × cosθ). Therefore, two independent equations associated with the given: Here, conducting an integration of Equation (14) over θ and t subjected t angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic s switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic close to the +z direction and applies a negative magnetic field -H0 (aligned alo induce a switch. A representation of the dependence of θ on t (in the unit of α × model described mathematically in Equation (14) is shown in Figure 7. As th progressed, different curves reflecting angle θ (given in the unit of radian) versu the distinct values of θ0. By the end of the switching process, the moment had direction, and the angle θ approached π.
where E SW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equation (6) is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent vectors appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with respect to the angles, one finally finds the following: dθ/dt = −sin −1 θ × ( By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθsinφ, −sinθ) and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordinate system, including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the form: where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equation (6) is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent vectors appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with respect to the angles, one finally finds the following: dϕ/dt = sin −1 θ × (əESW/əθ) -α × sin −2 θ × (əESW/əφ) We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in the context where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K × cos 2 θ − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle variables are given: Here, conducting an integration of Equation (14) over θ and t subjected to a provided initial angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching time. The switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment positioned close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) sufficient to induce a switch. A representation of the dependence of θ on t (in the unit of α × K) according to the model described mathematically in Equation (14) is shown in Figure 7. As the switching process progressed, different curves reflecting angle θ (given in the unit of radian) versus time, arose due to By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθsinφ, −sinθ) and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordinate system, including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the form: where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equation (6) is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent vectors appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with respect to the angles, one finally finds the following: dϕ/dt = sin −1 θ × (əESW/əθ) -α × sin −2 θ × (əESW/əφ) We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in the context where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K × cos 2 θ − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle variables are given: Here, conducting an integration of Equation (14) over θ and t subjected to a provided initial angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching time. The switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment positioned close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) sufficient to induce a switch. A representation of the dependence of θ on t (in the unit of α × K) according to the model described mathematically in Equation (14) is shown in Figure 7. As the switching process progressed, different curves reflecting angle θ (given in the unit of radian) versus time, arose due to ϕ) − α × ( By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθsinφ, − and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordinate sys including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the form: where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equatio is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent ve appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with re to the angles, one finally finds the following: We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in the co where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K × c − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle variable given: Here, conducting an integration of Equation (14) over θ and t subjected to a provided i angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching time. switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment positi close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) sufficie induce a switch. A representation of the dependence of θ on t (in the unit of α × K) according t model described mathematically in Equation (14) is shown in Figure 7. As the switching pr progressed, different curves reflecting angle θ (given in the unit of radian) versus time, arose d By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθs and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordin including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the f (dθ/dt) × eθ where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of E is approximated as dn/dt. Next, by equating the coefficients pertinent to both tang appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives w to the angles, one finally finds the following: We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle va given: Here, conducting an integration of Equation (14) over θ and t subjected to a prov angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) s induce a switch. A representation of the dependence of θ on t (in the unit of α × K) accor model described mathematically in Equation (14) is shown in Figure 7. As the switchi progressed, different curves reflecting angle θ (given in the unit of radian) versus time, a θ) (12) dφ/dt = sin −1 θ × ( By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθsinφ, −sinθ) and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordinate system, including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the form: where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equation (6) is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent vectors appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with respect to the angles, one finally finds the following: We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in the context where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K × cos 2 θ − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle variables are given: Here, conducting an integration of Equation (14) over θ and t subjected to a provided initial angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching time. The switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment positioned close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) sufficient to induce a switch. A representation of the dependence of θ on t (in the unit of α × K) according to the By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθsinφ, −sinθ) and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordinate system, including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the form: where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equation (6) is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent vectors appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with respect to the angles, one finally finds the following: dϕ/dt = sin −1 θ × (əESW/əθ) -α × sin −2 θ × (əESW/əφ) We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in the context where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K × cos 2 θ − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle variables are given: Here, conducting an integration of Equation (14) over θ and t subjected to a provided initial angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching time. The switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment positioned close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) sufficient to induce a switch. A representation of the dependence of θ on t (in the unit of α × K) according to the By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, cosθsin and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coordinate including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into the form where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS of Equ is approximated as dn/dt. Next, by equating the coefficients pertinent to both tangent appearing in the LHS of Equation (11) and rearranging the appropriate time derivatives with to the angles, one finally finds the following: We now investigated the switching behavior in the case of axial symmetry (Hx = 0) in the where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle varia given: Here, conducting an integration of Equation (14) over θ and t subjected to a provide angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching ti switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic moment po close to the +z direction and applies a negative magnetic field -H0 (aligned along −z) suff induce a switch. A representation of the dependence of θ on t (in the unit of α × K) accordin By differentiating the expressions n = (sinθcosφ, sinθsinφ, cosθ), eθ = (cosθcosφ, co and eφ = (−sinφ, cosφ, 0) with respect to time t while making use of the right-hand coor including a standard unit basis (e.g., eθ × eφ = n), one can transform Equation (6) into th where ESW is identified as the energy pertaining to an SW particle of interest, and the LHS is approximated as dn/dt. Next, by equating the coefficients pertinent to both ta appearing in the LHS of Equation (11) and rearranging the appropriate time derivativ to the angles, one finally finds the following: We now investigated the switching behavior in the case of axial symmetry (Hx = 0 where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = − − |M| × |H| × cosθ). Therefore, two independent equations associated with the angl given: Here, conducting an integration of Equation (14) over θ and t subjected to a p angle θ0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switch switching field hz is denoted as Hz × |M|. Suppose one initiates with a magnetic mom close to the +z direction and applies a negative magnetic field -H0 (aligned along − induce a switch. A representation of the dependence of θ on t (in the unit of α × K) ac ϕ) We now investigated the switching behavior in the case of axial symmetry (H x = 0) in the context where the magnetic energy was assumed to depend solely on the θ parameter (E(θ) = −1/2 × K × cos 2 θ − |M| × |H| × cosθ). Therefore, two independent equations associated with the angle variables are given: Here, conducting an integration of Equation (14) over θ and t subjected to a provided initial angle θ 0 ≡ θ(t = 0) was apparently responsible for resolving the magnetic switching time. The switching field h z is denoted as H z × |M|. Suppose one initiates with a magnetic moment positioned close to the +z direction and applies a negative magnetic field −H 0 (aligned along −z) sufficient to induce a switch. A representation of the dependence of θ on t (in the unit of α × K) according to the model described mathematically in Equation (14) is shown in Figure 7. As the switching process progressed, different curves reflecting angle θ (given in the unit of radian) versus time, arose due to the distinct values of θ 0 . By the end of the switching process, the moment had approached the −z direction, and the angle θ approached π. Two features of this graph are important for us to revisit. First, the switching behavior occurred much faster for larger values of θ0. Second, each trajectory path describing the time-dependent angle approached a final value that was sufficiently close to π in an exponential form. That is, the curve formally took a comparably long time to achieve the entire variation process (from θ < 0.5 to θ = π). To efficiently define a finite switching time, it is customary to select a cutoff angle θs and to calculate the corresponding time ts it takes to reach that angle. Figure 7 shows that the switching behavior confined in the angle range of 0 to π exists over a relatively short time interval. Prior to this interval, the angles inside different curves were close to θ0, and beyond this interval, the angles remained close to π. The midpoint for the time of this interval was roughly located at θs = π/2. This explains why it made sense to select π/2 as a cutoff angle. Such a critical angle value has been frequently used in the literature [32,33]. Additionally, it is noted that at the same time (t = ts), the projection moment Mz had reached zero.
Continuing with the exact analytic expression of Equations (12) and (14), we resolved the analytic solution to obtain the switching time of the magnetic particle in the context of axial symmetry by means of integration using separation of variables. By integrating Equation (14) over time parameter t as well as angle θ while imposing the condition θ(t = 0) = θ0, the analytic solution finally becomes the following: One crucial thing about the understanding of the initial angle θ0 requires emphasis. The switching time t (0) s(θ0) ultimately approached zero once θ0 was close to π/2. Additionally, this conclusion was certainly consistent with common sense. A similar expression concerning Equation (15) was derived in another report, in particular for the case where K = 0, and for an arbitrary K [34].

Switching Behavior of an SW particle in a Small Perpendicular Bias Field
In the presence of the switching field H0 pointing in the negative z-axis and the bias field Hperp pointing in the positive x-axis, the energy density of a Stoner-Wohlfarth particle is expressed as the following: Here, the switching field and the bias field were also recognized as hperp = Hperp × |M| and hz = H0 × |M|. If the bias field Hperp was significantly small (hperp = |hperp| << hz = |hz|), one can hope to find an analytic approximation to obtain the switching time by taking such a bias field into account as a Two features of this graph are important for us to revisit. First, the switching behavior occurred much faster for larger values of θ 0 . Second, each trajectory path describing the time-dependent angle approached a final value that was sufficiently close to π in an exponential form. That is, the curve formally took a comparably long time to achieve the entire variation process (from θ < 0.5 to θ = π). To efficiently define a finite switching time, it is customary to select a cutoff angle θ s and to calculate the corresponding time t s it takes to reach that angle. Figure 7 shows that the switching behavior confined in the angle range of 0 to π exists over a relatively short time interval. Prior to this interval, the angles inside different curves were close to θ 0 , and beyond this interval, the angles remained close to π. The midpoint for the time of this interval was roughly located at θ s = π/2. This explains why it made sense to select π/2 as a cutoff angle. Such a critical angle value has been frequently used in the literature [32,33]. Additionally, it is noted that at the same time (t = t s ), the projection moment M z had reached zero.
Continuing with the exact analytic expression of Equations (12) and (14), we resolved the analytic solution to obtain the switching time of the magnetic particle in the context of axial symmetry by means of integration using separation of variables. By integrating Equation (14) over time parameter t as well as angle θ while imposing the condition θ(t = 0) = θ 0 , the analytic solution finally becomes the following: One crucial thing about the understanding of the initial angle θ 0 requires emphasis. The switching time t (0) s (θ 0 ) ultimately approached zero once θ 0 was close to π/2. Additionally, this conclusion was certainly consistent with common sense. A similar expression concerning Equation (15) was derived in another report, in particular for the case where K = 0, and for an arbitrary K [34].

Switching Behavior of an SW particle in a Small Perpendicular Bias Field
In the presence of the switching field H 0 pointing in the negative z-axis and the bias field H perp pointing in the positive x-axis, the energy density of a Stoner-Wohlfarth particle is expressed as the following: Here, the switching field and the bias field were also recognized as h perp = H perp × |M| and h z = H 0 × |M|. If the bias field H perp was significantly small (h perp = |h perp | << h z = |h z |), one can hope to find an analytic approximation to obtain the switching time by taking such a bias field into account as a small perturbation. Similar reports addressing the precessional magnetization switching under a biased perpendicular anisotropy can be found in [35,36]. The result can be formulated as follows. Combining Equation (12) and (15), one notes that the magnetic switching time was identified from the following expression after integration over angle θ and time t, respectively: Treating the −h perp × sinφ term as a perturbation and depending on the previously obtained switching time appearing in Equation (15), the determination of t (0) s (θ 0 ) in the case of a perpendicular bias field was found by substituting an "effective initial angle" θ in for the initial angle θ 0 . The initial unit vector n(t = 0) was designated by the angles (θ 0 , ϕ 0 ). Accordingly, this effective initial angle was defined as the angle between the initial unit vector n(t = 0) and the position of the total energy maximum, which was characterized by the spherical angles θ max = θ * = |H perp | × (|H 0 | − K/|M|) −1 , φ max = π. In a word, the following approximation was effective as a method of computing the magnetic switching time as long as the angle θ * stayed typically small (e.g., θ * < 0.5), while the calculational accuracy turned out to be quadratic in terms in θ * : Furthermore, relying on the primitive LLG model, we numerically solved the magnetic switching time of an SW particle in the case of a small perpendicular bias field via the computational package Mathematica [37] while fulfilling the requirement that θ(t s ) = π/2. Performing this calculational process for different initial directions governed by n(t = 0), we determined under what circumstance the approximate expression (18) may lose its validity. First, we tested the case in which the field H perp was assumed to be constant, and the moment was set up at the initial position stated by (θ 0 , 0), while the switching field satisfied h 0 = h z = 1.5 × K. Using the analytic solution derived from Equation (18) along with the numerically solved LLG equation, we then plotted the t s (θ 0 ) relations in Figure 8, where time parameter was shown in the unit of α × K, and the magnitude of H perp was alternatively given by h perp = 0.02 × K. The initial angle θ 0 is expressed in the unit of radian. small perturbation. Similar reports addressing the precessional magnetization switching under a biased perpendicular anisotropy can be found in [35,36]. The result can be formulated as follows.
Combining Equation (12) and (15), one notes that the magnetic switching time was identified from the following expression after integration over angle θ and time t, respectively: Treating the −hperp × sinϕ term as a perturbation and depending on the previously obtained switching time appearing in Equation (15), the determination of t (0) s(θ0) in the case of a perpendicular bias field was found by substituting an "effective initial angle" θin for the initial angle θ0. The initial unit vector n(t = 0) was designated by the angles (θ0, φ0). Accordingly, this effective initial angle was defined as the angle between the initial unit vector n(t = 0) and the position of the total energy maximum, which was characterized by the spherical angles θmax = θ* = |Hperp| × (|H0| − K/|M|) −1 , ϕmax = π. In a word, the following approximation was effective as a method of computing the magnetic switching time as long as the angle θ* stayed typically small (e.g., θ* < 0.5), while the calculational accuracy turned out to be quadratic in terms in θ*: Furthermore, relying on the primitive LLG model, we numerically solved the magnetic switching time of an SW particle in the case of a small perpendicular bias field via the computational package Mathematica [37] while fulfilling the requirement that θ(ts) = π/2. Performing this calculational process for different initial directions governed by n(t = 0), we determined under what circumstance the approximate expression (18) may lose its validity. First, we tested the case in which the field Hperp was assumed to be constant, and the moment was set up at the initial position stated by (θ0, 0), while the switching field satisfied h0 = hz = 1.5 × K. Using the analytic solution derived from Equation (18) along with the numerically solved LLG equation, we then plotted the ts(θ0) relations in Figure 8, where time parameter was shown in the unit of α × K, and the magnitude of Hperp was alternatively given by hperp = 0.02 × K. The initial angle θ0 is expressed in the unit of radian. Figure 8. The magnetic switching time ts calculated using a numerical approach (points) and an analytical method (solid line) is shown as a function of the initial angle θ0 in the context of the switching field (h0 = 1.5 × K) supplied along the −z axis, in addition to the bias field (hperp = 0.02 × K) given in the (π/2, 0) direction.
In general, the analytical results obtained on the basis of Equation (12), (16) and (18) qualitatively agreed well with the numerical solutions (solid line) considering that the perturbation parameter θ* Figure 8. The magnetic switching time t s calculated using a numerical approach (points) and an analytical method (solid line) is shown as a function of the initial angle θ 0 in the context of the switching field (h 0 = 1.5 × K) supplied along the −z axis, in addition to the bias field (h perp = 0.02 × K) given in the (π/2, 0) direction.
In general, the anasof Equations (12), (16) and (18) qualitatively agreed well with the numerical solutions (solid line) considering that the pserturbation parameter θ * was kept small throughout the entire simulation process. Moreover, we observed from the graph that consistent and accurate correspondence was found among these two different calculational approaches for nearly all of the values of θ 0 . Additionally, the existence of the typically small bias field (h perp = 0.02 × K << h 0 = 1.5 × K) somehow significantly decreased the switching time from approximately 600 × α × K at θ 0 ≈ 0.01 to approximately 320 × α × K at θ 0 = 0.15 (corresponding to approximately 8.6 • ).
Second, we considered a situation where the vector n(t = 0) started its motion from the energy minimum state corresponding to the applied field H 0 , which was aligned along the positive z-axis. This typical starting point (initial unit vector) was parameterized by the angles θ 0 = h perp × (h 0 + K) −1 and φ 0 = 0. Figure 9 illustrates how analytic approximation coincides with the numerical results represented by discrete points. was kept small throughout the entire simulation process. Moreover, we observed from the graph that consistent and accurate correspondence was found among these two different calculational approaches for nearly all of the values of θ0. Additionally, the existence of the typically small bias field (hperp = 0.02 × K << h0 = 1.5 × K) somehow significantly decreased the switching time from approximately 600 × α × K at θ0 ≈ 0.01 to approximately 320 × α × K at θ0 = 0.15 (corresponding to approximately 8.6°). Second, we considered a situation where the vector n(t = 0) started its motion from the energy minimum state corresponding to the applied field H0, which was aligned along the positive z-axis. This typical starting point (initial unit vector) was parameterized by the angles θ0 = hperp × (h0 + K) −1 and ϕ0 = 0. Figure 9 illustrates how analytic approximation coincides with the numerical results represented by discrete points. Figure 9. Magnetic switching time ts calculated using the numerical approach (points) and analytical method (solid line) is shown as a function of the bias field hperp applied along the (π/2, 0) direction in the context of the switching field (h0 = 1.5 × K) supplied along the −z axis in addition to the initial angle determined to be hperp × (h0 + K) −1 .
Compared to Figure 8, the correspondence among the two different calculational approaches became slightly worse as the bias field hperp was increased beyond approximately 0.15 × K, as seen in Figure 9. On the other hand, with the increase of the bias field hperp lying in the range of 0.01 × K and 0.15 × K, the determined switching time arising from the analytic approximation decreased and agreed fairly well with the numerical data points, indicating that the remainder term (O(θ* 2 )) occurring in Equation (18) became non-trivial. Although the quality of the approximation begins to deteriorate slightly when hperp > 0.15 × K, the analytical approach in studying the magnetic switching is still believed to be acceptable because the actual bias field is normally set to be a few percent of the anisotropy field K.

Conclusions
The SW model appears to be quite an effective simulation of magnetic systems and exhibits extremely rich physical properties from, for instance, a dynamic point of view. The combined use of a uniaxial anisotropy and an externally applied field provided a thorough insight into the magnetic switching behavior tuned by a weak bias component hperp (or small initial angle θ0). The validity of the analytic approximation that treats a comparably small hperp (<<h0) as an effective perturbation is confirmed by the good agreement found between the numerical solution to the LLG equation and the results arising from the analytic method. Utilizing these two computational approaches, we consistently found that the calculated magnetic switching time during which the moment vector changed its angle by π/2 − θ0 demonstrated decreasing behavior in a somewhat exponential manner when hperp was increased from 0.01 × K to 0.15 × K (or θ0 was increased from 0.01 to 0.3). Furthermore, Figure 9. Magnetic switching time t s calculated using the numerical approach (points) and analytical method (solid line) is shown as a function of the bias field h perp applied along the (π/2, 0) direction in the context of the switching field (h 0 = 1.5 × K) supplied along the −z axis in addition to the initial angle determined to be h perp × (h 0 + K) −1 .
Compared to Figure 8, the correspondence among the two different calculational approaches became slightly worse as the bias field h perp was increased beyond approximately 0.15 × K, as seen in Figure 9. On the other hand, with the increase of the bias field h perp lying in the range of 0.01 × K and 0.15 × K, the determined switching time arising from the analytic approximation decreased and agreed fairly well with the numerical data points, indicating that the remainder term (O(θ * 2 )) occurring in Equation (18) became non-trivial. Although the quality of the approximation begins to deteriorate slightly when h perp > 0.15 × K, the analytical approach in studying the magnetic switching is still believed to be acceptable because the actual bias field is normally set to be a few percent of the anisotropy field K.

Conclusions
The SW model appears to be quite an effective simulation of magnetic systems and exhibits extremely rich physical properties from, for instance, a dynamic point of view. The combined use of a uniaxial anisotropy and an externally applied field provided a thorough insight into the magnetic switching behavior tuned by a weak bias component h perp (or small initial angle θ 0 ). The validity of the analytic approximation that treats a comparably small h perp (<<h 0 ) as an effective perturbation is confirmed by the good agreement found between the numerical solution to the LLG equation and the results arising from the analytic method. Utilizing these two computational approaches, we consistently found that the calculated magnetic switching time during which the moment vector changed its angle by π/2 − θ 0 demonstrated decreasing behavior in a somewhat exponential manner when h perp was increased from 0.01 × K to 0.15 × K (or θ 0 was increased from 0.01 to 0.3). Furthermore, a slight discrepancy existed between the numerical data points and the analytic results, given that h perp was tuned beyond 0.15 × K when t s was shown as a function of the perpendicular bias field. It is expected that the derived analytic expression, which treats h perp as small perturbation, can be extended to simplify the determinations for magnetic switching in granular magnetic media; however, the bias field can originate from the induced field of other grains or from the misalignment of the individual grain's axes with the applied field.
Author Contributions: D.X., the first author, undertook all of the (sometimes computationally intense) analytic calculations that are reported in the paper. He also wrote the majority of the first draft and participated in the editing of that draft which produced the final, submitted version. W.M., the second author, is a visiting Professor of Physics at Umeå University with beneficial suggestions and comments.