We investigate the ultraviolet (UV) spectral properties of faint Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.9<z<4.6 and provide material to prepare future observations of the faint Universe. We use data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median M$_{UV}$ of -18 and down to M$_{UV}$ of -16), low stellar mass (median value of $10^{8.4}$ and down to $10^{7}M_{\odot}$) LAEs at redshift z>3. 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.

We present strong-lensing analyses of three galaxy clusters, RXJ2129.4+0009 (z=0.235), MS0451.6-0305 (z=0.55), and MACSJ2129.4-0741 (z=0.589), using the powerful combination of Hubble Space Telescope (HST) multi-band observations, and Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy. We confirm the spectroscopic identification of 4 new multiple image systems, and 36 cluster members in RXJ2129. Our resulting mass model is constrained by 8 multiple image systems, and composed of 71 halos including one dark matter cluster-scale halo and 2 galaxy-scale halos optimized individually. For MS0451, we report the spectroscopic identification of 2 new systems of multiple images in the poorly-constrained Northern region, and 44 cluster members. Our mass model uses 16 multiple image systems, and 146 halos, including 2 large-scale halos, and 7 galaxy-scale halos independently optimized. For MACSJ2129, we confirm the spectroscopic identification of 4 new systems of multiple images, and 52 cluster members. We also report the spectroscopic identification of one new multiple image system at $z=4.41$. Our mass model is optimized by 14 multiple image systems, and composed of 151 halos, including 2 large-scale halos and 4 galaxy-scale halos independently optimized. Our best models have rms of 0.33",0.6", 0.73" in the image plane for RXJ2129, MS0451, and MACSJ2129 respectively. This analysis presents a detailed comparison with the existing literature showing excellent agreements, and discuss specific studies of lensed galaxies, e.g. a group of submilimeter galaxies at z=2.9 in MS0451, and a bright z=2.1472 red singly imaged galaxy in MACSJ2129.

Processing of raw data from modern astronomical instruments is nowadays often carried out using dedicated software, so-called "pipelines" which are largely run in automated operation. In this paper we describe the data reduction pipeline of the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph operated at ESO's Paranal observatory. This spectrograph is a complex machine: it records data of 1152 separate spatial elements on detectors in its 24 integral field units. Efficiently handling such data requires sophisticated software, a high degree of automation and parallelization. We describe the algorithms of all processing steps that operate on calibrations and science data in detail, and explain how the raw science data gets transformed into calibrated datacubes. We finally check the quality of selected procedures and output data products, and demonstrate that the pipeline provides datacubes ready for scientific analysis.

We measure the ionizing photon production efficiency ($\xi_{ion}$) of low-mass galaxies ($10^{7.8}$-$10^{9.8}$ $M_{\odot}$) at $1.4<z<2.7$, allowing us to better understand the contribution of dwarf galaxies to the ionizing background and cosmic reionization. We target galaxies that are magnified by the strong lensing galaxy clusters Abell 1689, MACS J0717, and MACS J1149. We utilize Keck/MOSFIRE spectra to measure optical nebular emission line fluxes and HST imaging to measure the rest-UV and rest-optical photometry. We present two methods of stacking. First, we take the average of the log(L$_{H\alpha}$ /L$_{UV}$) of galaxies in our sample to determine the typical log($\xi_{ion}$). Second, we take the logarithm of the total L$_{H\alpha}$ over the total L$_{UV}$. We prefer the latter as it provides the total ionizing UV luminosity density of galaxies when multiplied by the non-ionizing UV luminosity density from the UV luminosity function. log($\xi_{ion}$) calculated from the second method is $\sim$ 0.2 dex higher than the first method. We do not find any strong dependence between log($\xi_{ion}$) and stellar mass, M$_{UV}$ or UV spectral slope ($\beta$). We report a value of log($\xi_{ion}$) $\sim25.47\pm 0.09$ for our UV-complete sample ($-22<M_{UV}<-17.3$) and $\sim25.37\pm0.11$ for our mass-complete sample ($7.8<\log(M_*)<9.8)$. These values are consistent with measurements of more massive, more luminous galaxies in other high-redshift studies that use the same stacking technique. Our log($\xi_{ion}$) is $0.2-0.3$ dex higher than low-redshift galaxies of similar mass, indicating an evolution in the stellar properties, possibly due to metallicity, age, or the prevalence of binary stars. We also find a correlation between log($\xi_{ion}$) and the equivalent widths of H$\alpha$ and [OIII]$\lambda$5007 fluxes, confirming that these equivalent widths can be used to estimate $\xi_{ion}$.

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<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~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.

With the spatial resolution of the Atacama Large Millimetre Array (ALMA), dusty galaxies in the distant Universe typically appear as single, compact blobs of dust emission, with a median half-light radius, $\approx$ 1 kpc. Occasionally, strong gravitational lensing by foreground galaxies or galaxy clusters has probed spatial scales 1-2 orders of magnitude smaller, often revealing late-stage mergers, sometimes with tantalising hints of sub-structure. One lensed galaxy in particular, the Cosmic Eyelash at $z=$ 2.3, has been cited extensively as an example of where the interstellar medium exhibits obvious, pronounced clumps, on a spatial scale of $\approx$ 100 pc. Seven orders of magnitude more luminous than giant molecular clouds in the local Universe, these features are presented as circumstantial evidence that the blue clumps observed in many $z\sim$ 2-3 galaxies are important sites of ongoing star formation, with significant masses of gas and stars. Here, we present data from ALMA which reveal that the dust continuum of the Cosmic Eyelash is in fact smooth and can be reproduced using two S\'ersic profiles with effective radii, 1.2 and 4.4 kpc, with no evidence of significant star-forming clumps down to a spatial scale of $\approx$ 80 pc and a star-formation rate of $<$ 3 M$_\odot$ yr$^{-1}$.

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

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.

We present a new galaxy cluster lens modeling approach, hybrid-Lenstool, that is implemented in the publicly available modeling software Lenstool. hybrid-Lenstool combines a parametric approach to model the core of the cluster, and a non-parametric (free-form) approach to model the outskirts. hybrid-Lenstool optimizes both strong- and weak-lensing constraints simultaneously (Joint-Fit), providing a self-consistent reconstruction of the cluster mass distribution on all scales. In order to demonstrate the capabilities of the new algorithm, we tested it on a simulated cluster. hybrid-Lenstool yields more accurate reconstructed mass distributions than the former Sequential-Fit approach where the parametric and the non-parametric models are optimized successively. Indeed, we show with the simulated cluster that the mass density profile reconstructed with a Sequential-Fit deviates form the input by $2-3\sigma$ at all scales while the Joint-Fit gives a profile that is within $1-1.5\sigma$ of the true value. This gain in accuracy is consequential for recovering mass distributions exploiting cluster lensing and therefore for all applications of clusters as cosmological probes. Finally we found that the Joint-Fit approach yields shallower slope of the inner density profile than the Sequential-Fit approach, thus revealing possible biases in previous lensing studies.

[Abridged] We discovered in the Herschel Reference Survey an extremely bright IR source with $S_{500}$~120mJy (Red Virgo 4 - RV4). Based on IRAM/EMIR and IRAM/NOEMA detections of the CO(5-4), CO(4-3), and [CI] lines, RV4 is located at z=4.724, yielding a total observed L$_{IR}$ of 1.1+/-0.6x0$^{14}$L$_{\odot}$. At the position of the Herschel emission, three blobs are detected with the VLA at 10cm. The CO(5-4) line detection of each blob confirms that they are at the same redshift with the same line width, indicating that they are multiple images of the same source. In Spitzer and deep optical observations, two sources, High-z Lens 1 (HL1) West and HL1 East, are detected at the center of the three VLA/NOEMA blobs. These two sources are placed at z=1.48 with XSHOOTER spectra, suggesting that they could be merging and gravitationally lensing the emission of RV4. HL1 is the second most distant lens known to date in strong lensing systems. The Einstein radius of the lensing system is 2.2"+/-0.2 (20kpc). The high redshift of HL1 and the large Einstein radius are highly unusual for a strong lensing system. We present the ISM properties of the background source RV4. Different estimates of the gas depletion time yield low values suggesting that RV4 is a SB galaxy. Among all high-z SMGs, this source exhibits one of the lowest L$_{[CI]}$ to L$_{IR}$ ratios, 3.2+/-0.9x10$^{-6}$, suggesting an extremely short gas tdepl of only 14+/-5Myr. It also shows a relatively high L$_{[CI]}$ to L$_{CO(4-3)}$ ratio (0.7+/-0.2) and low L$_{CO(5-4)}$ to L$_{IR}$ ratio (only ~50% of the value expected for normal galaxies) hinting a low density of gas. Finally, we discuss that the short tdepl of RV4 can be explained by either a very high SFE, which is difficult to reconcile with major mergers simulations of high-z galaxies, or a rapid decrease of SF, which would bias the estimate of tdepl toward low value.

We present Keck/MOSFIRE and Keck/LRIS spectroscopy of A1689-217, a lensed (magnification $\sim7.9$), star-forming (SFR $\sim$ 16 $\rm M_\odot\ yr^{-1}$), dwarf (log($M_\ast/M_\odot$) = $8.07-8.59$), Ly$\alpha$-emitter ($\rm EW_0\sim138$ A) at $z$ = 2.5918. Dwarf galaxies similar to A1689-217 are common at high redshift and likely responsible for reionization, yet few have been studied with detailed spectroscopy. We report a 4.2$\sigma$ detection of the electron-temperature-sensitive [O III]$\lambda$4363 emission line and use this line to directly measure an oxygen abundance of 12+$\log$(O/H) = 8.06 $\pm$ 0.12 ($\sim1/4\ \rm Z_\odot$). A1689-217 is the lowest mass galaxy at $z>$ 2 with an [O III]$\lambda$4363 detection. Using the rest-optical emission lines, we measure A1689-217's other nebular conditions including electron temperature ($T_e$([O III]) $\sim$ 14,000 K), electron density ($n_e\sim$ 220 $\rm cm^{-3}$) and reddening ($E(B-V)\sim$ 0.39). We study relations between strong-line ratios and direct metallicities with A1689-217 and other galaxies with [O III]$\lambda$4363 detections at $z\sim0-3.1$, showing that the locally-calibrated, oxygen-based, strong-line relations are consistent from $z\sim0-3.1$. We also show additional evidence that the $\rm O_{32}$ vs. $\rm R_{23}$ excitation diagram can be utilized as a redshift-invariant, direct-metallicity-based, oxygen abundance diagnostic out to $z\sim$ 3.1. From this excitation diagram and the strong-line ratio $-$ metallicity plots, we observe that the ionization parameter at fixed O/H is consistent with no redshift evolution. Although A1689-217 is metal-rich for its $M_\ast$ and SFR, we find it to be consistent within the large scatter of the low-mass end of the Fundamental Metallicity Relation.

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.

The cold molecular gas in contemporary galaxies is structured in discrete cloud complexes. These giant molecular clouds (GMCs), with $10^4$-$10^7$ solar masses and radii of 5-100 parsecs, are the seeds of star formation. Highlighting the molecular gas structure at such small scales in distant galaxies is observationally challenging. Only a handful of molecular clouds were reported in two extreme submillimetre galaxies at high redshift. Here we search for GMCs in a typical Milky Way progenitor at z = 1.036. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we mapped the CO(4-3) emission of this gravitationally lensed galaxy at high resolution, reading down to 30 parsecs, which is comparable to the resolution of CO observations of nearby galaxies. We identify 17 molecular clouds, characterized by masses, surface densities and supersonic turbulence all of which are 10-100 times higher than present-day analogues. These properties question the universality of GMCs and suggest that GMCs inherit their properties from ambient interstellar medium. The measured cloud gas masses are similar to the masses of stellar clumps seen in the galaxy in comparable numbers. This corroborates the formation of molecular clouds by fragmentation of distant turbulent galactic gas disks, which then turn into stellar clumps ubiquitously observed in galaxies at cosmic noon.

We compare the molecular and ionized gas kinematics of two strongly lensed galaxies at $z\sim1$ based on observations from ALMA and MUSE. We derive the CO and [OII] rotation curves and dispersion profiles of these two galaxies. We find a difference between the observed molecular and ionized gas rotation curves for one of the galaxies, the Cosmic Snake, for which we obtain a spatial resolution of few hundred parsecs along the major axis. The rotation curve of the molecular gas is steeper than the rotation curve of the ionized gas. In the second galaxy, A521, the molecular and ionized gas rotation curves are consistent, but the spatial resolution is only of few kpc on the major axis. Using simulations, we investigate the effect of the thickness of the gas disk and effective radius on the observed rotation curves and find that a more extended and thicker disk smooths the curve. We also find that the presence of a strongly inclined thick disk (>1 kpc) can smooth the rotation curve because it degrades the spatial resolution along the line of sight. By building a model using a stellar disk and two gas disks, we reproduce the rotation curves of the Cosmic Snake with a molecular gas disk that is more massive and more radially and vertically concentrated than the ionized gas disk. Finally, we also obtain an intrinsic velocity dispersion in the Cosmic Snake of 18.5+-7 km/s and 19.5+-6 km/s for the molecular and ionized gas, respectively, which is consistent with a molecular disk with a smaller and thinner disk. For A521, the intrinsic velocity dispersion values are 11+-8 km/s and 54+-11 km/s, with a higher value for the ionized gas. This could indicate that the ionized gas disk is thicker and more turbulent in this galaxy. These results highlight the different spatial distribution of the molecular and ionized gas disks at $z\sim1$ and suggest the presence of thick ionized gas disks at this epoch.

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

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

Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on $H_0$ of a few-percent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the lens model. In this work, we present a detailed characterization of the environment of the lens system WFI2033-4723 ($z_{\rm src} = 1.662$, $z_{\rm lens}$ = 0.6575), one of the core targets of the H0LICOW project for which we present cosmological inferences in a companion paper (Rusu et al. 2019). We use the Gemini and ESO-Very Large telescopes to measure the spectroscopic redshifts of the brightest galaxies towards the lens, and use the ESO-MUSE integral field spectrograph to measure the velocity-dispersion of the lens ($\sigma_{\rm {los}}= 250^{+15}_{-21}$ km/s) and of several nearby galaxies. In addition, we measure photometric redshifts and stellar masses of all galaxies down to $i < 23$ mag, mainly based on Dark Energy Survey imaging (DR1). Our new catalog, complemented with literature data, more than doubles the number of known galaxy spectroscopic redshifts in the direct vicinity of the lens, expanding to 116 (64) the number of spectroscopic redshifts for galaxies separated by less than 3 arcmin (2 arcmin) from the lens. Using the flexion-shift as a measure of the amplitude of the gravitational perturbation, we identify 2 galaxy groups and 3 galaxies that require specific attention in the lens models. The ESO MUSE data enable us to measure the velocity-dispersions of three of these galaxies. These results are essential for the cosmological inference analysis presented in Rusu et al. (2019).

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.

We present the first measurements of the Lyman-continuum photon production efficiency $\xi_{\textrm{ion,0}}$ at $z\sim4$-5 for galaxies fainter than 0.2 $L^*$ ($-$19 mag). $\xi_{\textrm{ion,0}}$ quantifies the production rate of ionizing photons with respect to the UV luminosity density assuming a fiducial escape fraction of zero. Extending previous measurements of $\xi_{\textrm{ion,0}}$ to the faint population is important, as ultra-faint galaxies are expected to contribute the bulk of the ionizing emissivity. We probe $\xi_{\textrm{ion,0}}$ to such faint magnitudes by taking advantage of 200-hour depth Spitzer/IRAC observations from the GREATS program and $\approx$300 3<$z$<6 galaxies with spectroscopic redshifts from the MUSE GTO Deep + Wide programs. Stacked IRAC [3.6]$-$[4.5] colors are derived and used to infer the H$\alpha$ rest-frame equivalent widths, which range from 403\r{A} to 2818\r{A}. The derived $\xi_{\textrm{ion,0}}$ is $\log_{10}(\xi_{\textrm{ion,0}} / \textrm{Hz erg}^{-1}) = 25.36 \pm 0.08$ over $-$20.5 < M$_{\textrm{UV}}$ < $-$17.5, similar to those derived for brighter galaxy samples at the same redshift and therefore suggesting that $\xi_{\textrm{ion}}$ shows no strong dependence on $M_{UV}$. The $\xi_{\textrm{ion,0}}$ values found in our sample imply that the Lyman-continuum escape fraction for $M_{\textrm{UV}} \approx -19$ star-forming galaxies cannot exceed $\approx$8-20\% in the reionization era.

We present the results obtained with VLT/MUSE on the faint-end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The precise aim of the present study is to further constrain the abundance of Lyman-alpha emitters (LAEs) by taking advantage of the magnification provided by lensing clusters. We blindly selected a sample of 156 LAEs, with redshifts between $2.9 \le z \le 6.7$ and magnification-corrected luminosities in the range $ 39 \lesssim \log L_{Ly_{\alpha}}$ [erg s$^{-1}$] $\lesssim 43$. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure. To properly take into account the individual differences in detection conditions (including lensing configurations, spatial and spectral morphologies) when computing the LF, a new method based on the 1/Vmax approach was implemented. The LAE LF has been obtained in four different redshift bins with constraints down to $\log L_{Ly_{\alpha}} = 40.5$. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis including data from the literature to complete the present sample a steep faint-end slope was measured varying from $\alpha = -1.69^{+0.08}_{-0.08}$ to $\alpha = -1.87^{+0.12}_{-0.12}$ between the lowest and the highest redshift bins. The contribution of the LAE population to the star formation rate density at $z \sim 6$ is $\lesssim 50$% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. (abridged)

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

We recommend a conceptual design study for a spectroscopic facility in the southern hemisphere comprising a large diameter telescope, fiber system, and spectrographs collectively optimized for massively-multiplexed spectroscopy. As a baseline, we propose an 11.4-meter aperture, optical spectroscopic survey telescope with a five square degree field of view. Using current technologies, the facility could be equipped with 15,000 robotically-controlled fibers feeding spectrographs over 360<lambda<1330 nm with options for fiber-fed spectrographs at high resolution and a panoramic IFU at a separate focus. This would enable transformational progress via its ability to access a larger fraction of objects from Gaia, LSST, Euclid, and WFIRST than any currently funded or planned spectroscopic facility. An ESO-sponsored study (arXiv:1701.01976) discussed the scientific potential in ambitious new spectroscopic surveys in Galactic astronomy, extragalactic astronomy, and cosmology. The US community should establish links with European and other international communities to plan for such a powerful facility and maximize the potential of large aperture multi-object spectroscopy given the considerable investment in deep imaging surveys.

Using the Spitzer Space Telescope and Herschel Space Observatory, we have conducted a survey of infrared galaxies in the field of the galaxy cluster Abell S1063 (AS1063) at $z=0.347$, which is one of the most massive clusters known and a target of the HST CLASH and Frontier-Field surveys. The Spitzer/MIPS 24 $\mu$m and Herschel/PACS & SPIRE images revealed that the core of AS1063 is surprisingly devoid of infrared sources, showing only a few detectable sources within the central r$\sim1^{\prime}$. There is, however, one particularly bright source (2.3 mJy at 24 $\mu$m; 106 mJy at 160 $\mu$m), which corresponds to a background galaxy at $z=0.61$. The modest magnification factor (4.0$\times$) implies that this galaxy is intrinsically IR-luminous (L$_{\rm IR}=3.1\times10^{11}\ \rm L_{\odot}$). What is particularly interesting about this galaxy is that HST optical/near-infrared images show a remarkably bright and large (1 kpc) clump at one edge of the disk. Our follow-up optical/near-infrared spectroscopy shows Balmer (H$\alpha$-H8) and forbidden emission from this clump ([OII] $\lambda$3727, [OIII] $\lambda\lambda$4959,5007, [NII] $\lambda\lambda$6548,6583), indicating that it is a HII region. The HII region appears to have formed in-situ, as kinematically it is part of a rotating disk, and there is no evidence of nearby interacting galaxies. With an extinction correction of A$_{\rm V}=1.5$ mag, the star formation rate of this giant HII region is $\sim$10 M$_{\odot}$ yr$^{-1}$, which is exceptionally large, even for high redshift HII regions. Such a large and luminous HII region is often seen at $z\sim2$ but quite rare in the nearby Universe.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We present an updated strong-lensing analysis of the massive cluster Abell 370 (A370), continuing the work first presented in Lagattuta et al. (2017). In this new analysis, we take advantage of the deeper imaging data from the Hubble Space Telescope (HST) Frontier Fields program, as well as a large spectroscopic mosaic obtained with the Multi-Unit Spectroscopic Explorer (MUSE). Thanks to the extended coverage of this mosaic, we probe the full 3D distribution of galaxies in the field, giving us a unique picture of the extended structure of the cluster and its surroundings. Our final catalog contains 584 redshifts, representing the largest spectroscopic catalog of A370 to date. Constructing the model, we measure a total mass distribution that is quantitatively similar to our previous work -- though to ensure a low rms error in the model fit, we invoke a significantly large external shear term. Using the redshift catalog, we search for other bound groups of galaxies, which may give rise to a more physical interpretation of this shear. We identify three structures in narrow redshift ranges along the line of sight, highlighting possible infalling substructures into the main cluster halo. We also discover additional substructure candidates in low-resolution imaging at larger projected radii. More spectroscopic coverage of these regions (pushing close to the A370 virial radius) and more extended, high-resolution imaging will be required to investigate this possibility, further advancing the analysis of these interesting developments.

Groups and clusters of galaxies are a current focus of astronomical research owing to their role in determining the environmental effects on galaxies and the constraints they provide to cosmology. The eROSITA X-ray telescope on board the Spectrum Roentgen Gamma observatory will be launched in 2019 and will have completed eight scans of the full sky when 4MOST starts operating. The experiment will detect groups and clusters of galaxies through X-ray emission from the hot intergalactic medium. The purpose of the 4MOST eROSITA Galaxy Cluster Redshift Survey is to provide spectroscopic redshifts of the optical counterparts to the X-ray emission from 40,000 groups and clusters of galaxies so as to perform dynamical estimates of the total mass and to measure the properties of the member galaxies. The survey aims to obtain precise redshift measurements of the photometrically identified brightest cluster galaxies at redshift $z > 0.7$. At lower redshifts ($z < 0.7$), the programme aims to sample over 15 member galaxies per cluster and enable dynamical mass measurements to calibrate the clusters for cosmological experiments. At $z < 0.2$, eROSITA will also detect X-ray emission from galaxy groups and filaments. 4MOST spectroscopic data from the survey will be used for optical identification of galaxy groups down to eROSITA's mass detection limits of $10^{13} M_\odot$, as well as the detection of the largest filaments for pioneering studies of their X-ray emission.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission.

We present observations of two strong-lensing galaxy clusters located within the $90$ per cent credible sky localization maps released following LIGO-Virgo's discovery of the binary black hole (BH-BH) gravitational wave (GW) source GW170814. Our objectives were (1) to search for candidate electromagnetic (EM) counterparts to GW170814 under the hypothesis that it was strongly-lensed, and thus more distant and less massive than inferred by LIGO-Virgo, and (2) to demonstrate the feasibility of rapid target of opportunity observations to search for faint lensed transient point sources in crowded cluster cores located within GW sky localizations. Commencing $20$ hours after discovery, and continuing over $12$ nights, we observed Abell 3084 ($z=0.22$) and SMACSJ0304.3$-$4401 ($z=0.46$) with GMOS on the Gemini-South telescope, and Abell 3084 with MUSE on ESO's Very Large Telescope. We detect no candidate EM counterparts in these data. Calibration of our photometric analysis methods using simulations yield $5\sigma$ detection limits for transients in difference images of the cores of these clusters of $i=25$. This is the most sensitive photometric search to date for counterparts to GW sources, and rules out the possibility that GW170814 was lensed by these clusters with a kilonova-like EM counterpart. Based on the detector frame masses of the compact objects, and assuming that at least one Neutron Star (NS) is required in the merging system to produce a kilonova-like counterpart, implies that GW170814 was neither a NS-NS nor NS-BH merger at $z>8$ lensed by either of these clusters. Also, in the first ever emission line search for counterparts to GW sources, we detected no lines down to a $5\sigma$ detection limit of $5\times10^{-17}{\rm erg\,s^{-1}\,cm^{-2}}$.

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

Deep optical spectroscopic surveys of galaxies provide us a unique opportunity to investigate rest-frame ultra-violet (UV) emission line properties of galaxies at ${z \sim 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 HeII${\lambda}$1640 emitters at ${z\sim2}$. The deepest areas of our MUSE pointings reach a ${3\sigma}$ line flux limit of 3.1${\times}$10-19 erg/ s/ cm$^2$. After discarding broad line active galactic nuclei we find 13 HeII${\lambda}$1640 detections from MUSE with a median MUV = $-20.1$ and 21 tentative HeII${\lambda}$1640 detections from other public surveys. Excluding Ly${\alpha}$, all except two galaxies in our sample show at least one other rest-UV emission line, with CIII]${\lambda}$1907,${\lambda}$1909 being the most prominent. We use multi-wavelength data available in the Hubble legacy fields to derive basic galaxy properties of our sample via spectral energy distribution fitting techniques. Taking advantage of the high quality spectra obtained by MUSE (${\sim10 - 30}$h of exposure time per pointing), we use photo-ionisation models to study the rest-UV emission line diagnostics of the HeII${\lambda}$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 HeII${\lambda}$1640 equivalent widths (${\sim}$ 0.2 - 10 A), thus additional mechanisms are necessary in models to match the observed HeII${\lambda}$1640 properties.

We analyse maps of the spatially-resolved nebular emission of $\approx$1500 star-forming galaxies at $z\approx0.6$-$2.2$ from deep KMOS and MUSE observations to measure the average shape of their rotation curves. We use these to test claims for declining rotation curves at large radii in galaxies at $z\approx1$-$2$ that have been interpreted as evidence for an absence of dark matter. We show that the shape of the average rotation curves, and the extent to which they decline beyond their peak velocities, depends upon the normalisation prescription used to construct the average curve. Normalising in size by the galaxy stellar disk-scale length after accounting for seeing effects ($R_{\rm{d}}^{\prime}$), we construct stacked position-velocity diagrams that trace the average galaxy rotation curve out to $6R_{\rm{d}}^{\prime}$ ($\approx$13 kpc, on average). Combining these curves with average HI rotation curves for local systems, we investigate how the shapes of galaxy rotation curves evolve over $\approx$10 Gyr. The average rotation curve for galaxies binned in stellar mass, stellar surface mass density and/or redshift is approximately flat, or continues to rise, out to at least $6R_{\rm{d}}^{\prime}$. We find a trend between the outer slopes of galaxies' rotation curves and their stellar mass surface densities, with the higher surface density systems exhibiting flatter rotation curves. Drawing comparisons with hydrodynamical simulations, we show that the average shapes of the rotation curves for our sample of massive, star-forming galaxies at $z\approx0$-$2.2$ are consistent with those expected from $\Lambda$CDM theory and imply dark matter fractions within $6R_{\rm{d}}$ of at least $\approx60$ percent.

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 <100 pc scales, we find two correlations previously reported at kpc scales: a spatially resolved Mass-Metallicity Relation (RMZR) and a spatially resolved SFMS (RSFMS). We find no secondary dependency of the RMZR with the SFR density. We find that both resolved relations have a local origin, as they do not depend on the total stellar mass. The observational results of this paper are consistent with the inside-out scenario for the growth of galactic disks.

With a projected size of about 450 kpc at z~2.3, the Slug Ly-alpha nebula is a rare laboratory to study, in emission, the properties of the intergalactic gas in the Cosmic Web. Since its discovery, the Slug has been the subject of several spectroscopic follow-ups to constrain the properties of the emitting gas. Here we report the results of a deep MUSE integral-field spectroscopic search for non-resonant, extended HeII1640 and metal emission. Extended HeII radiation is detected on scales of about 100 kpc, but only in some regions associated with the bright Ly-alpha emission and a continuum-detected source, implying large and abrupt variations in the line ratios across adjacent regions in projected space. The recent detection of associated H-alpha emission and similar abrupt variations in the Ly-alpha kinematics, strongly suggest that the HeII/Ly-alpha gradient is due to large variations in the physical distances between the associated quasar and these regions. This implies that the overall length of the emitting structure could extend to physical Mpc scales and be mostly oriented along our line of sight. At the same time, the relatively low HeII/Ly-alpha values suggest that the emitting gas has a broad density distribution that - if expressed in terms of a lognormal - implies dispersions as high as those expected in the interstellar medium of galaxies. These results strengthen the possibility that the density distribution of intergalactic gas at high-redshift is extremely clumpy and multiphase on scales below our current observational spatial resolution of a few physical kpc.

We present a multi-wavelength analysis of the core of the massive galaxy cluster MACS\,J0417.5-1154 ($z = 0.441$; MACS\;J0417). Our analysis takes advantage of VLT/MUSE observations which allow the spectroscopic confirmation of three strongly-lensed systems. One of these, nick-named \emph{The Doughnut}, consists of three complete images of a complex ring galaxy at $z = 0.8718$ and a fourth, partial and radial image close to the Brightest Cluster Galaxy (BCG) only discernible thanks to its strong [OII] line emission. The best-fit mass model (rms of 0.38\arcsec) yields a two-dimensional enclosed mass of $M({\rm R < 200\,kpc}) = (1.77\pm0.03)\times10^{14}\,\msun$ and almost perfect alignment between the peaks of the BCG light and the dark matter of ($0.5\pm0.5$)\arcsec . Our finding that a significant misalignment results when the radial image of \emph{The Doughnut} is omitted serves as an important caveat for studies of BCG-dark matter offsets in galaxy clusters. Using \emph{Chandra} data to map the intra-cluster gas, we observe an offset between the gas and dark-matter peaks of ($1.7\pm0.5$)\arcsec, and excellent alignment of the X-ray peak with the location of optical emission line associated with the BCG. We interpret all observational evidence in the framework of on-going merger activity, noting specifically that the coincidence between the gas peak and the peak of blue light from the BCG may be evidence of dense, cold gas leading to direct star formation. We use the surface area $\sigma_{\mu}$ above a given magnification factor $\mu$ as a metric to estimate the lensing power of MACS\,J0417. We obtain $\sigma(\mu > 3) = 0.22$\,arcmin$^2$ confirming MACS\,J0417 as an efficient gravitational lens. Finally, we discuss the differences between our mass model and Mahler et al. (2018).

(Abridged) We investigate the Lyman $\alpha$ emitter luminosity function (LAE LF) within the redshift range $2.9 \leq z \leq 6$ from the first instalment of the blind integral field spectroscopic survey MUSE-Wide. This initial part of the survey probes a region of 22.2 arcmin$^2$ in the CANDELS/GOODS-S field. The dataset provided us with 237 LAEs from which we construct the LAE LF in the luminosity range $42.2 \leq \log L_\mathrm{Ly\alpha} [\mathrm{erg\,s}^{-1}]\leq 43.5$ within a volume of $2.3\times10^5$ Mpc$^3$. For the LF construction we utilise three different non-parametric estimators: The classical $1/V_\mathrm{max}$ 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 $\Phi(\log L_\mathrm{Ly\alpha} [\mathrm{erg\,s}^{-1}] = 43.5) \simeq 3\times 10^{-6}$ Mpc$^{-3}$ and $\Phi(\log L_\mathrm{Ly\alpha} [\mathrm{erg\,s}^{-1}] = 42.2) \simeq 2 \times 10^{-3}$ Mpc$^{-3}$ towards the bright- and faint-end of our survey, respectively. By employing a non-parametric statistical test, as well as by comparing the full sample to sub-samples in redshift bins, we find no supporting evidence for an evolving LAE LF over the probed redshift and luminosity range. We determine the best-fitting Schechter function parameters $\alpha = -1.84^{+0.42}_{-0.41}$ and $\log L^* [\mathrm{erg\,s}^{-1}] = 42.2^{+0.22}_{-0.16}$ with the corresponding normalisation $\log \phi^* [\mathrm{Mpc}^{-3}] = -2.71$. When correcting for completeness in the LAE LF determinations, we take into account that LAEs exhibit diffuse extended low surface-brightness haloes. We compare the resulting LF to one obtained where we apply 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.

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.

