Search

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides deep HST F275W and F435W imaging over four CANDELS fields (GOODS-N, GOODS-S, COSMOS, and EGS). We combine this newly acquired UV imaging with existing HST imaging from CANDELS as well as existing ancillary data to obtain robust photometric redshifts and reliable estimates for galaxy physical properties for over 150,000 galaxies in the $\sim$430 arcmin$^2$ UVCANDELS area. Here, we leverage the power of the new UV photometry to not only improve the photometric redshift measurements in these fields, but also constrain the full redshift probability distribution combining multiple redshift fitting tools. Furthermore, using the full UV-to-IR photometric dataset, we measure the galaxy physical properties by fitting templates from population synthesis models with two different parameterizations (flexible and fixed-form) of the star-formation histories (SFHs). Compared to the flexible SFH parametrization, we find that the fixed-form SFHs systematically underestimate the galaxy stellar masses, both at the low- ($\lesssim10^9 M_\odot$) and high- ($\gtrsim10^{10} M_\odot$) mass end, by as much as $\sim0.5$ dex. This underestimation is primarily due the limited ability of fixed-form SFH parameterization to simultaneously capture the chaotic nature of star-formation in these galaxies., Comment: 22 pages, 6 figures; accepted to ApJS; catalogs available via MAST

Galactic outflows provide important feedback effects to regulate the evolution of the host galaxies. Two primary diagnostics of galactic outflows are broad and/or blueshifted emission and absorption lines. Even though well-established methods exist to analyze these outflow signatures, connections between them are rarely studied and largely unknown. In this paper, we present the first detailed comparisons of the outflow properties measured independently from the two outflow diagnostics for a sample of 33 low-redshift star-forming galaxies. Their UV absorption lines are detected by the Hubble Space Telescope/Cosmic Origin Spectrograph, and optical emission lines are observed by the Keck/Echellette Spectrograph and Imager. We find that several outflow properties derived from emission and absorption lines are tightly correlated. These include outflow maximum velocity, line width, and sizes. Specifically, in a given galaxy, outflows seen in emission lines have smaller maximum velocities, narrower line widths, and smaller sizes than those measured from the absorption lines. These findings can be interpreted by the fact that emission line luminosity is weighted by density squared, while absorption line depth is weighted by density. We then test both spherical and bi-conical outflow models, and find the same outflow velocity and density distributions can explain the observed outflow features in emission and absorption lines for individual galaxies. These results provide novel calibration between galactic outflow properties measured from the two diagnostics and provide valuable insights for future models of galactic outflows by potentially doubling the number of observational constraints., Comment: 3 tables, 17 figures, 30 pages, submitted to ApJ

One of the fundamental questions of cosmology is the origin and mechanism(s) responsible for the reionization of the Universe beyond $z\sim6$. To address this question, many studies over the past decade have focused on local ($z\sim0.3$) galaxies which leak ionizing radiation (Lyman continuum or LyC). However, line-of-sight effects and data quality have prohibited deeper insight into the nature of LyC escape. To circumvent these limitations, we analyze stacks of a consolidated sample of {\it HST}/COS observations of the LyC in 89 galaxies at $z\sim0.3$. From fitting of the continuum, we obtain information about the underlying stellar populations and neutral ISM geometry. We find that most LyC non-detections are not leaking appreciable LyC ($f_{esc}^{\rm LyC}<1$\%) but also that exceptional cases point to spatial variations in the LyC escape fraction $f_{esc}^{\rm LyC}$. Stellar populations younger than 3 Myr lead to an increase in ionizing feedback, which in turn increases the isotropy of LyC escape. Moreover, mechanical feedback from supernovae in 8-10 Myr stellar populations is important for anisotropic gas distributions needed for LyC escape. While mechanical feedback is necessary for any LyC escape, high $f_{esc}^{\rm LyC}$ ($>5$\%) also requires a confluence of young stars and ionizing feedback. A two-stage burst of star formation could facilitate this optimal LyC escape scenario., Comment: Submitted to AAS publications, 42 pages, 25 figures

Feedback is widely recognized as an essential condition for Lyman continuum (LyC) escape in star-forming galaxies. However, the mechanisms by which galactic outflows clear neutral gas and dust remain unclear. In this paper, we model the Mg II 2796\r{A}, 2804\r{A} absorption + emission lines in 29 galaxies taken from the Low-z LyC Survey (LzLCS) to investigate the impact of (radiation + mechanical) feedback on LyC escape. Using constraints on Mg$^+$ and photoionization models, we map the outflows' neutral hydrogen content and predict $f_{esc}^{LyC}$ with a multiphase wind model. We measure mass, momentum, and energy loading factors for the neutral winds, which carry up to 10% of the momentum and 1% of the energy in SFR-based deposition rates. We use SED template fitting to determine the relative ages of stellar populations, allowing us to identify radiation feedback dominant systems. We then examine feedback related properties (stellar age, loading factors, etc.) under conditions that optimize feedback efficiency, specifically high star formation rate surface density and compact UV half-light radii. Our findings indicate that the strongest leakers are radiation feedback dominant, lack Mg II outflows, but have extended broad components in higher ionization lines like [O III] 5007\r{A}, as observed by Amor\'in et al. (2024). In contrast, galaxies experiencing supernovae feedback typically exhibit weaker $f_{esc}^{LyC}$ and show evidence of outflows in both Mg II and higher ionization lines. We attribute these findings to rapid or "catastrophic" cooling in the radiation-dominant systems, which, given the low metallicities in our sample, are likely experiencing delayed supernovae., Comment: 34 pages, 16 figures, 7 tables

We present morphological classifications based on Galaxy Zoo analysis of 71,052 galaxies with imaging from the United Kingdom Infrared Telescope Infrared Deep Sky Survey (UKIDSS). Galaxies were selected out of the Galaxy Zoo 2 (GZ2) sample, so also have gri imaging from the Sloan Digital Sky Survey. An identical classification tree, and vote weighting/aggregation was applied to both UKIDSS and GZ2 classifications enabling direct comparisons. With this Research Note we provide a public release of the GZ:UKIDSS morphologies and discuss some initial comparisons with GZ2., Comment: 3 pages, 1 figure

JWST is uncovering the properties of ever increasing numbers of galaxies at z>6, during the epoch of reionization. Connecting these observed populations to the process of reionization requires understanding how efficiently they produce Lyman continuum (LyC) photons and what fraction (fesc) of these photons escape into the intergalactic medium. By applying the Cox proportional hazards model, a survival analysis technique, to the Low-redshift Lyman Continuum Survey (LzLCS), we develop new, empirical, multivariate predictions for fesc. The models developed from the LzLCS reproduce the observed fesc for z~3 samples, which suggests that LyC emitters may share similar properties at low and high redshift. Our best-performing models for the z~3 galaxies include information about dust attenuation, ionization, and/or morphology. We then apply these models to z$\gtrsim$6 galaxies. For large photometric samples, we find a median predicted fesc=0.047-0.14. For smaller spectroscopic samples, which may include stronger emission line galaxies, we find that $\geq$33% of the galaxies have fesc >0.2, and we identify several candidate extreme leakers with fesc $\geq$0.5. The current samples show no strong trend between predicted fesc and UV magnitude, but limited spectroscopic information makes this result uncertain. Multivariate predictions can give significantly different results from single variable predictions, and the predicted fesc for high-redshift galaxies can differ significantly depending on whether star formation rate surface density or radius is used as a measure of galaxy morphology. We provide all parameters necessary to predict fesc for additional samples of high-redshift galaxies using these models., Comment: Accepted for publication in ApJ. 33 pages, 9 figures, 10 tables, plus appendix

