How to Weave a Star's Light in 14 Seconds ✛

Original: "The Payne Zero Project I: Stellar Spectra from Physical Models in Seconds"
· Yuan-Sen Ting, Elliot M. Kim
arXiv:2607.24141v1 · 2026-07-27 · CC BY 4.0 · 2 min
A new method shifts spectrum calculations onto graphics cards, replacing approximate estimates with a physical model of the star.

Every spectrum of a star is like a tapestry woven from thousands of fine absorption lines. These lines are left by atoms in the star's atmosphere, and from their pattern astronomers learn the temperature, composition, and motion of the star. But such a tapestry is hard to 'weave': old programs did it sequentially, thread by thread, taking nearly an hour for a single star.

Payne Zero moves this process onto a graphics processing unit (GPU), which can perform many operations simultaneously — like a loom with millions of needles. Now all spectral lines are calculated at once. The atmospheric model, which sets the tension and color of the threads, is accelerated on a multi-core CPU by a factor of 6–7, and a trained assistant picks initial conditions, reducing the number of iterations. The result: a full optical spectrum of the Sun is ready in 14 seconds — 54 times faster than before.

14 seconds — that's how long the new code takes to weave a complete portrait of a star, from the violet to the dark red edge.

This speed makes it possible, for the first time, to embed direct physical calculations into the study of millions of stars. In a test on 1600 stars from the APOGEE survey, the code measured the abundances of 12 elements — from carbon to nickel — in a single pass and clearly separated the stars into two chemical groups, as expected. No more relying on approximate surrogates: each star's spectrum is woven anew from atomic 'threads', calibrated using reference stars like the Sun. This approach promises to recalibrate the atomic data itself and reveal the detailed chemical evolution history of our Galaxy.

The project began as a hobby and made active use of large language models to translate code from Fortran to Python, but key decisions about data boundaries and parallelization remained with the human. Artificial intelligence translated the routine, but the architecture was guided by experience with the original programs.

\lambda_{i_\lambda} = \lambda_0 \left(1 + R^{-1}_{\rm grid}\right)^{i_\lambda}
Logarithmically uniform sampling, where R is the number of points per e-fold change in wavelength.
F_{\rm rad} + F_{\rm conv} = \sigma_{\rm SB} T_{\rm eff}^4
Energy balance in the atmosphere: the sum of radiative and convective fluxes equals the total blackbody emission.
Scientists
Christian DopplerD. B. McLaughlinDidier QuelozMichel MayorR. A. RossiterAlbert Einstein
Concepts
gpu parker solar probe pump probe spectroscopy single pixel imaging sdss square degree GPCP