# The X-shooter Spectral Library (XSL): Data release 2

A. Gonneau, M. Lyubenova, A. Lançon, S. Trager, R. Peletier, A. Arentsen, Y.-P. Chen, P. Coelho, M. Dries, J. Falcón-Barroso, et al.

► To cite this version:

A. Gonneau, M. Lyubenova, A. Lançon, S. Trager, R. Peletier, et al. The X-shooter Spectral Library (XSL): Data release 2. Astronomy and Astrophysics - A&A, 2020, 634, pp.A133. 10.1051/0004-6361/201936825. hal-03152017

HAL Id: hal-03152017

https://hal.science/hal-03152017

Submitted on 28 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
# The X-shooter Spectral Library (XSL): Data release 2 ★★★

A. Gonneau<sup>1,2,4</sup>, M. Lyubenova<sup>3,2</sup>, A. Lançon<sup>4</sup>, S. C. Trager<sup>2</sup>, R. F. Peletier<sup>2</sup>, A. Arentsen<sup>7,2</sup>, Y.-P. Chen<sup>6</sup>, P. R. T. Coelho<sup>11</sup>, M. Dries<sup>2</sup>, J. Falcón-Barroso<sup>8,9</sup>, P. Prugniel<sup>5</sup>, P. Sánchez-Blázquez<sup>10</sup>, A. Vazdekis<sup>8,9</sup>, and K. Verro<sup>2</sup>

¹ Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
e-mail: agonneau@ast.cam.ac.uk

² Kapteyn Astronomical Institute, University of Groningen, Landleven 12, 9747 AD Groningen, The Netherlands

³ ESO, Karl-Schwarzschild-Str. 2, 85748 Garching bei München, Germany

⁴ Observatoire Astronomique de Strasbourg, Université de Strasbourg, CNRS, UMR 7550, 11 rue de l'Université, 67000 Strasbourg, France

⁵ CRAL-Observatoire de Lyon, Université de Lyon, Lyon I, CNRS, UMR 5574, Lyon, France

⁶ New York University Abu Dhabi, PO Box 129188, Abu Dhabi, UAE

⁷ Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany

⁸ Instituto de Astrofísica de Canarias, Vía Láctea s/n, La Laguna, Tenerife, Spain

⁹ Departamento de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife, Spain

¹⁰ Departamento de Física de la Tierra y Astrofísica, UCM, 28040 Madrid, Spain

¹<fnref n="1" /> Universidade de São Paulo, Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Rua do Matão 1226, 05508-090 São Paulo, Brazil

Received 1 October 2019 / Accepted 29 December 2019

## ABSTRACT

We present the second data release (DR2) of the X-shooter Spectral Library (XSL), which contains all the spectra obtained over the six semesters of that program. This release supersedes our first data release from Chen et al. (2014, A&A, 565, A117), with a larger number of spectra (813 observations of 666 stars) and with a more extended wavelength coverage as the data from the near-infrared arm of the X-shooter spectrograph are now included. The DR2 spectra then consist of three segments that were observed simultaneously and, if combined, cover the range between ~300 nm and ~2.45 µm at a spectral resolving power close to R = 10 000. The spectra were corrected for instrument transmission and telluric absorption, and they were also corrected for wavelength-dependent flux-losses in 85% of the cases. On average, synthesized broad-band colors agree with those of the MILES library and of the combined IRTF and Extended IRTF libraries to within ~1%. The scatter in these comparisons indicates typical errors on individual colors in the XSL of 2–4%. The comparison with 2MASS point source photometry shows systematics of up to 5% in some colors, which we attribute mostly to zero-point or transmission curve errors and a scatter that is consistent with the above uncertainty estimates. The final spectra were corrected for radial velocity and are provided in the rest-frame (with wavelengths in air). The spectra cover a large range of spectral types and chemical compositions (with an emphasis on the red giant branch), which makes this library an asset when creating stellar population synthesis models or for the validation of near-ultraviolet to near-infrared theoretical stellar spectra across the Hertzsprung-Russell diagram.

**Key words.** Hertzsprung-Russell and C-M diagrams – catalogs

## 1. Introduction

Stellar spectral libraries are fundamental resources that shape our understanding of stellar astrophysics and allow us to study the stellar populations of galaxies across the Universe. These libraries come in two flavors: empirical, which are composed of a well-defined set of stars with certain stellar atmospheric parameters coverage, and theoretical where stellar spectra are computed for an arbitrarily large set of parameters and extensive wavelength coverage. The list of empirical and theoretical spectral libraries grows continuously due to their versatile usage in modern astrophysics and they are boosted by developments of

new instrumentation<fnref n="1" /> (e.g., Coelho 2009; Husser et al. 2013; Allende Prieto et al. 2018 for the theoretical side and Vazdekis et al. 2016; Villaume et al. 2017 for the empirical side).

For the purpose of building stellar population synthesis models, a good stellar library should have four key properties, as highlighted by Trager (2012): an exhaustive coverage of all stellar evolutionary phases and chemical compositions to represent as well as possible the integrated light of real stellar systems; a broad wavelength coverage since not all stellar phases contribute equally at all wavelengths; simultaneous observations at all wavelengths to avoid issues due to temporal stellar spectral variations; and good calibration of the individual stars in terms of flux and wavelength calibration as well as high-precision stellar atmospheric parameters. This latter point requires comparison with synthetic stellar spectra, which highlights another

<footnote>* Table C.1 and the spectra are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/634/A133</footnote>

<footnote>** Based on observations collected at ESO Paranal La Silla Observatory, Chile, Prog. IDs 084.B-0869, 085.B-0751, and 189.B-0925 (PI Trager).</footnote>

<footnote n="1">¹ For an extensive list of stellar spectral libraries see David Montes' web collection at https://webs.ucm.es/info/Astrof/invest/actividad/spectra.html and references therein.</footnote>
important application of libraries with the above properties: the validation and improvement of theoretical stellar models across the Hertzsprung-Russell (HR) diagram and across wavelengths.

These were the goals of the X-shooter Spectral Library (XSL), which is a moderate-resolution (R ~ 10 000) spectral library that was designed to cover most of the HR diagram. The observations were carried out with the X-shooter three-arm spectrograph on ESO's VLT (Vernet et al. 2011) in two phases, a pilot and an ESO Large Program, spanning six semesters in total. In our first data release (hereafter DR1, Chen et al. 2014), we present spectra of 237 unique stars, which were observed during the pilot program, for a wavelength range that was restricted to the two optical arms of X-shooter (300–1024 nm).

In the current data release (DR2), we present our full set of 813 observations of 666 stars, now also including data from the near-infrared (NIR) arm of X-shooter. This NIR extension is undoubtedly one of the main advantages of XSL over other empirical spectral libraries in the literature. Empirical NIR libraries for a relatively wide range of stellar parameters have been constructed in the past, first at very low spectral resolution (e.g., Johnson & Méndez 1970; Lançon & Rocca-Volmerange 1992), then progressively at intermediate resolution (e.g., Lançon & Wood 2000; Ivanov et al. 2004; Rayner et al. 2009; Villaume et al. 2017), or at higher spectral resolution but in restricted wavelength ranges (e.g., Cenarro et al. 2001; Majewski et al. 2017). We note that some of these libraries did not attempt to preserve the shapes of the stellar continua across the wavelengths observed. Only some have been combined with optical libraries for the purpose of calculating the spectra of synthetic stellar populations (e.g., Pickles 1998; Lançon et al. 1999, 2008; Vazdekis et al. 2003; Maraston 2005; Röck et al. 2016; Conroy et al. 2018). In these efforts, the need to merge optical and NIR observations of distinct samples of stars is an inevitable cause of systematic errors. The X-shooter instrument helps alleviate this issue, as it acquires simultaneous observations of the whole spectral range for every stellar target. In addition, the spectral resolution achieved is higher than those of any of the previous libraries that cover as broad spectral range.

The XSL DR2 data are homogeneously reduced and calibrated, and the spectra are made available in three spectral ranges, corresponding to the three arms of the spectrograph: UVB-300–556 nm, VIS-533–1020 nm, and NIR-994–2480 nm. Our data release contains, similar to DR1, repeated observations of several cool giants stars. Both data releases are available on our website<fnref n="2" />.

Our sample selection and observing strategy are described in Sect. 2. Section 3 gives details about the data reduction and calibration process. The final spectra are assessed in terms of spectral resolution and energy distribution in Sect. 4. In Sect. 5, we summarize the data products made available, before concluding in Sect. 6. Additional details about the input catalogs used to build this library can be found in Appendix A. Appendix B collects peculiarities of some of the stars in the program, and observational artifacts that affect some of the spectra. Finally, Appendix C provides the log of all the observations.

## 2. Sample selection and observing strategy

### 2.1. Selection criteria

The XSL target stars were selected to cover as much of the Hertzsprung-Russell diagram as possible in the allocated time,

<figure><img src="image_1.png" /><figcaption>Fig. 1. Positions of XSL stars in the sky (Aitoff projection).</figcaption></figure>

with a wide range of chemical compositions. Our original sample consists of 679 unique stars. The primary references for the construction of the XSL target list were existing spectral libraries or compilations of stellar parameter measurements. In particular, the XSL sample has a strong overlap with the MILES spectral library (142 stars, Sánchez-Blázquez et al. 2006; Cenarro et al. 2007) and the NGSL library (135 stars, Gregg et al. 2006). We completed the list with objects from the parameter compilation PASTEL (Soubiran et al. 2010, 2016) and from a variety of more specialised catalogs (see Table A.1 for details and corresponding references). In Fig. 1 we plot the distribution of the selected XSL stars on the sky.

More than half of the stars in the XSL library are red giants in the broad sense, which includes red supergiants or asymptotic giant branch stars. These stars provide strong (age-dependent) contributions to the NIR emission of galaxies (e.g., Lançon et al. 1999; Maraston 2005; Melbourne et al. 2012). The luminous red stars in XSL are located in star clusters, in the field, in the Galactic bulge and in the Magellanic Clouds (Table A.1). A fraction of the asymptotic giant branch stars are carbon stars, whose spectra have been studied in detail by Gonneau et al. (2016, 2017). We obtained repeated observations for a fraction of the XSL stars (~20% of our sample), mainly in the cases of luminous cool stars, known or suspected to vary in time.

Figure 2 presents the stellar atmospheric parameters coverage of the XSL sample as determined by Arentsen et al. (2019) and Gonneau et al. (2017). We refer the reader to Fig. 2 of Chen et al. (2014) where the HR diagram is presented using literature values. Arentsen et al. (2019) estimated atmospheric parameters in a homogeneous manner for 754 XSL spectra of 616 stars: they fit the ultraviolet and visible spectra with the ULySS package (Koleva et al. 2009) using the MILES spectral interpolator (Sharma et al. 2016) as a reference. Gonneau et al. (2017) determined atmospheric parameters of carbon-rich stars by comparing the observations to a set of high-resolution synthetic carbon star spectra, based on hydrostatic model atmospheres (Aringer et al. 2016)<fnref n="3" />.

The aim of the selection was to obtain representative spectra mapping the widest area of the parameter space defined by the effective temperature ($T_{eff}$), the surface gravity (log g), and the iron metallicity ([Fe/H]). Therefore, the sample does not reproduce the natural frequencies of various stellar types. At the time the target lists were prepared, large public catalogs of homogeneous stellar parameter estimates were still lacking, and we were unable to populate certain regions of parameters in a satisfactory

<footnote n="2">² http://xsl.astro.unistra.fr</footnote>

