Showing 1–50 of 495 results for author: Bell, E

GD-1 Stellar Stream and Cocoon in the DESI Early Data Release

Authors: Monica Valluri, Parker Fagrelius, Sergey. E. Koposov, Ting S. Li, Oleg Y. Gnedin, Eric F. Bell, Raymond G. Carlberg, Andrew P. Cooper, Jessia N. Aguilar, Carlos Allende Prieto, Vasily Belokurov, Leandro Beraldo e Silva, David Brooks, Amanda Byström, Todd Claybaugh, Kyle Dawson, Arjun Dey, Peter Doel, Jaime E. Forero-Romero, Enrique Gaztañaga, Satya Gontcho A Gontcho, Klaus Honscheid, T . Kisner, Anthony Kremin, A. Lambert , et al. (27 additional authors not shown)

Abstract: We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream compo… ▽ More We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream component (based on a median per star velocity precision of 2.7km/s). The cocoon extends over at least a ~ 20deg segment of the stream observed by DESI. The thin and cocoon components have similar mean values of [Fe/H]: -2.54+/- 0.04dex and -2.45+/-0.06dex suggestive of a common origin. The data are consistent with the following scenarios for the origin of the cocoon. The progenitor of the GD-1 stream was an accreted globular cluster (GC) and: (a) the cocoon was produced by pre-accretion tidal stripping of the GC while it was still inside its parent dwarf galaxy; (b) the cocoon is debris from the parent dwarf galaxy; (c) an initially thin GC tidal stream was heated by impacts from dark subhalos in the Milky Way; (d) an initially thin GC stream was heated by a massive Sagittarius dwarf galaxy; or a combination of some these. In the first two cases the velocity dispersion and mean metallicity are consistent with the parent dwarf galaxy having a halo mass of ~0^9\msun. Future DESI spectroscopy and detailed modeling may enable us to distinguish between these possible origins. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ, 23 pages, 13 figures 4 tables

The Size-Mass relation at Rest-Frame $1.5μ$m from JWST/NIRCam in the COSMOS-WEB and PRIMER-COSMOS fields

Authors: Marco Martorano, Arjen van der Wel, Maarten Baes, Eric F. Bell, Gabriel Brammer, Marijn Franx, Angelos Nersesian

Abstract: We present the galaxy stellar mass - size relation in the rest-frame near-IR ($1.5~μ{\text{m}}$) and its evolution with redshift up to $z=2.5$. Sérsic profiles are measured for $\sim$ $26\,000$ galaxies with stellar masses $M_\star > 10^9~{\text{M}}_\odot$ from JWST/NIRCam F277W and F444W imaging provided by the COSMOS-WEB and PRIMER surveys, using coordinates, redshifts, colors and stellar mass e… ▽ More We present the galaxy stellar mass - size relation in the rest-frame near-IR ($1.5~μ{\text{m}}$) and its evolution with redshift up to $z=2.5$. Sérsic profiles are measured for $\sim$ $26\,000$ galaxies with stellar masses $M_\star > 10^9~{\text{M}}_\odot$ from JWST/NIRCam F277W and F444W imaging provided by the COSMOS-WEB and PRIMER surveys, using coordinates, redshifts, colors and stellar mass estimates from the COSMOS2020 catalog. The new rest-frame near-IR effective radii are generally smaller than previously measured rest-frame optical sizes, on average by 0.14~dex, with no significant dependence on redshift. For quiescent galaxies this size offset does not depend on stellar mass, but for star-forming galaxies the offset increases from -0.1~dex at $M_\star = 10^{9.5}~{\text{M}}_\odot$ to -0.25~dex at $M_\star > 10^{11}~{\text{M}}_\odot$. That is, we find that the near-IR stellar mass - size relation for star-forming galaxies is flatter in the rest-frame near-IR than in the rest-frame optical at all redshifts $0.5<z<2.5$. The general pace of size evolution is the same in the near-IR as previously demonstrated in the optical, with slower evolution ($R_{\text{e}} \propto (1+z)^{-0.7}$) for $L^*$~star-forming galaxies and faster evolution ($R_{\text{e}} \propto (1+z)^{-1.3}$) for $L^*$~quiescent galaxies. Massive ($M_\star>10^{11}~{\text{M}}_\odot$) star-forming galaxies evolve in size almost as fast as quiescent galaxies. Low-mass ($M_\star<10^{10}~{\text{M}}_\odot$)~quiescent galaxies evolve as slow as star-forming galaxies. Our main conclusion is that the size evolution narrative as it has emerged over the past two decades does not radically change when accessing with JWST the rest-frame near-IR, a better proxy of the underlying stellar mass distribution. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Accepted for publication

First resolved stellar halo kinematics of a MW-mass galaxy outside the Local Group: A flat counter-rotating halo in NGC 4945

Authors: Camila Beltrand, Antonela Monachesi, Richard D'Souza, Eric F. Bell, Roelof S. de Jong, Facundo A. Gomez, Jeremy Bailin, In Sung Jang, Adam Smercina

Abstract: Stellar halos of galaxies, primarily formed through the accretion of smaller objects, are important to understand the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in-situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in-situ com… ▽ More Stellar halos of galaxies, primarily formed through the accretion of smaller objects, are important to understand the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in-situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in-situ component from the accreted stellar halos. The low surface brightness of stellar halos makes it inaccessible with traditional integrated light spectroscopy. In this work, using a novel technique, we study the kinematics of the stellar halo of the edge-on galaxy NGC 4945. We couple new deep Multi Unit Spectroscopic Explorer spectroscopic observations with existing Hubble Space Telescope imaging data to spectroscopically measure the line-of-sight (LOS) heliocentric velocity and velocity dispersion in two fields at a galactocentric distance of 12.2 kpc (outer disk field) and 34.6 kpc (stellar halo field) along NGC 4945 major axis, by stacking individual spectra of red giant branch and asymptotic giant branch stars. We obtain a LOS velocity and dispersion of 673+/-11 km/s and 73+/-14 km/s, respectively, for the outer disk field. This is consistent with the mean HI velocity of the disk at that distance. For the halo field we obtain a LOS velocity and dispersion of 519+/-12 km/s and 42+/-22 km/s. The halo fields' velocity measurement is within ~40 km/s from the systemic LOS velocity of NGC 4945, which is 563 km/s, suggesting that its stellar halo at 34.6 kpc along the major axis is counter-rotating and is of likely accretion origin. This provides the first ever kinematic measurement of the stellar halo of a Milky Way-mass galaxy outside the Local Group from its resolved stellar population, and establishes a powerful technique for measuring the velocity field of the stellar halos of nearby galaxies. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Accepted to A&A. 23 pages, 12 figures, 2 appendices

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). VI. The High-Mass Stellar Initial Mass Function of M33

Authors: Tobin M. Wainer, Benjamin F. Williams, L. Clifton Johnson, Daniel R. Weisz, Julianne J. Dalcanton, Anil C. Seth, Andrew Dolphin, Meredith J. Durbin, Eric F. Bell, Zhuo Chen, Puragra Guhathakurta, Eric W. Koch, Christina W. Lindberg, Erik Rosolowsky, Karin M. Sandstrom, Evan D. Skillman, Adam Smercina, Estephani E. TorresVillanueva

Abstract: We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $Γ$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). W… ▽ More We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $Γ$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). We then probabilistically model the distribution of MF slopes for a highly strict cluster sample of 9 clusters more massive than log(Mass/M$_{\odot}$)=3.6; above this mass, all clusters have well-populated main sequences of massive stars and should have accurate recovery of their MF slopes, based on extensive tests with artificial clusters. We find the ensemble IMF is best described by a mean high-mass slope of $\overlineΓ = 1.49\pm0.18$, with an intrinsic scatter of $σ^{2}_Γ = 0.02^{+0.16}_{0.00}$, consistent with a universal IMF. We find no dependence of the IMF on environmental impacts such as the local star formation rate or galactocentric radius within M33, which serves as a proxy for metallicity. This $\overlineΓ$ measurement is consistent with similar measurements in M31, despite M33 having a much higher star formation rate intensity. While this measurement is formally consistent with the canonical Kroupa ($Γ= 1.30$) IMF, as well as the Salpeter ($Γ= 1.35)$) value, it is the second Local Group cluster sample to show evidence for a somewhat steeper high-mass IMF slope. We explore the impacts a steeper IMF slope has on a number of astronomical sub-fields. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJ. 9 Figures, 1 Table

Less is less: photometry alone cannot predict the observed spectral indices of $z\sim1$ galaxies from the LEGA-C spectroscopic survey

Authors: Angelos Nersesian, Arjen van der Wel, Anna Gallazzi, Joel Leja, Rachel Bezanson, Eric F. Bell, Francesco D'Eugenio, Anna de Graaff, Yasha Kaushal, Marco Martorano, Michael Maseda, Stefano Zibetti

Abstract: We test whether we can predict optical spectra from deep-field photometry of distant galaxies. Our goal is to perform a comparison in data space, highlighting the differences between predicted and observed spectra. The Large Early Galaxy Astrophysics Census (LEGA-C) provides high-quality optical spectra of thousands of galaxies at redshift $0.6<z<1$. Broad-band photometry of the same galaxies, dra… ▽ More We test whether we can predict optical spectra from deep-field photometry of distant galaxies. Our goal is to perform a comparison in data space, highlighting the differences between predicted and observed spectra. The Large Early Galaxy Astrophysics Census (LEGA-C) provides high-quality optical spectra of thousands of galaxies at redshift $0.6<z<1$. Broad-band photometry of the same galaxies, drawn from the recent COSMOS2020 catalog, is used to predict the optical spectra with the spectral energy distribution (SED) fitting code Prospector and the MILES stellar library. The observed and predicted spectra are compared in terms of two age and metallicity-sensitive absorption features (H$δ_\mathrm{A}$ and Fe4383). The global bimodality of star-forming and quiescent galaxies in photometric space is recovered with the model spectra. But the presence of a systematic offset in the Fe4383 line strength and the weak correlation between the observed and modeled line strength imply that accurate age or metallicity determinations cannot be inferred from photometry alone. For now we caution that photometry-based estimates of stellar population properties are determined mostly by the modeling approach and not the physical properties of galaxies, even when using the highest-quality photometric datasets and state-of-the-art fitting techniques. When exploring a new physical parameter space (i.e. redshift or galaxy mass) high-quality spectroscopy is always needed to inform the analysis of photometry. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 13 pages, 8 figures, accepted 26 October 2023

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

CEERS: 7.7 $μ$m PAH Star Formation Rate Calibration with JWST MIRI

Authors: Kaila Ronayne, Casey Papovich, Guang Yang, Lu Shen, Mark Dickinson, Robert Kennicutt, Anahita Alavi, Pablo Arrabal Haro, Micaela Bagley, Denis Burgarella, Aurélien Le Bail, Eric Bell, Nikko Cleri, Justin Cole, Luca Costantin, Alexander de la Vega, Emanuele Daddi, David Elbaz, Steven Finkelstein, Norman Grogin, Benne Holwerda, Jeyhan Kartaltepe, Allison Kirkpatrick, Anton Koekemoer, Ray Lucas , et al. (11 additional authors not shown)