The aim of this paper is to introduce the WIRCam Ultra Deep Survey (WUDS), a near-IR photometric public survey carried out at the CFH Telescope in the field of the CFHTLS-D3 field (Groth Strip). WUDS includes four near-IR bands (Y, J, H and K_s) over a field of view of ~400 arcmin^2. The typical depth of WUDS data reaches between ~26.8 in Y and J, and ~26 in H and K_s (AB, 3 sigma in 1.3 arcsec aperture). The area and depth of this survey were specifically tailored to set strong constraints on the cosmic star formation rate and the luminosity function brighter or around L* in the z~6-10 redshift domain, although these data are also useful for a variety of extragalactic projects.This first paper is intended to present the properties of WUDS: catalog building, completeness and depth, number counts, photometric redshifts, and global properties of the galaxy population. We have also concentrated on the study of galaxy samples at z~[4.5-7] in this field. UV luminosity functions were derived at z~5 and z~6 taking advantage from the fact that WUDS covers a particularly interesting regime at intermediate luminosities, which allows a combined determination of M* and Phi* with increased accuracy. Our results on the luminosity function are consistent with a small evolution of both M* and Phi* between z=5 and z=6, irrespective of the method used to derive them, either photometric redshifts applied to blindly-selected dropout samples or the classical Lyman Break Galaxy color-preselected samples. Our results lend support to higher Phi* determinations at z=6 than usually reported. The selection and combined analysis of different galaxy samples at z>7 will be presented in a forthcoming paper. WUDS is intended to provide a robust database in the near-IR for the selection of targets for detailed spectroscopic studies, in particular for the EMIR/GTC GOYA Survey (Abridged)

Using an ultra-deep blind survey with the MUSE integral field spectrograph on the ESO Very Large Telescope, we obtain spectroscopic redshifts to a depth never explored before: galaxies with observed magnitudes m > 30 - 32. Specifically, we detect objects via Lyman-alpha emission at 2.9 < z < 6.7 without individual continuum counterparts in areas covered by the deepest optical/near-infrared imaging taken by the Hubble Space Telescope, the Hubble Ultra Deep Field. In total, we find 102 such objects in 9 square arcminutes at these redshifts. Detailed stacking analyses confirm the Lyman-alpha emission as well as the 1216 Angstrom-breaks and faint UV continua (M_UV ~ -15). This makes them the faintest spectroscopically-confirmed objects at these redshifts, similar to the sources believed to reionize the universe. A simple model for the expected fraction of detected/undetected Lyman-alpha emitters as a function of luminosity is consistent with these objects being the high-equivalent width tail of the normal Lyman-alpha-emitter population at these redshifts.

Observations have shown that massive star-forming clumps are present in the internal structure of high-redshift galaxies. One way to study these clumps in detail with a higher spatial resolution is by exploiting the power of strong gravitational lensing which stretches images on the sky. In this work, we present an analysis of the clumpy galaxy A68-HLS115 at $z=1.5858$, located behind the cluster Abell 68, but strongly lensed by a cluster galaxy member. Resolved observations with SINFONI/VLT in the near-infrared show Ha, Hb, [NII], and [OIII] emission lines. Combined with images covering the B band to the far-infrared and CO(2-1) observations, this makes this galaxy one of the only sources for which such multi-band observations are available and for which it is possible to study the properties of resolved star-forming clumps and to perform a detailed analysis of the integrated properties, kinematics, and metallicity. We obtain a stability of $\upsilon_{rot}/\sigma_0 = 2.73$ by modeling the kinematics, which means that the galaxy is dominated by rotation, but this ratio also indicates that the disk is marginally stable. We find a high intrinsic velocity dispersion of $80\pm10$ km s$^{-1}$ that could be explained by the high gas fraction of $f_{gas}=0.75\pm0.15$ observed in this galaxy. This high $f_{gas}$ and the observed sSFR of $\rm 3.12 \, Gyr^{-1}$ suggest that the disk turbulence and instabilities are mostly regulated by incoming gas. The direct measure of the Toomre stability criterion of $Q_{crit}=0.70$ could also indicate the presence of a quasi-stable thick disk. Finally, we identify three clumps in the Ha map which have similar velocity dispersions, metallicities, and seem to be embedded in the rotating disk. These three clumps contribute together to $40\%$ on the SFR(Ha) of the galaxy and show a SFR density about $100$ times higher than HII regions in the local Universe.

We present spatially resolved analysis of a lensed galaxy, SDSS1958+5950 at $z = 2.225$, from the Cambridge Sloan Survey of Wide Arcs in the Sky (CASSOWARY). We use our new high resolution imaging data to construct a robust lens model for the galaxy group at $z = 0.214$. We employ the updated lens model to combine the Integral Field Spectrographic observations on two highly distorted images of the lensed target. We adopt a forward-modeling approach to deconvolve the effects of point spread function from the combined source-plane reconstruction. The approach is adapted to the lens model magnification and enables a resolution of $\sim$170 pc in the galaxy-source plane. We propose an ongoing merger as the origin of the lensed system on the basis of its source-plane morphology, kinematics and rest-frame emission line ratios. Using our novel technique of adaptive coadded source plane reconstruction, we are able to detect different components in the velocity gradient that were not seen in previous studies of this object, plausibly belonging to different components in the merging system.

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.

We present a study of the [OIII]5007/[OII]3727 (O32) ratios of star-forming galaxies drawn from MUSE data spanning a redshift range 0.28<z<0.85. Recently discovered Lyman continuum (LyC) emitters have extremely high oxygen line ratios: O32>4. Here we aim to understand the properties and the occurrences of galaxies with such high line ratios. Combining data from several MUSE GTO programmes, we select a population of star-forming galaxies with bright emission lines, from which we draw 406 galaxies for our analysis based on their position in the z-dependent star formation rate (SFR) - stellar mass (M*) plane. Out of this sample 15 are identified as extreme oxygen emitters based on their O32 ratios (3.7%) and 104 galaxies have O32>1 (26%). Our analysis shows no significant correlation between M*, SFR, and the distance from the SFR-M* relation with O32. We find a decrease in the fraction of galaxies with O32>1 with increasing M*, however, this is most likely a result of the relationship between O32 and metallicity, rather than between O32 and M*. We draw a comparison sample of local analogues with <z>~0.03 from SDSS, and find similar incidence rates for this sample. In order to investigate the evolution in the fraction of high O32 emitters with redshift, we bin the sample into three redshift subsamples of equal number, but find no evidence for a dependence on redshift. Furthermore, we compare the observed line ratios with those predicted by nebular models with no LyC escape and find that most of the extreme oxygen emitters can be reproduced by low metallicity models. The remaining galaxies are likely LyC emitter candidates. Finally, based on a comparison between electron temperature estimates from the [OIII4363]/[OIII]5007 ratio of the extreme oxygen emitters and nebular models, we argue that the galaxies with the most extreme O32 ratios have young light-weighted ages.

We present adaptive optics (AO) assisted integral-field spectroscopy of the intermediate-age star cluster NGC 419 in the Small Magellanic Cloud. By investigating the cluster dynamics and the rotation properties of main sequence turn-off stars (MSTO), we demonstrate the power of AO-fed MUSE observations for this class of objects. Based on 1 049 radial velocity measurements, we determine a dynamical cluster mass of 1.4+/-0.2x10^5 M_sun and a dynamical mass-to-light ratio of 0.67+/-0.08, marginally higher than simple stellar population predictions for a Kroupa initial mass function. A stacking analysis of spectra at both sides of the extended MSTO reveals significant rotational broadening. Our analysis further provides tentative evidence that red MSTO stars rotate faster than their blue counterparts. We find average V sin i values of 87+/-16 km/s and 130+/-22 km/s for blue and red MSTO stars, respectively. Potential systematic effects due to the low spectral resolution of MUSE can reach 30 km/s but the difference in V sin i between the populations is unlikely to be affected.

The Maunakea Spectroscopic Explorer (MSE) Project is a planned replacement for the existing 3.6-m Canada France Hawaii Telescope (CFHT) into a 10-m class dedicated wide field highly multiplexed fibre fed spectroscopic facility. MSE seeks to tackle basic science questions ranging from the origin of stars and stellar systems, Galaxy archaeology at early times, galaxy evolution across cosmic time, to cosmology and the nature of dark matter and dark energy. MSE will be a primary follow-up facility for many key future photometric and astrometric surveys, as well as a major component in the study of the multi-wavelength Universe. The MSE is based on a prime focus telescope concept which illuminate 3200 fibres or more. These fibres are feeding a Low Moderate Resolution (LMR) spectrograph and a High Resolution (HR). The LMR will provide 2 resolution modes at R>2500 and R>5000 on a wavelength range of 360 to 950 nm and a resolution of R>3000 on the 950 nm to 1300 nm bandwidth. Possibly the H band will be also covered by a second NIR mode from ranging from 1450 to 1780 nm. The HR will have a resolution of R>39000 on the 360 to 600 nm wavelength range and R>20000 on the 600 to 900 nm bandwidth. This paper presents the LMR design after its Conceptual Design Review held in June 2017. It focuses on the general concept, optical and mechanical design of the instrument. It describes the associated preliminary expected performances especially concerning optical and thermal performances.

Galaxies at low-redshift typically possess negative gas-phase metallicity gradients (centres more metal-rich than their outskirts). Whereas, it is not uncommon to observe positive metallicity gradients in higher-redshift galaxies ($z \gtrsim 0.6$). Bridging these epochs, we present gas-phase metallicity gradients of 84 star-forming galaxies between $0.08 < z < 0.84$. Using the galaxies with reliably determined metallicity gradients, we measure the median metallicity gradient to be negative ($-0.039^{+0.007}_{-0.009}$ dex/kpc). Underlying this, however, is significant scatter: $(8\pm3)\%\ [7]$ of galaxies have significantly positive metallicity gradients, $(38 \pm 5)\%\ [32]$ have significantly negative gradients, $(31\pm5)\%\ [26]$ have gradients consistent with being flat. (The remaining $(23\pm5)\%\ [19]$ have unreliable gradient estimates.) We notice a slight trend for a more negative metallicity gradient with both increasing stellar mass and increasing star formation rate (SFR). However, given the potential redshift and size selection effects, we do not consider these trends to be significant. Indeed, once we normalize the SFR relative to that of the main sequence, we do not observe any trend between the metallicity gradient and the normalized SFR. This is contrary to recent studies of galaxies at similar and higher redshifts. We do, however, identify a novel trend between the metallicity gradient of a galaxy and its size. Small galaxies ($r_d < 3$ kpc) present a large spread in observed metallicity gradients (both negative and positive gradients). In contrast, we find no large galaxies ($r_d > 3$ kpc) with positive metallicity gradients, and overall there is less scatter in the metallicity gradient amongst the large galaxies. These large (well-evolved) galaxies may be analogues of present-day galaxies, which also show a common negative metallicity gradient.

Recent studies have revealed intense UV metal emission lines in a modest sample of z>7 Lyman-alpha emitters, indicating a hard ionizing spectrum is present. If such high ionization features are shown to be common, it may indicate that extreme radiation fields play a role in regulating the visibility of Lyman-alpha in the reionization era. Here we present deep near-infrared spectra of seven galaxies with Lyman-alpha emission at 5.4<z<8.7 (including a newly-confirmed lensed galaxy at z=6.031) and three bright z>7 photometric targets. In nine sources we do not detect UV metal lines. However in the z=8.683 galaxy EGSY8p7, we detect a 4.6 sigma emission line in the narrow spectral window expected for NV 1243. The feature is unresolved (FWHM<90 km/s) and is likely nebular in origin. A deep H-band spectrum of EGSY8p7 reveals non-detections of CIV, He II, and OIII]. The presence of NV requires a substantial flux of photons above 77 eV, pointing to a hard ionizing spectrum powered by an AGN or fast radiative shocks. Regardless of its origin, the intense radiation field of EGSY8p7 may aid the transmission of Lyman-alpha through what is likely a partially neutral IGM. With this new detection, five of thirteen known Lyman-alpha emitters at z>7 have now been shown to have intense UV line emission, suggesting that extreme radiation fields are commonplace among the Lyman-alpha population. Future observations with JWST will eventually clarify the origin of these features and explain their role in the visibility of Lyman-alpha in the reionization era.

We present a selection of 24 candidate galaxy-galaxy lensing (GGLs) identified from Hubble images in the outskirts of the massive galaxy clusters from the CLASH survey. These GGLs provide insights into the mass distributions at larger scales than the strong lensing region in the cluster cores. We built parametric mass models for three of these GGLs showing simple lensing configurations, in order to assess the properties of their lens and its environment. We show that the local shear estimated from the GGLs traces the gravitational potential of the clusters at 1-2 arcmin radial distance, allowing us to derive their velocity dispersion. We also find a good agreement between the strength of the shear measured at the GGL positions through strong-lensing modelling and the value derived independently from a weak-lensing analysis of the background sources. Overall, we show the advantages of using single GGL events in the outskirts of clusters to robustly constrain the local shear, even when only photometric redshift estimates are known for the source. We argue that the mass-luminosity scaling relation of cluster members can be tested by modelling the GGLs found around them, and show that the mass parameters can vary up to $\sim$30% between the cluster and GGL models assuming this scaling relation.

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.

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

With the Hubble Frontier Fields program, gravitational lensing has provided a powerful way to extend the study of the ultraviolet luminosity function (LF) of galaxies at $z \sim 6$ down to unprecedented magnitude limits. At the same time, significant discrepancies between different studies were found at the very faint end of the LF. In an attempt to understand such disagreements, we present a comprehensive assessment of the uncertainties associated with the lensing models and the size distribution of galaxies. We use end-to-end simulations from the source plane to the final LF that account for all lensing effects and systematic uncertainties by comparing several mass models. In addition to the size distribution, the choice of lens model leads to large differences at magnitudes fainter than $M_{UV} = -15~$ AB mag, where the magnification factor becomes highly uncertain. We perform MCMC simulations that include all these uncertainties at the individual galaxy level to compute the final LF, allowing, in particular, a crossover between magnitude bins. The best LF fit, using a modified Schechter function that allows for a turnover at faint magnitudes, gives a faint-end slope of $\alpha = -2.01_{-0.14}^{+0.12}$, a curvature parameter of $\beta = 0.48_{-0.25}^{+0.49}$, and a turnover magnitude of $M_{T} = -14.93_{-0.52}^{+0.61}$. Most importantly our procedure shows that robust constraints on the LF at magnitudes fainter than $M_{UV} = -15~$ AB remain unrealistic. More accurate lens modeling and future observations of lensing clusters with the James Webb Space Telescope can reliably extend the UV LF to fainter magnitudes.

A galaxy cluster acts as a cosmic telescope over background galaxies but also as a cosmic microscope of the lens imperfections. The diverging magnification of lensing caustics enhances the microlensing effect of substructure present within the lensing mass. Fine-scale structure can be accessed as a moving background source brightens and disappears when crossing these caustics. The recent recognition of a distant lensed star near the Einstein radius of the galaxy cluster MACSJ1149.5+2223 (Kelly et al. 2017) allows the rare opportunity to reach subsolar mass microlensing through a super-critical column of cluster matter. Here we compare these observations with high-resolution ray-tracing simulations that include stellar microlensing set by the observed intracluster starlight and also primordial black holes that may be responsible for the recently observed LIGO events. We explore different scenarios with microlenses from the intracluster medium and black holes, including primordial ones, and examine strategies to exploit these unique alignments. We find that the best constraints on the fraction of compact dark matter in the small-mass regime can be obtained in regions of the cluster where the intracluster medium plays a negligible role. This new lensing phenomenon should be widespread and can be detected within modest-redshift lensed galaxies so that the luminosity distance is not prohibitive for detecting individual magnified stars. Continuous {\it Hubble Space Telescope} monitoring of several such optimal arcs will be rewarded by an unprecedented mass spectrum of compact objects that can contribute to uncovering the nature of dark matter.

Discovery of strongly-lensed gravitational wave (GW) sources will unveil binary compact objects at higher redshifts and lower intrinsic luminosities than is possible without lensing. Such systems will yield unprecedented constraints on the mass distribution in galaxy clusters, measurements of the polarization of GWs, tests of General Relativity, and constraints on the Hubble parameter. Excited by these prospects, and intrigued by the presence of so-called "heavy black holes" in the early detections by LIGO-Virgo, we commenced a search for strongly-lensed GWs and possible electromagnetic counterparts in the latter stages of the second LIGO observing run (O2). Here, we summarise our calculation of the detection rate of strongly-lensed GWs, describe our review of BBH detections from O1, outline our observing strategy in O2, summarize our follow-up observations of GW170814, and discuss the future prospects of detection.

Strong lensing by massive galaxy clusters can provide magnification of the flux and even multiple images of the galaxies that lie behind them. This phenomenon facilitates observations of high-redshift supernovae (SNe), that would otherwise remain undetected. Type Ia supernovae (SNe Ia) detections are of particular interest because of their standard brightness, since they can be used to improve either cluster lensing models or cosmological parameter measurements. We present a ground-based, near-infrared search for lensed SNe behind the galaxy cluster Abell 370. Our survey was based on 15 epochs of J-band observations with the HAWK-I instrument on the Very Large Telescope (VLT). We use Hubble Space Telescope (HST) photometry to infer the global properties of the multiply-imaged galaxies. Using a recently published lensing model of Abell 370, we also present the predicted magnifications and time delays between the images. In our survey, we did not discover any live SNe from the 13 lensed galaxies with 47 multiple images behind Abell 370. This is consistent with the expectation of $0.09\pm0.02$ SNe calculated based on the measured star formation rate. We compare the expectations of discovering strongly lensed SNe in our survey and that performed with HST during the Hubble Frontier Fields (HFF) programme. We also show the expectations of search campaigns that can be conducted with future facilities, such as the James Webb Space Telescope (JWST) or the Wide-Field Infrared Survey Telescope (WFIRST). We show that the NIRCam instrument aboard the JWST will be sensitive to most SN multiple images in the strongly lensed galaxies and thus will be able to measure their time delays if observations are scheduled accordingly.

We analyse a sample of 8 highly magnified galaxies at redshift 0.6<z<1.5 observed with MUSE, exploring the resolved properties of these galaxies at sub-kiloparsec scales. Combining multi-band HST photometry and MUSE spectra, we derive the stellar mass, global star formation rates, extinction and metallicity from multiple nebular lines, concluding that our sample is representative of z~1 star-forming galaxies. We derive the 2D kinematics of these galaxies from the [OII] emission and model it with a new method that accounts for lensing effects and fits multiple images simultaneously. We use these models to calculate the 2D beam-smearing correction and derive intrinsic velocity dispersion maps. We find them to be fairly homogeneous, with relatively constant velocity dispersions between 15 - 80 km/s and Gini coefficient of <0.3. We do not find any evidence for higher (or lower) velocity dispersions at the positions of bright star-forming clumps. We derive resolved maps of dust attenuation and attenuation-corrected star formation rates from emission lines for two objects in the sample. We use this information to study the relation between resolved star formation rate and velocity dispersion. We find that these quantities are not correlated, and the high-velocity dispersions found for relatively low star-forming densities seems to indicate that, at sub-kiloparsec scales, turbulence in high-z discs is mainly dominated by gravitational instability rather than stellar feedback.

We report the discovery of eMACSJ1341-QG-1, a quiescent galaxy at $z=1.594$ located behind the massive galaxy cluster eMACSJ1341.9$-$2442 ($z=0.835$). The system was identified as a gravitationally lensed triple image in Hubble Space Telescope images obtained as part of a snapshot survey of the most X-ray luminous galaxy clusters at $z>0.5$ and spectroscopically confirmed in ground-based follow-up observations with the ESO/X-Shooter spectrograph. From the constraints provided by the triple image, we derive a first, crude model of the mass distribution of the cluster lens, which predicts a gravitational amplification of a factor of $\sim$30 for the primary image and a factor of $\sim$6 for the remaining two images of the source, making eMACSJ1341-QG-1 by far the most strongly amplified quiescent galaxy discovered to date. Our discovery underlines the power of SNAPshot observations of massive, X-ray selected galaxy clusters for lensing-assisted studies of faint background populations.

We present results from the KMOS lensing survey-KLENS which is exploiting gravitational lensing to study the kinematics of 24 star forming galaxies at $1.4<z<3.5$ with a median mass of $\rm log(M_\star/M_\odot)=9.6$ and median star formation rate (SFR) of $\rm 7.5\,M_\odot\,yr^{-1}$. We find that 25% of these low-mass/low-SFR galaxies are rotation dominated, while the majority of our sample shows no velocity gradient. When combining our data with other surveys, we find that the fraction of rotation dominated galaxies increases with the stellar mass, and decreases for galaxies with a positive offset from the main sequence. We also investigate the evolution of the intrinsic velocity dispersion, $\sigma_0$, as a function of the redshift, $z$, and stellar mass, $\rm M_\star$, assuming galaxies in quasi-equilibrium (Toomre Q parameter equal to 1). From the $z-\sigma_0$ relation, we find that the redshift evolution of the velocity dispersion is mostly expected for massive galaxies ($\rm log(M_\star/M_\odot)>10$). We derive a $\rm M_\star-\sigma_0$ relation, using the Tully-Fisher relation, which highlights that a different evolution of the velocity dispersion is expected depending on the stellar mass, with lower velocity dispersions for lower masses, and an increase for higher masses, stronger at higher redshift. The observed velocity dispersions from this work and from comparison samples spanning $0<z<3.5$ appear to follow this relation, except at higher redshift ($z>2$), where we observe higher velocity dispersions for low masses ($\rm log(M_\star/M_\odot)\sim 9.6$) and lower velocity dispersions for high masses ($\rm log(M_\star/M_\odot)\sim 10.9$) than expected. This discrepancy could, for instance, suggest that galaxies at high-$z$ do not satisfy the stability criterion, or that the adopted parametrisation of the specific star formation rate and molecular properties fail at high redshift.

Motivated by the preponderance of so-called "heavy black holes" in the binary black hole (BBH) gravitational wave (GW) detections to date, and the role that gravitational lensing continues to play in discovering new galaxy populations, we explore the possibility that the GWs are strongly-lensed by massive galaxy clusters. For example, if one of the GW sources were actually located at $z=1$, then the rest-frame mass of the associated BHs would be reduced by a factor $\sim2$. Based on the known populations of BBH GW sources and strong-lensing clusters, we estimate a conservative lower limit on the number of BBH mergers detected per detector year at LIGO/Virgo's current sensitivity that are multiply-imaged, of $R_{\rm detect}\simeq10^{-5}{\rm yr}^{-1}$. This is equivalent to rejecting the hypothesis that one of the BBH GWs detected to date was multiply-imaged at $<\sim4\sigma$. It is therefore unlikely but not impossible that one of the GWs is multiply-imaged. We identify three spectroscopically confirmed strong-lensing clusters with well constrained mass models within the $90\%$ credible sky localisations of the BBH GWs from LIGO's first observing run. In the event that one of these clusters multiply-imaged one of the BBH GWs, we predict that $20-60\%$ of the putative next appearances of the GWs would be detectable by LIGO, and that they would arrive at Earth within three years of first detection.

The "Antennae Galaxy" (NGC 4038/39) is the closest major interacting galaxy system and therefore often taken as merger prototype. We present the first comprehensive integral field spectroscopic dataset of this system, observed with the MUSE instrument at the ESO VLT. We cover the two regions in this system which exhibit recent star-formation: the central galaxy interaction and a region near the tip of the southern tidal tail. In these fields, we detect HII regions and diffuse ionized gas to unprecedented depth. About 15% of the ionized gas was undetected by previous observing campaigns. This newly detected faint ionized gas is visible everywhere around the central merger, and shows filamentary structure. We estimate diffuse gas fractions of about 60% in the central field and 10% in the southern region. We are able to show that the southern region contains a significantly different population of HII regions, showing fainter luminosities. By comparing HII region luminosities with the HST catalog of young star clusters in the central field, we estimate that there is enough Lyman-continuum leakage in the merger to explain the amount of diffuse ionized gas that we detect. We compare the Lyman-continuum escape fraction of each HII region against ionization-parameter sensitive emission line ratios. While we find no systematic trend between these properties, the most extreme line ratios seem to be strong indicators of density bounded ionization. Extrapolating the Lyman-continuum escape fractions to the southern region, we conclude that just from the comparison of the young stellar populations to the ionized gas there is no need to invoke other ionization mechanisms than Lyman-continuum leaking HII regions for the diffuse ionized gas in the Antennae.

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

Giant stellar clumps are ubiquitous in high-redshift galaxies. They are thought to play an important role in the build-up of galactic bulges and as diagnostics of star formation feedback in galactic discs. Hubble Space Telescope (HST) blank field imaging surveys have estimated that these clumps have masses up to 10$^{9.5}$ M$_{\odot}$ and linear sizes larger than ~1 kpc. Recently, gravitational lensing has also been used to get higher spatial resolution. However, both recent lensed observations and models suggest that the clumps properties may be overestimated by the limited resolution of standard imaging techniques. A definitive proof of this observational bias is nevertheless still missing. Here we investigate directly the effect of resolution on clump properties by analysing multiple gravitationally-lensed images of the same galaxy at different spatial resolutions, down to 30 pc. We show that the typical mass and size of giant clumps, generally observed at $\sim$1 kpc resolution in high-redshift galaxies, are systematically overestimated. The high spatial resolution data, only enabled by strong gravitational lensing using currently available facilities, support smaller scales of clump formation by fragmentation of the galactic gas disk via gravitational instabilities.

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.

We present rest-frame Ly$\alpha$ equivalent widths (EW) of 417 Ly$\alpha$ emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at $2.9 < z < 6.6$ in the Hubble Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured EW values taking into account extended Ly$\alpha$ emission and UV continuum slopes ($\beta$). Our LAEs reach unprecedented depths, both in Ly$\alpha$ luminosities and UV absolute magnitudes, from log($L_{\rm Ly\alpha}$/erg s$^{-1}$) $\sim$41.0 to 43.0 and from Muv $\sim$ -16 to -21 (0.01-1.0 $L^{*}_{\rm z=3}$). The EW values span the range of $\sim$ 5 to 240 \AA\ or larger, and their distribution can be well fitted by an exponential law $N = N_{\rm 0}$ exp($-$EW/$w_{\rm 0}$). Owing to the high dynamic range in Muv, we find that the scale factor, $w_{\rm 0}$, depends on Muv in the sense that including fainter Muv objects increases $w_{\rm 0}$, i.e., the Ando effect. The results indicate that selection functions affect the EW scale factor. Taking these effects into account, we find that our $w_{\rm 0}$ values are consistent with those in the literature within $1\sigma$ uncertainties at $2.9 < z < 6.6$ at a given threshold of Muv and $L_{\rm Ly\alpha}$. Interestingly, we find 12 objects with EW $>200$ \AA\ above $1\sigma$ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak CIV $\lambda 1549$ emission line. For the remaining 10 very large EW LAEs, we find that the EW values can be reproduced by young stellar ages ($< 100$ Myr) and low metallicities ($\lesssim 0.02$ $Z_{\rm \odot}$). Otherwise, at least part of the Ly$\alpha$ 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 present a measurement of the fraction of Lyman $\alpha$ (Ly$\alpha$) emitters ($X_{\rm{Ly} \alpha}$) amongst HST continuum-selected galaxies at $3<z<6$ with the Multi-Unit Spectroscopic Explorer (MUSE) on the VLT. Making use of the first 24 MUSE-Wide pointings in GOODS-South, each having an integration time of 1 hour, we detect 100 Ly$\alpha$ emitters and find $X_{\rm{Ly} \alpha}\gtrsim0.5$ for most of the redshift range covered, with 29 per cent of the Ly$\alpha$ sample exhibiting rest equivalent widths (rest-EWs) $\leq$ 15\AA. Adopting a range of rest-EW cuts (0 - 75\AA), we find no evidence of a dependence of $X_{\rm{Ly} \alpha}$ on either redshift or UV luminosity.

We report an integral-field spectroscopic (IFS) observation of a gravitationally lensed spiral galaxy A1689B11 at redshift $z=2.54$. It is the most ancient spiral galaxy discovered to date and the second kinematically confirmed spiral at $z\gtrsim2$. Thanks to gravitational lensing, this is also by far the deepest IFS observation with the highest spatial resolution ($\sim$ 400 pc) on a spiral galaxy at a cosmic time when the Hubble sequence is about to emerge. After correcting for a lensing magnification of 7.2 $\pm$ 0.8, this primitive spiral disk has an intrinsic star formation rate of 22 $\pm$ 2 $M_{\odot}$ yr$^{-1}$, a stellar mass of 10$^{9.8 \pm 0.3}$$M_{\odot}$ and a half-light radius of $r_{1/2}=2.6 \pm 0.7$ kpc, typical of a main-sequence star-forming (SF) galaxy at $z\sim2$. However, the H\alpha\ kinematics show a surprisingly tranquil velocity field with an ordered rotation ($V_{\rm c}$ = 200 $\pm$ 12 km/s) and uniformly small velocity dispersions ($V_{\rm \sigma, mean}$ = 23 $\pm$ 4 km/s and $V_{\rm \sigma, outer-disk}$ = 15 $\pm$ 2 km/s). The low gas velocity dispersion is similar to local spiral galaxies and is consistent with the classic density wave theory where spiral arms form in dynamically cold and thin disks. We speculate that A1689B11 belongs to a population of rare spiral galaxies at $z\gtrsim2$ that mark the formation epoch of thin disks. Future observations with JWST will greatly increase the sample of these rare galaxies and unveil the earliest onset of spiral arms.

We report the detection of extended Lyman-alpha (Lya) 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 at ESO-VLT. Our sample consists of continuum-faint (-15> M_{UV}> -22) Lya emitters (LAEs). Using a 2D, two-component decomposition of Lya emission assuming circular exponential distributions, we measure scale lengths and luminosities of Lya haloes. We find that 80% of our objects having reliable Lya halo measurements show Lya 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 Lya haloes in our sample is ~4.5 kpc. By comparing the maximal detected extent of the Lya emission with the predicted dark matter halo virial radii of simulated galaxies, we show that the detected Lya emission of our selected sample of LAEs probes a significant portion of the cold circum-galactic medium (CGM) of these galaxies (>50% in average). This result shows that there must be significant HI reservoirs in the CGM and reinforces the idea that Lya haloes are ubiquitous around high-redshift Lya emitting galaxies. Our characterization of the Lya haloes indicates that the majority of the Lya flux comes from the halo (~65%) and that their scale lengths seem to be linked to the UV properties of the galaxies. We do not observe a significant Lya halo size evolution with redshift. We also find that the Lya lines cover a large range of full width at half maximum (FWHM) from 118 to 512 km/s. While the FWHM does not seem to be correlated to the Lya scale length, most compact Lya haloes and those that are not detected with high significance tend to have narrower Lya profiles. Finally, we investigate the origin of the extended Lya emission but we conclude that our data do not allow us to disentangle the possible processes.

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 present spatially resolved stellar kinematic maps, for the first time, for a sample of 17 intermediate redshift galaxies (0.2 < z < 0.8). We used deep MUSE/VLT integral field spectroscopic observations in the Hubble Deep Field South (HDFS) and Hubble Ultra Deep Field (HUDF), resulting from ~30h integration time per field, each covering 1'x1' field of view, with ~0.65" spatial resolution. We selected all galaxies brighter than 25mag in the I band and for which the stellar continuum is detected over an area that is at least two times larger than the spatial resolution. The resulting sample contains mostly late-type disk, main-sequence star-forming galaxies with 10^8.5 - 10^10.5 Msun. Using a full-spectrum fitting technique, we derive two-dimensional maps of the stellar and gas kinematics, including the radial velocity V and velocity dispersion sigma. We find that most galaxies in the sample are consistent with having rotating stellar disks with roughly constant velocity dispersions and that the Vrms=sqrt{V^2+sigma^2} of the gas and stars, a scaling proxy for the galaxy gravitational potential, compare well to each other. These spatially resolved observations of intermediate redshift galaxies suggest that the regular stellar kinematics of disk galaxies that is observed in the local Universe was already in place 4 - 7 Gyr ago and that their gas kinematics traces the gravitational potential of the galaxy, thus is not dominated by shocks and turbulent motions. Finally, we build dynamical axisymmetric Jeans models constrained by the derived stellar kinematics for two specific galaxies and derive their dynamical masses. These are in good agreement (within 25%) with those derived from simple exponential disk models based on the gas kinematics. The obtained mass-to-light ratios hint towards dark matter dominated systems within a few effective radii.

We present the results of a multiwavelength investigation of the very X-ray luminous galaxy cluster MACSJ0553.4-3342 ($z = 0.4270$; hereafter MACSJ0553). Combining high-resolution data obtained with the Hubble Space Telescope and the Chandra X-ray Observatory with ground-based galaxy spectroscopy, our analysis establishes the system unambiguously as a binary, post-collision merger of massive clusters. Key characteristics include perfect alignment of luminous and dark matter for one component, a separation of almost 650 kpc (in projection) between the dark-matter peak of the other subcluster and the second X-ray peak, extremely hot gas (k$T > 15$ keV) at either end of the merger axis, a potential cold front in the east, an unusually low gas mass fraction of approximately 0.075 for the western component, a velocity dispersion of $1490_{-130}^{+104}$ km s$^{-1}$, and no indication of significant substructure along the line of sight. We propose that the MACSJ0553 merger proceeds not in the plane of the sky, but at a large inclination angle, is observed very close to turnaround, and that the eastern X-ray peak is the cool core of the slightly less massive western component that was fully stripped and captured by the eastern subcluster during the collision. If correct, this hypothesis would make MACSJ0553 a superb target for a competitive study of ram-pressure stripping and the collisional behaviour of luminous and dark matter during cluster formation.

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.

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.

