MOONS is the new Multi-Object Optical and Near-infrared Spectrograph currently under construction for the Very Large Telescope (VLT) at ESO. This remarkable instrument combines, for the first time, the collecting power of an 8-m telescope, 1000 fibres with individual robotic positioners, and both low- and high-resolution simultaneous spectral coverage across the 0.64–1.8 μm wavelength range. This facility will provide the astronomical community with a powerful, world-leading instrument able to serve a wide range of Galactic, extragalactic and cosmological studies. Construction is now proceeding full steam ahead and this overview article presents some of the science goals and the technical description of the MOONS instrument. More detailed information on the MOONS surveys is provided in the other dedicated articles in this Messenger issue.

ABSTRACT We have obtained data for 41 star forming galaxies in the MUSE Atlas of Discs (MAD) survey with VLT/MUSE. These data allow us, at high resolution of a few 100 pc, to extract ionized gas kinematics (V, σ) of the centres 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 six galaxies. We find that the diffuse ionized gas rotates on average slower than the star forming gas with lags of 0–10 km s−1 while also having higher velocity dispersion. The magnitude of these lags is on average 5 km s−1 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 disc.

The Multi Object Optical and Near-infrared Spectrograph (MOONS) instrument is the next generation multi-object spectrograph for the VLT. This powerful instrument will combine for the first time: the large collecting power of the VLT with a high multipexing capability offered by 1000 optical fibres moved with individual robotic positioners and a novel, very fast spectrograph able to provide both low- and high-resolution spectroscopy simultaneously across the wavelength range 0.64 mu m - 1.8 mu m. Such a facility will provide the astronomical community with a powerful, world-leading instrument able to serve a wide range of Galactic, Extragalactic and Cosmological studies. The final assembly, integration and verification phase of the instrument is now about to start performance testing.

Following the success of the Plateau de Bure high-z Blue Sequence Survey (PHIBSS), we present the PHIBSS2 legacy program, a survey of the molecular gas properties of star-forming galaxies on and around the star-formation main sequence (MS) at different redshifts using IRAM’s NOrthern Extended Millimeter Array (NOEMA). This survey significantly extends the existing sample of star-forming galaxies with CO molecular gas measurements, probing the peak epoch of star formation (z = 1 − 1.6) as well as its building-up (z = 2 − 3) and winding-down (z = 0.5 − 0.8) phases. The targets are drawn from the well-studied GOODS, COSMOS, and AEGIS cosmological deep fields and uniformly sample the MS in the stellar mass (M⋆) – star formation rate (SFR) plane with log(M⋆/M⊙) = 10 − 11.8 and SFR = 3.5 − 500 M⊙ yr−1 without morphological selection, thus providing a statistically meaningful census of star-forming galaxies at different epochs. We describe the survey strategy and sample selection before focusing on the results obtained at redshift z = 0.5 − 0.8, where we report 60 CO(2-1) detections out of 61 targets. We determine molecular gas masses between 2 × 109 and 5 × 1010 M⊙ and separately obtain disc sizes and bulge-to-total (B/T) luminosity ratios from HST I-band images. The median molecular gas-to-stellar mass ratio μgas∼ = 0.28 ± 0.04, gas fraction fgas∼ = 0.22 ± 0.02, and depletion time $ \widetilde{t_{\mathrm{depl}}} = 0.84 \pm 0.07\,\mathrm{Gyr} $ as well as their dependence with stellar mass and offset from the MS follow published scaling relations for a much larger sample of galaxies spanning a significantly wider range of redshifts, the cosmic evolution of the SFR being mainly driven by that of the molecular gas fraction. The galaxy-averaged molecular Kennicutt–Schmidt (KS) relation between molecular gas and SFR surface densities is strikingly linear, pointing towards similar star formation timescales within galaxies at any given epoch. In terms of morphology, the molecular gas content, the SFR, the disc stellar mass, and the disc molecular gas fraction do not seem to correlate with B/T and the stellar surface density, which suggests an ongoing supply of fresh molecular gas to compensate for the build-up of the bulge. Our measurements do not yield any significant variation of the depletion time with B/T and hence no strong evidence for morphological quenching within the scatter of the MS.

