Search

We test whether Lyman alpha emitters (LAEs) and Lyman-break galaxies (LBGs) can be good tracers of high-z large-scale structures, using the Horizon Run 5 cosmological hydrodynamical simulation. We identify LAEs using the Ly{\alpha} emission line luminosity and its equivalent width, and LBGs using the broad-band magnitudes at z~2.4, 3.1, and 4.5. We first compare the spatial distributions of LAEs, LBGs, all galaxies, and dark matter around the filamentary structures defined by dark matter. The comparison shows that both LAEs and LBGs are more concentrated toward the dark matter filaments than dark matter. We also find an empirical fitting formula for the vertical density profile of filaments as a binomial power-law relation of the distance to the filaments. We then compare the spatial distributions of the samples around the filaments defined by themselves. LAEs and LBGs are again more concentrated toward their filaments than dark matter. We also find the overall consistency between filamentary structures defined by LAEs, LBGs, and dark matter, with the median spatial offsets that are smaller than the mean separation of the sample. These results support the idea that the LAEs and LBGs could be good tracers of large-scale structures of dark matter at high redshifts., Comment: 20 pages, 15 figures, 2 tables, accepted for publication in ApJ

DESI is a groundbreaking international project to observe more than 40 million quasars and galaxies over a 5-year period to create a 3D map of the sky. This map will enable us to probe multiple aspects of cosmology, from dark energy to neutrino mass. We are focusing here on one type of object observed by DESI, the Lyman Break Galaxies (LBGs). The aim is to use their spectra to determine whether they are indeed LBGs, and if so, to determine their distance from the Earth using a phenomenon called redshift. This will enable us to place these galaxies on the DESI 3D map. The aim is therefore to develop a convolutional neural network (CNN) inspired by QuasarNET (See arXiv:1808.09955), performing simultaneously a classification (LBG type or not) and a regression task (determine the redshift of the LBGs). Initially, data augmentation techniques such as shifting the spectra in wavelengths, adding noise to the spectra, or adding synthetic spectra were used to increase the model training dataset from 3,019 data to over 66,000. In a second phase, modifications to the QuasarNET architecture, notably through transfer learning and hyperparameter tuning with Bayesian optimization, boosted model performance. Gains of up to 26% were achieved on the Purity/Efficiency curve, which is used to evaluate model performance, particularly in areas with interesting redshifts, at low (around 2) and high (around 4) redshifts. The best model obtained an average score of 94%, compared with 75% for the initial model.

We study the environment of the z=6.33 ultraluminous quasar SDSS J010013.02+280225.8 (J0100) to understand its association with large-scale structure. Theoretical models propose high-redshift quasars as markers of galaxy overdensities residing in the most massive dark matter halos (DMHs) in the early universe. J0100 is an ultraluminous quasar with the most massive black hole known at z>6, suggesting a high likelihood of residing in a massive DMH. We present wide-field ($\sim$522 square arcminute) imaging in the r-, i-, and z-bands from the Large Binocular Camera on the Large Binocular Telescope, with Y- and J-band imaging from the Wide-field Infrared Camera on the Canada-France-Hawaii Telescope, centered on J0100. Applying color selections, we identify 23 objects as i-droput Lyman Break Galaxy (LBG) candidates in the J0100 field. We use the deep photometric catalog in the 1.27 square degree COSMOS field to calculate the density of LBGs in a blank field, and to estimate the selection completeness and purity. The observed surface density of LBG candidates in the J0100 field corresponds to a galaxy overdensity of $\delta$=4 (at 8.4$\sigma$). This large-scale overdensity suggests that the $\sim$ 22 square arcminute overdensity found by Kashino et al. using JWST data extends out to much larger scales. We calculate the angular auto-correlation function of the candidates and find a positive correlation on $\lesssim$ 10 arcminute scales as well as evidence of asymmetries in their spatial distribution, further suggesting a direct detection of large-scale structure around the ultra-luminous quasar J0100., Comment: 21 pages, 11 figures, 3 tables, to be published in The Astrophysical Journal (ApJ)

We use data from the JWST Public Release IMaging for Extragalactic Research (PRIMER) survey to measure the size scaling relations of 1668 rest-frame UV-bright Lyman-break galaxies (LBGs) at $z=3-5$ with stellar masses $\mathrm{log}_{10}(M_{\star}/M_{\odot}) > 9$. The sample was selected from seeing-dominated ground-based data, presenting an unbiased sampling of the morphology and size distributions of luminous sources. We fit S\'ersic profiles to eight NIRCam bands and also measure a non-parametric half-light radius. We find that the size distributions with both measurements are well-fit by a log-normal distribution at all redshifts, consistent with disk formation models where size is governed by host dark-matter halo angular momentum. We find a size-redshift evolution of $R_{e} = 3.51(1+z)^{-0.60\pm0.22}$ kpc, in agreement with JWST studies. When considering the typical (modal) size over $z=3-5$, we find little evolution with bright LBGs remaining compact at $R_{e}\simeq0.7-0.9$ kpc. Simultaneously, we find evidence for a build-up of large ($R_{e} > 2$ kpc) galaxies by $z=3$. We find some evidence for a negatively sloped size-mass relation at $z=5$ when S\'ersic profiles are used to fit the data in F200W. The intrinsic scatter in our size-mass relations increases at higher redshifts. Additionally, measurements probing the rest-UV (F200W) show larger intrinsic scatter than those probing the rest-optical (F356W). Finally, we leverage rest-UV and rest-optical photometry to show that disky galaxies are well established by $z=5$, but are beginning to undergo dissipative processes, such as mergers, by $z=3$. The agreement of our size-mass and size-luminosity relations with simulations provides tentative evidence for centrally concentrated star formation at high-redshift., Comment: 15 pages, 9 figures, 3 tables. Accepted to MNRAS. 2 figures added in appendix

