Search

Lyman$-\alpha$ (Ly$\alpha$) emission lines are typically found to be redshifted with respect to the systemic redshifts of galaxies, likely due to resonant scattering of Ly$\alpha$ photons. Here we measure the average velocity offset for a sample of 96 $z\approx3.3$ Ly$\alpha$ emitters (LAEs) with a median Ly$\alpha$ flux (luminosity) of $\approx 10^{-17}~\rm erg~cm^{-2}~s^{-1}$ ($\approx10^{42}~\rm erg~s^{-1}$) and a median star formation rate (SFR) of $\approx1.3 \rm M_{\odot} yr^{-1}$ (not corrected for possible dust extinction), detected by the Multi-Unit Spectroscopic Explorer as part of our MUSEQuBES circumgalactic medium (CGM) survey. By postulating that the stacked CGM absorption profiles of these LAEs, probed by 8 background quasars, must be centered on the systemic redshift, we measure an average velocity offset, V$_{\rm offset} = 171 \pm 8$ $\rm km s^{-1}$, between the Ly$\alpha$ emission peak and the systemic redshift. The observed V$_{\rm offset}$ is lower by factors of $\approx1.4$ and $\approx2.6$ compared to the velocity offsets measured for narrow-band selected LAEs and Lyman break galaxies, respectively, which probe galaxies with higher masses and SFRs. Consistent with earlier studies based on direct measurements for individual objects, we find that the V$_{\rm offset}$ is correlated with the full width at half-maximum of the red peak of the Ly$\alpha$ line, and anti-correlated with the rest-frame equivalent width. Moreover, we find that $V_{\rm offset}$ is correlated with SFR with a sub-linear scaling relation, V$_{\rm offset}\propto \rm SFR^{0.16\pm0.03}$. Adopting the mass scaling for main sequence galaxies, such a relation suggests that V$_{\rm offset}$ scales with the circular velocity of the dark matter halos hosting the LAEs., Comment: Published on 15 May, 2020 along with an erratum (Ref: MNRAS 498, 4424, 2020, doi:10.1093/mnras/staa2668), since Fig. 2 has not been updated in the published version. There are other typos in the published version, introduced by the typesetter, which cannot be revised in the erratum. We feel that the astrp-ph version is more accurate than the published one

We study the physical properties of the ionized gas in local disks using the sample of 38 nearby $\sim10^{8.5-11.2}$M$_\odot$ 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 \AA\ to investigate the resolved ionized gas properties on $\sim$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 HII regions is on average $\sim$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, Comment: 21 pages + appendices. Accepted for publication in MNRAS

We present the first characterization of the circumgalactic medium of Ly α emitters (LAEs), using a sample of 96 z ≈ 3.3 LAEs detected with the VLT/MUSE in fields centred on eight bright background quasars. The LAEs have low Ly α luminosities ($\sim 10^{42}\, \text{erg}\, \text{s}^{-1}$) and star formation rates (SFRs) $\sim 1~\text{M}_\odot \, \text{yr}^{-1}$, which for main-sequence galaxies correspond to stellar masses of only $\sim 10^{8.6}\, \text{M}_\odot$. The median transverse distance between the LAEs and the quasar sightlines is 165 proper kpc (pkpc). We stacked the high-resolution quasar spectra and measured significant excess H i and C iv absorption near the LAEs out to 500 $\text{km}\, \text{s}^{-1}$ and at least ≈250 pkpc (corresponding to ≈7 virial radii). At $\lesssim 30~\text{km}\, \text{s}^{-1}$ from the galaxies, the median H i and C iv optical depths are enhanced by an order of magnitude. The absorption is significantly stronger around the ≈1/3 of our LAEs that are part of ‘groups’, which we attribute to the large-scale structures in which they are embedded. We do not detect any strong dependence of either the H i or C iv absorption on transverse distance (over the range ≈50–250 pkpc), redshift, or the properties of the Ly α emission line (luminosity, full width at half-maximum, or equivalent width). However, for H i, but not C iv, the absorption at $\lesssim 100\, \text{km}\, \text{s}^{-1}$ from the LAE does increase with the SFR. This suggests that LAEs surrounded by more H i tend to have higher SFRs.