To understand how galaxies reionized the universe, we must determine how the escape fraction of Lyman Continuum (LyC) photons (fesc) depends on galaxy properties. Using the z~0.3 Low-redshift Lyman Continuum Survey (LzLCS), we develop and analyze new multivariate predictors of fesc. These predictions use the Cox proportional hazards model, a survival analysis technique that incorporates both detections and upper limits. Our best model predicts the LzLCS fesc detections with a root-mean-square (RMS) scatter of 0.31 dex, better than single-variable correlations. According to ranking techniques, the most important predictors of fesc are the equivalent width (EW) of Lyman-series absorption lines and the UV dust attenuation, which track line-of-sight absorption due to HI and dust. The HI absorption EW is uniquely crucial for predicting fesc for the strongest LyC emitters, which show properties similar to weaker LyC emitters and whose high fesc may therefore result from favorable orientation. In the absence of HI information, star formation rate surface density ($\Sigma_{\rm SFR}$) and [O III]/[O II] ratio are the most predictive variables and highlight the connection between feedback and fesc. We generate a model suitable for z>6, which uses only the UV slope, $\Sigma_{\rm SFR}$, and [O III]/[O II]. We find that $\Sigma_{\rm SFR}$ is more important in predicting fesc at higher stellar masses, whereas [O III]/[O II] plays a greater role at lower masses. We also analyze predictions for other parameters, such as the ionizing-to-non ionizing flux ratio and Ly=alpha escape fraction. These multivariate models represent a promising tool for predicting fesc at high redshift., Comment: Accepted for publication in ApJ. 34 pages + appendix, 12 figures

The escape fraction of LyC ionizing radiation $f_{LyC}^{esc}$ is crucial for understanding reionization, yet impossible to measure at z $\gtrsim$ 5.3. Recently, studies have focused on calibrating indirect indicators of $f_{LyC}^{esc}$ at z $\sim$ 0.3, finding that Ly$\alpha$ is closely linked to it. What is still unclear is whether the LyC - Ly$\alpha$ relation evolves with redshift, and if Ly$\alpha$ is truly applicable as an $f_{LyC}^{esc}$ indicator during the reionization epoch. In this study, we investigate seven $-21 \lesssim M_{UV} \lesssim -19$ gravitationally lensed galaxies from the BELLS GALLERY survey at z $\sim$ 2.3. Our targets have rest-frame Ly$\alpha$ equivalent widths ranging from 10 \r{A} to 100 \r{A} and low dust content ($-2.5 \lesssim \beta \lesssim -1.9$), both indicative of high LyC escape. Surprisingly, direct estimates of $f_{LyC}^{esc}$ using Hubble Space Telescope imaging with F275W and F225W reveal that our targets are not LyC emitters, with absolute $f_{LyC}^{esc}$ $\lesssim$ 6.5% at 3$\sigma$ significance (with two sources having absolute $f_{LyC}^{esc}$(3$\sigma$) $\lesssim$ 10% and $\lesssim$ 16%). The low $f_{LyC}^{esc}$, coupled with the high Ly$\alpha$ escape fraction and equivalent width could potentially be attributed to the redshift evolution of the neutral hydrogen column density and dust content. Our results challenge previous studies based on local samples, suggesting that the extrapolation of z ~ 0 Ly$\alpha$-based LyC indirect estimators into the reionization epoch might not be correct., Comment: 19 Pages, 8 Figures. Submitted to ApJ

Low-mass galaxies can significantly contribute to reionization due to their potentially high Lyman continuum (LyC) escape fraction and relatively high space density. We present a constraint on the LyC escape fraction from low-mass galaxies at z = 1.3 - 3.0. We obtained rest-frame UV continuum imaging with the ACS/SBC and the WFC3/UVIS from the Hubble Space Telescope for eight strongly-lensed galaxies that were identified in the Sloan Giant Arc Survey (SGAS) and the Cluster Lensing And Supernova survey with Hubble (CLASH). The targeted galaxies were selected to be spectroscopically confirmed, highly magnified, and blue in their UV spectral shapes ($\beta<-1.7$). Our targets include intrinsically low luminosity galaxies down to a magnification-corrected absolute UV magnitude of $M_{\rm UV}\sim-14$. We perform custom-defined aperture photometry to place the most reliable upper limits of LyC escape from our sample. From our observations, we report no significant ($>$$2\sigma$) detections of LyC fluxes, placing 1$\sigma$ upper limits on the absolute LyC escape fractions of 3 - 15%. Our observations do not support the expected increased escape fractions of LyC photons from intrinsically UV faint sources. Considering the highly anisotropic geometry of LyC escape, increasing the sample size of faint galaxies in future LyC observations is crucial., Comment: Submitted to ApJ (16 pages, 6 figures, 5 tables)

The geometry of the neutral gas in and around galaxies is a key regulator of the escape of ionizing photons. We present the first statistical study aiming at linking the neutral and ionized gas distributions to the Lyman continuum (LyC) escape fraction (fesc(LyC)) in a sample of 22 confirmed LyC leakers and non-leakers at z~0.35 using the Keck Cosmic Web Imager (Keck/KCWI) and the Low Resolution Spectrograph 2 (HET/LRS2). Our integral field unit data enable the detection of neutral and low-ionization gas, as traced by Mg II, and ionized gas, as traced by [O II], extending beyond the stellar continuum for 7 and 10 objects, respectively. All but one object with extended Mg II emission also shows extended [O II] emission; in this case, Mg II emission is always more extended than [O II] by a factor 1.3 on average. Most of the galaxies with extended emission are non or weak LyC leakers (fesc(LyC) < 5%), but we find a large diversity of neutral gas configurations around these weakly LyC-emitting galaxies. Conversely, the strongest leakers (fesc(LyC) > 10%) appear uniformly compact in both Mg II and [O II] with exponential scale lengths <1 kpc. We also find a trend between fesc(LyC) and the spatial offsets of the nebular gas and the stellar continuum emission. Moreover, we find significant anti-correlations between the spatial extent of the neutral gas and the [O III]/[O II] ratio, and H$\beta$ equivalent width, as well as positive correlations with metallicity and UV size, suggesting that galaxies with more compact neutral gas sizes are more highly ionized. The observations suggest that strong LyC emitters do not have extended neutral gas halos and ionizing photons may be emitted in many directions. Combined with high ionization diagnostics, we propose the Mg II, and potentially [O II], spatial compactness are indirect indicators of LyC emitting galaxies at high-redshift., Comment: 31 pages, 24 figures, 4 tables, submitted to A&A

