Showing 1–39 of 39 results for author: Iyer, K G

When, Where, and How Star Formation Happens in a Galaxy Pair at Cosmic Noon Using CANUCS JWST/NIRISS Grism Spectroscopy

Authors: Vicente Estrada-Carpenter, Marcin Sawicki, Gabe Brammer, Guillaume Desprez, Roberto Abraham, Yoshihisa Asada, Maruša Bradač, Kartheik G. Iyer, Nicholas S. Martis, Jasleen Matharu, Lamiya Mowla, Adam Muzzin, Gaël Noirot, Ghassan T. E. Sarrouh, Victoria Strait, Chris J. Willott

Abstract: Spatially resolved studies are key to understanding when, where, and how stars form within galaxies. Using slitless grism spectra and broadband imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS) we study the spatially resolved properties of a strongly lensed ($μ$ = 5.4$\pm$1.8) z = 0.8718 galaxy pair consisting of a blue face-on galaxy (10.2 $\pm$ 0.2 log($M/M_\odot$)) with multiple… ▽ More Spatially resolved studies are key to understanding when, where, and how stars form within galaxies. Using slitless grism spectra and broadband imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS) we study the spatially resolved properties of a strongly lensed ($μ$ = 5.4$\pm$1.8) z = 0.8718 galaxy pair consisting of a blue face-on galaxy (10.2 $\pm$ 0.2 log($M/M_\odot$)) with multiple star-forming clumps and a dusty red edge-on galaxy (9.9 $\pm$ 0.3 log($M/M_\odot$)). We produce accurate H$α$ maps from JWST/NIRISS grism data using a new methodology that accurately models spatially varying continuum and emission line strengths. With spatially resolved indicators, we probe star formation on timescales of $\sim$ 10 Myr (NIRISS H$α$ emission line maps) and $\sim$ 100 Myr (UV imaging and broadband SED fits). Taking the ratio of the H$α$ to UV flux ($η$), we measure spatially resolved star formation burstiness. We find that in the face-on galaxy both H$α$ and broadband star formation rates (SFRs) drop at large galactocentric radii by a factor of $\sim$ 4.7 and 3.8 respectively, while SFR over the last $\sim$ 100 Myrs has increased by a factor of 1.6. Additionally, of the 20 clumps identified in the galaxy pair we find that 7 are experiencing bursty star formation, while 10 clumps are quenching and 3 are in equilibrium (either being in a state of steady star formation or post-burst). Our analysis reveals that the blue face-on galaxy disk is predominantly in a quenching or equilibrium phase. However, the most intense quenching within the galaxy is seen in the quenching clumps. This pilot study demonstrates what JWST/NIRISS data can reveal about spatially varying star formation in galaxies at Cosmic Noon. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: MNRAS Accepted, 15 pages, 12 figures

CANUCS: Constraining the MACS J0416.1-2403 Strong Lensing Model with JWST NIRISS, NIRSpec and NIRCam

Authors: Gregor Rihtaršič, Maruša Bradač, Guillaume Desprez, Anishya Harshan, Gaël Noirot, Vicente Estrada-Carpenter, Nicholas S. Martis, Roberto G. Abraham, Yoshihisa Asada, Gabriel Brammer, Kartheik G. Iyer, Jasleen Matharu, Lamiya Mowla, Adam Muzzin, Ghassan T. E. Sarrouh, Marcin Sawicki, Victoria Strait, Chris J. Willott, Rachel Gledhill, Vladan Markov, Roberta Tripodi

Abstract: Strong gravitational lensing in galaxy clusters has become an essential tool in astrophysics, allowing us to directly probe the dark matter distribution and study magnified background sources. The precision and reliability of strong lensing models rely heavily on the number and quality of multiple images of background sources with spectroscopic redshifts. We present an updated strong lensing model… ▽ More Strong gravitational lensing in galaxy clusters has become an essential tool in astrophysics, allowing us to directly probe the dark matter distribution and study magnified background sources. The precision and reliability of strong lensing models rely heavily on the number and quality of multiple images of background sources with spectroscopic redshifts. We present an updated strong lensing model of the galaxy cluster MACS J0416.1-2403 with the largest sample of multiple images with spectroscopic redshifts in a galaxy cluster field to date. Furthermore, we aim to demonstrate the effectiveness of JWST particularly its NIRISS camera, for strong lensing studies. We use the JWST 's NIRCam imaging and NIRSpec and NIRISS spectroscopy from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). The cluster mass model is constrained using Lenstool software. Our new dataset, used for constraining the lens model, comprises 303 secure multiple images from 111 background sources and includes systems with previously known MUSE redshift and systems for which we obtained spectroscopic redshift for the first time using NIRISS and NIRSpec spectroscopy. The total number of secure spectroscopic systems is >20% higher than in the previous strong lensing studies of this cluster. The derived strong lensing model can reproduce multiple images with the root-mean-square distance of 0.53''. We also provide a full catalogue with 415 multiple images, including less reliable candidates. We furthermore demonstrate the effectiveness of JWST particularly NIRISS, for strong lensing studies. As NIRISS F115W, F150W, and F200W grism spectroscopy captures at least two of the [OII] λ3727, [OIII] λλ4959, 5007, and Hα lines at 1<z<3 (a redshift range particularly relevant for strong lensing studies) without target pre-selection, it complements MUSE and NIRSpec observations extremely well. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: Submitted to A&A, 23 pages, 12 figures, 6 tables

Stochastic prior for non-parametric star-formation histories

Authors: Jenny T. Wan, Sandro Tacchella, Benjamin D. Johnson, Kartheik G. Iyer, Joshua S. Speagle, Roberto Maiolino

Abstract: The amount of power contained in the variations in galaxy star-formation histories (SFHs) across a range of timescales encodes key information about the physical processes which modulate star formation. Modelling the SFHs of galaxies as stochastic processes allows the relative importance of different timescales to be quantified via the power spectral density (PSD). In this paper, we build upon the… ▽ More The amount of power contained in the variations in galaxy star-formation histories (SFHs) across a range of timescales encodes key information about the physical processes which modulate star formation. Modelling the SFHs of galaxies as stochastic processes allows the relative importance of different timescales to be quantified via the power spectral density (PSD). In this paper, we build upon the PSD framework and develop a physically-motivated, "stochastic" prior for non-parametric SFHs in the spectral energy distribution (SED)-modelling code Prospector. We test this prior in two different regimes: 1) massive, $z = 0.7$ galaxies with both photometry and spectra, analogous to those observed with the LEGA-C survey, and 2) $z = 8$ galaxies with photometry only, analogous to those observed with NIRCam on JWST. We find that it is able to recover key galaxy parameters (e.g. stellar mass, stellar metallicity) to the same level of fidelity as the commonly-used continuity prior. Furthermore, the realistic variability information incorporated by the stochastic SFH model allows it to fit the SFHs of galaxies more accurately and precisely than traditional non-parametric models. In fact, the stochastic prior is $\gtrsim 2\times$ more accurate than the continuity prior in measuring the recent star-formation rates (log SFR$_{100}$ and log SFR$_{10}$) of both the $z = 0.7$ and $z = 8$ mock systems. While the PSD parameters of individual galaxies are difficult to constrain, the stochastic prior implementation presented in this work allows for the development hierarchical models in the future, i.e. simultaneous SED-modelling of an ensemble of galaxies to measure their underlying PSD. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: Comments welcome!

Katachi: Decoding the Imprints of Past Star Formation on Present Day Morphology in Galaxies with Interpretable CNNs

Authors: Juan Pablo Alfonzo, Kartheik G. Iyer, Masayuki Akiyama, Greg L. Bryan, Suchetha Cooray, Eric Ludwig, Lamiya Mowla, Kiyoaki C. Omori, Camilla Pacifici, Joshua S. Speagle, John F. Wu

Abstract: The physical processes responsible for shaping how galaxies form and quench over time leave imprints on both the spatial (galaxy morphology) and temporal (star formation history; SFH) tracers that we use to study galaxies. While the morphology-SFR connection is well studied, the correlation with past star formation activity is not as well understood. To quantify this we present Katachi, an interpr… ▽ More The physical processes responsible for shaping how galaxies form and quench over time leave imprints on both the spatial (galaxy morphology) and temporal (star formation history; SFH) tracers that we use to study galaxies. While the morphology-SFR connection is well studied, the correlation with past star formation activity is not as well understood. To quantify this we present Katachi, an interpretable convolutional neural network (CNN) framework that learns the connection between the factors regulating star formation in galaxies on different spatial and temporal scales. Katachi is trained on 9904 galaxies at 0.02$<$z$<$0.1 in the SDSS-IV MaNGA DR17 sample to predict stellar mass (M$_*$; RMSE 0.22 dex), current star formation rate (SFR; RMSE 0.31 dex) and half-mass time (t$_{50}$; RMSE 0.23 dex). This information allows us to reconstruct non-parametric SFHs for each galaxy from \textit{gri} imaging alone. To quantify the morphological features informing the SFH predictions we use SHAP (SHapley Additive exPlanations). We recover the expected trends of M$_*$ governed by the growth of galaxy bulges, and SFR correlating with spiral arms and other star-forming regions. We also find the SHAP maps of D4000 are more complex than those of M$_*$ and SFR, and that morphology is correlated with t$_{50}$ even at fixed mass and SFR. Katachi serves as a scalable public framework to predict galaxy properties from large imaging surveys including Rubin, Roman, and Euclid, with large datasets of high SNR imaging across limited photometric bands. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Comments: 35 pages, 25 figures. Accepted for publication in ApJ, comments are welcomed

Exposing Line Emission: A First Look At The Systematic Differences of Measuring Stellar Masses With JWST NIRCam Medium Versus Wide Band Photometry

Authors: Ghassan T. Sarrouh, Adam Muzzin, Kartheik G. Iyer, Lamiya Mowla, Roberto G. Abraham, Yoshihisa Asada, Marusa Bradac, Gabriel B. Brammer, Guillaume Desprez, Nicholas S. Martis, Jasleen Matharu, Gaël Noirot, Marcin Sawicki, Victoria Strait, Chris Willott, Johannes Zabl

