Micromagnetic simulations for magnetic multipoles

Resumo Técnico: Simulações Micromagnéticas para Multipolos Magnéticos

Definição do Problema
Os multipolos magnéticos de cluster servem como parâmetros de ordem críticos que caracterizam a simetria de spin em materiais magnéticos avançados, particularmente altermagnetos e antiferromagnetos não colineares (ex: Mn3Sn). Embora esses multipolos de ordem superior governem fenômenos elétricos e ópticos únicos, seu comportamento em escalas de comprimento mesoscópicas permanece amplamente inexplorado. Modelos analíticos existentes capturam certos aspectos do movimento de paredes de domínio, mas carecem da resolução espacial necessária para dinâmicas complexas, inerentemente duas e três dimensais. Por outro lado, modelos de spin atomísticos, embora detalhados, são computacionalmente proibitivos para simular texturas mesoscópicas em escalas de micrômetros. Existe atualmente uma falta de um arcabouço micromagnético geral e abrangente capaz de analisar quantitativamente sistemas de multipolos não uniformes, especificamente aqueles onde o estado de spin dentro de uma célula unitária é descrito por um único vetor representativo.

Metodologia
Os autores desenvolvem um formalismo micromagnético para sistemas magnéticos não triviais contendo multipolos magnéticos de cluster. A premissa central é que as orientações relativas dos momentos magnéticos dentro de uma célula unitária são preservadas rigidamente, permitindo que o estado de spin coletivo seja caracterizado por um campo de vetor representativo que varia suavemente, $\mathbf{g}(\mathbf{r}, t)$.

• Estrutura Teórica: Partindo de uma densidade de Lagrangiana e de uma densidade de função de dissipação de Rayleigh para os dipolos magnéticos individuais, os autores derivam uma equação de movimento para o vetor representativo $\mathbf{g}$. Esta equação serve como o contraparte de multipolo para a equação de Landau-Lifshitz para ferromagnetos.
• Discretização: Para implementação numérica, o contínuo é discretizado em volumes micromagnéticos ($V_c$), cada um contendo células unitárias suficientes para que o campo vetorial seja quase uniforme. Um vetor representativo $\mathbf{G}$ é definido como a média volumétrica de $\mathbf{g}$.
• Estudo de Caso (Mn3Sn): O arcabouço é aplicado ao Mn3Sn, um antiferromagneto não colinear com uma rede de kagome e ordem de octupolo magnético. O modelo incorpora energia Zeeman, anisotropia de seis polos ($K_6$) e rigidez de troca ($A_{ex}$), incluindo interações de Dzyaloshinskii-Moriya.
• Validação: O modelo micromagnético é validado simulando a dinâmica de um octupolo único e perfis de parede de domínio de 60° unidimensionais, comparando os resultados diretamente com cálculos de modelos de spin atomísticos.

Principais Contribuições

1. Formalismo Unificado: O artigo introduz um arcabouço micromagnético geral para multipolos magnéticos de cluster do tipo vetorial, preenchendo a lacuna entre a física de escala atômica e os fenômenos de escala de dispositivo para antiferromagnetos não colineares.
2. Validação: O modelo é rigorosamente validado contra simulações de spin atomístico, demonstrando resultados quase idênticos para a dinâmica de octupolo uniforme (dirigida por campo e por torque de spin) e perfis estáticos de parede de domínio de 60° (tanto do tipo Néel quanto do tipo Bloch).
3. Dinâmica Mesoscópica: O arcabouço permite a simulação de fenômenos além do alcance dos modelos atomísticos, especificamente o movimento de paredes de domínio de 180° complexas compostas por três paredes consecutivas de 60°.

Resultos

• Dinâmica Uniforme: O modelo micromagnético reproduz com precisão a evolução temporal do ângulo do octupolo ($\Phi_G$) e as mudanças na energia magnética total sob campos magnéticos externos e torques de spin-órbita, correspondendo aos resultados de spin atomístico.
• Paredes de Domínio de 60°: As simulações revelam que as paredes de domínio de 60° no Mn3Sn exibem pequenas diferenças estruturais dependendo da direção de propagação (tipo Néel vs. tipo Bloch), com larguras calculadas de aproximadamente 145 nm e 128 nm, respectivamente.
• Movimento de Parede de Domínio de 180°: Quando um campo magnético conduz uma parede de domínio de 180° (composta por três segmentos de 60°), a parede sofre uma deformação significativa de perfil. Como as paredes de 60° constituintes movem-se a velocidades diferentes (a parede central $[+30^\circ, -30^\circ]$ move-se aproximadamente duas vezes mais rápido que as outras), a parede de 180° comprime-se e deforma-se durante o movimento.
• Massa Inercial Efetiva: A deformação da parede de domínio aumenta sua energia magnética ($E_{def}$). Os autores descobrem que, em baixas velocidades, $E_{def}$ escala quadraticamente com a velocidade ($v_{dw}^2$). Isso permite a definição de uma massa inercial efetiva para a parede de domínio de octupolo.
  • Massa calculada para a parede de 180° do tipo Néel: $M \approx 8,8 \times 10^{-28}$ kg.
  • Massa calculada para a parede do tipo Bloch: $M \approx 1,7 \times 10^{-27}$ kg.
  • Estes valores são comparáveis à densidade de massa de paredes de domínio ferromagnéticas em Permalloy.

Significância
O artigo afirma que este trabalho fornece uma abordagem unificada para investigar a dinâmica mesoscópica de multipolos de cluster de alta ordem. Ao estabelecer um arcabouço micromagnético computacionalmente eficiente, os autores permitem a análise quantitativa e espacialmente resolvida de sistemas de multipolos não uniformes. Esta capacidade é apresentada como essencial para compreender e projetar as propriedades físicas de materiais magnéticos funcionais, como altermagnetos e antiferromagnetos não colineares, que são candidatos para tecnologias spintrônicas de próxima geração. A descoberta de massa inercial emergente nestes sistemas sugere que a dinâmica de spin inercial pode desempenhar um papel significativo na comutação e dinâmica de spin, de forma semelhante aos efeitos observados em antiferromagnetos colineares.