🔬 materials science

Dynamic Magnetic Pair-Density Function of a One-Dimensional Ferromagnet

This paper establishes the theoretical foundation for analyzing local spin dynamics in magnetic materials by deriving and validating the dynamic magnetic pair-density function (DymPDF) for a one-dimensional Heisenberg ferromagnet, demonstrating its ability to capture energy-dependent spin correlations and magnon-mode transitions through analytical derivation and SpinW simulations.

Original authors: Shin-ichi Shamoto

Published 2026-08-25

📖 3 min read☕ Coffee break read

CC BY 4.0

Original authors: Shin-ichi Shamoto

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

To understand how a solid material holds its shape or conducts electricity, scientists often look at the arrangement of its atoms. But in magnetic materials, there is a second, invisible layer of order: the way tiny atomic magnets, called spins, point in relation to one another. For decades, researchers have mapped these static arrangements, creating a kind of snapshot that shows where spins are located and how they align at a single moment in time. This technique, known as the magnetic pair distribution function, has been essential for studying disordered magnets where the spins do not line up in a perfect, repeating pattern. However, magnets are rarely still. Just as atoms vibrate and move, these spins fluctuate and change direction over time. To truly understand how a magnetic material behaves, scientists need to see not just where the spins are, but how they move and interact as energy flows through the system. This requires a method that can capture these fleeting changes in real space, bridging the gap between a static picture and the dynamic reality of a living magnetic system.

A researcher has now developed a new mathematical framework to do exactly this, creating what they call the dynamic magnetic pair-density function. By extending the existing methods for static magnets to include energy changes, they have created a tool that can track the motion of spins as they shift and oscillate. To test this new approach, the scientist applied it to a simple, one-dimensional chain of magnetic atoms, a system that acts as a clean laboratory for understanding more complex materials. They used advanced computer simulations to model how these spins behave when hit with a specific amount of energy, generating a theoretical map of the spin dynamics. The results were striking: the new method successfully predicted how the relationship between neighboring spins changes as the energy of the system varies. Most notably, the researcher found that the nature of the connection between these spins flips completely at a specific energy threshold. Below this point, the spins tend to align in one way, but as the energy increases past a certain halfway mark, the relationship reverses, causing the signal to change its sign. This reversal happens precisely when the energy reaches half of the maximum possible energy the spin waves in the chain can carry.

The researcher confirmed that their new mathematical description matches the computer simulations with high precision. They showed that the behavior of the spins is governed by two main factors: the density of available energy states for the spin waves and a specific geometric factor related to how the spins are arranged in space. By including realistic details such as the limits of the measuring equipment and the finite distance over which the spins influence each other, they were able to reproduce the complex patterns seen in the simulations. The study demonstrates that this new dynamic method can accurately describe the local behavior of spins, even in a system where the spins are constantly changing. This success suggests that the technique is robust enough to be used on more complicated and disordered magnetic materials, such as those found in nanomagnets or frustrated systems where spins struggle to find a stable arrangement. By providing a clear window into the real-space motion of spins, this work lays the groundwork for future experiments that could reveal how local dynamics drive the physical properties of magnetic materials, offering a new way to see the invisible machinery of magnetism in action.

Problem and Motivation
Fourier-transform analysis, specifically the Pair Distribution Function (PDF) method, is a well-established technique for studying local structures in liquids, amorphous materials, and disordered crystals using total neutron and X-ray scattering. While the magnetic counterpart, the magnetic Pair Distribution Function (mPDF), has been developed to investigate local magnetic structures in frustrated magnets, the dynamic extension of this formalism remains theoretically undeveloped. Although recent advances in high-intensity pulsed neutron and synchrotron X-ray sources have enabled the study of both static and dynamic correlations via Fourier transforms, a rigorous theoretical framework for the Dynamic Magnetic Pair-Density Function (DymPDF), denoted as $D_M(r, E)$ , has not yet been established. This paper addresses the need to derive the analytical expression for DymPDF by extending the static mPDF formalism to finite energy transfers, aiming to provide a theoretical foundation for investigating local spin dynamics in magnetic materials.

