Search

We present the design, rationale, properties and catalogs of the MusE Gas FLOw and Wind survey (MEGAFLOW), a survey of the cool gaseous halos of $z\sim1$ galaxies using low-ionization MgII absorption systems. The survey consists of 22 quasar fields selected from the Sloan Digital Sky Survey (SDSS) having multiple ($\geq3$) strong MgII absorption lines over the redshift range $0.31$, i.e. more than 60 galaxies per arcmin$^{2}$. They were identified using a dual detection algorithm based on both continuum and emission line objects. The achieved [OII] 50\%\ completeness is 3.7$\times 10^{-18}$ erg/s/cm$^2$ (corresponding to SFR$>0.01$ M$_\odot$ yr$^{-1}$ at $z=1$) using realistic mock [OII] emitters and the 50\%\ completeness is $m_{F775W}\approx26$ AB magnitudes for continuum sources. We find that (i) the fraction of [OII] emitting galaxies which have no continuum is $\sim15$\%; (ii) the success rate in identifying at least one galaxy within 500 km/s and 100 kpc is $\approx90$\%\ for MgII absorptions with $W_r^{2796}\gtrsim0.5$~\AA; (iii) the mean number of galaxies per MgII absorption is $2.9\pm1.6$ within the MUSE field-of-view; (iv) of the 80 MgII systems at $0.310^9$ M$_\odot$, and star-formation rates $>1$ M$_\odot$ yr$^{-1}$., Comment: 23 pages; 16 figures; in press for publication in A&A; accepted version

After over three decades of unsuccessful attempts, we report the first detection of molecular gas emission in Malin 1, the largest spiral galaxy observed to date, and one of the most iconic giant low surface brightness galaxies. Using ALMA, we detect significant $^{12}$CO(J=1-0) emission in the galaxy's central region and tentatively identify CO emission across three regions on the disc. These observations allow for a better estimate of the H$_2$ mass and molecular gas mass surface density, both of which are remarkably low given the galaxy's scale. By integrating data on its HI mass, we derive a very low molecular-to-atomic gas mass ratio. Overall, our results highlight the minimal presence of molecular gas in Malin 1, contrasting sharply with its extensive, homogeneous atomic gas reservoir. For the first time, we position Malin 1 on the Kennicutt-Schmidt (K-S) diagram, where it falls below the main sequence for normal spirals, consistent with previous upper limits but now with more accurate figures. These findings are crucial for constraining our understanding of star formation processes in environments characterized by extremely low molecular gas densities and for refining models of galaxy formation, thereby improving predictions concerning the formation, evolution, and distribution of these giant, elusive galaxies., Comment: Accepted for publication in The Astrophysical Journal Letters. Published

Melnick 34 (Mk 34) is one of the most massive binary systems known and is one of the brightest X-ray point sources in the 30 Doradus region. We investigated the impact of this massive system on the surrounding interstellar medium (ISM) using the optical spectroscopic capabilities of the narrow-field mode (NFM) of the Multi-Unit Spectroscopic Explorer (MUSE). MUSE-NFM spatially resolved the ISM in the vicinity of Mk 34 with a resolution comparable to that of the HST. The analysis of the [NII]$\lambda$6583 and [SII]$\lambda$6717 emission lines reveals a cone-like structure apparently originating from Mk 34 and extending southeast. Electron density maps and radial velocity measurements of the ISM lines further support an outflow scenario traced by these emissions. While no clear northwestern counterpart to this outflow was observed, we note increased extinction in that direction, towards the R136 cluster. The ISM material along the projected diagonal of the outflow on both sides of Mk 34 shows similar properties in terms of the emission line ratios seen in the Baldwin-Phillips-Terlevich diagram. These results are consistent across two observational epochs. Additionally, we examined the residual maps within a 0.5" radius of Mk 34 after modeling and subtracting the point spread function. The observed variations in the residuals could potentially be linked to Mk 34's known periodic behavior. However, further observations with appropriate cadence are needed to fully monitor the 155 day periodicity of Mk 34's X-ray emissions., Comment: 9 pages, 10 figures. Accepted in A&A

BlueMUSE is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph under development for the Very Large Telescope (VLT). With an optimised transmission down to 350 nm, spectral resolution of R$\sim$3500 on average across the wavelength range, and a large FoV (1 arcmin$^2$), BlueMUSE will open up a new range of galactic and extragalactic science cases facilitated by its specific capabilities. The BlueMUSE consortium includes 9 institutes located in 7 countries and is led by the Centre de Recherche Astrophysique de Lyon (CRAL). The BlueMUSE project development is currently in Phase A, with an expected first light at the VLT in 2031. We introduce here the Top Level Requirements (TLRs) derived from the main science cases, and then present an overview of the BlueMUSE system and its subsystems fulfilling these TLRs. We specifically emphasize the tradeoffs that are made and the key distinctions compared to the MUSE instrument, upon which the system architecture is built., Comment: 20 pages, 10 figures, proceedings of the SPIE astronomical telescopes and instrumentation conference, Yokohama, 16-21 June

BlueMUSE is a blue, medium spectral resolution, panoramic integral-field spectrograph under development for the Very Large Telescope (VLT). We demonstrate and discuss an early End-To-End simulation software for final BlueMUSE datacube products. Early access to such simulations is key to a number of aspects already in the development stage of a new major instrument. We outline the software design choices, including lessons learned from the MUSE instrument in operation at the VLT since 2014. The current simulation software package is utilized to evaluate some of the technical specifications of BlueMUSE as well as giving assistance in the assessment of certain trade offs regarding instrument capabilities, e.g., spatial and spectral resolution and sampling. By providing simulations of the end-user product including realistic environmental conditions such as sky contamination and seeing, BlueSi can be used to devise and prepare the science of the instrument by individual research teams., Comment: 20 pages, 13 figures

Aims. The central region of the Giant Low Surface Brightness galaxy Malin 1 has long been known to have a complex morphology with evidence of a bulge, disc, and potentially a bar hosting asymmetric star formation. In this work, we use VLT/MUSE data to resolve the central region of Malin 1 in order to determine its structure. Methods. We use careful light profile fitting in every image slice of the datacube to create wavelength-dependent models of each morphological component, from which we could cleanly extract their spectra. We then used the kinematics and emission line properties from these spectra to better understand the nature of each component extracted from our model fit. Results. We report the detection of a pair of distinct sources at the centre of this galaxy with a separation of ~1.05", which corresponds to a separation on sky of ~1.9 kpc. The radial velocity data of each object confirms that they both lie in the kinematic core of the galaxy, and analysis of the emission lines reveals that the central compact source is more consistent with being ionized by star formation and/or a LINER, while the off-centre compact source lies closer to the separation between star-forming galaxies and AGN. Conclusions. This evidence suggests that the centre of Malin 1 hosts either a bar with asymmetric star formation or two distinct components in which the off-centre compact source could either be a star-forming clump containing one or more star clusters that is in the process of falling into the core of the galaxy and which will eventually merge with the central NSC, or a clump of gas infalling into the centre of the galaxy from either outside or from the disc and triggering star formation there., Comment: 13 pages, 5 figures, accepted for publication in A&A

