Galaxy evolution in the William Herschel deep field

In this Thesis we investigate the evolutionary histories of faint field galaxies using extremely deep optical and near-infrared photometry. Our work is centred on a 50 arcmin(^2) region at high galactic latitude which we call "The William Herschel Deep Field" (WHDF). In this work we describe three new near-infrared surveys of this field. In considering both this infrared data and the existing optical data, our broad aims are to increase our understanding of both the growth of galaxy clustering in the Universe and also to determine the star-formation histories of the field galaxy population. We consider our observations primarily in the context of luminosity evolution models in low density universes, but alternative scenarios are considered. Near-infrared galaxy counts derived from our catalogues are consistent with the predictions of our models, without the need for a steep faint-end slope for the galaxy luminosity function. We find that optical-infrared colour distributions of infrared-selected galaxies in the WHDF are deficient in red, early-type galaxies. This is consistent with the predictions of evolutionary models in which these systems have a small amount of on-going star-formation. We measure the amplitude of galaxy clustering in the WHDF for galaxies selected in optical and near-infrared bandpasses using the projected two-point correlation function. By comparing our measured clustering amplitudes with the predictions of our models we find that in all bandpasses the growth of galaxy clustering is approximately fixed in proper co-ordinates, again assuming a low-density Universe. Finally, an analysis of errors on the correlation function measurements suggest that discrepancies between our work and those of other authors may be explained by an underestimation of statistical errors

Extragalactic Fields Optimized for Adaptive Optics

In this paper we present the coordinates of 67 55' x 55' patches of sky which have the rare combination of both high stellar surface density (>0.5 arcmin^{-2} with 13<R<16.5 mag) and low extinction (E(B-V)<0.1). These fields are ideal for adaptive-optics based follow-up of extragalactic targets. One region of sky, situated near Baade's Window, contains most of the patches we have identified. Our optimal field, centered at RA: 7h24m3s, Dec: -1deg27'15", has an additional advantage of being accessible from both hemispheres. We propose a figure of merit for quantifying real-world adaptive optics performance, and use this to analyze the performance of multi-conjugate adaptive optics in these fields. We also compare our results to those that would be obtained in existing deep fields. In some cases adaptive optics observations undertaken in the fields given in this paper would be orders of magnitude more efficient than equivalent observations undertaken in existing deep fields.Comment: 28 pages, 15 figures, 1 table; accepted for publication in PAS

The Next Generation Virgo Cluster Survey - Infrared (NGVS-IR): I. A new Near-UV/Optical/Near-IR Globular Cluster selection tool

The NGVS-IR project (Next Generation Virgo Survey - Infrared) is a contiguous near-infrared imaging survey of the Virgo cluster of galaxies. It complements the optical wide-field survey of Virgo (NGVS). The current state of NGVS-IR consists of Ks-band imaging of 4 deg^2 centered on M87, and J and Ks-band imaging of 16 deg^2 covering the region between M49 and M87. In this paper, we present the observations of the central 4 deg^2 centered on Virgo's core region. The data were acquired with WIRCam on the Canada-France-Hawaii Telescope and the total integration time was 41 hours distributed in 34 contiguous tiles. A survey-specific strategy was designed to account for extended galaxies while still measuring accurate sky brightness within the survey area. The average 5\sigma limiting magnitude is Ks=24.4 AB mag and the 50% completeness limit is Ks=23.75 AB mag for point source detections, when using only images with better than 0.7" seeing (median seeing 0.54"). Star clusters are marginally resolved in these image stacks, and Virgo galaxies with \mu_Ks=24.4 AB mag arcsec^-2 are detected. Combining the Ks data with optical and ultraviolet data, we build the uiK color-color diagram which allows a very clean color-based selection of globular clusters in Virgo. This diagnostic plot will provide reliable globular cluster candidates for spectroscopic follow-up campaigns needed to continue the exploration of Virgo's photometric and kinematic sub-structures, and will help the design of future searches for globular clusters in extragalactic systems. Equipped with this powerful new tool, future NGVS-IR investigations based on the uiK diagram will address the mapping and analysis of extended structures and compact stellar systems in and around Virgo galaxies.Comment: 23 pages, 18 figures. Accepted for publication in ApJ

Deep Spitzer 24 μm COSMOS Imaging. I. The Evolution of Luminous Dusty Galaxies—Confronting the Models