<footnote n="3">³ Only the stars less affected by pulsation (namely stars from Groups A, B and C as defined in Gonneau et al. 2017) are plotted in Fig. 2.</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 2.** Stellar atmospheric parameters, calculated by Arentsen et al. (2019, open symbols) and Gonneau et al. (2017, star symbols), for 769 XSL spectra.</figcaption></figure>

way. In particular, some of the Bulge stars selected on the basis of high estimated metallicities turned out to be less metal-rich than expected (Arentsen et al. 2019).

For applications requiring a higher level of completeness, we note that the X-shooter instrument archive has become a useful source of complementary data over the years. Low mass dwarfs from the archive were combined with XSL data for population synthesis applications by Dries (in prep.). The red supergiants of Davies et al. (2013) were observed with settings similar to XSL and can usefully complement the XSL library in that region. Finally, the ESO archive<fnref n="4" /> contains a very large number of observations of B-type stars, because these have been used over many years for the correction of telluric absorption. We do not discuss these external data in this article.

## 2.2. Observing strategy

We gathered the XSL observations in two phases: a pilot program (periods P84 and P85, October 2009–September 2010) and a Large Program over four semesters (P89-P92, April 2012–March 2014).

We observed each science target with two slit widths. We first took a spectrum with a narrow slit to achieve the required spectral resolution. The narrow-slit widths for the UVB, VIS and NIR arms of X-shooter were 0.5", 0.7", 0.6", respectively, leading to nominal resolutions $R = \lambda/\Delta\lambda \sim 9200$, ~11 000 and ~7770. With the narrow slits, a wavelength-dependent fraction of the stellar flux is lost. The loss-fraction depends on the seeing, but also on the centering of the star in the slit and on the performance of the atmospheric dispersion corrector of the instrument, both of which are not controlled precisely. The most reliable approach to correct these losses is to observe the targets with a wide slit positioned vertically on the sky. That observation can be calibrated into absolute flux when conditions are photometric and then can be used to correct the energy distribution of the narrow-slit observation of the same target. In non-photometric conditions, as tolerated in our observations ("thin cirrus" tolerance of the ESO VLT scheduling), absolute fluxes cannot be garanteed, but the shape of the energy distribution is properly recovered with this method under the assumption that the cirrus have a flat

transmission curve<fnref n="5" />. We took wide-slit spectra for every target immediately after the narrow-slit ones, using a 5" wide slit.

We performed our observations using the SLIT spectroscopy mode of X-shooter. In this mode three observing strategies are available: STARE, NODDING, and OFFSET. We used each of them depending on the desired outcome. In STARE mode the star was located at the center of the slit and the sky background was estimated from each side of the stellar signal on the observed frame. We used this mode for the wide-slit observations. The NODDING acquisition mode allows observations of the star at two positions (A and B) along the spectrograph slit. In this way we achieved an improved sky subtraction with a double-pass subtraction. Almost all our narrow-slit observations were taken in NODDING mode. An OFFSET acquisition mode alternates between the star and an empty sky region. We used this mode for 12 observations of very bright stars in the pilot program.

## 2.3. Spectro-photometric standard stars

In order to flux-calibrate our science spectra we used spectro-photometric standard stars (typically white dwarf stars of type DA). These stars were observed once per night using a 5" slit width.

In practice, we used only the five following stars, for which the spectral energy distribution and other parameters are given in the X-shooter manual<fnref n="6" />: BD+174708, EG274, Feige 110, LTT3218, and LTT7987. The other reference stars (CD-32 9927, EG21, GD 153, GD 71, LTT 1020 and LTT 4364) were discarded either because their spectral energy distribution was not reliably known or because the long exposure times of their default observations (300 or 600 seconds) resulted in saturated NIR K-band spectra.

# 3. Data reduction and calibration

We have reduced the full set of XSL observations over the pilot and the Large Programs in a uniform manner. This allows us to provide an updated calibration of the spectra initially released under DR1 and to ensure a homogeneity of the present DR2.

Our data reduction and calibration processes are illustrated in Fig. 3. We apply the first four steps (referred to as “common data reduction”) to all observations: science target stars and flux standards, wide and narrow slits, and all spectral arms. We use the X-shooter data reduction pipeline (XSH pipeline, Sect. 3.1, Modigliani et al. 2010), up to the creation of two-dimensional (2D) order spectra. Then we extract 1-dimensional (1D) spectra outside of the XSH pipeline to better control the rejection of bad pixels (see Sect. 3.2). We further apply corrections for the telluric absorption bands (Sect. 3.3) and for the continuous atmospheric extinction (Sect. 3.4).

After processing all observations in this uniform way, we derive response curves from the flux-standard spectra (Sect. 3.5), and we use these to process both the narrow and wide slit science spectra (Sect. 3.5). Finally, using the wide-slit spectrum corresponding to each narrow slit observation, we derive a wavelength dependent flux-loss correction and apply it to

<footnote n="4">⁴ http://archive.eso.org/wdb/wdb/eso/xshooter/form</footnote>

<footnote n="5">⁵ Using the extinction curve for normal cirrus of Lynch & Mazuk (2001), we find that the ratio of V-band to K-band transmission varies by less than 1% when transmission drops from 100% to 60% because of standard cirrus.</footnote>

<footnote n="6">⁶ http://www.eso.org/sci/facilities/paranal/instruments/xshooter/tools/specphot_list.html</footnote>
<figure><img src="image_3.png" /><figcaption>Fig. 3. Overview of the XSL data reduction process.</figcaption></figure>

achieve the final XSL spectrum for each stellar observation
(Sect. 3.6).

3.1. XSH pipeline

The XSH pipeline consists of a set of data processing mod-
ules that perform the individual tasks of the data reduction.
These “recipes” can be chained in a flexible way to fulfill the
requirement of a specific program. We use the XSH pipeline ver-
sion 2.6.8, in physical mode as per the pipeline recommendation.
In that mode, the wavelength and spatial scale calibrations are
performed by optimising a physical model of the instrument. To
achieve a more flexible automation of our process we run the
XSH pipeline via the Esorex environment<fnref n="7" /> (version 3.12).

For each science observation, the input to the XSH pipeline
is a set of raw data-frames that we identify and collect using the
associations of calibration files provided by the ESO archive.
When necessary, we modify the default associations to ensure
that the flux-standard frames are always reduced with the same
flat-field as the corresponding science frames, rather than with
the flat-field nearest-in-time (see Sect. 3.5). Saturation was a
known problem for several observations of the pilot program but
occurred rarely in later semesters. Here, we apply a scheme sim-
ilar to that of DR1 to flag strongly saturated observations auto-
matically and exclude these from further processing.

We used the default parameters for the creation of the
calibration frames (up to the *xsh_flexcomp* recipe). To trans-
form the science and flux-standard frames into flat-fielded,

rectified and wavelength-calibrated 2D order spectra, we used
recipes *xsh_scired_slit_stare* for wide-slit observations, and
*xsh_scired_slit_offset* for narrow-slit observations. In particular,
we reduced any frames observed in NODDING mode as if they
were taken in OFFSET mode, which gives us more flexibility
in the handling of bad pixels and improves the quality of our
final spectra. In STARE mode, we switch off the sky-subtraction
offered by the XSH pipeline because we apply our own algo-
rithm during the subsequent 1D-extraction step. Two key para-
meters of the above recipes control the sampling of the output
spectra along the wavelength axis and along the spatial axis. We
set them respectively to 0.015 nm and 0.16" for the UVB and
VIS arms, and to 0.06 nm and 0.21" for the NIR arm.

The final products of the reduction with the X-shooter
pipeline are 2D frames containing the spectral orders of each
arm. They include a map of pixel quality flags and a map of the
estimated variances of the pixel errors, which we use to propa-
gate errors through subsequent processing steps.

3.2. 1D extraction

A main driver for our choice of performing the extraction of
1D spectra outside of the XSH pipeline was the need for more
control over the rejection of bad pixels and the sky subtraction.
Our extraction procedure follows a prescription adapted from
Horne (1986). It implements a rejection of masked and out-
lier pixels, as well as a weighting scheme based on a smooth
throughput profile, that can follow residual distortions in the rec-
tified 2D-spectra, and on the local pixel-variance as estimated
by the earlier steps in the ESO pipeline. The extraction pro-
cedure is an updated version of the one presented in Gonneau
et al. (2016), and the description below focuses on the modified
elements.

3.2.1. Fraction of bad pixels

Before extracting a given order, we double-checked for satura-
tion and other invalid pixels by counting the percentage of pixels
flagged as bad along the spectrum. The count is restricted to a
box extending ±3/±2 pixels (respectively for the UVB/VIS arms
and the NIR arm) to either side of the peak of the spatial pro-
file of the spectrum. As a general rule, the extraction process is
abandoned if more than 30% of the pixels are determined to be
bad. Exceptions are implemented for pixels with XSH pipeline-
code 6 (unremoved cosmic ray), a warning that we discard to
improve the extraction of the observations taken in exception-
ally good seeing (see Sect. 3.2.5), and for XSH pipeline-code 21
(extrapolated flux in the NIR) in the longest wavelength order of
the NIR arm, because keeping those pixels was found to reduce
discontinuities between the two last orders of that arm signif-
icantly (Gonneau et al. 2016, this is an update with respect to
who describe the discontinuities).

3.2.2. Sky background removal

We carried out a sky subtraction for all spectra acquired in
STARE mode, as well as for any long-exposure VIS and NIR
frames observed in NODDING and OFFSET modes (science
frames with exposure times ≥600s in the VIS and exposure
time ≥200s in the NIR; flux standard frames with exposure time
≥190s in the VIS and all the NIR frames). The areas along the
slit from which the median sky levels are computed are selected
in each frame after fitting a Gaussian to the spatial profile of the
star, in order to account for the seeing. In STARE mode, these

<footnote n="7">⁷ https://www.eso.org/sci/software/cpl/esorex.html</footnote>
adjustments avoid artifacts such as the negative fluxes that sometimes result from over-subtraction in the standard XSH-pipeline process. In long-exposure spectra, the main purpose is to account for variations in the sky background between the two combined acquisitions.

### 3.2.3. Updated standard extraction

Horne's extraction scheme relies on obtaining first a rough guess for the 1D spectrum (that we refer to as the *standard* spectrum) and then optimizing this extraction with an inverse-variance weighted algorithm, that accounts for the spatial profile of the light at a given wavelength (we call the result the *final* spectrum).

To obtain the standard spectrum, we summed the (sky-subtracted) 2D spectrum along the spatial dimension over a limited range of pixels. We updated our procedure by setting the limits to ±4 or 8 pixels, respectively for NODDING and STARE/OFFSET frames, centered on the peak of the signal. The calculation of the final spectrum is as in Gonneau et al. (2016). Along the spatial axis, the ratio of the flux in a pixel to the normalized spatial profile is taken as an estimate of the stellar flux at the wavelength of interest and these estimates are averaged with an inverse-variance weight. A limited number of iterations are implemented to add outliers to the bad-pixel mask (maximum 3 at any wavelengths), and these pixels are discarded in the final sum.

### 3.2.4. Merging of the individual orders

Figure 4 shows the overlapping orders before merging for each arm of the spectrograph. The flat-fielding algorithm of the XSH pipeline normalizes flat-field frames only globally, hence the blaze function of the instrument is removed in the division by the flat-field, and the extracted orders align well. To merge them, we combined the spectra in the regions of overlap, using weights that combine the inverse variance and a coefficient, $\alpha$, set to vary linearly from 0 to 1 across those regions. This choice was made because we found that the single-order variance was frequently underestimated at the very end of orders. Hence the flux of each overlapping order, $F_{merging}$, is calculated as follows:

