Metastability and dynamics in remanent states of square artificial spin ice with long-range dipole interactions

This paper precisely determines the metastable remanent states and analyzes the oscillation modes of square artificial spin ice by incorporating all long-range dipole interactions, ultimately identifying stability limits based on the ratio of local anisotropy to dipolar coupling strength.

The Gist

Imagine a giant, microscopic city made of tiny, elongated magnets (like little bar magnets) arranged in a perfect grid. This is called Artificial Square Spin Ice.

In this city, every magnet wants to point in a specific direction, but they are also neighbors who don't want to bump into each other. This creates a "frustrated" situation, like a group of friends trying to decide where to sit at a round table where everyone has a different preference.

The Story of the "Remanent State"

Usually, scientists try to get these magnets to settle into their absolute perfect, lowest-energy arrangement (the "Ground State"). But it's really hard to get them there because they get stuck in traffic jams of their own making.

However, there's an easier way to get them to stop moving. Imagine you have a giant magnet that you slowly pull away from the city. As you pull it away, the little magnets in the city align with it. When you finally remove the big magnet completely, the little magnets don't go back to their perfect, chaotic ground state. Instead, they get stuck in a Remanent State.

Think of it like a crowd of people who were marching in a straight line. When the leader stops, the crowd doesn't instantly scatter into a random mess; they stay in a loose, organized line for a while. This "stuck" line is the Remanent State. It's not the most perfect arrangement possible, but it's stable enough to stay there for a long time.

The Big Question: Is it Stable?

The author of this paper, G. M. Wysin, asks a crucial question: "If we nudge these magnets slightly, will they snap back into place, or will the whole line collapse?"

To answer this, he treats the magnets not as rigid arrows that can only point North or South (like an on/off switch), but as flexible darts that can wiggle and tilt slightly. He then calculates the "vibrations" or "wiggles" of these magnets.

The Analogy of the Trampoline and the Springs

Imagine the Remanent State is a trampoline.

• The Ground State is a trampoline lying perfectly flat on the ground.
• The Remanent State is a trampoline that has been pulled up into a specific shape and held there by a few people (the magnetic field).

The paper analyzes what happens if you jump on this trampoline.

1. The Springs (Anisotropy): The magnets have their own internal "springs" that want to keep them pointing along their long axis. This is like the tension in the trampoline fabric.
2. The Neighbors (Dipole Interactions): The magnets also push and pull on each other from a distance. This is like people holding hands while jumping on the trampoline.

The Two Models: The "Neighborly" vs. The "Global" View

The paper compares two ways of looking at these interactions:

1. The "Neighborly" Model (Short-Range):
Imagine you only listen to the people standing immediately next to you.

• Result: If the internal springs (anisotropy) aren't strong enough, the trampoline collapses. The magnets need to be very "stubborn" (strong internal preference) to stay in this Remanent State. If they are too weak, the whole structure becomes unstable and falls apart.

2. The "Global" Model (Long-Range):
Now, imagine everyone in the city can feel the pull of everyone else, even those far away.

• Result: This changes everything! The paper finds that when you include these long-distance "whispers" between magnets, the Remanent State becomes much more stable.
• The Surprise: Even if the internal springs are very weak (the magnets are very "flexible"), the long-distance cooperation between all the magnets holds the structure together. It's like a crowd of people holding hands across a huge field; even if individual people are weak, the collective hold is strong.

The "Wiggle" Frequencies (Magnons)

The paper also calculates the specific "notes" or frequencies at which these magnets vibrate.

• Think of the Remanent State as a musical instrument.
• If the instrument is stable, it hums with a clear, high-pitched tone.
• If it's about to collapse (unstable), the tone drops to zero, and the instrument goes silent (or breaks).

The author found that including the long-range interactions changes the "music" significantly. It raises the pitch (frequency) of the vibrations and prevents the instrument from going silent, meaning the state is safer and more robust.

The Real-World Takeaway

Why does this matter?
Scientists build these artificial magnetic cities to study how information can be stored or processed. If a state is unstable, the data is lost.

This paper tells us: Don't worry too much about the magnets being perfectly rigid. Even if they are a bit wobbly, the fact that they "talk" to each other over long distances helps keep the data safe. The Remanent State is a sturdy, metastable home for magnetic information, and it's easier to maintain than we previously thought.

In summary: The paper uses math to prove that a specific, slightly "imperfect" arrangement of tiny magnets is actually very stable because the magnets help each other out from far away, acting like a giant, cooperative safety net.

Technical Summary

1. Problem Statement

Artificial square spin ice (SASI) systems exhibit geometric frustration, leading to complex ground states and metastable configurations. While protocols exist to drive these systems toward their ground states (satisfying the "two-in, two-out" ice rule), experimental protocols often leave the system in a remanent state (RS) after an applied magnetic field is slowly removed.

• The Challenge: These RS configurations are metastable local energy minima with non-zero residual magnetization and excess energy compared to the ground state.
• The Gap: Previous theoretical analyses often relied on nearest-neighbor (NN) approximations or Ising spin models (fixed direction). However, real magnetic islands have finite shape anisotropy and experience long-range dipole interactions (LRD). It is unclear how LRD affects the stability limits and the spin-wave spectrum (magnon modes) of these remanent states, particularly when the island dipoles are allowed to tilt (Heisenberg-like dynamics) rather than being strictly pinned.

