We address the origin of the observed bimodal rotational distribution of stars in massive young and intermediate age stellar clusters. This bimodality is seen as split main sequences at young ages and also has been recently directly observed in the $Vsini$ distribution of stars within massive young and intermediate age clusters. Previous models have invoked binary interactions as the origin of this bimodality, although these models are unable to reproduce all of the observational constraints on the problem. Here we suggest that such a bimodal rotational distribution is set up early within a cluster's life, i.e., within the first few Myr. Observations show that the period distribution of low-mass ($\la 2 M_\odot$) pre-main sequence (PMS) stars is bimodal in many young open clusters and we present a series of models to show that if such a bimodality exists for stars on the PMS that it is expected to manifest as a bimodal rotational velocity (at fixed mass/luminosity) on the main sequence for stars with masses in excess of $\sim1.5$~\msun. Such a bimodal period distribution of PMS stars may be caused by whether stars have lost (rapid rotators) or been able to retain (slow rotators) their circumstellar discs throughout their PMS lifetimes. We conclude with a series of predictions for observables based on our model.

Here we describe a simple, efficient, and most importantly fully operational point-spread-function(PSF)-reconstruction approach for laser-assisted ground layer adaptive optics (GLAO) in the frame of the Multi Unit Spectroscopic Explorer (MUSE) Wide Field Mode. Based on clear astrophysical requirements derived by the MUSE team and using the functionality of the current ESO Adaptive Optics Facility we aim to develop an operational PSF-reconstruction (PSFR) algorithm and test it both in simulations and using on-sky data. The PSFR approach is based on a Fourier description of the GLAO correction to which the specific instrumental effects of MUSE Wide Field Mode (pixel size, internal aberrations, etc.) have been added. It was first thoroughly validated with full end-to-end simulations. Sensitivity to the main atmospheric and AO system parameters was analysed and the code was re-optimised to account for the sensitivity found. Finally, the optimised algorithm was tested and commissioned using more than one year of on-sky MUSE data. We demonstrate with an on-sky data analysis that our algorithm meets all the requirements imposed by the MUSE scientists, namely an accuracy better than a few percent on the critical PSF parameters including full width at half maximum and global PSF shape through the kurtosis parameter of a Moffat function. The PSFR algorithm is publicly available and is used routinely to assess the MUSE image quality for each observation. It can be included in any post-processing activity which requires knowledge of the PSF.

Recent work has shown that the properties of multiple populations within massive stellar clusters (i.e., in the extent of their abundance variations as well as the fraction of stars that show the anomalous chemistry) depend on the mass as well as the age of the host cluster. Such correlations are largely unexpected in current models for the formation of multiple populations and hence provide essential insight into their origin. Here we extend or previous study into the presence or absence of multiple populations using integrated light spectroscopy of the $\sim600$~Myr, massive ($\sim10^7 - 10^8$~\msun) clusters, W3 and W30, in the galactic merger remnant, NGC 7252. Due to the extreme mass of both clusters, the expectation is that they should host rather extreme abundance spreads, manifested through, e.g., high mean [Na/Fe] abundances. However, we do not find evidence for a strong [Na/Fe] enhancement, with the observations being consistent with the solar value. This suggests that age is playing a key role, or alternatively that multiple populations only manifest below a certain stellar mass, as the integrated light at all ages above $\sim100$~Myr is dominated by stars near or above the main sequence turn-off.

The Sagittarius dwarf spheroidal galaxy (Sgr dSph) is in an advanced stage of disruption but still hosts its nuclear star cluster (NSC), M54, at its center. In this paper, we present a detailed kinematic characterization of the three stellar populations present in M54: young metal-rich (YMR); intermediate-age metal-rich (IMR); and old metal-poor (OMP), based on the spectra of $\sim6500$ individual M54 member stars extracted from a large MUSE/VLT dataset. We find that the OMP population is slightly flattened with a low amount of rotation ($\sim0.8$ km s$^{-1}$) and with a velocity dispersion that follows a Plummer profile. The YMR population displays a high amount of rotation ($\sim5$ km s$^{-1}$) and a high degree of flattening, with a lower and flat velocity dispersion profile. The IMR population shows a high but flat velocity dispersion profile, with some degree of rotation ($\sim2$ km s$^{-1}$). We complement our MUSE data with information from \textit{Gaia DR2} and confirm that the stars from the OMP and YMR populations are comoving in 3D space, suggesting that they are dynamically bound. While dynamical evolutionary effects (e.g. energy equipartition) are able to explain the differences in velocity dispersion between the stellar populations, the strong differences in rotation indicate different formation paths for the populations, as supported by an $N$-body simulation tailored to emulate the YMR-OMP system. This study provides additional evidence for the M54 formation scenario proposed in our previous work, where this NSC formed via GC accretion (OMP) and in situ formation from gas accretion in a rotationally supported disc (YMR).

We aim to provide spectroscopic evidence regarding the nature of the putative star cluster in Eridanus 2 and to place constraints on the mass and abundance of massive astrophysical compact halo objects (MACHOs) as a constituent of dark matter. Methods. We present spectroscopic observations of the central square arcminute of Eridanus 2 from MUSE-Faint, a survey of ultra-faint dwarf galaxies with the Multi Unit Spectroscopic Explorer on the Very Large Telescope. We derive line-of-sight velocities for possible member stars of the putative cluster and for stars in the centre of Eridanus 2. We discuss the existence of the cluster and determine new constraints for MACHOs using the Fokker-Planck diffusion approximation. Results. Out of 182 extracted spectra, we identify 26 member stars of Eridanus 2, seven of which are possible cluster members. We find intrinsic mean line-of-sight velocities of $79.7^{+3.1}_{-3.8}\,\mathrm{km}\,\mathrm{s}^{-1}$ and $76.0^{+3.2}_{-3.7}\,\mathrm{km}\,\mathrm{s}^{-1}$ for the cluster and the bulk of Eridanus 2, respectively, and intrinsic velocity dispersions of ${<}7.6\,\mathrm{km}\,\mathrm{s}^{-1}$ (68-$\%$ upper limit) and $10.3^{+3.9}_{-3.2}\,\mathrm{km}\,\mathrm{s}^{-1}$, respectively. This indicates the cluster most likely exists as a distinct dynamical population hosted by Eridanus 2, without surplus of dark matter over the background distribution. Among the member stars in the bulk of Eridanus 2, we find possible carbon stars, alluding to the existence of an intermediate-age population. We derive constraints on the fraction of dark matter that can consist of MACHOs with a given mass between $1$-$10^5\,M_\mathrm{sun}$. For dark matter consisting purely of MACHOs, the mass of the MACHOs must be less than ${\sim}7.6\,M_\mathrm{sun}$ and ${\sim}44\,M_\mathrm{sun}$ at a $68$- and $95$-$\%$ confidence level, respectively. (Abridged)

