Classical Heisenberg Spin Research Articles

Exact Numerical Solution of the Fully Connected Classical and Quantum Heisenberg Spin Glass.

We present the mean field solution of the quantum and classical Heisenberg spin glasses, using the combination of a high precision numerical solution of the Parisi full replica symmetry breaking equations and a continuous time quantum MonteCarlo algorithm. We characterize the spin glass order and its low-energy excitations down to zero temperature. The Heisenberg spin glass has a rougher energy landscape than its Ising analog, and exhibits a very slow temperature evolution of its dynamical properties. We extend our analysis to the doped, metallic Heisenberg spin glass, which displays unexpectedly slow spin dynamics, reflecting the proximity to the melting quantum critical point and its associated Sachdev-Ye-Kitaev Planckian dynamics.

Nematic ordering in the Heisenberg spin-glass system in three dimensions.

Nematic ordering, where the spins globally align along a spontaneously chosen axis irrespective of direction, occurs in spin-glass systems of classical Heisenberg spins in d=3. In this system where the nearest-neighbor interactions are quenched randomly ferromagnetic or antiferromagnetic, instead of the locally randomly ordered spin-glass phase, the system orders globally as a nematic phase. This nematic ordering of the Heisenberg spin-glass system is dramatically different from the spin-glass ordering of the Ising spin-glass system. The system is solved exactly on a hierarchical lattice and, equivalently, Migdal-Kadanoff approximately on a cubic lattice. The global phase diagram is calculated, exhibiting this nematic phase, and ferromagnetic, antiferromagnetic, disordered phases. The nematic phase of the classical Heisenberg spin-glass system is also found in other dimensions d>2: We calculate nematic transition temperatures in 24 different dimensions in 2<d≤4.

Exact classical spin dynamics of high spin nanoscale molecular magnetic clusters

A common feature of many organic-based molecular magnets is that inter-molecular magnetic interactions are extremely weak compared to intra-molecular ones. As a result a bulk sample can be described in terms of independent individual molecular magnets made up of a cluster of a small finite number of coupled magnetic spins. Because of the large value of individual spins of several molecular magnetic compounds it turns out that such systems can be described to very high accuracy, at sufficiently high temperatures, by a classical Heisenberg spin model. Only for sufficiently low temperatures need one incorporate the quantum character of the spins. In this work we focus on the exact classical spin dynamics of small nanoscale molecular magnetic clusters of classical Heisenberg spins. The present study permits us to obtain analytical solutions for the time-dependent spin-spin autocorrelation function for certain clusters of classical Heisenberg spins coupled through exchange interaction. This study underscores the potential utility of a classical spin treatment of nanoscale magnetic clusters and provides results useful to interpret data obtained from experiments.

Layer effects on the magnetic textures in magnets with local inversion asymmetry

Magnets with broken local inversion symmetries are interesting candidates for chiral magnetic textures such as skyrmions and spin spirals. The property of these magnets is that each subsequent layer can possess a different Dzyaloshniskii-Moriya interaction (DMI) originating from the local inversion symmetry breaking. Given that new candidates of such systems are emerging, with the Van der Waals crystals and magnetic multilayer systems, it is interesting to investigate how the chiral magnetic textures depend on the number of layers and the coupling between them. In this article, we model the magnetic layers with a classical Heisenberg spin model, where the sign of the DMI alternates for each consecutive layer. We use Monte Carlo simulations to examine chiral magnetic textures and show that the pitch of magnetic spirals is influenced by the interlayer coupling and the number of layers. We observe even-odd effects in the number of layers, where we observe a suppression of the spin spirals for even layer numbers. We give an explanation for our findings by proposing a net DMI in systems with strongly coupled layers. Our results can be used to determine the DMI in systems with a known number of layers, and for new technologies based on the tunability of the spiral wave length.

Robust nature of the chiral spin helix in CrNb3S6 nanostructures studied by off-axis electron holography

