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, 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 $\sim$300 nm and $\sim$2.45 $\mu$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 $\sim\!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.

The X-Shooter Spectral Library (XSL) contains more than 800 spectra of stars across the color-magnitude diagram, that extend from near-UV to near-IR wavelengths (320-2450 nm). We summarize properties of the spectra of O-rich Long Period Variables in XSL, such as phase-related features, and we confront the data with synthetic spectra based on static and dynamical stellar atmosphere models. We discuss successes and remaining discrepancies, keeping in mind the applications to population synthesis modeling that XSL is designed for.

The initial mass function (IMF) is a key ingredient in many studies of galaxy formation and evolution. Although the IMF is often assumed to be universal, there is continuing evidence that it is not universal. Spectroscopic studies that derive the IMF of the unresolved stellar populations of a galaxy often assume that this spectrum can be described by a single stellar population (SSP). To alleviate these limitations, in this paper we have developed a unique hierarchical Bayesian framework for modelling composite stellar populations (CSPs). Within this framework we use a parameterized IMF prior to regulate a direct inference of the IMF. We use this new framework to determine the number of SSPs that is required to fit a set of realistic CSP mock spectra. The CSP mock spectra that we use are based on semi-analytic models and have an IMF that varies as a function of stellar velocity dispersion of the galaxy. Our results suggest that using a single SSP biases the determination of the IMF slope to a higher value than the true slope, although the trend with stellar velocity dispersion is overall recovered. If we include more SSPs in the fit, the Bayesian evidence increases significantly and the inferred IMF slopes of our mock spectra converge, within the errors, to their true values. Most of the bias is already removed by using two SSPs instead of one. We show that we can reconstruct the variable IMF of our mock spectra for signal-to-noise ratios exceeding $\sim$75.

Recent studies based on the integrated light of distant galaxies suggest that the initial mass function (IMF) might not be universal. Variations of the IMF with galaxy type and/or formation time may have important consequences for our understanding of galaxy evolution. We have developed a new stellar population synthesis (SPS) code specifically designed to reconstruct the IMF. We implement a novel approach combining regularization with hierarchical Bayesian inference. Within this approach we use a parametrized IMF prior to regulate a direct inference of the IMF. This direct inference gives more freedom to the IMF and allows the model to deviate from parametrized models when demanded by the data. We use Markov Chain Monte Carlo sampling techniques to reconstruct the best parameters for the IMF prior, the age, and the metallicity of a single stellar population. We present our code and apply our model to a number of mock single stellar populations with different ages, metallicities, and IMFs. When systematic uncertainties are not significant, we are able to reconstruct the input parameters that were used to create the mock populations. Our results show that if systematic uncertainties do play a role, this may introduce a bias on the results. Therefore, it is important to objectively compare different ingredients of SPS models. Through its Bayesian framework, our model is well-suited for this.

We provide a new collection of spectra of 35 carbon stars obtained with the ESO/VLT X-shooter instrument as part of the X-shooter Spectral Library project. The spectra extend from 0.3$\mu$m to 2.4$\mu$m with a resolving power above $\sim$ 8000. The sample contains stars with a broad range of (J-K) color and pulsation properties located in the Milky Way and the Magellanic Clouds. We show that the distribution of spectral properties of carbon stars at a given (J-K) color becomes bimodal (in our sample) when (J-K) is larger than about 1.5. We describe the two families of spectra that emerge, characterized by the presence or absence of the absorption feature at 1.53$\mu$m, generally associated with HCN and C$_2$H$_2$. This feature appears essentially only in large-amplitude variables, though not in all observations. Associated spectral signatures that we interpret as the result of veiling by circumstellar matter, indicate that the 1.53$\mu$m feature might point to episodes of dust production in carbon-rich Miras.