Measuring the Variability of Hydroxyl Emissions in Infrared Sky Spectra using SPIRou: Dataset

Abstract

Subtracting the changing sky contribution from the near-infrared (NIR) spectra of faint astronomicalobjects is challenging and crucial to a wide range of science cases such as estimating the velocitydispersions of dwarf galaxies, studying the gas dynamics in faint galaxies, and accurate redshifts, andany spectroscopic studies of faint targets. Since the sky background varies with time and location, NIRspectral observations, especially those employing fiber spectrometers and targeting extended sources,require frequent sky-only observations for calibration. However, sky subtraction can be optimizedwith sufficient a priori knowledge of the sky’s variability. In this work, we explore how to optimize skysubtraction by analyzing 1075 high-resolution NIR spectra from the CFHT’s SPIRou on Maunakea,and we estimate the variability of 481 hydroxyl (OH) lines. These spectra were collected during twosets of three nights dedicated to obtaining sky observations every five and a half minutes. Duringthe first set, we observed how the Moon affects the NIR, which has not been accurately measured atthese wavelengths. We suggest that if one uses a principal component analysis reconstruction of thesky spectrum and attempts to observe targets at Y JHK mags fainter than ∼15 and attempts a skysubtraction better than 1%, then the Moon contribution must be accounted for at Moon separationdistances of at least 10◦. We also identified 126 spectral doublet, or OH lines that split into at least twocomponents, at SPIRou’s resolution. In addition, we used Lomb-Scargle Periodograms and Gaussianprocess regression to estimate most OH lines vary on similar timescales, which provides a valuableinput for IR spectroscopic survey strategies. The data and code developed for this study are publiclyavailable here.

Dataset description

In total, we collected 1075 sky observations, which spanned from July 28th, 2018 to January 10th, 2022, or approximately 3.5 years. These observations included two sets of three days dedicated to sky measurement where each day, a sky spectrum was observed approximately every 5.5 minutes for 12 hours. These days occurred on December 14th, 15th, and 16th of 2019 and January 22nd, 23rd, and 25th 2020. These 1D extracted and flat fielded sky spectra have a wavelength range of 0.965 − 2.500μm containing 285,377 wavelength bins resampled on a uniform wavelength grid with a step of 1 km/s/pixel. Since the pixels were constant in velocity, the change in wavelength increased from 3x10^−6 − 8x10^−6 μm per pixel. All observations were affected by a steep black body curve starting at 2.1μm, which was caused by thermal emission.

We also present a table of doublets identified during this study with the transition, the measured singlet line from Rousselot 2000 (mu0), the doublet lines (mu1 and mu2), and if the line was identified as a doublet (Y/N).