Magnetic soliton crystals with layered structures that host periodic chiral helimagnetic ordering are promising candidates for spintronic nanodevices. Among them, helimagnetic $\mathrm{Cr}{\mathrm{Nb}}_{3}{\mathrm{S}}_{6}$ is unique owing to its crystallographic chirality and monoaxial Dzyaloshinskii-Moriya interaction. It is crucial to explore its magnetic configurations and properties with respect to the temperature and thickness, especially in reduced dimensions. Here, the chiral helimagnetic ground state in $\mathrm{Cr}{\mathrm{Nb}}_{3}{\mathrm{S}}_{6}$ nanostructures is investigated using off-axis electron holography in the transmission electron microscope. The period of the helical state is found to be independent of both temperature and specimen thickness, while the temperature dependence of the saturation magnetization is shown to follow a classical Heisenberg spin model. Monte Carlo simulations based on a discrete classical Heisenberg model reproduce the experimental observations closely, confirming the applicability of a three-dimensional Heisenberg model even in a confined specimen geometry.

Grain boundary anisotropy on nano-polycrystalline magnetic thin films

Grain boundaries in polycrystalline thin films with crystallite sizes at nanoscale presents regions characterized by a high degree of local structural disorder. As a consequence, great values of the associated local anisotropies are expected. On this regard, a systematic investigation of the effect of the grain boundary anisotropy on the magnetic properties in such type of nanostructured systems is addressed. For developing this work, a standard Monte Carlo simulation in the framework of classical Heisenberg spins was carried out, with a Hamiltonian involving exchange couplings, dipolar interactions, Zeeman interaction and contributions of cubic magneto-crystalline anisotropy. A quantification of local structural disorder was considered. Results revealed that i) by keeping the same number of grains, different organizations give rise to different spontaneous magnetizations, ii) the critical exponent of the magnetization differs of pure models, which is attributed to the complexity of the lattice and consistent with a distribution of critical temperatures, iii) Boundary anisotropy varies with temperature and its strength are determinant factors for blocking temperatures, and iv) Boundary anisotropy inside in the hysteretic properties where coercive field variations are observed.

Invariant Measures for Integrable Spin Chains and an Integrable Discrete Nonlinear Schrödinger Equation

We consider discrete analogues of two well-known open problems regarding invariant measures for dispersive PDE, namely, the invariance of the Gibbs measure for the continuum (classical) Heisenberg model and the invariance of white noise under focusing cubic nonlinear Schrödinger equation. These continuum models are completely integrable and connected by the Hasimoto transform; correspondingly, we focus our attention on discretizations that are also completely integrable and also connected by a discrete Hasimoto transform. We consider these models on the infinite lattice ${\mathbb{Z}}$. Concretely, for a completely integrable variant of the classical Heisenberg spin chain model (introduced independently by Haldane, Ishimori, and Sklyanin) we prove the existence and uniqueness of solutions for initial data following a Gibbs law (which we show is unique) and show that the Gibbs measure is preserved under these dynamics. In the setting of the focusing Ablowitz--Ladik system, we prove invariance of a measure that we will show is the appropriate discrete analogue of white noise. We also include a thorough discussion of the Poisson geometry associated to the discrete Hasimoto transform introduced by Ishimori that connects the two models studied in this article.

Energy and spin diffusion in the one-dimensional classical Heisenberg spin chain at finite and infinite temperatures.

The energy and spin diffusion behaviors in the one-dimensional classical Heisenberg spin chain have been systematically investigated using the equilibrium diffusion method. The spatiotemporal autocorrelation functions for energy and spin are calculated at finite and infinite temperatures. As conserved quantities, the spreading of excess energy and spin can be used to determine their actual diffusion behaviors. At low temperatures, the energy diffusion shows almost ballistic behavior, and spin shows superdiffusion behavior for finite chain size. For energy diffusion, normal diffusion behavior can be obtained when the temperature is higher than 0.75. For spin diffusion, normal diffusion behavior is observed at infinite temperature.

Rotating magnetic field driven antiferromagnetic domain wall motion: Role of Dzyaloshinskii-Moriya interaction

In this work, we study the rotating magnetic field driven domain wall (DW) motion in antiferromagnetic nanowires, using the micromagnetic simulations of the classical Heisenberg spin model. We show that in low frequency region, the rotating field alone could efficiently drive the DW motion even in the absence of Dzyaloshinskii-Moriya interaction (DMI). In this case, the DW rotates synchronously with the magnetic field, and a stable precession torque is available and drives the DW motion with a steady velocity. In large frequency region, the DW only oscillates around its equilibrium position and cannot propagate. The dependences of the velocity and critical frequency differentiating the two motion modes on several parameters are investigated in details, and the direction of the DW motion can be controlled by modulating the initial phase of the field. Interestingly, a unidirectional DW motion is predicted attributing to the bulk DMI, and the nonzero velocity for high frequency is well explained. Thus, this work does provide useful information for future antiferromagnetic spintronics applications.