We present the results from a MUSE survey of twelve $z\simeq3.15$ quasars, which were selected to be much fainter (20, Comment: Accepted to MNRAS. 17 pages, 7 figures

This is a correction to: Monthly Notices of the Royal Astronomical Society, Volume 496, Issue 2, August 2020, Pages 1013–1022, https://doi.org/10.1093/mnras/staa1347

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time., Comment: 25 pages, 17 figures, 5 tables, accepted for publication in A&A

A prediction of the classic active galactic nuclei (AGN) unification model is the presence of ionisation cones with different orientations depending on the AGN type. Confirmations of this model exist for present times, but it is less clear in the early Universe. Here, we use the morphology of giant Ly$\alpha$ nebulae around AGNs at redshift z$\sim$3 to probe AGN emission and therefore the validity of the AGN unification model at this redshift. We compare the spatial morphology of 19 nebulae previously found around type I AGNs with a new sample of 4 Ly$\alpha$ nebulae detected around type II AGNs. Using two independent techniques, we find that nebulae around type II AGNs are more asymmetric than around type I, at least at radial distances $r>30$~physical kpc (pkpc) from the ionizing source. We conclude that the type I and type II AGNs in our sample show evidence of different surrounding ionising geometries. This suggests that the classical AGN unification model is also valid for high-redshift sources. Finally, we discuss how the lack of asymmetry in the inner parts (r$\lesssim$30 pkpc) and the associated high values of the HeII to Ly$\alpha$ ratios in these regions could indicate additional sources of (hard) ionizing radiation originating within or in proximity of the AGN host galaxies. This work demonstrates that the morphologies of giant Ly$\alpha$ nebulae can be used to understand and study the geometry of high redshift AGNs on circum-nuclear scales and it lays the foundation for future studies using much larger statistical samples., Comment: 15 pages, 13 figures, accepted for publication in MNRAS

Galactic outflows are thought to eject baryons back out to the circum-galactic medium (CGM). Studies based on metal absorption lines (MgII in particular) in the spectra of background quasars indicate that the gas is ejected anisotropically, with galactic winds likely leaving the host in a bi-conical flow perpendicular to the galaxy disk. In this paper, we present a detailed analysis of an outflow from a z = 0.7 "green-valley" galaxy (log($M_*$/$\mathrm{M}_\odot$) = 9.9; SFR = 0.5 $\mathrm{M}_\odot\,\mathrm{yr}^{-1}$) probed by two background sources part of the MUSE Gas Flow and Wind (MEGAFLOW) survey. Thanks to a fortuitous configuration with a background quasar (SDSSJ1358+1145) and a bright background galaxy at $z = 1.4$, both at impact parameters of $\approx 15\,\mathrm{kpc}$, we can - for the first time - probe both the receding and approaching components of a putative galactic outflow around a distant galaxy. We measure a significant velocity shift between the MgII absorption from the two sightlines ($84\pm17\,\mathrm{km}\,\mathrm{s}^{-1}$), which is consistent with the expectation from our simple fiducial wind model, possibly combined with an extended disk contribution., Comment: Submitted to MNRAS. After addressing referee's comments

CO observations allow estimations of the gas content of molecular clouds, which trace the reservoir of cold gas fuelling star formation, as well as to determine extinction via H$_2$ column density, N(H$_2$). Here, we studied millimetric and optical properties at 26 supernovae (SNe) locations of different types in a sample of 23 nearby galaxies by combining molecular $^{12}$C$^{16}$O (J = 1 $\rightarrow$ 0) resolved maps from the EDGE survey and optical Integral Field Spectroscopy from the CALIFA survey. We found an even clearer separation between type II and type Ibc SNe in terms of molecular gas than what we found in the optical using H$\alpha$ emission as a proxy for current SF rate, which reinforces the fact that SNe Ibc are more associated with SF-environments. While A$_V$ at SN locations is similar for SNe II and SNe Ibc, and higher compared to SNe Ia, N(H$_2$) is significantly higher for SNe Ibc than for SNe II and SNe Ia. When compared to alternative extinction estimations directly from SN photometry and spectroscopy, we find that our SNe Ibc have also redder color excess but showed standard Na I D absorption pseudo-equivalent widths ($\sim$1 \AA). In some cases we find no extinction when estimated from the environment, but high amounts of extinction when measured from SN observations, which suggests that circumstellar material or dust sublimation may be playing a role. This work serves as a benchmark for future studies combining last generation millimeter and optical IFS instruments to reveal the local environmental properties of extragalactic SNe., Comment: MNRAS accepted, 17 pages, 8 Figures, 4 Tables

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

We present an analysis of HI Lyman-alpha emission in deep VLT/MUSE observations of two highly magnified and extended galaxies at z=3.5 and 4.03, including a newly discovered, almost complete Einstein ring. While these Lyman-alpha haloes are intrinsically similar to the ones typically seen in other MUSE deep fields, the benefits of gravitational lensing allows us to construct exceptionally detailed maps of Lyman-alpha line properties at sub-kpc scales. By combining all multiple images, we are able to observe complex structures in the Lyman-alpha emission and uncover small (~ 120 km/s in Lyman-alpha peak shift), but significant at > 4 sigma, systematic variations in the shape of the Lyman-alpha line profile within each halo. Indeed, we observe a global trend for the line peak shift to become redder at large radii, together with a strong correlation between the peak wavelength and line width. This systematic intrahalo variation is markedly similar to the object-to-object variations obtained from the integrated properties of recent large samples. Regions of high surface brightness correspond to relatively small line shifts, which could indicate that Lyman-alpha emission escapes preferentially from regions where the line profile has been less severely affected by scattering of Lyman-alpha photons., 9 pages, 4 figures, accepted for publication in MNRAS

We present results from our on-going MusE GAs FLOw and Wind (MEGAFLOW) survey, which consists of 22 quasar lines-of-sight, each observed with the integral field unit (IFU) MUSE and the UVES spectrograph at the ESO Very Large Telescopes (VLT). The goals of this survey are to study the properties of the circum-galactic medium around $z\sim1$ star-forming galaxies. The absorption-line selected survey consists of 79 strong \MgII\ absorbers (with rest-frame equivalent width (REW)$\gtrsim$0.3\AA) and, currently, 86 associated galaxies within 100 projected~kpc of the quasar with stellar masses ($M_\star$) from $10^9$ to $10^{11}$ \msun. We find that the cool halo gas traced by \MgII\ is not isotropically distributed around these galaxies, as we show the strong bi-modal distribution in the azimuthal angle of the apparent location of the quasar with respect to the galaxy major-axis. This supports a scenario in which outflows are bi-conical in nature and co-exist with a coplanar gaseous structure extending at least up to 60 to 80 kpc. Assuming that absorbers near the minor axis probe outflows, the current MEGAFLOW sample allowed us to select 26 galaxy-quasar pairs suitable for studying winds. From this sample, using a simple geometrical model, we find that the outflow velocity only exceeds the escape velocity when $M_{\star}\lesssim 4\times10^9$~\msun, implying the cool material is likely to fall back except in the smallest halos. Finally, we find that the mass loading factor $\eta$, the ratio between the ejected mass rate and the star formation rate (SFR), appears to be roughly constant with respect to the galaxy mass., Comment: 15 pages, 7 figures (+26 in the Appendix), accepted for publication in MNRAS

This paper presents a study of the galactic environment of a chemically-pristine (, Comment: 13 pages, 4 figures, published in MNRAS

We investigate the Lyman α emitter (LAE) luminosity function (LF) within the redshift range 2.9 ≤ z ≤ 6 from the first instalment of the blind integral field spectroscopic MUSE-Wide survey. This initial part of the survey probes a region of 22.2 arcmin2 in the CANDELS/GOODS-S field (24 MUSE pointings with 1h integrations). The dataset provided us with 237 LAEs from which we construct the LAE LF in the luminosity range 42.2 ≤ log LLyα[erg s−1] ≤ 43.5 within a volume of 2.3 × 105 Mpc3. For the LF construction we utilise three different non-parametric estimators: the classical 1/Vmax method, the C− method, and an improved binned estimator for the differential LF. All three methods deliver consistent results, with the cumulative LAE LF being Φ(log L Lyα[erg s−1] = 43.5) ≃ 3 × 10−6 Mpc−3 and Φ(log L Lyα[erg s−1] = 42.2) ≃ 2 × 10−3 Mpc−3 towards the bright and faint end of our survey, respectively. By employing a non-parametric statistical test, and by comparing the full sample to subsamples in redshift bins, we find no supporting evidence for an evolving LAE LF over the probed redshift and luminosity range. Using a parametric maximum-likelihood technique we determine the best-fitting Schechter function parameters α = 1.84+04.2−0.41 and L∗[erg s−0.1] = 42.2−0.16+0.22 with the corresponding normalisation logϕ*[Mpc−3]= − 2.71. However, the dynamic range in Lyα luminosities probed by MUSE-Wide leads to a strong degeneracy between α and L*. Moreover, we find that a power-law parameterisation of the LF appears to be less consistent with the data compared to the Schechter function, even so when not excluding the X-Ray identified AGN from the sample. When correcting for completeness in the LAE LF determinations, we take into account that LAEs exhibit diffuse extended low surface brightness halos. We compare the resulting LF to one obtained by applying a correction assuming compact point-like emission. We find that the standard correction underestimates the LAE LF at the faint end of our survey by a factor of 2.5. Contrasting our results to the literature we find that at 42.2 ≤ log LLyα[erg s−1] ≲ 42.5 previous LAE LF determinations from narrow-band surveys appear to be affected by a similar bias.

Deep optical spectroscopic surveys of galaxies provide a unique opportunity to investigate rest-frame ultra-violet (UV) emission line properties of galaxies at z ∼ 2 − 4.5. Here we combine VLT/MUSE Guaranteed Time Observations of the Hubble Deep Field South, Ultra Deep Field, COSMOS, and several quasar fields with other publicly available data from VLT/VIMOS and VLT/FORS2 to construct a catalogue of He II λ1640 emitters at z ≳ 2. The deepest areas of our MUSE pointings reach a 3σ line flux limit of 3.1 × 10−19 erg s−1 cm−2. After discarding broad-line active galactic nuclei, we find 13 He II λ1640 detections from MUSE with a median MUV = −20.1 and 21 tentative He II λ1640 detections from other public surveys. Excluding Lyα, all except two galaxies in our sample show at least one other rest-UV emission line, with C III] λ1907, λ1909 being the most prominent. We use multi-wavelength data available in the Hubble legacy fields to derive basic galaxy properties of our sample through spectral energy distribution fitting techniques. Taking advantage of the high-quality spectra obtained by MUSE (∼10 − 30 h of exposure time per pointing), we use photo-ionisation models to study the rest-UV emission line diagnostics of the He II λ1640 emitters. Line ratios of our sample can be reproduced by moderately sub-solar photo-ionisation models, however, we find that including effects of binary stars lead to degeneracies in most free parameters. Even after considering extra ionising photons produced by extreme sub-solar metallicity binary stellar models, photo-ionisation models are unable to reproduce rest-frame He II λ1640 equivalent widths (∼0.2 − 10 Å), thus additional mechanisms are necessary in models to match the observed He II λ1640 properties.

The well-known quasar SDSS J095253.83+011421.9 (J0952+0114) at z=3.02 has one of the most peculiar spectra discovered so far, showing the presence of narrow Ly$\alpha$ and broad metal emission lines. Although recent studies have suggested that a Proximate Damped Ly$\alpha$ system (PDLA) causes this peculiar spectrum, the origin of the gas associated with the PDLA is unknown. Here we report the results of MUSE observations that reveal a new giant ($\approx$ 100 physical kpc) Lyman $\alpha$ nebula. The detailed analysis of the Ly$\alpha$ velocity, velocity dispersion, and surface brightness profiles suggests that the J0952+0114 Ly$\alpha$ nebula shares similar properties of other QSO nebulae previously detected with MUSE, implying that the PDLA in J0952+0144 is covering only a small fraction of the QSO emission solid angle. We also detected bright and spectrally narrow CIV$\lambda$1550 and HeII$\lambda$1640 extended emission around J0952+0114 with velocity centroids similar to the peak of the extended and central narrow Ly$\alpha$ emission. The presence of a peculiarly bright, unresolved, and relatively broad HeII$\lambda$1640 emission in the central region at exactly the same PDLA redshift hints at the possibility that the PDLA originates in a clumpy outflow with a bulk velocity of about 500 km/s. The smaller velocity dispersion of the large scale Ly$\alpha$ emission suggests that the high-speed outflow is confined to the central region. Lastly, the derived spatially resolved HeII/Ly$\alpha$ and CIV/Ly$\alpha$ maps show a positive gradient with the distance to the QSO hinting at a non-homogeneous ionization parameter distribution., Comment: 17 pages, 13 figures and 2 tables; Accepted for publication in ApJ

We use the MusE GAs FLOw and Wind (MEGAFLOW) survey to study the kinematics of extended disk-like structures of cold gas around $z\approx1$ star-forming galaxies. The combination of VLT/MUSE and VLT/UVES observations allows us to connect the kinematics of the gas measured through MgII quasar absorption spectroscopy to the kinematics and orientation of the associated galaxies constrained through integral field spectroscopy. Confirming previous results, we find that the galaxy-absorber pairs of the MEGAFLOW survey follow a strong bimodal distribution, consistent with a picture of MgII absorption being predominantly present in outflow cones and extended disk-like structures. This allows us to select a bona-fide sample of galaxy-absorber pairs probing these disks for impact parameters of 10-70 kpc. We test the hypothesis that the disk-like gas is co-rotating with the galaxy disks, and find that for 7 out of 9 pairs the absorption velocity shares the sign of the disk velocity, disfavouring random orbits. We further show that the data are roughly consistent with inflow velocities and angular momenta predicted by simulations, and that the corresponding mass accretion rates are sufficient to balance the star formation rates., 22 pages + Supplementary Appendix (15 pages); submitted to MNRAS; version after addressing referee's comments

We report the discovery of six spatially extended (10-100 kpc) line-emitting nebulae in the z=0.57 galaxy group hosting PKS0405-123, one of the most luminous quasars at z, 9 pages, 3 figures, 2 tables. Accepted to ApJ Letters

This article is an erratum for:[https://doi.org/10.1051/0004-6361/201834565]​​​​​​​

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

We present spatially resolved maps of six individually-detected Lyman alpha haloes (LAHs) as well as a first statistical analysis of the Lyman alpha (Lya) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies using MUSE. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lya line profile. Using a three-dimensional two-component model for the Lya emission, we measure the full width at half maximum (FWHM), the peak velocity shift and the asymmetry of the Lya line in the core and in the halo of 19 galaxies. We find that the Lya line shape is statistically different in the halo compared to the core for ~40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lya line both at the interstellar and circum-galactic scales. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies show broader, more redshifted and less asymmetric Lya lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lya lines. The generally broad and red line shapes found in the halo component suggests that the Lya haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lya line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs., Comment: 17+4 pages, 12+3 figures, 2+1 tables, accepted for publication in A&A

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., 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

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 MUSE (Multi Unit Spectroscopic Explorer) survey for dark galaxies fluorescently illuminated by quasars at $z>3$. Compared to previous studies which are based on deep narrow-band (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 (EW$_{0}$) distributions of the Ly$\alpha$ sources detected in proximity to the quasars and in control samples, we detect a clear correlation between the locations of high EW$_{0}$ objects and the quasars. This correlation is not seen in other properties such as Ly$\alpha$ luminosities or volume overdensities, suggesting the possible fluorescent nature of at least some of these objects. Among these, we find 6 sources without continuum counterparts and EW$_{0}$ limits larger than $240\,\mathrm{\AA}$ 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 previously detected candidates at $z\approx2.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., Comment: 28 pages, 13 figures, 5 tables, submitted to ApJ

International audience; Using an ultra-deep blind survey with the MUSE integral field spectrograph on the European Southern Observatory Very Large Telescope, we obtain spectroscopic redshifts to a depth never before explored: galaxies with observed magnitudes m AB 30-32. Specifically, we detect objects via Lyα emission at 2.9

We report the detection of extended Lyα haloes around 145 individual star-forming galaxies at redshifts 3 ≤ z ≤ 6 in the Hubble Ultra Deep Field observed with the Multi-Unit Spectroscopic Explorer (MUSE) at ESO-VLT. Our sample consists of continuum-faint (− 15 ≥ MUV ≥ −22) Lyα emitters (LAEs). Using a 2D, two-component (continuum-like and halo) decomposition of Lyα emission assuming circular exponential distributions, we measure scale lengths and luminosities of Lyα haloes. We find that 80% of our objects having reliable Lyα halo measurements show Lyα emission that is significantly more extended than the UV continuum detected by HST (by a factor ≈4 to >20). The median exponential scale length of the Lyα haloes in our sample is ≈4.5 kpc with a few haloes exceeding 10 kpc. By comparing the maximal detected extent of the Lyα emission with the predicted dark matter halo virial radii of simulated galaxies, we show that the detected Lyα emission of our selected sample of Lyα emitters probes a significant portion of the cold circum-galactic medium of these galaxies (>50% in average). This result therefore shows that there must be significant HI reservoirs in the circum-galactic medium and reinforces the idea that Lyα haloes are ubiquitous around high-redshift Lyα emitting galaxies. Our characterization of the Lyα haloes indicates that the majority of the Lyα flux comes from the halo (≈65%) and that their scale lengths seem to be linked to the UV properties of the galaxies (sizes and magnitudes). We do not observe a significant Lyα halo size evolution with redshift, although our sample for z> 5 is very small. We also explore the diversity of the Lyα line profiles in our sample and we find that the Lyα lines cover a large range of full width at half maximum (FWHM) from 118 to 512 km s-1. While the FWHM does not seem to be correlated to the Lyα scale length, most compact Lyα haloes and those that are not detected with high significance tend to have narrower Lyα profiles (-1). Finally, we investigate the origin of the extended Lyα emission but we conclude that our data do not allow us to disentangle the possible processes, i.e. scattering from star-forming regions, fluorescence, cooling radiation from cold gas accretion, and emission from satellite galaxies.

We present here a review of the latest results on the spatially-resolved analysis of the stellar populations and ionized gas of disk-dominated galaxies based on Calar Alto Legacy Integral Field Area (CALIFA) data. CALIFA is an ongoing integral field spectroscopy (IFS) survey of galaxies in the Local Universe (0.005 &lt, z &lt, 0.03) that has already obtained spectroscopic information up to (sim)2.5 (r_e) with a spatial resolution better than (sim)1 kpc for a total number of more than 600 galaxies of different morphological types, covering the color-magnitude diagram up to M(_{rm R}&lt, )18 mag. With nearly 2000 spectra obtained for each galaxy, CALIFA offers one of the best IFU datasets to study the star formation histories and chemical enrichment of galaxies. In this article, we focus on the main results from the analysis of the oxygen abundances based on the study of ionized gas in H II regions and individual spaxels and their relation to the global properties of galaxies, using an updated/revised dataset with more galaxies and ionized regions. In summary, we have confirmed previous published results indicating that: (1) the M-Z relation does not present a secondary relation to the star formation rate, when the abundance is measured at the effective radius, (2) the oxygen abundance presents a strong correlation with the stellar surface density (∑-Z relation), (3) the oxygen abundance profiles present three well-defined regimes: (i) an overall negative radial gradient between 0.5 and 2 (r_e), with a characteristic slope of (alpha_{O/H}) (sim)(-)0.1 dex/(r_e), (ii) a universal flattening beyond (&gt, )2 (r_e), and (iii) an inner drop at (&lt, )0.5 (r_e) that depends on mass, (4) the presence of bending in the surface brightness profile of disk galaxies is not clearly related to either the change in the shape of the oxygen abundance profile or the properties of the underlying stellar population. All of these results indicate that disk galaxies present an overall inside-out growth, with chemical enrichment and stellar mass growth tightly correlated and dominated by local processes and limited effects of radial mixing or global outflows. However, clear deviations are shown with respect to this simple scenario, which affect the abundance profiles in both the innermost and outermost regions of galaxies.

arXiv:1710.00826v1, We analyze the properties of a multiply imaged Lyα (Lyα) emitter at z = 5.75 identified through SHARDS Frontier Fields intermediate-band imaging of the Hubble Frontier Fields (HFF) cluster Abell 370. The source, A370-L57, has low intrinsic luminosity (M ∼ -16.5), steep UV spectral index (β = -2.4 ± 0.1), and extreme rest-frame equivalent width of Lyα. Two different gravitational lens models predict high magnification (μ ∼ 10-16) for the two detected counterimages, separated by 7″, while a predicted third counterimage (μ ∼ 3-4) is undetected. We find differences of ∼50% in magnification between the two lens models, quantifying our current systematic uncertainties. Integral field spectroscopy of A370-L57 with MUSE shows a narrow (FWHM = 204 ± 10 km s) and asymmetric Lyα profile with an integrated luminosity L(Lyα) ∼ 10 erg s. The morphology in the Hubble Space Telescope bands comprises a compact clump (r < 100 pc) that dominates the Lyα and continuum emission and several fainter clumps at projected distances ≲1 kpc that coincide with an extension of the Lyα emission in the SHARDS F823W17 and MUSE observations. The latter could be part of the same galaxy or an interacting companion. We find no evidence of a contribution from active galactic nuclei to the Lyα emission. Fitting of the spectral energy distribution with stellar population models favors a very young (t < 10 Myr), low-mass (), and metal-poor (Z ≲ 4 10) stellar population. Its modest star formation rate (SFR ∼ 1.0 yr) implies high specific SFR (sSFR ∼ 2.5 10 yr) and SFR density ( yr kpc). The properties of A370-L57 make it a good representative of the population of galaxies responsible for cosmic reionization., A.H.-C. and P.G.P.-G. acknowledge funding by the Spanish Ministry of Economy and Competitiveness (MINECO) under grants AYA2012-31277, AYA2015-70815-ERC, and AYA2015-63650-P. J.M.D. acknowledges support of the projects AYA2015-64508-P (MINECO/FEDER, UE), AYA2012-39475-C02-01, and the consolider project CSD2010-00064 funded by MINECO. D.J.L. and J.R. acknowledge support from the ERC starting grant 336736-CALENDS. R.A.M. acknowledges support by the Swiss National Science Foundation. A.A.-H. acknowledges funding by the MINECO grant AYA2015-64346-C2-1-P, which is partly funded by the FEDER program. H.D.S. acknowledges funding from the program ANR-T-ERC ASTROBRAIN. P.S. acknowledges funding from the Swiss National Science Foundation under grant P2GEP2_165426. J.M.R.E. acknowledges funding by the MINECO grant AYA2015-70498-C2-1-R. H.D. acknowledges financial support from the Spanish 2014 Ramón y Cajal program MINECO RYC-2014-15686.

Boxy/peanut bulges are considered to be part of the same stellar structure as bars and both could be linked through the buckling instability. The Milky Way is our closest example. The goal of this Letter is to determine if the mass assembly of the different components leaves an imprint in their stellar populations allowing the estimation the time of bar formation and its evolution. To this aim, we use integral field spectroscopy to derive the stellar age distributions, SADs, along the bar and disc of NGC 6032. The analysis clearly shows different SADs for the different bar areas. There is an underlying old (≥12 Gyr) stellar population for the whole galaxy. The bulge shows star formation happening at all times. The inner bar structure shows stars of ages older than 6 Gyr with a deficit of younger populations. The outer bar region presents an SAD similar to that of the disc. To interpret our results, we use a generic numerical simulation of a barred galaxy. Thus, we constrain, for the first time, the epoch of bar formation, the buckling instability period and the posterior growth from disc material. We establish that the bar of NGC 6032 is old, formed around 10 Gyr ago while the buckling phase possibly happened around 8 Gyr ago. All these results point towards bars being long-lasting even in the presence of gas., Monthly Notices of the Royal Astronomical Society, 470 (1), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

The radio-quiet quasar Q2059-360 at redshift $z=3.08$ is known to be close to a small Lyman $\alpha$ blob (LAB) and to be absorbed by a proximate damped Ly$\alpha$ (PDLA) system. Here, we present the Multi Unit Spectroscopic Explorer (MUSE) integral field spectroscopy follow-up of this quasi-stellar object (QSO). Our primary goal is to characterize this LAB in detail by mapping it both spatially and spectrally using the Ly$\alpha$ line, and by looking for high-ionization lines to constrain the emission mechanism. Combining the high sensitivity of the MUSE integral field spectrograph mounted on the Yepun telescope at ESO-VLT with the natural coronagraph provided by the PDLA, we map the LAB down to the QSO position, after robust subtraction of QSO light in the spectral domain. In addition to confirming earlier results for the small bright component of the LAB, we unveil a faint filamentary emission protruding to the south over about 80 pkpc (physical kpc); this results in a total size of about 120 pkpc. We derive the velocity field of the LAB (assuming no transfer effects) and map the Ly$\alpha$ line width. Upper limits are set to the flux of the N V $\lambda 1238-1242$, C IV $\lambda 1548-1551$, He II $\lambda 1640$, and C III] $\lambda 1548-1551$ lines. We have discovered two probable Ly$\alpha$ emitters at the same redshift as the LAB and at projected distances of 265 kpc and 207 kpc from the QSO; their Ly$\alpha$ luminosities might well be enhanced by the QSO radiation. We also find an emission line galaxy at $z=0.33$ near the line of sight to the QSO. This LAB shares the same general characteristics as the 17 others surrounding radio-quiet QSOs presented previously. However, there are indications that it may be centered on the PDLA galaxy rather than on the QSO., Comment: Accepted for publication in Astronomy & Astrophysics; 16 pages, 19 figures

We present a method to recover the gas-phase metallicity gradients from integral field spectroscopic (IFS) observations of barely resolved galaxies. We take a forward modelling approach and compare our models to the observed spatial distribution of emission line fluxes, accounting for the degrading effects of seeing and spatial binning. The method is flexible and is not limited to particular emission lines or instruments. We test the model through comparison to synthetic observations and use downgraded observations of nearby galaxies to validate this work. As a proof of concept we also apply the model to real IFS observations of high-redshift galaxies. From our testing we show that the inferred metallicity gradients and central metallicities are fairly insensitive to the assumptions made in the model and that they are reliably recovered for galaxies with sizes approximately equal to the half width at half maximum of the point-spread function. However, we also find that the presence of star forming clumps can significantly complicate the interpretation of metallicity gradients in moderately resolved high-redshift galaxies. Therefore we emphasize that care should be taken when comparing nearby well-resolved observations to high-redshift observations of partially resolved galaxies., Accepted for publication in MNRAS

We have conducted a two-layered spectroscopic survey (1'x1' ultra deep and 3'x3' deep regions) in the Hubble Ultra Deep Field (HUDF) with the Multi Unit Spectroscopic Explorer (MUSE). The combination of a large field of view, high sensitivity, and wide wavelength coverage provides an order of magnitude improvement in spectroscopically confirmed redshifts in the HUDF; i.e., 1206 secure spectroscopic redshifts for HST continuum selected objects, which corresponds to 15% of the total (7904). The redshift distribution extends well beyond z>3 and to HST/F775W magnitudes as faint as ~30 mag (AB, 1-sigma). In addition, 132 secure redshifts were obtained for sources with no HST counterparts that were discovered in the MUSE data cubes by a blind search for emission-line features. In total, we present 1338 high quality redshifts, which is a factor of eight increase compared with the previously known spectroscopic redshifts in the same field. We assessed redshifts mainly with the spectral features [OII] at z~20% up to 28-29 mag and ~27 mag, respectively. We used the determined redshifts to test continuum color selection (dropout) diagrams of high-z galaxies. The selection condition for F336W dropouts successfully captures ~80% of the targeted z~2.7 galaxies. However, for higher redshift selections (F435W, F606W, and F775W dropouts), the success rates decrease to ~20-40%. We empirically redefine the selection boundaries to make an attempt to improve them to ~60%. The revised boundaries allow bluer colors that capture Lya emitters with high Lya equivalent widths falling in the broadbands used for the color-color selection. Along with this paper, we release the redshift and line flux catalog., Comment: 26 pages, 26 figures, 4 tables, accepted for publication in A&A (MUSE UDF Series Paper II), the redshift catalogs will be available at the CDS

Galaxies in the local Universe are known to follow bimodal distributions in the global stellar population properties. We analyse the distribution of the local average stellar population ages of 654 053 sub-galactic regions resolved on ∼1 kpc scales in a volume-corrected sample of 394 galaxies, drawn from the Calar Alto Legacy Integral Field Area (CALIFA) DR3 integral-field-spectroscopy survey and complemented by Sloan Digital Sky Survey (SDSS) imaging. We find a bimodal local-age distribution, with an old and a young peak primarily due to regions in early-type galaxies and star-forming regions of spirals, respectively. Within spiral galaxies, the older ages of bulges and interarm regions relative to spiral arms support an internal age bimodality. Although regions of higher stellar mass surface density, μ*, are typically older, μ* alone does not determine the stellar population age and a bimodal distribution is found at any fixed μ*. We identify an ‘old ridge’ of regions of age ∼9 Gyr, independent of μ*, and a ‘young sequence’ of regions with age increasing with μ* from 1–1.5 to 4–5 Gyr. We interpret the former as regions containing only old stars, and the latter as regions where the relative contamination of old stellar populations by young stars decreases as μ* increases. The reason why this bimodal age distribution is not inconsistent with the unimodal shape of the cosmic-averaged star formation history is that (i) the dominating contribution by young stars biases the age low with respect to the average epoch of star formation, and (ii) the use of a single average age per region is unable to represent the full time extent of the star formation history of ‘young sequence’ regions., Monthly Notices of the Royal Astronomical Society, 468 (2), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We present the deepest study to date of the Lya luminosity function (LF) in a blank field using blind integral field spectroscopy from MUSE. We constructed a sample of 604 Lya emitters (LAEs) across the redshift range 2.91 < z < 6.64 using automatic detection software in the Hubble Ultra Deep Field. We calculate accurate total Lya fluxes capturing low surface brightness extended Lya emission now known to be a generic property of high-redshift star-forming galaxies. We simulated realistic extended LAEs to characterise the selection function of our samples, and performed flux-recovery experiments to test and correct for bias in our determination of total Lya fluxes. We find an accurate completeness correction accounting for extended emission reveals a very steep faint-end slope of the LF, alpha, down to luminosities of log10 L erg s^-1< 41.5, applying both the 1/Vmax and maximum likelihood estimators. Splitting the sample into three broad redshift bins, we see the faint-end slope increasing from -2.03+1.42-inf at z ~ 3.44 to -2.86+0.76-inf at z ~ 5.48, however no strong evolution is seen between the 68% confidence regions in L*-alpha parameter space. Using the Lya line flux as a proxy for star formation activity, and integrating the observed LFs, we find that LAEs' contribution to the cosmic SFRD rises with redshift until it is comparable to that from continuum-selected samples by z ~ 6. This implies that LAEs may contribute more to the star-formation activity of the early Universe than previously thought - any additional interglactic medium correction would act to further boost the Lya luminosities. Finally, assuming fiducial values for the escape of Lya and LyC radiation, and the clumpiness of the IGM, we integrated the maximum likelihood LF at 5.00 < z < 6.64 and find we require only a small extrapolation beyond the data (< 1 dex in L) for LAEs alone to maintain an ionised IGM at z ~ 6., Comment: 16 pages, 10 figures. Accepted for publication in A&A special issue: "The MUSE Hubble Ultra Deep Field Survey" (MUSE UDF Series Paper VI)

Emission signatures from galactic winds provide an opportunity to directly map the outflowing gas, but this is traditionally challenging because of the low surface brightness. Using deep observations (27 hours) of the Hubble Deep Field South from the Multi Unit Spectroscopic Explorer (MUSE) instrument, we identify signatures of an outflow in both emission and absorption from a spatially resolved galaxy at z = 1.29 with a stellar mass M* = 8 x 10^9 Msun, star formation rate SFR = 77 Msun/yr, and star formation rate surface brightness 1.6 Msun/kpc^2 within the [OII] half-light radius R_1/2,[OII] = 2.76 +- 0.17 kpc. From a component of the strong resonant MgII and FeII absorptions at -350 km/s, we infer a mass outflow rate that is comparable to the star formation rate. We detect non-resonant FeII* emission, at lambda 2626, 2612, 2396, and 2365, at 1.2-2.4-1.5-2.7 x 10^-18 egs s-1 cm-2 respectively. These flux ratios are consistent with the expectations for optically thick gas. By combining the four non-resonant FeII* emission lines, we spatially map the FeII* emission from an individual galaxy for the first time. The FeII* emission has an elliptical morphology that is roughly aligned with the galaxy minor kinematic axis, and its integrated half-light radius R_1/2,FeII* = 4.1 +- 0.4 kpc is 50% larger than the stellar continuum (R_1/2,* = 2.34 +- 0.17 kpc) or the [OII] nebular line. Moreover, the FeII* emission shows a blue wing extending up to -400 km/s, which is more pronounced along the galaxy minor kinematic axis and reveals a C-shaped pattern in a p-v diagram along that axis. These features are consistent with a bi-conical outflow., 13 pages, 6 figures, 4 tables, submitted to A&A - comments welcome

We present an updated version of the mass--metallicity relation (MZR) using integral field spectroscopy data obtained from 734 galaxies observed by the CALIFA survey. These unparalleled spatially resolved spectroscopic data allow us to determine the metallicity at the same physical scale ($\mathrm{R_{e}}$) for different calibrators. We obtain MZ relations with similar shapes for all calibrators, once the scale factors among them are taken into account. We do not find any significant secondary relation of the MZR with either the star formation rate (SFR) or the specific SFR for any of the calibrators used in this study, based on the analysis of the residuals of the best fitted relation. However we do see a hint for a (s)SFR-dependent deviation of the MZ-relation at low masses (M$, Comment: 22 pages, 9 figures, accepted for publishing in MNRAS

Understanding the nature of bulges in disc galaxies can provide important insights into the formation and evolution of galaxies. For instance, the presence of a classical bulge suggests a relatively violent history. In contrast, the presence of an inner disc instead (also referred to as a >pseudobulge>) indicates the occurrence of secular evolution processes in the main disc. However, we still lack criteria to effectively categorise bulges, limiting our ability to study their impact on the evolution of the host galaxies. Here we present a recipe to separate inner discs from classical bulges by combining four different parameters from photometric and kinematic analyses: the bulge Sérsic index n, the concentration index C, the Kormendy (1977, ApJ, 217, 406) relation and the inner slope of the radial velocity dispersion profile σ. With that recipe we provide a detailed bulge classification for a sample of 45 galaxies from the integral-field spectroscopic survey CALIFA. To aid in categorising bulges within these galaxies, we perform 2D image decomposition to determine bulge Sérsic index, bulge-to-total light ratio, surface brightness and effective radius of the bulge and use growth curve analysis to derive a new concentration index, C. We further extract the stellar kinematics from CALIFA data cubes and analyse the radial velocity dispersion profile. The results of the different approaches are in good agreement and allow a safe classification for approximately 95% of the galaxies. In particular, we show that our new >inner> concentration index performs considerably better than the traditionally used C when yielding the nature of bulges. We also found that a combined use of this index and the Kormendy relation gives a very robust indication of the physical nature of the bulge.© ESO, 2017., R.G.B. and R.G.D. acknowledge support from grants AYA2014-57490-P and JA-FQM-2828. R.A.M. acknowledges support by the Swiss National Science Foundation. I.M. acknowledges financial support from grants AYA2013-42227-P and AYA2016-76682-C3-1-P. SFS thanks the CONACYT-125180, DGAPA-IA100815 and DGAPA-IA101217 projects for providing him support in this study.

The C III] 1907,1909 emission doublet has been proposed as an alternative to Lyman-alpha in redshift confirmations of galaxies at z > 6 since it is not attenuated by the largely neutral intergalactic medium at these redshifts and is believed to be strong in the young, vigorously star-forming galaxies present at these early cosmic times. We present a statistical sample of 17 C III]-emitting galaxies beyond z~1.5 using 30 hour deep VLT/MUSE integral field spectroscopy covering 2 square arcminutes in the Hubble Deep Field South (HDFS) and Ultra Deep Field (UDF), achieving C III] sensitivities of ~2e-17 erg/s/cm^2 in the HDFS and ~7e-18 erg/s/cm^2 in the UDF. The rest-frame equivalent widths range from 2 to 19 Angstroms. These 17 galaxies represent ~3% of the total sample of galaxies found between 1.5 < z < 4. They also show elevated star formation rates, lower dust attenuation, and younger mass-weighted ages than the general population of galaxies at the same redshifts. Combined with deep slitless grism spectroscopy from the HST/WFC3 in the UDF, we can tie the rest-frame ultraviolet C III] emission to rest-frame optical emission lines, namely [O III] 5007, finding a strong correlation between the two. Down to the flux limits that we observe (~1e-18 erg/s/cm^2 with the grism data in the UDF), all objects with a rest-frame [O III] 4959,5007 equivalent width in excess of 250 Angstroms, the so-called Extreme Emission Line Galaxies, have detections of C III] in our MUSE data. More detailed studies of the C III]-emitting population at these intermediate redshifts will be crucial to understand the physical conditions in galaxies at early cosmic times and to determine the utility of C III] as a redshift tracer., 16 pages and 13 figures, accepted for publication in A&A

