Reading Between the Lines: How Much Information Remains After Fitting a Stellar Spectrum?
Abstract
What if we could model 100% of a stellar spectrum? Modern surveys infer stellar properties by comparing observed spectra with synthetic models, yet even the best fits leave structured residuals. We analyse 7087 GALAH survey spectra and ask what information remains after subtracting each best-fitting model. Across 13927 pixels, 4.5% of residuals exceed 3σ, far above 0.3% expected from Gaussian noise. We map the amplitude, significance, sign, and recurrence of these residuals and identify 380 significant regions. Comparison with stellar atlases, line lists, diffuse interstellar bands, and sky diagnostics links most regions to plausible causes: 62% are candidate line-strength or line-profile mismatches, while others trace missing or unidentified lines, interstellar absorption, sky emission, or broader reduction structure. We then train data-driven models on the residual spectra themselves. Residual-only models recover Teff to 80-150 K and [Fe/H] to 0.059-0.086 dex, depending on sample quality, and retain weaker, partly indirect information about elemental abundances. Extinction-related coefficients recover the strongest known interstellar features, whereas chromospheric activity leaves plausible Balmer-line signatures but is not recovered robustly out of sample. Most strikingly, the residual flux in 6% of wavelength pixels that were synthesised but excluded from the GALAH DR4 optimisation recovers stellar labels almost as well as the full residual spectrum and better than the remaining 94% alone. Residuals are therefore not merely imperfections to discard: they map missing line data, incomplete physics, unfitted astrophysical components, and observational systematics. Reading between the lines turns one spectroscopist's trash into another one's treasure, and provides a practical route towards extracting more of the information already present in stellar spectra.
Related
- FastSpec : Fast and versatile stellar spectrum synthesis code
- Baryonic feedback suppression of the matter power spectrum: a three-parameter fitting formula and its single-parameter reduction
- Unlocking new information in stellar circular polarisation with second-order least-squares deconvolution
- POSTELLAR: Posterior Stellar Spectrum Sampling - An Alternative to Approximate Stellar Spectra for Exoplanetary Analysis
- Planetesimal Formation Across the Stellar Mass Spectrum and its Influence on Exoplanet-Inherited Volatile Budgets
BibTeX
@misc{buder2026reading,
title = {{Reading Between the Lines: How Much Information Remains After Fitting a Stellar Spectrum?}},
author = {Sven Buder and Tomasz Różánski},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.25583},
eprint = {2610.25583},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.25583},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Reading Between the Lines: How Much Information Remains After Fitting a Stellar Spectrum?
AU - Sven Buder
AU - Tomasz Różánski
PY - 2026
DA - 2026/10/18
PB - arXiv
DO - 10.48550/arxiv.2610.25583
UR - https://arxiv.org/abs/2610.25583
AB - What if we could model 100% of a stellar spectrum? Modern surveys infer stellar properties by comparing observed spectra with synthetic models, yet even the best fits leave structured residuals. We analyse 7087 GALAH survey spectra and ask what information remains after subtracting each best-fitting model. Across 13927 pixels, 4.5% of residuals exceed 3σ, far above 0.3% expected from Gaussian noise. We map the amplitude, significance, sign, and recurrence of these residuals and identify 380 significant regions. Comparison with stellar atlases, line lists, diffuse interstellar bands, and sky diagnostics links most regions to plausible causes: 62% are candidate line-strength or line-profile mismatches, while others trace missing or unidentified lines, interstellar absorption, sky emission, or broader reduction structure. We then train data-driven models on the residual spectra themselves. Residual-only models recover T_eff to 80-150 K and [Fe/H] to 0.059-0.086 dex, depending on sample quality, and retain weaker, partly indirect information about elemental abundances. Extinction-related coefficients recover the strongest known interstellar features, whereas chromospheric activity leaves plausible Balmer-line signatures but is not recovered robustly out of sample. Most strikingly, the residual flux in 6% of wavelength pixels that were synthesised but excluded from the GALAH DR4 optimisation recovers stellar labels almost as well as the full residual spectrum and better than the remaining 94% alone. Residuals are therefore not merely imperfections to discard: they map missing line data, incomplete physics, unfitted astrophysical components, and observational systematics. Reading between the lines turns one spectroscopist's trash into another one's treasure, and provides a practical route towards extracting more of the information already present in stellar spectra.
ER -