Methodology
The author first reviews the formalism of the static magnetic pair-distribution function (mPDF) in polar coordinates, which relates the static magnetic structure factor $S_M(Q)$ to real-space spin correlations. They then extend this formalism to the dynamic regime by incorporating finite energy transfer $E$ . The dynamic magnetic structure factor $S_M(Q, E)$ is normalized by the Bose factor $\langle n(E) + 1 \rangle$ and Fourier-transformed to obtain the real-space DymPDF, $f(r, E)$ .

To validate this derived formalism, the author applies it to a representative system: a one-dimensional (1D) Heisenberg ferromagnet. The theoretical derivation involves:

Hamiltonian Definition: Defining a spin Hamiltonian with nearest-neighbor exchange interactions ( $J_{ij} = -1$ meV) and spin $S=1$ .
Dispersion Relation: Utilizing the linear spin-wave theory (Holstein-Primakoff approximation) to express the magnon dispersion $E(k)$ .
Correlation Functions: Deriving energy-resolved spin-correlation tensors $C_{ij}^{\alpha\beta}(E)$ to calculate the coefficients $A_{ij}(E)$ and $B_{ij}(E)$ that govern the DymPDF.
Simulation and Comparison: Generating powder-averaged dynamic structure factors and corresponding DymPDF maps using the SpinW software package. The analytical results derived from the 1D ferromagnet model are compared against these SpinW simulations.
Experimental Realism: The model is further refined to account for finite instrumental resolution (modeled via Gaussian broadening), magnetic correlation length, and Fourier-termination effects (due to finite $Q_{max}$ ).

Key Contributions and Results
The primary contribution is the derivation of the analytical expression for the DymPDF, $D_M(r, E)$ , which links the dynamic magnetic structure factor to real-space spin correlations at specific energy transfers.

For the 1D Heisenberg ferromagnet, the analysis yields several specific results:

Analytical Form: The DymPDF is shown to be governed by both the one-dimensional magnon density of states and a real-space mode factor. The expression includes terms dependent on the wave vector $q(E)$ , which is determined by the energy $E$ .
Sign Reversal: A critical finding is that the nearest-neighbor spin-pair correlation, $D_M(a, E)$ , changes sign at a specific energy transition. This transition occurs at $E_t = 2|J_a|S$ , which corresponds to one-half of the maximum magnon energy ( $E_{max}/2$ ).
Quantitative Agreement: The analytical expression for the nearest-neighbor correlation was fitted to the simulation data. The fitting yielded a transition energy $E_t = 1.99(2)$ meV and a width $\Gamma = 0.6(7)$ meV, which are in excellent agreement with the simulation parameters ( $E_t = 2$ meV).
Real-Space Reproduction: The real-space DymPDF at low energy (1 meV) was successfully reproduced using a model that incorporates finite instrumental resolution ( $\sigma$ ), magnetic correlation length ( $\xi$ ), and Fourier-termination effects ( $K(r)$ ). The fitted parameters ( $Q_{max} \approx 4.85$ Å $^{-1}$ , $\sigma \approx 0.81$ Å, $\xi \approx 16$ Å) were consistent with the simulation inputs.

Significance
The paper establishes the theoretical foundation for DymPDF analysis. By demonstrating that the DymPDF can be analytically derived and that it accurately reflects the underlying magnon dispersion and spin correlations, the author confirms the validity of this approach. The results suggest that DymPDF provides a robust real-space method for investigating local spin dynamics. The author posits that this formalism is applicable to a range of complex magnetic systems, including nanomagnets, clustered magnets, frustrated magnets, and amorphous magnetic systems, offering a new perspective on local dynamics beyond traditional reciprocal-space analysis.

More like this