skip to main content


Title: Data-driven Discovery of Diffuse Interstellar Bands with APOGEE Spectra
Abstract

Data-driven models of stellar spectra are useful tools to study nonstellar information, such as the diffuse interstellar bands (DIBs) caused by intervening interstellar material. Using ∼55,000 spectra of ∼17,000 red clump stars from the APOGEE DR16 data set, we create second-order polynomial models of the continuum-normalized flux as a function of stellar parameters (Teff,logg, [Fe/H], [α/Fe], and age). The model and data show good agreement within uncertainties across the APOGEE wavelength range, although many regions reveal residuals that are not in the stellar rest-frame. We show that many of these residual features—having average extrema at the level of ∼3% in stellar flux on average—can be attributed to incompletely removed spectral lines from the Earth’s atmosphere and DIBs from the interstellar medium (ISM). After removing most of the remaining contamination from Earth’s sky, we identify 84 absorption features not seen in unreddened sightlights that have <50% probability of being noise artifacts—with 25 of these features having <5% probability of being noise artifacts—including all 10 previously known DIBs in the APOGEE wavelength range. Because many of these features occur in the wavelength windows that APOGEE uses to measure chemical abundances, note that characterization and removal of this nonstellar contamination establish an important step in reaching the precision required for chemical tagging experiments. Proper characterization of these features will benefit Galactic ISM science and the currently ongoing Milky Way Mapper program of Sloan Digital Sky Survey V, which relies on the APOGEE spectrograph.

 
more » « less
NSF-PAR ID:
10500226
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
965
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 120
Size(s):
["Article No. 120"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) Stellar Library (MaStar) is a large collection of high-quality empirical stellar spectra designed to cover all spectral types and ideal for use in the stellar population analysis of galaxies observed in the MaNGA survey. The library contains 59,266 spectra of 24,130 unique stars with spectral resolutionR∼ 1800 and covering a wavelength range of 3622–10,354 Å. In this work, we derive five physical parameters for each spectrum in the library: effective temperature (Teff), surface gravity (logg), metallicity ([Fe/H]), microturbulent velocity (log(vmicro)), and alpha-element abundance ([α/Fe]). These parameters are derived with a flexible data-driven algorithm that uses a neural network model. We train a neural network using the subset of 1675 MaStar targets that have also been observed in the Apache Point Observatory Galactic Evolution Experiment (APOGEE), adopting the independently-derived APOGEE Stellar Parameter and Chemical Abundance Pipeline parameters for this reference set. For the regions of parameter space not well represented by the APOGEE training set (7000 ≤T≤ 30,000 K), we supplement with theoretical model spectra. We present our derived parameters along with an analysis of the uncertainties and comparisons to other analyses from the literature.

     
    more » « less
  2. Abstract

    The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work, we focus on the most massive galaxies in the nearby universelog(M/M)>11.2. We obtain high-quality Magellan/LDSS-3 long-slit spectroscopy with a wide wavelength coverage of 0.4–1.01μm for 41 early-type galaxies (ETGs) in the MASSIVE survey and derive high signal-to-noise spectra within an aperture ofRe/8. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter,αIMF= (M/L)/(M/L)MW, ranges from 1.1 to 3.1, with an average of 〈αIMF〉 = 1.84, suggesting that on average, the IMF is more bottom heavy than Salpeter. Comparing the estimated stellar masses with the dynamical masses, we find that most galaxies have stellar masses that are smaller than their dynamical masses within the 1σuncertainty. We complement our sample with lower-mass galaxies from the literature and confirm thatlog(αIMF)is positively correlated withlog(σ),log(M), andlog(Mdyn). From the combined sample, we show that the IMF in the centers of more massive ETGs is more bottom heavy. In addition, we find thatlog(αIMF)is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density ΣKroupamight be responsible for the variation ofαIMF. We conclude thatσ, [Mg/Fe], and [Z/H] are the primary drivers of the global stellar IMF variation.

     
    more » « less
  3. One of the top priorities in observational astronomy is the direct imaging and characterization of extrasolar planets (exoplanets) and planetary systems. Direct images of rocky exoplanets are of particular interest in the search for life beyond the Earth, but they tend to be rather challenging targets since they are orders-of-magnitude dimmer than their host stars and are separated by small angular distances that are comparable to the classicalλ<#comment/>/Ddiffraction limit, even for the coming generation of 30 m class telescopes. Current and planned efforts for ground-based direct imaging of exoplanets combine high-order adaptive optics (AO) with a stellar coronagraph observing at wavelengths ranging from the visible to the mid-IR. The primary barrier to achieving high contrast with current direct imaging methods is quasi-static speckles, caused largely by non-common path aberrations (NCPAs) in the coronagraph optical train. Recent work has demonstrated that millisecond imaging, which effectively “freezes” the atmosphere’s turbulent phase screens, should allow the wavefront sensor (WFS) telemetry to be used as a probe of the optical system to measure NCPAs. Starting with a realistic model of a telescope with an AO system and a stellar coronagraph, this paper provides simulations of several closely related regression models that take advantage of millisecond telemetry from the WFS and coronagraph’s science camera. The simplest regression model, called the naïve estimator, does not treat the noise and other sources of information loss in the WFS. Despite its flaws, in one of the simulations presented herein, the naïve estimator provides a useful estimate of an NCPA of∼<#comment/>0.5radian RMS (≈<#comment/>λ<#comment/>/13), with an accuracy of∼<#comment/>0.06radian RMS in 1 min of simulated sky time on a magnitude 8 star. Thebias-corrected estimatorgeneralizes the regression model to account for the noise and information loss in the WFS. A simulation of the bias-corrected estimator with 4 min of sky time included an NCPA of∼<#comment/>0.05radian RMS (≈<#comment/>λ<#comment/>/130) and an extended exoplanet scene. The joint regression of the bias-corrected estimator simultaneously achieved an NCPA estimate with an accuracy of∼<#comment/>5×<#comment/>10−<#comment/>3radian RMS and an estimate of the exoplanet scene that was free of the self-subtraction artifacts typically associated with differential imaging. The5σ<#comment/>contrast achieved by imaging of the exoplanet scene was∼<#comment/>1.7×<#comment/>10−<#comment/>4at a distance of3λ<#comment/>/Dfrom the star and∼<#comment/>2.1×<#comment/>10−<#comment/>5at10λ<#comment/>/D. These contrast values are comparable to the very best on-sky results obtained from multi-wavelength observations that employ both angular differential imaging (ADI) and spectral differential imaging (SDI). This comparable performance is despite the fact that our simulations are quasi-monochromatic, which makes SDI impossible, nor do they have diurnal field rotation, which makes ADI impossible. The error covariance matrix of the joint regression shows substantial correlations in the exoplanet and NCPA estimation errors, indicating that exoplanet intensity and NCPA need to be estimated self-consistently to achieve high contrast.

     
    more » « less
  4. Abstract

    We present rest-frame optical emission-line flux ratio measurements for fivez> 5 galaxies observed by the James Webb Space Telescope Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliablerelativeflux calibration of emission lines that are closely separated in wavelength, despite the uncertainabsolutespectrophotometry of the current version of the reductions. Compared toz∼ 3 galaxies in the literature, thez> 5 galaxies have similar [Oiii]λ5008/Hβratios, similar [Oiii]λ4364/Hγratios, and higher (∼0.5 dex) [NeIII]λ3870/[OII]λ3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]λ3870/[OII]λ3728, [Oiii]λ4364/Hγ, and [Oiii]λ5008/Hβemission-line ratios are consistent with an interstellar medium (ISM) that has very high ionization (log(Q)89, units of cm s−1), low metallicity (Z/Z≲ 0.2), and very high pressure (log(P/k)89, units of cm−3). The combination of [Oiii]λ4364/Hγand [Oiii]λ(4960 + 5008)/Hβline ratios indicate very high electron temperatures of4.1<log(Te/K)<4.4, further implying metallicities ofZ/Z≲ 0.2 with the application of low-redshift calibrations for “Te-based” metallicities. These observations represent a tantalizing new view of the physical conditions of the ISM in galaxies at cosmic dawn.

     
    more » « less
  5. Abstract

    We study the kinematics of the interstellar medium (ISM) viewed “down the barrel” in 20 gravitationally lensed galaxies during cosmic noon (z= 1.5–3.5). We use moderate-resolution spectra (R∼ 4000) from Keck’s Echellette Spectrograph and Imager and Magellan/MagE to spectrally resolve the ISM absorption in these galaxies into ∼10 independent elements and use double Gaussian fits to quantify the velocity structure of the gas. We find that the bulk motion of gas in this galaxy sample is outflowing, with average velocity centroidvcent=141km s−1(±111 km s−1scatter) measured with respect to the systemic redshift. A total of 16 out of the 20 galaxies exhibit a clear positive skewness, with a blueshifted tail extending to ∼ −500 km s−1. We examine scaling relations in outflow velocities with galaxy stellar mass and star formation rate, finding correlations consistent with a momentum-driven wind scenario. Our measured outflow velocities are also comparable to those reported for FIRE-2 and TNG50 cosmological simulations at similar redshift and galaxy properties. We also consider implications for interpreting results from lower-resolution spectra. We demonstrate that while velocity centroids are accurately recovered, the skewness, velocity width, and probes of high-velocity gas (e.g.,v95) are subject to large scatter and biases at lower resolution. We find thatR≳ 1700 is required for accurate results for the gas kinematics of our sample. This work represents the largest available sample of well-resolved outflow velocity structure atz> 2 and highlights the need for good spectral resolution to recover accurate properties.

     
    more » « less