We present the first estimate of the Ly{\alpha} luminosity function using blind spectroscopy from the Multi Unit Spectroscopic Explorer, MUSE, in the Hubble Deep Field South. Using automatic source-detection software, we assemble a homogeneously-detected sample of 59 Ly{\alpha} emitters covering a flux range of -18.0 < log10 (F) < -16.3 (erg s^-1 cm^-2), corresponding to luminosities of 41.4 < log10 (L) < 42.8 (erg s^-1). As recent studies have shown, Ly{\alpha} fluxes can be underestimated by a factor of two or more via traditional methods, and so we undertake a careful assessment of each object's Ly{\alpha} flux using a curve-of-growth analysis to account for extended emission. We describe our self-consistent method for determining the completeness of the sample, and present an estimate of the global Ly{\alpha} luminosity function between redshifts 2.91 < z < 6.64 using the 1/Vmax estimator. We find the luminosity function is higher than many number densities reported in the literature by a factor of 2 - 3, although our result is consistent at the 1{\sigma} level with most of these studies. Our observed luminosity function is also in good agreement with predictions from semi-analytic models, and shows no evidence for strong evolution between the high- and low-redshift halves of the data. We demonstrate that one's approach to Ly{\alpha} flux estimation does alter the observed luminosity function, and caution that accurate flux assessments will be crucial in measurements of the faint end slope. This is a pilot study for the Ly{\alpha} luminosity function in the MUSE deep-fields, to be built on with data from the Hubble Ultra Deep Field which will increase the size of our sample by almost a factor of 10., Comment: 13 pages, 8 figures, submitted to MNRAS

