Showing 1–50 of 304 results for author: Robotham, A

EPOCHS I. The Discovery and Star Forming Properties of Galaxies in the Epoch of Reionization at $6.5 < z < 18$ with PEARLS and Public JWST data

Authors: Christopher J. Conselice, Nathan Adams, Thomas Harvey, Duncan Austin, Leonardo Ferreira, Katherine Ormerod, Qiao Duan, James Trussler, Qiong Li, Ignas Juodzbalis, Lewi Westcott, Honor Harris, Louise T. C. Seeyave, Asa F. L. Bluck, Rogier A. Windhorst, Rachana Bhatawdekar, Dan Coe, Seth H. Cohen, Cheng Cheng, Simon P. Driver, Brenda Frye, Lukas J. Furtak, Norman A. Grogin, Nimish P. Hathi, Benne W. Holwerda , et al. (10 additional authors not shown)

Abstract: We present in this paper the discovery, properties, and a catalog of 1165 high redshift $6.5 < z < 18$ galaxies found in deep JWST NIRCam imaging from the GTO PEARLS survey combined with data from JWST public fields. We describe our bespoke homogeneous reduction process and our analysis of these areas including the NEP, CEERS, GLASS, NGDEEP, JADES, and ERO SMACS-0723 fields with over 214 arcmin… ▽ More We present in this paper the discovery, properties, and a catalog of 1165 high redshift $6.5 < z < 18$ galaxies found in deep JWST NIRCam imaging from the GTO PEARLS survey combined with data from JWST public fields. We describe our bespoke homogeneous reduction process and our analysis of these areas including the NEP, CEERS, GLASS, NGDEEP, JADES, and ERO SMACS-0723 fields with over 214 arcmin$^{2}$ imaged to depths of $\sim 30$ mag. We describe our rigorous methods for identifying these galaxies, involving the use of Lyman-break strength, detection significance criteria, visual inspection, and integrated photometric redshifts probability distributions predominately at high redshift. Our sample is a robust and highly pure collection of distant galaxies from which we also remove brown dwarf stars, and calculate completeness and contamination from simulations. We include a summary of the basic properties of these $z > 6.5$ galaxies, including their redshift distributions, UV absolute magnitudes, and star formation rates. Our study of these young galaxies reveals a wide range of stellar population properties as seen in their colors and SED fits which we compare to stellar population models, indicating a range of star formation histories, dust, AGN and/or nebular emission. We find a strong trend exists between stellar mass and $(U-V)$ color, as well as the existence of the `main-sequence' of star formation for galaxies as early as $z \sim 12$. This indicates that stellar mass, or an underlying variable correlating with stellar mass, is driving galaxy formation, in agreement with simulation predictions. We also discover ultra-high redshift candidates at $z > 12$ in our sample and describe their properties. Finally, we note a significant observed excess of galaxies compared to models at $z > 12$, revealing a tension between predictions and our observations. △ Less

Submitted 20 July, 2024; originally announced July 2024.

Comments: Submitted to AAS Journals, 37 pages

Milky Way-mass disc galaxies with low-mass stellar haloes have diverse merger histories

Authors: Katy L. Proctor, Aaron D. Ludlow, Claudia del P. Lagos, Aaron S. G. Robotham

Abstract: We use the Eagle simulation to study the relationship between the stellar haloes of Milky-Way-mass galaxies and their recent merger histories. The stellar mass ratio of the most massive merger that occurred since $z=1$ is strongly correlated with the $z=0$ fraction of ex situ stars in the galaxy, $f_{\rm ex\ situ}$, but is weakly correlated with stellar halo mass fraction, $f_{\rm SH}$, particular… ▽ More We use the Eagle simulation to study the relationship between the stellar haloes of Milky-Way-mass galaxies and their recent merger histories. The stellar mass ratio of the most massive merger that occurred since $z=1$ is strongly correlated with the $z=0$ fraction of ex situ stars in the galaxy, $f_{\rm ex\ situ}$, but is weakly correlated with stellar halo mass fraction, $f_{\rm SH}$, particularly for disc galaxies. Contrary to common belief, our results suggest that disc galaxies with low mass stellar haloes do not necessarily have quiescent merger histories; in fact, roughly one quarter have experienced a merger with stellar mass ratio $> 0.1$ since $z=1$. We demonstrate that the population of disc galaxies with low $f_{\rm SH}$ and active merger histories undergo mergers with satellites whose orbits are more circular than average and are approximately co-planar with the disc; instead of contributing significantly to the stellar halo, these mergers lead to discs that contain a substantial fraction of ex situ stellar mass and are thicker and more extended than those of quiescent galaxies. Such mergers also supply fuel that often incites a significant episode of in situ star formation in the disc. Our results suggest that seemingly quiescent disc galaxies with low-mass stellar haloes actually have diverse merger histories, which limits the extent to which stellar haloes alone can act as observable tracers of the hierarchical assembly of galaxies. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: 11 pages, 8 figures. Submitted to MNRAS

Galaxy Mergers in the Epoch of Reionization I: A JWST Study of Pair Fractions, Merger Rates, and Stellar Mass Accretion Rates at $z = 4.5-11.5$

Authors: Qiao Duan, Christopher J. Conselice, Qiong Li, Duncan Austin, Thomas Harvey, Nathan J. Adams, Kenneth J. Duncan, James Trussler, Leonardo Ferreira, Lewi Westcott, Honor Harris, Rogier A. Windhorst, Benne W. Holwerda, Thomas J. Broadhurst, Dan Coe, Seth H. Cohen, Simon P. Driver, Brenda Frye, Norman A. Grogin, Nimish P. Hathi, Rolf A. Jansen, Anton M. Koekemoer, Madeline A. Marshall, Mario Nonino, Rafael Ortiz III , et al. (7 additional authors not shown)

Abstract: We present a full analysis of galaxy major merger pair fractions, merger rates, and mass accretion rates, thus uncovering the role of mergers in galaxy formation at the earliest previously unexplored epoch of $4.5<z<11.5$. We target galaxies with masses $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) = 8.0 - 10.0$, utilizing data from eight JWST Cycle-1 fields (CEERS, JADES GOODS-S, NEP-TDF, NGDEEP, GLA… ▽ More We present a full analysis of galaxy major merger pair fractions, merger rates, and mass accretion rates, thus uncovering the role of mergers in galaxy formation at the earliest previously unexplored epoch of $4.5<z<11.5$. We target galaxies with masses $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) = 8.0 - 10.0$, utilizing data from eight JWST Cycle-1 fields (CEERS, JADES GOODS-S, NEP-TDF, NGDEEP, GLASS, El-Gordo, SMACS-0723, MACS-0416), covering an unmasked area of 189.36 $\mathrm{arcmin}^2$. We develop a new probabilistic pair-counting methodology that integrates full photometric redshift posteriors and corrects for detection incompleteness to quantify close pairs with physical projected separations between 20 and 50 kpc. Our analysis reveals an increase in pair fractions up to $z = 8$, reaching $0.211 \pm 0.065$, followed by a statistically flat evolution to $z = 11.5$. We find that the galaxy merger rate increases from the local Universe up to $z = 6$ and then stabilizes at a value of $\sim 6$ Gyr$^{-1}$ up to $z = 11.5$. We fit both a power-law and a power-law + exponential model to our pair fraction and merger rate redshift evolution, finding that the latter model describes the trends more accurately, particularly at $z = 8.0 - 11.5$. In addition, we measure that the average galaxy increases its stellar mass due to mergers by a factor of $2.77 \pm 0.99$ from redshift $z = 10.5$ to $z = 5.0$. Lastly, we investigate the impact of mergers on galaxy stellar mass growth, revealing that mergers contribute $71 \pm 25\%$ as much to galaxy stellar mass increases as star formation from gas. This indicates that mergers drive about half of galaxy assembly at high redshift. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 30 Pages, 15 Figures, Submitted to MNRAS

Predicting the Non-Thermal Pressure in Galaxy Clusters

Authors: Andrew Sullivan, Stanislav Shabala, Chris Power, Connor Bottrell, Aaron Robotham

Abstract: We investigate the relationship between a galaxy cluster's hydrostatic equilibrium state, the entropy profile, $K$, of the intracluster gas, and the system's non-thermal pressure (NTP), within an analytic model of cluster structures. When NTP is neglected from the cluster's hydrostatic state, we find that the gas' logarithmic entropy slope, $k\equiv \mathrm{d}\ln K/\mathrm{d}\ln r$, converges at l… ▽ More We investigate the relationship between a galaxy cluster's hydrostatic equilibrium state, the entropy profile, $K$, of the intracluster gas, and the system's non-thermal pressure (NTP), within an analytic model of cluster structures. When NTP is neglected from the cluster's hydrostatic state, we find that the gas' logarithmic entropy slope, $k\equiv \mathrm{d}\ln K/\mathrm{d}\ln r$, converges at large halocentric radius, $r$, to a value that is systematically higher than the value $k\simeq1.1$ that is found in observations and simulations. By applying a constraint on these `pristine equilibrium' slopes, $k_\mathrm{eq}$, we are able to predict the required NTP that must be introduced into the hydrostatic state of the cluster. We solve for the fraction, $\mathcal{F}\equiv p_\mathrm{nt}/p$, of NTP, $p_\mathrm{nt}$, to total pressure, $p$, of the cluster, and we find $\mathcal{F}(r)$ to be an increasing function of halocentric radius, $r$, that can be parameterised by its value in the cluster's core, $\mathcal{F}_0$, with this prediction able to be fit to the functional form proposed in numerical simulations. The minimum NTP fraction, as the solution with zero NTP in the core, $\mathcal{F}_0=0$, we find to be in excellent agreement with the mean NTP predicted in non-radiative simulations, beyond halocentric radii of $r\gtrsim0.7r_{500}$, and in tension with observational constraints derived at similar radii. For this minimum NTP profile, we predict $\mathcal{F}\simeq0.20$ at $r_{500}$, and $\mathcal{F}\simeq0.34$ at $2r_{500}$; this amount of NTP leads to a hydrostatic bias of $b\simeq0.12$ in the cluster mass $M_{500}$ when measured within $r_{500}$. Our results suggest that the NTP of galaxy clusters contributes a significant amount to their hydrostatic state near the virial radius, and must be accounted for when estimating the cluster's halo mass using hydrostatic equilibrium approaches. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 10 pages, 5 figures. Accepted for publication in PASA

Wide Area VISTA Extra-galactic Survey (WAVES): Unsupervised star-galaxy separation on the WAVES-Wide photometric input catalogue using UMAP and ${\rm{\scriptsize HDBSCAN}}$

Authors: Todd L. Cook, Behnood Bandi, Sam Philipsborn, Jon Loveday, Sabine Bellstedt, Simon P. Driver, Aaron S. G. Robotham, Maciej Bilicki, Gursharanjit Kaur, Elmo Tempel, Ivan Baldry, Daniel Gruen, Marcella Longhetti, Angela Iovino, Benne W. Holwerda, Ricardo Demarco

Abstract: Star-galaxy separation is a crucial step in creating target catalogues for extragalactic spectroscopic surveys. A classifier biased towards inclusivity risks including spurious stars, wasting fibre hours, while a more conservative classifier might overlook galaxies, compromising completeness and hence survey objectives. To avoid bias introduced by a training set in supervised methods, we employ an… ▽ More Star-galaxy separation is a crucial step in creating target catalogues for extragalactic spectroscopic surveys. A classifier biased towards inclusivity risks including spurious stars, wasting fibre hours, while a more conservative classifier might overlook galaxies, compromising completeness and hence survey objectives. To avoid bias introduced by a training set in supervised methods, we employ an unsupervised machine learning approach. Using photometry from the Wide Area VISTA Extragalactic Survey (WAVES)-Wide catalogue comprising 9-band $u-K_s$ data, we create a feature space with colours, fluxes, and apparent size information extracted by ${\rm P{\scriptsize RO} F{\scriptsize OUND}}$. We apply the non-linear dimensionality reduction method UMAP (Uniform Manifold Approximation and Projection) combined with the classifier ${\rm{\scriptsize HDBSCAN}}$ to classify stars and galaxies. Our method is verified against a baseline colour and morphological method using a truth catalogue from Gaia, SDSS, GAMA, and DESI. We correctly identify 99.72% of galaxies within the AB magnitude limit of $Z = 21.2$, with an F1 score of 0.9970 across the entire ground truth sample, compared to 0.9871 from the baseline method. Our method's higher purity (0.9966) compared to the baseline (0.9780) increases efficiency, identifying 11% fewer galaxy or ambiguous sources, saving approximately 70,000 fibre hours on the 4MOST instrument. We achieve reliable classification statistics for challenging sources including quasars, compact galaxies, and low surface brightness galaxies, retrieving 95.1%, 84.6%, and 99.5% of them respectively. Angular clustering analysis validates our classifications, showing consistency with expected galaxy clustering, regardless of the baseline classification. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS; 17 pages, 11 figures

JWST view of three infant galaxies at z=8.3 and implications for reionization

Authors: Zhiyuan Ma, Bangzheng Sun, Cheng Cheng, Haojing Yan, Fengwu Sun, Nicholas Foo, Eiichi Egami, Jose M. Diego, Seth H. Cohen, Rolf A. Jansen, Jake Summers, Rogier A. Windhorst, Jordan C. J. D'Silva, Anton M. Koekemoer, Dan Coe, Christopher J. Conselice, Simon P. Driver, Brenda Frye, Norman A. Grogin, Madeline A. Marshall, Mario Nonino, Rafael Ortiz III, Nor Pirzkal, Aaron Robotham, Russell E. Ryan, Jr. , et al. (12 additional authors not shown)

Abstract: New JWST/NIRCam wide-field slitless spectroscopy provides redshifts for two z > 8 galaxies located behind the lensing cluster MACS J0416.1-2403. Both galaxies are strong [O iii]λ5007 emitters. For one galaxy, "Y1", the existing redshift z = 8.31, based on ALMA measurements of [O iii] 88 μm and [C ii] 157.7 μm lines, is confirmed. JWST/NIRCam images resolve this galaxy into three components of simi… ▽ More New JWST/NIRCam wide-field slitless spectroscopy provides redshifts for two z > 8 galaxies located behind the lensing cluster MACS J0416.1-2403. Both galaxies are strong [O iii]λ5007 emitters. For one galaxy, "Y1", the existing redshift z = 8.31, based on ALMA measurements of [O iii] 88 μm and [C ii] 157.7 μm lines, is confirmed. JWST/NIRCam images resolve this galaxy into three components of similar colors, and the whole system extends over ~3.4 kpc. The other galaxy, "JD", is at z = 8.34 instead of the previously claimed z = 9.28. It has a companion, "JD-N", at the same redshift with projected separation ~2.3 kpc. All objects are only moderately magnified and have intrinsic MUV ranging from -19.66 to -20.85 mag. Their eight-band NIRCam spectral energy distributions show that the galaxies are all very young with ages $\lesssim$11 Myr and stellar masses about 108 $M_{\odot}$. These infant galaxies are actively forming stars at rates of a few tens to a couple of hundred $M_{\odot} yr^{-1}$, but only one of them (JD) has a blue rest-frame UV slope. This slope indicates a high Lyman-continuum photon escape fraction that could contribute significantly to the cosmic hydrogen-reionizing background. The other two systems have much flatter slopes largely because their dust extinction is twice as high as JD's albeit only AV ~ 0.90 mag. The much lower indicated escape fractions show that even very young, actively star-forming galaxies can have negligible contribution to reionization when they quickly form dust throughout their bodies. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 18 pages, 6 figures, submitted to ApJL