MUSE (Multi Unit Spectroscopic Explorer) is an integral-field spectrograph mounted on the Very Large Telescope (VLT) in Chile and made available to the European community since October 2014. The Centre de Recherche Astrophysique de Lyon has developed a dedicated software to help MUSE users analyze the reduced data. In this paper we introduce MPDAF, the MUSE Python Data Analysis Framework, based on several well-known Python libraries (Numpy, Scipy, Matplotlib, Astropy) which offers new tools to manipulate MUSE-specific data. We present different examples showing how this Python package may be useful for MUSE data analysis.

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<1.5 (473 objects) and Lya at 2.9<z<6.7 (692 objects). With respect to F775W magnitude, a 50% completeness is reached at 26.5 mag for ultra deep and 25.5 mag for deep fields, and the completeness remains >~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.

We present the MUSE Hubble Ultra Deep Survey, a mosaic of nine MUSE fields covering 90\% of the entire HUDF region with a 10-hour deep exposure time, plus a deeper 31-hour exposure in a single 1.15 arcmin2 field. The improved observing strategy and advanced data reduction results in datacubes with sub-arcsecond spatial resolution (0.65 arcsec at 7000 A) and accurate astrometry (0.07 arcsec rms). We compare the broadband photometric properties of the datacubes to HST photometry, finding a good agreement in zeropoint up to mAB=28 but with an increasing scatter for faint objects. We have investigated the noise properties and developed an empirical way to account for the impact of the correlation introduced by the 3D drizzle interpolation. The achieved 3 sigma emission line detection limit for a point source is 1.5 and 3.1 10-19 erg.s-1.cm-2 for the single ultra-deep datacube and the mosaic, respectively. We extracted 6288 sources using an optimal extraction scheme that takes the published HST source locations as prior. In parallel, we performed a blind search of emission line galaxies using an original method based on advanced test statistics and filter matching. The blind search results in 1251 emission line galaxy candidates in the mosaic and 306 in the ultradeep datacube, including 72 sources without HST counterparts (mAB>31). In addition 88 sources missed in the HST catalog but with clear HST counterparts were identified. This data set is the deepest spectroscopic survey ever performed. In just over 100 hours of integration time, it provides nearly an order of magnitude more spectroscopic redshifts compared to the data that has been accumulated on the UDF over the past decade. The depth and high quality of these datacubes enables new and detailed studies of the physical properties of the galaxy population and their environments over a large redshift range.

We analyze the properties of a multiply-imaged Lyman-alpha (Lya) 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_UV~-16.5), steep UV spectral index (\beta=-2.4+/-0.1), and extreme rest-frame equivalent width of Lya (EW(Lya)=420+180-120 \AA). Two different gravitational lens models predict high magnification (\mu~10--16) for the two detected counter-images, separated by 7", while a predicted third counter-image (\mu~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 Lya profile with an integrated luminosity L(Lya)~10^42 erg/s. The morphology in the HST bands comprises a compact clump (r_e<100 pc) that dominates the Lya and continuum emission and several fainter clumps at projected distances <1 kpc that coincide with an extension of the Lya 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 contribution from AGN to the Lya emission. Fitting of the spectral energy distribution with stellar population models favors a very young (t<10 Myr), low mass (M*~10^6.5 Msun), and metal poor (Z<4x10^-3) stellar population. Its modest star formation rate (SFR~1.0 Msun/yr) implies high specific SFR (sSFR~2.5x10^-7 yr^-1) and SFR density (Sigma_SFR ~ 7-35 Msun/yr/kpc^2). The properties of A370-L57 make it a good representative of the population of galaxies responsible for cosmic reionization.

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.

We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications.

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

We report on the galaxy MACSJ0032-arc at z=3.6314 discovered during the Herschel Lensing snapshot Survey of massive galaxy clusters, and strongly lensed by the cluster MACSJ0032.1+1808. The successful detections of its rest-frame UV, optical, FIR, millimeter, and radio continua, and of its CO emission enable us to characterize, for the first time at such a high redshift, the stellar, dust, and molecular gas properties of a compact star-forming galaxy with a size smaller than 2.5 kpc, a fairly low stellar mass of 4.8x10^9 Msun, and a moderate IR luminosity of 4.8x10^11 Lsun. We find that the bulk of the molecular gas mass and star formation seems to be spatially decoupled from the rest-frame UV emission. About 90% of the total star formation rate is undetected at rest-frame UV wavelengths because of severe obscuration by dust, but is seen through the thermal FIR dust emission and the radio synchrotron radiation. The observed CO(4-3) and CO(6-5) lines demonstrate that high-J transitions, at least up to J=6, remain excited in this galaxy, whose CO spectral line energy distribution resembles that of high-redshift submm galaxies, even though the IR luminosity of MACSJ0032-arc is ten times lower. This high CO excitation is possibly due to the compactness of the galaxy. We find evidence that this high CO excitation has to be considered in the balance when estimating the CO-to-H2 conversion factor. The inferred depletion time of the molecular gas in MACSJ0032-arc supports the decrease in the gas depletion timescale of galaxies with redshift, although to a lesser degree than predicted by galaxy evolution models. Instead, the measured molecular gas fraction as high as 60-79% in MACSJ0032-arc favors the continued increase in the gas fraction of galaxies with redshift as expected, despite the plateau observed between z~1.5 and z~2.5.

A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (HST). Here we describe two unusual transient events that appeared behind MACS0416 in a strongly lensed galaxy at redshift $z = 1.0054 \pm 0.0002$. These transients---designated HFF14Spo-NW and HFF14Spo-SE and collectively nicknamed "Spock"---were faster and fainter than any supernova (SN), but significantly more luminous than a classical nova. They reached peak luminosities of $\sim10^{41}$ erg s$^{-1}$ ($M_{\rm AB} < -14$ mag) in 5 rest-frame days, then faded below detectability in roughly the same time span. Models of the cluster lens suggest that these events may be spatially coincident at the source plane, but are most likely not temporally coincident. We find that HFF14Spo can be explained as a luminous blue variable (LBV), a recurrent nova (RN), or a pair of stellar microlensing events. To distinguish between these hypotheses will require a clarification of the positions of nearby critical curves, along with high-cadence monitoring of the field that could detect new transient episodes in the host galaxy.

We present a clustering analysis of a sample of 238 Ly{$\alpha$}-emitters at redshift 3<z<6 from the MUSE-Wide survey. This survey mosaics extragalactic legacy fields with 1h MUSE pointings to detect statistically relevant samples of emission line galaxies. We analysed the first year observations from MUSE-Wide making use of the clustering signal in the line-of-sight direction. This method relies on comparing pair-counts at close redshifts for a fixed transverse distance and thus exploits the full potential of the redshift range covered by our sample. A clear clustering signal with a correlation length of r0 = 2.9(+1.0/-1.1) Mpc (comoving) is detected. Whilst this result is based on only about a quarter of the full survey size, it already shows the immense potential of MUSE for efficiently observing and studying the clustering of Ly{$\alpha$}-emitters.

At redshift z = 2, when the Universe was just three billion years old, half of the most massive galaxies were extremely compact and had already exhausted their fuel for star formation(1-4). It is believed that they were formed in intense nuclear starbursts and that they ultimately grew into the most massive local elliptical galaxies seen today, through mergers with minor companions(5,6), but validating this picture requires higher-resolution observations of their centres than is currently possible. Magnification from gravitational lensing offers an opportunity to resolve the inner regions of galaxies(7). Here we report an analysis of the stellar populations and kinematics of a lensed z = 2.1478 compact galaxy, which surprisingly turns out to be a fast-spinning, rotationally supported disk galaxy. Its stars must have formed in a disk, rather than in a merger-driven nuclear starburst(8). The galaxy was probably fed by streams of cold gas, which were able to penetrate the hot halo gas until they were cut off by shock heating from the dark matter halo(9). This result confirms previous indirect indications(10-13) that the first galaxies to cease star formation must have gone through major changes not just in their structure, but also in their kinematics, to evolve into present-day elliptical galaxies.

Strong gravitational lensing by galaxy clusters is a fundamental tool to study dark matter and constrain the geometry of the Universe. Recently, the Hubble Space Telescope Frontier Fields programme has allowed a significant improvement of mass and magnification measurements but lensing models still have a residual root mean square between 0.2 arcsec and few arcsec- onds, not yet completely understood. Systematic errors have to be better understood and treated in order to use strong lensing clusters as reliable cosmological probes. We have analysed two simulated Hubble-Frontier-Fields-like clusters from the Hubble Frontier Fields Comparison Challenge, Ares and Hera. We use several estimators (relative bias on magnification, den- sity profiles, ellipticity and orientation) to quantify the goodness of our reconstructions by comparing our multiple models, optimized with the parametric software LENSTOOL , with the input models. We have quantified the impact of systematic errors arising, first, from the choice of different density profiles and configurations and, secondly, from the availability of con- straints (spectroscopic or photometric redshifts, redshift ranges of the background sources) in the parametric modelling of strong lensing galaxy clusters and therefore on the retrieval of cosmological parameters. We find that substructures in the outskirts have a significant im- pact on the position of the multiple images, yielding tighter cosmological contours. The need for wide-field imaging around massive clusters is thus reinforced. We show that competitive cosmological constraints can be obtained also with complex multimodal clusters and that photometric redshifts improve the constraints on cosmological parameters when considering a narrow range of (spectroscopic) redshifts for the sources.

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

We present IFS-RedEx, a spectrum and redshift extraction pipeline for integral-field spectrographs. A key feature of the tool is a wavelet-based spectrum cleaner. It identifies reliable spectral features, reconstructs their shapes, and suppresses the spectrum noise. This gives the technique an advantage over conventional methods like Gaussian filtering, which only smears out the signal. As a result, the wavelet-based cleaning allows the quick identification of true spectral features. We test the cleaning technique with degraded MUSE spectra and find that it can detect spectrum peaks down to S/N = 8 while reporting no fake detections. We apply IFS-RedEx to MUSE data of the strong lensing cluster MACSJ1931.8-2635 and extract 54 spectroscopic redshifts. We identify 29 cluster members and 22 background galaxies with z >= 0.4. IFS-RedEx is open source and publicly available.

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.

We map the lensing-inferred substructure in the first three clusters observed by the Hubble Space Telescope Frontier Fields Initiative (HSTFF): Abell 2744 (z = 0.308), MACSJ0416, (z = 0.396) and MACSJ1149 (z = 0.543). Statistically resolving dark-matter subhaloes down to ~10^{9.5} solar masses, we compare the derived subhalo mass functions (SHMFs) to theoretical predictions from analytical models and with numerical simulations in a Lambda Cold Dark Matter (LCDM) cosmology. Mimicking our observational cluster member selection criteria in the HSTFF, we report excellent agreement in both amplitude and shape of the SHMF over four decades in subhalo mass (10^{9-13} solar masses). Projection effects do not appear to introduce significant errors in the determination of SHMFs from simulations. We do not find evidence for a substructure crisis, analogous to the missing satellite problem in the Local Group, on cluster scales, but rather excellent agreement of the count-matched HSTFF SHMF down to M_{sub halo}/M_{halo} ~ 10^{-5}. However, we do find discrepancies in the radial distribution of sub haloes inferred from HSTFF cluster lenses compared to determinations from simulated clusters. This suggests that although the selected simulated clusters match the HSTFF sample in mass, they do not adequately capture the dynamical properties and complex merging morphologies of these observed cluster lenses. Therefore, HSTFF clusters are likely observed in a transient evolutionary stage that is presently insufficiently sampled in cosmological simulations. The abundance and mass function of dark matter substructure in cluster lenses continues to offer an important test of the LCDM paradigm, and at present we find no tension between model predictions and observations.

We present a VLT/MUSE survey of lensed high-redshift galaxies behind the z=0.77 cluster RCS0224-0002. We study the detailed internal properties of a highly magnified ({\mu}~29) z=4.88 galaxy seen through the cluster. We detect wide-spread nebular CIV{\lambda}{\lambda}1548,1551{\AA} emission from this galaxy as well as a bright Ly{\alpha} halo with a spatially-uniform wind and absorption profile across 12 kpc in the image plane. Blueshifted high- and low-ionisation interstellar absorption indicate the presence of a high-velocity outflow ({\Delta}v~300 km/s) from the galaxy. Unlike similar observations of galaxies at z=2-3, the Ly{\alpha} emission from the halo emerges close to the systemic velocity - an order of magnitude lower in velocity offset than predicted in "shell"-like outflow models. To explain these observations we favour a model of an outflow with a strong velocity gradient, which changes the effective column density seen by the Ly{\alpha} photons. We also search for high-redshift Ly{\alpha} emitters and identify 14 candidates between z=4.8-6.6, including an over-density at z=4.88, of which only one has a detected counterpart in HST/ACS+WFC3 imaging.

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.

We present a MUSE and KMOS dynamical study 405 star-forming galaxies at redshift z=0.28-1.65 (median redshift z=0.84). Our sample are representative of star-forming, main-sequence galaxies, with star-formation rates of SFR=0.1-30Mo/yr and stellar masses M=10^8-10^11Mo. For 49+/-4% of our sample, the dynamics suggest rotational support, 24+/-3% are unresolved systems and 5+/-2% appear to be early-stage major mergers with components on 8-30kpc scales. The remaining 22+/-5% appear to be dynamically complex, irregular (or face-on systems). For galaxies whose dynamics suggest rotational support, we derive inclination corrected rotational velocities and show these systems lie on a similar scaling between stellar mass and specific angular momentum as local spirals with j*=J/M*\propto M^(2/3) but with a redshift evolution that scales as j*\propto M^{2/3}(1+z)^(-1). We identify a correlation between specific angular momentum and disk stability such that galaxies with the highest specific angular momentum, log(j*/M^(2/3))>2.5, are the most stable, with Toomre Q=1.10+/-0.18, compared to Q=0.53+/-0.22 for galaxies with log(j*/M^(2/3))<2.5. At a fixed mass, the HST morphologies of galaxies with the highest specific angular momentum resemble spiral galaxies, whilst those with low specific angular momentum are morphologically complex and dominated by several bright star-forming regions. This suggests that angular momentum plays a major role in defining the stability of gas disks: at z~1, massive galaxies that have disks with low specific angular momentum, appear to be globally unstable, clumpy and turbulent systems. In contrast, galaxies with high specific angular have evolved in to stable disks with spiral structures.

(Abridged) We consider the scientific case for a large aperture (10-12m class) optical spectroscopic survey telescope with a field of view comparable to that of LSST. Such a facility could enable transformational progress in several areas of astrophysics, and may constitute an unmatched capability for decades. Deep imaging from LSST and Euclid will provide accurate photometry for spectroscopic targets beyond the reach of 4m class instruments. Such a facility would revolutionise our understanding of the assembly and enrichment history of the Milky Way and the role of dark matter through chemo-dynamical studies of tens of millions of stars in the Local Group. Emission and absorption line spectroscopy of redshift z=2-5 galaxies can be used to directly chart the evolution of the cosmic web and examine its connection with activity in galaxies. The facility will also have synergistic impact, e.g. in following up live and transpired transients found with LSST, as well as providing targets and the local environmental conditions for follow-up studies with E-ELT and future space missions. Although our study is exploratory, we highlight a specific telescope design with a 5 square degree field of view and an additional focus that could host a next-generation panoramic IFU. We discuss some technical challenges and operational models and recommend a conceptual design study aimed at completing a more rigorous science case in the context of a costed technical design.

The gravitationally-lensed galaxy A1689-zD1 is one of the most distant spectroscopically confirmed sources ($z=7.5$). It is the earliest known galaxy where the interstellar medium (ISM) has been detected; dust emission was detected with the Atacama Large Millimetre Array (ALMA). A1689-zD1 is also unusual among high-redshift dust emitters as it is a sub-L* galaxy and is therefore a good prospect for the detection of gaseous ISM in a more typical galaxy at this redshift. We observed A1689-zD1 with ALMA in bands 6 and 7 and with the Green Bank Telescope (GBT) in band $Q$. To study the structure of A1689-zD1, we map the mm thermal dust emission and find two spatial components with sizes about $0.4-1.7$\,kpc (lensing-corrected). The rough spatial morphology is similar to what is observed in the near-infrared with {\it HST} and points to a perturbed dynamical state, perhaps indicative of a major merger or a disc in early formation. The ALMA photometry is used to constrain the far-infrared spectral energy distribution, yielding a dust temperature ($T_{\rm dust} \sim 35$--$45$\,K for $\beta = 1.5-2$). We do not detect the CO(3-2) line in the GBT data with a 95\% upper limit of 0.3\,mJy observed. We find a slight excess emission in ALMA band~6 at 220.9\,GHz. If this excess is real, it is likely due to emission from the [CII] 158.8\,$\mu$m line at $z_{\rm [CII]} = 7.603$. The stringent upper limits on the [CII]/$L_{\rm FIR}$ luminosity ratio suggest a [CII] deficit similar to several bright quasars and massive starbursts.

We present the results of a 135-arcmin$^2$ search for high-redshift galaxies lensed by 29 clusters from the MAssive Cluster and extended MAssive Cluster Surveys (MACS and eMACS). We use relatively shallow images obtained with the Hubble Space Telescope in four passbands, namely, F606W, F814W, F110W, and F140W. We identify 130 F814W dropouts as candidates for galaxies at $z \le 6$. In order to fit the available broad-band photometry to galaxy spectral energy distribution (SED) templates, we develop a prior for the level of dust extinction at various redshifts. We also investigate the systematic biases incurred by the use of SED-fit software. The fits we obtain yield an estimate of 20 Lyman-break galaxies with photometric redshifts from $z \sim 7$ to 9. In addition, our survey has identified over 100 candidates with a significant probability of being lower-redshift ($z \sim 2$) interlopers. We conclude that even as few as four broad-band filters -- when combined with fitting the SEDs -- are capable of isolating promising objects. Such surveys thus allow one both to probe the bright end ($M_{1500} \le -19$) of the high-redshift UV luminosity function and to identify candidate massive evolved galaxies at lower redshifts.

Massive galaxy clusters are the most violent large scale structures undergoing merger events in the Universe. Based upon their morphological properties in X-rays, they are classified as un-relaxed and relaxed clusters and often host (a fraction of them) different types of non-thermal radio emitting components, viz., haloes, mini-haloes, relics and phoenix within their Intra Cluster Medium (ICM). The radio haloes show steep (alpha = -1.2) and ultra steep (alpha < -1.5) spectral properties at low radio frequencies, giving important insights on the merger (pre or post) state of the cluster. Ultra steep spectrum radio halo emissions are rare and expected to be the dominating population to be discovered via LOFAR and SKA in the future. Further, the distribution of matter (morphological information), alignment of hot X-ray emitting gas from the ICM with the total mass (dark + baryonic matter) and the bright cluster galaxy (BCG) is generally used to study the dynamical state of the cluster. We present here a multi wavelength study on 14 massive clusters from the CLASH survey and show the correlation between the state of their merger in X-ray and spectral properties (1.4 GHz - 150 MHz) at radio wavelengths. Using the optical data we also discuss about the gas-mass alignment, in order to understand the interplay between dark and baryonic matter in massive galaxy clusters.

We present a strong lensing analysis on the massive cluster Abell 370 (A370; z = 0.375), using a combination of deep multi-band Hubble Space Telescope (HST) imaging and Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy. From only two hours of MUSE data, we are able to measure 120 redshifts in the Southern BCG area, including several multiply-imaged lens systems. In total, we increase the number of multiply-imaged systems with a secure redshift from 4 to 15, nine of which are newly discovered. Of these, eight are located at z > 3, greatly extending the redshift range of spectroscopically-confirmed systems over previous work. Using these systems as constraints, we update a parametric lens model of A370, probing the mass distribution from cluster to galaxy scales. Overall, we find that a model with only two cluster- scale dark matter halos (one for each BCG) does a poor job of fitting these new image constraints. Instead, two additional mass clumps -- a central "bar" of mass located between the BCGs, and another clump located within a "crown" of galaxies in the Northern part of the cluster field -- provide significant improvements to the fit. Additional physical evidence suggests these clumps are indeed real features of the system, but with relatively few image constraints in the crown region, this claim is difficult to evaluate from a modeling perspective. Additional MUSE observations of A370 covering the entire strong-lensing region will greatly help these efforts, further improving our understanding of this intriguing cluster.

We present a ground-based near-infrared search for lensed supernovae behind the massive cluster Abell 1689 at z=0.18, one of the most powerful gravitational telescopes that nature provides. Our survey was based on multi-epoch $J$-band observations with the HAWK-I instrument on VLT, with supporting optical data from the Nordic Optical Telescope. Our search resulted in the discovery of five high-redshift, $0.671<z<1.703$, photometrically classified core-collapse supernovae with magnifications in the range $\Delta m$ = $-0.31$ to $-1.58$ mag, as calculated from lensing models in the literature. Thanks to the power of the lensing cluster, the survey had the sensitivity to detect supernovae up to very high-redshifts, $z$$\sim$$3$, albeit for a limited region of space. We present a study of the core-collapse supernova rates for $0.4\leq z< 2.9$, and find good agreement with both previous estimates, and the predictions from the star formation history. During our survey, we also discovered 2 type Ia supernovae in A1689 cluster members, which allowed us to determine the cluster Ia rate to be $0.14^{+0.19}_{-0.09}\pm0.01$ $\rm{SNuB}$$\,h^2$ (SNuB$\equiv 10^{-12} \,\rm{SNe} \, L^{-1}_{\odot,B} yr^{-1}$). The cluster rate normalized by the stellar mass is $0.10^{+0.13}_{-0.06}\pm0.02$ in $\rm SNuM$$\,h^2$ (SNuM$\equiv 10^{-12} \,\rm{SNe} \, M^{-1}_{\odot} yr^{-1}$). Furthermore, we explore the optimal future survey for improving the core-collapse supernova rate measurements at $z\gtrsim2$ using gravitational telescopes, as well as for the detections with multiply lensed images, and find that the planned WFIRST space mission has excellent prospects. Massive clusters can be used as gravitational telescopes to significantly expand the survey range of supernova searches, with important implications for the study of the high-$z$ transient universe.