Abstract: Photometrically derived stellar masses are known to suffer from systematic uncertainties, particularly due to nebular emission contributions to the spectral energy distribution. Using \emph{JWST} NIRCam imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS), we introduce a comparison study of photometrically-derived redshifts and stellar masses based on two photometric catalogs of the s… ▽ More Photometrically derived stellar masses are known to suffer from systematic uncertainties, particularly due to nebular emission contributions to the spectral energy distribution. Using \emph{JWST} NIRCam imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS), we introduce a comparison study of photometrically-derived redshifts and stellar masses based on two photometric catalogs of the same field spanning $\sim$0.4-4.5$μ$m: one consisting solely of wide band photometry, and another employing a combination of wide and medium band photometry. We find that \tilda70\% of galaxies have consistent photometric redshifts between both catalogs, with median stellar mass difference between the two catalogs of \lessthan\ 0.2 dex across all redshift bins. There are however a subset of galaxies (5\% at z\tilda2 up to 15\% at z\tilda6) where wide bands underestimate star formation rates and infer older stellar populations, leading to median stellar mass differences of \tilda0.7 dex. Examination of the SEDs for galaxies with inconsistent photometric redshifts shows this is caused by the inability of the wide bands to distinguish continuum emission from emission lines. Computing a stellar mass density with our sample we find that it is potentially underestimated using wide-band photometry by \tilda10-20\% at z \lessthan\ \ 4, and potentially overestimated by as much as a factor of 2-3 at z \greaterthan\ 5. These systematic differences caused by the poor spectral resolution of wide bands have implications for both ongoing and future planned observing programs which determine stellar mass and other physical properties of high redshift galaxies solely via wide band photometry. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 10 pages, 5 figures

MSC Class: 85-05

CEERS: Increasing Scatter along the Star-Forming Main Sequence Indicates Early Galaxies Form in Bursts

Authors: Justin W. Cole, Casey Papovich, Steven L. Finkelstein, Micaela B. Bagley, Mark Dickinson, Kartheik G. Iyer, L. Y. Aaron Yung, Laure Ciesla, Ricardo O. Amorin, Pablo Arrabal Haro, Rachana Bhatawdekar, Antonello Calabro, Nikko J. Cleri, Alexander de la Vega, Avishai Dekel, Ryan Endsley, Eric Gawiser, Mauro Giavalisco, Nimish P. Hathi, Michaela Hirschmann, Benne W. Holwerda, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Ray A. Lucas, Sara Mascia , et al. (7 additional authors not shown)

Abstract: We present the star-formation-rate -- stellar-mass (SFR-M$_\ast$) relation for galaxies in the CEERS survey at $4.5\leq z\leq 12$. We model the \jwst\ and \hst\ rest-UV and rest-optical photometry of galaxies with flexible star-formation histories (SFHs) using \bagpipes. We consider SFRs averaged from the SFHs over 10~Myr (\sfrten) and 100~Myr (\sfrcen), where the photometry probes SFRs on these t… ▽ More We present the star-formation-rate -- stellar-mass (SFR-M$_\ast$) relation for galaxies in the CEERS survey at $4.5\leq z\leq 12$. We model the \jwst\ and \hst\ rest-UV and rest-optical photometry of galaxies with flexible star-formation histories (SFHs) using \bagpipes. We consider SFRs averaged from the SFHs over 10~Myr (\sfrten) and 100~Myr (\sfrcen), where the photometry probes SFRs on these timescales, effectively tracing nebular emission lines in the rest-optical (on $\sim10$~Myr timescales) and the UV/optical continuum (on $\sim100$ Myr timescales). We measure the slope, normalization and intrinsic scatter of the SFR-M$_\ast$ relation, taking into account the uncertainty and the covariance of galaxy SFRs and $M_\ast$. From $z\sim 5-9$ there is larger scatter in the $\sfrten-M_\ast$ relation, with $σ(\log \sfrcen)=0.4$~dex, compared to the $\sfrcen-M_\ast$ relation, with $σ(\log \sfrten)=0.1$~dex. This scatter increases with redshift and increasing stellar mass, at least out to $z\sim 7$. These results can be explained if galaxies at higher redshift experience an increase in star-formation variability and form primarily in short, active periods, followed by a lull in star formation (i.e. ``napping'' phases). We see a significant trend in the ratio $R_\mathrm{SFR}=\log(\sfrten/\sfrcen)$ in which, on average, $R_\mathrm{SFR}$ decreases with increasing stellar mass and increasing redshift. This yields a star-formation ``duty cycle'' of $\sim40\%$ for galaxies with $\log M_\ast/M_\odot\geq 9.3$, at $z\sim5$, declining to $\sim20\%$ at $z\sim9$. Galaxies also experience longer lulls in star formation at higher redshift and at higher stellar mass, such that galaxies transition from periods of higher SFR variability at $z\gtrsim~6$ to smoother SFR evolution at $z\lesssim~4.5$. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 28 pages, 12 figures, 2 Appendix figures

Evidence for a Shallow Evolution in the Volume Densities of Massive Galaxies at $z=4$ to $8$ from CEERS

Authors: Katherine Chworowsky, Steven L. Finkelstein, Michael Boylan-Kolchin, Elizabeth J. McGrath, Kartheik G. Iyer, Casey Papovich, Mark Dickinson, Anthony J. Taylor, L. Y. Aaron Yung, Pablo Arrabal Haro, Micaela B. Bagley, Bren E. Backhaus, Rachana Bhatawdekar, Yingjie Cheng, Nikko J. Cleri, Justin W. Cole, M. C. Cooper, Luca Costantin, Avishai Dekel, Maximilien Franco, Seiji Fujimoto, Christopher C. Hayward, Benne W. Holwerda, Marc Huertas-Company, Michaela Hirschmann , et al. (14 additional authors not shown)

Abstract: We analyze the evolution of massive (log$_{10}$ [$M_\star/M_\odot$] $>10$) galaxies at $z \sim$ 4--8 selected from the JWST Cosmic Evolution Early Release Science (CEERS) survey. We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high redshift massive galaxies. Where available we inc… ▽ More We analyze the evolution of massive (log$_{10}$ [$M_\star/M_\odot$] $>10$) galaxies at $z \sim$ 4--8 selected from the JWST Cosmic Evolution Early Release Science (CEERS) survey. We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam wide-field slitless spectroscopy. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei (AGN). We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos for $z \sim$ 1--4. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 20 pages, 10 figures

A Steep Decline in the Galaxy Space Density Beyond Redshift 9 in the CANUCS UV Luminosity Function

Authors: Chris J. Willott, Guillaume Desprez, Yoshihisa Asada, Ghassan T. E. Sarrouh, Roberto Abraham, Maruša Bradač, Gabe Brammer, Vince Estrada-Carpenter, Kartheik G. Iyer, Nicholas S. Martis, Jasleen Matharu, Lamiya Mowla, Adam Muzzin, Gaël Noirot, Marcin Sawicki, Victoria Strait, Gregor Rihtaršič, Sunna Withers

Abstract: We present a new sample of 158 galaxies at redshift $z>7.5$ selected from deep \jwst\ NIRCam imaging of five widely-separated sightlines in the CANUCS survey. Two-thirds of the pointings and 80\% of the galaxies are covered by 12 to 14 NIRCam filters, including seven to nine medium bands, providing accurate photometric redshifts and robustness against low redshift interlopers. A sample of 28 galax… ▽ More We present a new sample of 158 galaxies at redshift $z>7.5$ selected from deep \jwst\ NIRCam imaging of five widely-separated sightlines in the CANUCS survey. Two-thirds of the pointings and 80\% of the galaxies are covered by 12 to 14 NIRCam filters, including seven to nine medium bands, providing accurate photometric redshifts and robustness against low redshift interlopers. A sample of 28 galaxies at $z>7.5$ with spectroscopic redshifts shows a low systematic offset and scatter in the difference between photometric and spectroscopic redshifts. We derive the galaxy UV luminosity function at redshifts 8 to 12, finding a slightly higher normalization than previously seen with \hst\ at redshifts 8 to 10. We observe a steeper decline in the galaxy space density from $z=8$ to $12$ than found by most \jwst\ Cycle 1 studies. In particular, we find only eight galaxies at $z>10$ and none at $z>12.5$, with no $z>10$ galaxies brighter than F277W AB=28 or $M_{\rm UV}=-20$ in our unmasked, delensed survey area of 53.4 square arcminutes. We attribute the lack of bright $z>10$ galaxies in CANUCS compared to GLASS and CEERS to intrinsic variance in the galaxy density along different sightlines. The evolution in the CANUCS luminosity function between $z=8$ and $12$ is comparable to that predicted by simulations that assume a standard star formation efficiency, without invoking any special adjustments. △ Less

Submitted 18 March, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Comments: 30 pages, 12 figures, ApJ, in press

A Milky Way-like barred spiral galaxy at a redshift of 3

Authors: Luca Costantin, Pablo G. Pérez-González, Yuchen Guo, Chiara Buttitta, Shardha Jogee, Micaela B. Bagley, Guillermo Barro, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Cristina Cabello, Enrico Maria Corsini, Jairo Méndez-Abreu, Alexander de la Vega, Kartheik G. Iyer, Laura Bisigello, Yingjie Cheng, Lorenzo Morelli, Pablo Arrabal Haro, Fernando Buitrago, M. C. Cooper, Avishai Dekel, Mark Dickinson, Steven L. Finkelstein, Mauro Giavalisco, Benne W. Holwerda , et al. (8 additional authors not shown)

Abstract: The majority of massive disk galaxies in the local Universe show a stellar barred structure in their central regions, including our Milky Way. Bars are supposed to develop in dynamically cold stellar disks at low redshift, as the strong gas turbulence typical of disk galaxies at high redshift suppresses or delays bar formation. Moreover, simulations predict bars to be almost absent beyond… ▽ More The majority of massive disk galaxies in the local Universe show a stellar barred structure in their central regions, including our Milky Way. Bars are supposed to develop in dynamically cold stellar disks at low redshift, as the strong gas turbulence typical of disk galaxies at high redshift suppresses or delays bar formation. Moreover, simulations predict bars to be almost absent beyond $z = 1.5$ in the progenitors of Milky Way-like galaxies. Here we report observations of ceers-2112, a barred spiral galaxy at redshift $z_{\rm phot} \sim 3$, which was already mature when the Universe was only 2 Gyr old. The stellar mass ($M_{\star} = 3.9 \times 10^9 M_{\odot}$) and barred morphology mean that ceers-2112 can be considered a progenitor of the Milky Way, in terms of both structure and mass-assembly history in the first 2 Gyr of the Universe, and was the closest in mass in the first 4 Gyr. We infer that baryons in galaxies could have already dominated over dark matter at $z \sim 3$, that high-redshift bars could form in approximately 400 Myr and that dynamically cold stellar disks could have been in place by redshift $z = 4-5$ (more than 12 Gyrs ago). △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: Published in Nature. Accepted on September 12