Abstract: We test the relationship between UV-derived star formation rates (SFRs) and the 7.7 $μ$m polycyclic aromatic hydrocarbon (PAH) luminosities from the integrated emission of galaxies at z ~ 0 - 2. We utilize multi-band photometry covering 0.2 - 160 $μ$m from HST, CFHT, JWST, Spitzer, and Herschel for galaxies in the Cosmic Evolution Early Release Science (CEERS) Survey. We perform spectral energy di… ▽ More We test the relationship between UV-derived star formation rates (SFRs) and the 7.7 $μ$m polycyclic aromatic hydrocarbon (PAH) luminosities from the integrated emission of galaxies at z ~ 0 - 2. We utilize multi-band photometry covering 0.2 - 160 $μ$m from HST, CFHT, JWST, Spitzer, and Herschel for galaxies in the Cosmic Evolution Early Release Science (CEERS) Survey. We perform spectral energy distribution (SED) modeling of these data to measure dust-corrected far-UV (FUV) luminosities, $L_{FUV}$, and UV-derived SFRs. We then fit SED models to the JWST/MIRI 7.7 - 21 $μ$m CEERS data to derive rest-frame 7.7 $μ$m luminosities, $L_{770}$, using the average flux density in the rest-frame MIRI F770W bandpass. We observe a correlation between $L_{770}$ and $L_{FUV}$, where log $L_{770}$ is proportional to (1.27+/-0.04) log $L_{FUV}$. $L_{770}$ diverges from this relation for galaxies at lower metallicities, lower dust obscuration, and for galaxies dominated by evolved stellar populations. We derive a "single-wavelength" SFR calibration for $L_{770}$ which has a scatter from model estimated SFRs (${σ_{ΔSFR}}$) of 0.24 dex. We derive a "multi-wavelength" calibration for the linear combination of the observed FUV luminosity (uncorrected for dust) and the rest-frame 7.7 $μ$m luminosity, which has a scatter of ${σ_{ΔSFR}}$ = 0.21 dex. The relatively small decrease in $σ$ suggests this is near the systematic accuracy of the total SFRs using either calibration. These results demonstrate that the rest-frame 7.7 $μ$m emission constrained by JWST/MIRI is a tracer of the SFR for distant galaxies to this accuracy, provided the galaxies are dominated by star-formation with moderate-to-high levels of attenuation and metallicity. △ Less

Submitted 13 October, 2023; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: 20 pages, 11 figures, 2 tables, submitted to ApJ

Black Hole Growth, Baryon Lifting, Star Formation, and IllustrisTNG

Authors: G. M. Voit, B. D. Oppenheimer, E. F. Bell, B. Terrazas, M. Donahue

Abstract: Quenching of star formation in the central galaxies of cosmological halos is thought to result from energy released as gas accretes onto a supermassive black hole. The same energy source also appears to lower the central density and raise the cooling time of baryonic atmospheres in massive halos, thereby limiting both star formation and black hole growth, by lifting the baryons in those halos to g… ▽ More Quenching of star formation in the central galaxies of cosmological halos is thought to result from energy released as gas accretes onto a supermassive black hole. The same energy source also appears to lower the central density and raise the cooling time of baryonic atmospheres in massive halos, thereby limiting both star formation and black hole growth, by lifting the baryons in those halos to greater altitudes. One predicted signature of that feedback mechanism is a nearly linear relationship between the central black hole's mass (MBH) and the original binding energy of the halo's baryons. We present the increasingly strong observational evidence supporting a such a relationship, showing that it extends up to halos of mass Mhalo ~10^14 MSun. We then compare current observational constraints on the MBH--Mhalo relation with numerical simulations, finding that black hole masses in IllustrisTNG appear to exceed those constraints at Mhalo < 10^13 MSun and that black hole masses in EAGLE fall short of observations at Mhalo ~ 10^14 MSun. A closer look at IllustrisTNG shows that quenching of star formation and suppression of black hole growth do indeed coincide with black hole energy input that lifts the halo's baryons. However, IllustrisTNG does not reproduce the observed MBH--Mhalo relation because its black holes gain mass primarily through accretion that does not contribute to baryon lifting. We suggest adjustments to some of the parameters in the IllustrisTNG feedback algorithm that may allow the resulting black hole masses to reflect the inherent links between black hole growth, baryon lifting, and star formation among the massive galaxies in those simulations. △ Less

Submitted 3 October, 2023; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: 20 pages, 10 figures, Accepted to ApJ (Figures 5 and 6 updated)

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). V. The Structure of M33 in Resolved Stellar Populations

Authors: Adam Smercina, Julianne J. Dalcanton, Benjamin F. Williams, Meredith J. Durbin, Margaret Lazzarini, Eric F. Bell, Yumi Choi, Andrew Dolphin, Karoline Gilbert, Puragra Guhathakurta, Eric W. Koch, Hans-Walter Rix, Erik Rosolowsky, Anil Seth, Evan D. Skillman, Daniel R. Weisz

Abstract: We present a detailed analysis of the the structure of the Local Group flocculent spiral galaxy M33, as measured using the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey. Leveraging the multiwavelength coverage of PHATTER, we find that the oldest populations are dominated by a smooth exponential disk with two distinct spiral arms and a classical central bar $-$… ▽ More We present a detailed analysis of the the structure of the Local Group flocculent spiral galaxy M33, as measured using the Panchromatic Hubble Andromeda Treasury Triangulum Extended Region (PHATTER) survey. Leveraging the multiwavelength coverage of PHATTER, we find that the oldest populations are dominated by a smooth exponential disk with two distinct spiral arms and a classical central bar $-$ completely distinct from what is seen in broadband optical imaging, and the first-ever confirmation of a bar in M33. We estimate a bar extent of $\sim$1 kpc. The two spiral arms are asymmetric in orientation and strength, and likely represent the innermost impact of the recent tidal interaction responsible for M33's warp at larger scales. The flocculent multi-armed morphology for which M33 is known is only visible in the young upper main sequence population, which closely tracks the morphology of the ISM. We investigate the stability of M33's disk, finding $Q{\sim}1$ over the majority of the disk. We fit multiple components to the old stellar density distribution and find that, when considering recent stellar kinematics, M33's bulk structure favors the inclusion of an accreted halo component, modeled as a broken power-law. The best-fit halo model has an outer power-law index of $-$3 and accurately describes observational evidence of M33's stellar halo from both resolved stellar spectroscopy in the disk and its stellar populations at large radius. Integrating this profile yields a total halo stellar mass of ${\sim}5{\times}10^8\ M_{\odot}$, giving a total stellar halo mass fraction of 16%, most of which resides in the innermost 2.5 kpc. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 31 pages, 17 figures, 5 tables, Accepted for publication in ApJ

Constraining the assembly time of the stellar haloes of nearby Milky Way-mass galaxies through AGB populations

Authors: Benjamin Harmsen, Eric F. Bell, Richard D'Souza, Antonela Monachesi, Roelof S. de Jong, Adam Smercina, In Sung Jang, Benne W. Holwerda

Abstract: The star formation histories (SFHs) of galactic stellar haloes offer crucial insights into the merger history of the galaxy and the effects of those mergers on their hosts. Such measurements have revealed that while the Milky Way's most important merger was 8-10 Gyr ago, M31's largest merger was more recent, within the last few Gyr. Unfortunately, the required halo SFH measurements are extremely o… ▽ More The star formation histories (SFHs) of galactic stellar haloes offer crucial insights into the merger history of the galaxy and the effects of those mergers on their hosts. Such measurements have revealed that while the Milky Way's most important merger was 8-10 Gyr ago, M31's largest merger was more recent, within the last few Gyr. Unfortunately, the required halo SFH measurements are extremely observationally expensive outside of the Local Group. Here we use asymptotic giant branch (AGB) stars brighter than the tip of the red giant branch (RGB) to constrain stellar halo SFHs. Both stellar population models and archival datasets show that the AGB/RGB ratio constrains the time before which 90% of the stars formed, $t_{90}$. We find AGB stars in the haloes of three highly-inclined roughly Milky Way-mass galaxies with resolved star measurements from the Hubble Space Telescope; this population is most prominent in the stellar haloes of NGC 253 and NGC 891, suggesting that their stellar haloes contain stars born at relatively late times, with inferred $t_{90}\sim 6\pm1.5$Gyr. This ratio also varies from region to region, tending towards higher values along the major axis and in tidal streams or shells. By combining our measurements with previous constraints, we find a tentative anticorrelation between halo age and stellar halo mass, a trend that exists in models of galaxy formation but has never been elucidated before, i.e, the largest stellar haloes of Milky-Way mass galaxies were assembled more recently. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 20 Pages, 10 Figures

Rest-Frame Near-Infrared Radial Light Profiles up to z=3 from JWST/NIRCam: Wavelength Dependence of the Sérsic Index

Authors: Marco Martorano, Arjen van der Wel, Eric F. Bell, Marijn Franx, Katherine E. Whitaker, Angelos Nersesian, Sedona H. Price, Maarten Baes, Katherine A. Suess, Erica J. Nelson, Tim B. Miller, Rachel Bezanson, Gabriel Brammer

Abstract: We examine the wavelength dependence of radial light profiles based on Sérsic index $n$ measurements of 1067 galaxies with M$_*\geq$ 10$^{9.5}$M$_\odot$ and in the redshift range $0.5 < z < 3$. The sample and rest-frame optical light profiles are drawn from CANDELS$+$3D-HST; rest-frame near-infrared light profiles are inferred from CEERS JWST/NIRCam imaging. $n$ shows only weak dependence on wavel… ▽ More We examine the wavelength dependence of radial light profiles based on Sérsic index $n$ measurements of 1067 galaxies with M$_*\geq$ 10$^{9.5}$M$_\odot$ and in the redshift range $0.5 < z < 3$. The sample and rest-frame optical light profiles are drawn from CANDELS$+$3D-HST; rest-frame near-infrared light profiles are inferred from CEERS JWST/NIRCam imaging. $n$ shows only weak dependence on wavelength, regardless of redshift, galaxy mass and type: on average, star-forming galaxies have $n = 1-1.5$ and quiescent galaxies have $n = 3-4$ in the rest-frame optical and near-infrared. The strong correlation at all wavelengths between $n$ and star-formation activity implies a physical connection between the radial stellar mass profile and star-formation activity. The main caveat is that the current sample is too small to discern trends for the most massive galaxies (M$_* > 10^{11}M_\odot$). △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 15 pages, 7 figures, Submitted to ApJ. Comments welcome

Journal ref: ApJ 957 46 (2023)

CEERS Key Paper VII: JWST/MIRI Reveals a Faint Population of Galaxies at Cosmic Noon Unseen by Spitzer

Authors: Allison Kirkpatrick, Guang Yang, Aurelien Le Bail, Greg Troiani, Eric F. Bell, Nikko J. Cleri, David Elbaz, Steven L. Finkelstein, Nimish P. Hathi, Michaela Hirschmann, Benne W. Holwerda, Dale D. Kocevski, Ray A. Lucas, Jed McKinney, Casey Papovich, Pablo G. Perez-Gonzalez, Alexander de la Vega, Micaela B. Bagley, Emanuele Daddi, Mark Dickinson, Henry C. Ferguson, Adriano Fontana, Andrea Grazian, Norman A. Grogin, Pablo Arrabal Haro , et al. (11 additional authors not shown)

Abstract: The Cosmic Evolution Early Release Science (CEERS) program observed the Extended Groth Strip with the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) in 2022. In this paper, we discuss the four MIRI pointings that observed with longer wavelength filters, including F770W, F1000W, F1280W, F1500W, F1800W, and F2100W. We compare the MIRI galaxies with the Spitzer/MIPS 24$μ$m po… ▽ More The Cosmic Evolution Early Release Science (CEERS) program observed the Extended Groth Strip with the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) in 2022. In this paper, we discuss the four MIRI pointings that observed with longer wavelength filters, including F770W, F1000W, F1280W, F1500W, F1800W, and F2100W. We compare the MIRI galaxies with the Spitzer/MIPS 24$μ$m population in the EGS field. We find that MIRI can observe an order of magnitude deeper than MIPS in significantly shorter integration times, attributable to JWST's much larger aperture and MIRI's improved sensitivity. MIRI is exceptionally good at finding faint ($L_{\rm IR}<10^{10} L_\odot$) galaxies at $z\sim1-2$. We find that a significant portion of MIRI galaxies are "mid-IR weak"--they have strong near-IR emission and relatively weaker mid-IR emission, and most of the star formation is unobscured. We present new IR templates that capture how the mid-IR to near-IR emission changes with increasing infrared luminosity. We present two color-color diagrams to separate mid-IR weak galaxies and active galactic nuclei (AGN) from dusty star-forming galaxies and find that these color diagrams are most effective when used in conjunction with each other. We present the first number counts of 10$μ$m sources and find that there are $\lesssim10$ IR AGN per MIRI pointing, possibly due to the difficulty of distinguishing AGN from intrinsically mid-IR weak galaxies (due to low metallicities or low dust content). We conclude that MIRI is most effective at observing moderate luminosity ($L_{\rm IR}=10^9-10^{10}L_\odot$) galaxies at $z=1-2$, and that photometry alone is not effective at identifying AGN within this faint population. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 21 pages, 10 figures. Resubmitted to ApJS after revision

