Longitudinal magnon transport properties in the easy-axis XXZ Heisenberg ferromagnet on the face-centered cubic lattice

Basic Information

• Paper ID: 2507.04082
• Title: Longitudinal magnon transport properties in the easy-axis XXZ Heisenberg ferromagnet on the face-centered cubic lattice
• Author: M. R. Parymuda (Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine)
• Classification: cond-mat.str-el (Condensed Matter Physics - Strongly Correlated Electron Systems)
• Publication Date: October 24, 2025 (arXiv v5)
• Paper Link: https://arxiv.org/abs/2507.04082

Abstract

This paper presents a detailed investigation of longitudinal magnon-mediated spin and thermal transport properties in the S=1/2 ferromagnetic XXZ model on the face-centered cubic (FCC) lattice with easy-axis exchange anisotropy (Δ>1) comprising four sublattices. The study employs linear spin-wave theory (LSWT) and the Kubo formalism to calculate longitudinal spin conductivity and thermal conductivity, both exhibiting activated temperature dependence in the low-temperature regime with determined magnetic field dependencies. The analysis demonstrates that the magnon energy gap is crucial for ensuring the convergence of these conductivities. Furthermore, by examining the ratio of thermal to spin conductivity, an analogue of the Wiedemann-Franz law for magnon transport is identified at low temperatures. Finally, the results are shown to be generalizable to systems with arbitrary spin.

Research Background and Motivation

1. Research Questions

This paper investigates magnon-mediated spin and thermal transport properties in easy-axis anisotropic ferromagnets on the FCC lattice, with particular focus on:

• Temperature dependence of longitudinal spin conductivity σ and thermal conductivity κ in the low-temperature regime
• Mechanisms of magnetic field regulation of transport properties
• Analogues of the Wiedemann-Franz law in magnon transport

2. Importance of the Problem

• Fundamental Scientific Significance: Spin-wave transport is an active research field in condensed matter physics, with important value for understanding the propagation mechanisms of spin excitations in magnetically ordered systems
• Technological Application Prospects: Magnetic insulators show significant advantages in magnon spintronics because the absence of free charge carriers reduces energy dissipation, allowing spin waves to propagate over long distances with small Gilbert damping parameters
• Material Relevance: Rare-earth chalcogenides such as EuO and EuS are FCC ferromagnets with low Gilbert damping, where spin-orbit coupling produces exchange anisotropy, making the XXZ model an appropriate theoretical description

3. Limitations of Existing Methods

• For isotropic Heisenberg models (Δ=1), gapless acoustic modes in the dispersion relation lead to divergence of transport coefficients
• In two-dimensional systems, even with finite gaps and small Gilbert damping, longitudinal transport coefficients diverge due to infrared singularities of the Bose function at zero frequency
• Lack of systematic studies of transport properties in anisotropic ferromagnets on three-dimensional FCC lattices

4. Research Motivation

• Explore how the anisotropy parameter Δ, lattice geometry, and external magnetic field collectively affect transport properties
• Establish a Wiedemann-Franz law for magnon transport, revealing universal relationships between spin and thermal conductivities
• Provide a theoretical prediction framework for transport properties in practical magnetic materials

Core Contributions

1. Systematic calculation of transport coefficients for easy-axis XXZ ferromagnets on FCC lattice: First detailed study of longitudinal spin conductivity and thermal conductivity in three-dimensional FCC structures containing four sublattices
2. Revelation of the critical role of magnon energy gap: Demonstrated the necessity of finite gap (Δ>1) for ensuring convergence of transport coefficients, with activated temperature dependence:
   • σ(T) ∝ T^(3/2)exp(-Δ_gp/k_BT)
   • κ(T) ∝ T^(5/2)exp(-Δ_gp/k_BT)
3. Establishment of Wiedemann-Franz law for magnon transport: Derived a universal Lorentz number L = 5/2(k_B/gμ_B)², independent of anisotropy parameter Δ and external magnetic field H
4. Analysis of magnetic field regulation of transport properties: Demonstrated how external magnetic field suppresses transport coefficients through gap adjustment and magnon density modulation
5. Provision of theoretical framework generalizable to arbitrary spin S: Proved that results can be extended to arbitrary spin systems through appropriate modification of Holstein-Primakoff transformation