UVCANDELS is a HST Cycle-26 Treasury Program awarded 164 orbits of primary ultraviolet (UV) F275W imaging and coordinated parallel optical F435W imaging in four CANDELS fields: GOODS-N, GOODS-S, EGS, and COSMOS, covering a total area of $\sim426$ arcmin$^2$. This is $\sim2.7$ times larger than the area covered by previous deep-field space UV data combined, reaching a depth of about 27 and 28 ABmag ($5\sigma$ in $0.2"$ apertures) for F275W and F435W, respectively. Along with the new photometric catalogs, we present an analysis of the rest-frame UV luminosity function (LF), relying on our UV-optimized aperture photometry method yielding a factor of $1.5\times$ increase than the H-isophot aperture photometry in the signal-to-noise ratios of galaxies in our F275W imaging. Using well tested photometric redshift measurements we identify 5810 galaxies at redshifts $0.6

Studying the galaxies responsible for reionization is often conducted through local reionization-era analogs; however, many of these local analogs are too massive to be representative of the low-mass star-forming galaxies that are thought to play a dominant role in reionization. The local, low-mass dwarf starburst galaxy Pox 186 is one such system with physical conditions representative of a reionization-era starburst galaxy. We present deep ultraviolet (UV) spectroscopy of Pox 186 to study its stellar population and ionization conditions and to compare these conditions to other local starburst galaxies. The new Cosmic Origins Spectrograph data are combined with archival observations to cover $\sim$1150-2000 A and allow for an assessment of Pox 186's stellar population, the relative enrichment of C and O, and the escape of ionizing photons. We detect significant Ly$\alpha$ and low-ionization state absorption features, indicative of previously undetected neutral gas in Pox 186. The C/O relative abundance, log(C/O) = -0.62$\pm$0.02, is consistent with other low-metallicity dwarf galaxies and suggests a comparable star formation history in these systems. We compare UV line ratios in Pox 186 to those of dwarf galaxies and photoionization models, and we find excellent agreement for the ratios utilizing the intense C III], O III], and double-peaked C IV lines. However, the UV and optical He II emission is faint and distinguishes Pox 186 from other local starburst dwarf galaxies. We explore mechanisms that could produce faint He II, which have implications for the low-mass reionization-era galaxies which may have similar ionization conditions., Comment: 22 pages, 9 figures, accepted for publication in The Astrophysical Journal

The UltraViolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) survey is a Hubble Space Telescope (HST) Cycle-26 Treasury Program, allocated in total 164 orbits of primary Wide-Field Camera 3 Ultraviolet and Visible light F275W imaging with coordinated parallel Advanced Camera for Surveys F435W imaging, on four of the five premier extragalactic survey fields: GOODS-N, GOODS-S, EGS, and COSMOS. We introduce this survey by presenting a thorough search for galaxies at $z\gtrsim2.4$ that leak significant Lyman continuum (LyC) radiation, as well as a stringent constraint on the LyC escape fraction ($f_{\rm esc}$) from stacking the UV images of a population of star-forming galaxies with secure redshifts. Our extensive search for LyC emission and stacking analysis benefit from the catalogs of high-quality spectroscopic redshifts compiled from archival ground-based data and HST slitless spectroscopy, carefully vetted by dedicated visual inspection efforts. We report a sample of five galaxies as individual LyC leaker candidates, showing $f_{\rm esc}^{\rm rel}\gtrsim60\%$ estimated using detailed Monte Carlo analysis of intergalactic medium attenuation. We develop a robust stacking method to apply to five samples of in total 85 non-detection galaxies in the redshift range of $z\in[2.4,3.7]$. Most stacks give tight 2-$\sigma$ upper limits below $f_{\rm esc}^{\rm rel}<6\%$. A stack for a subset of 32 emission-line galaxies shows tentative LyC leakage detected at 2.9-$\sigma$, indicating $f_{\rm esc}^{\rm rel}=5.7\%$ at $z\sim2.65$, supporting the key role of such galaxies in contributing to the cosmic reionization and maintaining the UV ionization background. These new F275W and F435W imaging mosaics from UVCANDELS have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes., Comment: 33 pages, 21 figures, and 5 tables. Resubmitted after addressing the referee report

Lyman-alpha line profiles are a powerful probe of ISM structure, outflow speed, and Lyman continuum escape fraction. In this paper, we present the Ly$\alpha$ line profiles of the COS Legacy Archive Spectroscopic SurveY, a sample rich in spectroscopic analogs of reionization-era galaxies. A large fraction of the spectra show a complex profile, consisting of a double-peaked Ly$\alpha$ emission profile in the bottom of a damped, Ly$\alpha$ absorption trough. Such profiles reveal an inhomogeneous interstellar medium (ISM). We successfully fit the damped Ly$\alpha$ absorption (DLA) and the Ly$\alpha$ emission profiles separately, but with complementary covering factors, a surprising result because this approach requires no Ly$\alpha$ exchange between high-$N_\mathrm{HI}$ and low-$N_\mathrm{HI}$ paths. The combined distribution of column densities is qualitatively similar to the bimodal distributions observed in numerical simulations. We find an inverse relation between Ly$\alpha$ peak separation and the [O III]/[O II] flux ratio, confirming that the covering fraction of Lyman-continuum-thin sightlines increases as the Ly$\alpha$ peak separation decreases. We combine measurements of Ly$\alpha$ peak separation and Ly$\alpha$ red peak asymmetry in a diagnostic diagram which identifies six Lyman continuum leakers in the CLASSY sample. We find a strong correlation between the Ly$\alpha$ trough velocity and the outflow velocity measured from interstellar absorption lines. We argue that greater vignetting of the blueshifted Ly$\alpha$ peak, relative to the redshifted peak, is the source of the well-known discrepancy between shell-model parameters and directly measured outflow properties. The CLASSY sample illustrates how scattering of Ly$\alpha$ photons outside the spectroscopic aperture reshapes Ly$\alpha$ profiles as the distances to these compact starbursts span a large range., Comment: 40 pages, 19 figures, 5 tables, accepted for publication in ApJ, comments welcome

In the current JWST era, rest-frame UV spectra play a crucial role in enhancing our understanding of the interstellar medium (ISM) and stellar properties of the first galaxies in the epoch of reionization (EoR, $z>6$). Here, we compare well-known and reliable optical diagrams sensitive to the main ionization source (i.e., star formation, SF; active galactic nuclei, AGN; shocks) to UV counterparts proposed in the literature - the so-called ``UV-BPT diagrams'' - using the HST COS Legacy Archive Spectroscopic SurveY (CLASSY), the largest high-quality, high-resolution and broad-wavelength range atlas of far-UV spectra for 45 local star-forming galaxies. In particular, we explore where CLASSY UV line ratios are located in the different UV diagnostic plots, taking into account state-of-the-art photoionization and shock models and, for the first time, the measured ISM and stellar properties (e.g., gas-phase metallicity, ionization parameter, carbon abundance, stellar age). We find that the combination of C III] $\lambda\lambda$1907,9 He II $\lambda1640$ and O III] $\lambda$1666 can be a powerful tool to separate between SF, shocks and AGN at sub-solar metallicities. We also confirm that alternative diagrams without O III] $\lambda$1666 still allow us to define a SF-locus with some caveats. Diagrams including C IV $\lambda\lambda$1548,51 should be taken with caution given the complexity of this doublet profile. Finally, we present a discussion detailing the ISM conditions required to detect UV emission lines, visible only in low gas-phase metallicity (12+log(O/H) $\lesssim8.3$) and high ionization parameter (log($U$) $\gtrsim-2.5$) environments. Overall, CLASSY and our UV toolkit will be crucial in interpreting the spectra of the earliest galaxies that JWST is currently revealing., Comment: 36 pages, 15 figures, accepted for publication in ApJ

Observations of low-ionization state (LIS) metal lines provide crucial insights into the interstellar medium of galaxies, yet, disentangling the physical processes responsible for the emerging line profiles is difficult. This work investigates how mock spectra generated using a single galaxy in a radiation-hydrodynamical simulation can help us interpret observations of a real galaxy. We create 22,500 C II and Si II spectra from the virtual galaxy at different times and through multiple lines of sight and compare them with the 45 observations of low-redshift star-forming galaxies from the COS Legacy Spectroscopic SurveY (CLASSY). We find that the mock profiles provide accurate replicates to the observations of 38 galaxies with a broad range of stellar masses ($10^6$ to $10^9$ $M_\odot$) and metallicities (0.02 to 0.55 $Z_\odot$). Additionally, we highlight that aperture losses explain the weakness of the fluorescent emission in several CLASSY spectra and must be accounted for when comparing simulations to observations. Overall, we show that the evolution of a single simulated galaxy can produce a large diversity of spectra whose properties are representative of galaxies of comparable or smaller masses. Building upon these results, we explore the origin of the continuum, residual flux, and fluorescent emission in the simulation. We find that these different spectral features all emerge from distinct regions in the galaxy's ISM, and their characteristics can vary as a function of the viewing angle. While these outcomes challenge simplified interpretations of down-the-barrel spectra, our results indicate that high-resolution simulations provide an optimal framework to interpret these observations., Comment: Accepted for publication in Apj