The Panchromatic Hubble Andromeda Treasury XXI. The Legacy Resolved Stellar Photometry Catalog

Authors: Benjamin F. Williams, Meredith Durbin, Dustin Lang, Julianne J. Dalcanton, Andrew E. Dolphin, Adam Smercina, Petia Yanchulova Merica-Jones, Daniel R. Weisz, Eric F. Bell, Karoline M. Gilbert, Leo Girardi, Karl Gordon, Puragra Guhathakurta, L. Clifton Johnson, Tod R. Lauer, Anil Seth, Evan Skillman

Abstract: We present the final legacy version of stellar photometry for the Panchromatic Hubble Andromeda Treasury (PHAT) survey. We have reprocessed all of the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) near ultraviolet (F275W, F336W), optical (F475W, F814W), and near infrared (F110W, F160W) imaging from the PHAT survey using an improved method that optimi… ▽ More We present the final legacy version of stellar photometry for the Panchromatic Hubble Andromeda Treasury (PHAT) survey. We have reprocessed all of the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) near ultraviolet (F275W, F336W), optical (F475W, F814W), and near infrared (F110W, F160W) imaging from the PHAT survey using an improved method that optimized the survey depth and chip gap coverage by including all overlapping exposures in all bands in the photometry. An additional improvement was gained through the use of charge transfer efficiency (CTE) corrected input images, which provide more complete star finding as well as more reliable photometry for the NUV bands, which had no CTE correction in the previous version of the PHAT photometry. While this method requires significantly more computing resources and time than earlier versions where the photometry was performed on individual pointings, it results in smaller systematic instrumental completeness variations as demonstrated by cleaner maps in stellar density, and it results in optimal constraints on stellar fluxes in all bands from the survey data. Our resulting catalog has 138 million stars, 18% more than the previous catalog, with lower density regions gaining as much as 40% more stars. The new catalog produces nearly seamless population maps which show relatively well-mixed distributions for populations associated with ages older than 1-2 Gyr, and highly structured distributions for the younger populations. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 28 pages, 4 tables, 18 figures, accepted for publication in The Astrophysical Journal Supplement Series

A census of star formation histories of massive galaxies at 0.6 < z < 1 from spectro-photometric modeling using Bagpipes and Prospector

Authors: Yasha Kaushal, Angelos Nersesian, Rachel Bezanson, Arjen van der Wel, Joel Leja, Adam Carnall, Stefano Zibetti, Gourav Khullar, Marijn Franx, Adam Muzzin, Anna De Graaff, Camilla Pacifici, Katherine E. Whitaker, Eric F. Bell, Marco Martorano

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

Submitted 9 November, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

Comments: 28 pages, 12 figures, 4 tables, Accepted to ApJ

Stellar Half-Mass Radii of $0.5<z<2.3$ Galaxies: Comparison with JWST/NIRCam Half-Light Radii

Authors: Arjen van der Wel, Marco Martorano, Boris Haussler, Kalina V. Nedkova, Tim B. Miller, Gabriel B. Brammer, Glenn van de Ven, Joel Leja, Rachel S. Bezanson, Adam Muzzin, Danilo Marchesini, Anna de Graaff, Mariska Kriek, Eric F. Bell, Marijn Franx

Abstract: We use CEERS JWST/NIRCam imaging to measure rest-frame near-IR light profiles of $>$500 $M_\star>10^{10}~M_\odot$ galaxies in the redshift range $0.5<z<2.3$. We compare the resulting rest-frame 1.5-2$μ$m half-light radii ($R_{\rm{NIR}}$) with stellar half-mass radii (\rmass) derived with multi-color light profiles from CANDELS HST imaging. In general agreement with previous work, we find that… ▽ More We use CEERS JWST/NIRCam imaging to measure rest-frame near-IR light profiles of $>$500 $M_\star>10^{10}~M_\odot$ galaxies in the redshift range $0.5<z<2.3$. We compare the resulting rest-frame 1.5-2$μ$m half-light radii ($R_{\rm{NIR}}$) with stellar half-mass radii (\rmass) derived with multi-color light profiles from CANDELS HST imaging. In general agreement with previous work, we find that $R_{\rm{NIR}}$ and \rmass~are up to 40\%~smaller than the rest-frame optical half-light radius $R_{\rm{opt}}$. The agreement between $R_{\rm{NIR}}$ and \rmass~is excellent, with negligible systematic offset ($<$0.03 dex) up to $z=2$ for quiescent galaxies and up to $z=1.5$ for star-forming galaxies. We also deproject the profiles to estimate \rmassd, the radius of a sphere containing 50\% of the stellar mass. We present the $R-M_\star$ distribution of galaxies at $0.5<z<1.5$, comparing $R_{\rm{opt}}$, \rmass~and \rmassd. The slope is significantly flatter for \rmass~and \rmassd~ compared to $R_{\rm{opt}}$, mostly due to downward shifts in size for massive star-forming galaxies, while \rmass~and \rmassd~do not show markedly different trends. Finally, we show rapid size evolution ($R\propto (1+z)^{-1.7\pm0.1}$) for massive ($M_\star>10^{11}~M_\odot$) quiescent galaxies between $z=0.5$ and $z=2.3$, again comparing $R_{\rm{opt}}$, \rmass~and \rmassd. We conclude that the main tenets of the size evolution narrative established over the past 20 years, based on rest-frame optical light profile analysis, still hold in the era of JWST/NIRCam observations in the rest-frame near-IR. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: Submitted to ApJ. Comments welcome

RomAndromeda: The Roman Survey of the Andromeda Halo

Authors: Arjun Dey, Joan Najita, Carrie Filion, Jiwon Jesse Han, Sarah Pearson, Rosemary Wyse, Adrien C. R. Thob, Borja Anguiano, Miranda Apfel, Magda Arnaboldi, Eric F. Bell, Leandro Beraldo e Silva, Gurtina Besla, Aparajito Bhattacharya, Souradeep Bhattacharya, Vedant Chandra, Yumi Choi, Michelle L. M. Collins, Emily C. Cunningham, Julianne J. Dalcanton, Ivanna Escala, Hayden R. Foote, Annette M. N. Ferguson, Benjamin J. Gibson, Oleg Y. Gnedin , et al. (28 additional authors not shown)

Abstract: As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the ha… ▽ More As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$σ$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

NANCY: Next-generation All-sky Near-infrared Community surveY

Authors: Jiwon Jesse Han, Arjun Dey, Adrian M. Price-Whelan, Joan Najita, Edward F. Schlafly, Andrew Saydjari, Risa H. Wechsler, Ana Bonaca, David J Schlegel, Charlie Conroy, Anand Raichoor, Alex Drlica-Wagner, Juna A. Kollmeier, Sergey E. Koposov, Gurtina Besla, Hans-Walter Rix, Alyssa Goodman, Douglas Finkbeiner, Abhijeet Anand, Matthew Ashby, Benedict Bahr-Kalus, Rachel Beaton, Jayashree Behera, Eric F. Bell, Eric C Bellm , et al. (184 additional authors not shown)

Abstract: The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GAL… ▽ More The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

CEERS: MIRI deciphers the spatial distribution of dust-obscured star formation in galaxies at $0.1<z<2.5$

Authors: Benjamin Magnelli, Carlos Gómez-Guijarro, David Elbaz, Emanuele Daddi, Casey Papovich, Lu Shen, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Véronique Buat, Luca Costantin, Mark Dickinson, Steven L. Finkelstein, Jonathan P. Gardner, Eric F. Jiménez-Andrade, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Yipeng Lyu, Pablo G. Pérez-González, Nor Pirzkal, Sandro Tacchella, Alexander de la Vega, Stijn Wuyts, Guang Yang, L. Y. Aaron Yung , et al. (1 additional authors not shown)