Direct Lyα imaging of intergalactic gas at $z\sim 2$ has recently revealed giant cosmological structures around quasars, e.g., the Slug Nebula. Despite their high luminosity, the detection rate of such systems in narrow-band and spectroscopic surveys is less than 10%, possibly encoding crucial information on the distribution of gas around quasars and the quasar emission properties. In this study, we use the MUSE integral-field instrument to perform a blind survey for giant $\mathrm{Ly}\alpha $ nebulae around 17 bright radio-quiet quasars at $3\lt z\lt 4$ that does not suffer from most of the limitations of previous surveys. After data reduction and analysis performed with specifically developed tools, we found that each quasar is surrounded by giant $\mathrm{Ly}\alpha $ nebulae with projected sizes larger than 100 physical kiloparsecs and, in some cases, extending up to 320 kpc. The circularly averaged surface brightness profiles of the nebulae appear to be very similar to each other despite their different morphologies and are consistent with power laws with slopes $\approx -1.8$. The similarity between the properties of all these nebulae and the Slug Nebula suggests a similar origin for all systems and that a large fraction of gas around bright quasars could be in a relatively “cold” (T ∼ 10(4 )K) and dense phase. In addition, our results imply that such gas is ubiquitous within at least 50 kpc from bright quasars at $3\lt z\lt 4$ independently of the quasar emission opening angle, or extending up to 200 kpc for quasar isotropic emission.