JWST's PEARLS: resolved study of the stellar and dust components in starburst galaxies at cosmic noon

Authors: M. Polletta, B. L. Frye, N. Garuda, S. P. Willner, S. Berta, R. Kneissl, H. Dole, R. A. Jansen, M. D. Lehnert, S. H. Cohen, J. Summers, R. A. Windhorst, J. C. J. D'Silva, A. M. Koekemoer, D. Coe, C. J. Conselice, S. P. Driver, N. A. Grogin, M. A. Marshall, M. Nonino, R. Ortiz III, N. Pirzkal, A. Robotham, R. E. Ryan, Jr., C. N. A. Willmer , et al. (13 additional authors not shown)

Abstract: Dusty star-forming galaxies (DSFGs) contribute significantly to the stellar buildup at cosmic noon. Major mergers and gas accretion are often invoked to explain DSFGs' prodigious star-formation rates (SFRs) and large stellar masses. We conducted a spatially-resolved morphological analysis of the rest-frame UV/NIR emission in three DSFGs at z~2.5. Initially discovered as CO emitters by NOEMA observ… ▽ More Dusty star-forming galaxies (DSFGs) contribute significantly to the stellar buildup at cosmic noon. Major mergers and gas accretion are often invoked to explain DSFGs' prodigious star-formation rates (SFRs) and large stellar masses. We conducted a spatially-resolved morphological analysis of the rest-frame UV/NIR emission in three DSFGs at z~2.5. Initially discovered as CO emitters by NOEMA observations of a bright Herschel source, we observed them with the JWST/NIRCam as part of the PEARLS program. The NIRCam data reveal the galaxies' stellar population and dust distribution on scales of 250 pc. Spatial variations in stellar mass, SFR, and dust extinction are determined in resolved maps obtained through pixel-based SED fitting. The CO emitters are massive, dusty starburst galaxies with SFRs ranging from 340 to 2500 Msun/yr, positioning them among the most active SFGs at 2<z<3. Notably, they belong to the ~1.5% of the entire JWST population with extremely red colors. Their morphologies are disk-like, with effective radii of 2.0-4.4 kpc, and exhibit sub-structures such as clumps and spiral arms. The galaxies have dust extinctions up to Av=5-7 mag with asymmetric distributions extending over several kpc and including off-center regions resembling bent spiral arms and clumps. The NIR dust-attenuation curve in these sources deviates from standard laws, implying different dust grain properties than commonly assumed in starburst galaxies. The proximity of galaxies with consistent redshifts, strong color gradients, overall disturbed appearance, asymmetric dust obscuration, and wide-spread star formation favor interactions (minor mergers and flybys) as the mechanism driving the CO galaxies' exceptional SFRs. Their large masses and rich environment hint at membership in two proto-structures, as initially inferred from their association with a Planck-selected high-z source. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 24 pages, 21 figures + appendix. Submitted to A&A. Comments welcome!

DEVILS/MIGHTEE/GAMA/DINGO: The Impact of SFR Timescales on the SFR-Radio Luminosity Correlation

Authors: Robin H. W. Cook, Luke J. M. Davies, Jonghwan Rhee, Catherine L. Hale, Sabine Bellstedt, Jessica E. Thorne, Ivan Delvecchio, Jordan D. Collier, Richard Dodson, Simon P. Driver, Benne W. Holwerda, Matt J. Jarvis, Kenda Knowles, Claudia Lagos, Natasha Maddox, Martin Meyer, Aaron S. G. Robotham, Sambit Roychowdhury, Kristof Rozgonyi, Nicholas Seymour, Malgorzata Siudek, Matthew Whiting, Imogen Whittam

Abstract: The tight relationship between infrared luminosity (L$_\mathrm{TIR}$) and 1.4 GHz radio continuum luminosity (L$_\mathrm{1.4GHz}$) has proven useful for understanding star formation free from dust obscuration. Infrared emission in star-forming galaxies typically arises from recently formed, dust-enshrouded stars, whereas radio synchrotron emission is expected from subsequent supernovae. By leverag… ▽ More The tight relationship between infrared luminosity (L$_\mathrm{TIR}$) and 1.4 GHz radio continuum luminosity (L$_\mathrm{1.4GHz}$) has proven useful for understanding star formation free from dust obscuration. Infrared emission in star-forming galaxies typically arises from recently formed, dust-enshrouded stars, whereas radio synchrotron emission is expected from subsequent supernovae. By leveraging the wealth of ancillary far-ultraviolet - far-infrared photometry from the Deep Extragalactic VIsible Legacy Survey (DEVILS) and Galaxy and Mass Assembly (GAMA) surveys, combined with 1.4 GHz observations from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey and Deep Investigation of Neutral Gas Origins (DINGO) projects, we investigate the impact of timescale differences between far-ultraviolet - far-infrared and radio-derived star formation rate (SFR) tracers. We examine how the SED-derived star formation histories (SFH) of galaxies can be used to explain discrepancies in these SFR tracers, which are sensitive to different timescales. Galaxies exhibiting an increasing SFH have systematically higher L$_\mathrm{TIR}$ and SED-derived SFRs than predicted from their 1.4 GHz radio luminosity. This indicates that insufficient time has passed for subsequent supernovae-driven radio emission to accumulate. We show that backtracking the SFR(t) of galaxies along their SED-derived SFHs to a time several hundred megayears prior to their observed epoch will both linearise the SFR-L$_\mathrm{1.4GHz}$ relation and reduce the overall scatter. The minimum scatter in the SFR(t)-L$_\mathrm{1.4GHz}$ is reached at 200 - 300 Myr prior, consistent with theoretical predictions for the timescales required to disperse the cosmic ray electrons responsible for the synchrotron emission. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 23 pages, 13 figures, 2 tables. Accepted for publication in MNRAS

EPOCHS III: Unbiased UV continuum slopes at 6.5<z<13 from combined PEARLS GTO and public JWST NIRCam imaging

Authors: Duncan Austin, Christopher J. Conselice, Nathan J. Adams, Thomas Harvey, Qiao Duan, James Trussler, Qiong Li, Ignas Juodzbalis, Katherine Ormerod, Leonardo Ferreira, Lewi Westcott, Honor Harris, Stephen M. Wilkins, Rachana Bhatawdekar, Joseph Caruana, Dan Coe, Seth H. Cohen, Simon P. Driver, Jordan C. J. D'Silva, Brenda Frye, Lukas J. Furtak, Norman A. Grogin, Nimish P. Hathi, Benne W. Holwerda, Rolf A. Jansen , et al. (12 additional authors not shown)

Abstract: We present an analysis of rest-frame UV continuum slopes, $β$, using a sample of 1011 galaxies at $6.5<z<13$ from the EPOCHS photometric sample collated from the GTO PEARLS and public ERS/GTO/GO (JADES, CEERS, NGDEEP, GLASS) JWST NIRCam imaging across $178.9~\mathrm{arcmin}^2$ of unmasked blank sky. We correct our UV slopes for the photometric error coupling bias using $200,000$ power law SEDs for… ▽ More We present an analysis of rest-frame UV continuum slopes, $β$, using a sample of 1011 galaxies at $6.5<z<13$ from the EPOCHS photometric sample collated from the GTO PEARLS and public ERS/GTO/GO (JADES, CEERS, NGDEEP, GLASS) JWST NIRCam imaging across $178.9~\mathrm{arcmin}^2$ of unmasked blank sky. We correct our UV slopes for the photometric error coupling bias using $200,000$ power law SEDs for each $β=\{-1,-1.5,-2,-2.5,-3\}$ in each field, finding biases as large as $Δβ\simeq-0.55$ for the lowest SNR galaxies in our sample. Additionally, we simulate the impact of rest-UV line emission (including Ly$α$) and damped Ly$α$ systems on our measured $β$, finding biases as large as $0.5-0.6$ for the most extreme systems. We find a decreasing trend with redshift of $β=-1.51\pm0.08-(0.097\pm0.010)\times z$, with potential evidence for Pop.~III stars or top-heavy initial mass functions (IMFs) in a subsample of 68 $β+σ_β<-2.8$ galaxies. At $z\simeq11.5$, we measure an extremely blue $β(M_{\mathrm{UV}}=-19)=-2.73\pm0.06$, deviating from simulations, indicative of low-metallicity galaxies with non-zero Lyman continuum escape fractions $f_{\mathrm{esc, LyC}}\gtrsim0$ and minimal dust content. The observed steepening of $\mathrm{d}β/\mathrm{d}\log_{10}(M_{\star}/\mathrm{M}_{\odot})$ from $0.22\pm0.02$ at $z=7$ to $0.81\pm0.13$ at $z=11.5$ implies that dust produced in core-collapse supernovae (SNe) at early times may be ejected via outflows from low mass galaxies. We also observe a flatter $\mathrm{d}β/\mathrm{d}M_{\mathrm{UV}}=0.03\pm0.02$ at $z=7$ and a shallower $\mathrm{d}β/\mathrm{d}\log_{10}(M_{\star} / \mathrm{M}_{\odot})$ at $z<11$ than seen by HST, unveiling a new population of low mass, faint, galaxies reddened by dust produced in the stellar winds of asymptotic giant branch (AGB) stars or carbon-rich Wolf-Rayet binaries. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 37 pages, 18 figures, 8 tables. Data available via the Github page at https://github.com/duncanaustin98/EPOCHS-III-data. Submitted to ApJ on 16/04/2024

PEARLS: Discovery of Point-Source Features Within Galaxies in the North Ecliptic Pole Time Domain Field

Authors: Rafael Ortiz III, Rogier A. Windhorst, Seth H. Cohen, S. P. Willner, Rolf A. Jansen, Timothy Carleton, Patrick S. Kamieneski, Michael J. Rutkowski, Brent Smith, Jake Summers, Tyler J. McCabe, Rosalia O'Brien, Jose M. Diego, Min S. Yun, Jordan C. J. D'Silva, Juno Li, Hansung B. Gim, Nimish P. Hathi, Benne W. Holwerda, Adi Zitrin, Cheng Cheng, Noah J. McLeod, Christopher J. Conselice, Simon P. Driver, Haojing Yan , et al. (9 additional authors not shown)