Following the success of the Plateau de Bure high-z Blue Sequence Survey (PHIBSS), we present the PHIBSS2 legacy program, a survey of the molecular gas properties of star-forming galaxies on and around the star formation main sequence (MS) at different redshifts using NOEMA. This survey significantly extends the existing sample of star-forming galaxies with CO molecular gas measurements, probing the peak epoch of star formation (z=1-1.6) as well as its building-up (z=2-3) and winding-down (z=0.5-0.8) phases. The targets are drawn from the GOODS, COSMOS, and AEGIS deep fields and uniformly sample the MS in the stellar mass (M*) - star formation rate (SFR) plane with log(M*/Msun) = 10-11.8. We describe the survey strategy and sample selection before focusing on the results obtained at z=0.5-0.8, where we report 60 CO(2-1) detections out of 61 targets. We determine their molecular gas masses and separately obtain disc sizes and bulge-to-total (B/T) luminosity ratios from HST I-band images. The median molecular gas-to-stellar mass ratio, gas fraction, and depletion time as well as their dependence with M* and offset from the MS follow published scaling relations for a much larger sample of galaxies spanning a wider range of redshifts. The galaxy-averaged Kennicutt-Schmidt (KS) relation between molecular gas and SFR surface densities is strikingly linear, pointing towards similar star formation timescales within galaxies at any given epoch. In terms of morphology, the molecular gas content, the SFR, the disc stellar mass, and the disc molecular gas fraction do not seem to correlate with B/T and the stellar surface density, which suggests an ongoing supply of fresh molecular gas to compensate for the build-up of the bulge. Our measurements do not yield any significant variation of the depletion time with B/T and hence no strong evidence for morphological quenching within the scatter of the MS.

We study the physical properties of the ionized gas in local discs using the sample of 38 nearby similar to 10(8.5-11.2) M-circle dot Star-Forming Main-Sequence (SFMS) galaxies observed so far as part of the MUSE Atlas of Disks (MAD). Specifically, we use all strong emission lines in the MUSE wavelength range 4650-9300 angstrom to investigate the resolved ionized gas properties on similar to 100 pc scales. This spatial resolution enables us to disentangle HII regions from the diffuse ionized gas (DIG) in the computation of gas metallicities and star formation rates (SFRs) of star-forming regions. The gas metallicities generally decrease with radius. The metallicity of the H II regions is on average similar to 0.1 dex higher than that of the DIG, but the metallicity radial gradient in both components is similar. The mean metallicities within the inner galaxy cores correlate with the total stellar mass of the galaxies. On our < 100 pc scales, we find two correlations previously reported at kpc scales: a spatially resolved mass-metallicity relation (RMZR) and a spatially resolved SFMS (RSFMS). We find no secondary dependence of the RMZR with the SFR density. We find that both resolved relations have a local origin, as they do not depend on the total stellar mass. The observational results of this paper are consistent with the inside-out scenario for the growth of galactic disks.

We present a measurement of the fraction of Lyman alpha (Ly alpha) emitters (XLy alpha) amongst HST continuum-selected galaxies at 3 < z < 6 with the Multi-Unit Spectroscopic Explorer (MUSE) on the VLT. Making use of the first 24 MUSE-Wide pointings in GOODS-South, each having an integration time of 1 h, we detect 100 Ly alpha emitters and find X-Ly alpha greater than or similar to 0.5 for most of the redshift range covered, with 29 per cent of the Ly alpha sample exhibiting rest equivalent widths (rest-EWs) <= 15 angstrom. Adopting a range of rest-EW cuts (0-75 angstrom), we find no evidence of a dependence of X-Ly alpha on either redshift or ultraviolet luminosity.

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(-1), 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 +/- 0.41 M circle dot. 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.