M. Lyubenova et. al., In this paper, we describe how we use stellar dynamics information to constrain the shape of the stellar initial mass function (IMF) in a sample of 27 early-type galaxies from the CALIFA survey. We obtain dynamical and stellar mass-to-light ratios, ¿dyn and ¿*, over a homogenous aperture of 0.5 Re. We use the constraint ¿dyn¿¿* to test two IMF shapes within the framework of the extended MILES stellar population models. We rule out a single power-law IMF shape for 75 per cent of the galaxies in our sample. Conversely, we find that a double power-law IMF shape with a varying high-mass end slope is compatible (within 1¿) with 95 per cent of the galaxies. We also show that dynamical and stellar IMF mismatch factors give consistent results for the systematic variation of the IMF in these galaxies. © 2016, Published by Oxford University Press on behalf of the Royal Astronomical Society., This Paper is based on data obtained by the CALIFA survey, funded by the Spanish Ministry of Science under grant ICTS-2009-10, and the CAHA. IMN and JFB acknowledge funding from grant AYA2013-48226-C3-1-P from the Spanish Ministry of Economy and Competitiveness (MINECO) and, together with and GvdV, from the FP7 Marie Curie Actions via the ITN DAGAL (grant 289313). CJW acknowledges support through the Marie Curie Career Integration Grant 303912. Support for LG is provided by the Ministry of Economy, Development, and Tourism's Millennium Science Initiative through grant IC120009 awarded to The Millennium Institute of Astrophysics (MAS), and CONICYT through FONDECYT grant 3140566. RGD acknowledges support from AyA2014-57490-P. JMA acknowledges support from the ERC Starting Grant (SEDmorph; P.I. V. Wild).