We present the first results obtained from the identification of ~30,000 sources in the Spitzer/24 μm observations of the COSMOS field at S_(24 μm) ≳ 80 μJy. Using accurate photometric redshifts (σ_ z ~ 0.12 at z ~ 2 for 24 μm sources with i ^+ ≳ 25 mag AB) and simple extrapolations of the number counts at faint fluxes, we resolve with unprecedented detail the buildup of the mid-infrared background across cosmic ages. We find that ~50% and ~80% of the 24 μm background intensity originate from galaxies at z ≳ 1 and z ≳ 2, respectively, supporting the scenario where highly obscured sources at very high redshifts (z ≳ 2) contribute only marginally to the cosmic infrared background. Assuming flux-limited selections at optical wavelengths, we also find that the fraction of i ^+-band sources with 24 μm detection strongly increases up to z ~ 2 as a consequence of the rapid evolution that star-forming galaxies have undergone with look-back time. Nonetheless, this rising trend shows a clear break at z ~ 1.3, probably due to k-correction effects implied by the complexity of spectral energy distributions in the mid-infrared. Finally, we compare our results with the predictions from different models of galaxy formation. We note that semianalytical formalisms currently fail to reproduce the redshift distributions observed at 24 μm. Furthermore, the simulated galaxies at S _(24 μm) > 80 μJy exhibit R–K colors much bluer than observed and the predicted K-band fluxes are systematically underestimated at z ≳ 0.5. Unless these discrepancies mainly result from an incorrect treatment of extinction in the models they may reflect an underestimate of the predicted density of high-redshift massive sources with strong ongoing star formation, which would point to more fundamental processes and/or parameters (e.g., initial mass function, critical density to form stars, feedback,...) that are still not fully controlled in the simulations. The most recent backward evolution scenarios reproduce reasonably well the flux/redshift distribution of 24 μm sources up to z ~ 3, although none of them is able to exactly match our results at all redshifts

The fmos-cosmos survey of star-forming galaxies at z ~ 1.6. I. Hα-based star formation rates and dust extinction

We present the first results from a near-IR spectroscopic survey of the COSMOS field, using the Fiber Multi-Object Spectrograph on the Subaru telescope, designed to characterize the star-forming galaxy population at 1.4 < z < 1.7. The high-resoluti

Tracking radiation damage of Euclid VIS detectors after one year in space

Due to the space radiation environment at L2, ESA’s Euclid mission will be subject to a large amount of highly energetic particles over its lifetime. These particles can cause damage to the detectors by creating defects in the silicon lattice. These defects degrade the returned image in several ways, one example being a degradation of the Charge Transfer Efficiency, which appears as readout trails in the image data. This can be problematic for the Euclid VIS instrument, which aims to measure the shapes of galaxies to a very high degree of accuracy. Using a special clocking technique called trap pumping, the single defects in the CCDs can be detected and characterised. Being the first instrument in space with this capability, it will provide novel insights into the creation and evolution of radiation-induced defects and give input to the radiation damage correction of the scientific d ata. We present the status of the radiation damage of the Euclid VIS CCDs and how it has evolved over the first year in space

Tracking radiation damage of Euclid VIS detectors after 1 year in space

Due to the space radiation environment at L2, ESA's Euclid mission will be subject to a large amount of highly energetic particles over its lifetime. These particles can cause damage to the detectors by creating defects in the silicon lattice. These defects degrade the returned image in several ways, one example being a degradation of the Charge Transfer Efficiency, which appears as readout trails in the image data. This can be problematic for the Euclid VIS instrument, which aims to measure the shapes of galaxies to a very high degree of accuracy. Using a special clocking technique called trap pumping, the single defects in the CCDs can be detected and characterised. Being the first instrument in space with this capability, it will provide novel insights into the creation and evolution of radiation-induced defects and give input to the radiation damage correction of the scientific data. We present the status of the radiation damage of the Euclid VIS CCDs and how it has evolved over the first year in space.Comment: Submitted to the Proceedings for the SPIE Space Telescopes and Instrumentation 2024: Optical, Infrared, and Millimeter Wave, paper no. 13092-24, 10 pages, 6 figures

Efficient PSF Modeling with ShOpt.jl: A PSF Benchmarking Study with JWST NIRCam Imaging