We present a study of a sample of 45 spectroscopically confirmed, UV luminous galaxies at $z\sim 6$. They were selected as bright Lyman-break galaxies (LBGs) using deep multi-band optical images in more than 2 deg$^2$ of the sky, and subsequently identified via their strong Ly$\alpha$ emission. The majority of these LBGs span an absolute UV magnitude range from $-22.0$ to $-20.5$ mag with Ly$\alpha$ equivalent width (EW) between $\sim$10 and $\sim$200 \AA, representing the most luminous galaxies at $z\sim 6$ in terms of both UV continuum emission and Ly$\alpha$ line emission. We model the SEDs of 10 LBGs that have deep infrared observations from HST, JWST, and/or Spitzer, and find that they have a wide range of stellar masses and ages. They also have high star-formation rates ranging from a few tens to a few hundreds of Solar mass per year. Five of the LBGs have JWST or HST images and four of them show compact morphology in these images, including one that is roughly consistent with a point source, suggesting that UV luminous galaxies at this redshift are generally compact. The fraction of our photometrically selected LBGs with strong Ly$\alpha$ emission ($\mathrm{EW}>25$ \AA) is about $0.2$, which is consistent with previous results and supports a moderate evolution of the IGM opacity at the end of cosmic reionization. Using deep X-ray images, we do not find evidence of strong AGN activity in these galaxies, but our constraint is loose and we are not able to rule out the possibility of any weak AGN activity., Comment: 19 pages, 11 figures, Accepted for publication in ApJ

Using integrated spectra for two gravitationally lensed galaxies from the JWST TEMPLATES Early Release Science program, we analyze faint auroral lines, which provide direct measurements of the gas-phase chemical abundance. For the brighter galaxy, SGAS1723$+$34 ($z = 1.3293$), we detect the [OIII]$\lambda4363$, [SIII]$\lambda6312$, and [OII]$\lambda\lambda$7320,7330 auroral emission lines, and set an upper limit for the [NII]$\lambda5755$ line. For the second galaxy, SGAS1226$+$21 ($z = 2.925$), we do not detect any auroral lines, and report upper limits. With these measurements and upper limits, we constrain the electron temperatures in different ionization zones within both of these galaxies. For SGAS1723$+$34, where auroral lines are detected, we calculate direct oxygen and nitrogen abundances, finding an N/O ratio consistent with observations of nearby ($z\sim 0$) galaxies. These observations highlight the potent combination of JWST and gravitational lensing to measure faint emission lines in individual distant galaxies and to directly study the chemical abundance patterns in those galaxies., Comment: Submitted to the Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program. 15 pages, 6 figures. Comments welcome

We report for the first time a relationship between galaxy kinematics and net Lyman-alpha equivalent width (net Lya EW) in star forming galaxies during the epoch of peak cosmic star formation. Building on the previously reported broadband imaging segregation of Lya-emitting and Lya-absorbing Lyman break galaxies (LBGs) at z~2 (Paper I in this series) and previously at z~3, we use the Lya spectral type classification method to study the relationship between net Lya EW and nebular emission-line kinematics determined from IFU spectroscopy. We show that z~2 and z~3 LBGs segregate in colour-magnitude space according to their kinematic properties and Lyman-alpha spectral type, and conclude that LBGs with Lya dominant in absorption are almost exclusively rotation-dominated (presumably disc-like) systems, and LBGs with Lya dominant in emission characteristically have dispersion-dominated kinematics. We quantify the relationship between the strength of rotational dynamic support and net Lya EW, and demonstrate the consistency of our result with other properties that scale with net Lya EW and kinematics. Based on these findings, we suggest a method by which large samples of rotation- and dispersion-dominated galaxies might be selected using broadband imaging in as few as three filters and/or net Lya EW alone. Application of this method will enable an understanding of galaxy kinematic behaviour over large scales in datasets from current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs in redshift ranges from z~2-6 across many hundreds to thousands of Mpc. Finally, we speculate that the combination of our result linking net Lya EW and nebular emission-line kinematics with the known large-scale clustering behaviour of Lya-absorbing and Lya-emitting LBGs is evocative of a nascent morphology-density relation at z~2-3., Comment: Published in PASA (Version 2: Typos corrected and minor clarifying edits in S3.1 & S3.2.1)

High-redshift Lyman break galaxies (LBGs) are efficiently selected in deep images using as few as three broadband filters, and have been shown to have multiple intrinsic and small- to large-scale environmental properties related to Lyman-alpha. In this paper we demonstrate a statistical relationship between net Lyman-alpha equivalent width (net Lya EW) and the optical broadband photometric properties of LBGs at z~2. We show that LBGs with the strongest net Lya EW in absorption (aLBGs) and strongest net Lya EW in emission (eLBGs) separate into overlapping but discrete distributions in $(U_n-R)$ colour and $R$-band magnitude space, and use this segregation behaviour to determine photometric criteria by which sub-samples with a desired Lya spectral type can be selected using data from as few as three broadband optical filters. We propose application of our result to current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs across many hundreds to thousands of Mpc, and for which spectroscopic follow-up to obtain Lya spectral information is prohibitive. To this end, we use spectrophotometry of composite spectra derived from a sample of 798 LBGs divided into quartiles on the basis of net Lya EW to calculate criteria for the selection of Lya absorbing and Lya emitting populations of z~3 LBGs using $ugri$ broadband photometric data from the Vera Rubin Observatory Legacy Survey of Space and Time (LSST)., Comment: Published in Publications of the Astronomical Society of Australia 16 pages, 6 figures, 8 tables

We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-infrared (FIR) properties which reveal a significant fraction ($f_{\rm obs} = 0.4$-$0.7$) of obscured star formation, consistent with previous studies. From measurements of the rest-frame UV slope, we find that the brightest LBGs at these redshifts show bluer ($\beta \simeq -2.2$) colours than expected from an extrapolation of the colour-magnitude relation found at fainter magnitudes. Assuming a modified blackbody spectral-energy distribution (SED) in the FIR (with dust temperature of $T_{\rm d} = 46\,{\rm K}$ and $\beta_{\rm d} = 2.0$), we find that the REBELS sources are in agreement with the local ''Calzetti-like'' starburst Infrared-excess (IRX)-$\beta$ relation. By reanalysing the data available for 108 galaxies at $z \simeq 4$-$6$ from the ALPINE ALMA large program using a consistent methodology and assumed FIR SED, we show that from $z \simeq 4$-$8$, massive galaxies selected in the rest-frame UV have no appreciable evolution in their derived IRX-$\beta$ relation. When comparing the IRX-$M_{\star}$ relation derived from the combined ALPINE and REBELS sample to relations established at $z < 4$, we find a deficit in the IRX, indicating that at $z > 4$ the proportion of obscured star formation is lower by a factor of $\gtrsim 3$ at a given a $M_{\star}$. Our IRX-$\beta$ results are in good agreement with the high-redshift predictions of simulations and semi-analytic models for $z \simeq 7$ galaxies with similar stellar masses and SFRs., Comment: 18 pages, 10 figures, 2 tables (plus 1 figure and 2 tables in the appendix). Submitted to MNRAS, comments welcome

