El método científico, desde Newton, se ha basado en buscar una única ecuación para un único fenómeno. Incluso la regresión simbólica moderna, que deduce fórmulas automáticamente a partir de datos, hereda en silencio ese supuesto: busca una única curva óptima. Pero los sistemas reales, desde la goma hasta las galaxias, cambian de mecanismo cuando varían las condiciones. Junto a este trabajo, en el archivo hay estudios donde redes neuronales aprenden a deducir leyes naturales o a construir mapas gravitatorios, pero todos ellos también asumen una sola ley. Este artículo da un paso al lado: no pregunta «qué ecuación es mejor», sino «cuántas existen a la vez».
Los autores desarrollaron un proceso en cuatro pasos. Primero, tres búsquedas paralelas de regresión simbólica (global, por condiciones y por bloques espaciales) reúnen un conjunto de ecuaciones con distinta topología. Luego, para cada celda del espacio de parámetros, se calcula el peso de dominancia de cada ecuación mediante la probabilidad a posteriori bayesiana. En el tercer paso, una regresión simbólica secundaria encuentra leyes analíticas para la evolución de esos pesos, con restricción de monotonicidad si es necesario para procesos irreversibles. Finalmente, los pesos se extrapolan a nuevas condiciones y la predicción se construye como suma ponderada de todas las ecuaciones. El número de mecanismos K* se determina automáticamente mediante el criterio de información BIC.
En los cuatro sistemas analizados (elastómeros, ebullición, galaxias y evaporación de gotas), el BIC eligió exactamente tres ecuaciones independientes. Para los elastómeros, son la ley potencial de endurecimiento por cristalización, la ley de raíz cuadrada de elasticidad entrópica y la ley lineal de Hooke; la predicción mejoró un 72% en comparación con el gradient boosting. En la ebullición, la restricción de monotonicidad añadió un 20% de precisión. Con 3269 mediciones de 163 galaxias, el sistema encontró por sí mismo las leyes newtoniana, de transición y de «campo constante», y la diferencia del factor 2,2 en la aceleración característica a0 entre galaxias pobres y ricas en gas excluye una curva MOND única con p < 10^-34. En la evaporación de gotas, se detectó automáticamente el momento de desprendimiento de la línea de contacto.
Este trabajo se distingue de los estudios vecinos porque no busca una única fórmula óptima, sino que cartografía todo un ecosistema de mecanismos. Los trabajos adyacentes usan redes neuronales para deducir leyes o construyen mapas gravitatorios, pero todos permanecen en el paradigma de K=1. Aquí, K* se convierte en la primera magnitud medible, lo que abre una dirección ortogonal: no mejorar las ecuaciones, sino describir su coexistencia. Esto convierte el debate sobre la universalidad de una ley en una prueba computacional, no filosófica.
El tema podría evolucionar hacia una teoría dinámica de la ecología de leyes: no un mapa estático, sino una descripción de la simbiosis, la competencia y la sucesión de ecuaciones a lo largo del tiempo. Es posible ampliar el conjunto de operadores para incluir funciones trigonométricas y especiales. Aplicarlo a la turbulencia, la fatiga de materiales y la evolución estelar permitiría comprobar la universalidad del enfoque. También sería prometedor crear un punto de referencia interdisciplinario para diagnosticar el número de mecanismos en distintos campos.
El método es aplicable en mecánica de materiales, termofísica, astrofísica y microfluídica. Podría cambiar la forma de verificar leyes universales en biología y climatología.
El siguiente paso es aplicar el marco al escalado metabólico y a la sensibilidad climática, donde hay debate sobre si existe una o varias leyes. También hay que ampliar el vocabulario de operadores y mejorar el algoritmo de selección de K.
En el ámbito de la computación científica sigue abierta la cuestión de cómo determinar automáticamente el número de regímenes en los datos sin suposiciones a priori. Este trabajo propone un criterio BIC concreto y un campo continuo de pesos en lugar de límites rígidos, lo que resuelve el antiguo problema de las transiciones suaves entre mecanismos.
Una vela, que vemos como una única llama, en realidad se rige por cuatro ecuaciones distintas: conducción de calor en la cera, combustión difusiva en la zona azul, radiación en la parte brillante y convección en la lengüeta sobre la llama.