Pressure-Driven Evolution of the Physical Properties of Layered Y2GaB and Y2InB MAB Phases: A Comparative First-Principles Study

Resumo Técnico: Evolução da Propriedade Física dos Fases MAB de Camadas Y2GaB e Y2InB sob Pressão

Definição do Problema
As fases MAB, boretos ternários em camadas, exibem uma combinação única de comportamento eletrônico metálico e rigidez cerâmica devido à sua estrutura nanolaminada, onde ligações fortes metal-boro no plano coexistem com interações intercamadas mais fracas. Embora estudos anteriores tenham caracterizado compostos MAB específicos sob pressão ambiente, há uma falta de compreensão unificada sobre como a pressão hidrostática e a substituição no sítio A (especificamente Ga vs. In) influenciam sua estabilidade estrutural, anisotropia mecânica, topologia eletrônica e propriedades termofísicas. Além disso, trabalhos computacionais anteriores sobre Y2GaB e Y2InB relataram modos de fônons de frequência negativa, levantando questões sobre sua estabilidade dinâmica e a sensibilidade dessas avaliações aos parâmetros computacionais (tratamentos de célula unitária vs. supercélula). Este estudo aborda essas lacunas comparando sistematicamente a evolução dependente da pressão de Y2GaB e Y2InB de 0 a 20 GPa para estabelecer janelas de pressão mecanicamente viáveis e esclarecer as relações estrutura-propriedade.

Metodologia
O estudo emprega cálculos de primeiros princípios baseados na teoria do funcional da densidade (DFT) de Kohn-Sham usando o código CASTEP com a aproximação de gradiente generalizado (GGA) de Perdew-Burke-Ernzerhof (PBE). As principais características metodológicas incluem:

• Otimização Estrutural: Estruturas hexagonais $P6_3/mmc$ foram otimizadas sob pressões hidrostáticas (0, 5, 10, 15 e 20 GPa) usando o algoritmo BFGS.
• Estabilidade Elástica: Constantes elásticas foram derivadas via o método de tensão-deformação finita. A estabilidade foi avaliada usando critérios corrigidos por pressão para simetria hexagonal, distinguindo entre estabilidade de deformação homogênea e estabilidade de fônons.
• Derivação de Propriedades: Módulos policristalinos (médias de Voigt-Reuss-Hill), velocidades acústicas, temperaturas de Debye e estimativas de condutividade térmica (modelos de Clarke e Slack) foram calculados para estados elasticamente estáveis.
• Análise Eletrônica e Óptica: Estruturas de bandas, densidades de estados (DOS) e superfícies de Fermi foram computadas com e sem acoplamento spin-órbita (SOC) usando o Quantum ESPRESSO. Propriedades ópticas foram resolvidas para polarizações de campo elétrico ao longo de [100] e [001].
• Ligação e Dinâmica: Mapas de diferença de densidade eletrônica, análise de carga de Mulliken/Hirshfeld e populações de sobreposição de ligação foram usados para caracterizar a ligação. Dispersões de fônons foram calculadas usando o método de deslocamento finito tanto para configurações de célula unitária quanto para supercélulas expandidas para avaliar a estabilidade dinâmica e a sensibilidade numérica.

Resultos Principais

• Evolução Estrutural: Ambos os compostos mantêm a estrutura $P6_3/mmc$ sob pressão. A compressão é anisotrópica, com o parâmetro de rede basal ($a$) contraindo-se mais significativamente que o parâmetro axial ($c$). Consequentemente, a razão $c/a$ aumenta para ambos, e o volume da célula unitária diminui monotonicamente (aprox. 21,8% para Y2GaB e 20,3% para Y2InB de 0–20 GPa).
• Estabilidade Elástica e Propriedades Mecânicas:
  • Janelas de Estabilidade: O Y2GaB satisfaz os critérios de estabilidade elástica corrigidos por pressão a 0, 5 e 10 GPa, mas torna-se instável a 15 e 20 GPa. O Y2InB é elasticamente estável apenas à pressão ambiente (0 GPa).
  • Resposta Mecânica: Dentro de seu intervalo estável, o Y2GaB exibe aumento nos módulos de bulk, cisalhamento e Young com a pressão. Ambos os compostos são previstos como dúcteis (razão de Pugh $G/B < 0,57$). O Y2GaB mostra alta anisotropia elástica a 0 GPa ($A^U = 5,25$), que diminui com a pressão, enquanto o Y2InB é menos anisotrópico à pressão ambiente.
  • Propriedades Termofísicas: Para o Y2GaB estável, a compressão aumenta as velocidades do som, a temperatura de Debye e as temperaturas de fusão estimadas. A condutividade térmica de rede estimada por Slack a 300 K sobe de 0,605 W m⁻¹ K⁻¹ (0 GPa) para 1,543 W m⁻¹ K⁻¹ (10 GPa). Ambos os compostos exibem condutividades térmicas mínimas muito baixas, sugerindo potencial para aplicações de barreira térmica em baixas temperaturas, limitadas por seus pontos de fusão relativamente baixos.