We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$\alpha$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MIRI playing a crucial role in robustly assessing the LAE's stellar mass and ages. Our LAEs are low-mass objects (log$_{10}$(M$_\star$[M$_\odot$]) ~ 7.5), with little or no dust extinction (E(B - V) ~ 0.1) and a blue UV continuum slope ($\beta$ ~ -2.2). While 75% of our LAEs are young (< 100 Myr), the remaining 25% have significantly older stellar populations (> 100 Myr). These old LAEs are statistically more massive, less extinct and have lower specific star formation rate (sSFR) compared to young LAEs. Besides, they populate the M$_\star$ - SFR plane along the main-sequence (MS) of star-forming galaxies, while young LAEs populate the starburst region. The comparison between the LAEs properties to those of a stellar-mass matched sample of LBGs shows no statistical difference between these objects, except for the LBGs redder UV continuum slope and marginally larger E(B - V) values. Interestingly, 48% of the LBGs have ages < 10 Myr and are classified as starbursts, but lack detectable Ly$\alpha$ emission. This is likely due to HI resonant scattering and/or selective dust extinction. Overall, we find that JWST observations are crucial in determining the properties of LAEs and shedding light on the properties and similarities between LAEs and LBGs., Comment: 22 pages, 15 figures, 5 tables. Accepted for publication in APJ

We report the identification of Lyman Break Galaxy (LBG) candidates around the most luminous Hot Dust-Obscured Galaxy (Hot DOG) known, WISE J224607.56$-$052634.9 (W2246$-$0526) at $z=4.601$, using deep \textit{r}-, \textit{i}-, and \textit{z}-band imaging from the Gemini Multi-Object Spectrograph South (GMOS-S). We use the surface density of LBGs to probe the Mpc-scale environment of W2246$-$0526 to characterize its richness and evolutionary state. We identify LBG candidates in the vicinity of W2246$-$0526 using the selection criteria developed by \cite{2004VOuchi} and \cite{2006Yoshida} in the Subaru Deep Field and in the Subaru XMM-Newton Deep Field, slightly modified to account for the difference between the filters used, and we find 37 and 55 LBG candidates, respectively. Matching to the $z$-band depths of those studies, this corresponds to $\delta = 5.8^{+2.4}_{-1.9}$ times the surface density of LBGs expected in the field. Interestingly, the Hot DOG itself, as well as a confirmed neighbor, do not satisfy either LBG selection criteria, suggesting we may be missing a large number of companion galaxies. Our analysis shows that we are most likely only finding those with higher-than-average IGM optical depth or moderately high dust obscuration. The number density of LBG candidates is not concentrated around W2246$-$0526, suggesting either an early evolutionary stage for the proto-cluster or that the Hot DOG may not be the most massive galaxy, or that the Hot DOG may be affecting the IGM transparency in its vicinity. The overdensity around W2246$-$0526 is comparable to overdensities found around other Hot DOGs and is somewhat higher than typically found for radio galaxies and luminous quasars at a similar redshift., Comment: 20 pages, 15 figures. The main results are in Figures 9 and 12. Accepted for publication in A&A

We present our new Atacama Large Millimeter/submillimeter Array (ALMA) observations targeting CO(6-5) emission from three luminous Lyman break galaxies (LBGs) at $z_{\rm spec} = 6.0293$-$6.2037$ found in the Subaru/Hyper Suprime-Cam survey, whose [OIII]$88\mu$m and [CII]$158\mu$m emission have been detected with ALMA. We find a marginal detection of the CO(6-5) line from one of our LBGs, J0235-0532, at the $\simeq 4 \sigma$ significance level and obtain upper limits for the other two LBGs, J1211-0118 and J0217-0208. Our $z=6$ luminous LBGs are consistent with the previously found correlation between the CO luminosity and the infrared luminosity. The unique ensemble of the multiple far-infrared emission lines and underlying continuum fed to a photodissociation region model reveal that J0235-0532 has a relatively high hydrogen nucleus density that is comparable to those of low-$z$ (U)LIRGs, quasars, and Galactic star-forming regions with high $n_{\rm H}$ values, while the other two LBGs have lower $n_{\rm H}$ consistent with local star-forming galaxies. By carefully taking account of various uncertainties, we obtain total gas mass and gas surface density constraints from their CO luminosity measurements. We find that J0235-0532 locates below the Kennicutt-Schmidt (KS) relation, comparable to the previously CO(2-1) detected $z=5.7$ LBG, HZ10. Combined with previous results for dusty starbursts at similar redshifts, the KS relation at $z=5$-$6$ is on average consistent with the local one., Comment: 33 pages, 12 figures, accepted for publication in ApJ

Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identified 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing stacking analyses, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift., Comment: 26 pages, 14 figures, 8 tables. Accepted for publication in A&A