Star formation in the early Universe has left its imprint on the chemistry of observable stars in galaxies. We derive elemental abundances and the slope of the low-mass end of the initial mass function (IMF) for a sample of 25 very massive galaxies, separated into brightest cluster galaxies (BCGs) and their massive satellites. The elemental abundances of BGCs and their satellites are similar, but for some elements, satellite galaxies show a correlation with the global velocity dispersion. Using a subset of derived elemental abundances, we model the star formation histories of these galaxies with chemical evolution models, and predict the high-mass end slope of the IMF and star formation timescales. The high-mass end IMF slope of the satellite galaxies correlates with the global velocity dispersion. The low- and the high-mass end IMF slopes are weakly correlated in a general sense that top heavy IMFs are paired with bottom heavy IMFs. Our results do not necessarily imply that the IMF was simultaneously bottom and top heavy. Instead, our findings can be considered consistent with a temporal variation in the IMF, where, for massive galaxies, the high-mass end IMF slope is representative of the very early age and the low-mass end slope of the later star formation. The small but noticeable differences between the BCGs and the satellites in terms of their elemental abundances and IMF slopes, together with their stellar kinematical properties, suggest somewhat different formation pathways, where BCGs experience more major, gas-free mergers., Comment: Accepted for publication in MNRAS, 27 pages, 18 figures, 5 tables (including Appendices and supplementary material)

The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate, Comment: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

$\omega$ Centauri is considered the most massive globular cluster of the Milky Way and likely the former nuclear star cluster of a galaxy accreted by the Milky Way. It is speculated to contain an intermediate-mass black hole (IMBH) from several dynamical models. However, uncertainties regarding the location of the cluster center or the retention of stellar remnants limit the robustness of the IMBH detections reported so far. In this paper, we derive and study the stellar kinematics from the highest-resolution spectroscopic data yet, using the Multi Unit Spectroscopic Explorer (MUSE) in the narrow field mode (NFM) and wide field mode (WFM). Our exceptional data near the center reveal for the first time that stars within the inner 20" ($\sim$0.5 pc) counter-rotate relative to the bulk rotation of the cluster. Using this dataset, we measure the rotation and line-of-sight velocity dispersion (LOSVD) profile out to 120$''$ with different centers proposed in the literature. We find that the velocity dispersion profiles using different centers match well with those previously published. Based on the counter--rotation, we determine a kinematic center and look for any signs of an IMBH using the high-velocity stars close to the center. We do not find any significant outliers $>$60 km/s within the central 20$''$, consistent with no IMBH being present at the center of $\omega$ Centauri. A detailed analysis of Jeans' modeling of the putative IMBH will be presented in the next paper of the series., Comment: Accepted for publication in MNRAS, 16 pages, 9 figures

We introduce the MUSE gAlaxy Groups in COSMOS (MAGIC) survey, which was built to study the impact of environment on galaxy evolution over the last 8 Gyr. It consists of 17 MUSE fields targeting 14 massive structures at intermediate redshift ($0.321.5$. The spectroscopic redshift completeness is high: in the redshift range of [OII] emitters ($0.25 \le z < 1.5$), where most of the groups are found, it globally reaches a maximum of 80% down to $z_{app}^{++}=25.9$, and locally decreases from $\sim 100$% to $\sim50$% in magnitude bins from $z_{app}^{++}=23-24$ to $z_{app}^{++}=25.5$. We find that the fraction of quiescent galaxies increases with local density and with the time spent in groups. A morphological dichotomy is also found between bulge-dominated quiescent and disk-dominated star-forming galaxies. As environment gets denser, the peak of the stellar mass distribution shifts towards $M_*>10^{10}~M_\odot$, and the fraction of galaxies with $M_*<10^9~M_\odot$ decreases significantly, even for star-forming galaxies. We also highlight peculiar features such as close groups, extended nebulae, and a gravitational arc. Our results suggest that galaxies are preprocessed in groups of increasing mass before entering rich groups and clusters. We publicly release two catalogs containing the properties of galaxies and groups, respectively., Comment: 27 pages, 22 figures, accepted for publication in A&A

A quarter of a century has passed since the observing technique of integral field spectroscopy (IFS) was first applied to planetary nebulae (PNe). Progress after the early experiments was relatively slow, mainly because of the limited field-of-view (FoV) of first generation instruments.With the advent of MUSE at the ESO Very Large Telescope, this situation has changed. MUSE is a wide field-of-view, high angular resolution, one-octave spanning optical integral field spectrograph with high throughput. Its major science mission has enabled an unprecedented sensitive search for Ly{\alpha} emitting galaxies at redshift up to z=6.5. This unique property can be utilized for faint objects at low redshift as well. It has been demonstrated that MUSE is an ideal instrument to detect and measure extragalactic PNe with high photometric accuracy down to very faint magnitudes out to distances of 30 Mpc, even within high surface brightness regions of their host galaxies. When coupled with a differential emission line filtering (DELF) technique, MUSE becomes far superior to conventional narrow-band imaging, and therefore MUSE is ideal for accurate Planetary Nebula Luminosity Function (PNLF) distance determinations. MUSE enables the PNLF to become a competitive tool for an independent measure of the Hubble constant, and stellar population studies of the host galaxies that present a sufficiently large number of PNe., Comment: 12 pages, 8 figures, to appear in proceedings of IAU Symposium 384 "Planetary Nebulae", Krakow, September 4-8, 2023

Giant low-surface brightness (GLSB) galaxies are an extreme class of objects with very faint and extended gas-rich disks. Malin 1 is the largest GLSB galaxy known to date, but its formation is still poorly understood. We use VLT/MUSE IFU spectroscopic observations of Malin 1 to reveal, for the first time, the presence of H$\alpha$ emission distributed across numerous regions along its disk, up to radial distances of $\sim$100 kpc. We made an estimate of the dust attenuation using the Balmer decrement and found that Malin 1 has a mean H$\alpha$ attenuation of 0.36 mag. We observe a steep decline in the star formation rate surface density ($\Sigma_{\rm SFR}$) within the inner 20 kpc, followed by a shallow decline in the extended disk. Similarly, the gas phase metallicity we estimated shows a steep gradient in the inner 20 kpc, followed by a flattening of the metallicity in the extended disk with a relatively high value of $\sim$0.6 $Z_{\odot}$. We found that the normalized abundance gradient of the inner disk is similar to values found in normal galaxies but with an extreme value in the extended disk. A comparison of the star formation rate surface density and gas surface density shows that, unlike normal disk galaxies or other LSBs, Malin 1 exhibits a very low star formation efficiency. Owing to the detection of emission lines over a large part of the disk of Malin 1, this work sheds light on the star formation processes in this unique galaxy, highlighting its extended star-forming disk, dust attenuation, almost flat metallicity distribution in the outer disk, and exceptionally low star-formation efficiency. Our findings contribute to a more detailed understanding of the formation of the giant disk of Malin 1 and also constrain possible proposed scenarios on the nature of GLSB galaxies in general., Comment: 12 pages, 9 figures, accepted for publication in A&A

Thanks to the MUSE integral field spectrograph on the VLT, extragalactic distance measurements with the [O III] 5007 A planetary nebula luminosity function (PNLF) are now possible out to approx. 40 Mpc. Here we analyze the VLT/MUSE data for 20 galaxies from the ESO public archive to identify the systems' planetary nebulae (PNe) and determine their PNLF distances. Three of the galaxies do not contain enough PNe for a robust measure of the PNLF, and the results for one other system are compromised by the galaxy's internal extinction. However, we obtain robust PNLF distances for the remaining 16 galaxies, two of which are isolated and beyond 30 Mpc in a relatively unperturbed Hubble flow. From these data, we derive a Hubble Constant of 74.2 +/- 7.2 (stat) +/-3.7 (sys) km/s/Mpc, a value that is very similar to that found from other quality indicators (e.g., Cepheids, the tip of the red giant branch, and surface brightness fluctuations). At present, the uncertainty is dominated by the small number of suitable galaxies in the ESO archival and their less than ideal observing conditions and calibrations. Based on our experience with these systems, we identify the observational requirements necessary for the PNLF to yield a competitive value for H0 that is independent of the SN Ia distance scale, and help resolve the current tension in the Hubble constant., Comment: 60 pages, 65 figures, published at ApJS 271, 40. DOI 10.3847/1538-4365/ad2166

Tidal dwarf galaxies (TDGs) are low-mass objects that form within tidal and/or collisional debris ejected from more massive interacting galaxies. We use CO($1-0$) observations from ALMA and integral-field spectroscopy from MUSE to study molecular and ionized gas in three TDGs: two around the collisional galaxy NGC 5291 and one in the late-stage merger NGC 7252. The CO and H$\alpha$ emission is more compact than the HI emission and displaced from the HI dynamical center, so these gas phases cannot be used to study the internal dynamics of TDGs. We use CO, HI, and H$\alpha$ data to measure the surface densities of molecular gas ($\Sigma_{\rm mol}$), atomic gas ($\Sigma_{\rm atom}$) and star-formation rate ($\Sigma_{\rm SFR}$), respectively. We confirm that TDGs follow the same spatially integrated $\Sigma_{\rm SFR}-\Sigma_{\rm gas}$ relation of regular galaxies, where $\Sigma_{\rm gas} = \Sigma_{\rm mol} + \Sigma_{\rm atom}$, even though they are HI dominated. We find a more complex behaviour in terms of the spatially resolved $\Sigma_{\rm SFR}-\Sigma_{\rm mol}$ relation on sub-kpc scales. The majority ($\sim$60$\%$) of SF regions in TDGs lie on the same $\Sigma_{\rm SFR}-\Sigma_{\rm mol}$ relation of normal spiral galaxies but show a higher dispersion around the mean. The remaining fraction of SF regions ($\sim$40$\%$) lie in the starburst region and are associated with the formation of massive super star clusters, as shown by Hubble Space Telescope images. We conclude that the local SF activity in TDGs proceeds in a hybrid fashion, with some regions comparable to normal spiral galaxies and others to extreme starbursts., Comment: 11 pages, 6 figures, 3 tables, Accepted for publication in MNRAS

Leo T ($M_V = -8.0$) is both the faintest and the least massive galaxy known to contain neutral gas and to display signs of recent star formation. We analyse photometry and stellar spectra to identify member stars and to better understand the overall dynamics and stellar content of the galaxy and to compare the properties of its young and old stars. We use data from the Multi Unit Spectroscopic Explorer (MUSE) on the VLT. We supplement this information with spectroscopic data from the literature and with Hubble Space Telescope (HST) photometry. Our analysis reveals two distinct populations of stars in Leo T. The first population, with an age of $\lesssim 500~\mathrm{Myr}$, includes three emission-line Be stars comprising 15% of the total number of young stars. The second population of stars is much older, with ages ranging from $>5~\mathrm{Gyr}$ to as high as $10~\mathrm{Gyr}$. We combine MUSE data with literature data to obtain an overall velocity dispersion of $\sigma_{v} = 7.07^{+1.29}_{-1.12}~\mathrm{km\ s^{-1}}$ for Leo T. When we divide the sample of stars into young and old populations, we find that they have distinct kinematics. Specifically, the young population has a velocity dispersion of $2.31^{+2.68}_{-1.65}\,\mathrm{km\ s^{-1}}$, contrasting with that of the old population, of $8.14^{+1.66}_{-1.38}\,\mathrm{km\ s^{-1}}$. The fact that the kinematics of the cold neutral gas is in good agreement with the kinematics of the young population suggests that the recent star formation in Leo T is linked with the cold neutral gas. We assess the existence of extended emission-line regions and find none to a surface brightness limit of~$< 1\times 10^{-20}\,\mathrm{erg}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}~\mathrm{arcsec}^{-2}$ which corresponds to an upper limit on star formation of $\sim 10^{-11}~\mathrm{M_\odot~yr^{-1}~pc^{-2}}$, implying that the star formation in Leo T has ended., Comment: Accepted for publication in A&A on 04-08-2023. 17 pages, 17 figures, 3 tables

