🔬 materials science

Development of a magnetic interatomic potential for cubic anti-ferromagnets: the case of NiO

Este artículo presenta una metodología novedosa para desarrollar potenciales interatómicos magnéticos para antiferromagnetos cúbicos mediante la integración de los modelos de intercambio de Heisenberg y de Néel, aplicándola con éxito al NiO para crear potenciales validados que permiten simulaciones de dinámica molecular a gran escala de fenómenos magnetoelásticos acoplados.

Autores originales: Ievgeniia Korniienko, Pablo Nieves, Jakub Sebesta, Roberto Iglesias, Dominik Legut

Publicado 2026-07-20

📖 4 min de lectura☕ Lectura para el café

CC BY 4.0

Autores originales: Ievgeniia Korniienko, Pablo Nieves, Jakub Sebesta, Roberto Iglesias, Dominik Legut

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Imagina el mundo invisible dentro de un objeto sólido, como una pequeña y bulliciosa ciudad hecha de átomos. En esta ciudad, los átomos no están simplemente quietos; están constantemente agitándose, chocando entre sí y tomándose de las manos de formas específicas. Pero en algunos materiales especiales, estos átomos tienen un superpoder secreto: actúan como diminutos imanes. Normalmente, cuando pensamos en imanes, imaginamos un imán de nevera que se pega al metal. Pero en un tipo diferente de material llamado "antiferromagnético", los vecinos son como un juego de tirar de la cuerda donde todos tiran en direcciones opuestas. Un átomo apunta hacia "arriba", su vecino apunta hacia "abajo", y se cancelan entre sí para que todo el bloque no actúe como un imán hacia el mundo exterior.

Los científicos están muy interesados en estos materiales porque pueden cambiar de forma cuando se aplica un campo magnético, o pueden cambiar su magnetismo cuando se les aprieta. Esto se llama "magnetoelasticidad". Es como si apretaras una pelota antiestrés y, de repente, decidiera apuntar al norte, o si acercaras un imán y esta se aplastara. Para entender cómo ocurre esto, los investigadores utilizan potentes simulaciones por computadora. Construyen una "receta" llamada potencial interatómico, que es básicamente un conjunto de instrucciones que le dice a la computadora cómo deben empujar, tirar y reaccionar los átomos a las fuerzas magnéticas. La parte difícil es que, para estos materiales de tirar de la cuerda, a las recetas les ha faltado un ingrediente crucial: las reglas magnéticas. Sin ellas, las simulaciones por computadora son como intentar predecir el clima sin saber sobre el viento.

Este artículo trata de escribir finalmente ese libro de recetas que faltaba para un material específico llamado Óxido de Níquel (NiO), que es el "modelo a seguir" de estos imanes de tirar de la cuerda. Los investigadores, un equipo de científicos de la República Checa y España, desarrollaron una nueva forma de mezclar las reglas del magnetismo con las reglas de cómo se mueven y estiran los átomos. Crearon dos versiones diferentes de esta receta. Una versión se basa en una idea clásica sobre cómo se comportan los sólidos iónicos (como la sal), y la otra es una receta nueva y personalizada, construida desde cero utilizando cálculos computacionales de alto nivel.

El equipo probó sus nuevas recetas simulando una gran ciudad de más de 85.000 átomos. Comprobaron si el material simulado se encogería o expandiría la cantidad justa cuando el orden magnético cambiaba, y si reaccionaría a la presión de la manera correcta. Los resultados fueron prometedores: ambas recetas predijeron el comportamiento del material con una excelente precisión, coincidiendo con los complejos cálculos realizados por supercomputadoras. Específicamente, encontraron que el material se encoge aproximadamente un 0,14% cuando ordena sus espines magnéticos, y recrearon con éxito la resistencia del material a ser aplastado en diferentes direcciones.