2. Methodology

The author employs a classical Heisenberg-like dipole model to analyze the system.

• Model System: A square lattice of magnetic islands with fixed dipole magnitude $\mu$ but variable direction. The islands possess uniaxial anisotropy ($K_1$) along their long axes and planar anisotropy ($K_3$) due to their thin geometry.
• Hamiltonian: The total energy includes:
  1. Dipole-Dipole Interactions: Summed over all pairs of islands, extending from nearest neighbors (NN) to infinite range.
  2. Anisotropy Energy: Quadratic approximations for uniaxial ($K_1$) and planar ($K_3$) terms.
• Analytical Approach:
  1. Equilibrium Determination: The remanent state is defined as the local energy minimum where sublattices A and B tilt inward toward the magnetization direction ($\phi_B = -\phi_A$).
  2. Linearization: Small deviations from equilibrium are analyzed using quadratic expansions of the Hamiltonian ($H^{(2)}$) for in-plane ($\phi$) and out-of-plane ($\theta$) fluctuations.
  3. Dynamics: The equations of motion are derived from the torque equation (Landau-Lifshitz without damping), leading to a matrix eigenvalue problem.
  4. Comparison: The study compares three models:
     • NN-only model.
     • Model including up to 2nd nearest neighbors.
     • Model including infinite-range dipole interactions (LRD).

3. Key Contributions

• Inclusion of Long-Range Interactions: The paper provides a rigorous analytical framework for summing dipole interactions over an infinite lattice, classifying bonds into intra-sublattice (AA/BB) and inter-sublattice (AB) types.
• Heisenberg Dynamics: Unlike Ising models, this work accounts for the continuous tilting of dipoles away from the island axes, which is crucial for accurately determining the equilibrium angles and stability limits.
• Stability Criteria Derivation: The author derives specific stability thresholds for the uniaxial anisotropy ($K_1$) required to maintain a metastable remanent state, comparing NN vs. LRD models.
• Dispersion Relations: The calculation of full magnon mode dispersion relations ($\omega$ vs. wave vector $q$) for the remanent state, identifying soft modes that signal instability.

4. Key Results

A. Equilibrium State and Tilting

• In the remanent state, the two sublattices tilt inward toward the magnetization direction (along the $[10]$ island axis).
• Effect of LRD: Including infinite-range interactions increases the inward tilting angle ($\phi_A$) compared to the NN model. This is primarily driven by 4th nearest-neighbor interactions (AB bonds).
• Energy: The inclusion of LRD lowers the energy of the remanent state, making it more favorable, though it remains metastable relative to the ground state.

B. Stability Limits

The stability of the remanent state depends on the ratio of uniaxial anisotropy to dipolar coupling ($K_1/D$).

• NN-Model: Stability requires $K_1 > K_{1,min} \approx 2.947 D$. Below this, in-plane modes at wave vectors $(\pi, 0)$ and $(0, \pi)$ go to zero frequency, causing instability.
• LRD-Model: The stability threshold is significantly reduced to $K_1 > 1.094 D$.
  • Significance: Long-range dipole interactions effectively "stiffen" the system, allowing the remanent state to remain stable even with much weaker uniaxial anisotropy.
  • The instability mechanism involves the softening of in-plane spin-wave modes to zero frequency at the Brillouin zone boundary.

C. Spin-Wave Spectrum (Magnon Modes)

• Mode Splitting: The spectrum consists of two branches ($\omega^+$ and $\omega^-$).
• Directional Dependence:
  • Along the magnetization direction ($[10]$), the modes can touch or cross, exhibiting linear dispersion near $q=0$.
  • Along the perpendicular direction ($[01]$) and diagonal ($[11]$), the modes remain widely separated.
• LRD Impact: Including LRD shifts the entire spectrum to higher frequencies and alters the dispersion shape. Specifically, the crossing point of modes along the $[10]$ direction shifts to lower wave vectors.
• Realistic Parameters: For a realistic Permalloy-based spin ice (with $K_1 \approx 38D$ and $K_3 \approx 84D$), the system is deeply stable, and the mode frequencies are high, confirming that thermal fluctuations at room temperature would not easily destroy the remanent state.

5. Significance

• Experimental Relevance: The calculated spin-wave spectra provide a "fingerprint" for identifying remanent states in experimental settings (e.g., via Brillouin Light Scattering or Ferromagnetic Resonance). The distinct frequencies and dispersion relations differentiate the RS from the ground state.
• Theoretical Advancement: The work demonstrates that neglecting long-range dipole interactions leads to an overestimation of the anisotropy required for stability. It validates the necessity of using Heisenberg-like models over Ising approximations for accurate dynamic predictions in artificial spin ice.
• Metastability Control: The findings suggest that long-range interactions play a stabilizing role, which is critical for designing spin-ice devices where specific metastable states are desired for information storage or logic operations.

In summary, this paper establishes that while square artificial spin ice remanent states are inherently metastable, long-range dipole interactions significantly enhance their stability and fundamentally alter their dynamic excitation spectra, a factor that must be accounted for in both theoretical modeling and experimental characterization.