Abstract: [Abridged] We combined HST images from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey with JWST images from the Cosmic Evolution Early Release Science (CEERS) survey to measure the stellar and dust-obscured star formation distributions of a mass-complete ($>10^{10}M_\odot$) sample of 69 star-forming galaxies (SFGs) at $0.1<z<2.5$. Rest-mid-infrared (rest-MIR) morphologies (size… ▽ More [Abridged] We combined HST images from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey with JWST images from the Cosmic Evolution Early Release Science (CEERS) survey to measure the stellar and dust-obscured star formation distributions of a mass-complete ($>10^{10}M_\odot$) sample of 69 star-forming galaxies (SFGs) at $0.1<z<2.5$. Rest-mid-infrared (rest-MIR) morphologies (sizes and Sérsic indices) were determined using their sharpest Mid-InfraRed Instrument (MIRI) images dominated by dust emission. Rest-MIR Sérsic indices were only measured for the brightest MIRI sources ($S/N>75$; 35 galaxies). At lower $S/N$, simulations show that simultaneous measurements of the size and Sérsic index become unreliable. We extended our study to fainter sources ($S/N>10$; 69 galaxies) by fixing their Sérsic index to unity. The Sérsic index of bright galaxies ($S/N>75$) has a median value of 0.7, which, together with their axis ratio distribution, suggests a disk-like morphology in the rest-MIR. Galaxies above the main sequence (MS; i.e., starbursts) have rest-MIR sizes that are a factor 2 smaller than their rest-optical sizes. The median rest-optical to rest-MIR size ratio of MS galaxies increases with stellar mass, from 1.1 at $10^{9.8}M_\odot$ to 1.6 at $10^{11}M_\odot$. This mass-dependent trend resembles the one found in the literature between the rest-optical and rest-near-infrared sizes of SFGs, suggesting that it is due to radial color gradients affecting rest-optical sizes and that the sizes of the stellar and star-forming components of SFGs are, on average, consistent at all masses. There is, however, a small population of SFGs (15%) with a compact star-forming component embedded in a larger stellar structure. This could be the missing link between galaxies with an extended stellar component and those with a compact stellar component, the so-called blue nuggets. △ Less

Submitted 16 October, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: Published in A&A. 18 pages, 14 figures

Journal ref: A&A 678, A83 (2023)

Origins of the Evil Eye: M64's Stellar Halo Reveals the Recent Accretion of an SMC-mass Satellite

Authors: Adam Smercina, Eric F. Bell, Paul A. Price, Jeremy Bailin, Julianne J. Dalcanton, Roelof S. de Jong, Richard D'Souza, Katya Gozman, In Sung Jang, Antonela Monachesi, David Nidever, Colin T. Slater

Abstract: M64, often called the "Evil Eye" galaxy, is unique among local galaxies. Beyond its dramatic, dusty nucleus, it also hosts an outer gas disk that counter-rotates relative to its stars. The mass of this outer disk is comparable to the gas content of the Small Magellanic Cloud (SMC), prompting the idea that it was likely accreted in a recent minor merger. Yet, detailed follow-up studies of M64's out… ▽ More M64, often called the "Evil Eye" galaxy, is unique among local galaxies. Beyond its dramatic, dusty nucleus, it also hosts an outer gas disk that counter-rotates relative to its stars. The mass of this outer disk is comparable to the gas content of the Small Magellanic Cloud (SMC), prompting the idea that it was likely accreted in a recent minor merger. Yet, detailed follow-up studies of M64's outer disk have shown no evidence of such an event, leading to other interpretations, such as a "flyby" interaction with the distant diffuse satellite Coma P. We present Subaru Hyper Suprime-Cam observations of M64's stellar halo, which resolve its stellar populations and reveal a spectacular radial shell feature, oriented $\sim$30$^{\circ}$ relative to the major axis and along the rotation axis of the outer gas disk. The shell is $\sim$45 kpc southeast of M64, while a similar but more diffuse plume to the northwest extends to $>$100 kpc. We estimate a stellar mass and metallicity for the southern shell of $M_{\star} {=} 1.80~{\pm}~0.54{\times}10^8~M_{\odot}$ and [M/H] $=$ $-$1.0, respectively, and a similar mass of $1.42~{\pm}~0.71{\times}10^8 M_{\odot}$ for the northern plume. Taking into account the accreted material in M64's inner disk, we estimate a total stellar mass for the progenitor satellite of $M_{\rm \star,prog}~{\simeq}~5{\times}10^8~M_{\odot}$. These results suggest that M64 is in the final stages of a minor merger with a gas-rich satellite strikingly similar to the SMC, in which M64's accreted counter-rotating gas originated, and which is responsible for the formation of its dusty inner star-forming disk. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 12 pages, 5 figures, accepted for publication in ApJ Letters

Clusters, Clouds, and Correlations: Relating Young Clusters to Giant Molecular Clouds in M33 and M31

Authors: Joshua Peltonen, Erik Rosolowsky, L. Clifton Johnson, Anil C. Seth, Julianne Dalcanton, Eric F. Bell, Jonathan Braine, Eric W. Koch, Margaret Lazzarini, Adam K. Leroy, Evan D. Skillman, Adam Smercina, Tobin Wainer, Benjamin F. Williams

Abstract: We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249… ▽ More We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249 clusters to 251 GMCs measured from a CARMA $^{12}$CO(1-0) survey with 20 pc resolution. Through two-point correlation analysis, we find that young clusters have a high probability of being near other young clusters, but correlation between GMCs is suppressed by the cloud identification algorithm. By comparing the positions, we find that younger clusters are closer to GMCs than older clusters. Through cross-correlation analysis of the M33 cluster data, we find that clusters are statistically associated when they are $\leq$10 Myr old. Utilizing the high precision ages of the clusters, we find that clusters older than $\approx 18$ Myr are uncorrelated with the molecular ISM. Using the spatial coincidence of the youngest clusters and GMCs in M33, we estimate that clusters spend $\approx$4-6 Myr inside their parent GMC. Through similar analysis, we find that the GMCs in M33 have a total lifetime of $\approx 11$-15 Myr. We also develop a drift model and show that the above correlations can be explained if the clusters in M33 have a 5-10 km s$^{-1}$ velocity dispersion relative to the molecular ISM. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 14 pages, 13 figures, 1 tables, accepted for publication at MNRAS

The Gas-Phase Mass--Metallicity Relation for Massive Galaxies at $z\sim0.7$ with the LEGA-C Survey

Authors: Zach J. Lewis, Brett H. Andrews, Rachel Bezanson, Michael Maseda, Eric F. Bell, Romeel Davé, Francesco D'Eugenio, Marijn Franx, Anna Gallazzi, Anna de Graaff, Yasha Kaushal, Angelos Nersesian, Jeffrey A. Newman, Arjen van der Wel, Po-Feng Wu

Abstract: The massive end of the gas-phase mass--metallicity relation (MZR) is a sensitive probe of active galactic nuclei (AGN) feedback that is a crucial but highly uncertain component of galaxy evolution models. In this paper, we extend the $z\sim0.7$ MZR by $\sim$0.5 dex up to log$(M_\star/\textrm{M}_\odot)\sim11.1$. We use extremely deep VLT VIMOS spectra from the Large Early Galaxy Astrophysics Census… ▽ More The massive end of the gas-phase mass--metallicity relation (MZR) is a sensitive probe of active galactic nuclei (AGN) feedback that is a crucial but highly uncertain component of galaxy evolution models. In this paper, we extend the $z\sim0.7$ MZR by $\sim$0.5 dex up to log$(M_\star/\textrm{M}_\odot)\sim11.1$. We use extremely deep VLT VIMOS spectra from the Large Early Galaxy Astrophysics Census (LEGA-C) survey to measure metallicities for 145 galaxies. The LEGA-C MZR matches the normalization of the $z\sim0.8$ DEEP2 MZR where they overlap, so we combine the two to create an MZR spanning from 9.3 to 11.1 log$(M_\star/\textrm{M}_\odot)$. The LEGA-C+DEEP2 MZR at $z\sim0.7$ is offset to slightly lower metallicities (0.05-0.13 dex) than the $z\sim0$ MZR, but it otherwise mirrors the established power law rise at low/intermediate stellar masses and asymptotic flattening at high stellar masses. We compare the LEGA-C+DEEP2 MZR to the MZR from two cosmological simulations (IllustrisTNG and SIMBA), which predict qualitatively different metallicity trends for high-mass galaxies. This comparison highlights that our extended MZR provides a crucial observational constraint for galaxy evolution models in a mass regime where the MZR is very sensitive to choices about the implementation of AGN feedback. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: 10 pages, 4 figures, 1 table

The Panchromatic Hubble Andromeda Treasury XX: The Disk of M31 is Thick

Authors: Julianne J. Dalcanton, Eric F. Bell, Yumi Choi, Andrew E. Dolphin, Morgan Fouesneau, Léo Girardi, David W. Hogg, Anil C. Seth, Benjamin F. Williams

Abstract: We present a new approach to measuring the thickness of a partially face-on stellar disk, using dust geometry. In a moderately-inclined disk galaxy, the fraction of reddened stars is expected to be 50% everywhere, assuming that dust lies in a thin midplane. In a thickened disk, however, a wide range of radii project onto the line of sight. Assuming stellar density declines with radius, this geomet… ▽ More We present a new approach to measuring the thickness of a partially face-on stellar disk, using dust geometry. In a moderately-inclined disk galaxy, the fraction of reddened stars is expected to be 50% everywhere, assuming that dust lies in a thin midplane. In a thickened disk, however, a wide range of radii project onto the line of sight. Assuming stellar density declines with radius, this geometrical projection leads to differences in the numbers of stars on the near and far sides of the thin dust layer. The fraction of reddened stars will thus differ from the 50% prediction, with a deviation that becomes larger for puffier disks. We map the fraction of reddened red giant branch (RGB) stars across M31, which shows prominent dust lanes on only one side of the major axis. The fraction of reddened stars varies systematically from 20% to 80%, which requires that these stars have an exponential scale height h_z that is 0.14+/-0.015 times the exponential scale length (h_r~5.5kpc). M31's RGB stars must therefore have h_z=770+/-80pc, which is far thicker than the Milky Way's thin disk, but comparable to its thick disk. The lack of a significant thin disk in M31 is unexpected, but consistent with its interaction history and high disk velocity dispersion. We suggest that asymmetric reddening be used as a generic criteria for identifying ``thick disk'' dominated systems, and discuss prospects for future 3-dimensional tomographic mapping of the gas and stars in M31. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 22 pages. Accepted to the Astrophysical Journal

Saying Hallo to M94's Stellar Halo: Investigating the Accretion History of the Largest Pseudobulge Host in the Local Universe

Authors: Katya Gozman, Eric F. Bell, Adam Smercina, Paul Price, Jeremy Bailin, Roelof S. de Jong, Richard D'Souza, In Sung Jang, Antonela Monachesi, Colin Slater

Abstract: It is not yet settled how the combination of secular processes and merging gives rise to the bulges and pseudobulges of galaxies. The nearby ($D\sim$ 4.2 Mpc) disk galaxy M94 (NGC 4736) has the largest pseudobulge in the local universe, and offers a unique opportunity for investigating the role of merging in the formation of its pseudobulge. We present a first ever look at M94's stellar halo, whic… ▽ More It is not yet settled how the combination of secular processes and merging gives rise to the bulges and pseudobulges of galaxies. The nearby ($D\sim$ 4.2 Mpc) disk galaxy M94 (NGC 4736) has the largest pseudobulge in the local universe, and offers a unique opportunity for investigating the role of merging in the formation of its pseudobulge. We present a first ever look at M94's stellar halo, which we expect to contain a fossil record of M94's past mergers. Using Subaru's Hyper Suprime-Cam, we resolve and identify red giant branch (RGB) stars in M94's halo, finding two distinct populations. After correcting for completeness through artificial star tests, we can measure the radial profile of each RGB population. The metal-rich RGB stars show an unbroken exponential profile to a radius of 30 kpc that is a clear continuation of M94's outer disk. M94's metal poor stellar halo is detectable over a wider area and clearly separates from its metal-rich disk. By integrating the halo density profile, we infer a total accreted stellar mass of $\sim 2.8 \times 10^8 M_\odot$, with a median metallicity of [M/H] $=-$1.4. This indicates that M94's most-massive past merger was with a galaxy similar to, or less massive than, the Small Magellanic Cloud. Few nearby galaxies have had such a low-mass dominant merger; therefore we suggest that M94's pseudobulge was not significantly impacted by merging. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 19 pages, 12 figures, 1 table; Published in ApJ on April 14, 2023

CEERS Key Paper VI: JWST/MIRI Uncovers a Large Population of Obscured AGN at High Redshifts

Authors: G. Yang, K. I. Caputi, C. Papovich, P. Arrabal Haro, M. B. Bagley, P. Behroozi, E. F. Bell, L. Bisigello, V. Buat, D. Burgarella, Y. Cheng, N. J. Cleri, R. Dave, M. Dickinson, D. Elbaz, H. C. Ferguson, S. L. Finkelstein, N. A. Grogin, N. P. Hathi, M. Hirschmann, B. W. Holwerda, M. Huertas-Company, T. Hutchison, E. Iani, J. S. Kartaltepe , et al. (13 additional authors not shown)

Abstract: Mid-infrared observations are powerful in identifying heavily obscured Active Galactic Nuclei (AGN) which have weak emission in other wavelengths. Data from the Mid-Infrared Instrument (MIRI) onboard JWST provides an excellent opportunity to perform such studies. We take advantage of the MIRI imaging data from the Cosmic Evolution Early Release Science Survey (CEERS) to investigate the AGN populat… ▽ More Mid-infrared observations are powerful in identifying heavily obscured Active Galactic Nuclei (AGN) which have weak emission in other wavelengths. Data from the Mid-Infrared Instrument (MIRI) onboard JWST provides an excellent opportunity to perform such studies. We take advantage of the MIRI imaging data from the Cosmic Evolution Early Release Science Survey (CEERS) to investigate the AGN population in the distant universe. We estimate the source properties of MIRI-selected objects by utilizing spectral energy distribution (SED) modelling, and classify them into star-forming galaxies (SF), SF-AGN mixed objects, and AGN. The source numbers of these types are 418, 111, and 31, respectively, from 4 MIRI pointings covering $\sim 9$ arcmin$^2$. The sample spans a redshift range of $\approx 0$--5. We derive the median SEDs for all three source types, respectively, and publicly release them. The median MIRI SED of AGN is similar to the typical SEDs of hot dust-obscured galaxies and Seyfert 2s, for which the mid-IR SEDs are dominantly from AGN-heated hot dust. Based on our SED-fit results, we estimate the black-hole accretion density (BHAD; i.e., total BH growth rate per comoving volume) as a function of redshift. At $z<3$, the resulting BHAD agrees with the X-ray measurements in general. At $z>3$, we identify a total of 27 AGN and SF-AGN mixed objects, leading to that our high-$z$ BHAD is substantially higher than the X-ray results ($\sim 0.5$ dex at $z \approx 3$--5). This difference indicates MIRI can identify a large population of heavily obscured AGN missed by X-ray surveys at high redshifts. △ Less

Submitted 15 May, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: ApJL in press, as part of the CEERS Focus Issue. 13 pages and 7 figures (median SEDs in Fig. 5 are available upon request)

Evolution in the orbital structure of quiescent galaxies from MAGPI, LEGA-C and SAMI surveys: direct evidence for merger-driven growth over the last 7 Gy

Authors: Francesco D'Eugenio, Arjen van der Wel, Joanna M. Piotrowska, Rachel Bezanson, Edward N. Taylor, Jesse van de Sande, William M. Baker, Eric F. Bell, Sabine Bellstedt, Joss Bland-Hawthorn, Asa F. L. Bluck, Sarah Brough, Julia J. Bryant, Matthew Colless, Luca Cortese, Scott M. Croom, Caro Derkenne, Pieter van Dokkum, Deanne Fisher, Caroline Foster, Anna Gallazzi, Anna de Graaff, Brent Groves, Josha van Houdt, Claudia del P. Lagos , et al. (15 additional authors not shown)

Abstract: We present the first study of spatially integrated higher-order stellar kinematics over cosmic time. We use deep rest-frame optical spectroscopy of quiescent galaxies at redshifts z=0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. Conservatively using a redshift-in… ▽ More We present the first study of spatially integrated higher-order stellar kinematics over cosmic time. We use deep rest-frame optical spectroscopy of quiescent galaxies at redshifts z=0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. Conservatively using a redshift-independent cut in stellar mass ($M_\star = 10^{11}\,{\rm M}_\odot$), and matching the stellar-mass distributions of our samples, we find 7 $σ$ evidence of $h_4$ increasing with cosmic time, from a median value of 0.019$\pm$0.002 at z=0.8 to 0.059$\pm$0.004 at z=0.06. Alternatively, we use a physically motivated sample selection, based on the mass distribution of the progenitors of local quiescent galaxies as inferred from numerical simulations; in this case, we find 10 $σ$ evidence. This evolution suggests that, over the last 7 Gyr, there has been a gradual decrease in the rotation-to-dispersion ratio and an increase in the radial anisotropy of the stellar velocity distribution, qualitatively consistent with accretion of gas-poor satellites. These findings demonstrate that massive galaxies continue to accrete mass and increase their dispersion support after becoming quiescent. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: 19 pages, 9 figures Accepted for publication in MNRAS

Different higher-order kinematics between star-forming and quiescent galaxies based on the SAMI, MAGPI and LEGA-C surveys

Authors: Francesco D'Eugenio, Arjen van der Wel, Caro Derkenne, Josha van Houdt, Rachel Bezanson, Edward N. Taylor, Jesse van de Sande, William M. Baker, Eric F. Bell, Joss Bland-Hawthorn, Asa F. L. Bluck, Sarah Brough, Julia J. Bryant, Matthew Colless, Luca Cortese, Scott M. Croom, Pieter van Dokkum, Deanne Fisher, Caroline Foster, Amelia Fraser-McKelvie, Anna Gallazzi, Anna de Graaff, Brent Groves, Claudia del P. Lagos, Tobias J. Looser , et al. (16 additional authors not shown)

Abstract: We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass $M_\star > 10^{10.5} {\rm M}_\odot$) at redshifts z = 0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys, to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the lat… ▽ More We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass $M_\star > 10^{10.5} {\rm M}_\odot$) at redshifts z = 0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys, to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. We find that at all redshifts where we have large enough samples, $h_4$ anti-correlates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anti-correlation with rotation-to-dispersion ratio, we also find a correlation between $h_4$ and $M_\star$, potentially connected to the assembly history of galaxies. In contrast, after controlling for mass, we find no evidence of independent correlation between $h_4$ and aperture velocity dispersion or galaxy size. These results hold for both star-forming and quiescent galaxies. For quiescent galaxies, $h_4$ also correlates with projected shape, even after controlling for the rotation-to-dispersion ratio. At any given redshift, star-forming galaxies have lower $h_4$ compared to quiescent galaxies, highlighting the link between kinematic structure and star-forming activity. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 26 pages, 15 figures Accepted for publication in MNRAS

