Search

Stellar velocity dispersion ($\sigma$) of massive elliptical galaxies is a key ingredient to breaking the mass-sheet degeneracy and obtaining precise and accurate cosmography from gravitational time delays. The relative uncertainty on the Hubble constant H$_0$ is double the relative error on $\sigma$. Therefore, time-delay cosmography imposes much more demanding requirements on the precision and accuracy of $\sigma$ than galaxy studies. While precision can be achieved with an adequate signal-to-noise ratio (SNR), accuracy depends critically on factors such as elemental abundances and temperature of stellar templates, flux calibration, and wavelength ranges. We carry out a detailed study of the problem using multiple sets of galaxy spectra of massive elliptical galaxies with SNR$\sim$30-160 \AA$^{-1}$, and state-of-the-art empirical and semi-empirical stellar libraries and stellar population synthesis templates. We show that the choice of stellar library is generally the dominant source of residual systematic error. We propose a general recipe to mitigate and account for residual uncertainties. We show that sub-percent accuracy can be achieved on individual spectra with our data quality. Covariance between velocity dispersions measured for a sample of spectra can also be reduced to sub-percent levels. We recommend this recipe for all applications that require high precision and accurate stellar kinematics, and make all software publicly available to facilitate its implementation. This recipe will be used in future TDCOSMO collaboration papers., Comment: Submitted to A&A

The electron density (${n_{\rm e}}$) of the interstellar medium (ISM) in star-forming galaxies is intimately linked to star formation and ionization condition. Using the high-resolution spectra obtained from the JWST NIRSpec micro shutter assembly (MSA) as part of the GLASS-JWST program, we have assembled the largest sample to date (34 galaxies) with individual ${n_{\rm e}}$ measurements derived from the [OII] $\lambda\lambda$3726,29 and/or [SII] $\lambda\lambda$6718,32 doublets at $0.7\lesssim z\lesssim 9.3$. The gravitational lensing magnification by the foreground Abell~2744 cluster allows us to probe ${n_{\rm e}}$ in galaxies with stellar masses ($M_{*}$) down to $\simeq 10^{7.5} M_\odot$ across the entire redshift range. Our analysis reveals that the [OII] flux ratios are marginally anti-correlated with specific star formation rate (sSFR) within a 1-$\sigma$ confidence interval, whereas the [SII] flux ratios show no significant correlation with sSFR. Despite clear correlation between sSFR and redshift within our sample, we find no apparent redshift evolution of ${n_{\rm e}}$ at $z \simeq 1-9$. Our dataset also includes 13 galaxies where ${n_{\rm e}}$ can be measured from both [OII] and [SII]. Contrary to findings at lower redshifts, we observe considerable scatter in ${n_{\rm e}}$ measurements from [OII] and [SII], indicating a complex gaseous environment with significant variations in ${n_{\rm e}}$ in high-redshift galaxies. This work highlights the unique capability of JWST NIRSpec/MSA high-resolution spectroscopy to characterize the detailed physical properties of the ISM in individual high-redshift galaxies., Comment: 15 pages, 7 figures, 2 tables; Accepted for publication in ApJL

We present measurements of the gas-phase Oxygen and Nitrogen abundances obtained by applying the direct method to JWST NIRspec $R\sim1000$ spectroscopy for 6 galaxies at redshift greater than 3. Our measurements are based on rest-frame optical Nitrogen [N II]$_{\lambda\lambda6548,6583}$ lines and are complemented by 6 additional objects from the literature at $3\leq z \leq 6$. We find that 9 out of 12 objects have values of log(N/O) that are compatible with those found for low-redshift, metal-poor, dwarf galaxies and for HII regions of more luminous local galaxies. However, 3 out of 12 objects have log(N/O) values that are overabundant compared to what is expected on the basis of their Oxygen abundance. We explore a few standard scenarios to explain the observations and conclude that, within the limited statistics available to us, none of them can be definitely excluded even though we prefer dilution by pristine gas infall in between star formation bursts as this is predicted by simulations to take place as a natural part of bursty star formation., Comment: Submitted to ApJ

We measure the spectral properties of a sample of 20 galaxies at z~0.35 selected for having surprisingly red JWST/NIRCAM F200W-F444W colors. 19 galaxies were observed with JWST/NIRSpec in the PRISM configuration, while one galaxy was observed with the high resolution gratings. 17/20 galaxies in our sample exhibit strong 3.3 $\mu m$ polycyclic aromatic hydrocarbon (PAH$_{3.3}$) emission (equivalent width EW(PAH$_{3.3}$)$>0.03\mu m$). In these galaxies, the strength of the color excess does not depend on environment and it correlates with EW(PAH$_{3.3}$). Nonetheless, the presence of the PAH$_{3.3}$ alone can not fully explain the color excess, as an equivalent width of ~0.1$\mu m$ is able to increase the color of galaxies by only 0.13 mag. A contribution from a hot dust component is required to explain the excess. Both the EW(PAH$_{3.3}$) and flux correlate with the H$\alpha$ equivalent width and flux, suggesting that they are produced by the same mechanism. 5/20 galaxies showing PAH would be classified as passive based on broad band rest frame colors ((B-V) and/or UVJ diagrams) and are hence "faux-passive". Of these, 3 galaxies have a significantly lower EW(PAH$_{3.3}$) given their color and also have low EW(H$\alpha$) and we tentatively conclude this behaviour is due to the presence of an AGN. The three galaxies with no PAH$_{3.3}$ in emission have passive spectra, as do the 8 galaxies in our sample with normal F200W-F444W colors. We therefore conclude that the PAH$_{3.3}$ feature is linked to dust-enshrouded star formation. The dust corrected SFR from PAH$_{3.3}$ is a factor of 3.5 higher than the SFR obtained from H$\alpha$, suggesting that these galaxies are characterized by significant amounts of dust., Comment: A&A in press