Abstract: $ $The first public 0.9-4.4 $μ$m NIRCam images of the North Ecliptic Pole (NEP) Time Domain Field (TDF) uncovered many galaxies that display point-source features in their cores as seen in the longer wavelength filters. We visually identified a sample of 66 galaxies ($\sim$1 galaxy per arcmin$^2… ▽ More $ $The first public 0.9-4.4 $μ$m NIRCam images of the North Ecliptic Pole (NEP) Time Domain Field (TDF) uncovered many galaxies that display point-source features in their cores as seen in the longer wavelength filters. We visually identified a sample of 66 galaxies ($\sim$1 galaxy per arcmin$^2$) with point-like cores and fit their spectral energy distributions (SED)s using EAZY and CIGALE to characterize the sample's active galactic nucleus (AGN) and host galaxy parameters. Single-template fitting best fits $70\%$ of the sample with a Seyfert-blended SED. With CIGALE we compute the median fractional AGN contribution to the 0.1-30.0 $μ$m flux to be $0.30\pm0.06$, and that $56\%$ of the 66 galaxies have star-formation rates in the starburst range whereas the remainder are near the star-formation main sequence. There are Very Large Array (VLA) 3 GHz detections for 24/66 galaxies, implying some combination of AGN emission and vigorous star formation. We present a novel sample selection procedure in tandem to the morphological sample selection based on objects' light profiles at 4.4 $μ$m. This procedure identifies a parameter space which probes point-source features and automatically recovers our visual sample with minimal contamination from both stars and brighter galaxies. This procedure may be used in other extant and future NIRCam images to streamline the search for galaxies with point-like cores. The morphological approach to recognizing AGN is being resurrected by the James Webb Space Telescope (JWST) by virtue of its superb angular resolution at infrared wavelengths. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 11 pages, 6 figures, 1 table

Birds of a Feather: Resolving Stellar Mass Assembly With JWST/NIRCam in a Pair of Kindred $z \sim 2$ Dusty Star-forming Galaxies Lensed by the PLCK G165.7+67.0 Cluster

Authors: Patrick S. Kamieneski, Brenda L. Frye, Rogier A. Windhorst, Kevin C. Harrington, Min S. Yun, Allison Noble, Massimo Pascale, Nicholas Foo, Seth H. Cohen, Rolf A. Jansen, Timothy Carleton, Anton M. Koekemoer, Christopher N. A. Willmer, Jake S. Summers, Nikhil Garuda, Reagen Leimbach, Benne W. Holwerda, Justin D. R. Pierel, Eric F. Jimenez-Andrade, S. P. Willner, Belen Alcalde Pampliega, Amit Vishwas, William C. Keel, Q. Daniel Wang, Cheng Cheng , et al. (16 additional authors not shown)

Abstract: We present a new parametric lens model for the G165.7+67.0 galaxy cluster, which was discovered with $Planck$ through its bright submillimeter flux, originating from a pair of extraordinary dusty star-forming galaxies (DSFGs) at $z\approx 2.2$. Using JWST and interferometric mm/radio observations, we characterize the intrinsic physical properties of the DSFGs, which are separated by only… ▽ More We present a new parametric lens model for the G165.7+67.0 galaxy cluster, which was discovered with $Planck$ through its bright submillimeter flux, originating from a pair of extraordinary dusty star-forming galaxies (DSFGs) at $z\approx 2.2$. Using JWST and interferometric mm/radio observations, we characterize the intrinsic physical properties of the DSFGs, which are separated by only $\sim 1^{\prime\prime}$ (8 kpc) and a velocity difference $ΔV \lesssim 600~{\rm km}~{\rm s}^{-1}$ in the source plane, and thus likely undergoing a major merger. Boasting intrinsic star formation rates ${\rm SFR}_{\rm IR} = 320 \pm 70$ and $400 \pm 80~ M_\odot~{\rm yr}^{-1}$, stellar masses ${\rm log}[M_\star/M_\odot] = 10.2 \pm 0.1$ and $10.3 \pm 0.1$, and dust attenuations $A_V = 1.5 \pm 0.3$ and $1.2 \pm 0.3$, they are remarkably similar objects. We perform spatially-resolved pixel-by-pixel SED fitting using rest-frame near-UV to near-IR imaging from JWST/NIRCam for both galaxies, resolving some stellar structures down to 100 pc scales. Based on their resolved specific SFRs and $UVJ$ colors, both DSFGs are experiencing significant galaxy-scale star formation events. If they are indeed interacting gravitationally, this strong starburst could be the hallmark of gas that has been disrupted by an initial close passage. In contrast, the host galaxy of the recently discovered triply-imaged SN H0pe has a much lower SFR than the DSFGs, and we present evidence for the onset of inside-out quenching and large column densities of dust even in regions of low specific SFR. Based on the intrinsic SFRs of the DSFGs inferred from UV through FIR SED modeling, this pair of objects alone is predicted to yield an observable $1.1 \pm 0.2~{\rm CCSNe~yr}^{-1}$, making this cluster field ripe for continued monitoring. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 47 pages, 21 figures, 5 tables. Submitted to ApJ, comments welcome!

Lensed Type Ia Supernova "Encore" at z=2: The First Instance of Two Multiply-Imaged Supernovae in the Same Host Galaxy

Authors: J. D. R. Pierel, A. B. Newman, S. Dhawan, M. Gu, B. A. Joshi, T. Li, S. Schuldt, L. G. Strolger, S. H. Suyu, G. B. Caminha, S. H. Cohen, J. M. Diego, J. C. J. Dsilva, S. Ertl, B. L. Frye, G. Granata, C. Grillo, A. M. Koekemoer, J. Li, A. Robotham, J. Summers, T. Treu, R. A. Windhorst, A. Zitrin, S. Agarwal , et al. (38 additional authors not shown)

Abstract: A bright ($m_{\rm F150W,AB}$=24 mag), $z=1.95$ supernova (SN) candidate was discovered in JWST/NIRCam imaging acquired on 2023 November 17. The SN is quintuply-imaged as a result of strong gravitational lensing by a foreground galaxy cluster, detected in three locations, and remarkably is the second lensed SN found in the same host galaxy. The previous lensed SN was called "Requiem", and therefore… ▽ More A bright ($m_{\rm F150W,AB}$=24 mag), $z=1.95$ supernova (SN) candidate was discovered in JWST/NIRCam imaging acquired on 2023 November 17. The SN is quintuply-imaged as a result of strong gravitational lensing by a foreground galaxy cluster, detected in three locations, and remarkably is the second lensed SN found in the same host galaxy. The previous lensed SN was called "Requiem", and therefore the new SN is named "Encore". This makes the MACS J0138.0$-$2155 cluster the first known system to produce more than one multiply-imaged SN. Moreover, both SN Requiem and SN Encore are Type Ia SNe (SNe Ia), making this the most distant case of a galaxy hosting two SNe Ia. Using parametric host fitting, we determine the probability of detecting two SNe Ia in this host galaxy over a $\sim10$ year window to be $\approx3\%$. These observations have the potential to yield a Hubble Constant ($H_0$) measurement with $\sim10\%$ precision, only the third lensed SN capable of such a result, using the three visible images of the SN. Both SN Requiem and SN Encore have a fourth image that is expected to appear within a few years of $\sim2030$, providing an unprecedented baseline for time-delay cosmography. △ Less

Submitted 22 July, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: Accepted, ApJL

JWST Photometric Time-Delay and Magnification Measurements for the Triply-Imaged Type Ia "Supernova H0pe" at z = 1.78

Authors: J. D. R. Pierel, B. L. Frye, M. Pascale, G. B. Caminha, W. Chen, S. Dhawan, D. Gilman, M. Grayling, S. Huber, P. Kelly, S. Thorp, N. Arendse, S. Birrer, M. Bronikowski, R. Canameras, D. Coe, S. H. Cohen, C. J. Conselice, S. P. Driver, J. C. J. Dsilva, M. Engesser, N. Foo, C. Gall, N. Garuda, C. Grillo , et al. (38 additional authors not shown)

Abstract: Supernova (SN) H0pe is a gravitationally lensed, triply-imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply-imaged SNe provide a rare opportunity to constrain the Hubble constant ($H_0$), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is locate… ▽ More Supernova (SN) H0pe is a gravitationally lensed, triply-imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply-imaged SNe provide a rare opportunity to constrain the Hubble constant ($H_0$), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is located at $z=1.783$, and is the first SN Ia with sufficient light curve sampling and long enough time delays for an $H_0$ inference. Here we present photometric time-delay measurements and SN properties of SN H0pe. Using JWST/NIRCam photometry we measure time delays of $Δt_{ab}=-116.6^{+10.8}_{-9.3}$ and $Δt_{cb}=-48.6^{+3.6}_{-4.0}$ observer-frame days relative to the last image to arrive (image 2b; all uncertainties are $1σ$), which corresponds to a $\sim5.6\%$ uncertainty contribution for $H_0$ assuming $70 \rm{km s^{-1} Mpc^{-1}}$. We also constrain the absolute magnification of each image to $μ_{a}=4.3^{+1.6}_{-1.8}$, $μ_{b}=7.6^{+3.6}_{-2.6}$, $μ_{c}=6.4^{+1.6}_{-1.5}$ by comparing the observed peak near-IR magnitude of SN H0pe to the non-lensed population of SNe Ia. △ Less

Submitted 22 July, 2024; v1 submitted 27 March, 2024; originally announced March 2024.

Comments: Accepted in ApJ

Predicting the Scaling Relations between the Dark Matter Halo Mass and Observables from Generalised Profiles II: Intracluster Gas Emission

Authors: Andrew Sullivan, Chris Power, Connor Bottrell, Aaron Robotham, Stas Shabala

Abstract: We investigate the connection between a cluster's structural configuration and observable measures of its gas emission that can be obtained in X-ray and Sunyaev-Zeldovich (SZ) surveys. We present an analytic model for the intracluster gas density profile: parameterised by the dark matter halo's inner logarithmic density slope, $α$, the concentration, $c$, the gas profile's inner logarithmic densit… ▽ More We investigate the connection between a cluster's structural configuration and observable measures of its gas emission that can be obtained in X-ray and Sunyaev-Zeldovich (SZ) surveys. We present an analytic model for the intracluster gas density profile: parameterised by the dark matter halo's inner logarithmic density slope, $α$, the concentration, $c$, the gas profile's inner logarithmic density slope, $\varepsilon$, the dilution, $d$, and the gas fraction, $η$, normalised to cosmological content. We predict four probes of the gas emission: the emission-weighted, $T_\mathrm{X}$, and mean gas mass-weighted, $T_\mathrm{m_g}$, temperatures, and the spherically, $Y_\mathrm{sph}$, and cylindrically, $Y_\mathrm{cyl}$, integrated Compton parameters. Over a parameter space of clusters, we constrain the X-ray temperature scaling relations, $M_{200} - T_\mathrm{X}$ and $M_{500} - T_\mathrm{X}$, within $57.3\%$ and $41.6\%$, and $M_{200} - T_\mathrm{m_g}$ and $M_{500} - T_\mathrm{m_g}$, within $25.7\%$ and $7.0\%$, all respectively. When excising the cluster's core, the $M_{200} - T_\mathrm{X}$ and $M_{500} - T_\mathrm{X}$ relations are further constrained, to within $31.3\%$ and $17.1\%$, respectively. Similarly, we constrain the SZ scaling relations, $M_{200} - Y_\mathrm{sph}$ and $M_{500} - Y_\mathrm{sph}$, within $31.1\%$ and $17.7\%$, and $M_{200} - Y_\mathrm{cyl}$ and $M_{500} - Y_\mathrm{cyl}$, within $25.2\%$ and $22.0\%$, all respectively. The temperature observable $T_\mathrm{m_g}$ places the strongest constraint on the halo mass, whilst $T_\mathrm{X}$ is more sensitive to the parameter space. The SZ constraints are sensitive to the gas fraction, whilst insensitive to the form of the gas profile itself. In all cases, the halo mass is recovered with an uncertainty that suggests the cluster's structural profiles only contribute a minor uncertainty in its scaling relations. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: This paper has been accepted for publication with PASA

EMU/GAMA: A Technique for Detecting Active Galactic Nuclei in Low Mass Systems

Authors: Jahang Prathap, Andrew M. Hopkins, Aaron S. G. Robotham, Sabine Bellstedt, José Afonso, Ummee T. Ahmed, Maciej Bilicki, Malcolm N. Bremer, Sarah Brough, Michael J. I. Brown, Yjan Gordon, Benne W. Holwerda, Denis Leahy, Ángel R. López-Sánchez, Joshua R. Marvil, Tamal Mukherjee, Isabella Prandoni, Stanislav S. Shabala, Tessa Vernstrom, Tayyaba Zafar

Abstract: We propose a new method for identifying active galactic nuclei (AGN) in low mass ($\rm M_*\leq10^{10}M_\odot$) galaxies. This method relies on spectral energy distribution (SED) fitting to identify galaxies whose radio flux density has an excess over that expected from star formation alone. Combining data in the Galaxy and Mass Assembly (GAMA) G23 region from GAMA, Evolutionary Map of the Universe… ▽ More We propose a new method for identifying active galactic nuclei (AGN) in low mass ($\rm M_*\leq10^{10}M_\odot$) galaxies. This method relies on spectral energy distribution (SED) fitting to identify galaxies whose radio flux density has an excess over that expected from star formation alone. Combining data in the Galaxy and Mass Assembly (GAMA) G23 region from GAMA, Evolutionary Map of the Universe (EMU) early science observations, and Wide-field Infrared Survey Explorer (WISE), we compare this technique with a selection of different AGN diagnostics to explore the similarities and differences in AGN classification. We find that diagnostics based on optical and near-infrared criteria (the standard BPT diagram, the WISE colour criterion, and the mass-excitation, or MEx diagram) tend to favour detection of AGN in high mass, high luminosity systems, while the ``ProSpect'' SED fitting tool can identify AGN efficiently in low mass systems. We investigate an explanation for this result in the context of proportionally lower mass black holes in lower mass galaxies compared to higher mass galaxies and differing proportions of emission from AGN and star formation dominating the light at optical and infrared wavelengths as a function of galaxy stellar mass. We conclude that SED-derived AGN classification is an efficient approach to identify low mass hosts with low radio luminosity AGN. △ Less

Submitted 18 February, 2024; originally announced February 2024.

Comments: 16 pages, 12 figures, 5 tables. Accepted for publication in PASA

TREASUREHUNT: Transients and Variability Discovered with HST in the JWST North Ecliptic Pole Time Domain Field

Authors: Rosalia O'Brien, Rolf A. Jansen, Norman A. Grogin, Seth H. Cohen, Brent M. Smith, Ross M. Silver, W. P. Maksym III, Rogier A. Windhorst, Timothy Carleton, Anton M. Koekemoer, Nimish P. Hathi, Christopher N. A. Willmer, Brenda L. Frye, M. Alpaslan, M. L. N. Ashby, T. A. Ashcraft, S. Bonoli, W. Brisken, N. Cappelluti, F. Civano, C. J. Conselice, V. S. Dhillon, S. P. Driver, K. J. Duncan, R. Dupke , et al. (34 additional authors not shown)

Abstract: The JWST North Ecliptic Pole (NEP) Time Domain Field (TDF) is a $>$14 arcmin diameter field optimized for multi-wavelength time-domain science with JWST. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over 3 cycles (the "TREASUREHUNT" program), deep images were obtained with… ▽ More The JWST North Ecliptic Pole (NEP) Time Domain Field (TDF) is a $>$14 arcmin diameter field optimized for multi-wavelength time-domain science with JWST. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over 3 cycles (the "TREASUREHUNT" program), deep images were obtained with ACS/WFC in F435W and F606W that cover almost the entire JWST NEP TDF. Many of the individual pointings of these programs partially overlap, allowing an initial assessment of the potential of this field for time-domain science with HST and JWST. The cumulative area of overlapping pointings is ~88 arcmin$^2$, with time intervals between individual epochs that range between 1 day and 4$+$ years. To a depth of $m_{AB}$ $\simeq$ 29.5 mag (F606W), we present the discovery of 12 transients and 190 variable candidates. For the variable candidates, we demonstrate that Gaussian statistics are applicable, and estimate that ~80 are false positives. The majority of the transients will be supernovae, although at least two are likely quasars. Most variable candidates are AGN, where we find 0.42% of the general $z$ $<$ 6 field galaxy population to vary at the $~3σ$ level. Based on a 5-year timeframe, this translates into a random supernova areal density of up to ~0.07 transients per arcmin$^2$ (~245 deg$^{-2}$) per epoch, and a variable AGN areal density of ~1.25 variables per arcmin$^2$ (~4500 deg$^{-2}$) to these depths. △ Less

Submitted 2 May, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: 32 pages, 11 figures, 5 tables, 1 Appendix

JWST's PEARLS: 119 multiply imaged galaxies behind MACS0416, lensing properties of caustic crossing galaxies, and the relation between halo mass and number of globular clusters at $z=0.4$

Authors: Jose M. Diego, Nathan J. Adams, Steven Willner, Tom Harvey, Tom Broadhurst, Seth H. Cohen, Rolf A. Jansen, Jake Summers, Rogier A. Windhorst, Jordan C. J. D'Silva, Anton M. Koekemoer, Dan Coe, Christopher J. Conselice, Simon P. Driver, Brenda Frye, Norman A. Grogin, Madeline A. Marshall, Mario Nonino, Rafael Ortiz III, Nor Pirzkal, Aaron Robotham, Russell E. Ryan, Jr., Christopher N. A. Willmer, Haojing Yan, Fengwu Sun , et al. (4 additional authors not shown)

Abstract: We present a new lens model for the $z=0.396$ galaxy cluster MACS J0416.1$-$2403 based on a previously known set of 77 spectroscopically confirmed, multiply imaged galaxies plus an additional set of 42 candidate multiply imaged galaxies from past HST and new JWST data. The new galaxies lack spectroscopic redshifts but have geometric and/or photometric redshift estimates that are presented here. Th… ▽ More We present a new lens model for the $z=0.396$ galaxy cluster MACS J0416.1$-$2403 based on a previously known set of 77 spectroscopically confirmed, multiply imaged galaxies plus an additional set of 42 candidate multiply imaged galaxies from past HST and new JWST data. The new galaxies lack spectroscopic redshifts but have geometric and/or photometric redshift estimates that are presented here. The new model predicts magnifications and time delays for all multiple images. The full set of constraints totals 343, constituting the largest sample of multiple images lensed by a single cluster to date. Caustic-crossing galaxies lensed by this cluster are especially interesting. Some of these galaxies show transient events, most of which are interpreted as micro-lensing of stars at cosmological distances. These caustic-crossing arcs are expected to show similar events in future, deeper JWST observations. We provide time delay and magnification models for all these arcs. The time delays and the magnifications for different arcs are generally anti-correlated, as expected from $N$-body simulations. In the major sub-halos of the cluster, the dark-matter mass from our lens model correlates well with the observed number of globular clusters. This confirms earlier results, derived at lower redshifts, which suggest that globular clusters can be used as powerful mass proxies for the halo masses when lensing constraints are scarce or not available. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 21 pages and 11 figures

JWST's PEARLS: Improved Flux Calibration for NIRCam

Authors: Zhiyuan Ma, Haojing Yan, Bangzheng Sun, Seth H. Cohen, Rolf A. Jansen, Jake Summers, Rogier A. Windhorst, Jordan C. J. D'Silva, Anton M. Koekemoer, Dan Coe, Christopher J. Conselice, Simon P. Driver, Brenda Frye, Norman A. Grogin, Madeline A. Marshall, Mario Nonino, Rafael Ortiz III, Nor Pirzkal, Aaron Robotham, Russell E. Ryan, Jr., Christopher N. A. Willmer, Heidi B. Hammel, Stefanie N. Milam, Nathan J. Adams, Cheng Cheng , et al. (1 additional authors not shown)

Abstract: The Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS), a JWST GTO program, obtained a set of unique NIRCam observations that have enabled us to significantly improve the default photometric calibration across both NIRCam modules. The observations consisted of three epochs of 4-band (F150W, F200W, F356W, and F444W) NIRCam imaging in the Spitzer IRAC Dark Field (IDF). The three… ▽ More The Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS), a JWST GTO program, obtained a set of unique NIRCam observations that have enabled us to significantly improve the default photometric calibration across both NIRCam modules. The observations consisted of three epochs of 4-band (F150W, F200W, F356W, and F444W) NIRCam imaging in the Spitzer IRAC Dark Field (IDF). The three epochs were six months apart and spanned the full duration of Cycle 1. As the IDF is in the JWST continuous viewing zone, we were able to design the observations such that the two modules of NIRCam, modules A and B, were flipped by 180 degrees and completely overlapped each other's footprints in alternate epochs. We were therefore able to directly compare the photometry of the same objects observed with different modules and detectors, and we found significant photometric residuals up to ~ 0.05 mag in some detectors and filters, for the default version of the calibration files that we used (jwst_1039.pmap). Moreover, there are multiplicative gradients present in the data obtained in the two long-wavelength bands. The problem is less severe in the data reduced using the latest pmap (jwst_1130.pmap as of September 2023), but it is still present, and is non-negligible. We provide a recipe to correct for this systematic effect to bring the two modules onto a more consistent calibration, to a photometric precision better than ~ 0.02 mag. △ Less

