Search

A crucial yet challenging task in galaxy evolution studies is the identification of distant merging galaxies, a task which suffers from a variety of issues ranging from telescope sensitivities and limitations to the inherently chaotic morphologies of young galaxies. In this paper, we use random forests and convolutional neural networks to identify high-redshift JWST CEERS galaxy mergers. We train these algorithms on simulated $3

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 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

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

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 present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifting bandpasses to sample the same rest-frame emission blueward and redward of the Balmer break. This approach allows us to identify LRDs over a wider redshift range and is less susceptible to contamination from galaxies with strong breaks that otherwise lack a rising red continuum. The redshift distribution of our sample increases at $z<8$ and then undergoes a rapid decline at $z\sim4.5$, which may tie the emergence, and obscuration, of these sources to the inside-out growth that galaxies experience during this epoch. We find that LRDs are 2-3 dex more numerous than bright quasars at $z\sim5-7$, but their number density is only 0.6-1 dex higher than X-ray and UV selected AGN at these redshifts. Within our sample, we have identified the first X-ray detected LRDs at $z=3.1$ and $z=4.66$. An X-ray spectral analysis confirms that these AGN are moderately obscured with $\log\,(N_{\rm H}/{\rm cm}^{2}$) of $23.3^{+0.4}_{-1.3}$ and $22.72^{+0.13}_{-0.16}$. Our analysis reveals that reddened AGN emission dominates their rest-optical light, while the rest-UV originates from their host galaxies. We also present NIRSpec follow-up spectroscopy of 17 LRDs that show broad emission lines consistent with AGN activity. The confirmed AGN fraction of our sample is $71\%$ for sources with F444W$<26.5$. In addition, we find three LRDs with narrow blue-shifted Balmer absorption features in their spectra, suggesting an outflow of high-density, low ionization gas from near the central engine of these faint, red AGN., Comment: 23 pages, 17 figures, submitted to ApJ

We utilize deep JWST NIRCam observations for the first direct constraints on the Galaxy Stellar Mass Function (GSMF) at z>10. Our EPOCHS v1 sample includes 1120 galaxy candidates at 6.5 9 than some studies predict, suggesting a rapid buildup of stellar mass in the early Universe. We find no incompatibility between our results and those of standard cosmological models, as suggested previously, although the most massive galaxies may require a high star formation efficiency. We find that the 'Little Red Dot' galaxies dominate the z=7 GSMF at high-masses, necessitating a better understanding of the relative contributions of AGN and stellar emission. We show that assuming a theoretically motivated top-heavy IMF reduces stellar mass by 0.5 dex without affecting fit quality, but our results remain consistent with existing cosmological models with a standard IMF., Comment: 40 pages, 15 figures. Submitted to ApJ. Data is avaible on Github: https://github.com/tHarvey303/EpochsIV. Comments to corresponding author welcome at thomas.harvey-3@manchester.ac.uk

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

Ultraviolet (UV; rest-frame $\sim1200-2000$ A) spectra provide a wealth of diagnostics to characterize fundamental galaxy properties, such as their chemical enrichment, the nature of their stellar populations, and their amount of Lyman-continuum (LyC) radiation. In this work, we leverage publicly released JWST data to construct the rest-frame UV-to-optical composite spectrum of a sample of 63 galaxies at $5.6

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

The wavelength-coverage and sensitivity of JWST now enables us to probe the rest-frame UV - optical spectral energy distributions (SEDs) of galaxies at high-redshift ($z>4$). From these SEDs it is, in principle, through SED fitting possible to infer key physical properties, including stellar masses, star formation rates, and dust attenuation. These in turn can be compared with the predictions of galaxy formation simulations allowing us to validate and refine the incorporated physics. However, the inference of physical properties, particularly from photometry alone, can lead to large uncertainties and potential biases. Instead, it is now possible, and common, for simulations to be \emph{forward-modelled} to yield synthetic observations that can be compared directly to real observations. In this work, we measure the JWST broadband fluxes and colours of a robust sample of $58$ the distributions differ somewhat, though our observed sample size is small and thus susceptible to statistical fluctuations. Likewise, the predicted and observed colour evolution show broad agreement, at least at $58$, though, again, the sample size is small here., Comment: 11 pages, 10 figures, submitted to MNRAS

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

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

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

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

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

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

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

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

