Depuis Newton, la méthode scientifique repose sur la recherche d'une équation unique pour un phénomène donné. Même la régression symbolique moderne, qui déduit automatiquement des formules à partir des données, hérite silencieusement de cette hypothèse : elle cherche une seule courbe optimale. Mais les systèmes réels, du caoutchouc aux galaxies, basculent entre plusieurs mécanismes lorsque les conditions changent. À côté de ce travail, des études dans les archives montrent des réseaux de neurones apprenant à déduire des lois de la nature ou à construire des cartes gravitationnelles, mais toutes supposent aussi une loi unique. Cet article fait un pas de côté : il ne demande pas « quelle équation est la meilleure », mais « combien existent simultanément ».
Les auteurs ont développé un pipeline en quatre étapes. D'abord, trois recherches parallèles de régression symbolique (globale, par conditions et par blocs spatiaux) rassemblent un ensemble d'équations de topologies différentes. Ensuite, pour chaque cellule de l'espace des paramètres, on calcule le poids de domination de chaque équation via la probabilité a posteriori bayésienne. À la troisième étape, une régression symbolique secondaire trouve les lois analytiques d'évolution de ces poids, avec éventuellement une contrainte de monotonie pour les processus irréversibles. Enfin, les poids sont extrapolés à de nouvelles conditions, et la prédiction est construite comme une somme pondérée de toutes les équations. Le nombre de mécanismes K* est déterminé automatiquement par le critère d'information BIC.
Dans les quatre systèmes testés (élastomères, ébullition, galaxies, évaporation de gouttes), le BIC a choisi exactement trois équations indépendantes. Pour les élastomères, il s'agit de la loi de puissance du durcissement par cristallisation, de la loi en racine de l'élasticité entropique et de la loi linéaire de Hooke ; la prédiction s'est améliorée de 72 % par rapport au gradient boosting. Dans l'ébullition, la contrainte de monotonie a apporté 20 % de précision supplémentaire. Sur 3269 mesures de 163 galaxies, le système a lui-même trouvé les lois newtonienne, de transition et de « champ constant », et la différence de l'accélération caractéristique a0 d'un facteur 2,2 entre les galaxies pauvres et riches en gaz exclut une courbe MOND unique avec p < 10^-34. Dans l'évaporation des gouttes, le moment du détachement de la ligne de contact a été automatiquement détecté.
Ce travail se distingue des études voisines en ce qu'il ne cherche pas une seule formule optimale, mais cartographie tout un écosystème de mécanismes. Les travaux voisins utilisent des réseaux de neurones pour déduire des lois ou construisent des cartes gravitationnelles, mais ils restent tous dans le paradigme K=1. Ici, K* devient la première grandeur mesurable, ce qui ouvre une direction orthogonale : non pas améliorer les équations, mais décrire leur coexistence. Cela transforme le débat sur l'universalité d'une loi d'une question philosophique en un test computationnel.
Le sujet pourrait évoluer vers une théorie dynamique de l'écologie des lois : non pas une carte statique, mais une description de la symbiose, de la compétition et de la succession des équations dans le temps. Il est possible d'étendre l'ensemble des opérateurs pour couvrir les fonctions trigonométriques et spéciales. L'application à la turbulence, à la fatigue des matériaux et à l'évolution stellaire permettrait de tester l'universalité de l'approche. La création d'un benchmark interdisciplinaire pour diagnostiquer le nombre de mécanismes dans différents domaines est également prometteuse.
La méthode s'applique à la mécanique des matériaux, à la thermophysique, à l'astrophysique et à la microfluidique. Elle pourrait changer la façon de tester les lois universelles en biologie et en climatologie.
La prochaine étape consiste à appliquer le cadre à la loi d'échelle métabolique et à la sensibilité climatique, où l'on débat de l'existence d'une ou de plusieurs lois. Il faut aussi élargir le vocabulaire des opérateurs et améliorer l'algorithme de sélection de K.
Dans le domaine du calcul scientifique, la question de savoir comment déterminer automatiquement le nombre de régimes dans les données sans hypothèses a priori reste ouverte. Ce travail propose un critère BIC concret et un champ continu de poids au lieu de frontières rigides, ce qui résout un problème ancien de transitions douces entre mécanismes.
Une bougie, que nous voyons comme une flamme unique, est en réalité régie par quatre équations différentes : la conduction thermique dans la cire, la combustion par diffusion dans la zone bleue, le rayonnement dans la partie brillante et la convection dans la langue au-dessus de la flamme.