Submitted 16 January, 2024; v1 submitted 22 November, 2023; originally announced November 2023.

Comments: 12 pages, 8 figures, 3 tables. Accepted to PASP

ProPane: Image Warping with Fire

Authors: A. S. G. Robotham, R. Tobar, S. Bellstedt, S. Casura, R. H. W. Cook, J. C. J. D'Silva, L. J. Davies, S. P. Driver, J. Li, L. P. Garate-Nuñez

Abstract: In this paper we introduce the software package ProPane, written for the R data analysis language. ProPane combines the full range of wcslib projections with the C++ image manipulation routines provided by the CImg library. ProPane offers routines for image warping and combining (including stacking), and various related tasks such as image alignment tweaking and pixel masking. It can stack an effe… ▽ More In this paper we introduce the software package ProPane, written for the R data analysis language. ProPane combines the full range of wcslib projections with the C++ image manipulation routines provided by the CImg library. ProPane offers routines for image warping and combining (including stacking), and various related tasks such as image alignment tweaking and pixel masking. It can stack an effectively unlimited number of target frames using multiple parallel cores, and offers threading for many lower level routines. It has been used for a number of current and upcoming large surveys, and we present a range of its capabilities and features. ProPane is already available under a permissive open-source LGPL-3 license at github.com/asgr/ProPane (DOI: 10.5281/zenodo.10057053). △ Less

Submitted 1 February, 2024; v1 submitted 3 November, 2023; originally announced November 2023.

Comments: 19 pages, 17 figures, 5 tables, accepted to MNRAS

JWST NIRCam Photometry: A Study of Globular Clusters Surrounding Bright Elliptical Galaxy VV 191a at z=0.0513

Authors: Jessica M. Berkheimer, Timothy Carleton, Rogier A. Windhorst, William C. Keel, Benne W. Holwerda, Mario Nonino, Seth H. Cohen, Rolf A. Jansen, Dan Coe, Christopher J. Conselice, Simon P. Driver, Brenda L. Frye, Norman A. Grogin, Anton M. Koekemoer, Ray Lucas, Madeline A. Marshall, Nor Pirzkal, Clayton Robertson, Aaron Robotham, Russell E. Ryan Jr., Brent M. Smith, Jake Summers, Scott Tompkins, Christopher N. A. Willmer, Haojing Yan

Abstract: James Webb Space Telescope NIRCam images have revealed 154 reliable globular cluster (GC) candidates around the $z = 0.0513$ elliptical galaxy VV~191a after subtracting 34 likely interlopers from background galaxies inside our search area. NIRCam broadband observations are made at 0.9-4.5 $μ$m using the F090W, F150W, F356W, and F444W filters. Using PSF-matched photometry, the data are analyzed to… ▽ More James Webb Space Telescope NIRCam images have revealed 154 reliable globular cluster (GC) candidates around the $z = 0.0513$ elliptical galaxy VV~191a after subtracting 34 likely interlopers from background galaxies inside our search area. NIRCam broadband observations are made at 0.9-4.5 $μ$m using the F090W, F150W, F356W, and F444W filters. Using PSF-matched photometry, the data are analyzed to present color-magnitude diagrams (CMDs) and color distributions that suggest a relatively uniform population of GCs, except for small fractions of reddest (5-8%) and bluest (2-4%) outliers. GC models in the F090W vs. (F090-F150W) diagram fit the NIRCam data well and show that the majority of GCs detected have a mass of approximately $\sim$$10^{6.5}$$M_{\odot}$, with metallicities [Fe/H] spanning the typical range expected for GCs (-2.5$\le$ [Fe/H]$\le$ 0.5). However, the models predict $\sim$0.3-0.4 mag bluer (F356W-F444W) colors than the NIRCam data for a reasonable range of GC ages, metallicities, and reddening. Although our data does not quite reach the luminosity function turnover, the measured luminosity function is consistent with previous measurements, suggesting an estimated peak at $m_{\rm AB}$$\sim$-9.4 mag, $\pm$0.2 mag in the F090W filter. △ Less

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

Comments: 11 pages, 7 figures

Report number: AAS50418

Star formation and AGN activity 500 Myr after the Big Bang: Insights from JWST

Authors: Jordan C. J. D'Silva, Simon P. Driver, Claudia D. P. Lagos, Aaron S. G. Robotham, Jake Summers, Rogier A. Windhorst

Abstract: We consider the effect of including an active galactic nuclei (AGN) component when fitting spectral energy distributions of 109 spectroscopically confirmed $z\approx 3.5-12.5$ galaxies with JWST. Remarkably, we find that the resulting cosmic star formation history is $\approx 0.4$ dex lower at $z\gtrsim 9.5$ when an AGN component is included in the fitting. This alleviates previously reported exce… ▽ More We consider the effect of including an active galactic nuclei (AGN) component when fitting spectral energy distributions of 109 spectroscopically confirmed $z\approx 3.5-12.5$ galaxies with JWST. Remarkably, we find that the resulting cosmic star formation history is $\approx 0.4$ dex lower at $z\gtrsim 9.5$ when an AGN component is included in the fitting. This alleviates previously reported excess star formation at $z\gtrsim 9.5$ compared to models based on typical baryon conversion efficiencies inside dark matter halos. We find that the individual stellar masses and star formation rates can be as much as $\approx 4$ dex lower when fitting with an AGN component. These results highlight the importance of considering both stellar mass assembly and supermassive black hole growth when interpreting the light distributions of among the first galaxies to ever exist. △ Less

Submitted 11 December, 2023; v1 submitted 4 October, 2023; originally announced October 2023.

Comments: 15 pages, 5 figures. Accepted for publication in ApJL after addressing reviewer comments

The Hyper Suprime-Cam extended Point Spread Functions and applications

Authors: L. P. Garate-Nuñez, A. S. G. Robotham, S. Bellstedt, L. J. M. Davies, C. Martínez-Lombilla

Abstract: We present extended point spread function (PSF) models for the Hyper Suprime-Cam Subaru Strategic Program Public Data Release 3 (HSC-SSP PDR3) in all $\textit{g,r,i,Z}$ and $\textit{Y}$-bands. Due to its 8.2m primary mirror and long exposure periods, HSC combines deep images with wide-field coverage. Both properties make HSC one of the most suitable observing facilities for low surface brightness… ▽ More We present extended point spread function (PSF) models for the Hyper Suprime-Cam Subaru Strategic Program Public Data Release 3 (HSC-SSP PDR3) in all $\textit{g,r,i,Z}$ and $\textit{Y}$-bands. Due to its 8.2m primary mirror and long exposure periods, HSC combines deep images with wide-field coverage. Both properties make HSC one of the most suitable observing facilities for low surface brightness (LSB) studies, which are particularly sensitive to the PSF. By applying a median stacking technique of point-like sources with different brightness, we show how to construct the HSC-SSP PDR3 PSF models to an extent of R $\sim$ 5.6 arcmin. These models are appropriate for the HSC-PDR3 intermediate-state data which do not have applied the final aggressive background subtraction. The intermediate-state data is especially stored for users interested in large extended objects, where our new PSFs provide them with a crucial tool to characterise LSB properties at large angles. We demonstrate that our HSC PSFs behave reasonably in two scenarios. In the first one, we generate 2-D models of a bright star, showing no evidence of residual structures across the five bands. In the second scenario, we recreate the PSF-scattered light on mock images with special consideration of the effect of this additional flux on LSB measurements. We find that, despite the well-behaved nature of the HSC-PDR3 PSFs, there is a non-negligible impact on the faint light present in the mock images. This impact could lead to incorrect LSB measurements if a proper star subtraction is not applied. △ Less

Submitted 16 May, 2024; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: 16 pages, 12 figures, accepted in MNRAS May 2024

PEARLS: A Potentially Isolated Quiescent Dwarf Galaxy with a TRGB Distance of 30 Mpc

Authors: Timothy Carleton, Timothy Ellsworth-Bowers, Rogier A. Windhorst, Seth H. Cohen, Christopher J. Conselice, Jose M. Diego, Adi Zitrin, Haylee N. Archer, Isabel McIntyre, Patrick Kamieneski, Rolf A. Jansen, Jake Summers, Jordan C. J. D'Silva, Anton M. Koekemoer, Dan Coe, Simon P. Driver, Brenda Frye, Norman A. Grogin, Madeline A. Marshall, Mario Nonino, Nor Pirzkal, Aaron Robotham, Russell E. Ryan, Jr., Rafael Ortiz III, Scott Tompkins , et al. (3 additional authors not shown)

Abstract: A wealth of observations have long suggested that the vast majority of isolated classical dwarf galaxies ($M_*=10^7$-$10^9$ M$_\odot$) are currently star-forming. However, recent observations of the large abundance of "Ultra-Diffuse Galaxies" beyond the reach of previous large spectroscopic surveys suggest that our understanding of the dwarf galaxy population may be incomplete. Here we report the… ▽ More A wealth of observations have long suggested that the vast majority of isolated classical dwarf galaxies ($M_*=10^7$-$10^9$ M$_\odot$) are currently star-forming. However, recent observations of the large abundance of "Ultra-Diffuse Galaxies" beyond the reach of previous large spectroscopic surveys suggest that our understanding of the dwarf galaxy population may be incomplete. Here we report the serendipitous discovery of an isolated quiescent dwarf galaxy in the nearby Universe, which was imaged as part of the PEARLS GTO program. Remarkably, individual red-giant branch stars are visible in this near-IR imaging, suggesting a distance of $30\pm4$ Mpc, and a wealth of archival photometry point to an sSFR of $2\times10^{-11}$ yr$^{-1}$ and SFR of $4\times10^{-4}$ M$_\odot$ yr$^{-1}$. Spectra obtained with the Lowell Discovery Telescope find a recessional velocity consistent with the Hubble Flow and ${>}1500$ km/s separated from the nearest massive galaxy in SDSS, suggesting that this galaxy was either quenched from internal mechanisms or had a very high-velocity ($>1000$ km/s) interaction with a nearby massive galaxy in the past. This analysis highlights the possibility that many nearby quiescent dwarf galaxies are waiting to be discovered and that JWST has the potential to resolve them. △ Less

Submitted 4 January, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: Accepted to ApJ Letters

PEARLS: JWST counterparts of micro-Jy radio sources in the Time Domain Field

Authors: S. P. Willner, H. B. Gim, M. del Carmen Polletta, S. H. Cohen, C. N. A. Willmer, X. Zhao, J. C. J. D'Silva, R. A. Jansen, A. M. Koekemoer, J. Summers, R. A. Windhorst, D. Coe, C. J. Conselice, S. P. Driver, B. Frye, N. A. Grogin, M. A. Marshall, M. Nonino, R. Ortiz III, N. Pirzkal, A. Robotham, M. J. Rutkowski, R. E. Ryan, Jr., S. Tompkins, H. Yan , et al. (16 additional authors not shown)

Abstract: The Time Domain Field (TDF) near the North Ecliptic Pole in JWST's continuous-viewing zone will become a premier "blank field" for extragalactic science. JWST/NIRCam data in a 16 arcmin$^2$ portion of the TDF identify 4.4 $μ$m counterparts for 62 of 63 3 GHz sources with S(3 GHz) > 5 μJy. The one unidentified radio source may be a lobe of a nearby Seyfert galaxy, or it may be an infrared-faint rad… ▽ More The Time Domain Field (TDF) near the North Ecliptic Pole in JWST's continuous-viewing zone will become a premier "blank field" for extragalactic science. JWST/NIRCam data in a 16 arcmin$^2$ portion of the TDF identify 4.4 $μ$m counterparts for 62 of 63 3 GHz sources with S(3 GHz) > 5 μJy. The one unidentified radio source may be a lobe of a nearby Seyfert galaxy, or it may be an infrared-faint radio source. The bulk properties of the radio-host galaxies are consistent with those found by previous work: redshifts range from 0.14 to 4.4 with a median redshift of 1.33. The radio emission arises primarily from star formation in $\sim 2/3$ of the sample and from an active galactic nucleus in $\sim 1/3$, but just over half the sample shows evidence for an AGN either in the spectral energy distribution or by radio excess. All but three counterparts are brighter than magnitude 23 AB at 4.4 $μ$m, and the exquisite resolution of JWST identifies correct counterparts for sources for which observations with lower angular resolution would mis-identify a nearby bright source as the counterpart when the correct one is faint and red. Up to 11% of counterparts might have been unidentified or misidentified absent NIRCam observations. △ Less

Submitted 26 September, 2023; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJ. V2 adds an author and some acknowledgments inadvertently omitted

The JWST Discovery of the Triply-imaged Type Ia "Supernova H0pe" and Observations of the Galaxy Cluster PLCK G165.7+67.0

Authors: Brenda L. Frye, Massimo Pascale, Justin Pierel, Wenlei Chen, Nicholas Foo, Reagen Leimbach, Nikhil Garuda, Seth Cohen, Patrick Kamieneski, Rogier Windhorst, Anton M. Koekemoer, Pat Kelly, Jake Summers, Michael Engesser, Daizhong Liu, Lukas Furtak, Maria Polletta, Kevin Harrington, Steve Willner, Jose M. Diego, Rolf Jansen, Dan Coe, Christopher J. Conselice, Liang Dai, Herve Dole , et al. (17 additional authors not shown)

