Researcher profile

Joel Leja

Joel Leja contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
24works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2026arXiv

The X-Ray Dot: Exotic Dust or a Late-Stage Little Red Dot?

JWST's "Little Red Dots" (LRDs) are increasingly interpreted as active galactic nuclei (AGN) obscured by dense thermalized gas rather than dust as evidenced by their X-ray weakness, blackbody-like continua, and Balmer line profiles. A key question is how LRDs connect to standard UV-luminous AGN and whether transitional phases exist and if they are observable. We present the "X-Ray Dot" (XRD), a compact source at $z=3.28$ observed by the NIRSpec WIDE GTO survey. The XRD exhibits LRD hallmarks: a blackbody-like ($T_{\rm eff} \simeq 6400\,$K) red continuum, a faint but blue rest-UV excess, falling mid-IR emission, and broad Balmer lines ($\rm FWHM \sim 2700-3200\,km\,s^{-1}$). Unlike LRDs, however, it is remarkably X-ray luminous ($L_\textrm{2$-$10$\,$keV} = 10^{44.18}\,$erg$\,$s$^{-1}$) and has a continuum inflection that is bluewards of the Balmer limit. We find that the red rest-optical and blue mid-IR continuum cannot be reproduced by standard dust-attenuated AGN models without invoking extremely steep extinction curves, nor can the weak mid-IR emission be reconciled with well-established X-ray--torus scaling relations. We therefore consider an alternative scenario: the XRD may be an LRD in transition, where the gas envelope dominates the optical continuum but optically thin sightlines allow X-rays to escape. The XRD may thus provide a physical link between LRDs and standard AGN, offering direct evidence that LRDs are powered by supermassive black holes and providing insight into their accretion properties.

preprint2023arXiv

Early JWST imaging reveals strong optical and NIR color gradients in galaxies at $z\sim2$ driven mostly by dust

Recent studies have shown that galaxies at cosmic noon are redder in the center and bluer in the outskirts, mirroring results in the local universe. These color gradients could be caused by either gradients in the stellar age or dust opacity; however, distinguishing between these two causes is impossible with rest-frame optical photometry alone. Here we investigate the underlying causes of the gradients from spatially-resolved rest-frame $U-V$ vs. $V-J$ color-color diagrams, measured from early observations with the James Webb Space Telescope. We use $1\, μm - 4\, μm$ NIRCam photometry from the CEERS survey of a sample of 54 galaxies with $M_* / M_\odot>10$ at redshifts $1.7<z<2.3$ selected from the 3D-HST catalog. We model the light profiles in the F115W, F200W and F356W NIRCam bands using \texttt{imcascade}, a Bayesian implementation of the Multi-Gaussian expansion (MGE) technique which flexibly represents galaxy profiles using a series of Gaussians. We construct resolved rest-frame $U-V$ and $V-J$ color profiles. The majority of star-forming galaxies have negative gradients (i.e. redder in the center, bluer in the outskirts) in both $U-V$ and $V-J$ colors consistent with radially decreasing dust attenuation. A smaller population (roughly 15\%) of star-forming galaxies have positive $U-V$ but negative $V-J$ gradients implying centrally concentrated star-formation. For quiescent galaxies we find a diversity of UVJ color profiles, with roughly one-third showing star-formation in their center. This study showcases the potential of JWST to study the resolved stellar populations of galaxies at cosmic noon.

preprint2023arXiv

Flexible Models for Galaxy Star Formation Histories Both Shift and Scramble the Optical Color-M/L Relationship

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 SED model with a more sophisticated 14-dimensional 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.

preprint2023arXiv

Molecular Gas Reservoirs in Massive Quiescent Galaxies at $\mathrm{z\sim0.7}$ Linked to Late Time Star Formation