We present the results of deep, ground based U-band imaging with the Large Binocular Telescope of the Cosmic Evolution Survey (COSMOS) field as part of the near-UV imaging program, UVCANDELS. We utilize a seeing sorted stacking method along with night-to-night relative transparency corrections to create optimal depth and optimal resolution mosaics in the U-band, which are capable of reaching point source magnitudes of AB 26.5 mag at 3 sigma. These ground based mosaics bridge the wavelength gap between the HST WFC3 F27W and ACS F435W images and are necessary to understand galaxy assembly in the last 9-10 Gyr. We use the depth of these mosaics to search for the presence of U-band intragroup light (IGrL) beyond the local Universe. Regardless of how groups are scaled and stacked, we do not detect any U-band IGrL to unprecedented U-band depths of 29.1-29.6 mag/arcsec2, which corresponds to an IGrL fraction of less than 1% of the total group light. This stringent upper limit suggests that IGrL does not contribute significantly to the Extragalactic Background Light at short wavelengths. Furthermore, the lack of UV IGrL observed in these stacks suggests that the atomic gas observed in the intragroup medium (IGrM) is likely not dense enough to trigger star formation on large scales. Future studies may detect IGrL by creating similar stacks at longer wavelengths or by pre-selecting groups which are older and/or more dynamically evolved similar to past IGrL observations of compact groups and loose groups with signs of gravitational interactions., Comment: Accepted to PASP

High-resolution imaging of galaxies in rest-frame UV has revealed the existence of giant star-forming clumps prevalent in high redshift galaxies. Studying these sub-structures provides important information about their formation and evolution and informs theoretical galaxy evolution models. We present a new method to identify clumps in galaxies' high-resolution rest-frame UV images. Using imaging data from CANDELS and UVCANDELS, we identify star-forming clumps in an HST/F160W$\leq 25$ AB mag sample of 6767 galaxies at $0.5\leq z\leq 3$ in four fields, GOODS-N, GOODS-S, EGS, and COSMOS. We use a low-pass band filter in Fourier space to reconstruct the background image of a galaxy and detect small-scale features (clumps) on the background-subtracted image. Clumpy galaxies are defined as those having at least one off-center clump that contributes a minimum of 10$\%$ of the galaxy's total rest-frame UV flux. We measure the fraction of clumpy galaxies ($\rm f_{clumpy}$) as a function of stellar mass, redshift, and galaxy environment. Our results indicate that $\rm f_{clumpy}$ increases with redshift, reaching $\sim 65\%$ at $z\sim 1.5$. We also find that $\rm f_{clumpy}$ in low-mass galaxies ($\rm 9.5\leq log(M_*/M_\odot)\leq 10$) is 10$\%$ higher compared to that of their high-mass counterparts ($\rm log(M_*/M_\odot)>10.5$). Moreover, we find no evidence of significant environmental dependence of $\rm f_{clumpy}$ for galaxies at the redshift range of this study. Our results suggest that the fragmentation of gas clouds under violent disk instability remains the primary driving mechanism for clump formation, and incidents common in dense environments, such as mergers, are not the dominant processes., Comment: 16 pages, 11 figures, 2 tables, accepted for publication in ApJ

We investigate the impact of field-to-field variation, deriving from cosmic variance, in measured Lyman-$\alpha$ emitter (LAE) luminosity functions (LFs) and this variation's impact on inferences of the neutral fraction of the intergalactic medium (IGM) during reionization. We post-process a z=7 IGM simulation to populate the dark matter halos with LAEs. These LAEs have realistic UV magnitudes, Ly$\alpha$ fluxes, and Ly$\alpha$ line profiles. We calculate the attenuation of Ly$\alpha$ emission in universes with varying IGM neutral fraction, $\bar{\rm{x}}_{\rm{HI}}$. In a $\bar{\rm{x}}_{\rm{HI}}=0.3$ simulation, we perform 100 realizations of a mock 2 square degree survey with a redshift window $\Delta z = 0.5$ and flux limit $\rm{f}_{Ly\alpha}>1\times10^{-17}\:\rm{ergs}\:\: \rm{s}^{-1} \: \rm{cm}^{-2}$; such a survey is typical in depth and volume of the largest LAE surveys conducted today. For each realization, we compute the LAE LF and use it to recover the input $\bar{\rm{x}}_{\rm{HI}}$. Comparing the inferred values of $\bar{\rm{x}}_{\rm{HI}}$ across the ensemble of the surveys, we find that cosmic variance, deriving from large-scale structure and variation in the neutral gas along the sightline, imposes a floor in the uncertainty of $\Delta \bar{\rm{x}}_{\rm{HI}} \sim 0.2$ when $\bar{\rm{x}}_{\rm{HI}}$ $=0.3$. We explore mitigation strategies to decrease this uncertainty, such as increasing the volume, decreasing the flux limit, or probing the volume with many independent fields. Increasing the area and/or depth of the survey does not mitigate the uncertainty, but composing a survey with many independent fields is effective. This finding highlights the best strategy for LAE surveys aiming at constraining $\bar{\rm{x}}_{\rm{HI}}$ of the universe during reionization., Comment: 17 pages, 13 figures

Lyman $\alpha$ (Ly$\alpha$) is now routinely used as a tool for studying high-redshift galaxies and its resonant nature means it can trace neutral hydrogen around star-forming galaxies. Integral field spectrograph measurements of high-redshift Ly$\alpha$ emitters indicate that significant extended Ly$\alpha$ halo emission is ubiquitous around such objects. We present a sample of redshift 0.23 to 0.31 galaxies observed with the Hubble Space Telescope selected to match the star formation properties of high-$z$ samples while optimizing the observations for detection of low surface brightness Ly$\alpha$ emission. The Ly$\alpha$ escape fractions range between 0.7\% and 37\%, and we detect extended Ly$\alpha$ emission around six out of seven targets. We find Ly$\alpha$ halo to UV scale length ratios around 6:1 which is marginally lower than high-redshift observations, and halo flux fractions between 60\% and 85\% -- consistent with high-redshift observations -- when using comparable methods. However, our targets show additional extended stellar UV emission: we parametrize this with a new double exponential model. We find that this parametrization does not strongly affect the observed Ly$\alpha$ halo fractions. We find that deeper H$\alpha$ data would be required to firmly determine the origin of Ly$\alpha$ halo emission, however, there are indications that H$\alpha$ is more extended than the central FUV profile, potentially indicating conditions favorable for the escape of ionizing radiation. We discuss our results in the context of high-redshift galaxies, cosmological simulations, evolutionary studies of the circumgalactic medium in emission, and the emission of ionizing radiation., Comment: 20 page, 14 figures, 6 tables. Accepted for publication in MNRAS

We study the cross-correlation between the stochastic gravitational-wave background (SGWB) generated by binary black hole (BBH) mergers across the universe and the distribution of galaxies across the sky. We use the anisotropic SGWB measurement obtained using data from the third observing run (O3) of Advanced LIGO detectors and galaxy over-density obtained from the Sloan Digital Sky Survey (SDSS) spectroscopic catalog. We compute, for the first time, the angular power spectrum of their cross-correlation. Instead of integrating the SGWB across frequencies, we analyze the cross-correlation in 10 Hz wide SGWB frequency bands to study the frequency dependence of the cross-correlation angular power spectrum. Finally, we compare the observed cross-correlation to the spectra predicted by astrophysical models. We apply a Bayesian formalism to explore the parameter space of the theoretical models, and we set constraints on a set of (effective) astrophysical parameters describing the galactic process of gravitational wave (GW) emission. Parameterizing with a Gaussian function the astrophysical kernel describing the local process of GW emission at galactic scales, we find the 95\% upper limit on kernel amplitude to be $2.7 \times 10^{-32}$ erg cm$^{-3}$s$^{-1/3}$ when ignoring the shot noise in the GW emission process, and $2.16 \times 10^{-32}$ erg cm$^{-3}$s$^{-1/3}$ when the shot noise is included in the analysis. As the sensitivity of the LIGO-Virgo-KAGRA network improves, we expect to be able to set more stringent bounds on this kernel function and constrain its parameters., Comment: 19 pages, 12 figures