Detailed Methodology

Task Definition

Calculate longitudinal spin conductivity σ(T,H) and thermal conductivity κ(T,H) for the S=1/2 easy-axis XXZ ferromagnet on three-dimensional FCC lattice, where:

• Input: Lattice structure, exchange interaction parameters (J₁, J₂), anisotropy parameter Δ, external magnetic field H, temperature T
• Output: Spin conductivity σ, thermal conductivity κ and their temperature and magnetic field dependencies
• Constraints: Low-temperature regime (k_BT ≪ Δ_gp), small Gilbert damping parameter α ≪ 1

Model Architecture

1. Lattice Structure and Hamiltonian

The FCC lattice contains four sublattices (α=1,2,3,4), with lattice positions defined as R_mα = R_m + r_α, where R_m = m_xe_x + m_ye_y + m_ze_z.

The XXZ exchange Hamiltonian:

H = J₁∑_{⟨mα;nβ⟩} S_mα·S_nβ + J₂∑_{⟨⟨mα;nβ⟩⟩} S_mα·S_nβ

where J₁, J₂ < 0 denote ferromagnetic exchange, with anisotropic scalar product defined as:

S_mα·S_nβ = (1/2)(S⁺_mαS⁻_nβ + S⁻_mαS⁺_nβ) + ΔS^z_mαS^z_nβ

External magnetic field term: -gμ_BH∑_m∑_α S^z_mα

2. Holstein-Primakoff Transformation

Maps spin operators to boson operators (leading order approximation):

S⁺_mα = √(2S)a_mα
S⁻_mα = √(2S)a†_mα  
S^z_mα = S - a†_mαa_mα

satisfying commutation relations: a_mα, a†_nβ = δ_mnδ_αβ

3. Momentum Space Diagonalization

Introduce Fourier transformation: a_mα = (1/√N)∑_q e^(iq·R_m)a_qα

The Hamiltonian in momentum space is written as a 4×4 matrix form:

H = H₀ + ∑_q (a†_q1, a†_q2, a†_q3, a†_q4) M (a_q1, a_q2, a_q3, a_q4)^T

Four dispersion relations are obtained by diagonalizing matrix M.

4. Magnon Dispersion Relations

In the long-wavelength limit (|q|→0):

• Two flat bands: ω_q1 = ω_q2 = 2|J₁|Δ + gμ_BH
• Two dispersive bands: ω_q± = Δ^±_gp + D±|q|²

where the gap:

Δ^±_gp = |J₂|Δ± + gμ_BH
Δ± = [Δ(3+8δ)-3]/2 ± (1/2)√C₁

Stiffness constants D± are given in Appendix A.

Technical Innovations

1. Critical Role of Energy Gap

Innovation: Explicitly demonstrated the necessity of magnon energy gap for convergence of transport coefficients

• For Δ=1 (isotropic Heisenberg model), gapless acoustic modes lead to σ,κ→∞
• For Δ>1, anisotropy automatically opens a gap, ensuring finite transport coefficients even at H=0

Physical Mechanism:

• Gap suppresses thermal excitation of low-energy magnons
• Bose distribution function n_B ≃ exp(-ω_qα/k_BT) provides exponential cutoff
• Magnon density N_mag ∝ T^(3/2)exp(-Δ_gp/k_BT)

2. Kubo Formalism Calculation

Method Advantages:

• Rigorously calculates transport coefficients through linear response theory
• Introduces Gilbert damping parameter α ensuring Green function convergence
• Utilizes analytic continuation from Matsubara to real frequencies

Transport Coefficient Expressions:

σ±(T) = L±₁₁
κ±(T) = (1/T)[L±₂₂ - (L±₁₂)²/L±₁₁]

where L±_αβ are transport coefficients defined through magnetization current and heat current correlation functions.

3. Low-Temperature Asymptotic Behavior

Key Technique: Utilize integral asymptotic expansion

∫^∞_Δ dx (dn_B/dx)[(x-Δ)^(μ/2)/x] ≃ -Γ(μ/2+1)(k_BT)^(μ/2)e^(-Δ/k_BT)/Δ

valid when Δ/k_BT ≫ 1.

Yields:

σ(T) ∝ T^(3/2)e^(-Δ_gp/k_BT)
κ(T) ∝ T^(5/2)e^(-Δ_gp/k_BT)

