14 segundos. Eso es lo que tarda el nuevo código en tejer el espectro completo de una estrella, desde el ultravioleta hasta el infrarrojo cercano. Este espectro no es un simple arcoíris, sino un detallado lienzo compuesto por miles de oscuras líneas de absorción, cada una dejada por un átomo específico en la atmósfera de la estrella. Los astrónomos trabajan con estos lienzos como tejedores: cada línea es un hilo, y el espectro completo — la tela acabada, que almacena información sobre la temperatura, composición química y movimiento de la estrella.
Durante décadas, calcular uno de estos lienzos tomaba decenas de minutos. En la era de los grandes sondeos, que recogen millones de espectros, esto es increíblemente lento. En su ayuda llegaron los emuladores — modelos simplificados entrenados con ejemplos precalculados — y las interpolaciones sobre rejillas. Pero como una copia de una copia, heredan errores ajenos y pierden precisión cuando es necesario considerar la influencia de múltiples elementos químicos. La simulación numérica directa seguía siendo el estándar, pero demasiado costosa para su aplicación masiva.
Payne Zero cambia las reglas del juego. El código, nombrado en honor a Cecilia Payne-Gaposchkin, reescribe la lógica de la síntesis: en lugar de procesar las longitudes de onda secuencialmente, todas se manejan simultáneamente en la GPU. Es como reemplazar un telar manual por un jacquard de múltiples lanzaderas — los hilos se entrelazan a la vez. Las iteraciones atmosféricas, que requieren pasos de presión y temperatura, se aceleran entre 6 y 7 veces en CPU multinúcleo, y un inicializador entrenado reduce a la mitad el número de estas iteraciones. Resultado: la región infrarroja de APOGEE (1.5–1.7 µm) se sintetiza en 1.4 segundos.
Con esta velocidad, el análisis estelar puede, por primera vez, incorporar el ajuste físico directo en el procesamiento de millones de objetos. En una demostración con 1600 estrellas reales del sondeo APOGEE, Payne Zero determinó simultáneamente 12 elementos (carbono, nitrógeno, oxígeno, magnesio, aluminio, silicio, azufre, calcio, titanio, manganeso, níquel) y el parámetro [Fe/H]. Las secuencias químicas separaron claramente las estrellas de α alta y baja, en concordancia con el catálogo. Y todo ello con un recálculo autoconsistente de la atmósfera al final.
El toque final — la calibración de los hilos. El método utiliza la diferenciación automática de toda la cadena de transporte radiativo para ajustar simultáneamente más de cien mil correcciones a los datos atómicos, utilizando los atlas de referencia del Sol y de Arturo. La discrepancia con las observaciones se redujo en más del 80%. Ahora cada transición atómica recibe un 'hilo' calibrado, y la tela del espectro queda uniforme, sin distorsiones.
El cálculo directo óptimo dentro del bucle de ajuste elimina la necesidad misma de emuladores. Se pueden introducir cualquier elemento sin reestructurar las rejillas, y cualquier discrepancia con los datos se puede verificar a nivel de constantes atómicas o de la estructura de la atmósfera. La metalicidad y las abundancias individuales ahora se miden de manera consistente, teniendo en cuenta la influencia mutua a través de la ecuación de estado. Así se teje un lienzo unificado de la historia galáctica: más rápido, más preciso y con menores errores sistemáticos.
El proyecto comenzó como un hobby y utilizó activamente grandes modelos de lenguaje para traducir código de Fortran a Python, pero las decisiones clave sobre límites de datos y paralelización siguieron siendo humanas. La inteligencia artificial se encargó de la rutina, pero la arquitectura la dictó la experiencia con los programas originales.