Physics-informed X-ray tomography of magnetization and material parameters

Resumen Técnico: Tomografía de Rayos X informada por la física de la magnetización y parámetros de materiales

Planteamiento del problema
Las texturas magnéticas tridimensionales, como las cuerdas de skyrmion, los hofniones y los anillos de vórtice, están estabilizadas por interacciones microscópicas dentro de los materiales anfitriones. Si bien la tomografía de vectores magnéticos por rayos X (utilizando el dicroísmo circular de rayos X, o XMCD) puede imágenes estas estructuras a escala nanométrica, la técnica enfrenta limitaciones significativas. El débil contraste magnético suele resultar en proyecciones ruidosas, y las restricciones prácticas sobre el tiempo de exposición en el sincrotrón, la dosis de radiación y los ángulos de medición accesibles conducen frecuentemente a datos incompletos (por ejemplo, cuñas o conos faltantes). En consecuencia, las reconstrucciones iterativas estándar suelen dejar las características a escala nanométrica mal restringidas o introducen artefactos, particularmente cuando la resolución experimental se aproxima a la escala de las características magnéticas. Además, los enfoques tradicionales que aplican modelos micromagnéticos solo después de la reconstrucción no pueden recuperar las características perdidas durante el proceso de imagen y dependen de un ajuste manual que consume mucho tiempo.

Metodología
Los autores proponen un marco de reconstrucción "informado por la física" que integra el funcional de energía micromagnética directamente en el proceso de optimización tomográfica. En lugar de tratar la reconstrucción como un problema inverso puramente matemático, el método minimiza una función objetivo conjunta $L(m)$:
$$L(m) = L^2_{\text{data}}(m) + \lambda_E E_{\text{total}}(m)$$
donde $L^2_{\text{data}}$ mide el desajuste entre las proyecciones simuladas y las medidas por XMCD, y $E_{\text{total}}$ representa la energía micromagnética total de la muestra. Este término de energía incluye la rigidez de intercambio ($A$), las interacciones de desmagnetización (escaladas por la magnetización de saturación $M_s$) y la interacción Dzyaloshinskii–Moriya (DMI, constante $D$).

La reconstrucción utiliza un optimizador basado en gradientes (Adam de Riemann) con diferenciación automática para actualizar el campo de vectores de magnetización $\mathbf{m}$ (donde $|\mathbf{m}|=1$). El modelo de propagación directa simula las integrales de línea de los componentes de la magnetización a lo largo de la dirección del haz para series de inclinación de doble eje. La energía micromagnética actúa como un prior físico, restringiendo la solución a estados físicamente plausibles sin desenfocar las paredes de dominio, ya que se preserva el funcional de energía completo (incluyendo las interacciones competitivas) en lugar de un simple término de suavizado.

Contribuciones clave y resultados
El método se aplicó a un disco de 1.1 µm de diámetro y 300 nm de espesor del imán quiral de temperatura ambiente $\text{Co}_8\text{Zn}_9\text{Mn}_3$, procesado mediante pticografía de rayos X blandos dicroicos.

1. Recuperación a partir de datos dispersos y de un solo eje:
   El estudio demuestra que el enfoque informado por la física puede recuperar la textura helicoidal 3D a partir de conjuntos de datos significativamente reducidos.

   • Muestreo disperso: Se lograron reconstrucciones precisas utilizando solo 9 proyecciones (cada 8ª proyección de un conjunto estándar de doble eje) en comparación con las 67 completas.
   • Tomografía de un solo eje: El método reconstruyó con éxito la magnetización 3D utilizando únicamente datos de un único eje de inclinación (ya sea la serie $xz$ o la $yz$). Aunque un solo eje no puede medir el componente de la magnetización paralelo al eje de rotación, la energía micromagnética acopla el componente no medido con los componentes medidos, permitiendo que el optimizador infiera la información faltante de manera consistente con la física del material.
   • Validación: Utilizando una verdad de terreno sintética generada a partir de los datos experimentales, los autores demostraron que, sin el prior físico, el error en los componentes en el plano era alto ($\epsilon_m = 0.74$). Con el prior físico equilibrado, el error cayó a $\epsilon_m = 0.11$, recuperando con precisión todos los tres componentes vectoriales.

2. Sensibilidad de los parámetros y lectura:
   Los autores investigaron si la reconstrucción podía determinar los parámetros del material ($A$, $D$, $M_s$) analizando el desajuste de los datos con pesos de alta energía ($\lambda_E$).

   • Sensibilidad: Se encontró que el desajuste de los datos es relativamente insensible a $M_s$, pero altamente sensible a la relación entre la rigidez de intercambio y la DMI ($A/D$).
   • Valor de lectura: Mediante la minimización del desajuste de los datos a través de una cuadrícula de parámetros, el método identificó una constante DMI óptima de $D = 0.66 \times 10^{-3} \text{ J m}^{-2}$ (asumiendo la $A_0$ de la literatura). Este valor es aproximadamente un 10% superior al valor de la literatura ($D_0 = 0.60 \times 10^{-3} \text{ J m}^{-2}$) derivado de la relación de la hélice en volumen ($L_D = 4\pi A/D$).
   • Explicación: La discrepancia se atribuye a la diferencia entre la relación de volumen (válida para cristales sin bordes) y la geometría confinada de la nanoestructura. El modelo informado por la física, que tiene en cuenta la energía completa del disco confinado, produce un conjunto de parámetros que reproduce mejor la textura medida que el valor derivado del volumen.

Significancia y afirmaciones
El artículo afirma que este enfoque ofrece un doble beneficio para la imagen magnética:

1. Eficiencia: Reduce el tiempo y el tiempo de exposición al haz requeridos para la tomografía vectorial al permitir reconstrucciones 3D robustas a partir de datos dispersos o adquisiciones de un solo eje, lo que facilita los estudios dependientes del campo y en el tiempo.
2. Caracterización de materiales: Extiende la utilidad de la tomografía más allá de la imagen hacia la "magnetometría impulsada por la textura", permitiendo la determinación de parámetros de interacción magnética fundamentales (específicamente la relación $A/D$) directamente a partir de la textura de la nanoestructura.

Los autores enfatizan que el método es modular; el prior físico puede combinarse con cualquier modelo de propagación directa diferenciable (por ejemplo, holografía electrónica o tomografía de Lorentz) y extenderse para incluir otras interacciones de materiales, siempre que el modelo de propagación directa sea diferenciable. El trabajo sugiere una estrategia donde una única medición bien muestreada de un estado de referencia puede validar el modelo de energía y fijar los parámetros, permitiendo que las condiciones subsecuentes se imaginen con adquisiciones dispersas.