The metallicity of a supernova progenitor, together with its mass, is one of the main parameters that can rule the progenitor's fate. We present the second study of nearby supernova (SN) host galaxies (0.005 10 dex) by targeted searches. We neither found evidence that the metallicity at the SN location differs from the average metallicity at the galactocentric distance of the SNe. By extending our SN sample with published metallicities at the SN location, we are able to study the metallicity distributions for all SN subtypes split into SN discovered in targeted and untargeted searches. We confirm a bias toward higher host masses and metallicities in the targeted searches. By combining data from targeted and untargeted searches, we found a sequence from higher to lower local metallicity: SN Ia, Ic, and II show the highest metallicity, which is significantly higher than those of SN Ib, IIb, and Ic-BL. Our results support the scenario according to which SN Ib result from binary progenitors. Additionally, at least part of the SN Ic are the result of single massive stars that were stripped of their outer layers by metallicity-driven winds. We studied several proxies of the local metallicity that are frequently used in the literature and found that the total host metallicity allows estimating the metallicity at the SN location with an accuracy better than 0.08 dex and very small bias. In addition, weak AGNs that cannot be seen in the total spectrum may weakly bias (by 0.04 dex) the metallicity estimate that is derived from the galaxy-integrated spectrum. © ESO, 2016., This work was partly funded by FCT with the research grant PTDC/CTE-AST/112582/2009. Support for L.G. is partially provided by FCT, by CONICYT through FONDECYT grant 3140566, and from the Ministry of Economy, Development, and Tourism's Millennium Science Initiative through grant IC12009, awarded to The Millennium Institute of Astrophysics (MAS). V.S. acknowledges financial support from Fundacao para a Ciencia e a Tecnologia (FCT) under program Ciencia 2008. C.J.W. acknowledges support through the Marie Curie Career Integration Grant 303912.