$$F_{merging} = \frac{(1 - \alpha) * F_L/V_L + \alpha * F_R/V_R}{((1 - \alpha)/V_L + \alpha/V_R)}, \quad (1)$$

where $F$ and $V$ stand for pixel fluxes and variances, and subscripts $L$ and $R$ refer to the spectral orders on the left and right sides of region of overlap.

For the frames observed in NODDING mode, the 1D extraction procedure produces two separate spectra that we averaged together. The end product of this procedure is a continuous 1D spectrum with two extensions: a flux spectrum and an error spectrum.

### 3.2.5. Chessboard effect

Some of our NIR frames are affected by an issue that we call the “chessboard effect”. When the seeing is exceptional (typically <2 pixels per FWHM of the spectral profile in parts of the raw data), the interpolation scheme used by the XSH pipeline for the geometric transformation into a rectified, wavelength calibrated frame does not perform well. The resulting images (known as ORDER2D frames in the pipeline) display a chessboard pattern as illustrated in the top panel of Fig. 5. Because the wavy patterns seen along single lines of the rectified 2D-spectrum vary

<figure><img src="image_4.png" /><figcaption>Fig. 4. Examples of order-by-order spectra for the three arms (XSL observation X0365). From top to bottom: UVB (12 orders), VIS (15 orders) and NIR (16 orders).</figcaption></figure>

abruptly between one line and the next, the crests of these waves are mistaken for cosmic ray hits.

To deal with this, we first identify pixels with bit code 6 in the quality-mask images produced by the XSH pipeline (flag 6 = unremoved cosmic rays, cf. the bad pixel code conventions in the X-shooter Pipeline Manual) and reset these flags to “good pixel”<fnref n="8" />. Then, to avoid any re-flagging by our own extraction procedure, we define a “safe zone” in the 2D frame in which any cosmic-ray identification is switched off. This rectangular zone extends to ±4 or 6 pixels on either side of the peak of the spatial profile of the spectrum (respectively for NODDING or for STARE/OFFSET frames). This scheme was also used in a few VIS frames. The bottom panels of Fig. 5 show an example of a 2D quality-mask before and after editing, and the extent of the safe zone (represented as a blue rectangle). Figure 6 shows a comparison of the different extraction methods, with and without the flagging and the safe zone, for the ORDER 16 of the NIR arm.

### 3.3. Telluric absorption correction

We corrected the 1D extracted spectra for telluric absorption with `molecfit` (Smette et al. 2015; Kausch et al. 2015). `Molecfit` uses a radiative transfer code and a molecular

<footnote n="8">⁸ A similar solution is now recommended by ESO on the X-shooter FAQ webpage: http://www.eso.org/sci/data-processing/faq/my-spectrum-of-a-bright-source-looks-much-noisier-than-i-expected-why.html</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 5.** Example of ORDER2D frame affected by the chessboard effect. *Top*: input 2D flux spectrum, before extraction. *Middle*: input 2D quality mask. *Bottom*: optimal 2D mask. For both 2D masks, the bad pixels are shown in black and gray.</figcaption></figure>

line database to calculate a synthetic spectrum of the Earth's atmosphere based on local weather conditions and standard atmospheric profiles. In addition, *molecfit* also models the instrumental line spread function. Once the characteristics of the instrument have been taken into account, the synthetic transmission spectrum can be used to correct an observed spectrum for telluric absorption.

Table 1 details the relevant parameters that we used for running *molecfit*. Most of these are taken from Kausch et al. (2015), who also applied *molecfit* to the processing of X-shooter spectra. In the UVB arm we only corrected for absorption by $O_3$, in the VIS arm we corrected for absorption by $H_2O$, $O_2$ and $O_3$, and in the NIR arm for $H_2O$, $CO_2$, $CO$, $CH_4$ and $O_2$. Except for water vapor, the abundances of the different molecular species are assumed to be constant and equal to the standard atmospheric values. The total amount of water vapor along the pointing direction of the observation is a free parameter of *molecfit* that can be fit together with the other free parameters. The wavelength regions used for the fitting procedure are specified in Table 2. For the NIR arm and the VIS arm, these regions are almost the same as those used in Kausch et al. (2015).

It is known that the wavelength calibration of X-shooter is not perfect: Moehler (2015) reports shifts with an average peak-to-peak amplitude of 0.28 pixels, corresponding to 0.0042 nm and 0.0168 nm for our VIS and NIR data, respectively. This results in small offsets between the wavelength scale of the spectrum and the synthetic telluric transmission spectrum

<figure><img src="image_6.png" /><figcaption>**Fig. 6.** Examples of extracted spectra with the chessboard effect (X0648). The red spectra correspond to the inverse-variance weighted extractions, while the black spectra are for the standard extractions. *Upper panel*: original extraction, and *lower panel*: updated extraction to deal with the chessboard effect (flagging and safe zone in the rejection scheme).</figcaption></figure>

<table><caption>**Table 1.** Parameters of molecfit used to correct XSL spectra for telluric absorption.</caption><thead><tr><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td>FTOL</td><td>0.01</td></tr><tr><td>XTOL</td><td>0.01</td></tr><tr><td>LIST_MOLEC</td><td>O3 (UVB)<br>H2O O2 O3 (VIS)<br>H2O CO2 CO CH4 O2 (NIR)</td></tr><tr><td>FIT_MOLEC</td><td>0 (UVB)<br>1 0 0 (VIS)<br>1 0 0 0 0 (NIR)</td></tr><tr><td>RELCOL</td><td>1.0 (UVB)<br>1.0 1.0 1.0 (VIS)<br>1.0 1.0 1.0 1.0 1.0 (NIR)</td></tr><tr><td>FIT_CONT</td><td>1</td></tr><tr><td>CONT_N</td><td>3</td></tr><tr><td>CONT_CONST</td><td>Average of spectrum</td></tr><tr><td>FIT_WLC</td><td>1</td></tr><tr><td>WLC_N</td><td>0</td></tr><tr><td>WLC_CONST</td><td>0.0</td></tr><tr><td>FIT_REShti</td><td>0</td></tr><tr><td>KERNMODE</td><td>0</td></tr><tr><td>FIT_RES_GAUSS</td><td>1</td></tr><tr><td>RES_GAUSS</td><td>1.0</td></tr><tr><td>VARKERN</td><td>1</td></tr></tbody></table>

**Notes.** See Noll et al. (2014) for an explanation of the meaning of the different parameters.

produced by molecfit. We alleviate this problem by applying molecfit in a slightly different way. First, we use the “classical” molecfit approach, as described in Smette et al. (2015), to derive the precipitable water vapor column (PWV). Then we cut the spectrum into different wavelength segments. The segments are chosen such that each one contains at least one telluric feature and, if possible, also a “clean” region which is (almost) free of telluric contamination. We fix the PWV parameter to the value found in the first iteration and apply molecfit to each of the wavelength segments. As a final step, the corrected wavelength segments are recombined. Across regions where wavelength
<table><caption>**Table 2.** Wavelength regions (air) used by *molecfit* for fitting atmospheric transmission spectrum.</caption><thead><tr><th>Arm</th><th>Wavelength region [nm]</th><th>Molecule</th></tr></thead><tbody><tr><td>UVB</td><td>302–350</td><td>O<sub>3</sub></td></tr><tr><td>VIS</td><td>685–694</td><td>O<sub>2</sub></td></tr><tr><td>VIS</td><td>758–777</td><td>O<sub>2</sub></td></tr><tr><td>VIS</td><td>929–945</td><td>H<sub>2</sub>O</td></tr><tr><td>NIR</td><td>1120–1130</td><td>H<sub>2</sub>O</td></tr><tr><td>NIR</td><td>1470–1480</td><td>H<sub>2</sub>O</td></tr><tr><td>NIR</td><td>1800–1810</td><td>H<sub>2</sub>O</td></tr><tr><td>NIR</td><td>2060–2070</td><td>CO<sub>2</sub></td></tr><tr><td>NIR</td><td>2350–2360</td><td>CH<sub>4</sub></td></tr></tbody></table>

**Notes.** For each region, the most important molecule of the transmission spectrum is also given.

segments overlap, we use linearly progressive weights to average the two sets of data.

Our “extended molecfit” approach, where the spectrum of a particular arm is divided into smaller wavelength segments, allows molecfit to find better local wavelength solutions and line spread functions. The final telluric-corrected spectra have a much smaller variance. We use the extended molecfit approach to correct spectra in the VIS and NIR arms. The boundaries of the different wavelength segments are given in Table 3. Figure 7 compares the two telluric correction approaches: the classical molecfit approach (in gray) and the extended molecfit approach (in red). Figure 8 illustrates the results for the full VIS and NIR wavelength ranges. For the UVB arm, we use the classical molecfit approach since for this arm the wavelength region that contains telluric lines is relatively small.

### 3.4. Atmospheric extinction correction

Before computing instrument response curves and applying them, the spectra of science and flux-standard stars were corrected for observation-specific effects. We corrected the frames for exposure time, gain and continuous atmospheric extinction as follows:

$$F_{\text{out}} = \frac{F_{\text{in}}}{\text{exptime} \times \text{gain}} \times 10^{(0.4 \cdot \text{airmass-extinct})}, \quad (2)$$

where $F_{\text{in}}$ is the flux of the observations as available after correction for telluric absorption bands. The airmass, exposure time (exptime) and gain were taken from the observation header.

The atmospheric extinction (extinct, expressed in magnitudes per unit airmass) combines aerosol scattering (Mie diffusion, $k_{\text{aero}}$) and Rayleigh scattering ($k_{\text{ray}}$), following Patat et al. (2011). Molecular absorption (including broad features that are sometimes dealt with like a contribution to the extinction continuum) were already accounted for by molecfit. We model the aerosol contribution as advocated by Moehler et al. (2014):

$$k_{\text{aero}} = 0.014 \times (\lambda[\mu\text{m}]^{-1.38}). \quad (3)$$

The Rayleigh scattering mostly depends on the observation airmass, and only weakly on the atmospheric conditions. In the zenith direction, we model it as:

$$k_{\text{ray}} = \frac{p}{1013.25} \times (0.00864 + 6.5 \times 10^{-6} H) \times \lambda^{-(3.916+0.074\lambda+0.050/\lambda)}, \quad (4)$$

<table><caption>**Table 3.** Wavelength segments (air) used for the extended molecfit approach, as described in the text.</caption><thead><tr><th>Arm</th><th>Wavelength segment [nm]</th><th>Arm</th><th>Wavelength segment [nm]</th></tr></thead><tbody><tr><td>VIS</td><td>λ<sub>start</sub>–637</td><td>NIR</td><td>λ<sub>start</sub>–1105</td></tr><tr><td>VIS</td><td>633–712</td><td>NIR</td><td>1095–1185</td></tr><tr><td>VIS</td><td>708–752</td><td>NIR</td><td>1175–1255</td></tr><tr><td>VIS</td><td>748–782</td><td>NIR</td><td>1245–1305</td></tr><tr><td>VIS</td><td>778–862</td><td>NIR</td><td>1295–1455</td></tr><tr><td>VIS</td><td>858–926</td><td>NIR</td><td>1445–1605</td></tr><tr><td rowspan="7">VIS</td><td rowspan="7">922–λ<sub>end</sub></td><td>NIR</td><td>1595–1755</td></tr><tr><td>NIR</td><td>1745–1880</td></tr><tr><td>NIR</td><td>1870–1985</td></tr><tr><td>NIR</td><td>1975–2040</td></tr><tr><td>NIR</td><td>2030–2085</td></tr><tr><td>NIR</td><td>2075–2285</td></tr><tr><td>NIR</td><td>2275–2385</td></tr><tr><td>NIR</td><td>2375–λ<sub>end</sub></td></tr></tbody></table>

