Microscopic Origin of the Ultralow Lattice Thermal Conductivity in Vacancy-Ordered Halide Double Perovskites Cs$_2BX_6$ ($B$ = Zr, Pd, Sn, Te, Hf, and Pt; $X$= Cl, Br, and I)

Resumo Técnico: Origem Microscópica da Condutividade Térmica de Rede Ultrabaixa em Perovskitas Duplas de Haleto com Ordem de Vacância

Definição do Problema
Perovskitas duplas de haleto sem chumbo com ordem de vacância, com a fórmula $Cs_2BX_6$ (onde B = Zr, Pd, Sn, Te, Hf, Pt e X = Cl, Br, I), surgiram como candidatas promissoras para aplicações optoeletrônicas e termoelétricas devido à sua intrínseca condutividade térmica de rede ($\kappa_L$) ultrabaixa. Embora estudos anteriores tenham identificado compostos específicos como $Cs_2SnI_6$ e $Cs_2SnBr_6$ como possuidores de um $\kappa_L$ excepcionalmente baixo, carece-se de uma compreensão sistemática de como o $\kappa_L$ evolui em toda a família $Cs_2BX_6$. Especificamente, os mecanismos microscópicos que regem esses valores ultrabaixos — se derivam do forte espalhamento de fônons devido a modos de "chocalho" (rattling) em grandes vazios estruturais ou de uma ligação química intrinsecamente fraca — ainda não foram plenamente esclarecidos. Além disso, as contribuições relativas dos processos de transporte de fônons de três, quatro e coerentes para a resistência térmica total exigem uma quantificação rigorosa.

Metodologia
Os autores empregaram cálculos de primeira ordem de última geração baseados na Teoria do Funcional da Densidade (DFT) utilizando o Vienna Ab initio Simulation Package (VASP). O estudo utilizou o funcional de troca-correlação PBEsol e o método de onda plana aumentada pelo projetor (PAW). As etapas computacionais principais incluíram:

• Constantes de Força: Constantes de força harmônicas foram calculadas usando o método de deslocamento finito (Phonopy) com supercélulas grandes (4×4×4). Constantes de força anarmônicas de ordem superior foram extraídas usando o método de dinâmica de rede por compressão (Compressive Sensing Lattice Dynamics - CSLD).
• Renormalização de Fônons: Para considerar os efeitos de temperatura, as constantes de força de segunda ordem foram renormalizadas usando o método de fônons autoconsistentes (SCPH) (pacote ALAMODE).
• Transporte Térmico: A condutividade térmica de rede foi calculada resolvendo iterativamente a equação de transporte de Peierls–Boltzmann. O estudo incorporou processos de espalhamento de três e quatro fônons (pacote FourPhonon) e avaliou a contribuição coerente ($\kappa^C_L$) usando uma teoria unificada de transporte térmico para cristais e vidros.
• Aprendizado de Máquina: Um modelo preditivo para $\kappa_L$ foi construído usando o algoritmo Sure Independence Screening and Sparsifying Operator (SISSO) para identificar os descritores físicos mais críticos que regem o transporte térmico.
• Análise de Ligação Química: A análise de População de Hamilton de Orbitais Cristalinos (COHP) foi realizada para quantificar as forças de ligação química e sua relação com as constantes de força.

Resultos Principais

