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.

NoiseChisel is a program to detect very low signal-to-noise ratio (S/N) features with minimal assumptions on their morphology. It was introduced in 2015 and released within a collection of data analysis programs and libraries known as GNU Astronomy Utilities (Gnuastro). Over the last ten stable releases of Gnuastro, NoiseChisel has significantly improved: detecting even fainter signal, enabling better user control over its inner workings, and many bug fixes. The most important change may be that NoiseChisel's segmentation features have been moved into a new program called Segment. Another major change is the final growth strategy of its true detections, for example NoiseChisel is able to detect the outer wings of M51 down to S/N of 0.25, or 28.27 mag/arcsec2 on a single-exposure SDSS image (r-band). Segment is also able to detect the localized HII regions as "clumps" much more successfully. Finally, to orchestrate a controlled analysis, the concept of a "reproducible paper" is discussed: this paper itself is exactly reproducible (snapshot v4-0-g8505cfd).

The Hubble Ultra Deep field (HUDF) is the deepest region ever observed with the Hubble Space Telescope. With the main objective of unveiling the nature of galaxies up to $z \sim 7-8$, the observing and reduction strategy have focused on the properties of small and unresolved objects, rather than the outskirts of the largest objects, which are usually over-subtracted. We aim to create a new set of WFC3/IR mosaics of the HUDF using novel techniques to preserve the properties of the low surface brightness regions. We created ABYSS: a pipeline that optimises the estimate and modelling of low-level systematic effects to obtain a robust background subtraction. We have improved four key points in the reduction: 1) creation of new absolute sky flat fields, 2) extended persistence models, 3) dedicated sky background subtraction and 4) robust co-adding. The new mosaics successfully recover the low surface brightness structure removed on the previous HUDF published reductions. The amount of light recovered with a mean surface brightness dimmer than $\overline{\mu}=26$ mar arcsec$^{-2}$ is equivalent to a m=19 mag source when compared to the XDF and a m=20 mag compared to the HUDF12. We present a set of techniques to reduce ultra-deep images ($\mu>32.5$ mag arcsec$^{-2}$, $3\sigma$ in $10\times10$ arcsec boxes), that successfully allow to detect the low surface brightness structure of extended sources on ultra deep surveys. The developed procedures are applicable to HST, JWST, EUCLID and many other space and ground-based observatories. We will make the final ABYSS WFC3/IR HUDF mosaics publicly available at http://www.iac.es/proyecto/abyss/.

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.

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.

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 study the UV luminosity functions (LFs) at $z\sim 4$, $5$, $6,$ and $7$ based on the deep large-area optical images taken by the Hyper Suprime-Cam (HSC) Subaru strategic program (SSP). On the 100 deg$^2$ sky of the HSC SSP data available to date, we make enormous samples consisting of a total of 579,565 dropout candidates at $z\sim 4-7$ by the standard color selection technique, 358 out of which are spectroscopically confirmed by our follow-up spectroscopy and other studies. We obtain UV LFs at $z \sim 4-7$ that span a very wide UV luminosity range of $\sim 0.002 - 100 \, L_{\rm UV}^\ast$ ($-26 < M_{\rm UV} < -14$ mag) by combining LFs from our program and the ultra-deep Hubble Space Telescope legacy surveys. We derive three parameters of the best-fit Schechter function, $\phi^\ast$, $M_{\rm UV}^\ast$, and $\alpha$, of the UV LFs in the magnitude range where the AGN contribution is negligible, and find that $\alpha$ and $\phi^\ast$ decrease from $z\sim 4$ to $7$ with no significant evolution of $M_{\rm UV}^\ast$. Because our HSC SSP data bridge the LFs of galaxies and AGNs with great statistical accuracy, we carefully investigate the bright end of the galaxy UV LFs that are estimated by the subtraction of the AGN contribution either aided with spectroscopy or the best-fit AGN UV LFs. We find that the bright end of the galaxy UV LFs cannot be explained by the Schechter function fits at $> 2 \sigma$ significance, and require either double power-law functions or modified Schechter functions that consider a magnification bias due to gravitational lensing.

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.