Abstract: A Type Ia supernova (SN) at $z=1.78$ was discovered in James Webb Space Telescope Near Infrared Camera imaging of the galaxy cluster PLCK G165.7+67.0 (G165; $z = 0.35$). The SN is situated 1.5-2 kpc from the host-galaxy nucleus and appears in three different locations as a result of gravitational lensing by G165. These data can yield a value for Hubble's constant using time delays from this multip… ▽ More A Type Ia supernova (SN) at $z=1.78$ was discovered in James Webb Space Telescope Near Infrared Camera imaging of the galaxy cluster PLCK G165.7+67.0 (G165; $z = 0.35$). The SN is situated 1.5-2 kpc from the host-galaxy nucleus and appears in three different locations as a result of gravitational lensing by G165. These data can yield a value for Hubble's constant using time delays from this multiply-imaged SN Ia that we call "SN H0pe." Over the cluster, we identified 21 image multiplicities, confirmed five of them using the Near-Infrared Spectrograph, and constructed a new lens model that gives a total mass within 600 kpc of ($2.6 \pm 0.3) \times 10^{14}$ $M_{\odot}$. The photometry uncovered a galaxy overdensity coincident with the SN host galaxy. NIRSpec confirmed six member galaxies, four of which surround the SN host galaxy with relative velocity $\lesssim$900 km s$^{-1}$ and projected physical extent $\lesssim$33 kpc. This compact galaxy group is dominated by the SN host galaxy, which has a stellar mass of $(5.0 \pm 0.1) \times 10^{11}$ $M_{\odot}$. The group members have specific star-formation rates of 2-260 Gyr$^{-1}$ derived from the H$α$-line fluxes corrected for stellar absorption, dust extinction, and slit losses. Another group centered on a strongly-lensed dusty star forming galaxy is at $z=2.24$. The total (unobscured and obscured) SFR of this second galaxy group is estimated to be ($\gtrsim$100 $M_{\odot}$ yr$^{-1}$), which translates to a supernova rate of $\sim$1 SNe yr$^{-1}$, suggesting that regular monitoring of this cluster may yield additional SNe. △ Less

Submitted 4 December, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: 29 pages, Accepted to ApJ on November 24, 2023

Quenching massive galaxies across cosmic time with the semi-analytic model SHARK v2.0

Authors: Claudia D. P. Lagos, Matias Bravo, Rodrigo Tobar, Danail Obreschkow, Chris Power, Aaron S. G. Robotham, Katy L. Proctor, Samuel Hansen, Angel Chandro-Gomez, Julian Carrivick

Abstract: We introduce version 2.0 of the SHARK semi-analytic model of galaxy formation after many improvements to the physics included. The most significant being: (i) a model describing the exchange of angular momentum (AM) between the interstellar medium and stars; (ii) a new active galactic nuclei feedback model which has two modes, a wind and a jet mode, with the jet mode tied to the jet energy product… ▽ More We introduce version 2.0 of the SHARK semi-analytic model of galaxy formation after many improvements to the physics included. The most significant being: (i) a model describing the exchange of angular momentum (AM) between the interstellar medium and stars; (ii) a new active galactic nuclei feedback model which has two modes, a wind and a jet mode, with the jet mode tied to the jet energy production; (iii) a model tracking the development of black hole (BH) spins; (iv) more sophisticated modelling of environmental effects on satellite galaxies; and (v) automatic parameter exploration using Particle Swarm Optimisation. We focus on two timely research topics: the structural properties of galaxies and the quenching of massive galaxies. For the former, SHARK v2.0 is capable of producing a more realistic stellar size-mass relation with a plateau marking the transition from disk- to bulge-dominated galaxies, and scaling relations between specific AM and mass that agree well with observations. For the quenching of massive galaxies, SHARK v2.0 produces massive galaxies that are more quenched than the previous version, reproducing well the observed relations between star formation rate (SFR) and stellar mass, and specific SFR and BH mass at $z=0$. SHARK v2.0 produces a number density of massive-quiescent galaxies >1dex higher than the previous version, in good agreement with JWST observations at $z\le 5$; predicts a stellar mass function of passive galaxies in reasonably good agreement with observations at $0.5<z<5$; and environmental quenching to already be effective at $z=5$. △ Less

Submitted 15 April, 2024; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS (change of terminology from radio, QSO mode to jet and wind mode, and a new section 5.3 discussing results in the context of other simulations). Supplementary material with additional comparisons with observations can be found here https://clagos.com/files/Shark_v2_SupplementaryMaterial.pdf

EPOCHS IX. When cosmic dawn breaks: Evidence for evolved stellar populations in $7 < z < 12$ galaxies from PEARLS GTO and public NIRCam imaging

Authors: James A. A. Trussler, Christopher J. Conselice, Nathan Adams, Duncan Austin, Leonardo Ferreira, Tom Harvey, Qiong Li, Aswin P. Vijayan, Stephen M. Wilkins, Rogier A. Windhorst, Rachana Bhatawdekar, Cheng Cheng, Dan Coe, Seth H. Cohen, Simon P. Driver, Brenda Frye, Norman A. Grogin, Nimish Hathi, Rolf A. Jansen, Anton Koekemoer, Madeline A. Marshall, Mario Nonino, Rafael Ortiz, Nor Pirzkal, Aaron Robotham , et al. (6 additional authors not shown)

Abstract: The presence of evolved stars in high-redshift galaxies can place valuable indirect constraints on the onset of star formation in the Universe. Thus we use PEARLS GTO and public NIRCam photometric data to search for Balmer-break candidate galaxies at $7 < z < 12$. We find that our Balmer-break candidates at $z \sim 10.5$ tend to be older (115 Myr), have lower inferred [O III] + H$β$ emission line… ▽ More The presence of evolved stars in high-redshift galaxies can place valuable indirect constraints on the onset of star formation in the Universe. Thus we use PEARLS GTO and public NIRCam photometric data to search for Balmer-break candidate galaxies at $7 < z < 12$. We find that our Balmer-break candidates at $z \sim 10.5$ tend to be older (115 Myr), have lower inferred [O III] + H$β$ emission line equivalent widths (120 Å), have lower specific star formation rates (6 Gyr$^{-1}$) and redder UV slopes ($β= -1.8$) than our control sample of galaxies. However, these trends all become less strong at $z \sim 8$, where the F444W filter now probes the strong rest-frame optical emission lines, thus providing additional constraints on the current star formation activity of these galaxies. Indeed, the bursty nature of Epoch of Reionisation galaxies can lead to a disconnect between their current SED profiles and their more extended star-formation histories. We discuss how strong emission lines, the cumulative effect of weak emission lines, dusty continua and AGN can all contribute to the photometric excess seen in the rest-frame optical, thus mimicking the signature of a Balmer break. Additional medium-band imaging will thus be essential to more robustly identify Balmer-break galaxies. However, the Balmer break alone cannot serve as a definitive proxy for the stellar age of galaxies, being complexly dependent on the star-formation history. Ultimately, deep NIRSpec continuum spectroscopy and MIRI imaging will provide the strongest indirect constraints on the formation era of the first galaxies in the Universe, thereby revealing when cosmic dawn breaks. △ Less

Submitted 7 March, 2024; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 25 pages, 16 figures, 2 tables. Updated to published version in MNRAS

Journal ref: MNRAS, 527, 11627 (2024)

A search for high-redshift direct-collapse black hole candidates in the PEARLS north ecliptic pole field

Authors: Armin Nabizadeh, Erik Zackrisson, Fabio Pacucci, Peter W. Maksym, Weihui Li, Francesca Civano, Seth H. Cohen, Jordan C. J. D'Silva, Anton M. Koekemoer, Jake Summers, Rogier A. Windhorst, Nathan Adams, Christopher J. Conselice, Dan Coe, Simon P. Driver, Brenda Frye, Norman A. Grogin, Rolf A. Jansen, Madeline A. Marshall, Mario Nonino, Nor Pirzkal, Aaron Robotham, Michael J. Rutkowski, Russell E. Ryan, Jr., Scott Tompkins , et al. (9 additional authors not shown)

Abstract: Direct-collapse black holes (DCBHs) of mass $\sim 10^4$-$10^5 {M}_\odot$ that form in HI-cooling halos in the early Universe are promising progenitors of the $\gtrsim 10^9 {M}_\odot$ supermassive black holes that fuel observed $z \gtrsim 7$ quasars. Efficient accretion of the surrounding gas onto such DCBH seeds may render them sufficiently bright for detection with the JWST up to $z\approx 20$. A… ▽ More Direct-collapse black holes (DCBHs) of mass $\sim 10^4$-$10^5 {M}_\odot$ that form in HI-cooling halos in the early Universe are promising progenitors of the $\gtrsim 10^9 {M}_\odot$ supermassive black holes that fuel observed $z \gtrsim 7$ quasars. Efficient accretion of the surrounding gas onto such DCBH seeds may render them sufficiently bright for detection with the JWST up to $z\approx 20$. Additionally, the very steep and red spectral slope predicted across the $\approx 1$-5 $μ$m wavelength range of the JWST/NIRSpec instrument during their initial growth phase should make them photometrically identifiable up to very high redshifts. In this work, we present a search for such DCBH candidates across the 34 arcmin$^{2}$ in the first two spokes of the JWST cycle-1 PEARLS survey of the north ecliptic pole time-domain field covering eight NIRCam filters down to a maximum depth of $\sim$ 29 AB mag. We identify two objects with spectral energy distributions consistent with the Pacucci et al. (2016) DCBH models. However, we also note that even with data in eight NIRCam filters, objects of this type remain degenerate with dusty galaxies and obscured active galactic nuclei over a wide range of redshifts. Follow-up spectroscopy would be required to pin down the nature of these objects. Based on our sample of DCBH candidates and assumptions on the typical duration of the DCBH steep-slope state, we set a conservative upper limit of $\lesssim 5\times 10^{-4}$ comoving Mpc$^{-3}$ (cMpc$^{-3}$) on the comoving density of host halos capable of hosting DCBHs with spectral energy distributions similar to the Pacucci et al. (2016) models at $z\approx 6$-14. △ Less

Submitted 19 January, 2024; v1 submitted 14 August, 2023; originally announced August 2023.

Comments: 10 pages, 6 figures, 1 table, accepted for publication in A&A

Identifying the discs, bulges, and intra-halo light of simulated galaxies through structural decomposition

Authors: Katy L. Proctor, Claudia del P. Lagos, Aaron D. Ludlow, Aaron S. G. Robotham

Abstract: We perform a structural decomposition of galaxies identified in three cosmological hydrodynamical simulations by applying Gaussian Mixture Models (GMMs) to the kinematics of their stellar particles. We study the resulting disc, bulge, and intra-halo light (IHL) components of galaxies whose host dark matter haloes have virial masses in the range $M_{200}=10^{11}$-- $10^{15}\,{\rm M_\odot}$. Our dec… ▽ More We perform a structural decomposition of galaxies identified in three cosmological hydrodynamical simulations by applying Gaussian Mixture Models (GMMs) to the kinematics of their stellar particles. We study the resulting disc, bulge, and intra-halo light (IHL) components of galaxies whose host dark matter haloes have virial masses in the range $M_{200}=10^{11}$-- $10^{15}\,{\rm M_\odot}$. Our decomposition technique isolates galactic discs whose mass fractions, $f_{\rm disc}$, correlate strongly with common alternative morphology indicators; for example, $f_{\rm disc}$ is approximately equal to $κ_{\rm co}$, the fraction of stellar kinetic energy in co-rotation. The primary aim of our study, however, is to characterise the IHL of galaxies in a consistent manner and over a broad mass range, and to analyse its properties from the scale of galactic stellar haloes up to the intra-cluster light. Our results imply that the IHL fraction, $f_{\rm IHL}$, has appreciable scatter and is strongly correlated with galaxy morphology: at fixed stellar mass, the IHL of disc galaxies is typically older and less massive than that of spheroids. Above $M_{200}\approx 10^{13}\,{\rm M_\odot}$, we find, on average, $f_{\rm IHL}\approx 0.45$, albeit with considerable scatter. The transition radius beyond which the IHL dominates the stellar mass of a galaxy is roughly $30\,{\rm kpc}$ for $M_{200}\lesssim 10^{12.8}\,{\rm M_\odot}$, but increases strongly towards higher masses. However, we find that no alternative IHL definitions -- whether based on the ex-situ stellar mass, or the stellar mass outside a spherical aperture -- reproduce our dynamically-defined IHL masses. △ Less

Submitted 6 December, 2023; v1 submitted 24 July, 2023; originally announced July 2023.

Comments: 16 pages, 11 figures. MNRAS published

JWST's PEARLS: Mothra, a new kaiju star at z=2.091 extremely magnified by MACS0416, and implications for dark matter models

Authors: J. M. Diego, Bangzheng Sun, Haojing Yan, Lukas J. Furtak, Erik Zackrisson, Liang Dai, Patrick Kelly, Mario Nonino, Nathan Adams, Ashish K. Meena, S. P. Willner, Adi Zitrin, Seth H. Cohen, Jordan C. J. D Silva, Rolf A. Jansen, Jake Summers, Rogier A. Windhorst, Dan Coe, Christopher J. Conselice, Simon P. Driver, Brenda Frye, Norman A. Grogin, Anton M. Koekemoer, Madeline A. Marshall, Nor Pirzkal , et al. (6 additional authors not shown)

Abstract: We report the discovery of Mothra, an extremely magnified monster star, likely a binary system of two supergiant stars, in one of the strongly lensed galaxies behind the galaxy cluster MACS0416. The star is in a galaxy with spectroscopic redshift $z=2.091$ in a portion of the galaxy that is parsecs away from the cluster caustic. The binary star is observed only on the side of the critical curve wi… ▽ More We report the discovery of Mothra, an extremely magnified monster star, likely a binary system of two supergiant stars, in one of the strongly lensed galaxies behind the galaxy cluster MACS0416. The star is in a galaxy with spectroscopic redshift $z=2.091$ in a portion of the galaxy that is parsecs away from the cluster caustic. The binary star is observed only on the side of the critical curve with negative parity but has been detectable for at least eight years, implying the presence of a small lensing perturber. Microlenses alone cannot explain the earlier observations of this object made with the Hubble Space Telescope. A larger perturber with a mass of at least $10^4$\,\Msun\ offers a more satisfactory explanation. Based on the lack of perturbation on other nearby sources in the same arc, the maximum mass of the perturber is $M< 2.5\times10^6$\,\Msun, making it the smallest substructure constrained by lensing above redshift 0.3. The existence of this millilens is fully consistent with the expectations from the standard cold dark matter model. On the other hand, the existence of such small substructure in a cluster environment has implications for other dark matter models. In particular, warm dark matter models with particle masses below 8.7\,keV are excluded by our observations. Similarly, axion dark matter models are consistent with the observations only if the axion mass is in the range $0.5\times10^{-22}\, {\rm eV} < m_a < 5\times10^{-22}\, {\rm eV}$. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 26 pages and 27 figures

Journal ref: A&A 679, A31 (2023)

JWST's PEARLS: Transients in the MACS J0416.1-2403 Field