We introduce the Bias-free Extragalactic Analysis for Cosmic Origins with NIRCam (BEACON) survey, a JWST Cycle2 program allocated up to 600 pure-parallel hours of observations. BEACON explores high-latitude areas of the sky with JWST/NIRCam over $\sim100$ independent sightlines, totaling $\sim0.3$deg$^2$, reaching a median F444W depth of $\approx28.2$AB mag (5$\sigma$). Based on existing JWST observations in legacy fields, we estimate that BEACON will photometrically identify 25--150 galaxies at $z>10$ and 500--1000 at $z\sim7$--10 uniquely enabled by an efficient multiple filter configuration spanning $0.9$--5.0$\mu$m. The expected sample size of $z>10$ galaxies will allow us to obtain robust number density estimates and to discriminate between different models of early star formation. In this paper, we present an overview of the survey design and initial results using the first 19 fields. We present 129 galaxy candidates at $z>7$ identified in those fields, including 11 galaxies at $z>10$ and several UV-luminous ($M_{\rm UV}<-21$mag) galaxies at $z\sim8$. The number densities of $z<13$ galaxies inferred from the initial fields are overall consistent with those in the literature. Despite reaching a considerably large volume ($\sim10^5$Mpc$^3$), however, we find no galaxy candidates at $z>13$, providing us with a complimentary insight into early galaxy evolution with minimal cosmic variance. We publish imaging and catalog data products for these initial fields. Upon survey completion, all BEACON data will be coherently processed and distributed to the community along with catalogs for redshift and other physical quantities., Comment: Submitted to ApJ; DR1 data release will be made on the team website (https://beacon-jwst.github.io); Fig. 8 has been updated

The origin of the tight scaling relation between the mass of supermassive black holes (SMBHs; $M_{\rm BH}$) and their host-galaxy properties remains unclear. Active galactic nuclei (AGNs) probe phases of ongoing SMBH growth and offer the only opportunity to measure $M_{\rm BH}$ beyond the local Universe. However, determining AGN host galaxy stellar velocity dispersion $\sigma_\star$, and their galaxy dynamical masses $M_{\rm dyn}$, is complicated by AGN contamination, aperture effects and different host galaxy morphologies. We select a sample of AGNs for which $M_{\rm BH}$ has been independently determined to high accuracy by state-of-the-art techniques: dynamical modeling of the reverberation signal and spatially resolving the broad-line region with VLTI/GRAVITY. Using IFU observations, we spatially map the host galaxy stellar kinematics across the galaxy and bulge effective radii. We find that that the dynamically hot component of galaxy disks correlates with $M_{\rm BH}$; however, the correlations are tightest for aperture-integrated $\sigma_\star$ measured across the bulge. Accounting for the different $M_{\rm BH}$ distributions, we demonstrate - for the first time - that AGNs follow the same $M_{\rm BH}$-$\sigma_\star$ and $M_{\rm BH}$-$M_{\rm bulge, dyn}$ relations as quiescent galaxies. We confirm that the classical approach of determining the virial factor as sample-average, yielding ${\rm log }f= 0.65 \pm 0.18$, is consistent with the average $f$ from individual measurements. The similarity between the underlying scaling relations of AGNs and quiescent galaxies implies that the current AGN phase is too short to have altered BH masses on a population level. These results strengthen the local calibration of $f$ for measuring single-epoch $M_{\rm BH}$ in the distant Universe., Comment: Submitted to ApJ. 40 pages, 19 figures, 6 tables. Reuploaded after acceptance for publication

We present JWST/NIRSpec PRISM spectroscopic characterization of GHZ9 at z= 10.145 $\pm$ 0.010, currently the most distant source detected by the Chandra X-ray Observatory. The spectrum reveals several UV high-ionization lines, including CII, SiIV, [NIV], CIV, HeII, OIII], NIII], and CIII]. The prominent rest-frame equivalent widths (EW(CIV)$\simeq$65A, EW(HeII)$\simeq$18A, EW(CIII])$\simeq$48A) show the presence of a hard radiation field, while the analysis of line ratio diagnostics suggest this galaxy hosts both AGN and star-formation activity. GHZ9 is nitrogen-enriched (6--9.5 times solar), carbon-poor (0.2--0.65 times solar), metal-poor (Z = 0.01--0.1 Z$_{\odot}$), and compact ($<$ 106 pc), similarly to GNz11, GHZ2, and recently discovered N-enhanced high redshift objects. We exploited the newly available JWST/NIRSpec and NIRCam dataset to perform an independent analysis of the Chandra data confirming that GHZ9 is the most likely JWST source associated to X-ray emission at 0.5-7 keV. Assuming a spectral index $\Gamma$ = 2.3 (1.8), we estimate a black hole (BH) mass of 1.60 $\pm$ 0.31 (0.48 $\pm$ 0.09) $\times$ 10$^8$M$_{\odot}$, which is consistent either with Eddington-accretion onto heavy ($\geq$ 10$^6$ M$_{\odot}$) BH seeds formed at z=18, or super-Eddington accretion onto a light seed of $\sim$ 10$^2-10^4$ M$_{\odot}$ at z = 25. The corresponding BH-to-stellar mass ratio M$_{BH}$/M$_{star}$= 0.33$\pm$0.22 (0.10$\pm$0.07), with a stringent limit $>$0.02, implies an accelerated growth of the BH mass with respect to the stellar mass. GHZ9 is the ideal target to constrain the early phases of AGN-galaxy coevolution with future multi-frequency observations., Comment: Submitted to ApJ

Strongly lensed quasars can be used to constrain cosmological parameters through time-delay cosmography. Models of the lens masses are a necessary component of this analysis. To enable time-delay cosmography from a sample of $\mathcal{O}(10^3)$ lenses, which will soon become available from surveys like the Rubin Observatory's Legacy Survey of Space and Time (LSST) and the Euclid Wide Survey, we require fast and standardizable modeling techniques. To address this need, we apply neural posterior estimation (NPE) for modeling galaxy-scale strongly lensed quasars from the Strong Lensing Insights into the Dark Energy Survey (STRIDES) sample. NPE brings two advantages: speed and the ability to implicitly marginalize over nuisance parameters. We extend this method by employing sequential NPE to increase precision of mass model posteriors. We then fold individual lens models into a hierarchical Bayesian inference to recover the population distribution of lens mass parameters, accounting for out-of-distribution shift. After verifying our method using simulated analogs of the STRIDES lens sample, we apply our method to 14 Hubble Space Telescope single-filter observations. We find the population mean of the power-law elliptical mass distribution slope, $\gamma_{\text{lens}}$, to be $\mathcal{M}_{\gamma_{\text{lens}}}=2.13 \pm 0.06$. Our result represents the first population-level constraint for these systems. This population-level inference from fully automated modeling is an important stepping stone towards cosmological inference with large samples of strongly lensed quasars.

Strongly lensed supernovae (SNe) are a rare class of transient that can offer tight cosmological constraints that are complementary to methods from other astronomical events. We present a follow-up study of one recently-discovered strongly lensed SN, the quadruply-imaged Type Ia SN 2022qmx (aka, "SN Zwicky") at z = 0.3544. We measure updated, template-subtracted photometry for SN Zwicky and derive improved time delays and magnifications. This is possible because SNe are transient, fading away after reaching their peak brightness. Specifically, we measure point spread function (PSF) photometry for all four images of SN Zwicky in three Hubble Space Telescope WFC3/UVIS passbands (F475W, F625W, F814W) and one WFC3/IR passband (F160W), with template images taken $\sim 11$ months after the epoch in which the SN images appear. We find consistency to within $2\sigma$ between lens model predicted time delays ($\lesssim1$ day), and measured time delays with HST colors ($\lesssim2$ days), including the uncertainty from chromatic microlensing that may arise from stars in the lensing galaxy. The standardizable nature of SNe Ia allows us to estimate absolute magnifications for the four images, with images A and C being elevated in magnification compared to lens model predictions by about $6\sigma$ and $3\sigma$ respectively, confirming previous work. We show that millilensing or differential dust extinction is unable to explain these discrepancies and find evidence for the existence of microlensing in images A, C, and potentially D, that may contribute to the anomalous magnification., Comment: Submitted to ApJ

