Showing 1–50 of 65 results for author: Vijayan, A

First Light And Reionisation Epoch Simulations (FLARES) XVI: Size Evolution of Massive Dusty Galaxies at Cosmic Dawn from UV to IR

Authors: Paurush Punyasheel, Aswin P. Vijayan, Thomas R. Greve, William J. Roper, Hiddo Algera, Steven Gillman, Bitten Gullberg, Dimitrios Irodotou, Christopher C. Lovell, Louise T. C. Seeyave, Peter A. Thomas, Stephen M. Wilkins

Abstract: We use the First Light And Reionisation Epoch Simulations (FLARES) to study the evolution of the rest-frame ultraviolet (UV) and far-infrared (FIR) sizes for a statistical sample of massive ($\gtrsim10^{9}$M$_{\odot}$) high redshift galaxies (z $\in$ [5,10]). Galaxies are post-processed using the SKIRT radiative transfer code, to self-consistently obtain the full spectral energy distribution and s… ▽ More We use the First Light And Reionisation Epoch Simulations (FLARES) to study the evolution of the rest-frame ultraviolet (UV) and far-infrared (FIR) sizes for a statistical sample of massive ($\gtrsim10^{9}$M$_{\odot}$) high redshift galaxies (z $\in$ [5,10]). Galaxies are post-processed using the SKIRT radiative transfer code, to self-consistently obtain the full spectral energy distribution and surface brightness distribution. We create mock observations of the galaxies for the Near Infrared Camera (NIRCam) to study the rest-frame UV 1500 $\unicode{xC5}$ morphology. We also generate mock rest-frame FIR (50 $μ$m) photometry and mock ALMA (158 $μ$m) (0.01"-0.03" and $\approx$0.3" angular resolution) observations to study the dust-continuum. We find the effect of dust on observed sizes reduces with increasing wavelength from the UV to optical ($\sim$0.6 times the UV at 0.4$μ$m), with no evolution in FIR sizes. Observed sizes vary within 0.4-1.2 times the intrinsic sizes at different signal to noise ratios (SNR = 5-20) across redshifts. The effect of PSF and noise makes bright structures prominent, whereas fainter regions blend with noise, leading to an underestimation (factor of 0.4-0.8) of sizes at SNR=5. At SNR=15-20, the underestimation reduces (factor of 0.6-0.9) at z=5-8 but due to PSF, at z=9-10, bright cores are dominant, resulting in an overestimation (factor of 1.0-1.2). For ALMA, low resolution sizes are effected by noise which acts as extended emission. The size evolution in UV broadly agrees with current observational samples and other simulations. This work is one of the first to analyse the panchromatic sizes of a statistically significant sample of simulated high-redshift galaxies, complementing a growing body of research highlighting the importance of conducting an equivalent comparison between observed galaxies and their simulated counterparts in the early Universe. △ Less

Submitted 20 August, 2024; originally announced August 2024.

Comments: 22 pages, 28 figures, data along with code to reproduce figures will be available on GitHub after publication

SCUBADive I: JWST+ALMA Analysis of 289 sub-millimeter galaxies in COSMOS-Web

Authors: Jed McKinney, Caitlin M. Casey, Arianna S. Long, Olivia R. Cooper, Sinclaire M. Manning, Maximilien Franco, Hollis Akin, Erini Lambrides, Elaine Gammon, Camila Silva, Fabrizio Gentile, Jorge A. Zavala, Aristeidis Amvrosiadis, Irma Andika, Malte Brinch, Jaclyn B. Champagne, Nima Chartab, Nicole E. Drakos, Andreas L. Faisst, Seiji Fujimoto, Steven Gillman, Ghassem Gozaliasl, Thomas R. Greve, Santosh Harish, Christopher C. Hayward , et al. (14 additional authors not shown)

Abstract: JWST has enabled detecting and spatially resolving the heavily dust-attenuated stellar populations of sub-millimeter galaxies, revealing detail that was previously inaccessible. In this work we construct a sample of 289 sub-millimeter galaxies with detailed joint ALMA and JWST constraints in the COSMOS field. Sources are originally selected using the SCUBA-2 instrument and have archival ALMA obser… ▽ More JWST has enabled detecting and spatially resolving the heavily dust-attenuated stellar populations of sub-millimeter galaxies, revealing detail that was previously inaccessible. In this work we construct a sample of 289 sub-millimeter galaxies with detailed joint ALMA and JWST constraints in the COSMOS field. Sources are originally selected using the SCUBA-2 instrument and have archival ALMA observations from various programs. Their JWST NIRCam imaging is from COSMOS-Web and PRIMER. We extract multi-wavelength photometry in a manner that leverages the unprecedented near-infrared spatial resolution of JWST, and fit the data with spectral energy distribution models to derive photometric redshifts, stellar masses, star-formation rates and optical attenuation. The sample has an average z=2.6, A_V=2.5, SFR=270 and log(M*)=11.1. There are 81 (30%) galaxies that have no previous optical/near-infrared detections, including 75% of the z>4 sub-sample (n=28). The faintest observed near-infrared sources have the highest redshifts and largest A_V=4. In a preliminary morphology analysis we find that ~10% of our sample exhibit spiral arms and 5% host stellar bars, with one candidate bar found at z>3. Finally, we find that the clustering of JWST galaxies within 10 arcseconds of a sub-mm galaxy is a factor of 2 greater than what is expected based on either random clustering or the distribution of sources around any red galaxy irrespective of a sub-mm detection. △ Less

Submitted 15 August, 2024; originally announced August 2024.

Comments: 37 pages (15 for RGBs + references), 14 figures, submitted to ApJ

Accurate Simultaneous Constraints on the Dust Mass, Temperature and Emissivity Index of a Galaxy at Redshift 7.31

Authors: Hiddo Algera, Hanae Inami, Ilse De Looze, Andrea Ferrara, Hiroyuki Hirashita, Manuel Aravena, Tom Bakx, Rychard Bouwens, Rebecca Bowler, Elisabete Da Cunha, Pratika Dayal, Yoshinobu Fudamoto, Jacqueline Hodge, Alexander Hygate, Ivana van Leeuwen, Themiya Nanayakkara, Marco Palla, Andrea Pallottini, Lucie Rowland, Renske Smit, Laura Sommovigo, Mauro Stefanon, Aswin Vijayan, Paul van der Werf

Abstract: We present new multi-frequency ALMA continuum observations of the massive [$\log_{10}(M_\star/M_\odot) = 10.3_{-0.2}^{+0.1}$], UV-luminous [$M_\mathrm{UV} = -21.7 \pm 0.2$] $z=7.31$ galaxy REBELS-25 in Bands 3, 4, 5, and 9. Combining the new observations with previously-taken data in Bands 6 and 8, we cover the dust continuum emission of the galaxy in six distinct bands -- spanning rest-frame… ▽ More We present new multi-frequency ALMA continuum observations of the massive [$\log_{10}(M_\star/M_\odot) = 10.3_{-0.2}^{+0.1}$], UV-luminous [$M_\mathrm{UV} = -21.7 \pm 0.2$] $z=7.31$ galaxy REBELS-25 in Bands 3, 4, 5, and 9. Combining the new observations with previously-taken data in Bands 6 and 8, we cover the dust continuum emission of the galaxy in six distinct bands -- spanning rest-frame $50-350\,μ$m -- enabling simultaneous constraints on its dust mass ($M_\mathrm{dust}$), temperature ($T_\mathrm{dust}$) and emissivity index ($β_\mathrm{IR}$) via modified blackbody fitting. Given a fiducial model of optically thin emission, we infer a cold dust temperature of $T_\mathrm{dust} = 32_{-6}^{+9}\,$K and a high dust mass of $\log_{10}(M_\mathrm{dust}/M_\odot) = 8.2_{-0.4}^{+0.6}$, and moderately optically thick dust does not significantly alter these estimates. If we assume dust production is solely through supernovae (SNe), the inferred dust yield would be high, $y = 0.7_{-0.4}^{+2.3}\,M_\odot$ per SN. Consequently, we argue grain growth in the interstellar medium of REBELS-25 also contributes to its dust build-up. This is supported by the steep dust emissivity index $β_\mathrm{IR} = 2.5 \pm 0.4$ we measure for REBELS-25, as well as by its high stellar mass, dense interstellar medium, and metal-rich nature. Our results suggest that constraining the dust emissivity indices of high-redshift galaxies is important not only to mitigate systematic uncertainties in their dust masses and obscured star formation rates, but also to assess if dust properties evolve across cosmic time. We present an efficient observing setup to do so with ALMA, combining observations of the peak and Rayleigh-Jeans tail of the dust emission. △ Less

Submitted 5 August, 2024; originally announced August 2024.

Comments: 14 pages main text, 6 figures + appendices; re-submitted to MNRAS after a positive referee report

The RAdio Galaxy Environment Reference Survey (RAGERS): Evidence of an anisotropic distribution of submillimeter galaxies in the 4C 23.56 protocluster at z=2.48

Authors: Dazhi Zhou, Thomas R. Greve, Bitten Gullberg, Minju M. Lee, Luca Di Mascolo, Simon R. Dicker, Charles E. Romero, Scott C. Chapman, Chian-Chou Chen, Thomas Cornish, Mark J. Devlin, Luis C. Ho, Kotaro Kohno, Claudia D. P. Lagos, Brian S. Mason, Tony Mroczkowski, Jeff F. W. Wagg, Q. Daniel Wang, Ran Wang, Malte. Brinch, Helmut Dannerbauer, Xue-Jian Jiang, Lynge R. B. Lauritsen, Aswin P. Vijayan, David Vizgan , et al. (19 additional authors not shown)

Abstract: High-redshift radio(-loud) galaxies (H$z$RGs) are massive galaxies with powerful radio-loud active galactic nuclei (AGNs) and serve as beacons for protocluster identification. However, the interplay between H$z$RGs and the large-scale environment remains unclear. To understand the connection between H$z$RGs and the surrounding obscured star formation, we investigated the overdensity and spatial di… ▽ More High-redshift radio(-loud) galaxies (H$z$RGs) are massive galaxies with powerful radio-loud active galactic nuclei (AGNs) and serve as beacons for protocluster identification. However, the interplay between H$z$RGs and the large-scale environment remains unclear. To understand the connection between H$z$RGs and the surrounding obscured star formation, we investigated the overdensity and spatial distribution of submillimeter-bright galaxies (SMGs) in the field of 4C\,23.56, a well-known H$z$RG at $z=2.48$. We used SCUBA-2 data ($σ\,{\sim}\,0.6$\,mJy) to estimate the $850\,{\rm μm}$ source number counts and examine the radial and azimuthal overdensities of the $850\,{\rm μm}$ sources in the vicinity of the H$z$RG. The angular distribution of SMGs is inhomogeneous around the H$z$RG 4C\,23.56, with fewer sources oriented along the radio jet. We also find a significant overdensity of bright SMGs (${\rm S}_{850\rm\,μm}\geq5\,$mJy). Faint and bright SMGs exhibit different spatial distributions. The former are concentrated in the core region, while the latter prefer the outskirts of the H$z$RG field. High-resolution observations show that the seven brightest SMGs in our sample are intrinsically bright, suggesting that the overdensity of bright SMGs is less likely due to the source multiplicity. △ Less

Submitted 4 August, 2024; originally announced August 2024.

Comments: 19 pages, 17 figures, 5 tables, accepted to A&A

The Radio Galaxy Environment Reference Survey (RAGERS): a submillimetre study of the environments of massive radio-quiet galaxies at $z = 1{\rm -}3$

Authors: Thomas M. Cornish, Julie L. Wardlow, Thomas R. Greve, Scott Chapman, Chian-Chou Chen, Helmut Dannerbauer, Tomotsugu Goto, Bitten Gullberg, Luis C. Ho, Xue-Jian Jiang, Claudia Lagos, Minju Lee, Stephen Serjeant, Hyunjin Shim, Daniel J. B. Smith, Aswin Vijayan, Jeff Wagg, Dazhi Zhou

Abstract: Measuring the environments of massive galaxies at high redshift is crucial to understanding galaxy evolution and the conditions that gave rise to the distribution of matter we see in the Universe today. While high-$z$ radio galaxies (H$z$RGs) and quasars tend to reside in protocluster-like systems, the environments of their radio-quiet counterparts are relatively unexplored, particularly in the su… ▽ More Measuring the environments of massive galaxies at high redshift is crucial to understanding galaxy evolution and the conditions that gave rise to the distribution of matter we see in the Universe today. While high-$z$ radio galaxies (H$z$RGs) and quasars tend to reside in protocluster-like systems, the environments of their radio-quiet counterparts are relatively unexplored, particularly in the submillimetre, which traces dust-obscured star formation. In this study we search for 850 $μ$m-selected submillimetre galaxies in the environments of massive ($M_{\star} > 10^{11} M_{\odot}$), radio-quiet ($L_{500 {\rm MHz}} \lesssim 10^{25}$ W Hz$^{-1}$) galaxies at $z \sim 1\text{--}3$ using S2COSMOS data. By constructing number counts in circular regions of radius 1--6 arcmin and comparing with blank-field measurements, we find no significant overdensities of SMGs around massive radio-quiet galaxies at any of these scales, despite being sensitive down to overdensities of $δ\sim 0.4$. To probe deeper than the catalogue we also examine the distribution of peaks in the SCUBA-2 SNR map, which reveals only tentative signs of any difference in the SMG densities of the radio-quiet galaxy environments compared to the blank field, and only on smaller scales (1$^{\prime}$ radii, corresponding to $\sim0.5$ Mpc) and higher SNR thresholds. We conclude that massive, radio-quiet galaxies at cosmic noon are typically in environments with $δ\lesssim0.4$, which are either consistent with the blank field or contain only weak overdensities spanning sub-Mpc scales. The contrast between our results and studies of H$z$RGs with similar stellar masses and redshifts implies an intrinsic link between the wide-field environment and radio AGN luminosity at high redshift. △ Less

Submitted 30 August, 2024; v1 submitted 30 July, 2024; originally announced July 2024.

Comments: 13 pages, 5 figures. Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Vol. 533, Issue 1 (2024) pp. 1032-1044

A photo-$z$ cautionary tale: Redshift confirmation of COSBO-7 at $z=2.625$

Authors: Shuowen Jin, Nikolaj B. Sillassen, Jacqueline Hodge, Georgios E. Magdis, Caitlin Casey, Francesca Rizzo, Anton M. Koekemoer, Francesco Valentino, Vasily Kokorev, Benjamin Magnelli, Raphael Gobat, Steven Gillman, Maximilien Franco, Andreas Faisst, Jeyhan Kartaltepe, Eva Schinnerer, Sune Toft, Hiddo S. B. Algera, Santosh Harish, Minju Lee, Daizhong Liu, Marko Shuntov, Margherita Talia, Aswin Vijayan

Abstract: Photometric redshifts are widely used in studies of dusty star-forming galaxies (DSFGs), but catastrophic photo-$z$ failure can undermine all redshift-dependent results. Here we report the spectroscopic redshift confirmation of COSBO-7, a strongly lensed DSFG in the COSMOS-PRIMER field. Recently, using 10 bands of JWST NIRCam and MIRI imaging data on COSBO-7, Ling et al. (2024) reported a photomet… ▽ More Photometric redshifts are widely used in studies of dusty star-forming galaxies (DSFGs), but catastrophic photo-$z$ failure can undermine all redshift-dependent results. Here we report the spectroscopic redshift confirmation of COSBO-7, a strongly lensed DSFG in the COSMOS-PRIMER field. Recently, using 10 bands of JWST NIRCam and MIRI imaging data on COSBO-7, Ling et al. (2024) reported a photometric redshift solution of $z\gtrsim7.0$, favored by four independent spectral energy distribution (SED) fitting codes, and therefore providing an appealing candidate of the most distant massive DSFG. This photo-$z$ solution was also supported by a single line detection in ALMA Band 3 consistent with CO(7-6) at $z=7.46$. However, our new ALMA observations robustly detect two lines in Band 6 identified as CO(7-6) and [CI](2-1) at $z_{\rm spec}=2.625$, and thus the Band 3 line as CO(3-2). The three robust line detections decidedly place COSBO-7 at $z=2.625$, refuting the photo-$z$ solution. We derive physical parameters by fitting NIR-to-mm photometry and lens modeling, revealing that COSBO-7 is a main-sequence galaxy. We examine possible reasons for this photo-$z$ failure and attribute it to (1) the likely underestimation of photometric uncertainty at 0.9$μ$m, and (2) the lack of photometry at wavelengths beyond 20$μ$m. Notably, we recover a bona-fide $z_{\rm phot}\sim 2.3$ by including the existing MIPS $24μ$m photometry, demonstrating the critical importance of mid-infrared data for photo-$z$ robustness. This work highlights a common challenge in modeling SEDs of DSFGs, cautioning against the reliability of photometric redshifts, as well as pseudo-spectroscopic redshifts based on single line detection. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Submitted to A&A Letter

The Structure of Massive Star-Forming Galaxies from JWST and ALMA: Dusty, High Redshift Disk Galaxies

Authors: Steven Gillman, Ian Smail, Bitten Gullberg, A. M. Swinbank, Aswin P. Vijayan, Minju Lee, Gabe Brammer, U. Dudzevičiūtė, Thomas R. Greve, Omar Almaini, Malte Brinch, Scott C. Chapman, Chian-Chou Chen, Soh Ikarashi, Yuichi Matsuda, Wei-Hao Wang, Fabian Walter, Paul P. van der Werf

Abstract: We present an analysis of the JWST NIRCam and MIRI morphological properties of 80 massive ($\log_{10}(M_\ast[M_{\odot}])$=11.2$\pm$0.1) dusty star-forming galaxies at $z$$=$2.7$^{+1.2}_{-0.7}$, identified as sub-millimetre galaxies (SMGs) by ALMA, that have been observed as part of the JWST PRIMER project. To compare the structure of these massive, active galaxies to more typical star-forming gala… ▽ More We present an analysis of the JWST NIRCam and MIRI morphological properties of 80 massive ($\log_{10}(M_\ast[M_{\odot}])$=11.2$\pm$0.1) dusty star-forming galaxies at $z$$=$2.7$^{+1.2}_{-0.7}$, identified as sub-millimetre galaxies (SMGs) by ALMA, that have been observed as part of the JWST PRIMER project. To compare the structure of these massive, active galaxies to more typical star-forming galaxies, we define a sample of 850 field galaxies with matched redshifts and specific star formation rates. We identify 20$\pm$5% of the SMGs as candidate late-stage major mergers, a further 40$\pm$10% as potential minor mergers and 40$\pm$10% which have comparatively undisturbed disk-like morphologies, with no obvious massive neighbours on $\lesssim$20-30kpc (projected) scales. These rates are comparable to those for the field sample and indicate that the majority of the sub-millimetre-detected galaxies are not late-stage major mergers, but have interaction rates similar to the less-active population at $z$$\sim$2-3. Through a multi-wavelength morphological analysis, we establish that SMGs have comparable near-infrared sizes to the less active population, but exhibit lower Sérsic indices, consistent with bulge-less disks and have more structured morphologies at 2$μ$m relative to 4$μ$m. We find evidence for dust reddening as the origin of the morphological differences between the populations, identifying a strong correlation between the F200W$-$F444W pixel colour and the 870$μ$m surface brightness. We conclude that SMGs and less active galaxies at the same epochs share a common disk-like structure, but the weaker bulge components of the SMGs results in a lower dynamical stability. Consequently, instabilities triggered either secularly or by minor external perturbations result in higher levels of activity (and dust content) in SMGs compared to typical star-forming galaxies. [Abridged] △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 29 pages, 13 figures including appendices. Submitted to A&A. Comments welcome

COSMOS-Web: The Role of Galaxy Interactions and Disk Instabilities in Producing Starbursts at z<4

Authors: A. L. Faisst, M. Brinch, C. M. Casey, N. Chartab, M. Dessauges-Zavadsky, N. E. Drakos, S. Gillman, G. Gonzaliasl, C. C. Hayward, O. Ilbert, P. Jablonka, J. S. Kartaltepe, A. M. Koekemoer, V. Kokorev, E. Lambrides, D. Liu, C. Maraston, C. L. Martin, A. Renzini, B. E. Robertson, D. B. Sanders, Z. Sattari, N. Scoville, C. M. Urry, A. P. Vijayan , et al. (27 additional authors not shown)

Abstract: We study of the role of galaxy-galaxy interactions and disk instabilities in producing starburst activity in galaxies out to z=4. For this, we use a sample of 387 galaxies with robust total star formation rate measurements from Herschel, gas masses from ALMA, stellar masses and redshifts from multi-band photometry, and JWST/NIRCam rest-frame optical imaging. Using mass-controlled samples, we find… ▽ More We study of the role of galaxy-galaxy interactions and disk instabilities in producing starburst activity in galaxies out to z=4. For this, we use a sample of 387 galaxies with robust total star formation rate measurements from Herschel, gas masses from ALMA, stellar masses and redshifts from multi-band photometry, and JWST/NIRCam rest-frame optical imaging. Using mass-controlled samples, we find an increased fraction of interacting galaxies in the starburst regime at all redshifts out to z=4. This increase correlates with star formation efficiency (SFE), but not with gas fraction. However, the correlation is weak (and only significant out to z=2), which could be explained by the short duration of SFE increase during interaction. In addition, we find that isolated disk galaxies make up a significant fraction of the starburst population. The fraction of such galaxies with star-forming clumps ("clumpy disks") is significantly increased compared to the main-sequence disk population. Furthermore, this fraction directly correlates with SFE. This is direct observational evidence for a long-term increase of SFE maintained due to disk instabilities, contributing to the majority of starburst galaxies in our sample and hence to substantial mass growth in these systems. This result could also be of importance for explaining the growth of the most massive galaxies at z>6. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 14 pages, 5 figures. Submitted to ApJ. Comments welcome

When are galactic winds molecular?

Authors: Aditi Vijayan, Mark R. Krumholz

Abstract: The molecular phase of supernova-driven outflows originates from the cold, molecular gas in the disc of a star-forming galaxy, and may carry a substantial fraction of the wind mass flux in some galaxies, but it remains poorly understood. Observations of this phase come mostly from very nearby galaxies due its low surface brightness and covering fraction, and simulations often lack the spatial reso… ▽ More The molecular phase of supernova-driven outflows originates from the cold, molecular gas in the disc of a star-forming galaxy, and may carry a substantial fraction of the wind mass flux in some galaxies, but it remains poorly understood. Observations of this phase come mostly from very nearby galaxies due its low surface brightness and covering fraction, and simulations often lack the spatial resolution necessary to resolve it. Here we analytically estimate the survivability of this phase in order to understand under what conditions an galactic wind can contain a significant molecular phase. We show that the molecular content of outflows is primarily determined by two dimensionless numbers: a generalised Eddington ratio describing the strength of the outflow and an ionisation parameter-like quantity describing the strength of the radiation field per baryon. We apply this model to a sample of galaxies and show that, while any molecules entrained in the winds of normal star-forming galaxies should be destroyed close to the galactic disc, the outflows of strong starburst should become increasingly dominated by molecules. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 8 pages, 5 figures, submitted to MNRAS. Comments welcome

Like a candle in the wind: The embers of once aflame, now smouldering galaxies at $5 < z < 8$

Authors: James Trussler, Christopher Conselice, Nathan Adams, Duncan Austin, Joseph Caruana, Tom Harvey, Qiong Li, Christopher Lovell, Louise Seeyave, Aswin Vijayan, Stephen Wilkins

Abstract: We develop a photometric search method for identifying smouldering galaxies at $5< z < 8$, which are defined to have weak emission lines and thus generally have low specific star formation rates and may even be in a state of (temporary) quiescence. The deep NIRCam imaging (${\sim}29.5$ AB mag, 5$σ$) from the JADES second data release is essential for finding these systems, as they are faint, relat… ▽ More We develop a photometric search method for identifying smouldering galaxies at $5< z < 8$, which are defined to have weak emission lines and thus generally have low specific star formation rates and may even be in a state of (temporary) quiescence. The deep NIRCam imaging (${\sim}29.5$ AB mag, 5$σ$) from the JADES second data release is essential for finding these systems, as they are faint, relatively quiescent dwarf galaxies ($M_* \sim 10^{8}$-$10^9$ $\mathrm{M}_\odot)$ in the Epoch of Reionisation (EoR). Moreover, medium-band imaging is key, enabling a clear identification of the lack of emission lines in these galaxies, thus betraying their dormant flame. Owing to the young age of the Universe, combined with the likely bursty star formation in these first dwarf galaxies, conventional colour-selection methods like the UVJ diagram likely miss a large fraction of the quiescent population in the EoR. Indeed, we find that smouldering galaxies constitute a considerable fraction (0.10-0.35) of the EoR dwarf galaxy population ($M_* \sim 10^{8}$-$10^{9}$ $\mathrm{M}_\odot$). As predicted by simulations, these first dwarf galaxies are fragile, the star formation in their shallow potential wells easily snuffed out by feedback-driven winds triggered by secular or merger-driven starbursts, with the smouldering fraction increasing with decreasing stellar mass. Finally, we provide observational constraints on the smouldering galaxy comoving number density (${\sim}10^{-4}$-$10^{-5}$ dex$^{-1}$ Mpc$^{-3}$), which, although hampered by incompleteness, should aid in our understanding of the primordial baryon cycle, as current simulations greatly disagree on whether these systems are rare (${\sim}1\%$) or common (${\sim}50\%$) in the EoR. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 22 pages, 17 figures, 2 tables

First Light and Reionization Epoch Simulations (FLARES) -- XV: The physical properties of super-massive black holes and their impact on galaxies in the early universe

Authors: Stephen M. Wilkins, Jussi K. Kuusisto, Dimitrios Irodotou, Shihong Liao, Christopher C. Lovell, Sonja Soininen, Sabrina C. Berger, Sophie L. Newman, William J. Roper, Louise T. C. Seeyave, Peter A. Thomas, Aswin P. Vijayan

Abstract: Understanding the co-evolution of super-massive black holes (SMBHs) and their host galaxies remains a key challenge of extragalactic astrophysics, particularly the earliest stages at high-redshift. However, studying SMBHs at high-redshift with cosmological simulations, is challenging due to the large volumes and high-resolution required. Through its innovative simulation strategy, the First Light… ▽ More Understanding the co-evolution of super-massive black holes (SMBHs) and their host galaxies remains a key challenge of extragalactic astrophysics, particularly the earliest stages at high-redshift. However, studying SMBHs at high-redshift with cosmological simulations, is challenging due to the large volumes and high-resolution required. Through its innovative simulation strategy, the First Light And Reionisation Epoch Simulations (FLARES) suite of cosmological hydrodynamical zoom simulations allows us to simulate a much wider range of environments which contain SMBHs with masses extending to $M_{\bullet}>10^{9}\ M_{\odot}$ at $z=5$. In this paper, we use FLARES to study the physical properties of SMBHs and their hosts in the early Universe ($5\le\, z \le10$). FLARES predicts a sharply declining density with increasing redshift, decreasing by a factor of 100 over the range $z=5\to 10$. Comparison between our predicted bolometric luminosity function and pre-\emph{JWST} observations yield a good match. However, recent \emph{JWST} observations appear to suggest a larger contribution of SMBHs than previously observed, or predicted by FLARES. Finally, by using a re-simulation with AGN feedback disabled, we explore the impact of AGN feedback on their host galaxies. This reveals that AGN feedback results in a reduction of star formation activity, even at $z>5$, but only in the most massive galaxies. A deeper analysis reveals that AGN are also the cause of suppressed star formation in passive galaxies but that the presence of an AGN doesn't necessarily result in the suppression of star formation. △ Less

Submitted 9 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 17 pages, 13 figures, to be submitted to the Open Journal of Astrophysics. The first two authors are joint primary authors. Updated with new observational constraints

The JWST-PRIMAL Legacy Survey. A JWST/NIRSpec reference sample for the physical properties and Lyman-$α$ absorption and emission of $\sim 500$ galaxies at $z=5.5-13.4$

Authors: K. E. Heintz, G. B. Brammer, D. Watson, P. A. Oesch, L. C. Keating, M. J. Hayes, Abdurro'uf, K. Z. Arellano-Córdova, A. C. Carnall, C. R. Christiansen, F. Cullen, R. Davé, P. Dayal, A. Ferrara, K. Finlator, J. P. U. Fynbo, S. R. Flury, V. Gelli, S. Gillman, R. Gottumukkala, K. Gould, T. R. Greve, S. E. Hardin, T. Y. -Y Hsiao, A. Hutter , et al. (23 additional authors not shown)

Abstract: One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$α$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$α$ absorption (DLA) wings from high column densities of neu… ▽ More One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$α$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$α$ absorption (DLA) wings from high column densities of neutral atomic hydrogen (HI), signifying major gas accretion events in the formation of these galaxies. To explore this new phenomenon systematically, we assemble the JWST/NIRSpec PRImordial gas Mass AssembLy (PRIMAL) legacy survey of 494 galaxies at $z=5.5-13.4$. We characterize this benchmark sample in full and spectroscopically derive the galaxy redshifts, metallicities, star-formation rates, and ultraviolet slopes. We define a new diagnostic, the Ly$α$ damping parameter $D_{\rm Lyα}$ to measure and quantify the Ly$α$ emission strength, HI fraction in the IGM, or local HI column density for each source. The JWST-PRIMAL survey is based on the spectroscopic DAWN JWST Archive (DJA-Spec). All the software, reduced spectra, and spectroscopically derived quantities and catalogs are made publicly available in dedicated repositories. The fraction of strong galaxy DLAs are found to be in the range $65-95\%$ at $z>5.5$. The fraction of strong Ly$α$ emitters (LAEs) is found to increase with decreasing redshift, in qualitative agreement with previous observational results, and are predominantly associated with low-metallicity and UV faint galaxies. By contrast, strong DLAs are observed in galaxies with a variety of intrinsic physical properties. Our results indicate that strong DLAs likely reflect a particular early assembly phase of reionization-era galaxies, at which point they are largely dominated by pristine HI gas accretion. [abridged] △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 20 pages, 18 figures, 3 tables. Submitted to A&A. Comments welcome! All data and catalogs are available through the DAWN JWST Archive (DJA): https://dawn-cph.github.io/dja/ and https://github.com/keheintz/jwst-primal

EPOCHS IV: SED Modelling Assumptions and their impact on the Stellar Mass Function at 6.5 < z < 13.5 using PEARLS and public JWST observations

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

Abstract: 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<z<13.5 taken from a consistent reduction and analysis of publicly available deep JWST NIRCam data covering the PEARLS, CEERS, GLASS, JADES GOOD-S, NGDEEP, and SMACS0723 surveys, totalling 187 arcmin2. We investigat… ▽ More 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<z<13.5 taken from a consistent reduction and analysis of publicly available deep JWST NIRCam data covering the PEARLS, CEERS, GLASS, JADES GOOD-S, NGDEEP, and SMACS0723 surveys, totalling 187 arcmin2. We investigate the impact of SED fitting methods, assumed star formation histories (SFH), dust laws, and priors on galaxy masses and the resultant GSMF. Whilst our fiducial GSMF agrees with the literature at z<13.5, we find that the assumed SFH model has a large impact on the GSMF and stellar mass density (SMD), finding a 0.75 dex increase in the SMD at z=10.5 between a flexible non-parametric and standard parametric SFH. Overall, we find a flatter SMD evolution at z > 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. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 40 pages, 15 figures. Submitted to ApJ. Data is avaible on Github: https://github.com/tHarvey303/EpochsIV. Comments to corresponding author welcome at [email protected]

Tracing the rise of supermassive black holes: A panchromatic search for faint, unobscured quasars at z > 6 with COSMOS-Web and other surveys

Authors: Irham T. Andika, Knud Jahnke, Masafusa Onoue, John D. Silverman, Itsna K. Fitriana, Angela Bongiorno, Malte Brinch, Caitlin M. Casey, Andreas Faisst, Steven Gillman, Ghassem Gozaliasl, Christopher C. Hayward, Michaela Hirschmann, Dale Kocevski, Anton M. Koekemoer, Vasily Kokorev, Erini Lambrides, Minju M. Lee, R. Michael Rich, Benny Trakhtenbrot, C. Megan Urry, Stephen M. Wilkins, Aswin P. Vijayan

Abstract: We report the identification of 64 new candidates of compact galaxies, potentially hosting faint quasars with bolometric luminosities of $L_\mathrm{bol} = 10^{43}$--10$^{46}$ erg s$^{-1}$, residing in the reionization epoch within the redshift range of $6 \lesssim z \lesssim 8$. These candidates were selected by harnessing the rich multiband datasets provided by the emerging JWST-driven extragalac… ▽ More We report the identification of 64 new candidates of compact galaxies, potentially hosting faint quasars with bolometric luminosities of $L_\mathrm{bol} = 10^{43}$--10$^{46}$ erg s$^{-1}$, residing in the reionization epoch within the redshift range of $6 \lesssim z \lesssim 8$. These candidates were selected by harnessing the rich multiband datasets provided by the emerging JWST-driven extragalactic surveys, focusing on COSMOS-Web, as well as JADES, UNCOVER, CEERS, and PRIMER. Our search strategy includes two stages: applying stringent photometric cuts to catalog-level data and detailed spectral energy distribution fitting. These techniques effectively isolate the quasar candidates while mitigating contamination from low-redshift interlopers, such as brown dwarfs and nearby galaxies. The selected candidates indicate physical traits compatible with low-luminosity active galactic nuclei, likely hosting $\approx10^5$--$10^7~M_\odot$ supermassive black holes (SMBHs) living in galaxies with stellar masses of $\approx10^8$--$10^{10}~M_\odot$. The SMBHs selected in this study, on average, exhibit an elevated mass compared to their hosts, with the mass ratio distribution slightly higher than those of galaxies in the local Universe. As with other high-$z$ studies, this is at least in part due to the selection method for these quasars. An extensive Monte Carlo analysis provides compelling evidence that heavy black hole seeds from the direct collapse scenario appear to be the preferred pathway to mature this specific subset of SMBHs by $z\approx7$. This work underscores the significance of further spectroscopic observations, as the quasar candidates presented here offer exceptional opportunities to delve into the nature of the earliest galaxies and SMBHs that formed during cosmic infancy. △ Less

Submitted 2 February, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astronomy & Astrophysics journal. 19 pages, 10 figures, and 4 tables. We welcome comments from the reader

Cosmic Evolution Early Release Science (CEERS) survey: The colour evolution of galaxies in the distant Universe

Authors: Stephen M. Wilkins, Jack C. Turner, Micaela B. Bagley, Steven L. Finkelstein, Ricardo O. Amorín, Adrien Aufan Stoffels D Hautefort, Peter Behroozi, Rachana Bhatawdekar, Avishai Dekel, James Donnellan, Nicole E. Drakos, Flaminia Fortuni, Nimish P. Hathi, Michaela Hirschmann, Benne W. Holwerda, Dimitrios Irodotou, Anton M. Koekemoer, Christopher C. Lovell, Emiliano Merlin, Will J. Roper, Louise T. C. Seeyave, Aswin P. Vijayan, L. Y. Aaron Yung

Abstract: 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 g… ▽ More 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 $5<z<10$ galaxies using the Cosmic Evolution Early Release Science (CEERS) Survey. We then analyse predictions from a variety of models using the same methodology and compare the NIRCam/F277W magnitude distribution and NIRCam colours with observations. We find that the predicted and observed magnitude distributions are similar, at least at $5<z<8$. At $z>8$ 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 $5<z<8$. There is however some disagreement between the observed and modelled strength of the strong line contribution. In particular all the models fails to reproduce the F410M-F444W colour at $z>8$, though, again, the sample size is small here. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 11 pages, 10 figures, submitted to MNRAS

Cosmic Vine: A z=3.44 large-scale structure hosting massive quiescent galaxies

Authors: Shuowen Jin, Nikolaj B. Sillassen, Georgios E. Magdis, Malte Brinch, Marko Shuntov, Gabriel Brammer, Raphael Gobat, Francesco Valentino, Adam C. Carnall, Minju Lee, Aswin P. Vijayan, Steven Gillman, Vasily Kokorev, Aurélien Le Bail, Thomas R. Greve, Bitten Gullberg, Katriona M. L. Gould, Sune Toft

Abstract: We report the discovery of a large-scale structure at z=3.44 revealed by JWST data in the Extended Groth Strip (EGS) field. This structure, called the Cosmic Vine, consists of 20 galaxies with spectroscopic redshifts at 3.43<z<3.45 and six galaxy overdensities ($4-7σ$) with consistent photometric redshifts, making up a vine-like structure extending over a ~4x0.2 pMpc^2 area. The two most massive g… ▽ More We report the discovery of a large-scale structure at z=3.44 revealed by JWST data in the Extended Groth Strip (EGS) field. This structure, called the Cosmic Vine, consists of 20 galaxies with spectroscopic redshifts at 3.43<z<3.45 and six galaxy overdensities ($4-7σ$) with consistent photometric redshifts, making up a vine-like structure extending over a ~4x0.2 pMpc^2 area. The two most massive galaxies ($M_*\approx10^{10.9}~M_\odot$) of the Cosmic Vine are found to be quiescent with bulge-dominated morphologies ($B/T>70\%$). Comparisons with simulations suggest that the Cosmic Vine would form a cluster with halo mass $M_{\rm halo}>10^{14}M_\odot$ at z=0, and the two massive galaxies are likely forming the brightest cluster galaxies (BCGs). The results unambiguously reveal that massive quiescent galaxies can form in growing large-scale structures at z>3, thus disfavoring the environmental quenching mechanisms that require a virialized cluster core. Instead, as suggested by the interacting and bulge-dominated morphologies, the two galaxies are likely quenched by merger-triggered starburst or active galactic nucleus (AGN) feedback before falling into a cluster core. Moreover, we found that the observed specific star formation rates of massive quiescent galaxies in z>3 dense environments are one to two orders of magnitude lower than that of the BCGs in the TNG300 simulation. This discrepancy potentially poses a challenge to the models of massive cluster galaxy formation. Future studies comparing a large sample with dedicated cluster simulations are required to solve the problem. △ Less

Submitted 18 February, 2024; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: A&A Letters in press, open access at https://doi.org/10.1051/0004-6361/202348540

DEIMOS spectroscopy of $z=6$ protocluster candidate in COSMOS -- A massive protocluster embedded in a large scale structure?

Authors: Malte Brinch, Thomas R. Greve, David B. Sanders, Conor J. R. McPartland, Nima Chartab, Steven Gillman, Aswin P. Vijayan, Minju M. Lee, Gabriel Brammer, Caitlin M. Casey, Olivier Ilbert, Shuowen Jin, Georgios Magdis, H. J. McCracken, Nikolaj B. Sillassen, Sune Toft, Jorge A. Zavala

Abstract: We present the results of our Keck/DEIMOS spectroscopic follow-up of candidate galaxies of i-band-dropout protocluster candidate galaxies at $z\sim6$ in the COSMOS field. We securely detect Lyman-$α$ emission lines in 14 of the 30 objects targeted, 10 of them being at $z=6$ with a signal-to-noise ratio of $5-20$, the remaining galaxies are either non-detections or interlopers with redshift too dif… ▽ More We present the results of our Keck/DEIMOS spectroscopic follow-up of candidate galaxies of i-band-dropout protocluster candidate galaxies at $z\sim6$ in the COSMOS field. We securely detect Lyman-$α$ emission lines in 14 of the 30 objects targeted, 10 of them being at $z=6$ with a signal-to-noise ratio of $5-20$, the remaining galaxies are either non-detections or interlopers with redshift too different from $z=6$ to be part of the protocluster. The 10 galaxies at $z\approx6$ make the protocluster one of the riches at $z>5$. The emission lines exhibit asymmetric profiles with high skewness values ranging from 2.87 to 31.75, with a median of 7.37. This asymmetry is consistent with them being Ly$α$, resulting in a redshift range of $z=5.85-6.08$. Using the spectroscopic redshifts, we re-calculate the overdensity map for the COSMOS field and find the galaxies to be in a significant overdensity at the $4σ$ level, with a peak overdensity of $δ=11.8$ (compared to the previous value of $δ=9.2$). The protocluster galaxies have stellar masses derived from Bagpipes SED fits of $10^{8.29}-10^{10.28} \rm \,M_{\rm \odot}$ and star formation rates of $2-39\,\rm M_{\rm \odot}\rm\,yr^{-1}$, placing them on the main sequence at this epoch. Using a stellar-to-halo-mass relationship, we estimate the dark matter halo mass of the most massive halo in the protocluster to be $\sim 10^{12}\rm M_{\rm \odot}$. By comparison with halo mass evolution tracks from simulations, the protocluster is expected to evolve into a Virgo- or Coma-like cluster in the present day. △ Less

Submitted 18 December, 2023; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 26 pages, 14 figues, 5 tables, main text is 16 pages, appendix is 10 pages, published in MNRAS

The impact of binary stars on the dust and metal evolution of galaxies

Authors: Robert M. Yates, David Hendriks, Aswin P. Vijayan, Robert G. Izzard, Peter A. Thomas, Payel Das

Abstract: We present detailed implementations of (a) binary stellar evolution (using binary_c) and (b) dust production and destruction into the cosmological semi-analytic galaxy evolution simulation, L-Galaxies. This new version of L-Galaxies is compared to a version assuming only single stars and to global and spatially-resolved observational data across a range of redshifts ($z$). We find that binaries ha… ▽ More We present detailed implementations of (a) binary stellar evolution (using binary_c) and (b) dust production and destruction into the cosmological semi-analytic galaxy evolution simulation, L-Galaxies. This new version of L-Galaxies is compared to a version assuming only single stars and to global and spatially-resolved observational data across a range of redshifts ($z$). We find that binaries have a negligible impact on the stellar masses, gas masses, and star formation rates of galaxies only if the total mass ejected by massive stars is unchanged. This is because massive stars determine the strength of supernova (SN) feedback, which in turn regulates galaxy growth. Binary effects, such as common envelope ejection and novae, affect carbon and nitrogen enrichment in galaxies, however heavier alpha elements are more affected by the choice of SN and wind yields. Unlike many other simulations, the new L-Galaxies reproduces observed dust-to-metal (DTM) and dust-to-gas (DTG) ratios at $z\sim{}0-4$. This is mainly due to shorter dust accretion timescales in dust-rich environments. However, dust masses are under-predicted at $z>4$, highlighting the need for enhanced dust production at early times in simulations, possibly accompanied by increased star formation. On sub-galactic scales, there is very good agreement between L-Galaxies and observed dust and metal radial profiles at $z=0$. A drop in DTM ratio is also found in diffuse, low-metallicity regions, contradicting the assumption of a universal value. We hope that this work serves as a useful template for binary stellar evolution implementations in other cosmological simulations in future. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 19 pages, 10 figures. Submitted to MNRAS. Comments welcome!

Industry 4.0 and Beyond: The Role of 5G, WiFi 7, and TSN in Enabling Smart Manufacturing

Authors: Jobish John, Md. Noor-A-Rahim, Aswathi Vijayan, H. Vincent Poor, Dirk Pesch

Abstract: This paper explores the role that 5G, WiFi-7, and Time-Sensitive Networking (TSN) can play in driving smart manufacturing as a fundamental part of the Industry 4.0 vision. The paper provides an in-depth analysis of each technology's application in industrial communications, with a focus on TSN and its key elements that enable reliable and secure communication in industrial networks. In addition, t… ▽ More This paper explores the role that 5G, WiFi-7, and Time-Sensitive Networking (TSN) can play in driving smart manufacturing as a fundamental part of the Industry 4.0 vision. The paper provides an in-depth analysis of each technology's application in industrial communications, with a focus on TSN and its key elements that enable reliable and secure communication in industrial networks. In addition, the paper includes a comparative study of these technologies, analyzing them based on a number of industrial use-cases, supported secondary applications, industry adoption, and current market trends. The paper concludes by highlighting the challenges and future directions for the adoption of these technologies in industrial networks and emphasizes their importance in realizing the Industry 4.0 vision within the context of smart manufacturing. △ Less

Submitted 3 October, 2023; originally announced October 2023.

The COSMOS-Web ring: in-depth characterization of an Einstein ring lensing system at z~2

Authors: W. Mercier, M. Shuntov, R. Gavazzi, J. W. Nightingale, R. Arango, O. Ilbert, A. Amvrosiadis, L. Ciesla, C. Casey, S. Jin, A. L. Faisst, I. T. Andika, N. E. Drakos, A. Enia, M. Franco, S. Gillman, G. Gozaliasl, C. C. Hayward, M. Huertas-Company, J. S. Kartaltepe, A. M. Koekemoer, C. Laigle, D. Le Borgne, G. Magdis, G. Mahler , et al. (12 additional authors not shown)

Abstract: Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z=2 that we serendipitously discovered in the COSMOS-Web survey and possibly the most distant lens discovered to date. Methods. We extract the visible and NIR photometry from more than 25 bands and we derive the photometric redshifts and physical properties of both the lens and the source with three different SED f… ▽ More Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z=2 that we serendipitously discovered in the COSMOS-Web survey and possibly the most distant lens discovered to date. Methods. We extract the visible and NIR photometry from more than 25 bands and we derive the photometric redshifts and physical properties of both the lens and the source with three different SED fitting codes. Using JWST/NIRCam images, we also produce two lens models to (i) recover the total mass of the lens, (ii) derive the magnification of the system, (iii) reconstruct the morphology of the lensed source, and (iv) measure the slope of the total mass density profile of the lens. Results. The lens is a very massive and quiescent (sSFR < 10^(-13) yr-1) elliptical galaxy at z = 2.02 \pm 0.02 with a total mass Mtot(<thetaE) = (3.66 \pm 0.36) x 10^11 Msun and a stellar mass M* = (1.37 \pm 0.14) x 10^11 Msun. Compared to SHMRs from the literature, we find that the total mass is consistent with the presence of a DM halo of mass Mh = 1.09^(+1.46)_(-0.57) x 10^13 Msun. In addition, the background source is a M* = (1.26 \pm 0.17) x 10^10 Msun star-forming galaxy (SFR=(78 \pm 15) Msun/yr) at z = 5.48 \pm 0.06. Its reconstructed morphology shows two components with different colors. Dust attenuation values from SED fitting and nearby detections in the FIR also suggest it could be partially dust-obscured. Conclusions. We find the lens at z=2. Its total, stellar, and DM halo masses are consistent within the Einstein ring, so we do not need any unexpected changes in our description of the lens (e.g. change its IMF or include a non-negligible gas contribution). The most likely solution for the lensed source is at z = 5.5. Its reconstructed morphology is complex and highly wavelength dependent, possibly because it is a merger or a main sequence galaxy with a heterogeneous dust distribution. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 16 pages, submitted to A&A

Journal ref: A&A 687, A61 (2024)

Quokka-based Understanding of Outflows (QED) I. Metal loading, phase structure, and convergence testing for Solar neighbourhood conditions

Authors: Aditi Vijayan, Mark R. Krumholz, Benjamin D. Wibking

Abstract: Multiphase galactic outflows, generated by supernova feedback, are likely to be more metal-rich than the interstellar media from which they are driven due to incomplete mixing between supernova ejecta and the ambient ISM. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to res… ▽ More Multiphase galactic outflows, generated by supernova feedback, are likely to be more metal-rich than the interstellar media from which they are driven due to incomplete mixing between supernova ejecta and the ambient ISM. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to resolve mass, momentum and energy exchanges between the different phases of the outflows. In this context, we present simulations of outflows, driven by SN feedback, conducted using \textsc{Quokka}, a new GPU-optimised AMR radiation-hydrodynamics code. This code allows us to reach combinations of resolution, simulation volume, and simulation duration larger than those that have previously been possible, and to resolve all gas phases from cold neutral medium, $T \sim 100$ K, to hot ionised gas, $T \gtrsim 10^7$ K. In this, a first of a series of papers exploring generation and evolution of multiphase outflows from a wide range of galactic environments and star formation rates, we quantify the extent of selective metal loading in Solar neighbourhood-like environments. We explain the selective metal loading we find as a result of the transport of metals within and between phases, a phenomenon we can study owing to the parsec-scale resolution that our simulations achieve. We also quantify the sensitivity of metal loading studies to numerical resolution, and present convergence criteria for future studies. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS. Full title="Quokka-based Understanding of Outflows Derived from Extensive, Repeated, Accurate, Thorough, Demanding, Expensive, Memory-consuming, Ongoing Numerical Simulations of Transport, Removal, Accretion, Nucleosynthesis, Deposition, and Uplifting of Metals (QUOD ERAT DEMONSTRANDUM, i.e, QED)"

COSMOS-Web: Intrinsically Luminous z$\gtrsim$10 Galaxy Candidates Test Early Stellar Mass Assembly

Authors: Caitlin M. Casey, Hollis B. Akins, Marko Shuntov, Olivier Ilbert, Louise Paquereau, Maximilien Franco, Christopher C. Hayward, Steven L. Finkelstein, Michael Boylan-Kolchin, Brant E. Robertson, Natalie Allen, Malte Brinch, Olivia R. Cooper, Xuheng Ding, Nicole E. Drakos, Andreas L. Faisst, Seiji Fujimoto, Steven Gillman, Santosh Harish, Michaela Hirschmann, Shuowen Jin, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Vasily Kokorev, Daizhong Liu , et al. (17 additional authors not shown)

Abstract: We report the discovery of 15 exceptionally luminous $10\lesssim z\lesssim14$ candidate galaxies discovered in the first 0.28 deg$^2$ of JWST/NIRCam imaging from the COSMOS-Web Survey. These sources span rest-frame UV magnitudes of $-20.5>M_{\rm UV}>-22$, and thus constitute the most intrinsically luminous $z\gtrsim10$ candidates identified by JWST to-date. Selected via NIRCam imaging with Hubble… ▽ More We report the discovery of 15 exceptionally luminous $10\lesssim z\lesssim14$ candidate galaxies discovered in the first 0.28 deg$^2$ of JWST/NIRCam imaging from the COSMOS-Web Survey. These sources span rest-frame UV magnitudes of $-20.5>M_{\rm UV}>-22$, and thus constitute the most intrinsically luminous $z\gtrsim10$ candidates identified by JWST to-date. Selected via NIRCam imaging with Hubble ACS/F814W, deep ground-based observations corroborate their detection and help significantly constrain their photometric redshifts. We analyze their spectral energy distributions using multiple open-source codes and evaluate the probability of low-redshift solutions; we conclude that 12/15 (80%) are likely genuine $z\gtrsim10$ sources and 3/15 (20%) likely low-redshift contaminants. Three of our $z\sim12$ candidates push the limits of early stellar mass assembly: they have estimated stellar masses $\sim5\times10^{9}\,M_\odot$, implying an effective stellar baryon fraction of $ε_{\star}\sim0.2-0.5$, where $ε_{\star}\equiv M_{\star}/(f_{b}M_{halo})$. The assembly of such stellar reservoirs is made possible due to rapid, burst-driven star formation on timescales $<$100\,Myr where the star-formation rate may far outpace the growth of the underlying dark matter halos. This is supported by the similar volume densities inferred for $M_\star\sim10^{10}\,M_\odot$ galaxies relative to $M_\star\sim10^{9}\,M_\odot$ -- both about $10^{-6}$ Mpc$^{-3}$ -- implying they live in halos of comparable mass. At such high redshifts, the duty cycle for starbursts would be of order unity, which could cause the observed change in the shape of the UVLF from a double powerlaw to Schechter at $z\approx8$. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understanding how, and if, such early massive galaxies push the limits of galaxy formation in $Λ$CDM. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 30 pages, 9 figures; ApJ submitted

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

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

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

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

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

Journal ref: MNRAS, 527, 11627 (2024)

Unveiling the distant Universe: Characterizing $z\ge9$ Galaxies in the first epoch of COSMOS-Web

Authors: Maximilien Franco, Hollis B. Akins, Caitlin M. Casey, Steven L. Finkelstein, Marko Shuntov, Katherine Chworowsky, Andreas L. Faisst, Seiji Fujimoto, Olivier Ilbert, Anton M. Koekemoer, Daizhong Liu, Christopher C. Lovell, Claudia Maraston, Henry Joy McCracken, Jed McKinney, Brant E. Robertson, Micaela B. Bagley, Jaclyn B. Champagne, Olivia R. Cooper, Xuheng Ding, Nicole E. Drakos, Andrea Enia, Steven Gillman, Christopher C. Hayward, Michaela Hirschmann , et al. (25 additional authors not shown)

Abstract: We report the identification of 15 galaxy candidates at $z\ge9$ using the initial COSMOS-Web JWST observations over 77 arcmin$^2$ through four NIRCam filters (F115W, F150W, F277W, F444W) with an overlap with MIRI (F770W) of 8.7 arcmin$^2$. We fit the sample using several publicly-available SED fitting and photometric redshift codes and determine their redshifts between $z=9.3$ and $z=10.9$ (… ▽ More We report the identification of 15 galaxy candidates at $z\ge9$ using the initial COSMOS-Web JWST observations over 77 arcmin$^2$ through four NIRCam filters (F115W, F150W, F277W, F444W) with an overlap with MIRI (F770W) of 8.7 arcmin$^2$. We fit the sample using several publicly-available SED fitting and photometric redshift codes and determine their redshifts between $z=9.3$ and $z=10.9$ ($\langle z\rangle=10.0$), UV-magnitudes between M$_{\rm UV}$ = $-$21.2 and $-$19.5 (with $\langle $M$_{\rm UV}\rangle=-20.2$) and rest-frame UV slopes ($\langle β\rangle=-2.4$). These galaxies are, on average, more luminous than most $z\ge9$ candidates discovered by JWST so far in the literature, while exhibiting similar blue colors in their rest-frame UV. The rest-frame UV slopes derived from SED-fitting are blue ($β\sim$[$-$2.0, $-$2.7]) without reaching extremely blue values as reported in other recent studies at these redshifts. The blue color is consistent with models that suggest the underlying stellar population is not yet fully enriched in metals like similarly luminous galaxies in the lower redshift Universe. The derived stellar masses with $\langle \log_{\rm 10} ($M$_\star/$M$_\odot)\rangle\approx8-9$ are not in tension with the standard $Λ$CDM model and our measurement of the volume density of such UV luminous galaxies aligns well with previously measured values presented in the literature at $z\sim9-10$. Our sample of galaxies, although compact, are significantly resolved. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: Submitted to ApJ

Size - Stellar Mass Relation and Morphology of Quiescent Galaxies at $z\geq3$ in Public $JWST$ Fields

Authors: Kei Ito, Francesco Valentino, Gabriel Brammer, Andreas L. Faisst, Steven Gillman, Carlos Gomez-Guijarro, Katriona M. L. Gould, Kasper E. Heintz, Olivier Ilbert, Christian Kragh Jespersen, Vasily Kokorev, Mariko Kubo, Georgios E. Magdis, Conor McPartland, Masato Onodera, Francesca Rizzo, Masayuki Tanaka, Sune Toft, Aswin P. Vijayan, John R. Weaver, Katherine E. Whitaker, Lillian Wright

Abstract: We present the results of a systematic study of the rest-frame optical morphology of quiescent galaxies at $z \geq 3$ using the Near-Infrared Camera (NIRCam) onboard $JWST$. Based on a sample selected by $UVJ$ color or $NUVUVJ$ color, we focus on 26 quiescent galaxies with $9.8<\log{(M_\star/M_\odot)}<11.4$ at $2.8<z_{\rm phot}<4.6$ with publicly available $JWST$ data. Their sizes are constrained… ▽ More We present the results of a systematic study of the rest-frame optical morphology of quiescent galaxies at $z \geq 3$ using the Near-Infrared Camera (NIRCam) onboard $JWST$. Based on a sample selected by $UVJ$ color or $NUVUVJ$ color, we focus on 26 quiescent galaxies with $9.8<\log{(M_\star/M_\odot)}<11.4$ at $2.8<z_{\rm phot}<4.6$ with publicly available $JWST$ data. Their sizes are constrained by fitting the Sérsic profile to all available NIRCam images. We see a negative correlation between the observed wavelength and the size in our sample and derive their size at the rest-frame $0.5\, {\rm μm}$ taking into account this trend. Our quiescent galaxies show a significant correlation between the rest-frame $0.5\, {\rm μm}$ size and the stellar mass at $z\geq3$. The analytical fit for them at $\log{(M_\star/M_\odot)}>10.3$ implies that our size - stellar mass relations are below those at lower redshifts, with the amplitude of $\sim0.6\, {\rm kpc}$ at $M_\star = 5\times 10^{10}\, M_\odot$. This value agrees with the extrapolation from the size evolution of quiescent galaxies at $z<3$ in the literature, implying that the size of quiescent galaxies increases monotonically from $z\sim3-5$. Our sample is mainly composed of galaxies with bulge-like structures according to their median Sérsic index and axis ratio of $n\sim3-4$ and $q\sim0.6-0.8$, respectively. On the other hand, there is a trend of increasing fraction of galaxies with low Sérsic index, suggesting $3<z<5$ might be the epoch of onset of morphological transformation with a fraction of very notable disky quenched galaxies. △ Less

Submitted 6 February, 2024; v1 submitted 13 July, 2023; originally announced July 2023.

Comments: 25 pages, 17 figures, 3 tables; accepted for publication in ApJ

First Light And Reionisation Epoch Simulations (FLARES) XIV: The Balmer/4000~Å Breaks of Distant Galaxies

Authors: Stephen M. Wilkins, Christopher C. Lovell, Dimitrios Irodotou, Aswin P. Vijayan, Anton Vikaeus, Erik Zackrisson, Joseph Caruana, Elizabeth R. Stanway, Christopher J. Conselice, Louise T. C. Seeyave, William J. Roper, Katherine Chworowsky, Steven L. Finkelstein

Abstract: With the successful launch and commissioning of JWST we are now able to routinely spectroscopically probe the rest-frame optical emission of galaxies at $z>6$ for the first time. Amongst the most useful spectral diagnostics used in the optical is the Balmer/4000~Å break; this is, in principle, a diagnostic of the mean ages of composite stellar populations. However, the Balmer break is also sensiti… ▽ More With the successful launch and commissioning of JWST we are now able to routinely spectroscopically probe the rest-frame optical emission of galaxies at $z>6$ for the first time. Amongst the most useful spectral diagnostics used in the optical is the Balmer/4000~Å break; this is, in principle, a diagnostic of the mean ages of composite stellar populations. However, the Balmer break is also sensitive to the shape of the star formation history, the stellar (and gas) metallicity, the presence of nebular continuum emission, and dust attenuation. In this work we explore the origin of the Balmer/4000~Å break using the SYNTHESIZER synthetic observations package. We then make predictions of the Balmer/4000~Å break using the First Light and Reionisation Epoch Simulations (FLARES) at $5<z<10$. We find that the average break strength weakly correlates with stellar mass and rest-frame far-UV luminosity, but that this is predominantly driven by dust attenuation. We also find that break strength provides a weak diagnostic of the age but performs better as a means to constrain star formation and stellar mass, alongside the UV and optical luminosity, respectively. △ Less

Submitted 29 May, 2023; originally announced May 2023.

Comments: 9 pages, 9 figures, submitted to MNRAS

First Light And Reionisation Epoch Simulations (FLARES) XIII: The Lyman-continuum emission of high-redshift galaxies

Authors: Louise T. C. Seeyave, Stephen M. Wilkins, Jussi K. Kuusisto, Christopher C. Lovell, Dimitrios Irodotou, Charlotte Simmonds, Aswin P. Vijayan, Peter A. Thomas, William J. Roper, Conor M. Byrne, Gareth T. Jones, Jack C. Turner, Christopher J. Conselice

Abstract: The history of reionisation is highly dependent on the ionising properties of high-redshift galaxies. It is therefore important to have a solid understanding of how the ionising properties of galaxies are linked to physical and observable quantities. In this paper, we use the First Light and Reionisation Epoch Simulations (FLARES) to study the Lyman-continuum (LyC, i.e. hydrogen-ionising) emission… ▽ More The history of reionisation is highly dependent on the ionising properties of high-redshift galaxies. It is therefore important to have a solid understanding of how the ionising properties of galaxies are linked to physical and observable quantities. In this paper, we use the First Light and Reionisation Epoch Simulations (FLARES) to study the Lyman-continuum (LyC, i.e. hydrogen-ionising) emission of massive ($M_*>10^8\,\mathrm{M_\odot}$) galaxies at redshifts $z=5-10$. We find that the specific ionising emissivity (i.e. intrinsic ionising emissivity per unit stellar mass) decreases as stellar mass increases, due to the combined effects of increasing age and metallicity. FLARES predicts a median ionising photon production efficiency (i.e. intrinsic ionising emissivity per unit intrinsic far-UV luminosity) of $\log_{10}(ξ_{\rm ion}\rm{/erg^{-1}Hz})=25.40^{+0.16}_{-0.17}$, with values spanning the range $\log_{10}(ξ_{\rm ion}\rm{/erg^{-1}Hz})=25-25.75$. This is within the range of many observational estimates, but below some of the extremes observed. We compare the production efficiency with observable properties, and find a weak negative correlation with the UV-continuum slope, and a positive correlation with the OIII equivalent width. We also consider the dust-attenuated production efficiency (i.e. intrinsic ionising emissivity per unit dust-attenuated far-UV luminosity), and find a median of $\log_{10}(ξ_{\rm ion}\rm{/erg^{-1}Hz})\sim25.5$. Within our sample of $M_*>10^8\,\mathrm{M_\odot}$ galaxies, it is the stellar populations in low mass galaxies that contribute the most to the total ionising emissivity. Active galactic nuclei (AGN) emission accounts for $10-20$ % of the total emissivity at a given redshift, and extends the LyC luminosity function by $\sim0.5$ dex. △ Less

Submitted 21 August, 2023; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: 19 pages, 17 figures, accepted for publication in MNRAS

Sub-Millimetre Galaxies with Webb: Near-Infrared Counterparts and Multi-wavelength Morphology

Authors: S. Gillman, B. Gullberg, G. Brammer, A. Vijayan, M. Lee, D. Blánquez, M. Brinch, T. Greve, I. Jermann, S. Jin, V. Kokorev, L. Liu, G. Magdis, F. Rizzo, F. Valentino

Abstract: We utilise the unprecedented depth and resolution of recent early-release science (ERS) JWST observations to define the near-infrared counterparts of sub-millimetre selected galaxies (SMGs). We identify 45 SCUBA-2 SMG positions within The Cosmic Evolution Early Release Science Survey (CEERS) JWST/NIRCam fields. Through an analysis of multi-wavelength $p$-values, NIRCam colours and predicted SCUBA-… ▽ More We utilise the unprecedented depth and resolution of recent early-release science (ERS) JWST observations to define the near-infrared counterparts of sub-millimetre selected galaxies (SMGs). We identify 45 SCUBA-2 SMG positions within The Cosmic Evolution Early Release Science Survey (CEERS) JWST/NIRCam fields. Through an analysis of multi-wavelength $p$-values, NIRCam colours and predicted SCUBA-2 fluxes, we define 43 JWST/NIRCam counterparts to the SCUBA-2 SMGs, finding a 63 per cent agreement with those identified in prior $HST$ studies. Using EaZy-py we fit the available HST and JWST observations to quantify the photometric redshifts of the NIRCam-SMGs, establishing a broad range of redshift from $z$$\approx$0.2$-$5.4 with a median of $z$$\approx$2.29, in agreement with other studies of SMGs. We analyse their rest-frame optical and near-infrared morphological properties (e.g. effective radius (R$_{\rm e}$), Sérsic index ($n$), CAS, Gini and M$_{20}$), finding, on average, late-type disc-like morphologies with large scatter into the intermediate and merger regions of the non-parametric parameter space. For the non-merging galaxies, we find a median rest-frame optical size and Sérsic index (and $1σ$ scatter) of R$_{\rm e}$=3.10$\pm$1.67kpc and $n$=0.96$\pm$0.66. Whilst in the rest-frame near-infrared we establish more compact, higher Sérsic index morphologies (R$_{\rm e}$=1.64$\pm$0.97, $n$=1.85$\pm$0.63). We further establish that both the rest-frame optical and near-infrared effective radii correlate negatively (at a 2$σ$ level) with redshift whilst the Sérsic index remains constant with cosmic time. Our results are consistent with the picture of inside-out galaxy evolution, with more centrally concentrated older stellar populations, and more extended, younger star-forming regions whose stellar emission is heavily attenuated in the central regions. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: Submitted for publication in A&A. Comments welcome. 18 pages

Journal ref: A&A 676, A26 (2023)

Galaxy Evolution in $\ddotμ$ based Cosmologies

Authors: Will J. Roper, Stephen M. Wilkins, Stephen Riggs, Jessica Pilling, Aswin P. Vijayan, Dimitrios Irodotou, Violetta Korbina, Jussi Kuusisto

Abstract: We present the first study of galaxy evolution in $\ddotμ$ based cosmologies. We find that recent JWST observations of massive galaxies at extremely high redshifts are consistent with such a cosmology. However, the low redshift Universe is entirely divergent from the $\ddotμ$ cosmic star formation rate density. We thus propose that our Universe was at one point dominated by a Primordial Bovine Her… ▽ More We present the first study of galaxy evolution in $\ddotμ$ based cosmologies. We find that recent JWST observations of massive galaxies at extremely high redshifts are consistent with such a cosmology. However, the low redshift Universe is entirely divergent from the $\ddotμ$ cosmic star formation rate density. We thus propose that our Universe was at one point dominated by a Primordial Bovine Herd (PBH) which later decayed producing dark energy. Note that we do not detail the mechanisms by which this decay process takes place. Despite its vanishingly small probability for existence, a $\ddotμ$ based cosmological model marries the disparate findings in the high and low redshift Universe. △ Less

Submitted 29 March, 2023; originally announced March 2023.

First Light And Reionisation Epoch Simulations (FLARES) XII: The consequences of star-dust geometry on galaxies in the EoR

Authors: Aswin P. Vijayan, Peter A. Thomas, Christopher C. Lovell, Stephen M. Wilkins, Thomas R. Greve, Dimitrios Irodotou, William J. Roper, Louise T. C. Seeyave

Abstract: Using the First Light And Reionisation Epoch Simulations (${\rm F{\small LARES}}$), a suite of hydrodynamical simulations we explore the consequences of a realistic model for star--dust geometry on the observed properties of galaxies. We find that the UV attenuation declines rapidly from the central regions of galaxies, and bright galaxies have spatially extended star formation that suffers less o… ▽ More Using the First Light And Reionisation Epoch Simulations (${\rm F{\small LARES}}$), a suite of hydrodynamical simulations we explore the consequences of a realistic model for star--dust geometry on the observed properties of galaxies. We find that the UV attenuation declines rapidly from the central regions of galaxies, and bright galaxies have spatially extended star formation that suffers less obscuration than their fainter counterparts, demonstrating a non-linear relationship between the UV luminosity and the UV attenuation, giving a double power-law shape to the UVLF. Spatially distinct stellar populations within galaxies experience a wide range of dust attenuation due to variations in the dust optical depth along their line-of-sight; which can range from completely dust obscured to being fully unobscured. The overall attenuation curve of a galaxy is then a complex combination of various lines-of-sight within the galaxy. We explore the manifestation of this effect to study the reliability of line ratios to infer galaxy properties, in particular the Balmer decrement and the BPT diagram. We find the Balmer decrement predicted Balmer line attenuation to be higher (factor of $1$ to $\gtrsim10$) than expected from commonly used attenuation curves. The observed BPT line ratios deviate from their intrinsic values (median difference of 0.08 (0.02) and standard deviation of 0.2 (0.05) for log$_{10}$([N${\rm {\small II}}$]$λ6585/$H$_α$) (log$_{10}$([O${\rm {\small III}}$]$λ5008/$H$_β$)). Finally, we explore the variation in observed properties (UV attenuation, UV slope and Balmer decrement) with viewing angle, finding average differences of $\sim0.3$ magnitudes in the UV attenuation. △ Less

Submitted 11 December, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

Comments: 17 pages, 13 figures. Published MNRAS. Data available at https://flaresimulations.github.io

Computational Orbital Mechanics of Marble Motion on a 3D Printed Surface -- 1. Formal Basis

Authors: Pooja Bhambhu, Preety, Paridhi Goel, Chinkey, Manisha Siwach, Ananya Kumari, Sudarshana, Sanjana Yadav, Shikha Yadav, Bharti, Poonam, Anshumali, Athira Vijayan, Divakar Pathak

Abstract: Simulating curvature due to gravity through warped surfaces is a common visualization aid in Physics education. We reprise a recent experiment exploring orbital trajectories on a precise 3D-printed surface to mimic Newtonian gravity, and elevate this analogy past the status of a mere visualization tool. We present a general analysis approach through which this straightforward experiment can be use… ▽ More Simulating curvature due to gravity through warped surfaces is a common visualization aid in Physics education. We reprise a recent experiment exploring orbital trajectories on a precise 3D-printed surface to mimic Newtonian gravity, and elevate this analogy past the status of a mere visualization tool. We present a general analysis approach through which this straightforward experiment can be used to create a reasonably advanced computational orbital mechanics lab at the undergraduate level, creating a convenient hands-on, computational pathway to various non-trivial nuances in this discipline, such as the mean, eccentric, and true anomalies and their computation, Laplace-Runge-Lenz vector conservation, characterization of general orbits, and the extraction of orbital parameters. We show that while the motion of a marble on such a surface does not truly represent a orbital trajectory under Newtonian gravity in a strict theoretical sense, but through a proposed projection procedure, the experimentally recorded trajectories closely resemble the Kepler orbits and approximately respect the known conservation laws for orbital motion. The latter fact is demonstrated through multiple experimentally-recorded elliptical trajectories with wide-ranging eccentricities and semi-major axes. In this first part of this two-part sequence, we lay down the formal basis of this exposition, describing the experiment, its calibration, critical assessment of the results, and the computational procedures for the transformation of raw experimental data into a form useful for orbital analysis. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: Article 1 of a two-part sequence. Article 2 (applications) to be released shortly. Comments are most welcome !

An Atlas of Color-selected Quiescent Galaxies at $z>3$ in Public $JWST$ Fields

Authors: Francesco Valentino, Gabriel Brammer, Katriona M. L. Gould, Vasily Kokorev, Seiji Fujimoto, Christian Kragh Jespersen, Aswin P. Vijayan, John R. Weaver, Kei Ito, Masayuki Tanaka, Olivier Ilbert, Georgios E. Magdis, Katherine E. Whitaker, Andreas L. Faisst, Anna Gallazzi, Steven Gillman, Clara Gimenez-Arteaga, Carlos Gomez-Guijarro, Mariko Kubo, Kasper E. Heintz, Michaela Hirschmann, Pascal Oesch, Masato Onodera, Francesca Rizzo, Minju Lee , et al. (2 additional authors not shown)

Abstract: We present the results of a systematic search for candidate quiescent galaxies in the distant Universe in eleven $JWST$ fields with publicly available observations collected during the first three months of operations and covering an effective sky area of $\sim145$ arcmin$^2$. We homogeneously reduce the new $JWST$ data and combine them with existing observations from the… ▽ More We present the results of a systematic search for candidate quiescent galaxies in the distant Universe in eleven $JWST$ fields with publicly available observations collected during the first three months of operations and covering an effective sky area of $\sim145$ arcmin$^2$. We homogeneously reduce the new $JWST$ data and combine them with existing observations from the $Hubble\,Space\,Telescope$. We select a robust sample of $\sim80$ candidate quiescent and quenching galaxies at $3 < z < 5$ using two methods: (1) based on their rest-frame $UVJ$ colors, and (2) a novel quantitative approach based on Gaussian Mixture Modeling of the $NUV-U$, $U-V$, and $V-J$ rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive ($M_\star\geq 10^{10.6} M_\odot$) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of $2-3$, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at $z>3$, as it could be expected for highly clustered massive systems. Importantly, $JWST$ enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 15 pages, 6 Figures + Appendix. Accepted for publication in ApJ on Feb, 9. Data release: - Reduced HST+JWST mosaics + photometric catalogs and Eazy-py modeling: https://erda.ku.dk/archives/7166d013c1ca1371aac3c57b9e73190d/published-archive.html - Supplementary material and tables: https://zenodo.org/record/7614908#.Y-4ZruzMLmE - MAST: https://doi.org/10.17909/g3nt-a370 - See also Gould et al. 2023

First Light And Reionisation Epoch Simulations (FLARES) XI: [OIII] emitting galaxies at $5<z<10$

Authors: Stephen M. Wilkins, Christopher C. Lovell, Aswin P. Vijayan, Dimitrios Irodotou, Nathan J. Adams, William J. Roper, Joseph Caruana, Jorryt Matthee, Louise T. C. Seeyave, Christopher J. Conselice, Pablo G. Pérez-González, Jack C. Turner, James M. S. Donnellan

Abstract: JWST has now made it possible to probe the rest-frame optical line emission of high-redshift galaxies extending to z~9, and potentially beyond. To aid in the interpretation of these emerging constraints, in this work we explore predictions for [OIII] emission in high-redshift galaxies using the First Light and Reionisation Epoch Simulations (FLARES). We produce predictions for the [OIII] luminosit… ▽ More JWST has now made it possible to probe the rest-frame optical line emission of high-redshift galaxies extending to z~9, and potentially beyond. To aid in the interpretation of these emerging constraints, in this work we explore predictions for [OIII] emission in high-redshift galaxies using the First Light and Reionisation Epoch Simulations (FLARES). We produce predictions for the [OIII] luminosity function, its correlation with the UV luminosity, and the distribution of equivalent widths (EWs). We also explore how the [OIII] EW correlates with physical properties including specific star formation rate, metallicity, and dust attenuation. Our predictions are largely consistent with recent observational constraints on the luminosity function, average equivalent widths, and line ratios. However, they fail to reproduce the observed tail of high-EW sources and the number density of extreme line emitters. Possibilities to explain these discrepancies include an additional source of ionising photons and/or greater stochasticity in star formation in the model or photometric scatter and/or bias in the observations. With JWST now rapidly building larger samples and a wider range of emission lines the answer to this remaining discrepancy should be available imminently. △ Less

Submitted 14 April, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 15 pages, accepted for publication in MNRAS, minor changes from original version

First Light and Reionisation Epoch Simulations (FLARES) X: Environmental Galaxy Bias and Survey Variance at High Redshift

Authors: Peter A. Thomas, Christopher C. Lovell, Maxwell G. A. Maltz, Aswin P. Vijayan, Stephen M. Wilkins, Dimitrios Irodotou, William J. Roper, Louise Seeyave

Abstract: Upcoming deep galaxy surveys with JWST will probe galaxy evolution during the epoch of reionisation (EoR, $5\leq z\leq10$) over relatively compact areas (e.g. $\sim$ 300\,arcmin$^2$ for the JADES GTO survey). It is therefore imperative that we understand the degree of survey variance, to evaluate how representative the galaxy populations in these studies will be. We use the First Light And Reionis… ▽ More Upcoming deep galaxy surveys with JWST will probe galaxy evolution during the epoch of reionisation (EoR, $5\leq z\leq10$) over relatively compact areas (e.g. $\sim$ 300\,arcmin$^2$ for the JADES GTO survey). It is therefore imperative that we understand the degree of survey variance, to evaluate how representative the galaxy populations in these studies will be. We use the First Light And Reionisation Epoch Simulations (FLARES) to measure the galaxy bias of various tracers over an unprecedentedly large range in overdensity for a hydrodynamic simulation, and use these relations to assess the impact of bias and clustering on survey variance in the EoR. Star formation is highly biased relative to the underlying dark matter distribution, with the mean ratio of the stellar to dark matter density varying by a factor of 100 between regions of low and high matter overdensity (smoothed on a scale of 14$\,h^{-1}$cMpc). This is reflected in the galaxy distribution --the most massive galaxies are found solely in regions of high overdensity. As a consequence of the above, galaxies in the EoR are highly clustered, which can lead to large variance in survey number counts. For mean number counts $N\lesssim 100$ (1000), in a unit redshift slice of angular area 300\,arcmin$^2$ (1.4\,deg$^2$), the 2-sigma range in $N$ is roughly a factor of four (two). We present relations between the expected variance and survey area for different survey geometries; these relations will be of use to observers wishing to understand the impact of survey variance on their results. △ Less

Submitted 14 June, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: Accepted version. Thanks to the referee contains more extensive description of testing, comparison to online Cosmic Variance Calculator, and greater range of survey areas than the original submission

FLARES IX: The Physical Mechanisms Driving Compact Galaxy Formation and Evolution

Authors: William J. Roper, Christopher C. Lovell, Aswin P. Vijayan, Dimitrios Irodotou, Jussi K. Kuusisto, Jasleen Matharu, Louise T. C. Seeyave, Peter A. Thomas, Stephen M. Wilkins

Abstract: In the FLARES (First Light And Reionisation Epoch Simulations) suite of hydrodynamical simulations, we find the high redshift ($z>5$) intrinsic size-luminosity relation is, surprisingly, negatively sloped. However, after including the effects of dust attenuation we find a positively sloped UV observed size-luminosity relation in good agreement with other simulated and observational studies. In thi… ▽ More In the FLARES (First Light And Reionisation Epoch Simulations) suite of hydrodynamical simulations, we find the high redshift ($z>5$) intrinsic size-luminosity relation is, surprisingly, negatively sloped. However, after including the effects of dust attenuation we find a positively sloped UV observed size-luminosity relation in good agreement with other simulated and observational studies. In this work, we extend this analysis to probe the underlying physical mechanisms driving the formation and evolution of the compact galaxies driving the negative size-mass/size-luminosity relation. We find the majority of compact galaxies ($R_{1/2, \star}< 1 \mathrm{pkpc}$), which drive the negative slope of the size-mass relation, have transitioned from extended to compact sizes via efficient centralised cooling, resulting in high specific star formation rates in their cores. These compact stellar systems are enshrouded by non-star forming gas distributions as much as $100\times$ larger than their stellar counterparts. By comparing with galaxies from the EAGLE simulation suite, we find that these extended gas distributions `turn on' and begin to form stars between $z=5$ and $z=0$ leading to increasing sizes, and thus the evolution of the size-mass relation from a negative to a positive slope. This explicitly demonstrates the process of inside-out galaxy formation in which compact bulges form earlier than the surrounding discs. △ Less

Submitted 6 February, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Massive galaxy formation caught in action at z~5 with JWST

Authors: Shuowen Jin, Nikolaj B. Sillassen, Georgios E. Magdis, Aswin P. Vijayan, Gabriel B. Brammer, Vasily Kokorev, John R. Weaver, Raphael Gobat, Clara Giménez-Arteaga, Francesco Valentino, Malte Brinch, Carlos Gómez-Guijarro, Marko Shuntov, Sune Toft, Thomas R. Greve, David Blanquez Sese

Abstract: We report the discovery of a compact group of galaxies, CGG-z5, at z~5.2 in the EGS field covered by the JWST/CEERS survey. CGG-z5 was selected as the highest overdensity of galaxies at z>2 in recent JWST public surveys and it consists of six candidate members lying within a projected area of $1.5"\times3"$ (10$\times$20~kpc$^2$). All group members are HST/F435W and HST/F606W dropouts while secure… ▽ More We report the discovery of a compact group of galaxies, CGG-z5, at z~5.2 in the EGS field covered by the JWST/CEERS survey. CGG-z5 was selected as the highest overdensity of galaxies at z>2 in recent JWST public surveys and it consists of six candidate members lying within a projected area of $1.5"\times3"$ (10$\times$20~kpc$^2$). All group members are HST/F435W and HST/F606W dropouts while securely detected in the JWST/NIRCam bands, yielding a narrow range of robust photometric redshifts $5.0<z<5.3$. The most massive galaxy in the group has a stellar mass log$(M_{*}/M_{\odot})\approx9.8$, while the rest are low-mass satellites (log$(M_{*}/M_{\odot})\approx8.4-9.2$). While several group members were already detected in the HST and IRAC bands, the low stellar masses and the compactness of the structure required the sensitivity and resolution of JWST for its identification. To assess the nature and evolutionary path of CGG-z5, we searched for similar compact structures in the \textsc{Eagle} simulations and followed their evolution with time. We find that all the identified structures merge into a single galaxy by z=3 and form a massive galaxy (log$(M_{*}/M_{\odot})>11$) at z~1. This implies that CGG-z5 could be a "proto-massive galaxy" captured during a short-lived phase of massive galaxy formation. △ Less

Submitted 9 November, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: Published in A&A Letter

Journal ref: A&A 670, L11 (2023)

The gas and stellar content of a metal-poor galaxy at $z=8.496$ as revealed by JWST and ALMA

Authors: K. E. Heintz, C. Giménez-Arteaga, S. Fujimoto, G. Brammer, D. Espada, S. Gillman, J. González-López, T. R. Greve, Y. Harikane, B. Hatsukade, K. K. Knudsen, A. M. Koekemoer, K. Kohno, V. Kokorev, M. M. Lee, G. E. Magdis, E. J. Nelson, F. Rizzo, R. L. Sanders, D. Schaerer, A. E. Shapley, V. B. Strait, S. Toft, F. Valentino, A. van der Wel , et al. (5 additional authors not shown)

Abstract: We present a joint analysis of the galaxy S04590 at $z=8.496$ based on NIRSpec, NIRCam, and NIRISS observations obtained through as part of Early Release Observations programme of the James Webb Space Telescope (JWST) and the far-infrared [CII]-$158μ$m emission line detected by dedicated Atacama Large Millimeter/submillimeter Array (ALMA) observations. We determine the physical properties of S0459… ▽ More We present a joint analysis of the galaxy S04590 at $z=8.496$ based on NIRSpec, NIRCam, and NIRISS observations obtained through as part of Early Release Observations programme of the James Webb Space Telescope (JWST) and the far-infrared [CII]-$158μ$m emission line detected by dedicated Atacama Large Millimeter/submillimeter Array (ALMA) observations. We determine the physical properties of S04590 from modelling of the spectral energy distribution (SED) and through the redshifted optical nebular emission lines detected with JWST/NIRSpec. The best-fit SED model reveals a low-mass ($M_\star = 10^{7.2}-10^{8}\,M_{\odot}$) galaxy with a low oxygen abundance of $12+\log{\rm (O/H)} = 7.16^{+0.10}_{-0.12}$ derived from the strong nebular and auroral emission lines. Assuming that [CII] effectively traces the interstellar medium (ISM), we estimate the total gas mass of the galaxy to be $M_{\rm gas} = (8.0\pm 4.0)\times 10^{8}\,M_\odot$ based on the luminosity and spatial extent of [CII]. This yields an exceptionally high gas fraction, $f_{\rm gas} = M_{\rm gas}/(M_{\rm gas} + M_\star) \gtrsim 90\%$, though still consistent within the range expected for its low metallicity. We further derive the metal mass of the galaxy based on the gas mass and gas-phase metallicity, which we find to be consistent with the expected metal production from Type II supernovae. Finally, we make the first constraints on the dust-to-gas (DTG) and dust-to-metals (DTM) ratios of galaxies in the epoch of reionization at $z\gtrsim 6$, showing overall low mass ratios of logDGT $<-3.8$ and logDTM $<-0.5$, though consistent with local scaling relations and in particular the local metal-poor galaxy I Zwicky 18. Our analysis highlights the synergy between ALMA and JWST in characterizing the gas, metal, and stellar content of the first generation of galaxies. △ Less

Submitted 6 February, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: Accepted for publication in ApJL

Dilution of chemical enrichment in galaxies 600 Myr after the Big Bang

Authors: Kasper E. Heintz, Gabriel B. Brammer, Clara Giménez-Arteaga, Victoria B. Strait, Claudia del P. Lagos, Aswin P. Vijayan, Jorryt Matthee, Darach Watson, Charlotte A. Mason, Anne Hutter, Sune Toft, Johan P. U. Fynbo, Pascal A. Oesch

Abstract: Galaxies throughout the last 12 Gyr of cosmic time follow a single, universal relation that connects their star-formation rates (SFRs), stellar masses ($M_\star$) and chemical abundances. Deviation from these fundamental scaling relations would imply a drastic change in the processes that regulate galaxy evolution. Observations have hinted at the possibility that this relation may be broken in the… ▽ More Galaxies throughout the last 12 Gyr of cosmic time follow a single, universal relation that connects their star-formation rates (SFRs), stellar masses ($M_\star$) and chemical abundances. Deviation from these fundamental scaling relations would imply a drastic change in the processes that regulate galaxy evolution. Observations have hinted at the possibility that this relation may be broken in the very early universe. However, until recently, chemical abundances of galaxies could be only measured reliably as far back as redshift $z=3.3$. With JWST, we can now characterize the SFR, $M_\star$, and chemical abundance of galaxies during the first few hundred million years after the Big Bang, at redshifts $z=7-10$. Here we show that galaxies at this epoch follow unique SFR-$M_\star$--main-sequence and mass-metallicity scaling relations, but their chemical abundance is a factor of three lower than expected from the fundamental-metallicity relation of later galaxies. These findings suggest that galaxies at this time are still intimately connected with the intergalactic medium and subject to continuous infall of pristine gas which effectively dilutes their metal abundances. △ Less

Submitted 23 February, 2023; v1 submitted 6 December, 2022; originally announced December 2022.

Comments: Resubmitted version, taking reviewers comments into account

First Light And Reionisation Epoch Simulations (FLARES) VIII. The Emergence of Passive Galaxies at $z \geqslant 5$

Authors: Christopher C. Lovell, Will Roper, Aswin P. Vijayan, Louise Seeyave, Dimitrios Irodotou, Stephen M. Wilkins, Christopher J. Conselice, Flaminia Fortuni, Jussi K. Kuusisto, Emiliano Merlin, Paola Santini, Peter Thomas

Abstract: Passive galaxies are ubiquitous in the local universe, and various physical channels have been proposed that lead to this passivity. To date, robust passive galaxy candidates have been detected up to $z \leqslant 5$, but it is still unknown if they exist at higher redshifts, what their relative abundances are, and what causes them to stop forming stars. We present predictions from the First Light… ▽ More Passive galaxies are ubiquitous in the local universe, and various physical channels have been proposed that lead to this passivity. To date, robust passive galaxy candidates have been detected up to $z \leqslant 5$, but it is still unknown if they exist at higher redshifts, what their relative abundances are, and what causes them to stop forming stars. We present predictions from the First Light And Reionisation Epoch Simulations (FLARES), a series of zoom simulations of a range of overdensities using the EAGLE code. Passive galaxies occur naturally in the EAGLE model at high redshift, and are in good agreement with number density estimates from HST and early JWST results at $3 \leqslant z \leqslant 5$. Due to the unique FLARES approach, we extend these predictions to higher redshifts, finding passive galaxy populations up to $z \sim 8$. Feedback from supermassive black holes is the main driver of passivity, leading to reduced gas fractions and star forming gas reservoirs. We find that passive galaxies at $z \geqslant 5$ are not identified in the typical UVJ selection space due to their still relatively young stellar populations, and present new rest--frame selection regions. We also present NIRCam and MIRI fluxes, and find that significant numbers of passive galaxies at $z \geqslant 5$ should be detectable in upcoming wide surveys with JWST. Finally, we present JWST colour distributions, with new selection regions in the observer--frame for identifying these early passive populations. △ Less

Submitted 21 August, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 21 pages, 20 figures. Accepted to MNRAS

X-ray Emission from the Interstellar and Circumgalactic Medium of Elliptical Galaxies based on MACER simulations

Authors: Aditi Vijayan, Bocheng Zhu, Miao Li, Feng Yuan, Luis C. Ho

Abstract: Interstellar (ISM) and circumgalactic mediums (CGM) around galaxies are linked to several physical processes that drive galaxy evolution. For example, the X-ray emission from the CGM gas around ellipticals has been linked to the AGN feedback occurring in the host. Upcoming telescopes, such as HUBS with ~2 eV resolution, can provide us with deep insights about the hot gas properties of such galaxie… ▽ More Interstellar (ISM) and circumgalactic mediums (CGM) around galaxies are linked to several physical processes that drive galaxy evolution. For example, the X-ray emission from the CGM gas around ellipticals has been linked to the AGN feedback occurring in the host. Upcoming telescopes, such as HUBS with ~2 eV resolution, can provide us with deep insights about the hot gas properties of such galaxies thus constrain these processes. In this project, we discuss X-ray emission of the ISM and CGM of elliptical galaxies simulated using MACER code. We generate X-ray emission data from the MACER simulations with various feedback models and produce mock observations for an instrument with high spectral resolution, which is a necessary step of selecting sources for the future observations with planned mission such as HUBS. More importantly, we establish connections between the physics of AGN and stellar feedback with the emission spectra from the ISM and CGM to investigate the possibility of using observations to constrain feedback models. We fit the X-ray spectra from these simulations with standard fitting procedures and compare the retrieved physical properties with their counterparts from the simulations to understand whether the future high-resolution observations can reliably reveal the properties of the gas in the galaxies. △ Less

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

Comments: 10 pages, 13 Figures, accepted in MNRAS

Unveiling the main sequence of galaxies at $z \geq 5$ with the James Webb Space Telescope: predictions from simulations

Authors: Jordan C. J. D'Silva, Claudia D. P. Lagos, Luke J. M. Davies, Christopher C. Lovell, Aswin P. Vijayan

Abstract: We use two independent, galaxy formation simulations, FLARES, a cosmological hydrodynamical simulation, and SHARK, a semi-analytic model, to explore how well the James Webb Space Telescope (JWST) will be able to uncover the existence and parameters of the star-forming main sequence (SFS) at $z=5\to10$, i.e. shape, scatter, normalisation. Using two independent simulations allows us to isolate predi… ▽ More We use two independent, galaxy formation simulations, FLARES, a cosmological hydrodynamical simulation, and SHARK, a semi-analytic model, to explore how well the James Webb Space Telescope (JWST) will be able to uncover the existence and parameters of the star-forming main sequence (SFS) at $z=5\to10$, i.e. shape, scatter, normalisation. Using two independent simulations allows us to isolate predictions (e.g., stellar mass, star formation rate, SFR, luminosity functions) that are robust to or highly dependent on the implementation of the physics of galaxy formation. Both simulations predict that JWST can observe $\ge 70-90\%$ (for SHARK and FLARES respectively) of galaxies up to $z\sim10$ (down to stellar masses of $\approx 10^{8.3}\,\rm M_{\odot}$ and SFRs of $\approx 10^{0.5}\,\rm M_{\odot}\, yr^{-1}$) in modest integration times and given current proposed survey areas (e.g. the Web COSMOS $0.6\,\rm deg^2$) to accurately constrain the parameters of the SFS. Although both simulations predict qualitatively similar distributions of stellar mass and SFR, there are important quantitative differences, such as the abundance of massive, star-forming galaxies, with FLARES predicting a higher abundance than SHARK; the early onset of quenching as a result of black hole growth in FLARES (at $z\approx 8$), not seen in SHARK until much lower redshifts; and the implementation of synthetic photometry, with FLARES predicting more JWST-detected galaxies ($\sim 90\%$) than SHARK ($\sim 70\%$) at $z=10$. JWST observations will distinguish between these models, leading to a significant improvement upon our understanding of the formation of the very first galaxies. △ Less

Submitted 5 October, 2022; v1 submitted 12 August, 2022; originally announced August 2022.

Comments: 21 pages, 13 figures including 1 in appendix, accepted to MNRAS

First Light And Reionisation Epoch Simulations (FLARES) VII: The Star Formation and Metal Enrichment Histories of Galaxies in the early Universe

Authors: Stephen M. Wilkins, Aswin P. Vijayan, Christopher C. Lovell, William J. Roper, Erik Zackrisson, Dimitrios Irodotou, Louise T. C. Seeyave, Jussi K. Kuusisto, Peter A. Thomas, Joseph Caruana, Christopher J. Conselice

Abstract: The star formation and metal enrichment histories of galaxies - at any epoch - constitute one of the key properties of galaxies, and their measurement is a core aim of observational extragalactic astronomy. The lack of deep rest-frame optical coverage at high-redshift has made robust constraints elusive, but this is now changing thanks to the \emph{James Webb Space Telescope (JWST)}. In preparatio… ▽ More The star formation and metal enrichment histories of galaxies - at any epoch - constitute one of the key properties of galaxies, and their measurement is a core aim of observational extragalactic astronomy. The lack of deep rest-frame optical coverage at high-redshift has made robust constraints elusive, but this is now changing thanks to the \emph{James Webb Space Telescope (JWST)}. In preparation for the constraints provided by \emph{JWST} we explore the star formation and metal enrichment histories of galaxies at $z=5-13$ using the First Light And Reionisation Epoch Simulations (FLARES) suite. Built on the EAGLE model, the unique strategy of FLARES allows us to simulate a wide range of stellar masses (and luminosities) and environments. While we predict significant redshift evolution of average ages and specific star formation rates our core result is a mostly flat relationship of age and specific star formation rate with stellar mass. We also find that galaxies in this epoch predominantly have strongly rising star formation histories, albeit with the magnitude dropping with redshift and stellar mass. In terms of chemical enrichment we predict a strong stellar mass - metallicity relation present at $z=10$ and beyond alongside significant $α$-enhancement. Finally, we find no environmental dependence of the relationship between age, specific star formation rate, or metallicity with stellar mass. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: To be submitted to MNRAS. Comments welcome

Seeing sharper and deeper: JWST's first glimpse of the photometric and spectroscopic properties of galaxies in the epoch of reionisation

Authors: James A. A. Trussler, Nathan J. Adams, Christopher J. Conselice, Leonardo Ferreira, Duncan Austin, Rachana Bhatawdekar, Joseph Caruana, Brenda L. Frye, Tom Harvey, Christopher C. Lovell, Massimo Pascale, William J. Roper, Aprajita Verma, Aswin P. Vijayan, Stephen M. Wilkins

Abstract: We analyse the photometric and spectroscopic properties of four galaxies in the epoch of reionisation (EoR) within the SMACS 0723 JWST Early Release Observations field. Given the known spectroscopic redshifts of these sources, we investigated the accuracy with which photometric redshifts can be derived using NIRCam photometry alone, finding that F115W imaging is essential to distinguish between z~… ▽ More We analyse the photometric and spectroscopic properties of four galaxies in the epoch of reionisation (EoR) within the SMACS 0723 JWST Early Release Observations field. Given the known spectroscopic redshifts of these sources, we investigated the accuracy with which photometric redshifts can be derived using NIRCam photometry alone, finding that F115W imaging is essential to distinguish between z~8 galaxies with high equivalent width (EW) [O III] λ5007 emission and z~10 Balmer break galaxies. We find that all four sources exhibit strong (\geq 0.6 mag) F356W-F444W colours, which sit at the extreme end of theoretical predictions from numerical simulations. We find that these galaxies deviate (by roughly 0.5 dex) from the local correlation between [O III] λ5007/Hβand [Ne III] λ3869/[O II], which is consistent with the predictions from simulations of high-redshift galaxies having elevated line excitation ratios. We measure the [O III] λ5007 rest-frame equivalent widths both directly from the spectroscopy, and indirectly as inferred from the strong F356W-F444W colours, finding large [O III] λ5007 EWs of 225-1740 Å. The [O III] λ5007 and HβEWs are consistent with those seen in extreme, intensely star-forming dwarf galaxies in the local Universe. Our structural analysis indicates that these galaxies are resolved, exhibiting irregular shapes with bright clumps. In line with the predictions from the FLARES hydrodynamic simulations, such intense star formation and extreme nebular conditions are likely the norm, rather than the exception, in the EoR. △ Less

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

Comments: 17 pages, 9 figures. Now accepted and published in MNRAS

Journal ref: MNRAS, 2023, 523, 3423

Discovery and properties of ultra-high redshift galaxies ($9<z<12$) in the JWST ERO SMACS 0723 Field

Authors: N. J. Adams, C. J. Conselice, L. Ferreira, D. Austin, J. Trussler, I. Juodžbalis, S. M. Wilkins, J. Caruana, P. Dayal, A. Verma, A. P. Vijayan

Abstract: We present a reduction and analysis of the \textit{James Webb Space Telescope} (JWST) SMACS~0723 field using new post-launch calibrations to conduct a search for ultra-high-redshift galaxies ($z > 9$) present within the Epoch of Reionisation. We conduct this search by modelling photometric redshifts in several ways for all sources and by applying conservative magnitude cuts ($m_{\rm F200W} < 28$)… ▽ More We present a reduction and analysis of the \textit{James Webb Space Telescope} (JWST) SMACS~0723 field using new post-launch calibrations to conduct a search for ultra-high-redshift galaxies ($z > 9$) present within the Epoch of Reionisation. We conduct this search by modelling photometric redshifts in several ways for all sources and by applying conservative magnitude cuts ($m_{\rm F200W} < 28$) to identify strong Lyman breaks greater than 1 magnitude. We find four $z > 9$ candidate galaxies which have not previously been identified, with one object at $z = 11.5$, and another which is possibly a close pair of galaxies. We measure redshifts for candidate galaxies from other studies and find the recovery rate to be only 23 per cent, with many being assigned lower redshift, dusty solutions in our work. Most of our $z > 9$ sample show evidence for Balmer-breaks, or extreme emission lines from H$β$ and [OIII], demonstrating that the stellar populations could be advanced in age or very young depending on the cause of the F444W excess. We discuss the resolved structures of these early galaxies and find that the Sérsic indices reveal a mixture of light concentration levels, but that the sizes of all our systems are exceptionally small ($< 0.5$~kpc). These systems have stellar masses M$_{*} \sim 10^{9.0}$ M$_{\odot}$, with our $z \sim 11.5$ candidate a dwarf galaxy with a stellar mass M$_{*} \sim 10^{7.8}$ -- $10^{8.2}$ M$_{\odot}$. These candidate ultra high-redshift galaxies are excellent targets for future NIRSpec observations aimed to better understand their physical nature. △ Less

Submitted 3 January, 2023; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: Updated to accepted mnras manuscript, 12 pages, 8 figures, 6 tables, added section 3.3 and Table 5

First Light And Reionisation Epoch Simulations (FLARES) VI: The colour evolution of galaxies $z=5-15$

Authors: Stephen M. Wilkins, Aswin P. Vijayan, Christopher C. Lovell, William J. Roper, Dimitrios Irodotou, Joseph Caruana, Louise T. C. Seeyave, Jussi K. Kuusisto, Peter A. Thomas

Abstract: With its exquisite sensitivity, wavelength coverage, and spatial and spectral resolution, the James Webb Space Telescope is poised to revolutionise our view of the distant, high-redshift ($z>5$) Universe. While Webb's spectroscopic observations will be transformative for the field, photometric observations play a key role in identifying distant objects and providing more comprehensive samples than… ▽ More With its exquisite sensitivity, wavelength coverage, and spatial and spectral resolution, the James Webb Space Telescope is poised to revolutionise our view of the distant, high-redshift ($z>5$) Universe. While Webb's spectroscopic observations will be transformative for the field, photometric observations play a key role in identifying distant objects and providing more comprehensive samples than accessible to spectroscopy alone. In addition to identifying objects, photometric observations can also be used to infer physical properties and thus be used to constrain galaxy formation models. However, inferred physical properties from broadband photometric observations, particularly in the absence of spectroscopic redshifts, often have large uncertainties. With the development of new tools for forward modelling simulations it is now routinely possible to predict observational quantities, enabling a direct comparison with observations. With this in mind, in this work, we make predictions for the colour evolution of galaxies at $z=5-15$ using the FLARES: First Light And Reionisation Epoch Simulations cosmological hydrodynamical simulation suite. We predict a complex evolution, driven predominantly by strong nebular line emission passing through individual bands. These predictions are in good agreement with existing constraints from Hubble and Spitzer as well as some of the first results from Webb. We also contrast our predictions with other models in the literature: while the general trends are similar we find key differences, particularly in the strength of features associated with strong nebular line emission. This suggests photometric observations alone should provide useful discriminating power between different models. △ Less

Submitted 6 September, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 10 pages, 9 figures, accepted for publication in MNRAS

The BlueTides Mock Image Catalogue: Simulated observations of high-redshift galaxies and predictions for JWST imaging surveys

Authors: Madeline A. Marshall, Katelyn Watts, Stephen Wilkins, Tiziana Di Matteo, Jussi K. Kuusisto, William J. Roper, Aswin P. Vijayan, Yueying Ni, Yu Feng, Rupert A. C. Croft

Abstract: We present a mock image catalogue of ~100,000 MUV=-22.5 to -19.6 mag galaxies at z=7-12 from the BlueTides cosmological simulation. We create mock images of each galaxy with the James Webb (JWST), Hubble, Roman, and Euclid Space Telescopes, as well as Subaru, and VISTA, with a range of near- and mid-infrared filters. We perform photometry on the mock images to estimate the success of these instrum… ▽ More We present a mock image catalogue of ~100,000 MUV=-22.5 to -19.6 mag galaxies at z=7-12 from the BlueTides cosmological simulation. We create mock images of each galaxy with the James Webb (JWST), Hubble, Roman, and Euclid Space Telescopes, as well as Subaru, and VISTA, with a range of near- and mid-infrared filters. We perform photometry on the mock images to estimate the success of these instruments for detecting high-z galaxies. We predict that JWST will have unprecedented power in detecting high-z galaxies, with a 95% completeness limit at least 2.5 magnitudes fainter than VISTA and Subaru, 1.1 magnitudes fainter than Hubble, and 0.9 magnitudes fainter than Roman, for the same wavelength and exposure time. Focusing on JWST, we consider a range of exposure times and filters, and find that the NIRCam F356W and F277W filters will detect the faintest galaxies, with 95% completeness at m=27.4 mag in 10ks exposures. We also predict the number of high-z galaxies that will be discovered by upcoming JWST imaging surveys. We predict that the COSMOS-Web survey will detect ~1000 MUV<-20.1 mag galaxies at 6.5<z<7.5, by virtue of its large survey area. JADES-Medium will detect almost 100% of MUV<-20 mag galaxies at z<8.5 due to its significant depth, however with its smaller survey area it will detect only ~100 of these galaxies at 6.5<z<7.5. Cosmic variance results in a large range in the number of predicted galaxies each survey will detect, which is more evident in smaller surveys such as CEERS and the PEARLS NEP and GOODS-S fields. △ Less

Submitted 25 August, 2022; v1 submitted 17 June, 2022; originally announced June 2022.

Comments: Published in MNRAS. The BlueTides Mock Image Catalogue is publicly available as a High Level Science Product (HLSP) via the Mikulski Archive for Space Telescopes (MAST), at https://doi.org/10.17909/er09-4527. Example codes for using the catalogue data are available at https://github.com/madelinemarshall/BlueTidesMockImageCatalogue

CellTypeGraph: A New Geometric Computer Vision Benchmark

Authors: Lorenzo Cerrone, Athul Vijayan, Tejasvinee Mody, Kay Schneitz, Fred A. Hamprecht

Abstract: Classifying all cells in an organ is a relevant and difficult problem from plant developmental biology. We here abstract the problem into a new benchmark for node classification in a geo-referenced graph. Solving it requires learning the spatial layout of the organ including symmetries. To allow the convenient testing of new geometrical learning methods, the benchmark of Arabidopsis thaliana ovule… ▽ More Classifying all cells in an organ is a relevant and difficult problem from plant developmental biology. We here abstract the problem into a new benchmark for node classification in a geo-referenced graph. Solving it requires learning the spatial layout of the organ including symmetries. To allow the convenient testing of new geometrical learning methods, the benchmark of Arabidopsis thaliana ovules is made available as a PyTorch data loader, along with a large number of precomputed features. Finally, we benchmark eight recent graph neural network architectures, finding that DeeperGCN currently works best on this problem. △ Less

Submitted 17 May, 2022; originally announced May 2022.

First Light And Reionisation Epoch Simulations (FLARES) V: The redshift frontier

Authors: Stephen M. Wilkins, Aswin P. Vijayan, Christopher C. Lovell, William J. Roper, Dimitrios Irodotou, Joseph Caruana, Louise T. C. Seeyave, Jussi K. Kuusisto, Peter A. Thomas, Shedeur A. K. Parris

Abstract: The James Webb Space Telescope (JWST) is set to transform many areas of astronomy, one of the most exciting is the expansion of the redshift frontier to $z>10$. In its first year alone JWST should discover hundreds of galaxies, dwarfing the handful currently known. To prepare for these powerful observational constraints, we use the First Light And Reionisation Epoch (FLARES) simulations to predict… ▽ More The James Webb Space Telescope (JWST) is set to transform many areas of astronomy, one of the most exciting is the expansion of the redshift frontier to $z>10$. In its first year alone JWST should discover hundreds of galaxies, dwarfing the handful currently known. To prepare for these powerful observational constraints, we use the First Light And Reionisation Epoch (FLARES) simulations to predict the physical and observational properties of the $z>10$ population of galaxies accessible to JWST. This is the first time such predictions have been made using a hydrodynamical model validated at low redshift. Our predictions at $z=10$ are broadly in agreement with current observational constraints on the far-UV luminosity function and UV continuum slope $β$, though the observational uncertainties are large. We note tension with recent constraints $z\sim 13$ from Harikane et al. 2022 - compared to these constraints, FLARES predicts objects with the same space density should have an order of magnitude lower luminosity, though this is mitigated slightly if dust attenuation is negligible in these systems. Our predictions suggest that in JWST's first cycle alone, around $600$ galaxies should be identified at $z>10$, with the first small samples available at $z>13$. △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: Submitted to MNRAS 20/04/2022

First Light And Reionisation Epoch Simulations (FLARES) IV: The size evolution of galaxies at $z\geq5$

Authors: William J. Roper, Christopher C. Lovell, Aswin P. Vijayan, Madeline A. Marshall, Dimitrios Irodotou, Jussi K. Kuusisto, Peter A. Thomas, Stephen M. Wilkins

Abstract: We present the intrinsic and observed sizes of galaxies at $z\geq5$ in the First Light And Reionisation Epoch Simulations (FLARES). We employ the large effective volume of FLARES to produce a sizeable sample of high redshift galaxies with intrinsic and observed luminosities and half light radii in a range of rest frame UV and visual photometric bands. This sample contains a significant number of i… ▽ More We present the intrinsic and observed sizes of galaxies at $z\geq5$ in the First Light And Reionisation Epoch Simulations (FLARES). We employ the large effective volume of FLARES to produce a sizeable sample of high redshift galaxies with intrinsic and observed luminosities and half light radii in a range of rest frame UV and visual photometric bands. This sample contains a significant number of intrinsically ultra-compact galaxies in the far-UV (1500 angstrom), leading to a negative intrinsic far-UV size-luminosity relation. However, after the inclusion of the effects of dust these same compact galaxies exhibit observed sizes that are as much as 50 times larger than those measured from the intrinsic emission, and broadly agree with a range of observational samples. This increase in size is driven by the concentration of dust in the core of galaxies, heavily attenuating the intrinsically brightest regions. At fixed luminosity we find a galaxy size redshift evolution with a slope of $m=1.21-1.87$ depending on the luminosity sample in question, and we demonstrate the wavelength dependence of the size-luminosity relation which will soon be probed by the Webb Space Telescope. △ Less

Submitted 23 March, 2022; originally announced March 2022.

The Impact of Dust on the Sizes of Galaxies in the Epoch of Reionization

Authors: Madeline A. Marshall, Stephen Wilkins, Tiziana Di Matteo, William J. Roper, Aswin P. Vijayan, Yueying Ni, Yu Feng, Rupert A. C. Croft

Abstract: We study the sizes of galaxies in the Epoch of Reionization using a sample of ~100,000 galaxies from the BlueTides cosmological hydrodynamical simulation from z=7 to 11. We measure the galaxy sizes from stellar mass and luminosity maps, defining the effective radius as the minimum radius which could enclose the pixels containing 50% of the total mass/light in the image. We find an inverse relation… ▽ More We study the sizes of galaxies in the Epoch of Reionization using a sample of ~100,000 galaxies from the BlueTides cosmological hydrodynamical simulation from z=7 to 11. We measure the galaxy sizes from stellar mass and luminosity maps, defining the effective radius as the minimum radius which could enclose the pixels containing 50% of the total mass/light in the image. We find an inverse relationship between stellar mass and effective half-mass radius, suggesting that the most massive galaxies are more compact and dense than lower mass galaxies, which have flatter mass distributions. We find a mildly negative relation between intrinsic far-ultraviolet luminosity and size, while we find a positive size-luminosity relation when measured from dust-attenuated images. This suggests that dust is the predominant cause of the observed positive size-luminosity relation, with dust preferentially attenuating bright sight lines resulting in a flatter emission profile and thus larger measured effective radii. We study the size-luminosity relation across the rest-frame ultraviolet and optical, and find that the slope decreases at longer wavelengths; this is a consequence of the relation being caused by dust, which produces less attenuation at longer wavelengths. We find that the far-ultraviolet size-luminosity relation shows mild evolution from z=7 to 11, and galaxy size evolves with redshift as $R\propto(1+z)^{-m}$, where $m=0.662\pm0.009$. Finally, we investigate the sizes of z=7 quasar host galaxies, and find that while the intrinsic sizes of quasar hosts are small relative to the overall galaxy sample, they have comparable sizes when measured from dust-attenuated images. △ Less

Submitted 9 February, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 12 pages, 10 figures. Accepted for publication in MNRAS