4. Wiedemann-Franz Law

Derivation:

L = κ/(σT) = (1/T²)[L₂₂/L₁₁ - (L₁₂/L₁₁)²]
  = [Γ(9/2)/Γ(5/2) - (Γ(7/2)/Γ(5/2))²](k_B/gμ_B)²
  = (5/2)(k_B/gμ_B)²

Universality: Lorentz number L is independent of Δ, H, temperature T, and material parameters, depending only on the g-factor.

Experimental Setup

Parameter Selection

This is a theoretical calculation study using the following model parameters:

• Exchange Interactions: |J₁| = |J| (energy scale), |J₂| = 0.8|J|, i.e., δ = |J₁|/|J₂| = 1.25
• Anisotropy Parameters: Δ = 1.0001, 1.00012, 1.00014, 1.00016, 1.00018 (slightly greater than 1 to open small gaps)
• Gilbert Damping: α = 0.01 (small damping approximation)
• g-factor: g = 2 (corresponding to S=1/2)
• Temperature Range: T ∈ 0, 2.5×10⁻⁴|J|/k_B (low-temperature regime)
• Magnetic Field Range: H ∈ 0, 7×10⁻⁴|J|/gμ_B

Computational Methods

1. Dispersion Relation Calculation: Diagonalize 4×4 dynamical matrix M to obtain four dispersion branches
2. Brillouin Zone Integration: Convert momentum summation to integration (1/N)∑_q → (1/(2π)³)∫d³q
3. Frequency Integration: Use delta function approximation to handle broadening introduced by Gilbert damping
4. Numerical Integration: Perform numerical integration for temperature and field-dependent integrals

Theoretical Framework Validity Conditions

• Low-Temperature Condition: k_BT ≪ Δ_gp, ensuring small magnon density
• Small Damping Approximation: α ≪ 1, LSWT neglecting magnon-magnon interactions is valid
• Long-Wavelength Approximation: Consider low-energy excitations near |q|→0

Experimental Results

Main Results

1. Temperature Dependence (H=0)

Spin Conductivity σ(T) (Figure 3(a)):

• σ=0 at T=0 (no thermal magnon excitations)
• Activated behavior in low-temperature regime: σ ∝ T^(3/2)exp(-Δ_gp/k_BT)
• σ decreases significantly with increasing Δ (larger gap)
• For Δ=1 (Heisenberg model), σ diverges

Quantitative Data (read from Figure 3(a), Δ=1.0001):

• T=1.0×10⁻⁴: σ ≈ 0.004 (units: 1/(aμ²_B))
• T=2.0×10⁻⁴: σ ≈ 0.008
• Shows approximately linear growth (exponential term dominates)

Thermal Conductivity κ(T) (Figure 3(b)):

• Similar activated behavior: κ ∝ T^(5/2)exp(-Δ_gp/k_BT)
• Higher power-law temperature dependence compared to σ(T)
• Also diverges for Δ=1

Quantitative Data (Δ=1.0001):

• T=1.0×10⁻⁴: κ ≈ 0.5×10⁻⁷ (units: 1/(ak_B|J|))
• T=1.6×10⁻⁴: κ ≈ 1.0×10⁻⁷

2. Magnetic Field Dependence (T=10⁻⁴ fixed)

Spin Conductivity σ(H) (Figure 4(a)):

• Maximum value at H=0
• Monotonically decreases with increasing H
• Approaches zero at high fields

Quantitative Data (Δ=1.0001):

• H=0: σ ≈ 0.9×10⁻⁴
• H=3×10⁻⁴: σ ≈ 0.45×10⁻⁴ (approximately halved)
• H=7×10⁻⁴: σ ≈ 0.1×10⁻⁴

Thermal Conductivity κ(H) (Figure 4(b)):

• Shows same qualitative behavior as σ(H)
• Numerically differs by factor (5/2)T

Quantitative Data (Δ=1.0001):

• H=0: κ ≈ 2.0×10⁸
• H=7×10⁻⁴: κ ≈ 0.2×10⁸ (approximately 10-fold reduction)

3. Dispersion Relations (Figure 2)