We obtain spatially resolved kinematics with the Keck Cosmic Web Imager (KCWI) integral-field spectrograph for a sample of 14 massive (11 < log10 M* < 12) lensing early-type galaxies (ETGs) at redshifts z=0.15-0.35 from the Sloan Lens ACS (SLACS) survey. We integrate within the galaxy effective radius and examine the rotational and dispersion velocities, showing that 11/14 are quantitatively classified as slow rotators in comparison with local galaxy surveys. Of key interest is the ability of this data to enable the precision required for cosmological inference with lensing time delays on scales of 1-2% uncertainty. The dataset is unprecedented for galaxy-scale lens galaxies, in terms of signal-to-noise ratio, sampling, and calibration. We test sources of systematic error and identify primary contributions from choice of stellar template library and wavelength range of the spectral fit. Systematics are quantified at the spatial bin level, resulting in systematic error at 3% and positive spatial covariance of 2%. We examine the effects of integration of the kinematic maps within circular apertures of different sizes and compare with SDSS single-aperture velocity dispersions. The most recent velocity dispersion estimates from SDSS spectra are found to be biased by a factor of 5.3% with respect to KCWI data, and to underestimate uncertainties. We examine correlations between scaling relations and show the correlations to agree with previous SLACS analysis with no statistically significant disagreement. A follow-up paper will present Jeans modeling and discuss the context of these observations within broader studies of galaxy evolution and cosmology., Comment: 33 pages, 14 figures, submitted for publication in ApJ

Understanding the rapid formation of supermassive black holes (SMBHs) in the early universe requires an understanding of how stellar mass grows in the host galaxies. Here, we perform an analysis of rest-frame optical spectra and imaging from JWST of two quasar host galaxies at z>6 which exhibit Balmer absorption lines. These features in the stellar continuum indicate a lack of young stars, similar to low-redshift post-starburst galaxies whose star formation was recently quenched. We find that the stellar mass (log(M_* / M_sun) > 10.6) of each quasar host grew in a starburst episode at redshift 7 or 8. One of the targets exhibits little ongoing star formation, as evidenced by the photometric signature of the Balmer break and a lack of spatially resolved H-alpha emission, placing it well below the star formation main sequence at z = 6. The other galaxy is transitioning to a quiescent phase; together, the two galaxies represent the most distant massive post-starburst galaxies known. The maturity of these two galaxies is further supported by the stellar velocity dispersions of their host galaxies, placing them slightly above the upper end of the local M_BH - sigma_* relation. The properties of our two post-starburst galaxies, each hosting an active SMBH with log(M_BH / M_sun) > 9, suggests that black holes played a major role in shaping the formation of the first massive galaxies in the Universe., Comment: 24 pages, 7 figures, submitted to a Nature journal

The radial gradient of gas-phase metallicity is a powerful probe of the chemical and structural evolution of star-forming galaxies, closely tied to disk formation and gas kinematics in the early universe. We present spatially resolved chemical and dynamical properties for a sample of 25 galaxies at $0.5 \lesssim z \lesssim 1.7$ from the \msasd survey. These innovative observations provide 3D spectroscopy of galaxies at a spatial resolution approaching JWST's diffraction limit and a high spectral resolution of $R\simeq2700$. The metallicity gradients measured in our galaxy sample range from $-$0.03 to 0.02 dex~kpc$^{-1}$. Most galaxies exhibit negative or flat radial gradients, indicating lower metallicity in the outskirts or uniform metallicity throughout the entire galaxy. We confirm a tight relationship between stellar mass and metallicity gradient at $z\sim1$ with small intrinsic scatter of 0.02 dex~kpc$^{-1}$. Our results indicate that metallicity gradients become increasingly negative as stellar mass increases, likely because the more massive galaxies tend to be more ``disky". This relationship is consistent with the predictions from cosmological hydrodynamic zoom-in simulations with strong stellar feedback. This work presents the effort to harness the multiplexing capability of JWST NIRSpec/MSA in slit-stepping mode to map the chemical and kinematic profiles of high-redshift galaxies in large samples and at high spatial and spectral resolution., Comment: 24 pages, 6 figures, 1 table; Accepted for publication in ApJL

We report the identification of two galaxy overdensities at $z\sim5.7$ in the sightline of the galaxy cluster Abell 2744. These overdensities consist of 25 and 17 member galaxies, spectroscopically confirmed with JWST NIRSpec/MSA and NIRCam/WFSS. Each overdensity has a total stellar mass of $\sim2\times10^{10} M_\odot$ and a star formation rate of $\sim200 M_\odot$/yr within a central region of radius $R=2$ Mpc (physical). The sensitive PRISM spectra allow us to identify six galaxies that show weak Ha+[NII] emissions within the overdensities ($27\pm6\%$), whereas the fraction of such galaxies is found significantly lower ($6\pm2\%$) in field samples of the equivalent redshift range. These weak emission line galaxies, dubbed as wELGs, exhibit a strong continuum break at $4000$AA rest-frame, a characteristic feature of evolved stellar populations. The high observed fraction of wELGs in the two overdensities is consistent with the idea that high-density environments are an ideal site where galaxies can accelerate their evolutionary pace compared to field analogs. Our study pinpoints an early onset of environmental effects, already important within one billion years after the Big Bang, and provides a complementary perspective on the emergence of quenched, massive galaxies at lower redshifts. Potential contributions from black hole accretion feedback to the reduction of star formation activity are discussed, but the connection to the local environments remains unclear., Comment: Submitted to ApJ; Accepted for publication

We present a new measurement of the dark and luminous matter distribution of massive elliptical galaxies, and their evolution with redshift, by combining strong lensing and dynamical observables. Our sample of 58 lens galaxies covers a redshift range of $0.090\leq z_{\rm l}\leq0.884$. By combining new Hubble Space Telescope imaging with previously observed velocity dispersion and line-of-sight measurements, we decompose the luminous matter profile from the dark matter profile and perform a Bayesian hierarchical analysis to constrain the population-level properties of both profiles. We find that the inner slope of the dark matter density profile ("cusp"; $\rho_{\rm DM}\propto r^{-\gamma_{\rm in}}$) is slightly steeper ($\mu_{\gamma_{\rm in}}=1.18^{+0.03}_{-0.03}$ at $z=0.35$ with $\leq0.16$ intrinsic scatter) than a standard Navarro$-$Frenk$-$White (NFW; $\gamma_{\rm in}=1$), with an appreciable evolution with redshift ($d\log(\gamma_{\rm in})/dz=-0.33\pm0.13$) and is consistent with NFW-like distributions at higher redshifts ($z\geq0.56$ for $\leq1\sigma$ consistency). Additionally, we find the stellar mass-to-light ratio at the population level consistent with that of a Salpeter initial mass function, a small stellar mass-to-light gradient ($\kappa_{*}(r)\propto r^{-\eta}$, with $\overline{\eta}\leq9.4\times10^{-3}$), and isotropic stellar orbits. Our averaged total mass density profile is consistent with a power-law profile within $0.25-4$ Einstein radii ($\overline{\gamma}=2.14\pm0.06$), with an internal mass-sheet transformation parameter $\overline{\lambda}=1.02\pm0.01$ consistent with no mass sheet. Our findings confirm the validity of the standard mass models used for time-delay cosmography. However, our results are in strong tension with predictions from hydrodynamical simulations such as IllustrisTNG, highlighting the need to better understand the formation of massive galaxies., Comment: 28 pages, 20 figures

