The standard cosmological model ($\Lambda$CDM) predicts the existence of the cosmic web: a distribution of matter into sheets and filaments connecting massive halos. However, observational evidence has been elusive due to the low surface brightness of the filaments. Recent deep MUSE/VLT data and upcoming observations offer a promising avenue for Ly$\alpha$ detection, motivating the development of modern theoretical predictions. We use hydrodynamical cosmological simulations run with the AREPO code to investigate the potential detectability of large-scale filaments, excluding contributions from the halos embedded in them. We focus on filaments connecting massive ($M_{200c}\sim(1-3)\times10^{12} M_\odot$) halos at z=3, and compare different simulation resolutions, feedback levels, and mock-image pixel sizes. We find increasing simulation resolution does not substantially improve detectability notwithstanding the intrinsic enhancement of internal filament structure. By contrast, for a MUSE integration of 31 hours, including feedback increases the detectable area by a factor of $\simeq$5.5 on average compared with simulations without feedback, implying that even the non-bound components of the filaments have substantial sensitivity to feedback. Degrading the image resolution from the native MUSE scale of (0.2")$^2$ per pixel to (5.3")$^2$ apertures has the strongest effect, increasing the detectable area by a median factor of $\simeq$200 and is most effective when the size of the pixel roughly matches the width of the filament. Finally, we find the majority of Ly$\alpha$ emission is due to electron impact collisional excitations, as opposed to radiative recombination.

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.

We present velocities of galactic outflows in seven star-forming galaxies at $z=$5-6 with stellar masses of $M_* \sim10^{10.1}\ \rm{M_\odot}$. Although it is challenging to observationally determine the outflow velocities, we overcome this by using ALMA [CII]158 $\mu$m emission lines for systemic velocities and deep Keck spectra with metal absorption lines for velocity profiles available to date. We construct a composite Keck spectrum of the galaxies at $z=$5-6 with the [CII]-systemic velocities, and fit outflow-line profiles to the SiII1260, CII1335, and SiIV1394,1403 absorption lines in the composite spectrum. We measure the maximum (90%) and central outflow velocities to be $v_\rm{max}=700^{+180}_{-110}\ \rm{km\ s^{-1}}$ and $v_\rm{out}= 400^{+100}_{-150}\ \rm{km\ s^{-1}}$ on average, respectively, showing no significant differences between the outflow velocities derived with the low to high-ionization absorption lines. For $M_* \sim10^{10.1}\ \rm{M_\odot}$, we find that the $v_\rm{max}$ value of our $z=$5-6 galaxies is 3 times higher than those of $z\sim0$ galaxies and comparable to $z\sim2$ galaxies. Estimating the halo circular velocity $v_\rm{cir}$ from the stellar masses and the abundance matching results, we investigate a $v_\rm{max}$-$v_\rm{cir}$ relation. Interestingly, $v_\rm{max}$ for galaxies with $M_*=10^{10.0-10.8}\ \rm{M_\odot}$ shows a clear positive correlation with $v_\rm{cir}$ and/or the galaxy star formation rate over $z=$0-6 with a small scatter of $\simeq \pm 0.1$ dex, which is in good agreement with theoretical predictions. This positive correlation suggests that the outflow velocity is physically related to the halo circular velocity, and that the redshift evolution of $v_\rm{max}$ at fixed $M_*$ is explained by the increase in $v_\rm{cir}$ toward high redshift.

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 ~ 3.3 Ly\alpha emitters (LAEs) with a median Ly\alpha flux (luminosity) of ~10^{-17} erg cm^{-2} s^{-1} (~10^{42} erg s^{-1}) and a median unobscured star formation rate (SFR) of ~1.3 Msun/yr, 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, Voffset = 169 \pm 10 km/s, between the Ly\alpha emission peak and the systemic redshift. The observed Voffset is lower by factors of ~1.4 and ~2.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 Voffset 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 Voffset is correlated with SFR with a sub-linear scaling relation, Voffset ~ SFR^(0.16 \pm 0.03). Adopting the mass scaling for main sequence galaxies, such a relation suggests that Voffset scales with the circular velocity of the dark matter halos hosting the LAEs.

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.