The Complete CEERS Early Universe Galaxy Sample: A Surprisingly Slow Evolution of the Space Density of Bright Galaxies at z ~ 8.5-14.5

Authors: Steven L. Finkelstein, Gene C. K. Leung, Micaela B. Bagley, Mark Dickinson, Henry C. Ferguson, Casey Papovich, Hollis B. Akins, Pablo Arrabal Haro, Romeel Dave, Avishai Dekel, Jeyhan S. Kartaltepe, Dale D. Kocevski, Anton M. Koekemoer, Norbert Pirzkal, Rachel S. Somerville, L. Y. Aaron Yung, Ricardo Amorin, Bren E. Backhaus, Peter Behroozi, Laura Bisigello, Volker Bromm, Caitlin M. Casey, Oscar A. Chavez Ortiz, Yingjie Cheng, Katherine Chworowsky , et al. (30 additional authors not shown)

Abstract: We present a sample of 88 candidate z~8.5-14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science (CEERS) survey. These data cover ~90 arcmin^2 (10 NIRCam pointings) in six broad-band and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than p… ▽ More We present a sample of 88 candidate z~8.5-14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science (CEERS) survey. These data cover ~90 arcmin^2 (10 NIRCam pointings) in six broad-band and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z~9, 11 and 14, and show that the space density of bright (M_UV=-20) galaxies changes only modestly from z~14 to z~9, compared to a steeper increase from z~8 to z~4. While our candidates are photometrically selected, spectroscopic followup has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z>10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models which naturally have enhanced star-formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic followup of a large sample of early galaxies can distinguish between these competing scenarios. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: Submitted to ApJL. Main paper is 33 pages, 13 figures, 5 tables. Two appendices with additional figures and tables

Evolution of the Size-Mass Relation of Star-forming Galaxies Since $z=5.5$ Revealed by CEERS

Authors: Ethan M. Ward, Alexander de la Vega, Bahram Mobasher, Elizabeth J. McGrath, Kartheik G. Iyer, Antonello Calabro, Luca Costantin, Mark Dickinson, Benne W. Holwerda, Marc Huertas-Company, Michaela Hirschmann, Ray A. Lucas, Viraj Pandya, Stephen M. Wilkins, L. Y. Aaron Yung, Pablo Arrabal Haro, Micaela B. Bagley, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Casey Papovich, Nor Pirzkal

Abstract: We combine deep imaging data from the CEERS early release JWST survey and HST imaging from CANDELS to examine the size-mass relation of star-forming galaxies and the morphology-quenching relation at stellar masses $\textrm{M}_{\star} \geq 10^{9.5} \ \textrm{M}_{\odot}$ over the redshift range $0.5 < z < 5.5$. In this study with a sample of 2,450 galaxies, we separate star-forming and quiescent gal… ▽ More We combine deep imaging data from the CEERS early release JWST survey and HST imaging from CANDELS to examine the size-mass relation of star-forming galaxies and the morphology-quenching relation at stellar masses $\textrm{M}_{\star} \geq 10^{9.5} \ \textrm{M}_{\odot}$ over the redshift range $0.5 < z < 5.5$. In this study with a sample of 2,450 galaxies, we separate star-forming and quiescent galaxies based on their star-formation activity and confirm that star-forming and quiescent galaxies have different morphologies out to $z=5.5$, extending the results of earlier studies out to higher redshifts. We find that star-forming and quiescent galaxies have typical Sérsic indices of $n\sim1.3$ and $n\sim4.3$, respectively. Focusing on star-forming galaxies, we find that the slope of the size-mass relation is nearly constant with redshift, as was found previously, but shows a modest increase at $z \sim 4.2$. The intercept in the size-mass relation declines out to $z=5.5$ at rates that are similar to what earlier studies found. The intrinsic scatter in the size-mass relation is relatively constant out to $z=5.5$. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 20 pages, 8 figures, 1 table, submitted to ApJ

Galaxies Going Bananas: Inferring the 3D Geometry of High-Redshift Galaxies with JWST-CEERS

Authors: Viraj Pandya, Haowen Zhang, Marc Huertas-Company, Kartheik G. Iyer, Elizabeth McGrath, Guillermo Barro, Steven L. Finkelstein, Martin Kuemmel, William G. Hartley, Henry C. Ferguson, Jeyhan S. Kartaltepe, Joel Primack, Avishai Dekel, Sandra M. Faber, David C. Koo, Greg L. Bryan, Rachel S. Somerville, Ricardo O. Amorin, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Emmanuel Bertin, Luca Costantin, Romeel Dave, Mark Dickinson , et al. (31 additional authors not shown)

Abstract: The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean e… ▽ More The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size and covariances with samples as small as $\sim50$ galaxies. We find high 3D ellipticities for all mass-redshift bins suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be $\sim1$ for $\log M_*/M_{\odot}=9.0-9.5$ dwarfs at $z>1$ (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a ``banana'' in the projected $b/a-\log a$ diagram with an excess of low $b/a$, large $\log a$ galaxies. The dwarf prolate fraction rises from $\sim25\%$ at $z=0.5-1.0$ to $\sim50-80\%$ at $z=3-8$. If these are disks, they cannot be axisymmetric but instead must be unusually oval (triaxial) unlike local circular disks. We simultaneously constrain the 3D size-mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar Sérsic indices ($n\sim1$), non-parametric morphological properties and specific star formation rates. Both tend to be visually classified as disks or irregular but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects and theoretical implications. △ Less

Submitted 15 January, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: Accepted version to appear in ApJ, main body is 36 pages of which ~half are full-page figures

$Λ$CDM not dead yet: massive high-z Balmer break galaxies are less common than previously reported

Authors: Guillaume Desprez, Nicholas S. Martis, Yoshihisa Asada, Marcin Sawicki, Chris J. Willott, Adam Muzzin, Roberto G. Abraham, Maruša Bradač, Gabe Brammer, Vicente Estrada-Carpenter, Kartheik G. Iyer, Jasleen Matharu, Lamiya Mowla, Gaël Noirot, Ghassan T. E. Sarrouh, Victoria Strait, Rachel Gledhill, Gregor Rihtaršič

Abstract: Early JWST observations that targeted so-called double-break sources (attributed to Lyman and Balmer breaks at $z>7$), reported a previously unknown population of very massive, evolved high-redshift galaxies. This surprising discovery led to a flurry of attempts to explain these objects' unexpected existence including invoking alternatives to the standard $Λ$CDM cosmological paradigm. To test thes… ▽ More Early JWST observations that targeted so-called double-break sources (attributed to Lyman and Balmer breaks at $z>7$), reported a previously unknown population of very massive, evolved high-redshift galaxies. This surprising discovery led to a flurry of attempts to explain these objects' unexpected existence including invoking alternatives to the standard $Λ$CDM cosmological paradigm. To test these early results, we adopted the same double-break candidate galaxy selection criteria to search for such objects in the JWST images of the CAnadian NIRISS Unbiased Cluster Survey (CANUCS), and found a sample of 19 sources over five independent CANUCS fields that cover a total effective area of $\sim60\,$arcmin$^2$ at $z\sim8$. However, (1) our SED fits do not yield exceptionally high stellar masses for our candidates, while (2) spectroscopy of five of the candidates shows that while all five are at high redshifts, their red colours are due to high-EW emission lines in star-forming galaxies rather than Balmer breaks in massive, evolved systems. Additionally, (3) field-to-field variance leads to differences of $\sim 1.5$ dex in the maximum stellar masses measured in the different fields, suggesting that the early single-field JWST observations may have suffered from cosmic variance and/or sample bias. Finally, (4) we show that the presence of even a single massive outlier can dominate conclusions from small samples such as those in early JWST observations. In conclusion, we find that the double-break sources in CANUCS are not sufficiently massive or numerous to warrant questioning the standard $Λ$CDM paradigm. △ Less

Submitted 19 April, 2024; v1 submitted 4 October, 2023; originally announced October 2023.

Comments: V3 Accepted for publication in MNRAS, 13 pages (+5 in Appendix), 7 figures (+4), 1 table (+2)

Star Formation at the Epoch of Reionization with CANUCS: The ages of stellar populations in MACS1149-JD1

Authors: Maruša Bradač, Victoria Strait, Lamiya Mowla, Kartheik G. Iyer, Gaël Noirot, Chris Willott, Gabe Brammer, Roberto Abraham, Yoshihisa Asada, Guillaume Desprez, Vince Estrada-Carpenter, Anishya Harshan, Nicholas S. Martis, Jasleen Matharu, Adam Muzzin, Gregor Rihtaršič, Ghassan T. E. Sarrouh, Marcin Sawicki

Abstract: We present measurements of stellar populations properties of a z = 9.1 gravitationally lensed galaxy MACS1149-JD1 using deep JWST NIRISS slitless spectroscopy as well as NIRISS and NIRCam imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). The galaxy is split into four components. Three magnified ($μ$ ~ 17) star-forming components are unresolved, giving intrinsic sizes of < 50pc. In… ▽ More We present measurements of stellar populations properties of a z = 9.1 gravitationally lensed galaxy MACS1149-JD1 using deep JWST NIRISS slitless spectroscopy as well as NIRISS and NIRCam imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). The galaxy is split into four components. Three magnified ($μ$ ~ 17) star-forming components are unresolved, giving intrinsic sizes of < 50pc. In addition, the underlying extended component contains the bulk of the stellar mass, formed the majority of its stars ~ 50Myr earlier than the other three components and is not the site of the most active star formation currently. The NIRISS and NIRCam resolved photometry does not confirm a strong Balmer break previously seen in Spitzer. The NIRISS grism spectrum has been extracted for the entire galaxy and shows a clear continuum and Lyman-break, with no Lyman-$α$ detected. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: This is the version of the article before peer review or editing, as submitted by an author to ApJLetters