Integral field spectroscopy (IFS) is a powerful tool for understanding the formation of galaxies across cosmic history. We present the observing strategy and first results of MSA-3D, a novel JWST program using multi-object spectroscopy in a slit-stepping strategy to produce IFS data cubes. The program observed 43 normal star-forming galaxies at redshifts $0.5 \lesssim z \lesssim 1.5$, corresponding to the epoch when spiral thin-disk galaxies of the modern Hubble sequence are thought to emerge, obtaining kpc-scale maps of rest-frame optical nebular emission lines with spectral resolution $R\simeq2700$. Here we describe the multiplexed slit-stepping method which is $>15$ times more efficient than the NIRSpec IFS mode for our program. As an example of the data quality, we present a case study of an individual galaxy at $z=1.104$ (stellar mass $M_{*} = 10^{10.3}~M_{\odot}$, star formation rate~$=3~M_{\odot}$ yr$^{-1}$) with prominent face-on spiral structure. We show that the galaxy exhibits a rotationally supported disk with moderate velocity dispersion ($\sigma = 36^{+5}_{-4}$~\kms), a negative radial metallicity gradient ($-0.020\pm0.002$~dex\,kpc$^{-1}$), a dust attenuation gradient, and an exponential star formation rate density profile which closely matches the stellar continuum. These properties are characteristic of local spirals, indicating that mature galaxies are in place at $z\sim1$. We also describe the customized data reduction and original cube-building software pipelines which we have developed to exploit the powerful slit-stepping technique. Our results demonstrate the ability of JWST slit-stepping to study galaxy populations at intermediate to high redshifts, with data quality similar to current surveys of the $z\sim0.1$ universe., Comment: 24 pages, 12 figures, 1 table; submitted to ApJ, comments welcome. Our custom designed data reduction pipeline (MSA3D) and example processed data-cube will be publicly released before the JWST Cycle 4 deadline. MSA3D package on GitHub: https://github.com/barisiciv/msa3d.git

We present new results on the rest-frame UV luminosity function (UVLF) and stellar mass-to-light (M/L) ratio of bright (M$_{\rm UV}\lesssim-20$ mag) spectroscopically-confirmed galaxies at $z=7-9$ derived from the BoRG-$JWST$ survey, a unique data set of NIRSpec prism follow up of $HST$-selected sources from random-pointing imaging. By selecting galaxies from over 200 independent sight lines, the survey minimizes cosmic variance ensuring a statistically robust sample of the bright-galaxy population during the epoch of reionization. The data is used to constrain, for the first time, the bright end of the UVLF at $z=7-9$ from spectroscopically-confirmed galaxies over eight independent fields. We find that the bright end of the UVLF is higher than found using imaging over $JWST$ legacy fields, suggesting the latter may be significantly affected by cosmic variance, and thus reducing the tension with recent findings from $JWST$ at $z>10$ and comparable to models invoking little dust attenuation and bursty star formation. Additionally, we use the galaxies' $JWST$ spectra to infer their stellar masses and M/L ratios relative to other $HST$ and $JWST$ studies. We show that the stellar mass scales almost linearly with UV luminosity (M$_* \propto L_{\rm UV}^{0.85\pm0.12}$), albeit with large ($\sim0.5$ dex) intrinsic scatter, consistent with stochastic bursts of star formation in early galaxy formation., Comment: 13 pages, 4 figures, 2 tables

We present the BoRG-JWST survey, a combination of two JWST Cycle 1 programs aimed at obtaining NIRSpec spectroscopy of representative, UV-bright $75$, highlighting the large abundance of high-redshift galaxies even in individual WFC3 pointings. The primary sample span an absolute magnitude range $-20.47$ sources currently known and comparable to the brightest sources at $z>10$. Prominent [O III]+H$\beta$ lines are found across the full sample, while a stack of sources reveals a plethora of other rest-optical lines and additional rest-UV C III]1909 \r{A} emission. Despite their luminosities, none of the low-resolution spectra display evidence for Type 1 AGN activity based on a search for broad-line emission. Lastly, we present a spectroscopic data release of 188 confirmed $0.5\lesssim z\lesssim5.0$ sources from filler MSA observations, highlighting the legacy value of the survey and a representative benchmark for comparisons to deep field observations., Comment: Submitted to ApJ, comments welcome

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038-4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain $D_{\Delta t} = 1.68^{+0.40}_{-0.38}$ Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat $\Lambda$CDM cosmology, we constrain the Hubble constant to be $H_{0} = 65^{+23}_{-14}\, \rm km\ s^{-1}\,Mpc^{-1}$. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory's Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future., Comment: 8 pages, 5 figures, 3 tables; accepted for publication in Astronomy & Astrophysics

This is the second in a series of papers in which we use JWST MIRI multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of 9 systems. The warm dust region is compact and sensitive to perturbations by populations of halos down to masses $\sim 10^6$ M$_{\odot}$. Using these warm dust flux-ratio measurements in combination with 5 previous narrow-line flux-ratio measurements, we constrain the halo mass function. In our model, we allow for complex deflector macromodels with flexible third and fourth-order multipole deviations from ellipticity, and we introduce an improved model of the tidal evolution of subhalos. We constrain a WDM model and find an upper limit on the half-mode mass of $10^{7.6} M_\odot$ at posterior odds of 10:1. This corresponds to a lower limit on a thermally produced dark matter particle mass of 6.1 keV. This is the strongest gravitational lensing constraint to date, and comparable to those from independent probes such as the Ly$\alpha$ forest and Milky Way satellite galaxies.

A quadruply lensed source, J125856.3-031944, has been discovered using the DELVE survey and WISE W1 - W2 colors. Followup direct imaging carried out with the MPIA 2.2 m and the Baade 6.5 m telescopes is analyzed, as is spectroscopy from the 2.5 m Nordic Optical Telescope. The lensed image configuration is kite-like, with the faintest image 2 magnitudes fainter than the other three. Redward of 6000 AA that image is badly blended with the lensing galaxy, which is elongated along the symmetry axis of the kite. Magellan direct imaging carried out in Sloan g permits better deblending. As the lensed image configuration is nearly circular, simple models give individual magnifications of at least 25 for the 3 brighter images. The source's narrow emission lines and low intrinsic luminosity qualify it as a type 2 AGN., Comment: 8 pages, 5 figures, submitted to AJ, comments welcome

Observations of caustic-crossing galaxies at redshift $0.72$ and the number density of events is greater around substructureand the number density of events is greater around substructures, and (ii) negative imaging regime where $\beta<2$. We study the particular case of seven microlensing events found by HST in the Dragon arc (at z=0.725). We find that a population of supergiant stars with a steep LF with $\beta=2.55$ fits the distribution of these events. We identify a small region of high density of microlensing events, and interpret it as evidence of a possible invisible substructure, for which we derive a mass of $\sim 1.3 \times 10^8\,\Msun$ (within its Einstein radius)., Comment: 23 pages, 16 figures

