Search

We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly-modeled sample to-date. Using the Prospector stellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift of $z=0.64^{+0.83}_{-0.32}$ ($68\%$ confidence), including 10 new or revised photometric redshifts at $z\gtrsim1$. We further find a median mass-weighted age of $t_m=0.8^{+2.71}_{-0.53}$Gyr, stellar mass of $\log(M_*/M_\odot)=9.69^{+0.75}_{-0.65}$, star formation rate of SFR=$1.44^{+9.37}_{-1.35}M_\odot$yr$^{-1}$, stellar metallicity of $\log(Z_*/Z_\odot)=-0.38^{+0.44}_{-0.42}$, and dust attenuation of $A_V=0.43^{+0.85}_{-0.36}$~mag (68\% confidence). Overall, the majority of short GRB hosts are star-forming ($\approx84\%$), with small fractions that are either transitioning ($\approx6\%$) or quiescent ($\approx10\%$); however, we observe a much larger fraction ($\approx40\%$) of quiescent and transitioning hosts at $z\lesssim0.25$, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal field galaxies, but do not include as many high-mass galaxies, implying that their binary neutron star (BNS) merger progenitors are dependent on a combination of host star formation and stellar mass. The distribution of ages and redshifts implies a broad delay-time distribution, with a fast-merging channel at $z>1$ and a decreased BNS formation efficiency at lower redshifts. If short GRB hosts are representative of BNS merger hosts within the horizon of current gravitational wave detectors, these results can inform future searches for electromagnetic counterparts. All of the data and modeling products are available on the BRIGHT website., 32 pages, 15 figures, 3 tables, published with ApJ

We present individual star-formation histories of $\sim3000$ massive galaxies (log($\mathrm{M_*/M_{\odot}}$) > 10.5) from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey at a lookback time of $\sim$7 billion years and quantify the population trends leveraging 20hr-deep integrated spectra of these $\sim$ 1800 star-forming and $\sim$ 1200 quiescent galaxies at 0.6 < $z$ < 1.0. Essentially all galaxies at this epoch contain stars of age < 3 Gyr, in contrast with older massive galaxies today, facilitating better recovery of previous generations of star formation at cosmic noon and earlier. We conduct spectro-photometric analysis using parametric and non-parametric Bayesian SPS modeling tools - Bagpipes and Prospector to constrain the median star-formation histories of this mass-complete sample and characterize population trends. A consistent picture arises for the late-time stellar mass growth when quantified as $t_{50}$ and $t_{90}$, corresponding to the age of the universe when galaxies formed 50\% and 90\% of their total stellar mass, although the two sets of models disagree at the earliest formation times (e.g. $t_{10}$). Our results reveal trends in both stellar mass and stellar velocity dispersion as in the local universe - low-mass galaxies with shallower potential wells grow their stellar masses later in cosmic history compared to high-mass galaxies. Unlike local quiescent galaxies, the median duration of late-time star-formation ($τ_{SF,late}$ = $t_{90}$ - $t_{50}$) does not consistently depend on the stellar mass. This census sets a benchmark for future deep spectro-photometric studies of the more distant universe., 25 pages, 12 figures, 4 tables, submitted to ApJ

We present a comprehensive catalog of observations and stellar population properties for 23 highly secure host galaxies of fast radio bursts (FRBs). Our sample comprises six repeating FRBs and 17 apparent non-repeaters. We present 82 new photometric and eight new spectroscopic observations of these hosts. Using stellar population synthesis modeling and employing non-parametric star formation histories (SFHs), we find that FRB hosts have a median stellar mass of $\approx 10^{9.9}\,M_{\odot}$, mass-weighted age $\approx 5.1$ Gyr, and ongoing star formation rate $\approx 1.3\,M_{\odot}$ yr$^{-1}$ but span wide ranges in all properties. Classifying the hosts by degree of star formation, we find that 87% (20/23 hosts) are star-forming, two are transitioning, and one is quiescent. The majority trace the star-forming main sequence of galaxies, but at least three FRBs in our sample originate in less active environments (two non-repeaters and one repeater). Across all modeled properties, we find no statistically significant distinction between the hosts of repeaters and non-repeaters. However, the hosts of repeating FRBs generally extend to lower stellar masses, and the hosts of non-repeaters arise in more optically luminous galaxies. While four of the galaxies with the most clear and prolonged rises in their SFHs all host repeating FRBs, demonstrating heightened star formation activity in the last $\lesssim 100$ Myr, one non-repeating host shows this SFH as well. Our results support progenitor models with short delay channels (i.e., magnetars formed via core-collapse supernova) for most FRBs, but the presence of some FRBs in less active environments suggests a fraction form through more delayed channels., Comment: 52 pages, 32 figures, 6 tables, submitted

The first JWST data revealed an unexpected population of red galaxies that appear to have redshifts of $z\sim 7-9$ and high masses of $M_*$ $\sim$ 10$^{10}$ M$_{\odot}$ (Labbé et al. 2023). Here we fit Sérsic profiles to the F200W NIRCam images of the 13 massive galaxy candidates of Labbé et al., to determine their structural parameters. Satisfactory fits were obtained for nine galaxies. We find that their effective radii are extremely small, ranging from $r_{\rm e}\sim 80$ pc to $r_{\rm e} \sim 300$ pc, with a mean of $\langle r_{\rm e}\rangle \approx 150$ pc. For their apparent stellar masses, the galaxies are smaller than any other galaxy population that has been observed at any other redshift. We use the fits to derive circularized three-dimensional stellar mass profiles of the galaxies, and compare these to the mass profiles of massive quiescent galaxies at $z\sim$2.3 and nearby elliptical galaxies. We find that, despite the high redshift galaxies having $10-20$ times smaller half-light radii, the central stellar densities are comparable to those of their putative descendants at later times. The most straightforward interpretation is that the dense compact inner regions of the most massive ellipticals today were already in place $\sim 600$ Myr after the Big Bang. We caution that the redshifts and masses of the galaxies remain to be confirmed, and that the complex NIRCam point spread function is not yet fully characterized., Submitted to ApJ Letters

We present high-resolution 1.5--6 GHz Karl G. Jansky Very Large Array (VLA) and $\textit{Hubble Space Telescope}$ ($\textit{HST}$) optical and infrared observations of the extremely active repeating fast radio burst (FRB) FRB$\,$20201124A and its barred spiral host galaxy. We constrain the location and morphology of star formation in the host and search for a persistent radio source (PRS) coincident with FRB$\,$20201124A. We resolve the morphology of the radio emission across all frequency bands and measure a star formation rate SFR $\approx 8.9\,M_{\odot}$ yr$^{-1}$, a factor of $\approx 4-6$ larger than optically-inferred SFRs, demonstrating dust-obscured star formation throughout the host. Compared to a sample of all known FRB hosts with radio emission, the host of FRB$\,$20201124A has the most significant obscured star formation. While ${\it HST}$ observations show the FRB to be offset from the bar or spiral arms, the radio emission extends to the FRB location. We propose that the FRB progenitor could have formed $\textit{in situ}$ (e.g., a magnetar central engine born from the explosion of a massive star). It is still plausible, although less likely, that the progenitor of FRB$\,$20201124A migrated from the central bar of the host, e.g., via a runaway massive star. We further place a limit on the luminosity of a putative PRS at the FRB position of $L_{\rm 6.0 \ GHz}$ $\lesssim$ 2.6 $\times$ $10^{27}$ erg s$^{-1}$ Hz$^{-1}$, two orders of magnitude below any PRS known to date. However, this limit is still broadly consistent with both magnetar nebulae and hypernebulae models assuming a constant energy injection rate of the magnetar and an age of $\gtrsim 10^{5}$ yr in each model, respectively., Comment: 21 pages, 6 figures, 3 tables, Submitted