• Four dispersion branches along high-symmetry path Γ-M-A-R-Γ
• Upper and lower dispersive bands have energies approximately 4|J| and 2|J| at Γ point, respectively
• Two flat bands (do not contribute to transport)

Key Findings

1. Dual Role of Energy Gap

• Necessity: Ensures convergence of transport coefficients (Δ>1)
• Trade-off: Optimal transport requires Δ_gp ≪ 1 (gap as small as possible but nonzero)

2. Low-Energy Branch Dominance

At low temperatures, only the lower dispersive band (ω_q-) contributes to transport because:

• Δ⁺_gp ≫ Δ⁻_gp
• Upper branch Bose occupation is exponentially suppressed

3. Flat Bands Do Not Contribute

Although flat bands have high density of states, their group velocity v_g=0, thus contributing nothing to σ and κ.

4. Wiedemann-Franz Law

κ/(σT) = L = (5/2)(k_B/gμ_B)² ≈ 3.125×10⁻⁵ eV²/K² (for g=2)

This ratio is independent of Δ, H, and T, exhibiting universality.

5. Magnetic Field Regulation Mechanism

External magnetic field suppresses transport through two mechanisms:

• Increases gap: Δ_gp → Δ_gp + gμ_BH
• Reduces magnon density: N_mag ∝ exp(-gμ_BH/k_BT)

Comparison with Literature

Consistency with Ref.25:

• Simple cubic lattice antiferromagnetic and ferromagnetic Heisenberg models also yield L=5/2
• Confirms universality of WF law (independent of lattice structure)

Differences from Two-Dimensional Systems:

• Two-dimensional systems diverge even with finite gaps due to infrared divergence
• Three-dimensional systems have convergent phase space integral ∫q²dq

Related Work

1. Magnon Spintronics

• YIG Materials 6-8: Most famous low-damping magnetic insulators, spin waves propagate long distances at room temperature
• Rare-Earth Compounds 9,10: EuO and EuS are FCC ferromagnets, directly relevant to this paper's model
• Other Materials 11-13: MnFe₂O₄, Cr₂O₃, Cu₂OSeO₃ etc. exhibit rich magnon transport phenomena

2. Theoretical Methods

• Linear Spin-Wave Theory 14-19: Holstein-Primakoff transformation is standard tool for low-temperature magnetism
• Kubo Formalism 15,16: Rigorous framework of linear response theory for calculating transport coefficients
• Boltzmann Equation 25: Alternative method for deriving WF law, results consistent with Kubo formalism

3. Transport Property Studies

• Frequency Dependence 18: Frequency response of spin conductivity on FCC lattice
• Transverse Transport 34-36: Anomalous transport in systems with broken inversion symmetry
• Two-Dimensional Systems 34,35: Topological magnon transport on checkerboard lattice

4. Wiedemann-Franz Law

• Electronic Systems: Classical WF law L=(π²/3)(k_B/e)²
• Magnon Systems 24-27:
  • Ref.24 first proposed magnon WF law
  • Ref.25 derived L=5/2 using Boltzmann equation
  • Ref.27 studied breakdown in strongly nonlinear regime

5. Positioning of This Work

• Supplements Existing Research: Provides systematic analysis of longitudinal transport on FCC lattice
• Extends to Anisotropy: Studies Δ>1 case, not limited to Δ=1
• Emphasizes Gap Role: Explicitly points out necessity of gap for convergence
• Three-Dimensional vs. Two-Dimensional: Avoids infrared divergence problem in 2D systems

Conclusions and Discussion

Main Conclusions

1. Activated Transport Behavior: In easy-axis anisotropic (Δ>1) FCC ferromagnets, spin conductivity and thermal conductivity exhibit low-temperature behavior σ∝T^(3/2)exp(-Δ_gp/k_BT) and κ∝T^(5/2)exp(-Δ_gp/k_BT)
2. Critical Role of Energy Gap: Magnon energy gap is necessary condition for ensuring convergence of transport coefficients; for isotropic Heisenberg model (Δ=1), transport coefficients diverge
3. Universal Wiedemann-Franz Law: Established WF law for magnon transport with Lorentz number L=(5/2)(k_B/gμ_B)², independent of Δ, H, and T
4. Magnetic Field Regulation: External magnetic field suppresses transport by increasing gap and reducing magnon density, providing microscopic mechanism for magnetothermal effects
5. Generalizability: Theoretical framework can be extended to arbitrary spin S systems through modified HP transformation