Aims. We use the spectra of more than 30,000 red giant branch (RGB) stars in 25 globular clusters (GC), obtained within the MUSE survey of Galactic globular clusters, to calibrate the Ca II triplet (CaT) metallicity relation and derive metallicities for all individual stars. We investigate the overall metallicity distributions as well as those of the different populations within each cluster. Methods. The Ca II triplet in the near-infrared at 8498, 8542, and 8662 AA is visible in stars with spectral types between F and M and can be used to determine their metallicities. In this work, we calibrate the relation using average cluster metallicities from literature and MUSE spectra, and extend it below the horizontal branch - a cutoff that has traditionally been made to avoid a non-linear relation - using a quadratic function. In addition to the classic relation based on V-VHB we also present calibrations based on absolute magnitude and luminosity. The obtained relations are then used to calculate metallicities for all the stars in the sample and to derive metallicity distributions for different populations within a cluster, which have been separated using so-called "chromosome maps" based on HST photometry. Results. We show that, despite the relatively low spectral resolution of MUSE (R=1900-3700) we can derive single star metallicities with a mean statistical intra-cluster uncertainty of ~0.12 dex. We present metallicity distributions for the RGB stars in 25 GCs, and investigate the different metallicities of the populations P3 (and higher) in so-called metal-complex or Type II clusters, finding metallicity variations in all of them. We also detected unexpected metallicity variations in the Type I cluster NGC 2808 and confirm the Type II status of NGC 7078.

We present a detailed study of stellar rotation in the massive 1.5 Gyr old cluster NGC 1846 in the Large Magellanic Cloud. Similar to other clusters at this age, NGC 1846 shows an extended main sequence turn-off (eMSTO), and previous photometric studies have suggested it could be bimodal. In this study, we use MUSE integral-field spectroscopy to measure the projected rotational velocities (vsini) of around 1400 stars across the eMSTO and along the upper main sequence of NGC 1846. We measure vsini values up to ~250 km/s and find a clear relation between the vsini of a star and its location across the eMSTO. Closer inspection of the distribution of rotation rates reveals evidence for a bimodal distribution, with the fast rotators centred around vsini = 140 km/s and the slow rotators centred around vsini = 60 km/s. We further observe a lack of fast rotating stars along the photometric binary sequence of NGC 1846, confirming results from the field that suggest that tidal interactions in binary systems can spin down stars. However, we do not detect a significant difference in the binary fractions of the fast and slowly rotating sub-populations. Finally, we report on the serendipitous discovery of a planetary nebula associated with NGC 1846.

Observed mass-to-light ratios (M/L) of metal-rich globular clusters (GCs) disagree with theoretical predictions. This discrepancy is of fundamental importance since stellar population models provide the stellar masses that underpin most of extragalactic astronomy, near and far. We have derived radial velocities for 1,622 stars located in the centres of 59 Milky Way GCs - twelve of which have no previous kinematic information - using integral-field unit data from the WAGGS project. Using N-body models, we then determine dynamical masses and M/L ratios for the studied clusters. Our sample includes NGC 6528 and NGC 6553, which extend the metallicity range of GCs with measured M/L up to [Fe/H] ~ -0.1 dex. We find that metal-rich clusters have M/L more than two times lower than what is predicted by simple stellar population models. This confirms that the discrepant M/L-[Fe/H] relation remains a serious concern. We explore how our findings relate to previous observations, and the potential causes for the divergence, which we conclude is most likely due to dynamical effects.

We combine MUSE spectroscopy and Hubble Space Telescope ultraviolet (UV) photometry to perform a study of the chemistry and dynamics of the Galactic globular cluster Messier 80 (M80, NGC 6093). Previous studies have revealed three stellar populations that not only vary in their light-element abundances, but also in their radial distributions, with concentration decreasing with increasing nitrogen enrichment. This remarkable trend, which sets M80 apart from the other Galactic globular clusters, points towards a complex formation and evolutionary history. To better understand how M80 formed and evolved, revealing its internal kinematics is key. We find that the most N-enriched population rotates faster than the other two populations at a 2 sigma confidence level. While our data further suggest that the intermediate population shows the least amount of rotation, this trend is rather marginal (1 - 2 sigma). Using axisymmetric Jeans models, we show that these findings can be explained from the radial distributions of the populations if they possess different angular momenta. Our findings suggest that the populations formed with primordial kinematical differences.

We investigate the binary content of the two stellar populations that coexist in the globular cluster NGC 3201. Previous studies of binary stars in globular clusters have reported higher binary fractions in their first populations (P1, having field-like abundances) compared to their second populations (P2, having anomalous abundances). This is interpreted as evidence for the latter forming more centrally concentrated. In contrast to previous studies, our analysis focuses on the cluster centre, where comparable binary fractions between the populations are predicted because of the short relaxation times. However, we find that even in the centre of NGC 3201, the observed binary fraction of P1 is higher, (23.1 +/- 6.2)% compared to (8.2 +/- 3.5)% in P2. Our results are difficult to reconcile with a scenario where the populations only differ in their initial concentrations, but instead suggests that the populations also formed with different fractions of binary stars.

We investigate the morphology of the colour-magnitude diagram (CMD) of the open cluster NGC 2509 in comparison with other Galactic open clusters of similar age using Gaia photometry. At $\sim900$ Myr Galactic open clusters in our sample all show an extended main sequence turn off (eMSTO) with the exception of NGC 2509, which presents an exceptionally narrow CMD. Our analysis of the Gaia data rules out differential extinction, stellar density, and binaries as a cause for the singular MSTO morphology in this cluster. We interpret this feature as a consequence of the stellar rotation distribution within the cluster and present the analysis with MIST stellar evolution models that include the effect of stellar rotation on which we based our conclusion. In particular, these models point to an unusually narrow range of stellar rotation rates ($\Omega/\Omega_{\rm{crit,ZAMS}} = [0.4, 0.6]$) within the cluster as the cause of this singular feature in the CMD of NGC 2509. Interestingly, models that do not include rotation are not as good at reproducing the morphology of the observed CMD in this cluster.

Context. Moderately metal-poor inner bulge globular clusters are relics of a generation of long-lived stars that formed in the early Galaxy. Terzan 9, projected at 4d 12 from the Galactic center, is among the most central globular clusters in the Milky Way, showing an orbit which remains confined to the inner 1 kpc. Aims. Our aim is the derivation of the cluster's metallicity, together with an accurate measurement of the mean radial velocity. In the literature, metallicities in the range between have been estimated for Terzan 9 based on color-magnitude diagrams and CaII triplet (CaT) lines. Aims. Our aim is the derivation of the cluster's metallicity, together with an accurate measurement of the mean radial velocity. In the literature, metallicities in the range between -2.0 and -1.0 have been estimated for Terzan 9 based on color-magnitude diagrams and CaII triplet (CaT) lines. Methods. Given its compactness, Terzan 9 was observed using the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope. The extraction of spectra from several hundreds of individual stars allowed us to derive their radial velocities, metallicities, and [Mg/Fe]. The spectra obtained with MUSE were analysed through full spectrum fitting using the ETOILE code. Results. We obtained a mean metallicity of [Fe/H] -1.10 0.15, a heliocentric radial velocity of vhr = 58.1 1.1 km/s , and a magnesium-to-iron [Mg/Fe] = 0.27 0.03. The metallicity-derived character of Terzan 9 sets it among the family of the moderately metal-poor Blue Horizontal Branch clusters HP 1, NGC 6558, and NGC 6522.

The center of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) hosts a nuclear star cluster, M54, which is the only galaxy nucleus that can be resolved into individual stars at optical wavelengths. It is thus a key target for understanding the formation of nuclear star clusters and their relation to globular clusters. We present a large Multi-Unit Spectroscopic Explorer (MUSE) data set that covers M54 out to $\sim$2.5 half-light radius, from which we extracted the spectra of $\sim$6600 cluster member stars. We use these data in combination with HST photometry to derive age and metallicity for each star. The stellar populations show a well defined age-metallicity relation, implying an extended formation history for the central region of Sgr dSph. We classify these populations into three groups, all with the same systemic velocity: young metal-rich (YMR; 2.2\,Gyr, \mbox{[Fe/H]$=-0.04$}); intermediate-age metal-rich (IMR; 4.3\,Gyr, \mbox{[Fe/H]$=-0.29$}); and old metal-poor (OMP; 12.2\,Gyr, \mbox{[Fe/H]$=-1.41$}). The YMR and OMP populations are more centrally concentrated than the IMR population, which are likely stars of the Sgr dSph. We suggest the OMP population is the result of accretion and merging of two or more old and metal-poor globular clusters dragged to the center by dynamical friction. The YMR is consistent with being formed by in situ star formation in the nucleus. The ages of the YMR population suggest that it may have been triggered into forming when the Sgr dSph began losing its gas during the most recent interaction with the Milky Way, $\sim$3\,Gyr ago.

Galactic globular clusters (GCs) are now known to host multiple populations displaying particular abundance variations. The different populations within a GC can be well distinguished following their position in the pseudo two-colors diagrams, also referred to as "chromosome maps". These maps are constructed using optical and near-UV photometry available from the Hubble Space Telescope (HST) UV survey of GCs. However, the chemical tagging of the various populations in the chromosome maps is hampered by the fact that HST photometry and elemental abundances are both available only for a limited number of stars. The spectra collected as part of the MUSE survey of globular clusters provide a spectroscopic counterpart to the HST photometric catalogs covering the central regions of GCs. In this paper, we use the MUSE spectra of 1115 red giant branch (RGB) stars in NGC 2808 to characterize the abundance variations seen in the multiple populations of this cluster. We use the chromosome map of NGC 2808 to divide the RGB stars into their respective populations. We then combine the spectra of all stars belonging to a given population, resulting in one high signal-to-noise ratio spectrum representative of each population. Variations in the spectral lines of O, Na, Mg, and Al are clearly detected among four of the populations. In order to quantify these variations, we measured equivalent width differences and created synthetic populations spectra that were used to determine abundance variations with respect to the primordial population of the cluster. Our results are in good agreement with the values expected from previous studies based on high-resolution spectroscopy. We do not see any significant variations in the spectral lines of Ca, K, and Ba. We also do not detect abundance variations among the stars belonging to the primordial population of NGC 2808.

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

Globular clusters produce many exotic stars due to a much higher frequency of dynamical interactions in their dense stellar environments. Some of these objects were observed together with several hundred thousands other stars in our MUSE survey of 26 Galactic globular clusters. Assuming that at least a few exotic stars have exotic spectra, that means spectra that contain emission lines, we can use this large spectroscopic data set of over a million stellar spectra as a blind survey to detect stellar exotica in globular clusters. To detect emission lines in each spectrum, we model the expected shape of an emission line as a Gaussian curve. This template is used for matched filtering on the differences between each observed 1D spectrum and its fitted spectral model. The spectra with the most significant detections of H$\alpha$ emission are checked visually and cross-matched with published catalogues. We find 156 stars with H$\alpha$ emission, including several known cataclysmic variables (CV) and two new CVs, pulsating variable stars, eclipsing binary stars, the optical counterpart of a known black hole, several probable sub-subgiants and red stragglers, and 21 background emission-line galaxies. We find possible optical counterparts to 39 X-ray sources, as we detect H$\alpha$ emission in several spectra of stars that are close to known positions of Chandra X-ray sources. This spectral catalogue can be used to supplement existing or future X-ray or radio observations with spectra of potential optical counterparts to classify the sources.

We utilize multi-epoch MUSE spectroscopy to study binaries in the core of NGC 3201. Our sample consists of 3553 stars with 54883 spectra in total comprising 3200 main-sequence stars up to 4 magnitudes below the turn-off. Each star in our sample has between 3 and 63 (with a median of 14) reliable radial velocity (RV) measurements within five years of observations. We introduce a statistical method to determine the probability of a star showing RV variations based on the whole inhomogeneous RV sample. Using HST photometry and an advanced dynamical MOCCA simulation of this specific GC we overcome observational biases that previous spectroscopic studies had to deal with. This allows us to infer a binary frequency in the MUSE FoV and enables us to deduce the underlying true binary frequency of (6.75+-0.72) % in NGC 3201. The comparison of the MUSE observations with the MOCCA simulation suggests a significant fraction of primordial binaries. We can also confirm a radial increase of the binary fraction towards the GC centre due to mass segregation. We discovered that in our sample at least (57.5+-7.9) % of blue straggler stars (BSS) are in a binary system. For the first time in a study of GCs, we were able to fit Keplerian orbits to a significant sample of 95 binaries. We present the binary system properties of eleven BSS and show evidence that two BSS formation scenarios, the mass transfer in binary (or triple) star systems and the coalescence due to binary-binary interactions, are present in our data. We also describe the binary and spectroscopic properties of four sub-subgiant (or red straggler) stars. Furthermore, we discovered two new black hole (BH) candidates with minimum masses (Msini) of (7.68+-0.50) M_sun, (4.4+-2.8) M_sun, and refine the minimum mass estimate on the already published BH to (4.53+-0.21) M_sun. These BHs are consistent with an extensive BH subsystem hosted by NGC 3201.

The presence of star-to-star light-element abundance variations (a.k.a. multiple populations, MPs) appears to be ubiquitous within old and massive clusters in the Milky Way and all studied nearby galaxies. Most previous studies have focussed on resolved images or spectroscopy of individual stars, although there has been significant effort in the past few years to look for multiple population signatures in integrated light spectroscopy. If proven feasible, integrated light studies offer a potential way to vastly open parameter space, as clusters out to tens of Mpc can be studied. We use the NaD lines in the integrated spectra of two clusters with similar ages ($2-3$ Gyr) but very different masses, NGC 1978 ($\sim3\times10^5$ Msun) in the LMC and G114 ($1.7\times10^7$ Msun) in NGC 1316. For NGC 1978, our findings agree with resolved studies of individual stars which did not find evidence for Na spreads. However, for G114, we find clear evidence for the presence of multiple populations. The fact that the same anomalous abundance patterns are found in both the intermediate age and ancient GCs lends further support to the notion that young massive clusters are effectively the same as the ancient globular clusters, only separated in age.

We compare the results of a large grid of N-body simulations with the surface brightness and velocity dispersion profiles of the globular clusters $\omega$ Cen and NGC 6624. Our models include clusters with varying stellar-mass black hole retention fractions and varying masses of a central intermediate-mass black hole (IMBH). We find that an $\sim 45,000$ M$_\odot$ IMBH, whose presence has been suggested based on the measured velocity dispersion profile of $\omega$ Cen, predicts the existence of about 20 fast-moving, $m>0.5$ M$_\odot$ main-sequence stars with a (1D) velocity $v>60$ km/sec in the central 20 arcsec of $\omega$ Cen. However no such star is present in the HST/ACS proper motion catalogue of Bellini et al. (2017), strongly ruling out the presence of a massive IMBH in the core of $\omega$ Cen. Instead, we find that all available data can be fitted by a model that contains 4.6% of the mass of $\omega$ Cen in a centrally concentrated cluster of stellar-mass black holes. We show that this mass fraction in stellar-mass BHs is compatible with the predictions of stellar evolution models of massive stars. We also compare our grid of $N$-body simulations with NGC 6624, a cluster recently claimed to harbor a 20,000 M$_\odot$ black hole based on timing observations of millisecond pulsars. However, we find that models with $M_{IMBH}>1,000$ M$_\odot$ IMBHs are incompatible with the observed velocity dispersion and surface brightness profile of NGC 6624,ruling out the presence of a massive IMBH in this cluster. Models without an IMBH provide again an excellent fit to NGC 6624.

A nova is a cataclysmic event on the surface of a white dwarf in a binary system that increases the overall brightness by several orders of magnitude. Although binary systems with a white dwarf are expected to be overabundant in globular clusters (GCs) compared to the Galaxy, only two novae from Galactic globular clusters have been observed. We present the discovery of an emission nebula in the Galactic globular cluster M 22 (NGC 6656) in observations made with the integral-field spectrograph MUSE. We extract the spectrum of the nebula and use the radial velocity determined from the emission lines to confirm that the nebula is part of NGC 6656. Emission-line ratios are used to determine the electron temperature and density. It is estimated to have a mass of 1 to $17 \times 10^{-5}$ solar masses. This mass and the emission-line ratios indicate that the nebula is a nova remnant. Its position coincides with the reported location of a 'guest star', an ancient Chinese term for transients, observed in May 48 BCE. With this discovery, this nova may be one of the oldest confirmed extrasolar events recorded in human history.

We present the first stellar spectroscopy in the low-luminosity (M_V ~-9.3 mag), dwarf galaxy Leo P. Its significantly low oxygen abundance (3% solar) and relative proximity (~1.6 Mpc) make it a unique galaxy to investigate the properties of massive stars with near-primordial compositions akin to those in the early Universe. From our VLT-MUSE spectroscopy we find the first direct evidence for an O-type star in the prominent HII region, providing an important test case to investigate the potential environmental dependence of the upper end of the initial mass function in the dwarf galaxy regime. We classify 14 further sources as massive stars (and 17 more as candidate massive stars), most likely B-type objects. From comparisons with published evolutionary models we argue that the absolute visual magnitudes of massive stars in very metal-poor systems such as Leo P and I Zw 18 may be fainter by ~0.5 mag compared to Galactic stars. We also present spectroscopy of two carbon stars identified previously as candidate asymptotic-giant-branch stars. Two of three further candidate asymptotic-giant-branch stars display CaII absorption, confirming them as cool, evolved stars; we also recover CaII absorption in the stacked data of the next brightest 16 stars in the upper red giant branch. These discoveries will provide targets for future observations to investigate the physical properties of these objects and to calibrate evolutionary models of luminous stars at such low metallicity. The MUSE data also reveal two 100pc-scale ring structures in Halpha emission, with the HII region located on the northern edge of the southern ring. Lastly, we report serendipitous observations of 20 galaxies, with redshifts ranging from z=0.39, to a close pair of star-forming galaxies at z=2.5.

The physics governing the formation of star clusters is still not entirely understood. One open question concerns the amount of angular momentum that newly formed clusters possess after emerging from their parent gas clouds. Recent results suggest an alignment of stellar spins and binary orbital spins in star clusters, which support a scenario in which clusters are born with net angular momentum cascading down to stellar scales. In this paper, we combine Gaia data and published line of sight velocities to explore if NGC6791 and NGC6819, two of the clusters for which an alignment of stellar spins has been reported, rotate in the same plane as their stars. We find evidence for rotation in NGC6791 using both proper motions and line of sight velocities. Our estimate of the inclination angle is broadly consistent with the mean inclination that has been determined for its stars, but the uncertainties are still substantial. Our results identify NGC6791 as a promising follow-up candidate to investigate the link between cluster and stellar rotation. We find no evidence for rotation in NGC6819.

We use Hubble Space Telescope multicolour photometry of the globular cluster 47 Tucanae to uncover a population of 24 objects with no previous classification that are outliers from the single-star model tracks in the colour-magnitude diagram and yet are likely cluster members. By comparing those sources with evolutionary models and X-ray source catalogues, we were able to show that the majority of those sources are likely binary systems that do not have any X-ray source detected nearby, most possibly formed by a white dwarf and a main-sequence star and a small number of possible double-degenerate systems.

Using data from the WiFeS Atlas of Galactic Globular cluster Spectra we study the behaviour of the calcium triplet (CaT), a popular metallicity indicator in extragalactic stellar population studies. A major caveat of these studies is that the potential sensitivity to other stellar population parameters such as age, calcium abundance and the initial mass function has not yet been empirically evaluated. Here we present measurements of the strength of the CaT feature for 113 globular clusters in the Milky Way and its satellite galaxies. We derive empirical calibrations between the CaT index and both the iron abundance ([Fe/H]) and calcium abundance ([Ca/H]), finding a tighter relationship for [Ca/H] than for [Fe/H]. For stellar populations 3 Gyr and older the CaT can be used to reliably measure [Ca/H] at the 0.1 dex level but becomes less reliable for ages of $\sim 2$ Gyr and younger. We find that the CaT is relatively insensitive to the horizontal branch morphology. The stellar mass function however affects the CaT strengths significantly only at low metallicities. Using our newly derived empirical calibration, we convert our measured CaT indices into [Ca/H] values for the globular clusters in our sample.

Westerlund 2 (Wd2) is the central ionizing star cluster of the \ion{H}{2} region RCW~49 and the second most massive young star cluster (${\rm M} = (3.6 \pm 0.3)\times 10^4\,{\rm M}_\odot$) in the Milky Way. Its young age ($\sim2\,$Myr) and close proximity to the Sun ($\sim 4\,$kpc) makes it a perfect target to study stars emerging from their parental gas cloud, the large number of OB-stars and their feedback onto the gas, and the gas dynamics. We combine high-resolution multi-band photometry obtained in the optical and near-infrared with the \textit{Hubble} Space Telescope (HST), and VLT/MUSE integral field spectroscopy to study the gas, the stars, and their interactions, simultaneously. In this paper we focus on a small, $64\times64\,{\rm arcsec}^2$ region North of the main cluster center, which we call the Northern Bubble (NB), a circular cavity carved into the gas of the cluster region. Using MUSE data, we determined the spectral types of 17 stars in the NB from G9III to O7.5. With the estimation of these spectral types we add 2 O and 5 B-type stars to the previously published census of 37 OB-stars in Wd2. To measure radial velocities we extracted 72 stellar spectra throughout Wd2, including the 17 of the NB, and show that the cluster member stars follow a bimodal velocity distribution centered around $(8.10 \pm 1.53)\,{\rm km}\,{\rm s}^{-1}$ and $(25.41 \pm 1.57)\,{\rm km}\,{\rm s}^{-1}$ with a dispersion of $(4.52 \pm 1.78)\,{\rm km}\,{\rm s}^{-1}$ and $(3.46 \pm 1.29)\,{\rm km}\,{\rm s}^{-1}$, respectively. These are in agreement with CO($J=1$-2) studies of RCW~49 leaving cloud-cloud collision as a viable option for the formation scenario of Wd2. The bimodal distribution is also detected in the Gaia DR2 proper motions.

We present an analysis of the relatively low mass ($\sim2400$~M$_{\odot}$), $\sim800$~Myr, Galactic open cluster, NGC~2818, using Gaia DR2 results combined with VLT/FLAMES spectroscopy. Using Gaia DR2 proper motions and parallax measurements we are able to select a clean sample of cluster members. This cluster displays a clear extended main sequence turn-off (eMSTO), a phenomenon previously studied mainly in young and intermediate age massive clusters in the Magellanic clouds. The main sequence of NGC~2818 is extremely narrow, with a width of $\sim0.01$ magnitudes (G$_{\rm BP} - $ G$_{\rm RP}$), suggesting very low levels of differential extinction. Using VLT/FLAMES spectroscopy of 60 cluster members to measure the rotational velocity of the stars (Vsini) we find that stars on the red side of the eMSTO have high Vsini ($>160$~km/s) while stars on the blue side have low Vsini ($<160$~km/s), in agreement with model predictions. The cluster also follows the previously discovered trend between the age of the cluster and the extent of the eMSTO. We conclude that stellar rotation is the likely cause of the eMSTO phenomenon.

We present adaptive optics (AO) assisted integral-field spectroscopy of the intermediate-age star cluster NGC 419 in the Small Magellanic Cloud. By investigating the cluster dynamics and the rotation properties of main sequence turn-off stars (MSTO), we demonstrate the power of AO-fed MUSE observations for this class of objects. Based on 1 049 radial velocity measurements, we determine a dynamical cluster mass of 1.4+/-0.2x10^5 M_sun and a dynamical mass-to-light ratio of 0.67+/-0.08, marginally higher than simple stellar population predictions for a Kroupa initial mass function. A stacking analysis of spectra at both sides of the extended MSTO reveals significant rotational broadening. Our analysis further provides tentative evidence that red MSTO stars rotate faster than their blue counterparts. We find average V sin i values of 87+/-16 km/s and 130+/-22 km/s for blue and red MSTO stars, respectively. Potential systematic effects due to the low spectral resolution of MUSE can reach 30 km/s but the difference in V sin i between the populations is unlikely to be affected.

Aims. As a new approach to the study of resolved stellar populations in nearby galaxies, our goal is to demonstrate in NGC300 that integral field spectroscopy with high spatial resolution and excellent seeing conditions reaches an unprecedented depth in severely crowded fields. Methods. MUSE observations with seven pointings in NGC300 have resulted in datacubes that are analyzed in four ways: (1) PSF-fitting 3D spectroscopy with PampelMUSE yields deblended spectra of individually distinguishable stars. The technique also provides samples of planetary nebulae that are complete down to m5007=28. (2) pseudo-monochromatic images, created at the wavelengths of the most important emission lines and corrected for continuum light by using the P3D visualization tool, provide maps of HII regions, SNR, and the diffuse ISM at a high level of sensitivity, allowing for the discovery of planetary nebulae, WR stars etc. (3) The use of the P3D line-fitting tool yields emission line fluxes, surface brightness, and kinematic information for gaseous objects, corrected for absorption line profiles of the underlying stellar population. (4) Visual inspection of the datacubes is demonstrated to be effcient for data mining and the discovery of background galaxies and unusual objects. Results. We present a catalogue of luminous stars, rare stars such as WR and other emission line stars, carbon stars, symbiotic star candidates, planetary nebulae, HII regions, supernova remnants, giant shells, peculiar diffuse and filamentary emission line objects, and background galaxies, along with their spectra. Conclusions. The technique of crowded-field 3D spectroscopy is capable of deblending individual bright stars, the unresolved background of faint stars, gaseous nebulae, and the diffuse component of the interstellar medium, resulting in unprecedented legacy value for observations of nearby galaxies with MUSE.

We present a detailed study of the nuclear star clusters (NSCs) and massive black holes (BHs) of four of the nearest low-mass early-type galaxies: M32, NGC205, NGC5012, and NGC5206. We measure dynamical masses of both the BHs and NSCs in these galaxies using Gemini/NIFS or VLT/SINFONI stellar kinematics, Hubble Space Telescope (HST) imaging, and Jeans Anisotropic Models. We detect massive BHs in M32, NGC5102, and NGC5206, while in NGC205, we find only an upper limit. These BH mass estimates are consistent with previous measurements in M32 and NGC205, while those in NGC5102 and NGC5206 are estimated for the first time, and both found to be $<$$10^6~M_{\odot}$. This adds to just a handful of galaxies with dynamically measured sub-million $M_{\odot}$ central BHs. Combining these BH detections with our recent work on NGC404's BH, we find that 80\% (4/5) of nearby, low-mass ($10^9-10^{10}~M_{\odot}$; $\sigma_{\star}\sim20-70$ km/s) early-type galaxies host BHs. Such a high occupation fraction suggests the BH seeds formed in the early epoch of cosmic assembly likely resulted in abundant seeds, favoring a low-mass seed mechanism of the remnants, most likely from the first generation of massive stars. We find dynamical masses of the NSCs ranging from $2-80\times10^6~M_{\odot}$ and compare these masses to scaling relations for NSCs based primarily on photometric mass estimates. Color gradients suggest younger stellar populations lie at the centers of the NSCs in three of the four galaxies (NGC205, NGC5102, and NGC5206), while the morphology of two are complex and are best-fit with multiple morphological components (NGC5102 and NGC5206). The NSC kinematics show they are rotating, especially in M32 and NGC5102 ($V/\sigma_{\star}\sim0.7$).

As part of our massive spectroscopic survey of 25 Galactic globular clusters with MUSE, we performed multiple epoch observations of NGC 3201 with the aim of constraining the binary fraction. In this cluster, we found one curious star at the main-sequence turn-off with radial velocity variations of the order of 100 km/s, indicating the membership to a binary system with an unseen component since no other variations appear in the spectra. Using an adapted variant of the generalized Lomb-Scargle periodogram, we could calculate the orbital parameters and found the companion to be a detached stellar-mass black hole with a minimum mass of 4.36 $\pm$ 0.41 solar masses. The result is an important constraint for binary and black hole evolution models in globular clusters as well as in the context of gravitational wave sources.

The "Antennae Galaxy" (NGC 4038/39) is the closest major interacting galaxy system and therefore often taken as merger prototype. We present the first comprehensive integral field spectroscopic dataset of this system, observed with the MUSE instrument at the ESO VLT. We cover the two regions in this system which exhibit recent star-formation: the central galaxy interaction and a region near the tip of the southern tidal tail. In these fields, we detect HII regions and diffuse ionized gas to unprecedented depth. About 15% of the ionized gas was undetected by previous observing campaigns. This newly detected faint ionized gas is visible everywhere around the central merger, and shows filamentary structure. We estimate diffuse gas fractions of about 60% in the central field and 10% in the southern region. We are able to show that the southern region contains a significantly different population of HII regions, showing fainter luminosities. By comparing HII region luminosities with the HST catalog of young star clusters in the central field, we estimate that there is enough Lyman-continuum leakage in the merger to explain the amount of diffuse ionized gas that we detect. We compare the Lyman-continuum escape fraction of each HII region against ionization-parameter sensitive emission line ratios. While we find no systematic trend between these properties, the most extreme line ratios seem to be strong indicators of density bounded ionization. Extrapolating the Lyman-continuum escape fractions to the southern region, we conclude that just from the comparison of the young stellar populations to the ionized gas there is no need to invoke other ionization mechanisms than Lyman-continuum leaking HII regions for the diffuse ionized gas in the Antennae.

This is the first of a series of papers presenting the results from our survey of 25 Galactic globular clusters with the MUSE integral-field spectrograph. In combination with our dedicated algorithm for source deblending, MUSE provides unique multiplex capabilities in crowded stellar fields and allows us to acquire samples of up to 20 000 stars within the half-light radius of each cluster. The present paper focuses on the analysis of the internal dynamics of 22 out of the 25 clusters, using about 500 000 spectra of 200 000 individual stars. Thanks to the large stellar samples per cluster, we are able to perform a detailed analysis of the central rotation and dispersion fields using both radial profiles and two-dimensional maps. The velocity dispersion profiles we derive show a good general agreement with existing radial velocity studies but typically reach closer to the cluster centres. By comparison with proper motion data we derive or update the dynamical distance estimates to 14 clusters. Compared to previous dynamical distance estimates for 47 Tuc, our value is in much better agreement with other methods. We further find significant (>3sigma) rotation in the majority (13/22) of our clusters. Our analysis seems to confirm earlier findings of a link between rotation and the ellipticities of globular clusters. In addition, we find a correlation between the strengths of internal rotation and the relaxation times of the clusters, suggesting that the central rotation fields are relics of the cluster formation that are gradually dissipated via two-body relaxation.

We map the interstellar medium (ISM) including the diffuse interstellar bands (DIBs) in absorption toward the globular cluster NGC 6397 using VLT/MUSE. This pilot study demonstrates the power of MUSE for mapping the local ISM on very small scales which provides a new window for ISM observations. Assuming the absorbers are located at the rim of the Local Bubble we trace small-scale variations in NaI and KI as well as in several DIBs structures on the order of mpc (milliparsec, a few thousand AU). The sightlines defined by binned stellar spectra are separated by only a few arcseconds and we probe the absorption within a physically connected region. This analysis utilized the fitting residuals of individual stellar spectra of NGC 6397 member stars and analyzed lines from neutral species and several DIBs in Voronoi-binned composite spectra with high signal-to-noise ratio (S/N). We verify the suitability of the MUSE 3D spectrograph for such measurements and gain new insights by probing a single physical absorber with multiple sight lines.

Sub-subgiants are stars observed to be redder than normal main-sequence stars and fainter than normal subgiant (and giant) stars in an optical color-magnitude diagram. The red straggler stars, which lie redward of the red giant branch, may be related and are often grouped together with the sub-subgiants in the literature. These stars defy our standard theory of single-star evolution, and are important tests for binary evolution and stellar collision models. In total, we identify 65 sub-subgiants and red stragglers in 16 open and globular star clusters from the literature; 50 of these, including 43 sub-subgiants, pass our strict membership selection criteria (though the remaining sources may also be cluster members). In addition to their unique location on the color-magnitude diagram, we find that at least 58% (25/43) of sub-subgiants in this sample are X-ray sources with typical 0.5-2.5 keV luminosities of order 10^30 - 10^31 erg/s. Their X-ray luminosities and optical-to-X-ray flux ratios are similar to those of RS CVn active binaries. At least 65% (28/43) of the sub-subgiants in our sample are variables, 21 of which are known to be radial-velocity binaries. Typical variability periods are <15 days. At least 33% (14/43) of the sub-subgiants are H-alpha emitters. These observational demographics provide strong evidence that binarity is important for sub-subgiant formation. Finally, we find that the number of sub-subgiants per unit mass increases toward lower-mass clusters, such that the open clusters in our sample have the highest specific frequencies of sub-subgiants.

We describe the dynamical evolution of a unique type of dark star cluster model in which the majority of the cluster mass at Hubble time is dominated by an intermediate-mass black hole (IMBH). We analyzed results from about 2000 star cluster models (Survey Database I) simulated using the Monte-Carlo code MOCCA and identified these dark star cluster models. Taking one of these models, we apply the method of simulating realistic "mock observations" by utilizing the COCOA and SISCO codes to obtain the photometric and kinematic observational properties of the dark star cluster model at 12 Gyr. We find that the perplexing Galactic globular cluster NGC 6535 closely matches the observational photometric and kinematic properties of the dark star cluster model presented in this paper. Based on our analysis and currently observed properties of NGC 6535, we suggest that this globular cluster could potentially harbour an IMBH. If it exists, the presence of this IMBH can be detected robustly with proposed kinematic observations of NGC 6535.

Aims. We demonstrate the high multiplex advantage of crowded field 3D spectroscopy using the new integral field spectrograph MUSE by means of a spectroscopic analysis of more than 12,000 individual stars in the globular cluster NGC 6397. Methods. The stars are deblended with a PSF fitting technique, using a photometric reference catalogue from HST as prior, including relative positions and brightnesses. This catalogue is also used for a first analysis of the extracted spectra, followed by an automatic in-depth analysis using a full-spectrum fitting method based on a large grid of PHOENIX spectra. Results. With 18,932 spectra from 12,307 stars in NGC 6397 we have analysed the largest sample so far available for a single globular cluster. We derived a mean radial velocity of vrad=17.84+-0.07 km/s and a mean metallicity of [Fe/H]=-2.120+-0.002, with the latter seemingly varying with temperature for stars on the RGB. We determine T_eff and [Fe/H] from the spectra, and log g from HST photometry. This is the first very comprehensive HRD for a globular cluster based on the analysis of several thousands of stellar spectra, ranging from the main sequence to the tip of the RGB. Furthermore, two interesting objects were identified with one being a post-AGB star and the other a possible millisecond-pulsar companion.

We present a detailed analysis of the kinematics of the Galactic globular cluster NGC 6397 based on more than ~18,000 spectra obtained with the novel integral field spectrograph MUSE. While NGC 6397 is often considered a core collapse cluster, our analysis suggests a flattening of the surface brightness profile at the smallest radii. Although it is among the nearest globular clusters, the low velocity dispersion of NGC 6397 of <5km/s imposes heavy demands on the quality of the kinematical data. We show that despite its limited spectral resolution, MUSE reaches an accuracy of 1km/s in the analysis of stellar spectra. We find slight evidence for a rotational component in the cluster and the velocity dispersion profile that we obtain shows a mild central cusp. To investigate the nature of this feature, we calculate spherical Jeans models and compare these models to our kinematical data. This comparison shows that if a constant mass-to-light ratio is assumed, the addition of an intermediate-mass black hole with a mass of 600M_sun brings the model predictions into agreement with our data, and therefore could be at the origin of the velocity dispersion profile. We further investigate cases with varying mass-to-light ratios and find that a compact dark stellar component can also explain our observations. However, such a component would closely resemble the black hole from the constant mass-to-light ratio models as this component must be confined to the central ~5arcsec of the cluster and must have a similar mass. Independent constraints on the distribution of stellar remnants in the cluster or kinematic measurements at the highest possible spatial resolution should be able to distinguish the two alternatives.

Since 2014, MUSE, the Multi-Unit Spectroscopic Explorer, is in operation at the ESO-VLT. It combines a superb spatial sampling with a large wavelength coverage. By design, MUSE is an integral-field instrument, but its field-of-view and large multiplex make it a powerful tool for multi-object spectroscopy too. Every data-cube consists of 90,000 image-sliced spectra and 3700 monochromatic images. In autumn 2014, the observing programs with MUSE have commenced, with targets ranging from distant galaxies in the Hubble Deep Field to local stellar populations, star formation regions and globular clusters. This paper provides a brief summary of the key features of the MUSE instrument and its complex data reduction software. Some selected examples are given, how multi-object spectroscopy for hundreds of continuum and emission-line objects can be obtained in wide, deep and crowded fields with MUSE, without the classical need for any target pre-selection.

We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3.

We present a new integral-field spectroscopic dataset of the central part of the Orion Nebula (M 42), observed with the MUSE instrument at the ESO VLT. We reduced the data with the public MUSE pipeline. The output products are two FITS cubes with a spatial size of ~5.9'x4.9' (corresponding to ~0.76 pc x 0.63 pc) and a contiguous wavelength coverage of 4595...9366 Angstrom, spatially sampled at 0.2". We provide two versions with a sampling of 1.25 Angstrom and 0.85 Angstrom in dispersion direction. Together with variance cubes these files have a size of 75 and 110 GiB on disk. They represent one of the largest integral field mosaics to date in terms of information content. We make them available for use in the community. To validate this dataset, we compare world coordinates, reconstructed magnitudes, velocities, and absolute and relative emission line fluxes to the literature and find excellent agreement. We derive a two-dimensional map of extinction and present de-reddened flux maps of several individual emission lines and of diagnostic line ratios. We estimate physical properties of the Orion Nebula, using the emission line ratios [N II] and [S III] (for the electron temperature $T_e$) and [S II] and [Cl III] (for the electron density $N_e$), and show two-dimensional images of the velocity measured from several bright emission lines.

MUSE observations of NGC5813 reveal a complex structure in the velocity dispersion map, previously hinted by SAURON observations. The structure is reminiscent of velocity dispersion maps of galaxies comprising two counter-rotating discs, and may explain the existence of the kinematically distinct core (KDC). Further evidence for two counter-rotating components comes from the analysis of the higher moments of the stellar line-of-sight velocity distributions and fitting MUSE spectra with two separate Gaussian line-of-sight velocity distributions. The emission-line kinematics show evidence of being linked to the present cooling flows and the buoyant cavities seen in X-rays. We detect ionised gas in a nuclear disc-like structure, oriented like the KDC, which is, however, not directly related to the KDC. We build an axisymmetric Schwarzschild dynamical model, which shows that the MUSE kinematics can be reproduced well with two counter-rotating orbit families, characterised by relatively low angular momentum components, but clearly separated in integral phase space and with radially varying contributions. The model indicates that the counter-rotating components in NGC5813 are not thin discs, but dynamically hot structures. Our findings give further evidence that KDCs in massive galaxies should not necessarily be considered as structurally or dynamically decoupled regions, but as the outcomes of the mixing of different orbital families, where the balance in the distribution of mass of the orbital families is crucial. We discuss the formation of the KDC in NGC5813 within the framework of gas accretion, binary mergers and formation of turbulent thick discs from cold streams at high redshift.

Diffuse Interstellar Bands (DIBs) are non-stellar weak absorption features of unknown origin found in the spectra of stars viewed through one or several clouds of Interstellar Medium (ISM). Research of DIBs outside the Milky Way is currently very limited. Specifically spatially resolved investigations of DIBs outside of the Local Group is, to our knowledge, inexistent. Here, we explore the capability of the high sensitivity Integral Field Spectrograph, MUSE, as a tool to map diffuse interstellar bands at distances larger than 100 Mpc. We use MUSE commissioning data for AM 1353-272 B, the member with highest extinction of the "The Dentist's Chair", an interacting system of two spiral galaxies. High signal-to-noise spectra were created by co-adding the signal of many spatial elements distributed in a geometry of concentric elliptical half-rings. We derived decreasing radial profiles for the equivalent width of the $\lambda$5780.5 DIB both in the receding and approaching side of the companion galaxy up to distances of $\sim$4.6 kpc from the center of the galaxy. Likewise, interstellar extinction, as derived from the Halpha/Hbeta line ratio displays a similar trend, with decreasing values towards the external parts. This translates into an intrinsic correlation between the strength of the DIB and the extinction within AM 1353-272 B consistent with the current existing global trend between these quantities when using measurements for both Galactic and extragalactic sight lines. Mapping of DIB strength in the Local Universe as up to now only done for the Milky Way seems feasible. This offers a new approach to study the relationship between DIBs and other characteristics and species of the ISM in different conditions as those found in our Galaxy to the use of galaxies in the Local Group and/or single sightlines towards supernovae, quasars and galaxies outside the Local Group.

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

We used the PMAS integral field spectrograph to obtain large sets of radial velocities in the central regions of three northern Galactic globular clusters: M3, M13, and M92. By applying the novel technique of crowded field 3D spectroscopy, we measured radial velocities for about 80 stars within the central ~ 10 arcsec of each cluster. These are by far the largest spectroscopic datasets obtained in the innermost parts of these clusters up to now. To obtain kinematical data across the whole extent of the clusters, we complement our data with measurements available in the literature. We combine our velocity measurements with surface brightness profiles to analyse the internal dynamics of each cluster using spherical Jeans models, and investigate whether our data provide evidence for an intermediate-mass black hole in any of the clusters. The surface brightness profiles reveal that all three clusters are consistent with a core profile, although shallow cusps cannot be excluded. We find that spherical Jeans models with a constant mass-to-light ratio provide a good overall representation of the kinematical data. A massive black hole is required in none of the three clusters to explain the observed kinematics. Our 1sigma (3sigma) upper limits are 5300 M_sun (12000 M_sun) for M3, 8600 M_sun (13000 M_sun) for M13, and 980 M_sun (2700 M_sun) for M92. A puzzling circumstance is the existence of several potential high velocity stars in M3 and M13, as their presence can account for the majority of the discrepancies that we find in our mass limits compared to M92.

(Abridged) We describe a new method to extract spectra of stars from observations of crowded stellar fields with integral field spectroscopy (IFS). Our approach extends the well-established concept of crowded field photometry in images into the domain of 3-dimensional spectroscopic datacubes. The main features of our algorithm are: (1) We assume that a high-fidelity input source catalogue already exists and that it is not needed to perform sophisticated source detection in the IFS data. (2) Source positions and properties of the point spread function (PSF) vary smoothly between spectral layers of the datacube, and these variations can be described by simple fitting functions. (3) The shape of the PSF can be adequately described by an analytical function. Even without isolated PSF calibrator stars we can therefore estimate the PSF by a model fit to the full ensemble of stars visible within the field of view. (4) By using sparse matrices to describe the sources, the problem of extracting the spectra of many stars simultaneously becomes computationally tractable. We present extensive performance and validation tests of our algorithm using realistic simulated datacubes that closely reproduce actual IFS observations of the central regions of Galactic globular clusters. We investigate the quality of the extracted spectra under the effects of crowding. The main effect of blending between two nearby stars is a decrease in the S/N in their spectra. The effect increases with the crowding in the field in a way that the maximum number of stars with useful spectra is always ~0.2 per spatial resolution element. This balance breaks down when exceeding a total source density of ~1 significantly detected star per resolution element. We close with an outlook by applying our method to a simulated globular cluster observation with the upcoming MUSE instrument at the ESO-VLT.

[Abridged] Detecting cosmic ray hits (cosmics) in fiber-fed IFS data of single exposures is a challenging task, because of the complex signal recorded by IFS instruments. Existing detection algorithms are commonly found to be unreliable in the case of IFS data and the optimal parameter settings are usually unknown a-priori for a given dataset. The CALIFA survey generates hundreds of IFS datasets for which a reliable and robust detection algorithm for cosmics is required as an important part of the fully automatic CALIFA data reduction pipeline. We developed a novel algorithm, PyCosmic, which combines the edge-detection algorithm of L.A.Cosmic with a point-spread function convolution scheme. We generated mock data to compute the efficiency of different algorithms for a wide range of characteristic fibre-fed IFS datasets using the PMAS and VIMOS IFS instruments as representative cases. PyCosmic is the only algorithm that achieves an acceptable detection performance for CALIFA data. We find that PyCosmic is the most robust tool with a detection rate of >~90% and a false detection rate <5% for any of the tested IFS data. It has one less free parameter than the L.A.Cosmic algorithm. Only for strongly undersampled IFS data does L.A.Cosmic exceed the performance of PyCosmic by a few per cent. DCR never reaches the efficiency of the other two algorithms and should only be used if computational speed is a concern. Thus, PyCosmic appears to be the most versatile cosmics detection algorithm for IFS data. It is implemented in the new CALIFA data reduction pipeline as well as in recent versions of the multi-instrument IFS pipeline P3D.

We present a new library of synthetic spectra based on the stellar atmosphere code PHOENIX. It covers the wavelength range from 500{\AA} to 55000{\AA} with a resolution of R=500000 in the optical and near IR, R=100000 in the IR and {\Delta}{\lambda}=0.1{\AA} in the UV. The parameter space covers 2300K<=Teff<=8000K, 0.0<=log(g)<=6.0, -4.0<=[Fe/H]<=+1.0 and -0.3<=[{\alpha}/Fe]<=+0.8. The library is work-in-progress and going to be extended to at least Teff=25000K. We use a new self-consistent way of describing the microturbulence for our model atmospheres. The entire library of synthetic spectra will be available for download. Futhermore we present a method for fitting spectra, especially designed to work with the new 2nd generation VLT instrument MUSE. We show that we can determine stellar parameters (Teff, log(g), [Fe/H] and [{\alpha}/Fe]) and even single element abundances.