GHZ2/GLASS-z12 has been recently observed by JWST with both NIRSpec and MIRI spectrographs, making it the most distant galaxy ($z_{spec}=12.34$) with complete spectroscopic coverage from rest-frame UV to optical. It is identified as a strong CIV$_{1549}$ emitter with many detected emission lines (NIV], HeII, OIII], NIII], CIII], [OII], [NeIII], [OIII], and H$\alpha$), including a remarkable OIII$_{1333}$ Bowen fluorescence line. We analyze in this paper the joint NIRSpec+MIRI spectral data set. Combining six optical diagnostics (R2, R3, R23, O32, Ne3O2, and Ne3O2Hd), we find extreme ionization conditions, with O32 $=1.39 \pm 0.19$ and Ne3O2 $=0.37 \pm 0.18$ in stark excess compared to typical values in the ISM at lower redshifts. These line properties are compatible either with an AGN or with a compact, dense star-forming environment ($\Sigma_{\rm SFR}$ $\sim 10^2$-$10^3$ Msun/yr/kpc$^2$), with a high ionization parameter ($\log_{10}$(U) $=-1.75 \pm 0.16$), a high ionizing photon production efficiency $\log(\xi_{\rm ion}) = 25.7_{-0.1}^{+0.2}$, and a low, although not pristine, metal content between $5\%$ and $11\%$ Z$_\odot$ (confirmed by the T$_e$ method), indicating a rapid metal enrichment in the last few Myrs. These properties also suggest that a substantial amount of ionizing photons ($\sim 10\%$) are leaking outside. The general lessons learned from GHZ2 are the following: (i) the UV to optical combined nebular indicators are broadly in agreement with UV-only or optical-only indicators. (ii) UV+optical diagnostics fail to discriminate between an AGN and star-formation in a low metallicity, high density, and extreme ionization environment. (iii) comparing the nebular line ratios with local analogs may be approaching its limits at $z \gtrsim 10$, as this approach is potentially challenged by the unique conditions of star formation experienced by galaxies at these extreme redshifts., Comment: 14 pages, 3 figures, and 1 table. Resubmitted to ApJ. A section (3.4) has been added to the previous version

We present an analysis of the galaxy merger rate in the redshift range $4.0

We spectroscopically confirm the $M_{\rm UV} = -20.5$ mag galaxy GHZ2/GLASS-z12 to be at redshift $z=12.34$. The source was selected via NIRCam photometry in GLASS-JWST ERS data, providing the first evidence of a surprising abundance of bright galaxies at $z \gtrsim 10$. The NIRSpec PRISM spectrum shows detections of N IV, C IV, He II, O III, C III, O II, and Ne III lines, and the first detection at high-redshift of the O III Bowen fluorescence line at 3133 {\AA} rest-frame. The prominent C IV line with rest-frame equivalent width (EW) $\approx 46$ {\AA} puts GHZ2 in the category of extreme C IV emitters. GHZ2 displays UV lines with EWs that are only found in active galactic nuclei (AGNs) or composite objects at low/intermediate redshifts. The UV line-intensity ratios are compatible both with AGNs and star formation in a low-metallicity environment, with the low limit on the [Ne IV]/[N IV] ratio favoring a stellar origin of the ionizing photons. We discuss a possible scenario in which the high ionizing output is due to low metallicity stars forming in a dense environment. We estimate a metallicity $\lesssim 0.1 Z/{\rm Z}_{\odot}$, a high ionization parameter logU $> -2$, a N/O abundance 4--5 times the solar value, and a subsolar C/O ratio similar to the recently discovered class of nitrogen-enhanced objects. Considering its abundance patterns and the high stellar mass density ($10^4$~M$_{\odot}$~pc$^{-2}$), GHZ2 is an ideal formation site for the progenitors of today's globular clusters. The remarkable brightness of GHZ2 makes it a ``Rosetta stone'' for understanding the physics of galaxy formation within just 360 Myr after the Big Bang., Comment: Accepted for publication in The Astrophysical Journal. 20 pages, 7 figures, 2 tables

The internal structure and abundance of dark matter halos and subhalos are powerful probes of the nature of dark matter. In order to compare observations with dark matter models, accurate theoretical predictions of these quantities are needed. We present a fast and accurate method to describe the tidal evolution of subhalos within their parent halo, based on a semi-analytic approach. We first consider idealized N-body simulations of subhalos within their host halo, using a generalized mass density profile that describes their properties in a variety of dark matter models at infall, including popular warm, cold, and self-interacting ones. Using these simulations we construct tidal "tracks" for the evolution of subhalos based on their conditions at infall. Second, we use the results of these simulations to build semi-analytic models for tidal effects, including stripping and heating and implement them within the code GALACTICUS. Our semi-analytic models can accurately predict the tidal evolution of both cored and cuspy subhalos, including the bound mass and density profiles, providing a powerful and efficient tool for studying the post-infall properties of subhalos in different dark matter models., Comment: 25 pages, 50 figures. Version 2: matches published version. NOTE that the definition of dynamical timescale is slightly changed in version 2. The best-fit tidal stripping parameter is updated accordingly

Historically, metallicity profiles of galaxies have been modelled using a radially symmetric, two-parameter linear model, which reveals that most galaxies are more metal-rich in their central regions than their outskirts. However, this model is known to yield inaccurate results when the point-spread function (PSF) of a telescope is large. Furthermore, a radially symmetric model cannot capture asymmetric structures within a galaxy. In this work, we present an extension of the popular forward-modelling python package LENSTRONOMY, which allows the user to overcome both of these obstacles. We demonstrate the new features of this code base through two illustrative examples on simulated data. First, we show that through forward modelling, LENSTRONOMY is able to recover accurately the metallicity gradients of galaxies, even when the PSF is comparable to the size of a galaxy, as long as the data is observed with a sufficient number of pixels. Additionally, we demonstrate how LENSTRONOMY is able to fit irregular metallicity profiles to galaxies that are not well-described by a simple surface brightness profile. This opens up pathways for detailed investigations into the connections between morphology and chemical structure for galaxies at cosmological distances using the transformative capabilities of JWST. Our code is publicly available and open source, and can also be used to model spatial distributions of other galaxy properties that are traced by its surface brightness profile., Comment: 12 pages, 6 figures, including one appendix. Accepted for publication in Royal Astronomical Society Techniques & Instruments

The exceptional spectra of the most luminous $z>10$ sources observed so far have challenged our understanding of early galaxy evolution, requiring a new observational benchmark for meaningful interpretation. As such, we construct spectroscopic templates representative of high-redshift, star-forming populations, using 482 confirmed sources at $z=5.0-12.9$ with JWST/NIRSpec prism observations, and report on their average properties. We find $z=5-11$ galaxies are dominated by blue UV continuum slopes ($\beta=-2.3$ to $-2.7$) and reduced Balmer indices, characteristic of dust-poor and young systems, with a shift towards bluer slopes and younger ages with redshift. The evolution is mirrored by ubiquitous CIII] detections across all redshifts (rest-frame equivalent widths of $=5-14$ \r{A}), which increase in strength towards early times. Rest-frame optical lines reveal elevated ratios ($O32=7-31$, $R23=5-8$, and $Ne3O2=1-2$) and subsolar metallicities (log O/H$=7.3-7.9$), typical of ionization conditions and metallicities rarely observed in $z\sim0$ populations. Within our sample, we identify 57 Ly$\alpha$ emitters, which we stack and compare to a matched sample of nonemitters. The former are characterized by more extreme ionizing conditions with enhanced CIII], CIV, and HeII+[OIII] line emission, younger stellar populations from Balmer jumps, and a more pristine interstellar medium seen through bluer UV slopes and elevated rest-frame optical line ratios. The novel comparison illustrates important intrinsic differences between the two populations, with implications for Ly$\alpha$ visibility. The spectral templates derived here represent a new observational benchmark with which to interpret high-redshift sources, lifting our constraints on their global properties to unprecedented heights and extending out to the earliest of cosmic times., Comment: Accepted for publication in ApJ (with minor changes). Replaced to match accepted version. All composite spectra available online or upon request