In the coming era of massive surveys (e.g. LSST, SKA), the role of the database designers and the algorithms they choose to adopt becomes the decisive factor in scientific progress. Systems that allow/encourage users/scientists to be more creative with the reduction/analysis algorithms can greatly enhance scientific productivity. The separation/modularity of the detection processes and catalog production is one proposal for achieving `Reduction/analysis algorithms for large databases and vice versa' (a key theme for the 26th ADASS). With the new noise-based detection paradigm, non-parametric detection is now possible for astronomical objects to very low surface brightness limits. In our implementation, one software (NoiseChisel) is in charge of detection and another (MakeCatalog) is in charge of catalog production. This modularity has many advantages for pipeline developers, and more importantly, it empowers scientific curiosity and creativity.

A noise-based non-parametric technique for detecting nebulous objects, for example, irregular or clumpy galaxies, and their structure in noise is introduced. "Noise-based" and "non-parametric" imply that this technique imposes negligible constraints on the properties of the targets and that it employs no regression analysis or fittings. The sub-sky detection threshold is defined and initial detections are found, independently of the sky value. False detections are then estimated and removed using the ambient noise as a reference. This results in a purity level of 0.89 for the final detections as compared to 0.29 for SExtractor when a completeness of 1 is desired for a sample of extremely faint and diffuse mock galaxy profiles. The difference in the mean of the undetected pixels with the known background of mock images is decreased by 4.6 times depending on the diffuseness of the test profiles, quantifying the success in their detection. A non-parametric approach to defining substructure over a detected region is also introduced. NoiseChisel is our software implementation of this new technique. Contrary to the existing signal-based approach to detection, in its various implementations, signal-related parameters such as the image point spread function or known object shapes and models are irrelevant here. Such features make this technique very useful in astrophysical applications such as detection, photometry, or morphological analysis of nebulous objects buried in noise, for example, galaxies that do not generically have a known shape when imaged.

Despite expanding research activity in gravitational lens modeling, there is no particular software which is considered a standard. Much of the gravitational lens modeling software is written by individual investigators for their own use. Some gravitational lens modeling software is freely available for download but is widely variable with regard to ease of use and quality of documentation. This review of 13 software packages was undertaken to provide a single source of information. Gravitational lens models are classified as parametric models or non-parametric models, and can be further divided into research and educational software. Software used in research includes the GRAVLENS package (with both gravlens and lensmodel), Lenstool, LensPerfect, glafic, PixeLens, SimpLens, Lensview, and GRALE. In this review, GravLensHD, G-Lens, Gravitational Lensing, lens and MOWGLI are categorized as educational programs that are useful for demonstrating various aspects of lensing. Each of the 13 software packages is reviewed with regard to software features (installation, documentation, files provided, etc.) and lensing features (type of model, input data, output data, etc.) as well as a brief review of studies where they have been used. Recent studies have demonstrated the utility of strong gravitational lensing data for mass mapping, and suggest increased use of these techniques in the future. Coupled with the advent of greatly improved imaging, new approaches to modeling of strong gravitational lens systems are needed. This is the first systematic review of strong gravitational lens modeling software, providing investigators with a starting point for future software development to further advance gravitational lens modeling research.

Using the very deep Subaru images of the GOODS-N region, from the MOIRCS Deep Survey and images from the HST/ACS, we have measured the Luminosity Ratio (LR) of the outer to the central regions of massive (M>10^{10.5}M_{Sun}) galaxies at fixed radii in a single rest-frame for z<3.5 as a new approach to the problem of size evolution. We didn't observe any evolution in the median LR. Had a significant size growth occurred, the outer to central luminosity ratios would have demonstrated a corresponding increase with a decrease in redshift.

We present scaling relations between stellar-mass (Mstar) and the size of galaxies at 0.3 < z < 3 for half- (R_50) and 90 percent-light (R_90) radii, using a deep K-band selected catalogue taken with the Subaru Telescope and MOIRCS in the GOODS-North region. The logarithmic slope R \propto Mstar^{0.1-0.2} is independent of redshift in a wide mass range of Mstar ~ 10^8-10^11 Msun, irrespective of galaxy populations (star-forming, quiescent). The offset change is < 50 percent. Provided that optical light in the rest frame traces the stellar mass of galaxies, the universal relation demonstrates that the stellar mass was built up in galaxies over their cosmic histories in a similar manner on average irrelevant to galaxy mass. The small offset in each stellar mass bin from the universal relation shows weak size evolution at a given mass. There is a moderate increase of 30-50 percent for R_50 and R_90 for less massive galaxies (Mstar < 10^10 Msun) from z~3 to z~1, while the sizes remains unchanged or slightly decrease towards z~0.3. For massive galaxies (Mstar > 10^11 Msun), the evolution is ~70-80 % increase in R_90 from z~3 to z~0.3, though that in R_50 is weaker. The evolution of compactness factor, R_50/R_90, which becomes smaller at lower redshift, is suggestive of minor merging effect in the outer envelope of massive galaxies.