We use CEERS JWST/NIRCam imaging to measure rest-frame near-IR light profiles of $>$500 $M_\star>10^{10}~M_\odot$ galaxies in the redshift range $0.510^{11}~M_\odot$) quiescent galaxies between $z=0.5$ and $z=2.3$, again comparing $R_{\rm{opt}}$, \rmass~and \rmassd. We conclude that the main tenets of the size evolution narrative established over the past 20 years, based on rest-frame optical light profile analysis, still hold in the era of JWST/NIRCam observations in the rest-frame near-IR., Submitted to ApJ. Comments welcome

The James Webb Space Telescope (JWST) is revolutionizing our knowledge of $z>5$ galaxies and their actively accreting black holes. Using the JWST Cycle 1 Treasury program Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) in the lensing field Abell 2744, we report the identification of a sample of little red dots at $3 < z_{\rm{phot}} < 7$ that likely contain highly-reddened accreting supermassive black holes. Using a NIRCam-only selection to F444W$, submitted to ApJ

We present the results of a search for high-redshift ($z>9$) galaxy candidates in the JWST UNCOVER survey, using deep NIRCam and NIRISS imaging in 7 bands over $\sim45$ arcmin$^2$ and ancillary HST observations. The NIRCam observations reach a $5-σ$ limiting magnitude of $\sim 29.2$ AB. The identification of high$-z$ candidates relies on a combination of a dropout selection and photometric redshifts. We find 16 candidates at $99$. A few galaxies at $z\sim9-10$ appear to show a clear Balmer break between the F356W and F444W/F410M bands, which helps constrain their stellar mass. We estimate blue UV continuum slopes between $β=-1.8$ and $-2.3$, typical for early galaxies at $z>9$ but not as extreme as the bluest recently discovered sources. We also find evidence for a rapid redshift-evolution of the mass-luminosity relation and a redshift-evolution of the UV continuum slope for a given range of intrinsic magnitude, in line with theoretical predictions. These findings suggest that deeper JWST observations are needed to reach the fainter galaxy population at those early epochs, and follow-up spectroscopy will help better constrain the physical properties and star formation histories of a larger sample of galaxies., Submitted to MNRAS

Flagship near-future surveys targeting $10^8-10^9$ galaxies across cosmic time will soon reveal the processes of galaxy assembly in unprecedented resolution. This creates an immediate computational challenge on effective analyses of the full data-set. With simulation-based inference (SBI), it is possible to attain complex posterior distributions with the accuracy of traditional methods but with a $>10^4$ increase in speed. However, it comes with a major limitation. Standard SBI requires the simulated data to have identical characteristics to the observed data, which is often violated in astronomical surveys due to inhomogeneous coverage and/or fluctuating sky and telescope conditions. In this work, we present a complete SBI-based methodology, ``SBI$^{++}$,'' for treating out-of-distribution measurement errors and missing data. We show that out-of-distribution errors can be approximated by using standard SBI evaluations and that missing data can be marginalized over using SBI evaluations over nearby data realizations in the training set. In addition to the validation set, we apply SBI$^{++}$ to galaxies identified in extragalactic images acquired by the James Webb Space Telescope, and show that SBI$^{++}$ can infer photometric redshifts at least as accurately as traditional sampling methods and crucially, better than the original SBI algorithm using training data with a wide range of observational errors. SBI$^{++}$ retains the fast inference speed of $\sim$1 sec for objects in the observational training set distribution, and additionally permits parameter inference outside of the trained noise and data at $\sim$1 min per object. This expanded regime has broad implications for future applications to astronomical surveys., Accepted for publication in ApJL; 13 pages, 5 figures. Code and a Jupyter tutorial are made publicly available at https://github.com/wangbingjie/sbi_pp