On the nature of disks at high redshift seen by JWST/CEERS with contrastive learning and cosmological simulations

Authors: J. Vega-Ferrero, M. Huertas-Company, L. Costantin, P. G. Pérez-González, R. Sarmiento, J. S. Kartaltepe, A. Pillepich, M. B. Bagley, S. L. Finkelstein, E. J. McGrath, J. H. Knapen, P. Arrabal Haro, E. F. Bell, F. Buitrago, A. Calabrò, A. Dekel, M. Dickinson, H. Domínguez Sánchez, D. Elbaz, H. C. Ferguson, M. Giavalisco, B. W. Holwerda, D. D. Kocesvski, A. M. Koekemoer, V. Pandya , et al. (4 additional authors not shown)

Abstract: Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to… ▽ More Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to be robust to noise and to correlate well with the physical properties of the simulated galaxies, including their 3D structure. We apply the method simultaneously to F200W and F356W galaxy images of a mass-complete sample ($M_*/M_\odot>10^9$) at $ 3 \leq z \leq 6$ from the first JWST/NIRCam CEERS data release. We find that the simulated and observed galaxies do not exactly populate the same manifold in the representation space from contrastive learning. We also find that half the galaxies classified as disks -- either CNN-based or visually -- populate a similar region of the representation space as TNG50 galaxies with low stellar specific angular momentum and non-oblate structure. Although our data-driven study does not allow us to firmly conclude on the true nature of these galaxies, it suggests that the disk fraction at $z \geq 3$ remains uncertain and possibly overestimated by traditional supervised classifications. Deeper imaging and spectroscopic follow-ups as well as comparisons with other simulations will help to unambiguously determine the true nature of these galaxies, and establish more robust constraints on the emergence of disks at very high redshift. △ Less

Submitted 25 October, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

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

Hidden Little Monsters: Spectroscopic Identification of Low-Mass, Broad-Line AGN at $z>5$ with CEERS

Authors: Dale D. Kocevski, Masafusa Onoue, Kohei Inayoshi, Jonathan R. Trump, Pablo Arrabal Haro, Andrea Grazian, Mark Dickinson, Steven L. Finkelstein, Jeyhan S. Kartaltepe, Michaela Hirschmann, Seiji Fujimoto, Stephanie Juneau, Ricardo O. Amorin, Micaela B. Bagley, Guillermo Barro, Eric F. Bell, Laura Bisigello, Antonello Calabro, Nikko J. Cleri, M. C. Cooper, Xuheng Ding, Norman A. Grogin, Luis C. Ho, Akio K. Inoue, Linhua Jiang , et al. (12 additional authors not shown)

Abstract: We report on the discovery of two low-luminosity, broad-line AGN at $z>5$ identified using JWST NIRSpec spectroscopy from the CEERS Survey. We detect broad H$α$ emission from both sources, with FWHM of $2038\pm286$ and $1807\pm207$ km s$^{-1}$, resulting in black hole (BH) masses that are 1-2 dex below that of existing samples of luminous quasars at $z>5$. The first source, CEERS 1670 at… ▽ More We report on the discovery of two low-luminosity, broad-line AGN at $z>5$ identified using JWST NIRSpec spectroscopy from the CEERS Survey. We detect broad H$α$ emission from both sources, with FWHM of $2038\pm286$ and $1807\pm207$ km s$^{-1}$, resulting in black hole (BH) masses that are 1-2 dex below that of existing samples of luminous quasars at $z>5$. The first source, CEERS 1670 at $z=5.242$, is 2-3 dex fainter than known quasars at similar redshifts and was previously identified as a candidate low-luminosity AGN based on its rest-frame optical SED. We measure a BH mass of $M_{\rm BH}=1.3\pm0.4\times 10^{7}~M_{\odot}$, confirming that this AGN is powered by the least-massive BH known in the universe at the end of cosmic reionization. The second source, CEERS 3210 at $z=5.624$, is inferred to be a heavily obscured, broad-line AGN caught in a transition phase between a dust-obscured starburst and an unobscured quasar. We estimate its BH mass to be $M_{\rm BH}\simeq 0.9-4.7 \times 10^{7}~M_{\odot}$, depending on the level of dust obscuration assumed. We derive host stellar masses, $M_\star$, allowing us to place constraints on the BH-galaxy mass relationship in the lowest mass range yet probed in the early universe. The $M_{\rm BH}/M_\star$ ratio for CEERS 1670, in particular, is consistent with or higher than the empirical relationship seen in massive galaxies at $z=0$. We examine the emission-line ratios of both sources and find that their location on the BPT and OHNO diagrams is consistent with model predictions for low-metallicity AGN with $Z/Z_\odot \simeq 0.2-0.4$. The spectroscopic identification of low-luminosity, broad-line AGN at $z>5$ with $M_{\rm BH}\simeq 10^{7}~M_{\odot}$ demonstrates the capability of JWST to push BH masses closer to the range predicted for the BH seed population and provides a unique opportunity to study the early stages of BH-galaxy assembly. △ Less

Submitted 31 January, 2023; originally announced February 2023.

Comments: 16 pages, 8 figures, submitted to ApJL

CEERS: Spatially Resolved UV and mid-IR Star Formation in Galaxies at 0.2 < z < 2.5: The Picture from the Hubble and James Webb Space Telescopes

Authors: Lu Shen, Casey Papovich, Guang Yang, Jasleen Matharu, Xin Wang, Benjamin Magnelli, David Elbaz, Shardha Jogee, Anahita Alavi, Pablo Arrabal Haro, Bren E. Backhaus, Micaela B. Bagley, Eric F. Bell, Laura Bisigello, Antonello Calabrò, M. C. Cooper, Luca Costantin, Emanuele Daddi, Mark Dickinson, Steven L. Finkelstein, Seiji Fujimoto, Mauro Giavalisco, Norman A. Grogin, Yuchen Guo, Benne W. Holwerda , et al. (16 additional authors not shown)

Abstract: We present the mid-IR (MIR) morphologies for 64 star-forming galaxies at $0.2<z<2.5$ with stellar mass $\rm{M_*>10^{9}~M_\odot}$ using JWST MIRI observations from the Cosmic Evolution Early Release Science survey (CEERS). The MIRI bands span the MIR (7.7--21~$μ$m), enabling us to measure the effective radii ($R_{\rm{eff}}$) and Sérsic indexes of these SFGs at rest-frame 6.2 and 7.7 $μ$m, which con… ▽ More We present the mid-IR (MIR) morphologies for 64 star-forming galaxies at $0.2<z<2.5$ with stellar mass $\rm{M_*>10^{9}~M_\odot}$ using JWST MIRI observations from the Cosmic Evolution Early Release Science survey (CEERS). The MIRI bands span the MIR (7.7--21~$μ$m), enabling us to measure the effective radii ($R_{\rm{eff}}$) and Sérsic indexes of these SFGs at rest-frame 6.2 and 7.7 $μ$m, which contains strong emission from Polycyclic aromatic hydrocarbon (PAH) features, a well-established tracer of star formation in galaxies. We define a ``PAH-band'' as the MIRI bandpass that contains these features at the redshift of the galaxy. We then compare the galaxy morphologies in the PAH-bands to those in rest-frame Near-UV (NUV) using HST ACS/F435W or ACS/F606W and optical/near-IR using HST WFC3/F160W imaging from UVCANDELS and CANDELS, where the NUV-band and F160W trace the profile of (unobscured) massive stars and the stellar continuum, respectively. The $R_{\rm{eff}}$ of galaxies in the PAH-band are slightly smaller ($\sim$10\%) than those in F160W for galaxies with $\rm{M_*\gtrsim10^{9.5}~M_\odot}$ at $z\leq1.2$, but the PAH-band and F160W have a similar fractions of light within 1 kpc. In contrast, the $R_{\rm{eff}}$ of galaxies in the NUV-band are larger, with lower fractions of light within 1 kpc compared to F160W for galaxies at $z\leq1.2$. Using the MIRI data to estimate the $\rm{SFR_{\rm{IR}}}$ surface density, we find the correlation between the $\rm{SFR_{\rm{IR}}}$ surface density and stellar mass has a steeper slope than that of the $\rm{SFR_{\rm{UV}}}$ surface density and stellar mass, suggesting more massive galaxies having increasing amounts of obscured fraction of star formation in their inner regions. This paper demonstrates how the high-angular resolution data from JWST/MIRI can reveal new information about the morphology of obscured-star formation. △ Less

Submitted 2 April, 2023; v1 submitted 13 January, 2023; originally announced January 2023.

Comments: 28 pages, 11 figures, Accepted by ApJ

CEERS Key Paper IV: Galaxies at $4 < z < 9$ are Bluer than They Appear -- Characterizing Galaxy Stellar Populations from Rest-Frame $\sim 1$ micron Imaging

Authors: Casey Papovich, Justin Cole, Guang Yang, Steven L. Finkelstein, Guillermo Barro, Véronique Buat, Denis Burgarella, Pablo G. Pérez-González, Paola Santini, Lise-Marie Seillé, Lu Shen, Pablo Arrabal Haro, Micaela B. Bagley, Eric F. Bell, Laura Bisigello, Antonello Calabrò, Caitlin M. Casey, Marco Castellano, Katherine Chworowsky, Nikko J. Cleri, M. C. Cooper, Luca Costantin, Mark Dickinson, Henry C. Ferguson, Adriano Fontana , et al. (24 additional authors not shown)