We use James Webb Space Telescope Near-Infrared Camera Wide Field Slitless Spectroscopy (NIRCam WFSS) and Near-Infrared spectrograph (NIRSpec) in the Cosmic Evolution Early Release survey (CEERS) to measure rest-frame optical emission-line of 155 galaxies at z>2. The blind NIRCam grism observations include a sample of galaxies with bright emission lines that were not observed on the NIRSpec masks. We study the changes of the Ha, [OIII]/Hb, and [NeIII]/[OII] emission lines in terms of redshift by comparing to lower redshift SDSS and CLEAR samples. We find a significant (>3$\sigma$) correlation between [OIII]/Hb with redshift, while [NeIII]/[OII] has a marginal (2$\sigma$) correlation with redshift. We compare [OIII]/Hb and [NeIII]/[OII] to stellar mass and Hb SFR. We find that both emission-line ratios have a correlation with Hb SFR and an anti-correlation with stellar mass across the redshifts 0

We investigate the morphology and resolved physical properties of a sample of 22 IR-selected DSFG at cosmic noon using the JWST/NIRCam images obtained in the EGS field for the CEERS survey. The resolution of the NIRCam images allowed to spatially resolve these galaxies up to 4.4um and identify their bulge even when extinguished by dust. The goal of this study is to obtain a better understanding of the formation and evolution of FIR-bright galaxies by spatially resolving their properties using JWST in order to look through the dust and bridge the gap between the compact FIR sources and the larger optical SFG. Based on RGB images from the NIRCam filters, we divided each galaxy into several uniformly colored regions, fitted their respective SEDs, and measured physical properties. After classifying each region as SF or quiescent, we assigned galaxies to three classes, depending on whether active SF is located in the core, in the disk or in both. We find (i) that galaxies at a higher z tend to have a fragmented disk with a low core mass fraction. They are at an early stage of bulge formation. When moving toward a lower z, the core mass fraction increases, and the bulge growth is associated with a stabilization of the disk: the NIRCam data clearly point toward bulge formation in preexisting disks. (ii) Lopsidedness is a common feature of DSFGs. It could have a major impact on their evolution; (iii) 23% of galaxies have a SF core embedded in a quiescent disk. They seem to be undergoing outside-in quenching, often facilitated by their strong lopsidedness inducing instabilities. (iv) We show that half of our galaxies with SF concentrated in their core are good SMG counterpart candidates, demonstrating that compact SMGs are usually surrounded by a larger, less obscured disk. (v) Finally, we found surprising evidence for clump-like substructures being quiescent or residing in quiescent regions., Comment: 23 pages, 24 figures, Accepted by A&A

We present a spectroscopic survey of Ly$\alpha$ emitters in the Extended Groth Strip (EGS) field, targeting the regime near the Epoch of Reionization. Using Keck/DEIMOS, we observed 947 high-$z$ candidates with photometric redshifts from 3 $< z_\text{phot} <$ 7 and down to an $H$-band (HST/WFC3 F160W) magnitude limit of < 27.5. Observations were taken over the course of 8 nights, with integration times ranging from 4 to 7.8 hours. Our survey secured 137 unique redshifts, 126 of which are Ly$\alpha$ emitters at 2.8 $< z <$ 6.5 with a mean redshift of $\overline{z} = 4.3$. We provide a comprehensive redshift catalog for our targets, as well as the reduced one- and two- dimensional spectra for each object. These observations will provide an important auxiliary dataset for the JWST Directors Discretionary Early Release Science (DD-ERS) program the Cosmic Evolution Early Release Science Survey (CEERS), which recently completed near- and mid-IR imaging and spectroscopy of galaxies in the EGS field., Comment: 8 pages, 7 figures, 2 tables, Accepted in MNRAS

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 study a new population of extremely red objects (EROs) recently discovered by JWST based on their NIRCam colors F277W$-$F444W $>1.5$ mag. We find 37 EROs in the CEERS field with F444W $<28$ mag and photometric redshifts between $57$ by up to a factor $\sim$60. Similarly, if they are OSOs with luminosities in the L$_{\rm bol}>10^{46-47}$ erg s$^{-1}$ range, their number would exceed that of bright blue QSOs by more than two orders of magnitude. Additional photometry at mid-IR wavelengths will reveal the true nature of the red continuum emission in these EROs and will place this puzzling population in the right context of galaxy evolution., Comment: 26 pages, 10 figures, submitted to ApJ

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{JWST}$/NIRCam+MIRI imaging over the first 20 arcmin$^2$ of publicly-available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W$-$F444W ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED fitting with a variety of codes (including BAGPIPES, PROSPECTOR, BEAGLE, and CIGALE), and find a $>95\%$ probability that these indeed lie at $z>7$. Both sources are compact ($R_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn., Comment: 23 pages, 8 figures, submitted to ApJL, comments welcome

In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. In particular, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on the question of why OFGs and similar galaxies are so faint at optical/near-IR wavelengths. We find that stellar mass is the primary proxy for dust attenuation, among the properties studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation being close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$); OFGs exhibit a 30% smaller effective radius than the average SFG at the same stellar mass and redshift. On the contrary, OFGs do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show that stellar mass is the primary proxy for dust attenuation and compact stellar light profiles behind the thick dust columns obscuring typical massive SFGs., Comment: Accepted by A&A. 13 pages, 9 figures

