Search

A large fraction of the accreting supermassive black hole population is shrouded by copious amounts of gas and dust, particularly in the distant ($z\gtrsim1$) Universe. While much of the obscuration is attributed to a parsec-scale torus, there is a known contribution from the larger-scale host galaxy. Using JWST/NIRCam imaging from the COSMOS-Web survey, we probe the galaxy-wide dust distribution in X-ray selected AGN up to $z\sim2$. Here, we focus on a sample of three AGNs with their host galaxies exhibiting prominent dust lanes, potentially due to their edge-on alignment. These represent 27% (3 out of 11 with early NIRCam data) of the heavily obscured ($N_H>10^{23}$ cm$^{-2}$) AGN population. With limited signs of a central AGN in the optical and near-infrared, the NIRCam images are used to produce reddening maps $E(B-V)$ of the host galaxies. We compare the mean central value of $E(B-V)$ to the X-ray obscuring column density along the line-of-sight to the AGN ($N_H\sim10^{23-23.5}$ cm$^{-2}$). We find that the extinction due to the host galaxy is present ($0.6\lesssim E(B-V) \lesssim 0.9$; $1.9 \lesssim A_V \lesssim 2.8$) and significantly contributes to the X-ray obscuration at a level of $N_H\sim10^{22.5}$ cm$^{-2}$ assuming an SMC gas-to-dust ratio which amounts to $\lesssim$30% of the total obscuring column density. These early results, including three additional cases from CEERS, demonstrate the ability to resolve such dust structures with JWST and separate the different circumnuclear and galaxy-scale obscuring structures., 12 pages, 8 figures, Accepted for publication in ApJL

