The thermodynamics of CaSiO3 in Earth's lower mantle

Definição do Problema
O manto inferior da Terra, caracterizado por pressões de 24–127 GPa e temperaturas de 1900–2600 K, permanece inacessível à observação direta. O silicato de cálcio (CaSiO₃) constitui aproximadamente 10% da massa do manto inferior, contudo, sua caracterização mecânica e termodinâmica precisa sob estas condições extremas está incompleta. Embora experimentos e teoria confirmem a existência de uma estrutura de Ca-perovskita acima de 12 GPa, o diagrama de fases permanece controverso. Especificamente, existem dados conflitantes em relação à simetria cristalina em altas pressões e temperaturas: se a fase estável é cúbica ($Pm\bar{3}m$) ou tetragonal ($I4/mcm$). Além disso, as velocidades de propagação das ondas sísmicas (longitudinais e transversais) e a condutividade térmica de rede ($\kappa$) são mal restringidas. Abordagens teóricas padrão, como a aproximação quase-harmônica (QHA), falham na CaSiO₃ cúbica porque os fônons harmônicos são instáveis nesta fase. Adicionalmente, o papel do transporte coerente do tipo onda na condutividade térmica em altas temperaturas permanece não investigado neste sistema.

Metodologia
Para enfrentar esses desafios, os autores empregam simulações de primeiros princípios combinando a Teoria do Funcional da Densidade (DFT) ao nível PBEsol com a Aproximação Harmônica Auto-Consistente Estocástica (SSCHA). O método SSCHA é utilizado para capturar plenamente as flutuações iônicas anarmônicas e o movimento zero-ponto quântico, que são críticos nas condições do manto inferior.

• Estabilidade de Fase: A estabilidade da fase cúbica de 5 átomos e da fase tetragonal é avaliada através do cálculo dos autovalores do Hessiano de energia livre anarmônico e das diferenças de energia livre. O método minimiza a energia livre quântica usando um ansatz de matriz de densidade Gaussiana para os graus de liberdade iônicos.
• Propriedades Elásticas: As velocidades do som são derivadas da inclinação dos fônons físicos acústicos em temperatura finita, obtidos a partir do Hessiano de energia livre. O estudo distingue entre velocidades de som isotérmicas e isentrópicas utilizando relações termodinâmicas padrão.
• Transporte Térmico: A condutividade térmica de rede é avaliada utilizando o formalismo de Wigner, recentemente introduzido. Esta abordagem considera tanto o transporte do tipo partícula (Peierls-Boltzmann) quanto o transporte do tipo onda por tunelamento entre bandas de fônons não degeneradas, permitindo a avaliação de efeitos de transporte coerente em um sistema com forte anarmonicidade.

Principais Contribuições e Resultados

1. Diagrama de Fases e Mecanismo de Transição:

   • O estudo estabelece que a fase cúbica ($Pm\bar{3}m$) é o estado estável da CaSiO₃ sob as condições termodinâmicas do manto inferior.
   • A fronteira de fase entre as fases cúbica e tetragonal é determinada como sendo do tipo de primeira ordem, aumentando linearmente de 300 K para 1000 K entre 12 GPa e 100 GPa.
   • A temperatura de transição é identificada em 800 K (a 100 GPa). Os autores descartam uma transição displaciva de segunda ordem, observando que a fase tetragonal possui uma energia livre menor no ponto onde a fase cúbica se estabiliza.
   • O mecanismo envolve o endurecimento dos modos de rotação dos octaedros de oxigênio (modos instáveis M-R) com a temperatura, enquanto bandas de fônons de alta energia amolecem. A transição de fase é impulsionada pela diferença de energia livre, e não pelo amolecimento de um único modo de fônon para energia zero.

2. Propriedades Elásticas e Velocidade do Som:

   • As simulações calculam os módulos de compressibilidade e de cisalhamento, bem como as velocidades de som longitudinal e transversal em função da profundidade.
   • Os resultados mostram boa concordância com dados experimentais para a velocidade do som longitudinal.
   • Em relação à velocidade do som transversal, o estudo minimiza o papel das rotações octaédricas, uma hipótese anteriormente defendida para explicar discrepâncias entre teoria e experimento. Os autores demonstram que a deformação de cisalhamento não estabiliza as inclinações octaédricas em altas temperaturas; em vez disso, o cisalhamento aumenta a barreira de energia livre para a rotação.
   • Embora a velocidade transversal simulada concorde com as previsões de Dinâmica Molecular de Primeiros Princípios (FPMD), ela permanece ligeiramente superestimada em comparação com os experimentos, o que os autores atribuem às limitações do funcional de troca-correlação (PBEsol/GGA).

3. Condutividade Térmica:

   • Utilizando o formalismo de Wigner, o estudo avalia a condutividade térmica de rede da CaSiO₃.
   • Os resultados indicam uma predominância de transporte do tipo partícula sob as condições do manto inferior (alta pressão).
   • O transporte do tipo onda (coerente), que se torna relevante em pressão ambiente e altas temperaturas, é suprimido pelas altas pressões encontradas no manto inferior.
   • Consequentemente, a equação de transporte de Boltzmann padrão (BTE) é justificada para uso neste sistema, apesar da forte anarmonicidade iônica.

Significância
O artigo fornece uma caracterização termodinâmica rigorosa da CaSiO₃, um constituinte importante do manto inferior da Terra, superando as limitações das aproximações harmônicas através do método SSCHA. Ao resolver o diagrama de fases como sendo de primeira ordem e identificar a fase cúbica como estável, o trabalho esclarece a base estrutural para a modelagem sísmica do interior profundo da Terra. As descobertas sobre as velocidades do som oferecem uma base teórica refinada para interpretar dados sísmicos, enquanto a análise da condutividade térmica valida o uso de teorias de transporte padrão no regime de alta pressão. O estudo ressalta a necessidade de considerar as flutuações quânticas e anarmônicas para modelar com precisão a dinâmica do interior da Terra, particularmente na fronteira núcleo-manto, onde a regulação do fluxo de calor é crítica.