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

Résumé Technique : Évolution sous pression des propriétés physiques des phases MAB de type Y2GaB et Y2InB

Énoncé du Problème
Les borures ternaires stratifiés, connus sous le nom de phases MAB, présentent une combinaison unique de comportement électronique métallique et de rigidité de type céramique en raison de leur structure nanolaminée, où des liaisons métal-bore fortes dans le plan coexistent avec des interactions interfoliaires plus faibles. Bien que des études antérieures aient caractérisé des composés MAB spécifiques à la pression ambiante, il manque une compréhension unifiée de la manière dont la pression hydrostatique et la substitution sur le site A (spécifiquement Ga vs In) influencent leur stabilité structurelle, leur anisotropie mécanique, leur topologie électronique et leurs propriétés thermophysiques. De plus, des travaux de calcul antérieurs sur Y2GaB et Y2InB ont rapporté des modes de phonons à fréquence négative, soulevant des questions quant à leur stabilité dynamique et à la sensibilité de ces évaluations aux paramètres de calcul (traitements par cellule élémentaire vs supercellule). Cette étude aborde ces lacunes en comparant systématiquement l'évolution dépendante de la pression de Y2GaB et Y2InB de 0 à 20 GPa afin d'établir des fenêtres de pression mécaniquement viables et de clarifier les relations structure-propriétés.

Méthodologie
L'étude emploie des calculs de premier principe basés sur la théorie de la fonctionnelle de la densité (DFT) de Kohn-Sham à l'aide du code CASTEP utilisant l'approximation du gradient généralisé (GGA) de Perdew-Burke-Ernzerhof (PBE). Les caractéristiques méthodologiques clés incluent :

• Optimisation Structurelle : Les structures hexagonales $P6_3/mmc$ ont été optimisées sous pressions hydrostatiques (0, 5, 10, 15 et 20 GPa) en utilisant l'algorithme BFGS.
• Stabilité Élastique : Les constantes élastiques ont été dérivées via la méthode de contrainte-déformation finie. La stabilité a été évaluée à l'aide de critères corrigés par la pression pour la symétrie hexagonale, distinguant la stabilité de déformation homogène de la stabilité de phonon.
• Dérivation des Propriétés : Les modules polycristallins (moyennes de Voigt-Reuss-Hill), les vitesses acoustiques, les températures de Debye et les estimations de la conductivité thermique (modèles de Clarke et de Slack) ont été calculés pour les états élastiquement stables.
• Analyse Électronique et Optique : Les structures de bandes, les densités d'états (DOS) et les surfaces de Fermi ont été calculées avec et sans couplage spin-orbite (SOC) via Quantum ESPRESSO. Les propriétés optiques ont été résolues pour les polarisations de champ électrique selon [100] et [001].
• Liaison et Dynamique : Des cartes de différence de densité électronique, des analyses de charges de Mulliken/Hirshfeld et des populations de recouvrement de liaison ont été utilisées pour caractériser les liaisons. Les dispersions de phonons ont été calculées en utilisant la méthode des déplacements finis pour des configurations de cellule élémentaire et de supercellule étendue afin d'évaluer la stabilité dynamique et la sensibilité numérique.

Résultats Clés

• Évolution Structurelle : Les deux composés conservent la structure $P6_3/mmc$ sous pression. La compression est anisotrope, le paramètre de réseau basal ($a$) se contractant plus significativement que le paramètre axial ($c$). Par conséquent, le rapport $c/a$ augmente pour les deux, et le volume de la cellule unité diminue de manière monotone (environ 21,8 % pour Y2GaB et 20,3 % pour Y2InB sur 0–20 GPa).
• Stabilité Élastique et Propriétés Mécaniques :
  • Fenêtres de Stabilité : Y2GaB satisfait les critères de stabilité élastique corrigés par la pression à 0, 5 et 10 GPa, mais devient instable à 15 et 20 GPa. Y2InB est élastiquement stable uniquement à la pression ambiante (0 GPa).
  • Réponse Mécanique : Dans sa plage de stabilité, Y2GaB présente une augmentation des modules de compressibilité, de cisaillement et de Young avec la pression. Les deux composés sont prédits comme ductiles (rapport de Pugh $G/B < 0,57$). Y2GaB présente une forte anisotropie élastique à 0 GPa ($A^U = 5,25$), qui diminue avec la pression, tandis que Y2InB est moins anisotrope à la pression ambiante.
  • Propriétés Thermophysiques : Pour Y2GaB stable, la compression augmente les vitesses du son, la température de Debye et les températures de fusion estimées. La conductivité thermique de réseau estimée par Slack à 300 K passe de 0,605 W m⁻¹ K⁻¹ (0 GPa) à 1,543 W m⁻¹ K⁻¹ (10 GPa). Les deux composés présentent de très faibles conductivités thermiques minimales, suggérant un potentiel pour des applications de barrière thermique à basse température, limité par leurs points de fusion relativement bas.