We use a newly assembled sample of 3545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z = 1.5–2.5 to constrain the relationship between UV slope (β) and dust attenuation (LIR/LUV ≡ IRX). Our sample significantly extends the range of LUV and β probed in previous UV-selected samples, including those as faint as M1600 = −17.4 ( ) and −2.6 ≲ β ≲ 0.0. IRX is measured using stacks of deep Herschel data, and the results are compared with predictions of the IRX−β relation for different assumptions of the stellar population model and obscuration curve. We find that z = 1.5–2.5 galaxies have an IRX−β relation that is consistent with the predictions for an SMC curve if we invoke subsolar-metallicity models currently favored for high-redshift galaxies, while the commonly assumed starburst curve overpredicts the IRX at a given β by a factor of ≳3. IRX is roughly constant with LUV for LUV ≳ 3 × 109 L⊙. Thus, the commonly observed trend of fainter galaxies having bluer β may simply reflect bluer intrinsic slopes for such galaxies, rather than lower obscurations. The IRX−β relation for young/low-mass galaxies at z ≳ 2 implies a dust curve that is steeper than the SMC. The lower attenuations and higher ionizing photon output for low-metallicity stellar populations point to Lyman continuum production efficiencies, ξion, that may be elevated by a factor of ≈2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z ∼ 2.

We present an analysis of MUSE observations obtained on the massive Frontier Fields cluster Abell 2744. This new dataset covers the entire multiply-imaged region around the cluster core. We measure spectroscopic redshifts for HST-selected continuum sources together with line emitters blindly detected in the datacube. The combined catalog consists of 514 spectroscopic redshifts (with 414 new identifications), including 156 cluster members and 326 magnified background sources. We use this redshift information to perform a strong-lensing analysis of all multiple images previously found in the deep Frontier Field images, and add three new MUSE-detected multiply-imaged systems with no obvious HST counterpart. The combined strong lensing constraints include a total of 60 systems producing 188 images altogether, out of which 29 systems and 83 images are spectroscopically confirmed, making Abell 2744 one of the most well-constrained clusters to date. A parametric mass model including two cluster-scale components in the core and several group-scale substructures at larger radii accurately reproduces all the spectroscopic multiple systems, reaching an rms of 0.67" in the image plane. Overall, the large number of spectroscopic redshifts gives us a robust model and we estimate the systematics on the mass density and magnification within the cluster core to be typically 9