Galaxy Morphology from $z\sim6$ through the eyes of JWST

Authors: M. Huertas-Company, K. G. Iyer, E. Angeloudi, M. B. Bagley, S. L. Finkelstein, J. Kartaltepe, R. Sarmiento, J. Vega-Ferrero, P. Arrabal Haro, P. Behroozi, F. Buitrago, Y. Cheng, L. Costantin, A. Dekel, M. Dickinson, D. Elbaz, N. A. Grogin, N. P. Hathi, B. W. Holwerda, A. M. Koekemoer, R. A. Lucas, C. Papovich, P. G. Pérez-González, N. Pirzkal, L-M. Seillé , et al. (4 additional authors not shown)

Abstract: We analyze the Near Infrared ($\sim0.8-1μ$m) rest-frame morphologies of galaxies with $\log M_*/M_\odot>9$ in the redshift range $0<z<6$, compare with previous HST-based results and release the first JWST-based morphological catalog of $\sim20,000$ galaxies in the CEERS survey. Galaxies are classified into four main broad classes -- spheroid, disk+spheroid, disk, and disturbed -- based on imaging… ▽ More We analyze the Near Infrared ($\sim0.8-1μ$m) rest-frame morphologies of galaxies with $\log M_*/M_\odot>9$ in the redshift range $0<z<6$, compare with previous HST-based results and release the first JWST-based morphological catalog of $\sim20,000$ galaxies in the CEERS survey. Galaxies are classified into four main broad classes -- spheroid, disk+spheroid, disk, and disturbed -- based on imaging with four filters -- $F150W$, $F200W$, $F356W$, and $F444W$ -- using Convolutional Neural Networks trained on HST/WFC3 labeled images and domain-adapted to JWST/NIRCam. We find that $\sim90\%$ and $\sim75\%$ of galaxies at $z<3$ have the same early/late and regular/irregular classification, respectively, in JWST and HST imaging when considering similar wavelengths. For small (large) and faint objects, JWST-based classifications tend to systematically present less bulge-dominated systems (peculiar galaxies) than HST-based ones, but the impact on the reported evolution of morphological fractions is less than $\sim10\%$. Using JWST-based morphologies at the same rest-frame wavelength ($\sim0.8-1μ$m), we confirm an increase in peculiar galaxies and a decrease in bulge-dominated galaxies with redshift, as reported in previous HST-based works, suggesting that the stellar mass distribution, in addition to light distribution, is more disturbed in the early universe. However, we find that undisturbed disk-like systems already dominate the high-mass end of the late-type galaxy population ($\log M_*/M_\odot>10.5$) at $z\sim5$, and bulge-dominated galaxies also exist at these early epochs, confirming a rich and evolved morphological diversity of galaxies $\sim1$ Gyr after the Big Bang. Finally, we find that the morphology-quenching relation is already in place for massive galaxies at $z>3$, with massive quiescent galaxies ($\log M_*/M_\odot>10.5$) being predominantly bulge-dominated. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: Submitted to A&A, comments welcome

An extremely compact, low-mass post-starburst galaxy at $z=5.2$

Authors: Victoria Strait, Gabriel Brammer, Adam Muzzin, Guillaume Dezprez, Yoshihisi Asada, Roberto Abraham, Maruša Bradač, Kartheik G. Iyer, Nicholas Martis, Lamiya Mowla, Gaël Noirot, Ghassan Sarrouh, Marcin Sawicki, Chris Willott, Katriona Gould, Tess Grindlay, Jasleen Matharu, Gregor Rihtaršič

Abstract: We report the discovery of a low-mass $z=5.200\pm 0.002$ galaxy that is in the process of ceasing its star formation. The galaxy, MACS0417-z5PSB, is multiply imaged with magnification factors $\sim40$ by the galaxy cluster MACS J0417.5-1154, observed as part of the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). Using observations of MACS0417-z5PSB with a JWST/NIRSpec Prism spectrum and NIRCam i… ▽ More We report the discovery of a low-mass $z=5.200\pm 0.002$ galaxy that is in the process of ceasing its star formation. The galaxy, MACS0417-z5PSB, is multiply imaged with magnification factors $\sim40$ by the galaxy cluster MACS J0417.5-1154, observed as part of the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). Using observations of MACS0417-z5PSB with a JWST/NIRSpec Prism spectrum and NIRCam imaging, we investigate the mechanism responsible for the cessation of star formation of the galaxy, and speculate about possibilities for its future. Using spectrophotometric fitting, we find a remarkably low stellar mass of $\rm{M_*}=4.3\pm^{0.9}_{0.8} \times 10^{7} \rm{M_{\odot}}$, less than 1% of the characteristic stellar mass at $z\sim5$. We measure a de-lensed rest-UV half-light radius in the source plane of $30\pm^{7}_{5}$ pc, and measure a star formation rate from H$α$ of $0.14\pm^{0.17}_{0.12}$ $\rm{M_{\odot}/yr}$. We find that under the assumption of a double power law star formation history, MACS0417-z5PSB has seen a recent rise in star formation, peaking $\sim10-30$ Myr ago and declining precipitously since then. Together, these measurements reveal a low-mass, extremely compact galaxy which is in the process of ceasing star formation. We investigate the possibilities of mechanisms that have led to the cessation of star formation in MACS0417-z5PSB, considering stellar and AGN feedback, and environmental processes. We can likely rule out an AGN and most environmental processes, but leave open the possibility that MACS0417-z5PSB could be a star forming galaxy in the lull of a bursty star formation history. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 9 pages, 4 figures, 1 table

CEERS Spectroscopic Confirmation of NIRCam-Selected z > 8 Galaxy Candidates with JWST/NIRSpec: Initial Characterization of their Properties

Authors: Seiji Fujimoto, Pablo Arrabal Haro, Mark Dickinson, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Rebecca L. Larson, Denis Burgarella, Micaela B. Bagley, Peter Behroozi, Katherine Chworowsky, Michaela Hirschmann, Jonathan R. Trump, Stephen M. Wilkins, L. Y. Aaron Yung, Anton M. Koekemoer, Casey Papovich, Nor Pirzkal, Henry C. Ferguson, Adriano Fontana, Norman A. Grogin, Andrea Grazian, Lisa J. Kewley, Dale D. Kocevski, Jennifer M. Lotz, Laura Pentericci , et al. (16 additional authors not shown)

Abstract: We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of $z\simeq9-13$ and $M_{\rm\,UV} \in[-21,-18]$ newly identified in NIRCam images in the Cosmic Evolution Early Release Science (CEERS) Survey. We confirm emission line redshifts for 7 galaxies at $z=7.762-8.998$ using spectra at $\sim1-5μ$m either with the NIRSpec prism or its three medium resolution gratings… ▽ More We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of $z\simeq9-13$ and $M_{\rm\,UV} \in[-21,-18]$ newly identified in NIRCam images in the Cosmic Evolution Early Release Science (CEERS) Survey. We confirm emission line redshifts for 7 galaxies at $z=7.762-8.998$ using spectra at $\sim1-5μ$m either with the NIRSpec prism or its three medium resolution gratings. For $z\simeq9$ photometric candidates, we achieve a high confirmation rate of $\simeq$90\%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at $z>10$, where the strong [OIII] and H$β$ lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Lyman-$α$ continuum break is not detected with the current sensitivity. Compared with HST-selected bright galaxies ($M_{\rm\,UV}\simeq-22$) that are similarly spectroscopically confirmed at $z\gtrsim8$, these NIRCam-selected galaxies are characterized by lower star formation rates (SFR$\simeq4\,M_{\odot}$~yr$^{-1}$) and lower stellar masses ($\simeq10^{8}\,M_{\odot}$), but with higher [OIII]+H$β$ equivalent widths ($\simeq$1100$Å$), and elevated production efficiency of ionizing photons ($\log(ξ_{\rm\,ion}/{\rm\,Hz\,erg}^{-1})\simeq25.8$) induced by young stellar populations ($<10$~Myrs) accounting for $\simeq20\%$ of the galaxy mass, highlighting the key contribution of faint galaxies to cosmic reionization. Taking advantage of the homogeneous selection and sensitivity, we also investigate metallicity and ISM conditions with empirical calibrations using the [OIII]/H$β$ ratio. We find that galaxies at $z\sim8-9$ have higher SFRs and lower metallicities than galaxies at similar stellar masses at $z\sim2-6$, which is generally consistent with the current galaxy formation and evolution models. △ Less

Submitted 26 May, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: 24 pages, 12 figures, 2 tables. ApJL Focus Issue in press

Spatially Resolved Stellar Populations of $0.3<z<6.0$ Galaxies in WHL0137-08 and MACS0647+70 Clusters as Revealed by JWST: How do Galaxies Grow and Quench Over Cosmic Time?

Authors: Abdurro'uf, Dan Coe, Intae Jung, Henry C. Ferguson, Gabriel Brammer, Kartheik G. Iyer, Larry D. Bradley, Pratika Dayal, Rogier A. Windhorst, Adi Zitrin, Ashish Kumar Meena, Masamune Oguri, Jose M. Diego, Vasily Kokorev, Paola Dimauro, Angela Adamo, Christopher J. Conselice, Brian Welch, Eros Vanzella, Tiger Yu-Yang Hsiao, Xinfeng Xu, Namrata Roy, Celia R. Mulcahey