In November 2022, the James Webb Space Telescope (JWST) returned deep near-infrared images of Abell~2744 -- a powerful lensing cluster capable of magnifying distant, incipient galaxies beyond it. Together with the existing Hubble Space Telescope (HST) imaging, this publicly available dataset opens a fundamentally new discovery space to understand the remaining mysteries of the formation and evolution of galaxies across cosmic time. In this work, we detect and measure some 50,000 objects across the 45 arcmin$^2$ JWST footprint down to a $5\,\sigma$ limiting magnitude of $\sim$29.9\,mag in 0.32" apertures. Photometry is performed using circular apertures on images matched to the point spread function of the reddest NIRCam band, F444W, and cleaned of bright cluster galaxies and the related intra-cluster light. To give an impression of the photometric performance, we measure photometric redshifts and achieve a $\sigma_{\rm NMAD}\approx0.03$ based on known, but relatively small, spectroscopic samples. With this paper, we publicly release HST and JWST PSF-matched photometric catalogs optimized for bright and extended sources (0.7" apertures) and compact and faint sources (0.32" apertures) along with basic photometric redshifts, rest-frame colors, and individual magnification estimates. These catalogs will set the stage for efficient and deep spectroscopic follow-up of the first JWST-selected samples in Summer 2023., Comment: 27 pages, 19 figures, submitted to ApJS. Comments welcome!

Star formation in half of massive galaxies was quenched by the time the Universe was three billion years old. Very low amounts of molecular gas appear responsible for this, at least in some cases, though morphological gas stabilization, shock heating, or activity associated with accretion onto a central supermassive black hole is invoked in other cases. Recent studies of quenching by gas depletion have been based upon upper limits that are insufficiently sensitive to determine this robustly, or stacked emission with its problems of averaging. Here we report 1.3mm observations of dust emission from six strongly lensed galaxies where star formation has been quenched, with magnifications of up to a factor of 30. Four of the six galaxies are undetected in dust emission, with an estimated upper limit on the dust mass of 0.0001 times the stellar mass, and by proxy (assuming a Milky Way molecular gas-to-dust ratio) 0.01 times the stellar mass in molecular gas. This is two orders of magnitude less molecular gas per unit stellar mass than seen in star forming galaxies at similar redshifts. It remains difficult to extrapolate from these small samples, but these observations establish that gas depletion is responsible for a cessation of star formation in some fraction of high-redshift galaxies., 17 pages, 3 figures. Authors' version. Published online by Nature on September 22, 2021

We explore how the presence of detectable molecular gas depends on the inferred star formation histories (SFHs) in eight massive, quiescent galaxies at z ∼ 0.7. Half of the sample have clear detections of molecular gas, traced by CO(2–1). We find that the molecular gas content is unrelated to the rate of star formation decline prior to the most recent 1 Gyr, suggesting that the gas reservoirs are not left over from their primary star formation epoch. However, the recent SFHs of CO-detected galaxies demonstrate evidence for secondary bursts of star formation in their last Gyr. The fraction of stellar mass formed in these secondary bursts ranges from f burst ≈ 0.3%–6% and ended between t end-burst ≈ 0–330 Myr ago. The CO-detected galaxies form a higher fraction of mass in the last Gyr ( f M 1 Gyr = 2.6 % ± 1.8 % ) compared to the CO-undetected galaxies ( f M 1 Gyr = 0.2 % ± 0.1 % ). The galaxies with gas reservoirs have enhanced late-time star formation, highlighting this as a contributing factor to the observed heterogeneity in the gas reservoirs in high-redshift quiescent galaxies. We find that the amount of gas and star formation driven by these secondary bursts are inconsistent with that expected from dry minor mergers, and instead are likely driven by recently accreted gas, i.e., gas-rich minor mergers. This conclusion would not have been made based on SFRUV+IR measurements alone, highlighting the power of detailed SFH modeling in the interpretation of gas reservoirs. Larger samples are needed to understand the frequency of low-level rejuvenation among quiescent galaxies at intermediate redshifts, and to what extent this drives the diversity of molecular gas reservoirs.

The first few 100 Myr at z > 10 mark the last major uncharted epoch in the history of the universe, where only a single galaxy (GN-z11 at z ≈ 11) is currently spectroscopically confirmed. Here we present a search for luminous z > 10 galaxies with JWST/NIRCam photometry spanning ≈1–5 μm and covering 49 arcmin2 from the public JWST Early Release Science programs (CEERS and GLASS). Our most secure candidates are two MUV ≈ −21 systems: GLASS-z12 and GLASS-z10. These galaxies display abrupt ≳1.8 mag breaks in their spectral energy distributions (SEDs), consistent with complete absorption of flux bluewards of Lyα that is redshifted to =+z12.40.3 0.1and=+z10.40.5 0.4. Lower redshift interlopers such as quiescent galaxies with strong Balmer breaks would be comfortably detected at >5σ in multiple bands where instead we find no flux. From SED modeling we infer that these galaxies have already built up ∼109 solar masses in stars over the ≲300–400 Myr after the Big Bang. The brightness of these sources enable morphological constraints. Tantalizingly, GLASS-z10 shows a clearly extended exponential light profile, potentially consistent with a disk galaxy of r50 ≈ 0.7 kpc. These sources, if confirmed, join GN-z11 in defying number density forecasts for luminous galaxies based on Schechter UV luminosity functions, which require a survey area >10× larger than we have studied here to find such luminous sources at such high redshifts. They extend evidence from lower redshifts for little or no evolution in the bright end of the UV luminosity function into the cosmic dawn epoch, with implications for just how early these galaxies began forming. This, in turn, suggests that future deep JWST observations may identify relatively bright galaxies to much earlier epochs than might have been anticipated., R.P.N. acknowledges funding from JWST programs GO-1933 and GO-2279. Support for this work was provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51515.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. We acknowledge support from: the Swiss National Science Foundation through project grant 200020_207349 (P.A.O., A.W.). The Cosmic Dawn Center (DAWN) is funded by the Danish National Research Foundation under grant No. 140. R.J.B. and M.S. acknowledge support from NWO grant TOP1.16.057

With just a month of data, JWST is already transforming our view of the universe, revealing and resolving starlight in unprecedented populations of galaxies. Although “HST-dark” galaxies have previously been detected at long wavelengths, these observations generally suffer from a lack of spatial resolution, which limits our ability to characterize their sizes and morphologies. Here we report on a first view of starlight from a subset of the HST-dark population that is bright with JWST/NIRCam (4.4 μm < 24.5 mag) and very faint or even invisible with HST (μm). In this Letter we focus on a dramatic and unanticipated population of physically extended galaxies (≳0.″25). These 12 galaxies have photometric redshifts 2 < z < 6, high stellar masses M ⋆ ≳ 1010 M ⊙, and significant dust-attenuated star formation. Surprisingly, the galaxies have elongated projected axis ratios at 4.4 μm, suggesting that the population is disk dominated or prolate and we hence refer to them as ultrared flattened objects. Most of the galaxies appear red at all radii, suggesting significant dust attenuation throughout. With R e (F444W) ∼ 1–2 kpc, the galaxies are similar in size to compact massive galaxies at z ∼ 2 and the cores of massive galaxies and S0s at z ∼ 0. The stellar masses, sizes, and morphologies of the sample suggest that some could be progenitors of lenticular or fast-rotating galaxies in the local universe. The existence of this population suggests that our previous censuses of the universe may have missed massive, dusty edge-on disks, in addition to dust-obscured starbursts.

We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $��_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($��_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langle��_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and "late bloomers" are uncommon. In summary, our analysis supports the "grow & quench" framework and is consistent with a wide and continuously-populated diversity of quenching timescales., 30 pages, 12 figures + appendix, accepted by ApJ

The variety of star formation histories (SFHs) of $z\gtrsim6$ galaxies provides important insights into early star formation, but has been difficult to systematically quantify. Some observations suggest that many $z\sim6-9$ galaxies are dominated by $\gtrsim200$ Myr stellar populations, implying significant star formation at $z\gtrsim9$, while others find that most reionization era galaxies are $\lesssim10$ Myr, consistent with little $z\gtrsim9$ star formation. Here, we quantify the distribution of ages of UV-bright ($-22.5\lesssim M_{UV}\lesssim-21$) galaxies colour-selected to lie at $z\simeq6.6-6.9$, an ideal redshift range to systematically study the SFHs of reionization era galaxies with ground-based observatories and Spitzer. We infer galaxy properties with two SED modelling codes and compare results, finding that stellar masses are largely insensitive to the model, but the inferred ages can vary by an order of magnitude. We infer a distribution of ages assuming a simple, parametric SFH model, finding a median age of $\sim30-70$ Myr depending on SED model. We quantify the fractions of $\leq10$ Myr and $\geq250$ Myr galaxies, finding that these systems comprise $\sim15-30$ percent and $\sim20-25$ percent of the population, respectively. With a flexible SFH model, the shapes of the SFHs are consistent with those implied by the simple model (e.g. young galaxies have rapidly rising SFHs). However, stellar masses can differ significantly, with those of young systems sometimes being more than an order of magnitude larger with the flexible SFH. We quantify the implications of these results for $z\gtrsim9$ stellar mass assembly and discuss improvements expected from JWST., Comment: 24 pages, 20 figures, 2 appendices, accepted to MNRAS

Dust plays a pivotal role in determining the observed spectral energy distribution (SED) of galaxies. Yet our understanding of dust attenuation is limited and our observations suffer from the dust-metallicity-age degeneracy in SED fitting (single galaxies), large individual variances (ensemble measurements), and the difficulty in properly dealing with uncertainties (statistical considerations). In this study, we create a population Bayesian model to rigorously account for correlated variables and non-Gaussian error distributions and demonstrate the improvement over a simple Bayesian model. We employ a flexible 5-D linear interpolation model for the parameters that control dust attenuation curves as a function of stellar mass, star formation rate (SFR), metallicity, redshift, and inclination. Our setup allows us to determine the complex relationships between dust attenuation and these galaxy properties simultaneously. Using Prospector fits of nearly 30,000 3D-HST galaxies, we find that the attenuation slope ($n$) flattens with increasing optical depth ($\tau$), though less so than in previous studies. $\tau$ increases strongly with SFR, though when $\log~{\rm SFR}\lesssim 0$, $\tau$ remains roughly constant over a wide range of stellar masses. Edge-on galaxies tend to have larger $\tau$ than face-on galaxies, but only for $\log~M_*\gtrsim 10$, reflecting the lack of triaxiality for low-mass galaxies. Redshift evolution of dust attenuation is strongest for low-mass, low-SFR galaxies, with higher optical depths but flatter curves at high redshift. Finally, $n$ has a complex relationship with stellar mass, highlighting the intricacies of the star-dust geometry. We have publicly released software (https://github.com/Astropianist/DustE) for users to access our population model., Comment: 27 pages, 14 figures, 1 table. Submitted to ApJ

With just a month of data, JWST is already transforming our view of the Universe, revealing and resolving starlight in unprecedented populations of galaxies. Although ``HST-dark" galaxies have previously been detected at long wavelengths, these observations generally suffer from a lack of spatial resolution which limits our ability to characterize their sizes and morphologies. Here we report on a first view of starlight from a subset of the HST-dark population that are bright with JWST/NIRCam (4.4$μ$m

We present BVRI and unfiltered Clear light curves of 70 stripped-envelope supernovae (SESNe), observed between 2003 and 2020, from the Lick Observatory Supernova Search (LOSS) follow-up program. Our SESN sample consists of 19 spectroscopically normal SNe~Ib, two peculiar SNe Ib, six SN Ibn, 14 normal SNe Ic, one peculiar SN Ic, ten SNe Ic-BL, 15 SNe IIb, one ambiguous SN IIb/Ib/c, and two superluminous SNe. Our follow-up photometry has (on a per-SN basis) a mean coverage of 81 photometric points (median of 58 points) and a mean cadence of 3.6d (median of 1.2d). From our full sample, a subset of 38 SNe have pre-maximum coverage in at least one passband, allowing for the peak brightness of each SN in this subset to be quantitatively determined. We describe our data collection and processing techniques, with emphasis toward our automated photometry pipeline, from which we derive publicly available data products to enable and encourage further study by the community. Using these data products, we derive host-galaxy extinction values through the empirical colour evolution relationship and, for the first time, produce accurate rise-time measurements for a large sample of SESNe in both optical and infrared passbands. By modeling multiband light curves, we find that SNe Ic tend to have lower ejecta masses and lower ejecta velocities than SNe~Ib and IIb, but higher $^{56}$Ni masses., Comment: Accepted by MNRAS

$JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5\mu$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $\Lambda$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $2\mu$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($\Delta z\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical., Comment: Submitted to ApJL. Fig. 1 shows the "superposed" state, Fig. 2 presents the potential z~5 protocluster, Fig. 3 shows some atypical z~5 SEDs, Fig. 4 studies the morphology, Fig. 5 summarizes the extraordinary situation at z~17, and Fig. 6 argues for why such interlopers may not be a major issue for z>10 searches. Comments warmly welcomed

W-CDF-S, ELAIS-S1, and XMM-LSS will be three Deep-Drilling Fields (DDFs) of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), but their extensive multi-wavelength data have not been fully utilized as done in the COSMOS field, another LSST DDF. To prepare for future science, we fit source spectral energy distributions (SEDs) from X-ray to far-infrared in these three fields mainly to derive galaxy stellar masses and star-formation rates. We use CIGALE v2022.0, a code that has been regularly developed and evaluated, for the SED fitting. Our catalog includes 0.8 million sources covering $4.9~\mathrm{deg^2}$ in W-CDF-S, 0.8 million sources covering $3.4~\mathrm{deg^2}$ in ELAIS-S1, and 1.2 million sources covering $4.9~\mathrm{deg^2}$ in XMM-LSS. Besides fitting normal galaxies, we also select candidates that may host active galactic nuclei (AGNs) or are experiencing recent star-formation variations and use models specifically designed for these sources to fit their SEDs; this increases the utility of our catalog for various projects in the future. We calibrate our measurements by comparison with those in well-studied smaller regions and briefly discuss the implications of our results. We also perform detailed tests of the completeness and purity of SED-selected AGNs. Our data can be retrieved from a public website., Comment: 49 pages, 39 figures, 9 tables, accepted for publication in ApJS

Rejuvenating galaxies are unusual galaxies that fully quench and then subsequently experience a "rejuvenation" event to become star-forming once more. Rejuvenation rates vary substantially in models of galaxy formation: 10%-70% of massive galaxies are expected to experience rejuvenation by z = 0. Measuring the rate of rejuvenation is therefore important for calibrating the strength of star formation feedback mechanisms. However, these observations are challenging because rejuvenating systems blend in with normal star-forming galaxies in broadband photometry. In this paper, we use the galaxy spectral energy distribution (SED)-fitting code Prospector to search for observational markers that distinguish normal star-forming galaxies from rejuvenating galaxies. We find that rejuvenating galaxies have smaller Balmer absorption line equivalent widths (EWs) than star-forming galaxies. This is analogous to the well-known "K + A" or post-starburst galaxies, which have strong Balmer absorption due to A-stars dominating the light: in this case, rejuvenating systems have a lack of A-stars, instead resembling "O - A" systems. We find star-forming galaxies that have H$\beta$, H$\gamma$, and/or H$\delta$ absorption EWs $\lesssim 3${\AA} corresponds to a highly pure selection of rejuvenating systems. Interestingly, while this technique is highly effective at identifying mild rejuvenation, "strongly" rejuvenating systems remain nearly indistinguishable from star-forming galaxies due to the well-known stellar outshining effect. We conclude that measuring Balmer absorption line EWs in star-forming galaxy populations is an efficient method to identify rejuvenating populations, and discuss several techniques to either remove or resolve the nebular emission which typically lies on top of these absorption lines., Comment: Published in ApJ, 16 pages, 9 figures

We present a new rest-frame color-color selection method using "synthetic $u_s-g_s$ and $g_s-i_s$'', $(ugi)_s$ colors to identify star-forming and quiescent galaxies. Our method is similar to the widely-used $U-V$ versus $V-J$ ($UVJ$) diagram. However, $UVJ$ suffers known systematics. Spectroscopic campaigns have shown that $UVJ$-selected quiescent samples at $z \gtrsim 3$ include $\sim 10-30\%$ contamination from galaxies with dust-obscured star formation and strong emission lines. Moreover, at $z>3$, $UVJ$ colors are extrapolated because the rest-frame J-band shifts beyond the coverage of the deepest bandpasses at $< 5~\mu m$ (typically $Spitzer$/IRAC 4.5 $\mu m$ or future $JWST$/NIRCam observations). We demonstrate that $(ugi)_s$ offers improvements to $UVJ$ at $z>3$, and can be applied to galaxies in the $JWST$ era. We apply $(ugi)_s$ selection to galaxies at $0.5, Comment: 33 pages, 20 figures, accepted for publication in ApJ

REQUIEM-2D (REsolving QUIEscent Magnified galaxies with 2D grism spectroscopy) is comprised of a sample of 8 massive ($\log M_*/M_\odot > 10.6$) strongly lensed quiescent galaxies at $z\sim2$. REQUIEM-2D combines the natural magnification from strong gravitational lensing with the high spatial-resolution grism spectroscopy of \emph{Hubble Space Telescope} through a spectrophotometric fit to study spatially resolved stellar populations. We show that quiescent galaxies in the REQUIEM-2D survey have diverse formation histories manifesting as a gradient in stellar ages, including examples of (1) a younger central region supporting outside-in formation, (2) flat age gradients that show evidence for both spatially-uniform early formation or inside-out quenching, and (3) regions at a fixed radial distance having different ages (such asymmetries cannot be recovered when averaging stellar population measurements azimuthally). The typical dust attenuation curve for the REQUIEM-2D galaxies is constrained to be steeper than Calzetti's law in the UV and generally consistent with $A_V, Submitted to ApJ

The advent of the James Webb Space Telescope (JWST) signals a new era in exploring galaxies in the high-$z$ universe. Current and upcoming JWST imaging will potentially detect galaxies out to $z \sim 20$, creating a new urgency in the quest to infer accurate photometric redshifts (photo-$z$) for individual galaxies from their spectral energy distributions, as well as masses, ages and star formation rates. Here we illustrate the utility of informed priors encoding previous observations of galaxies across cosmic time in achieving these goals. We construct three joint priors encoding empirical constraints of redshifts, masses, and star formation histories in the galaxy population within the \prospector\ Bayesian inference framework. In contrast with uniform priors, our model breaks an age-mass-redshift degeneracy, and thus reduces the mean bias error in masses from 0.3 to 0.1 dex, and in ages from 0.6 to 0.2 dex in tests done on mock JWST observations. Notably, our model recovers redshifts at least as accurately as the state-of-the-art photo-$z$ code \eazy\ in deep JWST fields, but with two advantages: tailoring a model based on a particular survey renders mostly unnecessary given well-motivated priors; obtaining joint posteriors describing stellar, active galactic nuclei, gas, and dust contributions becomes possible. We can now confidently use the joint distribution to propagate full non-Gaussian redshift uncertainties into inferred properties of the galaxy population. This model, ``\prospector-$\beta$'', is intended for fitting galaxy photometry where the redshift is unknown, and will be instrumental in ensuring the maximum science return from forthcoming photometric surveys with JWST. The code is made publicly available online as a part of \prospector., Comment: Accepted for publication in ApJL. 13 pages, 6 figures, 2 tables. The code is made publicly available online as a part of Prospector; the version used in this work corresponds to the state of the Git repository at commit 820ad72

The remarkably tight relationship between galaxy optical color and stellar mass-to-light ratio (M */L) is widely used for efficient stellar mass estimates. However, it remains unclear whether this low scatter comes from a natural order in the galaxy population, or whether it is driven by simple relationships in the models used to describe them. In this work, we investigate the origins of the relationship by contrasting the derived relationship from a simple 4 parameter spectral energy distribution (SED) model with a more sophisticated 14D Prospector–α model including nonparametric star formation histories (SFHs). We apply these models to 63,430 galaxies at 0.5 < z < 3 and fit a hierarchical Bayesian model (HBM) to the population distribution in the (g − r)– log ( M / L g ) plane. We find that Prospector–α infers systematically higher M */L by 0.12 dex, a result of nonparametric SFHs producing older ages, and also systematically redder rest frame (g − r) by 0.06 mag owing to the contribution from nebular emission. Surprisingly, the combined effects of the M */L and (g − r) offsets produce a similar average relationship for the two models; though Prospector–α produces a higher scatter of 0.28 dex compared to the simple model of 0.12 dex. Also, unlike the simple model, the Prospector–α relationship shows substantial redshift evolution due to stellar aging. These expected and testable effects produce overall older and redder galaxies; though the color–M */L relationship is measured only at 0.5 < z < 3. Finally, we demonstrate that the HBM produces substantial shrinkage in the individual posteriors of faint galaxies, an important first step toward using the observed galaxy population directly to inform the SED-fitting priors.

Modelling reionization often requires significant assumptions about the properties of ionizing sources. Here, we infer the total output of hydrogen-ionizing photons (the ionizing emissivity, $\dot{N}_\textrm{ion}$) at $z=4-14$ from current reionization constraints, being maximally agnostic to the properties of ionizing sources. We use a Bayesian analysis to fit for a non-parametric form of $\dot{N}_\textrm{ion}$, allowing us to flexibly explore the entire prior volume. We infer a declining $\dot{N}_\textrm{ion}$ with redshift at $z>6$, which can be used as a benchmark for reionization models. Model-independent reionization constraints from the CMB optical depth and Ly$\alpha$ and Ly$\beta$ forest dark pixel fraction produce $\dot{N}_\textrm{ion}$ evolution ($d\log_{10}\dot{N}_\textrm{ion}/dz|_{z=6\rightarrow8} = -0.31\pm0.35$ dex) consistent with the declining UV luminosity density of galaxies, assuming constant ionizing photon escape fraction and efficiency. Including measurements from Ly$\alpha$ damping of galaxies and quasars produces a more rapid decline: $d\log_{10}\dot{N}_\textrm{ion}/dz|_{z=6\rightarrow8} =-0.44\pm0.22$ dex, steeper than the declining galaxy luminosity density (if extrapolated beyond $M_\mathrm{UV} \lesssim -13$), and constrains the mid-point of reionization to $z = 6.93\pm0.14$., Comment: submitted to MNRAS

Galaxies are complicated physical systems which obey complex scaling relationships; as a result, properties measured from broadband photometry are often highly correlated, degenerate, or both. Therefore, the accuracy of basic properties like stellar masses and star formation rates (SFRs) depend on the accuracy of many second-order galaxy properties, including star formation histories (SFHs), stellar metallicities, dust properties, and many others. Here, we re-assess measurements of galaxy stellar masses and SFRs using a 14-parameter physical model built in the Prospector Bayesian inference framework. We find that galaxies are ∼0.2 dex more massive and have ∼0.2 dex lower star formation rates than classic measurements. These measurements lower the observed cosmic star formation rate density and increase the observed buildup of stellar mass, finally bringing these two metrics into agreement at the factor-of-two level at 0.5 < z < 2.5.

Accurate models of the star formation histories (SFHs) of recently-quenched galaxies can provide constraints on when and how galaxies shut down their star formation. The recent development of "non-parametric" SFH models promises the flexibility required to make these measurements. However, model and prior choices significantly affect derived SFHs, particularly for post-starburst galaxies (PSBs) which have sharp changes in their recent SFH. In this paper, we create mock PSBs, then use the Prospector SED fitting software to test how well four different SFH models recover key properties. We find that a two-component parametric model performs well for our simple mock galaxies, but is sensitive to model mismatches. The fixed- and flexible-bin non-parametric models included in Prospector are able to rapidly quench a major burst of star formation, but systematically underestimate the post-burst age by up to 200 Myr. We develop a custom SFH model that allows for additional flexibility in the recent SFH. Our flexible non-parametric model is able to constrain post-burst ages with no significant offset and just ~90 Myr of scatter. Our results suggest that while standard non-parametric models are able to recover first-order quantities of the SFH (mass, SFR, average age), accurately recovering higher-order quantities (burst fraction, quenching time) requires careful consideration of model flexibility. These mock recovery tests are a critical part of future SFH studies. Finally, we show that our new, public SFH model is able to accurately recover the properties of mock star-forming and quiescent galaxies and is suitable for broader use in the SED fitting community., 24 pages, 10 figures. Posted after referee report; comments welcome. Key results in Figures 8 & 9. New SFH model publicly available as "continuity_psb_sfh" in Prospector: https://github.com/bd-j/prospector

We use the panchromatic SED-fitting code Prospector to measure the galaxy logM$^*$-logSFR relationship (the `star-forming sequence') across $0.2 < z < 3.0$ using the COSMOS-2015 and 3D-HST UV-IR photometric catalogs. We demonstrate that the chosen method of identifying star-forming galaxies introduces a systematic uncertainty in the inferred normalization and width of the star-forming sequence, peaking for massive galaxies at $\sim 0.5$ dex and $\sim0.2$ dex respectively. To avoid this systematic, we instead parameterize the density of the full galaxy population in the logM$^*$-logSFR-redshift plane using a flexible neural network known as a normalizing flow. The resulting star-forming sequence has a low-mass slope near unity and a much flatter slope at higher masses, with a normalization $0.2-0.5$ dex lower than typical inferences in the literature. We show this difference is due to the sophistication of the Prospector stellar populations modeling: the nonparametric star formation histories naturally produce higher masses while the combination of individualized metallicity, dust, and star formation history constraints produce lower star formation rates than typical UV+IR formulae. We introduce a simple formalism to understand the difference between SFRs inferred from spectral energy distribution fitting and standard template-based approaches such as UV+IR SFRs. Finally, we demonstrate the inferred star-forming sequence is consistent with predictions from theoretical models of galaxy formation, resolving a long-standing $\sim0.2-0.5$ dex offset with observations at $0.5, Accepted to ApJ. Trained flow available at https://github.com/jrleja/sfs_leja_trained_flow

We present the complete set of {\it Hubble Space Telescope} imaging of the binary neutron star merger GW170817 and its optical counterpart AT 2017gfo. Including deep template imaging in F814W, F110W, F140W, and F160W at 3.4 years post-merger, we re-analyze the full light curve of AT 2017gfo across 12 bands from 5--1273 rest-frame days after merger. We obtain four new detections of the short $\gamma$-ray burst (GRB) 170817A afterglow from 109--170 rest-frame days post-merger. These detections are consistent with the previously observed $\beta=-0.6$ spectral index in the afterglow light curve with no evidence for spectral evolution. We also analyze our limits in the case of novel late-time optical and IR emission signatures, such as a kilonova afterglow or infrared dust echo, but find our limits are not constraining in these contexts. We use the new data to construct deep optical and infrared stacks, reaching limits of $M=-6.3$ to $-4.6$ mag, to analyze the local environment around AT 2017gfo and low surface brightness features in its host galaxy NGC 4993. We rule out the presence of any globular cluster at the position of AT 2017gfo to $2.3 \times 10^{4} L_{\odot}$, including those with the reddest $V-H$ colors. Finally, we analyze the substructure of NGC 4993 in deep residual imaging, and find shell features which extend up to 71.8\arcsec\ (14.2 kpc) from the center of the galaxy. We find that the shells have a cumulative stellar mass of $6.3\times10^{8} M_{\odot}$, roughly 2% the total stellar mass of NGC 4993, and mass-weighted ages of $>$3 Gyr. We conclude that it was unlikely the GW170817 progenitor system formed in the galaxy merger, which based on dynamical signatures and the stellar population in the shells mostly likely occurred 220--685 Myr ago., Comment: 21 pages, 9 figures, submitted to ApJ

Using deep rest-frame optical spectroscopy from the Large Early Galaxy Astrophysical Census (LEGA-C) survey, conducted using VIMOS on the ESO Very Large Telescope, we systematically search for low-ionization [OII] 3726,3729 emission in the spectra of a mass-complete sample of z~0.85 galaxies. Intriguingly, we find that 59 percent of UVJ-quiescent (i.e. non star-forming) galaxies in the sample have ionized gas, as traced by [OII] emission, detected above our completeness limit of 1.5 Angstroms. The median stacked spectrum of the lowest equivalent width quiescent galaxies also shows [OII] emission. The overall fraction of sources with [OII] above our equivalent width limit is comparable to what we find in the low-redshift Universe from GAMA and MASSIVE, except perhaps at the highest stellar masses (log Mstar/Msol > 11.5). However, stacked spectra for the individual low-equivalent width systems uniquely indicates ubiquitous [OII] emission in the higher-z LEGA-C sample, with typical [OII] luminosities per unit stellar mass that are a factor of 3 larger than the lower-z GAMA sample. Star formation in these otherwise quiescent galaxies could play a role in producing the [OII] emission at higher-z, although it is unlikely to provide the bulk of the ionizing photons. More work is required to fully quantify the contributions of evolved stellar populations or active galactic nuclei to the observed spectra., ApJ accepted: 16 pages, 11 figures

Calcium-rich (Ca-rich) transients are a new class of supernovae (SNe) that are known for their comparatively rapid evolution, modest peak luminosities, and strong nebular calcium emission lines. Currently, the progenitor systems of Ca-rich transients remain unknown. Although they exhibit spectroscopic properties not unlike core-collapse Type Ib/c SNe, nearly half are found in the outskirts of their host galaxies that are predominantly elliptical, suggesting a closer connection to the older stellar populations of SNe Ia. In this paper, we present a compilation of publicly available multiwavelength observations of all known and/or suspected host galaxies of Ca-rich transients ranging from FUV to IR, and use these data to characterize their stellar populations with prospector. We estimate several galaxy parameters including integrated star formation rate, stellar mass, metallicity, and age. For nine host galaxies, the observations are sensitive enough to obtain nonparametric star formation histories, from which we recover SN rates and estimate probabilities that the Ca-rich transients in each of these host galaxies originated from a core-collapse vs. Type Ia-like explosion. Our work supports the notion that the population of Ca-rich transients do not come exclusively from core-collapse explosions, and must either be only from white dwarf stars or a mixed population of white dwarf stars with other channels, potentially including massive star explosions. Additional photometry and explosion site spectroscopy of larger samples of Ca-rich host galaxies will improve these estimates and better constrain the ratio of white dwarf vs. massive star progenitors of Ca-rich transients., Comment: 22 pages, 7 figures, 5 tables, accepted to ApJ

When the Universe was merely three billion years old, about half of massive galaxies had already formed the bulk of their stars and new star formation plummeted [1]. How galaxies quench at such early times remains a puzzle, as their dark matter halos contain large gas reservoirs [2-4]. This gas should cool efficiently, sustaining star formation over long periods [5,6]. Here we present sensitive 1.3mm wavelength observations of cold dust in six quenched galaxies in the redshift range z=1.6 to z=3.2 with stellar masses ranging from 2.5x1010M⊙ to 5x1011M⊙, which are magnified by foreground galaxy clusters. Even with factors of up to 30 in magnification, four of the six galaxies are undetected at this wavelength. We show that these quenched galaxies have extremely little dust at early times, and by proxy very little cold molecular gas. The median dust mass is

Star formation in half of massive galaxies was quenched by the time the Universe was 3 billion years old

We present SPECULATOR - a fast, accurate, and flexible framework for emulating stellar population synthesis (SPS) models for predicting galaxy spectra and photometry. For emulating spectra, we use principal component analysis to construct a set of basis functions, and neural networks to learn the basis coefficients as a function of the SPS model parameters. For photometry, we parameterize the magnitudes (for the filters of interest) as a function of SPS parameters by a neural network. The resulting emulators are able to predict spectra and photometry under both simple and complicated SPS model parameterizations to percent-level accuracy, giving a factor of $10^3$-$10^4$ speed up over direct SPS computation. They have readily-computable derivatives, making them amenable to gradient-based inference and optimization methods. The emulators are also straightforward to call from a GPU, giving an additional order-of-magnitude speed-up. Rapid SPS computations delivered by emulation offers a massive reduction in the computational resources required to infer the physical properties of galaxies from observed spectra or photometry and simulate galaxy populations under SPS models, whilst maintaining the accuracy required for a range of applications., Comment: 15 pages, 9 figures, accepted by ApJS April 2020

Using a local reference sample of 21 galaxies, we compare observations of the $��$2.16 $��$m Brackett-$��$ (Br$��$) hydrogen recombination line with predictions from the Prospector Bayesian inference framework, which was used to fit the broadband photometry of these systems. This is a clean test of the spectral-energy-distribution-derived star formation rates (SFRs), as dust is expected to be optically thin at this wavelength in nearly all galaxies; thus, the internal conversion of SFR to predicted line luminosity does not depend strongly on the adopted dust model and posterior dust parameters, as is the case for shorter-wavelength lines such as H$��$. We find that Prospector predicts Br$��$ luminosities and equivalent widths with small offsets ($\sim$0.05 dex), and scatter ($\sim$0.2 dex), consistent with measurement uncertainties, though we caution that the derived offset is dependent on the choice of stellar isochrones. We demonstrate that even when the Prospector-derived dust attenuation does not well describe, e.g., H$��$ line properties or observed reddening between H$��$ and Br$��$, the underlying SFRs are accurate, as verified by the dust-free Br$��$ comparison. Finally, we discuss in what ways Br$��$ might be able to help constrain model parameters when treated as an input to the model, and comment on its potential as an accurate monochromatic SFR indicator in the era of JWST multiobject near-IR spectroscopy., 18 pages, 9 figures, 2 tables. Version accepted in ApJ

We present optical observations of the Swift short-duration gamma-ray burst (GRB) GRB 161104A and its host galaxy at $z=0.793 \pm 0.003$. We model the multiband photometry and spectroscopy with the stellar population inference code Prospector, and explore the posterior using nested sampling. We find that the mass-weighted age $t_m = 2.12^{+0.23}_{-0.21}$~Gyr, stellar mass $\log{(M/M_\odot)} = 10.21 \pm 0.04$, metallicity $\log{(Z/Z_\odot)} = 0.08^{+0.05}_{-0.06}$, dust extinction $A_V = 0.08^{+0.08}_{-0.05}$ mag, and the star formation rate $\text{SFR} = 9.9 \times 10^{-2} M_\odot$~yr$^{-1}$. These properties, along with a prominent 4000 Angstrom break and optical absorption lines classify this host as an early-type, quiescent galaxy. Using Dark Energy Survey galaxy catalogues, we demonstrate that the host of GRB 161104A resides on the outskirts of a galaxy cluster at $z\approx 0.8$, situated $\approx 1$ Mpc from the likely brightest cluster galaxy. We also present new modeling for 20 additional short GRB hosts ($\approx33\%$ of which are early-type galaxies), finding population medians of $\log(M/M_\odot) = 9.94^{+0.88}_{-0.98}$ and $t_m = 1.07^{+1.98}_{-0.67}$~Gyr ($68\%$ confidence). We further find that the host of GRB 161104A is more distant, less massive, and younger than the four other short GRB hosts known to be associated with galaxy clusters. Cluster short GRBs have faint afterglows, in the lower $\approx 11\%$ ($\approx 30\%$) of observed X-ray (optical) luminosities. We place a lower limit on the fraction of short GRBs in galaxy clusters versus those in the field of $\approx 5-13\%$, consistent with the fraction of stellar mass $\approx 10-20\%$ in galaxy clusters at redshifts $0.1 \leq z \leq 0.8$., Comment: 20 pages, 9 figures, ApJ: Vol. 904, No. 1

We present a novel Bayesian methodology to jointly model photometry and deep Hubble Space Telescope 2D grism spectroscopy of high-redshift galaxies. Our requiem2d code measures both unresolved and resolved stellar populations, ages, and star formation histories (SFHs) for the ongoing REsolving QUIEscent Magnified (REQUIEM) Galaxies Survey, which targets strong gravitationally lensed quiescent galaxies a z similar to 2. We test the accuracy of requiem2d using a simulated sample of massive galaxies at z similar to 2 from the Illustris cosmological simulation and find that we recover the general trends in SFH and median stellar ages. We further present a pilot study for the REQUIEM Galaxies Survey: MRG-S0851, a quintuply imaged, massive (logM/M-circle dot = 11.02 +/- 0.04) red galaxy at z = 1.883.+/-.0.001. With an estimated gravitational magnification of mu = 5.7(-0.2)(+0.4), we sample the stellar populations on 0.6 kpc physical size bins. The global mass-weighted median age is constrained to be 1.8(-0.2)(+0.3) Gyr, and our spatially resolved analysis reveals that MRG-S0851 has a flat age gradient in the inner 3 kpc core after taking into account the subtle effects of dust and metallicity on age measurements, favoring an early formation scenario. The analysis for the full REQUIEM-2D sample will be presented in a forthcoming paper with a beta release of the requiem2d code.

We present BVRI and unfiltered light curves of 93 Type Ia supernovae (SNe Ia) from the Lick Observatory Supernova Search (LOSS) follow-up program conducted between 2005 and 2018. Our sample consists of 78 spectroscopically normal SNe Ia, with the remainder divided between distinct subclasses (three SN 1991bg-like, three SN 1991T-like, four SNe Iax, two peculiar, and three super-Chandrasekhar events), and has a median redshift of 0.0192. The SNe in our sample have a median coverage of 16 photometric epochs at a cadence of 5.4 days, and the median first observed epoch is ~4.6 days before maximum B-band light. We describe how the SNe in our sample are discovered, observed, and processed, and we compare the results from our newly developed automated photometry pipeline to those from the previous processing pipeline used by LOSS. After investigating potential biases, we derive a final systematic uncertainty of 0.03 mag in BVRI for our dataset. We perform an analysis of our light curves with particular focus on using template fitting to measure the parameters that are useful in standardising SNe Ia as distance indicators. All of the data are available to the community, and we encourage future studies to incorporate our light curves in their analyses., Comment: 29 pages, 13 figures, accepted for publication in MNRAS

There has been a long-standing factor-of-two tension between the observed star formation rate density and the observed stellar mass buildup after $z\sim2$. Recently we have proposed that sophisticated panchromatic SED models can resolve this tension, as these methods infer systematically higher masses and lower star formation rates than standard approaches. In a series of papers we now extend this analysis and present a complete, self-consistent census of galaxy formation over 0.2 < z < 3 inferred with the Prospector galaxy SED-fitting code. In this work, Paper I, we present the evolution of the galaxy stellar mass function using new mass measurements of ~10^5 galaxies in the 3D-HST and COSMOS-2015 surveys. We employ a new methodology to infer the mass function from the observed stellar masses: instead of fitting independent mass functions in a series of fixed redshift intervals, we construct a continuity model that directly fits for the redshift evolution of the mass function. This approach ensures a smooth picture of galaxy assembly and makes use of the full, non-Gaussian uncertainty contours in our stellar mass inferences. The resulting mass function has higher number densities at a fixed stellar mass than almost any other measurement in the literature, largely owing to the older stellar ages inferred by Prospector. The stellar mass density is ~50% higher than previous measurements, with the offset peaking at z~1. This talk describes a paper submitted for publication in the Astrophysical Journal.

We analyze the distribution of rest-frame U-V and V-J colors for star-forming galaxies at 0.5 < z < 2.5. Using stellar population synthesis, stochastic star formation histories, and a simple prescription for the dust attenuation that accounts for the shape and inclination of galaxies, we construct a model for the distribution of galaxy colors. With only two free parameters, this model is able to reproduce the observed galaxy colors as a function of redshift and stellar mass remarkably well. Our analysis suggests that the wide range of dust attenuation values measured for star-forming galaxies at a given redshift and stellar mass is almost entirely due to the effect of inclination; if all galaxies were observed edge-on, they would show very similar dust attenuation. This result has important implications for the interpretation of dust attenuation measurements, the treatment of UV and IR luminosity, and the comparison between numerical simulations and observations., 7 pages, 5 figures

Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a "normal" molecular-gas to dust mass ratio, $\delta_{\mathrm{GDR}}$, typically of order one hundred. Recent null detections of quiescent galaxies in deep dust continuum observations support a picture where the cold gas and dust has been rapidly depleted or expelled. In this work, we present another viable explanation: a significant fraction of galaxies with low star formation per unit stellar mass are predicted to have extreme $\delta_{\mathrm{GDR}}$ ratios. We show that simulated massive quiescent galaxies at $0 < z < 3$ in the \textsc{simba} cosmological simulations have $\delta_{\mathrm{GDR}}$ values that extend $>$4 orders of magnitude. The dust in most simulated quiescent galaxies is destroyed significantly more rapidly than the molecular gas depletes, and cannot be replenished. The transition from star-forming to quiescent halts dust formation via star formation processes, with dust subsequently destroyed by supernova shocks and thermal sputtering of dust grains embedded in hot plasma. After this point, the dust growth rate in the models is not sufficient to overcome the loss of $>$3 orders of magnitude in dust mass to return to normal values of $\delta_{\mathrm{GDR}}$ despite having high metallicity. Our results indicate that it is not straight forward to use a single observational indicator to robustly pre-select exotic versus normal ratios. These simulations make strong predictions that can be tested with millimeter facilities., Comment: 9 pages, 5 figures, version accepted for publication in the Astrophysical Journal Letters

Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST., 21 Pages

The UVJ color-color diagram is a popular and efficient method to distinguish between quiescent and star-forming galaxies through their rest-frame $U-V$ vs. $V-J$ colors. Here we explore the information content of this color-color space using the Bayesian inference machine Prospector. We fit the same physical model to two datasets: (i) UVJ fluxes alone, and (ii) full UV-mid IR (MIR) broadband SEDs from the 3D-HST survey. Notably this model uses both nonparametric SFHs and a flexible dust attenuation curve, both of which have the potential to `break' the typical correlations observed in UVJ color-color space. Instead, these fits confirm observed trends between UVJ colors and observed galaxy properties, including specific star formation rate (sSFR), dust attenuation, stellar age, and stellar metallicity. They also demonstrate that UVJ colors do not, on their own, constrain stellar age or metallicity; the observed trends in the UVJ diagram are instead driven by galaxy scaling relationships and thus will evolve with cosmological time. We also show that UVJ colors 'saturate' below $\log(\mathrm{sSFR/yr}^{-1})\lesssim -10.5$, i.e. changing sSFR no longer produces substantial changes in UVJ colors. We show that far-UV and/or MIR fluxes continue to correlate with sSFR down to low sSFRs and can be used in color-color diagrams to efficiently target galaxies with much lower levels of ongoing star formation. We provide selection criteria in these new color-color spaces as a function of desired sample sSFR., Comment: update to published version

Modelling reionization often requires significant assumptions about the properties of ionizing sources. Here, we infer the total output of hydrogen-ionizing photons (the ionizing emissivity, $\dot{N}_\textrm{ion}$) at $z=4-14$ from current reionization constraints, being maximally agnostic to the properties of ionizing sources. We use a Bayesian analysis to fit for a non-parametric form of $\dot{N}_\textrm{ion}$, allowing us to flexibly explore the entire prior volume. We infer a declining $\dot{N}_\textrm{ion}$ with redshift at $z>6$, which can be used as a benchmark for reionization models. Model-independent reionization constraints from the CMB optical depth and Ly$��$ and Ly$��$ forest dark pixel fraction produce $\dot{N}_\textrm{ion}$ evolution ($d\log_{10}\dot{N}_\textrm{ion}/dz|_{z=6\rightarrow8} = -0.31\pm0.35$ dex) consistent with the declining UV luminosity density of galaxies, assuming constant ionizing photon escape fraction and efficiency. Including measurements from Ly$��$ damping of galaxies and quasars produces a more rapid decline: $d\log_{10}\dot{N}_\textrm{ion}/dz|_{z=6\rightarrow8} =-0.44\pm0.22$ dex, steeper than the declining galaxy luminosity density (if extrapolated beyond $M_\mathrm{UV} \lesssim -13$), and constrains the mid-point of reionization to $z = 6.93\pm0.14$., submitted to MNRAS

We present optical photometry and spectroscopy of SN 2016iet, an unprecedented Type I supernova (SN) at $z=0.0676$ with no obvious analog in the existing literature. The peculiar light curve has two roughly equal brightness peaks ($\approx -19$ mag) separated by 100 days, and a subsequent slow decline by 5 mag in 650 rest-frame days. The spectra are dominated by emission lines of calcium and oxygen, with a width of only $3400$ km s$^{-1}$, superposed on a strong blue continuum in the first year, and with a large ratio of $L_{\rm [Ca\,II]}/L_{\rm [O\,I]}\approx 4$ at late times. There is no clear evidence for hydrogen or helium associated with the SN at any phase. We model the light curves with several potential energy sources: radioactive decay, central engine, and circumstellar medium (CSM) interaction. Regardless of the model, the inferred progenitor mass near the end of its life (i.e., CO core mass) is $\gtrsim 55$ M$_\odot$ and up to $120$ M$_\odot$, placing the event in the regime of pulsational pair instability supernovae (PPISNe) or pair instability supernovae (PISNe). The models of CSM interaction provide the most consistent explanation for the light curves and spectra, and require a CSM mass of $\approx 35$ M$_\odot$ ejected in the final decade before explosion. We further find that SN 2016iet is located at an unusually large offset ($16.5$ kpc) from its low metallicity dwarf host galaxy ($Z\approx 0.1$ Z$_\odot$, $M\approx 10^{8.5}$ M$_\odot$), supporting the PPISN/PISN interpretation. In the final spectrum, we detect narrow H$\alpha$ emission at the SN location, likely due to a dim underlying galaxy host or an H II region. Despite the overall consistency of the SN and its unusual environment with PPISNe and PISNe, we find that the inferred properties of SN\,2016iet challenge existing models of such events., Comment: 26 Pages, 17 Figures, Submitted to ApJ

There has been a long-standing factor-of-two tension between the observed star formation rate density and the observed stellar mass buildup after $z\sim2$. Recently we have proposed that sophisticated panchromatic SED models can resolve this tension, as these methods infer systematically higher masses and lower star formation rates than standard approaches. In a series of papers we now extend this analysis and present a complete, self-consistent census of galaxy formation over $0.2 < z < 3$ inferred with the \texttt{Prospector} galaxy SED-fitting code. In this work, Paper I, we present the evolution of the galaxy stellar mass function using new mass measurements of $\sim$10$^5$ galaxies in the 3D-HST and COSMOS-2015 surveys. We employ a new methodology to infer the mass function from the observed stellar masses: instead of fitting independent mass functions in a series of fixed redshift intervals, we construct a continuity model that directly fits for the redshift evolution of the mass function. This approach ensures a smooth picture of galaxy assembly and makes use of the full, non-Gaussian uncertainty contours in our stellar mass inferences. The resulting mass function has higher number densities at a fixed stellar mass than almost any other measurement in the literature, largely owing to the older stellar ages inferred by \texttt{Prospector}. The stellar mass density is $\sim$50% higher than previous measurements, with the offset peaking at $z\sim1$. The next two papers in this series will present the new measurements of star-forming main sequence and the cosmic star formation rate density, respectively., 22 pages, 11 figures