Sick, the spectroscopic inference crank

Abstract

There exists an inordinate amount of spectral data in both public and private astronomical archives which remain severely under-utilised. The lack of reliable open-source tools for analysing large volumes of spectra contributes to this situation, which is poised to worsen as large surveys successively release orders of magnitude more spectra. In this Article I introduce sick, the spectroscopic inference crank, a flexible and fast Bayesian tool for inferring astrophysical parameters from spectra. sick can be used to provide a nearest-neighbour estimate of model parameters, a numerically optimised point estimate, or full Markov Chain Monte Carlo sampling of the posterior probability distributions. This generality empowers any astronomer to capitalise on the plethora of published synthetic and observed spectra, and make precise inferences for a host of astrophysical (and nuisance) quantities. Model intensities can be reliably approximated from existing grids of synthetic or observed spectra using linear multi-dimensional interpolation, or a Cannon-based model. Additional phenomena that transform the data (e.g., redshift, rotational broadening, continuum, spectral resolution) are incorporated as free parameters and can be marginalised away. Outlier pixels (e.g., cosmic rays or poorly modelled regimes) can be treated with a Gaussian mixture model, and a noise model is included to account for systematically underestimated variance. Combining these phenomena into a scalar-justified, quantitative model permits precise inferences with credible uncertainties on noisy data. Using a forward model on low-resolution, high S/N spectra of M67 stars reveals atomic diffusion processes on the order of 0.05 dex, previously only measurable with differential analysis techniques in high-resolution spectra. [abridged]

arXiv commentAccepted to ApJS

References in corpus (30)

• emcee: The MCMC Hammer
• PARSEC: stellar tracks and isochrones with the PAdova and TRieste Stellar Evolution Code
• A new extensive library of PHOENIX stellar atmospheres and synthetic spectra
• SEGUE: A Spectroscopic Survey of 240,000 stars with g=14-20
• The Radial Velocity Experiment (RAVE): first data release
• The GALAH Survey: Scientific Motivation
• Na-O Anticorrelation and HB. VII. The chemical composition of first and second-generation stars in 15 globular clusters from GIRAFFE spectra
• The WiggleZ Dark Energy Survey: Survey Design and First Data Release
• Determining the Type, Redshift, and Age of a Supernova Spectrum
• The SEGUE Stellar Parameter Pipeline. I. Description and Initial Validation Tests
• The RAdial Velocity Experiment (RAVE): Fourth data release
• Evolutionary synthesis of galaxies at high spectral resolution with the code PEGASE-HR
• Improved spectroscopic parameters for transiting planet hosts
• DAOSPEC: an automatic code for measuring equivalent widths in high-resolution stellar spectra
• The SEGUE Stellar Parameter Pipeline. III. Comparison with High-Resolution Spectroscopy of SDSS/SEGUE Field Stars
• Gaia-ESO Survey: The analysis of high-resolution UVES spectra of FGK-type stars
• Elemental Abundance Ratios in Stars of the Outer Galactic Disk. I. Open Clusters
• POLLUX : a database of synthetic stellar spectra
• Automated derivation of stellar atmospheric parameters and chemical abundances: the MATISSE algorithm
• The AMBRE project: A new synthetic grid of high-resolution FGKM stellar spectra
• Solar twins in M67
• Abundances and possible diffusion of elements in M67 stars
• A High Relative-Precision Color-Magnitude Diagram of M67
• Abundances of four open clusters from solar stars
• M67-1194, an unusually Sun-like solar twin in M67
• Deep 2MASS Photometry of M67 and Calibration of the Main Sequence J-Ks Color Difference as an Age Indicator
• A High-Resolution Spectroscopic Search for the Remaining Donor for Tycho's Supernova
• Close binary stars in the solar-age Galactic open cluster M67
• Tracing mixing in stars: new beryllium observations of the open clusters NGC 2516, Hyades, and M67
• An Effective temperature calibration for solar type stars using equivalent width ratios - A fast and easy spectroscopic temperature estimation

Cited by in corpus (9)

• Modern stellar spectroscopy caveats
• The EMBLA Survey -- Metal-poor stars in the Galactic bulge
• The RAVE-on catalog of stellar atmospheric parameters and chemical abundances for chemo-dynamic studies in the Gaia era
• The Gaia-ESO Survey: Evidence of atomic diffusion in M67?
• Using the Multi-Object Adaptive Optics demonstrator RAVEN to observe metal-poor stars in and towards the Galactic Centre
• The Most Metal-poor Stars in the Inner Bulge
• Mapping the distribution of OB stars and associations in Auriga
• An Interpretable Machine Learning Framework for Modeling High-Resolution Spectroscopic Data
• Analytic marginalization of absorption line continua