Micromagnetic simulations for magnetic multipoles

Résumé technique : Simulations micromagnétiques pour les multipôles magnétiques

Énoncé du problème
Les multipôles magnétiques de grappes (clusters) servent de paramètres d'ordre critiques caractérisant la symétrie de spin dans les matériaux magnétiques avancés, particulièrement les altermagnets et les antiferromagnétiques non colinéaires (par exemple, Mn3Sn). Bien que ces multipôles d'ordre élevé régissent des phénomènes électriques et optiques uniques, leur comportement aux échelles de longueur mésoscopiques reste largement inexploré. Les modèles analytiques existants capturent certains aspects du mouvement des parois de domaines, mais manquent de la résolution spatiale requise pour les dynamiques complexes, intrinsèquement bidimensionnelles et tridimensionnelles. Inversement, les modèles de spins atomistiques, bien que détaillés, sont trop coûteux en termes de calcul pour simuler des textures mésoscopiques à l'échelle du micromètre. Il existe actuellement un manque de cadre micromagnétique général et complet capable d'analyser quantitativement les systèmes de multipôles non uniformes, spécifiquement ceux où l'état de spin au sein d'une maille élémentaire est décrit par un vecteur représentatif unique.

Méthodologie
Les auteurs développent un formalisme micromagnétique pour les systèmes magnétiques non triviaux contenant des multipôles magnétiques de grappes. L'hypothèse centrale est que les orientations relatives des moments magnétiques au sein d'une maille élémentaire sont rigidement préservées, permettant de caractériser l'état de spin collectif par un champ de vecteurs représentatifs variant de manière lisse, $\mathbf{g}(\mathbf{r}, t)$.

• Cadre théorique : En partant d'une densité lagrangienne et d'une densité de fonction de dissipation de Rayleigh pour les dipôles magnétiques individuels, les auteurs dérivent une équation du mouvement pour le vecteur représentatif $\mathbf{g}$. Cette équation sert de contrepartie multipolaire à l'équation de Landau-Lifshitz pour les ferromagnétiques.
• Discrétisation : Pour la mise en œuvre numérique, le continuum est discrétisé en volumes micromagnétiques ($V_c$), chacun contenant suffisamment de mailles élémentaires pour que le champ vectoriel soit presque uniforme. Un vecteur représentatif $\mathbf{G}$ est défini comme la moyenne volumique de $\mathbf{g}$.
• Étude de cas (Mn3Sn) : Le cadre est appliqué au Mn3Sn, un antiferromagnétique non colinéaire possédant un réseau de kagomé et un ordre octupolaire magnétique. Le modèle incorpore l'énergie de Zeeman, une anisotropie à sixité ($K_6$) et une rigidité d'échange ($A_{ex}$), incluant les interactions de Dzyaloshinskii-Moriya.
• Validation : Le modèle micromagnétique est validé par la simulation de la dynamique d'un octupôle unique et de profils de parois de domaines à 60° unidimensionnels, en comparant directement les résultats avec des calculs de modèles de spins atomistiques.

Contributions clés

1. Formalisme unifié : L'article introduit un cadre micromagnétique général pour les multipôles magnétiques de grappes de type vectoriel, comblant le fossé entre la physique à l'échelle atomique et les phénomènes à l'échelle des dispositifs pour les antiferromagnétiques non colinéaires.
2. Validation : Le modèle est rigoureusement validé par des simulations de spins atomistiques, démontrant des résultats quasi identiques pour la dynamique d'octupôle uniforme (pilotée par champ ou par couple de spin) et pour les profils statiques de parois de domaines à 60° (de type Néel et de type Bloch).
3. Dynamique mésoscopique : Le cadre permet la simulation de phénomènes hors de portée des modèles atomistiques, spécifiquement le mouvement piloté par champ de parois de domaines à 180° complexes composées de trois segments de 60° consécutifs.

Résultats

• Dynamique uniforme : Le modèle micromagnétique reproduit avec précision l'évolution temporelle de l'angle de l'octupôle ($\Phi_G$) et les changements d'énergie magnétique totale sous champs magnétiques externes et couples de spin-orbite, correspondant aux résultats des spins atomistiques.
• Parois de domaines à 60° : Les simulations révèlent que les parois de domaines à 60° dans le Mn3Sn présentent de légères différences structurelles selon la direction de propagation (type Néel vs type Bloch), avec des largeurs calculées d'environ 145 nm et 128 nm, respectivement.
• Mouvement de paroi de domaine à 180° : Lorsqu'un champ magnétique pilote une paroi de domaine à 180° (composée de trois segments de 60°), la paroi subit une déformation de profil significative. Parce que les parois de 60° constituantes se déplacent à des vitesses différentes (la paroi centrale $[+30^\circ, -30^\circ]$ se déplace environ deux fois plus vite que les autres), la paroi de 180° se comprime et se déforme pendant le mouvement.
• Masse inertielle effective : La déformation de la paroi de domaine augmente son énergie magnétique ($E_{def}$). Les auteurs constatent qu'à basses vitesses, $E_{def}$ suit une loi quadratique par rapport à la vitesse ($v_{dw}^2$). Cela permet de définir une masse inertielle effective pour la paroi de domaine d'octupôle.
  • Masse calculée pour la paroi de 180° de type Néel : $M \approx 8,8 \times 10^{-28}$ kg.
  • Masse calculée pour la paroi de type Bloch : $M \approx 1,7 \times 10^{-27}$ kg.
  • Ces valeurs sont comparables à la densité de masse des parois de domaines ferromagnétiques dans le Permalloy.

Signification
L'article affirme que ce travail fournit une approche unifiée pour étudier la dynamique mésoscopique des multipôles de grappes de haute intensité. En établissant un cadre micromagnétique efficace sur le plan computationnel, les auteurs permettent l'analyse quantitative et spatialement résolue des systèmes de multipôles non uniformes. Cette capacité est présentée comme essentielle pour comprendre et concevoir les propriétés physiques des matériaux magnétiques fonctionnels, tels que les altermagnets et les antiferromagnétiques non colinéaires, qui sont des candidats pour les technologies spintroniques de prochaine génération. La découverte d'une masse inertielle émergente dans ces systèmes suggère que la dynamique de spin inertielle peut jouer un rôle important dans le basculement (switching) et la dynamique de spin, de manière similaire aux effets observés dans les antiferromagnétiques colinéaires.