Sin embargo, los autores toman nota cuidadosamente de que estos resultados son actualmente válidos para simulaciones a temperatura cero (un estado teórico donde los átomos dejan de agitarse). Si bien los modelos funcionan maravillosamente en la computadora, los científicos sugieren que se necesita más trabajo para ver cómo se comportan estos materiales en el mundo real y cálido. También señalan que, aunque su nueva receta personalizada (llamada RF-MEAM) es mejor para predecir cómo el material resiste la torsión que la receta antigua, todavía hay cierto misterio sobre qué tan fuerte es el estiramiento magnético del material en comparación con otros estudios.

En resumen, este artículo no pretende haber resuelto todos los misterios de los materiales magnéticos, pero ha construido un puente muy fuerte y confiable para futuros exploradores. Al combinar con éxito las reglas magnéticas con las reglas de movimiento para el Óxido de Níquel, han abierto la puerta a simulaciones mucho más grandes y detalladas. Esto significa que, en el futuro, los científicos podrían ser capaces de diseñar mejores sensores, actuadores o incluso tecnologías cuánticas que utilicen estos ingeniosos imanes que cambian de forma, todo porque finalmente tienen un buen mapa de cómo danzan juntos los átomos.

Planteamiento del Problema
Los potenciales interatómicos son fundamentales para las simulaciones de dinámica molecular (MD) de materiales magnéticos; sin embargo, la incorporación de características magnéticas en los potenciales para antiferromagnetos (AFM) complejos sigue siendo un desafío significativo. Si bien se han formulado y validado con éxito modelos de espín-red que combinan espín clásico y MD (SD-MD) para ferromagnetos (FM) cúbicos, extenderlos a sistemas AFM es difícil. La dificultad central radica en el requisito de que el modelo reproduzca simultáneamente la dinámica magnética bajo distancias interatómicas variables y las propiedades elásticas gobernadas por el potencial interatómico en presencia de interacciones magnéticas. Actualmente, los potenciales interatómicos disponibles para el material AFM cúbico prototípico, el óxido de níquel (NiO), pueden describir respuestas de carga y de campo eléctrico, pero carecen de interacciones magnéticas. Esta limitación obstaculiza el estudio de fenómenos acoplados de campos eléctricos y magnéticos, como la magnetoelasticidad y los efectos magnetoeléctricos, en los óxidos metálicos.

Metodología
Los autores proponen una metodología para integrar las propiedades magnéticas en los potenciales interatómicos para AFMs cúbicos mediante la adición de un hamiltoniano magnético al potencial interatómico estándar. El hamiltoniano total ( $H_{sl}$ ) incluye la energía cinética, un potencial interatómico no magnético ( $V$ ) y un potencial magnético ( $V_{mag}$ ).

Hamiltoniano Magnético: El término magnético incluye la interacción de intercambio de Heisenberg y la interacción de Néel ( $H_{N\acute{e}el}$ ). La interacción de Néel da cuenta de la anisotropía magnética mediante términos de dipolo ( $l$ ) y cuadrupolo ( $q$ ), que describen la magnetostricción anisotrópica y la anisotropía magnetocristalina (MCA), respectivamente.
Parametrización de Bethe-Slater: La dependencia espacial de los parámetros de intercambio ( $J$ ), dipolo ( $l$ ) y cuadrupolo ( $q$ ) se modela utilizando curvas de Bethe-Slater. Para el caso AFM, los autores derivan expresiones específicas para estos parámetros basadas en el modelo de Néel para una estructura cúbica simple (SC) con espines antiparalelos.
Especificidades del NiO: Dada la complejidad de la celda unitaria del NiO (subred de Ni cúbica centrada en las caras con vecinos más cercanos tanto paralelos como antiparalelos), el modelo simplifica las interacciones centrándose exclusivamente en los segundos vecinos más cercanos (NNN) de los átomos de Ni. Esta elección se justifica porque la interacción de intercambio de los NNN es más fuerte que las interacciones de los vecinos más cercanos en el NiO, y la subred de Ni de los NNN forma una estructura cúbica simple donde todos los espines son antiparalelos, lo que simplifica la derivación de los parámetros.
Construcción del Potencial: Se construyeron dos potenciales interatómicos distintos para validar la metodología:
- SD-MD 1: Basado en el modelo de Born de sólidos iónicos (Fisher y Matsubara), utilizando un potencial de Buckingham para interacciones de corto alcance.
- SD-MD 2: Basado en un método de átomo embebido modificado libre de referencia (RF-MEAM), ajustado a un conjunto de datos de 45 estructuras cristalinas de NiO distorsionadas calculadas mediante la Teoría del Funcional de la Densidad (DFT) con corrección de Hubbard.
Parámetros de Entrada: Los parámetros clave para las curvas de Bethe-Slater (temperatura de Néel, constantes elásticas, magnetostricción de volumen espontánea, constantes magnetoelásticas anisotrópicas y constantes de MCA) se derivaron de una combinación de datos experimentales y cálculos ab initio (usando los códigos VASP y MAELAS).