[ABRIDGED] Stellar feedback in high-redshift galaxies plays an important role in the re-ionization epoch of the Universe. Green Pea galaxies (GPs) postulate as favorite local laboratories. However, at their typical redshift of $z\sim0.2$, the most intimate interaction between stars and surrounding ISM cannot be disentangled. Detailed studies of Blue Compact Dwarf galaxies (BCDs) are necessary to anchor our investigations on them. We present here a study in detail UM 462, a BCD with similar properties to GPs uisng high quality optical IFS data with MUSE. Total oxygen abundance by means of the direct method is 12+$\log$(O/H)$\sim$8.02 and homogenous all over the galaxy, in stark contrast with the metallicities derived from several strong line methods. The velocity field for the ionised gas presents a velocity stratification in the area towards the north with redder velocities in the high ionisation lines and bluer velocities in the low ionisation lines. This is the only area with velocity dispersions clearly above the MUSE instrumental width, and it is surrounded by two $\sim$1 kpc-long structures nicknamed \emph{the horns}. We interpret the observational evidence in that area as a fragmented super-bubble fruit of the stellar feedback and it may constitute a preferred channel for LyC photons from the youngest generation of stars to escape. The most recent SF seems to propagate from the outer to the inner parts of the galaxy, and then from east to west. We identified a supernova remnant and Wolf-Rayet stars - as traced by the red bump - that support this picture. The direction of the propagation implies the presence of younger Wolf-Rayet stars at the maximum in H$\alpha$. The ensemble of results exemplifies the potential of 2D detailed spectroscopic studies of dwarf star-forming galaxies at high spatial resolution as key reference for similar studies on primeval galaxies., Comment: 29 pages, 25 figures; accepted by A&A, 03/05/2023

The formation process of multiple populations in globular clusters is still up for debate. Kinematic differences between the populations are particularly interesting in this respect, because they allow us to distinguish between single-epoch formation scenarios and multi-epoch formation scenarios. We analyze the kinematics of 25 globular clusters and aim to find kinematic differences between multiple populations to constrain their formation process. We split red-giant branch (RGB) stars in each cluster into three populations (P1, P2, P3) for the type-II clusters and two populations (P1 and P2) otherwise using Hubble photometry. We derive the rotation and dispersion profiles for each cluster and its populations by using all stars with radial velocity measurements obtained from MUSE spectroscopy. Based on these profiles, we calculate the rotation strength in terms of ordered-over-random motion $\left(v/\sigma\right)_\mathrm{HL}$ evaluated at the half-light radius of the cluster. We detect rotation in all but four clusters. For NGC~104, NGC~1851, NGC~2808, NGC~5286, NGC~5904, NGC~6093, NGC~6388, NGC~6541, NGC~7078 and NGC~7089 we also detect rotation for P1 and/or P2 stars. For NGC~2808, NGC~6093 and NGC~7078 we find differences in $\left(v/\sigma\right)_\mathrm{HL}$ between P1 and P2 that are larger than $1\sigma$. Whereas we find that P2 rotates faster than P1 for NGC~6093 and NGC~7078, the opposite is true for NGC~2808. However, even for these three clusters, the differences are still of low significance. We find that the strength of rotation of a cluster generally scales with its median relaxation time. For P1 and P2, the corresponding relation is very weak at best. We observe no correlation between the difference in rotation strength between P1 and P2 and cluster relaxation time. The MUSE stellar radial velocities that this analysis is based on are made publicly available.