• Estrutura Eletrônica: Ambos os materiais são metais multibanda. O Y2GaB possui nove folhas de superfície de Fermi, e o Y2InB possui oito à pressão ambiente. O acoplamento spin-órbita (SOC) modifica as dispersões de banda locais, mas não abre um gap nem altera o número de cruzamentos da superfície de Fermi. Nenhuma transição de Lifshitz (reconstrução topológica da superfície de Fermi) é observada entre 0 e 10 GPa para o Y2GaB, embora o bolso central encolha significativamente. Os estados eletrônicos próximos ao nível de Fermi são dominados por orbitais Y-4d, com contribuições secundárias de orbitais B-2p e orbitais p do sítio A.
• Caráter de Ligação: A análise revela ligação mista iônico-covalente entre Y e B, caracterizada pelo acúmulo de carga ao redor do B e depleção ao redor do Y. Os elementos do sítio A (Ga/In) exibem ligação metálica mais difusa.
• Propriedades Ópticas: Ambos os compostos exibem respostas Drude metálicas e forte anisotropia óptica entre as polarizações [100] e [001]. Eles atuam como refletores eficientes na região visível e absorvedores no ultravioleta. A pressão induz deslocamentos não uniformes nas características ópticas, geralmente elevando as energias características para o Y2GaB.
• Dinâmica de Fônons: Cálculos de fônons de célula unitária mostram frequências positivas para ambos os compostos a 0 e 10 GPa, indicando estabilidade dinâmica dentro deste tratamento restrito. No entanto, cálculos de supercélula expandida revelam modos imaginários, destacando a sensibilidade da estabilidade de fônons de baixa frequência ao tamanho da célula e ao tratamento da constante de força. Notavelmente, o Y2InB permanece elasticamente instável a 10 GPa, apesar de mostrar frequências positivas no espectro de fônons de célula unitária, ressaltando a distinção entre critérios de estabilidade elástica e dinâmica.

Significância e Alegações
O artigo estabelece um quadro comparativo para entender como a substituição do sítio A (Ga vs. In) e a pressão hidrostática modulam as propriedades das fases MAB. Suas principais contribuições são:

1. Definição de Janelas de Estabilidade: Identifica intervalos de pressão específicos onde o Y2GaB é mecanicamente viável (0–10 GPa) e esclarece que o Y2InB é estável apenas à pressão ambiente sob os critérios adotados.
2. Desacoplamento de Tipos de Estabilidade: O estudo demonstra explicitamente que a estabilidade elástica (resistência à deformação homogênea) e a estabilidade dinâmica (resistência a deslocamentos atômicos) são distintas e devem ser avaliadas separadamente, notando particularmente a sensibilidade dos cálculos de fônons ao tamanho da supercélula.
3. Relações Estrutura-Propriedade: Detalha como a pressão aumenta sistematicamente a rigidez, o transporte térmico e a resposta óptica no Y2GaB, ao mesmo tempo que reduz a anisotropia elástica, fornecendo uma base para futuras validações experimentais.
4. Escopo Modesto: Os autores afirmam explicitamente que não existem derivados MABene experimentalmente estabelecidos desses compostos; o estudo foca nas fases parentais como uma base necessária para qualquer avaliação futura de derivados de dimensão reduzida. Os valores de condutividade térmica são apresentados como estimativas de modelo (Slack/Clarke) em vez de previsões precisas de transporte, reconhecendo a ausência de cálculos explícitos de espalhamento elétron-fônon e anarmônico.