1. Condutividade Térmica Ultrabaixa: Todos os compostos $Cs_2BX_6$ investigados exibem valores de $\kappa_L$ intrinsecamente ultrabaixos abaixo de 1,0 W m$^{-1}$ K$^{-1}$ a 300 K. Os valores desviam-se significativamente da dependência de temperatura convencional de $T^{-1}$, mostrando uma dependência de temperatura muito mais fraca ($\kappa_L \propto T^{-\alpha}$ com $\alpha \ll 1$).
2. Mecanismo Dominante: Contrariando a expectativa de que grandes vazios estruturais levam ao forte espalhamento de fônons por "chocalho", a análise revela que o principal fator que suprime o transporte de calor na maioria dos compostos é a baixa velocidade de grupo de fônons ($\nu_g$), variando de 1100 a 1600 m s$^{-1}$. Esta baixa velocidade decorre de ligações químicas intrinsecamente fracas dentro da estrutura de haleto.
3. Contribuições de Espalhamento:
   • O espalhamento de três fônons é o mecanismo de resistência dominante na temperatura ambiente para a maioria dos compostos.
   • O espalhamento de quatro fônons reduz o $\kappa_L$ em aproximadamente 20–25% na maioria dos casos, embora isso varie com a composição (por exemplo, até 35% em $Cs_2ZrCl_6$).
   • O transporte coerente contribui com cerca de 15% para o $\kappa_L$ total na maioria dos sistemas, subindo para 21% em $Cs_2SnI_6$.
4. Anomalias Específicas:
   • $Cs_2SnI_6$: Exibe o $\kappa_L$ mais baixo (0,28 W m$^{-1}$ K$^{-1}$) entre os compostos estudados. Diferente da tendência geral, sua condutividade ultrabaixa é impulsionada por um espalhamento de três fônons excepcionalmente forte na faixa de baixa frequência (0,5–1,0 THz) devido a um grande espaço de fase de espalhamento associado a ramos ópticos planos dominados por vibrações de Iodo.
   • $Cs_2TeX_6$: Mostra um comportamento anômalo onde o $\kappa_L$ não diminui estritamente com o aumento da massa do haleto (I > Br > Cl), atribuído a um espalhamento de três fônons mais forte nos variantes de cloreto e brometo em comparação ao iodeto.
5. Insights de Aprendizado de Máquina: O modelo SISSO, que superou os modelos tradicionais de Slack e o modelo de Slack modificado (alcançando um $R^2$ de 0,96), identificou que o $\kappa_L$ escala linearmente com a velocidade do som ($\nu_g$) em vez de $\nu_g^3$. O modelo destacou que a combinação da velocidade do som, módulo de compressibilidade e o deslocamento médio quadrático do átomo do sítio-B é crucial para a predição, reforçando o papel da ligação química fraca.
6. Tendências de Ligação Química: A análise COHP confirmou que a força de ligação B–X diminui de Cl para I. A configuração eletrônica do cátion do sítio-B influencia significativamente a ligação; por exemplo, compostos baseados em Pd mostram ligações mais fracas (menor COHP integrado) do que os baseados em Zr devido à ocupação de estados antibondantes $t_{2g}$.

Significância e Alegações
O artigo afirma fornecer uma compreensão microscópica unificada do transporte térmico em perovskitas duplas com ordem de vacância. Ele desafia a suposição prevalecente de que os modos de "chocalho" em grandes vazios são a causa primária do $\kappa_L$ ultrabaixo nesta família. Em vez disso, os autores afirmam que a ligação química intrinsecamente fraca, levando a baixas velocidades de som, é o mecanismo dominante para a maioria dos compostos $Cs_2BX_6$.

O estudo esclarece os papos distintos da química do sítio-B e do sítio-X:

• A espécie de haleto (X) desempenha um papel dominante em determinar o $\kappa_L$ nos compostos baseados em cloreto.
• A influência do cátion do sítio-B é mais pronunciada nos cloretos, mas menos significativa nos brometos e iodetos, com a exceção notável dos compostos baseados em Telúrio.

Ao estabelecer que a ligação fraca e as baixas velocidades de som são as causas raízes da condutividade térmica de rede ultrabaixa, o trabalho oferece princípios orientadores para o design racional de materiais à base de haletos com propriedades de transporte térmico ajustadas para aplicações termoelétricas e de gestão térmica. As descobertas também destacam o $Cs_2SnI_6$ como um caso único onde o forte espalhamento anarmônico, e não apenas a baixa velocidade, impulsiona a condutividade térmica recorde.