The thermodynamics of CaSiO3 in Earth's lower mantle

Planteamiento del Problema
El manto inferior de la Tierra, caracterizado por presiones de 24–127 GPa y temperaturas de 1900–2600 K, permanece inaccesible para la observación directa. El silicato de calcio (CaSiO₃) constituye aproximadamente el 10% de la masa del manto inferior, sin embargo, su caracterización mecánica y termodinámica precisa bajo estas condiciones extremas está incompleta. Aunque los experimentos y la teoría confirman la existencia de una estructura de Ca-perovskita por encima de los 12 GPa, el diagrama de fases sigue siendo controvertido. Específicamente, existen datos conflictivos con respecto a la simetría cristalina a altas presiones y temperaturas: si la fase estable es cúbica ($Pm\bar{3}m$) o tetragonal ($I4/mcm$). Además, las velocidades de propagación de las ondas sísmicas (longitudinales y transversales) y la conductividad térmica de red ($\kappa$) están insuficientemente delimitadas. Los enfoques teóricos estándar, como la aproximación cuasi-armónica (QHA), fallan en el CaSiO₃ cúbico porque los fonones armónicos son inestables en esta fase. Asimismo, el papel del transporte coherente de tipo ondulatorio en la conductividad térmica a altas temperaturas permanece sin investigar en este sistema.

Metodología
Para abordar estos desafíos, los autores emplean simulaciones de primeros principios combinando la Teoría del Funcional de la Densidad (DFT) al nivel PBEsol con la Aproximación Armónica Auto-Consistente Estocástica (SSCHA). El método SSCHA se utiliza para capturar plenamente las fluctuaciones iónicas anharmónicas y el movimiento cero cuántico, que son críticos en las condiciones del manto inferior.

• Estabilidad de Fase: La estabilidad de la fase cúbica de 5 átomos y la fase tetragonal se evalúa mediante el cálculo de los autovalores del Hessiano de la energía libre anharmónica y las diferencias de energía libre. El método minimiza la energía libre cuántica utilizando un ansatz de matriz de densidad gaussiana para los grados de libertad iónicos.
• Propiedades Elásticas: Las velocidades del sonido se derivan de la pendiente de los fonones físicos acústicos a temperaturas finitas, obtenidos a partir del Hessiano de la energía libre. El estudio distingue entre velocidades de sonido isotérmicas e isentrópicas utilizando relaciones termodinámicas estándar.
• Transporte Térmico: La conductividad térmica de la red se evalúa utilizando el recientemente introducido formalismo de Wigner. Este enfoque tiene en cuenta tanto el transporte de tipo partícula (Peierls-Boltzmann) como el transporte ondulatorio por tunelamiento entre bandas de fonones no degeneradas, lo que permite la evaluación de los efectos de transporte coherente en un sistema con fuerte anharmonía.

Contribuciones Clave y Resultados

1. Diagrama de Fases y Mecanismo de Transición:

   • El estudio establece que la fase cúbica ($Pm\bar{3}m$) es el estado estable del CaSiO₃ bajo las condiciones termodinámicas del manto inferior.
   • Se determina que la frontera de fase entre las fases cúbica y tetragonal es de tipo primer orden, aumentando linealmente de 300 K a 1000 K entre 12 GPa y 100 GPa.
   • La temperatura de transición se identifica en 800 K (a 100 GPa). Los autores descartan una transición displaciva de segundo orden, señalando que la fase tetragonal tiene una energía libre menor en el punto donde la fase cúbica se estabiliza.
   • El mecanismo implica el endurecimiento de los modos de rotación de los octaedros de oxígeno (modos inestables M-R) con la temperatura, mientras que las bandas de fonones de alta energía se suavizan. La transición de fase es impulsada por la diferencia de energía libre más que por el ablandamiento de un solo modo de fonón hacia una energía cero.

2. Propiedades Elásticas y Velocidad del Sonido:

   • Las simulaciones calculan los módulos de compresibilidad y de cizalla, así como las velocidades de sonido longitudinal y transversal en función de la profundidad.
   • Los resultados muestran un buen acuerdo con los datos experimentales para la velocidad de sonido longitudinal.
   • Respecto a la velocidad de sonido transversal, el estudio resta importancia al papel de las rotaciones octaédricas, una hipótesis defendida previamente para explicar las discrepancias entre la teoría y el experimento. Los autores demuestran que la deformación de cizalla no estabiliza las inclinaciones octaédricas a altas temperaturas; en su lugar, la cizalla aumenta la barrera de energía libre para la rotación.
   • Aunque la velocidad transversal simulada concuerda con las predicciones de la Dinámica Molecular de Primeros Principios (FPMD), sigue estando ligeramente sobreestimada en comparación con los experimentos, lo que los autores atribuyen a las limitaciones del funcional de intercambio-correlación (PBEsol/GGA).

3. Conductividad Térmica:

   • Utilizando el formalismo de Wigner, el estudio evalúa la conductividad térmica de la red del CaSiO₃.
   • Los resultados indican una predominancia del transporte de tipo partícula bajo las condiciones del manto inferior (alta presión).
   • El transporte ondulatorio (coherente), que se vuelve relevante a presión ambiente y altas temperaturas, se ve suprimido por las altas presiones presentes en el manto inferior.
   • En consecuencia, se justifica el uso de la ecuación de transporte de Boltzmann estándar para este sistema, a pesar de la fuerte anharmonía iónica.

Significancia
El artículo proporciona una caracterización termodinámica rigurosa del CaSiO₃, un constituyente mayoritario del manto inferior de la Tierra, superando las limitaciones de las aproximaciones armónicas mediante el método SSCHA. Al resolver el diagrama de fases como de primer orden e identificar la fase cúbica como estable, el trabajo esclarece la base estructural para el modelado sísmico en la Tierra profunda. Los hallazgos sobre las velocidades del sonido ofrecen una base teórica refinada para la interpretación de datos sísmicos, mientras que el análisis de la conductividad térmica valida el uso de las teorías de transporte estándar en el régimen de alta presión. El estudio subraya la necesidad de tener en cuenta las fluctuaciones cuánticas y anharmónicas para modelar con precisión la dinámica del interior de la Tierra, particularmente en la frontera núcleo-manto, donde la regulación del flujo de calor es crítica.