Contribuciones Clave

Marco Metodológico: El artículo establece un procedimiento general para el desarrollo de potenciales interatómicos para AFMs cúbicos que incluyen tanto el intercambio de Heisenberg como las interacciones de Néel, abordando específicamente el desafío de modelar configuraciones de espín antiparalelos.
Validación de Doble Potencial: Los autores aplicaron con éxito este marco al NiO, creando dos modelos distintos (SD-MD 1 y SD-MD 2) que incorporan interacciones magnetoelásticas.
Derivación de Parámetros: El trabajo proporciona una derivación detallada de los parámetros de Bethe-Slater para sistemas AFM, vinculándolos con propiedades macroscópicas como la temperatura de Néel, la magnetostricción de volumen espontánea y las constantes magnetoelásticas anisotrópicas.
Desarrollo de Potencial Personalizado: El estudio demuestra la creación de un potencial RF-MEAM personalizado para el NiO, el cual ofrece una mejor concordancia con las constantes elásticas calculadas por DFT (específicamente $C_{44}$ ) en comparación con el potencial tradicional del modelo de Born.

Resultados
Los modelos desarrollados fueron validados frente a cálculos de DFT y datos experimentales a temperatura cero:

Propiedades Mecánicas y Magnéticas: Ambos modelos SD-MD mostraron una excelente concordancia con los valores de DFT y experimentales para las propiedades mecánicas (constantes elásticas $C_{11}, C_{12}, C_{44}$ , módulo de compresibilidad) y las propiedades magnéticas (temperatura de Néel, magnetostricción de volumen espontánea $\omega_s$ y la constante de anisotropía magnetocristalina $K_1$ ).
Magnetostricción: Los modelos reprodujeron con precisión la magnetostricción de volumen de entrada ( $\omega_s \approx -0.0014$ ) y las constantes magnetostrictivas anisotrópicas ( $\lambda_{100}, \lambda_{111}$ ).
Dependencia de la Presión: Los modelos capturaron correctamente la dependencia lineal de la constante de anisotropía magnetocristalina $K_1$ respecto a la presión hidrostática, coincidiendo con la pendiente derivada de los cálculos de DFT.
Comparación de Potenciales: El modelo basado en RF-MEAM (SD-MD 2) demostró un rendimiento superior en la reproducción de la constante elástica de cizalla $C_{44}$ y los coeficientes magnetostrictivos resultantes en comparación con el modelo de Born (SD-MD 1), abordando las limitaciones conocidas de este último para el NiO.

Significancia
El artículo presenta un "primer intento" de desarrollar potenciales interatómicos que incluyan propiedades magnéticas específicamente para antiferromagnetos de simetría cúbica. Los autores afirman que este trabajo sirve como prueba de concepto de la metodología a temperatura cero. La principal significancia radica en permitir simulaciones de dinámica molecular a gran escala de fenómenos magnetoelásticos en antiferromagnetos, una capacidad previamente limitada por la falta de potenciales adecuados. Los autores sugieren que estos modelos abren vías para estudiar los efectos acoplados de campos eléctricos y magnéticos en óxidos metálicos, el acoplamiento magnetoelástico y magnetoeléctrico, los efectos de forma en nanoestructuras magnéticas y la dinámica de magnones-fonones. El trabajo es modesto en sus afirmaciones, señalando que los resultados deben interpretarse con cautela respecto a los estudios de temperatura finita, los cuales requieren una investigación más detallada.

Más como este