Authors: Haojing Yan, Zhiyuan Ma, Bangzheng Sun, Lifan Wang, Patrick Kelly, Jose M. Diego, Seth H. Cohen, Rogier A. Windhorst, Rolf A. Jansen, Norman A. Grogin, John F. Beacom, Christopher J. Conselice, Simon P. Driver, Brenda Frye, Dan Coe, Madeline A. Marshall, Anton Koekemoer, Christopher N. A. Willmer, Aaron Robotham, Jordan C. J. D'Silva, Jake Summers, Mario Nonino, Nor Pirzkal, Russell E. Ryan, Jr., Rafael Ortiz III , et al. (5 additional authors not shown)

Abstract: With its unprecedented sensitivity and spatial resolution, the James Webb Space Telescope (JWST) has opened a new window for time-domain discoveries in the infrared. Here we report observations in the only field that has received four epochs (spanning 126 days) of JWST NIRCam observations in Cycle 1. This field is towards MACS J0416.1-2403, which is a rich galaxy cluster at redshift z=0.4 and is o… ▽ More With its unprecedented sensitivity and spatial resolution, the James Webb Space Telescope (JWST) has opened a new window for time-domain discoveries in the infrared. Here we report observations in the only field that has received four epochs (spanning 126 days) of JWST NIRCam observations in Cycle 1. This field is towards MACS J0416.1-2403, which is a rich galaxy cluster at redshift z=0.4 and is one of the Hubble Frontier Fields. We have discovered 14 transients from these data. Twelve of these transients happened in three galaxies (with z=0.94, 1.01, and 2.091) crossing a lensing caustic of the cluster,and these transients are highly magnified by gravitational lensing. These 12 transients are likely of similar nature to those previously reported based on the Hubble Space Telescope (HST) data in this field, i.e., individual stars in the highly magnified arcs. However, these twelve could not have been found by HST because they are too red and too faint. The other two transients are associated with background galaxies (z=2.205 and 0.7093) that are only moderately magnified, and they are likely supernovae. They indicate a de-magnified supernova surface density, when monitored at a time cadence of a few months to a ~3--4 micron survey limit of AB ~ 28.5 mag, of ~0.5 per sq. arcmin integrated to z ~ 2. This survey depth is beyond the capability of HST but can be easily reached by JWST. △ Less

Submitted 11 October, 2023; v1 submitted 14 July, 2023; originally announced July 2023.

Comments: ApJS accepted

EPOCHS VII: Discovery of high redshift ($6.5 < z < 12$) AGN candidates in JWST ERO and PEARLS data

Authors: Ignas Juodžbalis, Christopher J. Conselice, Maitrayee Singh, Nathan Adams, Katherine Ormerod, Thomas Harvey, Duncan Austin, Marta Volonteri, Seth H. Cohen, Rolf A. Jansen, Jake Summers, Rogier A. Windhorst, Jordan C. J. D'Silva, Anton M. Koekemoer, Dan Coe, Simon P. Driver, Brenda Frye, Norman A. Grogin, Madeline A. Marshall, Mario Nonino, Nor Pirzkal, Aaron Robotham, Russell E. Ryan, Jr., Rafael Ortiz III, Scott Tompkins , et al. (2 additional authors not shown)

Abstract: We present an analysis of a sample of robust high redshift galaxies selected photometrically from the `blank' fields of the Prime Extragalactic Areas for Reionization Science (PEARLS) survey and Early Release Observations (ERO) data of the James Webb Space Telescope (JWST) with the aim of selecting candidate high redshift active galactic nuclei (AGN). Sources were identified from the parent sample… ▽ More We present an analysis of a sample of robust high redshift galaxies selected photometrically from the `blank' fields of the Prime Extragalactic Areas for Reionization Science (PEARLS) survey and Early Release Observations (ERO) data of the James Webb Space Telescope (JWST) with the aim of selecting candidate high redshift active galactic nuclei (AGN). Sources were identified from the parent sample using a threefold selection procedure, which includes spectral energy distribution (SED) fitting to identify sources that are best fitted by AGN SED templates, a further selection based on the relative performance of AGN and non-AGN models, and finally morphological fitting to identify compact sources of emission, resulting in a purity-oriented procedure. Using this procedure, we identify a sample of nine AGN candidates at $6.5 < z < 12$, from which we constrain their physical properties as well as measure a lower bound on the AGN fraction in this redshift range of $5 \pm 1$\%. As this is an extreme lower limit due to our focus on purity and our SEDs being calibrated for unobscured Type 1 AGN, this demonstrates that AGN are perhaps quite common at this early epoch. The rest-frame UV colors of our candidate objects suggest that these systems are potentially candidate obese black hole galaxies (OBG), or AGN with very little galaxy component. We also investigate emission from our sample sources from fields overlapping with Chandra and VLA surveys, allowing us to place X-ray and 3 GHz radio detection limits on our candidates. Of note is a $z = 11.9$ candidate source exhibiting an abrupt morphological shift in the reddest band as compared to the bluer bands, indicating a potential merger or an unusually strong outflow. △ Less

Submitted 3 August, 2023; v1 submitted 14 July, 2023; originally announced July 2023.

Comments: Accepted by MNRAS, 12 pages, 11 figures, updated to the final accepted version

Resolving cosmic star formation histories of present-day bulges, disks, and spheroids with ProFuse

Authors: Sabine Bellstedt, Aaron S. G. Robotham, Simon P. Driver, Claudia del P. Lagos, Luke J. M. Davies, Robin H. W. Cook

Abstract: We present the first look at star formation histories of galaxy components using ProFuse, a new technique to model the 2D distribution of light across multiple wavelengths using simultaneous spectral and spatial fitting of purely imaging data. We present a number of methods to classify galaxies structurally/morphologically, showing the similarities and discrepancies between these schemes. We show… ▽ More We present the first look at star formation histories of galaxy components using ProFuse, a new technique to model the 2D distribution of light across multiple wavelengths using simultaneous spectral and spatial fitting of purely imaging data. We present a number of methods to classify galaxies structurally/morphologically, showing the similarities and discrepancies between these schemes. We show the variation in component-wise mass functions that can occur simply due to the use of a different classification method, which is most dramatic in separating bulges and spheroids. Rather than identifying the best-performing scheme, we use the spread of classifications to quantify uncertainty in our results. We study the cosmic star formation history (CSFH), forensically derived using ProFuse with a sample of ~7,000 galaxies from the Galaxy And Mass Assembly (GAMA) survey. Remarkably, the forensic CSFH recovered via both our method (ProFuse) and traditional SED fitting (ProSpect) are not only exactly consistent with each other over the past 8 Gyr, but also with the in-situ CSFH measured using ProSpect. Furthermore, we separate the CSFH by contributions from spheroids, bulges and disks. While the vast majority (70%) of present-day star formation takes place in the disk population, we show that 50% of the stars that formed at cosmic noon (8-12 Gyr ago) now reside in spheroids, and present-day bulges are composed of stars that were primarily formed in the very early Universe, with half their stars already formed ~12 Gyr ago. △ Less

Submitted 18 June, 2024; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 27 pages, 20 figures, accepted for publication in MNRAS

SimSpin v2.6.0 -- Constructing synthetic spectral IFU cubes for comparison with observational surveys

Authors: K. E. Harborne, A. Serene, E. J. A. Davies, C. Derkenne, S. Vaughan, A. I. Burdon, C. del P. Lagos, R. McDermid, S. O'Toole, C. Power, A. S. G. Robotham, G. Santucci, R. Tobar

Abstract: In this work, we present a methodology and a corresponding code-base for constructing mock integral field spectrograph (IFS) observations of simulated galaxies in a consistent and reproducible way. Such methods are necessary to improve the collaboration and comparison of observation and theory results, and accelerate our understanding of how the kinematics of galaxies evolve over time. This code,… ▽ More In this work, we present a methodology and a corresponding code-base for constructing mock integral field spectrograph (IFS) observations of simulated galaxies in a consistent and reproducible way. Such methods are necessary to improve the collaboration and comparison of observation and theory results, and accelerate our understanding of how the kinematics of galaxies evolve over time. This code, SimSpin, is an open-source package written in R, but also with an API interface such that the code can be interacted with in any coding language. Documentation and individual examples can be found at the open-source website connected to the online repository. SimSpin is already being utilised by international IFS collaborations, including SAMI and MAGPI, for generating comparable data sets from a diverse suite of cosmological hydrodynamical simulations. △ Less

Submitted 30 August, 2023; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 20 pages, 15 figures, 2 tables. Accepted for publication in PASA. 30/08/23

GAMA/DEVILS: Cosmic star formation and AGN activity over 12.5 billion years

Authors: Jordan C. J. D'Silva, Simon P. Driver, Claudia D. P. Lagos, Aaron S. G. Robotham, Sabine Bellstedt, Luke J. M. Davies, Jessica E. Thorne, Joss Bland-Hawthorn, Matias Bravo, Benne Holwerda, Steven Phillipps, Nick Seymour, Malgorzata Siudek, Rogier A. Windhorst

Abstract: We use the Galaxy and Mass Assembly (GAMA) and the Deep Extragalactic Visible Legacy Survey (DEVILS) observational data sets to calculate the cosmic star formation rate (SFR) and active galactic nuclei (AGN) bolometric luminosity history (CSFH/CAGNH) over the last 12.5 billion years. SFRs and AGN bolometric luminosities were derived using the spectral energy distribution fitting code ProSpect, whi… ▽ More We use the Galaxy and Mass Assembly (GAMA) and the Deep Extragalactic Visible Legacy Survey (DEVILS) observational data sets to calculate the cosmic star formation rate (SFR) and active galactic nuclei (AGN) bolometric luminosity history (CSFH/CAGNH) over the last 12.5 billion years. SFRs and AGN bolometric luminosities were derived using the spectral energy distribution fitting code ProSpect, which includes an AGN prescription to self consistently model the contribution from both AGN and stellar emission to the observed rest-frame ultra-violet to far-infrared photometry. We find that both the CSFH and CAGNH evolve similarly, rising in the early Universe up to a peak at look-back time $\approx 10$~Gyr ($z \approx 2$), before declining toward the present day. The key result of this work is that we find the ratio of CAGNH to CSFH has been flat ($\approx 10^{42.5}\mathrm{erg \, s^{-1}M_{\odot}^{-1}yr}$) for $11$~Gyr up to the present day, indicating that star formation and AGN activity have been coeval over this time period. We find that the stellar masses of the galaxies that contribute most to the CSFH and CAGNH are similar, implying a common cause, which is likely gas inflow. The depletion of the gas supply suppresses cosmic star formation and AGN activity equivalently to ensure that they have experienced similar declines over the last 10 Gyr. These results are an important milestone for reconciling the role of star formation and AGN activity in the life cycle of galaxies. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 16 pages, 10 figures. Figures 9 and 10 are the main results. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Hidden giants in JWST's PEARLS: An ultra-massive z=4.26 sub-millimeter galaxy that is invisible to HST

Authors: Ian Smail, Ugne Dudzeviciute, Mark Gurwell, Giovanni G. Fazio, S. P. Willner, A. M. Swinbank, Vinodiran Arumugam, Jake Summers, Seth H. Cohen, Rolf A. Jansen, Rogier A. Windhorst, Ashish Meena, Adi Zitrin, William C. Keel, Dan Coe, Christopher J. Conselice, Jordan C. J. D'Silva, Simon P. Driver, Brenda Frye, Norman A. Grogin, Anton M. Koekemoer, Madeline A. Marshall, Mario Nonino, Nor Pirzkal, Aaron Robotham , et al. (10 additional authors not shown)

Abstract: We present a multi-wavelength analysis using SMA, JCMT, NOEMA, JWST, HST, and SST of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A1489. These SMA-located sources both lie at z=4.26 and have bright dust continuum emission, but 850.2 is a UV-detected Lyman-break galaxy, while 850.1 is undetected at <2um, even with deep JWST/NIRCam observations. We… ▽ More We present a multi-wavelength analysis using SMA, JCMT, NOEMA, JWST, HST, and SST of two dusty strongly star-forming galaxies, 850.1 and 850.2, seen through the massive cluster lens A1489. These SMA-located sources both lie at z=4.26 and have bright dust continuum emission, but 850.2 is a UV-detected Lyman-break galaxy, while 850.1 is undetected at <2um, even with deep JWST/NIRCam observations. We investigate their stellar, ISM, and dynamical properties, including a pixel-level SED analysis to derive sub-kpc-resolution stellar-mass and Av maps. We find that 850.1 is one of the most massive and highly obscured, Av~5, galaxies known at z>4 with M*~10^11.8 Mo (likely forming at z>6), and 850.2 is one of the least massive and least obscured, Av~1, members of the z>4 dusty star-forming population. The diversity of these two dust-mass-selected galaxies illustrates the incompleteness of galaxy surveys at z>3-4 based on imaging at <2um, the longest wavelengths feasible from HST or the ground. The resolved mass map of 850.1 shows a compact stellar mass distribution, Re(mass)~1kpc, but its expected evolution to z~1.5 and then z~0 matches both the properties of massive, quiescent galaxies at z~1.5 and ultra-massive early-type galaxies at z~0. We suggest that 850.1 is the central galaxy of a group in which 850.2 is a satellite that will likely merge in the near future. The stellar morphology of 850.1 shows arms and a linear bar feature which we link to the active dynamical environment it resides within. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: Submitted to ApJ, comments welcome!

Magellanic System Stars Identified in SMACS J0723.3-7327 James Webb Space Telescope Early Release Observations Images

Authors: Jake Summers, Rogier A. Windhorst, Seth H. Cohen, Rolf A. Jansen, Timothy Carleton, Patrick S. Kamieneski, Benne W. Holwerda, Christopher J. Conselice, Nathan J. Adams, Brenda Frye, Jose M. Diego, Christopher N. A. Willmer, Rafael Ortiz III, Cheng Cheng, Alex Pigarelli, Aaron Robotham, Jordan C. J. D'Silva, Scott Tompkins, Simon P. Driver, Haojing Yan, Dan Coe, Norman Grogin, Anton Koekemoer, Madeline A. Marshall, Nor Pirzkal , et al. (1 additional authors not shown)

Abstract: We identify 71 distant stars in JWST/NIRCam ERO images of the field of galaxy cluster SMACS J0723.3-7327 (SMACS 0723). Given the relatively small ($\sim$$10^{\circ}$) angular separation between SMACS 0723 and the Large Magellanic Cloud, it is likely that these stars are associated with the LMC outskirts or Leading Arm. This is further bolstered by a spectral energy distribution analysis, which sug… ▽ More We identify 71 distant stars in JWST/NIRCam ERO images of the field of galaxy cluster SMACS J0723.3-7327 (SMACS 0723). Given the relatively small ($\sim$$10^{\circ}$) angular separation between SMACS 0723 and the Large Magellanic Cloud, it is likely that these stars are associated with the LMC outskirts or Leading Arm. This is further bolstered by a spectral energy distribution analysis, which suggests an excess of stars at a physical distance of $40-100$ kpc, consistent with being associated with or located behind the Magellanic system. In particular, we find that the overall surface density of stars brighter than 27.0 mag in the field of SMACS 0723 is $\sim$2.3 times that of stars in a blank field with similar galactic latitude (the North Ecliptic Pole Time Domain Field), and that the density of stars in the SMACS 0723 field with SED-derived distances consistent with the Magellanic system is $\sim$6.1 times larger than that of the blank field. The candidate stars at these distances are consistent with a stellar population at the same distance modulus with [Fe/H] $= -1.0$ and an age of $\sim$$5.0$ Gyr. On the assumption that all of the 71 stars are associated with the LMC, then the stellar density of the LMC at the location of the SMACS 0723 field is $\sim$$740$ stars kpc$^{-3}$, which helps trace the density of stars in the LMC outskirts. △ Less

