Non-equilibrium formulation of helicity-dependent thermal field for ultrafast magnetization dynamics

This paper proposes a novel non-equilibrium thermal field model based on atomic spin flip probabilities that successfully reproduces ultrafast demagnetization across various grid sizes, overcoming the limitations of traditional micromagnetic models.

The Gist

Imagine you have a giant crowd of people (atoms) in a stadium, all holding flags (spins) pointing in the same direction. This represents a magnet. Usually, if you want to change the direction of all those flags, you have to push them one by one, which takes time.

But what happens if you hit that stadium with a super-fast, intense flash of light (a femtosecond laser)? The light is so fast that it doesn't just push the flags; it makes the people spin, stumble, and flip their flags chaotically in a split second. This is called ultrafast demagnetization.

Scientists have been trying to simulate this on computers for years, but they hit a wall. Here is the problem:

• The Atomistic View: If you simulate every single person in the stadium, the computer works perfectly but takes forever to run.
• The Micromagnetic View: To save time, scientists group people into "blocks" (cells) and treat each block as one giant person. This is fast, but it fails to capture the chaos of the laser flash. It's like trying to predict a riot by only looking at the average mood of a neighborhood; you miss the individual punches.

The New Solution: A "Chaos Calculator"

This paper introduces a clever new way to simulate these fast flips without needing to track every single atom. The author, Ezio Iacocca, proposes a "Non-Equilibrium Thermal Field."

Here is how it works, using some everyday analogies:

1. The "Coin Flip" Probability

Imagine each block of atoms isn't just a solid rock, but a bag of coins. When the laser hits, it doesn't just heat the bag; it changes the odds of the coins flipping.

• If the laser is "right-handed" (circularly polarized), it makes it more likely for a "Heads" coin to flip to "Tails."
• The paper calculates the probability of these flips happening inside a block.

2. The "Weather Forecast" vs. The "Storm"

Traditional simulations treat heat like a gentle, random breeze (thermal noise) that blows equally in all directions. It's like a calm day where the wind speed is the same everywhere.

This new method realizes that during a laser flash, the "wind" is actually a storm.

• Because the laser forces specific flips, the "wind" has a direction (it pushes the flags toward a new state).
• It also has a variability (some blocks flip more than others).
• The paper creates a mathematical "storm model" that tells the computer: "At this exact moment, the wind is pushing this hard in this specific direction, with this much chaos."

3. The "Energy Budget"

The author calculates how much energy is released when these spins flip. Think of it like a bank account.

• Every time a spin flips, it spends a tiny bit of "angular momentum" currency.
• The paper sums up all these tiny transactions in a block to figure out the total "energy bill."
• This energy bill is then converted into a massive, temporary "temperature" (thousands of degrees!) that drives the simulation.

Why is this a Big Deal?

It's Grid-Independent:
In the old "Micromagnetic" method, if you changed the size of your simulation blocks (the grid), the results would change completely. It was like measuring a room with a ruler that stretched or shrank depending on how you held it.
This new method works regardless of the block size. Whether you simulate a tiny speck or a large chunk of material, the physics remains consistent.

It Bridges the Gap:
This approach allows scientists to simulate huge magnetic structures (like the ones used in real hard drives) with the accuracy of a tiny, atom-by-atom simulation. It's like being able to predict the traffic flow of an entire city using a model that understands the behavior of a single driver.

The Bottom Line

The paper says: "We found a way to mathematically describe the chaotic, super-fast flipping of spins caused by a laser. By treating these flips as a probability game with a specific energy cost, we can simulate ultrafast magnetism on a computer without needing to track every single atom, making our models faster, more accurate, and ready for real-world applications."

In short, they built a universal translator that lets fast computers understand the chaotic language of ultrafast lasers.

Technical Summary

Here is a detailed technical summary of the paper "Non-equilibrium formulation of helicity-dependent thermal field for ultrafast magnetization dynamics" by Ezio Iacocca.

1. Problem Statement

The paper addresses a critical limitation in modeling ultrafast magnetization dynamics (occurring on femtosecond to picosecond timescales) using micromagnetic simulations.

• The Gap: Traditional micromagnetic models discretize magnetic materials into cells (typically hundreds of atoms wide) and apply a standard thermal field derived from the fluctuation-dissipation theorem. This approach assumes thermal equilibrium and uncorrelated noise.
• The Failure: While atomistic spin dynamics (ASD) can accurately reproduce ultrafast demagnetization driven by optical lasers or electric currents, standard micromagnetic models fail to do so in a grid-independent manner. As cell sizes increase, the impact of thermal fluctuations diminishes artificially, preventing the simulation of phenomena like domain patterns or solitons under ultrafast excitation without resorting to computationally expensive atomistic scales.
• The Need: A method is required to bridge the gap between atomistic physics and continuum micromagnetics that accounts for non-equilibrium effects (specifically spin flips induced by helicity-dependent lasers) in a way that is independent of the simulation grid size.