We present the first search for $z\geqslant7$, continuum-confirmed sources with NIRISS/WFS spectroscopy over the Abell 2744 Frontier Fields cluster, as part of the GLASS-JWST ERS survey. With $\sim15$ hrs of pre-imaging and multi-angle grism exposures in the F115W, F150W, and F200W filters, we describe the general data handling (i.e., reduction, cleaning, modeling, and extraction processes) and analysis for the GLASS-JWST survey. We showcase the power of JWST to peer deep into reionization, when most intergalactic hydrogen is neutral, by confirming two galaxies at $z=8.04\pm0.15$ and $z=7.90\pm0.13$ by means of their Lyman breaks. Fainter continuum spectra are observed in both the F150W and F200W bands, indicative of blue ($-1.69$ and $-1.33$) UV slopes and moderately-bright absolute magnitudes ($-20.37$ and $-19.68$ mag). We do not detect strong Ly$\alpha$ in either galaxy, but do observe tentative ($\sim2.7-3.8\sigma$) HeII$\lambda$1640 A, OIII]$\lambda\lambda$1661,1666 A, and NIII]$\lambda\lambda$1747,1749 A line emission in one, suggestive of low metallicity, star-forming systems with possible non-thermal contributions. These novel observations provide a first look at the extraordinary potential of JWST/NIRISS for confirming representative samples of bright $z\geqslant7$ sources in the absence of strong emission lines, and gain unprecedented insight into their contributions towards cosmic reionization., Comment: Accepted for publication in ApJL on 1st September 2022, after minor modifications based on the referee report. 10 pages, 3 figures, 1 table

We measure, for the first time, the median equivalent width (EW) of H$\alpha$+[N II] in star-forming galaxies at $z\sim8$. Our estimate leverages the unique photometric depth of the Spitzer/IRAC $5.8\mu$m-band mosaics (probing $\approx 5500 - 7100$ A at $z\sim8$) of the GOODS Reionization Era Wide Area Treasury from Spitzer (GREATS) program. We median stacked the stamps of $102$ Lyman-break galaxies in the $3.6, 4.5, 5.8$ and $8.0\mu$m bands, after carefully removing potential contamination from neighbouring sources. We infer an extreme rest-frame EW$_0$(H$\alpha$+[N II])$=2328^{+1326}_{-1127}$ A from the measured red $[3.6]-[5.8]=0.82\pm0.27$ mag, consistent with young ($\lesssim10^7$ yr) average stellar population ages at $z\sim8$. This implies an ionizing photon production efficiency of $\log(\xi_{\mathrm{ion},0}/\mathrm{erg\ Hz}^{-1})=25.97^{+0.18}_{-0.28}$. Such a high value for photo production, similar to the highest values found at $z\lesssim4$, indicates that only modest escape fractions $f_\mathrm{esc}\lesssim0.3$ (at $2\sigma$) are sufficient for galaxies brighter than $M_\mathrm{UV}<-18$ mag to re-ionize the neutral Hydrogen at $z\sim8$. This requirement is relaxed even more to $f_\mathrm{esc}\le 0.1$ when considering galaxies brighter than $M_\mathrm{UV}\approx -13$ mag, consistent with recent luminosity functions and as typically assumed in studies addressing re-ionization. These exceptional results clearly indicate that galaxies can be the dominant source of reionizing photons, and provide us with an exciting glimpse into what we might soon learn about the early universe, and particularly about the Reionization Epoch, from forthcoming JWST/MIRI and NIRCam programs., Comment: Submitted to ApJ. 12 pages, 5 figures (main) + appendix. Comments welcome

We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $\tau_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (\psi_{\rm UV}) = -0.05 ~[log (\psi)]^{2} + 0.86 ~log(\psi) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $\tau_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions., Comment: Accepted to MNRAS

We present a Bayesian inference on the neutral hydrogen fraction of the intergalactic medium (IGM), $x_{\textrm{HI}}$, at $z \sim$ 6-8 using the properties of Lyman break galaxies during the Epoch of Reionization. We use large samples of LBG candidates at $5.5 \leq z \leq 8.2$ with spectroscopy from Keck/DEIMOS and Keck/MOSFIRE. For each galaxy, we incorporate either the Lyman-$\alpha$ equivalent width (EW) for detections or the EW limit spectrum for nondetections to parameterize the EW distribution at various ultraviolet brightnesses for a given redshift. Using our reference sample of galaxy candidates from the ionized universe at $z$ $\sim$ 6.0, we are able to infer $x_{\textrm{HI}}$ at two redshifts: $z$ $\sim$ 6.7 and $z$ $\sim$ 7.6. This work includes intrinsically faint, gravitationally lensed galaxies at $z$ $\sim$ 6.0 in order to constrain the intrinsic faint-end Ly$\alpha$ EW distribution and provide a comparable population of galaxies to counterparts in our sample that are at higher redshift. The inclusion of faint galaxy candidates, in addition to a more sophisticated modelling framework, allows us to better isolate effects of the interstellar medium and circumgalactic medium on the observed Lyman-$\alpha$ distribution from those of the IGM. We infer an upper limit of $x_{\textrm{HI}}$ $\leq$ 0.25 at $z$ = 6.7 $\pm$ 0.2 and a neutral fraction of $x_{\textrm{HI}}$ = $0.83^{+0.08}_{-0.11}$ at $z$ = 7.6$\pm$ 0.6, both within 1$\sigma$ uncertainty, results which favor a moderately late and rapid reionization., Comment: 11 pages, 7 figures

Some studies suggest that the dust temperatures ($T_\mathrm{d}$) in high-redshift ($z\gtrsim 5$) Lyman break galaxies (LBGs) are high. However, possible observational bias in $T_\mathrm{d}$ is yet to be understood. Thus, we perform a simple test using random realizations of LBGs with various stellar masses, dust temperatures, and dust-to-stellar mass ratios, and examine how the sample detected by ALMA is biased in terms of $T_\mathrm{d}$. We show that ALMA tends to miss high-$T_\mathrm{d}$ objects even at total dust luminosity. LBGs are, however, basically selected by the stellar UV luminosity. The dust-temperature bias in a UV-selected sample is complicated because of the competing effects between high $T_\mathrm{d}$ and low dust abundance. For ALMA Band 6, there is no tendency of high-$T_\mathrm{d}$ LBGs being more easily detected in our experiment. Thus, we suggest that the observed trend of high $T_\mathrm{d}$ in $z\gtrsim 5$ LBGs is real. We also propose that the 450 $\mu$m band is useful in further clarifying the dust temperatures. To overcome the current shallowness of 450 $\mu$m observations, we examine a future Antarctic 30-m class telescope with a suitable atmospheric condition for wavelengths $\lesssim 450~\mu$m, where the detection is not confusion-limited. We find that, with this telescope, an $L_\mathrm{IR}$-selected sample with $\log(L_\mathrm{IR}/\mathrm{L}_{\odot})>11$ is constructed for $z\gtrsim 5$, and detection in the intermediate-$M_\star$ (stellar mass) range [$9<\log (M_\star /M_{\odot})<9.5$] is much improved, especially at high $T_\mathrm{d}$., Comment: 11 pages, 6 figures, accepted for publication in MNRAS