**Notes.** The extended molecfit approach is only used for the VIS arm and the NIR arm.

with pressure $p = 744$ hPa, height $H = 2.64$ km and wavelength $\lambda$ in $\mu\text{m}$ (values from Noll et al. 2012).

Figure 9 compares our extinction curve to the XSH pipeline one. From the bottom panel of that figure, we can see that the main differences are in the ozone Huggins bands, which are not deep enough in the reference curve.

### 3.5. Response curve and first flux calibration

In the XSH pipeline the 2D flat-field spectra are not normalized across the surface of each grating order before they are divided into the program data. Hence the division corrects both the pixel-to-pixel variations of the sensitivity of the detector and major transmissions variations, such as those produced by the blaze functions of the gratings. This choice ensures that extracted single-order spectra within one arm are easily connectable. In exchange for this benefit, the arm-spectra (as illustrated in Figs. 4 and 8) contain an imprint of the (inverted) spectra of the lamps used to illuminate the flat-field exposures. Unfortunately, the X-shooter flat-field lamps are less stable than other elements of the acquisition chain, in particular in the NIR arm. Hence, when pairing observations of a science target and of a spectro-photometric standard star, it is recommended to use the same flat-field images for the two data sets<fnref n="9" />.

To determine a response curve from the 1D spectrum of a flux standard, we performed a $\lambda^2$-minimization between that observation on one hand, and the product of its theoretical spectrum and the unknown response curve on the other. The response curve is represented with a high-order spline polynomial. The typical number of spline nodes is 35 for the UVB and NIR arms and 60 for the VIS. These spline nodes are not all regularly spaced; they are placed sparsely where spurious oscillations should be avoided (e.g., across the telluric water bands in the NIR arm) and more tightly where small-scale instrumental effects must be accounted for (e.g., features at 365 nm due to the combination of two different flat-fields in the UVB

<footnote n="9">⁹ See http://www.eso.org/observing/dfo/quality/XSHOOTER/qc/problems/problems_xshooter.html#NIR FF; current recommendations extend to the VIS and UVB arms.</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig. 7.** Comparison between telluric correction with the classical *molecfit* approach (MF TC) and with the extended *molecfit* approach (extended MF TC). The upper star corresponds to a flux-standard star and the lower star to X0073 (HD 18769). The figure shows different telluric regions in the VIS arm and NIR arm. Within each of these regions, the corrected spectrum is normalized to unit mean.</figcaption></figure>

arm<fnref n="10" />, the dichroic features at 560-620 nm<fnref n="11" />, a flat-field bump at 2100 nm; see Moehler et al. 2014 for more details). Furthermore, we masked the Balmer lines in the UVB and VIS arms, the O₂ lines in the VIS and the deepest telluric regions in the NIR.

We used the response curves for a first flux-calibration of the science spectra both for the narrow and wide slit observations. The resulting intensities are physical flux densities, in erg s⁻¹ cm⁻² Å⁻¹. These spectra are cleaned from the atmospheric and instrumental effects, *except* for the slit losses. The narrow-slit spectra, in particular, still do not reflect the absolute distribution of the stellar flux.

### 3.6. Flux-loss correction and final energy distributions

For 85% of the spectra, wide-slit observations of sufficient quality are available, and thus the narrow-slit spectra are corrected for wavelength-dependent flux-losses. To this end, we computed a smoothed version of the ratio of the wide-slit over the narrow-slit spectrum, and then fit a low-order polynomial (typically of order 1 for UVB/VIS frames, and 3 for NIR frames). Finally, we multiplied our narrow-slit spectra by this polynomial. In a few cases (3 observations of 2 a priori non-variable stars, namely X0194 [HD 193281] and X0306/X0311 [HD 29391]), when only one of several observations of the same star had an associated wide-slit spectrum or an otherwise more reliable continuum,

<figure><img src="image_8.png" /><figcaption>**Fig. 8.** Before (black) and after (red) the telluric absorption correction process for X0365 (VIS and NIR arms).</figcaption></figure>

that observation was used as a reference for the other observation(s).

In summary, our data reduction procedure uses the ESO pipeline until the 2D spectra are rectified and wavelength-calibrated. Subsequent steps were adjusted, mainly to improve the usage of pixel-flags, the sky subtraction, the telluric correction

<footnote n="10">¹⁰ https://www.eso.org/observing/dfo/quality/XSHOOTER/pipeline/pipe_problems.html#uvb flux</footnote>

<footnote n="11">¹¹ https://www.eso.org/sci/facilities/paranal/instruments/xshooter/knownpb.html</footnote>
<figure><img src="image_9.png" /><figcaption>**Fig. 9.** Top: UVB flux-standard spectrum (related to X0365) corrected for O<sub>3</sub> with Molecfit (red spectrum), compared to an uncorrected one (black spectrum). Middle: our own transmission curve (red curve) – with the details of its three components – compared to the pipeline one (black curve). Bottom: ratios of the transmissions.</figcaption></figure>

and the combined usage of narrow-slit and broad-slit observation for flux calibration. 85% of the DR2 spectra are flux calibrated in the sense that they are corrected for all stable transmission factors of the acquisition chain, and for wavelength-dependent losses due to the narrow widths of the slits that provide the desired spectral resolution. Residual flux calibration errors may come from changes in telescope+instrument+sky transmission-curves between the observations of the science targets and of the spectrophotometric standard star, and from residual slit losses in the broad-slit observations. They are assessed via comparisons with external data sets in Sect. 4.5. A few examples of final DR2 spectra can be found in Fig. 10.

# 4. Properties of the DR2 spectra

In this section we present an assessment of the overall properties of the spectra, including the quality of the wavelength and flux calibration, and properties of the line-spread function. We describe the quality flags associated with each published spectrum.

## 4.1. Typical signal-to-noise ratios

Figure 11 shows the distribution of the median signal-to-noise ratios (S/Ns) per pixel for our observations. Except for the coolest stars, typical S/N values are about 70, 90 and 96 for the

UVB, VIS and NIR arms respectively. In fact, the S/N for cool stars ($T_{\text{eff}} < 5000$ K) varies drastically with wavelength across the UVB arm. S/N similar to those of warmer stars are found in a temperature-dependent range of wavelength toward the red end of that arm.

Correlations between the errors in neighboring pixels are expected as a result of the redistribution of pixel fluxes when the original 2D-images are transformed into rectified order-images. They depend on the observation mode, on the seeing and on wavelength in a complex fashion. We have not attempted to characterize them in detail. In general however, they would extend over a few pixels, on a wavelength range similar to the (arm-dependent) resolution element of the spectra.

## 4.2. Spectral resolution and line spread function

We have determined the spectral resolution and the line spread function for each arm with the Penalized PiXel-Fitting method (pPXF, Cappellari & Emsellem 2004; Cappellari 2017), using a subset of the PHOENIX theoretical spectral library (Husser et al. 2013) as templates. The available PHOENIX library samples effective temperatures with a step of 200 K (or of 500 K above 12 000 K), surface gravities with a step of 0.5 dex, and [Fe/H] with a step of 0.5 dex (or of 1 dex where [Fe/H] < -1). The subset retained contains the 388 theoretical stars that are closest to the XSL library stars in the stellar parameter space, based on the stellar parameters derived from XSL observations by Arentsen et al. (2019).

All XSL UVB, VIS and NIR spectra were fit one 200 Å-wide wavelength region at a time. In the NIR, the telluric-contaminated regions (13 500–14 200 Å and 17 700–21 000 Å) were not analyzed.

### 4.2.1. The spectral resolution

The spectral resolution is described by the velocity dispersion measurements of the pPXF fit, as the template spectrum has a much higher resolution ($R \sim 500 000$, Husser et al. 2013) than the X-shooter spectra ($R \sim 10 000$). As the observations span multiple observing periods, the measurements for stars from different observing periods were analyzed separately and then combined. In each observing period, the velocity dispersion results from the pPXF fitting were sigma-clipped with a 5-$\sigma$ threshold. Sigma-clipping was iterated until convergence was achieved.

Figure 12 shows the wavelength dependence of the resolution for each observing period and X-shooter arm. The 1 $\sigma$ error bars belong to period P90 measurements, but are similar for all observing periods. The scatter around the median value comes from stars with lower S/N or a low number of spectral lines in each wavelength region. The red points show the weighted averages of the medians of all observing periods. The resolution is stable throughout the observing periods, with only P84 observations showing slightly lower resolution. The resolution $R = \lambda/\Delta\lambda$ is also relatively constant with wavelength, although Chen et al. (2014) initially showed otherwise for the UVB arm. However, they only used the F, G, K stars of observing periods P84 and P85. The adopted spectral resolution for each arm is the average value over all wavelengths: $13 \pm 1 \text{ km s}^{-1}$ for UVB, $11 \pm 1 \text{ km s}^{-1}$ for VIS, and $16 \pm 1 \text{ km s}^{-1}$ for NIR. These values are shown with the black dashed lines and the $\pm 1\sigma$ gray shaded areas in Fig. 12.

The spectral resolution measures are converted from $\sigma$ to FWHM assuming Gaussian line profiles in velocity space ($FWHM = 2.35\sigma$) and used to re-evaluate the resolving power
<figure><img src="image_10.png" /><figcaption>**Fig. 10.** Typical XSL spectra in the OBAFGKMC sequence in log-scale of wavelength. For each star, the corresponding XSL observation number and spectral type are indicated in brackets. The gaps indicate the overlap regions between consecutive arms. Gray bands mask the deepest telluric regions in the NIR arm, as well as a few other local artifacts.</figcaption></figure>

$R = c/FWHM$ (c being the speed of light). Our adopted values are close to the nominal ones<fnref n="12" />, as can be seen in Table 4.

### 4.2.2. The shape of the line spread function

The deviations of the line spread function (LSF) from a Gaussian are described with Gauss-Hermite velocity moments $h_3$ and $h_4$, which are related to skewness and kurtosis respectively. These moments can be obtained with pXF only when the

<footnote n="12">¹² https://www.eso.org/sci/facilities/paranal/instruments/xshooter/inst.html</footnote>
<figure><img src="image_11.png" /><figcaption>Fig. 11. Distribution of the signal-to-noise ratio per arm.</figcaption></figure>

signal-to-noise ratio is excellent. Therefore, $h_3$ and $h_4$ measurements from all observing periods were analyzed together. We computed the S/N in each wavelength bin (window of 200 Å each) and rejected measurements from spectra with $S/N < 200$. Additionally, measurements from the very blue wavelengths, $\lambda < 3600$ Å, were not included due to low signal. The Gauss-Hermite moments were finally computed from 6510 highest-quality pXF measurements. Table 5 shows the median and the standard deviation of the results per arm. To within one standard deviation, the line spread function in all arms is a Gaussian.

### 4.3. Rest-frame velocity correction

The DR2 data are corrected for radial velocity and provided in the rest-frame (with wavelengths in air). The barycentric radial velocities, $cz$, were estimated separately in each of the three arms, because differences between arms are common and often in excess of $10 \text{ km s}^{-1}$ (see Moehler 2015, for a discussion). The $cz$ values are available in the header of the final FITS files, together with the barycentric corrections we applied$^{13}$.

In the UVB and VIS arms, the final estimates were obtained with the ULySS package (Koleva et al. 2009) and the ELODIE 3.2 spectral interpolator (Wu et al. 2011), as also used in the study of stellar parameters of Arentsen et al. (2019). This full-spectrum fitting approach requires an approximate a priori knowledge of the radial velocity for the method to converge, because it performs a local minimization using a Levenberg-Marquardt algorithm. The accuracy of the $cz$ guess should be better than a few times the FWHM of the spectral lines. At the

