Comprehensive calculations of Fe-Si system via deep learning force field: the core structure of Mercury

Resumen Técnico: Cálculos Exhaustivos del Sistema Fe-Si mediante Campos de Fuerza de Aprendizaje Profundo

Planteamiento del Problema
La estructura interna y la evolución térmica del núcleo de Mercurio siguen siendo poco comprendidas debido a la falta de restricciones observacionales y a las limitaciones de los experimentos de alta presión tradicionales, que a menudo proporcionan solo cantidades físicas aisladas. Los modelos existentes suelen tratar las propiedades termoelásticas como valores fijos, ignorando sus variaciones con la temperatura, la presión y la composición. Además, los enfoques teóricos enfrentan un compromiso: la dinámica molecular ab initio (AIMD) ofrece una alta precisión pero es computacionalmente prohibitiva para sistemas grandes y escalas de tiempo largas, mientras que la dinámica molecular (MD) clásica carece de la transferencia y precisión necesarias para entornos químicos complejos. Existe una necesidad crítica de un marco unificado que pueda extraer sistemáticamente los patrones de variación de todas las cantidades físicas relevantes (ecuación de estado, propiedades termoelásticas, curvas de fusión y propiedades de transporte) como funciones de las variables de estado para modelar con precisión el núcleo de Mercurio.

Metodología
Los autores proponen un marco híbrido que combina cálculos ab initio con redes neuronales profundas para lograr la precisión de ab initio a un costo computacional reducido.

• Campo de Fuerza de Aprendizaje Profundo (DLFF): Se construyó una superficie de energía potencial profunda utilizando el software DeePMD-kit. El modelo fue entrenado con datos generados a partir de simulaciones AIMD (código VASP) utilizando el funcional PBE, potenciales PAW y un modelo de Hubbard U (U ≈ 4.3 eV) para contabilizar los efectos magnéticos en líquidos de Fe-Si. El conjunto de datos de entrenamiento incluyó 15,000 configuraciones de 47 trayectorias AIMD isotérmicas-isobáricas, con un error cuadrático medio (RMSE) de ~1 meV para la energía y ~50 meV/Å para las fuerzas.
• Simulaciones de Dinámica Molecular: Se realizaron simulaciones de MD a gran escala (hasta 2,000 átomos) utilizando el DLFF entrenado dentro del paquete LAMMPS. Estas simulaciones cubrieron sistemas de Fe-Si con concentraciones de Si que oscilan entre el 5 y el 25 at.% bajo presiones de 6–40 GPa y temperaturas de 1600–3000 K.
• Cálculo de Propiedades:
  • Ecuación de Estado (EOS): Se aplicó un ajuste polinómico multivariante a los datos P-ρ-T y E-ρ-T para derivar las propiedades termoelásticas (módulo de compresibilidad, expansión térmica, parámetro de Grüneisen, velocidad del sonido).
  • Curvas de Fusión: Se emplearon métodos de integración termodinámica (TI) para calcular la diferencia de energía libre de Gibbs entre las fases líquida y sólida (bcc), identificando los puntos de fusión donde $G_{líquido} = G_{sólido}$.
  • Propiedades de Transporte: La conductividad térmica electrónica se calculó mediante la fórmula de Kubo-Greenwood utilizando instantáneas de AIMD. La conductividad térmica de la red, la viscosidad y los coeficientes de difusión se derivaron de simulaciones de MD utilizando fórmulas de Green-Kubo y relaciones de Einstein.
  • Modelado del Núcleo: Se utilizó un modelo termodinámico de promedio esférico 1D para resolver los perfiles de temperatura y presión adiabáticos, incorporando la EOS calculada y las condiciones de contorno en la frontera núcleo-manto (CMB).

Contribuciones Clave y Resultados

• Cobertura Unificada de Propiedades Físicas: El estudio generó con éxito un conjunto de datos exhaustivo para el sistema Fe-Si (<25 at.% Si) que cubre la ecuación de estado, propiedades termoelásticas, curvas adiabáticas, curvas de fusión y propiedades de transporte con precisión de ab initio.
• Validación del DLFF: El campo de fuerza de aprendizaje profundo reprodujo con éxito la precisión de nivel AIMD para la presión y la energía, con discrepancias que caen dentro del error estadístico de AIMD (1%). Esto valida el uso de DLFF para el cálculo de propiedades complejas como la viscosidad y la difusión en sistemas grandes.
• Curvas de Fusión y Adiabáticas: Las curvas de fusión calculadas para Fe-10 at.% Si y Fe-25 at.% Si muestran que una parte significativa del núcleo de Mercurio es probablemente sólida. Las curvas adiabáticas son relativamente planas en comparación con las curvas de fusión.
• Propiedades de Transporte:
  • Viscosidad y Difusión: A bajas temperaturas (1600 K), el sistema exhibe una alta viscosidad (45 mPa s) y bajos coeficientes de difusión (0.71 × 10⁻⁹ m² s⁻¹), lo que indica una fase sólida. A temperaturas más altas (2000 K), la viscosidad cae a ~6.4 mPa s y la difusión aumenta a ~4.0 × 10⁻⁹ m² s⁻¹, lo cual es consistente con un comportamiento líquido.
  • Conductividad Térmica: La conductividad térmica total se estima en 48.0 W/m/K (electrónica: 46.6 W/m/K; de red: ~1.4 W/m/K). Este valor es significativamente mayor que el de FeS (4 W/m/K) pero menor que el del líquido de Fe no magnético en el núcleo de la Tierra.
• Estructura del Núcleo de Mercurio:
  • Núcleo Interno Sólido: La intersección de las curvas adiabáticas y de fusión sugiere la existencia de un núcleo interno sólido, que comprende aproximadamente la mitad a dos tercios del radio del núcleo. Se infiere que la concentración de Si en el núcleo interno sólido disminuye con el radio.
  • Rechazo de la Nieve de Hierro: Las pendientes de las curvas adiabáticas y de fusión calculadas ($dP_{ad}/dT_{ad} > dP_{melt}/dT_{melt}$) no respaldan los mecanismos de cristalización de arriba hacia abajo (nieve de hierro).
  • Estratificación Térmica: La alta conductividad térmica (48.0 W/m/K) resulta en un flujo de calor adiabático (~20 mW m⁻²) que excede las estimaciones para el flujo de calor de la frontera núcleo-manto (10–15 mW m⁻²). Esta condición ($Q_{ad} > Q_{cmb}$) sugiere la posibilidad de una estratificación térmica sobre el núcleo.

Significancia y Reivindicaciones
El artículo afirma proporcionar la "imagen física más posible" de la estructura y la evolución térmica del núcleo de Mercurio basada en las propiedades calculadas. Los autores concluyen que el núcleo de Mercurio probablemente contiene un núcleo interno sólido sobre el cual se asienta una capa térmicamente estratificada, mientras que descartan explícitamente el mecanismo de "nieve de hierro" (cristalización de arriba hacia abajo) bajo las condiciones modeladas actuales.

El estudio enfatiza que estas conclusiones se derivan específicamente para sistemas de Fe-Si con concentraciones de Si por debajo del 25 at.%. Los autores señalan modestamente que sus cálculos de la curva de fusión se centran en la transición líquido-bcc y no contemplan la cinética de transiciones de fase sólida complejas u otros elementos potenciales (p. ej., S) que podrían conducir a la inmiscibilidad o a diferentes conductividades térmicas. El trabajo sirve como un conjunto de datos fundacional para futuros modelos de balance de energía y entropía de la evolución térmica de Mercurio.