We present deep ALMA dust continuum observations for a sample of luminous ($M_{\rm UV} < -22$) star-forming galaxies at $z \simeq 7$. We detect five of the six sources in the far-infrared (FIR), providing key constraints on the obscured star-formation rate (SFR) and the infrared-excess-$\beta$ (IRX-$\beta$) relation without the need for stacking. Despite the galaxies showing blue rest-frame UV slopes ($\beta \simeq -2$) we find that 35-75 percent of the total SFR is obscured. We find the IRX-$\beta$ relation derived for these $z \simeq 7$ sources is consistent with that found for local star-burst galaxies. Using our relatively high-resolution (FWHM $\simeq 0.7\,{\rm arcsec}$) observations we identify a diversity of dust morphologies in the sample. We find both compact emission that appears offset relative to the unobscured components and extended dust emission that is co-spatial with the rest-frame UV light. In the majority of the sources we detect strong rest-frame UV colour gradients (with up to $\Delta \beta \simeq 0.7$-$1.4$) as probed by the multi-band UltraVISTA ground-based data. The observed redder colours are spatially correlated with the location of the FIR detection. Our results show that even in bright Lyman-break galaxies at $z \simeq 7$ the peak of the star-formation is typically hosted by the fainter, redder, regions in the rest-frame UV, which have an obscured fraction of $f_{\rm obs} \ge 0.8$. As well as demonstrating the importance of dust obscured star-formation within the Epoch of Reionization, these observations provide an exciting taster of the rich spatially resolved datasets that will be obtained from JWST and high-resolution ALMA follow-up at these redshifts., Comment: 14 pages, 6 figures, 3 tables. Updated to MNRAS accepted version

We present two bright galaxy candidates at z~12-13 identified in our H-dropout Lyman break selection with 2.3 deg2 near-infrared deep imaging data. These galaxy candidates, selected after careful screening of foreground interlopers, have spectral energy distributions showing a sharp discontinuity around 1.7 um, a flat continuum at 2-5 um, and non-detections at <1.2 um in the available photometric datasets, all of which are consistent with z>12 galaxy. An ALMA program targeting one of the candidates shows a tentative 4sigma [OIII]88um line at z=13.27, in agreement with its photometric redshift estimate. The number density of the z~12-13 candidates is comparable to that of bright z~10 galaxies, and is consistent with a recently proposed double power-law luminosity function rather than the Schechter function, indicating little evolution in the abundance of bright galaxies from z~4 to 13. Comparisons with theoretical models show that the models cannot reproduce the bright end of rest-frame ultraviolet luminosity functions at z~10-13. Combined with recent studies reporting similarly bright galaxies at z~9-11 and mature stellar populations at z~6-9, our results indicate the existence of a number of star-forming galaxies at z>10, which will be detected with upcoming space missions such as James Webb Space Telescope, Nancy Grace Roman Space Telescope, and GREX-PLUS., Comment: 20 pages, 13 figures, accepted for publication in ApJ

We make use of ALMA continuum observations of $15$ luminous Lyman-break galaxies at $z$$\sim$$7$$-$$8$ to probe their dust-obscured star-formation. These observations are sensitive enough to probe to obscured SFRs of $20$ $M_{\odot}$$/$$yr$ ($3\sigma$). Six of the targeted galaxies show significant ($\geq$$3$$\sigma$) dust continuum detections, more than doubling the number of known dust-detected galaxies at $z$$>$$6.5$. Their IR luminosities range from $2.7$$\times$$10^{11}$ $L_{\odot}$ to $1.1$$\times$$10^{12}$ $L_{\odot}$, equivalent to obscured SFRs of $20$ to $105$ $M_{\odot}$$/$$yr$. We use our results to quantify the correlation of the infrared excess IRX on the UV-continuum slope $\beta_{UV}$ and stellar mass. Our results are most consistent with an SMC attenuation curve for intrinsic $UV$-slopes $\beta_{UV,intr}$ of $-2.63$ and most consistent with an attenuation curve in-between SMC and Calzetti for $\beta_{UV,intr}$ slopes of $-2.23$, assuming a dust temperature $T_d$ of $50$ K. Our fiducial IRX-stellar mass results at $z$$\sim$$7$$-$$8$ are consistent with marginal evolution from $z$$\sim$$0$. We then show how both results depend on $T_d$. For our six dust-detected sources, we estimate their dust masses and find that they are consistent with dust production from SNe if the dust destruction is low ($<$$90$%). Finally we determine the contribution of dust-obscured star formation to the star formation rate density for $UV$ luminous ($<$$-$$21.5$ mag: $\gtrsim$$1.7$$L_{UV} ^*$) $z$$\sim$$7$$-$$8$ galaxies, finding that the total SFR density at $z$$\sim$$7$ and $z$$\sim$$8$ from bright galaxies is $0.18_{-0.10}^{+0.08}$ dex and $0.20_{-0.09}^{+0.05}$ dex higher, respectively, i.e. $\sim$$\frac{1}{3}$ of the star formation in $\gtrsim$$1.7$$L_{UV} ^*$ galaxies at $z$$\sim$$7$$-$$8$ is obscured by dust., Comment: 21 pages, 12 figures, 2 tables. Submitted to ApJ