Abstract: We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4<z<9$ using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 $μ$m data extend the coverage of the rest-frame… ▽ More We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4<z<9$ using imaging data from the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) combined with data from the Hubble Space Telescope and the Spitzer Space Telescope. The JWST/MIRI 5.6 and 7.7 $μ$m data extend the coverage of the rest-frame spectral-energy distribution (SED) to nearly 1 micron for galaxies in this redshift range. By modeling the galaxies' SEDs the MIRI data show that the galaxies have, on average, rest-frame UV (1600 Å) $-$ $I$-band colors 0.4 mag bluer than derived when using photometry that lacks MIRI. Therefore, the galaxies have lower (stellar)-mass-to-light ratios. The MIRI data reduce the stellar masses by $\langle Δ\log M_\ast\rangle=0.25$ dex at $4<z<6$ (a factor of 1.8) and 0.37 dex at $6<z<9$ (a factor of 2.3). This also reduces the star-formation rates (SFRs) by $\langle Δ\log\mathrm{SFR} \rangle=0.14$ dex at $4<z<6$ and 0.27 dex at $6<z<9$. The MIRI data also improve constraints on the allowable stellar mass formed in early star-formation. We model this using a star-formation history that includes both a "burst' at $z_f=100$ and a slowly varying ("delayed-$τ$") model. The MIRI data reduce the allowable stellar mass by 0.6 dex at $4<z< 6$ and by $\approx$1 dex at $6<z<9$. Applying these results globally, this reduces the cosmic stellar-mass density by an order of magnitude in the early universe ($z\approx9$). Therefore, observations of rest-frame $\gtrsim$1 $μ$m are paramount for constraining the stellar-mass build-up in galaxies at very high-redshifts. △ Less

Submitted 25 March, 2023; v1 submitted 30 December, 2022; originally announced January 2023.

Comments: Updated with accepted ApJ version. Part of the CEERS Focus Issue. 27 pages, many figures (4 Figure Sets, available upon reasonable request)

Rapid radio brightening of GRB 210702A

Authors: G. E. Anderson, T. D. Russell, H. M. Fausey, A. J. van der Horst, P. J. Hancock, A. Bahramian, M. E. Bell, J. C. A. Miller-Jones, G. Rowell, M. W. Sammons, R. A. M. J. Wijers, T. J. Galvin, A. J. Goodwin, R. Konno, A. Rowlinson, S. D. Ryder, F. Schussler, S. J. Wagner, S. J. Zhu

Abstract: We observed the rapid radio brightening of GRB 210702A with the Australian Telescope Compact Array (ATCA) just 11hr post-burst, tracking early-time radio variability over a 5hr period on ~15min timescales at 9.0, 16.7, and 21.2GHz. A broken power-law fit to the 9.0GHz light curve showed that the 5hr flare peaked at a flux density of 0.4+/-0.1mJy at ~13hr post-burst with a steep rise and decline. T… ▽ More We observed the rapid radio brightening of GRB 210702A with the Australian Telescope Compact Array (ATCA) just 11hr post-burst, tracking early-time radio variability over a 5hr period on ~15min timescales at 9.0, 16.7, and 21.2GHz. A broken power-law fit to the 9.0GHz light curve showed that the 5hr flare peaked at a flux density of 0.4+/-0.1mJy at ~13hr post-burst with a steep rise and decline. The observed temporal and spectral evolution are not expected in the standard internal-external shock model, where forward and reverse shock radio emission evolves on much longer timescales. The early-time (<1day) optical and X-ray light curves from the Neil Gehrels Swift Observatory demonstrated typical afterglow forward shock behaviour, allowing us to use blast wave physics to determine a likely homogeneous circumburst medium and an emitting electron population power-law index of p=2.9+/-0.1. We suggest the early-time radio flare is likely due to weak interstellar scintillation (ISS), which boosted the radio afterglow emission above the ATCA sensitivity limit on minute timescales. Using weak ISS relations, we were able to place an upper limit on the size of the blast wave of $\leq6 \times 10^{16}$cm in the plane of the sky, which is consistent with the theoretical forward shock size prediction of $8\times10^{16}$cm for GRB 210702A at ~13h post-burst. This represents the earliest ISS size constraint on a GRB blast wave to date, demonstrating the importance of rapid (<1day) radio follow-up of GRBs using several-hour integrations to capture the early afterglow evolution, and to track scintillation over a broad frequency range. △ Less

Submitted 22 June, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: Submitted to MNRAS 13 pages, 6 figures, 2 tables

CEERS Key Paper V: A triality on the nature of HST-dark galaxies

Authors: Pablo G. Pérez-González, Guillermo Barro, Marianna Annunziatella, Luca Costantin, Ángela García-Argumánez, Elizabeth J. McGrath, Rosa M. Mérida, Jorge A. Zavala, Pablo Arrabal Haro, Micaela B. Bagley, Bren E. Backhaus, Peter Behroozi, Eric F. Bell, Laura Bisigello, Véronique Buat, Antonello Calabrò, Caitlin M. Casey, Nikko J. Cleri, Rosemary T. Coogan, M. C. Cooper, Asantha R. Cooray, Avishai Dekel, Mark Dickinson, David Elbaz, Henry C. Ferguson , et al. (30 additional authors not shown)

Abstract: The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometr… ▽ More The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometric redshifts in 2 dimensions and stellar populations properties in a pixel-by-pixel basis. We select 138 galaxies with F150W-F356W>1.5 mag, F356W<27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies at 2<z<6 with masses 9<log M/M_sun<11 and a variety of specific SFRs (<1 to >100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3<z<5, masses log M/M_sun~10 and post-starburst stellar mass-weighted ages (0.5-1 Gyr); and (3) 11% are strong young starbursts with indications of high-EW emission lines (typically, [OIII]+Hbeta) at 6<z<7 and log M/M_sun~9.5. The sample is dominated by disk-like galaxies with a remarkable compactness for XELG-z6 (effective radii smaller than 0.4 kpc). Large attenuations in SFGs, 2<A(V)<5 mag, are found within 1.5 times the effective radius, approximately 2 kpc, while QGs present A(V)~0.2 mag. Our SED-fitting technique reproduces the expected dust emission luminosities of IR-bright and sub-millimeter galaxies. This study implies high levels of star formation activity between z~20 and z~10, where virtually 100% of our galaxies had already formed 10^8 M_sun of their stellar content, 60% of them had assembled 10^9 M_sun, and 10% up to 10^10 M_sun (in situ or ex situ). (abridged) △ Less

Submitted 3 April, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Published in CEERS ApJL Focus Issue, ApJL 946, L16

Journal ref: ApJL 946, L16 https://iopscience.iop.org/article/10.3847/2041-8213/acb3a5

First Look at z > 1 Bars in the Rest-Frame Near-Infrared with JWST Early CEERS Imaging

Authors: Yuchen Guo, Shardha Jogee, Steven L. Finkelstein, Zilei Chen, Eden Wise, Micaela B. Bagley, Guillermo Barro, Stijn Wuyts, Dale D. Kocevski, Jeyhan S. Kartaltepe, Elizabeth J. McGrath, Henry C. Ferguson, Bahram Mobasher, Mauro Giavalisco, Ray A. Lucas, Jorge A. Zavala, Jennifer M. Lotz, Norman A. Grogin, Marc Huertas-Company, Jesús Vega-Ferrero, Nimish P. Hathi, Pablo Arrabal Haro, Mark Dickinson, Anton M. Koekemoer, Casey Papovich , et al. (23 additional authors not shown)

Abstract: Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of {\it{JWST}} CEERS NIRCam images to present the first quantitative identification and characterization of stell… ▽ More Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of {\it{JWST}} CEERS NIRCam images to present the first quantitative identification and characterization of stellar bars at $z>1$ based on rest-frame NIR F444W images of high resolution (~1.3 kpc at z ~ 1-3). We identify stellar bars in these images using quantitative criteria based on ellipse fits. For this pilot study, we present six examples of robustly identified bars at $z>1$ with spectroscopic redshifts, including the two highest redshift bars at ~2.136 and 2.312 quantitatively identified and characterized to date. The stellar bars at $z$ ~ 1.1-2.3 presented in our study have projected semi-major axes of ~2.9-4.3 kpc and projected ellipticities of ~0.41-0.53 in the rest-frame NIR. The barred host galaxies have stellar masses ~ $ 1 \times 10^{10}$ to $2 \times 10^{11}$ $M_{\odot}$, star formation rates of ~ 21-295 $M_{\odot}$ yr$^{-1}$, and several have potential nearby companions. Our finding of bars at $z$ ~1.1-2.3 demonstrates the early onset of such instabilities and supports simulations where bars form early in massive dynamically cold disks. It also suggests that if these bars at lookback times of 8-10 Gyr survive out to present epochs, bar-driven secular processes may operate over a long time and have a significant impact on some galaxies by z ~ 0. △ Less

Submitted 11 December, 2022; v1 submitted 16 October, 2022; originally announced October 2022.

Comments: 16 pages, 5 figures. Accepted for Publication in Astrophysical Journal Letters

The Velocity Dispersion Function for Massive Quiescent and Star-Forming Galaxies at 0.6 $<$ z $\leq$ 1.0

Authors: Lance Taylor, Rachel Bezanson, Arjen van der Wel, Alan Pearl, Eric F. Bell, Francesco D'Eugenio, Marijn Franx, Michael V. Maseda, Adam Muzzin, David Sobral, Caroline Straatman, Katherine E. Whitaker, Po-Feng Wu

Abstract: We present the first direct spectroscopic measurement of the stellar velocity dispersion function (VDF) for massive quiescent and star-forming galaxies at $0.6 < z \leq 1.0$. For this analysis we use individual measurements of stellar velocity dispersion from high-S/N spectra from the public Large Early Galaxy Astrophysics Census (LEGA-C) survey. We report a remarkable stability of the VDF for bot… ▽ More We present the first direct spectroscopic measurement of the stellar velocity dispersion function (VDF) for massive quiescent and star-forming galaxies at $0.6 < z \leq 1.0$. For this analysis we use individual measurements of stellar velocity dispersion from high-S/N spectra from the public Large Early Galaxy Astrophysics Census (LEGA-C) survey. We report a remarkable stability of the VDF for both quiescent and star-forming galaxies within this redshift range, though we note the presence of weak evolution in the number densities of star-forming galaxies. We compare both VDFs with previous direct and inferred measurements at local and intermediate redshifts, with the caveat that previous measurements of the VDF for star-forming galaxies are poorly constrained at all epochs. We emphasize that this work is the first to directly push to low-stellar velocity dispersion ($σ_\star > 100$ km s$^{-1}$) and extend to star-forming galaxies. We are largely consistent with the high-sigma tail measured from BOSS, and we find that the VDF remains constant from the median redshift of LEGA-C, $z\sim0.8$, to the present day. △ Less

Submitted 30 September, 2022; originally announced October 2022.

Comments: 11 pages, 3 figures, accepted for publication in ApJ

Ultrafaint Dwarf Galaxy Candidates in the M81 Group: Signatures of Group Accretion

Authors: Eric F. Bell, Adam Smercina, Paul A. Price, Richard D'Souza, Jeremy Bailin, Roelof S. de Jong, Katya Gozman, In Sung Jang, Antonela Monachesi, Oleg Y. Gnedin, Colin T. Slater

Abstract: The faint and ultrafaint dwarf galaxies in the Local Group form the observational bedrock upon which our understanding of small-scale cosmology rests. In order to understand whether this insight generalizes, it is imperative to use resolved-star techniques to discover similarly faint satellites in nearby galaxy groups. We describe our search for ultrafaint galaxies in the M81 group using deep grou… ▽ More The faint and ultrafaint dwarf galaxies in the Local Group form the observational bedrock upon which our understanding of small-scale cosmology rests. In order to understand whether this insight generalizes, it is imperative to use resolved-star techniques to discover similarly faint satellites in nearby galaxy groups. We describe our search for ultrafaint galaxies in the M81 group using deep ground-based resolved-star data sets from Subaru's Hyper Suprime-Cam. We present one new ultrafaint dwarf galaxy in the M81 group and identify five additional extremely low surface brightness candidate ultrafaint dwarfs that reach deep into the ultrafaint regime to $M_V \sim -6$ (similar to current limits for Andromeda satellites). These candidates' luminosities and sizes are similar to known Local Group dwarf galaxies Tucana B, Canes Venatici I, Hercules, and Boötes I. Most of these candidates are likely to be real, based on tests of our techniques on blank fields. Intriguingly, all of these candidates are spatially clustered around NGC 3077, which is itself an M81 group satellite in an advanced state of tidal disruption. This is somewhat surprising, as M81 itself and its largest satellite M82 are both substantially more massive than NGC 3077 and by virtue of their greater masses, would have been expected to host as many or more ultrafaint candidates. These results lend considerable support to the idea that satellites of satellites are an important contribution to the growth of satellite populations around Milky Way-mass galaxies. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: The Astrophysical Journal Letters; in press. 11 pages, 4 figures, 1 table

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

The Mass Scale of High-Redshift Galaxies: Virial Mass Estimates Calibrated with Stellar Dynamical Models from LEGA-C

Authors: Arjen van der Wel, Josha van Houdt, Rachel Bezanson, Marijn Franx, Francesco D'Eugenio, Caroline Straatman, Eric F. Bell, Adam Muzzin, David Sobral, Michael V. Maseda, Anna de Graaff, Bradford P. Holden

Abstract: Dynamical models for $673$ galaxies at $z=0.6-1.0$ with spatially resolved (long-slit) stellar kinematic data from LEGA-C are used to calibrate virial mass estimates defined as $M_{\rm{vir}}=K σ'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $σ'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a Sérsic p… ▽ More Dynamical models for $673$ galaxies at $z=0.6-1.0$ with spatially resolved (long-slit) stellar kinematic data from LEGA-C are used to calibrate virial mass estimates defined as $M_{\rm{vir}}=K σ'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $σ'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a Sérsic profile fit to HST imaging. The sample is representative for $M_{\star}>3\times10^{10}~M_{\odot}$ and includes all types of galaxies, irrespective of morphology and color. We demonstrate that using $R=R_{\rm{sma}}$~(the semi-major axis length of the ellipse that encloses 50\% of the light) in combination with an inclination correction on $σ'_{\star,\rm{int}}$~produces an unbiased $M_{\rm{vir}}$. We confirm the importance of projection effects on $σ'_{\star,\rm{int}}$ by showing the existence of a similar residual trend between virial mass estimates and inclination for the nearby early-type galaxies in the ATLAS$^{\rm{3D}}$~survey. Also, as previously shown, when using a Sérsic profile-based $R$ estimate, then a Sérsic index-dependent correction to account for non-homology in the radial profiles is required. With respect to analogous dynamical models for low-redshift galaxies from the ATLAS$^{\rm{3D}}$~survey we find a systematic offset of 0.1 dex in the calibrated virial constant for LEGA-C, which may be due to physical differences between the galaxy samples or an unknown systematic error. Either way, with our work we establish a common mass scale for galaxies across 8 Gyr of cosmic time with a systematic uncertainty of at most 0.1 dex. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: Published in ApJ

Journal ref: The Astrophysical Journal, 2022, 936, 1

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). IV. Star Cluster Catalog