We present a study on stellar properties of Lyman-alpha (Ly$\alpha$) emitters at 5 $< z <$ 8.2. We use 247 photometrically-selected, lensed, high-redshift, low luminosity galaxy candidates with spectroscopic follow-up. Of these, 38 are confirmed spectroscopically to be between 5 $< z <$ 8.2 via detection of Ly$\alpha$. For each galaxy and candidate, we estimate stellar mass, star formation rate, specific star formation rate, and mass-weighted age with spectral energy distribution fitting. We also measure the UV $\beta$ slope and luminosity using values from photometry. We find no strong correlation between Ly$\alpha$ equivalent width and any of these properties, as well as no significant difference between the physical properties of Ly$\alpha$ emitters and candidates without Ly$\alpha$ detected. This lack of expected trends may be explained by a combination of the evolving opacity of the IGM at these redshifts as well as the unique phase space probed by our lensed sample. Via tests on other galaxy samples which show varying strengths of correlations, we conclude that if there exist any relationships between Ly$\alpha$ EW and physical properties in the underlying population of faint galaxies, they are weak correlations. We also present the results of a spectroscopic search for CIII] emission in confirmed Ly$\alpha$ emitters at $z \sim 7$, finding no CIII] detections, but putting constraints on strong AGN activity and extreme nebular emission., Comment: v2, 14 pages, 8 figures, accepted to MNRAS after minor updates

We present deep JWST NIRSpec observations in the sightline of MACS J1149.5+2223, a massive cluster of galaxies at $z=0.54$. We report the spectroscopic redshift of 28 sources at $37$, i.e. comparable to the most metal poor galaxies in the local Universe. The search of primordial galaxies may be accomplished by extending toward a lower mass and/or by investigating inhomogeneities at smaller spatial scales. Lastly, we investigate potential systematics caused by the limitation of JWST's MSA observations. Caution is warranted when the target exceeds the slit size, as this situation could allow an overestimation of "global" metallicity, especially under the presence of strong negative metallicity gradient., Comment: Accepted for publication on May 29, 2024 to ApJ

We present a measurement of the mass-metallicity relation (MZR) at cosmic noon, using the JWST near-infrared wide-field slitless spectroscopy obtained by the GLASS-JWST Early Release Science program. By combining the power of JWST and the lensing magnification by the foreground cluster A2744, we extend the measurements of the MZR to the dwarf mass regime at high redshifts. A sample of 50 galaxies with several emission lines is identified across two wide redshift ranges of $z=1.8-2.3$ and $2.6-3.4$ in the stellar mass range of $\log{(M_*/M_\odot)}\in [6.9, 10.0]$. The observed slope of MZR is $0.223 \pm 0.017$ and $0.294 \pm 0.010$ at these two redshift ranges, respectively, consistent with the slopes measured in field galaxies with higher masses. In addition, we assess the impact of the morphological broadening on emission line measurement by comparing two methods of using 2D forward modeling and line profile fitting to 1D extracted spectra. We show that ignoring the morphological broadening effect when deriving line fluxes from grism spectra results in a systematic reduction of flux by $\sim30\%$ on average. This discrepancy appears to affect all the lines and thus does not lead to significant changes in flux ratio and metallicity measurements. This assessment of the morphological broadening effect using JWST data presents, for the first time, an important guideline for future work deriving galaxy line fluxes from wide-field slitless spectroscopy, such as Euclid, Roman, and the Chinese Space Station Telescope., Comment: 19 pages, 6 figures, accepted for publication in ApJL

The mass distribution in massive elliptical galaxies encodes their evolutionary history, thus providing an avenue to constrain the baryonic astrophysics in their evolution. The power-law assumption for the radial mass profile in ellipticals has been sufficient to describe several observables to the noise level, including strong lensing and stellar dynamics. In this paper, we quantitatively constrained any deviation, or the lack thereof, from the power-law mass profile in massive ellipticals through joint lensing-dynamics analysis of a large statistical sample with 77 galaxy-galaxy lens systems. We performed an improved and uniform lens modelling of these systems from archival Hubble Space Telescope imaging using the automated lens modelling pipeline dolphin. We combined the lens model posteriors with the stellar dynamics to constrain the deviation from the power law after accounting for the line-of-sight lensing effects, a first for analyses on galaxy-galaxy lenses. We find that the Sloan Lens ACS Survey (SLACS) lens galaxies with a mean redshift of 0.2 are consistent with the power-law profile within 1.1$\sigma$ (2.8$\sigma$) and the Strong Lensing Legacy Survey (SL2S) lens galaxies with a mean redshift of 0.6 are consistent within 0.8$\sigma$ (2.1$\sigma$), for a spatially constant (Osipkov-Merritt) stellar anisotropy profile. We adopted the spatially constant anisotropy profile as our baseline choice based on previous dynamical observables of local ellipticals. However, spatially resolved stellar kinematics of lens galaxies are necessary to differentiate between the two anisotropy models. Future studies will use our lens models to constrain the mass distribution individually in the dark matter and baryonic components., Comment: 31 pages, 20 figures, 9 tables. This version: accepted to MNRAS

The fraction of local dwarf galaxies that hosts massive black holes is arguably the cleanest diagnostic of the dominant seed formation mechanism of today's supermassive black holes. A 5 per cent constraint on this quantity can be achieved with AXIS observations of 3300 galaxies across the mass spectrum through a combination of serendipitous extra-galactic fields plus a dedicated 1 Msec GO program., Comment: This White Paper is part of a series commissioned for the AXIS Probe Concept Mission

The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ($z>\sim 10$). We study low-mass galaxy formation including this "stream velocity" using high resolution AREPO hydrodynamics simulations, and present theoretical predictions of the UV luminosity function (UVLF) and galaxy stellar mass function (GSMF) down to extremely faint and low mass galaxies ($M_{UV}>\sim-15$, $10^4M_\odot<=M_*<=10^8 M_\odot)$. We show that, although the stream velocity suppresses early star formation overall, it induces a short period of rapid star formation in some larger dwarfs, leading to an enhancement in the faint-end of the UVLF at $z=12$. We demonstrate that JWST observations are close to this enhanced regime, and propose that the UVLF may constitute an important probe of the stream velocity at high redshift for JWST and future observatories., Comment: 13 pages, 7 figures. Accepted ApJL