We show that the James Webb Space Telescope will be able to detect dwarf satellites of high-$z$ Lyman Break Galaxies (LBGs). To this aim, we use cosmological simulations following the evolution of a typical $M_\star\simeq10^{10}\rm M_\odot$ LBG up to $z\simeq6$, and analyse the observational properties of its five satellite dwarf galaxies ($10^7{\rm M_\odot}

Several evidences indicate that Lyman Break Galaxies (LBG) in the Epoch of Reionization (redshift $z>6$) might host massive black holes (MBH). We address this question by using a merger-tree model combined with tight constraints from the 7 Ms Chandra survey, and the known high-$z$ super-MBH population. We find that a typical LBG with $M_{\rm UV}=-22$ residing in a $M_h\approx 10^{12} M_\odot$ halo at $z=6$ host a MBH with mass $M_\bullet \approx 2\times 10^8 M_\odot$. Depending on the fraction, $f_{\rm seed}$, of early halos planted with a direct collapse black hole seed ($M_{\rm seed}=10^5 M_\odot$), the model suggests two possible scenarios: (a) if $f_{\rm seed}=1$, MBH in LBGs mostly grow by merging, and must accrete at a low ($\lambda_E\simeq 10^{-3}$) Eddington ratio not to exceed the experimental X-ray luminosity upper bound $L_X^* = 10^{42.5} {\rm erg\, s}^{-1}$; (b) if $f_{\rm seed}=0.05$ accretion dominates ($\lambda_E\simeq 0.22$), and MBH emission in LBGs must be heavily obscured. In both scenarios the UV luminosity function is largely dominated by stellar emission up to very bright mag, $M_{\rm UV} > -23$, with BH emission playing a subdominant role. Scenario (a) poses extremely challenging, and possibly unphysical, requirements on DCBH formation. Scenario (b) entails testable implications on the physical properties of LBGs involving the FIR luminosity, emission lines, and presence of outflows.

We examine the absorption line spectra of a sample of 31 luminous (M_UV=-23) Lyman break galaxies at redshift z~6 using data taken with the FOCAS and OSIRIS spectrographs on the Subaru and GTC telescopes. For two of these sources we present longer exposure data taken at higher spectral resolution from ESO's X-shooter spectrograph. Using these data, we demonstrate the practicality of stacking our lower resolution data to measure the depth of various interstellar and stellar absorption lines to probe the covering fraction of low ionization gas and the gas-phase and stellar metallicities near the end of the era of cosmic reionization. From maximum absorption line depths of SiII1260 and CII1334, we infer a mean covering fraction of >0.85+/-0.16 for our sample. This is larger than that determined using similar methods for lower luminosity galaxies at slightly lower redshifts, suggesting that the most luminous galaxies appear to have a lower escape fraction than fainter galaxies, and therefore may not play a prominent role in concluding reionization. Using various interstellar absorption lines we deduce gas-phase metallicities close to solar indicative of substantial early enrichment. Using selected stellar absorption lines, we model our spectra with a range of metallicities using techniques successfully employed at lower redshift and deduce a stellar metallicity of 0.4 +0.3/-0.1 solar, consistent with the stellar mass - stellar metallicity relation recently found at z~3-5. We discuss the implications of these metallicity estimates for the typical ages of our luminous galaxies and conclude our results imply initial star formation at redshifts z~10, consistent with independent analyses of earlier objects., Comment: 13 pages, 10 figures, 3 tables, Accepted for Publication in ApJ

We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves., Comment: 17 pages, 12 figures. Resubmitted to MNRAS after referee report

The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works., Comment: Accepted for publication in A&A, 28 pages and 20 figures. Full versions of Tables C.1, C.2, C.3, C.4, and C.5 are available for download at http://astro.puc.cl/~rcarvaja/docs/tables_ALMA_FF_V_RCP.tar.gz

Characterizing high-z quasar environments is key to understanding the co-evolution of quasars and the surrounding galaxies. To restrict their global picture, we statistically examine the g-dropout galaxy overdensity distribution around 570 faint quasar candidates at z ~ 4, based on the Hyper Suprime-Cam Subaru Strategic Program survey. We compare the overdensity significances of g-dropout galaxies around the quasars with those around g-dropout galaxies, and find no significant difference between their distributions. A total of 4 (22) out of the 570 faint quasars, 0.7_{-0.4}^{+0.4} (3.9_{-0.8}^{+0.8}) %, are found to be associated with the > 4 sigma overdense regions within an angular separation of 1.8 (3.0) arcmin, which is the typical size of protoclusters at this epoch. This is similar to the fraction of g-dropout galaxies associated with the > 4 sigma overdense regions. This result is consistent with our previous work that 1.3_{-0.9}^{+0.9} % and 2.0_{-1.1}^{+1.1} % of luminous quasars detected in the Sloan Digital Sky Survey exist in the > 4 sigma overdense regions within 1.8 and 3.0 arcmin separations, respectively. Therefore, we suggest that the galaxy number densities around quasars are independent of their luminosity, and most quasars do not preferentially appear in the richest protocluster regions at z ~ 4. The lack of an apparent positive correlation between the quasars and the protoclusters implies that: i) the gas-rich major merger rate is relatively low in the protocluster regions, ii) most high-z quasars may appear through secular processes, or iii) some dust-obscured quasars exist in the protocluster regions., Comment: 12 pages, 4 figures

This work explores from a statistical point of view the rest-frame Far-ultraviolet (FUV) to Far-infrared (FIR) emission of a population of Lyman break galaxies (LBGs) at $z\sim3$ that cannot be individually detected from current FIR observations. We perform a stacking analysis over a sample of $\sim$17000 LBGs at redshift $2.5 10.25$)., Comment: 20 pages, 12 figures, Accepted for publication in A&A

We analyse a unique sample of 517 bright ($L>L^{*}$) LBGs at redshift z$\sim$3 in order to characterise the distribution of their UV slopes $\beta$ and infer their dust extinction under standard assumptions. We exploited multi-band observations over 750 arcmin$^2$ of the COSMOS field that were acquired with three different ground-based facilities: the Large Binocular Camera (LBC) on the Large Binocular Telescope (LBT), the Suprime-Cam on the SUBARU telescope, and the VIRCAM on the VISTA telescope (ULTRAVISTA DR2). Our multi-band photometric catalogue is based on a new method that is designed to maximise the signal-to-noise ratio in the estimate of accurate galaxy colours from images with different point spread functions (PSF). We adopted an improved selection criterion based on deep Y-band data to isolate a sample of galaxies at $z\sim 3 $ to minimise selection biases. We measured the UV slopes ($\beta$) of the objects in our sample and then recovered the intrinsic probability density function of $\beta$ values (PDF($\beta$)), taking into account the effect of observational uncertainties through detailed simulations. The galaxies in our sample are characterised by mildly red UV slopes with $<\beta>\simeq -1.70$ throughout the enitre luminosity range that is probed by our data ($-24\lesssim M_{1600}\lesssim -21$). The resulting dust-corrected star formation rate density (SFRD) is $log(SFRD)\simeq-1.6 M_{\odot}/yr/Mpc^{3}$, corresponding to a contribution of about 25% to the total SFRD at z$\sim$3 under standard assumptions. Ultra-bright LBGs at $z \sim 3$ match the known trends, with UV slopes being redder at decreasing redshifts, and brighter galaxies being more highly dust extinct and more frequently star-forming than fainter galaxies. [abridged], Comment: Matched to journal version. 11 pages, 13 figures, Astronomy & Astrophysics in press