Using infrared data from the Herschel Space Observatory and Karl G. Jansky Very Large Array 3 GHz observations in the COSMOS field, we investigate the redshift evolution of the infrared radio correlation (IRRC) for star-forming galaxies (SFGs) we classify as either spheroid- or disc-dominated based on their morphology. The sample predominantly consists of disc galaxies with stellar mass greater than or similar to 10(10)M(circle dot), and residing on the star-forming main sequence (MS). After the removal of AGN using standard approaches. we observe a significant difference between the redshift evolution of the median IR/radio ratio (q) over bar (TIR) of (i) a sample of ellipticals, plus discs with a substantial bulge component ("spheroid-dominated' SFGs) and, (ii) virtually pure discs and irregular systems ('disc-dominated' SFGs). The spheroid-dominated population follows a declining R versus z trend similar to that measured in recent evolutionary studies of the IRRC. However, for disc-dominated galaxies, where radio and IR emission should be linked to star formation in the most straightforward way, we measure very little change in (q) over bar (TIR). This suggests that low-redshift calibrations of radio emission as a star formation rate (SFR) tracer may remain valid out to at least z similar or equal to 1-1.5 for pure star-forming systems. We find that the different redshift evolution of (RIR for the spheroid- and disc-dominated sample is mainly due to an increasing radio excess for spheroid-dominated galaxies at z greater than or similar to 0.8, hinting at some residual AGN activity in these systems. This finding demonstrates that in the absence of AGN, the IRRC is independent of redshift, and that radio observations can therefore be used to estimate SFRs at all redshifts for genuinely star-forming galaxies.

We present a structural and morphological catalogue for 45 million objects selected from the first year data of the Dark Energy Survey (DES). Single Sersic fits and non-parametric measurements are produced for g, r, and i filters. The parameters from the best-fitting Sersic model (total magnitude, half-light radius, Sersic index, axis ratio, and position angle) are measured with GALFIT; the non-parametric coefficients (concentration, asymmetry, dumpiness, Gini, M20) are provided using the Zurich Estimator of Structural Types (ZEST+). To study the statistical uncertainties, we consider a sample of state-of-the-art image simulations with a realistic distribution in the input parameter space and then process and analyse them as we do with real data: this enables us to quantify the observational biases due to PSF blurring and magnitude effects and correct the measurements as a function of magnitude, galaxy size, Sersic index (concentration for the analysis of the non-parametric measurements) and ellipticity. We present the largest structural catalogue to date: we find that accurate and complete measurements for all the structural parameters are typically obtained for galaxies with SEXTRACTOR MAG_AUTO_I <= 21. Indeed, the parameters in the filters i and r can be overall well recovered up to MAG_AUTO <= 21.5, corresponding to a fitting completeness of similar to 90 per cent below this threshold, for a total of 25 million galaxies. The combination of parametric and non-parametric structural measurements makes this catalogue an important instrument to explore and understand how galaxies form and evolve. The catalogue described in this paper will be publicly released alongside the DES collaboration Y1 cosmology data products at the following URL: https://des.ncsa.illinois.edu/releases.

This paper provides an update of our previous scaling relations (Genzel et al.2015) between galaxy integrated molecular gas masses, stellar masses and star formation rates, in the framework of the star formation main-sequence (MS), with the main goal to test for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ~1mm dust photometry, in a large sample of 1444 star forming galaxies (SFGs) between z=0 and 4. The sample covers the stellar mass range log(M*/M_solar)=9.0-11.8, and star formation rates relative to that on the MS, delta_MS=SFR/SFR(MS), from 10^{-1.3} to 10^{2.2}. Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t_depl, defined as the ratio of molecular gas mass to star formation rate, scales as (1+z)^{-0.6}x(delta_MS)^{-0.44}, and is only weakly dependent on stellar mass. The ratio of molecular-to-stellar mass mu_gas depends on (1+z)^{2.5}x (delta_MS)^{0.52}x(M*)^{-0.36}, which tracks the evolution of the specific star formation rate. The redshift dependence of mu_gas requires a curvature term, as may the mass-dependences of t_depl and mu_gas. We find no or only weak correlations of t_depl and mu_gas with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high-z.

After completion of its final-design review last year, it is full steam ahead for the construction of the MOONS instrument - the next generation multi-object spectrograph for the VLT. This remarkable instrument will combine for the first time: the 8 m collecting power of the VLT, 1000 optical fibres with individual robotic positioners and both medium- and high-resolution spectral coverage acreoss the wavelength range 0.65μm - 1.8 μm. Such a facility will allow a veritable host of Galactic, Extragalactic and Cosmological questions to be addressed. In this paper we will report on the current status of the instrument, details of the early testing of key components and the major milestones towards its delivery to the telescope.

We report the discovery of a 10^4 kpc^2 gaseous structure detected in [OII] in an over-dense region of the COSMOS-Gr30 galaxy group at z~0.725 thanks to deep MUSE Guaranteed Time Observations. We estimate the total amount of diffuse ionised gas to be of the order of (~5+-3)x10^10 Msun and explore its physical properties to understand its origin and the source(s) of the ionisation. The MUSE data allow the identification of a dozen of group members embedded in this structure from emission and absorption lines. We extracted spectra from small apertures defined for both the diffuse ionised gas and the galaxies. We investigated the kinematics and ionisation properties of the various galaxies and extended gas regions thanks to line diagnostics (R23, O32 and [OIII]/H\beta) available within the MUSE wavelength range. We compared these diagnostics to photo-ionisation models and shock models. The structure is divided in two kinematically distinct sub-structures. The most extended sub-structure of ionised gas is likely rotating around a massive galaxy and displays filamentary patterns linking some galaxies. The second sub-structure links another massive galaxy hosting an Active Galactic Nucleus to a low mass galaxy but also extends orthogonally to the AGN host disk over ~35 kpc. This extent is likely ionised by the AGN itself. The location of small diffuse regions in the R23 vs. O32 diagram is compatible with photo-ionisation. However, the location of three of these regions in this diagram (low O32, high R23) can also be explained by shocks, which is supported by their large velocity dispersions. One edge-on galaxy shares the same properties and may be a source of shocks. Whatever the hypothesis, the extended gas seems to be non primordial. We favour a scenario where the gas has been extracted from galaxies by tidal forces and AGN triggered by interactions between at least the two sub-structures.

Recent theoretical models suggest that the early phase of galaxy formation could involve an epoch when galaxies are gas rich but inefficient at forming stars: a “dark galaxy” phase. Here, we report the results of our Multi-Unit Spectroscopic Explorer (MUSE) survey for dark galaxies fluorescently illuminated by quasars at z > 3. Compared to previous studies which are based on deep narrowband (NB) imaging, our integral field survey provides a nearly uniform sensitivity coverage over a large volume in redshift space around the quasars as well as full spectral information at each location. Thanks to these unique features, we are able to build control samples at large redshift distances from the quasars using the same data taken under the same conditions. By comparing the rest-frame equivalent width (EW0) distributions of the Lyα sources detected in proximity to the quasars and in control samples, we detect a clear correlation between the locations of high-EW0 objects and the quasars. This correlation is not seen in other properties, such as Lyα luminosities or volume overdensities, suggesting the possible fluorescent nature of at least some of these objects. Among these, we find six sources without continuum counterparts and EW0 limits larger than 240 Å that are the best candidates for dark galaxies in our survey at z > 3.5. The volume densities and properties, including inferred gas masses and star formation efficiencies, of these dark galaxy candidates are similar to those of previously detected candidates at z ≈ 2.4 in NB surveys. Moreover, if the most distant of these are fluorescently illuminated by the quasar, our results also provide a lower limit of t = 60 Myr on the quasar lifetime.

Cosmological simulations suggest that most of the matter in the Universe is distributed along filaments connecting galaxies. Illuminated by the cosmic UV background (UVB), these structures are expected to glow in fluorescent Ly alpha emission with a surface brightness (SB) that is well below current limits for individual detections. Here, we perform a stacking analysis of the deepest MUSE/VLT data using three-dimensional regions (subcubes) with orientations determined by the position of neighbouring Ly alpha galaxies at 3 < z < 4. Our method increase the probability of detecting filamentary Ly alpha emission, provided that these structures are Lyman-limit systems (LLSs). By stacking 390 oriented subcubes we reach a 2 sigma sensitivity level of SB approximate to 0.44 x 10(-20) erg s(-1) cm(-2) arcsec(-2) in an aperture of 1 arcsec(2) x 6.25 angstrom, three times below the expected fluorescent Ly alpha signal from the Haardt & Madan UVB at z similar to 3.5. No detectable emission is found on intergalactic scales, implying that at least two thirds of our subcubes do not contain oriented LLSs. On the other hand, significant emission is detected in the circumgalactic medium (CGM) in the direction of the neighbours. The signal is stronger for galaxies with a larger number of neighbours and appears to be independent of any other galaxy properties. We estimate that preferentially oriented satellite galaxies cannot contribute significantly to this signal, suggesting instead that gas densities in the CGM are typically larger in the direction of neighbouring galaxies on cosmological scales.