The large-scale structure (LSS) of the Universe is comprised of galaxy filaments, tendrils, and voids. The majority of the Universe's volume is taken up by these voids, which exist as underdense, but not empty, regions. The galaxies found inside these voids are expected to be some of the most isolated objects in the Universe. This study, using the Galaxy and Mass Assembly (GAMA) and Galaxy Zoo surveys, aims to investigate basic physical properties and morphology of void galaxies versus field (filament and tendril) galaxies. We use void galaxies with stellar masses of $9.35 < log(M/M_\odot) < 11.25$, and this sample is split by identifying two redshift-limited regions, 0 < z < 0.075, and, $0.075 < z < 0.15$. To find comparable objects in the sample of field galaxies from GAMA and Galaxy Zoo, we identify "twins" of void galaxies as field galaxies within $\pm$0.05 dex and $\pm$0.15 dex of M and specific star formation rate. We determine the statistical significance of our results using the Kolmogorov-Smirnov (KS) test. We see that void galaxies, in contrast with field galaxies, seem to be disk-dominated and have predominantly round bulges (with > 50 percent of the Galaxy Zoo citizen scientists agreeing that bulges are present)., Comment: 13 pages, 16 figures, 3 tables, accepted by MNRAS

We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from the Cosmic Evolution Early Release Science Survey (CEERS) NIRCam imaging with photometric redshifts z_phot>8. We measure emission line redshifts of z=7.65 and 8.64 for two galaxies, and z=9.77(+0.37,-0.29) and 10.01(+0.14,-0.19) for two others via the detection of continuum breaks consistent with Lyman-alpha opacity from a mostly neutral intergalactic medium. The presence (absense) of strong breaks (strong emission lines) give high confidence that these two galaxies are at z>9.6, but the break-derived redshifts have large uncertainties given the low spectral resolution and relatively low signal-to-noise of the CEERS NIRSpec prism data. The two z~10 sources are relatively luminous (M_UV<-20), with blue continua (-2.38 candidates with CEERS NIRSpec spectroscopy do not have secure redshifts, but the absence of emission lines in their spectra is consistent with redshifts z>9.6. We find that z>8 photometric redshifts are generally in agreement (within uncertainties) with the spectroscopic values. However, the photometric redshifts tend to be slightly overestimated (average Delta(z)=0.50+/-0.12), suggesting that current templates do not fully describe the spectra of very high-z sources. Overall, our results solidifies photometric evidence for a high space density of bright galaxies at z>8 compared to theoretical model predictions, and further disfavors an accelerated decline in the integrated UV luminosity density at z>8., Comment: Published in ApJL. 19 pages, 9 figures, 7 tables. File with Table 6 included in source .tar file

We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$\alpha$-break by Hubble and with a Ly$\alpha$ redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we observed this source with JWST/NIRSpec spectroscopy, MIRI and NIRCam imaging, and NIRCam/WFSS slitless spectroscopy. The NIRSpec spectra uncover many emission lines, and the strong [O III] emission line confirms the ground-based Ly$\alpha$ redshift. We detect a significant broad (FWHM~1200 km/s) component in the H$\beta$ emission line, which we conclude originates in the broad-line region of an active galactic nucleus (AGN), as the lack of a broad component in the forbidden lines rejects an outflow origin. This hypothesis is supported by the presence of high-ionization lines, as well as a spatial point-source component embedded within a smoother surface brightness profile. The mass of the black hole is log($M_{BH}/M_{\odot})=6.95{\pm}0.37$, and we estimate that it is accreting at 1.2 ($\pm$0.5) x the Eddington limit. The 1-8 $\mu$m photometric spectral energy distribution (SED) from NIRCam and MIRI shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M$_{\odot}$~9.5) and highly star-forming (SFR~30 M$_{\odot}$ yr$^{-1}$). Ratios of the strong emission lines show that the gas in this galaxy is metal-poor (Z/Z$_{\odot}$~0.1), dense (n$_{e}$~10$^{3}$ cm$^{-3}$), and highly ionized (log U~-2.1), consistent with the general galaxy population observed with JWST at high redshifts. We use this presently highest-redshift AGN discovery to place constraints on black hole seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from massive black hole seeds is required to form this object by the observed epoch., Comment: 26 pages, 14 figures, 6 tables, published in ApJL