Quasars (QSOs) hosting supermassive black holes are believed to reside in massive halos harboring galaxy overdensities. However, many observations revealed average or low galaxy densities around $z\gtrsim6$ QSOs. This could be partly because they measured galaxy densities in only tens of arcmin$^2$ around QSOs and might have overlooked potential larger scale galaxy overdensities. Some previous studies also observed only Lyman break galaxies (LBGs, massive older galaxies) and missed low mass young galaxies like Ly$\alpha$ emitters (LAEs) around QSOs. Here we present observations of LAE and LBG candidates in $\sim700$ arcmin$^2$ around a $z=6.61$ luminous QSO using Subaru Telescope Suprime-Cam with narrow/broadbands. We compare their sky distributions, number densities and angular correlation functions with those of LAEs/LBGs detected in the same manner and comparable data quality in our control blank field. In the QSO field, LAEs and LBGs are clustering in 4-20 comoving Mpc angular scales, but LAEs show mostly underdensity over the field while LBGs are forming $30\times60$ comoving Mpc$^2$ large scale structure containing 3-$7\sigma$ high density clumps. The highest density clump includes a bright (23.78 mag in the narrowband) extended ($\gtrsim 16$ kpc) Ly$\alpha$ blob candidate, indicative of a dense environment. The QSO could be part of the structure but is not located exactly at any of the high density peaks. Near the QSO, LAEs show underdensity while LBGs average to $4\sigma$ excess densities compared to the control field. If these environments reflect halo mass, the QSO may not be in the most massive halo, but still in a moderately massive one., Comment: 35 pages, 15 figures, 2 tables, accepted for publication in The Astrophysical Journal

We present Atacama Large Millimeter/Submillimeter Array observations of the rest-frame far-infrared (FIR) dust continuum emission of six bright Lyman-break galaxies (LBGs) at $z \simeq 7$. One LBG is detected ($5.2\sigma$ at peak emission), while the others remain individually undetected at the $3\sigma$ level. The average FIR luminosity of the sample is found to be $L_{\rm FIR} \simeq 2 \times 10^{11}\,{\rm L}_{\odot}$, corresponding to an obscured star-formation rate (SFR) that is comparable to that inferred from the unobscured UV emission. In comparison to the infrared excess (IRX$\,=L_{\rm FIR}/L_{\rm UV}$)-$\beta$ relation, our results are consistent with a Calzetti-like attenuation law (assuming a dust temperature of T = 40-50 K). We find a physical offset of 3 kpc between the dust continuum emission and the rest-frame UV light probed by Hubble Space Telescope imaging for galaxy ID65666 at $z = 7.17^{+0.09}_{-0.06}$. The offset is suggestive of an inhomogeneous dust distribution, where 75% of the total star formation activity (SFR$ \,\simeq 70\,{\rm M}_{\odot}/{\rm yr}$) of the galaxy is completely obscured. Our results provide direct evidence that dust obscuration plays a key role in shaping the bright-end of the observed rest-frame UV luminosity function at $z \simeq 7$, in agreement with cosmological galaxy formation simulations. The existence of a heavily-obscured component of galaxy ID65666 indicates that dusty star-forming regions, or even entire galaxies, that are "UV-dark" are significant even in the $z \simeq 7$ galaxy population., Comment: 15 pages, 8 figures, 3 tables. Updated to match MNRAS accepted version

We use a sample of 4178 Lyman break galaxies (LBGs) at z = 3, 4 and 5 in the UKIRT Infrared Deep Sky Survey (UKIDSS) Ultra Deep Survey (UDS) field to investigate the relationship between the observed slope of the stellar continuum emission in the ultraviolet, {\beta}, and the thermal dust emission, as quantified via the so-called 'infrared excess' (IRX = LIR/LUV). Through a stacking analysis we directly measure the 850-{\mu}m flux density of LBGs in our deep (0.9mJy) James Clerk Maxwell Telescope (JCMT) SCUBA-2 850-{\mu}m map, as well as deep public Herschel/SPIRE 250-, 350- and 500-{\mu}m imaging. We establish functional forms for the IRX-{\beta} relation to z ~ 5, confirming that there is no significant redshift evolution of the relation and that the resulting average IRX-{\beta} curve is consistent with a Calzetti-like attenuation law. We compare our results with recent work in the literature, finding that discrepancies in the slope of the IRX-{\beta} relation are driven by biases in the methodology used to determine the ultraviolet slopes. Consistent results are found when IRX-{\beta} is evaluated by stacking in bins of stellar mass, M, and we argue that the near-linear IRX-M relationship is a better proxy for correcting observed UV luminosities to total star formation rates, provided an accurate handle on M can be had, and also gives clues as to the physical driver of the role of dust-obscured star formation in high-redshift galaxies., Comment: 12 pages, 10 figures, submitted to MNRAS