We describe the results of a new reverberation mapping program focused on the nearby Seyfert galaxy NGC 3227. Photometric and spectroscopic monitoring were carried out from 2022 December to 2023 June with the Las Cumbres Observatory network of telescopes. We detected time delays in several optical broad emission lines, with H$\beta$ having the longest delay at $\tau_{\rm cent}=4.0^{+0.9}_{-0.9}$ days and He II having the shortest delay with $\tau_{\rm cent}=0.9^{+1.1}_{-0.8}$ days. We also detect velocity-resolved behavior of the H$\beta$ emission line, with different line-of-sight velocities corresponding to different observed time delays. Combining the integrated H$\beta$ time delay with the width of the variable component of the emission line and a standard scale factor suggests a black hole mass of $M_{\rm BH}=1.1^{+0.2}_{-0.3} \times 10^7 M_{\odot}$. Modeling of the full velocity-resolved response of the H$\beta$ emission line with the phenomenological code CARAMEL finds a similar mass of $M_{\rm BH}=1.2^{+1.5}_{-0.7} \times 10^7 M_{\odot}$, and suggests that the H$\beta$-emitting broad line region (BLR) may be represented by a biconical or flared disk structure that we are viewing at an inclination angle of $\theta_i \approx 33^{\circ}$ and with gas motions that are dominated by rotation. The new photoionization-based BLR modeling tool BELMAC finds general agreement with the observations when assuming the best-fit CARAMEL results, however BELMAC prefers a thick disk geometry and kinematics that are equally comprised of rotation and inflow. Both codes infer a radially extended and flattened BLR that is not outflowing., Comment: 16 pages, 9 figures, 4 tables; submitted to ApJ

We analyze new JWST NIRCam and NIRSpec data on the redshift 9.11 galaxy MACS1149-JD1. Our NIRCam imaging data reveal that JD1 comprises three spatially distinct components. Our spectroscopic data indicate that JD1 appears dust-free but is already enriched, $12 + \log {\rm (O/H) } = 7.90^{+0.04}_{-0.05}$. We also find that the Carbon and Neon abundances in JD1 are below the solar abundance ratio. Particularly the Carbon under-abundance is suggestive of recent star formation where Type~II supernovae have already enriched the ISM in Oxygen but intermediate mass stars have not yet enriched the ISM in Carbon. A recent burst of star formation is also revealed by the star formation history derived from NIRCam photometry. Our data do not reveal the presence of a significant amount of old populations, resulting in a factor of $\sim7\times$ smaller stellar mass than previous estimates. Thus, our data support the view that JD1 is a young galaxy., Comment: Accepted for publication

We present a comprehensive search and analysis of high-redshift galaxies in a suite of nine public JWST extragalactic fields taken in Cycle 1, covering a total effective search area of $\sim358{\rm arcmin^2}$. Through conservative ($8\sigma$) photometric selection, we identify 341 galaxies at $51.5\,\sigma$ below the derived size-mass slope. These compact sources exhibit a high star formation surface density $\Sigma_{\rm SFR}>10\,M_\odot\,{\rm yr^{-1}\,kpc^{-2}}$, a range in which only $<0.01\,\%$ of the local star-forming galaxy sample is found. For those with available NIRSpec data, no evidence of ongoing supermassive black hole accretion is observed. A potential explanation for the observed high [OIII]-to-Hbeta ratios could be high shock velocities, likely originating within intense star-forming regions characterized by high $\Sigma_{\rm SFR}$. Lastly, we find that the rest-frame UV and optical sizes of our sample are comparable. Our results are consistent with these early galaxies building up their structures inside-out and yet to exhibit the strong color gradient seen at lower redshift., Comment: Accepted for publication in ApJ

An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 months observing campaign. One-dimensional echo-mapping models fit these variations if a slowly varying background is included for each emission line. These variations are more evident in the CIV light curve, which is the line least affected by intrinsic absorption in Mrk817 and least blended with neighboring emission lines. We identify five temporal windows with distinct emission line response, and measure their corresponding time delays, which range from 2 to 13 days. These temporal windows are plausibly linked to changes in the UV and X-ray obscuration occurring during these same intervals. The shortest time lags occur during periods with diminishing obscuration, whereas the longest lags occur during periods with rising obscuration. We propose that the obscuring outflow shields the ultraviolet broad lines from the ionizing continuum. The resulting change in the spectral energy distribution of the ionizing continuum, as seen by clouds at a range of distances from the nucleus, is responsible for the changes in the line response., Comment: 20 pages, 8 figures, Accepted for publication in ApJ

Time delays from strong gravitational lensing provide a one-step absolute distance measurement. Thus, they measure $H_0$ independently of all other probes. We first review the foundations and history of time-delay cosmography. Then, we illustrate the current state of the art by means of two recent case studies that have been real breakthroughs: i) the quadruply imaged quasar lensed by a galaxy-scale deflector RXJ1131$-$1231, for which spatially resolved stellar kinematics is available; ii) the multiply imaged supernova "Refsdal", the first with measured time delays, lensed by cluster MACS1149.5$+$2223. We conclude by discussing the exciting future prospects of time-delay cosmography in the coming decade., Comment: Invited chapter for the edited book "Hubble Constant Tension" (Eds. E. Di Valentino and D. Brout, Springer Singapore, expected in 2024)

Continuum reverberation mapping probes the sizescale of the optical continuum-emitting region in active galactic nuclei (AGN). Through 3 years of multiwavelength photometric monitoring in the optical with robotic observatories, we perform continuum reverberation mapping on Mrk~876. All wavebands show large amplitude variability and are well correlated. Slow variations in the light curves broaden the cross-correlation function (CCF) significantly, requiring detrending in order to robustly recover interband lags. We measure consistent interband lags using three techniques (CCF, JAVELIN, PyROA), with a lag of around 13~days from $u$ to $z$. These lags are longer than the expected radius of 12~days for the self-gravitating radius of the disk. The lags increase with wavelength roughly following $\lambda^{4/3}$, as would be expected from thin disk theory, but the lag normalization is approximately a factor of 3 longer than expected, as has also been observed in other AGN. The lag in the $i$ band shows an excess which we attribute to variable H$\alpha$ broad-line emission. A flux-flux analysis shows a variable spectrum that follows $f_\nu \propto \lambda^{-1/3}$ as expected for a disk, and an excess in the $i$ band that also points to strong variable H$\alpha$ emission in that band., Comment: 14 pages, 6 figures, 6 tables, accepted for publication in ApJ

The broad line region (BLR) size-luminosity relation has paramount importance for estimating the mass of black holes in active galactic nuclei (AGNs). Traditionally, the size of the H$\beta$ BLR is often estimated from the optical continuum luminosity at 5100\angstrom{} , while the size of the H$\alpha$ BLR and its correlation with the luminosity is much less constrained. As a part of the Seoul National University AGN Monitoring Project (SAMP) which provides six-year photometric and spectroscopic monitoring data, we present our measurements of the H$\alpha$ lags of 6 high-luminosity AGNs. Combined with the measurements for 42 AGNs from the literature, we derive the size-luminosity relations of H$\alpha$ BLR against broad H$\alpha$ and 5100\angstrom{} continuum luminosities. We find the slope of the relations to be $0.61\pm0.04$ and $0.59\pm0.04$, respectively, which are consistent with the \hb{} size-luminosity relation. Moreover, we find a linear relation between the 5100\angstrom{} continuum luminosity and the broad H$\alpha$ luminosity across 7 orders of magnitude. Using these results, we propose a new virial mass estimator based on the H$\alpha$ broad emission line, finding that the previous mass estimates based on the scaling relations in the literature are overestimated by up to 0.7 dex at masses lower than $10^7$~M$_{\odot}$., Comment: Accepted for publication in ApJ (Jun. 25th, 2023). 21 pages, 12 figures