It is still a challenge to assess the merger fraction of galaxies at different cosmic epochs in order to probe the evolution of their mass assembly. Using the Illustris cosmological simulations, we investigate the relation between the separation of galaxies in a pair, both in velocity and projected spatial separation space, and the probability that these interacting galaxies will merge in the future. From this analysis, we propose a new set of criteria to select close pairs of galaxies along with a new corrective term to be applied to the computation of the galaxy merger fraction. We then probe the evolution of the major and minor merger fraction using the latest MUSE deep observations over the HUDF, HDFS, COSMOS-Gr30 and Abell 2744 regions. From a parent sample of 2483 galaxies with spectroscopic redshifts, we identify 366 close pairs spread over a large range of redshifts ($0.2<z<6$) and stellar masses ($10^7-10^{11}M_{\odot}$). Using the stellar mass ratio between the secondary and primary galaxy as a proxy to split the sample into major, minor and very minor mergers, we found a total of 183 major, 142 minor and 47 very minor close pairs corresponding to a mass ratio range of 1:1-1:6, 1:6-1:100 and lower than 1:100, respectively. Due to completeness issues, we do not consider the very minor pairs in the analysis. Overall, the major merger fraction increases up to $z\approx 2-3$ reaching 25% for pairs with the most massive galaxy with a stellar mass $M^*\geq 10^{9.5}M_{\odot}$. Beyond this redshift, the fraction decreases down to $\sim 5$% at $z\approx 6$. The evolution of the minor merger fraction is roughly constant with cosmic time, with a fraction of 20% at $z<3$ and a slow decrease between $3\leq z \leq6$ to 8-13%.

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.

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.

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 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<1. The discovery is enabled by the Multi Unit Spectroscopic Explorer (MUSE) and provides tantalizing evidence connecting large-scale gas streams with nuclear activity on scales of <10 proper kpc (pkpc). One of the nebulae exhibits a narrow, filamentary morphology extending over 50 pkpc toward the quasar with narrow internal velocity dispersion (50 km/s) and is not associated with any detected galaxies, consistent with a cool intragroup medium (IGrM) filament. Two of the nebulae are 10 pkpc North and South of the quasar with tidal arm like morphologies. These two nebulae, along with a continuum emitting arm extending 60 pkpc from the quasar are signatures of interactions which are expected to redistribute angular momentum in the host interstellar medium (ISM) to facilitate star formation and quasar fueling in the nucleus. The three remaining nebulae are among the largest and most luminous [O III] emitting `blobs' known (1400-2400 pkpc^2) and correspond both kinematically and morphologically with interacting galaxy pairs in the quasar host group, consistent with arising from stripped ISM rather than large-scale quasar outflows. The presence of these large- and small-scale nebulae in the vicinity of a luminous quasar bears significantly on the effect of large-scale environment on galaxy and black hole fueling, providing a natural explanation for the previously known correlation between quasar luminosity and cool circumgalactic medium (CGM).

Star-forming galaxies have been found to follow a relatively tight relation between stellar mass ($M_{*}$) and star formation rate (SFR), dubbed the `star formation sequence'. A turnover in the sequence has been observed, where galaxies with $M_{*} < 10^{10} {\rm M}_{\odot}$ follow a steeper relation than their higher mass counterparts, suggesting that the low-mass slope is (nearly) linear. In this paper, we characterise the properties of the low-mass end of the star formation sequence between $7 \leq \log M_{*}[{\rm M}_{\odot}] \leq 10.5$ at redshift $0.11 < z < 0.91$. We use the deepest MUSE observations of the Hubble Ultra Deep Field and the Hubble Deep Field South to construct a sample of 179 star-forming galaxies with high signal-to-noise emission lines. Dust-corrected SFRs are determined from H$\beta$ $\lambda 4861$ and H$\alpha$ $\lambda 6563$. We model the star formation sequence with a Gaussian distribution around a hyperplane between $\log M_{*}$, $\log {\rm SFR}$, and $\log (1+z)$, to simultaneously constrain the slope, redshift evolution, and intrinsic scatter. We find a sub-linear slope for the low-mass regime where $\log {\rm SFR}[{\rm M}_{\odot}/{\rm yr}] = 0.83^{+0.07}_{-0.06} \log M_{*}[{\rm M}_{\odot}] + 1.74^{+0.66}_{-0.68} \log (1+z)$, increasing with redshift. We recover an intrinsic scatter in the relation of $\sigma_{\rm intr} = 0.44^{+0.05}_{-0.04}$ dex, larger than typically found at higher masses. As both hydrodynamical simulations and (semi-)analytical models typically favour a steeper slope in the low-mass regime, our results provide new constraints on the feedback processes which operate preferentially in low-mass halos.

The physical origin of the near-ultraviolet MgII emission remains an under-explored domain, contrary to more typical emission lines detected in the spectra of star-forming galaxies. We explore the nebular and physical properties for a sample of 381 galaxies between 0.70 < z < 2.34 drawn from the MUSE Hubble Ultra Deep Survey. The spectra of these galaxies show a wide variety of profiles of the MgII 2796,2803 resonant doublet, from absorption to emission. We present a study on the main drivers for the detection of MgII emission in galaxy spectra. By exploiting photoionization models we verified that the emission-line ratios observed in galaxies with MgII in emission are consistent with nebular emission from HII regions. From a simultaneous analysis of MUSE spectra and ancillary HST information via spectral energy distribution (SED) fitting, we find that galaxies with MgII in emission have lower stellar masses, smaller sizes, bluer spectral slopes and lower optical depth than those with absorption. This leads us to suggest that MgII emission is a potential tracer of physical conditions not merely related to those of the ionized gas. We show that these differences in MgII emission/absorption can be explained in terms of a higher dust and neutral gas content in the interstellar medium (ISM) of galaxies showing MgII in absorption, confirming the extreme sensitivity of MgII to the presence of the neutral ISM. We conclude with an analogy between the MgII doublet and the Ly-alpha line, due to their resonant nature. Further investigations with current and future facilities, including JWST, are promising as the detection of MgII emission and its potential connection with Ly-alpha could provide new insights on the ISM content in the early Universe.