resolution of XSL, the FWHM of the lines is on the order of 30 km s$^{-1}$, and since the barycentric velocities of the Galactic and LMC/SMC XSL stars are expected to span a range of up to about 400 km s$^{-1}$, we repeated the fit starting with a series of $cz$ guess in the interval from -300 to +300 km s$^{-1}$ with a step of 100 km s$^{-1}$. For each value of the guess, the atmospheric parameters $T_{\text{eff}}$, log g and [Fe/H] as well as $cz$ and the broadening were re-adjusted. The solution with the lowest $\chi^2$ was adopted.

In the NIR, the adopted values of $cz$ were estimated by cross-correlation with the nearest model in the library of synthetic spectra of Husser et al. (2013). Starting from the estimated parameters of Arentsen et al. (2019), we used $\chi^2$-minimization to identify an adequate NIR representation of each NIR spectrum. The spectral ranges used for the cross-correlation were adjusted based on the effective temperature and spectral type, in order to include a sufficient number of spectral features. In the coolest stars, we avoided regions dominated by a single steep molecular bandhead, as the resulting $cz$ then depend on the way the spectra are renormalized to remove the dominant slope, before the cross-correlation.

For some spectra of cool stars, it was not possible to determine $cz$ in the UVB, or even in the VIS, because of an insufficient flux. When a measurement in the NIR was available, we then simply assigned the latter to the UVB and VIS.

During the preparation of DR2, several methods other than those finally adopted were employed to re-estimate the radial velocities of various data subsets. The distributions of velocity differences between methods typically display an approximately Gaussian peak, with a standard deviation of 2 to 3 km s$^{-1}$; non-Gaussian tails of the distributions typically contain ~13%, ~18% and ~10% of the stars in the UVB, the VIS and the NIR arms. Hence the value of 2 to 3 km s$^{-1}$ is representative of the uncertainties on the rest-frame wavelengths of the DR2 data.

### 4.4. Quality assurance

After the data reduction, we assigned a number of quality flags to our data based on visual inspection and on experience gained by comparing the observations with synthetic or other empirical spectra (e.g., Gonneau et al. 2017; Lançon et al. 2018; Arentsen et al. 2019). Table 8 (item “Quality flags (per arm)”) describes our flags. The flags were added to the headers of the final FITS files, and peculiarities specific to individual spectra are listed in Table B.1. We note that this list is not exhaustive, as some artifacts might have been overlooked.

More generic caveats concerning the data are the following. In the UVB, the spectra tend to be wavy near 0.35 µm (the combination of two different flat-fields D2 and QTH in this range leads to errors, see footnote$^{10}$) and between 0.545–0.55 µm (known instabilities of the transmission of the dichroic plate that feeds the UVB and VIS arms of X-shooter, see footnote$^{11}$).

For the VIS arm, frequent artifacts are seen in the spectral shape between 0.54–0.59 µm (dichroic issue, see footnote$^{11}$). The continuum at the red end, beyond 0.97 µm, can also have some artifacts (second dichroic region; effect >5%). Molecfit corrects only for molecular absorption lines, hence telluric absorption by metals is not removed (an example is the frequent contamination of the stellar NaI doublet by telluric Na, between 588 and 592 nm).

In the NIR, a step sometimes arises between the two reddest orders (near 2.25 µm). The noise spectrum does not capture the 1.1 µm telluric band sufficiently strongly. The slope of the spectrum is sometimes poorly determined below 1.05 µm (dichroic region; effect >5%). The K-band is usually of poor quality for

<footnote n="13">13 The barycentric radial velocity is the sum of the topocentric velocity and the barycentric correction.</footnote>
<figure><img src="image_12.png" /><figcaption>**Fig. 12.** Spectral resolution wavelength dependence for each observing periods. Gray, black, green, blue, dark blue and indigo markers show the medians of measurements of data from observing periods "P84", "P85", "P89", "P90", "P91" and "P92" respectively. The red point shows the weighted average of all observing periods. The error bars belong to period P90 measurements and are shown for illustrative purposes.</figcaption></figure>

<table><caption>**Table 4.** Adopted spectral resolution values per arm.</caption><thead><tr><th>Arm</th><th>σ<br>[km s<sup>-1</sup>]</th><th>FWHM<br>[km s<sup>-1</sup>]</th><th>R<br>(our values)</th><th>R<br>(nominal values)</th></tr></thead><tbody><tr><td>UVB</td><td>13</td><td>30.61</td><td>9793</td><td>9700</td></tr><tr><td>VIS</td><td>11</td><td>25.90</td><td>11 573</td><td>11 400</td></tr><tr><td>NIR</td><td>16</td><td>37.68</td><td>7956</td><td>8100</td></tr></tbody></table>

hot stars (K-band continuum wavy, lower S/N). Some of these stars are flagged with HAIR_NIR (meaning that there are some narrow spikes in the NIR spectrum, cf Table 8).

## 4.5. Photometric comparison

Chen et al. (2014) compared the DR1 (UVB and VIS) spectra with literature spectra taken from NGSL (Gregg et al. 2006), the calcium-triplet library CaT (Cenarro et al. 2001), and also with higher-resolution spectral libraries such as UVES POP (Bagnulo et al. 2003) and ELODIE (Prugniel & Soubiran 2001). They found a good agreement in the line shapes and depths between XSL and these libraries.

Here we focus on the broad-band spectral energy distribution across the single-arm spectra of DR2. Based on synthetic photometry, we compare colors within the UVB and VIS arms to those of the MILES library (Sánchez-Blázquez et al. 2006); in the NIR, we confront the XSL colors with those of the 2MASS survey (Cohen et al. 2003; Skrutskie et al. 2006), and with synthetic colors of the IRTF/E-IRTF spectra (Rayner et al. 2009; Villaume et al. 2017).

<table><caption>**Table 5.** Line-spread function parameters $h_3$ and $h_4$.</caption><thead><tr><th></th><th>UVB</th><th>VIS</th><th>NIR</th></tr></thead><tbody><tr><th>$h_3$</th><td>0.00 (0.06)</td><td>-0.03 (0.06)</td><td>0.03 (0.05)</td></tr><tr><th>$h_4$</th><td>-0.03 (0.09)</td><td>-0.01 (0.07)</td><td>-0.01 (0.06)</td></tr></tbody></table>

<table><caption>**Table 6.** Artificial filters used in Sect. 4.5.1.</caption><thead><tr><th>Reference set</th><th>Filter</th><th>λ<sub>min</sub> [Å]</th><th>λ<sub>max</sub> [Å]</th></tr></thead><tbody><tr><td rowspan="4">MILES</td><td>box1</td><td>3900</td><td>4500</td></tr><tr><td>box2</td><td>4600</td><td>5200</td></tr><tr><td>box3</td><td>5400</td><td>6200</td></tr><tr><td>box4</td><td>6300</td><td>7100</td></tr></tbody></table>

### 4.5.1. Comparison with MILES

Many of the XSL stars were initially selected from the MILES collection (Sánchez-Blázquez et al. 2006). Forty MILES stars were present in XSL DR1, and Chen et al. (2014) found a good agreement between these data in the UVB and VIS arms. In the present release, 173 spectra of 142 stars have a MILES counterpart, of which 164 are corrected for slit-losses. To avoid complications due to the edges of the X-shooter arms, we define artificial broad-band filters to compare the synthetic photometry of the two data sets (Table 6 and Fig. 13).

Figure 14 shows how the MILES and XSL synthetic colors compare. We find excellent agreement in general, with a mean offset of -0.005/-0.01 mag and an rms scatter of
<figure><img src="image_13.png" /><figcaption>**Fig. 13.** Filters designed for the comparison with MILES spectra, overlaid on the XSL spectrum of a cool star. The filters avoid the spectral regions most strongly affected by the unstable dichroic between the UVB and VIS arms (plotted in gray).</figcaption></figure>

0.056/0.032 mag, respectively for colors (box1-box2) and (box3-box4). We have examined outliers with color differences larger than 0.1 mag individually, using plots of the spectra and the extensive comparison of XSL spectra with the theoretical spectra of Husser et al. (2013) that is the subject of another article (in prep.). The clear outlier at the red extreme in the bottom panel is BD+19 5116B (XSL observation X0440, spectral type M4.0Ve), a flare star. We observed this star at a different phase than MILES. The star is also an outlier in (box1-box2). Six other outliers at the red end of (box1-box2) are cluster stars. Five of these are in NGC 6838. We trust their XSL energy distributions to better than 10% because the comparison of these spectra with reddened theoretical models raises no alarms and leads to an extinction estimate that agrees with the values listed in Di Cecco et al. (2015) for the cluster (article in preparation). One object is in NGC 1904, and again best-fit models are not alarming and agree with the very small reddening toward that cluster found in the literature. Still in (box1-box2), the bluest outlier is a hot star, which is also represented very well with models. Removing those outliers, the rms scatter in (box1-box2) drops from 0.056 to 0.030 mag. The two outliers in (box3-box4) not yet discussed are the hot star just mentioned, and an F-type star for which the fit with model spectra excludes large error on the slope in the VIS arm. Removing these, the rms scatter for (box3-box4) drops from 0.032 to 0.022 mag. Table 7 summarizes these dispersion measurements.

### 4.5.2. Comparison with 2MASS

The 2MASS filters *J*, *H* and *K*<sub>s</sub> fit entirely into the wavelength range of the NIR arm of X-shooter, as seen in Fig. 15. We compute synthetic photometry compatible with 2MASS using the relative spectral response curves and reference fluxes of Cohen et al. (2003). The 2MASS filters *J*, *H* and *K*<sub>s</sub> fit entirely into the wavelength range of the NIR arm of X-shooter, as seen in Fig. 15).

Figure 16 shows the results of the comparison. Overall there is good agreement. Small systematic offsets indicate small zero point errors in the exercise; they are in opposite directions for (*J* − *H*) and (*H* − *K*<sub>s</sub>), suggesting a normalization difference in the *H* band (as also mentioned in Villaume et al. 2017). We find an rms scatter around the mean of 0.050 mag for (*J* − *H*) and 0.053 mag for (*H* − *K*<sub>s</sub>). After removing the contribution explained by 2MASS errors in quadrature, the residual scatter

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Comparison of synthetic colors between XSL and MILES. Only XSL-spectra with slit-loss corrections are used.</figcaption></figure>

<table><caption>**Table 7.** Standard deviations of the color-differences between XSL and external data sets.</caption><thead><tr><th>External reference</th><th>Color</th><th>Std. dev. σ [mag]</th><th>Comments</th></tr></thead><tbody><tr><td>MILES</td><td>(box1-box2)</td><td>0.030</td><td>0.056 incl. 7 outliers</td></tr><tr><td>MILES</td><td>(box3-box4)</td><td>0.022</td><td>0.032 incl. 3 outliers</td></tr><tr><td>2MASS</td><td>(J - H)</td><td>0.032</td><td>0.050 incl. 2MASS errors</td></tr><tr><td>2MASS</td><td>(H - K<sub>s</sub>)</td><td>0.038</td><td>0.056 incl. 2MASS errors</td></tr><tr><td>IRTF</td><td>(J - H)</td><td>0.023</td><td>0.033 incl. 6 outliers</td></tr><tr><td>IRTF</td><td>(H - K<sub>s</sub>)</td><td>0.019</td><td>0.029 incl. 4 outliers</td></tr></tbody></table>

**Notes.** The outliers are not the same for all colors. The 2MASS standard deviations listed each include six 4-σ outliers. See text for details.