2. Methodology

The author proposes a non-equilibrium thermal field ($H_{n-e}$) formulation based on statistical spin-flip probabilities rather than equilibrium thermodynamics.

• Statistical Model of Spin Flips:
  • The model considers a micromagnetic cell containing $N$ spins.
  • It assumes a left circularly polarized laser pulse with helicity $\sigma$ creates a differential probability for spin-up ($P$) and spin-down ($B$) switching, where $P > B$.
  • The final magnetization state is treated as a probability distribution resulting from $n$ up-spin flips and $m$ down-spin flips.
• Derivation of Field Statistics:
  • Using the Central Limit Theorem (for large $N$), the distribution of the final magnetization is approximated as Gaussian.
  • The non-equilibrium field direction is defined by a unit vector $\hat{h}_{n-e}$ where the $z$-component ($\eta$) follows a Gaussian distribution with a mean ($\mu_{PB}$) and standard deviation ($\sigma_{PB}$) derived from the spin-flip probabilities.
  • Unlike traditional thermal fields (uniformly distributed on a sphere), this field has a specific mean direction and variance determined by the laser helicity and intensity.
• Energy and Magnitude Calculation:
  • The magnitude of the field is linked to the magnon energy ($\Delta E$) generated by the net change in spin orientation during a time step $\delta t$.
  • $\Delta E$ is calculated based on the difference in the number of flipped spins relative to the background.
  • The field magnitude is derived via $H_{n-e} = \sqrt{\frac{2\alpha \Delta E}{\gamma \mu_0 M_s V \delta t}}$.
  • Key Insight: For small cells and short times, this results in effective temperatures on the order of 10,000 K and fields on the order of 10 Tesla, which are necessary to drive ultrafast dynamics.
• Numerical Implementation:
  • The approach is implemented within a Pseudospectral Landau-Lifshitz (PS-LL) model.
  • The effective field includes local anisotropy, exchange interactions (via Fourier transform), the traditional thermal field ($H_{th}$), and the new non-equilibrium field ($H_{n-e}$).
  • To handle computational complexity (large factorials in probability calculations), Stirling's approximation is used for $N > 100$.

3. Key Contributions

• Grid-Independent Formulation: The proposed non-equilibrium field allows micromagnetic simulations to reproduce ultrafast demagnetization rates consistently across different cell sizes (from atomic resolution up to $\sim 2$ nm), overcoming the "grid dependence" failure of standard thermal models.
• Quantification of Magnon Energy: The paper provides a method to quantify magnon energy as a function of cell size, enabling a consistent determination of the non-equilibrium field magnitude.
• Helicity Dependence: The model explicitly incorporates the effect of laser helicity (circular polarization) by defining differential switching probabilities for up and down spins, leading to a directional bias in the thermal field.
• Multiscale Bridging: It offers a pathway to simulate ultrafast dynamics in macroscopic sample dimensions (using micromagnetics) that match the results of atomistic spin dynamics, enabling the study of domain patterns and solitons under ultrafast excitation.

4. Results

The model was verified using FePt material parameters with various cell sizes (0.4 nm to 2 nm) and laser pulse profiles:

• Demagnetization Consistency: Simulations showed that with the $H_{n-e}$ correction, the rate of demagnetization was comparable across all cell sizes (0.5 nm, 1 nm, 2 nm). In contrast, standard thermal models showed negligible demagnetization for larger cells (e.g., 2 nm).
• Atomic Limit: At the atomic limit ($N=1$), the correction vanishes ($\Delta E = 0$), correctly reverting to the standard atomistic behavior where demagnetization is driven purely by uncorrelated noise.
• Laser Temperature and Pulse Width:
  • Increasing the laser peak temperature (and thus switching probability) increased the demagnetization depth.
  • Increasing the pulse duration (standard deviation $t_\sigma$) allowed the system more time to react, leading to full demagnetization and the formation of labyrinth domain patterns, consistent with experimental observations.
• Effective Temperatures: The model yielded equivalent temperatures in the range of thousands of Kelvin during the pulse, consistent with the energy scales required for ultrafast switching.

5. Significance

This work represents a significant step forward in computational micromagnetism:

• Efficiency: It enables the simulation of ultrafast phenomena in large, experimentally relevant volumes without the prohibitive computational cost of full atomistic simulations.
• Physical Accuracy: It resolves the discrepancy where micromagnetic models failed to capture ultrafast dynamics, providing a physically grounded mechanism (non-equilibrium spin flips) rather than an ad-hoc adjustment.
• Future Applications: The framework is extendable to ferrimagnets, multi-pulse switching scenarios, and the investigation of complex magnetic textures (solitons, domain walls) under non-equilibrium conditions. It paves the way for designing and analyzing next-generation ultrafast magnetic storage and logic devices.