The Lyman alpha (lya) line of Hydrogen is a prominent feature in the spectra of star-forming galaxies, usually redshifted by a few hundreds of km/s compared to the systemic redshift. This large offset hampers follow-up surveys, galaxy pair statistics and correlations with quasar absorption lines when only lya is available. We propose diagnostics that can be used to recover the systemic redshift directly from the properties of the lya line profile. We use spectroscopic observations of Lyman-Alpha Emitters (LAEs) for which a precise measurement of the systemic redshift is available. Our sample contains 13 sources detected between z~3 and z~6 as part of various Multi Unit Spectroscopic Explorer (MUSE) Guaranteed Time Observations (GTO). We also include a compilation of spectroscopic lya data from the literature spanning a wide redshift range (z~0-8). First, restricting our analysis to double-peaked lya spectra, we find a tight correlation between the velocity offset of the red peak with respect to the systemic redshift, Vpeak, and the separation of the peaks. Secondly, we find a correlation between Vpeak and the full width at half maximum of the lya line. Fitting formulas, to estimate systemic redshifts of galaxies with an accuracy of +-100 km/s when only the lya emission line is available, are given for the two methods.

We present the "SINS/zC-SINF AO survey" of 35 star-forming galaxies, the largest sample with deep adaptive optics-assisted (AO) near-infrared integral field spectroscopy at z~2. The observations, taken with SINFONI at the Very Large Telescope, resolve the Ha and [NII] line emission and kinematics on scales of ~1.5 kpc. In stellar mass, star formation rate, rest-optical colors and size, the AO sample is representative of its parent seeing-limited sample and probes the massive (M* ~ 2x10^9 - 3x10^11 Msun), actively star-forming (SFR ~ 10-600 Msun/yr) part of the z~2 galaxy population over a wide range in colors ((U-V)_rest ~ 0.15-1.5 mag) and half-light radii (R_e,H ~ 1-8.5 kpc). The sample overlaps largely with the "main sequence" of star-forming galaxies in the same redshift range to a similar K_AB = 23 magnitude limit; it has ~0.3 dex higher median specific SFR, ~0.1 mag bluer median (U-V)_rest color, and ~10% larger median rest-optical size. We describe the observations, data reduction, and extraction of basic flux and kinematic properties. With typically 3-4 times higher resolution and 4-5 times longer integrations (up to 23hr) than the seeing-limited datasets of the same objects, the AO data reveal much more detail in morphology and kinematics. The now complete AO observations confirm the majority of kinematically-classified disks and the typically elevated disk velocity dispersions previously reported based on subsets of the data. We derive typically flat or slightly negative radial [NII]/Ha gradients, with no significant trend with global galaxy properties, kinematic nature, or the presence of an AGN. Azimuthal variations in [NII]/Ha are seen in several sources and are associated with ionized gas outflows, and possible more metal-poor star-forming clumps or small companions. [Abridged]

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.

We provide, for the first time, robust observational constraints on the galaxy major merger fraction up to $z\approx 6$ using spectroscopic close pair counts. Deep Multi Unit Spectroscopic Explorer (MUSE) observations in the Hubble Ultra Deep Field (HUDF) and Hubble Deep Field South (HDF-S) are used to identify 113 secure close pairs of galaxies among a parent sample of 1801 galaxies spread over a large redshift range ($0.2<z<6$) and stellar masses ($10^7-10^{11} M_\odot$), thus probing about 12 Gyr of galaxy evolution. Stellar masses are estimated from spectral energy distribution (SED) fitting over the extensive UV-to-NIR HST photometry available in these deep Hubble fields, adding Spitzer IRAC bands to better constrain masses for high-redshift ($z\geqslant 3$) galaxies. These stellar masses are used to isolate a sample of 54 major close pairs with a galaxy mass ratio limit of 1:6. Among this sample, 23 pairs are identified at high redshift ($z\geqslant 3$) through their Ly$\alpha$ emission. The sample of major close pairs is divided into five redshift intervals in order to probe the evolution of the merger fraction with cosmic time. Our estimates are in very good agreement with previous close pair counts with a constant increase of the merger fraction up to $z\approx 3$ where it reaches a maximum of 20%. At higher redshift, we show that the fraction slowly decreases down to about 10% at $z\approx6$. The sample is further divided into two ranges of stellar masses using either a constant separation limit of $10^{9.5} M_\odot$ or the median value of stellar mass computed in each redshift bin. Overall, the major close pair fraction for low-mass and massive galaxies follows the same trend. These new, homogeneous, and robust estimates of the major merger fraction since $z\approx6$ are in good agreement with recent predictions of cosmological numerical simulations.

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 <1 Msun/year have MgII emission without accompanying FeII* emission, whereas galaxies with masses above 10^10 Msun and SFRs >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.

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)<0.05 down to F775W=30 and spectroscopic incompleteness is unlikely to strongly modify this statement. We have, however, identified clear systematic biases in the determination of photometric redshifts: in the 0.4<z<1.5 range, photometric redshifts are systematically biased low by as much as Dzn=-0.04 in the median, and at 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.

Star-forming galaxies (SFGs) are forming stars at a regular pace, forming the so-called main sequence (MS). However, all studies of their gas content show that their gas reservoir ought to be depleted in 0.5-2 Gyr. Thus, SFGs are thought to be fed by the continuous accretion of intergalactic gas in order to sustain their star-formation activity. However, direct observational evidence for this accretion phenomenon has been elusive. Theoretically, the accreted gas coming from the intergalactic medium is expected to orbit about the halo, delivering not just fuel for star-formation but also angular momentum to the galaxy. This accreting material is thus expected to form a gaseous structure that should be co-rotating with the host once at $r<0.3\;R_{\rm vir}$ or $r<10-30$ kpc. Because of the rough alignment between the star-forming disk and this extended gaseous structure, the accreting material can be most easily detected with the combination of background quasars and integral field units (IFUs). In this chapter, accretion studies using this technique are reviewed.

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.