The core-collapse supernova of a massive star rapidly brightens when a shock, produced following the collapse of its core, reaches the stellar surface. As the shock-heated star subsequently expands and cools, its early-time light curve should have a simple dependence on the progenitor's size and therefore final evolutionary state. Measurements of the progenitor's radius from early light curves exist for only a small sample of nearby supernovae, and almost all lack constraining ultraviolet observations within a day of explosion. The several-day time delays and magnifying ability of galaxy-scale gravitational lenses, however, should provide a powerful tool for measuring the early light curves of distant supernovae, and thereby studying massive stellar populations at high redshift. Here we analyse individual rest-frame ultraviolet-through-optical exposures taken with the Hubble Space Telescope that simultaneously capture, in three separate gravitationally lensed images, the early phases of a supernova at redshift $z \approx 3$ beginning within $5.8\pm 3.1$ hr of explosion. The supernova, seen at a lookback time of $\sim11.5$ billion years, is strongly lensed by an early-type galaxy in the Abell 370 cluster. We constrain the pre-explosion radius to be $533^{+154}_{-119}$ solar radii, consistent with a red supergiant. Highly confined and massive circumstellar material at the same radius can also reproduce the light curve, but is unlikely since no similar low-redshift examples are known., Comment: 69 pages, 12 figures/tables (4 main text, 8 extended data). Published in Nature

In late 2014, four images of Supernova (SN) "Refsdal," the first known example of a strongly lensed SN with multiple resolved images, were detected in the MACS J1149 galaxy-cluster field. Following the images' discovery, the SN was predicted to reappear within hundreds of days at a new position ~8 arcseconds away in the field. The observed reappearance in late 2015 makes it possible to carry out Refsdal's (1964) original proposal to use a multiply imaged SN to measure the Hubble constant H0, since the time delay between appearances should vary inversely with H0. Moreover, the position, brightness, and timing of the reappearance enable a novel test of the blind predictions of galaxy-cluster models, which are typically constrained only by the positions of multiply imaged galaxies. We have developed a new photometry pipeline that uses DOLPHOT to measure the fluxes of the five images of SN Refsdal from difference images. We apply four separate techniques to perform a blind measurement of the relative time delays and magnification ratios (mu_i/mu_1) between the last image SX and the earlier images S1-S4. We measure the relative time delay of SX-S1 to be 376.0+5.6-5.5 days and the relative magnification to be 0.30+0.05-0.03. This corresponds to a 1.5% precision on the time delay and 17% precision for the magnification ratios, and includes uncertainties due to millilensing and microlensing. In an accompanying paper, we place initial and blind constraints on the value of the Hubble constant., Comment: Published in ApJ. Companion paper presenting H0 constraints published in Science (DOI: 10.1126/science.abh1322)

The gravitationally lensed Supernova Refsdal appeared in multiple images, produced through gravitational lensing by a massive foreground galaxy cluster. After the supernova appeared in 2014, lens models of the galaxy cluster predicted an additional image of the supernova would appear in 2015, which was subsequently observed. We use the time delays between the images to perform a blinded measurement of the expansion rate of the Universe, quantified by the Hubble constant (H0). Using eight cluster lens models, we infer H0 = 64.8 +4.4-4.3 km / s / Mpc, where Mpc is the megaparsec. Using the two models most consistent with the observations, we find H0 = 66.6 +4.1-3.3 km / s / Mpc. The observations are best reproduced by models that assign dark-matter halos to individual galaxies and the overall cluster., Comment: Published in Science on May 11, 2023; this version updated to reflect minor edits to galley proofs. Companion paper presenting time-delay and relative magnification measurements published in ApJ (DOI: 10.3847/1538-4357/ac4ccb)

We combine our dynamical modeling black hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors, f, used in traditional reverberation mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have utilized our methods. Aiming to improve the precision of black hole mass estimates, as well as uncover any regularities in the behavior of the broad-line region (BLR), we search for correlations between f and other AGN/BLR parameters. We find (i) evidence for a correlation between the virial coefficient log10(fmean,{\sigma}) and black hole mass, (ii) marginal evidence for a similar correlation between log10(frms,{\sigma}) and black hole mass, (iii) marginal evidence for an anti-correlation of BLR disk thickness with log10(fmean,FWHM)and log10(frms,FWHM), and (iv) marginal evidence for an anti-correlation of inclination angle with log10(fmean,FWHM), log10(frms,{\sigma}), and log10(fmean,{\sigma}). Lastly, we find marginal evidence for a correlation between line-profile shape, when using the root-meansquare spectrum, log10(FWHM/{\sigma})rms, and the virial coefficient, log10(frms,{\sigma}), and investigate how BLR properties might be related to line-profile shape using caramel models., Comment: 17 pages, 10 figures

Dark matter subhalos with extended profiles and density cores, and globular stars clusters of mass $10^6-10^8 M_\odot$, that live near the critical curves in galaxy cluster lenses can potentially be detected through their lensing magnification of stars in background galaxies. In this work we study the effect such subhalos have on lensed images, and compare to the case of more well studied microlensing by stars and black holes near critical curves. We find that the cluster density gradient and the extended mass distribution of subhalos are important in determining image properties. Both lead to an asymmetry between the image properties on the positive and negative parity sides of the cluster that is more pronounced than in the case of microlensing. For example, on the negative parity side, subhalos with cores larger than about $50\,$pc do not generate any images with magnification above $\sim 100$ outside of the immediate vicinity of the cluster critical curve. We discuss these factors using analytical and numerical analysis, and exploit them to identify observable signatures of subhalos: subhalos create pixel-to-pixel flux variations of $\gtrsim 0.1$ magnitudes, on the positive parity side of clusters. These pixels tend to cluster around (otherwise invisible) subhalos. Unlike in the case of microlensing, signatures of subhalo lensing can be found up to $1''$ away from the critical curves of massive clusters., Comment: ApJ, accepted, 22 pages, 17 figures

Strong lensed quasi-stellar objects (QSOs) are valuable probes of the universe in numerous aspects. Two of these applications, reverberation mapping and measuring time delays for determining cosmological parameters, require the source QSOs to be variable with sufficient amplitude. In this paper, we forecast the number of strong lensed QSOs with sufficient variability to be detected by the Vera C. Rubin Telescope Legacy Survey of Space and Time (LSST). The damped random walk model is employed to model the variability amplitude of lensed QSOs taken from a mock catalog by Oguri & Marshall (2010). We expect 30--40% of the mock lensed QSO sample, which corresponds to $\sim$1000, to exhibit variability detectable with LSST. A smaller subsample of 250 lensed QSOs will show larger variability of $>0.15$~mag for bright lensed images with $i<21$ mag, allowing for monitoring with smaller telescopes. We discuss systematic uncertainties in the prediction by considering alternative prescriptions for variability and mock lens catalog with respect to our fiducial model. Our study shows that a large-scale survey of lensed QSOs can be conducted for reverberation mapping and time delay measurements following up on LSST., Comment: 11 pages, 5 figures, 1 table; submitted to MNRAS

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