We explore the properties of galaxies in the proximity (within a $\sim$2 Mpc radius sphere) of Type I quasars at 0.1

Galaxy morphology in atomic hydrogen (HI) and in the ultra-violet (UV) are closely linked. This has motivated their combined use to quantify morphology over the full H i disk for both H i and UV imaging. We apply galaxy morphometrics: Concentration, Asymmetry, Gini, M20 and Multimode-Intensity-Deviation statistics to the first moment-0 maps of the WALLABY survey of galaxies in the Hydra cluster center. Taking advantage of this new HI survey, we apply the same morphometrics over the full HI extent on archival GALEX FUV and NUV data to explore how well HI truncated, extended ultraviolet disk (XUV) and other morphological phenomena can be captured using pipeline WALLABY data products. Extended HI and UV disks can be identified relatively straightforward from their respective concentration. Combined with WALLABY HI, even the shallowest GALEX data is sufficient to identify XUV disks. Our second goal is to isolate galaxies undergoing ram-pressure stripping in the H i morphometric space. We employ four different machine learning techniques, a decision tree, a k-nearest neighbour, a support-vector machine, and a random forest. Up to 80% precision and recall are possible with the Random Forest giving the most robust results., Comment: 17 figures, 12 figures, 7 tables, accepted by MNRAS

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

We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4

