Article · IA & agents

Accélérer les Hessiennes des potentiels interatomiques sans tout remplir

La différentiation automatique éparse tronquée exploite le graphe de messages d’un potentiel interatomique pour réduire le coût des dérivées secondes.

Un graphe atomique épars guide une différentiation tronquée des Hessiennes d’un potentiel interatomique.

Le coût caché des dérivées secondes

Marcel F. Langer, Adrian Hill et Michele Ceriotti étudient les Hessiennes de potentiels interatomiques fondés sur des réseaux de messages. Une Hessienne complète demande beaucoup plus de produits Hessienne-vecteur que le calcul d’une force ; la structure de voisinage du graphe peut pourtant révéler des blocs nuls ou négligeables.

Colorier le graphe des atomes

La méthode reconstruit la structure éparse à partir de puissances de la matrice d’adjacence. Une version tronquée conserve seulement les couplages jusqu’à un nombre donné de sauts. Elle gagne des couleurs pour les produits Hessienne-vecteur, au prix d’une approximation des termes lointains.

Des gains dépendants de l’observable

Sur MACE et deux modèles Pet évalués sur des matériaux poreux, les couplages décroissent rapidement avec la distance de graphe. La différentiation éparse exacte reste proche de la dense mais n’apporte que des accélérations médianes d’environ 1,7× pour MACE, 1,4× pour Pet-XS et 1,2× pour Pet-S. Avec un seuil de 0,1 % sur la capacité thermique, les versions tronquées atteignent des accélérations médianes annoncées de 11×, 1,6× et 13× selon le modèle.

Interprétation éditoriale : tronquer avec un contrat d’erreur

Pour un agent scientifique, « plus rapide » doit être associé à une observable et à un seuil. Une troncature acceptable pour une capacité thermique peut être insuffisante pour une grandeur qui amplifie les couplages lointains. Le nombre de sauts est donc un hyperparamètre de validation, pas une constante universelle.

Limites et statut de la source

Les tests portent sur des potentiels, structures et configurations matérielles spécifiques, avec calculs en simple précision sur un GPU H100. Les zéolithes convergent plus lentement et les Hessiennes exactes restent peu éparses lorsque la portée d’interaction est grande. Les auteurs laissent ouverte l’analyse théorique de l’erreur et préviennent que les dérivées d’ordre supérieur ne se comporteront pas forcément aussi bien.

Source

Marcel F. Langer, Adrian Hill et Michele Ceriotti, « Truncated automatic sparse differentiation for machine learning interatomic potentials », ArXiv 2609.20510v1 : résumé et preprint, PDF.