Abstract: We study the spatially resolved stellar populations of 444 galaxies at $0.3<z<6.0$ in two clusters (WHL0137-08 and MACS0647+70) and a blank field, combining imaging data from HST and JWST to perform spatially resolved spectral energy distribution (SED) modeling using pixedfit. The high spatial resolution of the imaging data combined with magnification from gravitational lensing in the cluster fiel… ▽ More We study the spatially resolved stellar populations of 444 galaxies at $0.3<z<6.0$ in two clusters (WHL0137-08 and MACS0647+70) and a blank field, combining imaging data from HST and JWST to perform spatially resolved spectral energy distribution (SED) modeling using pixedfit. The high spatial resolution of the imaging data combined with magnification from gravitational lensing in the cluster fields allows us to resolve some galaxies to sub-kpc scales (for 109 of our galaxies). At redshifts around cosmic noon and higher ($2.5\lesssim z\lesssim 6.0$), we find mass doubling times to be independent of radius, inferred from flat specific star formation rate (sSFR) radial profiles and similarities between the half-mass and half-SFR radii. At lower redshifts ($1.5\lesssim z\lesssim 2.5$), a significant fraction of our star-forming galaxies show evidence for nuclear starbursts, inferred from centrally elevated sSFR, and a much smaller half-SFR radius compared to the half-mass radius. At later epochs, we find more galaxies suppress star formation in their center but are still actively forming stars in the disk. Overall, these trends point toward a picture of inside-out galaxy growth consistent with theoretical models and simulations. We also observe a tight relationship between the central mass surface density and global stellar mass with $\sim 0.38$ dex scatter. Our analysis demonstrates the potential of spatially resolved SED analysis with JWST data. Future analysis with larger samples will be able to further explore the assembly of galaxy mass and the growth of their structures △ Less

Submitted 7 February, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 31 pages, 18 figures, accepted for publication in ApJ. Some examples and tutorials of spatially resolved SED analysis will be available at https://github.com/aabdurrouf/JWST-HST_resolvedSEDfits

The Art of Measuring Physical Parameters in Galaxies: A Critical Assessment of Spectral Energy Distribution Fitting Techniques

Authors: Camilla Pacifici, Kartheik G. Iyer, Bahram Mobasher, Elisabete da Cunha, Viviana Acquaviva, Denis Burgarella, Gabriela Calistro Rivera, Adam C. Carnall, Yu-Yen Chang, Nima Chartab, Kevin C. Cooke, Ciaran Fairhurst, Jeyhan Kartaltepe, Joel Leja, Katarzyna Malek, Brett Salmon, Marianna Torelli, Alba Vidal-Garcia, Mederic Boquien, Gabriel G. Brammer, Michael J. I. Brown, Peter L. Capak, Jacopo Chevallard, Chiara Circosta, Darren Croton , et al. (30 additional authors not shown)

Abstract: The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have signifi… ▽ More The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply fourteen of the most commonly used SED fitting codes on samples from the CANDELS photometric catalogs at z~1 and z~3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust, and AGN models) on the derived stellar masses, SFRs, and A_V values. We then assess the difference among the codes on the SFR-stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (~0.1dex), SFR (~0.3dex), and dust attenuation (~0.3mag). Finally, we present some resources summarizing best practices in SED fitting. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 25 pages, 11 figures. Accepted for publication in ApJ

AGN Feedback in SDSS-IV MaNGA: AGNs have Suppressed Central Star Formation Rates

Authors: Caleb Lammers, Kartheik G. Iyer, Hector Ibarra-Medel, Camilla Pacifici, Sebastián F. Sánchez, Sandro Tacchella, Joanna Woo

Abstract: Despite the importance of feedback from active galactic nuclei (AGNs) in models of galaxy evolution, observational constraints on the influence of AGN feedback on star formation remain weak. To this end, we have compared the star formation trends of 279 low-redshift AGN galaxies with 558 inactive control galaxies using integral field unit spectroscopy from the SDSS-IV MaNGA survey. With a Gaussian… ▽ More Despite the importance of feedback from active galactic nuclei (AGNs) in models of galaxy evolution, observational constraints on the influence of AGN feedback on star formation remain weak. To this end, we have compared the star formation trends of 279 low-redshift AGN galaxies with 558 inactive control galaxies using integral field unit spectroscopy from the SDSS-IV MaNGA survey. With a Gaussian process-based methodology, we reconstruct nonparametric star formation histories in spatially resolved spaxels covering the face of each galaxy. Based on galaxy-wide star formation rates (SFRs) alone, we find no obvious signatures of AGN feedback. However, the AGN galaxies have significantly suppressed central (kiloparsec-scale) SFRs, lying up to a factor of $2$ below those of the control galaxies, providing direct observational evidence of AGN feedback suppressing star formation. The suppression of central SFRs in the AGN galaxies began in the central regions ${\sim} 6$ Gyr ago (redshift $z {\sim} 0.7$), taking place over a few gigayears. A small subset of the AGN galaxies were rapidly driven to quiescence shortly before being observed (in the last $500$ Myr), potentially indicating instances of AGN-driven feedback. More frequently, however, star formation continues in the AGN galaxies, with suppression primarily in the central regions. This is suggestive of a picture in which integrated (Gyr-timescale) AGN feedback can significantly affect central star formation, but may be inefficient in driving galaxy-wide quenching in low-redshift galaxies, instead leaving them in the green valley. △ Less

Submitted 28 June, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 22 pages, 15 figures. Accepted for publication in ApJ

A spatially resolved analysis of star-formation burstiness by comparing UV and H$α$ in galaxies at z$\sim$1 with UVCANDELS

Authors: Vihang Mehta, Harry I. Teplitz, Claudia Scarlata, Xin Wang, Anahita Alavi, James Colbert, Marc Rafelski, Norman Grogin, Anton Koekemoer, Laura Prichard, Rogier Windhorst, Justin M. Barber, Christopher J. Conselice, Y. Sophia Dai, Jonathan P. Gardner, Eric Gawiser, Yicheng Guo, Nimish Hathi, Pablo Arrabal Haro, Matthew Hayes, Kartheik G. Iyer, Rolf A. Jansen, Zhiyuan Ji, Peter Kurczynski, Maxwell Kuschel , et al. (8 additional authors not shown)

Abstract: The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides HST/UVIS F275W imaging for four CANDELS fields. We combine this UV imaging with existing HST/near-IR grism spectroscopy from 3D-HST$+$AGHAST to directly compare the resolved rest-frame UV and H$α$ emission for a sample of 979 galaxies at $0.7<z<1.5$ spanning a range in… ▽ More The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides HST/UVIS F275W imaging for four CANDELS fields. We combine this UV imaging with existing HST/near-IR grism spectroscopy from 3D-HST$+$AGHAST to directly compare the resolved rest-frame UV and H$α$ emission for a sample of 979 galaxies at $0.7<z<1.5$ spanning a range in stellar mass of $10^{8-11.5}$ M$_\odot$. Using a stacking analysis, we perform a resolved comparison between homogenized maps of rest-UV and H$α$ to compute the average UV-to-H$α$ luminosity ratio (an indicator of burstiness in star-formation) as a function of galactocentric radius. We find that galaxies below stellar mass of $\sim$10$^{9.5}$ M$_\odot$, at all radii, have a UV-to-H$α$ ratio higher than the equilibrium value expected from constant star-formation, indicating a significant contribution from bursty star-formation. Even for galaxies with stellar mass $\gtrsim$10$^{9.5}$ M$_\odot$, the UV-to-H$α$ ratio is elevated towards in their outskirts ($R/R_{eff}>1.5$), suggesting that bursty star-formation is likely prevalent in the outskirts of even the most massive galaxies but is likely over-shadowed by their brighter cores. Furthermore, we present the UV-to-H$α$ ratio as a function of galaxy surface brightness, a proxy for stellar mass surface density, and find that regions below $\sim$10$^{7.5}$ M$_\odot$ kpc$^{-2}$ are consistent with bursty star-formation, regardless of their galaxy stellar mass, potentially suggesting that local star-formation is independent of global galaxy properties at the smallest scales. Lastly, we find galaxies at $z>1.1$ to have bursty star-formation regardless of radius or surface brightness. △ Less

Submitted 15 June, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 22 pages, 10 figures; accepted to ApJ

Star Formation Variability as a Probe for the Baryon Cycle within Galaxies

Authors: Eun-jin Shin, Sandro Tacchella, Ji-hoon Kim, Kartheik G. Iyer, Vadim A. Semenov

Abstract: We investigate the connection of the regulation of star formation and the cycling of baryons within and in and out of galaxies. We use idealized numerical simulations of Milky Way-mass galaxies, in which we systemically vary the galaxy morphology (bulge-to-total mass ratio) and stellar feedback strength (total eight setups with 80 simulations). By following individual gas parcels through the disk,… ▽ More We investigate the connection of the regulation of star formation and the cycling of baryons within and in and out of galaxies. We use idealized numerical simulations of Milky Way-mass galaxies, in which we systemically vary the galaxy morphology (bulge-to-total mass ratio) and stellar feedback strength (total eight setups with 80 simulations). By following individual gas parcels through the disk, spiral arms, and massive star-forming clumps, we quantify how gas moves and oscillates through the different phases of the interstellar medium (ISM) and forms stars. We show that the residence time of gas in the dense ISM phase ($τ_{\rm SF}$), the nature of spiral arms (strength, number), and the clump properties (number, mass function, and young star fraction) depend on both the galaxy morphology and stellar feedback. Based on these results, we quantify signatures of the baryon cycle within galaxies using the temporal and spatial power spectrum density (PSD) of the star formation history (SFH). Stronger stellar feedback leads to more bursty star formation while the correlation timescale of the SFH is longer, because stronger feedback dissolves the dense, star-forming ISM phase, leading to a more homogeneous ISM and a decrease in $τ_{\rm SF}$. The bulge strength has a similar effect: the deep gravitational potential in a bulge-dominant galaxy imposes a strong shear force that effectively breaks apart gas clumps in the ISM; this subsequently inhibits the fragmentation of cool gas and therefore the star formation in the disk, leading to a decrease in the spatial power on scales of $\sim$ 1 kpc. We conclude that measurements of the temporal and spatial PSD of the SFH can provide constraints on the baryon cycle and the star formation process. △ Less

Submitted 3 November, 2022; originally announced November 2022.

Comments: 22 pages including 17 figures. Submitted to ApJ

Journal ref: 2023 ApJ 947 61

CEERS Key Paper III: The Diversity of Galaxy Structure and Morphology at z=3-9 with JWST

Authors: Jeyhan S. Kartaltepe, Caitlin Rose, Brittany N. Vanderhoof, Elizabeth J. McGrath, Luca Costantin, Isabella G. Cox, L. Y. Aaron Yung, Dale D. Kocevski, Stijn Wuyts, Henry C. Ferguson Brett H. Andrews, Micaela B. Bagley, Steven L. Finkelstein, Ricardo O. Amorin, Pablo Arrabal Haro, Bren E. Backhaus, Peter Behroozi, Laura Bisigello, Antonello Calabro, Caitlin M. Casey, Rosemary T. Coogan, Darren Croton, Alexander de la Vega, Mark Dickinson, M. C. Cooper, Adriano Fontana , et al. (36 additional authors not shown)