In gravitational lensing, central images in quads can serve as a powerful probe of the inner regions of lens galaxies. The presence of an offset central supermassive black hole (SMBH) has the potential to distort the time-delay surface in a way such that 3 central images form: a strongly de-magnified image near the SMBH, and two less de-magnified (and potentially observable) images at a central maximum and saddle point. Using a quad lens macro model, we simulate the constraints that could be placed on various lens galaxy parameters based on their central images probability of detection or non-detection. Informed by a recent low-redshift distribution of off-nucleus AGN, we utilize Bayesian inference to constrain the mean SMBH off-nucleus distance and galactic core radius for a sample of 6 quads. In general, we find that a detection of the central image in any quad would favor larger SMBH off-nucleus distances and galaxy core sizes. Assuming a linear relationship between core radii and velocity dispersion $r_c = b\sigma$, these results similarly imply strong constraints on $b$, where the likely case of a central image non-detection in each quad constraining $b$ to $3.11^{+2.72}_{-2.26} \times 10^{-4}$ kpc km$^{-1}$ s. Our results show that tight constraints on lens galaxy parameters can be made regardless of a detection or non-detection of a central image. Therefore, we recommend observational searches for the central image, possibly using our suggested novel detection technique in UV filters, to formalize stronger constraints on lens galaxy parameters., Comment: MNRAS, Accepted, 20 pages, 15 figures

We present the Texas Euclid Survey for Lyman-Alpha (TESLA), a spectroscopic survey in the 10 square degree of the Euclid North Ecliptic Pole (NEP) field. Using TESLA, we study how the physical properties of Lyman-alpha emitters (LAEs) correlate with Lyman-alpha emission to understand the escape of Lyman alpha from galaxies at redshifts 2 -- 3.5. We present an analysis of 43 LAEs performed in the NEP field using early data from the TESLA survey. We use Subaru Hyper Suprime-Cam imaging in the grizy-bands, Spitzer/IRAC channels 1 and 2 from the Hawaii 20 square degree (H20) survey and spectra acquired by the Visible Integral-Field Replicable Unit Spectrograph (VIRUS) on the Hobby-Eberly Telescope. We perform spectral energy distribution (SED) fitting to compute the galaxy properties of 43 LAEs, and study correlations between stellar mass, star formation rate (SFR), and dust, to the Lyman-alpha rest-frame equivalent widths (EW). We uncover marginal (1 sigma significance) correlations between stellar mass and Lyman-alpha EW, and star formation rate (SFR) and Lyman-alpha EW, with a Spearman correlation coefficient of -0.$34_{-.14}^{+.17}$ and -0.$37_{-.14}^{+.16}$ respectively. We show that the Lyman-alpha distribution of the 43 LAEs is consistent with being drawn from an exponential distribution with an e-folding scale of 150 Angstrom. Once complete the TESLA survey will enable the study of ~ thousands of LAEs to explore correlations between galaxy properties and Lyman-alpha EW. The large sample size will allow the construction of a predictive model for the Lyman-alpha EW as a function of SED-derived galaxy properties, which could be used to improve Lyman-alpha based constraints on reionization.

Lyman-alpha (Ly$\alpha$) emission from galaxies can be used to trace neutral hydrogen in the epoch of reionization, however, there is a degeneracy between the attenuation of Ly$\alpha$ in the intergalactic medium (IGM) and the line profile emitted from the galaxy. Large shifts of Ly$\alpha$ redward of systemic due to scattering in the interstellar medium can boost Ly$\alpha$ transmission in the IGM during reionization. The relationship between Ly$\alpha$ velocity offset from systemic and other galaxy properties is not well-established at high-redshift or low luminosities, due to the difficulty of observing emission lines which trace systemic redshift. Rest-frame optical spectroscopy with JWST/NIRSpec has opened a new window into understanding of Ly$\alpha$ at z>3. We present a sample of 12 UV-faint galaxies ($-20 \lesssim$ MUV $\lesssim -16$) at $3 \lesssim z \lesssim 6$, with Ly$\alpha$ velocity offsets, $\Delta v_{\mathrm{Ly}\alpha}$, measured from VLT/MUSE and JWST/NIRSpec from the GLASS-JWST Early Release Program. We find median $\Delta v_{\mathrm{Ly}\alpha}$ of 205 km s$^{-1}$ and standard deviation 75 km s$^{-1}$, compared to 320 and 170km s$^{-1}$ for MUV < -20 galaxies in the literature. Our new sample demonstrates the previously observed trend of decreasing Ly$\alpha$ velocity offset with decreasing UV luminosity and optical line velocity dispersion, extends to MUV $\gtrsim$ -20, consistent with a picture where the Ly$\alpha$ profile is shaped by gas close to the systemic redshift. Our results imply that during reionization Ly$\alpha$ from UV-faint galaxies will be preferentially attenuated, but that detecting Ly$\alpha$ with low $\Delta v_{\mathrm{Ly}\alpha}$ can be an indicator of large ionized bubbles., Comment: 13 pages, 4 figures

Ly$\alpha$ emission is possibly the best indirect diagnostic of Lyman continuum (LyC) escape since the conditions that favor the escape of Ly$\alpha$ photons are often the same that allows for the escape of LyC photons. In this work, we present the rest UV-optical spectral characteristics of 11 Ly$\alpha$ emitting galaxies at 3 $<$ z $<$ 6 - the optimal redshift range chosen to avoid the extreme IGM attenuation while simultaneously studying galaxies close enough to the epoch of reionization. From a combined analysis of JWST/NIRSpec and MUSE data, we present the Ly$\alpha$ escape fraction and study their correlations with other physical properties of galaxies that might facilitate Ly$\alpha$ escape. We find that our galaxies have low masses (80\% of the sample with $\rm log_{10} \ M_{\star} < 9.5\ M_{\odot}$), compact sizes (median $\rm R_e \sim 0.7 \ kpc $), low dust content, moderate [OIII]/[OII] flux ratios (mean $\sim$ 6.8 $\pm$ 1.2), and moderate Ly$\alpha$ escape fraction (mean $\rm f_{esc}^{Ly\alpha} \ \sim$ 0.11). Our sample show characteristics that are broadly consistent with the low redshift galaxies with Ly$\alpha$ emission, which are termed as "analogs" of high redshift population. We predict the Lyman continuum escape fraction in our sample to be low (0.03-0.07), although larger samples in the post-reionization epoch are needed to confirm these trends., Comment: 15 pages, 4 figures, 1 table; Submitted to ApJL