Limitations

1. Theoretical Approximation Constraints

• Low-Temperature Limitation: LSWT valid only when k_BT≪Δ_gp; magnon-magnon interactions become important at elevated temperatures
• Small Damping Assumption: Condition α≪1 may not hold at high temperatures; actual α should have temperature dependence α(T)∝T² 37,38
• Leading-Order Approximation: Neglects higher-order terms in HP transformation, unable to describe magnon interaction effects

2. Model Simplifications

• Anisotropy Parameters: Assumes Δ₁=Δ₂=Δ (nearest and next-nearest neighbor identical), actual materials may be more complex
• Exchange Interactions: Considers only nearest and next-nearest neighbors, neglects longer-range interactions
• Topologically Trivial: Considers systems preserving inversion symmetry, does not involve topologically nontrivial phases

3. Missing Experimental Verification

• Paper is purely theoretical work lacking direct comparison with experimental data
• Predictions need verification in actual materials like EuO, EuS

4. Scope of Applicability

• Lattice Structure: Specific to FCC lattice; other lattice structures require recalculation
• Magnetic Order Type: Applies only to ferromagnets; antiferromagnets require different treatment
• Dimensional Effects: Three-dimensional results cannot be directly applied to two-dimensional systems

Future Directions

1. Theoretical Extensions

• Include Magnon Interactions: Consider higher-order HP terms to study transport at finite temperatures
• Temperature-Dependent Damping: Introduce α(T) for more accurate description of high-temperature behavior
• Nonlinear Transport: Study nonlinear effects under large temperature or magnetic field gradients

2. Material Applications

• Actual Material Parameters: Use experimental parameters from EuO, EuS etc. for quantitative predictions
• Experimental Comparison: Compare with spin Seebeck effect and magnetothermal effect experiments
• Material Design: Optimize Δ and J for achieving optimal transport properties

3. Extension to Other Systems

• Antiferromagnets: Study transport properties of antiferromagnetic XXZ model
• Other Lattices: Simple cubic, body-centered cubic, honeycomb lattices etc.
• Topological Systems: Introduce topological nontrivial phases through broken inversion symmetry

4. Many-Body Effects

• Magnon-Phonon Coupling: Study effects of lattice vibrations on magnon transport
• Magnon-Magnon Scattering: Many-body calculations beyond LSWT
• Quantum Fluctuations: Quantum effects in low-dimensional systems

In-Depth Evaluation

Strengths

1. Rigor of Methodology

• Complete Theoretical Framework: Derives transport coefficients from microscopic Hamiltonian through LSWT and Kubo formalism
• Detailed Mathematical Derivation: Appendices provide complete analytical expressions and integral calculations
• Clear Approximation Conditions: Explicitly states low-temperature, small-damping and other validity conditions

2. Deep Physical Insights

• Clarification of Gap Role: Explicitly shows necessity of gap for convergence and explains physical mechanism
• Universality of WF Law: Demonstrates Lorentz number independence from multiple parameters, revealing deep physical connections
• Low-Energy Branch Dominance: Explains why only lower dispersive band contributes to transport

3. Systematic Nature of Results

• Multi-Parameter Study: Examines effects of Δ, T, H
• Qualitative and Quantitative Combination: Provides both analytical asymptotic behavior and numerical results
• Clear Figures: Dispersion relations, temperature and field dependencies clearly presented

4. Academic Contributions

• Fills Research Gap: First systematic study of longitudinal transport on FCC lattice with easy-axis XXZ model
• Consistency with Literature: WF law results consistent with other methods (Boltzmann equation), enhancing credibility
• Strong Generalizability: Framework applicable to other lattices and arbitrary spins

Weaknesses

1. Theoretical Limitations

• Narrow Temperature Range: Only applicable to extremely low temperatures (T~10⁻⁴|J|/k_B), limiting practical applications
• Neglects Interactions: Leading-order LSWT approximation fails at higher magnon densities
• Simplified Damping Treatment: Constant α cannot describe temperature-dependent dissipation mechanisms

2. Weak Experimental Connection

• Lacks Material Parameters: Does not use experimental data from specific materials (e.g., EuO)
• No Experimental Verification: Pure theoretical predictions lacking experimental comparison
• Measurement Feasibility Not Discussed: Does not discuss how to experimentally measure predicted transport coefficients