Authors: L. Clifton Johnson, Tobin M. Wainer, Estephani E. TorresVillanueva, Anil C. Seth, Benjamin F. Williams, Meredith J. Durbin, Julianne J. Dalcanton, Daniel R. Weisz, Eric F. Bell, Puragra Guhathakurta, Evan Skillman, Adam Smercina

Abstract: We construct a catalog of star clusters from Hubble Space Telescope images of the inner disk of the Triangulum Galaxy (M33) using image classifications collected by the Local Group Cluster Search, a citizen science project hosted on the Zooniverse platform. We identify 1214 star clusters within the Hubble Space Telescope imaging footprint of the Panchromatic Hubble Andromeda Treasury: Triangulum E… ▽ More We construct a catalog of star clusters from Hubble Space Telescope images of the inner disk of the Triangulum Galaxy (M33) using image classifications collected by the Local Group Cluster Search, a citizen science project hosted on the Zooniverse platform. We identify 1214 star clusters within the Hubble Space Telescope imaging footprint of the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. Comparing this catalog to existing compilations in the literature, 68% of the clusters are newly identified. The final catalog includes multi-band aperture photometry and fits for cluster properties via integrated light SED fitting. The cluster catalog's 50% completeness limit is ~1500 solar masses at an age of 100 Myr, as derived from comprehensive synthetic cluster tests. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 24 pages, 15 figures, 9 tables, Accepted to ApJ

DESI Observations of the Andromeda Galaxy: Revealing the Immigration History of our Nearest Neighbor

Authors: Arjun Dey, Joan R. Najita, S. E. Koposov, J. Josephy-Zack, Gabriel Maxemin, Eric F. Bell, C. Poppett, E. Patel, L. Beraldo e Silva, A. Raichoor, D. Schlegel, D. Lang, A. Meisner, Adam D. Myers, J. Aguilar, S. Ahlen, C. Allende Prieto, D. Brooks, A. P. Cooper, K. S. Dawson, A. de la Macorra, P. Doel, A. Font-Ribera, Juan Garcia-Bellido, S. Gontcho A Gontcho , et al. (23 additional authors not shown)

Abstract: We present DESI observations of the inner halo of M31, which reveal the kinematics of a recent merger - a galactic immigration event - in exquisite detail. Of the 11,416 sources studied in 3.75 hour of on-sky exposure time, 7,438 are M31 sources with well measured radial velocities. The observations reveal intricate coherent kinematic structure in the positions and velocities of individual stars:… ▽ More We present DESI observations of the inner halo of M31, which reveal the kinematics of a recent merger - a galactic immigration event - in exquisite detail. Of the 11,416 sources studied in 3.75 hour of on-sky exposure time, 7,438 are M31 sources with well measured radial velocities. The observations reveal intricate coherent kinematic structure in the positions and velocities of individual stars: streams, wedges, and chevrons. While hints of coherent structures have been previously detected in M31, this is the first time they have been seen with such detail and clarity in a galaxy beyond the Milky Way. We find clear kinematic evidence for shell structures in the Giant Stellar Stream, the Northeast Shelf and Western Shelf regions. The kinematics are remarkably similar to the predictions of dynamical models constructed to explain the spatial morphology of the inner halo. The results are consistent with the interpretation that much of the substructure in the inner halo of M31 is produced by a single galactic immigration event 1 - 2 Gyr ago. Significant numbers of metal-rich stars ([Fe/H]$>-0.5$) are present in all of the detected substructures, suggesting that the immigrating galaxy had an extended star formation history. We also investigate the ability of the shells and Giant Stellar Stream to constrain the gravitational potential of M31, and estimate the mass within a projected radius of 125 kpc to be ${\rm log_{10}}\, M_{\rm NFW}(<125\,{\rm kpc})/M_\odot = 11.80_{-0.10}^{+0.12}$. The results herald a new era in our ability to study stars on a galactic scale and the immigration histories of galaxies. △ Less

Submitted 20 January, 2023; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: 45 pages, 22 figures, 8 tables; Astrophysical Journal in press; Data at https://zenodo.org/record/6977494

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 common origin for the Fundamental Plane of quiescent and star-forming galaxies in the EAGLE simulations

Authors: Anna de Graaff, Marijn Franx, Eric F. Bell, Rachel Bezanson, Matthieu Schaller, Joop Schaye, Arjen van der Wel

Abstract: We use the EAGLE cosmological simulations to perform a comprehensive and systematic analysis of the $z=0.1$ Fundamental Plane (FP), the tight relation between galaxy size, mass and velocity dispersion. We first measure the total mass and velocity dispersion (including both random and rotational motions) within the effective radius to show that simulated galaxies obey a total mass FP that is very c… ▽ More We use the EAGLE cosmological simulations to perform a comprehensive and systematic analysis of the $z=0.1$ Fundamental Plane (FP), the tight relation between galaxy size, mass and velocity dispersion. We first measure the total mass and velocity dispersion (including both random and rotational motions) within the effective radius to show that simulated galaxies obey a total mass FP that is very close to the virial relation ($<10\%$ deviation), indicating that the effects of non-homology are weak. When we instead use the stellar mass, we find a strong deviation from the virial plane, which is driven by variations in the dark matter content. The dark matter fraction is a smooth function of the size and stellar mass, and thereby sets the coefficients of the stellar mass FP without substantially increasing the scatter. Hence, both star-forming and quiescent galaxies obey the same FP, with equally low scatter (0.02 dex). We employ simulations with a variable stellar initial mass function (IMF) to show that IMF variations have a modest additional effect on this FP. Moreover, when we use luminosity-weighted mock observations of the size and spatially-integrated velocity dispersion, the inferred FP changes only slightly. However, the scatter increases significantly, due to the luminosity-weighting and line-of-sight projection of the velocity dispersions, and measurement uncertainties on the half-light radii. Importantly, we find significant differences between the simulated FP and observations, which likely reflects a systematic difference in the stellar mass distributions. Therefore, we suggest the stellar mass FP offers a simple test for cosmological simulations, requiring minimal post-processing of simulation data. △ Less

Submitted 8 November, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: accepted for publication in MNRAS; resolution of some figures reduced for arxiv version

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

walter: A Tool for Predicting Resolved Stellar Population Observations with Applications to the Roman Space Telescope

Authors: Lachlan Lancaster, Sarah Pearson, Benjamin F. Williams, Kathryn V. Johnston, Tjitske K. Starkenburg, Erin Kado-Fong, Anil C. Seth, Eric F. Bell

Abstract: Studies of resolved stellar populations in the Milky Way and nearby galaxies reveal an amazingly detailed and clear picture of galaxy evolution. Within the Local Group, the ability to probe the stellar populations of small and large galaxies opens up the possibility of exploring key questions such as the nature of dark matter, the detailed formation history of different galaxy components, and the… ▽ More Studies of resolved stellar populations in the Milky Way and nearby galaxies reveal an amazingly detailed and clear picture of galaxy evolution. Within the Local Group, the ability to probe the stellar populations of small and large galaxies opens up the possibility of exploring key questions such as the nature of dark matter, the detailed formation history of different galaxy components, and the role of accretion in galactic formation. Upcoming wide-field surveys promise to extend this ability to all galaxies within 10~Mpc, drastically increasing our capability to decipher galaxy evolution and enabling statistical studies of galaxies' stellar populations. To facilitate the optimum use of these upcoming capabilities we develop a simple formalism to predict the density of resolved stars for an observation of a stellar population at fixed surface brightness and population parameters. We provide an interface to calculate all quantities of interest to this formalism via a public release of the code: \texttt{walter}. This code enables calculation of (i) the expected number density of detected stars, (ii) the exposure time needed to reach certain population features, such as the horizontal branch, and (iii) an estimate of the crowding limit, among other features. These calculations will be very useful for planning surveys with NASA's upcoming Nancy Grace Roman Space Telescope (Roman, formerly WFIRST), which we use for example calculations throughout this work. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 16 pages, 6 figures, submitted to AJ. Comments welcome

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) II. The Spatially Resolved Recent Star Formation History of M33

Authors: Margaret Lazzarini, Benjamin F. Williams, Meredith J. Durbin, Julianne J. Dalcanton, Adam Smercina, Eric F. Bell, Yumi Choi, Andrew Dolphin, Karoline Gilbert, Puragra Guhathakurta, Erik Rosolowsky, Evan Skillman, O. Grace Telford, Daniel Weisz

