Constructing Polynomial Spectral Models for Stars
Abstract
Stellar spectra depend on the stellar parameters and on dozens of photospheric elemental abundances. Simultaneous fitting of these model labels to observed spectra has been deemed unfeasible, because the number of ab initio spectral model grid calculations scales exponentially with . We suggest instead the construction of a polynomial spectral model (PSM) of order for the model flux at each wavelength. Building this approximation requires a minimum of only calculations: e.g. a quadratic spectral model () to fit labels simultaneously, can be constructed from as few as ab initio spectral model calculations; in practice, a somewhat larger number () of randomly chosen models lead to a better performing PSM. Such a PSM can be a good approximation only over a portion of label space, which will vary case by case. Yet, taking the APOGEE survey as an example, a single quadratic PSM provides a remarkably good approximation to the exact ab initio spectral models across much of this survey: for random labels within that survey the PSM approximates the flux to within , and recovers the abundances to within dex rms of the exact models. This enormous speed-up enables the simultaneous many-label fitting of spectra with computationally expensive ab initio models for stellar spectra, such as non-LTE models. A PSM also enables the simultaneous fitting of observational parameters, such as the spectrum's continuum or line-spread function.
Cite
@article{arxiv.1603.06574,
title = {Constructing Polynomial Spectral Models for Stars},
author = {Hans-Walter Rix and Yuan-Sen Ting and Charlie Conroy and David W. Hogg},
journal= {arXiv preprint arXiv:1603.06574},
year = {2016}
}
Comments
4 pages, 2 figures, ApJL (Accepted for publication- 2016 May 9)