We perform a deep survey of planetary nebulae (PNe) in the spiral galaxy NGC 300 to construct its planetary nebula luminosity function (PNLF). We aim to derive the distance using the PNLF and to probe the characteristics of the most luminous PNe. We analyse 44 fields observed with MUSE at the VLT, covering a total area of $\sim11$ kpc$^2$. We find [OIII]5007 sources using the differential emission line filter (DELF) technique. We identify PNe through spectral classification using the aid of the BPT-diagram. The PNLF distance is derived using the maximum likelihood estimation technique. For the more luminous PNe, we also measure their extinction using the Balmer decrement. We estimate the luminosity and effective temperature of the central stars of the luminous PNe, based on estimates of the excitation class and the assumption of optically thick nebulae. We identify 107 PNe and derive a most-likely distance modulus $(m-M)_0 = 26.48^{+0.11}_{-0.26}$ ($d = 1.98^{+0.10}_{-0.23}$ Mpc). We find that the PNe at the PNLF cut-off exhibit relatively low extinction, with some high extinction cases caused by local dust lanes. We present the lower limit luminosities and effective temperatures of the central stars for some of the brighter PNe. We also identify a few Type I PNe that come from a young population with progenitor masses $>2.5 \, M_\odot$, however do not populate the PNLF cut-off. The spatial resolution and spectral information of MUSE allow precise PN classification and photometry. These capabilities also enable us to resolve possible contamination by diffuse gas and dust, improving the accuracy of the PNLF distance to NGC 300., Comment: 20 pages, 11 figures. Accepted for publication in A&A

Stellar winds of cool and pulsating asymptotic giant branch (AGB) stars enrich the interstellar medium with large amounts of processed elements and various types of dust. We present the first study on the influence of gas-to-dust drift on ab initio simulations of stellar winds of M-type stars driven by radiation pressure on forsterite particles. Our study is based on our radiation hydrodynamic model code T-800 that includes frequency-dependent radiative transfer, dust extinction based on Mie scattering, grain growth and ablation, gas-to-dust drift using one mean grain size, a piston that simulates stellar pulsations, and an accurate high spatial resolution numerical scheme. To enable this study, we calculated new gas opacities based on the ExoMol database, and we extended the model code to handle the formation of minerals that may form in M-type stars. We determine the effects of drift by comparing drift models to our new and extant non-drift models. Three out of four new drift models show high drift velocities, 87-310 km/s. Our new drift model mass-loss rates are 1.7-13 per cent of the corresponding values of our non-drift models, but compared to the results of two extant non-drift models that use the same stellar parameters, these same values are 0.33-1.5 per cent. Meanwhile, a comparison of other properties such as the expansion velocity and grain size show similar values. Our results, which are based on single-component forsterite particles, show that the inclusion of gas-to-drift is of fundamental importance in stellar wind models driven by such transparent grains. Assuming that the drift velocity is insignificant, properties such as the mass-loss rate may be off from more realistic values by a factor of 50 or more., Comment: 15 pages, 6 figures, accepted and in press