Photoionization models of H ii regions require as input a description of the ionizing spectral energy distribution (SED) and of the gas distribution, in terms of ionization parameter U and chemical abundances (e.g., O/H and N/O).A strong degeneracy exists between the hardness of the SED and U, which in turn leads to high uncertainties in the determination of the other parameters, including abundances. One way to resolve the degeneracy is to fix one of the parameters using additional information. For each of the ~20 000 sources of the CALIFA H ii regions catalog, a grid of photoionization models is computed assuming the ionizing SED to be described by the underlying stellar population obtained from spectral synthesis modeling. The ionizing SED is then defined as the sum of various stellar bursts of different ages and metallicities. This solves the degeneracy between the shape of the ionizing SED and U. The nebular metallicity (associated with O/H) is defined using the classical strong line method O3N2 (which gives our models the status of >hybrids>). The remaining free parameters are the abundance ratio N/O and the ionization parameter U, which are determined by looking for the model fitting [N ii]/Hα and [O iii]/Hβ. The models are also selected to fit [O ii]/Hβ. This process leads to a set of ~3200 models that reproduce the three observations simultaneously. We find that the regions associated with young stellar bursts (i.e., ionized by OB stars) are affected by leaking of ionizing photons,the proportion of escaping photons having a median of 80%. The set of photoionization models satisfactorily reproduces the electron temperature derived from the [O iii]λ4363/5007 line ratio. We determine new relations between the nebular parameters, like the ionization parameter U and the [O ii]/[O iii] or [S ii]/[S iii] line ratios. A new relation between N/O and O/H is obtained, mostly compatible with previous empirical determinations (and not with previous results obtained using photoionization models). A new relation between U and O/H is also determined. All the models are publicly available on the Mexican Million Models database 3MdB. © 2016 ESO., Part of the results presented here have been obtained using computers from the CONACyT project CB2010/153985 and the UNAM-PAPIIT project 107215. The computers >Tychos> (Posgrado en Astrofisica-UNAM, Instituto de Astronomia-UNAM and PNPC-CONACyT) have been used for this research. G.D.I. gratefully acknowledges support from the Mexican CONACYT grant CB-2014-241732. S.F.S. thanks the CONACYT-125180 and DGAPA-IA100815 projects for providing support during this study. Support for L.G. is provided by the Ministry of Economy, Development, and Tourism's Millennium Science Initiative through grant IC120009, awarded to The Millennium Institute of Astrophysics, MAS. L.G. acknowledges support by CONICYT through FONDECYT grant 3140566.

Non-resonant FeII* 2365, 2396, 2612, 2626 emission can potentially trace galactic winds in emission and provide useful constraints to wind models. From the 3'x3' mosaic of the Hubble Ultra Deep Field (UDF) obtained with the VLT/MUSE integral field spectrograph, we identify a statistical sample of 40 FeII* emitters and 50 MgII 2796, 2803 emitters from a sample of 271 [OII] 3726, 3729 emitters with reliable redshifts from z = 0.85 - 1.5 down to 2E-18 (3 sigma) ergs/s/cm^2 (for [OII]), covering the stellar mass range 10^8 - 10^11 Msun. The FeII* and MgII emitters follow the galaxy main sequence, but with a clear dichotomy. Galaxies with masses below 10^9 Msun and star formation rates (SFRs) of 10 Msun/year have FeII* emission without accompanying MgII emission. Between these two regimes, galaxies have both MgII and FeII* emission, typically with MgII P-Cygni profiles. Indeed, the MgII profile shows a progression along the main sequence from pure emission to P-Cygni profiles to strong absorption, due to resonant trapping. Combining the deep MUSE data with HST ancillary information, we find that galaxies with pure MgII emission profiles have lower star formation rate surface densities than those with either MgII P-Cygni profiles or FeII* emission. These spectral signatures produced through continuum scattering and fluorescence, MgII P-Cygni profiles and FeII* emission, are better candidates for tracing galactic outflows than pure MgII emission, which may originate from HII regions. We compare the absorption and emission rest-frame equivalent widths for pairs of FeII transitions to predictions from outflow models and find that the observations consistently have less total re-emission than absorption, suggesting either dust extinction or non-isotropic outflow geometries., 16 pages, 12 figures, 3 tables, accepted for publication in A&A (14 October 2017)