The physical properties of galactic winds are one of the keys to understand galaxy formation and evolution. These properties can be constrained thanks to background quasar lines of sight (LOS) passing near star-forming galaxies (SFGs). We present the first results of the MusE GAs FLOw and Wind (MEGAFLOW) survey obtained of 2 quasar fields which have 8 MgII absorbers of which 3 have rest-equivalent width greater than 0.8 \AA. With the new Multi Unit Spectroscopic Explorer (MUSE) spectrograph on the Very Large Telescope (VLT), we detect 6 (75$\%$) MgII host galaxy candidates withing a radius of 30 arcsec from the quasar LOS. Out of these 6 galaxy--quasar pairs, from geometrical arguments, one is likely probing galactic outflows, two are classified as "ambiguous", two are likely probing extended gaseous disks and one pair seems to be a merger. We focus on the wind$-$pair and constrain the outflow using a high resolution quasar spectra from Ultraviolet and Visual Echelle Spectrograph (UVES). Assuming the metal absorption to be due to gas flowing out of the detected galaxy through a cone along the minor axis, we find outflow velocities of the order of $\approx$ 150 km/s (i.e. smaller than the escape velocity) with a loading factor, $\eta =\dot M_{\rm out}/$SFR, of $\approx$ 0.7. We see evidence for an open conical flow, with a low-density inner core. In the future, MUSE will provide us with about 80 multiple galaxy$-$quasar pairs in two dozen fields.

(Abridged) We make use of the deepest VLT/MUSE observations performed so far on the Hubble Deep Field South (HDFS) to characterize the low-mass (< $10^{10}$M$_\odot$) galaxy population at intermediate redshift. We identify a sample of 28 spatially-resolved emission-line galaxies in the deep (27h integration time) MUSE data cube, spread over a redshift interval of 0.2 < z < 1.4. The public HST images and multi-band photometry over the HDFS are used to constrain the stellar mass and star formation rate (SFR) of the galaxies and to perform a morphological analysis. We derive the resolved ionized gas properties of these galaxies from the MUSE data and model the disk (both in 2D and with GalPaK$^{\rm 3D}$) to retrieve their intrinsic gas kinematics. We build a sample of resolved emission-line galaxies of much lower stellar mass and SFR (by $\sim$1-2 orders of magnitude) than previous 3D spectroscopic surveys. Most of the spatially-resolved MUSE-HDFS galaxies have gas kinematics consistent with disk-like rotation, but about 20% have velocity dispersions larger than the rotation velocities, and 30% are part of a close pair and/or show clear signs of recent gravitational interactions. In the high-mass regime, the MUSE-HDFS galaxies follow the Tully-Fisher relation defined from previous surveys in a similar redshift range. This scaling relation extends also to lower masses/velocities but with a higher dispersion. The MUSE-HDFS galaxies follow the scaling relations defined in the local universe between the specific angular momentum and the stellar mass. However, we find that intermediate-redshift star-forming galaxies fill a continuum transition from the spiral to elliptical local scaling relations, according to the dynamical state of the gas. This indicates that some galaxies may lose their angular momentum and become dispersion-dominated prior to becoming passive.

Direct Ly $\alpha$ imaging of intergalactic gas at $z\sim2$ has recently revealed giant cosmological structures around quasars, e.g. the Slug Nebula (Cantalupo et al. 2014). 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 Ly $\alpha$ nebulae around 17 bright radio-quiet quasars at $3<z<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 Ly $\alpha$ nebulae with projected sizes larger than 100 physical kpc and, in some cases, extending up to 320 kpc. The circularly averaged surface brightness profiles of the nebulae appear 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$\sim$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<z<4$ independently of the quasar emission opening angle, or extending up to 200 kpc for quasar isotropic emission.