We report radio imaging and monitoring observations in the frequency range 0.235 - 2.7 GHz during the flaring mode of PKS 2155-304, one of the brightest BL Lac objects. The high sensitivity GMRT observations not only reveal extended kpc-scale jet and FRI type lobe morphology in this erstwhile `extended-core' blazar but also delineate the morphological details, thanks to its arcsec scale resolution. The radio light curve during the end phase of the outburst measured in 2008 shows high variability (8.5%) in the jet emission in the GHz range, compared to the lower core variability (3.2%) seen at the lowest frequencies. The excess of flux density with a very steep spectral index in the MHz range supports the presence of extra diffuse emission at low frequencies. The analysis of multi wavelength (radio/ optical/ gamma-ray) light curves at different radio frequencies confirms the variability of the core region and agrees with the scenario of high energy emission in gamma-rays due to inverse Compton emission from a collimated relativistic plasma jet followed by synchrotron emission in radio. Clearly, these results give an interesting insight of the jet emission mechanisms in blazars and highlight the importance of studying such objects with low frequency radio interferometers like LOFAR and the SKA and its precursor instruments.

[Abridged] We present a robust measurement of the rest-frame UV luminosity function (LF) and its evolution during the peak epoch of cosmic star formation at 1<z<3. We use our deep near ultraviolet imaging from WFC3/UVIS on the Hubble Space Telescope (HST) and existing ACS/WFC and WFC3/IR imaging of three lensing galaxy clusters, Abell 2744 and MACSJ0717 from the Hubble Frontier Field survey and Abell 1689. We use photometric redshifts to identify 780 ultra-faint galaxies with $M_{UV}$<-12.5 AB mag at 1<z<3. From these samples, we identified 5 new, faint, multiply imaged systems in A1689. We compute the rest-frame UV LF and find the best-fit faint-end slopes of $\alpha=-1.56\pm0.04$, $\alpha=-1.72\pm0.04$ and $\alpha=-1.94\pm0.06$ at 1.0<z<1.6, 1.6<z<2.2 and 2.2<z<3.0, respectively. Our results demonstrate that the UV LF becomes steeper from z\sim1.3 to z\sim2.6 with no sign of a turnover down to $M_{UV}=-14$ AB mag. We further derive the UV LFs using the Lyman break "dropout" selection and confirm the robustness of our conclusions against different selection methodologies. Because the sample sizes are so large, and extend to such faint luminosities, the statistical uncertainties are quite small, and systematic uncertainties (due to the assumed size distribution, for example), likely dominate. If we restrict our analysis to galaxies and volumes above > 50% completeness in order to minimize these systematics, we still find that the faint-end slope is steep and getting steeper with redshift, though with slightly shallower (less negative) values ($\alpha=-1.55\pm0.06$, $-1.69\pm0.07$ and $-1.79\pm0.08$ for $z\sim1.3$, 1.9 and 2.6, respectively). Finally, we conclude that the faint star-forming galaxies with UV magnitudes of $-18.5<M_{UV}<-12.5$ covered in this study, produce the majority (55%-60%) of the unobscured UV luminosity density at 1<z<3.

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.

We describe a revised collection of the number of sunspot groups from 1610 to the present. This new collection is based on the work of Hoyt and Schatten (Solar Phys. 179, 189, 1998). The main changes are the elimination of a considerable number of observations during the Maunder Minimum (hereafter, MM) and the inclusion of several long series of observations. Numerous minor changes are also described. Moreover, we have calculated the active-day percentage during the MM from this new collection as a reliable index of the solar activity. Thus, the level of solar activity obtained in this work is greater than the level obtained using the original Hoyt and Schatten data, although it remains compatible with a grand minimum of solar activity. The new collection is available in digital format.

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.

We present a new strong lensing mass reconstruction of the Bullet cluster (1E 0657-56) at z=0.296, based on WFC3 and ACS HST imaging and VLT/FORS2 spectroscopy. The strong lensing constraints underwent substantial revision compared to previously published analysis, there are now 14 (six new and eight previously known) multiply-imaged systems, of which three have spectroscopically confirmed redshifts (including one newly measured from this work). The reconstructed mass distribution explicitly included the combination of three mass components: i) the intra-cluster gas mass derived from X-ray observation, ii) the cluster galaxies modeled by their fundamental plane scaling relations and iii) dark matter. The model that includes the intra-cluster gas is the one with the best Bayesian evidence. This model has a total RMS value of 0.158" between the predicted and measured image positions for the 14 multiple images considered. The proximity of the total RMS to resolution of HST/WFC3 and ACS (0.07-0.15" FWHM) demonstrates the excellent precision of our mass model. The derived mass model confirms the spatial offset between the X-ray gas and dark matter peaks. The fraction of the galaxy halos mass to total mass is found to be f_s=11+/-5% for a total mass of 2.5+/-0.1 x 10^14 solar mass within a 250 kpc radial aperture.

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.

Emission line galaxies are one of the main tracers of the large-scale structure to be targeted by the next-generation dark energy surveys. To provide a better understanding of the properties and statistics of these galaxies, we have collected spectroscopic data from the VVDS and DEEP2 deep surveys and estimated the galaxy luminosity functions (LFs) of three distinct emission lines, [OII] ($0.5 < z < 1.3$), H{\beta} ($0.3 < z < 0.8$) and [OIII] ($0.3 < z < 0.8$). Our measurements are based on 35,639 emission line galaxies and cover a volume of $\sim10^7$Mpc$^3$. We present the first measurement of the H{\beta} LF at these redshifts. We have also compiled LFs from the literature that were based on independent data or covered different redshift ranges, and we fit the entire set over the whole redshift range with analytic Schechter and Saunders models, assuming a natural redshift dependence of the parameters. We find that the characteristic luminosity ($L_*$) and density ($\phi_*$) of all LFs increase with redshift. Using the Schechter model over the redshift ranges considered, we find that, for [OII] emitters, the characteristic luminosity $L_*(z=0.5)=3.2\times10^{41}$ erg/s increases by a factor of $2.7 \pm 0.2$ from z=0.5 to 1.3; for H{\beta} emitters $L_*(z=0.3)=1.3\times10^{41}$ erg/s increases by a factor of $2.0 \pm 0.2$ from z=0.3 to 0.8; and for [OIII] emitters $L_*(z=0.3)=7.3\times10^{41}$ erg/s increases by a factor of $3.5 \pm 0.4$ from z=0.3 to 0.8.

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

The far-infrared fine-structure line [CII] at 1900.5\,GHz is known to be one of the brightest cooling lines in local galaxies, and therefore it has been suggested to be an efficient tracer for star-formation in very high-redshift galaxies. However, recent results for galaxies at $z>6$ have yielded numerous non-detections in star-forming galaxies, except for quasars and submillimeter galaxies. We report the results of ALMA observations of two lensed, star-forming galaxies at $z = 6.029$ and $z=6.703$. The galaxy A383-5.1 (star formation rate [SFR] of 3.2 M$_\odot$ yr$^{-1}$ and magnification of $\mu = 11.4\pm1.9$) shows a line detection with $L_{\rm [CII]} = 8.9\times10^{6}$ L$_\odot$, making it the lowest $L_{\rm [CII]}$ detection at $z>6$. For MS0451-H (SFR = 0.4 M$_\odot$ yr$^{-1}$ and $\mu = 100\pm20$) we provide an upper limit of $L_{\rm [CII]} < 3\times10^{5}$ L$_\odot$, which is 1\,dex below the local SFR-$L_{\rm [CII]}$ relations. The results are consistent with predictions for low-metallicity galaxies at $z>6$, however, other effects could also play a role in terms of decreasing $L_{\rm [CII]}$. The detection of A383-5.1 is encouraging and suggests that detections are possible, but much fainter than initially predicted.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The Frontier Fields are a director's discretionary time campaign with HST and the Spitzer Space Telescope to see deeper into the universe than ever before. The Frontier Fields combine the power of HST and Spitzer with the natural gravitational telescopes of massive high-magnification clusters of galaxies to produce the deepest observations of clusters and their lensed galaxies ever obtained. Six clusters - Abell 2744, MACSJ0416.1-2403, MACSJ0717.5+3745, MACSJ1149.5+2223, Abell S1063, and Abell 370 - were selected based on their lensing strength, sky darkness, Galactic extinction, parallel field suitability, accessibility to ground-based facilities, HST, Spitzer and JWST observability, and pre-existing ancillary data. These clusters have been targeted by the HST ACS/WFC and WFC3/IR with coordinated parallels of adjacent blank fields for over 840 HST orbits. The Spitzer Space Telescope has dedicated > 1000 hours of director's discretionary time to obtain IRAC 3.6 and 4.5 micron imaging to ~26.5, 26.0 ABmag 5-sigma point-source depths in the six cluster and six parallel Frontier Fields. The Frontier Field parallel fields are the second-deepest observations thus far by HST with ~29th ABmag 5-sigma point source depths in seven optical - near-infrared bandpasses. Galaxies behind the Frontier Field cluster lenses experience typical magnification factors of a few, with small regions near the critical curves magnified by factors 10-100. Therefore, the Frontier Field cluster HST images achieve intrinsic depths of ~30-33 magnitudes over very small volumes. Early studies of the Frontier Fields have probed galaxies fainter than any seen before during the epoch of reionization 6 < z < 10, mapped out the cluster dark matter to unprecedented resolution, and followed lensed transient events.

We present a complete census of all Herschel-detected sources within the six massive lensing clusters of the HST Frontier Fields (HFF). We provide a robust legacy catalogue of 263 sources with Herschel fluxes, primarily based on imaging from the Herschel Lensing Survey (HLS) and PEP/HerMES Key Programmes. We optimally combine Herschel, Spitzer and WISE infrared (IR) photometry with data from HST, VLA and ground-based observatories, identifying counterparts to gain source redshifts. For each Herschel-detected source we also present magnification factor (mu), intrinsic IR luminosity and characteristic dust temperature, providing a comprehensive view of dust-obscured star formation within the HFF. We demonstrate the utility of our catalogues through an exploratory overview of the magnified population, including more than 20 background sub-LIRGs unreachable by Herschel without the assistance gravitational lensing.

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 present the results obtained with MUSE on the core of the lensing cluster A1689. Integral-field observations with MUSE provide a unique view of the central region, allowing us to conduct a complete census on both cluster galaxies and lensed background sources, identified based on their spectral features without preselection. We investigate the multiple-image configuration for all known sources in the field. Previous to our survey, 28 different lensed galaxies displaying 46 multiple images were known in the MUSE field of view, most of them based on photometric redshifts and lensing considerations. Among them, we spectroscopically confirm 12 images based on their emission-lines, corresponding to 7 different lensed galaxies between z = 0.95 and 5.0. In addition, 14 new galaxies have been spectroscopically identified in this area, with redshifts ranging between 0.8 and 6.2. All background sources within the MUSE field of view correspond to multiple-imaged systems lensed by A1689. 17 sources in total are found at z > 3 based on their Lyman-alpha emission, with Lyman-alpha luminosities ranging between 40.5 < log(Ly{\alpha}) < 42.5 after correction for magnification. This sample is particularly sensitive to the slope of the LF toward the faintest-end. The density of sources obtained in this survey is consistent with a steep value of {\alpha} < -1.5, although this result still needs further investigation. These results illustrate the efficiency of MUSE in the characterization of lensing clusters on one hand, and the study of faint and distant populations of galaxies on the other hand. In particular, our current survey of lensing clusters should provide a unique census of sources responsible for the reionization in a representative volume at z ~ 4-7.

We discuss spatially resolved emission line spectroscopy secured for a total sample of 15 gravitationally lensed star-forming galaxies at a mean redshift of $z\simeq2$ based on Keck laser-assisted adaptive optics observations undertaken with the recently-improved OSIRIS integral field unit (IFU) spectrograph. By exploiting gravitationally lensed sources drawn primarily from the CASSOWARY survey, we sample these sub-L$^{\ast}$ galaxies with source-plane resolutions of a few hundred parsecs ensuring well-sampled 2-D velocity data and resolved variations in the gas-phase metallicity. Such high spatial resolution data offers a critical check on the structural properties of larger samples derived with coarser sampling using multiple-IFU instruments. We demonstrate how kinematic complexities essential to understanding the maturity of an early star-forming galaxy can often only be revealed with better sampled data. Although we include four sources from our earlier work, the present study provides a more representative sample unbiased with respect to emission line strength. Contrary to earlier suggestions, our data indicates a more diverse range of kinematic and metal gradient behavior inconsistent with a simple picture of well-ordered rotation developing concurrently with established steep metal gradients in all but merging systems. Comparing our observations with the predictions of hydrodynamical simulations suggests that gas and metals have been mixed by outflows or other strong feedback processes, flattening the metal gradients in early star-forming galaxies.

We introduce Zurich Atmosphere Purge (ZAP), an approach to sky subtraction based on principal component analysis (PCA) that we have developed for the Multi Unit Spectrographic Explorer (MUSE) integral field spectrograph. ZAP employs filtering and data segmentation to enhance the inherent capabilities of PCA for sky subtraction. Extensive testing shows that ZAP reduces sky emission residuals while robustly preserving the flux and line shapes of astronomical sources. The method works in a variety of observational situations from sparse fields with a low density of sources to filled fields in which the target source fills the field of view. With the inclusion of both of these situations the method is generally applicable to many different science cases and should also be useful for other instrumentation. ZAP is available for download at http://muse-vlt.eu/science/tools.

[abridged] We present a strong-lensing analysis of MACSJ0717.5+3745, based on the full depth of the Hubble Frontier Field (HFF) observations, which brings the number of multiply imaged systems to 61, ten of which are spectroscopically confirmed. The total number of images comprised in these systems rises to 165. Our analysis uses a parametric mass reconstruction technique, as implemented in the Lenstool software, to constrain a mass distribution composed of four large-scale mass components + galaxy-scale perturbers. We find a superposition of cored isothermal mass components to provide a good fit to the observational constraints, resulting in a very shallow mass distribution for the smooth (large-scale) component. Given the implications of such a flat mass profile, we investigate whether a model composed of "peaky" non-cored mass components can also reproduce the observational constraints. We find that such a non-cored mass model reproduces the observational constraints equally well. Although the total mass distributions of both models are consistent, as well as the integrated two dimensional mass profiles, we find that the smooth and the galaxy-scale components are very different. We conclude that, even in the HFF era, the generic degeneracy between smooth and galaxy-scale components is not broken, in particular in such a complex galaxy cluster. Consequently, insights into the mass distribution of MACS J0717 remain limited, underlining the need for additional probes beyond strong lensing. Our findings also have implications for estimates of the lensing magnification: we show that the amplification difference between the two models is larger than the error associated with either model. This uncertainty decreases the area of the image plane where we can reliably study the high-redshift Universe by 50 to 70%.

We present a high-precision mass model of the galaxy cluster MACSJ1149.6+2223, based on a strong-gravitational-lensing analysis of Hubble Space Telescope Frontier Fields (HFF) imaging data and spectroscopic follow-up with Gemini/GMOS and VLT/MUSE. Our model includes 12 new multiply imaged galaxies, bringing the total to 22, comprised of 65 individual lensed images. Unlike the first two HFF clusters, Abell 2744 and MACSJ0416.1-2403, MACSJ1149 does not reveal as many multiple images in the HFF data. Using the Lenstool software package and the new sets of multiple images, we model the cluster with several cluster-scale dark-matter halos and additional galaxy-scale halos for the cluster members. Consistent with previous analyses, we find the system to be complex, composed of five cluster-scale halos. Their spatial distribution and lower mass, however, makes MACSJ1149 a less powerful lens. Our best-fit model predicts image positions with an RMS of 0.91". We measure the total projected mass inside a 200~kpc aperture as ($1.840\pm 0.006$)$\times 10^{14}$M$_{\odot}$, thus reaching again 1% precision, following our previous HFF analyses of MACSJ0416.1-2403 and Abell 2744. In light of the discovery of the first resolved quadruply lensed supernova, SN Refsdal, in one of the multiply imaged galaxies identified in MACSJ1149, we use our revised mass model to investigate the time delays and predict the rise of the next image between November 2015 and January 2016.

Observations of the cosmic microwave background indicate that baryons account for 5% of the Universe's total energy content. In the local Universe, the census of all observed baryons falls short of this estimate by a factor of two. Cosmological simulations indicate that the missing baryons might not have condensed into virialized haloes, but reside throughout the filaments of the cosmic web (where matter density is larger than average) as a low-density plasma at temperatures of $10^5-10^7$ kelvin, known as the warm-hot intergalactic medium. There have been previous claims of the detection of warm baryons along the line of sight to distant blazars and of hot gas between interacting clusters. These observations were, however, unable to trace the large-scale filamentary structure, or to estimate the total amount of warm baryons in a representative volume of the Universe. Here we report X-ray observations of filamentary structures of gas at $10^7$ kelvin associated with the galaxy cluster Abell 2744. Previous observations of this cluster were unable to resolve and remove coincidental X-ray point sources. After subtracting these, we reveal hot gas structures that are coherent over scales of 8 mergaparsecs. The filaments coincide with over-densities of galaxies and dark matter, with 5-10% of their mass in baryonic gas. This gas has been heated up by the cluster's gravitational pull and is now feeding its core. Our findings strengthen evidence for a picture of the Universe in which a large fraction of the missing baryons reside in the filaments of the cosmic web.

Spatially resolved studies of high redshift galaxies, an essential insight into galaxy formation processes, have been mostly limited to stacking or unusually bright objects. We present here the study of a typical (L$^{*}$, M$_\star$ = 6 $\times 10^9$ $M_\odot$) young lensed galaxy at $z=3.5$, observed with MUSE, for which we obtain 2D resolved spatial information of Ly$\alpha$ and, for the first time, of CIII] emission. The exceptional signal-to-noise of the data reveals UV emission and absorption lines rarely seen at these redshifts, allowing us to derive important physical properties (T$_e\sim$15600 K, n$_e\sim$300 cm$^{-3}$, covering fraction f$_c\sim0.4$) using multiple diagnostics. Inferred stellar and gas-phase metallicities point towards a low metallicity object (Z$_{\mathrm{stellar}}$ = $\sim$ 0.07 Z$_\odot$ and Z$_{\mathrm{ISM}}$ $<$ 0.16 Z$_\odot$). The Ly$\alpha$ emission extends over $\sim$10 kpc across the galaxy and presents a very uniform spectral profile, showing only a small velocity shift which is unrelated to the intrinsic kinematics of the nebular emission. The Ly$\alpha$ extension is $\sim$4 times larger than the continuum emission, and makes this object comparable to low-mass LAEs at low redshift, and more compact than the Lyman-break galaxies and Ly$\alpha$ emitters usually studied at high redshift. We model the Ly$\alpha$ line and surface brightness profile using a radiative transfer code in an expanding gas shell, finding that this model provides a good description of both observables.

We present a pilot study of the z=2.923 radio galaxy MRC0943-242, where we for the first time combine information from ALMA and MUSE data cubes. Even with modest integration times, we disentangle an AGN and a starburst dominated set of components. These data reveal a highly complex morphology, as the AGN, starburst, and molecular gas components show up as widely separated sources in dust continuum, optical continuum and CO line emission observations. CO(1-0) and CO(8-7) line emission suggest that there is a molecular gas reservoir offset from both the dust and the optical continuum that is located ~90kpc from the AGN. The UV line emission has a complex structure in emission and absorption. The line emission is mostly due to i) a large scale ionisation cone energised by the AGN, ii) a Ly-alpha emitting bridge of gas between the radio galaxy and a heavily star-forming set of components. Strangely, the ionisation cone has no Ly-alpha emission. We find this is due to an optically thick layer of neutral gas with unity covering fraction spread out over a region of at least ~100kpc from the AGN. Other, less thick absorption components are associated with Ly-alpha emitting gas within a few tens of kpc from the radio galaxy and are connected by a bridge of emission. We speculate that this linear structure of dust, Ly-alpha and CO emission, and the redshifted absorption seen in the circum-nuclear region may represent an accretion flow feeding gas into this massive AGN host galaxy.

We use deep Hubble Space Telescope imaging of the Frontier Fields to accurately measure the galaxy rest-frame ultraviolet luminosity function (UV LF) in the redshift range $z \sim 6-8$. We combine observations in three lensing clusters A2744, MACS0416, MACS0717 and their associated parallels fields to select high-redshift dropout candidates. We use the latest lensing models to estimate the flux magnification and the effective survey volume in combination with completeness simulations performed in the source plane. We report the detection of 227 galaxy candidates at $z=6-7$ and 25 candidates at $z \sim 8$. While the total survey area is about 4 arcmin$^{2}$ in each parallel field, it drops to about 0.6 to 1 arcmin$^{2}$ in the cluster core fields because of the strong lensing. We compute the UV luminosity function at $z \sim 7$ using the combined galaxy sample and perform Monte Carlo simulations to determine the best fit Schechter parameters. We are able to reliably constrain the LF down to an absolute magnitude of $M_{UV}=-15.25$, which corresponds to 0.005$L^{\star}$. More importantly, we find that the faint-end slope remains steep down to this magnitude limit with $\alpha=-2.04_{-0.17}^{+0.13}$. Our results confirm the most recent results in deep blank fields but extend the LF measurements more than two magnitudes deeper. The UV LF at $z \sim 8$ is not very well constrained below $M_{UV}=-18$ due to the small number statistics and incompleteness uncertainties. To assess the contribution of galaxies to cosmic reionization we derive the UV luminosity density at $z\sim7$ by integrating the UV LF down to an observationally constrained limit of $M_{UV} = -15$. We show that our determination of Log($\rho_{UV}$)=$26.2\pm0.13$ (erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$) can be sufficient to maintain the IGM ionized.

SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 lens models, 13 of which were generated before the SN magnification was known, qualifying as pure "blind tests". The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with the measured mu to within 1-sigma. However, there is a subtle systematic bias: the significant disagreements all involve models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or methodological explanation to account for it. We do find that model accuracy can be improved to some extent with stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. In addition to testing model accuracies as we have done here, Type Ia SN magnifications could also be used as inputs for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional strong- and weak-lensing information is unavailable.

We present a new integral-field spectroscopic dataset of the central part of the Orion Nebula (M 42), observed with the MUSE instrument at the ESO VLT. We reduced the data with the public MUSE pipeline. The output products are two FITS cubes with a spatial size of ~5.9'x4.9' (corresponding to ~0.76 pc x 0.63 pc) and a contiguous wavelength coverage of 4595...9366 Angstrom, spatially sampled at 0.2". We provide two versions with a sampling of 1.25 Angstrom and 0.85 Angstrom in dispersion direction. Together with variance cubes these files have a size of 75 and 110 GiB on disk. They represent one of the largest integral field mosaics to date in terms of information content. We make them available for use in the community. To validate this dataset, we compare world coordinates, reconstructed magnitudes, velocities, and absolute and relative emission line fluxes to the literature and find excellent agreement. We derive a two-dimensional map of extinction and present de-reddened flux maps of several individual emission lines and of diagnostic line ratios. We estimate physical properties of the Orion Nebula, using the emission line ratios [N II] and [S III] (for the electron temperature $T_e$) and [S II] and [Cl III] (for the electron density $N_e$), and show two-dimensional images of the velocity measured from several bright emission lines.

We present a high-precision mass model of galaxy cluster Abell 2744, based on a strong-gravitational-lensing analysis of the \emph{Hubble Space Telescope Frontier Fields} (HFF) imaging data, which now include both \emph{Advanced Camera for Surveys} and \emph{Wide-Field Camera 3} observations to the final depth. Taking advantage of the unprecedented depth of the visible and near-infrared data, we identify 34 new multiply imaged galaxies, bringing the total to 61, comprising 181 individual lensed images. In the process, we correct previous erroneous identifications and positions of multiple systems in the northern part of the cluster core. With the \textsc{Lenstool} software and the new sets of multiple images, we model the cluster using two cluster-scale dark matter halos plus galaxy-scale halos for the cluster members. Our best-fit model predicts image positions with an \emph{RMS} error of 0.69$\arcsec$, which constitutes an improvement by almost a factor of two over previous parametric models of this cluster. We measure the total projected mass inside a 200~kpc aperture as ($2.162\pm 0.005$)$\times 10^{14}M_{\odot}$, thus reaching 1\% level precision for the second time, following the recent HFF measurement of MACSJ0416.1-2403. Importantly, the higher quality of the mass model translates into an overall improvement by a factor of 4 of the derived magnification factor. % for the high-redshift lensed background galaxies. Together with our previous HFF gravitational lensing analysis, this work demonstrates that the HFF data enables high-precision mass measurements for massive galaxy clusters and the derivation of robust magnification maps to probe the early Universe.

We present Keck/MOSFIRE observations of UV metal lines in four bright gravitationally-lensed z~6-8 galaxies behind the cluster Abell 1703. The spectrum of A1703-zd6, a highly-magnified star forming galaxy with a Lyman-alpha redshift of z=7.045, reveals a confident detection of the nebular CIV emission line (unresolved with FWHM < 125 km/s). UV metal lines are not detected in the three other galaxies. At z~2-3, nebular CIV emission is observed in just 1% of UV-selected galaxies. The presence of strong CIV emission in one of the small sample of galaxies targeted in this paper may indicate hard ionizing spectra are more common at z~7. The total estimated equivalent width of the CIV doublet (38 A) and CIV/Lyman-alpha flux ratio (0.3) are comparable to measurements of narrow-lined AGNs. Photoionization models show that the nebular CIV line can also be reproduced by a young stellar population, with very hot metal poor stars dominating the photon flux responsible for triply ionizing carbon. Regardless of the origin of the CIV, we show that the ionizing spectrum of A1703-zd6 is different from that of typical galaxies at z~2, producing more H ionizing photons per unit 1500A luminosity and a larger flux density at 30-50 eV. If such extreme radiation fields are typical in UV-selected systems at z>7, it would indicate that reionization-era galaxies are more efficient ionizing agents than previously thought. Alternatively, we suggest that the small sample of Lyman-alpha emitters at z>7 may trace a rare population with intense radiation fields capable of ionizing their surrounding hydrogen distribution. Additional constraints on high ionization emission lines in galaxies with and without Lyman-alpha detections will help clarify whether hard ionizing spectra are common in the reionization era.

Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system threaded through the cluster core. We find that each of the central galaxies retains a dark matter halo, but that (at least) one of these is spatially offset from its stars. The best-constrained offset is 1.62+/-0.48kpc, where the 68% confidence limit includes both statistical error and systematic biases in mass modelling. Such offsets are not seen in field galaxies, but are predicted during the long infall to a cluster, if dark matter self-interactions generate an extra drag force. With such a small physical separation, it is difficult to definitively rule out astrophysical effects operating exclusively in dense cluster core environments - but if interpreted solely as evidence for self-interacting dark matter, this offset implies a cross-section sigma/m=(1.7+/-0.7)x10^{-4}cm^2/g x (t/10^9yrs)^{-2}, where t is the infall duration.

We present adaptive optics-assisted integral field spectroscopy around the Ha or Hb lines of 12 gravitationally lensed galaxies obtained with VLT/SINFONI, Keck/OSIRIS and Gemini/NIFS. We combine these data with previous observations and investigate the dynamics and star formation properties of 17 lensed galaxies at z = 1-4. Thanks to gravitational magnification of 1.4 - 90x by foreground clusters, effective spatial resolutions of 40 - 700 pc are achieved. The magnification also allows us to probe lower star formation rates and stellar masses than unlensed samples; our target galaxies feature dust-corrected SFRs derived from Ha or Hb emission of 0.8 - 40Msol/yr, and stellar masses M* ~ 4e8 - 6e10 Msol. All of the galaxies have velocity gradients, with 59% consistent with being rotating discs and a likely merger fraction of 29%, with the remaining 12% classed as 'undetermined.' We extract 50 star-forming clumps with sizes in the range 60pc - 1kpc from the Ha (or Hb) maps, and find that their surface brightnesses and their characteristic luminosities evolve to higher luminosities with redshift. We show that this evolution can be described by fragmentation on larger scales in gas-rich discs, and is likely to be driven by evolving gas fractions.

We present deep MUSE integral-field unit (IFU) spectroscopic observations of the giant (~150 x 80 kpc) Ly-alpha halo around the z=4.1 radio galaxy TNJ J1338-1942. This 9-hr observation maps the two-dimensional kinematics of the Ly-alpha emission across the halo. We identify two HI absorbers which are seen against the Ly-alpha emission, both of which cover the full 150 x 80 kpc extent of the halo and so have covering fractions ~1. The stronger and more blue-shifted absorber (dv~1200 km/s) has dynamics that mirror that of the underlying halo emission and we suggest that this high column material (n(HI) ~ 10^19.4 /cm^2), which is also seen in CIV absorption, represents an out-flowing shell that has been driven by the AGN (or star formation) within the galaxy. The weaker (n(HI)~10^14 /cm^2) and less blue shifted (dv~500 km/s) absorber most likely represents material in the cavity between the out-flowing shell and the Ly-alpha halo. We estimate that the mass in the shell must be of order 10^10 Msol -- a significant fraction of the ISM from a galaxy at z=4. The large scales of these coherent structures illustrate the potentially powerful influence of AGN feedback on the distribution and energetics of material in their surroundings. Indeed, the discovery of high-velocity (~1000 km/s), group-halo-scale (i.e. >150 kpc) and mass-loaded winds in the vicinity of the central radio source are broadly in agreement with the requirements of models that invoke AGN-driven outflows to regulate star formation and black-hole growth in massive galaxies at early times.

Candidates for the modest galaxies that formed most of the stars in the early universe, at redshifts $z > 7$, have been found in large numbers with extremely deep restframe-UV imaging. But it has proved difficult for existing spectrographs to characterise them in the UV. The detailed properties of these galaxies could be measured from dust and cool gas emission at far-infrared wavelengths if the galaxies have become sufficiently enriched in dust and metals. So far, however, the most distant UV-selected galaxy detected in dust emission is only at $z = 3.25$, and recent results have cast doubt on whether dust and molecules can be found in typical galaxies at this early epoch. Here we report thermal dust emission from an archetypal early universe star-forming galaxy, A1689-zD1. We detect its stellar continuum in spectroscopy and determine its redshift to be $z = 7.5\pm0.2$ from a spectroscopic detection of the Ly{\alpha} break. A1689-zD1 is representative of the star-forming population during reionisation, with a total star-formation rate of about 12M$_\odot$ yr$^{-1}$. The galaxy is highly evolved: it has a large stellar mass, and is heavily enriched in dust, with a dust-to-gas ratio close to that of the Milky Way. Dusty, evolved galaxies are thus present among the fainter star-forming population at $z > 7$, in spite of the very short time since they first appeared.

Our objectives are to determine the properties of the interstellar medium (ISM) and of star-formation in typical star-forming galaxies at high redshift. Following up on our previous multi-wavelength observations with HST, Spitzer, Herschel, and the Plateau de Bure Interferometer (PdBI), we have studied a strongly lensed z=2.013 galaxy, the arc behind the galaxy cluster MACS J0451+0006, with ALMA to measure the [CII] 158 micron emission line, one of the main coolants of the ISM. [CII] emission from the southern part of this galaxy is detected at 10 $\sigma$. Taking into account strong gravitational lensing, which provides a magnification of $\mu=49$, the intrinsic lensing-corrected [CII]158 micron luminosity is $L(CII)=1.2 \times 10^8 L_\odot$. The observed ratio of [CII]-to-IR emission, $L(CII)/L(FIR) \approx (1.2-2.4) \times 10^{-3}$, is found to be similar to that in nearby galaxies. The same also holds for the observed ratio $L(CII)/L(CO)=2.3 \times 10^3$, which is comparable to that of star-forming galaxies and active galaxy nuclei (AGN) at low redshift. We utilize strong gravitational lensing to extend diagnostic studies of the cold ISM to an order of magnitude lower luminosity ($L(IR) \sim (1.1-1.3) \times 10^{11} L_\odot$) and SFR than previous work at high redshift. While larger samples are needed, our results provide evidence that the cold ISM of typical high redshift galaxies has physical characteristics similar to normal star forming galaxies in the local Universe.

To extend the molecular gas measurements to typical star-forming galaxies (SFGs) with SFR < 40 Msun yr^{-1} and M* < 2.5x10^{10} Msun at z~1.5-3, we have observed CO emission with the IRAM Plateau de Bure Interferometer and 30m telescope for five strongly-lensed galaxies selected from the Herschel Lensing Survey. These observations are combined with a compilation of CO measurements from the literature. We infer the luminosity correction factors r2,1 = 0.81+/-0.20 and r3,1 = 0.57+/-0.15 for the J=2 and J=3 CO transitions, respectively, valid for SFGs at z>1. The combined sample of CO-detected SFGs at z>1 shows a large spread in star formation efficiency (SFE), such that SFE extend beyond the low values of local spirals and overlap the distribution of z>1 sub-mm galaxies. We find that the spread in SFE (or equivalently in molecular gas depletion timescale) is due to primarily the specific star formation rate, but also stellar mass and redshift. Correlations of SFE with the offset from the main-sequence and the compactness of the starburst are less clear. The increase of the molecular gas depletion timescale with M* now revealed by low M* SFGs at z>1 and observed at z=0 is in contrast to the admitted constant molecular gas depletion timescale and the linear Kennicutt-Schmidt relation. We confirm an increase of the molecular gas fraction (fgas) from z~0.2 to z~1.2, followed by a very mild increase toward higher redshifts. At each redshift fgas shows a large dispersion due to the dependence of fgas on M*, producing a gradient of increasing fgas with decreasing M*. We provide the first measure of fgas of z>1 SFGs at the low-M* end (10^{9.4} < M*/Msun < 10^{9.9}), reaching a mean fgas = 0.69+/-0.18, which shows a clear fgas upturn. Finally, we find evidence for a non-universal dust-to-gas ratio among high-redshift SFGs and sub-mm galaxies, local spirals and ULIRGs with near-solar metallicities.

We present the [OII] luminosity function measured in the redshift range 0.1<z<1.65 with unprecedented depth and accuracy. Our measurements are based on medium resolution flux-calibrated spectra of emission line galaxies with the FORS2 instrument at VLT and with the SDSS-III/BOSS spectrograph. The FORS2 spectra and the corresponding catalog containing redshifts and line fluxes are released along with this paper. In this work we use a novel method to combine the aforementioned surveys with GAMA, zCOSMOS and VVDS, which have different target selection, producing a consistent weighting scheme to derive the [OII] luminosity function. The measured luminosity function is in good agreement with previous independent estimates. The comparison with two state-of-the-art semi-analytical models is good, which is encouraging for the production of mock catalogs of [OII] flux limited surveys. We observe the bright end evolution over 8.5 Gyr: we measure the decrease of log L* from 42.4 erg/s at redshift 1.44 to 41.2 at redshift 0.165 and we find that the faint end slope flattens when redshift decreases. This measurement confirms the feasibility of the target selection of future baryonic acoustic oscillation surveys aiming at observing [OII] flux limited samples.

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 a new analysis of HST, Spitzer telescope imaging and VLT imaging and spectroscopic data of a bright lensed galaxy at $z$=1.0334 in the lensing cluster Abell~2667. Using this high-resolution imaging we present an updated lens model that allows us to fully understand the lensing geometry and reconstruct the lensed galaxy in the source plane. This giant arc gives a unique opportunity to peer into the structure of a high-redshift disk galaxy. We find that the lensed galaxy of Abell 2667 is a typical spiral galaxy with morphology similar to the structure of its counterparts at higher redshift $z\sim 2$. The surface brightness of the reconstructed source galaxy in the z$_{850}$ band reveals the central surface brightness $I(0)=20.28\pm0.22$ mag arcsec$^{-2}$ and the characteristic radius $r_s=2.01\pm0.16$ kpc at redshift $z \sim 1$. The morphological reconstruction in different bands shows obvious negative radial color gradients for this galaxy. Moreover, the redder central bulge tends to contain a metal-rich stellar population, rather than being heavily reddened by dust due to high and patchy obscuration. We analyze the VIMOS/IFU spectroscopic data and find that, in the given wavelength range ($\sim 1800-3200$ \AA), the combined arc spectrum of the source galaxy is characterized by a strong continuum emission with strong UV absorption lines (FeII and MgII) and shows the features of a typical starburst Wolf-Rayet galaxy NGC5253. More specifically, we have measured the EWs of FeII and MgII lines in the Abell 2667 spectrum, and obtained similar values for the same wavelength interval of the NGC5253 spectrum. Marginal evidence for CIII] 1909 emission at the edge of the grism range further confirms our expectation.

We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios.

Exploiting the power of gravitational lensing, the Hubble Frontier Fields (HFF) program aims at observing six massive galaxy clusters to explore the distant Universe far beyond the depth limits of blank field surveys. Using the complete Hubble Space Telescope observations of the first HFF cluster Abell 2744, we report the detection of 50 galaxy candidates at $z \sim 7$ and eight candidates at $z \sim 8$ in a total survey area of 0.96 arcmin$^{2}$ in the source plane. Three of these galaxies are multiply-imaged by the lensing cluster. Using an updated model of the mass distribution in the cluster we were able to calculate the magnification factor and the effective survey volume for each galaxy in order to compute the ultraviolet galaxy luminosity function at both redshifts 7 and 8. Our new measurements extend the $z \sim 7$ UV LF down to an absolute magnitude of $M_{UV} \sim -15.5$. We find a characteristic magnitude of $M^{\star}_{UV}=-20.63^{+0.69}_{-0.56}$ mag and a faint-end slope $\alpha = -1.88^{+0.17}_{-0.20}$ close to previous determinations in blank fields. We show here for the first time that this slope remains steep down to very faint luminosities of 0.01$L^{\star}$. Although prone to large uncertainties, our results at $z \sim 8$ also seem to confirm a steep faint-end slope below 0.1$L^{\star}$. The HFF program is therefore providing an extremely efficient way to study the faintest galaxy populations at $z > 7$ that would otherwise be inaccessible with current instrumentation. The full sample of six galaxy clusters will provide yet better constraints on the build-up of galaxies at early epochs and their contribution to cosmic reionization.

We present new observations of the lensing cluster SMACSJ2031.8-4036 obtained with the MUSE integral field spectrograph as part of its commissioning on the Very Large Telescope. By providing medium-resolution spectroscopy over the full 4750-9350 Angstroms domain and a 1x1 arcmin2 field of view, MUSE is ideally suited for identifying lensed galaxies in the cluster core, in particular multiple-imaged systems. We perform a redshift analysis of all sources in the datacube and identify a total of 12 systems ranging from $z=1.46$ to $z=6.4$, with all images of each system confirmed by a spectroscopic redshift. This allows us to accurately constrain the cluster mass profile in this region. We foresee that future MUSE observations of cluster cores should help us discover very faint Lyman-alpha emitters thanks to the strong magnification and the high sensitivity of this instrument.

Deep spectroscopic observations of z~6.5 galaxies have revealed a marked decline with increasing redshift in the detectability of Lyman-alpha emission. While this may offer valuable insight into the end of the reionisation process, it presents a fundamental challenge to the detailed spectroscopic study of the many hundreds of photometrically-selected distant sources now being found via deep HST imaging, and particularly those bright sources viewed through foreground lensing clusters. In this paper we demonstrate the validity of a new way forward via the convincing detection of an alternative diagnostic line, CIII]1909, seen in spectroscopic exposures of two star forming galaxies at z=6.029 and 7.213. The former detection is based on a 3.5 hour X-shooter spectrum of a bright (J=25.2) gravitationally-lensed galaxy behind the cluster Abell 383. The latter detection is based on a 4.2 hour MOSFIRE spectra of one of the most distant spectroscopically confirmed galaxies, GN-108036, with J=25.2. Both targets were chosen for their continuum brightness and previously-known redshift (based on Lyman-alpha), ensuring that any CIII] emission would be located in a favorable portion of the near-infrared sky spectrum. We compare our CIII] and Lyman-alpha equivalent widths in the context of those found at z~2 from earlier work and discuss the motivation for using lines other than Lyman-alpha to study galaxies in the reionisation era.

We present deep spectroscopy of 17 very low mass (M* ~ 2.0x10^6 Msun to 1.4x10^9 Msun) and low luminosity (M_UV ~ -13.7 to -19.9) gravitationally lensed galaxies in the redshift range z~1.5-3.0. Deep rest-frame ultraviolet spectra reveal large equivalent width emission from numerous lines (NIV], OIII], CIV, Si III], CIII]) which are rarely seen in individual spectra of more massive star forming galaxies. CIII] is detected in 16 of 17 low mass star forming systems with rest-frame equivalent widths as large as 13.5 Angstroms. Nebular CIV emission is present in the most extreme CIII] emitters, requiring an ionizing source capable of producing a substantial component of photons with energies in excess of 47.9 eV. Photoionization models support a picture whereby the large equivalent widths are driven by the increased electron temperature and enhanced ionizing output arising from metal poor gas and stars, young stellar populations, and large ionization parameters. The young ages implied by the emission lines and continuum SEDs indicate that the extreme line emitters in our sample are in the midst of a significant upturn in their star formation activity. The low stellar masses, blue UV colors, and large sSFRs of our sample are similar to those of typical z>6 galaxies. Given the strong attenuation of Ly-alpha in z>6 galaxies we suggest that CIII] is likely to provide our best probe of early star forming galaxies with ground-based spectrographs and one of the most efficient means of confirming z>10 galaxies with the James Webb Space Telescope.

Extending over three Hubble Space Telescope (HST) cycles, the Hubble Frontier Fields (HFF) initiative constitutes the largest commitment ever of HST time to the exploration of the distant Universe via gravitational lensing by massive galaxy clusters. We here present models of the mass distribution in the six HFF cluster lenses, derived from a joint strong- and weak-lensing analysis anchored by a total of 88 multiple-image systems identified in existing HST data. The resulting maps of the projected mass distribution and of the gravitational magnification effectively calibrate the HFF clusters as gravitational telescopes. Allowing the computation of search areas in the source plane, these maps are provided to the community to facilitate the exploitation of forthcoming HFF data for quantitative studies of the gravitationally lensed population of background galaxies. Our models of the gravitational magnification afforded by the HFF clusters allow us to quantify the lensing-induced boost in sensitivity over blank-field observations and predict that galaxies at $z>10$ and as faint as m(AB)=32 will be detectable, up to 2 magnitudes fainter than the limit of the Hubble Ultra Deep Field.

We present a high-precision mass model of the galaxy cluster MACSJ0416.1-2403, based on a strong-gravitational-lensing analysis of the recently acquired Hubble Space Telescope Frontier Fields (HFF) imaging data. Taking advantage of the unprecedented depth provided by HST/ACS observations in three passbands, we identify 51 new multiply imaged galaxies, quadrupling the previous census and bringing the grand total to 68, comprising 194 individual lensed images. Having selected a subset of the 57 most securely identified multiply imaged galaxies, we use the Lenstool software package to constrain a lens model comprised of two cluster-scale dark-matter halos and 98 galaxy-scale halos. Our best-fit model predicts image positions with an $RMS$ error of 0.68'', which constitutes an improvement of almost a factor of two over previous, pre-HFF models of this cluster. We find the total projected mass inside a 200~kpc aperture to be $(1.60\pm0.01)\times 10^{14}\ M_\odot$, a measurement that offers a three-fold improvement in precision, reaching the percent level for the first time in any cluster. Finally, we quantify the increase in precision of the derived gravitational magnification of high-redshift galaxies and find an improvement by a factor of $\sim$2.5 in the statistical uncertainty. Our findings impressively confirm that HFF imaging has indeed opened the domain of high-precision mass measurements for massive clusters of galaxies.

We present a comprehensive study of star-forming (SF) galaxies in the HST Frontier Field recent cluster merger A2744 (z=0.308). Wide-field, ultraviolet-infrared (UV-IR) imaging enables a direct constraint of the total star formation rate (SFR) for 53 cluster galaxies, with SFR{UV+IR}=343+/-10 Msun/yr. Within the central 4 arcmin (1.1 Mpc) radius, the integrated SFR is complete, yielding a total SFR{UV+IR}=201+/-9 Msun/yr. Focussing on obscured star formation, this core region exhibits a total SFR{IR}=138+/-8 Msun/yr, a mass-normalised SFR{IR} of Sigma{SFR}=11.2+/-0.7 Msun/yr per 10^14 Msun and a fraction of IR-detected SF galaxies f{SF}=0.080(+0.010,-0.037). Overall, the cluster population at z~0.3 exhibits significant intrinsic scatter in IR properties (total SFR{IR}, Tdust distribution) apparently unrelated to the dynamical state: A2744 is noticeably different to the merging Bullet cluster, but similar to several relaxed clusters. However, in A2744 we identify a trail of SF sources including jellyfish galaxies with substantial unobscured SF due to extreme stripping (SFR{UV}/SFR{IR} up to 3.3). The orientation of the trail, and of material stripped from constituent galaxies, indicates that the passing shock front of the cluster merger was the trigger. Constraints on star formation from both IR and UV are crucial for understanding galaxy evolution within the densest environments.

Using three magnified Type Ia supernovae (SNe Ia) detected behind CLASH clusters, we perform a first pilot study to see whether standardizable candles can be used to calibrate cluster mass maps created from strong lensing observations. Such calibrations will be crucial when next generation HST cluster surveys (e.g. FRONTIER) provide magnification maps that will, in turn, form the basis for the exploration of the high redshift Universe. We classify SNe using combined photometric and spectroscopic observations, finding two of the three to be clearly of type SN Ia and the third probable. The SNe exhibit significant amplification, up to a factor of 1.7 at $\sim5\sigma$ significance (SN-L2). We conducted this as a blind study to avoid fine tuning of parameters, finding a mean amplification difference between SNe and the cluster lensing models of $0.09 \pm 0.09^{stat} \pm 0.05^{sys}$ mag. This impressive agreement suggests no tension between cluster mass models and high redshift standardized SNe Ia. However, the measured statistical dispersion of $\sigma_{\mu}=0.21$ mag appeared large compared to the dispersion expected based on statistical uncertainties ($0.14$). Further work with the supernova and cluster lensing models, post unblinding, reduced the measured dispersion to $\sigma_{\mu}=0.12$. An explicit choice should thus be made as to whether SNe are used unblinded to improve the model, or blinded to test the model. As the lensed SN samples grow larger, this technique will allow improved constraints on assumptions regarding e.g. the structure of the dark matter halo.

The Hubble Frontier Fields (HFF) program combines the capabilities of the Hubble Space Telescope (HST) with the gravitational lensing of massive galaxy clusters to probe the distant Universe to an unprecedented depth. Here, we present the results of the first combined HST and Spitzer observations of the cluster Abell 2744. We combine the full near-infrared data with ancillary optical images to search for gravitationally lensed high-redshift (z > 6) galaxies. We report the detection of 15 I814-dropout candidates at z ~ 6-7 and one Y105-dropout at z ~ 8 in a total survey area of 1.43 arcmin^2 in the source plane. The predictions of our lens model allow us to also identify five multiply-imaged systems lying at redshifts between z ~ 6 and z ~ 8. Thanks to constraints from the mass distribution in the cluster, we were able to estimate the effective survey volume corrected for completeness and magnification effects. This was in turn used to estimate the rest-frame ultraviolet luminosity function (LF) at z ~ 6-8. Our LF results are generally in agreement with the most recent blank field estimates, confirming the feasibility of surveys through lensing clusters. Although based on a shallower observations than what will be achieved in the final dataset including the full ACS observations, the LF presented here extends down to Muv ~ -18.5 at z ~ 7 with one identified object at Muv ~ -15 thanks to the highly-magnified survey areas. This early study forecasts the power of using massive galaxy clusters as cosmic telescopes and its complementarity to blank fields.

We present in this letter the first analysis of a z~8 galaxy candidate found in the Hubble and Spitzer imaging data of Abell 2744, as part of the Hubble Frontier Fields legacy program. We applied the most commonly-used methods to select exceptionally high-z galaxies by combining non-detection and color-criteria using seven HST bands. We used GALFIT on IRAC images for fitting and subtracting contamination of bright nearby sources. The physical properties have been inferred from SED-fitting using templates with and without nebular emission. This letter is focussed on the brightest candidate we found (m$_{F160W}$=26.2) over the 4.9 arcmin$^2$ field of view covered by the WFC3. It shows a non-detection in the ACS bands and at 3.6{\mu}m whereas it is clearly detected at 4.5{\mu}m with rather similar depths. This break in the IRAC data could be explained by strong [OIII]+H{\beta} lines at z~8 which contribute to the 4.5{\mu}m photometry. The best photo-z is found at z~8.0$^{+0.2}_{-0.5}$, although solutions at low-redshift (z~1.9) cannot be completely excluded, but they are strongly disfavoured by the SED-fitting work. The amplification factor is relatively small at {\mu}=1.49$\pm$0.02. The Star Formation Rate in this object is ranging from 8 to 60 Mo/yr, the stellar mass is in the order of M$_{\star}$=(2.5-10) x 10$^{9}$Mo and the size is r~0.35$\pm$0.15 kpc. This object is one of the first z~8 LBG candidates showing a clear break between 3.6{\mu}m and 4.5{\mu}m which is consistent with the IRAC properties of the first spectroscopically confirmed galaxy at a similar redshift. Due to its brightness, the redshift of this object could potentially be confirmed by near infrared spectroscopy with current 8-10m telescopes. The nature of this candidate will be revealed in the coming months with the arrival of new ACS and Spitzer data, increasing the depth at optical and near-IR wavelengths.

We present Submillimeter Array (SMA) [CII] 158um and Jansky Very Large Array (JVLA) $^{12}$CO(1-0) line emission maps for the bright, lensed, submillimeter source at $z=5.2430$ behind Abell 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple $^{12}$CO rotational transitions, [CI], water and [NII], providing some of the best constraints on the properties of the interstellar medium (ISM) in a galaxy at $z>5$. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8-1000um) = (1.6$\pm$0.1)x10^14 L_sun/mu, where the total magnification mu_total = 8.9$\pm$1.9, via a new lens model from the [CII] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components which correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source plane separation <1 kpc. The reddest dominates the continuum map (de-magnified L_FIR = (1.1$\pm$0.2)x10^13 L_sun), and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (dV ~ 500 km/s) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region, and photo-dissociation region (PDR) analysis suggests that it is lower density, cooler and forming stars less vigorously than the other components. Strikingly, it has very strong [NII] emission which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.

MS$\,$0451.6$-$0305 is a rich galaxy cluster whose strong lensing is particularly prominent at submm wavelengths. We combine new SCUBA-2 data with imaging from Herschel SPIRE and PACS and HST in order to try to understand the nature of the sources being lensed. In the region of the "giant submm arc," we uncover seven multiply imaged galaxies (up from the previously known three), of which six are found to be at a redshift of $z\sim2.9$, and possibly constitute an interacting system. Using a novel forward-modelling approach, we are able to simultaneously deblend and fit SEDs to the individual galaxies that contribute to the giant submm arc, constraining their dust temperatures, far infrared luminosities and star formation rates. The submm arc first identified by SCUBA can now be seen to be composed of at least five distinct sources, four of these within the galaxy group at $z\sim2.9$. The total unlensed luminosity for this galaxy group is $(3.1\pm0.3) \times 10^{12}\,\mathrm{L}_\odot$, which gives an unlensed star formation rate of $(450\pm50)$ M$_\odot$ yr$^{-1}$. From the properties of this system, we see no evidence of evolution towards lower temperatures in the dust temperature versus far-infrared luminosity relation for high redshift galaxies.

We identified the z~2 Lyman break galaxies using deep HST ultraviolet (F275W/F336W) imaging of Abell 1689. Because of the imaging depth and the large magnification provided by the cluster, we detect galaxies 100x fainter (-19.5< M_1500 <-13) than previous surveys at this redshift. We are able to calculate the intrinsic sensitivity of the observations as a function of source plane position, allowing determinations of effective volume as a function of luminosity. We fit the faint-end slope of the luminosity function to be alpha = -1.74 +/-0.08, consistent with the values obtained for 2.5 < z < 6. There is no turnover in the luminosity function down to MUV = -13. The trend of increasingly redder UV spectral slopes with luminosity at higher redshifts is observed in our sample, but with redder slopes at all luminosities and average reddening of < E(B - V) >= 0.15. We assume the stars in these galaxies are metal poor (0.2Z_sun) compared to their brighter counterparts (Z_sun), resulting in bluer assumed intrinsic UV slopes and larger derived dust extinction. The total UV luminosity density at z ~ 2 is 4.31x10^26 erg/s/Hz/Mpc^3, more than 70% of which is emitted by galaxies in the luminosity range of our sample. We determine the star formation rate density at z ~ 2 (assuming constant dust extinction correction of 4.2 over all luminosities and a Kroupa IMF) of 0.148 M/yr/Mpc^3, significantly higher than previous determinations because of the additional population of fainter galaxies and the larger dust correction factors.[abridged]

Multi-wavelength, optical to IR/sub-mm observations of 5 strongly lensed galaxies identified by the Herschel Lensing Survey, plus two well-studied lensed galaxies, MS1512-cB58 and the Cosmic Eye, for which we also provide updated Herschel measurements, are used to determine the physical properties of z~1.5-3 star-forming galaxies close to or below the detection limits of blank fields. We constrain their stellar and dust content, determine star formation rates and histories, dust attenuation and extinction laws, and other related properties. We perform SED-fits of the full photometry of each object as well for the optical and infrared parts separately, exploring various parameters, including nebular emission. The IR observations and emission line measurements, where available, are used a posteriori constraints on the models. Besides the various stellar population models we explore, we use the observed IR/UV ratio to estimate the extinction and create "energy conserving models", that constrain most accurately the physical properties of our sources. Our sample has a median lensing-corrected IR luminosity ~ 3e11 Lsun, stellar masses between 2e9 and 2e11 Msun, and IR/UV luminosity ratios spanning a wide range. The dust masses of our galaxies are in the range 2 to 17e7 Msun, extending previous studies at the same redshift down to lower masses. We do not find any particular trend of the dust temperature Tdust with IR luminosity, suggesting an overall warmer dust regime at our redshift regardless of luminosity. Lensing enables us to study the detailed physical properties of individual IR-detected z~1.5-3 galaxies up to a factor ~10 fainter than achieved with deep blank field observations. We demonstrate that multi-wavelength observations combining stellar and dust emission can constrain star formation histories and extinction laws of star-forming galaxies.

In the course of our 870um APEX/LABOCA follow up of the Herschel Lensing Survey we have detected a source in AS1063 (RXC J2248.7-4431), that has no counterparts in any of the Herschel PACS/SPIRE bands, it is a Herschel 'drop-out' with S_870/S_500>0.5. The 870um emission is extended and centered on the brightest cluster galaxy suggesting either a multiply imaged background source or substructure in the Sunyaev-Zel'dovich (SZ) increment due to inhomogeneities in the hot cluster gas of this merging cluster. We discuss both interpretations with emphasis on the putative lensed source. Based on the observed properties and on our lens model we find that this source could be the first SMG with a moderate far infrared luminosity (L_FIR<10^12 L_sol) detected so far at z>4. In deep HST observations we identified a multiply imaged z~6 source and we measured its spectroscopic redshift z=6.107 with VLT/FORS. This source could be associated with the putative SMG but it is most likely offset spatially by 10-30kpc and they could be interacting galaxies. With a FIR luminosity in the range [5-15]x10^{11} L_sol corresponding to a star formation rate in the range [80-260]M_sol/yr, this SMG would be more representative than the extreme starbursts usually detected at z>4. With a total magnification of ~25 it would open a unique window to the 'normal' dusty galaxies at the end of the epoch of reionization.

Quiescent massive galaxies at z~2 are thought to be the progenitors of present-day massive ellipticals. Observations revealed them to be extraordinarily compact. The determination of stellar ages, star formation rates and dust properties via spectroscopic measurements has up to now only been feasible for the most luminous and massive specimens (~3x M*). Here we present a spectroscopic study of two near-infrared selected galaxies which are close to the characteristic stellar mass M* (~0.9x M* and ~1.3x M*) and whose observed brightness has been boosted by the gravitational lensing effect. We measure the redshifts of the two galaxies to be z=1.71\pm0.02 and z=2.15\pm0.01. By fitting stellar population synthesis models to their spectro-photometric SEDs we determine their ages to be 2.4^{+0.8}_{-0.6} Gyr and 1.7\pm0.3 Gyr, respectively, which implies that the two galaxies have higher mass-to-light ratios than most quiescent z~2 galaxies in other studies. We find no direct evidence for active star-formation or AGN activity in either of the two galaxies, based on the non-detection of emission lines. Based on the derived redshifts and stellar ages we estimate the formation redshifts to be z=4.3^{+3.4}_{-1.2} and z=4.3^{+1.0}_{-0.6}, respectively. We use the increased spatial resolution due to the gravitational lensing to derive constraints on the morphology. Fitting Sersic profiles to the de-lensed images of the two galaxies confirms their compactness, with one of them being spheroid-like, and the other providing the first confirmation of a passive lenticular galaxy at a spectroscopically derived redshift z~2.

We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2<z<4. We construct a deep K<24 sample of 220000 galaxies selected using the UltraVISTA DR1 data release. Our analysis is based on precise 30-band photometric redshifts. By comparing these photometric redshifts with 10800 spectroscopic redshifts from the zCOSMOS bright and faint surveys, we find a precision of sigma(dz/(1+z))=0.008 at i<22.5 and sigma(dz/(1+zs))=0.03 at 1.5<z<4. We derive the stellar mass function and correct for the Eddington bias. We find a mass-dependent evolution of the global and star-forming populations. This mass-dependent evolution is a direct consequence of the star formation being quenched in galaxies more massive than M>10^10.7Msun. For the mass function of the quiescent galaxies, we do not find any significant evolution of the high-mass end at z<1; however we observe a clear flattening of the faint-end slope. From z~3 to z~1, the density of quiescent galaxies increases over the entire mass range. Their comoving stellar mass density increases by 1.6 dex between z~3 and z~1 and by less than 0.2dex at z<1. We infer the star formation history from the mass density evolution and we find an excellent agreement with instantaneous star formation rate measurements at z<1.5, while we find differences of 0.2dex at z>1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 10^10Msun galaxies increases continuously in the redshift range 1<z<4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO J=4-3, and double that in CO J=1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks -- a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have M(dyn) of several x 10^11 Msun, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star-formation rates place them ~>5x above the main sequence, which supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile with a simple volumetric star-formation law. N-body and dark matter simulations suggest this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster.

We present rest-frame optical spectra for a sample of 9 low-mass star-forming galaxies in the redshift range 1.5 < z < 3 which are gravitationally lensed by foreground clusters. We used Triplespec, an echelle spectrograph at the Palomar 200-inch telescope that is very effective for this purpose as it samples the entire near-infrared spectrum simultaneously. By measuring the flux of nebular emission lines we derive gas phase metallicities and star formation rates, and by fitting the optical to infrared spectral energy distributions we obtain stellar masses. Taking advantage of the high magnification due to strong lensing we are able to probe the physical properties of galaxies with stellar masses in the range 7.8 < log M/Msun < 9.4 whose star formation rates are similar to those of typical star-forming galaxies in the local universe. We compare our results with the locally determined relation between stellar mass, gas metallicity and star formation rate. Our data are in excellent agreement with this relation, with an average offset <Delta log O/H> = 0.01 +/- 0.08, suggesting a universal relationship. Remarkably, the scatter around the fundamental metallicity relation is only 0.24 dex, smaller than that observed locally at the same stellar masses, which may provide an important additional constraint for galaxy evolution models.

Clusters of galaxies are excellent locations to probe the distribution of baryons and dark matter (DM) over a wide range of scales. We study a sample of seven massive, relaxed galaxy clusters with centrally-located brightest cluster galaxies (BCGs) at z=0.2-0.3. Using the observational tools of strong and weak gravitational lensing, combined with resolved stellar kinematics within the BCG, we measure the total radial density profile, comprising both dark and baryonic matter, over scales of ~3-3000 kpc. Lensing-derived mass profiles typically agree with independent X-ray estimates within ~15%, suggesting that departures from hydrostatic equilibrium are small and that the clusters in our sample (except A383) are not strongly elongated along the line of sight. The inner logarithmic slope gamma_tot of the total density profile measured over r/r200=0.003-0.03, where rho_tot ~ r^(-gamma_tot), is found to be nearly universal, with a mean <gamma_tot> = 1.16 +- 0.05 (random) +0.05-0.07 (systematic) and an intrinsic scatter of < 0.13 (68% confidence). This is further supported by the very homogeneous shape of the observed velocity dispersion profiles, obtained via Keck spectroscopy, which are mutually consistent after a simple scaling. Remarkably, this slope agrees closely with numerical simulations that contain only dark matter, despite the significant contribution of stellar mass on the scales we probe. The Navarro-Frenk-White profile characteristic of collisionless cold dark matter is a better description of the total mass density at radii >~ 5-10 kpc than that of dark matter alone. Hydrodynamical simulations that include baryons, cooling, and feedback currently provide a poorer match. We discuss the significance of our findings for understanding the assembly of BCGs and cluster cores, particularly the influence of baryons on the inner DM halo. [abridged]

The RCS 2319+00 supercluster is a massive supercluster at z=0.9 comprising three optically selected, spectroscopically confirmed clusters separated by <3 Mpc on the plane of the sky. This supercluster is one of a few known examples of the progenitors of present-day massive clusters (10^{15} Msun by z~0.5). We present an extensive spectroscopic campaign carried out on the supercluster field resulting, in conjunction with previously published data, in 1961 high confidence galaxy redshifts. We find 302 structure members spanning three distinct redshift walls separated from one another by ~65 Mpc. The component clusters have spectroscopic redshifts of 0.901, 0.905 and 0.905. The velocity dispersions are consistent with those predicted from X-ray data, giving estimated cluster masses of ~10^{14.5} - 10^{14.9} Msun. The Dressler-Shectman test finds evidence of substructure in the supercluster field and a friends-of-friends analysis identified 5 groups in the supercluster, including a filamentary structure stretching between two cluster cores previously identified in the infrared by Coppin et al. (2012). The galaxy colors further show this filamentary structure to be a unique region of activity within the supercluster, comprised mainly of blue galaxies compared to the ~43-77% red-sequence galaxies present in the other groups and cluster cores. Richness estimates from stacked luminosity function fits results in average group mass estimates consistent with ~10^{13} Msun halos. Currently, 22% of our confirmed members reside in >~10^{13} Msun groups/clusters destined to merge onto the most massive cluster, in agreement with the massive halo galaxy fractions important in cluster galaxy pre-processing in N-body simulation merger tree studies.

We present VLT spectroscopic observations of 7 discovered galaxy groups between 0.3<z<0.7. The groups were selected from the Strong Lensing Legacy Survey (SL2S), a survey that consists in a systematic search for strong lensing systems in the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). We give details about the target selection, spectroscopic observations and data reduction for the first release of confirmed SL2S groups. The dynamical analysis of the systems reveals that they are gravitationally bound structures, with at least 4 confirmed members and velocity dispersions between 300 and 800 km/s. Their virial masses are between 10^13 and 10^14 M_sun, and so can be classified as groups or low mass clusters. Most of the systems are isolated groups, except two of them that show evidence of an ongoing merger of two sub-structures. We find a good agreement between the velocity dispersions estimated from the analysis of the kinematics of group galaxies and the weak lensing measurements, and conclude that the dynamics of baryonic matter is a good tracer of the total mass content in galaxy groups.

We present a comprehensive observational study of the gas phase metallicity of star-forming galaxies from z ~ 0 -> 3. We combine our new sample of gravitationally lensed galaxies with existing lensed and non-lensed samples to conduct a large investigation into the mass-metallicity (MZ) relation at z > 1. We apply a self-consistent metallicity calibration scheme to investigate the metallicity evolution of star-forming galaxies as a function of redshift. The lensing magnification ensures that our sample spans an unprecedented range of stellar mass (3*10^{7}-6*10^{10} M_sun). We find that at the median redshift of z=2.07, the median metallicity of the lensed sample is 0.35 dex lower than the local SDSS star-forming galaxies and 0.18 dex lower than the z ~ 0.8 DEEP2 galaxies. We also present the z ~ 2 MZ relation using 19 lensed galaxies. A more rapid evolution is seen between z ~ 1->3 than z ~ 0 -> 1 for the high-mass galaxies (10^{9.5-11} M_sun), with almost twice as much enrichment between z ~ 1 -> 3 than between z ~ 1 -> 0. We compare this evolution with the most recent cosmological hydrodynamic simulations with momentum driven winds. We find that the model metallicity is consistent with the observed metallicity within the observational error for the low mass bins. However, for higher masses, the model over-predicts the metallicity at all redshifts. The over-prediction is most significant in the highest mass bin of 10^{10-11} M_sun.

We present results of a combined X-ray/optical analysis of the dynamics of the massive cluster MACS J0358.8-2955 (z=0.428) based on observations with the Chandra X-ray Observatory, the Hubble Space Telescope, and the Keck-I telescope on Mauna Kea. MACS J0358.8-2955 is found to be one of the most X-ray luminous clusters known at z>0.3, featuring L_X(<r_500) = 4.24*10^45 erg/s, kT = (9.55 +0.58/-0.37) keV, M^{3D}_{gas}(<r_500) = (9.18+/-1.45)*10^13 M_sun, and M^{3D}_{tot}(<r_500) = (1.12+/-0.18)*10^15 M_sun. The system's high velocity dispersion of (1440 +130/-110) km/s (890 km/s when the correct relativistic equation is used), however, is inflated by infall along the line of sight, as the result of a complex merger of at least three sub-clusters. One collision proceeds close to head-on, while the second features a significant impact parameter. The temperature variations in the intra-cluster gas, two tentative cold fronts, the radial velocities measured for cluster galaxies, and the small offsets between collisional and non-collisional cluster components all suggest that both merger events are observed close to core passage and along axes that are greatly inclined with respect to the plane of the sky. A strong-lensing analysis of the system anchored upon three triple-image systems (two of which have spectroscopic redshifts) yields independent constraints on the mass distribution. For a gas fraction of 8.2%, the resulting strong-lensing mass profile is in good agreement with our X-ray estimates, and the details of the mass distribution are fully consistent with our interpretation of the three-dimensional merger history of this complex system.

Time delays of gravitationally lensed sources can be used to constrain the mass model of a deflector and determine cosmological parameters. We here present an analysis of the time-delay distribution of multiply imaged sources behind 17 strong lensing galaxy clusters with well-calibrated mass models. We find that for time delays less than 1000 days, at z=3.0, their logarithmic probability distribution functions are well represented by P (log \Delta t)=5.3 x 10^-4 \Delta t^\beta M_250^-2\beta, with \beta=0.77, where M_250 is the projected cluster mass inside 250 kpc (in 10^14 M_sun), and \beta is the power-law slope of the distribution. The resultant probability distribution function enables us to estimate the time-delay distribution in a lensing cluster of known mass. For a cluster with M_250=2 x 10^14 M_sun, the fraction of time delays less than 1000 days is approximately 3%. Taking Abell 1689 as an example, its dark halo and brightest galaxies, with central velocity dispersions larger than 500 km/s, mainly produce large time delays, while galaxy-scale mass clumps are responsible for generating smaller time delays. We estimate the probability of observing multiple images of a supernova in the known images of Abell 1689. A two-component model of estimating the supernova rate is applied in this work. For a magnitude threshold of m_AB=26.5, the yearly rate of Type Ia (core-collapse) supernovae with time delays less than 1000 days is 0.004 +- 0.002 (0.029 +- 0.001). If the magnitude threshold is lowered to m_AB ~ 27.0, the rate of core-collapse supernovae suitable for time delay observation is 0.044 +- 0.015 per year.

We present the results of a search for molecular gas emission from a star-forming galaxy at z = 4.9. The galaxy benefits from magnification of 22 +/- 5x due to strong gravitational lensing by the foreground cluster MS1358+62. We target the CO(5-4) emission at a known position and redshift from existing Hubble Space Telescope/ACS imaging and Gemini/NIFS [OII]3727 imaging spectroscopy, and obtain a tentative detection at the 4.3sigma level with a flux of 0.104 +/- 0.024Jkm/s. From the CO line luminosity and assuming a CO-to-H2 conversion factor alpha=2, we derive a gas mass M_gas ~ 1^{+1}_{-0.6} x 10^9 M_sun. Combined with the existing data, we derive a gas fraction Mgas/(Mgas + M*) = 0.59^{+0.11}_{-0.06}. The faint line flux of this galaxy highlights the difficulty of observing molecular gas in representative galaxies at this epoch, and suggests that routine detections of similar galaxies in the absence of gravitational lensing will remain challenging even with ALMA in full science operations.

MUSE, the Multi Unit Spectroscopic Explorer, is a 2nd generation integral-field spectrograph under final assembly to see first light at the Very Large Telescope in 2013. By capturing ~ 90000 optical spectra in a single exposure, MUSE represents a challenge for data reduction and analysis. We summarise here the main features of the Data Reduction System, as well as some of the tools under development by the MUSE consortium and the DAHLIA team to handle the large MUSE datacubes (about 4x?10^8 pixels) to recover the original astrophysical signal.

We present HST/WFC3 narrowband imaging of the H-alpha emission in a sample of eight gravitationally-lensed galaxies at z = 1 - 1.5. The magnification caused by the foreground clusters enables us to obtain a median source plane spatial resolution of 360pc, as well as providing magnifications in flux ranging from ~10x to ~50x. This enables us to identify resolved star-forming HII regions at this epoch and therefore study their H-alpha luminosity distributions for comparisons with equivalent samples at z ~ 2 and in the local Universe. We find evolution in the both luminosity and surface brightness of HII regions with redshift. The distribution of clump properties can be quantified with an HII region luminosity function, which can be fit by a power law with an exponential break at some cut-off, and we find that the cut-off evolves with redshift. We therefore conclude that `clumpy' galaxies are seen at high redshift because of the evolution of the cut-off mass; the galaxies themselves follow similar scaling relations to those at z = 0, but their HII regions are larger and brighter and thus appear as clumps which dominate the morphology of the galaxy. A simple theoretical argument based on gas collapsing on scales of the Jeans mass in a marginally unstable disk shows that the clumpy morphologies of high-z galaxies are driven by the competing effects of higher gas fractions causing perturbations on larger scales, partially compensated by higher epicyclic frequencies which stabilise the disk.