Relaxation to equilibrium in models of classical spins with long-range interactions

For a model long-range interacting system of classical Heisenberg spins, we study how fluctuations, such as those arising from having a finite system size or through interaction with the environment, affect the dynamical process of relaxation to Boltzmann–Gibbs equilibrium. Under deterministic spin precessional dynamics, we unveil the full range of quasistationary behavior observed during relaxation to equilibrium, whereby the system is trapped in nonequilibrium states for times that diverge with the system size. The corresponding stochastic dynamics, modeling interaction with the environment and constructed in the spirit of the stochastic Landau–Lifshitz–Gilbert equation, however shows a fast relaxation to equilibrium on a size-independent timescale and no signature of quasistationarity, provided the noise is strong enough. Similar fast relaxation is also seen in Glauber Monte Carlo dynamics of the model, thus establishing the ubiquity of what has been reported earlier in particle dynamics (hence distinct from the spin dynamics considered here) of long-range interacting systems, that quasistationarity observed in deterministic dynamics is washed away by fluctuations induced through contact with the environment.

Analytical and numerical investigations of noncollinear magnetic ordering in the frustrated delafossite CuCrO2

The magnetic propagation vector in delafossite CuCrO2 with classical Heisenberg spins is calculated analytically as a function of exchange interactions up to fourth-nearest neighbors. Exchange inte ...

Secular dynamics of long-range interacting particles on a sphere in the axisymmetric limit.

We investigate the secular dynamics of long-range interacting particles moving on a sphere, in the limit of an axisymmetric mean-field potential. We show that this system can be described by the general kinetic equation, the inhomogeneous Balescu-Lenard equation. We use this approach to compute long-term diffusion coefficients, that are compared with direct simulations. Finally, we show how the scaling of the system's relaxation rate with the number of particles fundamentally depends on the underlying frequency profile. This clarifies why systems with a monotonic profile undergo a kinetic blocking and cannot relax as a whole under 1/N resonant effects. Because of its general form, this framework can describe the dynamics of globally coupled classical Heisenberg spins, long-range couplings in liquid crystals, or the orbital inclination evolution of stars in nearly Keplerian systems.

Microwave fields driven domain wall motions in antiferromagnetic nanowires

In this work, we study the microwave field driven domain wall (DW) motion in an antiferromagnetic nanowire, using the numerical calculations based on a classical Heisenberg spin model with the biaxial magnetic anisotropy. We show that a proper combination of a static magnetic field plus an oscillating field perpendicular to the nanowire axis is sufficient to drive the DW propagation along the nanowire. More importantly, the drift velocity at the resonance frequency is comparable to that induced by temperature gradients, suggesting that microwave field can be a very promising tool to control DW motions in antiferromagnetic nanostructures. The dependences of resonance frequency and drift velocity on the static and oscillating fields, the axial anisotropy, and the damping constant are discussed in details. Furthermore, the optimal orientations of the field are also numerically determined and explained. This work provides useful information for the spin dynamics in antiferromagnetic nanostructures for spintronics applications.

Parity-time symmetry breaking in spin chains

We investigate nonequilibrium phase transitions in classical Heisenberg spin chains associated with spontaneous breaking of parity-time ($\mathcal{PT}$) symmetry of the system under the action of Slonczewski spin-transfer torque (STT) modeled by an applied \textit{imaginary} magnetic field. We reveal the STT-driven $\mathcal{PT}$ symmetry-breaking phase transition between the regimes of precessional and exponentially damped spin dynamics and show that its several properties can be derived from the distribution of zeros of the system's partition function, the approach first introduced by Yang and Lee for studying equilibrium phase transitions in Ising spin chains. The physical interpretation of imaginary magnetic field as describing the action of non-conservative forces opens the possibility of direct observations of Lee-Yang zeros in nonequilibrium physical systems.

Excitations of breathers and rogue wave in the Heisenberg spin chain

We study the excitations of breathers and rogue wave in a classical Heisenberg spin chain with twist interaction, which is governed by a fourth-order integrable nonlinear Schrödinger equation. The dynamics of these waves have been extracted from an exact solution. In particular, the corresponding existence conditions based on the parameters of perturbation wave number K, magnon number N, background wave vector ks and amplitude c are presented explicitly. Furthermore, the characteristics of magnetic moment distribution corresponding to these nonlinear waves are also investigated in detail. Finally, we discussed the state transition of three types nonlinear localized waves under the different excitation conditions.