We present the second data release of the MUSE Hubble UDF surveys, which includes the deepest spectroscopic survey ever performed. The MUSE data, with their 3D content, amazing depth, wide spectral range, and excellent spatial and medium spectral resolution, are rich in information. This update of the first release incorporates a new 141-hour adaptive-optics-assisted MXDF field (1' diameter FoV) in addition to the reprocessed 10-hour mosaic (3'x3') and the single 31-hour deep field (1'x1'). We have securely identified and measured the redshift of 2221 sources, an increase of 41% compared to the first release. With the exception of 8 stars, the collected sample consists of 25 nearby galaxies (z < 0.25), 677 OII emitters (z=0.25-1.5), 201 galaxies in the MUSE redshift desert range (z=1.5-2.8), and 1308 LAEs (z=2.8-6.7). This represents an order of magnitude more redshifts than the collection of all spectroscopic redshifts obtained before MUSE in the Hubble UDF area (2221 vs 292). At z > 3, the difference is even more striking, with a factor of 65 increase (1308 vs 20). We compared the measured redshifts against three published photometric redshift catalogs and find the photo-z accuracy to be lower than the constraints provided by photo-z fitting codes. 80% of the galaxies have an HST counterpart. They are on average faint, with a median magnitude of 25.7 and 28.7 for the OII and Ly-alpha emitters, respectively. SED fits show that these galaxies tend to be low-mass star-forming galaxies, with a median stellar mass of 6.2 10**8 M and a median SFR of 0.4 M/yr. 20% of our catalog, or 424 galaxies, have no HST counterpart. The vast majority of these new sources are high EQW z>2.8 LAEs that are detected by MUSE thanks to their bright and asymmetric broad Ly-alpha line. We release advanced data products, specific software, and a web interface to select and download data sets., Comment: 46 pages, 48 figures

There are known differences between the physical properties of HII and diffuse ionized gas (DIG), but most of the studied regions in the literature are relatively bright. We compiled a faint sample of 390 HII regions with median $\log_{10}H\alpha$=34.7 in the spiral galaxy NGC300, derived their physical properties in terms of metallicity, density, extinction, and kinematics, and performed a comparative analysis of the properties of the DIG. We used MUSE data of nine fields in NGC300, covering a galactocentric distance of zero to ~450 arcsec (~4 projected kpc), including spiral arm and inter-arm regions. We binned the data in dendrogram leaves and extracted all strong nebular emission lines. We identified HII and DIG regions and compared their electron densities, metallicity, extinction, and kinematic properties. We also tested the effectiveness of unsupervised machine-learning algorithms in distinguishing between the HII and DIG regions. The gas density in the HII and DIG regions is close to the low-density limit in all fields. The average velocity dispersion in the DIG is higher than in the HII regions, which can be explained by the DIG being 1.8 kK hotter than HII gas. The DIG manifests a lower ionization parameter than HII gas, and the DIG fractions vary between 15-77%, with strong evidence of a contribution by hot low-mass evolved stars and shocks to the DIG ionization. Most of the DIG is consistent with no extinction and an oxygen metallicity that is indistinguishable from that of the HII gas.We observe a flat metallicity profile in the central region, without a sign of a gradient. The differences between extremely faint HII and DIG regions follow the same trends and correlations as their much brighter cousins. HII and DIG are so heterogeneous, however, that the differences within each class are larger than the differences between the two classes., Comment: Accepted in A&A

BlueMUSE is an integral field spectrograph in an early development stage for the ESO VLT. For our design of the data reduction software for this instrument, we are first reviewing capabilities and issues of the pipeline of the existing MUSE instrument. MUSE has been in operation at the VLT since 2014 and led to discoveries published in more than 600 refereed scientific papers. While BlueMUSE and MUSE have many common properties we briefly point out a few key differences between both instruments. We outline a first version of the flowchart for the science reduction, and discuss the necessary changes due to the blue wavelength range covered by BlueMUSE. We also detail specific new features, for example, how the pipeline and subsequent analysis will benefit from improved handling of the data covariance, and a more integrated approach to the line-spread function, as well as improvements regarding the wavelength calibration which is of extra importance in the blue optical range. We finally discuss how simulations of BlueMUSE datacubes are being implemented and how they will be used to prepare the science of the instrument., Comment: 10 pages, 4 figures. Presented at SPIE "Astronomical Telescopes and Instrumentation", July 2022, in Montr\'eal, as part of "Software and Cyberinfrastructure for Astronomy VII"

We investigate the warm ionized gas in the blue compact galaxy (BCG) Haro 14 by means of integral field spectroscopic observations taken with the MUSE/VLT. The large FoV of MUSE and its unprecedented sensitivity enable observations of the galaxy nebular emission up to large galactocentric distances. This allowed us to trace the ionized gas morphology and ionization structure up to kiloparsec scales and, for the first time, to accurately investigate the excitation mechanism operating in the outskirts of a typical BCG. The intensity and diagnostic maps reveal at least two highly distinct components of ionized gas: the bright central regions, mostly made of individual clumps, and a faint component which extends up to kiloparsec scales and consists of widespread diffuse emission, well-delineated filamentary structures, and faint knots. Noteworthy are the two curvilinear filaments extending up to 2 and 2.3 kpc southwest, which likely trace the edges of supergiant expanding bubbles driven by galactic outflows. We find that while the central clumps in Haro 14 are HII-region complexes, the morphology and line ratios of the whole low-surface-brightness component are not compatible with star formation photoionization. In the spatially resolved emission-line-ratio diagnostic diagrams, spaxels above the maximum starburst line form the majority. Moreover, our findings suggest that more than one alternative mechanism is ionizing the outer galaxy regions. The properties of the diffuse component are consistent with ionization by diluted radiation and the large filaments and shells are most probably shocked areas at the edge of bubbles. The mechanism responsible for the ionization of the faint individual clumps observed in the galaxy periphery is more difficult to assess. These clumps could be the shocked debris of fragmented shells or regions where star formation is proceeding under extreme conditions., Comment: 27 pages, 13 figures. Accepted in A&A

The UV bump is a broad absorption feature centered at 2175{\AA} that is seen in the attenuation/extinction curve of some galaxies, but its origin is not well known. Here, we use a sample of 86 star-forming galaxies at z=1.7-2.7 with deep rest-frame UV spectroscopy from the MUSE HUDF Survey to study the connection between the strength of the observed UV 2175{\AA} bump and the Spitzer/MIPS 24 micron photometry, which at the redshift range of our sample probes mid-IR polycyclic aromatic hydrocarbon (PAH) emission at ~6-8 micron. The sample has robust spectroscopic redshifts and consists of typical main-sequence galaxies with a wide range in stellar mass (log(Mstar/Msun) ~ 8.5-10.7) and star formation rates (SFRs; SFR ~ 1-100 Msun/yr). Galaxies with MIPS detections have strong UV bumps, except for those with mass-weighted ages younger than ~150 Myr. We find that the UV bump amplitude does not change with SFR at fixed stellar mass but increases with mass at fixed SFR. The UV bump amplitude and the PAH strength (defined as mid-IR emission normalized by SFR) are highly correlated and both also correlate strongly with stellar mass. We interpret these correlations as the result of the mass-metallicity relationship, such that at low metallicities PAH emission is weak due to a lower abundance of PAH molecules. The weak or complete absence of the 2175{\AA} bump feature on top of the underlying smooth attenuation curve at low mass/metallicities is then expected if the PAH carriers are the main source of the additional UV absorption., Comment: 9 pages, 7 figures, accepted for publication in mnras

Supermassive black hole masses (MBH) can dynamically be estimated with various methods and using different kinematic tracers. Different methods have only been cross-checked for a small number of galaxies and often show discrepancies. To understand these discrepancies, detailed cross-comparisons of additional galaxies are needed. We present the first part of our cross-comparison between stellar- and gas-based MBH estimates in the nearby fast-rotating early-type galaxy NGC 6958. The measurements presented here are based on ground-layer adaptive optics-assisted Multi-Unit Spectroscopic Explorer (MUSE) science verification data at around 0.6 arcsec spatial resolution. The spatial resolution is a key ingredient for the measurement and we provide a Gaussian parametrisation of the adaptive optics-assisted point spread function (PSF) for various wavelengths. From the MUSE data, we extracted the stellar kinematics and constructed dynamical models. Using an axisymmetric Schwarzschild technique, we measured an MBH of (3.6+2.7-2.4)\times 10^8 Msun at 3\sigma significance taking kinematical and dynamical systematics (e.g.,radially-varying mass-to-light ratio) into account. We also added a dark halo, but our data does not allow to constrain the dark matter fraction. Adding dark matter with an abundance matching prior results in a 25 per cent more massive black hole. Jeans anisotropic models return MBH of (4.6+2.5-2.7) \times 10^8 Msun and (8.6+0.8-0.8) \times 10^8 Msun at 3\sigma confidence for spherical and cylindrical alignment of the velocity ellipsoid, respectively. In a follow-up study, we will compare the stellar-based MBH with those from cold and warm gas tracers, which will provide additional constraints for the MBH for NGC 6958, and insights into assumptions that lead to potential systematic uncertainty., Comment: Accepted for publication in MNRAS, 21 pages, 18 Figures (Including Appendix). Supplementary information is added at the end of the file (3 pages)

We use spectra observed with the integral-field spectrograph MUSE to reveal the central kinematics of the Galactic globular cluster Messier 80 (M80, NGC 6093). Using observations obtained with the recently commissioned narrow-field mode of MUSE, we are able to analyse 932 stars in the central 7.5 arcsec by 7.5 arcsec of the cluster for which no useful spectra previously existed. Mean radial velocities of individual stars derived from the spectra are compared to predictions from axisymmetric Jeans models, resulting in radial profiles of the velocity dispersion, the rotation amplitude, and the mass-to-light ratio. The new data allow us to search for an intermediate-mass black hole (IMBH) in the centre of the cluster. Our Jeans model finds two similarly probable solutions around different dynamical cluster centres. The first solution has a centre close to the photometric estimates available in the literature and does not need an IMBH to fit the observed kinematics. The second solution contains a location of the cluster centre that is offset by about 2.4 arcsec from the first one and it needs an IMBH mass of $4600^{+1700}_{-1400}$ solar masses. N-body models support the existence of an IMBH in this cluster with a mass of up to 6000 solar masses in this cluster, although models without an IMBH provide a better fit to the observed surface brightness profile. They further indicate that the cluster has lost nearly all stellar-mass black holes. We further discuss the detection of two potential high-velocity stars with radial velocities of 80 to 90 km/s relative to the cluster mean., Comment: accepted for publication in MNRAS; data can be found here: https://doi.org/10.25625/VCNHOR

(Abriged) Blue compact galaxies (BCGs) are low-luminosity, metal-poor, gas-rich objects that form stars at high rates, excellent analogs to the high-redshift star-forming galaxy population. Being low-mass starbursts, they also constitute ideal laboratories for investigating star formation and massive stellar feedback. This work presents results from integral field spectroscopic observations of the BCG Haro 14 taken with the Multi Unit Spectroscopic Explorer (MUSE). The large MUSE field of view enables simultaneous observations of the starburst and the host galaxy. We built galaxy maps in continuum and in emission lines and generated synthetic VRI images, from which we produced color index maps and surface brightness profiles. We detected numerous clumps spread throughout the galaxy, both in continuum and in emission lines, and produced a catalog with their position, size, and photometry. This analysis allowed us to study the morphology and stellar populations of Haro 14 in detail. The stellar distribution shows a pronounced asymmetry; the intensity peak in continuum is not centered with respect to the stellar host but is displaced by about 500 pc southwest. At the position of the continuum peak we find a bright stellar cluster that with M$_{V}=-12.18$ appears as a strong super stellar cluster candidate. We also find a highly asymmetric, blue, but nonionizing stellar component that occupies almost the whole eastern part of the galaxy. We conclude that there are at least three different stellar populations in Haro 14: the current starburst of about 6 Myr; an intermediate-age component of between ten and several hundred million years; and a red and regular host of several gigayears. The pronounced lopsidedness in the continuum and also in the color maps, and the presence of numerous stellar clusters, are consistent with a scenario of mergers or interactions acting in Haro 14., Comment: 20 pages, 6 figures. Accepted in A&A

The [O III ] 5007 Planetary Nebula Luminosity Function (PNLF) is an established distance indicator that has been used for more than 30 years to measure the distances of galaxies out to ~15 Mpc. With the advent of the Multi-Unit Spectroscopic Explorer on the Very Large Telescope (MUSE) as an efficient wide-field integral field spectrograph, the PNLF method is due for a renaissance, as the spatial and spectral information contained in the instrument's datacubes provides many advantages over classical narrow-band imaging. Here we use archival MUSE data to explore the potential of a novel differential emission-line filter (DELF) technique to produce spectrophotometry that is more accurate and more sensitive than other methods. We show that DELF analyses are superior to classical techniques in high surface brightness regions of galaxies and we validate the method both through simulations and via the analysis of data from two early-type galaxies (NGC 1380 and NGC 474) and one late-type spiral (NGC 628). We demonstrate that with adaptive optics support or under excellent seeing conditions, the technique is capable of producing precision (< 0.05 mag) [O III ] photometry out to distances of 40 Mpc while providing discrimination between planetary nebulae and other emission-line objects such as H II regions, supernova remnants, and background galaxies. These capabilities enable us to use MUSE to measure precise PNLF distances beyond the reach of Cepheids and the tip of the red giant branch method, and become an additional tool for constraining the local value of the Hubble constant., Comment: 59 pages, 36 figures. Accepted for publication at ApJ

We report the discovery of diffuse extended Ly-alpha emission from redshift 3.1 to 4.5, tracing cosmic web filaments on scales of 2.5-4 comoving Mpc. These structures have been observed in overdensities of Ly-alpha emitters in the MUSE Extremely Deep Field, a 140 hour deep MUSE observation located in the Hubble Ultra Deep Field. Among the 22 overdense regions identified, 5 are likely to harbor very extended Ly-alpha emission at high significance with an average surface brightness of $\mathrm{5 \times 10^{-20} erg s^{-1} cm^{-2} arcsec^{-2}}$. Remarkably, 70% of the total Ly-alpha luminosity from these filaments comes from beyond the circumgalactic medium of any identified Ly-alpha emitters. Fluorescent Ly-alpha emission powered by the cosmic UV background can only account for less than 34% of this emission at z$\approx$3 and for not more than 10% at higher redshift. We find that the bulk of this diffuse emission can be reproduced by the unresolved Ly-alpha emission of a large population of ultra low luminosity Ly-alpha emitters ($\mathrm{<10^{40} erg s^{-1}}$), provided that the faint end of the Ly-alpha luminosity function is steep ($\alpha \lessapprox -1.8$), it extends down to luminosities lower than $\mathrm{10^{38} - 10^{37} erg s^{-1}}$ and the clustering of these Ly-alpha emitters is significant (filling factor $< 1/6$). If these Ly-alpha emitters are powered by star formation, then this implies their luminosity function needs to extend down to star formation rates $\mathrm{< 10^{-4} M_\odot yr^{-1}}$. These observations provide the first detection of the cosmic web in Ly-alpha emission in typical filamentary environments and the first observational clue for the existence of a large population of ultra low luminosity Ly-alpha emitters at high redshift., Comment: 28 pages, 19 figures

The evolution of the most massive stars is a puzzle with many missing pieces. Statistical analyses are the key to provide anchors to calibrate theory, however performing these studies is an arduous job. The state-of-the-art integral field spectrograph MUSE has stirred up stellar astrophysicists who are excited about the capability to take spectra of up to a thousand stars in a single exposure. The excitement was even higher with the commissioning of the MUSE narrow-field-mode (NFM) that has demonstrated angular resolutions akin to the Hubble Space Telescope. We present the first mapping of the dense stellar core R136 in the Tarantula nebula based on a MUSE-NFM mosaic. We aim to deliver the first homogeneous analysis of the most massive stars in the local Universe and to explore the impact of these peculiar objects to the interstellar medium., Comment: 10 pages, 5 figures, to appear in edited form in ESO Messenger 182 | 2021

To better characterize the upper end of the galaxy stellar mass range, the MUSE Most Massive Galaxies (M3G) Survey targeted the most massive galaxies (M$>10^{12}$ M$_{\odot}$) found in the densest known clusters of galaxies at $z\sim0.046$. The sample is composed by 25 early-type galaxies: 14 BCGs, of which 3 are in the densest region of the Shapley Super Cluster (SSC), and 11 massive satellites in the SSC. In this work we aim at deriving the spatial distribution and kinematics of the gas, and discussing its ionisation mechanism and origin in the optical wavelength range with MUSE data. We fit the continuum of the spectra using an extensive library of single stellar population models and model the emission lines employing up to three Gaussian functions. In the M3G sample, ionized-gas was detected in 5 BCGs, of which one is in the densest region of the SSC, and 6 massive satellites in the SSC. Among these objects, [OI] and [NI] were detected in 3 BCGs and one satellite. The gas is centrally concentrated in almost all objects, except for 2 BCGs that show filaments and 2 massive satellites with extended emission. The emission line profiles of 3 BCGs present red/blueshifted components. The presence of dust was revealed by analysing Balmer line ratios obtaining a mean $E(B-V)$ of 0.2-0.3. The emission-line diagnostic diagrams show predominately LINER line ratios with little contamination from star formation. The gas was detected in 80% of fast rotators and 35% of slow rotators. The orientations of stellar and gaseous rotations are aligned with respect to each other for 60% of satellites and 25% of BCGs. The presence of misalignments points to an external origin of the gas for 3 BCGs and 2 satellites. On the other hand, some of these systems are characterized by triaxial and prolate-like stellar rotation that could support an internal origin of the gas even in case of misalignments., Comment: 32 pages, 25 figures, A&A accepted

Spectroscopic surveys of massive galaxy clusters reveal the properties of faint background galaxies, thanks to the magnification provided by strong gravitational lensing. We present a systematic analysis of integral-field-spectroscopy observations of 12 massive clusters, conducted with the Multi Unit Spectroscopic Explorer (MUSE). All data were taken under very good seeing conditions (0.6") in effective exposure times between two and 15 hrs per pointing, for a total of 125 hrs. Our observations cover a total solid angle of ~23 arcmin$^2$ in the direction of clusters, many of which were previously studied by the MACS, Frontier Fields, GLASS and CLASH programs. The achieved emission line detection limit at 5$\sigma$ for a point source varies between (0.77--1.5)$\times$10$^{-18}$ erg\,s$^{-1}$\,cm$^{-2}$ at 7000\AA. We present our developed strategy to reduce these observational data, detect sources and determine their redshifts. We construct robust mass models for each cluster to further confirm our redshift measurements using strong-lensing constraints, and identify a total of 312 strongly lensed sources producing 939 multiple images. The final redshift catalogs contain more than 3300 robust redshifts, of which 40\% are for cluster members and $\sim$30\% for lensed Lyman-$\alpha$ emitters. 14\% of all sources are line emitters not seen in the available HST images, even at the depth of the FFs ($\sim29$ AB). We find that the magnification distribution of the lensed sources in the high-magnification regime ($\mu{=}$ 2--25) follows the theoretical expectation of $N(z)\propto\mu^{-2}$. The quality of this dataset, number of lensed sources, and number of strong-lensing constraints enables detailed studies of the physical properties of both the lensing cluster and the background galaxies. The full data products from this work are made available to the community. [abridged], Comment: 47 pages, 13 figures, 18 tables. Accepted for publication in A&A. Data Release products available on CDS and at https://cral-perso.univ-lyon1.fr/labo/perso/johan.richard/MUSE_data_release

We investigate the ultraviolet (UV) spectral properties of faint Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.93. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$\alpha$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$\alpha$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$\alpha$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope., Comment: 26 pages, 18 figures. Accepted for publication in A&A

Processing of raw data from modern astronomical instruments is nowadays often carried out using dedicated software, so-called "pipelines" which are largely run in automated operation. In this paper we describe the data reduction pipeline of the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph operated at ESO's Paranal observatory. This spectrograph is a complex machine: it records data of 1152 separate spatial elements on detectors in its 24 integral field units. Efficiently handling such data requires sophisticated software, a high degree of automation and parallelization. We describe the algorithms of all processing steps that operate on calibrations and science data in detail, and explain how the raw science data gets transformed into calibrated datacubes. We finally check the quality of selected procedures and output data products, and demonstrate that the pipeline provides datacubes ready for scientific analysis., Comment: 31 pages, 19 figures, includes the short appendix. Paper accepted by A&A. Software available from ESO (https://www.eso.org/sci/software/pipelines/muse/muse-pipe-recipes.html) and AIP (https://data.aip.de/projects/musepipeline.html)

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., Comment: 28 pages, 29 figures, 7 tables, accepted for publication in A&A

We have modelled the stellar and nebular continua and emission-line intensity ratios of massive stellar populations in the Antennae galaxy using high resolution and self-consistent libraries of model HII regions around central clusters of aging stars. The model libraries are constructed using the stellar population synthesis code, Starburst99, and photoionisation model, Cloudy. The Geneva and PARSEC stellar evolutionary models are plugged into Starburst99 to allow comparison between the two models. Using a spectrum-fitting methodology that allows the spectral features in the stellar and nebular continua (e.g. Wolf-Rayet features, Paschen jump), and emission-line diagnostics to constrain the models, we apply the libraries to the high-resolution MUSE spectra of the starbursting regions in the Antennae galaxy. Through this approach, we were able to model the continuum emission from Wolf-Rayet stars and extract stellar and gas metallicities, ages, electron temperatures and densities of starbursts by exploiting the full spectrum. From the application to the Antennae galaxy, we find that (1) the starbursts in the Antennae galaxy are characterised by stellar and gas metallicities of around solar, (2) the star-forming gas in starbursts in the Western loop of NGC 4038 appear to be more enriched, albeit slightly, than the rest of galaxy, (3) the youngest starbursts are found across the overlap region and over parts of the western-loop, though in comparison, the regions in the western-loop appear to be at a slightly later stage in star-formation than the overlap region, and (4) the results obtained from fitting the Geneva and Parsec models are largely consistent., Comment: Submitted to MNRAS - 27 figures, 34 pages including 2 appendices

High-performance instruments at large ground-based telescopes have made integral field spectroscopy (IFS) a powerful tool for the study of extended objects such as galaxies, nebulae, or even larger survey fields on the sky. Here we discuss the capabilities of IFS for the study of resolved stellar populations, using the new method of PSF-fitting crowded field IFS, analogous to the well-established technique of crowded field photometry with image sensors. We review early pioneering work with first generation integral field spectrographs, the breakthrough achieved with the MUSE instrument at the ESO Very Large Telescope, the remarkable progress accomplished with MUSE in the study of globular clusters, and first results on nearby galaxies. We discuss the synergy of integral field spectrographs at 8-10m class telescopes with future facilities such as the Extremely Large Telescope (ELT)., Comment: 6 pages, 3 figures, accepted for AN special issue for STARS2019/SMFNS2019

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., Comment: Resubmitted to A&A after favourable referee report. 4 pages, 2 figures

We present our ongoing work of using two independent tracers to estimate the supermassive black hole mass in the nearby early-type galaxy NGC 6958; namely integrated stellar and molecular gas kinematics. We used data from the Atacama Large Millimeter/submillimeter Array (ALMA), and the adaptive-optics assisted Multi-Unit Spectroscopic Explorer (MUSE) and constructed state-of-the-art dynamical models. The different methods provide black hole masses of $(2.89\pm 2.05) \times 10^8M_{\odot}$ from stellar kinematics and $(1.35\pm 0.09) \times 10^8M_{\odot}$ from molecular gas kinematics which are consistent within their $3\sigma$ uncertainties. Compared to recent M$_{\rm BH}$ - $\sigma_{\rm e}$ scaling relations, we derive a slightly over-massive black hole. Our results also confirm previous findings that gas-based methods tend to provide lower black hole masses than stellar-based methods. More black hole mass measurements and an extensive analysis of the method-dependent systematics are needed in the future to understand this noticeable discrepancy., Comment: 5 pages, 3 figures, Proceeding of IAU Symposium 353, Galactic Dynamics in the Era of Large Surveys, ed. M. Valluri & J. A. Sellwood, Cambridge Univ. Press, in press; submitted September, 2019

We have obtained data for 41 star forming galaxies in the MUSE Atlas of Disks (MAD) survey with VLT/MUSE. These data allow us, at high resolution of a few 100 pc, to extract ionized gas kinematics ($V, \sigma$) of the centers of nearby star forming galaxies spanning 3 dex in stellar mass. This paper outlines the methodology for measuring the ionized gas kinematics, which we will use in subsequent papers of this survey. We also show how the maps can be used to study the kinematics of diffuse ionized gas for galaxies of various inclinations and masses. Using two different methods to identify the diffuse ionized gas, we measure rotation velocities of this gas for a subsample of 6 galaxies. We find that the diffuse ionized gas rotates on average slower than the star forming gas with lags of 0-10 km/s while also having higher velocity dispersion. The magnitude of these lags is on average 5 km/s lower than observed velocity lags between ionized and molecular gas. Using Jeans models to interpret the lags in rotation velocity and the increase in velocity dispersion we show that most of the diffuse ionized gas kinematics are consistent with its emission originating from a somewhat thicker layer than the star forming gas, with a scale height that is lower than that of the stellar disk., Comment: 22 pages, accepted by MNRAS

NGC6240 is a well-studied nearby galaxy system in the process of merging. Based on optical, X-ray, and radio observations, it is thought to harbor two active nuclei We carried out a detailed optical 3D spectroscopic study to investigate the inner region of this system in connection with existing MERLIN and VLBA data. We observed NGC6240 with very high spatial resolution using the MUSE instrument in the Narrow-Field Mode with the four-laser GALACSI adaptive optics system on the ESO VLT under seeing conditions of 0.49 arcsec. Our 3D spectra cover the wavelength range from 4725 to 9350 AA at a spatial resolution of ~75 mas. We report the discovery of three nuclei in the final state of merging within a region of only 1 kpc in the NGC6240 system.Thanks to MUSE we are able to show that the formerly unresolved southern component actually consists of two distinct nuclei separated by only 198 pc. In combination with Gaia data we reach an absolute positional accuracy of only 30 mas that is essential to compare optical spectra with MERLIN and VLBA radio positions. The verification and detailed study of a system with three nuclei, two of which are active and each with a mass in excess of $9\times10^{7} M_{\odot}$, is of great importance for the understanding of hierarchical galaxy formation via merging processes since multiple mergers lead to a faster evolution of massive galaxies in comparison to binary mergers. So far it has been suggested that the formation of galactic nuclei with multiple supermassive black holes (SMBHs) is expected to be rare in the local universe.Triple massive black hole systems might be of fundamental importance for the coalescence of massive black hole binaries in less than a Hubble time leading to the loudest sources of gravitational waves in the millihertz regime., Comment: 15 pages, 13 figures, typos corrected, Astronomy & Astrophysics in press

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., Comment: accepted for publication in Astronomy & Astrophysics

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., Comment: Accepted for publication in A&A

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 existence of a spatially resolved Star-Forming Main Sequence (rSFMS) and a spatially resolved Mass-Metallicity Relation (rMZR) is now well established for local galaxies. Moreover, gradients with metallicity decreasing with radius seem to be common in local disc galaxies. These observations suggest that galaxy formation is a self-regulating process, and provide constraints for galaxy evolution models. Studying the evolution of these relations at higher redshifts is still however very challenging. In this paper, we analyse three gravitationally lensed galaxies at z = 0.6, 0.7 and 1, observed with MUSE and SINFONI. These galaxies are highly magnified by galaxy clusters, which allow us to observe resolved scaling relations and metallicity gradients on physical scales of a couple of hundred parsecs, comparable to studies of local galaxies. We confirm that the rSFMS is already in place at these redshifts on sub-kpc scales, and establish, for the first time, the existence of the rMZR at higher redshifts. We develop a forward-modelling approach to fit 2D metallicity gradients of multiply imaged lensed galaxies in the image plane, and derive gradients of -0.027+/-0.003, -0.019+/-0.003 and -0.039+/-0.060 dex/kpc. Despite the fact that these are clumpy galaxies, typical of high redshift discs, the metallicity variations in the galaxies are well described by global linear gradients, and we do not see any difference in metallicity associated with the star-forming clumps., Comment: 12 pages, 10 figures, Accepted for publication in MNRAS

Observations of Neptune, made in 2018 using the new Narrow Field Adaptive Optics mode of the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope (VLT) from 0.48 - 0.93 micron, are analysed here to determine the latitudinal and vertical distribution of cloud opacity and methane abundance in Neptune's observable troposphere (0.1 - ~3 bar). Previous observations at these wavelengths in 2003 by HST/STIS (Karkoschka and Tomasko 2011, Icarus 205, 674-694) found that the mole fraction of methane above the cloud tops (at ~2 bar) varied from ~4% at equatorial latitudes to ~2% at southern polar latitudes, by comparing the observed reflectivity at wavelengths near 825 nm controlled primarily by either methane absorption or H2-H2/H2-He collision-induced absorption. We find a similar variation in cloud-top methane abundance in 2018, which suggests that this depletion of methane towards Neptune's pole is potentially a long-lived feature, indicative of long-term upwelling at mid-equatorial latitudes and subsidence near the poles. By analysing these MUSE observations along the central meridian with a retrieval model, we demonstrate that a broad boundary between the nominal and depleted methane abundances occurs at between 20 - 40S. We also find a small depletion of methane near the equator, perhaps indicating subsidence there, and a local enhancement near 60 - 70S, which we suggest may be associated with South Polar Features (SPFs) seen in Neptune's atmosphere at these latitudes. Finally, by the use of both a reflectivity analysis and a principal component analysis, we demonstrate that this depletion of methane towards the pole is apparent at all locations on Neptune's disc, and not just along its central meridian.

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., Comment: 7 pages, 3 figures; accepted for publication in Astronomy & Astrophysics

The formation of globular clusters remains an open debate. Dwarf starburst galaxies are efficient at forming young massive clusters with similar masses as globular clusters and may hold the key to understanding their formation. We study star cluster formation in a tidal debris - including the vicinity of three tidal dwarf galaxies - in a massive gas dominated collisional ring around NGC~5291. These dwarfs have physical parameters which differ significantly from local starbursting dwarfs. They are gas-rich, highly turbulent, have a gas metallicity already enriched up to half-solar, and are expected to be free of dark matter. The aim is to study massive star cluster formation in this as yet unexplored type of environment. We use imaging from the Hubble Space Telescope using broadband filters covering the wavelength range from the near-ultraviolet to the near-infrared. We determine the masses and ages of the cluster candidates by using the spectral energy distribution-fitting code CIGALE, carefully considering age-extinction degeneracy effects on the estimation of the physical parameters. We find that the tidal dwarf galaxies in the ring of NGC 5291 are forming star clusters with an average efficiency of $\sim40\%$, comparable to blue compact dwarf galaxies. We also find massive star clusters for which the photometry suggests that they were formed at the very birth of the tidal dwarf galaxies and have survived for several hundred million years. Therefore our study shows that extended tidal dwarf galaxies and compact clusters may be formed simultaneously. In the specific case observed here, the young star clusters are not massive enough to survive for a Hubble time. However one may speculate that similar objects at higher redshift, with higher star formation rate, might form some of the long lived globular clusters., Comment: Accepted for publication in A&A. 19 pages, 17 figures

We present the MUSE-Wide survey, a blind, 3D spectroscopic survey in the CANDELS/GOODS-S and CANDELS/COSMOS regions. Each MUSE-Wide pointing has a depth of 1 hour and hence targets more extreme and more luminous objects over 10 times the area of the MUSE-Deep fields (Bacon et al. 2017). The legacy value of MUSE-Wide lies in providing "spectroscopy of everything" without photometric pre-selection. We describe the data reduction, post-processing and PSF characterization of the first 44 CANDELS/GOODS-S MUSE-Wide pointings released with this publication. Using a 3D matched filtering approach we detected 1,602 emission line sources, including 479 Lyman-$\alpha$ (Lya) emitting galaxies with redshifts $2.9 \lesssim z \lesssim 6.3$. We cross-match the emission line sources to existing photometric catalogs, finding almost complete agreement in redshifts and stellar masses for our low redshift (z < 1.5) emitters. At high redshift, we only find ~55% matches to photometric catalogs. We encounter a higher outlier rate and a systematic offset of $\Delta$z$\simeq$0.2 when comparing our MUSE redshifts with photometric redshifts. Cross-matching the emission line sources with X-ray catalogs from the Chandra Deep Field South, we find 127 matches, including 10 objects with no prior spectroscopic identification. Stacking X-ray images centered on our Lya emitters yielded no signal; the Lya population is not dominated by even low luminosity AGN. A total of 9,205 photometrically selected objects from the CANDELS survey lie in the MUSE-Wide footprint, which we provide optimally extracted 1D spectra of. We are able to determine the spectroscopic redshift of 98% of 772 photometrically selected galaxies brighter than 24th F775W magnitude. All the data in the first data release - datacubes, catalogs, extracted spectra, maps - are available on the website https://musewide.aip.de. [abridged], Comment: 25 pages 15+1 figures. Accepted, A&A. Comments welcome

Galaxies are surrounded by large reservoirs of gas, mostly hydrogen, fed by inflows from the intergalactic medium and by outflows due to galactic winds. Absorption-line measurements along the sightlines to bright and rare background quasars indicate that this circumgalactic medium pervades far beyond the extent of starlight in galaxies, but very little is known about the spatial distribution of this gas. A new window into circumgalactic environments was recently opened with the discovery of ubiquitous extended Lyman-alpha emission from hydrogen around high-redshift galaxies, facilitated by the extraordinary sensitivity of the MUSE instrument at the ESO Very Large Telescope. Due to the faintness of this emission, such measurements were previously limited to especially favourable systems or to massive statistical averaging. Here we demonstrate that low surface brightness Lyman-alpha emission surrounding faint galaxies at redshifts between 3 and 6 adds up to a projected sky coverage of nearly 100%. The corresponding rate of incidence (the mean number of Lyman-alpha emitters penetrated by any arbitrary line of sight) is well above unity and similar to the incidence rate of high column density absorbers frequently detected in the spectra of distant quasars. This similarity suggests that most circumgalactic atomic hydrogen at these redshifts has now been detected also in emission., Comment: published in Nature, online on 01 Oct 2018, in print on 11 Oct 2018. This is the Authors' version, for the published version see https://www.nature.com/articles/s41586-018-0564-6 ; free viewing access via https://rdcu.be/8eCx