Understanding the star formation rate (SFR) variability and how it depends on physical properties of galaxies is important for developing and testing the theory of galaxy formation. We investigate how statistical measurements of the extragalactic background light (EBL) can shed light on this topic and complement traditional methods based on observations of individual galaxies. Using semi-empirical models of galaxy evolution and SFR indicators sensitive to different star formation timescales (e.g., H$\alpha$ and UV continuum luminosities), we show that the SFR variability, quantified by the joint probability distribution of the SFR indicators (i.e., the bivariate conditional luminosity function), can be characterized as a function of galaxy mass and redshift through the cross-correlation between deep, near-infrared maps of the EBL and galaxy distributions. As an example, we consider combining upcoming SPHEREx maps of the EBL with galaxy samples from Rubin/LSST. We demonstrate that their cross-correlation over a sky fraction of $f_\mathrm{sky}\sim0.5$ can constrain the joint SFR indicator distribution at high significance up to $z\sim2.5$ for mass-complete samples of galaxies down to $M_{*}\sim10^9\,M_{\odot}$. These constraints not only allow models of different SFR variability to be distinguished, but also provide unique opportunities to investigate physical mechanisms that require large number statistics such as environmental effects. The cross-correlations investigated illustrate the power of combining cosmological surveys to extract information inaccessible from each data set alone, while the large galaxy populations probed capture ensemble-averaged properties beyond the reach of targeted observations towards individual galaxies., Comment: 12 pages, 7 figures, MNRAS accepted

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{JWST}$/NIRCam+MIRI imaging over the first 20 arcmin$^2$ of publicly-available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W$-$F444W ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED fitting with a variety of codes (including BAGPIPES, PROSPECTOR, BEAGLE, and CIGALE), and find a $>95\%$ probability that these indeed lie at $z>7$. Both sources are compact ($R_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn., Comment: 23 pages, 8 figures, submitted to ApJL, comments welcome

A growing number of far-infrared bright sources completely invisible in deep extragalactic optical surveys hint at an elusive population of z>4 dusty, star-forming galaxies. Cycle 1 JWST surveys are now detecting their rest-frame optical light, which provides key insight into their stellar properties and statistical constraints on the population as a whole. This work presents the JWST/NIRCam counterpart from the COSMOS-Web survey to a far-infrared SCUBA-2 and ALMA source, AzTECC71, which was previously undetected at wavelengths shorter than 850 microns. AzTECC71, amongst the reddest galaxies in COSMOS-Web with F277W - F444W~0.9, is undetected in NIRCam/F150W and F115W and fainter in F444W than other sub-millimeter galaxies identified in COSMOS-Web by 2-4 magnitudes. This is consistent with the system having both a lower stellar mass and higher redshift than the median dusty, star-forming galaxy. With deep ground- and space-based upper limits combined with detections in F277W, F444W and the far-IR including ALMA Band 6, we find a high probability (99%) that AzTECC71 is at z>4 with z_phot=5.7(+0.8,-0.7). This galaxy is massive (logM*/Msun~10.7) and IR-luminous (logLIR/Lsun~12.7), comparable to other optically-undetected but far-IR bright dusty, star-forming galaxies at z>4. This population of luminous, infrared galaxies at z>4 is largely unconstrained but comprises an important bridge between the most extreme dust-obscured galaxies and more typical high-redshift star-forming galaxies. If further far-IR-selected galaxies that drop out of the F150W filter in COSMOS-Web have redshifts z>4 like AzTECC71, then the volume density of such sources may be ~3-10x greater than previously estimated., 13 pages, 7 figures, submitted to ApJ

Utilizing the James Webb Space Telescope Early Release NIRCam Observations, we perform a weak-lensing analysis of the massive galaxy cluster SMACS J0723.3-7327 ($z=0.39$). We investigate the spatial variation of the PSF from the stars in the mosaic image. Our measurements show that the PSF for both modules has very small spatial and temporal variation with average complex ellipticity components of $e_1=0.007\pm0.001$ and $e_2=0.029\pm0.001$ in the observed north-up reference frame. We create PSF models through a principal component analysis of the stars and show that they properly account for the ellipticity of the PSF with residual shapes of $e_1=(0.3\pm3.5)\times10^{-4}$ and $e_2=(1.8\pm4.0)\times10^{-4}$. We select background galaxies by their photometric redshift and measure galaxy shapes by model fitting. Our weak-lensing source catalog achieves 215 galaxies arcmin$^{-2}$. We map the projected mass density of SMACSJ0723 and detect the cluster with a peak significance of $12.2\sigma$. The mass distribution is found to elongate in the east-west direction with an extension to the northeast edge of the field of view where a candidate substructure is found in the Chandra X-ray imaging. We fit the tangential shear with a Navarro-Frenk-White model and estimate the mass of the cluster to be $M_{500}=7.9\pm1.1\times10^{14}$ M$_{\odot}$ ($M_{200}=11.4\pm1.5\times10^{14}$ M$_\odot$ ), which agrees with existing mass estimates. Combining the multiwavelength evidence from literature with our weak-lensing analysis, we hypothesize that SMACSJ0723 is observed near first pericenter passage and we identify candidate radio relics., Comment: 14 pages, 11 figures, submitted to ApJ

We present the results of a systematic search for candidate quiescent galaxies in the distant Universe in eleven $JWST$ fields with publicly available observations collected during the first three months of operations and covering an effective sky area of $\sim145$ arcmin$^2$. We homogeneously reduce the new $JWST$ data and combine them with existing observations from the $Hubble\,Space\,Telescope$. We select a robust sample of $\sim80$ candidate quiescent and quenching galaxies at $3 < z < 5$ using two methods: (1) based on their rest-frame $UVJ$ colors, and (2) a novel quantitative approach based on Gaussian Mixture Modeling of the $NUV-U$, $U-V$, and $V-J$ rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive ($M_\star\geq 10^{10.6} M_\odot$) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of $2-3$, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at $z>3$, as it could be expected for highly clustered massive systems. Importantly, $JWST$ enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available., 15 pages, 6 Figures + Appendix. Accepted for publication in ApJ on Feb, 9. Data release: - Reduced HST+JWST mosaics + photometric catalogs and Eazy-py modeling: https://erda.ku.dk/archives/7166d013c1ca1371aac3c57b9e73190d/published-archive.html - Supplementary material and tables: https://zenodo.org/record/7614908#.Y-4ZruzMLmE - MAST: https://doi.org/10.17909/g3nt-a370 - See also Gould et al. 2023

We present the results of a systematic study of the rest-frame optical morphology of quiescent galaxies at $z \geq 3$ using the Near-Infrared Camera (NIRCam) onboard $JWST$. Based on a sample selected by $UVJ$ color or $NUVUVJ$ color, we focus on 26 quiescent galaxies with $9.810.3$ implies that our size - stellar mass relations are below those at lower redshifts, with the amplitude of $\sim0.6\, {\rm kpc}$ at $M_\star = 5\times 10^{10}\, M_\odot$. This value agrees with the extrapolation from the size evolution of quiescent galaxies at $z, Comment: 23 pages, 16 figures, 2 tables; submitted to ApJ

We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$\sigma$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W) reaching 5$\sigma$ point source depths of $\sim$25.3-26.0 magnitudes. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization ($64$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.5$ and measure its variance with galaxies' star formation histories and morphologies. In addition, we anticipate COSMOS-Web's legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool sub-dwarf stars in the Galactic halo, and possibly the identification of $z>10$ pair-instability supernovae. In this paper we provide an overview of the survey's key measurements, specifications, goals, and prospects for new discovery., Comment: 46 pages, 16 figures, ApJ accepted

We present the radio properties of 66 spectroscopically-confirmed normal star-forming galaxies (SFGs) at $4.4, Comment: Accepted for publication in Astrophysical Journal, 22 pages, 5 figure

ALMA observations of the long wavelength dust continuum are used to estimate the gas masses in a sample of 708 star-forming (SF) galaxies at z = 0.3 to 4.5. We determine the dependence of gas masses and star formation efficiencies (SFE=SFR per unit gass mass). We find that 70 percent of the increase in SFRs of the MS is due to the increased gas masses at earlier epochs while 30 percent is due to increased efficiency of SF. For galaxies above the MS this is reversed with 70 percent of the increased SFR relative to the MS being due to elevated SFEs. Thus, the major evolution of star formation activity at early epochs is driven by increased gas masses, while the starburst activity taking galaxies above the MS is due to enhanced triggering of star formation (likely due to galactic merging). The interstellar gas peaks at z = 2 and dominates the stellar mass down to z = 1.2. Accretion rates needed to maintain continuity of the MS evolution exceed 100 Msun per yr at z > 2. The galactic gas contents are likely the driving determinant for both the rise in SF and AGN activity from z = 5 to their peak at z = 2 and subsequent fall to lower z. We suggest that for self-gravitating clouds with supersonic turbulence, cloud collisions and the filamentary structure of the clouds regulate the star formation activity., AAS Russel Lecture June 2022 (N. Scoville), 18 pages. arXiv admin note: text overlap with arXiv:1702.04729

ISCEA (Infrared Satellite for Cosmic Evolution Astrophysics) is a small astrophysics mission whose Science Goal is to discover how galaxies evolved in the cosmic web of dark matter at cosmic noon. Its Science Objective is to determine the history of star formation and its quenching in galaxies as a function of local density and stellar mass when the Universe was 3-5 Gyrs old (1.2 10 Mpc in each of 50 protocluster (cluster and cosmic web) fields at 1.2 < z < 2.1. ISCEA measures the star formation quenching factor in those fields, and galaxy kinematics with a precision < 50 km/s to deduce the 3D spatial distribution in each field. ISCEA will transform our understanding of galaxy evolution at cosmic noon. ISCEA is a small satellite observatory with a 30cm equivalent diameter aperture telescope with a FoV of 0.32 deg^2, and a multi-object spectrograph with a digital micro-mirror device (DMD) as its programmable slit mask. ISCEA will obtain spectra of 1000 galaxies simultaneously at an effective resolving power of R=1000, with 2.8"x2.8" slits, over the NIR wavelength range of 1.1 to 2.0 microns, a regime not accessible from the ground without large gaps in coverage. ISCEA will achieve a pointing accuracy of, Comment: 40 pages, 31 figures

The faint-end slope of the quasar luminosity function at z~6 and its implication on the role of quasars in reionizing the intergalactic medium at early times has been an outstanding problem for some time. The identification of faint high-redshift quasars with luminosities of 25mag, with 30% of sources having their flux contaminated by foreground objects when the seeing resolution is ~0.7". We mitigate these issues by performing a pixel-level joint processing of ground and space-based data from Subaru/HSC and HST/ACS. We create a deconfused catalog over the 1.64 square-degrees of the COSMOS field, after accounting for spatial varying PSFs and astrometric differences between the two datasets. We identify twelve low-luminosity (M_UV ~ -21 mag) z>6 quasar candidates through (i) their red color measured between ACS/F814W and HSC/i-band and (ii) their compactness in the space-based data. Non-detections of our candidates in Hubble DASH data argues against contamination from late-type stars. Our constraints on the faint end of the quasar luminosity function at z~6.4 suggests a negligibly small contribution to reionization compared to the star-forming galaxy population. The confirmation of our candidates and the evolution of number density with redshift could provide better insights into how supermassive galaxies grew in the first billion years of cosmic time., 24 pages (incl. Appendix), 21 figures, 1 table; accepted for publication in ApJ

The redshift range z=4-6 marks a transition phase between primordial and mature galaxy formation in which galaxies considerably increase their stellar mass, metallicity, and dust content. The study of galaxies in this redshift range is therefore important to understand early galaxy formation and the fate of galaxies at later times. Here, we investigate the burstiness of the recent star-formation history (SFH) of 221 $z\sim4.5$ main-sequence galaxies at log(M) > 9.7 by comparing their ultra-violet (UV) continuum, H$��$ luminosity, and H$��$ equivalent-width (EW). The H$��$ properties are derived from the Spitzer [3.6$��$m]-[4.5$��$m] broad-band color, thereby properly taking into account model and photometric uncertainties. We find a significant scatter between H$��$ and UV-derived luminosities and star-formation rates (SFRs). About half of the galaxies show a significant excess in H$��$ compared to expectations from a constant smooth SFH. We also find a tentative anti-correlation between H$��$ EW and stellar mass, ranging from 1000$��$ at log(M) < 10 to below 100$��$ at log(M) > 11. Consulting models suggests that most $z\sim4.5$ galaxies had a burst of star-formation within the last 50 Myrs, increasing their SFRs by a factor of > 5. The most massive galaxies on the other hand might decrease their SFRs, and may be transitioning to a quiescent stage by z=4. We identify differential dust attenuation (f) between stars and nebular regions as the main contributor to the uncertainty. With local galaxies selected by increasing H$��$ EW (reaching values similar to high-z galaxies), we predict that f approaches unity at $z>4$ consistent with the extrapolation of measurements out to z=2., 31 pages (including appendix), 21 figures, accepted for publication in ApJ

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years., Comment: 30 pages, 10 figures, accepted in ApJ

We present an empirical parameterization of the [NII]/H$��$ flux ratio as a function of stellar mass and redshift valid at 0 < z < 2.7 and 8.5 < log(M) < 11.0. This description can easily be applied to (i) simulations for modeling [NII]$\lambda6584$ line emission, (ii) deblend [NII] and H$��$ in current low-resolution grism and narrow-band observations to derive intrinsic H$��$ fluxes, and (iii) to reliably forecast the number counts of H$��$ emission-line galaxies for future surveys, such as those planned for Euclid and the Wide Field Infrared Survey Telescope (WFIRST). Our model combines the evolution of the locus on the Baldwin, Phillips & Terlevich (BPT) diagram measured in spectroscopic data out to z = 2.5 with the strong dependence of [NII]/H$��$ on stellar mass and [OIII]/H$��$ observed in local galaxy samples. We find large variations in the [NII]/H$��$ flux ratio at a fixed redshift due to its dependency on stellar mass; hence, the assumption of a constant [NII] flux contamination fraction can lead to a significant under- or overestimate of H$��$ luminosities. Specifically, measurements of the intrinsic H$��$ luminosity function derived from current low-resolution grism spectroscopy assuming a constant 29% contamination of [NII] can be overestimated by factors of ~8 at log(L) > 43.0 for galaxies at redshifts z = 1.5. This has implications for the prediction of H$��$ emitters for Euclid and WFIRST. We also study the impact of blended H$��$ and [NII] on the accuracy of measured spectroscopic redshifts., 18 pages, 11 figures, 1 table, accepted for publication in ApJ

We trace the specific star formation rate (sSFR) of massive star-forming galaxies ($\gtrsim\!10^{10}\,\mathcal{M}_\odot$) from $z\sim2$ to 7. Our method is substantially different from previous analyses, as it does not rely on direct estimates of star formation rate, but on the differential evolution of the galaxy stellar mass function (SMF). We show the reliability of this approach by means of semi-analytical and hydrodynamical cosmological simulations. We then apply it to real data, using the SMFs derived in the COSMOS and CANDELS fields. We find that the sSFR is proportional to $(1+z)^{1.1\pm0.2}$ at $z>2$, in agreement with other observations but in tension with the steeper evolution predicted by simulations from $z\sim4$ to 2. We investigate the impact of several sources of observational bias, which however cannot account for this discrepancy. Although the SMF of high-redshift galaxies is still affected by significant errors, we show that future large-area surveys will substantially reduce them, making our method an effective tool to probe the massive end of the main sequence of star-forming galaxies., ApJ accepted

Recent studies have found a significant evolution and scatter in the IRX-$��$ relation at z > 4, suggesting different dust properties of these galaxies. The total far-infrared (FIR) luminosity is key for this analysis but poorly constrained in normal (main-sequence) star-forming z > 5 galaxies where often only one single FIR point is available. To better inform estimates of the FIR luminosity, we construct a sample of local galaxies and three low-redshift analogs of z > 5 systems. The trends in this sample suggest that normal high-redshift galaxies have a warmer infrared (IR) SED compared to average z < 4 galaxies that are used as prior in these studies. The blue-shifted peak and mid-IR excess emission could be explained by a combination of a larger fraction of the metal-poor inter-stellar medium (ISM) being optically thin to ultra-violet (UV) light and a stronger UV radiation field due to high star formation densities. Assuming a maximally warm IR SED suggests 0.6 dex increased total FIR luminosities, which removes some tension between dust attenuation models and observations of the IRX-$��$ relation at z > 5. Despite this, some galaxies still fall below the minimum IRX-$��$ relation derived with standard dust cloud models. We propose that radiation pressure in these highly star-forming galaxies causes a spatial offset between dust clouds and young star-forming regions within the lifetime of O/B stars. These offsets change the radiation balance and create viewing-angle effects that can change UV colors at fixed IRX. We provide a modified model that can explain the location of these galaxies on the IRX-$��$ diagram., 34 pages including Appendix, 18 figures (4 in Appendix), 5 tables (1 in Appendix). Accepted for publication in ApJ

Knowledge of the number density of H$��$ emitting galaxies is vital for assessing the scientific impact of the Euclid and WFIRST missions. In this work we present predictions from a galaxy formation model, Galacticus, for the cumulative number counts of H$��$-emitting galaxies. We couple Galacticus to three different dust attenuation methods and examine the counts using each method. A $��^2$ minimisation approach is used to compare the model predictions to observed galaxy counts and calibrate the dust parameters. We find that weak dust attenuation is required for the Galacticus counts to be broadly consistent with the observations, though the optimum dust parameters return large values for $��^2$, suggesting that further calibration of Galacticus is necessary. The model predictions are also consistent with observed estimates for the optical depth and the H$��$ luminosity function. Finally we present forecasts for the redshift distributions and number counts for two Euclid-like and one WFIRST-like survey. For a Euclid-like survey with redshift range $0.9\leqslant z\leqslant 1.8$ and H$��+{\rm [NII]}$ blended flux limit of $2\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 3900--4800 galaxies per square degree. For a WFIRST-like survey with redshift range $1\leqslant z\leqslant 2$ and blended flux limit of $1\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 10400--15200 galaxies per square degree., 21 pages (including appendix), 12 figures, 6 tables. Accepted by

We study the effects of local environment and stellar mass on galaxy properties using a mass complete sample of quiescent and star-forming systems in the COSMOS field at $z\lesssim$ 3. We show that at $z\lesssim$ 1, the median star-formation rate (SFR) and specific SFR (sSFR) of all galaxies depend on environment, but they become independent of environment at $z\gtrsim$ 1. However, we find that only for \textit{star-forming} galaxies, the median SFR and sSFR are similar in different environments, regardless of redshift and stellar mass. We find that the quiescent fraction depends on environment at $z\lesssim$ 1, and on stellar mass out to $z\sim$ 3. We show that at $z\lesssim$ 1, galaxies become quiescent faster in denser environments and that the overall environmental quenching efficiency increases with cosmic time. Environmental and mass quenching processes depend on each other. At $z\lesssim$ 1, denser environments more efficiently quench galaxies with higher masses (log($M/M_{\odot}$)$\gtrsim$ 10.7), possibly due to a higher merger rate of massive galaxies in denser environments, and that mass quenching is more efficient in denser regions. We show that the overall mass quenching efficiency ($\epsilon_{mass}$) for more massive galaxies (log($M/M_{\odot}$)$\gtrsim$ 10.2) rises with cosmic time until $z\sim$ 1 and flattens out since then. However, for less massive galaxies, the rise in $\epsilon_{mass}$ continues to the present time. Our results suggest that environmental quenching is only relevant at $z\lesssim$ 1, likely a fast process, whereas mass quenching is the dominant mechanism at $z\gtrsim$ 1, with a possible stellar feedback physics., Comment: 17 pages, 10 figures, Accepted for publication in the ApJ

The intrinsic escape fraction of ionizing Lyman continuum photons ($f_{esc}$) is crucial to understand whether galaxies are capable of reionizing the neutral hydrogen in the early universe at z>6. Unfortunately, it is not possible to access $f_{esc}$ at z>4 with direct observations and the handful of measurements from low redshift galaxies consistently find $f_{esc}$ < 10%, while at least $f_{esc}$ ~ 10% is necessary for galaxies dominate reionization. Here, we present the first empirical prediction of $f_{esc}$ at z>6 by combining the (sparsely populated) relation between [OIII]/[OII] and $f_{esc}$ with the redshift evolution of [OIII]/[OII] as predicted from local high-z analogs selected by their H$��$ equivalent-width. We find $f_{esc}$ = $5.7_{-3.3}^{+8.3}$% at z=6 and $f_{esc}$ = $10.4_{-6.3}^{+15.5}$% at z=9 for galaxies with log(M/M$_{sun}$) ~ 9.0 (errors given as 1$��$). However, there is a negative correlation with stellar mass and we find up to 50% larger $f_{esc}$ per 0.5 dex decrease in stellar mass. The population averaged escape fraction increases according to $f_{esc}$ = $f_{esc,0} ((1+z)/3)^a$, with $f_{esc,0} = 2.3 \pm 0.05$% and $a=1.17 \pm 0.02$ at z > 2 for log(M/M$_{sun}$) ~ 9.0. With our empirical prediction of $f_{esc}$ (thus fixing an important previously unknown variable) and further reasonable assumption on clumping factor and the production efficiency of Lyman continuum photons, we conclude that the average population of galaxies is just capable to reionize the universe by z ~ 6., 10 pages, 4 figure, 1 table. Accepted by ApJ

We use the Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) strong lensing image catalog of the merging galaxy cluster Abell 370 to obtain a mass model using the free-form lens inversion algorithm GRALE. The improvement of the strong lensing data quality results in a lens plane rms of only 0.45 arcsec, about a factor of two lower than that of our existing HFF v4 reconstruction. We attribute the improvement to spectroscopic data and use of the full reprocessed HST mosaics. In our reconstructed mass model, we found indications of three distinct mass features in Abell 370: (i) a $\sim 35$ kpc offset between the northern BCG and the nearest mass peak, (ii) a $\sim 100$ kpc mass concentration of roughly critical density $\sim 250$ kpc east of the main cluster, and (iii) a probable filament-like structure passing N-S through the cluster. While (i) is present in some form in most publicly available reconstructions spanning the range of modeling techniques: parametric, hybrid, and free-form, (ii) and (iii) are recovered by only about half of the reconstructions. We tested our hypothesis on the presence of the filament-like structure by creating a synthetic cluster - Irtysh IIIc - mocking the situation of a cluster with external mass. We also computed the source plane magnification distributions. Using them we estimated the probabilities of magnifications in the source plane, and scrutinized their redshift dependence. Finally, we explored the lensing effects of Abell 370 on the luminosity functions of sources at $z_s=9.0$, finding it consistent with published results., Comment: 16 pages, 13 figures