We present a new analysis of the rest-frame UV and optical spectra of a sample of three $z>8$ galaxies discovered behind the gravitational lensing cluster RX\,J2129.4+0009. We combine these observations with $z>7.5$ galaxies from the literature, for which similar measurements are available. As already pointed out in other studies, the high [\oiii]$\lambda$5007/[\oii]$\lambda$3727 ratios ($O_{32}$) and steep UV continuum slopes ($\beta$) are consistent with the values observed for low redshift Lyman continuum emitters, suggesting that such galaxies contribute to the ionizing budget of the intergalactic medium. We construct a logistic regression model to estimate the probability of a galaxy being a Lyman continuum emitter based on the measured \MUV, $\beta$, and $O_{32}$. Using this probability and the UV luminosity function, we construct an empirical model that estimates the contribution of high redshift galaxies to reionization. The preferred scenario in our analysis shows that at $z\sim8$, the average escape fraction of the galaxy population (i.e., including both LyC emitters and non-emitters) varies with \MUV, with intermediate UV luminosity ($-19

Recent observations of GN-z11 with JWST have revealed a Ly$\alpha$ emission line with an equivalent width of 18$\pm 2$ angstroms. At z=10.6, this galaxy is expected to lie in the heart of reionization. We use a series of inhomogeneous reionization simulations to derive the distribution of the Ly$\alpha$ EW after traveling through the neutral intergalactic medium with varying average neutral gas fraction, $x_{HI}$. We use these distribution to place an upper limit of $x_{HI} < $ 0.88 at z=10.6 at 95% confidence level. We compare our upper limit to different reionization history models, which include the recently identified enhancement at the bright end of the luminosity function at z>8. We find that models in which faint galaxies have higher escape fraction compared to bright galaxies are favored by the new data., Comment: 7 pages, 3 figures

We develop a new method to determine the distance between a high-redshift galaxy and a foreground screen of atomic hydrogen. In a partially neutral universe, and assuming spherical symmetry, this equates to the radius of a ionized 'bubble' (R_B) surrounding the galaxy. The method requires an observed Lya equivalent width, its velocity offset from systemic, and an input Lya profile for which we adopt scaled versions of the profiles observed in low-z galaxies. We demonstrate the technique in a sample of 23 galaxies at z > 6, including eight at z = 7.2 - 10.6 recently observed with JWST. Our model estimates the emergent Lya properties, and the foreground distance to the absorbing IGM. We find that galaxies at z > 7.5 occupy smaller bubbles (~0.5 - 1 pMpc) than those at lower-z. With a relationship that is secure at 99% confidence, we empirically demonstrate the growth of ionized regions during the reionization epoch for the first time. We independently estimate the upper limit on the Str\"omgren radii (R_S), and derive the escape fraction of ionizing photons budget necessary for reionization., Comment: 8 pages, 7 figures. Submitted to ApJL after addressing referee comments

Exploiting broad- and narrow-band images of the Hubble Space Telescope from the near-UV to I-band rest frame, we study the star-forming clumps of six galaxies of the GASP sample undergoing strong ram pressure stripping. Clumps are detected in H alpha and near-UV, tracing star formation on different timescales. We consider clumps located in galaxy disks and stripped tails and formed in stripped gas but still close to the disk, called extraplanar. We detect 2406 H alpha-selected clumps (1708 in disks, 375 in extraplanar regions, and 323 in tails) and 3745 UV-selected clumps (2021 disk, 825 extraplanar, and 899 tail clumps). Only 15 per cent of star-forming clumps are spatially resolved, meaning that most are smaller than 140 pc. We study the luminosity and size distribution functions (LDFs and SDFs, respectively) and the luminosity-size relation. The average LDF slope is 1.79 +/- 0.09, while the average SDF slope is 3.1 +/- 0.5. The results suggest that the star formation is turbulence-driven and scale-free, as in main-sequence alaxies. All of the clumps, whether they are in the disks or tails, have an enhanced H alpha luminosity at a given size, compared to the clumps in main-sequence galaxies. Indeed, their H alpha luminosity is closer to that of clumps in starburst galaxies, indicating that ram pressure is able to enhance the luminosity. No striking differences are found among disk and tail clumps, suggesting that the different environments in which they are embedded play a minor role in influencing the star formation., Comment: 34 pages, 18 figures, 10 tables. Published on ApJ

We demonstrate how the stellar and nebular conditions in star-forming galaxies modulate the emission and spectral profile of HI Lya emission line. We examine the net Lya output, kinematics, and in particular emission of blue-shifted Lya radiation, using spectroscopy from with the Cosmic Origins Spectrograph on HST, giving a sample of 87 galaxies at redshift z=0.05-0.44. We contrast the Lya spectral measurements with properties of the ionized gas (from optical spectra) and stars (from stellar modeling). We demonstrate correlations of unprecedented strength between the Lya escape fraction (and equivalent width) and the ionization parameter (p~10^-15). The relative contribution of blue-shifted emission to the total Lya also increases from ~0 to ~40% over the range of O_32 ratios (p~10^-6). We also find particularly strong correlations with estimators of stellar age and nebular abundance, and weaker correlations regarding thermodynamic variables. Low ionization stage absorption lines suggest the Lya emission and line profile are predominantly governed by the column of absorbing gas near zero velocity. Simultaneous multi-parametric analysis over many variables shows we can predict 80% of the variance on Lya luminosity, and ~50% on the EW. We determine the most crucial predictive variables, finding that for tracers of the ionization state and Hb luminosity dominate the luminosity prediction whereas the Lya EW is best predicted by Hb EW and the Ha/Hb ratio. We discuss our results with reference to high redshift observations, focussing upon the use of Lya to probe the nebular conditions in high-z galaxies and cosmic reionization., Comment: Accepted by MNRAS. 27 pages and 19 figures. ~90 online-only Figures are available here: https://ttt.astro.su.se/~mhaye/t/onlineonly_resub.pdf

Galaxy formation and evolution are regulated by the feedback from galactic winds. Absorption lines provide the most widely available probe of winds. However, since most data only provide information integrated along the line-of-sight, they do not directly constrain the radial structure of the outflows. In this paper, we present a method to directly measure the gas electron density in outflows (ne), which in turn yields estimates of outflow cloud properties (e.g., density, volume filling-factor, and sizes/masses). We also estimate the distance (r) from the starburst at which the observed densities are found. We focus on 22 local star-forming galaxies primarily from the COS Legacy Archive Spectroscopic SurveY (CLASSY). In half of them, we detect absorption lines from fine structure excited transitions of Si II (i.e., Si II*). We determine ne from relative column densities of Si II and Si II*, given Si II* originates from collisional excitation by free electrons. We find that the derived ne correlates well with the galaxy's star-formation rate per unit area. From photoionization models or assuming the outflow is in pressure equilibrium with the wind fluid, we get r ~ 1 to 2 * rstar or ~ 5 * rstar, respectively, where rstar is the starburst radius. Based on comparisons to theoretical models of multi-phase outflows, nearly all of the outflows have cloud sizes large enough for the clouds to survive their interaction with the hot wind fluid. Most of these measurements are the first-ever for galactic winds detected in absorption lines and, thus, will provide important constraints for future models of galactic winds., Comment: 15 pages, 8 figures, 3 tables, submitted to ApJ

Star-forming, H$\alpha$-emitting clumps are found embedded in the gaseous tails of galaxies undergoing intense ram-pressure stripping in galaxy clusters, so-called jellyfish galaxies. These clumps offer a unique opportunity to study star formation under extreme conditions, in the absence of an underlying disk and embedded within the hot intracluster medium. Yet, a comprehensive, high spatial resolution study of these systems is missing. We obtained UVIS/HST data to observe the first statistical sample of clumps in the tails and disks of six jellyfish galaxies from the GASP survey; we used a combination of broad-band filters and a narrow-band H{\alpha} filter. HST observations are needed to study the sizes, stellar masses and ages of the clumps and their clustering hierarchy. These observations will be used to study the clump scaling relations, the universality of the star formation process and verify whether a disk is irrelevant, as hinted by jellyfish galaxy results. This paper presents the observations, data reduction strategy, and some general results based on the preliminary data analysis: the UVIS high spatial resolution gives an unprecedented sharp view of the complex structure of the inner regions of the galaxies and of the substructures in the galaxy disks; we found clear signatures of stripping in regions very close in projection to the galactic disk; the star-forming regions in the stripped tails are extremely bright and compact while we did not detect a significant number of star-forming clumps outside those detected by MUSE. The paper finally presents the development plan for the project., Comment: 17 pages, 7 figures. Accepted for publication in ApJ

A tight positive correlation between the stellar mass and the gas-phase metallicity of galaxies has been observed at low redshifts. The redshift evolution of this correlation can strongly constrain theories of galaxy evolution. The advent of JWST allows probing the mass-metallicity relation at redshifts far beyond what was previously accessible. Here we report the discovery of two emission-line galaxies at redshifts 8.15 and 8.16 in JWST NIRCam imaging and NIRSpec spectroscopy of targets gravitationally lensed by the cluster RXJ2129.4$+$0005. We measure their metallicities and stellar masses along with nine additional galaxies at $7.2 < z_{\rm spec} < 9.5$ to report the first quantitative statistical inference of the mass-metallicity relation at $z\approx8$. We measure $\sim 0.9$ dex evolution in the normalization of the mass-metallicity relation from $z \approx 8$ to the local Universe; at fixed stellar mass, galaxies are 8 times less metal enriched at $z \approx 8$ compared to the present day. Our inferred normalization is in agreement with the predictions of the FIRE simulations. Our inferred slope of the mass-metallicity relation is similar to or slightly shallower than that predicted by FIRE or observed at lower redshifts. We compare the $z \approx 8$ galaxies to extremely low metallicity analog candidates in the local Universe, finding that they are generally distinct from extreme emission-line galaxies or "green peas" but are similar in strong emission-line ratios and metallicities to "blueberry galaxies". Despite this similarity, at fixed stellar mass, the $z \approx 8$ galaxies have systematically lower metallicities compared to blueberry galaxies., Comment: Published in ApJ

Once only accessible in nearby galaxies, we can now study individual stars across much of the observable universe aided by galaxy-cluster gravitational lenses. When a star, compact object, or multiple such objects in the foreground galaxy-cluster lens become aligned, they can magnify a background individual star, and the timescale of a magnification peak can limit its size to tens of AU. The number and frequency of microlensing events therefore opens a window into the population of stars and compact objects, as well as high-redshift stars. To assemble the first statistical sample of stars in order to constrain the initial mass function (IMF) of massive stars at redshift z=0.7-1.5, the abundance of primordial black holes in galaxy-cluster dark matter, and the IMF of the stars making up the intracluster light, we are carrying out a 192-orbit program with the Hubble Space Telescope called "Flashlights," which is now two-thirds complete owing to scheduling challenges. We use the ultrawide F200LP and F350LP long-pass WFC3 UVIS filters and conduct two 16-orbit visits separated by one year. Having an identical roll angle during both visits, while difficult to schedule, yields extremely clean subtraction. Here we report the discovery of more than a dozen bright microlensing events, including multiple examples in the famous "Dragon Arc" discovered in the 1980s, as well as the "Spocks" and "Warhol" arcs that have hosted already known supergiants. The ultradeep observer-frame ultraviolet-through-optical imaging is sensitive to hot stars, which will complement deep James Webb Space Telescope infrared imaging. We are also acquiring Large Binocular Telescope LUCI and Keck-I MOSFIRE near-infrared spectra of the highly magnified arcs to constrain their recent star-formation histories.

We present a new selection of 358 blue compact dwarf galaxies (BCDs) from 5,000 square degrees in the Dark Energy Survey (DES), and the spectroscopic follow-up of a subsample of 68 objects. For the subsample of 34 objects with deep spectra, we measure the metallicity via the direct T$_e$ method using the auroral [\oiii]$\lambda$ 4363 emission line. These BCDs have an average oxygen abundance of 12+log(O/H)= 7.8, with stellar masses between 10$^7$ to 10$^8$ M$_\odot$ and specific star formation rates between $\sim$ 10$^{-9}$ to 10$^{-7}$ yr$^{-1}$. We compare the position of our BCDs with the Mass-metallicity (M-Z) and Luminosity-metallicity (L-Z) relation derived from the Local Volume Legacy sample. We find the scatter about the M-Z relation is smaller than the scatter about the L-Z relation. We identify a correlation between the offsets from the M-Z and L-Z relation that we suggest is due to the contribution of metal-poor inflows. Finally, we explore the validity of the mass-metallicity-SFR fundamental plane in the mass range probed by our galaxies. We find that BCDs with stellar masses smaller than $10^{8}$M$_{\odot}$ do not follow the extrapolation of the fundamental plane. This result suggests that mechanisms other than the balance between inflows and outflows may be at play in regulating the position of low mass galaxies in the M-Z-SFR space., Comment: 19 pages, 10 figures, accepted by ApJ

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides HST/UVIS F275W imaging for four CANDELS fields. We combine this UV imaging with existing HST/near-IR grism spectroscopy from 3D-HST$+$AGHAST to directly compare the resolved rest-frame UV and H$\alpha$ emission for a sample of 979 galaxies at $0.71.5$), suggesting that bursty star-formation is likely prevalent in the outskirts of even the most massive galaxies but is likely over-shadowed by their brighter cores. Furthermore, we present the UV-to-H$\alpha$ ratio as a function of galaxy surface brightness, a proxy for stellar mass surface density, and find that regions below $\sim$10$^{7.5}$ M$_\odot$ kpc$^{-2}$ are consistent with bursty star-formation, regardless of their galaxy stellar mass, potentially suggesting that local star-formation is independent of global galaxy properties at the smallest scales. Lastly, we find galaxies at $z>1.1$ to have bursty star-formation regardless of radius or surface brightness., Comment: 22 pages, 10 figures; accepted to ApJ

