The thermodynamics of CaSiO3 in Earth's lower mantle

Problematica
Il mantello inferiore della Terra, caratterizzato da pressioni di 24–127 GPa e temperature di 1900–2600 K, rimane inaccessibile all'osservazione diretta. Il silicato di calcio (CaSiO₃) costituisce circa il 10% della massa del mantello inferiore, eppure la sua precisa caratterizzazione meccanica e termodinamica in queste condizioni estreme è incompleta. Sebbene gli esperimenti e la teoria confermino l'esistenza di una struttura a perovskite di Ca sopra i 12 GPa, il diagramma di fase rimane controverso. Nello specifico, esistono dati contrastanti riguardo alla simmetria cristallina ad alte pressioni e temperature: se la fase stabile sia cubica ($Pm\bar{3}m$) o tetragonale ($I4/mcm$). Inoltre, le velocità di propagazione delle onde sismiche (longitudinali e trasversali) e la conducibilità termica reticolare ($\kappa$) sono scarsamente vincolate. Gli approcci teorici standard, come l'approssimazione quasi-armonica (QHA), falliscono nel CaSiO₃ cubico perché i fononi armonici sono instabili in questa fase. Inoltre, il ruolo del trasporto coerente di tipo ondoso nella conducibilità termica ad alte temperature rimane non investigato in questo sistema.

Metodologia
Per affrontare queste sfide, gli autori impiegano simulazioni basate sui primi principi combinando la Teoria del Funzionale della Densità (DFT) al livello PBEsol con l'Approssimazione Armonica Auto-Coerente Stocastica (SSCHA). Il metodo SSCHA è utilizzato per catturare completamente le fluttuazioni ioniche anarmoniche e il moto zero quantistico, che sono critici nelle condizioni del mantello inferiore.

• Stabilità di Fase: La stabilità della fase cubica a 5 atomi e della fase tetragonale viene valutata calcolando gli autovalori dell'essenzialità (Hessian) dell'energia libera anarmonica e le differenze di energia libera. Il metodo minimizza l'energia libera quantistica utilizzando un ansatz di matrice di densità gaussiana per i gradi di libertà ionici.
• Proprietà Elastiche: Le velocità del suono sono derivate dalla pendenza dei fononi fisici acustici a temperature finite, ottenuti dall'essenzialità dell'energia libera. Lo studio distingue tra velocità del suono isoterme e isentropiche utilizzando le relazioni termodinamiche standard.
• Trasporto Termico: La conducibilità termica reticolare è valutata utilizzando il recentemente introdotto formalismo di Wigner. Questo approccio tiene conto sia del trasporto di tipo particella (Peierls-Boltzmann) sia del trasporto di tipo onda (tunneling) tra bande fononiche non degenerate, permettendo la valutazione degli effetti di trasporto coerente in un sistema con forte anarmonicità.

Contributi Chiave e Risultati

1. Diagramma di Fase e Meccanismo di Transizione:

   • Lo studio stabilisce che la fase cubica ($Pm\bar{3}m$) è lo stato stabile di CaSiO₃ nelle condizioni termodinamiche del mantello inferiore.
   • Il confine di fase tra la fase cubica e quella tetragonale è determinato essere di tipo primo ordine, aumentando linearmente da 300 K a 1000 K tra 12 GPa e 100 GPa.
   • La temperatura di transizione è identificata a 800 K (a 100 GPa). Gli autori escludono una transizione displaciva del secondo ordine, notando che la fase tetragonale ha un'energia libera inferiore nel punto in cui la fase cubica si stabilizza.
   • Il meccanismo prevede l'irrigidimento dei modi di rotazione degli ottaedri di ossigeno (modi instabili M-R) con la temperatura, mentre le bande fononiche ad alta energia si ammorsevoliscono. La transizione di fase è guidata dalla differenza di energia libera piuttosto che da un ammorbidimento di un singolo modo fononico verso l'energia zero.

2. Proprietà Elastiche e Velocità del Suono:

   • Le simulazioni calcolano i moduli di bulk e di taglio, nonché le velocità del suono longitudinale e trasversale in funzione della profondità.
   • I risultati mostrano un buon accordo con i dati sperimentali per la velocità del suono longitudinale.
   • Per quanto riguarda la velocità del suono trasversale, lo studio ridimensiona il ruolo delle rotazioni ottaedriche, un'ipotesi precedentemente sostenuta per spiegare le discrepanze tra teoria ed esperimento. Gli autori dimostrano che la deformazione di taglio non stabilizza le inclinazioni ottaedriche ad alte temperature; al contrario, il taglio aumenta la barriera di energia libera per la rotazione.
   • Sebbene la velocità trasversale simulata sia in accordo con le previsioni della Dinamica Molecolare basata sui Primi Principi (FPMD), essa rimane leggermente sovrastimata rispetto agli esperimenti, il che gli autori attribuiscono ai limiti del funzionale di scambio-correlazione (PBEsol/GGA).

3. Conducibilità Termica:

   • Utilizzando il formalismo di Wigner, lo studio valuta la conducibilità termica reticolare di CaSiO₃.
   • I risultati indicano una predominanza del trasporto di tipo particella nelle condizioni del mantello inferiore (alta pressione).
   • Il trasporto di tipo onda (coerente), che diventa rilevante alla pressione ambiente e ad alte temperature, è soppresso dalle alte pressioni presenti nel mantello inferiore.
   • Di conseguenza, l'equazione di trasporto di Boltzmann standard (BTE) è giustificata per l'uso in questo sistema, nonostante la forte anarmonicità ionica.

Significato
Il lavoro fornisce una rigorosa caratterizzazione termodinamica di CaSiO₃, un componente principale del mantello inferiore della Terra, superando i limiti delle approssimazioni armoniche attraverso il metodo SSCHA. Risolvendo il diagramma di fase come primo ordine e identificando la fase cubica come stabile, il lavoro chiarisce la base strutturale per la modellazione sismica della Terra profonda. I risultati riguardanti le velocità del suono offrono una base teorica raffinata per l'interpretazione dei dati sismici, mentre l'analisi della conducibilità termica valida l'uso delle teorie di trasporto standard nel regime di alta pressione. Lo studio sottolinea la necessità di tenere conto delle fluttuazioni quantistiche e anarmoniche per modellare accuratamente la dinamica dell'interno terrestre, particolarmente al confine mantello-nucleo dove la regolazione del flusso di calore è critica.