We present rest-frame Lyα equivalent widths (E W 0 ) of 417 Lyα emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at 2.9 Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured E W 0 values taking into account extended Lyα emission and UV continuum slopes (β ). Our LAEs reach unprecedented depths, both in Lyα luminosities and UV absolute magnitudes, from log (L Lyα /erg s-1 ) ~ 41.0 to 43.0 and from M UV ~ −16 to −21 (0.01−1.0 L * z =3 ). The E W 0 values span the range of ~ 5 to 240 A or larger, and their distribution can be well fitted by an exponential law N = N 0 exp(−E W 0 /w 0 ). Owing to the high dynamic range in M UV , we find that the scale factor, w 0 , depends on M UV in the sense that including fainter M UV objects increases w 0 , i.e., the Ando effect. The results indicate that selection functions affect the E W 0 scale factor. Taking these effects into account, we find that our w 0 values are consistent with those in the literature within 1σ uncertainties at 2.9 UV and L Lyα . Interestingly, we find 12 objects with E W 0 > 200 A above 1σ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak Civλ 1549 emission line. For the remaining 10 very large E W 0 LAEs, we find that the E W 0 values can be reproduced by young stellar ages ( ⊙ ). Otherwise, at least part of the Lyα emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling.

We tested the performance of photometric redshifts for galaxies in the Hubble Ultra Deep field down to 30th magnitude. We compared photometric redshift estimates from three spectral fitting codes from the literature (EAZY, BPZ and BEAGLE) to high quality redshifts for 1227 galaxies from the MUSE integral field spectrograph. All these codes can return photometric redshifts with bias |Dzn|=|z-z_phot|/(1+z)3 they are systematically biased high by up to Dzn = 0.05, an offset that can in part be explained by adjusting the amount of intergalactic absorption applied. In agreement with previous studies we find little difference in the performance of the different codes, but in contrast to those we find that adding extensive ground-based and IRAC photometry actually can worsen photo-z performance for faint galaxies. We find an outlier fraction, defined through |Dzn|>0.15, of 8% for BPZ and 10% for EAZY and BEAGLE, and show explicitly that this is a strong function of magnitude. While this outlier fraction is high relative to numbers presented in the literature for brighter galaxies, they are very comparable to literature results when the depth of the data is taken into account. Finally, we demonstrate that while a redshift might be of high confidence, the association of a spectrum to the photometric object can be very uncertain and lead to a contamination of a few percent in spectroscopic training samples that do not show up as catastrophic outliers, a problem that must be tackled in order to have sufficiently accurate photometric redshifts for future cosmological surveys., Comment: 24 pages and 25 figures, accepted for A&A

Astronomy and Astrophysics 559 (2013): A114 reproduced with permission from Astronomy and Astrophysics, The use of integral field spectroscopy is since recently allowing to measure the emission line fluxes of an increasingly large number of star-forming galaxies, both locally and at high redshift. Many studies have used these fluxes to derive the gas-phase metallicity of the galaxies by applying the so-called strong-line methods. However, the metallicity indicators that these datasets use were empirically calibrated using few direct abundance data points (Te-based measurements). Furthermore, a precise determination of the prediction intervals of these indicators is commonly lacking in these calibrations. Such limitations might lead to systematic errors in determining the gas-phase metallicity, especially at high redshift, which might have a strong impact on our understanding of the chemical evolution of the Universe. The main goal of this study is to review the most widely used empirical oxygen calibrations, O3N2 and N2, by using newdirect abundance measurements. We pay special attention to (1) the expected uncertainty of these calibrations as a function of the index value or abundance derived and (2) the presence of possible systematic offsets. This is possible thanks to the analysis of the most ambitious compilation of Te-based H ii regions to date. This new dataset compiles the Te-based abundances of 603 H ii regions extracted from the literature but also includes new measurements from the CALIFA survey. Besides providing new and improved empirical calibrations for the gas abundance, we also present a comparison between our revisited calibrations with a total of 3423 additional CALIFA H ii complexes with abundances derived using the ONS calibration from the literature. The combined analysis of T e-based and ONS abundances allows us to derive their most accurate calibration to date for both the O3N2 and N2 single-ratio indicators, in terms of all statistical significance, quality, and coverage of the parameters space. In particular, we infer that these indicators show shallower abundance dependencies and statistically significant offsets compared to others'. The O3N2 and N2 indicators can be empirically applied to derive oxygen abundances calibrations from either direct abundance determinations with random errors of 0.18 and 0.16, respectively, or from indirect ones (but based on a large amount of data), reaching an average precision of 0.08 and 0.09 dex (random) and 0.02 and 0.08 dex (systematic; compared to the direct estimations), respectively, R.A. Marino is funded by the Spanish program of International Campus of Excellence Moncloa (CEI). D. Mast thank the Plan Nacional de Investigación y Desarrollo funding programs, AYA2012-31935 of the Spanish Ministerio de Economía y Competitividad, for the support given to this project. S.F.S thanks the the Ramón y Cajal project RyC-2011-07590 of the spanish Ministerio de Economía y Competitividad, for the support giving to this project. F.F.R.O. acknowledges the Mexican National Council for Science and Technology (CONACYT) for financial support under the program Estancias Postdoctorales y Sabáticas al Extranjero para la Consolidación de Grupos de Investigación, 2010-2012. We acknowledge financial support for the ESTALLIDOS collaboration by the Spanish Ministerio de Ciencia e Innovación under grant AYA2010- 21887-C04-03. BG-L also acknowledges support from the Spanish Ministerio de Economía y Competitividad (MINECO) under grant AYA2012- 39408-C02-02. J.F.-B. acknowledges financial support from the Ramón y Cajal Program and grant AYA2010-21322-C03-02 from the Spanish Ministry of Economy and Competitiveness (MINECO), as well as to the DAGAL network from the People’s Program (Marie Curie Actions) of the European Union’s Seventh Framework Program FP7/2007-2013/ under REA grant agreement number PITN-GA-2011-289313. CK has been funded by project AYA2010-21887 from the Spanish PNAYA. P.P. acknowledges support by the Fundação para a Ciência e a Tecnologia (FCT) under project FCOMP-01-0124-FEDER-029170 (Reference FCT PTDC/FIS-AST/3214/2012), funded by FCT-MEC (PIDDAC) and FEDER (COMPETE). R.M.G.D. and R.G.B. also acknowledge support from the Spanish Ministerio de Economía y Competitividad (MINECO) under grant AyA2010-15081. V.S., L.G., and A.M.M. acknowledge financial support from the Fundação para a Ciência e a Tecnologia (FCT) under program Ciência 2008 and the research grant PTDC/CTE-AST/112582/2009

The velocity function is a fundamental observable statistic of the galaxy population, similarly impor- tant as the luminosity function, but much more difficult to measure. In this work we present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km/s for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population (60 < v_circ < 170 km/s), we use the HIPASS velocity function. For the massive galaxy population (170 < v_circ < 320 km/s), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume- corrected galaxy distribution functions. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations.

Truncations in the stellar population at the edges of disk galaxies are thought to be a common morphological feature (e.g., Erwin et al. 2005; and more recently Marino et al. 2016). In fact, using imaging data from the SDSS, Pohlen & Trujillo (2006) showed that only ~ 10% of face-on to intermediate inclined, nearby, late-type (Sb-Sdm) spiral galaxies have a normal/standard purely exponential disk down to the noise limit. In situations like these, the simultaneous fit of two lines, joined or not at an intermediate point (the break radius), constitutes a natural step towards the modelling of radial variation in surface brightness, metallicity, or any other relevant parameter. This work shows the results of simple simulations in which the simultaneous fit to two joined lines is compared to the simultaneous fit of two independent lines (i.e., two lines that do not necessarily coincide at an intermediate point), and also to the traditional single ordinary least squares fit. These simulations reveal some biases that should be taken into account when facing these kind of fitting procedures.