Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra: Application to Chemical Homogeneity in the Open Cluster M67

Abstract

High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archaeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of $18,933$ giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires $\approx$10 functional principal components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: $\mathrm{Fe}\lesssim 0.02$ dex; $\mathrm{C}\lesssim 0.03$ dex; $\mathrm{O},\mathrm{Mg},\mathrm{Si},\mathrm{Ni}\lesssim 0.04$ dex; $\mathrm{Ca}\lesssim 0.05$ dex; $\mathrm{N},\mathrm{Al}\lesssim 0.07$ dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way.

30 pages, 16 figures, 5 tables, Accepted to ApJ, Code available at https://github.com/aaryapatil/specdims

References in corpus (26)

• Array Programming with NumPy
• The 2.5 m Telescope of the Sloan Digital Sky Survey
• The GALAH Survey: Scientific Motivation
• The Apache Point Observatory Galactic Evolution Experiment (APOGEE) Spectrographs
• 4MOST: Project overview and information for the First Call for Proposals
• Target Selection for the SDSS-IV APOGEE-2 Survey
• Atomic Diffusion and Mixing in Old Stars I. VLT/FLAMES-UVES Observations of Stars in NGC 6397
• Fast likelihood-free cosmology with neural density estimators and active learning
• Inferring the eccentricity distribution
• Exploring Anticorrelations and Light Element Variations in Northern Globular Clusters Observed by the APOGEE Survey
• Non-LTE aluminium abundances in late-type stars
• Chemical Homogeneity in Collinder 261 and Implications for Chemical Tagging
• Chemically Tagging the HR1614 Moving Group
• SDSS-V: Pioneering Panoptic Spectroscopy
• Chemical Abundances of Main-Sequence, Turnoff, Subgiant and red giant Stars from APOGEE spectra II: Atomic Diffusion in M67 Stars
• Early turbulent mixing as the origin of chemical homogeneity in open star clusters
• A Hierarchical Bayesian SED Model for Type Ia Supernovae in the Optical to Near-Infrared
• Mapping the Interstellar Medium with Near-Infrared Diffuse Interstellar Bands
• Chemical (in)homogeneity and atomic diffusion in the open cluster M67
• Atomic diffusion and mixing in old stars II. Observations of stars in the globular cluster NGC 6397 with VLT/FLAMES-GIRAFFE
• How Many Elements Matter?
• Open Cluster Chemical Homogeneity Throughout the Milky Way
• Blind chemical tagging with DBSCAN: prospects for spectroscopic surveys
• The GALAH survey: Chemical homogeneity of the Orion complex
• Disentangled Representation Learning for Astronomical Chemical Tagging
• Neural Density Estimation and Likelihood-free Inference