is of 0.032 mag in (J − *H*) and 0.038 mag in (H − *K*<sub>s</sub>). These are still conservative estimates of the errors in the NIR XSL colors, which are clarified in Sect. 4.5.3. In both panels, the larger dispersion toward redder objects is due to cool stars (M and C stars,
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Response curves of the 2MASS filters, overlaid on an XSL spectrum. The transmission curves are taken from Cohen et al. (2003).</figcaption></figure>

and stars with $T_{\text{eff}} < 3500$ K), of which many are known or sus- pected variables.

### 4.5.3. Comparison with IRTF

Empirical NIR spectral libraries are still rare in the literature. We choose to compare XSL with the IRTF spectral library (Rayner et al. 2009) and its later extension (hereafter E-IRTF, Villaume et al. 2017) as their wavelength range overlaps with ours (0.7–2.5 µm) and the spectra are flux-calibrated. Their resolution is $R \sim 2000$, and their S/Ns are on the order of 100. IRTF contains mainly solar-metallicity stars with spectral types between F and M, plus some AGB stars and L dwarfs. E-IRTF extends the metallicity coverage of IRTF for late-type stars from mainly solar to $-1.7 < [\text{Fe/H}] < 0.3$. Our sample counts 19 stars in common with IRTF (in practice 25 XSL spectra) and 48 stars in common with E-IRTF (61 XSL spectra). The cross-match between XSL and IRTF then corresponds to 86 spectra in total, and 77 when we restrict the sample to XSL spectra with slit-loss corrections.

Figure 17 shows the comparison. We find a good agreement with insignificant offsets ($\le 0.007$ mag) and an initial rms scatter of 0.033 and 0.029 mag, respectively for ($J-H$) and ($H-K_s$). We checked all the outliers individually (respectively 6 and 4 cases for these two colors), by overplotting the IRTF and XSL spectra on one hand, and by inspecting the comparison of the XSL spectra with best-fit stellar models on the other. One outlier is a Mira that indeed varied in spectral type between the observations. Another is a star listed as abnormal in Table B.1 (X0878). For all other cases, the IRTF SEDs show signs of a discontinuity between the *J* and the *H* band, while the comparison between the XSL and the theoretical spectra essentially excludes errors larger than a few percent. Removing the 6 suspicious cases from the sample, the scatter drops to 0.023 mag in (*J* − *H*) and 0.019 mag in (*H* − *K*<sub>s</sub>). Hence a representative value of our internal uncertainties in these NIR colors, at least for this subsample, is 2%.

## 5. Final data products and availability

The final data products released as DR2 consist of single-arms spectra, with linearly sampled wavelengths in the rest-frame. The spectra are provided as FITS binary tables, each containing three columns: the wavelength (in nm), the flux spectrum (in erg s⁻¹ cm⁻² Å⁻¹), and an error spectrum (in the same units). A set of keywords in the headers of the FITS files keep track

<figure><img src="image_17.png" /><figcaption>**Fig. 16.** Comparison of synthetic colors between XSL and 2MASS. Only XSL-spectra with slit-loss corrections are used. The brown points correspond to cool stars ($T_{\text{eff}} < 3500$ K or spectral type M and C). The legends provide the average offset between the colors of the two data sets, the rms scatter around the mean, as well as the part of the rms scatter that is explained by the 2MASS error bars shown in the bottom panels.</figcaption></figure>

of various steps of the data-reduction and calibration processes, and of results of the quality assessments. A (non-exhaustive) list of these keywords is given in Table 8.

The DR2 spectra are made available as two separate sets. The first contains those spectra for which a correction for flux-losses across the spectrograph slit was possible, as described in Sect. 3.6). It represents 85% of the observations. We choose to also deliver the remaining 15% of the spectra, for which the flux-losses could not be estimated and corrected. They are identified via header keyword LOSS_COR, which is then set to F. The data can be downloaded from the XSL website, from the CDS interface or from the ESO Phase 3 interface.

## 6. Concluding remarks

In Chen et al. (2014), we release 246 near-ultraviolet and visible spectra for 237 stars. In this paper, we present the second data release of the X-shooter Spectral Library, which consists of
<figure><img src="image_18.png" /><figcaption>**Fig. 17.** Comparison of synthetic colors between XSL and IRTF (open circles for IRTF and filled circles for E-IRTF). The brown symbols correspond to cool stars (M/C stars, and stars with no spectral type and $T_{\text{eff}} < 3500$ K).</figcaption></figure>

all the spectra observed as part of this Large Program, in all the three wavelength ranges of the instrument. The major addition of this release is the NIR data, which extends out to 2.5 µm. The simultaneous optical and NIR observations in this unique data set eliminate a major risk of inconsistencies that has prevailed in previous empirical multi-wavelength libraries based on independent samples in the two wavelength regions. The library contains 2388 spectra of 666 unique stars, at a resolution $R \sim 10\,000$. This corresponds to 791 UVB spectra, 811 VIS spectra and 786 NIR spectra of these stars.

All the spectra, including those previously released, have been reduced in a homogeneous way with updated algorithms. They were corrected for atmospheric extinction and telluric absorption, which maximizes the wavelength ranges suitable for scientific exploitation. In particular, the wings of the water bands of the atmospheres of cool stars in the NIR arm can be measured, and synthetic photometry can be exploited across wavelength regions where it would otherwise be close to meaningless.

The observational setup and the data reduction were tailored to achieve a good relative flux calibration of the spectra, prioritizing the energy distributions over absolute fluxes. This effort was driven by two main scientific applications foreseen for XSL. In population synthesis calculations, the energy distributions are important in determining the relative fluxes of various stars at different wavelengths. They are also relevant if the data are used to test how well current theoretical stellar models are able to reproduce observed colors and spectral features simultaneously. Flux calibration is notoriously difficult to achieve with slit-fed echelle spectrographs and X-shooter is no exception. 85% of the DR2 spectra are flux calibrated in the sense that they are corrected for all stable transmission factors of the acquisition chain, and for wavelength-dependent losses due to the narrow widths of the slits that ensured the desired spectral resolution. For these 85%, residual flux calibration errors may come from changes in telescope+instrument+sky transmission curves between the observations of the science targets and of the spectrophotometric standard star, and from residual slit losses in broad-slit observations. The remaining 15% of the DR2 spectra are not corrected for slit-losss, because the matching broad-slit observations were faulty. Systematic differences of the energy distributions with those of previous libraries, if any, are small: the average differences between the synthetic colors measured on XSL spectra (after slit-loss correction) and on spectra from Miles or (E-)IRTF are below 1%; the dispersion in these comparisons indicates typical individual errors on these broad-band colors, outside the main telluric regions, of 2 to 4% in XSL (with some outliers). Anticipating future publications, we note that comparisons with theoretical spectral libraries across the HR-diagram show clear systematic differences that depend on location in that diagram, with excellent matches in some regions and offsets in certain colors in others; this demonstrates that the current flux calibration allows us to address the scientific questions that were planned. We note that it may be possible to improve the spectrophotometry of XSL spectra in the future, at least at optical wavelengths, pending on additional observations with other instruments (e.g., the integral-field spectrograph MUSE, as suggested by Ivanov et al. 2019).

The residual errors in the relative flux scales of the spectra of the three arms of X-shooter are generally smaller than a few percent. However, there are exceptions, and merging the arms for all the DR2 spectra is a non-trivial process. As previously mentioned, simply fitting the observations to (adequately reddened) theoretical spectra does not provide immediately useful re-scaling factors, as available collections of models do not match the observed energy distributions sufficiently well in all parts of the HR diagram (again, the systematics in these trends show that the issues are not artifacts of flux calibration procedures). Work on merging the arm-spectra is ongoing (Verro, in prep.). This work is tied to the development of a spectral interpolator for XSL, which is needed to predict XSL-based spectra along the isochrones of artificial stellar populations. The comparison between interpolated and original library spectra helps identifying and correcting merged spectra for which the initial relative scaling is inadequate.

The DR2 spectra are made available in the restframe, with radial velocity information in the headers. While the restframe calibration is accurate to typically better than 2 km s⁻¹, the radial velocity estimates are affected by larger instrumental instabilities and should be regarded with caution.

Because XSL extends into the NIR, it emphasises cool and infrared-luminous stars which otherwise tend to be scarce in stellar libraries. Studying the AGB stars and red supergiants by
<table><caption>**Table 8.** XSL keywords dictionary (primary header).</caption><thead><tr><th>General topic</th><th>Step</th><th>Keyword <sup>(1)</sup></th><th>Description</th><th>Default value</th></tr></thead><tbody><tr><td rowspan="5">General information</td><td></td><td>OBJECT</td><td>Machine-readable name of the star</td><td></td></tr><tr><td></td><td>HNAME</td><td>Human-readable name of the star</td><td></td></tr><tr><td></td><td>XSL_ID</td><td>X-shooter Spectral Library unique identifier</td><td></td></tr><tr><td></td><td>XSL_NAME</td><td>Internal name for the science spectrum</td><td></td></tr><tr><td></td><td>PROVi</td><td>Originating raw science file(s)</td><td></td></tr><tr><td rowspan="10">Data reduction</td><td>1D extraction</td><td>EXT_AVG</td><td>Two 1D-spectra averaged</td><td>F</td></tr><tr><td rowspan="3">Response curve</td><td>CAL_NAME</td><td>Internal name for the response curve used</td><td></td></tr><tr><td>CAL_STAR</td><td>Spectro-photometric standard star target</td><td></td></tr><tr><td>CAL_RAWi</td><td>Originating raw flux-standard file(s)</td><td></td></tr><tr><td>First flux calibration</td><td>FLUX_COR</td><td>First flux calibration applied</td><td>F</td></tr><tr><td rowspan="3">Flux-loss correction</td><td>LOSS_COR</td><td>Correction for slit-losses applied</td><td>F</td></tr><tr><td>STA_NAME</td><td>Internal name for the wide-slit spectrum used</td><td></td></tr><tr><td>STA_RAWi</td><td>Originating raw wide-slit file(s)</td><td></td></tr><tr><td>Rest-frame correction</td><td>REST_COR</td><td>Correction to rest-frame applied</td><td>F</td></tr><tr><td></td><td>REST_VAL</td><td>cz value <sup>(2)</sup> [in km s<sup>-1</sup>]</td><td></td></tr><tr><td>Barycor values</td><td>BARY_COR</td><td>Barycentric radial velocity correction value</td><td></td></tr><tr><td rowspan="10">Characterization and quality assurance</td><td rowspan="3">Quality flags</td><td>SNR</td><td>Median signal-to-noise ratio</td><td></td></tr><tr><td>WAVY_UVB</td><td>Wavy spectrum between 460 and 520 nm (suspected residuals of the blaze function)</td><td>F</td></tr><tr><td>WMIN_UVB</td><td>Shortest useful wavelength [in nm]</td><td>300</td></tr><tr><td>HAIR_VIS</td><td>Some narrow spikes in the VIS spectrum</td><td>F</td></tr><tr><td>NOIS_VIS</td><td>Noisy VIS spectrum</td><td>F</td></tr><tr><td>WAVY_VIS</td><td>Wavy in some part of the VIS spectrum</td><td>F</td></tr><tr><td></td><td>WMIN_VIS</td><td>Shortest useful wavelength [in nm]</td><td>540</td></tr><tr><td>HAIR_NIR</td><td>Some narrow spikes in the NIR spectrum</td><td>F</td></tr><tr><td>WAVY_NIR</td><td>Wavy in some part of the NIR spectrum</td><td>F</td></tr></tbody></table>