3. Analysis Depth

• Insufficient Flat Band Discussion: Simply states they don't contribute, lacks deeper exploration of physical significance
• Missing Topological Properties: Does not consider possible topological phases and edge states
• Finite-Size Effects Absent: Does not discuss limitations from actual sample sizes

4. Technical Details

• Numerical Methods Underspecified: Specific numerical algorithms for integration not detailed
• Missing Error Analysis: No precision assessment for numerical results
• Parameter Selection Justification: Why choose |J₂|=0.8|J| not sufficiently explained

Impact

1. Academic Contributions

• Theoretical Framework: Provides standard theoretical description for transport in FCC magnetic materials
• Universal Laws: WF law establishment has fundamental scientific value
• Method Demonstration: LSWT+Kubo formalism combination applicable to other systems

2. Practical Value

• Material Screening: Theory predictions guide search for materials with optimal transport properties
• Device Design: Magnetic field regulation mechanism useful for spintronics devices
• Magnetothermal Applications: Provides microscopic understanding for magnetic cooling technology

3. Reproducibility

• High Reproducibility:
  • Clear derivation steps
  • Explicit parameter choices
  • Complete formulas in appendices
• Code Implementation: Independent verification possible from provided formulas
• Easy Generalization: Framework directly applicable to other parameters and lattices

4. Follow-Up Research

• Citation Potential: Serves as standard reference for FCC magnetic material transport
• Research Directions: Stimulates studies on other lattices, antiferromagnets, topological systems
• Experimental Motivation: Theory predictions encourage related experimental measurements

Applicable Scenarios

1. Material Types

• Rare-Earth Chalcogenides: EuO, EuS and other FCC ferromagnetic insulators
• Low-Damping Materials: Systems with Gilbert damping α<0.1
• Easy-Axis Anisotropy: Materials with Δ>1 (from spin-orbit coupling)

2. Temperature Conditions

• Ultra-Low Temperatures: T < 0.1Δ_gp/k_B, typically requiring T < 1K
• Dilute Magnon Gas: Magnon density ≪ lattice site density

3. Experimental Techniques

• Spin Seebeck Effect Measurement: Temperature gradient-driven spin current
• Magnetothermal Effects: Magnetic field-controlled entropy change and temperature variation
• Nonlocal Transport: Long-distance spin current transmission

4. Theoretical Research

• Benchmark Comparison: Serves as benchmark for more complex theories (including interactions)
• Parameter Fitting: Extract exchange parameters and anisotropy from experimental data
• Generalization Studies: Extension to other lattices, spins, dimensions

5. Inapplicable Scenarios

• High-Temperature Region: T~J/k_B where magnon interactions become important
• Two-Dimensional Systems: Suffer from infrared divergence
• High-Damping Materials: α>0.1 where LSWT fails
• Topologically Nontrivial Phases: Systems with broken inversion symmetry

Key References

14 T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 (1940) - Original HP transformation paper

15 K. Aoyama and H. Kawamura, Phys. Rev. B 100, 144416 (2019) - Transport properties of antiferromagnets on square lattice

16 K. Aoyama, Phys. Rev. B 106, 224407 (2022) - Critical behavior of three-dimensional antiferromagnets

24 K. Nakata et al., Phys. Rev. B 92, 134425 (2015) - First proposal of magnon WF law

25 K. Nakata et al., Phys. Rev. B 96, 224414 (2017) - Derivation of WF law using Boltzmann equation

37 S. Tyc and B. I. Halperin, Phys. Rev. B 42, 2096 (1990) - Spin-wave damping in two-dimensional antiferromagnets

Overall Assessment: This is a rigorous theoretical physics paper systematically investigating magnon transport properties in easy-axis XXZ ferromagnets on FCC lattice. Main contributions include establishing a complete theoretical framework, revealing the critical role of energy gap, and proving universality of WF law. The methodology combines LSWT and Kubo formalism with detailed derivations and reliable results. Limitations include applicability only to ultra-low temperatures and small damping, plus lack of direct experimental comparison. The paper provides important theoretical foundation for understanding transport mechanisms in magnetic insulators with potential impact on magnon spintronics field.