We explore how the presence of detectable molecular gas depends on the inferred star formation histories (SFHs) in 8 massive, quiescent galaxies at $\mathrm{z\sim0.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 leftover 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 $\mathrm{f_{burst}\approx0.3-6\%}$, and ended between $\mathrm{t_{end\mbox{-}burst}\approx0-330~Myr}$ ago. The CO-detected galaxies form a higher fraction of mass in the last Gyr ($\mathrm{f_{M_{1Gyr}}=2.6\pm1.8\%}$) compared to the CO-undetected galaxies ($\mathrm{f_{M_{1Gyr}}=0.2\pm0.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 $\mathrm{SFR_{UV+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.

preprint2023arXiv

Star formation histories of UV-luminous galaxies at $z \simeq 6.8$: implications for stellar mass assembly at early cosmic times

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.

preprint2023arXiv

UNCOVER: Illuminating the Early Universe -- JWST/NIRSpec Confirmation of $z > 12$ Galaxies

Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection at $z_{\rm spec} = 12.393^{+0.004}_{-0.001}$, and a plausible candidate at $z_{\rm spec} = 13.079^{+0.013}_{-0.001}$. The galaxies are discovered in JWST/NIRCam imaging and their distances are inferred with JWST/NIRSpec spectroscopy, all from the JWST Cycle 1 UNCOVER Treasury survey. Detailed stellar population modeling using JWST NIRCam and NIRSpec data corroborates the primeval characteristics of these galaxies: low mass ($\sim 10^8~{\rm M_\odot}$), young, rapidly-assembling, metal-poor, and star-forming. Interestingly, both galaxies are spatially resolved, having lensing-corrected rest-UV effective radii on the order of 300-400 pc, which are notably larger than other spectroscopically confirmed systems at similar redshifts. The observed dynamic range of $z \gtrsim 10$ sizes spans over 1 order of magnitude, implying a significant scatter in the size-mass relation at early times. Deep into the epoch of reionization, these discoveries elucidate the emergence of the first galaxies.

preprint2022arXiv

A Bayesian Population Model for the Observed Dust Attenuation in Galaxies

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 ($τ$), though less so than in previous studies. $τ$ increases strongly with SFR, though when $\log~{\rm SFR}\lesssim 0$, $τ$ remains roughly constant over a wide range of stellar masses. Edge-on galaxies tend to have larger $τ$ 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.

preprint2022arXiv

Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8

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 &#34;late bloomers&#34; are uncommon. In summary, our analysis supports the &#34;grow & quench&#34; framework and is consistent with a wide and continuously-populated diversity of quenching timescales.

preprint2022arXiv

Hierarchical Bayesian inference of photometric redshifts with stellar population synthesis models

We present a Bayesian hierarchical framework to analyze photometric galaxy survey data with stellar population synthesis (SPS) models. Our method couples robust modeling of spectral energy distributions with a population model and a noise model to characterize the statistical properties of the galaxy populations and real observations, respectively. By self-consistently inferring all model parameters, from high-level hyper-parameters to SPS parameters of individual galaxies, one can separate sources of bias and uncertainty in the data.We demonstrate the strengths and flexibility of this approach by deriving accurate photometric redshifts for a sample of spectroscopically-confirmed galaxies in the COSMOS field, all with 26-band photometry and spectroscopic redshifts. We achieve a performance competitive with publicly-released photometric redshift catalogs based on the same data. Prior to this work, this approach was computationally intractable in practice due to the heavy computational load of SPS model calls; we overcome this challenge using with neural emulators. We find that the largest photometric residuals are associated with poor calibration for emission line luminosities and thus build a framework to mitigate these effects. This combination of physics-based modeling accelerated with machine learning paves the path towards meeting the stringent requirements on the accuracy of photometric redshift estimation imposed by upcoming cosmological surveys. The approach also has the potential to create new links between cosmology and galaxy evolution through the analysis of photometric datasets.

preprint2022arXiv

How Well Can We Measure Galaxy Dust Attenuation Curves? The Impact of the Assumed Star-Dust Geometry Model in SED Fitting

One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the Simba cosmological simulation with the Prospector SED fitting code. We consider the impact of the commonly-assumed uniform screen model and introduce a new non-uniform screen model parameterized by the fraction of unobscured stellar light. This non-uniform screen model allows for a non-zero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a non-uniform screen model, with the median offset in UV attenuation decreasing from $-0.30$ dex with a uniform screen model to $-0.17$ dex with the non-uniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the non-uniform screen model, decreasing the SFR offset on average by $0.12$ dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations.

preprint2022arXiv

Recovering the star formation histories of recently-quenched galaxies: the impact of model and prior choices

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 &#34;non-parametric&#34; 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.

preprint2022arXiv

Schrodinger&#39;s Galaxy Candidate: Puzzlingly Luminous at $z\approx17$, or Dusty/Quenched at $z\approx5$?

$JWST$&#39;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μ$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 $Λ$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy&#39;s nearest neighbors at $<2.5$&#34; have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749&#39;s photometry. The emission lines at $z\approx5$ conspire to boost the $>2μ$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 ($Δ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.

preprint2022arXiv

Spectral Energy Distributions in Three Deep-Drilling Fields of the Vera C. Rubin Observatory Legacy Survey of Space and Time: Source Classification and Galaxy Properties

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.

preprint2022arXiv

The Lick Observatory Supernova Search follow-up program: photometry data release of 70 stripped-envelope supernovae

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.

preprint2021arXiv

Deep Hubble Space Telescope Observations of GW170817: Complete Light Curves and the Properties of the Galaxy Merger of NGC 4993

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 $γ$-ray burst (GRB) 170817A afterglow from 109--170 rest-frame days post-merger. These detections are consistent with the previously observed $β=-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.

preprint2021arXiv

Diagnosing DASH: A Catalog of Structural Properties for the COSMOS-DASH Survey

We present the $H_{160}$ morphological catalogs for the COSMOS-DASH survey, the largest area near-IR survey using HST-WFC3 to date. Utilizing the &#34;Drift And SHift&#34; observing technique for HST-WFC3 imaging, the COSMOS-DASH survey imaged approximately 0.5 deg$^2$ of the UltraVISTA deep stripes (0.7 deg$^2$ when combined with archival data). Global structural parameters are measured for 51,586 galaxies within COSMOS-DASH using GALFIT (excluding the CANDELS area) with detection using a deep multi-band HST image. We recover consistent results with those from the deeper 3D-HST morphological catalogs, finding that, in general, sizes and Sérsic indices of typical galaxies are accurate to limiting magnitudes of $H_{160}<23$ and $H_{160}<22$ ABmag, respectively. In size-mass parameter space, galaxies in COSMOS-DASH demonstrate robust morphological measurements out to $z\sim2$ and down to $\log(M_\star/M_\odot)\sim9$. With the advantage of the larger area of COSMOS-DASH, we measure a flattening of the quiescent size-mass relation below $\log(M_\star/M_\odot)\sim10.5$ that persists out to $z\sim2$. We show that environment is not the primary driver of this flattening, at least out to $z=1.2$, whereas internal physical processes may instead govern the structural evolution.

preprint2021arXiv

SQuIGGLE: Studying Quenching in Intermediate-z Galaxies -- Gas, AnguLar Momentum, and Evolution

We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.5<z<~0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be useful analogs for high-redshift quenching galaxies, but low enough redshift that multi-wavelength follow-up observations are feasible with modest telescope investments. We use the Prospector code to infer the stellar population properties and non-parametric star formation histories of all galaxies in the sample. We find that SQuIGGLE galaxies are both very massive (M* ~ 10^11.25 Msun) and quenched, with inferred star formation rates <~1Msun/yr, more than an order of magnitude below the star-forming main sequence. The best-fit star formation histories confirm that these galaxies recently quenched a major burst of star formation: >75% of SQuIGGLE galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200Myr before observation. We find that SQuIGGLE galaxies are on average younger and more burst-dominated than most other z<~1 post-starburst samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full SQuIGGLE survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, AGN incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process.

preprint2021arXiv

The Diverse Molecular Gas Content of Massive Galaxies Undergoing Quenching at z~1

We present a detailed study of the molecular gas content and stellar population properties of three massive galaxies at 1 < z < 1.3 that are in different stages of quenching. The galaxies were selected to have a quiescent optical/near-infrared spectral energy distribution and a relatively bright emission at 24 micron, and show remarkably diverse properties. CO emission from each of the three galaxies is detected in deep NOEMA observations, allowing us to derive molecular gas fractions Mgas/Mstar of 13-23%. We also reconstruct the star formation histories by fitting models to the observed photometry and optical spectroscopy, finding evidence for recent rejuvenation in one object, slow quenching in another, and rapid quenching in the third system. To better constrain the quenching mechanism we explore the depletion times for our sample and other similar samples at z~0.7 from the literature. We find that the depletion times are highly dependent on the method adopted to measure the star formation rate: using the UV+IR luminosity we obtain depletion times about 6 times shorter than those derived using dust-corrected [OII] emission. When adopting the star formation rates from spectral fitting, which are arguably more robust, we find that recently quenched galaxies and star-forming galaxies have similar depletion times, while older quiescent systems have longer depletion times. These results offer new, important constraints for physical models of galaxy quenching.

preprint2021arXiv

Ubiquitous [OII] Emission in Quiescent Galaxies at z ~ 0.85

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.

preprint2020arXiv

A New Census of the 0.2 < z < 3.0 Universe, Part I: The Stellar Mass Function

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.

preprint2020arXiv

Brackett-$γ$ as a Gold-standard Test of Star Formation Rates Derived from SED Fitting

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.

preprint2020arXiv

Recent Star-formation in a Massive Slowly-Quenched Lensed Quiescent Galaxy at z=1.88

In this letter, we reconstruct the formation pathway of MRG-S0851, a massive, $\log M_*/M_\odot=11.02\pm0.04$, strongly lensed, red, galaxy at $z=1.883\pm0.001$. While the global photometry and spatially-resolved outskirts of MRG-S0851 imply an early-formation scenario with a slowly decreasing or constant star-formation history, a joint fit of 2D grism spectroscopy and photometry reveals a more complex scenario: MRG-S0851 is likely to be experiencing a centrally-concentrated rejuvenation in the inner $\sim$1 kpc in the last $\sim$100 Myr of evolution. We estimate $0.5\pm0.1\%$ of the total stellar mass is formed in this phase. Rejuvenation episodes are suggested to be infrequent for massive galaxies at $z\sim2$, but as our analyses indicate, more examples of complex star-formation histories may yet be hidden within existing data. By adding a FUV color criterion to the standard U-V/V-J diagnostic, thereby heightening our sensitivity to recent star formation, we show that we can select populations of galaxies with similar spectral energy distributions to that of MRG-S0851, but note that deep follow-up spectroscopic observations and/or spatially resolved analyses are necessary to robustly confirm the rejuvenation of these candidates. Using our criteria with MRG-S0851 as a prototype, we estimate that $\sim$1\% of massive quiescent galaxies at $1<z<2$ are potentially rejuvenating.

preprint2020arXiv

SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry

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.

preprint2020arXiv

The GOGREEN survey: Post-infall environmental quenching fails to predict the observed age difference between quiescent field and cluster galaxies at z>1

We study the star formation histories (SFHs) and mass-weighted ages of 331 UVJ-selected quiescent galaxies in 11 galaxy clusters and in the field at 1<z<1.5 from the Gemini Observations of Galaxies in Rich Early ENvironments (GOGREEN) survey. We determine the SFHs of individual galaxies by simultaneously fitting rest-frame optical spectroscopy and broadband photometry to stellar population models. We confirm that the SFHs are consistent with more massive galaxies having on average earlier formation times. Comparing galaxies found in massive clusters with those in the field, we find galaxies with $M_\ast<10^{11.3}$ M$_{\odot}$ in the field have more extended SFHs. From the SFHs we calculate the mass-weighted ages, and compare age distributions of galaxies between the two environments, at fixed mass. We constrain the difference in mass-weighted ages between field and cluster galaxies to $0.31_{^{-0.33}}^{_{+0.51}}$ Gyr, in the sense that cluster galaxies are older. We place this result in the context of two simple quenching models and show that neither environmental quenching based on time since infall (without pre-processing) nor a difference in formation times alone can reproduce both the average age difference and relative quenched fractions. This is distinctly different from local clusters, for which the majority of the quenched population is consistent with having been environmentally quenched upon infall. Our results suggest that quenched population in galaxy clusters at z>1 has been driven by different physical processes than those at play at z=0.