Ultraviolet light from early galaxies is thought to have ionized gas in the intergalactic medium. However, there are few observational constraints on this epoch because of the faintness of those galaxies and the redshift of their optical light into the infrared. We report the observation, in JWST imaging, of a distant galaxy that is magnified by gravitational lensing. JWST spectroscopy of the galaxy, at rest-frame optical wavelengths, detects strong nebular emission lines that are attributable to oxygen and hydrogen. The measured redshift is z = 9.51 +- 0.01, corresponding to 510 million years after the Big Bang. The galaxy has a radius of 16.2+4.6-7.2 parsecs, which is substantially more compact than galaxies with equivalent luminosity at z = 6 to 8, leading to a high star formation rate surface density.

Green pea galaxies are starbursting, low-mass galaxies that are good analogues to star-forming galaxies in the early universe. We perform a survey of 23 Green Peas using the MUSE Integral Field Unit spectrograph on the VLT to search for companion galaxies. The survey reaches an average point-source depth of $\sim 10^{-18}$ erg cm$^{-2}$ s$^{-1}$ for emission lines. The MUSE field of view allows us to probe a 1$\times$1 arcmin$^2$ field around these galaxies and to search their surroundings for faint companions that could have interacted with them and induced their starburst episodes. We search for companions using a variety of methods including template matching to emission and absorption line spectra. When restricting the search to the same physical area (R = 78 kpc) for all galaxies, we find that the fraction of green pea galaxies with companions is $0.11_{-0.05}^{+0.07}$. We define a control sample of star-forming galaxies with the same stellar masses and redshifts as the green peas, but consistent with the star-formation main sequence. We find that green pea galaxies are as likely to have companions as the control sample; for which the fraction of objects with companions is $0.08_{-0.03}^{+0.05}$. Given that we do not find statistical evidence for an elevated companion fraction in the green peas in this study, we argue that the ``companions" are likely unrelated to the bursts in these galaxies., Comment: Accepted for publication in the Astrophysical Journal