We present the spatially resolved emission line ratio properties of a ~10^10 M_sun star-forming galaxy at redshift z=1.03. This galaxy is gravitationally lensed as a triple-image giant arc behind the massive lensing cluster Abell 2667. The main image of the galaxy has magnification factors of 14+/-2.1 in flux and ~ 2 by 7 in area, yielding an intrinsic spatial resolution of 115-405 pc after AO correction with OSIRIS at KECK II. The HST morphology shows a clumpy structure and the H\alpha\ kinematics indicates a large velocity dispersion with V_{max} sin(i)/\sigma ~ 0.73, consistent with high redshift disk galaxies of similar masses. From the [NII]/H\alpha\ line ratios, we find that the central 350 parsec of the galaxy is dominated by star formation. The [NII]/H\alpha\ line ratios are higher in the outer-disk than in the central regions. Most noticeably, we find a blue-shifted region of strong [NII]/H\alpha\ emission in the outer disk. Applying our recent HII region and slow-shock models, we propose that this elevated [NII]/H\alpha\ ratio region is contaminated by a significant fraction of shock excitation due to galactic outflows. Our analysis suggests that shocked regions may mimic flat or inverted metallicity gradients at high redshift.

Strong gravitational lensing magnifies the flux from distant galaxies, allowing us to detect emission lines that would otherwise fall below the detection threshold for medium-resolution spectroscopy. Here we present the detection of temperature-sensitive oxygen emission lines from three galaxies at 2<z<3.5, which enables us to directly determine the oxygen abundances and thereby double the number of galaxies at z>2 for which this has been possible. The three galaxies have ~10% solar oxygen abundances in agreement with strong emission line diagnostics. Carbon and nitrogen ratios relative to oxygen are sub-solar as expected for young metal-poor galaxies. Two of the galaxies are Lya emitters with rest-frame equivalent widths of 20 A and 40 A, respectively, and their high magnification factors allow us for the first time to gain insight into the physical characteristics of high-redshift Lya emitters. Using constraints from the physical properties of the galaxies, we accurately reproduce their line profiles with radiative transfer models. The models show a relatively small outflow in agreement with the observed small velocity offsets between nebular emission and interstellar absorption lines.

We present VLT/X-shooter spectra of 13 galaxies in the redshift range 1< z < 6, which are strongly lensed by massive galaxy clusters. Spectroscopic redshifts are measured for nine galaxies, while three sources have redshifts determined from continuum breaks in their spectra. The stellar masses of the galaxies span four orders of magnitude between 10^7 and 10^11 M_sun and have luminosities at 1500 A rest-frame between 0.004 and 9 L^* after correcting for the magnification. This allows us to probe a variety of galaxy types from young, low-mass starburst galaxies to massive evolved galaxies. The lensed galaxies with stellar masses less than 10^10 M_sun have a large scatter compared to the fundamental relation between stellar mass, star formation rates and oxygen abundances. We provide a modified fit to the fundamental relation for low-mass, low-metallicity galaxies with a weaker dependence of the metallicity on either the star formation rate or stellar mass compared to low-redshift, high-mass and high-metallicity SDSS galaxies.

We report the first weak-lensing detection of a large-scale filament funneling matter onto the core of the massive galaxy cluster MACSJ0717.5+3745. Our analysis is based on a mosaic of 18 multi-passband images obtained with ACS aboard the HST, covering an area of \sim 10x20 arcmin^2. We use a weak-lensing pipeline developed for the COSMOS survey, modified for the analysis of galaxy clusters, to produce a weak-lensing catalogue. A mass map is then computed by applying a weak-gravitational-lensing multi-scale reconstruction technique designed to describe irregular mass distributions such as the one investigated here. We test the resulting mass map by comparing the mass distribution inferred for the cluster core with the one derived from strong-lensing constraints and find excellent agreement. The filament is detected within the 3 sigma detection contour of the lensing mass reconstruction, and underlines the importance of filaments for theoretical and numerical models of the mass distribution in the Cosmic Web. We measure the filament's projected length as \sim 4.5 h_{74}^{-1} Mpc, and its mean density as (2.92 \pm 0.66)10^8 h_{74} M_{\odot} kpc^{-2}. Combined with the redshift distribution of galaxies obtained after an extensive spectroscopic follow-up in the area, we can rule out any projection effect resulting from the chance alignment on the sky of unrelated galaxy group-scale structures. Assuming plausible constraints concerning the structure's geometry based on its galaxy velocity field, we construct a 3D model of the large-scale filament. Within this framework, we derive the three-dimensional length of the filament to be 18 h_{74}^{-1} Mpc, and a deprojected density in terms of the critical density of the Universe of (206 \pm 46) \rho_{crit}, a value that lies at the very high end of the range predicted by numerical simulations.

We investigate the physical properties, geometry and dynamics of the massive cluster merger MACS J0140.0-0555 (z=0.451) using X-ray and optical diagnostics. Featuring two galaxy overdensities separated by about 250 kpc in projection on the sky, and a single peak in the X-ray surface brightness distribution located between them, MACS J0140.0-0555 shows the tell-tale X-ray/optical morphology of a binary, post-collision merger. Our spectral analysis of the X-ray emission, as measured by our Chandra ACIS-I observation of the system, finds the intra-cluster medium to be close to isothermal (~8.5 keV) with no clear signs of cool cores or shock fronts. Spectroscopic follow-up of galaxies in the field of MACS J0140.0-0555 yields a velocity dispersion of 875 (+70/-100) km/s (n_z=66) and no significant evidence of bimodality or substructure along the line of sight. In addition, the difference in radial velocity between the brightest cluster galaxies of the two sub-clusters of 144+/-25 km/s is small compared to typical collision velocities of several 1000 km/s. A strongly lensed background galaxy at z=0.873 (which features variable X-ray emission from an active nucleus) provides the main constraint on the mass distribution of the system. We measure M(<75 kpc) = (5.6+/- 0.5)*10^13 M_sun for the north-western cluster component and a much less certain estimate of (1.5-3)*10^13 M_sun for the south-eastern subcluster. These values are in good agreement with our X-ray mass estimates which yield a total mass of MACS J0140.0-0555 of M(<r_500) ~ (6.8-9.1)*10^14 M_sun. ......

