Search

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., Comment: Submitted to AAS Journals, 37 pages

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., Comment: 11 pages, 8 figures. Submitted to MNRAS

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

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., Comment: 10 pages, 5 figures. Accepted for publication in PASA

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., Comment: Submitted to MNRAS; 17 pages, 11 figures

New JWST/NIRCam wide-field slitless spectroscopy provides redshifts for four z>8 galaxies located behind the lensing cluster MACS J0416.1-2403. Two of them, "Y1" and "JD", have previously reported spectroscopic redshifts based on ALMA measurements of [OIII] 88 $\mu$m and/or [CII] 157.7 $\mu$m lines. Y1 is a merging system of three components, and the existing redshift z=8.31 is confirmed. However, JD is at z=8.34 instead of the previously claimed z=9.28. JD's close companion, "JD-N", which was a previously discovered z>8 candidate, is now identified at the same redshift as JD. JD and JD-N form an interacting pair. A new candidate at z>8, "f090d_018", is also confirmed and is at z=8.49. These four objects are likely part of an overdensity that signposts a large structure extending ~165 kpc in projected distance and ~48.7 Mpc in radial distance. They are magnified by less than one magnitude and have intrinsic $M_{UV}$ ranging from -19.57 to -20.83 mag. Their spectral energy distributions show that the galaxies are all very young with ages ~ 4-18 Myr and stellar masses about $10^{7-8}$ ${\rm M_\odot}$. These infant galaxies have very different star formation rates ranging from a few to over a hundred $\rm{M_\odot}$ yr$^{-1}$, but only two of them (JD and f090d_018) have blue rest-frame UV slopes $\beta<-2.0$ indicative of a high Lyman-continuum photon escape fraction that could contribute significantly to the cosmic hydrogen-reionizing background. Interestingly, these two galaxies are the least massive and least active ones among the four. The other two systems have much flatter UV slopes largely because of their high dust extinction ($A_{\rm V}$=0.9-1.0 mag). Their 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., Comment: 18 pages, 6 figures, after addressing the referee report

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., Comment: 23 pages, 13 figures, 2 tables. Accepted for publication in MNRAS

We present an analysis of rest-frame UV continuum slopes, $\beta$, using a sample of 1011 galaxies at $6.5