With their high angular resolutions of 30--100 mas, large fields of view, and complex optical systems, imagers on next-generation optical/near-infrared space observatories, such as the Near-Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST), present both new opportunities for science and also new challenges for empirical point spread function (PSF) characterization. In this context, we introduce ShOpt, a new PSF fitting tool developed in Julia and designed to bridge the advanced features of PIFF (PSFs in the Full Field of View) with the computational efficiency of PSFEx (PSF Extractor). Along with ShOpt, we propose a suite of non-parametric statistics suitable for evaluating PSF fit quality in space-based imaging. Our study benchmarks ShOpt against the established PSF fitters PSFEx and PIFF using real and simulated COSMOS-Web Survey imaging. We assess their respective PSF model fidelity with our proposed diagnostic statistics and investigate their computational efficiencies, focusing on their processing speed relative to the complexity and size of the PSF models. We find that ShOpt can already achieve PSF model fidelity comparable to PSFEx and PIFF while maintaining competitive processing speeds, constructing PSF models for large NIRCam mosaics within minutes.Comment: 53 pages, 27 figures, submitted to Astronomical Journal. Same as v2 with a corrected typo, which contained revisions from v

JWST and ALMA discern the assembly of structural and obscured components in a high-redshift starburst galaxy

We present observations and analysis of the starburst, PACS-819, at z=1.45 ($M_*=10^{10.7}$ M$_{ \odot}$), using high-resolution ($0^{\prime \prime}.1$; 0.8 kpc) ALMA and multi-wavelength JWST images from the COSMOS-Web program. Dissimilar to HST/ACS images in the rest-frame UV, the redder NIRCam and MIRI images reveal a smooth central mass concentration and spiral-like features, atypical for such an intense starburst. Through dynamical modeling of the CO J=5--4 emission with ALMA, PACS-819 is rotation-dominated thus has a disk-like nature. However, kinematic anomalies in CO and asymmetric features in the bluer JWST bands (e.g., F150W) support a more disturbed nature likely due to interactions. The JWST imaging further enables us to map the distribution of stellar mass and dust attenuation, thus clarifying the relationships between different structural components, not discernable in the previous HST images. The CO J = 5 -- 4 and FIR dust continuum emission are co-spatial with a heavily-obscured starbursting core (<1 kpc) which is partially surrounded by much less obscured star-forming structures including a prominent arc, possibly a tidally-distorted dwarf galaxy, and a clump, either a sign of an ongoing violent disk instability or a recently accreted low-mass satellite. With spatially-resolved maps, we find a high molecular gas fraction in the central area reaching $\sim3$ ($M_{\text{gas}}$/$M_*$) and short depletion times ($M_{\text{gas}}/SFR\sim$ 120 Myrs) across the entire system. These observations provide insights into the complex nature of starbursts in the distant universe and underscore the wealth of complementary information from high-resolution observations with both ALMA and JWST.Comment: 18 pages, 12 figures, Submitted to Ap

COSMOS-Web: The over-abundance and physical nature of "little red dots"--Implications for early galaxy and SMBH assembly

JWST has revealed a population of compact and extremely red galaxies at $z>4$, which likely host active galactic nuclei (AGN). We present a sample of 434 ``little red dots'' (LRDs), selected from the 0.54 deg$^2$ COSMOS-Web survey. We fit galaxy and AGN SED models to derive redshifts and physical properties; the sample spans $z\sim5$-$9$ after removing brown dwarf contaminants. We consider two extreme physical scenarios: either LRDs are all AGN, and their continuum emission is dominated by the accretion disk, or they are all compact star-forming galaxies, and their continuum is dominated by stars. If LRDs are AGN-dominated, our sample exhibits bolometric luminosities $\sim10^{45-47}$ erg\,s$^{-1}$, spanning the gap between JWST AGN in the literature and bright, rare quasars. We derive a bolometric luminosity function (LF) $\sim100$ times the (UV-selected) quasar LF, implying a non-evolving black hole accretion density of $\sim10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ from $z\sim2$-$9$. By contrast, if LRDs are dominated by star formation, we derive stellar masses $\sim10^{8.5-10}\,M_\odot$. MIRI/F770W is key to deriving accurate stellar masses; without it, we derive a mass function inconsistent with $\Lambda$CDM. The median stellar mass profile is broadly consistent with the maximal stellar mass surface densities seen in the nearby universe, though the most massive $\sim50$\% of objects exceed this limit, requiring substantial AGN contribution to the continuum. Nevertheless, stacking all available X-ray, mid-IR, far-IR/sub-mm, and radio data yields non-detections. Whether dominated by dusty AGN, compact star-formation, or both, the high masses/luminosities and remarkable abundance of LRDs implies a dominant mode of early galaxy/SMBH growth.Comment: 27 pages, 13 figures. Submitted to ApJ, comments welcome! Data access at https://github.com/hollisakins/akins24_c