**Notes.** <sup>(1)</sup>The values assigned to the keywords in the headers of the DR2 data files may depend on the arm (UVB, VIS or NIR), even when this is not explicit in the keyword name. <sup>(2)</sup>The *cz* values are barycentric. They are affected by uncertainties directly related to those on the standard wavelength calibration of the Xshooter data (cf. Sect. 4.3), as illustrated by differences between arms.

themselves will help reduce current discrepancies with synthetic spectra based on static and dynamical stellar atmosphere models (e.g., Lançon et al. 2018), and will contribute to our understanding of young and intermediate age stellar populations. However, XSL remains a general-purpose library that covers a wide range of stellar types and chemical compositions. Fundamental stellar parameters have already been determined and published (Arentsen et al. 2019). We expect XSL will serve a diversity of projects, and in particular improve stellar population synthesis models at a time where many observing facilities concentrate on offering NIR instrumentation.

**Acknowledgements.** We are grateful to the ESO astronomers who obtained the data presented in this paper in Service Mode operations at Paranal Observatory. This research has made use of the SIMBAD database and the VizieR catalog, both operated at CDS, Strasbourg, France. AG is supported by the European

Union FP7 programme through ERC grant number 320360. AG is grateful to the IAC where part of this laborious work has been done. AG and PSB thank the International Space Science Institute in Beijing (ISSI-BJ) for supporting and hosting the meeting of the International Team on “Stellar Libraries of 2020”, during which some discussions contributing to this publication were held. AG and AL thank Cécile Loup for useful discussions during this work. AL, PP, AG gratefully acknowledge repeated support by the Programme National de Physique Stellaire (PNPS) of CNRS/INSU, and by the Programme National Cosmology et Galaxies (PNCG) of CNRS/INSU with INP and IN2P3, co-funded by CEA and CNES. PC acknowledges support from Fundação de Amparo á Pesquisa do Estado de São Paulo (FAPESP) through project 2018/05392-8, and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) through project 310041/2018-0. JFB and AV acknowledge support from grant AYA2016-77237-C3-1-P from the Spanish Ministry of Economy and Competitiveness (MINECO). PSB acknowledges financial support from the Spanish Ministry of Science and Innovation under grant AYA2016-77237-C3-2P. AA gratefully acknowledges funding by the Emmy Noether program from the Deutsche Forschungsgemeinschaft (DFG). We thank the referee for the helpful comments on the text and the website.
References

Aaronson, M., Blanco, V. M., Cook, K. H., & Schechter, P. L. 1989, ApJS, 70, 637

Allende Prieto, C., Koesterke, L., Hubeny, I., et al. 2018, A&A, 618, A25

Arentsen, A., Prugniel, P., Gonneau, A., et al. 2019, A&A, 627, A138

Aringer, B., Girardi, L., Nowotny, W., Marigo, P., & Bressan, A. 2016, MNRAS, 457, 3611

Arp, H. 1965, ApJ, 141, 43

Bagnulo, S., Jehin, E., Ledoux, C., et al. 2003, The Messenger, 114, 10

Bergeat, J., Knapik, A., & Rutily, B. 2002, A&A, 390, 967

Blanco, V. M. 1986, AJ, 91, 290

Blanco, V. M., McCarthy, M. F., & Blanco, B. M. 1984, AJ, 89, 636

Cappellari, M. 2017, MNRAS, 466, 798

Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138

Cenarro, A. J., Cardiel, N., Gorgas, J., et al. 2001, MNRAS, 326, 959

Cenarro, A. J., Peletier, R. F., Sánchez-Blázquez, P., et al. 2007, MNRAS, 374, 664

Chen, Y., Trager, S., Peletier, R., et al. 2014, A&A, 565, A117

Christlieb, N., Green, P. J., Wisotzki, L., & Reimers, D. 2001, A&A, 375, 366

Cioni, M.-R. L., Blommaert, J. A. D. L., Groenewegen, M. A. T., et al. 2003, A&A, 406, 51

Coelho, P. 2009, in American Institute of Physics Conference Series, eds. G. Giobbi, A. Tornambe, G. Raimondo, et al., 1111, 67

Cohen, M., Wheaton, W. A., & Megeath, S. T. 2003, AJ, 126, 1090

Conroy, C., Villaume, A., van Dokkum, P. G., & Lind, K. 2018, ApJ, 854, 139

Cushing, M. C., Rayner, J. T., & Vacca, W. D. 2005, ApJ, 623, 1115

Davies, B., Kudritzki, R.-P., Plez, B., et al. 2013, ApJ, 767, 3

Di Cecco, A., Bono, G., Prada Moroni, P. G., et al. 2015, AJ, 150, 51

Gonneau, A., Lançon, A., Trager, S. C., et al. 2016, A&A, 589, A36

Gonneau, A., Lançon, A., Trager, S. C., et al. 2017, A&A, 601, A141

Gray, R. O., Corbally, C. J., Garrison, R. F., et al. 2006, AJ, 132, 161

Gregg, M. D., Silva, D., Rayner, J., et al. 2006, in The 2005 HST Calibration Workshop: Hubble After the Transition to Two-Gyro Mode, eds. A. M. Koekemoer, P. Goudfrooij, & L. L. Dressel, 209

Groenewegen, M. A. T., & Blommaert, J. A. D. L. 2005, A&A, 443, 143

Horne, K. 1986, PASP, 98, 609

Hughes, S. M. G., & Wood, P. R. 1990, AJ, 99, 784

Husser, T.-O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6

Ivanov, V. D., Rieke, M. J., Engelbracht, C. W., et al. 2004, ApJS, 151, 387

Ivanov, V. D., Coccato, L., Neeser, M. J., et al. 2019, A&A, 629, A100

Johnson, H. L., & Méndez, M. E. 1970, AJ, 75, 785

Kausch, W., Noll, S., Smette, A., et al. 2015, A&A, 576, A78

Koleva, M., Prugniel, P., Bouchard, A., & Wu, Y. 2009, A&A, 501, 1269

Lançon, A., & Rocca-Volmerange, B. 1992, A&AS, 96, 593

Lançon, A., & Wood, P. R. 2000, A&AS, 146, 217

Lançon, A., Mouhcine, M., Fioc, M., & Silva, D. 1999, A&A, 344, L21

Lançon, A., Gallagher, J. S. I., Mouhcine, M., et al. 2008, A&A, 486, 165

Lançon, A., Gonneau, A., Trager, S., et al. 2018, in Why Galaxies Care about AGB Stars, eds. T. Lebzelter, & H. Olofsson, IAU Symp., 343, 0

Levesque, E. M., Massey, P., Olsen, K. A. G., et al. 2005, ApJ, 628, 973

Lynch, D. K., & Mazuk, S. M. 2001, Wavelength Dependence of Cirrus Optical Depth, The Aerospace Corporation, El Segundo, California

Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, AJ, 154, 94

Maraston, C. 2005, MNRAS, 362, 799

Massey, P. 2002, ApJS, 141, 81

McDonald, I., Zijlstra, A. A., & Boyer, M. L. 2012, MNRAS, 427, 343

Melbourne, J., Williams, B. F., Dalcanton, J. J., et al. 2012, ApJ, 748, 47

Modigliani, A., Goldoni, P., Royer, F., et al. 2010, SPIE Conf. Ser., 7737, 773728

Moehler, S. 2015, ESO document XS_wlc_shift_150615.pdf

Moehler, S., Modigliani, A., Freudling, W., et al. 2014, A&A, 568, A9

Nassau, J. J., & Blanco, V. M. 1957, ApJ, 125, 195

Noll, S., Kausch, W., Barden, M., et al. 2012, A&A, 543, A92

Noll, S., Kausch, W., Barden, M., et al. 2014, User Manual for Molecfit VLT-MAN-ESO-19550-5772

Patat, F., Moehler, S., O'Brien, K., et al. 2011, A&A, 527, A91

Pickles, A. J. 1998, PASP, 110, 863

Pompéia, L., Hill, V., Spite, M., et al. 2008, A&A, 480, 379

Prugniel, P., & Soubiran, C. 2001, A&A, 369, 1048

Rayner, J. T., Cushing, M. C., & Vacca, W. D. 2009, ApJS, 185, 289

Reid, N., & Mould, J. 1990, ApJ, 360, 490

Röck, B., Vazdekis, A., Ricciardelli, E., et al. 2016, A&A, 589, A73

Samus, N. N., Kazarovets, E. V., Durlevich, O. V., Kireeva, N. N., & Pastukhova, E. N. 2017, Astron. Rep., 61, 80

Sánchez-Blázquez, P., Peletier, R. F., Jiménez-Vicente, J., et al. 2006, MNRAS, 371, 703

Sharma, K., Prugniel, P., & Singh, H. P. 2016, A&A, 585, A64

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Smette, A., Sana, H., Noll, S., et al. 2015, A&A, 576, A77

Soubiran, C., Le Campion, J.-F., de Cayrel Strobel, G., & Caillo, A. 2010, A&A, 515, A111

Soubiran, C., Le Campion, J.-F., Brouillet, N., & Chemin, L. 2016, A&A, 591, A118

Trager, S. C. 2012, eds. P. Prugniel, & H. P. Singh, Astronomical Society of India Conference Series, 6, 1

Vazdekis, A., Cenarro, A. J., Gorgas, J., Cardiel, N., & Peletier, R. F. 2003, MNRAS, 340, 1317

Vazdekis, A., Koleva, M., Ricciardelli, E., Röck, B., & Falcón-Barroso, J. 2016, MNRAS, 463, 3409

Vernet, J., Dekker, H., D'Odorico, S., et al. 2011, A&A, 536, A105

Villaume, A., Conroy, C., Johnson, B., et al. 2017, ApJS, 230, 23

Whitelock, P. A., Feast, M. W., Marang, F., & Groenewegen, M. A. T. 2006, MNRAS, 369, 751

Wood, P. R., Bessell, M. S., & Fox, M. W. 1983, ApJ, 272, 99

Wu, Y., Singh, H. P., Prugniel, P., Gupta, R., & Koleva, M. 2011, A&A, 525, A71

Yong, D., Norris, J. E., Bessell, M. S., et al. 2013, ApJ, 762, 26

Zoccali, M., Hill, V., Lecureur, A., et al. 2008, A&A, 486, 177
Appendix A: Input catalogs

Table A.1 provides a list of references other than Elodie, MILES, NGSL and PASTEL, that were used to select XSL targets.