The first public 0.9-4.4{\mu}m NIRCam images of the North Ecliptic Pole (NEP) Time Domain Field (TDF) uncovered galaxies displaying point-source features in their cores as seen in the longer wavelength filters. We visually identified a sample of 66 galaxies (~1 galaxy per arcmin2) with point-like cores and have modeled their two-dimensional light profiles with GalFit, identifying 16 galactic nuclei with measurable point-source components. GalFit suggests the visual sample is a mix of both compact stellar bulge and point-source galaxy cores. This core classification is complemented by spectral energy distribution (SED) modeling to infer the sample's active galactic nucleus (AGN) and host-galaxy parameters. For galaxies with measurable point-source components, the median fractional AGN contribution to their 0.1-30.0{\mu}m flux is 0.44, and 14/16 are color-classified AGN. We conclude that near-infrared point-source galaxy cores are signatures of AGN. In addition, we define an automated sample-selection criterion to identify these point-source features. These criteria can be used in other extant and future NIRCam images to streamline the search for galaxies with unresolved IR-luminous AGN. The James Webb Space Telescope's superb angular resolution and sensitivity at infrared wavelengths is resurrecting the morphological identification of AGN., Comment: 17 pages, 9 figures, 2 tables

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 $\Delta 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., Comment: 47 pages, 21 figures, 5 tables. Submitted to ApJ, comments welcome!

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, $\alpha$, the concentration, $c$, the gas profile's inner logarithmic density slope, $\varepsilon$, the dilution, $d$, and the gas fraction, $\eta$, 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., Comment: This paper has been accepted for publication with PASA

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., Comment: 16 pages, 12 figures, 5 tables. Accepted for publication in PASA

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., Comment: 21 pages and 11 figures. Accepted for publication in A&A

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., Comment: 12 pages, 8 figures, 3 tables. Accepted to PASP

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 $\mu$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., Comment: 11 pages, 7 figures

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., Comment: 15 pages, 5 figures. Accepted for publication in ApJL after addressing reviewer comments

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., Comment: Accepted to ApJ Letters

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$\alpha$-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., Comment: 29 pages, Accepted to ApJ on November 24, 2023

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

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$\beta$ emission line equivalent widths (120 \r{A}), have lower specific star formation rates (6 Gyr$^{-1}$) and redder UV slopes ($\beta = -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., Comment: 25 pages, 16 figures, 2 tables. Updated to published version in MNRAS

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 $\mu$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., Comment: 10 pages, 6 figures, 1 table, accepted for publication in A&A

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 $\kappa_{{\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., Comment: 16 pages, 11 figures. MNRAS published

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}$., Comment: 26 pages and 27 figures

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., Comment: ApJS accepted

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., Comment: Accepted by MNRAS, 12 pages, 11 figures, updated to the final accepted version

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., Comment: 27 pages, 20 figures, accepted for publication in MNRAS

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., Comment: Submitted to ApJ, comments welcome!

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., Comment: 16 pages, 10 figures. Figures 9 and 10 are the main results. Accepted for publication in Monthly Notices of the Royal Astronomical Society

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., Comment: Published in ApJ

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., Comment: Accepted in MNRAS; 18 pages, 12 figures

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)., Comment: Accepted for publication in MNRAS

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

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~\mu$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., Comment: 29 pages, 11 figures, 5 tables. Accepted for publication in ApJ

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$., Comment: Accepted to ApJ

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., Comment: 24 pages, accepted by ApJ

The timescales on which galaxies move out of the blue cloud to the red sequence ($\tau^{}_\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 $\tau^{}_\mathrm{Q}$ and gain physical insight into the colour evolution of galaxies. In this work, we measure $\tau^{}_\mathrm{Q}$ in a consistent approach for both observations and simulations. After accounting for selection bias, we find evidence for an increase in $\tau^{}_\mathrm{Q}$ in GAMA as a function of cosmic time (from $\tau^{}_\mathrm{Q}\sim1$ Gyr to $\tau^{}_\mathrm{Q}\sim2$ Gyr in the lapse of $\sim4$ Gyr), but not in SHARK ($\tau^{}_\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 $\tau^{}_\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., Comment: 17 pages, 14 figures. Submitted to MNRAS. Updated to reflect changes addressing the referee's comments

We present a revised measurement of the extra-galactic background light (EBL) at radio frequencies based on a near complete compendium of radio source counts. We present the radio-EBL at 150 MHz, 325 MHz, 610 MHz, 1.4 GHz, 3 GHz, 5 GHz, and 8.4 GHz. In all cases the contribution to the radio-EBL, per decade of flux, exhibits a two-humped distribution well matched to the AGN and star-forming galaxy (SFG) populations, and with each population contributing roughly equal energy. Only at 3 GHz are the source count contributions to the EBL fully convergent, and hence we report empirical lower limits to the radio-EBL in the remaining bands. Adopting predictions from the SHARK semi-analytic model for the form of the SFG population, we can fit the fainter source counts providing measurements of the total contribution to the radio-EBL for the SFG and the AGN populations separately. This constitutes an empirically constrained model-dependent measurement for the SFG contribution, but a fully empirical measurement of the AGN contribution. Using the {\sc ProSpect} spectral energy distribution code we can model the UV-optical-infrared-mm-radio SFG EBL at all frequencies from the cosmic star-formation history and the adoption of a Chabrier initial mass function. However, significant discrepancy remains ($5\times$) between our source-count estimates of the radio-EBL and the direct measurements reported from the ARCADE-2 experiment. We can rule out a significant missing discrete source radio population and suggest that the cause of the high ARCADE-2 radio-EBL values may need to be sought either in the foreground subtraction or as a yet unknown diffuse component in the radio sky.

We present the first JWST observations of the $z=4.11$ luminous radio galaxy TN J1338-1942, obtained as part of the ``Prime Extragalactic Areas for Reionization and Lensing Science'' (``PEARLS'') project. Our NIRCam observations, designed to probe the key rest-frame optical continuum and emission line features at this redshift, enable resolved spectral energy distribution modelling that incorporates both a range of stellar population assumptions and radiative shock models. With an estimated stellar mass of $\log_{10}(M/\text{M}_{\odot}) \sim 10.9$, TN J1338--1942 is confirmed to be one of the most massive galaxies known at this epoch. Our observations also reveal extremely high equivalent-width nebular emission coincident with the luminous AGN jets that is best fit by radiative shocks surrounded by extensive recent star-formation. We estimate the total star-formation rate (SFR) could be as high as $\sim1600\,\text{M}_{\odot}\,\text{yr}^{-1}$, with the SFR that we attribute to the jet induced burst conservatively $\gtrsim500\,\text{M}_{\odot}\,\text{yr}^{-1}$. The mass-weighted age of the star-formation, $t_{\text{mass}} <4$ Myr, is consistent with the likely age of the jets responsible for the triggered activity and significantly younger than that measured in the core of the host galaxy. The extreme scale of the potential jet-triggered star-formation activity indicates the potential importance of positive AGN feedback in the earliest stages of massive galaxy formation, with our observations also illustrating the extraordinary prospects for detailed studies of high-redshift galaxies with JWST., Comment: 17 pages, 15 figures, Accepted for publication in MNRAS after minor revisions

We present a pilot study to assess the potential of Hyper Suprime-Cam Public Data Release 2 (HSC-PDR2) images for the analysis of extended faint structures within groups of galaxies. We examine the intra-group light (IGL) of the group 400138 ($M_{\rm{dyn}}= 1.3 \pm 0.5 \times 10^{13} $M$_{\odot}$, $z\sim 0.2$) from the Galaxy And Mass Assembly (GAMA) survey using Hyper-Suprime Cam Subaru Strategic Program Public Data Release 2 (HSC-PDR2) images in $g$, $r$, and $i$ bands. We present the most extended IGL measurement to date, reaching down to $\mu_{g}^{\rm{lim}}=30.76$ mag arcsec$^{-2}$ ($3 \sigma$; $10 \times 10$ arcsec$^{2}$) at a semi-major axis of 275 kpc. The IGL shows mean colour values of $g-i=0.92$, $g-r=0.60$, and $r-i=0.32$ ($\pm$0.01). The IGL stellar populations are younger ($2-2.5$ Gyr) and less metal-rich ([Fe/H] $ \sim -$0.4) than those of the host group galaxies. We find a range of IGL fractions as a function of total group luminosity of $\sim 2-36 \%$ depending on the definition of IGL, with larger fractions the bluer the observation wavelength. The early-type to late-type galaxy ratio suggests that 400138 is a more evolved group, dominated by ETGs, and the IGL fraction agrees with that of other similarly evolved groups. These results are consistent with tidal stripping of the outer parts of Milky Way-like galaxies as the main driver of the IGL build-up. This is supported by the detection of substructure in the IGL towards the galaxy member 1660615 suggesting a recent interaction ($<1$ Gyr ago) of that galaxy with the core of the group., Comment: 19 pages, 8 figures, 2 tables, accepted for publication in MNRAS

Gas-phase metallicities of galaxies are typically measured through auroral or nebular emission lines, but metallicity also leaves an imprint on the overall spectral energy distribution (SED) of a galaxy and can be estimated through SED fitting. We use the ProSpect SED fitting code with a flexible parametric star formation history and an evolving metallicity history to self-consistently measure metallicities, stellar mass, and other galaxy properties for $\sim90\,000$ galaxies from the Deep Extragalactic VIsible Legacy Survey (DEVILS) and Galaxy and Mass Assembly (GAMA) survey. We use these to trace the evolution of the mass-metallicity relation (MZR) and show that the MZR only evolves in normalisation by $\sim0.1\,$dex at stellar mass $M_\star = 10^{10.5}\,M_\odot$. We find no difference in the MZR between galaxies with and without SED evidence of active galactic nuclei emission at low redshifts ($z<0.3$). Our results suggest an anti-correlation between metallicity and star formation activity at fixed stellar mass for galaxies with $M_\star > 10^{10.5}\,M_\odot$ for $z<0.3$. Using the star formation histories extracted using ProSpect we explore higher-order correlations of the MZR with properties of the star formation history including age, width, and shape. We find that at a given stellar mass, galaxies with higher metallicities formed most of their mass over shorter timescales, and before their peak star formation rate. This work highlights the value of exploring the connection of a galaxy's current gas-phase metallicity to its star formation history in order to understand the physical processes shaping the MZR., Comment: 25 pages, 20 figures. Accepted for publication in MNRAS

We present early science results from Deep Investigation of Neutral Gas Origins (DINGO), an HI survey using the Australian Square Kilometre Array Pathfinder (ASKAP). Using ASKAP sub-arrays available during its commissioning phase, DINGO early science data were taken over $\sim$ 60 deg$^{2}$ of the Galaxy And Mass Assembly (GAMA) 23 h region with 35.5 hr integration time. We make direct detections of six known and one new sources at $z < 0.01$. Using HI spectral stacking, we investigate the HI gas content of galaxies at $0.04 < z< 0.09$ for different galaxy colours. The results show that galaxy morphology based on optical colour is strongly linked to HI gas properties. To examine environmental impacts on the HI gas content of galaxies, three sub-samples are made based on the GAMA group catalogue. The average HI mass of group central galaxies is larger than those of satellite and isolated galaxies, but with a lower HI gas fraction. We derive a variety of HI scaling relations for physical properties of our sample, including stellar mass, stellar mass surface density, $NUV-r$ colour, specific star formation rate, and halo mass. We find that the derived HI scaling relations are comparable to other published results, with consistent trends also observed to $\sim$0.5 dex lower limits in stellar mass and stellar surface density. The cosmic HI densities derived from our data are consistent with other published values at similar redshifts. DINGO early science highlights the power of HI spectral stacking techniques with ASKAP., Comment: 27 pages, 25 figures, 10 tables, accepted for publication in MNRAS

We report the results of James Webb Space Telescope/NIRCam observations of 19 (sub)millimeter (submm/mm) sources detected by the Atacama Large Millimeter Array (ALMA). The accurate ALMA positions allowed unambiguous identifications of their NIRCam counterparts. Taking gravitational lensing into account, these represent 16 distinct galaxies in three fields and constitute the largest sample of its kind to date. The counterparts' spectral energy distributions from rest-frame ultraviolet to near infrared provide photometric redshifts ($110^{10.5}$ Msol), which are similar to sub-millimeter galaxy (SMG) hosts studied previously. However, our sample is fainter in submm/mm than the classic SMG samples are, and our sources exhibit a wider range of properties. They have dust-embedded star-formation rates as low as 10 Msol yr$^{-1}$, and the sources populate both the star-forming main sequence and the quiescent categories. The deep NIRCam data allow us to study the rest-frame near-IR morphologies. Excluding two multiply imaged systems and one quasar, the majority of the remaining sources are disk-like and show either little or no disturbance. This suggests that secular growth is a potential route for the assembly of high-mass disk galaxies. While a few hosts have large disks, the majority have small disks (median half-mass radius of 1.6 kpc). At this time, it is unclear whether this is due to the prevalence of small disks at these redshifts or some unknown selection effects of deep ALMA observations. A larger sample of ALMA sources with NIRCam observations will be able to address this question., Comment: 19 pages, 11 figures, accepted by ApJL

The diffuse, unresolved sky provides most of the photons that the Hubble Space Telescope (HST) receives, yet remains poorly understood. HST Archival Legacy program SKYSURF aims to measure the 0.2-1.6 $\mu$m sky surface brightness (sky-SB) from over 140,000 HST images. We describe a sky-SB measurement algorithm designed for SKYSURF that is able to recover the input sky-SB from simulated images to within 1% uncertainty. We present our sky-SB measurements estimated using this algorithm on the entire SKYSURF database. Comparing our sky-SB spectral energy distribution (SED) to measurements from the literature shows general agreements. Our SKYSURF SED also reveals a possible dependence on Sun angle, indicating either non-isotropic scattering of solar photons off interplanetary dust or an additional component to Zodiacal Light. Finally, we update Diffuse Light limits in the near-IR based on the methods from Carleton et al. (2022), with values of 0.009 MJy sr$^{-1}$ (22 nW m$^{-2}$ sr$^{-1}$) at 1.25 $\mu$m, 0.015 MJy sr$^{-1}$ (32 nW m$^{-2}$ sr$^{-1}$) at 1.4 $\mu$m, and 0.013 MJy sr$^{-1}$ (25 nW m$^{-2}$ sr$^{-1}$) at 1.6 $\mu$m. These estimates provide the most stringent all-sky constraints to date in this wavelength range. SKYSURF sky-SB measurements are made public on the official SKYSURF website and will be used to constrain Diffuse Light in future papers., Comment: Revised based on helpful comments from the reviewer, and accepted to AJ on April 12th, 2023. Main paper: 18 pages, 9 figures, 4 tables. Appendices: 16 pages, 10 figures, 1 table. Main results shown in Figure 7 and Table 4

The first JWST data on the massive colliding cluster El Gordo confirm 23 known families of multiply lensed images and identify 8 new members of these families. Based on these families, which have been confirmed spectroscopically by MUSE, we derived an initial lens model. This model guided the identification of 37 additional families of multiply lensed galaxies, among which 28 are entirely new systems, and 9 were previously known. The initial lens model determined geometric redshifts for the 37 new systems. The geometric redshifts agree reasonably well with spectroscopic or photometric redshifts when those are available. The geometric redshifts enable two additional models that include all 60 families of multiply lensed galaxies spanning a redshift range $20.8$ and has an estimated virial mass close the maximum mass allowed by standard cosmological models. The JWST images also reveal the presence of small-mass perturbers that produce small lensing distortions. The smallest of these is consistent with being a dwarf galaxy at $z=0.87$ and has an estimated mass of $3.8\times10^9$~\Msol, making it the smallest substructure found at $z>0.5$. The JWST images also show several candidate caustic-crossing events. One of them is detected at high significance at the expected position of the critical curve and is likely a red supergiant star at $z=2.1878$. This would be the first red supergiant found at cosmological distances. The cluster lensing should magnify background objects at $z>6$, making more of them visible than in blank fields of similar size, but there appears to be a deficiency of such objects., Comment: 27 pages, 21 figures

We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST "Prime Extragalactic Areas for Reionization and Lensing Science" ("PEARLS") project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift proto-clusters; and the iconic backlit VV 191 galaxy system to map its dust attenuation. PEARLS also includes NIRISS spectra for one of the NEP fields and NIRSpec spectra of two high-redshift quasars. The main goal of PEARLS is to study the epoch of galaxy assembly, AGN growth, and First Light. Five fields, the JWST NEP Time-Domain Field (TDF), IRAC Dark Field (IDF), and three lensing clusters, will be observed in up to four epochs over a year. The cadence and sensitivity of the imaging data are ideally suited to find faint variable objects such as weak AGN, high-redshift supernovae, and cluster caustic transits. Both NEP fields have sightlines through our Galaxy, providing significant numbers of very faint brown dwarfs whose proper motions can be studied. Observations from the first spoke in the NEP TDF are public. This paper presents our first PEARLS observations, their NIRCam data reduction and analysis, our first object catalogs, the 0.9-4.5 $\mu$m galaxy counts and Integrated Galaxy Light. We assess the JWST sky brightness in 13 NIRCam filters, yielding our first constraints to diffuse light at 0.9-4.5 {\mu}m. PEARLS is designed to be of lasting benefit to the community., Comment: Accepted to AJ, comments welcome. We ask anyone who uses our public PEARLS (NEP TDF) data to refer to this overview paper

Using the first epoch of four-band NIRCam observations obtained by the James Webb Space Telescope (JWST) Prime Extragalactic Areas for Reionization and Lensing Science Program in the Spitzer IRAC Dark Field, we search for F150W and F200W dropouts. In 14.2 arcmin^2, we have found eight F150W dropouts and eight F200W dropouts, all brighter than 27.5 mag (the brightest being ~24 mag) in the band to the red side of the break. As they are detected in multiple bands, these must be real objects. Their nature, however, is unclear, and characterizing their properties is important for realizing the full potential of JWST. If the observed color decrements are due to the Lyman break, these objects should be at z >~ 11.7 and z >~ 15.4, respectively. The color diagnostics show that at least four F150W dropouts are far away from the usual contaminators encountered in dropout searches (red galaxies at much lower redshifts or brown dwarf stars). While the diagnostics of the F200W dropouts are less certain due to the limited number of passbands, at least one of them is likely not a known type of contaminant, and the rest are consistent with either high-redshift galaxies with evolved stellar populations or old galaxies at z ~ 3 to 8. If a significant fraction of our dropouts are indeed at z ~ 12, we have to face the severe problem of explaining their high luminosities and number densities. Spectroscopic identifications of such objects are urgently needed., Comment: Accepted for publication in ApJL

We derive the spatial and wavelength behavior of dust attenuation in the multiple-armed spiral galaxy VV191b using backlighting by the superimposed elliptical system VV191a in a pair with an exceptionally favorable geometry for this measurement. Imaging using JWST and HST spans the wavelength range 0.3-4.5 microns with high angular resolution, tracing the dust in detail from 0.6 to 1.5 microns. Distinct dust lanes continue well beyond the bright spiral arms, and trace a complex web, with a very sharp radial cutoff near 1.7 Petrosian radii. We present attenuation profiles and coverage statistics in each band at radii 14-21 kpc. We derive the attenuation law with wavelength; the data both within and between the dust lanes clearly favor a stronger reddening behavior (R ~ 2.0 between 0.6 and 0.9 microns, approaching unity by 1.5 microns) than found for starbursts and star-forming regions of galaxies. Power-law extinction behavior lambda^(-beta) gives beta=2.1 from 0.6-0.9 microns. R decreases at increasing wavelengths (R~1.1 between 0.9 and 1.5 microns), while beta steepens to 2.5. Mixing regions of different column density flattens the wavelength behavior, so these results suggest a different grain population than in our vicinity. The NIRCam images reveal a lens arc and counterimage from a background galaxy at z~1, spanning 90 degrees azimuthally at 2.8" from the foreground elliptical galaxy nucleus, and an additional weakly-lensed galaxy. The lens model and imaging data give a mass/light ratio 7.6 in solar units within the Einstein radius 2.0 kpc., Comment: Accepted by Astron. J. Analysis redone since submission, using updated JWST calibrations. Dust reddening behavior is steeper with wavelength and lensed galaxy redshift lower than we first derived

We derive single S\'ersic fits and bulge-disk decompositions for 13096 galaxies at redshifts z < 0.08 in the GAMA II equatorial survey regions in the Kilo-Degree Survey (KiDS) g, r and i bands. The surface brightness fitting is performed using the Bayesian two-dimensional profile fitting code ProFit. We fit three models to each galaxy in each band independently with a fully automated Markov-chain Monte Carlo analysis: a single S\'ersic model, a S\'ersic plus exponential and a point source plus exponential. After fitting the galaxies, we perform model selection and flag galaxies for which none of our models are appropriate (mainly mergers/Irregular galaxies). The fit quality is assessed by visual inspections, comparison to previous works, comparison of independent fits of galaxies in the overlap regions between KiDS tiles and bespoke simulations. The latter two are also used for a detailed investigation of systematic error sources. We find that our fit results are robust across various galaxy types and image qualities with minimal biases. Errors given by the MCMC underestimate the true errors typically by factors 2-3. Automated model selection criteria are accurate to > 90 % as calibrated by visual inspection of a subsample of galaxies. We also present g-r component colours and the corresponding colour-magnitude diagram, consistent with previous works despite our increased fit flexibility. Such reliable structural parameters for the components of a diverse sample of galaxies across multiple bands will be integral to various studies of galaxy properties and evolution. All results are integrated into the GAMA database., Comment: 36 pages, 33 figures. Accepted for publication in MNRAS