Submitted 1 February, 2024; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: Published in ApJ

Journal ref: ApJ 958 108 (2023)

Spectroscopy of the Supernova H0pe Host Galaxy at Redshift 1.78

Authors: M. Polletta, M. Nonino, B. Frye, A. Gargiulo, S. Bisogni, N. Garuda, D. Thompson, M. Lehnert, M. Pascale, S. P. Willner, P. Kamieneski, R. Leimbach, C. Cheng, D. Coe, S. H. Cohen, C. J. Conselice, L. Dai, J. Diego, H. Dole, S. P. Driver, J. C. J. D'Silva, A. Fontana, N. Foo, L. J. Furtak, N. A. Grogin , et al. (20 additional authors not shown)

Abstract: Supernova (SN) H0pe was discovered as a new transient in James Webb Space Telescope (JWST) NIRCam images of the galaxy cluster PLCK G165.7+67.0 taken as part of the "Prime Extragalactic Areas for Reionization and Lensing Science" (PEARLS) JWST GTO program (# 1176) on 2023 March 30 (AstroNote 2023-96; Frye et al. 2023). The transient is a compact source associated with a background galaxy that is s… ▽ More Supernova (SN) H0pe was discovered as a new transient in James Webb Space Telescope (JWST) NIRCam images of the galaxy cluster PLCK G165.7+67.0 taken as part of the "Prime Extragalactic Areas for Reionization and Lensing Science" (PEARLS) JWST GTO program (# 1176) on 2023 March 30 (AstroNote 2023-96; Frye et al. 2023). The transient is a compact source associated with a background galaxy that is stretched and triply-imaged by the cluster's strong gravitational lensing. This paper reports spectra in the 950-1370 nm observer frame of two of the galaxy's images obtained with Large Binocular Telescope (LBT) Utility Camera in the Infrared (LUCI) in longslit mode two weeks after the \JWST\ observations. The individual average spectra show the [OII] doublet and the Balmer and 4000 Angstrom breaks at redshift z=1.783+/-0.002. The CIGALE best-fit model of the spectral energy distribution indicates that SN H0pe's host galaxy is massive (Mstar~6x10^10 Msun after correcting for a magnification factor ~7) with a predominant intermediate age (~2 Gyr) stellar population, moderate extinction, and a magnification-corrected star formation rate ~13 Msun/yr, consistent with being below the main sequence of star formation. These properties suggest that H0pe might be a type Ia SN. Additional observations of SN H0pe and its host recently carried out with JWST (JWST-DD-4446; PI: B. Frye) will be able to both determine the SN classification and confirm its association with the galaxy analyzed in this work. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 6 pages, 4 figures, Letter accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 675, L4 (2023)

Galaxy And Mass Assembly (GAMA): The group HI mass as a function of halo mass

Authors: Ajay Dev, Simon P. Driver, Martin Meyer, Sambit Roychowdhury, Jonghwan Rhee, Adam R. H. Stevens, Claudia del P. Lagos, Joss Bland-Hawthorn, Barbara Catinella, A. M. Hopkins, Jonathan Loveday, Danail Obreschkow, Steven Phillipps, Aaron S. G. Robotham

Abstract: We determine the atomic hydrogen (HI) to halo mass relation (HIHM) using Arecibo Legacy Fast ALFA survey HI data at the location of optically selected groups from the Galaxy and Mass Assembly (GAMA) survey. We make direct HI detections for 37 GAMA groups. Using HI group spectral stacking of 345 groups, we study the group HI content as function of halo mass across a halo mass range of… ▽ More We determine the atomic hydrogen (HI) to halo mass relation (HIHM) using Arecibo Legacy Fast ALFA survey HI data at the location of optically selected groups from the Galaxy and Mass Assembly (GAMA) survey. We make direct HI detections for 37 GAMA groups. Using HI group spectral stacking of 345 groups, we study the group HI content as function of halo mass across a halo mass range of $10^{11} - 10^{14.7}\text{ M}_\odot$. We also correct our results for Eddington bias. We find that the group HI mass generally rises as a function of halo mass from $1.3\%$ of the halo mass at $10^{11.6} \text{M}_\odot$ to $0.4\%$ at $10^{13.7} \text{M}_\odot$ with some indication of flattening towards the high-mass end. Despite the differences in optical survey limits, group catalogues, and halo mass estimation methods, our results are consistent with previous group HI-stacking studies. Our results are also consistent with mock observations from SHARK and IllustrisTNG. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: Accepted in MNRAS; 18 pages, 12 figures

The Long and the Short of It: The Benefits and Leverage of Ultraviolet-Radio Galaxy Fitting

Authors: Jessica E. Thorne, Aaron S. G. Robotham, Sabine Bellstedt, Luke J. M. Davies

Abstract: Traditionally, the far ultraviolet (FUV) to far-infrared (FIR) and radio spectral energy distributions (SEDs) of galaxies have been considered separately despite the common physical process shaping them. In this work, we explore the utility of simultaneously fitting FUV-radio SEDs using an extended version of the ProSpect SED fitting code considering contributions from both free-free and synchrotr… ▽ More Traditionally, the far ultraviolet (FUV) to far-infrared (FIR) and radio spectral energy distributions (SEDs) of galaxies have been considered separately despite the common physical process shaping them. In this work, we explore the utility of simultaneously fitting FUV-radio SEDs using an extended version of the ProSpect SED fitting code considering contributions from both free-free and synchrotron emission. We use a small sample of galaxies from the Deep Extragalactic VIsible Legacy Survey (DEVILS) and the Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel (KINGFISH) where high-quality and robust FUV-radio data are available to provide an ideal sample for testing a radio extension of ProSpect. As the parameterisation of the radio extension links the radio continuum to the FIR emission, we explore the benefit of using radio continuum measurements as a constraint on the energy balance between dust attenuation and emission. We find that for situations where MIR-FIR photometry is unavailable, including a 1.4 GHz continuum measurement allows for improved accuracy in recovered star formation rates and dust luminosities of galaxies reducing the median uncertainty by 0.1 and 0.2 dex respectively. We also demonstrate that incorporating 3 and 10 GHz measurements allows for further constraint on the energy balance and therefore the star formation rate and dust luminosity. This demonstrates the advantage of extending FUV-FIR SED fitting techniques to radio frequencies, especially as we move into an era where FIR surveys will remain limited and radio data become abundant (i.e. with the SKA and precursors). △ Less

Submitted 14 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS

Dynamic Wisp Removal in JWST NIRCam Images

Authors: A. S. G. Robotham, J. C. J. D'Silva, R. A. Windhorst, R. A. Jansen, J. Summers, S. P. Driver, C. N. A. Wilmer, S. Bellstedt

Abstract: The James Webb Space Telescope (JWST) near-infrared camera (NIRCam) has been found to exhibit serious wisp-like structures in four of its eight short-wavelength detectors. The exact structure and strength of these wisps is highly variable with the position and orientation of JWST, so the use of static templates is non-optimal. Here we investigate a dynamic strategy to mitigate these wisps using lo… ▽ More The James Webb Space Telescope (JWST) near-infrared camera (NIRCam) has been found to exhibit serious wisp-like structures in four of its eight short-wavelength detectors. The exact structure and strength of these wisps is highly variable with the position and orientation of JWST, so the use of static templates is non-optimal. Here we investigate a dynamic strategy to mitigate these wisps using long-wavelength reference images. Based on a suite of experiments where we embed a worst-case scenario median-stacked wisp into wisp-free images, we define suitable parameters for our wisp removal strategy. Using this setup we re-process wisp-affected public Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS) data in the North Ecliptic Pole Time Domain Field (NEP-TDF), resulting in significant visual improvement in our detector frames and reduced noise in the final stacked images. △ Less

Submitted 13 July, 2023; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: 20 pages, 13 figures, accepted to PASP

Hydra II: Characterisation of Aegean, Caesar, ProFound, PyBDSF, and Selavy source finders

Authors: M. M. Boyce, A. M. Hopkins, S. Riggi, L. Rudnick, M. Ramsay, C. L. Hale, J. Marvil, M. Whiting, P. Venkataraman, C. P. O'Dea, S. A. Baum, Y. A. Gordon, A. N. Vantyghem, M. Dionyssiou, H. Andernach, J. D. Collier, J. English, B. S. Koribalski, D. Leahy, M. J. Michałowski, S. Safi-Harb, M. Vaccari, E. Alexander, M. Cowley, A. D. Kapinska , et al. (2 additional authors not shown)

Abstract: We present a comparison between the performance of a selection of source finders using a new software tool called Hydra. The companion paper, Paper~I, introduced the Hydra tool and demonstrated its performance using simulated data. Here we apply Hydra to assess the performance of different source finders by analysing real observational data taken from the Evolutionary Map of the Universe (EMU) Pil… ▽ More We present a comparison between the performance of a selection of source finders using a new software tool called Hydra. The companion paper, Paper~I, introduced the Hydra tool and demonstrated its performance using simulated data. Here we apply Hydra to assess the performance of different source finders by analysing real observational data taken from the Evolutionary Map of the Universe (EMU) Pilot Survey. EMU is a wide-field radio continuum survey whose primary goal is to make a deep ($20μ$Jy/beam RMS noise), intermediate angular resolution ($15^{\prime\prime}$), 1\,GHz survey of the entire sky south of $+30^{\circ}$ declination, and expecting to detect and catalogue up to 40 million sources. With the main EMU survey expected to begin in 2022 it is highly desirable to understand the performance of radio image source finder software and to identify an approach that optimises source detection capabilities. Hydra has been developed to refine this process, as well as to deliver a range of metrics and source finding data products from multiple source finders. We present the performance of the five source finders tested here in terms of their completeness and reliability statistics, their flux density and source size measurements, and an exploration of case studies to highlight finder-specific limitations. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: Paper accepted for publication in PASA

ACM Class: D.2.8

Hydra I: An extensible multi-source-finder comparison and cataloguing tool

Authors: M. M. Boyce, A. M. Hopkins, S. Riggi, L. Rudnick, M. Ramsay, C. L. Hale, J. Marvil, M. Whiting, P. Venkataraman, C. P. O'Dea, S. A. Baum, Y. A. Gordon, A. N. Vantyghem, M. Dionyssiou, H. Andernach, J. D. Collier, J. English, B. S. Koribalski, D. Leahy, M. J. Michałowski, S. Safi-Harb, M. Vaccari, E. Alexander, M. Cowley, A. D. Kapinska , et al. (2 additional authors not shown)

Abstract: The latest generation of radio surveys are now producing sky survey images containing many millions of radio sources. In this context it is highly desirable to understand the performance of radio image source finder (SF) software and to identify an approach that optimises source detection capabilities. We have created Hydra to be an extensible multi-SF and cataloguing tool that can be used to comp… ▽ More The latest generation of radio surveys are now producing sky survey images containing many millions of radio sources. In this context it is highly desirable to understand the performance of radio image source finder (SF) software and to identify an approach that optimises source detection capabilities. We have created Hydra to be an extensible multi-SF and cataloguing tool that can be used to compare and evaluate different SFs. Hydra, which currently includes the SFs Aegean, Caesar, ProFound, PyBDSF, and Selavy, provides for the addition of new SFs through containerisation and configuration files. The SF input RMS noise and island parameters are optimised to a 90\% ''percentage real detections'' threshold (calculated from the difference between detections in the real and inverted images), to enable comparison between SFs. Hydra provides completeness and reliability diagnostics through observed-deep ($\mathcal{D}$) and generated-shallow ($\mathcal{S}$) images, as well as other statistics. In addition, it has a visual inspection tool for comparing residual images through various selection filters, such as S/N bins in completeness or reliability. The tool allows the user to easily compare and evaluate different SFs in order to choose their desired SF, or a combination thereof. This paper is part one of a two part series. In this paper we introduce the Hydra software suite and validate its $\mathcal{D/S}$ metrics using simulated data. The companion paper demonstrates the utility of Hydra by comparing the performance of SFs using both simulated and real images. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: Paper accepted for publication in PASA

ACM Class: D.2.4; D.2.6; D.2.8; D.2.10; D.2.11; D.2.13

EPOCHS Paper II: The Ultraviolet Luminosity Function from $7.5<z<13.5$ using 180 square arcminutes of deep, blank-fields from the PEARLS Survey and Public JWST data

Authors: Nathan J. Adams, Christopher J. Conselice, Duncan Austin, Thomas Harvey, Leonardo Ferreira, James Trussler, Ignas Juodzbalis, Qiong Li, Rogier Windhorst, Seth H. Cohen, Rolf Jansen, Jake Summers, Scott Tompkins, Simon P. Driver, Aaron Robotham, Jordan C. J. D'Silva, Haojing Yan, Dan Coe, Brenda Frye, Norman A. Grogin, Anton M. Koekemoer, Madeline A. Marshall, Nor Pirzkal, Russell E. Ryan, Jr., W. Peter Maksym , et al. (12 additional authors not shown)

Abstract: We present an analysis of the ultraviolet luminosity function (UV LF) and star formation rate density of distant galaxies ($7.5 < z < 13.5$) in the `blank' fields of the Prime Extragalactic Areas for Reionization Science (PEARLS) survey combined with Early Release Science (ERS) data from the CEERS, GLASS, NGDEEP surveys/fields and the first data release of JADES. We use strict quality cuts on EAZY… ▽ More We present an analysis of the ultraviolet luminosity function (UV LF) and star formation rate density of distant galaxies ($7.5 < z < 13.5$) in the `blank' fields of the Prime Extragalactic Areas for Reionization Science (PEARLS) survey combined with Early Release Science (ERS) data from the CEERS, GLASS, NGDEEP surveys/fields and the first data release of JADES. We use strict quality cuts on EAZY photometric redshifts to obtain a reliable selection and characterisation of high-redshift ($z>6.5$) galaxies from a consistently processed set of deep, near-infrared imaging. Within an area of 180 arcmin$^{2}$, we identify 1046 candidate galaxies at redshifts $z>6.5$ and we use this sample to study the ultraviolet luminosity function (UV LF) in four redshift bins between $7.5<z<13.5$. The measured number density of galaxies at $z=8$ and $z=9$ match those of past observations undertaken by the {\em Hubble Space Telescope} (HST). Our $z=10.5$ measurements lie between early JWST results and past HST results, indicating cosmic variance may be the cause of previous high density measurements. However, number densities of UV luminous galaxies at $z=12.5$ are high compared to predictions from simulations. When examining the star formation rate density of galaxies at this time period, our observations are still largely consistent with a constant star formation efficiency, are slightly lower than previous early estimations using JWST and support galaxy driven reionization at $z\leq8$. △ Less

Submitted 6 March, 2024; v1 submitted 26 April, 2023; originally announced April 2023.