We perform a ground-based near-infrared spectroscopic survey using the Keck/MOSFIRE spectrograph to target Ly$\alpha$ emission at $7.0$4$\sigma$ level in eight $z>7$ galaxies, which include additional members of the known $z\sim7.7$ Ly$\alpha$-emitter (LAE) cluster (Tilvi et al. 2020). With the addition of these newly-discovered $z\sim7.7$ LAEs, this is currently the largest measured LAE cluster at $z>7$. The unusually-high Ly$\alpha$ detection rate at $z\sim7.7$ in this field suggests significantly stronger Ly$\alpha$ emission from the clustered LAEs than from the rest of our targets. We estimate the ionized bubble sizes around these LAEs and conclude that the LAEs are clustered within an extended ionized structure created by overlapping ionized bubbles which allow the easier escape of Ly$\alpha$ from galaxies. It is remarkable that the brightest object in the cluster has the lowest measured redshift of the Ly$\alpha$ line, being placed in front of the other LAEs in the line-of-sight direction. This suggests that we are witnessing the enhanced IGM transmission of Ly$\alpha$ from galaxies on the rear side of an ionized area. This could be a consequence of Ly$\alpha$ radiative transfer: Ly$\alpha$ close to the central velocity is substantially scattered by the IGM while Ly$\alpha$ from the rear-side galaxies is significantly redshifted to where it has a clear path., Comment: Submitted to ApJ (17 pages, 8 figures, 4 tables)

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 examine the HI gas kinematics of galaxy pairs in two clusters and a group using Australian Square Kilometre Array Pathfinder (ASKAP) WALLABY pilot survey observations. We compare the HI properties of galaxy pair candidates in the Hydra I and Norma clusters, and the NGC 4636 group, with those of non-paired control galaxies selected in the same fields. We perform HI profile decomposition of the sample galaxies using a tool, {\sc baygaud} which allows us to de-blend a line-of-sight velocity profile with an optimal number of Gaussian components. We construct HI super-profiles of the sample galaxies via stacking of their line profiles after aligning the central velocities. We fit a double Gaussian model to the super-profiles and classify them as kinematically narrow and broad components with respect to their velocity dispersions. Additionally, we investigate the gravitational instability of HI gas disks of the sample galaxies using Toomre Q parameters and HI morphological disturbances. We investigate the effect of the cluster environment on the HI properties of galaxy pairs by dividing the cluster environment into three subcluster regions (i.e., outskirts, infalling and central regions). We find that the denser cluster environment (i.e., infalling and central regions) is likely to impact the HI gas properties of galaxies in a way of decreasing the amplitude of the kinematically narrow HI gas ($M_{\rm{narrow}}^{\rm{HI}}$/$M_{\rm{total}}^{\rm{HI}}$), and increasing the Toomre Q values of the infalling and central galaxies. This tendency is likely to be more enhanced for galaxy pairs in the cluster environment., Comment: Accepted for publication in MNRAS, 25 pages, 20 figures. Supplementary material included (20 pages)

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

The selection of high-redshift galaxies often involves spectral energy distribution (SED) fitting to photometric data, an expectation for contamination levels, and measurement of sample completeness -- all vetted through comparison to spectroscopic redshift measurements of a sub-sample. The first JWST data is now being taken over several extragalactic fields, to different depths and across various areas, which will be ideal for the discovery and classification of galaxies out to distances previously uncharted. As spectroscopic redshift measurements for sources in this epoch will not be initially available to compare with the first photometric measurements of z > 8 galaxies, robust photometric redshifts are of the utmost importance. Galaxies at z > 8 are expected to have bluer rest-frame ultraviolet (UV) colors than typically-used model SED templates, which could lead to catastrophic photometric redshift failures. We use a combination of BPASS and Cloudy models to create a supporting set of templates that match the predicted rest-UV colors of z > 8 the simulated galaxies in a mock catalog (Yung et al. 2022), which mimics expected field depths and areas of the Cosmic Evolution Early Release Science Survey (CEERS: m$_{5\sigma}$ ~ 28.6 over ~100 arcmin$^2$; Finkelstein et al. 2022a, Bagley et al. 2022). We use EAZY to highlight the improvements in redshift recovery with the inclusion of our new template set and suggest criteria for selecting galaxies at 8 < z < 10 with JWST, providing an important test case for observers venturing into this new era of astronomy., Comment: 14 pages, 10 figures, submitted to ApJ. SED templates published at: ceers.github.io/LarsonSEDTemplates

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions., Comment: Replaced with published version

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

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

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

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

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

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

Identifying merging galaxies is an important - but difficult - step in galaxy evolution studies. We present random forest classifications of galaxy mergers from simulated JWST images based on various standard morphological parameters. We describe (a) constructing the simulated images from IllustrisTNG and the Santa Cruz SAM, and modifying them to mimic future CEERS observations as well as nearly noiseless observations, (b) measuring morphological parameters from these images, and (c) constructing and training the random forests using the merger history information for the simulated galaxies available from IllustrisTNG. The random forests correctly classify $\sim60\%$ of non-merging and merging galaxies across $0.5 < z < 4.0$. Rest-frame asymmetry parameters appear more important for lower redshift merger classifications, while rest-frame bulge and clump parameters appear more important for higher redshift classifications. Adjusting the classification probability threshold does not improve the performance of the forests. Finally, the shape and slope of the resulting merger fraction and merger rate derived from the random forest classifications match with theoretical Illustris predictions, but are underestimated by a factor of $\sim 0.5$., Comment: 21 pages, 14 figures, submitted to ApJ