<table><caption>**Table A.1.** Input catalogs for stars not in MILES, NGSL or PASTEL.</caption><thead><tr><th>Object name <sup>(1)</sup><br>in Table C.1</th><th>Reference</th><th>Object category</th></tr></thead><tbody><tr><td>HIP*</td><td>Gray et al. (2006)</td><td>Main sequence stars in the MW</td></tr><tr><td>HE*</td><td>Yong et al. (2013)</td><td>Very metal poor stars in the MW</td></tr><tr><td>CD*, CPD*</td><td>Levesque et al. (2005)</td><td>Red supergiants in the MW</td></tr><tr><td>HD*, HE*</td><td>Bergeat et al. (2002), Christlieb et al. (2001)</td><td>C stars in the MW</td></tr><tr><td>IRAS* or GCVS names <sup>(2)</sup></td><td>Whitelock et al. (2006), or GCVS</td><td>LPVs in the MW (mostly C)</td></tr><tr><td>HD 1638</td><td>McDonald et al. (2012)</td><td>K giant in MW (1 object)</td></tr><tr><td>Miscellaneous</td><td>IRTF, Cushing et al. (2005), Rayner et al. (2009)</td><td>Miscellaneous stars in the MW</td></tr><tr><td>Misc. or GCVS names <sup>(2)</sup></td><td>Lançon &amp; Wood (2000)</td><td>Red giants, supergiants and LPVs</td></tr><tr><td>OGLEII DIA BUL*</td><td>Groenewegen &amp; Blommaert (2005)</td><td>LPVs in Bulge fields</td></tr><tr><td>BMB* or 2MASSJ*</td><td>Blanco et al. (1984)</td><td>LPVs near NGC 6522 (near the Bulge)</td></tr><tr><td>CL NGC6522 Arp* or B86*</td><td>Arp (1965), Blanco (1986)</td><td>LPVs near NGC 6522 (near the Bulge)</td></tr><tr><td>2MASSJ*</td><td>Zoccali et al. (2008)</td><td>Giants toward the Bulge</td></tr><tr><td>PHS*</td><td>Pompéia et al. (2008)</td><td>Metal-poor giants in the LMC (3 objects)</td></tr><tr><td>SHV*</td><td>Hughes &amp; Wood (1990)</td><td>LPVs in the LMC</td></tr><tr><td>M2002*</td><td>Massey (2002)</td><td>Red supergiants in the LMC and SMC</td></tr><tr><td>HV*</td><td>Wood et al. (1983)</td><td>LPVs in the LMC and SMC</td></tr><tr><td>CL NGC121*, CL NGC371*</td><td>Wood et al. (1983)</td><td>LPVs in/near SMC clusters</td></tr><tr><td>ISO MCMSJ*</td><td>Cioni et al. (2003)</td><td>LPVs in the SMC</td></tr><tr><td>BBB*</td><td>Reid &amp; Mould (1990)</td><td>SMC supergiants (2 objects)</td></tr></tbody></table>

**Notes.** <sup>(1)</sup>This column provides generic names or examples from the listed references, but in a non-exclusive way. For instance, while some of the objects with HD* names are from Bergeat et al. (2002), most HD* objects are from our main reference catalogs (Elodie, MILES, NGSL, PASTEL). <sup>(2)</sup>By GCVS names, we refer to the label+constellation used in the General Catalog of Variable Stars (e.g., Samus et al. 2017), such as TU Car, V354 Cen, etc.
## Appendix B: Comments about individual objects or spectra

In this DR2 of XSL, we have chosen to include spectra of peculiar objects that made it into the collection for a variety of

reasons, and some spectra with known observational artifacts, because they may be of interest to some future users. We list the most obvious cases in Table B.1. Many of these spectra should not to use in standard population synthesis applications.

<table><caption>**Table B.1.** Peculiar stars, abnormal spectra or deprecated objects.</caption><thead><tr><th>Category</th><th>Spectrum (XSL_ID)</th><th>Star (HNAME)</th><th>Comment</th></tr></thead><tbody><tr><td rowspan="10">Peculiar</td><td>X0085, X0133</td><td>HD 96446</td><td>He-rich</td></tr><tr><td>X0116</td><td>HD 57060</td><td>P Cygni emission in Hα & HeI (1.083 µm); known eclipsing binary</td></tr><tr><td>X0190</td><td>OGLEII DIA BUL-SC19 2332</td><td>O-rich mira, possible TiO emission near 1.25 µm (Lançon et al. 2018)</td></tr><tr><td>X0214, X0248</td><td>HD 190073</td><td>Herbig Ae/Be star, numerous emission lines, continuum contribution from cool component in the NIR arm</td></tr><tr><td>X0357</td><td>CD-69 1618</td><td>He-rich</td></tr><tr><td>X0420</td><td>SV* HV 11223</td><td>S-star signatures confirmed (ZrO, 940 nm; cf. Wood et al. 1983)</td></tr><tr><td>X0424</td><td>CL* NGC 6522 ARP 3213</td><td>Superimposed stars with very different velocities: One dominant in the UVB arm (cz ≈ −198 km s<sup>-1</sup>), the other in the NIR (cz ≈ +30 km s<sup>-1</sup>). The UVB velocity is used for the barycentric correction in all arms.</td></tr><tr><td>X0513</td><td>CL* NGC 419 LE 35</td><td>C-star with contamination by nearby hotter star below 450 nm</td></tr><tr><td>X0653</td><td>HD 172488</td><td>Possibly He-rich</td></tr><tr><td>X0765</td><td>SV* HV 11366</td><td>S-star signatures confirmed (ZrO, 940 nm; cf. Wood et al. 1983)</td></tr><tr><td>X0798</td><td>SHV 0527122-695006</td><td>M giant with contamination by nearby hot star below 450 nm</td></tr><tr><td rowspan="10">Spectral type discussion</td><td>X0478, X0675</td><td>V874 Aql</td><td>O-rich LPV (not C-star as was stated in Nassau & Blanco 1957)</td></tr><tr><td>X0587</td><td>IRAS 10151-6008?</td><td>O-rich star, (not C-star as assumed in Whitelock et al. 2006). We suspect Whitelock et al. observed 2MASS 10165028-6023549, while we observed 2MASS 10165173-6023466.</td></tr><tr><td>X0660</td><td>OGLE 204664c4?</td><td>G8V/K0V star. This name was taken from Zoccali et al. (2008), but does not directly relate to the OGLE survey catalogs. Prefer Gaia DR2 4049054206130689024.</td></tr><tr><td>X0876</td><td>TU Car</td><td>O-rich LPV (as already suggested in Aaronson et al. 1989), not C-star</td></tr><tr><td>X0061</td><td>SHV 0525478-690944</td><td>C star. NIR poor below 1.3 µm.</td></tr><tr><td>X0085</td><td>HD 96446</td><td>Wavy NIR continuum (prefer X0133)</td></tr><tr><td>X0144</td><td>OGLEII DIA BUL-SC26 0532</td><td>A few artificial waves in the VIS arm (residuals from stripes in raw images)</td></tr><tr><td>X0196</td><td>HD 194453</td><td>K-band continuum wavy</td></tr><tr><td>X0208</td><td>HD 179821</td><td>K-band continuum wavy</td></tr><tr><td>X0214</td><td>HD 190073</td><td>K-band continuum wavy (prefer X0248)</td></tr><tr><td>X0221</td><td>HD 175640</td><td>K-band continuum wavy</td></tr><tr><td>X0237</td><td>SHV 0510004-692755</td><td>A few artificial waves in the VIS arm, as for X0144.</td></tr><tr><td>X0269</td><td>HD 163346</td><td>A few artificial waves in the VIS arm, as for X0144.</td></tr><tr><td>X0304</td><td>HD 37828</td><td>Continuum poor between 1 and 1.4 µm (prefer X0097)</td></tr><tr><td>X0338</td><td>R Cha</td><td>Slightly affected by saturation in the H-band [1600–1700 nm] and K-band [2130–2230 nm] (NIR arm)</td></tr><tr><td>X0450</td><td>HD 166991</td><td>Artifacts in K-band (NIR arm)</td></tr><tr><td>X0758, X0766, X0790</td><td>Feige 110</td><td>S/N poor in NIR, use average</td></tr><tr><td>X0778</td><td>HD 11397</td><td>Artifacts in K-band (NIR arm)</td></tr><tr><td>X0834</td><td>HD 42143</td><td>Bad telluric correction (O<sub>2</sub>) in A-band (VIS arm)</td></tr><tr><td>X0878</td><td>HD 83632</td><td>Artifacts in H-band (NIR arm)</td></tr><tr><td rowspan="16">Abnormal</td><td>X0008</td><td>WX Psc</td><td>VIS: lack of signal in this wavelength range.</td></tr><tr><td>X0020</td><td>ISO-MCMSJ005714.4-730121</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0029</td><td>ISO-MCMS J005304.7-730409</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0100</td><td>SHV 0528537-695119</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0130, X0131</td><td>LHS 2065</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0154</td><td>BMB 286</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0253</td><td>OGLEII DIA BUL-SC22 1319</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0296</td><td>OGLEII DIA BUL-SC13 0324</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0353, X0605, X0606, X0635</td><td>[ABC89] Cir18</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0429</td><td>V CrA</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0457</td><td>OGLE-82717</td><td>All three arms: too many different issues and identification uncertain.</td></tr><tr><td>X0544</td><td>SHV 0527072-701238</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0574</td><td>[ABC89] Pup42</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0653</td><td>HD 172488</td><td>UVB: continuum shape unreliable.</td></tr><tr><td>X0759</td><td>HV 12149</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0833</td><td>HD 39970</td><td>All three arms: many issues (including waviness in the UVB, overcorrection of the VIS O<sub>2</sub>).</td></tr><tr><td rowspan="16">Deprecated</td><td>X0008</td><td>WX Psc</td><td>VIS: lack of signal in this wavelength range.</td></tr><tr><td>X0020</td><td>ISO-MCMSJ005714.4-730121</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0029</td><td>ISO-MCMS J005304.7-730409</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0100</td><td>SHV 0528537-695119</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0130, X0131</td><td>LHS 2065</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0154</td><td>BMB 286</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0253</td><td>OGLEII DIA BUL-SC22 1319</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0296</td><td>OGLEII DIA BUL-SC13 0324</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0353, X0605, X0606, X0635</td><td>[ABC89] Cir18</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0429</td><td>V CrA</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0457</td><td>OGLE-82717</td><td>All three arms: too many different issues and identification uncertain.</td></tr><tr><td>X0544</td><td>SHV 0527072-701238</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0574</td><td>[ABC89] Pup42</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0653</td><td>HD 172488</td><td>UVB: continuum shape unreliable.</td></tr><tr><td>X0759</td><td>HV 12149</td><td>UVB: lack of signal in this wavelength range.</td></tr><tr><td>X0833</td><td>HD 39970</td><td>All three arms: many issues (including waviness in the UVB, overcorrection of the VIS O<sub>2</sub>).</td></tr></tbody></table>
## Appendix C: List of XSL stars

Table C.1 details our list of stars and some of their properties.

<table><caption>**Table C.1.** List of XSL observations (ordered by increasing RA).</caption><thead><tr><th>Star</th><th>RA<br>(J2000)</th><th>Dec<br>(J2000)</th><th>XSL_ID</th><th>MJD</th><th>Arm missing<br>or deprecated</th></tr></thead><tbody><tr><td>HD 224926</td><td>00:01:49.45</td><td>-03:01:39.02</td><td>X0022</td><td>55112.9969</td><td></td></tr><tr><td>HD 224926</td><td>00:01:49.45</td><td>-03:01:39.02</td><td>X0224</td><td>55410.3310</td><td></td></tr><tr><td>HD 224926</td><td>00:01:49.45</td><td>-03:01:39.02</td><td>X0280</td><td>55457.2688</td><td></td></tr><tr><td>HD 358</td><td>00:08:23.26</td><td>+29:05:25.55</td><td>X0025</td><td>55113.0508</td><td></td></tr><tr><td>HD 1638</td><td>00:20:28.31</td><td>-30:33:19.88</td><td>X0023</td><td>55113.0104</td><td></td></tr><tr><td>HD 1638</td><td>00:20:28.31</td><td>-30:33:19.88</td><td>X0180</td><td>55407.3051</td><td></td></tr><tr><td>HD 1638</td><td>00:20:28.31</td><td>-30:33:19.88</td><td>X0258</td><td>55439.2983</td><td></td></tr><tr><td>HD 1638</td><td>00:20:28.31</td><td>-30:33:19.88</td><td>X0286</td><td>55461.3079</td><td></td></tr><tr><td>CL* NGC 121 T V8</td><td>00:26:48.53</td><td>-71:32:50.47</td><td>X0385</td><td>56090.3967</td><td></td></tr><tr><td>CL* NGC 121 T V1</td><td>00:26:56.45</td><td>-71:32:07.65</td><td>X0411</td><td>56111.3201</td><td></td></tr><tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr></tbody></table>

**Notes.** The full Table is available at the CDS.