Thermal properties of spin-S Kitaev-Heisenberg model on a honeycomb lattice

Temperature (T) dependence of heat capacity C(T) in the S=1/2 Kitaev honeycomb model shows a double-peak structure resulting from fractionalization of spins into two kinds of Majorana fermions. Recently it has been discussed that the double-peak structure in C(T) is also observed in magnetic ordered phases of the S=1/2 Kitaev-Heisenberg (KH) model on a honeycomb lattice when the system is located in the vicinity of the Kitaev's spin liquid phase. In addition to the S=1/2 spin case, similar double-peak structure has been confirmed in the KH honeycomb model for classical Heisenberg spins, where spin S is regarded as S→∞. We investigate spin-S dependence of C(T) for the KH honeycomb models by using thermal pure quantum state. We also perform classical Monte Carlo calculations to obtain C(T) for the classical KH model. From obtained results, we find that the origin of the high-temperature peak is different between the quantum spin case with small Ss and the classical Heisenberg spin case. Furthermore, the high-temperature peak in the quantum spin case, which is one of the clues for fractionalization of spins, disappears for S>1.

Monte Carlo simulation of dynamic phase transitions and frequency dispersions of hysteresis curves in core/shell ferrimagnetic cubic nanoparticle

By means of Monte Carlo simulation method with Metropolis algorithm, we elucidate the thermal and magnetic phase transition behaviors of a ferrimagnetic core/shell nanocubic system driven by a time dependent magnetic field. The particle core is composed of ferromagnetic spins, and it is surrounded by an antiferromagnetic shell. At the interface of the core/shell particle, we use antiferromagnetic spin–spin coupling. We simulate the nanoparticle using classical Heisenberg spins. After a detailed analysis, our Monte Carlo simulation results suggest that present system exhibits unusual and interesting magnetic behaviors. For example, at the relatively lower temperature regions, an increment in the amplitude of the external field destroys the antiferromagnetism in the shell part of the nanoparticle, leading to a ground state with ferromagnetic character. Moreover, particular attention has been dedicated to the hysteresis behaviors of the system. For the first time, we show that frequency dispersions can be categorized into three groups for a fixed temperature for finite core/shell systems, as in the case of the conventional bulk systems under the influence of an oscillating magnetic field.

Computation of disordered system from the first principles of classical mechanics and ℕℙ hard problem

We study the classical 1D Heisenberg spin glasses in the framework of nearest-neighboring model. Based on the Hamilton equations we obtained the system of recurrence equations which allows to perform node-by-node calculations of a spin-chain. It is shown that calculations from the first principles of classical mechanics lead to ℕℙ hard problem, that however in the limit of the statistical equilibrium can be calculated by ℙ algorithm. For the partition function of the ensemble a new representation is offered in the form of one-dimensional integral of spin-chains’ energy distribution.

Fractionalized Z_{2} Classical Heisenberg Spin Liquids.

Quantum spin systems are by now known to exhibit a large number of different classes of spin liquid phases. By contrast, for classical Heisenberg models, only one kind of fractionalized spin liquid phase, the so-called Coulomb or U(1) spin liquid, has until recently been identified: This exhibits algebraic spin correlations and impurity moments, "orphan spins," whose size is a fraction of that of the underlying microscopic degrees of freedom. Here, we present two Heisenberg models exhibiting fractionalization in combination with exponentially decaying correlations. These can be thought of as a classical continuous spin version of a Z_{2} spin liquid. Our work suggests a systematic search and classification of classical spin liquids as a worthwhile endeavor.

Energy repartition in the nonequilibrium steady state

The concept of temperature in nonequilibrium thermodynamics is an outstanding theoretical issue. We propose an energy repartition principle that leads to a spectral (mode-dependent) temperature in steady-state nonequilibrium systems. The general concepts are illustrated by analytic solutions of the classical Heisenberg spin chain connected to Langevin heat reservoirs with arbitrary temperature profiles. Gradients of external magnetic fields are shown to localize spin waves in a Wannier-Zeemann fashion, while magnon interactions renormalize the spectral temperature. Our generic results are applicable to other thermodynamic systems such as Newtonian liquids, elastic solids, and Josephson junctions.