Constructing a flexible likelihood function for spectroscopic inference

Ian Czekala, Sean M. Andrews, Kaisey S. Mandel, David W. Hogg, Gregory M. Green

    Research output: Contribution to journalArticlepeer-review


    We present a modular, extensible likelihood framework for spectroscopic inference based on synthetic model spectra. The subtraction of an imperfect model from a continuously sampled spectrum introduces covariance between adjacent datapoints (pixels) into the residual spectrum. For the high signal-to-noise data with large spectral range that is commonly employed in stellar astrophysics, that covariant structure can lead to dramatically underestimated parameter uncertainties (and, in some cases, biases). We construct a likelihood function that accounts for the structure of the covariance matrix, utilizing the machinery of Gaussian process kernels. This framework specifically addresses the common problem of mismatches in model spectral line strengths (with respect to data) due to intrinsic model imperfections (e.g., in the atomic/molecular databases or opacity prescriptions) by developing a novel local covariance kernel formalism that identifies and self-consistently downweights pathological spectral line "outliers." By fitting many spectra in a hierarchical manner, these local kernels provide a mechanism to learn about and build data-driven corrections to synthetic spectral libraries. An open-source software implementation of this approach is available at, including a sophisticated probabilistic scheme for spectral interpolation when using model libraries that are sparsely sampled in the stellar parameters. We demonstrate some salient features of the framework by fitting the high-resolution V-band spectrum of WASP-14, an F5 dwarf with a transiting exoplanet, and the moderate-resolution K-band spectrum of Gliese 51, an M5 field dwarf.

    Original languageEnglish (US)
    Article number128
    JournalAstrophysical Journal
    Issue number2
    StatePublished - Oct 20 2015


    • methods: data analysis
    • methods: statistical
    • stars: fundamental parameters
    • stars: late-type
    • stars: statistics
    • techniques: spectroscopic

    ASJC Scopus subject areas

    • Astronomy and Astrophysics
    • Space and Planetary Science


    Dive into the research topics of 'Constructing a flexible likelihood function for spectroscopic inference'. Together they form a unique fingerprint.

    Cite this