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.

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.

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.

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.

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 present a study of four different formation scenarios that may be able to explain the origin of the recently announced planet (`Proxima b') orbiting Proxima Centauri. The aim is to examine how the formation scenarios differ in their predictions for the multiplicity of the Proxima system, the water/volatile content of Proxima b and its orbital eccentricity, so that these can be tested by future observations. A scenario of in situ formation via giant impacts from a locally enhanced disc of planetary embryos/planetesimals predicts that Proxima b will be in a multiplanet system with a measurably finite eccentricity. Assuming that the local solid enhancement needed to form a Proxima b analogue arises because of the inwards drift of solids in the form of small planetesimals/boulders, this scenario also results in Proxima b analogues being only moderately endowed with water/volatiles. A scenario in which multiple embryos form, migrate and mutually collide within a gas disc results in Proxima b being a member of a multiple system, possibly displaying mean motion resonances, but where the constituent members are Ocean planets due to accretion occurring mainly outside of the snowline. A scenario in which a single accreting embryo forms outside the snowline, and migrates inwards while accreting planetesimals/pebbles results in Proxima b being an isolated Ocean planet on a circular orbit. A scenario in which Proxima b formed via pebble accretion interior to the snowline produces a dry planet on a circular orbit. Future observations that characterise the physical and orbital properties of Proxima b, and any additional planets in the system, will provide valuable insights into the formation history of this neighbouring planetary system.

At a distance of 1.295 parsecs, the red-dwarf Proxima Centauri ($\alpha$ Centauri C, GL 551, HIP 70890, or simply Proxima) is the Sun's closest stellar neighbor and one of the best studied low-mass stars. It has an effective temperature of only $\sim$ 3050 K, a luminosity of $\sim$0.1 per cent solar, a measured radius of 0.14 R$_\odot$ and a mass of about 12 per cent the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is $\sim$ 83 days, and its quiescent activity levels and X-ray luminosity are comparable to the Sun's. New observations reveal the presence of a small planet orbiting Proxima with a minimum mass of 1.3~Earth masses and an orbital period of $\sim$11.2 days. Its orbital semi-major axis is $\sim0.05$ AU, with an equilibrium temperature in the range where water could be liquid on its surface.