• Structure Électronique : Les deux matériaux sont des métaux multibandes. Y2GaB possède neuf feuillets de surface de Fermi, et Y2InB en possède huit à la pression ambiante. Le couplage spin-orbite (SOC) modifie les dispersions de bandes locales mais n'ouvre pas de gap ni ne change le nombre de croisements de la surface de Fermi. Aucune transition de Lifshitz (reconstruction topologique de la surface de Fermi) n'est observée entre 0 et 10 GPa pour Y2GaB, bien que la poche centrale rétrécisse significativement. Les états électroniques près du niveau de Fermi sont dominés par les orbitales Y-4d, avec des contributions secondaires des orbitales B-2p et des orbitales p du site A.
• Caractère de Liaison : L'analyse révèle une liaison mixte iono-covalente entre Y et B, caractérisée par une accumulation de charge autour de B et un appauvrissement autour de Y. Les éléments du site A (Ga/In) présentent une liaison métallique plus diffuse.
• Propriétés Optiques : Les deux composés affichent des réponses de Drude métalliques et une forte anisotropie optique entre les polarisations [100] et [001]. Ils agissent comme des réflecteurs efficaces dans le visible et des absorbeurs dans l'ultraviolet. La pression induit des décalages non uniformes des caractéristiques optiques, augmentant généralement les énergies caractéristiques pour Y2GaB.
• Dynamique des Phonons : Les calculs de phonons de la cellule élémentaire montrent des fréquences positives pour les deux composés à 0 et 10 GPa, indiquant une stabilité dynamique au sein de ce traitement restreint. Cependant, les calculs de supercellule étendue révèlent des modes imaginaires, soulignant une sensibilité de la stabilité des phonons de basse fréquence à la taille de la cellule et au traitement des constantes de force. Notamment, Y2InB reste élastiquement instable à 1 o kPa malgré l'affichage de fréquences positives dans le spectre de phonons de la cellule élémentaire, soulignant la distinction entre les critères de stabilité élastique et dynamique.

Signification et Revendications
L'article établit un cadre comparatif pour comprendre comment la substitution du site A (Ga vs In) et la pression hydrostatique modulent les propriétés des phases MAB. Ses principales contributions sont :

1. Définition des Fenêtres de Stabilité : Il identifie les plages de pression spécifiques où Y2GaB est mécaniquement viable (0–10 GPa) et clarifie que Y2InB n'est stable qu'à la pression ambiante selon les critères adoptés.
2. Découplage des Types de Stabilité : L'étude démontre explicitement que la stabilité élastique (résistance à la déformation homogène) et la stabilité dynamique (résistance aux déplacements atomiques) sont distinctes et doivent être évaluées séparément, notant particulièrement la sensibilité des calculs de phonons à la taille de la supercellule.
3. Relations Structure-Propriété : Elle détaille comment la pression augmente systématiquement la rigidité, le transport thermique et la réponse optique pour Y2GaB tout en réduisant l'anisotropie élastique, fournissant une base pour toute validation expérimentale future.
4. Portée Modeste : Les auteurs déclarent explicitement qu'aucun dérivé de type MABene de ces composés n'est établi expérimentalement ; l'étude se concentre sur les phases parentes comme une fondation nécessaire pour toute évaluation future de dérivés de dimension réduite. Les valeurs de conductivité thermique sont présentées comme des estimations de modèles (Slack/Clarke) plutôt que des prédictions de transport précises, reconnaissant l'absence de calculs explicites de diffusion électron-phonon et d'anharmonicité.