We present and discuss measurements of the gas-phase metallicity gradient in gravitationally lensed galaxies at z=2.0-2.4 based on adaptive optics-assisted imaging spectroscopy with the Keck II telescope. Through deep exposures we have secured high signal to noise data for four galaxies with well-understood kinematic properties. Three galaxies with well-ordered rotation reveal metallicity gradients in the sense of having lower gas-phase metallicities at larger galactocentric radii. Two of these display gradients much steeper than found locally, while a third has one similar to that seen in local disk galaxies. The fourth galaxy exhibits complex kinematics indicative of an ongoing merger and reveals an "inverted" gradient with lower metallicity in the central regions. By comparing our sample to similar data in the literature for lower redshift galaxies, we determine that, on average, metallicity gradients must flatten by a factor of 2.6 +/- 0.9 between z=2.2 and the present epoch. This factor is in rough agreement with the size growth of massive galaxies suggesting that inside-out growth can account for the evolution of metallicity gradients. Since the addition of our new data provides the first indication of a coherent picture of this evolution, we develop a simple model of chemical evolution to explain the collective data. We find that metallicity gradients and their evolution can be explained by the inward radial migration of gas together with a radial variation in the mass loading factor governing the ratio of outflowing gas to the local star formation rate. Average mass loading factors of \lsim 2 are inferred from our model in good agreement with direct measurements of outflowing gas in z \simeq 2 galaxies.