Abstract: We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z=3-9 using early JWST CEERS NIRCam observations. Our sample consists of 850 galaxies at z>3 detected in both CANDELS HST imaging and JWST CEERS NIRCam images to enable a comparison of HST and JWST morphologies. Our team conducted a set of visual classifications, wit… ▽ More We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z=3-9 using early JWST CEERS NIRCam observations. Our sample consists of 850 galaxies at z>3 detected in both CANDELS HST imaging and JWST CEERS NIRCam images to enable a comparison of HST and JWST morphologies. Our team conducted a set of visual classifications, with each galaxy in the sample classified by three different individuals. We also measure quantitative morphologies using the publicly available codes across all seven NIRCam filters. Using these measurements, we present the fraction of galaxies of each morphological type as a function of redshift. Overall, we find that galaxies at z>3 have a wide diversity of morphologies. Galaxies with disks make up a total of 60\% of galaxies at z=3 and this fraction drops to ~30% at z=6-9, while galaxies with spheroids make up ~30-40% across the whole redshift range and pure spheroids with no evidence for disks or irregular features make up ~20%. The fraction of galaxies with irregular features is roughly constant at all redshifts (~40-50%), while those that are purely irregular increases from ~12% to ~20% at z>4.5. We note that these are apparent fractions as many selection effects impact the visibility of morphological features at high redshift. The distributions of Sérsic index, size, and axis ratios show significant differences between the morphological groups. Spheroid Only galaxies have a higher Sérsic index, smaller size, and higher axis ratio than Disk/Irregular galaxies. Across all redshifts, smaller spheroid and disk galaxies tend to be rounder. Overall, these trends suggest that galaxies with established disks and spheroids exist across the full redshift range of this study and further work with large samples at higher redshift is needed to quantify when these features first formed. △ Less

Submitted 13 January, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: Accepted for publication to ApJL, 24 pages, 14 figures

CEERS Key Paper III: The Resolved Host Properties of AGN at 3 < z < 5 with JWST

Authors: Dale D. Kocevski, Guillermo Barro, Elizabeth J. McGrath, Steven L. Finkelstein, Micaela B. Bagley, Henry C. Ferguson, Shardha Jogee, Guang Yang, Mark Dickinson, Nimish P. Hathi, Bren E. Backhaus, Eric F. Bell, Laura Bisigello, Véronique Buat, Denis Burgarella, Caitlin M. Casey, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Darren Croton, Emanuele Daddi, Adriano Fontana, Seiji Fujimoto, Jonathan P. Gardner, Eric Gawiser , et al. (34 additional authors not shown)

Abstract: We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshif… ▽ More We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host galaxy by fitting its spectral energy distribution using a combination of galaxy and AGN templates. The AGN hosts have an average stellar mass of ${\rm log}(M_{*}/{\rm M_{\odot}} )= 11.0$, making them among the most massive galaxies detected at this redshift range in the current CEERS pointings, even after accounting for nuclear light from the AGN. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk and one host is dominated by point-like emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or post-starburst stellar population. The presence of AGN in passively evolving galaxies at $z>3$ is significant because a rapid feedback mechanism is required in most semi-analytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings are in general agreement with this picture and show that AGN can continue to inject energy into these systems after their star formation is curtailed, possibly helping to maintain their quiescent state. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: Submitted to ApJ; 12 Pages, 6 Figures

Stochastic Modelling of Star Formation Histories III. Constraints from Physically-Motivated Gaussian Processes

Authors: Kartheik G. Iyer, Joshua S. Speagle, Neven Caplar, John C. Forbes, Eric Gawiser, Joel Leja, Sandro Tacchella