Comments: 28 Pages, 4 Tables, 10 Figures, V2 updated paper: Accepted in ApJ, analysis now includes completed PEARLS NEP-TDF Field and JADES DR1, updated image calibrations for all fields, more detailed completeness analysis conducted. V3 updated paper: corrected minor issues in acknowledgements section

Are JWST/NIRCam color gradients in the lensed z=2.3 dusty star-forming galaxy El Anzuelo due to central dust attenuation or inside-out galaxy growth?

Authors: Patrick S. Kamieneski, Brenda L. Frye, Massimo Pascale, Seth H. Cohen, Rogier A. Windhorst, Rolf A. Jansen, Min S. Yun, Cheng Cheng, Jake S. Summers, Timothy Carleton, Kevin C. Harrington, Jose M. Diego, Haojing Yan, Anton M. Koekemoer, Christopher N. A. Willmer, Andreea Petric, Lukas J. Furtak, Nicholas Foo, Christopher J. Conselice, Dan Coe, Simon P. Driver, Norman A. Grogin, Madeline A. Marshall, Nor Pirzkal, Aaron S. G. Robotham , et al. (2 additional authors not shown)

Abstract: Gradients in the mass-to-light ratio of distant galaxies impede our ability to characterize their size and compactness. The long-wavelength filters of $JWST$'s NIRCam offer a significant step forward. For galaxies at Cosmic Noon ($z\sim2$), this regime corresponds to the rest-frame near-infrared, which is less biased towards young stars and captures emission from the bulk of a galaxy's stellar pop… ▽ More Gradients in the mass-to-light ratio of distant galaxies impede our ability to characterize their size and compactness. The long-wavelength filters of $JWST$'s NIRCam offer a significant step forward. For galaxies at Cosmic Noon ($z\sim2$), this regime corresponds to the rest-frame near-infrared, which is less biased towards young stars and captures emission from the bulk of a galaxy's stellar population. We present an initial analysis of an extraordinary lensed dusty star-forming galaxy (DSFG) at $z=2.3$ behind the $El~Gordo$ cluster ($z=0.87$), named $El~Anzuelo$ ("The Fishhook") after its partial Einstein-ring morphology. The FUV-NIR SED suggests an intrinsic star formation rate of $81^{+7}_{-2}~M_\odot~{\rm yr}^{-1}$ and dust attenuation $A_V\approx 1.6$, in line with other DSFGs on the star-forming main sequence. We develop a parametric lens model to reconstruct the source-plane structure of dust imaged by the Atacama Large Millimeter/submillimeter Array, far-UV to optical light from $Hubble$, and near-IR imaging with 8 filters of $JWST$/NIRCam, as part of the Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS) program. The source-plane half-light radius is remarkably consistent from $\sim 1-4.5~μ$m, despite a clear color gradient where the inferred galaxy center is redder than the outskirts. We interpret this to be the result of both a radially-decreasing gradient in attenuation and substantial spatial offsets between UV- and IR-emitting components. A spatial decomposition of the SED reveals modestly suppressed star formation in the inner kiloparsec, which suggests that we are witnessing the early stages of inside-out quenching. △ Less

Submitted 23 July, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: 29 pages, 11 figures, 5 tables. Accepted for publication in ApJ

PEARLS: Low Stellar Density Galaxies in the El Gordo Cluster Observed with JWST

Authors: Timothy Carleton, Seth H. Cohen, Brenda Frye, Alex Pigarelli, Jiashuo Zhang, Rogier A. Windhorst, Jose M. Diego, Christopher J. Conselice, Cheng Cheng, Simon P. Driver, Nicholas Foo, Rachana A. Bhatawdekar, Patrick Kamieneski, Rolf A. Jansen, Haojing Yan, Jake Summers, Aaron Robotham, Christopher N. A. Willmer, Anton Koekemoer, Scott Tompkins, Dan Coe, Norman Grogin, Madeline A. Marshall, Mario Nonino, Nor Pirzkal , et al. (1 additional authors not shown)

Abstract: A full understanding of how unusually large "Ultra Diffuse Galaxies" (UDGs) fit into our conventional understanding of dwarf galaxies remains elusive, despite the large number of objects identified locally. A natural extension of UDG research is the study of similar galaxies at higher redshift to establish how their properties may evolve over time. However, this has been a challenging task given h… ▽ More A full understanding of how unusually large "Ultra Diffuse Galaxies" (UDGs) fit into our conventional understanding of dwarf galaxies remains elusive, despite the large number of objects identified locally. A natural extension of UDG research is the study of similar galaxies at higher redshift to establish how their properties may evolve over time. However, this has been a challenging task given how severely systematic effects and cosmological surface brightness dimming inhibit our ability to study low-surface brightness galaxies at high-$z$. Here, we present an identification of low stellar surface density galaxies (LDGs), likely the progenitors of local UDGs, at moderate redshift with deep near-IR observations of the El Gordo cluster at $z = 0.87$ with JWST. By stacking 8 NIRCAM filters, we are able to achieve an apparent surface brightness sensitivity of $24.59$ mag arcsec$^{-2}$, faint enough to be complete to the bright end of the LDG population. Our analysis identifies significant differences between this population and local UDGs, such as their color and size distributions, which suggest that UDG progenitors are bluer and more extended at high-$z$ than at $z = 0$. This suggests that multiple mechanisms are responsible for UDG formation and that prolonged transformation of cluster dwarfs is not a primary UDG formation mechanism at high-$z$. Furthermore, we find a slight overabundance of LDGs in El Gordo, and, in contrast to findings in local clusters, our analysis does not show a deficit of LDGs in the center of El Gordo, implying that tidal destruction of LDGs is significant between $z = 0.87$ and $z = 0$. △ Less

Submitted 17 August, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

Comments: Accepted to ApJ

Paper 1: The JWST PEARLS View of the El Gordo Galaxy Cluster and of the Structure It Magnifies

Authors: Brenda L. Frye, Massimo Pascale, Nicholas Foo, Reagen Leimbach, Nikhil Garuda, Paulina Soto Robles, Jake Summers, Carlos Diaz, Patrick Kamieneski, Lukas Furtak, Seth Cohen, Jose Diego, Benjamin Beauchesne, Rogier Windhorst, Steve Willner, Anton M. Koekemoer, Adi Zitrin, Gabriel Caminha, Karina Caputi, Dan Coe, Christopher J. Conselice, Liang Dai, Herve Dole, Simon Driver, Norman Grogin , et al. (19 additional authors not shown)

Abstract: The massive galaxy cluster El Gordo (z=0.87) imprints multitudes of gravitationally lensed arcs onto James Webb Space Telescope (JWST) Near-Infrared Camera (NIRCam) images. Eight bands of NIRCam imaging were obtained in the ``Prime Extragalactic Areas for Reionization and Lensing Science'' (``PEARLS'') program. PSF-matched photometry across Hubble Space Telescope (HST) and NIRCam filters supplies… ▽ More The massive galaxy cluster El Gordo (z=0.87) imprints multitudes of gravitationally lensed arcs onto James Webb Space Telescope (JWST) Near-Infrared Camera (NIRCam) images. Eight bands of NIRCam imaging were obtained in the ``Prime Extragalactic Areas for Reionization and Lensing Science'' (``PEARLS'') program. PSF-matched photometry across Hubble Space Telescope (HST) and NIRCam filters supplies new photometric redshifts. A new light-traces-mass lens model based on 56 image multiplicities identifies the two mass peaks and yields a mass estimate within 500 kpc of ~(7.0 +/- 0.30) x 10^14 Msun. A search for substructure in the 140 cluster members with spectroscopic redshifts confirms the two main mass components. The southeastern mass peak that contains the BCG is more tightly bound than the northwestern one. The virial mass within 1.7 Mpc is (5.1 +/- 0.60) x 10^14 Msun, lower than the lensing mass. A significant transverse velocity component could mean the virial mass is underestimated. We contribute one new member to the previously known z=4.32 galaxy group. Intrinsic (delensed) positions of the five secure group members span a physical extent of ~60 kpc. Thirteen additional candidates selected by spectroscopic/photometric constraints are small and faint with a mean intrinsic luminosity ~2.2 mag fainter than L*. NIRCam imaging admits a fairly wide range of brightnesses and morphologies for the group members, suggesting a more diverse galaxy population in this galaxy overdensity. △ Less

Submitted 3 June, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 24 pages, accepted by ApJ

The X-ray invisible Universe. A look into the halos undetected by eROSITA

Authors: P. Popesso, A. Biviano, E. Bulbul, A. Merloni, J. Comparat, N. Clerc, Z. Igo, A. Liu, S. Driver, M. Salvato, M. Brusa, Y. E. Bahar, N. Malavasi, V. Ghirardini, G. Ponti, A. Robotham, J. Liske, S. Grandis

Abstract: The paper presents the analysis of optically selected GAMA groups and clusters in the SRG/eROSITA X-ray map of eFEDS (eROSITA Final Equatorial Depth Survey), in the halo mass range $10^{13}-5{\times}10^{14}$ $M_{\odot}$ and at $z < 0.2$. All X-ray detections have a clear GAMA counterpart, but most of the GAMA groups in the halo mass range $10^{13}-10^{14}$ $M_{\odot}$ remain undetected. We compare… ▽ More The paper presents the analysis of optically selected GAMA groups and clusters in the SRG/eROSITA X-ray map of eFEDS (eROSITA Final Equatorial Depth Survey), in the halo mass range $10^{13}-5{\times}10^{14}$ $M_{\odot}$ and at $z < 0.2$. All X-ray detections have a clear GAMA counterpart, but most of the GAMA groups in the halo mass range $10^{13}-10^{14}$ $M_{\odot}$ remain undetected. We compare the X-ray surface brightness profiles of the eROSITA detected groups with the mean stacked profile of the undetected low-mass halos at fixed halo mass. Overall, we find that the undetected groups exhibit less concentrated X-ray surface brightness, dark matter, and galaxy distributions with respect to the X-ray detected halos. The mean gas mass fraction profiles are consistent in the two samples within 1.5$σ$, indicating that the gas follows the dark matter profile. The low mass concentration and the magnitude gap indicate that these systems are young. They reside with a higher probability in filaments while X-ray detected groups favor the nodes of the Cosmic Web. Because of the lower central emission, the undetected systems tend to be X-ray under-luminous at fixed halo mass and to lie below the $L_X-M_{halo}$ relation. Interestingly, the X-ray detected systems inhabiting the nodes scatter the less around the relation, while those in filaments tend to lie below it. We do not observe any strong relationship between the system X-ray appearance and the AGN activity. We cannot exclude the role of the past AGN feedback in affecting the gas distribution over the halo lifetime. However, the data suggests that the observed differences might be related to the halo assembly bias. △ Less

Submitted 20 October, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: 16 pages, 14 figures, Accepted for publication in MNRAS

A New WISE Calibration of Stellar Mass

Authors: T. H. Jarrett, M. E. Cluver, Edward N. Taylor, Sabine Bellstedt, A. S. G Robotham, H. F. M. Yao

Abstract: We derive new empirical scaling relations between WISE mid-infrared galaxy photometry and well-determined stellar masses from SED modeling of a suite of optical-infrared photometry provided by the DR4 Catalogue of the GAMA-KiDS-VIKING survey of the southern G23 field. The mid-infrared source extraction and characterization are drawn from the WISE Extended Source Catalogue (WXSC) and the archival A… ▽ More We derive new empirical scaling relations between WISE mid-infrared galaxy photometry and well-determined stellar masses from SED modeling of a suite of optical-infrared photometry provided by the DR4 Catalogue of the GAMA-KiDS-VIKING survey of the southern G23 field. The mid-infrared source extraction and characterization are drawn from the WISE Extended Source Catalogue (WXSC) and the archival ALLWISE catalog, combining both resolved and compact galaxies in the G23 sample to a redshift of 0.15. Three scaling relations are derived: W1 3.4 micron luminosity versus stellar mass, and WISE W1-W2, W1-W3 colors versus mass-to-light ratio (sensitive to a variety of galaxy types from passive to star-forming). For each galaxy in the sample, we then derive the combined stellar mass from these scaling relations, producing Mstellar estimates with better than $\sim$25-30% accuracy for galaxies with $>$10$^{9}$ Msolar and $<$40 - 50% for lower luminosity dwarf galaxies. We also provide simple prescriptions for rest-frame corrections and estimating stellar masses using only the W1 flux and the W1-W2 color, making stellar masses more accessible to users of the WISE data. Given a redshift or distance, these new scaling relations will enable stellar mass estimates for any galaxy in the sky detected by WISE with high fidelity across a range of mass-to-light. △ Less

Submitted 14 January, 2023; originally announced January 2023.

Comments: Accepted for publication in the Astrophysical Journal (ApJ)

Galaxy quenching timescales from a forensic reconstruction of their colour evolution

Authors: Matías Bravo, Aaron S. G. Robotham, Claudia del P. Lagos, Luke J. M. Davies, Sabine Bellstedt, Jessica E. Thorne

Abstract: The timescales on which galaxies move out of the blue cloud to the red sequence ($τ^{}_\mathrm{Q}$) provide insight into the mechanisms driving quenching. Here, we build upon previous work, where we showcased a method to reconstruct the colour evolution of observed low-redshift galaxies from the Galaxy And Mass Assembly (GAMA) survey based on spectral energy distribution (SED) fitting with ProSpec… ▽ More The timescales on which galaxies move out of the blue cloud to the red sequence ($τ^{}_\mathrm{Q}$) provide insight into the mechanisms driving quenching. Here, we build upon previous work, where we showcased a method to reconstruct the colour evolution of observed low-redshift galaxies from the Galaxy And Mass Assembly (GAMA) survey based on spectral energy distribution (SED) fitting with ProSpect, together with a statistically-driven definition for the blue and red populations. We also use the predicted colour evolution from the SHARK semi-analytic model, combined with SED fits of our simulated galaxy sample, to study the accuracy of the measured $τ^{}_\mathrm{Q}$ and gain physical insight into the colour evolution of galaxies. In this work, we measure $τ^{}_\mathrm{Q}$ in a consistent approach for both observations and simulations. After accounting for selection bias, we find evidence for an increase in $τ^{}_\mathrm{Q}$ in GAMA as a function of cosmic time (from $τ^{}_\mathrm{Q}\sim1$ Gyr to $τ^{}_\mathrm{Q}\sim2$ Gyr in the lapse of $\sim4$ Gyr), but not in SHARK ($τ^{}_\mathrm{Q}\lesssim1$ Gyr). Our observations and simulations disagree on the effect of stellar mass, with GAMA showing massive galaxies transitioning faster, but is the opposite in SHARK. We find that environment only impacts galaxies below $\sim10^{10}$ M$_\odot$ in GAMA, with satellites having shorter $τ^{}_\mathrm{Q}$ than centrals by $\sim0.4$ Gyr, with SHARK only in qualitative agreement. Finally, we compare to previous literature, finding consistency with timescales in the order of couple Gyr, but with several differences that we discuss. △ Less

Submitted 23 March, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 17 pages, 14 figures. Submitted to MNRAS. Updated to reflect changes addressing the referee's comments