We present results of a strong-lensing analysis of MACS J0717.5+3745 (hereafter MACS J0717), an extremely X-ray luminous galaxy cluster at z=0.55. Observations at different wavelengths reveal a complex and dynamically very active cluster, whose core is connected to a large scale filament extended over several Mpc. Using multi-passband imaging data obtained with the Hubble Space Telescope Advanced Camera for Surveys (ACS), we identify 15 multiply imaged systems across the full field of view of ACS, five of which we confirmed spectroscopically in groundbased follow-up observations with the Keck telescope. We use these multiply imaged systems to constrain a parametric model of the mass distribution in the cluster core, employing a new parallelized version of the Lenstool software. The main result is that the most probable description of the mass distribution comprises four cluster-scale dark matter haloes. The total mass distribution follows the light distribution but strongly deviates from the distribution of the intra-cluster gas as traced by the X-ray surface brightness. This confirms the complex morphology proposed by previous studies. We interpret this segregation of collisional and collisionless matter as strong evidence of multiple mergers and ongoing dynamical activity. MACS J0717 thus constitutes one of the most disturbed clusters presently known and, featuring a projected mass within the ACS field of view (R=150"=960 kpc) of 2.11 +- 0.23 x 10^{15} M_{sun}, the system is also one of the most massive known.

We report new constraints on the galaxy luminosity function at z~9 based on observations carried out with ESO/VLT FORS2, HAWK-I and X-Shooter around the lensing cluster A2667, as part of our project aimed at selecting z~7-10 candidates accessible to spectroscopy. Only one J-dropout source was selected in this field fulfilling the color and magnitude criteria. This source was recently confirmed as a mid-z interloper based on X-Shooter spectroscopy. The depth and the area covered by our survey are well suited to set strong constraints on the bright-end of the galaxy luminosity function and hence on the star formation history at very high redshift. The non-detection of reliable J-dropout sources over the ~36arcmin2 field of view towards A2667 was used to carefully determine the lens-corrected effective volume and the corresponding upper-limit on the density of sources. The strongest limit is obtained for Phi(M_{1500}=-21.4+/-0.50)<6.70x10^{-6}Mpc^{-3}mag^{-1} at z~9. A maximum-likelihood fit of the luminosity function using all available data points including the present new result yields M*>-19.7 with fixed alpha=-1.74 and Phi*=1.10x10^{-3}Mpc^{-3}. The corresponding star formation rate density should be rho_{SFR}<5.97x10^{-3}M_{solar}/yr/Mpc^{3} at z~9. These results are in good agreement with the most recent estimates already published in this range of redshift and for this luminosity domain. This new result confirms the decrease in the density of luminous galaxies at very high-redshift, hence providing strong constraints for the design of future surveys aiming to explore the very high-redshift Universe.

Measuring the intrinsic shape and orientation of dark matter (DM) and intracluster (IC) gas in galaxy clusters is crucial to constraining their formation and evolution, and for enhancing the use of clusters as more precise cosmological probes. Extending our previous works, we present for the first time results from a triaxial joint analysis of the galaxy cluster Abell 1835, by means of X-ray, strong lensing (SL) and Sunyaev Zel'dovich (SZ) data. We parametrically reconstruct the full three-dimensional structure (triaxial shape and principal axis orientation) of both the DM and the IC gas, and the level of non-thermal pressure of the IC gas. We find that the intermediate-major and minor-major axis ratios of the DM are 0.71+/-0.08 and 0.59+/-0.05, respectively, and the major axis of the DM halo is inclined with respect to the line of sight at 18.3+/-5.2 deg. We present the first observational measurement of the non-thermal pressure out to R_{200}, which has been evaluated to be a few percent of the total energy budget in the internal regions, while reaching approximately 20% in the outer volumes. We discuss the implications of our method for the viability of the CDM scenario, focusing on the concentration parameter C and the inner slope of the DM gamma in order to test the cold dark matter (CDM) paradigm for structure formation: we measure gamma=1.01+/-0.06 and C=4.32+/-0.44, values which are close to the predictions of the CDM model. The combination of X-ray/SL data at high spatial resolution, capable of resolving the cluster core, with the SZ data, which are more sensitive to the cluster outer volume, allows us to characterize the level and the gradient of the gas entropy distribution and non-thermal pressure out to R_{200}, breaking the degeneracy among the physical models describing the thermal history of the ICM.

We searched for z=7.3 Lya emitters (LAEs) behind two lensing clusters, Abell 2390 and CL 0024, with the Subaru Telescope Suprime-Cam and a narrowband NB1006 (FWHM ~ 21 nm centered at 1005 nm). We investigated if there exist objects consistent with the color of z=7.3 LAEs behind the clusters but could not detect any LAEs to the unlensed line limit F(Lya) ~ 6.9 x 10^{-18} erg/s/cm^2. Using several z=7 Lya luminosity functions (LFs) from the literature, we estimated and compared the expected detection numbers of z ~ 7 LAEs in lensing and blank field surveys in the case of using an 8m class ground based telescope. Given the steep bright-end slope of the LFs, when the detector field-of view (FOV) is comparable to the angular extent of a massive lensing cluster, imaging cluster(s) is more efficient in detecting z ~ 7 LAEs than imaging a blank field. However, the gain is expected to be modest, a factor of two at most and likely much less depending on the adopted LFs. The main advantage of lensing-cluster survey, therefore, remains to be the gain in depth and not necessarily in detection efficiency. For much larger detectors, the lensing effect becomes negligible and the efficiency of LAE detection is proportional to the instrumental FOV. We also inspected NB1006 images of three z ~ 7 z-dropouts previously detected in Abell 2390 and found that none of them are detected in NB1006. Two of them are consistent with predictions from the previous studies that they would be at lower redshifts. The other one has a photometric redshift of z ~ 7.3, and if it is at z=7.3, its unlensed Lya line flux would be very faint: F(Lya) < 4.4 x 10^{-18} erg/s/cm^2 (1 sigma upper limit) or rest frame equivalent width of W(Lya) < 26A. Its Lya emission might be attenuated by neutral hydrogen, as recent studies show that the fraction of Lyman break galaxies displaying strong Lya emission is lower at z ~ 7 than at z <~ 6.

We present secure [NII[ detections in two mm-bright, strongly lensed objects at high redshift, APM08279+5255 (z=3.911) and MM18423+5938 (z=3.930), using the IRAM Plateau de Bure Interferometer. Due to its ionization energy [NII] is a good tracer of the ionized gas phase in the interstellar medium. The measured fluxes are S([NII])=(4.8+/-0.8) Jy km/s and (7.4+/-0.5) Jy km/s respectively, yielding line luminosities of L([NII]) =(1.8+/-0.3) x 10^9 \mu^{-1} Lsun for APM08279+5255 and L([NII]) =(2.8+/-0.2) x 10^9 \mu^{-1} Lsun for MM18423+5938. Our high-resolution map of the [NII] and 1 mm continuum emission in MM18423+5938 clearly resolves an Einstein ring in this source, and reveals a velocity gradient in the dynamics of the ionized gas. A comparison of these maps with high-resolution EVLA CO observations enables us to perform the first spatially-resolved study of the dust continuum-to-molecular gas surface brightness (Sigma_{FIR} Sigma_{CO}^N, which can be interpreted as the star formation law) in a high-redshift object. We find a steep relation (N=1.4+/-0.2), consistent with a starbursting environment. We measure a [NII]/FIR luminosity ratio in APM0828+5255 and MM18423+5938 of 9.0 x 10^{-6} and 5.8 x 10^{-6}, respectively. This is in agreement with the decrease of the [NII]/FIR ratio at high FIR luminosities observed in local galaxies.

We present the analysis of new NIR, intermediate-resolution spectra of the gravitationally lensed galaxy "the 8 o'clock arc" at z_sys = 2.7350 obtained with VLT/X-shooter. These rest-frame optical data, combined with HST and Spitzer images, provide very valuable information, which nicely complement our previous detailed rest-frame UV spectral analysis. From high-resolution HST images, we reconstruct the morphology of the arc in the source plane, and identify that the source is formed of two majors parts, the main galaxy component and a smaller blob separated by 1.2 kpc in projected distance. The blob, with a twice larger magnification factor, is resolved in the spectra. The multi-Gaussian fitting of detected nebular emission lines and the spectral energy distribution modeling of the available multi-wavelength photometry provide the census of gaseous and stellar dust extinctions, gas-phase metallicities, star-formation rates (SFRs), and stellar, gas, and dynamical masses for both the main galaxy and the blob. As a result, the 8 o'clock arc shows a marginal trend for a more attenuated ionized gas than stars, and supports a dependence of the dust properties on the SFR. With a high specific star-formation rate, SSFR = 33+/-19 Gyr^{-1}, this lensed Lyman-break galaxy deviates from the mass-SFR relation, and is characterized by a young age of 40^{+25}_{-20} Myr and a high gas fraction of about 72%. The 8 o'clock arc satisfies the fundamental mass, SFR, and metallicity relation, and favors that it holds up beyond z~2.5. We believe that the blob, with a gas mass M_gas = (2.2+/-0.9)x10^{9} Msun (one order of magnitude lower than the mass of the galaxy), a half-light radius r1/2 = 0.53+/-0.05 kpc, a star-formation rate SFR_Halpha = 33+/-19 Msun yr^{-1}, and in rotation around the main core of the galaxy, is one of these star-forming clumps commonly observed in z>1 star-forming galaxies. (Abridged)

During our Herschel Lensing Survey (HLS) of massive galaxy clusters, we have discovered an exceptionally bright source behind the z=0.22 cluster Abell 773, which appears to be a strongly lensed submillimeter galaxy (SMG) at z=5.2429. This source is unusual compared to most other lensed sources discovered by Herschel so far, because of its higher submm flux (\sim 200mJy at 500\micron) and its high redshift. The dominant lens is a foreground z=0.63 galaxy, not the cluster itself. The source has a far-infrared (FIR) luminosity of L_FIR= 1.1 10^{14}/\mu Lo, where \mu is the magnification factor, likely \sim 11. We report here the redshift identification through CO lines with the IRAM-30m, and the analysis of the gas excitation, based on CO(7-6), CO(6-5), CO(5-4) detected at IRAM and the CO(2-1) at the EVLA. All lines decompose into a wide and strong red component, and a narrower and weaker blue component, 540\kms apart. Assuming the ultraluminous galaxy (ULIRG) CO-to-H2 conversion ratio, the H2 mass is 5.8 10^{11}/\mu Mo, of which one third is in a cool component. From the CI line we derive a CI/H2 number abundance of 6 10^{-5} similar to that in other ULIRGs. The H2O line is strong only in the red velocity component, with an intensity ratio I(H_2O)/I(CO) \sim 0.5, suggesting a strong local FIR radiation field, possibly from an active nucleus (AGN) component. We detect the [NII]205\mics line for the first time at high-z. It shows comparable blue and red components, with a strikingly broad blue one, suggesting strong ionized gas flows.

Our ability to study the most remote supernova explosions, crucial for the understanding of the evolution of the high-redshift universe and its expansion rate, is limited by the light collection capabilities of telescopes. However, nature offers unique opportunities to look beyond the range within reach of our unaided instruments thanks to the light-focusing power of massive galaxy clusters. Here we report on the discovery of one of the most distant supernovae ever found, at redshift, z=1.703. Due to a lensing magnification factor of 4.3\pm0.3, we are able to measure a lightcurve of the supernova, as well as spectroscopic features of the host galaxy with a precision comparable to what will otherwise only be possible with future generation telescopes.

We have used the IRAM Plateau de Bure Interferometer and the Expanded Very Large Array to obtain a high resolution map of the CO(6-5) and CO(1-0) emission in the lensed, star-forming galaxy SMMJ2135-0102 at z=2.32. The kinematics of the gas are well described by a model of a rotationally-supported disk with an inclination-corrected rotation speed, v_rot = 320+/-25km/s, a ratio of rotational- to dispersion- support of v/sigma=3.5+/-0.2 and a dynamical mass of 6.0+/-0.5x10^10Mo within a radius of 2.5kpc. The disk has a Toomre parameter, Q=0.50+/-0.15, suggesting the gas will rapidly fragment into massive clumps on scales of L_J ~ 400pc. We identify star-forming regions on these scales and show that they are 10x denser than those in quiescent environments in local galaxies, and significantly offset from the local molecular cloud scaling relations (Larson's relations). The large offset compared to local molecular cloud linewidth-size scaling relations imply that supersonic turbulence should remain dominant on scales ~100x smaller than in the kinematically quiescent ISM of the Milky Way, while the molecular gas in SMMJ2135 is expected to be ~50x denser than that in the Milky Way on all scales. This is most likely due to the high external hydrostatic pressure we measure for the interstellar medium (ISM), P_tot/kB ~ (2+/-1)x10^7K/cm3. In such highly turbulent ISM, the subsonic regions of gravitational collapse (and star-formation) will be characterised by much higher critical densities, n_crit>=10^8/cm3, a factor ~1000x more than the quiescent ISM of the Milky Way.

Using deep Keck spectroscopy of Lyman break galaxies selected from infrared imaging data taken with WFC3/IR onboard the Hubble Space Telescope, we present new evidence for a reversal in the redshift-dependent fraction of star forming galaxies with detectable Lyman alpha emission in the redshift range 6.3 < z < 8.8. Our earlier surveys with the DEIMOS spectrograph demonstrated a significant increase with redshift in the fraction of line emitting galaxies over the interval 4 < z < 6, particularly for intrinsically faint systems which dominate the luminosity density. Using the longer wavelength sensitivities of LRIS and NIRSPEC, we have targeted 19 Lyman break galaxies selected using recent WFC3/IR data whose photometric redshifts are in the range 6.3 < z < 8.8 and which span a wide range of intrinsic luminosities. Our spectroscopic exposures typically reach a 5-sigma sensitivity of < 50 A for the rest-frame equivalent width (EW) of Lyman alpha emission. Despite the high fraction of emitters seen only a few hundred million years later, we find only 2 convincing and 1 possible line emitter in our more distant sample. Combining with published data on a further 7 sources obtained using FORS2 on the ESO VLT, and assuming continuity in the trends found at lower redshift, we discuss the significance of this apparent reversal in the redshift-dependent Lyman alpha fraction in the context of our range in continuum luminosity. Assuming all the targeted sources are at their photometric redshift and our assumptions about the Lyman alpha EW distribution are correct, we would expect to find so few emitters in less than 1% of the realizations drawn from our lower redshift samples. Our new results provide further support for the suggestion that, at the redshifts now being probed spectroscopically, we are entering the era where the intergalactic medium is partially neutral.

Massive galaxy clusters at intermediate redshifts act as gravitational lenses that can magnify supernovae (SNe) occurring in background galaxies. We assess the possibility to use lensed SNe to put constraints on the mass models of galaxy clusters and the Hubble parameter at high redshift. Due to the standard candle nature of Type Ia supernovae (SNe Ia), observational information on the lensing magnification from an intervening galaxy cluster can be used to constrain the model for the cluster mass distribution. A statistical analysis using parametric cluster models was performed to investigate the possible improvements from lensed SNe Ia for the accurately modeled galaxy cluster A1689 and the less well constrained cluster A2204. Time delay measurements obtained from SNe lensed by accurately modeled galaxy clusters can be used to measure the Hubble parameter. For a survey of A1689 we estimate the expected rate of detectable SNe Ia and of multiply imaged SNe. The velocity dispersion and core radius of the main cluster potential show strong correlations with the predicted magnifications and can therefore be constrained by observations of SNe Ia in background galaxies. This technique proves especially powerful for galaxy clusters with only few known multiple image systems. The main uncertainty for measurements of the Hubble parameter from the time delay of strongly lensed SNe is due to cluster model uncertainties. For the extremely well modeled cluster A1689, a single time delay measurement could be used to determine the Hubble parameter with a precision of ~ 10%. We conclude that observations of SNe Ia behind galaxy clusters can be used to improve the mass modeling of the large scale component of galaxy clusters and thus the distribution of dark matter. Time delays from SNe strongly lensed by accurately modeled galaxy clusters can be used to measure the Hubble constant at high redshifts.

The spectral energy distributions (SED) of dusty galaxies at intermediate redshift may look similar to very high redshift galaxies in the optical/near infrared (NIR) domain. This can lead to the contamination of high redshift galaxy searches based on broad band optical/NIR photometry by lower redshift dusty galaxies as both kind of galaxies cannot be distinguished. The contamination rate could be as high as 50%. {This work shows how the far infrared (FIR) domain can help to recognize likely low-z interlopers in an optical/NIR search for high-z galaxies.} We analyse the FIR SEDs of two galaxies proposed as very high redshift ($z>7$) dropout candidates based on deep Hawk-I/VLT observations. The FIR SEDs are sampled with PACS/Herschel at 100 and 160\,$\mu$m, with SPIRE/Herschel at 250, 350 and 500\,$\mu$m and with LABOCA/APEX at 870\,$\mu$m. We find that redshifts $>7$ would imply extreme FIR SEDs (with dust temperatures $>100$\,K and FIR luminosities $>10^{13}$\,$L_{\odot}$). At z$\sim$2, instead, the SEDs of both sources would be compatible with that of typical ULIRGs/SMGs. Considering all the data available for these sources from visible to FIR we re-estimate the redshifts and we find $z\sim$1.6--2.5. Due to the strong spectral breaks observed in these galaxies, standard templates from the literature fail to reproduce the visible-near IR part of the SEDs even when additional extinction is included. These sources resemble strongly dust obscured galaxies selected in Spitzer observations with extreme visible-to-FIR colors, and the galaxy GN10 at $z=4$. Galaxies with similar SEDs could contaminate other high redshift surveys.

We investigate the nature and the physical properties of z, Y and J-dropout galaxies selected behind the lensing cluster A2667. This field is part of our project aimed at identifying z~7-10 candidates accessible to spectroscopic studies, based on deep photometry with ESO/VLT HAWK-I and FORS2 (zYJH and Ks-band images, AB(3 sigma)~26-27) on a sample of lensing clusters extracted from our multi-wavelength combined surveys with SPITZER, HST, and Herschel. In this paper we focus on the complete Y and J-dropout sample, as well as the bright z-dropouts fulfilling the selection criteria by Capak et al. (2011). 10 candidates are selected within the common field of ~33 arcmin2 (effective area once corrected for contamination and lensing dilution). All of them are detected in H and Ks bands in addition to J and/or IRAC 3.6/4.5, with H(AB)~23.4 to 25.2, and have modest magnification factors. Although best-fit photometric redshifts place all these candidates at high-z, the contamination by low-z interlopers is estimated at 50-75% level based on previous studies, and the comparison with the blank-field WIRCAM Ultra-Deep Survey (WUDS). The same result is obtained when photometric redshifts include a luminosity prior, allowing us to remove half of the original sample as likely z~1.7-3 interlopers with young stellar pulations and strong extinction. Two additional sources among the remaining sample could be identified at low-z based on a detection at 24 microns and on the HST z_850 band. These low-z interlopers are not well described by current templates given the large break, and cannot be easily identified based solely on optical and near-IR photometry. Given the estimated dust extinction and high SFRs, some of them could be also detected in the IR or sub-mm bands. After correction for likely contaminants, the observed counts at z>7.5 seem to be in agreement with an evolving LF. (abridged)

We report the discovery of a unique $z=6.027$ galaxy, multiply imaged by the cluster Abell 383 and detected in new Hubble Space Telescope ACS and WFC3 imaging, as well as in Warm Spitzer observations. This galaxy was selected as a pair of i-dropouts; its suspected high redshift was confirmed by the measurement of a strong Lyman-alpha line in both images using Keck/DEIMOS. Combining Hubble and Spitzer photometry after correcting for contamination by line emission (estimated to be a small effect), we identify a strong Balmer break of 1.5 magnitudes. Taking into account the magnification factor of 11.4+/-1.9 (2.65+/-0.17 mag) for the brightest image, the unlensed AB magnitude for the source is 27.2+/-0.05 in the H band, corresponding to a 0.4 L* galaxy, and 25.7+/-0.08 at 3.6 um. The UV slope is consistent with beta~2.0, and from the rest-frame UV continuum we measure a current star formation rate of 2.4+/-1.1 Msol/yr. The unlensed half-light radius is measured to be 300 pc, from which we deduce a star-forming surface density of ~10 Msol/yr/kpc2. The Lyman-alpha emission is found to be extended over ~3" along the slit, corresponding to ~5 kpc in the source plane. This can be explained by the presence of a much larger envelope of neutral hydrogen around the star-forming region. Finally, fitting the spectral energy distribution using 7 photometric data points with simple SED models, we derive the following properties: very little reddening, an inferred stellar mass of M*=6e9 Msol, and an inferred age of ~800 Myrs (corresponding to a redshift of formation of ~18). The star-formation rate of this object was likely much stronger in the past than at the time of observation, suggesting that we may be missing a fraction of galaxies at z~6 which have already faded in rest-frame UV wavelengths.

We present the first metallicity gradient measurement for a grand-design face-on spiral galaxy at z~1.5. This galaxy has been magnified by a factor of 22$\times$ by a massive, X-ray luminous galaxy cluster MACS\,J1149.5+2223 at z=0.544. Using the Laser Guide Star Adaptive Optics aided integral field spectrograph OSIRIS on KECK II, we target the Halpha emission and achieve a spatial resolution of 0.1", corresponding to a source plane resolution of 170 pc. The galaxy has well-developed spiral arms and the nebular emission line dynamics clearly indicate a rotationally supported disk with V_{rot}/\sigma~4. The best-fit disk velocity field model yields a maximum rotation of V_{rot} sin{i}=150$\pm$15 km s^{-1}, and a dynamical mass of M_{dyn}=1.3$\pm0.2\times10^{10}csc^2(i) M_{\odot} (within 2.5\,kpc), where the inclination angle i=45$\pm10^{\circ}$. Based on the [NII] and Halpha ratios, we measured the radial chemical abundance gradient from the inner hundreds of parsecs out to ~5 kpc. The slope of the gradient is -0.16$\pm$0.02 dex kpc$^{-1}$, significantly steeper than the gradient of late-type or early-type galaxies in the local universe. If representative of disk galaxies at z~1.5, our results support an "inside-out" disk formation scenario in which early infall/collapse in the galaxy center builds a chemically enriched nucleus, followed by slow enrichment of the disk over the next 9 Gyr.

We present the results of our analyses of the X-ray emission and of the strong lensing systems in the relaxed galaxy cluster Abell 611 (z=0.288). We infer the X-ray mass estimate deriving the density and temperature profiles of the intra-cluster medium within the radius r ~ 700 kpc through a non-parametric approach; assuming that the cluster is in hydrostatic equilibrium and adopting a matter density profile, we can recover the total mass distribution of Abell 611 via the X-ray data. We derive the total projected mass in the central region of Abell 611 performing a parametric analysis of its strong lensing features through the publicly available analysis software Lenstool. As a final step we compare the results obtained with both methods. We derive a good agreement between the X-ray and strong lensing total mass estimates in the central region (i.e. within the radius r ~100 kpc), while a marginal disagreement is found between the two mass estimates when extrapolating the strong lensing results in the outer spatial range. We suggest that in this case the X-ray/strong lensing mass disagreement can be explained by an incorrect estimate of the relative contributions of the baryonic component and of the dark matter, caused by the intrinsic degeneracy between the different mass components in the strong lensing analysis. We discuss the effect of some possible systematic errors that influence both mass estimates. We find a slight dependence of the measurements of the X-ray temperatures (and therefore of the X-ray total masses) on the background adopted in the spectral analysis, with total deviations on the value of M_200 of the order of the 1-sigma statistical error. The strong lensing mass results are instead sensitive to the parameterisation of the galactic halo mass in the central regions, in particular to the modelling of the Brightest Cluster Galaxy (BCG) baryonic component.

We aim to probe the mass of SL2S\,J02140-0535, a galaxy group at $z$ = 0.44 from the Strong Lensing Legacy Survey (SL2S). We combine strong lensing modeling and dynamical constraints. The strong lensing analysis is based on multi-band HST/ACS observations exhibiting strong lensing features that we have followed-up spectroscopically with VLT/FORS2. To constrain the scale radius of an NFW mass profile that cannot be constrained by strong lensing, we propose a new method by taking advantage of the large-scale dynamical information provided by VLT/FORS2 and KECK/LRIS spectroscopy of group members. In constrast to other authors, we show that the observed lensing features in SL2S\,J02140-0535 belong to different background sources: one at $z$ = 1.7 $\pm$ 0.1 produces three images, while the other at $z$ = 1.023 $\pm$ 0.001 has only a single image. Our unimodal NFW mass model reproduces these images very well. It is characterized by a concentration parameter $c_{200}$ = 6.0 $\pm$ 0.6, which is slightly greater than the value expected from $\Lambda$CDM simulations for a mass of M$_{200}$ $\approx$ 1 $\times$ 10$^{14}$ M$_{\sun}$. The spectroscopic analysis of group members also reveals a unimodal structure that exhibits no evidence of merging. We compare our dynamic mass estimate with an independent weak-lensing based mass estimate finding that both are consistent. Our combined lensing and dynamical analysis of SL2S\,J02140-0535 demonstrates the importance of spectroscopic information in reliably identifying the lensing features. Our findings argue that the system is a relaxed, massive galaxy group where mass is traced by light. This work shows a potentially useful method for constraining large-scale properties inaccessible to strong lensing, such as the scale radius of the NFW profile.

We present and analyse near-infrared spectroscopy for a sample of 28 gravitationally- lensed star-forming galaxies in the redshift range 1.5 < z < 5, observed mostly with the Keck II telescope. With typical magnifications of ~1.5-4 magnitudes, our survey provides a valuable census of star formation rates, gas-phase metallicities and dynamical masses for a representative sample of low luminosity galaxies seen at a formative period in cosmic history. We find less evolution in the mass-metallicity relation compared to earlier work that focused on more luminous systems with z - 2-3, especially in the low mass (- 10^9 Msol) where our sample is - 0.25 dex more metal-rich. We interpret this offset as a result of the lower star formation rates (typically a factor of -10 lower) for a given stellar mass in our sub-luminous systems. Taking this effect into account, we conclude our objects are consistent with a fundamental metallicity relation recently proposed from unlensed observations.

We present near-infrared imaging spectroscopy of the strongly-lensed z=2.00 galaxy SDSS J120601.69+514227.8 (`the Clone arc'). Using OSIRIS on the Keck 2 telescope with laser guide star adaptive optics, we achieve resolved spectroscopy with 0.20 arcsecond FWHM resolution in the diagnostic emission lines [O III], Halpha, and [N II]. The lensing magnification allows us to map the velocity and star formation from Halpha emission at a physical resolution of ~300 pc in the galaxy source plane. With an integrated star formation rate of ~50 Msun/yr, the galaxy is typical of sources similarly studied at this epoch. It is dispersion-dominated with a velocity gradient of +/- 80 km/s and average dispersion sigma = 85 km/s; the dynamical mass is 2.4 \times 10^{10} Msun within a half-light radius of 2.9 kpc. Robust detection of [N II] emission across the entire OSIRIS field of view enables us to trace the gas-phase metallicity distribution with 500 pc resolution. We find a strong radial gradient in both the [N II]/Halpha and [O III]/Halpha ratios indicating a metallicity gradient of -0.27 +/- 0.05 dex/kpc with central metallicity close to solar. We demonstrate that the gradient is seen independently in two multiple images. While the physical gradient is considerably steeper than that observed in local galaxies, in terms of the effective radius at that epoch, the gradient is similar. This suggests that subsequent growth occurs in an inside-out manner with the inner metallicity gradient diminished over time due to radial mixing and enrichment from star formation.

Precise modelling of strong lensing systems can be affected by external mass distributions, e.g. the group or cluster within which the lens is embedded. In this article, we propose to turn this limitation to our advantage and to use precise strong lensing modelling to probe external mass distributions surrounding the lens. We consider SL2S J08544-0121, a galaxy group at z=0.35 that contains a strong lensing system. A simple elliptical isothermal potential cannot reproduce satisfactorily the strong lensing constraints. We include an external mass perturbation corresponding to the group within which the lens is embedded. The lensing properties of this perturbation are parametrised by its total mass M and a smoothing scale s that quantifies the characteristic scale over which M is distributed. For a range of these parameters, we are able to reproduce accurately the observations. This suggests that light is a good tracer of mass. Interestingly, this also shows that a localised strong lensing analysis (on scales of ~10") allows us to constrain global properties of the group as a whole (on scales of ~100). Indeed, we constrain the group mass-to-light ratio to be M/L=98+-27 (i band, solar units, not corrected for evolution) and s=20" +- 9 (2sigma confidence level). We demonstrate that these strong lensing only constraints are due to the perturbed strong lensing configuration, where the main arc is located at ~5" from the galaxy, whereas its counter-image is found at ~8". To test independently our resulting strong lensing model, we pursue an independent weak lensing analysis of the group and find a mass-to-light ratio in the range 66-146 (1sigma confidence level).

Current efforts in observational cosmology are focused on characterizing the mass-energy content of the Universe. We present results from a geometric test based on strong lensing in galaxy clusters. Based on Hubble Space Telescope images and extensive ground-based spectroscopic follow-up of the massive galaxy cluster Abell 1689, we used a parametric model to simultaneously constrain the cluster mass distribution and dark energy equation of state. Combining our cosmological constraints with those from X-ray clusters and the Wilkinson Microwave Anisotropy Probe 5-year data gives {\Omega}m = 0.25 +/- 0.05 and wx = -0.97 +/- 0.07 which are consistent with results from other methods. Inclusion of our method with all other techniques available brings down the current 2{\sigma} contours on the dark energy equation of state parameter wx by about 30%.

Gravitational lensing by massive galaxy clusters allows study of the population of intrinsically faint infrared galaxies that lie below the sensitivity and confusion limits of current infrared and submillimeter telescopes. We present ultra-deep PACS 100 and 160 microns observations toward the cluster lens Abell 2218, to penetrate the Herschel confusion limit. We derive source counts down to a flux density of 1 mJy at 100 microns and 2 mJy at 160 microns, aided by strong gravitational lensing. At these levels, source densities are 20 and 10 beams/source in the two bands, approaching source density confusion at 160 microns. The slope of the counts below the turnover of the Euclidean-normalized differential curve is constrained in both bands and is consistent with most of the recent backwards evolutionary models. By integrating number counts over the flux range accessed by Abell 2218 lensing (0.94-35 mJy at 100 microns and 1.47-35 mJy at 160 microns, we retrieve a cosmic infrared background (CIB) surface brightness of ~8.0 and ~9.9 nW m^-2 sr^-1, in the respective bands. These values correspond to 55% (+/- 24%) and 77% (+/- 31%) of DIRBE direct measurements. Combining Abell 2218 results with wider/shallower fields, these figures increase to 62% (+/- 25%) and 88% (+/- 32%) CIB total fractions, resolved at 100 and 160 microns, disregarding the high uncertainties of DIRBE absolute values.

The Herschel Lensing Survey (HLS) will conduct deep PACS and SPIRE imaging of ~40 massive clusters of galaxies. The strong gravitational lensing power of these clusters will enable us to penetrate through the confusion noise, which sets the ultimate limit on our ability to probe the Universe with Herschel. Here, we present an overview of our survey and a summary of the major results from our Science Demonstration Phase (SDP) observations of the Bullet Cluster (z=0.297). The SDP data are rich, allowing us to study not only the background high-redshift galaxies (e.g., strongly lensed and distorted galaxies at z=2.8 and 3.2) but also the properties of cluster-member galaxies. Our preliminary analysis shows a great diversity of far-infrared/submillimeter spectral energy distributions (SEDs), indicating that we have much to learn with Herschel about the properties of galaxy SEDs. We have also detected the Sunyaev-Zel'dovich (SZ) effect increment with the SPIRE data. The success of this SDP program demonstrates the great potential of the Herschel Lensing Survey to produce exciting results in a variety of science areas.

The Sunyaev-Zel'dovich (SZ) effect is a spectral distortion of the cosmic microwave background as observed through the hot plasma in galaxy clusters. This distortion is a decrement in the CMB intensity for lambda > 1.3 mm, an increment at shorter wavelengths, and small again by lambda ~250 um. As part of the Herschel Lensing Survey (HLS) we have mapped 1E0657-56 (the Bullet cluster) with SPIRE with bands centered at 250, 350 and 500 um and have detected the SZ effect at the two longest wavelengths. The measured SZ effect increment central intensities are Delta I_{0} = 0.097 +- 0.019 MJy sr^{-1} at 350 um and Delta I_{0} = 0.268 +- 0.031 MJy sr^{-1} at 500 um, consistent with the SZ effect spectrum derived from previous measurements at 2 mm. No other diffuse emission is detected. The presence of the finite temperature SZ effect correction is preferred by the SPIRE data at a significance of 2.1 sigma, opening the possibility that the relativistic SZ effect correction can be constrained by SPIRE in a sample of clusters. The results presented here have important ramifications for both sub-mm measurements of galaxy clusters and blank field surveys with SPIRE.

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 um sources as priors. Our catalogs are able to reliably (5 sigma) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 um channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.

We use deep, five band (100-500um) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFR_FIR, of galaxies in the Bullet cluster (z=0.296), and a smaller background system (z=0.35) in the same field. Herschel detects 23 Bullet cluster members with a total SFR_FIR = 144 +/- 14 M_sun yr^-1. On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher SFR_FIR (21 galaxies; 207 +/- 9 M_sun yr^-1). SFRs extrapolated from 24um flux via recent templates (SFR_24) agree well with SFR_FIR for ~60% of the cluster galaxies. In the remaining ~40%, SFR_24 underestimates SFR_FIR due to a significant excess in observed S_100/S_24 (rest frame S_75/S_18) compared to templates of the same FIR luminosity.

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 micron bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 micron and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 micron selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities.This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

We review progress in understanding dark matter by astrophysics, and particularly via the effect of gravitational lensing. Evidence from many different directions now all imply that five sixths of the material content of the universe is in this mysterious form, separate from and beyond the ordinary "baryonic" particles in the standard model of particle physics. Dark matter appears not to interact via the electromagnetic force, and therefore neither emits nor reflects light. However, it definitely does interact via gravity, and has played the most important role in shaping the Universe on large scales. The most successful technique with which to investigate it has so far been the effect of gravitational lensing. The curvature of space-time near any gravitating mass (including dark matter) deflects passing rays of light - observably shifting, distorting and magnifying the images of background galaxies. Measurements of such effects currently provide constraints on the mean density of dark matter, and its density relative to baryonic matter; the size and mass of individual dark matter particles; and its cross section under various fundamental forces.

Massive galaxies in the early Universe have been shown to be forming stars at surprisingly high rates. Prominent examples are dust-obscured galaxies which are luminous when observed at sub-millimeter (sub-mm) wavelengths and which may be forming stars at rates upto 1,000Mo/yr. These intense bursts of star formation are believed to be driven by mergers between gas rich galaxies. However, probing the properties of individual star-forming regions within these galaxies is beyond the spatial resolution and sensitivity of even the largest telescopes at present. Here, we report observations of the sub-mm galaxy SMMJ2135-0102 at redshift z=2.3259 which has been gravitationally magnified by a factor of 32 by a massive foreground galaxy cluster lens. This cosmic magnification, when combined with high-resolution sub-mm imaging, resolves the star-forming regions at a linear scale of just ~100 parsecs. We find that the luminosity densities of these star-forming regions are comparable to the dense cores of giant molecular clouds in the local Universe, but they are ~100x larger and 10^7 times more luminous. Although vigorously star-forming, the underlying physics of the star formation processes at z~2 appears to be similar to that seen in local galaxies even though the energetics are unlike anything found in the present-day Universe.

We present results of a gravitational-lensing and optical study of MACS ,J1423.8+2404 (z=0.545, MACS, J1423). Our analysis uses high-resolution images taken with the Hubble Space Telescope in the F555W and F814W passbands, ground based imaging in eight optical and near-infrared filters obtained with Subaru and CFHT, as well as extensive spectroscopic data gathered with the Keck telescopes. At optical wavelengths the cluster exhibits no sign of substructure and is dominated by a cD galaxy that is 2.1 magnitudes (K-band) brighter than the second brightest cluster member, suggesting that MACS, J1423 is close to be fully virialized. Analysis of the redshift distribution of 140 cluster members reveals a Gaussian distribution, mildly disturbed by the presence of a loose galaxy group that may be falling into the cluster along the line of sight. Combining strong-lensing constraints from two spectroscopically confirmed multiple-image systems near the cluster core with a weak-lensing measurement of the gravitational shear on larger scales, we derive a parametric mass model for the mass distribution. All constraints can be satisfied by a uni-modal mass distribution centred on the cD galaxy and exhibiting very little substructure. The derived projected mass of M(<65\arcsec [415 kpc])=(4.3\pm0.6)\times 10^{14} M_sun is about 30% higher than the one derived from X-ray analyses assuming spherical symmetry, suggesting a slightly prolate mass distribution consistent with the optical indication of residual line-of-sight structure. The similarity in shape and excellent alignment of the centroids of the total mass, K-band light, and intra-cluster gas distributions add to the picture of a highly evolved system [ABRIDGED]

We present the discovery of a large-scale structure of emission-line galaxies at redshift z=4.86 behind a massive cluster of galaxies, A1689. Previous spectroscopic observations of a galaxy, A1689-7.1 at z=4.87, near this structure, revealed a possible overdense region of inter-galactic medium (IGM) around the galaxy, which extends at least ~80 comoving Mpc along the line of sight. In order to investigate whether this z~5 IGM overdense region contains a galaxy overdensity, we undertook narrow- and broad-band imaging observations around A1689-7.1 with Subaru/Suprime-Cam. We detected 51 candidate Ly-alpha emitters at redshift z=4.86+-0.03 in the 32 x 24 arcmin^2 field of view. After correction for lensing by the foreground cluster, we found a large-scale (~20 x 60 comoving Mpc) overdense region of galaxies around A1689-7.1 in the source plane at z=4.86. The densest peak in this region has an overdensity of delta~4, suggesting that this structure is probably a good candidate for a protocluster. A1689-7.1 is located at the edge of this region, where the local galaxy density is ~1.6 times the mean density and is close to the density contrast in the IGM along the line of sight to A1689-7.1 estimated from the optical depth. The overdensities of galaxies we have found may suggest that at least some parts of the IGM overdense region have already started to form galaxies and moreover they relate to the formation of a protocluster. Although we lack information on the three dimensional distributions of both IGM and galaxy overdense regions, the similarity of the scales of both regions may suggest that the two are parts of a single large-scale structure, which would be an large edge-on sheet along the line of sight with a size of ~20 x 60 x 80 comoving Mpc.

We present new multi-wavelength observations of the first submm-selected galaxy, SMM J02399-0136 at z=2.8. These observations include mapping of the CO J=1-0 emission using elements of the Expanded VLA, as well as high-resolution 1.4-GHz imaging and optical/IR data from the VLA, HST, Spitzer and Keck. Together these new data provide fundamental insights into the mass and distribution of stars, gas and dust within this archetypal SMG. The CO J=1-0 emission, with its minimal excitation and density requirements, traces the bulk of the metal-rich molecular gas, and reveals a mass of ~10^11 M(sun), extending over ~5" (~25 kpc in the source plane), although there is tentative evidence that it may be significantly larger. Our data suggest that three or more distinct structures are encompassed by this molecular gas reservoir, including the BAL quasar from which the redshift of the SMG was initially determined. In particular, the new rest-frame near-IR observations identify a massive, obscured, starburst coincident with a previously known Ly-alpha cloud. This starburst dominates the far-IR emission from the system and requires a re-assessment of previous claims that the gas reservoir resides in a massive, extended disk around the BAL QSO. Instead it appears that SMM J02399-0136 comprises a merger between a far-IR-luminous, but highly obscured starburst, the BAL QSO host and a faint third component. Our findings suggest that this SMG and its immediate environment mark a vast and complex galactic nursery and that detailed studies of other SMGs are likely to uncover a similarly rich diversity of properties.

We present a statistical analysis of a sample of 20 strong lensing clusters drawn from the Local Cluster Substructure Survey (LoCuSS), based on high resolution Hubble Space Telescope imaging of the cluster cores and follow-up spectroscopic observations using the Keck-I telescope. We use detailed parameterized models of the mass distribution in the cluster cores, to measure the total cluster mass and fraction of that mass associated with substructures within R<250kpc.These measurements are compared with the distribution of baryons in the cores, as traced by the old stellar populations and the X-ray emitting intracluster medium. Our main results include: (i) the distribution of Einstein radii is log-normal, with a peak and 1sigma width of <log(RE(z=2))>=1.16+/-0.28; (ii) we detect an X-ray/lensing mass discrepancy of <M_SL/M_X>=1.3 at 3 sigma significance -- clusters with larger substructure fractions displaying greater mass discrepancies, and thus greater departures from hydrostatic equilibrium; (iii) cluster substructure fraction is also correlated with the slope of the gas density profile on small scales, implying a connection between cluster-cluster mergers and gas cooling. Overall our results are consistent with the view that cluster-cluster mergers play a prominent role in shaping the properties of cluster cores, in particular causing departures from hydrostatic equilibrium, and possibly disturbing cool cores. Our results do not support recent claims that large Einstein radius clusters present a challenge to the CDM paradigm.

We present a strong lensing analysis of the galaxy cluster Abell 370 (z=0.375) based on the recent multicolor ACS images obtained as part of the Early Release Observation (ERO) that followed the Hubble Service Mission #4. Back in 1987, the giant gravitational arc (z=0.725) in Abell 370 was one of the first pieces of evidence that massive clusters are dense enough to act as strong gravitational lenses. The new observations reveal in detail its disklike morphology, and we show that it can be interpreted as a complex five-image configuration, with a total magnification factor of 32+/-4. Moreover, the high resolution multicolor information allowed us to identify 10 multiply imaged background galaxies. We derive a mean Einstein radius of RE=39+/-2" for a source redshift at z=2, corresponding to a mass of M(<RE) = 2.82+/-0.15 1e14 Msol and M(<250 kpc)=3.8+/-0.2 1e14 Msol, in good agreement with Subaru weak-lensing measurements. The typical mass model error is smaller than 5%, a factor 3 of improvement compared to the previous lensing analysis. Abell 370 mass distribution is confirmed to be bi-modal with very small offset between the dark matter, the X-ray gas and the stellar mass. Combining this information with the velocity distribution reveals that Abell 370 is likely the merging of two equally massive clusters along the line of sight, explaining the very high mass density necessary to efficiently produce strong lensing. These new observations stress the importance of multicolor imaging for the identification of multiple images which is key to determining an accurate mass model. The very large Einstein radius makes Abell 370 one of the best clusters to search for high redshift galaxies through strong magnification in the central region.

We present the identification of a bright submm source, SMMJ163555.5+661300, detected in the lensing cluster Abell2218, for which we have accurately determined the position using observations from the Submillimeter Array (SMA). The identified optical counterpart has a spectroscopic redshift of z=4.044+-0.001 if we attribute the single emission line detected at lambda=6140AA to Lyman-alpha. This redshift identification is in good agreement with the optical/near-infrared photometric redshift as well as the submm flux ratio S_450/S_850~1.6, the radio-submm flux ratio S_1.4/S_850 < 0.004, and the 24um to 850um flux ratio S_24/S_850 < 0.005. Correcting for the gravitational lensing amplification of ~5.5, we find that the source has a far-infrared luminosity of 1.3x10^12 Lsun, which implies a star formation rate of 230 Msun/yr. This makes it the lowest-luminosity SMG known at z>4 to date. Previous CO(4-3) emission line obserations yielded a non-detection, for which we derived an upper limit of the CO line luminosity of L'_CO = 0.3x10^10 K km/s/pc^2, which is not inconsistent with the L'_CO - L_FIR relation for starburst galaxies. The best fit model to the optical and near-infrared photometry give a stellar population with an age of 1.4 Gyr and a stellar mass of 1.6x10^10 Msun. The optical morphology is compact and in the source plane the galaxy has an extent of ~6kpc x 3kpc with individual star forming knots of <500 pc in size. J163556 is not resolved in the SMA data and we place a strict upper limit on the size of the starburst region of 8kpc x 3kpc, which implies a lower limit on the star formation rate surface density of 12 Msun/yr/kpc^2. The redshift of J163556 extends the redshift distribution of faint, lensed SMGs, and we find no evidence that these have a different redshift distribution than bright SMGs.

We present Advanced Camera for Surveys observations of MACSJ1149.5+2223, an X-ray luminous galaxy cluster at z=0.544 discovered by the Massive Cluster Survey. The data reveal at least seven multiply-imaged galaxies, three of which we have confirmed spectroscopically. One of these is a spectacular face-on spiral galaxy at z=1.491, the four images of which are gravitationally magnified by ~8<mu<~23. We identify this as an L* (M_B=-20.7), disk-dominated (B/T<~0.5) galaxy, forming stars at ~6Msol/yr. We use a robust sample of multiply-imaged galaxies to constrain a parameterized model of the cluster mass distribution. In addition to the main cluster dark matter halo and the bright cluster galaxies, our best model includes three galaxy-group-sized halos. The relative probability of this model is P(N_halo=4)/P(N_halo<4)>=10^12 where N_halo is the number of cluster/group-scale halos. In terms of sheer number of merging cluster/group-scale components, this is the most complex strong-lensing cluster core studied to date. The total cluster mass and fraction of that mass associated with substructures within R<=500kpc, are measured to be M_tot=(6.7+/-0.4)x10^14Msol and f_sub=0.25+/-0.12 respectively. Our model also rules out recent claims of a flat density profile at >~7sigma confidence, thus highlighting the critical importance of spectroscopic redshifts of multiply-imaged galaxies when modeling strong lensing clusters. Overall our results attest to the efficiency of X-ray selection in finding the most powerful cluster lenses, including complicated merging systems.

The identification and study of the first galaxies remains one of the most exciting topics in observational cosmology. The determination of the best possible observing strategies is a very important choice in order to build up a representative sample of spectroscopically confirmed sources at high-z (z>7), beyond the limits of present-day observations. This paper is intended to precisely adress the relative efficiency of lensing and blank fields in the identification and study of galaxies at 6<z<12. The detection efficiency and field-to-field variance are estimated from direct simulations of both blank and lensing fields observations. The presence of a strong-lensing cluster along the line of sight has a dramatic effect on the number of observed sources, with a positive magnification bias in typical ground-based ``shallow'' surveys (AB<~25.5). The positive magnification bias increases with the redshift of sources and decreases with both depth of the survey and the size of the surveyed area. The maximum efficiency is reached for lensing clusters at z~0.1-0.3. Observing blank fields in shallow surveys is particularly inefficient as compared to lensing fields if the UV LF for LBGs is strongly evolving at z>~7. Also in this case, the number of z>8 sources expected at the typical depth of JWST (AB~28-29) is much higher in lensing than in blank fields (e.g. a factor of ~10 for AB<~28). Blank field surveys with a large field of view are needed to prove the bright end of the LF at z>6-7, whereas lensing clusters are particularly useful for exploring the mid to faint end of the LF.

We present spatially-resolved dynamics for six strongly lensed star-forming galaxies at z=1.7-3.1, each enlarged by a linear magnification factor ~8. Using the Keck laser guide star AO system and the OSIRIS integral field unit spectrograph we resolve kinematic and morphological detail in our sample with an unprecedented fidelity, in some cases achieving spatial resolutions of ~100 pc. With one exception our sources have diameters ranging from 1-7 kpc, star formation rates of 2-40 Msun/yr (uncorrected for extinction) and dynamical masses of 10^(9.7-10.3) Msun. With this exquisite resolution we find that four of the six galaxies display coherent velocity fields consistent with a simple rotating disk model, which can only be recovered with the considerably improved spatial resolution and sampling from the combination of adaptive optics and strong gravitational lensing. Our model fits imply ratios for the systemic to random motion, V sin(i)/sigma, ranging from 0.5-1.3 and Toomre disk parameters Q<1. The large fraction of well-ordered velocity fields in our sample is consistent with data analyzed for larger, more luminous sources at this redshift. Our high resolution data further reveal that all six galaxies contain multiple giant star-forming HII regions whose resolved diameters are in the range 300 pc - 1.0 kpc, consistent with the Jeans length expected in the case of dispersion support. The density of star formation in these regions is ~100 times higher than observed in local spirals; such high values are only seen in the most luminous local starbursts. The global dynamics and demographics of star formation in these HII regions suggest that vigorous star formation is primarily governed by gravitational instability in primitive rotating disks.

We present a detailed analysis of the baryonic and dark matter distribution in the lensing cluster Abell 611 (z=0.288), with the goal of determining the dark matter profile over an unprecedented range of cluster-centric distance. By combining three complementary probes of the mass distribution, weak lensing from deep multi-color imaging, strong lensing constraints based on the identification of multiply-imaged sources, and resolved stellar velocity dispersion measures for the brightest cluster galaxy (BCG), we extend the methodology for separating the dark and baryonic mass components introduced by Sand et al. (2008). Our resulting dark matter profile samples the cluster from ~3 kpc to 3.25 Mpc, thereby providing an excellent basis for comparisons with recent numerical models. We demonstrate that only by combining our three observational techniques can degeneracies in constraining the form of the dark matter profile be broken on scales crucial for detailed comparisons with numerical simulations. Our analysis reveals that a simple Navarro, Frenk, and White (NFW) profile is an unacceptable fit to our data. We confirm earlier claims that the inner profile of the dark matter profile deviates significantly from the NFW form and find a inner logarithmic slope \beta flatter than 0.3 (68%; where rho_DM ~ r^{-\beta} at small radii). In order to reconcile our data with cluster formation in a LambdaCDM cosmology, we speculate that it may be necessary to revise our understanding of the nature of baryon--dark matter interactions in cluster cores. Comprehensive weak and strong lensing data, when coupled with kinematic information on the brightest cluster galaxy, can readily be applied to a larger sample of clusters to test the universality of these results.

Massive galaxy clusters at intermediate redshift can magnify the flux of distant background sources by several magnitudes and we exploit this effect to search for lensed distant supernovae that may otherwise be too faint to be detected. A supernova search was conducted at near infrared wavelengths using the ISAAC instrument at the VLT. The galaxy clusters Abell 1689, Abell 1835 and AC114 were observed at multiple epochs of 2 hours of exposure time, separated by a month. Image-subtraction techniques were used to search for transient objects with light curve properties consistent with supernovae, both in our new and archival ISAAC/VLT data. The limiting magnitude of the individual epochs was estimated by adding artificial stars to the subtracted images. Most of the epochs reach 90% detection efficiency at SZ(J) ~= 23.8-24.0 mag (Vega). Two transient objects, both in archival images of Abell 1689 and AC114, were detected. The transient in AC114 coincides - within the position uncertainty - with an X-ray source and is likely to be a variable AGN at the cluster redshift. The transient in Abell 1689 was found at SZ=23.24 mag, ~0.5 arcsec away from a galaxy with photometric redshift z=0.6 +/-0.15. The light curves and the colors of the transient are consistent with a reddened Type IIP supernova at redshift z=0.59 +/- 0.05. The lensing model of Abell 1689 predicts ~1.4 mag of magnification at the position of the transient, making it the most magnified supernova ever found and only the second supernova found behind a galaxy cluster. Our pilot survey has demonstrated the feasibility to find distant gravitationally magnified supernovae behind massive galaxy clusters. One likely supernova was found behind Abell 1689, in accordance with the expectations for this survey, as shown in an accompanying analysis paper.

We present a detailed study of the spatially resolved kinematics, star-formation and stellar mass in a highly amplified galaxy at z=4.92 behind the lensing cluster MS1358+62. We use the observed optical, near- and mid-infrared imaging from HST ACS & NICMOS and Spitzer IRAC to derive the stellar mass and the Gemini/NIFS IFU to investigate the velocity structure of the galaxy from the nebular [OII] emission. Using a detailed gravitational lens model, we account for lensing amplification factor 12.+/-2.0 and find that this intrinsically L* galaxy has a stellar mass of M*=7+/-2x10^8Mo, a dynamical mass of Mdyn=3+/-1x10^9csc^2(i)Mo (within of 2kpc) and a star-formation rate of 42+/-8Mo/yr. The source-plane UV/optical morphology of this galaxy is dominated by five discrete star-forming regions. Exploiting the dynamical information we derive masses for individual star-forming regions of Mcl~10^(8-9)Mo with sizes of ~200pc. We find that, at a fixed size, the star-formation rate density within these HII regions is approximately two orders of magnitude greater than those observed in local spiral/starburst galaxies, but consistent with the most massive HII regions in the local Universe such as 30Dor. Finally, we compare the spatially resolved nebular emission line velocity with the Ly-alpha and UV-ISM lines and find that this galaxy is surrounded by a galactic scale outflow in which the Ly-alpha appears redshifted by ~150km/s and the UV-ISM lines blue-shifted by -200km/s from the (systemic) nebular emission. The velocity structure of the outflow mirrors that of the nebular emission suggesting the outflow is young (~15Myr), and has yet to burst out of the system. Taken together, these results suggest that this young galaxy is undergoing its first major epoch of mass assembly.

Powerful gravitational telescopes in the form of massive galaxy clusters can be used to enhance the light collecting power over a limited field of view by about an order of magnitude in flux. This effect is exploited here to increase the depth of a survey for lensed supernovae at near-IR wavelengths. A pilot SN search program conducted with the ISAAC camera at VLT is presented. Lensed galaxies behind the massive clusters A1689, A1835 and AC114 were observed for a total of 20 hours split into 2, 3 and 4 epochs respectively, separated by approximately one month to a limiting magnitude J<24 (Vega). Image subtractions including another 20 hours worth of archival ISAAC/VLT data were used to search for transients with lightcurve properties consistent with redshifted supernovae, both in the new and reference data. The feasibility of finding lensed supernovae in our survey was investigated using synthetic lightcurves of supernovae and several models of the volumetric Type Ia and core-collapse supernova rates as a function of redshift. We also estimate the number of supernova discoveries expected from the inferred star formation rate in the observed galaxies. The methods consistently predict a Poisson mean value for the expected number of SNe in the survey between N_SN=0.8 and 1.6 for all supernova types, evenly distributed between core collapse and Type Ia SN. One transient object was found behind A1689, 0.5" from a galaxy with photometric redshift z_gal=0.6 +- 0.15. The lightcurve and colors of the transient are consistent with being a reddened Type IIP SN at z_SN=0.59. The lensing model predicts 1.4 magnitudes of magnification at the location of the transient, without which this object would not have been detected in the near-IR ground based search described in this paper (unlensed magnitude J~25). (abridged)

We present the first measurement of the relationship between the Sunyaev-Zel'dovich effect signal and the mass of galaxy clusters that uses gravitational lensing to measure cluster mass, based on 14 X-ray luminous clusters at z~0.2 from the Local Cluster Substructure Survey. We measure the integrated Compton y-parameter, Y, and total projected mass of the clusters (M_GL) within a projected clustercentric radius of 350 kpc, corresponding to mean overdensities of 4000-8000 relative to the critical density. We find self-similar scaling between M_GL and Y, with a scatter in mass at fixed Y of 32%. This scatter exceeds that predicted from numerical cluster simulations, however, it is smaller than comparable measurements of the scatter in mass at fixed T_X. We also find no evidence of segregation in Y between disturbed and undisturbed clusters, as had been seen with T_X on the same physical scales. We compare our scaling relation to the Bonamente et al. relation based on mass measurements that assume hydrostatic equilibrium, finding no evidence for a hydrostatic mass bias in cluster cores (M_GL = 0.98+/-0.13 M_HSE), consistent with both predictions from numerical simulations and lensing/X-ray-based measurements of mass-observable scaling relations at larger radii. Overall our results suggest that the Sunyaev-Zel'dovich effect may be less sensitive than X-ray observations to the details of cluster physics in cluster cores.

The existence of strong lensing systems with Einstein radii (Re) covering the full mass spectrum, from ~1-2" (produced by galaxy scale dark matter haloes) to >10" (produced by galaxy cluster scale haloes) have long been predicted. Many lenses with Re around 1-2" and above 10" have been reported but very few in between. In this article, we present a sample of 13 strong lensing systems with Re in the range 3"- 8", i.e. systems produced by galaxy group scale dark matter haloes, spanning a redshift range from 0.3 to 0.8. This opens a new window of exploration in the mass spectrum, around 10^{13}- 10^{14} M_{sun}, which is a crucial range for understanding the transition between galaxies and galaxy clusters. Our analysis is based on multi-colour CFHTLS images complemented with HST imaging and ground based spectroscopy. Large scale properties are derived from both the light distribution of the elliptical galaxies group members and weak lensing of the faint background galaxy population. On small scales, the strong lensing analysis yields Einstein radii between 2.5" and 8". On larger scales, the strong lenses coincide with the peak of the light distribution, suggesting that mass is traced by light. Most of the luminosity maps have complicated shapes, indicating that these intermediate mass structures are dynamically young. Fitting the reduced shear with a Singular Isothermal Sphere, we find sigma ~ 500 km/s and an upper limit of ~900 km/s for the whole sample. The mass to light ratio for the sample is found to be M/L_i ~ 250 (solar units, corrected for evolution), with an upper limit of 500. This can be compared to mass to light ratios of small groups (with sigma ~ 300 km/s and galaxy clusters with sigma > 1000 km/s, thus bridging the gap between these mass scales.

We report the results of a Spitzer infrared study of the Cosmic Eye, a strongly lensed, L*_UV Lyman Break Galaxy (LBG) at z=3.074. We obtained Spitzer IRS spectroscopy as well as MIPS 24 and 70 micron photometry. The Eye is detected with high significance at both 24 and 70 microns and, when including a flux limit at 3.5 mm, we estimate an infrared luminosity of L_IR = 8.3 (+4.7-4.4) x10^11 L_sun assuming a magnification of 28+-3. This L_IR is eight times lower than that predicted from the rest-frame UV properties assuming a Calzetti reddening law. This has also been observed in other young LBGs, and indicates that the dust reddening law may be steeper in these galaxies. The mid-IR spectrum shows strong PAH emission at 6.2 and 7.7 microns, with equivalent widths near the maximum values observed in star-forming galaxies at any redshift. The L_PAH-to-L_IR ratio lies close to the relation measured in local starbursts. Therefore, L_PAH or L_MIR may be used to estimate L_IR and thus, star formation rate, of LBGs, whose fluxes at longer wavelengths are typically below current confusion limits. We also report the highest redshift detection of the 3.3 micron PAH emission feature. The PAH ratio, L_6.2/L_3.3=5.1+- 2.7, and the PAH-to-L_IR ratio, L_3.3/L_IR = 8.5 +- 4.7 x10^-4, are both in agreement with measurements in local starbursts and ULIRGs, suggesting that this line may serve as a good proxy for L_PAH or L_IR at z > 3 with the James Webb Space Telescope.

Strong gravitational lensing is a unique tool to model with great accuracy the inner mass distribution of massive galaxy clusters. In particular, clusters with large Einstein radii provide a wealth of multiply imaged systems in the cluster core allowing to determine precisely the shape of the central dark matter profile. This paper presents a spectroscopic survey in the massive cluster Abell 1703, displaying a large Einstein radius (28" at z=2.8) and a high number of strongly-lensed systems including a central ring-like configuration. We used LRIS on Keck to target multiple images and lensed galaxy candidates, and use the measured redshifts to constrain the mass distribution of the cluster using a parametric model. The data enable us to measure accurate redshifts in good agreement with their photometric redshifts, and to update the identification of multiply imaged systems by discovering 3 new systems and a radial counter image. We also report the discovery of a remarkably bright ~3.6 L* i-band dropout at z=5.827 in our mask which is only moderately magnified by the cluster (~3.0+/-0.08). The improved parametric mass model, including 16 multiple systems with 10 spectroscopic redshifts, further constrains the cluster-scale mass distribution with a generalized NFW profile of best-fit logarithmic slope alpha=0.92+/-0.04, concentration c200=4.72+/-0.40 and scale radius rs=476+/-45 kpc. Our strong-lensing model predicts a large scale shear signal consistent with Subaru weak-lensing measurements out to 4 Mpc h^-1. Together with the fact that the strong-lensing modeling requires a single dark matter clump, this suggests that Abell 1703 is a relaxed, unimodal cluster. This unique cluster could be probed further using deep X-ray, SZ and dynamics analysis, for a detailed study of the physics in a relaxed cluster. (abridged)

Recent studies of galaxies ~2-3 Gyr after the Big Bang have revealed large, turbulent rotating systems. The existence of well-ordered rotation in galaxies during this peak epoch of cosmic star formation may suggest that gas accretion through cold streams is likely to be the dominant mode by which most star-forming galaxies at high redshift since major mergers can completely disrupt the observed velocity fields. However poor spatial resolution and sensitivity have hampered this interpretation, limiting the study to the largest and most luminous galaxies, which may have fundamentally different modes of assembly than more typical star forming galaxies. Here we report observations of a typical star forming galaxy at z=3.07 with a linear resolution of ~100 parsec. This spatial sampling is made possible by the combination of gravitational lensing and laser guide star adaptive optics. We find a well-ordered compact source in which molecular gas is being converted efficiently into stars, likely assembling a spheroidal bulge and disk similar to those seen in spiral galaxies at the present day.

We analyze a new gravitational lens, OAC-GL J1223-1239, serendipitously found in a deep I-band image of the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). The lens is a L_*, edge-on S0 galaxy at z=0.4656. The gravitational arc has a radius of 0.42 arcsec. We have determined the total mass and the dark matter (DM) fraction within the Einstein radius as a function of the lensed source redshift, which is presently unknown. For z ~ 1.3, which is in the middle of the redshift range plausible for the source according to some external constraints, we find the central velocity dispersion to be ~180 km/s. With this value, close to that obtained by means of the Faber-Jackson relation at the lens redshift, we compute a 30% DM fraction within the Einstein radius (given the uncertainty in the source redshift, the allowed range for the DM fraction is 25-35 % in our lensing model). When compared with the galaxies in the local Universe, the lensing galaxy, OAC-GL J1223-1239 seems to fall in the transition regime between massive DM dominated galaxies and lower-mass, DM deficient systems.

We present the results of a systematic search for gravitationally-lensed continuum Lyman break `drop-outs' beyond a redshift 7 conducted via very deep imaging through six foreground clusters undertaken with the Hubble and Spitzer Space Telescopes. The survey has yielded 10 z-band and 2 J-band drop-out candidates to photometric limits of J_110~=26.2 AB (5sigma). Taking into account the magnifications afforded by our clusters (1-4 magnitudes), we probe the presence of z>7 sources to unlensed limits of J_{110}~=30 AB, fainter than those charted in the Hubble Ultradeep Field. To verify the fidelity of our candidates we conduct a number of tests for instrumental effects which would lead to spurious detections, and carefully evaluate the likelihood of foreground contamination by considering photometric uncertainties in the drop-out signature, the upper limits from stacked IRAC data and the statistics of multiply-imaged sources. Overall, we conclude that we can expect about half of our sample of z-band drop-outs are likely to be at high redshift. An ambitious infrared spectroscopic campaign undertaken with the NIRSPEC spectrograph at the WM Keck Observatory for seven of the most promising candidates failed to detect any Lyman-alpha emission highlighting the challenge of making further progress in this field. While the volume density of high redshift sources will likely remain uncertain until more powerful facilities are available, our data provides the first potentially interesting constraints on the UV luminosity function at z~=7.5 at intrinsically faint limits. We discuss the implications of our results in the context of the hypothesis that the bulk of the reionizing photons in the era 7<z<12 arise in low luminosity galaxies undetected by conventional surveys.

In this article, we apply strong lensing techniques in Abell 1703, a massive X-ray luminous galaxy cluster at z=0.28. Our analysis is based on imaging data both from space and ground in 8 bands, complemented with a spectroscopic survey. Abell 1703 looks rather circular from the general shape of its multiply imaged systems and present a dominant giant elliptical cD galaxy in its centre. This cluster exhibits a remarkable bright 'central ring' formed by 4 bright images at z_{spec}=0.888 located very close to the cD galaxy, providing observational constraints that are potentially very interesting to probe the central mass distribution. The stellar contribution from the cD galaxy (~1.25 10^{12} M_{sun} within 7") is accounted for in our parametric mass modelling, and the underlying smooth dark matter component distribution is described using a generalized NFW profile parametrized with a central logarithmic slope \alpha. We find that within the range where observational constraints are present (from ~5" to ~50"), the slope of the dark matter distribution in Abell1703 is equal to 1.09^{+0.05}_{-0.11} (3\sigma confidence level). The concentration parameter is equal to c_{200} ~ 3.5, and the scale radius is constrained to be larger than the region where observational constraints are available. Within this radius, the 2D mass is equal to M(50")=2.4 10^{14} M_{\sun}. We cannot draw any conclusions on cosmological models at this point since we lack results from realistic numerical simulations containing baryons to make a proper comparison. We advocate the need for a sample of observed and simulated unimodal relaxed galaxy clusters in order to make reliable comparisons, and potentially provide a test of cosmological models.

The strong lensing modelling of gravitational ``rings'' formed around massive galaxies is sensitive to the amplitude of the external shear and convergence produced by nearby mass condensations. In current wide field surveys, it is now possible to find out a large number of rings, typically 10 gravitational rings per square degree. We propose here, to systematically study gravitational rings around galaxy clusters to probe the cluster mass profile beyond the cluster strong lensing regions. For cluster of galaxies with multiple arc systems, we show that rings found at various distances from the cluster centre can improve the modelling by constraining the slope of the cluster mass profile. We outline the principle of the method with simple numerical simulations and we apply it to 3 rings discovered recently in Abell~1689. In particular, the lens modelling of the 3 rings confirms that the cluster is bimodal, and favours a slope of the mass profile steeper than isothermal at a cluster radius $\sim 300 \kpc$. These results are compared with previous lens modelling of Abell~1689 including weak lensing analysis. Because of the difficulty arising from the complex mass distribution in Abell~1689, we argue that the ring method will be better implemented on simpler and relaxed clusters.

We present a reconstruction of the mass distribution of galaxy cluster Abell 1689 at z = 0.18 using detected strong lensing features from deep HST/ACS observations and extensive ground based spectroscopy. Earlier analyses have reported up to 32 multiply imaged systems in this cluster, of which only 3 were spectroscopically confirmed. In this work, we present a parametric strong lensing mass reconstruction using 24 multiply imaged systems with newly determined spectroscopic redshifts, which is a major step forward in building a robust mass model. In turn, the new spectroscopic data allows a more secure identification of multiply imaged systems. The resultant mass model enables us to reliably predict the redshifts of additional multiply imaged systems for which no spectra are currently available, and to use the location of these systems to further constrain the mass model. In particular, we have detected 5 strong galaxy-galaxy lensing systems just outside the Einstein ring region, further constraining the mass profile. Our strong lensing mass model is consistent with that inferred from our large scale weak lensing analysis derived using CFH12k wide field images. Thanks to a new method for reliably selecting a well defined background lensed galaxy population, we resolve the discrepancy found between the strong and weak lensing mass models reported in earlier work. [ABRIDGED]

We present Spitzer/IRS mid-infrared spectra for 15 gravitationally lensed, 24 micron--selected galaxies, and combine the results with 4 additional very faint galaxies with IRS spectra in the literature. The median intrinsic 24 micron flux density of the sample is 130 microJy, enabling a systematic survey of the spectral properties of the very faint 24 micron sources that dominate the number counts of Spitzer cosmological surveys. Six of the 19 galaxy spectra (32%) show the strong mid-IR continuua expected of AGN; X-ray detections confirm the presence of AGN in three of these cases, and reveal AGNs in two other galaxies. These results suggest that nuclear accretion may contribute more flux to faint 24 micron--selected samples than previously assumed. Almost all the spectra show some aromatic (PAH) emission features; the measured aromatic flux ratios do not show evolution from z=0. In particular, the high S/N mid-IR spectrum of SMM J163554.2+661225 agrees remarkably well with low--redshift, lower--luminosity templates. We compare the rest-frame 8 micron and total infrared luminosities of star--forming galaxies, and find that the behavior of this ratio with total IR luminosity has evolved modestly from z=2 to z=0. Since the high aromatic--to--continuum flux ratios in these galaxies rule out a dominant contribution by AGN, this finding implies systematic evolution in the structure and/or metallicity of infrared sources with redshift. It also has implications for the estimates of star forming rates inferred from 24 micron measurements, in the sense that at z ~2, a given observed frame 24 micron luminosity corresponds to a lower bolometric luminosity than would be inferred from low-redshift templates of similar luminosity at the corresponding rest wavelength.

We present a parametric strong lensing model of the cluster Abell 2218 based on HST ACS data. We constrain the lens model using 8 previously known multiply imaged systems, of which 7 have spectroscopically confirmed redshifts. In addition, we propose five candidate multiply imaged systems and estimate their redshifts using our mass model. The model parameters are optimized in the source plane by a bayesian Monte Carlo Markov Chain as implemented in the the publicly available software Lenstool. We find rms=0."12 for the scatter of the sources in the source plane, which translates into rms=1."49 between the predicted and measured image positions in the image plane. We find that the projected mass distribution of Abell 2218 is bimodal, which is supported by an analysis of the light distribution. We also find evidence for two structures in velocity space, separated by ~1000 km/s, corresponding to the two large scale dark matter clumps. We find that the lensing constraints can not be well reproduced using only dark matter halos associated with the cluster galaxies, but that the dark matter is required to be smoothly distributed in large scale halos. At 100" the enclosed projected mass is 3.8e14 solar masses. At that radius, the large scale halos contribute ~85% of the mass. We find that the model is not very sensitive to the fainter galaxy sized halos, unless they locally perturb a given multiply imaged system. Therefore, dark galaxy sized substructure can be reliably constrained only if it locally perturbs one of the systems. In an appendix we give a self-contained description of the parametric profile we use, the dual pseudo isothermal elliptical mass distribution (dPIE). This profile is a two component pseudo isothermal mass distribution (PIEMD) with both a core radii and a scale radii. (Abridged)

We study the properties and nature of extremely red galaxies (ERO, R-K>5.6) found behind two lensing clusters and compare them with other known galaxy populations. New HST/ACS observations, Spitzer IRAC and MIPS, and Chandra/ACIS observations of the two lensing clusters Abell 1835 and AC114 contemplate our earlier optical and near-IR observations and have been used to study extremely red objects (EROs) in these deep fields. We have found 6 and 9 EROs in Abell 1835 and AC114. Several (7) of these objects are undetected up to the I and/or z band, and are hence ``optical'' drop-out sources. The photometric redshifts of most of our sources (80%) are z~0.7-1.5. According to simple colour-colour diagrams the majority of our objects would be classified as hosting old stellar populations. However, there are clear signs of dusty starbursts for several among them. These objects correspond to the most extreme ones in R-K colour. We estimate a surface density of (0.97+-0.31) arcmin-2 for EROs with (R-K>5.6) at K<20.5. Among our 15 EROs 6 (40 %) also classify as distant red galaxies (DRGs). 11 of 13 EROs with available IRAC photometry also fulfil the selection criteria for IRAC selected EROs (IEROs) of Yan et al. (2004). SED modelling shows that ~ 36 % of the IEROs in our sample are luminous or ultra-luminous infrared galaxies ((U)LIRG). Some very red DRGs are found to be very dusty starbursts, even (U)LIRGs, as also supported by their mid-IR photometry. No indication for AGNs is found, although faint activity cannot be excluded for all objects. From mid-IR and X-ray data 5 objects are clearly classified as starbursts. The derived properties are quite similar to those of DRGs and IEROs, except for 5 extreme objects in terms of colours, for which a very high extinction (Av>3) is found.

We report the detection of CO(3-2) emission from a bright, gravitationally lensed Lyman Break Galaxy, LBGJ213512.73-010143 (the 'Cosmic Eye'), at z=3.07 using the Plateau de Bure Interferometer. This is only the second detection of molecular gas emission from an LBG and yields an intrinsic molecular gas mass of (2.4+/-0.4)x10^9 Mo. The lens reconstruction of the UV morphology of the LBG indicates that it comprises two components separated by ~2 kpc. The CO emission is unresolved, and appears to be centered on the intrinsically fainter (and also less highly magnified) of the two UV components. The width of the CO line indicates a dynamical mass of (8+/-2)x10^9csc(i)^2 Mo within the central 2 kpc. Employing mid-infrared observations from Spitzer we derive a stellar mass of ~(6+/-2)x10^9 Mo and a star-formation rate of ~60 Mo/yr, indicating that the molecular gas will be consumed in ~40 Myr. The gas fractions, star-formation efficiencies and line widths suggests that LBGJ213512 is a high-redshift, gas-rich analog of a local luminous infrared galaxy. This galaxy has a similar gas-to-dynamical mass fraction as observed in the submillimeter-selected population, although the gas surface density and star-formation efficiency is a factor of 3x less, suggesting less vigorous activity. We discuss the uncertainties in our conclusions arising from adopting a CO-to-H2 conversion factor appropriate for either the Milky Way or local luminous infrared galaxies. These observations demonstrate that current facilities, when aided by fortuitous gravitational magnification, can study 'ordinary' galaxies at high-redshift and so act as pathfinders for ALMA.

We determine the nature, redshift, stellar populations and dust properties of optically faint or non-detected extremely red objects (ERO) found from our survey of the lensing clusters A1835 and AC114. We perform SED fitting and use deep optical, HST, VLT, Spitzer data, and for some objects and sub-mm data. For most of the lensed EROs we find photometric redshifts showing a strong degeneracy between "low-z" (z~1-3) and high-z (z~6-7). Although formally best fits are often found at high-z, their resulting bright absolute magnitudes, the number density of these objects, and in some cases Spitzer photometry or longer wavelength observations, suggest strongly that all of these objects are at "low-z". The majority of these objects are best fitted with young (<~ 0.5-0.7 Gyr) and dusty starbursts. Indications for strong extinction, with A_V~2.4-4, are found in some objects. For dusty objects star formation rates (SFR) have been estimated from the bolometric luminosity determined after fitting of semi-empirical starburst, ERO, and ULIRG templates. Typically we find SFR~(1-18) Msun/yr. Again, SMMJ14009+0252 stands out as a LIRG with SFR 1000 Msun/yr. Concerning the comparison objects, we argue that the massive post-starburst z~6.5 galaxy candidate HUDF-J2 showing observed properties very similar to our EROs, is more likely a dusty starburst at z~2.3-2.6. This interpretation also naturally explains the observed 24 micron emission. Both empirically and from our SED fits we find that the IRAC selectec EROs from Yan et al (2004) show very similar properties to our lensed EROs. Reasonable fits are found for most of them with relatively young and dusty stellar populations.

We describe two peculiar galaxies falling into the massive galaxy clusters Abell 1689 (z~0.18) and 2667 (z~0.23) respectively. Hubble Space Telescope images show extraordinary trails composed of bright blue knots (-16.5<M<-11.5 mag) and stellar streams associated with each of these systems. Combining optical, near and mid-infrared and radio observations we prove that while both galaxies show similar extended trails of star-forming knots, their recent star formation histories are different. One (~L*) is experiencing a strong burst of star formation, appearing as a rare example of a luminous infrared cluster galaxy. In comparison, the other (~ 0.1 L*) has recently ceased its star formation activity. Our model suggests that the morphologies and star formation in these galaxies have been influenced by the combined action of tidal interaction (likely with the cluster potential) and of ram pressure with the intracluster medium. These results can be used to gain more insights to the origin of S0s, dwarf and ultra-compact dwarf (UCD) cluster galaxies.

We have carried out an extensive spectroscopic survey with the Keck and VLT telescopes, targeting lensed galaxies in the background of the massive cluster Abell 68. Spectroscopic measurements are obtained for 26 lensed images, including a distant galaxy at z=5.4 . Redshifts have been determined for 5 out of 7 multiply-image systems. Through a careful modeling of the mass distribution in the strongly-lensed regime, we derive a mass estimate of 5.3 x 10^14 Msun within 500 kpc. Our mass model is then used to constrain the redshift distribution of the remaining multiply-imaged and singly-imaged sources. This enables us to examine the physical properties for a subsample of 7 Lyman-alpha emitters at 1.7 < z < 5.5, whose unlensed luminosities of ~ 10^41 ergs/s are fainter than similar objects found in blank fields. Of particular interest is an extended Lyman-alpha emission region surrounding a highly magnified source at z=2.6, detected in VIMOS Integral Field Spectroscopy data. The physical scale of the most distant lensed source at z=5.4 is very small (<300 pc), similar to the lensed z ~ 5.6 emitter reported by Ellis et al. (2001) in Abell 2218. New photometric data available for Abell 2218 allow for a direct comparison between these two unique objects. Our survey illustrates the practicality of using lensing clusters to probe the faint end of the z ~ 2-5 Lyman-alpha luminosity function in a manner that is complementary to blank field narrow-band surveys.

We discuss new observational constraints on the abundance of faint high redshift Lyman-alpha emitters secured from a deep Keck near-infrared spectroscopic survey which utilizes the strong magnification provided by lensing galaxy clusters at intermediate redshift (z=0.2-0.5). In each of 9 clusters we have undertaken a systematic `blind' search for line emission with NIRSPEC in the J-band within carefully-selected regions which offer very high magnifications (>10-50x) for background sources with redshifts z~10. The high magnification enables the detection of emission at unprecedented limits (10^41 - 10^42 erg/s), much fainter than those of conventional narrow band imaging and other spectroscopic searches. As the comoving volumes probed are small, our survey is designed to address the important question of whether intrinsically feeble star forming galaxies could provide the dominant ionizing flux at z~10. Our survey has yielded six promising (>5 sigma) candidate Lyman-alpha emitters which lie between z=8.7 and z=10.2. While the predicted reionization photon budget depends upon a large number of physical assumptions, our first glimpse at the z~10 Universe suggests that low luminosity star-forming galaxies contribute a significant proportion of the UV photons necessary for cosmic reionization. [Abridged]

We summarise the current status of our project to identify and study z~6-10 galaxies thanks to strong gravitational lensing. Building on the detailed work from Richard et al. (2006), we present results from new follow-up observations (imaging) undertaken with ACS/HST and the Spitzer Space Telescope and compare our results with findings from the Hubble Ultra-Deep Field (UDF). These new observations are in agreement with the high-z nature for the vast majority of the candidates presented in Richard et al. (2006). We also discuss the properties of other optical dropout sources found in our searches and related objects (EROs, sub-mm galaxies,...) from other surveys.

We report the discovery using Hubble Space Telescope imaging and Keck spectroscopy of a very bright, highly magnified (~30*) Lyman Break Galaxy (LBG) at z=3.07 in the field of the massive z=0.33 cluster MACSJ2135.2-0102. The system comprises two high-surface brightness arcs with a maximum extent of 3", bracketing a central object which we identify as a massive early-type galaxy at z=0.73. We construct a lens model which reproduces the main features of the system using a combination of a galaxy-scale lens and the foreground cluster. We show that the morphological, spectral and photometric properties of the arcs are consistent with them arising from the lensing of a single ~L* LBG. The most important feature of this system is that the lensing magnification results in an apparent magnitude of r=20.3, making this one of the brightest LBGs known. Such a high magnification provides the opportunity of obtaining very high signal to noise (and potentially spatially resolved) spectroscopy of a high redshift galaxy to study its physical properties. We present initial imaging and spectroscopy demonstrating the basic properties of the system and discuss the opportunities for future observations.

Observing the first galaxies formed during the reionisation epoch, i.e. approximately within the first billion years after the Big Bang, remains one of the challenges of contemporary astrophysics. Several efforts are being undertaken to search for such remote objects. Combining the near-IR imaging power of the VLT and the natural effect of strong gravitational lensing our pilot program has allowed us to identify several galaxy candidates at redshift 6 <~ z <~ 10. The properties of these objects and the resulting constraints on the star formation rate density at high redshift are discussed. Finally we present the status of follow-up observations (ISAAC spectroscopy, HST and Spitzer imaging) and discuss future developments.

(abridged) We present the first results of our deep survey of lensing clusters aimed at constraining the abundance of star-forming galaxies at z~6-10. Deep near-IR photometry of two lensing clusters (A1835 and AC114) was obtained with ISAAC/VLT. These images, combined with existing data in the optical bands, including HST images, were used to select very high redshift candidates at z>~6 among the optical-dropouts. We have identified 18(8) first and second-category optical dropouts in A1835 (AC114), detected in more than one filter up to H(Vega)~23.8 (AB~25.2,uncorrected for lensing). Among them, 8(5) exhibit homogeneous SEDs compatible with star-forming galaxies at z>~6, and 5(1) are more likely intermediate-redshift EROs based on luminosity considerations. We have also identified a number of fainter sources in these fields fulfilling our photometric selection and located around the critical lines. We use all these data to make a first attempt at constraining the density of star-forming galaxies present at 6<z<10 using lensing clusters. Taken at face value, it appears that the number of candidates found seems to be higher than the one achieved in blank fields with similar photometric depth in the near-IR. The luminosity function derived for z>6 candidates appears compatible with that of LBGs at z~3. The turnover observed by Bouwens et al. (2005) towards the bright end relative to the z~3 LF is not observed in this sample. Also the upper limit for the UV SFR density at z~6-10 is compatible with the usual values derived at z~5-6, but higher than the estimates obtained in the NICMOS Ultra Deep Field (UDF). Increasing the number of blank and lensing fields with ultra-deep near-IR photometry is essential to get more accurate constraints on the abundance of z>6 galaxies.

We present the on-going observational program of a VIMOS Integral Field Unit survey of the central regions of massive, gravitational lensing galaxy clusters at redshift z~0.2. We have observed six clusters using the low-resolution blue grism (R about 200), and the spectroscopic survey is complemented by a wealth of photometric data, including Hubble Space Telescope optical data and near infrared VLT data. The principal scientific aims of this project are: the study of the high-z lensed galaxies, the transformation and evolution of galaxies in cluster cores and the use of multiple images to constrain cosmography. We briefly report here on the first results from this project on the clusters Abell 2667 and Abell 68.

We present the gas-phase oxygen abundance (O/H) for a sample of 131 star-forming galaxies at intermediate redshifts ($0.2<z<1.0$). The sample selection, the spectroscopic observations (mainly with VLT/FORS) and associated data reduction, the photometric properties, the emission-line measurements, and the spectral classification are fully described in a companion paper (Paper I). We use two methods to estimate the O/H abundance ratio: the {} "standard" $R\_{\mathrm{23}}$ method which is based on empirical calibrations, and the CL01 method which is based on grids of photo-ionization models and on the fitting of emission lines. For most galaxies, we have been able to solve the problem of the metallicity degeneracy between the high- and low-metallicity branches of the O/H vs. $R\_{\mathrm{23}}$ relationship using various secondary indicators. The luminosity -- metallicity ($L-Z$) relation has been derived in the $B$- and $R$-bands, with metallicities derived with the two methods ($R\_{\mathrm{23}}$ and CL01). In the analysis, we first consider our sample alone and then a larger one which includes other samples of intermediate-redshift galaxies drawn from the literature. The derived $L-Z$ relations at intermediate redshifts are very similar (same slope) to the $L-Z$ relation obtained for the local universe. Our sample alone only shows a small evolution of the $L-Z$ relation with redshift up to $z\sim1.0$. We only find statistical variations consistent with the uncertainty in the derived parameters. Including other samples of intermediate-redshift galaxies, we find however that galaxies at $z\sim1$ appear to be metal-deficient by a factor of $\sim3$ compared with galaxies in the local universe. For a given luminosity, they contain on average about one third of the metals locked in local galaxies.

We present the spectrophotometric properties of a sample of 141 emission-line galaxies at redshifts in the range $0.2<z<1.0$ with a peak around $z\in[0.2,0.4]$. The analysis is based on medium resolution ($R\_{\mathrm{s}}=500-600$), optical spectra obtained at VLT and Keck. The targets are mostly {} "Canada-France Redshift Survey" emission-line galaxies, with the addition of field galaxies randomly selected behind lensing clusters. We complement this sample with galaxy spectra from the {} "Gemini Deep Deep Survey" public data release. We have computed absolute magnitudes of the galaxies and measured the line fluxes and equivalent widths of the main emission/absorption lines. The last two have been measured after careful subtraction of the fitted stellar continuum using the \texttt{platefit} software originally developed for the SDSS and adapted to our data. We present a careful comparison of this software with the results of manual measurements. The pipeline has also been tested on lower resolution spectra, typical of the {} "VIMOS/VLT Deep Survey" ($R\_{\mathrm{s}}=250$), by resampling our medium resolution spectra. We show that we can successfully deblend the most important strong emission lines. These data are primarily used to perform a spectral classification of the galaxies in order to distinguish star-forming galaxies from AGNs. Among the initial sample of 141 emission-line galaxies, we find 7 Seyfert 2 (narrow-line AGN), 115 star-forming galaxies and 16 {} "candidate" star-forming galaxies. Scientific analysis of these data, in terms of chemical abundances, stellar populations, etc, will be presented in subsequent papers of this serie.

(abridged) We present extensive multi-color imaging and low resolution VIMOS Integral Field Unit spectroscopic observations of the X-ray luminous cluster Abell 2667 (z=0.233). An extremely bright giant gravitational arc (z=1.0334) is easily identified as part of a triple image system and other fainter multiple images are also revealed by the HST-WFPC2 images. The VIMOS-IFU observations cover a field of view of 54'' x 54'' and enable us to determine the redshift of all galaxies down to V=22.5. Furthermore, redshifts could be identified for some sources down to V=23.2. In particular we identify 21 cluster members in the cluster inner region, from which we derive a velocity dispersion of \sigma=960 km/s, corresponding to a total mass of 7.1 x 10^{13} solar masses within a 110 kpc radius. Using the multiple images constraints and priors on the mass distribution of cluster galaxy halos we construct a detailed lensing mass model leading to a total mass of 2.9 x 10^{13} solar masses within the Einstein radius (16 arcsec). The lensing mass and dynamical mass are in good agreement although the dynamical one is much less accurate. Comparing these measurements with published X-ray analysis, is however less conclusive. Although the X-ray temperature matches the dynamical and lensing estimates, the published NFW mass model derived from the X-ray measurement with its small concentration of c ~3 can not account for the large Einstein radius observed in this cluster. A larger concentration of ~6 would however match the strong lensing measurements. These results are likely reflecting the complex structure of the cluster mass distribution, underlying the importance of panchromatic studies from small to large scale in order to better understand cluster physics.

We summarize the main results obtained recently by our group on the identification and study of very high-z galaxies (z>7) using lensing clusters as natural gravitational telescopes. A description of our pilot survey with ISAAC/VLT is presented, aimed at the spectroscopic confirmation of z>7 candidate galaxies photometrically selected from deep near-IR, HST and optical ground-based imaging. The first results issued from this survey are discussed, in particular the global photometric properties of our high-z candidates, and the implications for the global star formation rate at very high-z.

We report the detection of a z~7 galaxy strongly lensed by the massive galaxy cluster Abell 2218 (z=0.175) at 3.6 and 4.5 um using the Spitzer Observatory and at 1.1 um using the Hubble Space Telescope. The new data indicate a refined photometric redshift in the range of 6.6-6.8 depending on the presence of Ly-alpha emission. The spectral energy distribution is consistent with having a significant Balmer break, suggesting that the galaxy is in the poststarburst stage with an age of at least ~50 Myr and quite possibly a few hundred Myr. This suggests the possibility that a mature stellar population is already in place at such a high redshift. Compared with typical Lyman break galaxies at z~3-4, the stellar mass is an order of magnitude smaller (~10^{9} Msun), but the specific star formation rate (star formation rate/M_{star}) is similarly large (> 10^{-9} yr^{-1}), indicating equally vigorous star-forming activity.

We review in this paper the main results recently obtained on the identification and study of very high-z galaxies usinglensing clusters as natural gravitational telescopes. We present in detail our pilot survey with ISAAC/VLT, aimed at the detection of z>7 sources. Evolutionary synthesis models for extremely metal-poor and PopIII starbursts have been used to derive the observational properties expected for these high-z galaxies, such as expected magnitudes and colors, line fluxes for the main emission lines, etc. These models have allowed to define fairly robust selection criteria to find z~7-10 galaxies based on broad-band near-IR photometry in combination with the traditional Lyman drop-out technique. The first results issued from our photometric and spectroscopic survey are discussed, in particular the preliminary confirmation rate, and the global properties of our high-z candidates, including the latest results on the possible z=10.0 candidate A1835-1916. The search efficiency should be significantly improved by the future near-IR multi-object ground-based and space facilities. However, strong lensing clusters remain a factor of ~5-10 more efficient than blank fields in this redshift domain, within the FOV of a few arcminutes around the cluster core, for the typical depth required for this survey project.

In a recent posting of a submitted paper Weatherley et al. re-analyse the spectroscopic data from our observations of Abell 1835 IR 1916, a galaxy we claimed to be at redshift z=10.0 (Pello et al., 2004, A&A, 416, L35). They conclude that the emission line is not detected in their analysis. Although their reduction procedure is accurate and robust, there are several differences with ours. The most important one concerns the stacking of the spectral frames. In the case of compact target sources [here: seeing-limited object,~3 pixels FWHM in the spectral direction, on a 1" (~7 pixels) wide slit] spectral shifts occur between individual frames. Our reduction corrects for this effect using a nearby strong emission line from a compact reference object, as described in Pello et al. However, Weatherley et al. neglect this effect. This difference turns out to be crucial. E.g. adopting their combination scheme as closely as possible but correcting for the spectral drifts we recover the line detection. We conclude that more work is needed to fully understand the origin of our differences. However, their non-detection claim needs further investigation and cannot be uphold based on the arguments presented in their original version of the paper. More information is beeing posted on http://webast.ast.obs-mip.fr/people/roser/z10_spectroscopy_discussion. With this material we also invite other researchers to re-analyse independently the ISAAC spectroscopic observations.

We present the first results obtained from our pilot ultra-deep near-IR survey with ISAAC/VLT, aimed at the detection of z>7 sources using lensing clusters as natural Gravitational Telescopes. Evolutionary synthesis models of PopIII and extremely metal-poor starbursts have been used to derive observational properties expected for these high-z galaxies, such as expected magnitudes and colors, line fluxes for the main emission lines, etc. These models have allowed us to define fairly robust selection criteria to find z~ 7-10 galaxies based on broad-band near-IR photometry in combination with the traditional Lyman drop-out technique. The magnification in the core of lensing clusters improves the search efficiency and subsequent spectroscopic follow up. The method used to identify sources and the first results obtained are presented, in particular the first spectroscopic confirmation of a z=10 candidate. The research efficiency will be significantly improved by the future near-IR multi-object facilities such as EMIR/GTC and KMOS/VLT-2, and later by the JWST.

We report the first likely spectroscopic confirmation of a z 10.0 galaxy from our ongoing search for distant galaxies with ISAAC/VLT. Galaxy candidates at z >~ 7 are selected from ultra-deep JHKs images in the core of gravitational lensing clusters for which deep optical imaging is also available, including HST data. The object reported here, found behind Abell 1835, exhibits a faint emission line detected in the J band, leading to z=10.0 when identified as Ly-a, in excellent agreement with the photometric redshift determination. Redshifts z < 7 are very unlikely for various reasons we discuss. The object is located on the critical lines corresponding to z=9 to 11. The magnification factor \mu ranges from 25 to 100. For this object we estimate SFR(Ly-a) (0.8-2.2) Msun/yr and SFR(UV) (47-75) Msun/yr, both uncorrected for lensing. The steep UV slope indicates a young object with negligible dust extinction. SED fits with young low-metallicity stellar population models yield (adopting mu=25) a lensing corrected stellar mass of M*~8.e+6 Msun, and luminosities of 2.e+10 Lsun, corresponding to a dark matter halo of a mass of typically M_tot>~ 5.e+8 Msun. In general our observations show that under excellent conditions and using strong gravitational lensing direct observations of galaxies close to the ``dark ages'' are feasible with ground-based 8-10m class telescopes.

We discuss the observational properties of a remarkably faint triply-imaged galaxy revealed in a deep z-band ACS observation of the lensing cluster A2218. A well-constrained mass model for the cluster, which incorporates the outcome of recent Keck spectroscopic campaigns, suggests that the triple system arises via a z>6 source viewed at high magnification (~x25). Optical and IR photometry from HST and Keck confirms the lensing hypothesis and suggests a significant discontinuity occurs in the spectral energy distribution within 9250-9850A. If this break is associated with Gunn-Peterson absorption from neutral hydrogen,6.6<z<7.1 is inferred. Deep Keck spectroscopy conducted using both optical and IR spectrographs fails to reveal any prominent emission lines in this region. However, an IR stellar continuum is detected whose decline below 9800A suggests a spectroscopic redshift of z~7. Regardless of the precise redshift, the source is remarkably compact (<1 kpc) and faint (z=28.0) yet is undergoing vigorous star formation at a rate ~2.6 Mo/yr. An intriguing property is the steep slope of the UV continuum implied by the photometry which may suggest that the source is representative of an early population of galaxies responsible for cosmic reionization. Independent verification of these results is highly desirable but our attempts highlight the difficulty of studying such sources with present facilities and the challenges faced in pushing back the frontiers of the observable universe beyond z~6.5.

We report the discovery of an unusual emission line galaxy at redshift z=1.68 found from near-IR imaging and spectroscopic observations with ISAAC/VLT of a faint gravitationally lensed R-band drop-out behind the lensing cluster Abell 1835. From the emission lines of [O_III], and H_beta this galaxy shows a moderate to high excitation and line equivalent widths typical of relatively metal-poor HII galaxies. Its apparent J magnitude translates to an absolute B-band magnitude M_B ~ -16.4 once corrected for a gravitational magnification of 1.8 magnitudes. This makes it one of the faintest intermediate redshift galaxiesknown so far. From the presence of emission lines and the available VRIJHK photometry we determine constraints on its dominant stellar population. The only viable fit found is for a fairly young (~ 6-9 Myr) burst suffering from a considerable extinction (A_v ~ 1.2-1.8 mag). We find that this object resembles strongly HII galaxies and intermediate redshift compact emission line galaxies, albeit suffering from a larger extinction than usual. We also discuss the possible contamination introduced by such R-band drop-out galaxies in searches for z > 5 galaxies.

We derive the luminosity function of high-redshift Lyman alpha emitting sources from a deep, blind, spectroscopic survey that utilized strong-lensing magnification by intermediate-redshift clusters of galaxies. We observed carefully selected regions near 9 clusters, consistent with magnification factors generally greater than 10 for the redshift range 4.5<z<6.7. Eleven emission-line candidates were located in the range 2.2<z<5.6 whose identification we justify as Lyman alpha, in most cases via further spectroscopic observations. The selection function we constructed for our survey takes into account our varying intrinsic Lyman alpha line sensitivity as a function of wavelength and sky position. By virtue of the strong magnification factor, we provide constraints on the Lyman alpha luminosity function to unprecedented limits of 10^40 erg/s, corresponding to a star-formation rate of 0.01 Msun/yr. Our cumulative z=5 Lyman alpha luminosity function is consistent with a power law form, n(>L) proportional to L^-1 over 10^41 to 10^42.5 erg/s. When combined with the results of other surveys, limited at higher luminosities, our results suggest evidence for the suppression of star formation in low-mass halos, as predicted in popular models of galaxy formation.

We present different methods used to identify high redshift (z>5) objects in the high-magnification regions of lensing galaxy clusters, taking advantage of very well constrained lensing models. The research procedures are explained and discussed. The detection of emission lines in the optical/NIR spectra, such as Lyman-alpha, allows us to determine the redshift of these sources. Thanks to the lensing magnification, it is possible to identify and to study more distant or intrinsically fainter objects with respect to standard field surveys.

This paper reviews the most recent developments related to the use of lensing clusters of galaxies as Gravitational Telescopes in deep Universe studies. We summarize the state of the art and the most recent results aiming at studying the physical properties of distant galaxies beyond the limits of conventional spectroscopy. The application of photometric redshift techniques in the context of gravitational lensing is emphasized for the study of both lensing structures and the background population of lensed galaxies. A presently ongoing search for the first building blocks of galaxies behind lensing clusters is presented and discussed.