Abstract: Galaxy formation and evolution involves a variety of effectively stochastic processes that operate over different timescales. The Extended Regulator model provides an analytic framework for the resulting variability (or `burstiness') in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow,… ▽ More Galaxy formation and evolution involves a variety of effectively stochastic processes that operate over different timescales. The Extended Regulator model provides an analytic framework for the resulting variability (or `burstiness') in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow, equilibrium, and dynamical processes influencing GMC creation and destruction using the power spectral density (PSD) formalism. We use the connection between the PSD and auto-covariance function (ACF) for general stochastic processes to reformulate this model as an auto-covariance function, which we use to model variability in galaxy star formation histories (SFHs) using physically-motivated Gaussian Processes in log SFR space. Using stellar population synthesis models, we then explore how changes in model stochasticity can affect spectral signatures across galaxy populations with properties similar to the Milky Way and present-day dwarfs as well as at higher redshifts. We find that, even at fixed scatter, perturbations to the stochasticity model (changing timescales vs overall variability) leave unique spectral signatures across both idealized and more realistic galaxy populations. Distributions of spectral features including H$α$ and UV-based SFR indicators, H$δ$ and Ca-H,K absorption line strengths, D$_n$(4000) and broadband colors provide testable predictions for galaxy populations from present and upcoming surveys with Hubble, Webb \& Roman. The Gaussian process SFH framework provides a fast, flexible implementation of physical covariance models for the next generation of SED modeling tools. Code to reproduce our results can be found at https://github.com/kartheikiyer/GP-SFH △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 33 pages, 18 figures. Submitted to ApJ, comments are welcome! Code available at https://github.com/kartheikiyer/GP-SFH

The Sparkler: Evolved High-Redshift Globular Clusters Captured by JWST

Authors: Lamiya A. Mowla, Kartheik G. Iyer, Guillaume Desprez, Vicente Estrada-Carpenter, Nicholas S. Martis, Gaël Noirot, Ghassan T. Sarrouh, Victoria Strait, Yoshihisa Asada, Roberto G. Abraham, Gabriel Brammer, Marcin Sawicki, Chris J. Willott, Marusa Bradac, René Doyon, Kate Gould, Adam Muzzin, Camilla Pacifici, Swara Ravindranath, Johannes Zabl

Abstract: Using data from JWST, we analyze the compact sources ("sparkles") located around a remarkable $z_{\rm spec}=1.378$ galaxy (the "Sparkler") that is strongly gravitationally lensed by the $z=0.39$ galaxy cluster SMACS J0723.3-7327. Several of these compact sources can be cross-identified in multiple images, making it clear that they are associated with the host galaxy. Combining data from JWST's {\e… ▽ More Using data from JWST, we analyze the compact sources ("sparkles") located around a remarkable $z_{\rm spec}=1.378$ galaxy (the "Sparkler") that is strongly gravitationally lensed by the $z=0.39$ galaxy cluster SMACS J0723.3-7327. Several of these compact sources can be cross-identified in multiple images, making it clear that they are associated with the host galaxy. Combining data from JWST's {\em Near-Infrared Camera} (NIRCam) with archival data from the {\em Hubble Space Telescope} (HST), we perform 0.4-4.4$μ$m photometry on these objects, finding several of them to be very red and consistent with the colors of quenched, old stellar systems. Morphological fits confirm that these red sources are spatially unresolved even in strongly magnified JWST/NIRCam images, while JWST/NIRISS spectra show [OIII]5007 emission in the body of the Sparkler but no indication of star formation in the red compact sparkles. The most natural interpretation of these compact red companions to the Sparkler is that they are evolved globular clusters seen at $z=1.378$. Applying \textsc{Dense Basis} SED-fitting to the sample, we infer formation redshifts of $z_{form} \sim 7-11$ for these globular cluster candidates, corresponding to ages of $\sim 3.9-4.1$ Gyr at the epoch of observation and a formation time just $\sim$0.5~Gyr after the Big Bang. If confirmed with additional spectroscopy, these red, compact "sparkles" represent the first evolved globular clusters found at high redshift, could be amongst the earliest observed objects to have quenched their star formation in the Universe, and may open a new window into understanding globular cluster formation. Data and code to reproduce our results will be made available at \faGithub\href{https://niriss.github.io/sparkler.html}{http://canucs-jwst.com/sparkler.html}. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: Submitted to ApJL. Comments are welcome. Data and code to reproduce our results will be made available at niriss.github.io/sparkler.html

Dusty Starbursts Masquerading as Ultra-high Redshift Galaxies in JWST CEERS Observations

Authors: Jorge A. Zavala, Veronique Buat, Caitlin M. Casey, Denis Burgarella, Steven L. Finkelstein, Micaela B. Bagley, Laure Ciesla, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Maximilien Franco, E. F. Jim'enez-Andrade, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Aurélien Le Bail, E. J. Murphy, Casey Papovich, Sandro Tacchella, Stephen M. Wilkins, Itziar Aretxaga, Peter Behroozi, Jaclyn B. Champagne, Adriano Fontana, Mauro Giavalisco, Andrea Grazian , et al. (99 additional authors not shown)

Abstract: Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may als… ▽ More Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations. △ Less

Submitted 30 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: Published in The Astrophysical Journal Letters (updated to match the published version)

Journal ref: ApJL 943 L9 (2023)

A Long Time Ago in a Galaxy Far, Far Away: A Candidate z ~ 12 Galaxy in Early JWST CEERS Imaging

Authors: Steven L. Finkelstein, Micaela B. Bagley, Pablo Arrabal Haro, Mark Dickinson, Henry C. Ferguson, Jeyhan S. Kartaltepe, Casey Papovich, Denis Burgarella, Dale D. Kocevski, Marc Huertas-Company, Kartheik G. Iyer, Rebecca L. Larson, Pablo G. Pérez-González, Caitlin Rose, Sandro Tacchella, Stephen M. Wilkins, Katherine Chworowsky, Aubrey Medrano, Alexa M. Morales, Rachel S. Somerville, L. Y. Aaron Yung, Adriano Fontana, Mauro Giavalisco, Andrea Grazian, Norman A. Grogin , et al. (95 additional authors not shown)

Abstract: We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging f… ▽ More We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang. △ Less

Submitted 7 September, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

3D-DASH: The Widest Near-Infrared Hubble Space Telescope Survey

Authors: Lamiya A. Mowla, Sam E. Cutler, Gabriel B. Brammer, Ivelina G. Momcheva, Katherine E. Whitaker, Pieter G. van Dokkum, Rachel S. Bezanson, Natascha M. Forster Schreiber, Marijn Franx, Kartheik G. Iyer, Danilo Marchesini, Adam Muzzin, Erica J. Nelson, Rosalind E. Skelton, Gregory F. Snyder, David A. Wake, Stijn Wuyts, Arjen van der Wel

Abstract: The 3D-Drift And SHift (3D-DASH) program is a \textit{Hubble Space Telescope} WFC3 F160W imaging and G141 grism survey of the equatorial COSMOS field. 3D-DASH extends the legacy of HST near-infrared imaging and spectroscopy to degree-scale swaths of the sky, enabling the identification and study of distant galaxies ($z>2$) that are rare or in short-lived phases of galaxy evolution at rest-frame op… ▽ More The 3D-Drift And SHift (3D-DASH) program is a \textit{Hubble Space Telescope} WFC3 F160W imaging and G141 grism survey of the equatorial COSMOS field. 3D-DASH extends the legacy of HST near-infrared imaging and spectroscopy to degree-scale swaths of the sky, enabling the identification and study of distant galaxies ($z>2$) that are rare or in short-lived phases of galaxy evolution at rest-frame optical wavelengths. Furthermore, when combined with existing ACS/F814W imaging, the program facilitates spatially-resolved studies of the stellar populations and dust content of intermediate-redshift ($0.5<z<2$) galaxies. Here we present the reduced F160W imaging mosaic available to the community. Observed with the efficient DASH technique, the mosaic comprises 1256 individual WFC3 pointings, corresponding to an area of 1.35 deg$^2$ (1.43 deg$^2$ in 1912 when including archival data). The median $5σ$ point-source limit in $H_{160}$ is 24.74 mag. We also provide tools to determine the local point spread function (PSF), create cutouts, and explore the image at any location within the 3D-DASH footprint. 3D-DASH is the widest \textit{HST}/WFC3 imaging survey in the F160W filter to date, increasing the existing extragalactic survey area in the near-infrared at HST resolution by an order of magnitude. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: The 3D-DASH mosaic, PSF and cutout generator tool, and the image explorer are available at archive.stsci.edu/hlsp/3d-dash and at www.lamiyamowla.com/3d-dash. Accepted for publication in Astrophysical Journal. arXiv admin note: text overlap with arXiv:1808.04379

The Near Infrared Imager and Slitless Spectrograph for the James Webb Space Telescope -- II. Wide Field Slitless Spectroscopy

Authors: Chris J. Willott, René Doyon, Loic Albert, Gabriel B. Brammer, William V. Dixon, Koraljka Muzic, Swara Ravindranath, Aleks Scholz, Roberto Abraham, Étienne Artigau, Maruša Bradač, Paul Goudfrooij, John B. Hutchings, Kartheik G. Iyer, Ray Jayawardhana, Stephanie LaMassa, Nicholas Martis, Michael R. Meyer, Takahiro Morishita, Lamiya Mowla, Adam Muzzin, Gaël Noirot, Camilla Pacifici, Neil Rowlands, Ghassan Sarrouh , et al. (4 additional authors not shown)

Abstract: We present the wide field slitless spectroscopy mode of the NIRISS instrument on the James Webb Space Telescope. This mode employs two orthogonal low-resolution (resolving power $\approx 150$) grisms in combination with a set of six blocking filters in the wavelength range 0.8 to $2.3\,μ$m to provide a spectrum of almost every source across the field-of-view. When combined with the low background,… ▽ More We present the wide field slitless spectroscopy mode of the NIRISS instrument on the James Webb Space Telescope. This mode employs two orthogonal low-resolution (resolving power $\approx 150$) grisms in combination with a set of six blocking filters in the wavelength range 0.8 to $2.3\,μ$m to provide a spectrum of almost every source across the field-of-view. When combined with the low background, high sensitivity and high spatial resolution afforded by the telescope, this mode will enable unprecedented studies of the structure and evolution of distant galaxies. We describe the performance of the as-built hardware relevant to this mode and expected imaging and spectroscopic sensitivity. We discuss operational and calibration procedures to obtain the highest quality data. As examples of the observing mode usage, we present details of two planned Guaranteed Time Observations programs: The Canadian NIRISS Unbiased Cluster Survey (CANUCS) and The NIRISS Survey for Young Brown Dwarfs and Rogue Planets. △ Less

Submitted 3 February, 2022; originally announced February 2022.

Comments: 19 pages, 12 figures, accepted for publication in PASP

IQ Collaboratory III: The Empirical Dust Attenuation Framework -- Taking Hydrodynamical Simulations with a Grain of Dust

Authors: ChangHoon Hahn, Tjitske K. Starkenburg, Daniel Anglés-Alcázar, Ena Choi, Romeel Davé, Claire Dickey, Kartheik G. Iyer, Ariyeh H. Maller, Rachel S. Somerville, Jeremy L. Tinker, L. Y. Aaron Yung

Abstract: We present the Empirical Dust Attenuation (EDA) framework -- a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color-magnitud… ▽ More We present the Empirical Dust Attenuation (EDA) framework -- a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color-magnitude relations, $(g-r) - M_r$ and $(FUV-NUV)-M_r$, of the simulations to a $M_r < -20$ and UV complete SDSS galaxy sample using likelihood-free inference. Without dust, none of the simulations match observations, as expected. With the EDA, however, we can reproduce the observed color-magnitude with all three simulations. Furthermore, the attenuation curves predicted by our dust prescription are in good agreement with the observed attenuation-slope relations and attenuation curves of star-forming galaxies. However, the EDA does not predict star-forming galaxies with low $A_V$ since simulated star-forming galaxies are intrinsically much brighter than observations. Additionally, the EDA provides, for the first time, predictions on the attenuation curves of quiescent galaxies, which are challenging to measure observationally. Simulated quiescent galaxies require shallower attenuation curves with lower amplitude than star-forming galaxies. The EDA, combined with forward modeling, provides an effective approach for shedding light on dust in galaxies and probing hydrodynamical simulations. This work also illustrates a major limitation in comparing galaxy formation models: by adjusting dust attenuation, simulations that predict significantly different galaxy populations can reproduce the same UV and optical observations. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: 26 pages, 15 figures

SED Analysis of 13 Spectroscopically Confirmed Galaxies at z$\simeq$6 to Constrain UV-Slope, Model Dust Attenuation and Escape Fractions

Authors: Junehyoung Jeon, Rogier Windhorst, Seth H. Cohen, Rolf A. Jansen, Brent M. Smith, Timothy Carleton, Eiichi Egami, Kristian Finlator, Linhua Jiang, Kartheik G. Iyer, Vihang Mehta

Abstract: The reionization of the hydrogen in the Universe is thought to have completed by redshift $z\simeq5.5-6$. To probe this era, galaxy observations in the Subaru Deep Field (SDF) have identified more than 100 galaxies at $z\simeq6$, many spectroscopically confirmed through follow-up observations. We model the spectral energy distributions (SEDs) of 13 SDF galaxies with the CIGALE and Dense Basis code… ▽ More The reionization of the hydrogen in the Universe is thought to have completed by redshift $z\simeq5.5-6$. To probe this era, galaxy observations in the Subaru Deep Field (SDF) have identified more than 100 galaxies at $z\simeq6$, many spectroscopically confirmed through follow-up observations. We model the spectral energy distributions (SEDs) of 13 SDF galaxies with the CIGALE and Dense Basis codes using available optical/IR data. Modeling deep IR photometry has the potential to constrain the galaxy's Lyman continuum (LyC) escape fraction (\fesc). We use the modeled nebular emission lines and find that the implied escape fractions ranges from \textbf{0 to 0.8 with a median of $\sim$0.35 for Dense Basis and $\sim$0.55 for CIGALE.} Significant uncertainties in the data exist, so that fitting results in a large range of \fesc\ for individual objects. The implied median \fesc-values may be high enough for galaxies to finish reionization by $z\sim6$. Furthermore, we find no strong trends between the UV-slope $β$ or \EBminV with model \fesc. If true, the lack of trends suggest that other factors besides nebular emission or dust extinction could have led to LyC escaping, such as the presence of holes in the ISM with sufficiently wide opening angles from outflows of supernovae and/or weak AGN, resulting in a range of implied \fesc-values depending on the viewing angle of each galaxy. The current \textit{HST, Spitzer} and ground-based photometric and model errors for the galaxies remain large, so IR spectroscopic observations with the \textit{James Webb Space Telescope} are needed to constrain this possibility. △ Less

Submitted 18 March, 2022; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: 24 pages, 5 figures, 3 tables

IQ Collaboratory II: The Quiescent Fraction of Isolated, Low Mass Galaxies Across Simulations and Observations

Authors: Claire M Dickey, Tjitske K Starkenburg, Marla Geha, ChangHoon Hahn, Daniel Anglés-Alcázar, Ena Choi, Romeel Davé, Shy Genel, Kartheik G Iyer, Ariyeh H Maller, Nir Mandelker, Rachel S Somerville, L Y Aaron Yung

Abstract: We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each s… ▽ More We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each simulation, adding realistic noise and observational limits. All three simulations show a decrease in the number of quiescent, isolated galaxies in the mass range $\mathrm{M}_* = 10^{9-10} \ \mathrm{M}_\odot$, in broad agreement with observations. However, even after accounting for observational and selection biases, none of the simulations reproduce the observed absence of quiescent field galaxies below $\mathrm{M}_*=10^{9} \ \mathrm{M}_\odot$. We find that the low mass quiescent populations selected via synthetic observations have consistent quenching timescales, despite apparent variation in the late time star formation histories. The effect of increased numerical resolution is not uniform across simulations and cannot fully mitigate the differences between the simulations and the observations. The framework presented here demonstrates a path towards more robust and accurate comparisons between theoretical simulations and galaxy survey observations, while the quenching threshold serves as a sensitive probe of feedback implementations. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: 19 pages, 8 figures. Figure 4 presents the main result. Code used in this work may be accessed at github.com/IQcollaboratory/orchard. Submitted to ApJ

The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution

Authors: Kartheik G. Iyer, Sandro Tacchella, Shy Genel, Christopher C. Hayward, Lars Hernquist, Alyson M. Brooks, Neven Caplar, Romeel Davé, Benedikt Diemer, John C. Forbes, Eric Gawiser, Rachel S. Somerville, Tjitske K. Starkenburg

Abstract: Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIR… ▽ More Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability. △ Less

Submitted 15 July, 2020; originally announced July 2020.

Comments: 31 pages, 17 figures (+ appendix). Resubmitted to MNRAS after responding to referee's comments. Comments are welcome!

Evaluation of probabilistic photometric redshift estimation approaches for The Rubin Observatory Legacy Survey of Space and Time (LSST)

Authors: S. J. Schmidt, A. I. Malz, J. Y. H. Soo, I. A. Almosallam, M. Brescia, S. Cavuoti, J. Cohen-Tanugi, A. J. Connolly, J. DeRose, P. E. Freeman, M. L. Graham, K. G. Iyer, M. J. Jarvis, J. B. Kalmbach, E. Kovacs, A. B. Lee, G. Longo, C. B. Morrison, J. A. Newman, E. Nourbakhsh, E. Nuss, T. Pospisil, H. Tranin, R. H. Wechsler, R. Zhou , et al. (2 additional authors not shown)

Abstract: Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive exper… ▽ More Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced for The Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. We report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performancemetrics. △ Less

Submitted 31 July, 2021; v1 submitted 10 January, 2020; originally announced January 2020.

Journal ref: MNRAS 499 2 1587 (2020)

Learning the Relationship between Galaxies Spectra and their Star Formation Histories using Convolutional Neural Networks and Cosmological Simulations

Authors: Christopher C. Lovell, Viviana Acquaviva, Peter A. Thomas, Kartheik G. Iyer, Eric Gawiser, Stephen M. Wilkins

Abstract: We present a new method for inferring galaxy star formation histories (SFH) using machine learning methods coupled with two cosmological hydrodynamic simulations. We train Convolutional Neural Networks to learn the relationship between synthetic galaxy spectra and high resolution SFHs from the EAGLE and Illustris models. To evaluate our SFH reconstruction we use Symmetric Mean Absolute Percentage… ▽ More We present a new method for inferring galaxy star formation histories (SFH) using machine learning methods coupled with two cosmological hydrodynamic simulations. We train Convolutional Neural Networks to learn the relationship between synthetic galaxy spectra and high resolution SFHs from the EAGLE and Illustris models. To evaluate our SFH reconstruction we use Symmetric Mean Absolute Percentage Error (SMAPE), which acts as a true percentage error in the low-error regime. On dust-attenuated spectra we achieve high test accuracy (median SMAPE $= 10.5\%$). Including the effects of simulated observational noise increases the error ($12.5\%$), however this is alleviated by including multiple realisations of the noise, which increases the training set size and reduces overfitting ($10.9\%$). We also make estimates for the observational and modelling errors. To further evaluate the generalisation properties we apply models trained on one simulation to spectra from the other, which leads to only a small increase in the error (median SMAPE $\sim 15\%$). We apply each trained model to SDSS DR7 spectra, and find smoother histories than in the VESPA catalogue. This new approach complements the results of existing SED fitting techniques, providing star formation histories directly motivated by the results of the latest cosmological simulations. △ Less

Submitted 9 October, 2019; v1 submitted 25 March, 2019; originally announced March 2019.

Comments: 20 pages, 21 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

Non-parametric Star Formation History Reconstruction with Gaussian Processes I: Counting Major Episodes of Star Formation

Authors: Kartheik G. Iyer, Eric Gawiser, Sandra M. Faber, Henry C. Ferguson, Anton M. Koekemoer, Camilla Pacifici, Rachel Somerville

Abstract: The star formation histories (SFHs) of galaxies contain imprints of the physical processes responsible for regulating star formation during galaxy growth and quenching. We improve the Dense Basis SFH reconstruction method of Iyer & Gawiser (2017), introducing a nonparametric description of the SFH based on the lookback times at which a galaxy assembles certain quantiles of its stellar mass. The me… ▽ More The star formation histories (SFHs) of galaxies contain imprints of the physical processes responsible for regulating star formation during galaxy growth and quenching. We improve the Dense Basis SFH reconstruction method of Iyer & Gawiser (2017), introducing a nonparametric description of the SFH based on the lookback times at which a galaxy assembles certain quantiles of its stellar mass. The method uses Gaussian Processes to create smooth SFHs that are independent of any functional form, with a flexible number of parameters that is adjusted to extract the maximum possible amount of SFH information from the SEDs being fit. We apply the method to reconstruct the SFHs of 48,791 galaxies with $H<25$ at $0.5 < z < 3.0$ across the five CANDELS fields. Using these SFHs, we study the evolution of galaxies as they grow more massive over cosmic time. We quantify the fraction of galaxies that show multiple major episodes of star formation, finding that the median time between two peaks of star formation is $\sim 0.42_{-0.10}^{+0.15}t_{univ}$ Gyr, where $t_{univ}$ is the age of the universe at a given redshift and remains roughly constant with stellar mass. Correlating SFHs with morphology, we find that studying the median SFHs of galaxies at $0.5<z<1.0$ at the same mass ($10^{10}< M_* < 10^{10.5}M_\odot$) allows us to compare the timescales on which the SFHs decline for different morphological classifications, ranging from $0.60^{-0.54}_{+1.54}$ Gyr for galaxies with spiral arms to $2.50^{-1.50}_{+2.25}$ Gyr for spheroids. The Gaussian Process-based SFH description provides a general approach to reconstruct smooth, nonparametric SFH posteriors for galaxies with a flexible number of parameters that can be incorporated into Bayesian SED fitting codes to minimize the bias in estimating physical parameters due to SFH parametrization. △ Less

Submitted 9 January, 2019; originally announced January 2019.

Comments: 26 pages, 15 figures. Submitted to ApJ

The SFR-M$_*$ Correlation Extends to Low Mass at High Redshift

Authors: Kartheik G. Iyer, Eric Gawiser, Romeel Davé, Philip Davis, Steven L. Finkelstein, Dritan Kodra, Anton M. Koekemoer, Peter Kurczynski, Jeffery A. Newman, Camilla Pacifici, Rachel Somerville

Abstract: To achieve a fuller understanding of galaxy evolution, SED fitting can be used to recover quantities beyond stellar masses (M$_*$) and star formation rates (SFRs). We use Star Formation Histories (SFHs) reconstructed via the Dense Basis method of Iyer \& Gawiser (2017) for a sample of $17,873$ galaxies at $0.5<z<6$ in the CANDELS GOODS-S field to study the nature and evolution of the SFR-M$_*$ cor… ▽ More To achieve a fuller understanding of galaxy evolution, SED fitting can be used to recover quantities beyond stellar masses (M$_*$) and star formation rates (SFRs). We use Star Formation Histories (SFHs) reconstructed via the Dense Basis method of Iyer \& Gawiser (2017) for a sample of $17,873$ galaxies at $0.5<z<6$ in the CANDELS GOODS-S field to study the nature and evolution of the SFR-M$_*$ correlation. The reconstructed SFHs represent trajectories in SFR-M$_*$ space, enabling us to study galaxies at epochs earlier than observed by propagating them backwards in time along these trajectories. We study the SFR-M$_*$ correlation at $z=1,2,3,4,5,6$ using both direct fits to galaxies observed at those epochs and SFR-M$_*$ trajectories of galaxies observed at lower redshifts. The SFR-M$_*$ correlations obtained using the two approaches are found to be consistent with each other through a KS test. Validation tests using SFHs from semi-analytic models and cosmological hydrodynamical simulations confirm the sensitivity of the method to changes in the slope, normalization and shape of the SFR-M$_*$ correlation. This technique allows us to further probe the low-mass regime of the correlation at high-z by $\sim 1$ dex and over an effective volume of $\sim 10\times$ larger than possible with just direct fits. We find that the SFR-M$_*$ correlation is consistent with being linear down to M$_*\sim 10^7 M_\odot$ at $z>4$. The evolution of the correlation is well described by $\log SFR= (0.80\pm 0.029 - 0.017\pm 0.010\times t_{univ})\log M_*$ $- (6.487\pm 0.282-0.039\pm 0.008\times t_{univ})$, where $t_{univ}$ is the age of the universe in Gyr. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 22 pages, 10 figures. Accepted for publication in ApJ

Reconstruction of Galaxy Star Formation Histories through SED Fitting: The Dense Basis Approach

Authors: Kartheik G. Iyer, Eric Gawiser

Abstract: We introduce the Dense Basis method for Spectral Energy Distribution (SED) fitting. It accurately recovers traditional SED parameters, including M$_*$, SFR and dust attenuation, and reveals previously inaccessible information about the number and duration of star formation episodes and the timing of stellar mass assembly, as well as uncertainties in these quantities. This is done using basis Star… ▽ More We introduce the Dense Basis method for Spectral Energy Distribution (SED) fitting. It accurately recovers traditional SED parameters, including M$_*$, SFR and dust attenuation, and reveals previously inaccessible information about the number and duration of star formation episodes and the timing of stellar mass assembly, as well as uncertainties in these quantities. This is done using basis Star Formation Histories (SFHs) chosen by comparing the goodness-of-fit of mock galaxy SEDs to the goodness-of-reconstruction of their SFHs. We train and validate the method using a sample of realistic SFHs at $z =1$ drawn from stochastic realisations, semi-analytic models, and a cosmological hydrodynamical galaxy formation simulation. The method is then applied to a sample of 1100 CANDELS GOODS-S galaxies at $1<z<1.5$ to illustrate its capabilities at moderate S/N with 15 photometric bands. Of the six parametrizations of SFHs considered, we adopt linear-exponential, bessel-exponential, lognormal and gaussian SFHs and reject the traditional parametrizations of constant (Top-Hat) and exponential SFHs. We quantify the bias and scatter of each parametrization. $15\%$ of galaxies in our CANDELS sample exhibit multiple episodes of star formation, with this fraction decreasing above $M_*>10^{9.5}M_\odot$. About $40\%$ of the CANDELS galaxies have SFHs whose maximum occurs at or near the epoch of observation. The Dense Basis method is scalable and offers a general approach to a broad class of data-science problems. △ Less

Submitted 14 February, 2017; originally announced February 2017.

Comments: 27 pages including 20 color figures, revised in response to referee's report