Abstract: We measure the spatially resolved recent star formation history (SFH) of M33 using optical images taken with the Hubble Space Telescope as part of the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. The area covered by the observations used in this analysis covers a de-projected area of $\sim$38 kpc$^{2}$ and extends to $\sim$3.5 and $\sim$2 kpc from the center… ▽ More We measure the spatially resolved recent star formation history (SFH) of M33 using optical images taken with the Hubble Space Telescope as part of the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. The area covered by the observations used in this analysis covers a de-projected area of $\sim$38 kpc$^{2}$ and extends to $\sim$3.5 and $\sim$2 kpc from the center of M33 along the major and semi-major axes, respectively. We divide the PHATTER optical survey into 2005 regions that measure 24 arcsec, $\sim$100 pc, on a side and fit color magnitude diagrams for each region individually to measure the spatially resolved SFH of M33 within the PHATTER footprint. There are significant fluctuations in the SFH on small spatial scales and also galaxy-wide scales that we measure back to about 630 Myr ago. We observe a more flocculent spiral structure in stellar populations younger than about 80 Myr, while the structure of the older stellar populations is dominated by two spiral arms. We also observe a bar in the center of M33, which dominates at ages older than about 80 Myr. Finally, we find that the mean star formation rate (SFR) over the last 100 Myr within the PHATTER footprint is 0.32$\pm$0.02 M$_{\odot}$ yr$^{-1}$. We measure a current SFR (over the last 10 Myr) of 0.20$\pm$0.03 M$_{\odot}$ yr$^{-1}$. This SFR is slightly higher than previous measurements from broadband estimates, when scaled to account for the fraction of the D25 area covered by the PHATTER survey footprint. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: accepted for publication in ApJ, 37 pages, 18 figures, 6 tables

New Globular Cluster Candidates in the M81 group

Authors: Jiaming Pan, Eric F. Bell, Adam Smercina, Paul Price, Colin T. Slater, Jeremy Bailin, Roelof S. de Jong, Richard D'Souza, In Sung Jang, Antonela Monachesi

Abstract: The study of outer halo globular cluster (GC) populations can give insight into galaxy merging, globular cluster accretion and the origin of GCs. We use archival Subaru Hyper Suprime-Cam (HSC) data in concert with space-based GALEX, IRAC and Gaia EDR3 data to select candidate Globular clusters (GCs) in the outer halo of the M81 group for confirmation and future study. We use a small sample of prev… ▽ More The study of outer halo globular cluster (GC) populations can give insight into galaxy merging, globular cluster accretion and the origin of GCs. We use archival Subaru Hyper Suprime-Cam (HSC) data in concert with space-based GALEX, IRAC and Gaia EDR3 data to select candidate Globular clusters (GCs) in the outer halo of the M81 group for confirmation and future study. We use a small sample of previously-discovered GCs to tune our selection criteria, finding that bright already-known GCs in the M81 group have sizes that are typically slightly larger than the Subaru PSF in our fields. In the optical bands, GCs appear to have colours that are only slightly different from stars. The inclusion of archival IRAC data yields dramatic improvements in colour separation, as the long wavelength baseline aids somewhat in the separation from stars and clearly separates GCs from many compact background galaxies. We show that some previously spectroscopically-identified GCs in the M81 group are instead foreground stars or background galaxies. GCs close to M82 have radial velocities suggesting that they fell into the M81 group along with M82. The overall M81 GC luminosity function is similar to the Milky Way and M31. M81's outer halo GCs are similar to the Milky Way in their metallicities and numbers, and much less numerous than M31's more metal-rich outer halo GC population. These properties reflect differences in the three galaxies' merger histories, highlighting the possibility of using outer halo GCs to trace merger history in larger samples of galaxies. △ Less

Submitted 14 June, 2022; originally announced June 2022.

High time resolution search for prompt radio emission from the long GRB 210419A with the Murchison Widefield Array

Authors: J. Tian, G. E. Anderson, P. J. Hancock, J. C. A. Miller-Jones, M. Sokolowski, N. A. Swainston, A. Rowlinson, A. Williams, D. L. Kaplan, N. Hurley-Walker, J. Morgan, N. D. R. Bhat, D. Ung, S. Tingay, K. W. Bannister, M. E. Bell, B. W. Meyers, M. Walker

Abstract: We present a low-frequency (170\textendash200\,MHz) search for prompt radio emission associated with the long GRB 210419A using the rapid-response mode of the Murchison Widefield Array (MWA), triggering observations with the Voltage Capture System (VCS) for the first time. The MWA began observing GRB 210419A within 89\,s of its detection by \textit{Swift}, enabling us to capture any dispersion del… ▽ More We present a low-frequency (170\textendash200\,MHz) search for prompt radio emission associated with the long GRB 210419A using the rapid-response mode of the Murchison Widefield Array (MWA), triggering observations with the Voltage Capture System (VCS) for the first time. The MWA began observing GRB 210419A within 89\,s of its detection by \textit{Swift}, enabling us to capture any dispersion delayed signal emitted by this GRB for a typical range of redshifts. We conducted a standard single pulse search with a temporal and spectral resolution of $100\,\upmu$s and 10\,kHz over a broad range of dispersion measures from 1 to $5000\,\text{pc}\,\text{cm}^{-3}$, but none were detected. However, fluence upper limits of $77\text{--}224$\,Jy\,ms derived over a pulse width of $0.5\text{--}10$\,ms and a redshift of $0.6<z<4$ are some of the most stringent at low radio frequencies. We compared these fluence limits to the GRB jet-interstellar medium (ISM) interaction model, placing constraints on the fraction of magnetic energy ($ε_{\text{B}}\lesssim[0.05 \text{--} 0.1]$). We also searched for signals during the X-ray flaring activity of GRB 210419A on minute timescales in the image domain and found no emission, resulting in an intensity upper limit of $0.57\,\text{Jy}\,\text{beam}^{-1}$, corresponding to a constraint of $ε_{\text{B}}\lesssim10^{-3}$. Our non-detection could imply that GRB 210419A was at a high redshift, there was not enough magnetic energy for low-frequency emission, or that the radio waves did not escape from the GRB environment. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Comments: Accepted for publication in MNRAS

The DECam Local Volume Exploration Survey Data Release 2

Authors: A. Drlica-Wagner, P. S. Ferguson, M. Adamów, M. Aguena, F. Andrade-Oliveira, D. Bacon, K. Bechtol, E. F. Bell, E. Bertin, P. Bilaji, S. Bocquet, C. R. Bom, D. Brooks, D. L. Burke, J. A. Carballo-Bello, J. L. Carlin, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, F. J. Castander, W. Cerny, C. Chang, Y. Choi, C. Conselice, M. Costanzi , et al. (99 additional authors not shown)

Abstract: We present the second public data release (DR2) from the DECam Local Volume Exploration survey (DELVE). DELVE DR2 combines new DECam observations with archival DECam data from the Dark Energy Survey, the DECam Legacy Survey, and other DECam community programs. DELVE DR2 consists of ~160,000 exposures that cover >21,000 deg^2 of the high Galactic latitude (|b| > 10 deg) sky in four broadband optica… ▽ More We present the second public data release (DR2) from the DECam Local Volume Exploration survey (DELVE). DELVE DR2 combines new DECam observations with archival DECam data from the Dark Energy Survey, the DECam Legacy Survey, and other DECam community programs. DELVE DR2 consists of ~160,000 exposures that cover >21,000 deg^2 of the high Galactic latitude (|b| > 10 deg) sky in four broadband optical/near-infrared filters (g, r, i, z). DELVE DR2 provides point-source and automatic aperture photometry for ~2.5 billion astronomical sources with a median 5σ point-source depth of g=24.3, r=23.9, i=23.5, and z=22.8 mag. A region of ~17,000 deg^2 has been imaged in all four filters, providing four-band photometric measurements for ~618 million astronomical sources. DELVE DR2 covers more than four times the area of the previous DELVE data release and contains roughly five times as many astronomical objects. DELVE DR2 is publicly available via the NOIRLab Astro Data Lab science platform. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 23 pages, 10 figures, 5 tables; to be submitted to AAS Journals; public data release at https://datalab.noirlab.edu/delve/. arXiv admin note: substantial text overlap with arXiv:2103.07476

Report number: FERMILAB-PUB-22-209-LDRD-PPD

LEGA-C: analysis of dynamical masses from ionized gas and stellar kinematics at z~0.8

Authors: Caroline M. S. Straatman, Arjen van der Wel, Josha van Houdt, Rachel Bezanson, Eric F. Bell, Pieter van Dokkum, Francesco D'Eugenio, Marijn Franx, Anna Gallazzi, Anna de Graaff, Michael Maseda, Sharon E. Meidt, Adam Muzzin, David Sobral, Po-Feng Wu

Abstract: We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset provi… ▽ More We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset providing sufficiently deep measurements of a $K_s-$band limited sample. For $\mathrm{M_{dyn,*}}$ we use Jeans Anisotropic Multi-Gaussian Expansion models. For $\mathrm{M_{dyn,eml}}$ we first fit a custom model of a rotating exponential disk with uniform dispersion, whose light is projected through a slit and corrected for beam smearing. We then apply an asymmetric drift correction based on assumptions common in the literature to the fitted kinematic components to obtain the circular velocity, assuming hydrostatic equilibrium. Within the half-light radius, $\mathrm{M_{dyn,eml}}$ is on average lower than $\mathrm{M_{dyn,*}}$, with a mean offset of $-0.15\pm0.016$ dex and galaxy-to-galaxy scatter of $0.19$ dex, reflecting the combined random uncertainty. While data of higher spatial resolution are needed to understand this small offset, it supports the assumption that the galaxy-wide ionized gas kinematics do not predominantly originate from disruptive events such as star formation driven outflows. However, a similar agreement can be obtained without modeling from the integrated emission line dispersions for axis ratios $q<0.8$. This suggests that our current understanding of gas kinematics is not sufficient to efficiently apply asymmetric drift corrections to improve dynamical mass estimates compared to observations lacking the $S/N$ required for spatially extended dynamics. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: 19 pages, 7 figures, accepted for publication by ApJ, in press

R2-D2: Roman and Rubin -- From Data to Discovery

Authors: Suvi Gezari, Misty Bentz, Kishalay De, K. Decker French, Aaron Meisner, Michelle Ntampaka, Robert Jedicke, Ekta Patel, Daniel Perley, Robyn Sanderson, Christian Aganze, Igor Andreoni, Eric F. Bell, Edo Berger, Ian Dell'Antonio, Ryan Foley, Henry Hsieh, Mansi Kasliwal, Joel Kastner, Charles D. Kilpatrick, J. Davy Kirkpatrick, Casey Lam, Karen Meech, Dante Minniti, Ethan O. Nadler , et al. (6 additional authors not shown)

Abstract: The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to… ▽ More The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 29 pages, 12 figures, Table of Implementation Recommendations, Appendix of Community Science Pitches, AURA-commissioned whitepaper submitted to the Director of STScI (Ken Sembach) and the Director of NOIRLab (Pat McCarthy)

Uncertainties associated with the backward integration of dwarf satellites using simple parametric potentials

Authors: Richard D'Souza, Eric F. Bell

Abstract: In order to backward integrate the orbits of Milky Way (MW) dwarf galaxies, much effort has been invested in recent years to constrain their initial phase-space coordinates. Yet equally important are the assumptions on the potential that the dwarf galaxies experience over time, especially given the fact that the MW is currently accreting the Large Magellanic Cloud (LMC). In this work, using a dark… ▽ More In order to backward integrate the orbits of Milky Way (MW) dwarf galaxies, much effort has been invested in recent years to constrain their initial phase-space coordinates. Yet equally important are the assumptions on the potential that the dwarf galaxies experience over time, especially given the fact that the MW is currently accreting the Large Magellanic Cloud (LMC). In this work, using a dark matter-only zoom-in simulation, we test whether the use of common parametric forms of the potential is adequate to successfully backward integrate the orbits of the subhaloes from their present-day positions. We parametrise the recovered orbits and compare them with those from the simulations. We find that simple symmetric parametric forms of the potential fail to capture the complexities and the inhomogeneities of the true potential experienced by the subhaloes. More specifically, modelling a recent massive accretion like that of the LMC as a sum of two spherical parametric potentials leads to substantial errors in the recovered parameters of the orbits. These errors rival those caused due to a) a 30\% uncertainty in the virial mass of the MW and b) not modelling the potential of the recently accreted massive satellite. Our work suggests that i) the uncertainties in the parameters of the recovered orbits of some MW dwarfs may be under-estimated and that ii) researchers should characterise the uncertainties inherent to their choice of integration techniques and assumptions of the potential against cosmological zoom-in simulations of the MW, which include a recently-accreted LMC. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS