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.

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.

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.

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.