We present the first comprehensive evolutionary analysis of the rest-frame UV spectroscopic properties of star-forming galaxies at z ~ 2-4. We match samples at different redshifts in UV luminosity and stellar mass, and perform systematic measurements of spectral features and stellar population modeling. By creating composite spectra grouped according to Ly$\alpha$ equivalent width (EW), and various galaxy properties, we study the evolutionary trends among Ly$\alpha$, low- and high-ionization interstellar (LIS and HIS) absorption features, and integrated galaxy properties. We also examine the redshift evolution of Ly$\alpha$ and LIS absorption kinematics, and fine-structure emission EWs. The connections among the strengths of Ly$\alpha$, LIS lines, and dust extinction are redshift-independent, as is the decoupling of Ly$\alpha$ and HIS line strengths, and the bulk outflow kinematics as traced by LIS lines. Stronger Ly$\alpha$ emission is observed at higher redshift at fixed UV luminosity, stellar mass, SFR, and age. Much of this variation in average Ly$\alpha$ strength with redshift, and the variation in Ly$\alpha$ strength at fixed redshift, can be explained in terms of variations in neutral gas covering fraction and/or dust content in the ISM and CGM. However, based on the connection between Ly$\alpha$ and CIII] emission strengths, we additionally find evidence for variations in the intrinsic production rate of Ly$\alpha$ photons at the highest Ly$\alpha$ EWs. The challenge now is to understand the observed evolution in neutral gas covering fraction and dust extinction within a coherent model for galaxy formation, and make robust predictions for the escape of ionizing radiation at z > 6., Comment: 27 pages, 14 figures, ApJ, accepted

We combined deep U-band imaging from the KPNO-4m/MOSAIC camera with very deep multi-waveband data from the optical to infrared, to select Lyman Break Galaxies (LBGs) at z~3 using U-V and V-R colors in the Subaru Deep Field. With the resulting sample of 5161 LBGs, we construct the UV luminosity function down to $M_{UV} = -18$ and find a steep faint-end slope of $\alpha=-1.78 \pm 0.05$. We analyze rest-frame UV-to-IR spectral energy distributions generated from the median optical photometry and photometry on median-stacked IR images. In the stacks of faint LBGs, we find a background depression centered on the galaxy. This deficit results from the systematic difficulty of SExtractor in finding faint galaxies in regions with higher-than-average surface densities of foreground galaxies. We corrected our stacked magnitudes for this. Best-fit stellar population templates for the stacked LBG SEDs indicate stellar masses and star-formation rates of log M*/Msun = 10 and 50 M$_\odot$/yr at i' = 24, down to log M*/Msun = 8 and = 3 M$_\odot$/yr at i' = 27. For the faint stacked LBGs there is a 1-mag excess over the expected stellar continuum in the K-band, which we attribute to redshifted [OIII]4959+5007 and H$\beta$ lines. Their implied equivalent widths increase with decreasing mass, reaching $\rm{EW_0([O III]4959,5007+H\beta)}$ =1500A in the faintest bin. Such strong [OIII] emission is seen only in a miniscule fraction of the most extreme local emission-line galaxies, but it probably universal in the faint galaxies that reionized the universe. Finally, we analyze clustering by computing the angular correlation function and performing halo occupation distribution (HOD) analysis. We find a mean dark halo mass of log(Mhalo/h) Msun = 11.29$\pm 0.12$ for the full sample of LBGs, and log(Mhalo/h) Msun = 11.49$\pm 0.1$ for the brightest half., Comment: scheduled for publication in The Astrophysical Journal, 2017, 21 pages in preprint format (19 journal pages)

In the standard picture of structure formation, the first massive galaxies are expected to form at the highest peaks of the density field, which constitute the cores of massive proto-clusters. Luminous quasars (QSOs) at z~4 are the most strongly clustered population known, and should thus reside in massive dark matter halos surrounded by large overdensities of galaxies, implying a strong QSO-galaxy cross-correlation function. We observed six z~4 QSO fields with VLT/FORS exploiting a novel set of narrow band filters custom designed to select Lyman Break Galaxies (LBGs) in a thin redshift slice of Delta_z~0.3, mitigating the projection effects that have limited the sensitivity of previous searches for galaxies around z>~4 QSOs. We find that LBGs are strongly clustered around QSOs, and present the first measurement of the QSO-LBG cross-correlation function at z~4, on scales of 0.1<~R<~9 Mpc/h (comoving). Assuming a power law form for the cross-correlation function xi=(r/r0_QG)^gamma, we measure r0_QG=8.83^{+1.39}_{-1.51} Mpc/h for a fixed slope of gamma=2.0. This result is in agreement with the expected cross-correlation length deduced from measurements of the QSO and LBG auto-correlation function, and assuming a linear bias model. We also measure a strong auto-correlation of LBGs in our QSO fields finding r0_GG=21.59^{+1.72}_{-1.69} Mpc/h for a fixed slope of gamma=1.5, which is ~4 times larger than the LBG auto-correlation length in random fields, providing further evidence that QSOs reside in overdensities of LBGs. Our results qualitatively support a picture where luminous QSOs inhabit exceptionally massive (M_halo>10^12 M_sun) dark matter halos at z~4., Comment: 25 pages, 22 figures, submitted to the ApJ

We present the results of clustering analyses of Lyman break galaxies (LBGs) at $z\sim3$, $4$, and $5$ using the final data release of the Canada--France--Hawaii Telescope Legacy Survey (CFHTLS). Deep- and wide-field images of the CFHTLS Deep Survey enable us to obtain sufficiently accurate two-point angular correlation functions to apply a halo occupation distribution analysis. Mean halo masses, calculated as $\langle M_{h} \rangle=10^{11.7}-10^{12.8}h^{-1}M_{\odot}$, increase with stellar-mass limit of LBGs. The threshold halo mass to have a central galaxy follows the same increasing trend with the low-$z$ results, whereas the threshold halo mass to have a satellite galaxy shows higher values at $z=3-5$ than $z=0.5-1.5$ over the entire stellar mass range. Satellite fractions of dropout galaxies, even at less massive haloes, are found to drop sharply from $z=2$ down to less than $0.04$ at $z=3-5$. These results suggest that satellite galaxies form inefficiently within dark haloes at $z=3-5$ even for less massive satellites with $M_{\star}<10^{10}M_{\odot}$. We compute stellar-to-halo mass ratios (SHMRs) assuming a main sequence of galaxies, which is found to provide consistent SHMRs with those derived from a spectral energy distribution fitting method. The observed SHMRs are in good agreement with the model predictions based on the abundance-matching method within $1\sigma$ confidence intervals. We derive observationally, for the first time, $M_{{\rm h}}^{{\rm pivot}}$, which is the halo mass at a peak in the star-formation efficiency, at $33$. In addition, $M_{{\rm h}}^{{\rm pivot}}$ and its normalization are found to be almost unchanged during $0