Galaxy Zoo: Clump Scout is a web-based citizen science project designed to identify and spatially locate giant star forming clumps in galaxies that were imaged by the Sloan Digital Sky Survey Legacy Survey. We present a statistically driven software framework that is designed to aggregate two-dimensional annotations of clump locations provided by multiple independent Galaxy Zoo: Clump Scout volunteers and generate a consensus label that identifies the locations of probable clumps within each galaxy. The statistical model our framework is based on allows us to assign false-positive probabilities to each of the clumps we identify, to estimate the skill levels of each of the volunteers who contribute to Galaxy Zoo: Clump Scout and also to quantitatively assess the reliability of the consensus labels that are derived for each subject. We apply our framework to a dataset containing 3,561,454 two-dimensional points, which constitute 1,739,259 annotations of 85,286 distinct subjects provided by 20,999 volunteers. Using this dataset, we identify 128,100 potential clumps distributed among 44,126 galaxies. This dataset can be used to study the prevalence and demographics of giant star forming clumps in low-redshift galaxies. The code for our aggregation software framework is publicly available at: https://github.com/ou-astrophysics/BoxAggregator, Comment: 31 pages, 22 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We present calculations of ultraviolet spectra resulting from the scattering of photons by gas in-falling onto an isotropically emitting source of radiation. The model is based on an adaptation of the semi-analytical line transfer (SALT) code of Scarlata & Panagia (2015), and designed to interpret the inverse P-Cygni profiles observed in the spectra of partially ionized galactic inflows. In addition to presenting the model, we explore the parameter space of the inflowing SALT model and recreate various physically motivated scenarios including spherical inflows, inflows with covering fractions less than unity, and galactic fountains (i.e., galactic systems with both an inflowing and outflowing component). The resulting spectra from inflowing gas show spectral features that could be misinterpreted as ISM features in low resolution spectroscopy ($\sigma \approx 120$ $\rm{km }$ $\rm{s}^{-1}$), suggesting that the total number of galactic systems with inflows is undercounted. Our models suggest that observations at medium resolution ($R = 6000$ or $\sigma \approx 50$ $\rm{km }$ $\rm{s}^{-1}$) that can be obtained with 8m-class telescopes will be able to resolve the characteristic inverse P Cygni profiles necessary to identify inflows., Comment: Accepted by ApJ

Absorption line spectroscopy offers one of the best opportunities to constrain the properties of galactic outflows and the environment of the circumgalactic medium. Extracting physical information from line profiles is difficult, however, for the physics governing the underlying radiation transfer is complicated and depends on many different parameters. Idealized analytical models are necessary to constrain the large parameter spaces efficiently, but are typically plagued by model degeneracy and systematic errors. Comparison tests with idealized numerical radiation transfer codes offer an excellent opportunity to confront both of these issues. In this paper, we present a detailed comparison between SALT, an analytical radiation transfer model for predicting UV spectra of galactic outflows, with the numerical radiation transfer software, RASCAS. Our analysis has lead to upgrades to both models including an improved derivation of SALT and a customizable adaptive mesh refinement routine for RASCAS. We explore how well SALT, when paired with a Monte Carlo fitting procedure, can recover flow parameters from non-turbulent and turbulent flows. When the velocity and density gradients are excluded, we find that flow parameters are well recovered from high resolution (20 $\rm{km}$ $\rm{s}^{-1}$) data and moderately well from medium resolution (100 $\rm{km}$ $\rm{s}^{-1}$) data without turbulence at a S/N = 10, while derived quantities (e.g., mass outflow rates, column density, etc.) are well recovered at all resolutions. In the turbulent case, biased errors emerge in the recovery of individual parameters, but derived quantities are still well recovered.

We obtained 838 Sloan r-band images (~28 hrs) of the GOODS-North field with the Large Binocular Camera (LBC) on the Large Binocular Telescope in order to study the presence of extended, low surface brightness features in galaxies and investigate the trade-off between image depth and resolution. The individual images were sorted by effective seeing, which allowed for optimal resolution and optimal depth mosaics to be created with all images with seeing FWHM < 0.9" and FWHM < 2.0", respectively. Examining bright galaxies and their substructure as well as accurately deblending overlapping objects requires the optimal resolution mosaic, while detecting the faintest objects possible (to a limiting magnitude of $m_{AB}$ ~ 29.2 mag) requires the optimal depth mosaic. The better surface brightness sensitivity resulting from the larger LBC pixels, compared to those of extant WFC3/UVIS and ACS/WFC cameras aboard the Hubble Space Telescope (HST) allows for unambiguous detection of both diffuse flux and very faint tidal tails. We created azimuthally-averaged radial surface brightness profiles for the 360 brightest galaxies in the mosaics. We find little difference in the majority of the light profiles from the optimal resolution and optimal depth mosaics. However, $\lesssim$ 15% of the profiles show excess flux in the galaxy outskirts down to surface brightness levels of $\mu^{AB}_{r} $ $\simeq$ 31 mag arcsec $^{-2}$. This is relevant to Extragalactic Background Light (EBL) studies as diffuse light in the outer regions of galaxies are thought to be a major contribution to the EBL. While some additional diffuse light exists in the optimal depth profiles compared to the shallower, optimal resolution profiles, we find that diffuse light in galaxy outskirts is a minor contribution to the EBL overall in the r-band., Comment: 30 pages, 25 figures, submitted to PASP

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

Abell 1689 is a well studied cluster of galaxies and one of the largest gravitational lens systems ever observed. We have obtained a reconstruction of the cluster Abell 1689 using Grale, a free-form lens inversion method that relies exclusively on the multiple image data. Non-inclusion of any data related to cluster member galaxies ensures an unbiased measure of the mass distribution, which is the most notable feature of free-form methods like Grale. We used two different sets of multiple image systems from the available strong lensing data - one containing only the secure systems (107 images), and the other containing all available systems, only excluding some very non-secure systems (151 images). For the very well-constrained central $\sim$100 kpc region of the cluster we made detailed comparison of the Grale reconstructed lensing mass and stellar mass retrieved by the Spectral Energy Distribution (SED) fitting software FAST++. We found a light-unaccompanied mass peak in this region, whose existence, while tentative, is favored by the distribution of nearby images that are local maxima in the Fermat potential. However, further tests, using different methodologies are needed to confirm the reality of this feature. If it shown to be real, this light-unaccompanied mass peak is consistent with dark matter self-interaction cross-section $\sigma \lesssim 1$cm$^2$/g, while being in tension with larger cross-sections., Comment: 14 pages, 11 figures

The COS Legacy Archive Spectroscopic SurveY (CLASSY) is designed to provide the community with a spectral atlas of 45 nearby star-forming galaxies which were chosen to cover similar properties as those seen at high-z (z>6). The prime high level science product of CLASSY is accurately coadded UV spectra, ranging from ~1000-2000A, derived from a combination of archival and new data obtained with HST's Cosmic Origins Spectrograph (COS). This paper details the multi-stage technical processes of creating this prime data product, and the methodologies involved in extracting, reducing, aligning, and coadding far-ultraviolet (FUV) and near-ultraviolet (NUV) spectra. We provide guidelines on how to successfully utilize COS observations of extended sources, despite COS being optimized for point sources, and best-practice recommendations for the coaddition of UV spectra in general. Moreover, we discuss the effects of our reduction and coaddition techniques in the scientific application of the CLASSY data. In particular, we find that accurately accounting for flux calibration offsets can affect the derived properties of the stellar populations, while customized extractions of NUV spectra for extended sources are essential for correctly diagnosing the metallicity of galaxies via CIII] nebular emission. Despite changes in spectral resolution of up to ~25% between individual datasets (due to changes in the COS line spread function), no adverse affects were observed on the difference in velocity width and outflow velocities of isolated absorption lines when measured in the final combined data products, owing in-part to our signal-to-noise regime of S/N<20., Comment: 31 pages, 17 figures, accepted for publication in ApJS

We explore how the fraction of quenched galaxies changes in groups of galaxies with respect to the distance to the center of the group, redshift, and stellar mass to determine the dominant process of environmental quenching in $0.2 < z < 0.8$ groups. We use new UV data from the UVCANDELS project in addition to existing multiband photometry to derive new galaxy physical properties of the group galaxies from the zCOSMOS 20k Group Catalog. Limiting our analysis to a complete sample of log$(M_*/M_{\odot})>10.56$ group galaxies we find that the probability of being quenched increases slowly with decreasing redshift, diverging from the stagnant field galaxy population. A corresponding analysis on how the probability of being quenched increases with time within groups suggests that the dominant environmental quenching process is characterized by slow ($\sim$Gyr) timescales. We find a quenching time of approximately $4.91^{+0.91}_{-1.47} $Gyrs, consistent with the slow processes of strangulation (Larson et al. 1980) and delayed-then-rapid quenching (Wetzel et al. 2013 arXiv:1206.3571v2 [astro-ph.CO]), although more data are needed to confirm this result.