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 (
Diffuse interstellar bands (DIBs) are broad absorption features associated with interstellar dust and can serve as chemical and kinematic tracers. Conventional measurements of DIBs in stellar spectra are complicated by residuals between observations and best-fit stellar models. To overcome this, we simultaneously model the spectrum as a combination of stellar, dust, and residual components, with full posteriors on the joint distribution of the components. This decomposition is obtained by modeling each component as a draw from a high-dimensional Gaussian distribution in the data space (the observed spectrum)—a method we call “Marginalized Analytic Data-space Gaussian Inference for Component Separation” (MADGICS). We use a data-driven prior for the stellar component, which avoids missing stellar features not well modeled by synthetic spectra. This technique provides statistically rigorous uncertainties and detection thresholds, which are required to work in the low signal-to-noise regime that is commonplace for dusty lines of sight. We reprocess all public Gaia DR3 RVS spectra and present an improved 8621 Å DIB catalog, free of detectable stellar line contamination. We constrain the rest-frame wavelength to 8623.14 ± 0.087 Å (vacuum), find no significant evidence for DIBs in the Local Bubble from the 1/6th of RVS spectra that are public, and show unprecedented correlation with kinematic substructure in Galactic CO maps. We validate the catalog, its reported uncertainties, and biases using synthetic injection tests. We believe MADGICS provides a viable path forward for large-scale spectral line measurements in the presence of complex spectral contamination.
more » « less- NSF-PAR ID:
- 10451599
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 954
- Issue:
- 2
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 141
- Size(s):
- ["Article No. 141"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract T eff,α /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
Abstract Precise radial velocity (PRV) surveys are important for the search for Earth analogs around nearby bright stars, which induce a small stellar reflex motion with an RV amplitude of ∼10 cm s−1. Detecting such a small RV signal poses challenges to instrumentation, data analysis, and the precision of astrophysical models to mitigate stellar jitter. In this work, we investigate an important component in the PRV error budget—the spectral contamination from the Earth’s atmosphere (tellurics). We characterize the effects of telluric absorption on the RV precision and quantify its contribution to the RV error budget over time and across a wavelength range of 350 nm–2.5
μ m. We use simulated solar spectra with telluric contamination injected, and we extract the RVs using two commonly adopted algorithms: dividing out a telluric model before performing cross-correlation or forward modeling the observed spectrum incorporating a telluric model. We assume various degrees of cleanness in removing the tellurics. We conclude that the RV errors caused by telluric absorption can be suppressed to close to or even below 1–10 cm s−1in the blue optical region. At red through near-infrared wavelengths, however, the residuals of tellurics can induce an RV error on the meter-per-second level even under the most favorable assumptions for telluric removal, leading to significant systematic noise in the RV time series and periodograms. If the red-optical or near-infrared becomes critical in the mitigation of stellar activity, systematic errors from tellurics can be eliminated with a space mission such as EarthFinder. -
more » « less
We present Atacama Large Millimeter/submillimeter Array (ALMA) sub-kiloparsec- to kiloparsec-scale resolution observations of the [C
II ], CO (9–8), and OH+(11–01) lines along with their dust continuum emission toward the far-infrared (FIR) luminous quasar SDSS J231038.88+185519.7 atz = 6.0031, to study the interstellar medium distribution, the gas kinematics, and the quasar-host system dynamics. We decompose the intensity maps of the [CII ] and CO (9–8) lines and the dust continuum with two-dimensional elliptical Sérsic models. The [CII ] brightness follows a flat distribution with a Sérsic index of 0.59. The CO (9–8) line and the dust continuum can be fit with an unresolved nuclear component and an extended Sérsic component with a Sérsic index of ∼1, which may correspond to the emission from an active galactic nucleus dusty molecular torus and a quasar host galaxy, respectively. The different [CII ] spatial distribution may be due to the effect of the high dust opacity, which increases the FIR background radiation on the [CII ] line, especially in the galaxy center, significantly suppressing the [CII ] emission profile. The dust temperature drops with distance from the center. The effective radius of the dust continuum is smaller than that of the line emission and the dust mass surface density, but is consistent with that of the star formation rate surface density. This may indicate that the dust emission is a less robust tracer of the dust and gas distribution but is a decent tracer of the obscured star formation activity. The OH+(11–01) line shows a P-Cygni profile with an absorption at ∼–400 km s−1, which may indicate an outflow with a neutral gas mass of (6.2 ± 1.2)×108M ⊙along the line of sight. We employed a three-dimensional tilted ring model to fit the [CII ] and CO (9–8) data cubes. The two lines are both rotation dominated and trace identical disk geometries and gas motions. This suggest that the [CII ] and CO (9–8) gas are coplanar and corotating in this quasar host galaxy. The consistent circular velocities measured with [CII ] and CO (9–8) lines indicate that these two lines trace a similar gravitational potential. We decompose the circular rotation curve measured from the kinematic model fit to the [CII ] line into four matter components (black hole, stars, gas, and dark matter). The quasar-starburst system is dominated by baryonic matter inside the central few kiloparsecs. We constrain the black hole mass to be 2.97+0.51-0.77 × 109M ⊙; this is the first time that the dynamical mass of a black hole has been measured atz ∼ 6. This mass is consistent with that determined using the scaling relations from quasar emission lines. A massive stellar component (on the order of 109M ⊙) may have already existed when the Universe was only ∼0.93 Gyr old. The relations between the black hole mass and the baryonic mass of this quasar indicate that the central supermassive black hole may have formed before its host galaxy. -
Abstract Photoionization modeling of active galactic nuclei (AGN) predicts that diffuse continuum (DC) emission from the broad-line region makes a substantial contribution to the total continuum emission from ultraviolet through near-infrared wavelengths. Evidence for this DC component is present in the strong Balmer jump feature in AGN spectra, and possibly from reverberation measurements that find longer lags than expected from disk emission alone. However, the Balmer jump region contains numerous blended emission features, making it difficult to isolate the DC emission strength. In contrast, the Paschen jump region near 8200 Å is relatively uncontaminated by other strong emission features. Here, we examine whether the Paschen jump can aid in constraining the DC contribution, using Hubble Space Telescope Space Telescope Imaging Spectrograph spectra of six nearby Seyfert 1 nuclei. The spectra appear smooth across the Paschen edge, and we find no evidence of a Paschen spectral break or jump in total flux. We fit multicomponent spectral models over the range 6800–9700 Å and find that the spectra can still be compatible with a significant DC contribution if the DC Paschen jump is offset by an opposite spectral break resulting from blended high-order Paschen emission lines. The fits imply DC contributions ranging from ∼10% to 50% at 8000 Å, but the fitting results are highly dependent on assumptions made about other model components. These degeneracies can potentially be alleviated by carrying out fits over a broader wavelength range, provided that models can accurately represent the disk continuum shape, Fe ii emission, high-order Balmer line emission, and other components.more » « less
-
Context. Cool stars, such as M giants, can only be analyzed in the near-infrared (NIR) regime due to the ubiquitous titanium oxide features in optical spectra of stars with T eff < 4000 K. In dust-obscured regions, the inner bulge and Galactic center region, the intrinsically bright M giants observed in the NIR are an optimal option for studying stellar abundances and the chemical evolution of stellar populations. Because of the uncertainties in photometric methods, a method for determining the stellar parameters for M giants from the NIR spectra themselves is needed. Aims. We develop a method for determining the stellar parameters for M giants from the NIR spectra. We validate the method by deriving the stellar parameters for nearby well-studied M giants with spectra from the spectral library of the Immersion GRating INfrared Spectrograph (IGRINS). We demonstrate the accuracy and precision of our method by determining the stellar parameters and α -element trends versus metallicity for solar neighborhood M giants. Methods. We carried out new observations of 44 M giant stars with IGRINS mounted on the Gemini South telescope. We also obtained the full H and K band IGRINS spectra of six nearby well-studied M giants at a spectral resolving power of R = 45 000 from the IGRINS spectral library. We used the tool called spectroscopy made easy in combination with one-dimensional (1D) model atmospheres in a radiative and convective scheme (MARCS) stellar atmosphere models to model the synthetic spectrum that fits the observed spectrum best. Results. The effective temperatures that we derive from our new method (tested for 3400 ≲ T eff ≲ 4000 K here) agree excellently with those of the six nearby well-studied M giants, which indicates that the accuracy is indeed high. For the 43 solar neighborhood M giants, our T eff , log g , [Fe/H], ξ micro , [C/Fe], [N/Fe], and [O/Fe] agree with APOGEE with mean differences and a scatter (our method – APOGEE) of −67±33 K, −0.31±0.15 dex, 0.02±0.05 dex, 0.22±0.13 km s −1 , −0.05±0.06 dex, 0.06±0.06 dex, and 0.02±0.09 dex, respectively. Furthermore, the tight offset with a small dispersion compared to the APOGEE T eff indicates a high precision in our derived temperatures and those derived from the APOGEE pipeline. The typical uncertainties in the stellar parameters are found to be ±100 K in T eff , ±0.2 dex in log g , ±0.1 dex in [Fe/H], and ±0.1 km s −1 in ξ micro . The α -element trends versus metallicity for Mg, Si, Ca, and Ti are consistent with the APOGEE DR17 trends for the same stars and with the GILD optical trends. We also find a clear enhancement in the abundances for thick-disk stars.more » « less
