Publié le 16 février 2026 21:07:00. Des chercheurs ont franchi une étape importante dans la modélisation des matériaux en intégrant avec succès les effets relativistes dans une méthode de calcul quantique de pointe, ouvrant la voie à la conception de nouveaux matériaux aux propriétés inédites.
- Une nouvelle méthode permet de simuler avec précision le comportement des matériaux contenant des éléments lourds, en tenant compte des interactions complexes entre le spin des électrons et leur mouvement orbital (couplage spin-orbite).
- Cette avancée repose sur l’utilisation de pseudopotentiels entièrement relativistes dérivés de l’équation de Dirac, combinés à une technique de Monte Carlo quantique à champ auxiliaire sans phase (pw-AFQMC).
- Les calculs réalisés valident la méthode, notamment en prédisant avec précision l’énergie de dissociation de l’iode, l’énergie de cohésion du plomb et la transition de phase du phosphure d’indium.
Une difficulté de longue date en science des matériaux résidait dans l’intégration des effets relativistes, et plus particulièrement du couplage spin-orbite, dans les simulations numériques. Ces effets, bien que cruciaux pour comprendre les propriétés de nombreux matériaux, notamment ceux contenant des éléments lourds, complexifient considérablement les calculs et nécessitent des ressources informatiques importantes.
L’équipe de recherche, issue d’une collaboration entre le Centre d’études quantiques avancées et l’École de physique et d’astronomie de l’Université normale de Pékin, ainsi que le Centre de physique quantique computationnelle du Flatiron Institute, a développé une approche innovante pour surmonter ces obstacles. Zheng Liu, Shiwei Zhang et Fengjie Ma ont réussi à intégrer le couplage spin-orbite dans la méthode pw-AFQMC en utilisant des pseudopotentiels entièrement relativistes. Ces pseudopotentiels, basés sur l’équation de Dirac – un pilier de la mécanique quantique relativiste – permettent de représenter de manière plus efficace les interactions entre les électrons et les noyaux atomiques.
La méthode pw-AFQMC modifiée utilise un hamiltonien à deux composants dans une base de spineurs, ce qui double la taille de l’espace de calcul pour tenir compte du spin des électrons. Pour valider leur approche, les chercheurs ont calculé l’énergie de dissociation de la molécule d’iode (I2), obtenant une valeur de 7,637 eV, une amélioration significative par rapport aux calculs précédents utilisant des pseudopotentiels scalaires-relativistes (7,513 eV). Ils ont également calculé l’énergie de cohésion du plomb (Pb), obtenant une valeur de 2,341 eV, en accord étroit avec les données expérimentales (2,365 eV) et les calculs de la théorie fonctionnelle de la densité (DFT) (2,378 eV).
L’équipe a ensuite appliqué sa méthode à la prédiction de la pression de transition du phosphure d’indium (InP), un semi-conducteur, lors de son passage d’une structure cristalline de type zinc-blende à une structure cristalline de type sel gemme. En construisant et en analysant des équations d’état, ils ont déterminé avec précision la pression à laquelle cette transition de phase se produit, estimée à 14,2 GPa.
Cette avancée ouvre de nouvelles perspectives pour la découverte et la conception de matériaux innovants dans des domaines tels que la spintronique, l’informatique quantique topologique et le stockage magnétique avancé. Une modélisation précise du couplage spin-orbite est en effet essentielle pour le développement de ces technologies de nouvelle génération. La nouvelle méthode permet de combler le fossé entre la physique quantique fondamentale et les applications pratiques en science des matériaux.
L’utilisation de pseudopotentiels optimisés, conservant les normes sur plusieurs projecteurs, est une étape cruciale pour l’intégration du couplage spin-orbite dans la méthode pw-AFQMC. Ces pseudopotentiels remplacent efficacement les interactions complexes entre les électrons de valence et les noyaux atomiques par un potentiel plus simple à traiter par ordinateur. L’approximation sans phase, une innovation méthodologique clé, permet d’éviter le « problème de signe » qui affecte de nombreuses simulations quantiques de Monte Carlo, en limitant la fonction d’onde d’essai à des valeurs non négatives, améliorant ainsi l’efficacité et la stabilité des calculs.