We use a background quasar to detect the presence of circum-galactic gas around a $z=0.91$ low-mass star forming galaxy. Data from the new Multi Unit Spectroscopic Explorer (MUSE) on the VLT show that the host galaxy has a dust-corrected star-formation rate (SFR) of 4.7$\pm$0.2 Msun/yr, with no companion down to 0.22 Msun/yr (5 $\sigma$) within 240 kpc (30"). Using a high-resolution spectrum (UVES) of the background quasar, which is fortuitously aligned with the galaxy major axis (with an azimuth angle $\alpha$ of only $15^\circ$), we find, in the gas kinematics traced by low-ionization lines, distinct signatures consistent with those expected for a "cold flow disk" extending at least 12 kpc ($3\times R_{1/2}$). We estimate the mass accretion rate $\dot M_{\rm in}$ to be at least two to three times larger than the SFR, using the geometric constraints from the IFU data and the HI column density of $\log N_{\rm HI} \simeq 20.4$ obtained from a {\it HST}/COS NUV spectrum. From a detailed analysis of the low-ionization lines (e.g. ZnII, CrII, TiII, MnII, SiII), the accreting material appears to be enriched to about 0.4 $Z_\odot$ (albeit with large uncertainties: $\log Z/Z_\odot=-0.4~\pm~0.4$), which is comparable to the galaxy metallicity ($12+\log \rm O/H=8.7\pm0.2$), implying a large recycling fraction from past outflows. Blue-shifted MgII and FeII absorptions in the galaxy spectrum from the MUSE data reveal the presence of an outflow. The MgII and FeII doublet ratios indicate emission infilling due to scattering processes, but the MUSE data do not show any signs of fluorescent FeII* emission.

We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3.

We present a method to constrain galaxy parameters directly from three-dimensional data cubes. The algorithm compares directly the data with a parametric model mapped in $x,y,\lambda$ coordinates. It uses the spectral lines-spread function (LSF) and the spatial point-spread function (PSF) to generate a three-dimensional kernel whose characteristics are instrument specific or user generated. The algorithm returns the intrinsic modeled properties along with both an `intrinsic' model data cube and the modeled galaxy convolved with the 3D-kernel. The algorithm uses a Markov Chain Monte Carlo (MCMC) approach with a nontraditional proposal distribution in order to efficiently probe the parameter space. We demonstrate the robustness of the algorithm using 1728 mock galaxies and galaxies generated from hydrodynamical simulations in various seeing conditions from 0.6" to 1.2". We find that the algorithm can recover the morphological parameters (inclination, position angle) to within 10% and the kinematic parameters (maximum rotation velocity) to within 20%, irrespectively of the PSF in seeing (up to 1.2") provided that the maximum signal-to-noise ratio (SNR) is greater than $\sim3$ pixel$^{-1}$ and that the ratio of the galaxy half-light radius to seeing radius is greater than about 1.5. One can use such an algorithm to constrain simultaneously the kinematics and morphological parameters of (nonmerging) galaxies observed in nonoptimal seeing conditions. The algorithm can also be used on adaptive-optics (AO) data or on high-quality, high-SNR data to look for nonaxisymmetric structures in the residuals.

We present the results of a survey for intervening HI 21-cm absorbers at intermediate and low redshift (0<z<1.2). For our total sample of 24 systems, we obtained high quality data for 17 systems, the other seven being severely affected by radio frequency interference (RFI). Five of our targets are low redshift (z<0.17) optical galaxies with small impact parameters (<20 kpc) toward radio-bright background sources. Two of these were detected in 21-cm absorption, showing narrow, high optical depth absorption profiles, the narrowest having a velocity dispersion of only 1.5 km/s, which puts an upper limit on the kinetic temperature of T_k<270 K. Combining our observations with results from the literature, we measure a weak anti-correlation between impact parameter and integral optical depth in local (z<0.5) 21-cm absorbers. Of eleven CaII and MgII systems searched, two were detected in 21-cm absorption, and six were affected by RFI to a level that precludes a detection. For these two systems at z~0.6 we measure spin temperatures of T_s=(65+/-17) K and T_s>180 K. A subset of our systems were also searched for OH absorption, but no detections were made.

We present Karl G. Jansky Very Large Array observations of the CO J=1-0 transition in a sample of four $z\sim2$ main sequence galaxies. These galaxies are in the blue sequence of star-forming galaxies at their redshift, and are part of the IRAM Plateau de Bure HIgh-$z$ Blue Sequence Survey (PHIBSS) which imaged them in CO J=3-2. Two galaxies are imaged here at high signal-to-noise, allowing determinations of their disk sizes, line profiles, molecular surface densities, and excitation. Using these and published measurements, we show that the CO and optical disks have similar sizes in main-sequence galaxies, and in the galaxy where we can compare CO J=1-0 and J=3-2 sizes we find these are also very similar. Assuming a Galactic CO-to-H$_2$ conversion, we measure surface densities of $\Sigma_{mol}\sim1200$ M$_\odot$pc$^{-2}$ in projection and estimate $\Sigma_{mol}\sim500-900$ M$_\odot$pc$^{-2}$ deprojected. Finally, our data yields velocity-integrated Rayleigh-Jeans brightness temperature line ratios $r_{31}$ that are approximately unity. In addition to the similar disk sizes, the very similar line profiles in J=1-0 and J=3-2 indicate that both transitions sample the same kinematics, implying that their emission is coextensive. We conclude that in these two main sequence galaxies there is no evidence for significant excitation gradients or a large molecular reservoir that is diffuse or cold and not involved in active star-formation. We suggest that $r_{31}$ in very actively star-forming galaxies is likely an indicator of how well mixed the star formation activity and the molecular reservoir are.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

The physical properties of galactic winds are of paramount importance for our understanding of galaxy formation. Fortunately, they can be constrained using background quasars passing near star-forming galaxies (SFGs). From the 14 quasar$-$galaxy pairs in our VLT/SINFONI Mgii Program for Line Emitters (SIMPLE) sample, we reobserved the 10 brightest galaxies in H$_{\alpha}$ with the VLT/SINFONI with 0.7" seeing and the corresponding quasar with the VLT/UVES spectrograph. Applying geometrical arguments to these ten pairs, we find that four are likely probing galactic outflows, three are likely probing extended gaseous disks, and the remaining three are not classifiable because they are viewed face-on. In this paper we present a detailed comparison between the line-of-sight kinematics and the host galaxy emission kinematics for the pairs suitable for wind studies. We find that the kinematic profile shapes (asymmetries) can be well reproduced by a purely geometrical wind model with a constant wind speed, except for one pair (towards J2357$-$2736) that has the smallest impact parameter b = 6 kpc and requires an accelerated wind flow. Globally, the outflow speeds are $\sim$ 100 km/s and the mass ejection rates (or $\dot M _{\rm out}$) in the gas traced by the low-ionization species are similar to the star formation rate (SFR), meaning that the mass loading factor, $\eta$ = $\dot M _{\rm out}$/SFR, is $\sim$1.0. The outflow speeds are also smaller than the local escape velocity, which implies that the outflows do not escape the galaxy halo and are likely to fall back into the interstellar medium.

We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (tdepl) and gas to stellar mass ratio (Mmolgas/M*) of SFGs near the star formation main-sequence with redshift, specific star formation rate (sSFR) and stellar mass (M*). The CO- and dust-based scaling relations agree remarkably well. This suggests that the CO-H2 mass conversion factor varies little within 0.6dex of the main sequence (sSFR(ms,z,M*)), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that tdepl scales as (1+z)^-0.3 *(sSFR/sSFR(ms,z,M*))^-0.5, with little dependence on M*. The resulting steep redshift dependence of Mmolgas/M* ~(1+z)^3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M* are driven by the flattening of the SFR-M* relation. Throughout the redshift range probed a larger sSFR at constant M* is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the Mmolgas-SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine Mmolgas with an accuracy of 0.1dex in relative terms, and 0.2dex including systematic uncertainties.

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

We present the first observation of a galaxy (z=0.2) that exhibits metal-line absorption back-illuminated by the galaxy ("down-the-barrel") and transversely by a background quasar at a projected distance of 58 kpc. Both absorption systems, traced by MgII, are blueshifted relative to the galaxy systemic velocity. The quasar sight-line, which resides almost directly along the projected minor axis of the galaxy, probes MgI and MgII absorption obtained from Keck/LRIS and Lya, SiII and SiIII absorption obtained from HST/COS. For the first time, we combine two independent models used to quantify the outflow properties for down-the-barrel and transverse absorption. We find that the modeled down-the-barrel deprojected outflow velocities range between $V_{dtb}=45-255$ km/s. The transverse bi-conical outflow model, assuming constant-velocity flows perpendicular to the disk, requires wind velocities $V_{outflow}=40-80$ km/s to reproduce the transverse MgII absorption kinematics, which is consistent with the range of $V_{dtb}$. The galaxy has a metallicity, derived from H$\alpha$ and NII, of $[{\rm O/H}]=-0.21\pm0.08$, whereas the transverse absorption has $[{\rm X/H}]=-1.12\pm0.02$. The galaxy star-formation rate is constrained between $4.6-15$ M$_{\odot}$/yr while the estimated outflow rate ranges between $1.6-4.2$ M$_{\odot}$/yr and yields a wind loading factor ranging between $0.1-0.9$. The galaxy and gas metallicities, the galaxy-quasar sight-line geometry, and the down-the-barrel and transverse modeled outflow velocities collectively suggest that the transverse gas originates from ongoing outflowing material from the galaxy. The $\sim$1 dex decrease in metallicity from the base of the outflow to the outer halo suggests metal dilution of the gas by the time it reached 58 kpc.

Over the past 18 months we have revisited the science requirements for a multi-object spectrograph (MOS) for the European Extremely Large Telescope (E-ELT). These efforts span the full range of E-ELT science and include input from a broad cross-section of astronomers across the ESO partner countries. In this contribution we summarise the key cases relating to studies of high-redshift galaxies, galaxy evolution, and stellar populations, with a more expansive presentation of a new case relating to detection of exoplanets in stellar clusters. A general requirement is the need for two observational modes to best exploit the large (>40 sq. arcmin) patrol field of the E-ELT. The first mode ('high multiplex') requires integrated-light (or coarsely resolved) optical/near-IR spectroscopy of >100 objects simultaneously. The second ('high definition'), enabled by wide-field adaptive optics, requires spatially-resolved, near-IR of >10 objects/sub-fields. Within the context of the conceptual study for an ELT-MOS called MOSAIC, we summarise the top-level requirements from each case and introduce the next steps in the design process.

Galaxies are thought to be fed by the continuous accretion of intergalactic gas, but direct observational evidence has been elusive. The accreted gas is expected to orbit about the galaxy's halo, delivering not just fuel for star-formation but also angular momentum to the galaxy, leading to distinct kinematic signatures. Here we report observations showing these distinct signatures near a typical distant star-forming galaxy where the gas is detected using a background quasar passing 26 kpc from the host. Our observations indicate that gas accretion plays a major role in galaxy growth since the estimated accretion rate is comparable to the star-formation rate.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

We study the properties of luminous stellar clumps identified in deep, high resolution HST/NIC2 F160W imaging at 1.6um of six z~2 star-forming galaxies with existing near-IR integral field spectroscopy from SINFONI at the VLT. Individual clumps contribute ~0.5%-15% of the galaxy-integrated rest-frame ~5000A emission, with median of about 2%; the total contribution of clump light ranges from 10%-25%. The median intrinsic clump size and stellar mass are ~1kpc and log(Mstar[Msun])~9, in the ranges for clumps identified in rest-UV or line emission in other studies. The clump sizes and masses in the subset of disks are broadly consistent with expectations for clump formation via gravitational instabilities in gas-rich, turbulent disks given the host galaxies' global properties. By combining the NIC2 data with ACS/F814W imaging available for one source, and AO-assisted SINFONI Halpha data for another, we infer modest color, M/L, and stellar age variations within each galaxy. In these two objects, sets of clumps identified at different wavelengths do not fully overlap; NIC2-identified clumps tend to be redder/older than ACS- or Halpha-identified clumps without rest-frame optical counterparts. There is evidence for a systematic trend of older ages at smaller galactocentric radii among the clumps, consistent with scenarios where inward migration of clumps transports material towards the central regions. From constraints on a bulge-like component at radii <1-3kpc, none of the five disks in our sample appears to contain a compact massive stellar core, and we do not discern a trend of bulge stellar mass fraction with stellar age of the galaxy. Further observations are necessary to probe the build-up of stellar bulges and the role of clumps in this process.

We present the initial imaging and spectroscopic data acquired as part of the VLT VIMOS Lyman-break galaxy Survey. UBR (or UBVI) imaging covers five 36'x36' fields centred on bright z>3 QSOs, allowing ~21,000 2<z<3.5 galaxy candidates to be selected using the Lyman-break technique. We performed spectroscopic follow-up using VIMOS, measuring redshifts for 1020 z>2 LBGs and 10 z>2 QSOs from a total of 19 VIMOS pointings. From the galaxy spectra, we observe a 625+-510 km/s velocity offset between the ISM absorption and Ly-alpha emission line redshifts. Using the photometric and spectroscopic catalogues, we have analysed the galaxy clustering at z~3. In the photometric case, the angular correlation function, w(theta), is well fit by a double power-law with clustering scale-length, r_0 = 3.19+0.32-0.54 Mpc/h for r < 1 Mpc/h and r_0 = 4.59+0.31-0.33 Mpc/h at larger scales. Using the redshift sample we estimate the semi-projected correlation function, w_p(sigma) and find r_0 = 3.67+0.23-0.24 Mpc/h for the VLT sample and r_0 = 3.98+0.14-0.15 Mpc/h for a combined VLT+Keck sample. From the z-space correlation functions and assuming the above xi(r) models, we find that the combined VLT and Keck surveys require a galaxy velocity dispersion, <w_z^2>^1/2 ~ 700 km/s, higher than the ~400 km/s found by previous authors. We also measure a value for the gravitational growth rate parameter of beta(z=3) = 0.48+-0.17, implying a low value for the bias of b = 2.06+1.1-0.5. This value is consistent with the galaxy clustering amplitude which gives b = 2.22+-0.16, assuming the standard cosmology, implying that the evolution of the gravitational growth rate is also consistent with Einstein gravity. We have compared our LBG clustering amplitudes with lower redshift measurements and find that the clustering strength is not inconsistent with that of low-redshift L* spirals for simple 'long-lived' galaxy models.

We have studied the properties of giant star forming clumps in five z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. The clumps reside in disk regions where the Toomre Q-parameter is below unity, consistent with their being bound and having formed from gravitational instability. Broad H{\alpha}/[NII] line wings demonstrate that the clumps are launching sites of powerful outflows. The inferred outflow rates are comparable to or exceed the star formation rates, in one case by a factor of eight. Typical clumps may lose a fraction of their original gas by feedback in a few hundred million years, allowing them to migrate into the center. The most active clumps may lose much of their mass and disrupt in the disk. The clumps leave a modest imprint on the gas kinematics. Velocity gradients across the clumps are 10-40 km/s/kpc, similar to the galactic rotation gradients. Given beam smearing and clump sizes, these gradients may be consistent with significant rotational support in typical clumps. Extreme clumps may not be rotationally supported; either they are not virialized, or they are predominantly pressure supported. The velocity dispersion is spatially rather constant and increases only weakly with star formation surface density. The large velocity dispersions may be driven by the release of gravitational energy, either at the outer disk/accreting streams interface, and/or by the clump migration within the disk. Spatial variations in the inferred gas phase oxygen abundance are broadly consistent with inside-out growing disks, and/or with inward migration of the clumps.

We present deep and high-resolution HST/NIC2 F160W imaging at 1.6micron of six z~2 star-forming galaxies with existing near-IR integral field spectroscopy from SINFONI at the VLT. The unique combination of rest-frame optical imaging and nebular emission-line maps provides simultaneous insight into morphologies and dynamical properties. The overall rest-frame optical emission of the galaxies is characterized by shallow profiles in general (Sersic index n<1), with median effective radii of ~5kpc. The morphologies are significantly clumpy and irregular, which we quantify through a non-parametric morphological approach, estimating the Gini (G), Multiplicity (Psi), and M_20 coefficients. The strength of the rest-frame optical emission lines in the F160W bandpass indicates that the observed structure is not dominated by the morphology of line-emitting gas, and must reflect the underlying stellar mass distribution of the galaxies. The sizes and structural parameters in the rest-frame optical continuum and Halpha emission reveal no significant differences, suggesting similar global distributions of the on-going star formation and more evolved stellar population. While no strong correlations are observed between stellar population parameters and morphology within the NIC2/SINFONI sample itself, a consideration of the sample in the context of a broader range of z~2 galaxy types indicates that these galaxies probe the high specific star formation rate and low stellar mass surface density part of the massive z~2 galaxy population, with correspondingly large effective radii, low Sersic indices, low G, and high Psi and M_20. The combined NIC2 and SINFONI dataset yields insights of unprecedented detail into the nature of mass accretion at high redshift. [Abridged]

We study the growth of dark matter halos in the concordance LCDM cosmology using several N-body simulations of large cosmological volumes. We build merger trees from the Millennium and Millennium-II simulations, covering a range 10^9-10^15 Msun in halo mass and 1-10^5 in merger mass ratio. Our algorithm takes special care of halo fragmentation and ensures that the mass contribution of each merger to halo growth is only counted once. This way the integrated merger rate converges and we can consistently determine the contribution of mergers of different mass ratios to halo growth. We find that all resolved mergers, up to mass ratios of 10^5, contribute only ~60% of the total halo mass growth, while major mergers are subdominant, e.g. mergers with mass ratios smaller than 3:1 (10:1) contribute only ~20% (~30%). This is verified with an analysis of two additional simulation boxes, where we follow all particles individually throughout cosmic time. Our results are also robust against using several halo definitions. Under the assumption that the power-law behaviour of the merger rate at large mass ratios can be extrapolated to arbitrarily large mass ratios, it is found that, independently of halo mass, ~40% of the mass in halos comes from genuinely smooth accretion of dark matter that was never bound in smaller halos. We discuss possible implications of our findings for galaxy formation. One implication, assuming as is standard that the pristine intergalactic medium is heated and photoionized by UV photons, is that all halos accrete >40% of their baryons in smooth "cold" T>~10^4K gas, rather than as warm, enriched or clumpy gas or as stars.

High signal-to-noise, representative spectra of star-forming galaxies at z~2, obtained via stacking, reveal a high-velocity component underneath the narrow H-alpha and [NII] emission lines. When modeled as a single Gaussian, this broad component has FWHM > 1500 km/s; when modeled as broad wings on the H-alpha and [NII] features, it has FWHM > 500 km/s. This feature is preferentially found in the more massive and more rapidly star-forming systems, which also tend to be older and larger galaxies. We interpret this emission as evidence of either powerful starburst-driven galactic winds or active supermassive black holes. If galactic winds are responsible for the broad emission, the observed luminosity and velocity of this gas imply mass outflow rates comparable to the star formation rate. On the other hand, if the broad line regions of active black holes account for the broad feature, the corresponding black holes masses are estimated to be an order of magnitude lower than those predicted by local scaling relations, suggesting a delayed assembly of supermassive black holes with respect to their host bulges.

We developed a new method to extract halo merger rates from the Millennium Simulation. First, by removing superfluous mergers that are artifacts of the FOF halo finder, we find a lower merger rate compared to previous work. The reductions are more significant at lower redshifts, lower halo masses, and minor mergers. Our approach agrees better with predictions from the EPS model. Second, we find that the FOF halo finder overestimates the halo mass by up to 50% for halos that are about to merge, which leads to an additional ~20% overestimate of the merger rate. Therefore, we define halo masses by including only gravitationally bound particles. We provide new best-fitting parameters for a global formula to account for these improvements. In addition, we extract the merger rate per progenitor halo, as well as per descendant halo. The former is the quantity that is related to observed galaxy merger fractions when they are measured via pair counting. At low mass/redshift the merger rate increases with mass and redshift. At high mass/redshift (for halos with masses a few times the "knee" of the mass function) these trends break down, and the merger rate per progenitor halo decreases with mass and increases only moderately with redshift. Defining the merger rate per progenitor halo also allows us to quantify the rate at which halos are being accreted onto larger halos. We provide an analytic formula to convert given merger rates per descendant halo into merger rates per progenitor halo. Finally, we perform a comparison between observed merger fractions and the fraction of halos that have undergone a major merger during the recent dynamical friction time, and find a fair agreement, within the large uncertainties of the observations. Our new halo merger trees are available at http://www.mpe.mpg.de/ir/MillenniumMergerTrees/.

In these proceedings, we summarize recent results from our "SINS" VLT/SINFONI integral-field survey, focusing on the 52 detected UV/optically-selected star-forming galaxies at z~2. Our H-alpha emission-line imaging and kinematic data of these systems illustrates that a substantial fraction (> 1/3) of these galaxies are large, rotating disks and that these disks are clumpy, thick, and forming stars rapidly. We compare these systems to local disk scaling relations and find that the backbones of these relations are already in place at z~2. Detailed analysis of the large disks in our sample provides strong evidence that this population cannot result from a merger-dominated formation history and instead must be assembled by the smooth but rapid inflow of gas along filaments. These systems will then secularly evolve from clump-dominated disks to bulge-dominated disks on short timescales, a phenomenon that is observed in our SINS observations and is consistent with predictions from numerical simulations. These results provide new and exciting insights into the formation of bulge-dominated galaxies in the local Universe.

Recent observations of UV-/optically selected, massive star forming galaxies at z~2 indicate that the baryonic mass assembly and star formation history is dominated by continuous rapid accretion of gas and internal secular evolution, rather than by major mergers. We use the Millennium Simulation to build new halo merger trees, and extract halo merger fractions and mass accretion rates. We find that even for halos not undergoing major mergers the mass accretion rates are plausibly sufficient to account for the high star formation rates observed in z~2 disks. On the other hand, the fraction of major mergers in the Millennium Simulation is sufficient to account for the number counts of submillimeter galaxies (SMGs), in support of observational evidence that these are major mergers. When following the fate of these two populations in the Millennium Simulation to z=0, we find that subsequent mergers are not frequent enough to convert all z~2 turbulent disks into elliptical galaxies at z=0. Similarly, mergers cannot transform the compact SMGs/red sequence galaxies at z~2 into observed massive cluster ellipticals at z=0. We argue therefore, that secular and internal evolution must play an important role in the evolution of a significant fraction of z~2 UV-/optically and submillimeter selected galaxy populations.