Micromagnetic simulations for magnetic multipoles

Sintesi Tecnica: Simulazioni Micromagnetiche per Multipoli Magnetici

Definizione del Problema
I multipoli magnetici a cluster fungono da parametri d'ordine critici che caratterizzano la simmetria di spin in materiali magnetici avanzati, in particolare negli altermagneti e negli antiferromagneti non collineari (es. Mn3Sn). Sebbene questi multipoli di ordine superiore governino fenomeni elettrici e ottici unici, il loro comportamento alle scale di lunghezza mesoscopiche rimane ampiamente inesplorato. I modelli analitici esistenti catturano alcuni aspetti del moto delle pareti di dominio, ma mancano della risoluzione spaziale necessaria per dinamiche intrinsecamente bidimensionali e tridimensionali complesse. Al contrario, i modelli di spin atomistico, pur essendo dettagliati, sono computazionalmente proibitivi per simulare texture mesoscopiche su scala micrometrica. Esiste attualmente una mancanza di un quadro micromagnetico generale e completo capace di analizzare quantitativamente sistemi multipolari non uniformi, specificamente quelli in cui lo stato di spin all'interno di una cella unitaria è descritto da un singolo vettore rappresentativo.

Metodologia
Gli autori sviluppano un formalismo micromagnetico per sistemi magnetici non banali contenenti cluster magnetici multipolari. Il presupposto fondamentale è che le orientazioni relative dei momenti magnetici all'interno di una cella unitaria siano preservate rigidamente, permettendo allo stato di spin collettivo di essere caratterizzato da un campo vettoriale rappresentativo che varia dolcemente, $\mathbf{g}(\mathbf{r}, t)$.

• Framework Teorico: Partendo da una densità lagrangiana e da una densità di funzione di dissipazione di Rayleigh per i singoli dipoli magnetici, gli autori derivano un'equazione del moto per il vettore rappresentativo $\mathbf{g}$. Questa equazione funge da controparte multipolare dell'equazione di Landau-Lifshitz per i ferromagneti.
• Discretizzazione: Per l'implementazione numerica, il continuo viene discretizzato in volumi micromagnetici ($V_c$), ciascuno contenente un numero sufficiente di celle unitarie affinché il campo vettoriale sia quasi uniforme. Un vettore rappresentativo $\mathbf{G}$ è definito come la media volumetrica di $\mathbf{g}$.
• Caso di Studio (Mn3Sn): Il framework viene applicato al Mn3Sn, un antiferromagnete non collineare con un reticolo di kagome e ordine ottupolare magnetico. Il modello incorpora l'energia Zeeman, l'anisotropia esafolia ($K_6$) e la rigidità di scambio ($A_{ex}$), inclusi le interazioni di Dzyaloshinskii-Moriya.
• Validazione: Il modello micromagnetico viene validato simulando la dinamica di un singolo ottupolo e profili di parete di dominio a 60° monodimensionali, confrontando i risultati direttamente con i calcoli del modello di spin atomistico.

Contributi Chiave

1. Formalismo Unificato: Il lavoro introduce un framework micromagnetico generale per i cluster magnetici multipolari vettoriali, colmando il divario tra la fisica su scala atomica e i fenomeni su scala di dispositivo per gli antiferromagneti non collineari.
2. Validazione: Il modello è rigorosamente validato rispetto alle simulazioni di spin atomistico, dimostrando risultati quasi identici per la dinamica ottupolare uniforme (indotta da campo e da spin-torque) e per i profili statici di parete di dominio a 60° (sia di tipo Néel che di tipo Bloch).
3. Dinamica Mesoscopica: Il framework consente la simulazione di fenomeni oltre la portata dei modelli atomistici, specificamente il moto indotto da campo di pareti di dominio a 180° complesse composte da tre pareti consecutive a 60°.

Risultati

• Dinamica Uniforme: Il modello micromagnetico riproduce accuratamente l'evoluzione temporale dell'angolo ottupolare ($\Phi_G$) e le variazioni dell'energia magnetica totale sotto campi magnetici esterni e spin-orbit torques, corrispondendo ai risultati dello spin atomistico.
• Pareti di Dominio a 60°: Le simulazioni rivelano che le pareti di dominio a 60° nel Mn3Sn mostrano lievi differenze strutturali a seconda della direzione di propagazione (tipo Néel vs tipo Bloch), con larghezze calcolate di circa 145 nm e 128 nm, rispettivamente.
• Moto della Parete di Dominio a 180°: Quando un campo magnetico guida una parete di dominio a 180° (composta da tre segmenti a 60°), la parete subisce una significativa deformazione del profilo. Poiché le pareti a 60° costituenti si muovono a velocità differenti (la parete centrale $[+30^\circ, -30^\circ]$ si muove circa due volte più velocemente delle altre), la parete a 180° si comprime e si deforma durante il moto.
• Massa Inerziale Effettiva: La deformazione della parete di dominio aumenta la sua energia magnetica ($E_{def}$). Gli autori riscontrano che, a basse velocità, $E_{def}$ scala quadraticamente con la velocità ($v_{dw}^2$). Ciò permette la definizione di una massa inerziale effettiva per la parete di dominio ottupolare.
  • Massa calcolata per la parete a 180° di tipo Néel: $M \approx 8.8 \times 10^{-28}$ kg.
  • Massa calcolata per la parete di tipo Bloch: $M \approx 1.7 \times 10^{-27}$ kg.
  • Questi valori sono comparabili alla densità di massa delle pareti di dominio ferromagnetiche nel Permalloy.

Significatività
Il lavoro sostiene che questa ricerca fornisce un approccio unificato per investigare la dinamica mesoscopica degli alti ordini di multipoli a cluster. Stabilendo un framework micromagnetico computazionalmente efficiente, gli autori consentono l'analisi quantitativa e spazialmente risolta di sistemi multipolari non uniformi. Questa capacità è presentata come essenziale per comprendere e ingegnerizzare le proprietà fisiche di materiali magnetici funzionali, come gli altermagneti e gli antiferromagneti non collineari, che sono candidati per le tecnologie spintroniche di prossima generazione. La scoperta di una massa inerziale emergente in questi sistemi suggerisce che la dinamica di spin inerziale possa giocare un ruolo significativo nello switching e nella dinamica di spin, in modo simile agli effetti osservati nei sistemi antiferromagneti collineari.