Showing 1–50 of 364 results for author: Davies, J

Euclid. III. The NISP Instrument

Authors: Euclid Collaboration, K. Jahnke, W. Gillard, M. Schirmer, A. Ealet, T. Maciaszek, E. Prieto, R. Barbier, C. Bonoli, L. Corcione, S. Dusini, F. Grupp, F. Hormuth, S. Ligori, L. Martin, G. Morgante, C. Padilla, R. Toledo-Moreo, M. Trifoglio, L. Valenziano, R. Bender, F. J. Castander, B. Garilli, P. B. Lilje, H. -W. Rix , et al. (412 additional authors not shown)

Abstract: The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm wavelength range. In this reference article we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the proc… ▽ More The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm wavelength range. In this reference article we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the processes needed to understand how NISP operates and is calibrated, and its technical potentials and limitations. Links to articles providing more details and technical background are included. NISP's 16 HAWAII-2RG (H2RG) detectors with a plate scale of 0.3" pix^-1 deliver a field-of-view of 0.57deg^2. In photo mode, NISP reaches a limiting magnitude of ~24.5AB mag in three photometric exposures of about 100s exposure time, for point sources and with a signal-to-noise ratio (SNR) of 5. For spectroscopy, NISP's point-source sensitivity is a SNR = 3.5 detection of an emission line with flux ~2x10^-16erg/s/cm^2 integrated over two resolution elements of 13.4A, in 3x560s grism exposures at 1.6 mu (redshifted Ha). Our calibration includes on-ground and in-flight characterisation and monitoring of detector baseline, dark current, non-linearity, and sensitivity, to guarantee a relative photometric accuracy of better than 1.5%, and relative spectrophotometry to better than 0.7%. The wavelength calibration must be better than 5A. NISP is the state-of-the-art instrument in the NIR for all science beyond small areas available from HST and JWST - and an enormous advance due to its combination of field size and high throughput of telescope and instrument. During Euclid's 6-year survey covering 14000 deg^2 of extragalactic sky, NISP will be the backbone for determining distances of more than a billion galaxies. Its NIR data will become a rich reference imaging and spectroscopy data set for the coming decades. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue 'Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Modelling the galaxy radio continuum from star formation and active galactic nuclei in the Shark semi-analytic model

Authors: Samuel P. Hansen, Claudia D. P. Lagos, Matteo Bonato, Robin H. W. Cook, Luke J. M. Davies, Ivan Delvecchio, Scott A. Tompkins

Abstract: We present a model of radio continuum emission associated with star formation (SF) and active galactic nuclei (AGN) implemented in the Shark semi-analytic model of galaxy formation. SF emission includes free-free and synchrotron emission, which depend on the free-electron density and the rate of core-collapse supernovae with a minor contribution from supernova remnants, respectively. AGN emission… ▽ More We present a model of radio continuum emission associated with star formation (SF) and active galactic nuclei (AGN) implemented in the Shark semi-analytic model of galaxy formation. SF emission includes free-free and synchrotron emission, which depend on the free-electron density and the rate of core-collapse supernovae with a minor contribution from supernova remnants, respectively. AGN emission is modelled based on the jet production rate, which depends on the black hole mass, accretion rate and spin, and includes synchrotron self-absorption. Shark reproduces radio luminosity functions (RLFs) at 1.4 GHz and 150 MHz for 0 $\leq$ z $\leq$ 4, and scaling relations between radio luminosity, star formation rate and infrared luminosity of galaxies in the local and distant universe in good agreement with observations. The model also reproduces observed number counts of radio sources from 150 MHz to 8.4 GHz to within a factor of two on average, though larger discrepancies are seen at the very bright fluxes at higher frequencies. We use this model to understand how the radio continuum emission from radio-quiet AGNs can affect the measured RLFs of galaxies. We find current methods to exclude AGNs from observational samples result in large fractions of radio-quiet AGNs contaminating the "star-forming galaxies" selection and a brighter end to the resulting RLFs. We investigate how this effects the infrared-radio correlation (IRRC) and show that AGN contamination can lead to evolution of the IRRC with redshift. Without this contamination our model predicts a redshift- and stellar mass-independent IRRC, except at the dwarf-galaxy regime. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS. 17 pages, 7 figures

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

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

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

Submitted 1 May, 2024; originally announced May 2024.

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

Euclid preparation: TBD. The pre-launch Science Ground Segment simulation framework

Authors: Euclid Collaboration, S. Serrano, P. Hudelot, G. Seidel, J. E. Pollack, E. Jullo, F. Torradeflot, D. Benielli, R. Fahed, T. Auphan, J. Carretero, H. Aussel, P. Casenove, F. J. Castander, J. E. Davies, N. Fourmanoit, S. Huot, A. Kara, E. Keihanen, S. Kermiche, K. Okumura, J. Zoubian, A. Ealet, A. Boucaud, H. Bretonniere , et al. (251 additional authors not shown)

Abstract: The European Space Agency's Euclid mission is one of the upcoming generation of large-scale cosmology surveys, which will map the large-scale structure in the Universe with unprecedented precision. The development and validation of the SGS pipeline requires state-of-the-art simulations with a high level of complexity and accuracy that include subtle instrumental features not accounted for previous… ▽ More The European Space Agency's Euclid mission is one of the upcoming generation of large-scale cosmology surveys, which will map the large-scale structure in the Universe with unprecedented precision. The development and validation of the SGS pipeline requires state-of-the-art simulations with a high level of complexity and accuracy that include subtle instrumental features not accounted for previously as well as faster algorithms for the large-scale production of the expected Euclid data products. In this paper, we present the Euclid SGS simulation framework as applied in a large-scale end-to-end simulation exercise named Science Challenge 8. Our simulation pipeline enables the swift production of detailed image simulations for the construction and validation of the Euclid mission during its qualification phase and will serve as a reference throughout operations. Our end-to-end simulation framework starts with the production of a large cosmological N-body & mock galaxy catalogue simulation. We perform a selection of galaxies down to I_E=26 and 28 mag, respectively, for a Euclid Wide Survey spanning 165 deg^2 and a 1 deg^2 Euclid Deep Survey. We build realistic stellar density catalogues containing Milky Way-like stars down to H<26. Using the latest instrumental models for both the Euclid instruments and spacecraft as well as Euclid-like observing sequences, we emulate with high fidelity Euclid satellite imaging throughout the mission's lifetime. We present the SC8 data set consisting of overlapping visible and near-infrared Euclid Wide Survey and Euclid Deep Survey imaging and low-resolution spectroscopy along with ground-based. This extensive data set enables end-to-end testing of the entire ground segment data reduction and science analysis pipeline as well as the Euclid mission infrastructure, paving the way to future scientific and technical developments and enhancements. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 38 pages, 25 figures, A&A submitted

Galaxy Spectra neural Network (GaSNet). II. Using Deep Learning for Spectral Classification and Redshift Predictions

Authors: Fucheng Zhong, Nicola R. Napolitano, Caroline Heneka, Rui Li, Franz Erik Bauer, Nicolas Bouche, Johan Comparat, Young-Lo Kim, Jens-Kristian Krogager, Marcella Longhetti, Jonathan Loveday, Boudewijn F. Roukema, Benedict L. Rouse, Mara Salvato, Crescenzo Tortora, Roberto J. Assef, Letizia P. Cassarà, Luca Costantin, Scott Croom, Luke J M Davies, Alexander Fritz, Guillaume Guiglion, Andrew Humphrey, Emanuela Pompei, Claudio Ricci , et al. (3 additional authors not shown)

Abstract: Large sky spectroscopic surveys have reached the scale of photometric surveys in terms of sample sizes and data complexity. These huge datasets require efficient, accurate, and flexible automated tools for data analysis and science exploitation. We present the Galaxy Spectra Network/GaSNet-II, a supervised multi-network deep learning tool for spectra classification and redshift prediction. GaSNet-… ▽ More Large sky spectroscopic surveys have reached the scale of photometric surveys in terms of sample sizes and data complexity. These huge datasets require efficient, accurate, and flexible automated tools for data analysis and science exploitation. We present the Galaxy Spectra Network/GaSNet-II, a supervised multi-network deep learning tool for spectra classification and redshift prediction. GaSNet-II can be trained to identify a customized number of classes and optimize the redshift predictions for classified objects in each of them. It also provides redshift errors, using a network-of-networks that reproduces a Monte Carlo test on each spectrum, by randomizing their weight initialization. As a demonstration of the capability of the deep learning pipeline, we use 260k Sloan Digital Sky Survey spectra from Data Release 16, separated into 13 classes including 140k galactic, and 120k extragalactic objects. GaSNet-II achieves 92.4% average classification accuracy over the 13 classes (larger than 90% for the majority of them), and an average redshift error of approximately 0.23% for galaxies and 2.1% for quasars. We further train/test the same pipeline to classify spectra and predict redshifts for a sample of 200k 4MOST mock spectra and 21k publicly released DESI spectra. On 4MOST mock data, we reach 93.4% accuracy in 10-class classification and an average redshift error of 0.55% for galaxies and 0.3% for active galactic nuclei. On DESI data, we reach 96% accuracy in (star/galaxy/quasar only) classification and an average redshift error of 2.8% for galaxies and 4.8% for quasars, despite the small sample size available. GaSNet-II can process ~40k spectra in less than one minute, on a normal Desktop GPU. This makes the pipeline particularly suitable for real-time analyses of Stage-IV survey observations and an ideal tool for feedback loops aimed at night-by-night survey strategy optimization. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 23 pages and 31 figures. The draft has been submitted to MNRAS

ProPane: Image Warping with Fire

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

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

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

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

The Hyper Suprime-Cam extended Point Spread Functions and applications

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Submitted 28 June, 2023; originally announced June 2023.

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

Efficient Reionization in a Large Hydrodynamic Galaxy Formation Simulation

Authors: James E. Davies, Simeon Bird, Simon Mutch, Yueying Ni, Yu Feng, Rupert Croft, Tiziana Di Matteo, J. Stuart B. Wyithe

Abstract: Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models th… ▽ More Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models that are pre-computed via the initial density field, or post-processed where feedback effects are ignored. We introduce Astrid-ES, a resimulation of the Astrid epoch of reionisation $20 > z > 5.5$ which includes an on-the-fly excursion-set reionization algorithm. Astrid-ES produces more accurate reionization histories without significantly impacting the computational time. This model directly utilises the star particles produced in the simulation to calculate the EoR history and includes a UV background which heats the gas particles after their reionization. We contrast the reionization topology and statistics in Astrid-ES with the previously employed parametric reionisation model, finding that in Astrid-ES, ionised regions are more correlated with galaxies, and the 21cm power-spectrum shows an increase in large scale power. We calculate the relation between the size of HII regions and the UV luminosity of the brightest galaxy within them. Prior to the overlap phase, we find a power-law fit of $\mathrm{log} (R) = -0.314 M_\mathrm{UV} - 2.550 \mathrm{log}(1+z) + 7.408$ with a standard deviation $σしぐま_R < 0.15 \mathrm{dex}$ across all mass bins. We also examine the properties of halos throughout reionization, finding that while the properties of halos in the simulation are correlated with the redshift of reionisation, they are not greatly affected by reionisation itself. △ Less

Submitted 8 August, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 12 pages, 11 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

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

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

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

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

Comments: Accepted for publication in MNRAS

L1 and off Sun-Earth line visible-light imaging of Earth-directed CMEs: An analysis of inconsistent observations

Authors: Richard A. Harrison, Jackie A. Davies, David Barnes, Christian Möstl

Abstract: The efficacy of coronal mass ejection (CME) observations as a key input to space weather forecasting is explored by comparing on and off Sun-Earth line observations from the ESA/NASA SOHO and NASA STEREO spacecraft. A comparison is made of CME catalogues based on L1 coronagraph imagery and off Sun-Earth line coronagraph and heliospheric imager (HI) observations, for the year 2011. Analysis reveals… ▽ More The efficacy of coronal mass ejection (CME) observations as a key input to space weather forecasting is explored by comparing on and off Sun-Earth line observations from the ESA/NASA SOHO and NASA STEREO spacecraft. A comparison is made of CME catalogues based on L1 coronagraph imagery and off Sun-Earth line coronagraph and heliospheric imager (HI) observations, for the year 2011. Analysis reveals inconsistencies in the identification of a number of potentially Earth-directed CMEs. The catalogues reflect our ability to identify and characterise CMEs, so any discrepancies can impact our prediction of Earth-directed CMEs. We show that 15 CMEs, which were observed by STEREO, that had estimated directions compatible with Earth-directed events, had no identified halo/partial halo counterpart listed in the L1 coronagraph CME catalogue. In-situ data confirms that for 9 of these there was a consistent L1 Interplanetary CME (ICME). The number of such "discrepant" events is significant compared to the number of ICMEs recorded at L1 in 2011, stressing the need to address space weather monitoring capabilities, particularly with the inclusion of off Sun-Earth line observation. While the study provides evidence that some halo CMEs are simply not visible in near-Earth coronagraph imagery, there is evidence that some halo CMEs viewed from L1 are compromised by preceding CME remnants or the presence of multiple-CMEs. This underlines (1) the value of multiple vantage point CME observation, and (2) the benefit of off Sun-Earth line platform heliospheric imaging, and coronagraph imaging, for the efficient identification and tracking of Earth-directed events. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: 36 pages, 6 figures, in press at AGU Space Weather, 2023

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Are the fates of supermassive black holes and galaxies determined by individual mergers, or by the properties of their host haloes?

Authors: Jonathan J. Davies, Andrew Pontzen, Robert A. Crain

Abstract: The fates of massive galaxies are tied to the evolution of their central supermassive black holes (BHs), due to the influence of AGN feedback. Correlations within simulated galaxy populations suggest that the masses of BHs are governed by properties of their host dark matter haloes, such as the binding energy and assembly time, at a given halo mass. However, the full picture must be more complex a… ▽ More The fates of massive galaxies are tied to the evolution of their central supermassive black holes (BHs), due to the influence of AGN feedback. Correlations within simulated galaxy populations suggest that the masses of BHs are governed by properties of their host dark matter haloes, such as the binding energy and assembly time, at a given halo mass. However, the full picture must be more complex as galaxy mergers have also been shown to influence the growth of BHs and the impact of AGN. In this study, we investigate this problem through a controlled experiment, using the genetic modification technique to adjust the assembly history of a Milky Way-like galaxy simulated with the EAGLE model. We change the halo assembly time (and hence the binding energy) in the absence of any disruptive merger events, and find little change in the integrated growth of the BH. We attribute this to the angular momentum support provided by a galaxy disc, which reduces the inflow of gas towards the BH and effectively decouples the BH's growth from the halo's properties. Introducing major mergers into the assembly history disrupts the disc, causing the BH to grow $\approx 4\times$ more massive and inject feedback that reduces the halo baryon fraction by a factor of $\approx 2$ and quenches star formation. Merger events appear essential to the diversity in BH masses in EAGLE, and we also show that they increase the halo binding energy; correlations between these quantities may therefore be the result of merger events. △ Less

Submitted 6 November, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

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

Galaxy quenching timescales from a forensic reconstruction of their colour evolution

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

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

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

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

Constraints on axionlike particles from a combined analysis of three flaring $\textit{Fermi}$ flat-spectrum radio quasars

Authors: James Davies, Manuel Meyer, Garret Cotter

Abstract: Many theories beyond the Standard Model of particle physics predict the existence of axionlike particles (ALPs) that mix with photons in the presence of a magnetic field. Searching for the effects of ALP-photon mixing in gamma-ray observations of blazars has provided some of the strongest constraints on ALP parameter space so far. Previously, only individual sources have been analysed. We perform… ▽ More Many theories beyond the Standard Model of particle physics predict the existence of axionlike particles (ALPs) that mix with photons in the presence of a magnetic field. Searching for the effects of ALP-photon mixing in gamma-ray observations of blazars has provided some of the strongest constraints on ALP parameter space so far. Previously, only individual sources have been analysed. We perform a combined analysis on $\textit{Fermi}$ Large Area Telescope data of three bright, flaring flat-spectrum radio quasars, with the blazar jets themselves as the dominant mixing region. For the first time, we include a full treatment of photon-photon dispersion within the jet, and account for the uncertainty in our B-field model by leaving the field strength free in the fitting. Overall, we find no evidence for ALPs, but are able to exclude the ALP parameters $m_a\lesssim200$ neV and $g_{aγがんま}\gtrsim 5 \times 10^{-12}$ GeV$^{-1}$ with 95\% confidence. △ Less

Submitted 21 April, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: 18 pages, 14 figures; published Physical Review D; matches published version

Journal ref: Phys. Rev. D 107, 083027 (2023)

DEVILS: Cosmic evolution of SED-derived metallicities and their connection to star-formation histories

Authors: Jessica E. Thorne, Aaron S. G. Robotham, Sabine Bellstedt, Luke J. M. Davies, Robin H. W. Cook, Luca Cortese, Benne Holwerda, Steven Phillipps, Malgorzata Siudek

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

Submitted 24 October, 2022; originally announced October 2022.

Comments: 25 pages, 20 figures. Accepted for publication 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

MIGHTEE-HI: Evolution of HI scaling relations of star-forming galaxies at $z<0.5$

Authors: Francesco Sinigaglia, Giulia Rodighiero, Ed Elson, Mattia Vaccari, Natasha Maddox, Bradley S. Frank, Matt J. Jarvis, Tom Oosterloo, Romeel Davé, Mara Salvato, Maarten Baes, Sabine Bellstedt, Laura Bisigello, Jordan D. Collier, Robin H. W. Cook, Luke J. M. Davies, Jacinta Delhaize, Simon P. Driver, Caroline Foster, Sushma Kurapati, Claudia del P. Lagos, Christopher Lidman, Pavel E. Mancera Piña, Martin J. Meyer, K. Moses Mogotsi , et al. (11 additional authors not shown)

Abstract: We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with… ▽ More We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with ${\rm sSFR}\equiv M_*/{\rm SFR}$), obtaining $\gtrsim 5σしぐま$ detections in most cases, the strongest HI-stacking detections to date in this redshift range. With these detections, we are able to measure scaling relations in the probed redshift interval, finding evidence for a moderate evolution from the median redshift of our sample $z_{\rm med}\sim 0.37$ to $z\sim 0$. In particular, low-$M_*$ galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 9$) experience a strong HI depletion ($\sim 0.5$ dex in $\log_{10}(M_{\rm HI}/{\rm M}_\odot)$), while massive galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 11$) keep their HI mass nearly unchanged. When looking at the star formation activity, highly star-forming galaxies evolve significantly in $M_{\rm HI}$ ($f_{\rm HI}$, where $f_{\rm HI}\equiv M_{\rm}/M_*$) at fixed SFR (sSFR), while at the lowest probed SFR (sSFR) the scaling relations show no evolution. These findings suggest a scenario in which low-$M_*$ galaxies have experienced a strong HI depletion during the last $\sim4$ Gyr, while massive galaxies have undergone a significant HI replenishment through some accretion mechanism, possibly minor mergers. Interestingly, our results are in good agreement with the predictions of the SIMBA simulation. We conclude that this work sets novel important observational constraints on galaxy scaling relations. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: 12 pages, 5 figures. Accepted for publication in ApJ Letters on 30 July 2022

Environmental effects in Herschel observations of the ionized carbon content of star forming dwarf galaxies in the Virgo cluster

Authors: Robert Minchin, Dario Fadda, Rhys Taylor, Boris Deshev, Jonathan Davies

Abstract: We use archival Herschel data to examine the singly ionized carbon ([CII]) content of 14 star forming dwarf galaxies in the Virgo cluster. We use spectral energy distribution (SED) fits to far infrared, mid infrared, near infrared, optical and ultraviolet data to derive the total infrared continuum (TIR) for these galaxies. We compare the [CII]/TIR ratio for dwarf galaxies in the central part of V… ▽ More We use archival Herschel data to examine the singly ionized carbon ([CII]) content of 14 star forming dwarf galaxies in the Virgo cluster. We use spectral energy distribution (SED) fits to far infrared, mid infrared, near infrared, optical and ultraviolet data to derive the total infrared continuum (TIR) for these galaxies. We compare the [CII]/TIR ratio for dwarf galaxies in the central part of Virgo to those in the southern part of the cluster and to galaxies with similar TIR luminosities and metallicities in the Herschel Dwarf Galaxy Survey (DGS) sample of field dwarf galaxies to look for signs of [CII] formation independent of star formation. Our analysis indicates that the sample of Virgo dwarfs in the central part of the cluster has significantly higher values of [CII]/TIR than the sample from the southern part of the cluster and the sample from the DGS, while the southern sample is consistent with the DGS. This [CII]/TIR excess implies that a significant fraction of the [CII] in the dwarf galaxies in the cluster center has an origin other than star formation and is likely to be due to environmental processes in the central part of the cluster. We also find a surprisingly strong correlation between [CII]/TIR and the local ram pressure felt by the dwarf galaxies in the cluster. In this respect, we claim that the excess [CII] we see in these galaxies is likely to be due to formation in ram pressure shocks. △ Less

Submitted 31 May, 2022; originally announced June 2022.

Comments: 20 pages, accepted for publication in the Astronomical Journal

SKYSURF: Constraints on Zodiacal Light and Extragalactic Background Light through Panchromatic HST All-Sky Surface-Brightness Measurements: II. First Limits on Diffuse Light at 1.25, 1.4, and 1.6 microns

Authors: Timothy Carleton, Rogier A. Windhorst, Rosalia O'Brien, Seth H. Cohen, Delondrae Carter, Rolf Jansen, Scott Tompkins, Richard G. Arendt, Sarah Caddy, Norman Grogin, Scott J. Kenyon, Anton Koekemoer, John MacKenty, Stefano Casertano, Luke J. M. Davies, Simon P. Driver, Eli Dwek, Alexander Kashlinsky, Nathan Miles, Rushabh Pawnikar, Nor Pirzkal, Aaron Robotham, Russell Ryan, Haley Abate, Hanga Andras-Letanovszky , et al. (6 additional authors not shown)

Abstract: We present the first results from the HST Archival Legacy project "SKYSURF." As described in Windhorst et al. 2022, SKYSURF utilizes the large HST archive to study the diffuse UV, optical, and near-IR backgrounds and foregrounds in detail. Here we utilize SKYSURF's first sky-surface brightness measurements to constrain the level of near-IR diffuse Extragalactic Background Light (EBL). Our sky-surf… ▽ More We present the first results from the HST Archival Legacy project "SKYSURF." As described in Windhorst et al. 2022, SKYSURF utilizes the large HST archive to study the diffuse UV, optical, and near-IR backgrounds and foregrounds in detail. Here we utilize SKYSURF's first sky-surface brightness measurements to constrain the level of near-IR diffuse Extragalactic Background Light (EBL). Our sky-surface brightness measurements have been verified to an accuracy of better than 1%, which when combined with systematic errors associated with HST, results in sky brightness uncertainties of $\sim$2-4% $\simeq$ 0.005 MJy/sr in each image. We put limits on the amount of diffuse EBL in three near-IR filters (F125W, F140W, and F160W) by comparing our preliminary sky measurements of $> 30,000$ images to Zodiacal light models, carefully selecting the darkest images to avoid contamination from stray light. In addition, we investigate the impact that instrumental thermal emission has on our measurements, finding that it has a limited impact on F125W and F140W measurements, whereas uncertainties in the exact thermal state of HST results in significant uncertainties in the level of astrophysical diffuse light in F160W images. When compared to the Kelsall et al. (1998) Zodiacal model, an isotropic diffuse background of $30$ nW m$^{-2}$ sr$^{-1}$ remains, whereas using the Wright (1998) Zodiacal model results in no discernible diffuse background. Based primarily on uncertainties in the foreground model subtraction, we present limits on the amount of diffuse EBL of 29 nW m$^{-2}$ sr$^{-1}$, 40 nW m$^{-2}$ sr$^{-1}$, and 29 nW m$^{-2}$ sr$^{-1}$ for F125W, F140W, and F160W respectively. While this light is generally isotropic, our modeling at this point does not distinguish between a cosmological origin or a Solar System origin (such as a dim, diffuse, spherical cloud of cometary dust). △ Less

Submitted 23 August, 2022; v1 submitted 12 May, 2022; originally announced May 2022.

Comments: Accepted with Windhorst et al. 2022 to AJ. Main figures are Fig. 10 and 11. Comments welcome!

SKYSURF: Constraints on Zodiacal Light and Extragalactic Background Light through Panchromatic HST All-Sky Surface-Brightness Measurements: I. Survey Overview and Methods

Authors: Rogier A. Windhorst, Timothy Carleton, Rosalia O'Brien, Seth H. Cohen, Delondrae Carter, Rolf Jansen, Scott Tompkins, Richard G. Arendt, Sarah Caddy, Norman Grogin, Anton Koekemoer, John MacKenty, Stefano Casertano, Luke J. M. Davies, Simon P. Driver, Eli Dwek, Alexander Kashlinsky, Scott J. Kenyon, Nathan Miles, Nor Pirzkal, Aaron Robotham, Russell Ryan, Haley Abate, Hanga Andras-Letanovszky, Jessica Berkheimer , et al. (13 additional authors not shown)

Abstract: We give an overview and describe the rationale, methods, and testing of the Hubble Space Telescope (HST) Archival Legacy project "SKYSURF." SKYSURF uses HST's unique capability as an absolute photometer to measure the ~0.2-1.7 $μみゅー$m sky surface brightness (SB) from 249,861 WFPC2, ACS, and WFC3 exposures in ~1400 independent HST fields. SKYSURF's panchromatic dataset is designed to constrain the dis… ▽ More We give an overview and describe the rationale, methods, and testing of the Hubble Space Telescope (HST) Archival Legacy project "SKYSURF." SKYSURF uses HST's unique capability as an absolute photometer to measure the ~0.2-1.7 $μみゅー$m sky surface brightness (SB) from 249,861 WFPC2, ACS, and WFC3 exposures in ~1400 independent HST fields. SKYSURF's panchromatic dataset is designed to constrain the discrete and diffuse UV to near-IR sky components: Zodiacal Light (ZL; inner Solar System), Kuiper Belt Objects (KBOs; outer Solar System), Diffuse Galactic Light (DGL), and the discrete plus diffuse Extragalactic Background Light (EBL). We outline SKYSURF's methods to: (1) measure sky-SB levels between its detected objects; (2) measure the integrated discrete EBL, most of which comes from AB$\simeq$17-22 mag galaxies; and (3) estimate how much diffuse light may exist in addition to the extrapolated discrete galaxy counts. Simulations of HST WFC3/IR images with known sky-values and gradients, realistic cosmic ray (CR) distributions, and star plus galaxy counts were processed with nine different algorithms to measure the "Lowest Estimated Sky-SB" (LES) in each image between the discrete objects. The best algorithms recover the inserted LES values within 0.2% when there are no image gradients, and within 0.2-0.4% when there are 5-10% gradients. SKYSURF requires non-standard re-processing of these HST images that includes restoring the lowest sky-level from each visit into each drizzled image. We provide a proof of concept of our methods from the WFC3/IR F125W images, where any residual diffuse light that HST sees in excess of the Kelsall et al. (1998) Zodiacal model prediction does not depend on the total object flux that each image contains. This enables us to present our first SKYSURF results on diffuse light in Carleton et al. (2022). △ Less

Submitted 25 August, 2022; v1 submitted 12 May, 2022; originally announced May 2022.

Comments: Accepted to AJ; see accompanying paper Carleton et al. 2022: arXiv:2205.06347. Comments welcome!

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

Authors: Simon P. Driver, Sabine Bellstedt, Aaron S. G. Robotham, Ivan K. Baldry, Luke J. Davies, Jochen Liske, Danail Obreschkow, Edward N. Taylor, Angus H. Wright, Mehmet Alpaslan, Steven P. Bamford, Amanda E. Bauer, Joss Bland-Hawthorn, Maciej Bilicki, Matias Bravo, Sarah Brough, Sarah Casura, Michelle E. Cluver, Matthew Colless, Christopher J. Conselice, Scott M. Croom, Jelte de Jong, Franceso D'Eugenio, Roberto De Propris, Burak Dogruel , et al. (45 additional authors not shown)

Abstract: In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19… ▽ More In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS. GAMA Data Release 4 is available at: http://www.gama-survey.org/dr4/

Galaxy mergers can initiate quenching by unlocking an AGN-driven transformation of the baryon cycle

Authors: Jonathan J. Davies, Andrew Pontzen, Robert A. Crain

Abstract: We use zoom simulations to show how merger-driven disruption of the gas disc in a galaxy provides its central active galactic nucleus (AGN) with fuel to drive outflows that entrain and expel a significant fraction of the circumgalactic medium (CGM). This in turn suppresses replenishment of the interstellar medium, causing the galaxy to quench up to several Gyr after the merger. We start by perform… ▽ More We use zoom simulations to show how merger-driven disruption of the gas disc in a galaxy provides its central active galactic nucleus (AGN) with fuel to drive outflows that entrain and expel a significant fraction of the circumgalactic medium (CGM). This in turn suppresses replenishment of the interstellar medium, causing the galaxy to quench up to several Gyr after the merger. We start by performing a zoom simulation of a present-day star-forming disc galaxy with the EAGLE galaxy formation model. Then, we re-simulate the galaxy with controlled changes to its initial conditions, using the genetic modification technique. These modifications either increase or decrease the stellar mass ratio of the galaxy's last significant merger, which occurs at $z\approx 0.74$. The halo reaches the same present-day mass in all cases, but changing the mass ratio of the merger yields markedly different galaxy and CGM properties. We find that a merger can unlock rapid growth of the central supermassive black hole if it disrupts the co-rotational motion of gas in the black hole's vicinity. Conversely, if a less disruptive merger occurs and gas close to the black hole is not disturbed, the AGN does not strongly affect the CGM, and consequently the galaxy continues to form stars. Our result illustrates how a unified view of AGN feedback, the baryon cycle and the interstellar medium is required to understand how mergers and quenching are connected over long timescales. △ Less

Submitted 21 June, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: 14 pages, 8 figures. Accepted for publication in MNRAS

Deep Extragalactic VIsible Legacy Survey (DEVILS): The emergence of bulges and decline of disk growth since $z = 1$

Authors: Abdolhosein Hashemizadeh, Simon P. Driver, Luke J. M. Davies, Aaron S. G. Robotham, Sabine Bellstedt, Rogier A. Windhorst, Matt Jarvis, Benne W. Holwerda, Malgorzata Siudek, Caroline Foster, Steven Phillipps, Jessica E. Thorne, Christian Wolf

Abstract: We present a complete structural analysis of the ellipticals (E), diffuse bulges (dB), compact bulges (cB), and disks (D) within a redshift range $0 < z < 1$, and stellar mass $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) \geq 9.5$ volume-limited sample drawn from the combined DEVILS and HST-COSMOS region. We use the {\sc ProFit} code to profile over $\sim35,000$ galaxies for which visual classificati… ▽ More We present a complete structural analysis of the ellipticals (E), diffuse bulges (dB), compact bulges (cB), and disks (D) within a redshift range $0 < z < 1$, and stellar mass $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) \geq 9.5$ volume-limited sample drawn from the combined DEVILS and HST-COSMOS region. We use the {\sc ProFit} code to profile over $\sim35,000$ galaxies for which visual classification into single or double-component was predefined in Paper-I. Over this redshift range, we see a growth in the total stellar mass density (SMD) of a factor of 1.5. At all epochs we find that the dominant structure, contributing to the total SMD, is the disk, and holds a fairly constant share of $\sim60\%$ of the total SMD from $z = 0.8$ to $z = 0.2$, dropping to $\sim30\%$ at $z = 0.0$ (representing $\sim33\%$ decline in the total disk SMD). Other classes (E, dB, and cB) show steady growth in their numbers and integrated stellar mass densities. By number, the most dramatic change across the full mass range is in the growth of diffuse bulges. In terms of total SMD, the biggest gain is an increase in massive elliptical systems, rising from 20\% at $z = 0.8$ to equal that of disks at $z = 0.0$ (30\%) representing an absolute mass growth of a factor of 2.5. Overall we see a clear picture of the emergence and growth of all three classes of spheroids over the past 8 Gyrs, and infer that in the later half of the Universe's timeline spheroid forming-processes and pathways (secular evolution, mass-accretion, and mergers) appear to dominate mass transformation over quiescent disk growth. △ Less

Submitted 6 May, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

A low-threshold ultrahigh-energy neutrino search with the Askaryan Radio Array

Authors: P. Allison, S. Archambault, J. J. Beatty, D. Z. Besson, A. Bishop, C. C. Chen, C. H. Chen, P. Chen, Y. C. Chen, B. A. Clark, W. Clay, A. Connolly, L. Cremonesi, P. Dasgupta, J. Davies, S. de Kockere, K. D. de Vries, C. Deaconu, M. A. DuVernois, J. Flaherty, E. Friedman, R. Gaior, J. Hanson, N. Harty, B. Hendricks , et al. (55 additional authors not shown)

Abstract: In the pursuit of the measurement of the still-elusive ultrahigh-energy (UHE) neutrino flux at energies of order EeV, detectors using the in-ice Askaryan radio technique have increasingly targeted lower trigger thresholds. This has led to improved trigger-level sensitivity to UHE neutrinos. Working with data collected by the Askaryan Radio Array (ARA), we search for neutrino candidates at the lowe… ▽ More In the pursuit of the measurement of the still-elusive ultrahigh-energy (UHE) neutrino flux at energies of order EeV, detectors using the in-ice Askaryan radio technique have increasingly targeted lower trigger thresholds. This has led to improved trigger-level sensitivity to UHE neutrinos. Working with data collected by the Askaryan Radio Array (ARA), we search for neutrino candidates at the lowest threshold achieved to date, leading to improved analysis-level sensitivities. A neutrino search on a data set with 208.7~days of livetime from the reduced-threshold fifth ARA station is performed, achieving a 68\% analysis efficiency over all energies on a simulated mixed-composition neutrino flux with an expected background of $0.10_{-0.04}^{+0.06}$ events passing the analysis. We observe one event passing our analysis and proceed to set a neutrino flux limit using a Feldman-Cousins construction. We show that the improved trigger-level sensitivity can be carried through an analysis, motivating the Phased Array triggering technique for use in future radio-detection experiments. We also include a projection using all available data from this detector. Finally, we find that future analyses will benefit from studies of events near the surface to fully understand the background expected for a large-scale detector. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 15 pages, 8 figures

The variation of the gas content of galaxy groups and pairs compared to isolated galaxies

Authors: Sambit Roychowdhury, Martin J. Meyer, Jonghwan Rhee, Martin A. Zwaan, Garima Chauhan, Luke J. M. Davies, Sabine Bellstedt, Simon P. Driver, Claudia del P. Lagos, Aaron S. G. Robotham, Joss Bland-Hawthorn, Richard Dodson, Benne W. Holwerda, Andrew M. Hopkins, Maritza A. Lara-Lopez, Angel R. Lopez-Sanchez, Danail Obreschkow, Kristof Rozgonyi, Matthew T. Whiting, Angus H. Wright

Abstract: We measure how the atomic gas (HI) fraction ($f_{HI}={\rm \frac{M_{HI}}{M_{*}}}$) of groups and pairs taken as single units vary with average stellar mass ($\langle {\rm M_*} \rangle$) and average star-formation rate ($\langle {\rm SFR} \rangle$), compared to isolated galaxies. The HI 21 cm emission observation are from (i) archival ALFALFA survey data covering three fields from the GAMA survey (p… ▽ More We measure how the atomic gas (HI) fraction ($f_{HI}={\rm \frac{M_{HI}}{M_{*}}}$) of groups and pairs taken as single units vary with average stellar mass ($\langle {\rm M_*} \rangle$) and average star-formation rate ($\langle {\rm SFR} \rangle$), compared to isolated galaxies. The HI 21 cm emission observation are from (i) archival ALFALFA survey data covering three fields from the GAMA survey (provides environmental and galaxy properties), and (ii) DINGO pilot survey data of one of those fields. The mean $f_{HI}$ for different units (groups/pairs/isolated galaxies) are measured in regions of the log($\langle {\rm M_*} \rangle$) -- log($\langle {\rm SFR} \rangle$) plane, relative to the z $\sim 0$ star-forming main sequence (SFMS) of individual galaxies, by stacking $f_{HI}$ spectra of individual units. For ALFALFA, $f_{HI}$ spectra of units are measured by extracting HI spectra over the full groups/pair areas and dividing by the total stellar mass of member galaxies. For DINGO, $f_{HI}$ spectra of units are measured by co-adding HI spectra of individual member galaxies, followed by division by their total stellar mass. For all units the mean $f_{HI}$ decreases as we move to higher $\langle {\rm M_*} \rangle$ along the SFMS, and as we move from above the SFMS to below it at any $\langle {\rm M_*} \rangle$. From the DINGO-based study, mean $f_{HI}$ in groups appears to be lower compared to isolated galaxies for all $\langle {\rm M_*} \rangle$ along the SFMS. From the ALFALFA-based study we find substantially higher mean $f_{HI}$ in groups compared to isolated galaxies (values for pairs being intermediate) for ${\langle{\rm M_*}\rangle}\lesssim10^{9.5}~{\rm M_{\odot}}$, indicating the presence of substantial amounts of HI not associated with cataloged member galaxies in low mass groups. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: Accepted for publication in ApJ. Main text: 26 pages, 16 figures, 6 tables

Predictions for the X-ray circumgalactic medium of edge-on discs and spheroids

Authors: Anna Nica, Benjamin D. Oppenheimer, Robert A. Crain, Ákos Bogdán, Jonathan J. Davies, William R. Forman, Ralph P. Kraft, John A. ZuHone

Abstract: We investigate how the X-ray circumgalactic medium (CGM) of present-day galaxies depends on galaxy morphology and azimuthal angle using mock observations generated from the EAGLE cosmological hydrodynamic simulation. By creating mock stacks of {\it eROSITA}-observed galaxies oriented to be edge-on, we make several observationally-testable predictions for galaxies in the stellar mass range… ▽ More We investigate how the X-ray circumgalactic medium (CGM) of present-day galaxies depends on galaxy morphology and azimuthal angle using mock observations generated from the EAGLE cosmological hydrodynamic simulation. By creating mock stacks of {\it eROSITA}-observed galaxies oriented to be edge-on, we make several observationally-testable predictions for galaxies in the stellar mass range $M_\star=10^{10.7-11.2}\;$M$_{\odot}$. The soft X-ray CGM of disc galaxies is between 60 and 100\% brighter along the semi-major axis compared to the semi-minor axis, between 10-30 kpc. This azimuthal dependence is a consequence of the hot ($T>10^6$ K) CGM being non-spherical: specifically it is flattened along the minor axis such that denser and more luminous gas resides in the disc plane and co-rotates with the galaxy. Outflows enrich and heat the CGM preferentially perpendicular to the disc, but we do not find an observationally-detectable signature along the semi-minor axis. Spheroidal galaxies have hotter CGMs than disc galaxies related to spheroids residing at higher halos masses, which may be measurable through hardness ratios spanning the $0.2-1.5$ keV band. While spheroids appear to have brighter CGMs than discs for the selected fixed $M_\star$ bin, this owes to spheroids having higher stellar and halo masses within that $M_\star$ bin, and obscures the fact that both simulated populations have similar total CGM luminosities at the exact same $M_\star$. Discs have brighter emission inside 20 kpc and more steeply declining profiles with radius than spheroids. We predict that the {\it eROSITA} 4-year all-sky survey should detect many of the signatures we predict here, although targeted follow-up observations of highly inclined nearby discs after the survey may be necessary to observe some of our azimuthally-dependent predictions. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: 12 pages, 11 figures, 1 table. Submitted to MNRAS. Comments welcome

Relevance of photon-photon dispersion within the jet for blazar axionlike particle searches

Authors: James Davies, Manuel Meyer, Garret Cotter

Abstract: Axionlike particles (ALPs) could mix with photons in the presence of astrophysical magnetic fields. Searching for this effect in gamma-ray observations of blazars has provided some of the strongest constraints on ALP parameter space so far. Previously, photon-photon dispersion of gamma-rays off of the CMB has been shown to be important for these calculations, and is universally included in ALP-pho… ▽ More Axionlike particles (ALPs) could mix with photons in the presence of astrophysical magnetic fields. Searching for this effect in gamma-ray observations of blazars has provided some of the strongest constraints on ALP parameter space so far. Previously, photon-photon dispersion of gamma-rays off of the CMB has been shown to be important for these calculations, and is universally included in ALP-photon mixing models. Here, we assess the effects of dispersion off of other photon fields within the blazar (produced by the accretion disk, the broad line region, the dust torus, starlight, and the synchrotron field) by modelling the jet and fields of the flat spectrum radio quasar 3C454.3 and propagating ALPs through the model both with and without the full dispersion calculation. We find that the full dispersion calculation can strongly affect the mixing, particularly at energies above 100 GeV -- often reducing the ALP-photon conversion probability. This could have implications for future searches planned with, e.g., the Cherenkov Telescope Array, particularly those looking for a reduced opacity of the universe at the highest energies. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 16 pages, 12 figures, accepted for publication in Physical Review D

Deep Extragalactic VIsible Legacy Survey (DEVILS): Identification of AGN through SED Fitting and the Evolution of the Bolometric AGN Luminosity Function

Authors: Jessica E. Thorne, Aaron S. G. Robotham, Luke J. M. Davies, Sabine Bellstedt, Michael J. I. Brown, Scott M. Croom, Ivan Delvecchio, Brent Groves, Matt J. Jarvis, Stanislav S. Shabala, Nick Seymour, Imogen H. Whittam, Matias Bravo, Robin H. W. Cook, Simon P. Driver, Benne Holwerda, Steven Phillipps, Malgorzata Siudek

Abstract: Active galactic nuclei (AGN) are typically identified through radio, mid-infrared, or X-ray emission or through the presence of broad and/or narrow emission lines. AGN can also leave an imprint on a galaxy's spectral energy distribution (SED) through the re-processing of photons by the dusty torus. Using the SED fitting code ProSpect with an incorporated AGN component, we fit the far ultraviolet t… ▽ More Active galactic nuclei (AGN) are typically identified through radio, mid-infrared, or X-ray emission or through the presence of broad and/or narrow emission lines. AGN can also leave an imprint on a galaxy's spectral energy distribution (SED) through the re-processing of photons by the dusty torus. Using the SED fitting code ProSpect with an incorporated AGN component, we fit the far ultraviolet to far-infrared SEDs of $\sim$494,00 galaxies in the D10-COSMOS field and $\sim$230,000 galaxies from the GAMA survey. By combining an AGN component with a flexible star formation and metallicity implementation, we obtain estimates for the AGN luminosities, stellar masses, star formation histories, and metallicity histories for each of our galaxies. We find that ProSpect can identify AGN components in 91 per cent of galaxies pre-selected as containing AGN through narrow-emission line ratios and the presence of broad lines. Our ProSpect-derived AGN luminosities show close agreement with luminosities derived for X-ray selected AGN using both the X-ray flux and previous SED fitting results. We show that incorporating the flexibility of an AGN component when fitting the SEDs of galaxies with no AGN has no significant impact on the derived galaxy properties. However, in order to obtain accurate estimates of the stellar properties of AGN host galaxies, it is crucial to include an AGN component in the SED fitting process. We use our derived AGN luminosities to map the evolution of the AGN luminosity function for $0<z<2$ and find good agreement with previous measurements and predictions from theoretical models. △ Less

Submitted 12 December, 2021; originally announced December 2021.

Comments: Published in MNRAS

Deep Extragalactic VIsible Legacy Survey (DEVILS): Evolution of the $σしぐま_{\mathrm{SFR}}$-M$_{\star}$ relation and implications for self-regulated star formation

Authors: L. J. M. Davies, J. E. Thorne, S. Bellstedt, M. Bravo, A. S. G. Robotham, S. P. Driver, R. H. W. Cook, L. Cortese, J. D'Silva, M. W. Grootes, B. W. Holwerda, A. M. Hopkins, M. J. Jarvis, C. Lidman, S. Phillipps, M. Siudek

Abstract: We present the evolution of the star-formation dispersion - stellar mass relation ($σしぐま_{SFR}$-M$_{\star}$) in the DEVILS D10 region using new measurements derived using the ProSpect spectral energy distribution fitting code. We find that $σしぐま_{SFR}$-M$_{\star}$ shows the characteristic 'U-shape' at intermediate stellar masses from 0.1<z<0.7 for a number of metrics, including using the deconvolved int… ▽ More We present the evolution of the star-formation dispersion - stellar mass relation ($σしぐま_{SFR}$-M$_{\star}$) in the DEVILS D10 region using new measurements derived using the ProSpect spectral energy distribution fitting code. We find that $σしぐま_{SFR}$-M$_{\star}$ shows the characteristic 'U-shape' at intermediate stellar masses from 0.1<z<0.7 for a number of metrics, including using the deconvolved intrinsic dispersion. A physical interpretation of this relation is the combination of stochastic star-formation and stellar feedback causing large scatter at low stellar masses and AGN feedback causing asymmetric scatter at high stellar masses. As such, the shape of this distribution and its evolution encodes detailed information about the astrophysical processes affecting star-formation, feedback and the lifecycle of galaxies. We find that the stellar mass that the minimum $σしぐま_{SFR}$ occurs evolves linearly with redshift, moving to higher stellar masses with increasing lookback time and traces the turnover in the star-forming sequence. This minimum $σしぐま_{SFR}$ point is also found to occur at a fixed specific star-formation rate (sSFR) at all epochs (sSFR~10$^{-9.6}$yr$^{-1}$). The physical interpretation of this is that there exists a maximum sSFR at which galaxies can internally self-regulate on the tight sequence of star-formation. At higher sSFRs, stochastic stellar processes begin to cause galaxies to be pushed both above and below the star-forming sequence leading to increased SFR dispersion. As the Universe evolves, a higher fraction of galaxies will drop below this sSFR threshold, causing the dispersion of the low-stellar mass end of the star-forming sequence to decrease with time. △ Less

Submitted 12 December, 2021; originally announced December 2021.

Comments: 22 pages, 14 figures, Accepted MNRAS

Quantifying the uncertainty in CME kinematics derived from geometric modelling of Heliospheric Imager data

Authors: Luke Barnard, Mathew Owens, Christopher J. Scott, Mike Lockwood, Curt A. de Koning, Tanja Amerstorfer, Jürgen Hinterreiter, Christian Möstl, Jackie Davies, Pete Riley

Abstract: Geometric modelling of Coronal Mass Ejections (CMEs) is a widely used tool for assessing their kinematic evolution. Furthermore, techniques based on geometric modelling, such as ELEvoHI, are being developed into forecast tools for space weather prediction. These models assume that solar wind structure does not affect the evolution of the CME, which is an unquantified source of uncertainty. We use… ▽ More Geometric modelling of Coronal Mass Ejections (CMEs) is a widely used tool for assessing their kinematic evolution. Furthermore, techniques based on geometric modelling, such as ELEvoHI, are being developed into forecast tools for space weather prediction. These models assume that solar wind structure does not affect the evolution of the CME, which is an unquantified source of uncertainty. We use a large number of Cone CME simulations with the HUXt solar wind model to quantify the scale of uncertainty introduced into geometric modelling and the ELEvoHI CME arrival times by solar wind structure. We produce a database of simulations, representing an average, a fast, and an extreme CME scenario, each independently propagating through 100 different ambient solar wind environments. Synthetic heliospheric imager observations of these simulations are then used with a range of geometric models to estimate the CME kinematics. The errors of geometric modelling depend on the location of the observer, but do not seem to depend on the CME scenario. In general, geometric models are biased towards predicting CME apex distances that are larger than the true value. For these CME scenarios, geometric modelling errors are minimised for an observer in the L5 region. Furthermore, geometric modelling errors increase with the level of solar wind structure in the path of the CME. The ELEvoHI arrival time errors are minimised for an observer in the L5 region, with mean absolute arrival time errors of $8.2\pm1.2$~h, $8.3\pm1.0$~h, and $5.8\pm0.9$~h for the average, fast, and extreme CME scenarios. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Comments: 37 pages, 12 figures. Accepted for publication in Space Weather

Deep Extragalactic VIsible Legacy Survey (DEVILS): DR1 Blended Spectra Search for Candidate Strong Gravitational Lenses

Authors: B. W. Holwerda, S. Knabel, J. E Thorne, S. Bellstedt, M. Siudek, L. J. M. Davies

Abstract: Here, we present a catalog of blended spectra in Data Release 1 of the Deep Extragalactic VIsible Legacy Survey (DEVILS) on the Anglo-Australian Telescope (AAT). Of the 23197 spectra, 181 showed signs of a blend of redshifts and spectral templates. We examine these blends in detail for signs of either a candidate strong lensing galaxy or a useful overlapping galaxy pair. One of the three DEVILS… ▽ More Here, we present a catalog of blended spectra in Data Release 1 of the Deep Extragalactic VIsible Legacy Survey (DEVILS) on the Anglo-Australian Telescope (AAT). Of the 23197 spectra, 181 showed signs of a blend of redshifts and spectral templates. We examine these blends in detail for signs of either a candidate strong lensing galaxy or a useful overlapping galaxy pair. One of the three DEVILS target fields, COSMOS (D10) is close to complete and it is fully imaged with Hubble Space Telescope Advanced Camera for Surveys (HST/ACS) and we visually examine the 57 blended spectra in this field in the F814W postage stamps. Nine are classical strong lensing candidates with an elliptical as the lens, out to higher redshifts than any previous search with spectroscopic surveys such as SDSS or GAMA. The gravitational lens candidate success rate similar to earlier such searches (0.1%). Strong gravitational lenses identified with blended spectroscopy have typically shown a high success rate (>70%) which make these interesting targets for future higher resolution lensing studies, monitoring for supernovae cosmography, or searches for magnified atomic hydrogen signal. △ Less

Submitted 23 November, 2021; v1 submitted 19 November, 2021; originally announced November 2021.

Comments: 25 pages, 10 figures, 6 tables, v2 fixed figure 1 and 2 mislabeling

Searching for TeV gamma-ray emission from SGR\,1935+2154 during its 2020 X-ray and radio bursting phase

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Anguner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlohr, B. Bi, M. Bottcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose , et al. (230 additional authors not shown)

Abstract: Magnetar hyperflares are the most plausible explanation for fast radio bursts (FRB) -- enigmatic powerful radio pulses with durations of several milliseconds and high brightness temperatures. The first observational evidence for this scenario was obtained in 2020 April when a FRB was detected from the direction of the Galactic magnetar and soft gamma-ray repeater SGR\,1935+2154. The FRB was preced… ▽ More Magnetar hyperflares are the most plausible explanation for fast radio bursts (FRB) -- enigmatic powerful radio pulses with durations of several milliseconds and high brightness temperatures. The first observational evidence for this scenario was obtained in 2020 April when a FRB was detected from the direction of the Galactic magnetar and soft gamma-ray repeater SGR\,1935+2154. The FRB was preceded by two gamma-ray outburst alerts by the BAT instrument aboard the Swift satellite, which triggered follow-up observations by the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. has observed SGR\,1935+2154 for 2 hr on 2020 April 28. The observations are coincident with X-ray bursts from the magnetar detected by INTEGRAL and Fermi-GBM, thus providing the first very high energy (VHE) gamma-ray observations of a magnetar in a flaring state. High-quality data acquired during these follow-up observations allow us to perform a search for short-time transients. No significant signal at energies $E>0.6$~TeV is found and upper limits on the persistent and transient emission are derived. We here present the analysis of these observations and discuss the obtained results and prospects of the H.E.S.S. follow-up program for soft gamma-ray repeaters. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: Received 2021 June 13; accepted 2021 June 28; published 2021 September 29

Multipoint interplanetary coronal mass ejections observed with Solar Orbiter, BepiColombo, Parker Solar Probe, Wind and STEREO-A

Authors: C. Möstl, A. J. Weiss, M. A. Reiss, T. Amerstorfer, R. L. Bailey, J. Hinterreiter, M. Bauer, D. Barnes, J. A. Davies, R. A. Harrison, J. L. Freiherr von Forstner, E. E. Davies, D. Heyner, T. Horbury, S. D. Bale

Abstract: We report the result of the first search for multipoint in situ and imaging observations of interplanetary coronal mass ejections (ICMEs) starting with the first Solar Orbiter (SolO) data in 2020 April - 2021 April. A data exploration analysis is performed including visualizations of the magnetic field and plasma observations made by the five spacecraft SolO, BepiColombo, Parker Solar Probe (PSP),… ▽ More We report the result of the first search for multipoint in situ and imaging observations of interplanetary coronal mass ejections (ICMEs) starting with the first Solar Orbiter (SolO) data in 2020 April - 2021 April. A data exploration analysis is performed including visualizations of the magnetic field and plasma observations made by the five spacecraft SolO, BepiColombo, Parker Solar Probe (PSP), Wind and STEREO-A, in connection with coronagraph and heliospheric imaging observations from STEREO-A/SECCHI and SOHO/LASCO. We identify ICME events that could be unambiguously followed with the STEREO-A heliospheric imagers during their interplanetary propagation to their impact at the aforementioned spacecraft, and look for events where the same ICME is seen in situ by widely separated spacecraft. We highlight two events: (1) a small streamer blowout CME on 2020 June 23 observed with a triple lineup by PSP, BepiColombo and Wind, guided by imaging with STEREO-A, and (2) the first fast CME of solar cycle 25 ($ \approx 1600$ km s$^{-1}$) on 2020 November 29 observed in situ by PSP and STEREO-A. These results are useful for modeling the magnetic structure of ICMEs and the interplanetary evolution and global shape of their flux ropes and shocks, and for studying the propagation of solar energetic particles. The combined data from these missions are already turning out to be a treasure trove for space weather research and are expected to become even more valuable with an increasing number of ICME events expected during the rise and maximum of solar cycle 25. △ Less

Submitted 4 January, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: in press at ApJ Letters (submitted 2021 September 15, revised 2021 November 24, accepted 2021 December 9). 11 pages, 3 figures, 1 table

Predicting CMEs using ELEvoHI with STEREO-HI beacon data

Authors: Maike Bauer, Tanja Amerstorfer, Jürgen Hinterreiter, Andreas J. Weiss, Jackie A. Davies, Christian Möstl, Ute V. Amerstorfer, Martin A. Reiss, Richard A. Harrison

Abstract: Being able to accurately predict the arrival of coronal mass ejections (CMEs) at Earth has been a long-standing problem in space weather research and operations. In this study, we use the ELlipse Evolution model based on Heliospheric Images (ELEvoHI) to predict the arrival time and speed of 10 CME events that were observed by HI on the STEREO-A spacecraft between 2010 and 2020. Additionally, we in… ▽ More Being able to accurately predict the arrival of coronal mass ejections (CMEs) at Earth has been a long-standing problem in space weather research and operations. In this study, we use the ELlipse Evolution model based on Heliospheric Images (ELEvoHI) to predict the arrival time and speed of 10 CME events that were observed by HI on the STEREO-A spacecraft between 2010 and 2020. Additionally, we introduce a Python tool for downloading and preparing STEREO-HI data, as well as tracking CMEs. In contrast to most previous studies, we use not only science data, which has a relatively high spatial and temporal resolution, but also low-quality beacon data, which is - in contrast to science data - provided in real-time by the STEREO-A spacecraft. We do not use data from the STEREO-B spacecraft. We get a mean absolute error of 8.81 $\pm$ 3.18 h / 59 $\pm$ 31 kms$^{-1}$ for arrival time/speed predictions using science data and 11.36 $\pm$ 8.69 h / 106 $\pm$ 61 kms$^{-1}$ for beacon data. We find that using science data generally leads to more accurate predictions, but using beacon data with the ELEvoHI model is certainly a viable choice in the absence of higher resolution real-time data. We propose that these differences could be minimized if not eliminated altogether if higher quality real-time data was available, either by enhancing the quality of the already available data or coming from a new mission carrying a HI instrument on-board. △ Less

Submitted 7 December, 2021; v1 submitted 18 August, 2021; originally announced August 2021.

Comments: 25 pages, 9 figures, submitted to AGU Space Weather 2021 August 9, published in Space Weather 2021 December 3

Journal ref: 2021, Space Weather, Volume 19, Issue 12

gammaALPs: An open-source python package for computing photon-axion-like-particle oscillations in astrophysical environments

Authors: Manuel Meyer, James Davies, Julian Kuhlmann

Abstract: Axions and axion-like particles (ALPs) are hypothetical particles that occur in extensions of the Standard Model and are candidates for cold dark matter. They could be detected through their oscillations into photons in the presence of external electromagnetic fields. gammaALPs is an open-source python framework that computes the oscillation probability between photons and axions/ALPs. In addition… ▽ More Axions and axion-like particles (ALPs) are hypothetical particles that occur in extensions of the Standard Model and are candidates for cold dark matter. They could be detected through their oscillations into photons in the presence of external electromagnetic fields. gammaALPs is an open-source python framework that computes the oscillation probability between photons and axions/ALPs. In addition to solving the photon-ALP equations of motion, gammaALPs includes models for magnetic fields in different astrophysical environments such as jets of active galactic nuclei, intra-cluster and intergalactic media, and the Milky Way. Users are also able to easily incorporate their own custom magnetic-field models. We review the basic functionality and features of gammaALPs, which is heavily based on other open-source scientific packages such as numpy and scipy. Although focused on gamma-ray energies, gammaALPs can be easily extended to arbitrary photon energies. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: Proceedings of the 37th International Cosmic Ray Conference (ICRC 2021). The online documentation is available at https://gammaalps.readthedocs.io/

LMC N132D: A mature supernova remnant with a power-law gamma-ray spectrum extending beyond 8 TeV

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun , et al. (212 additional authors not shown)

Abstract: We analyzed 252 hours of High Energy Stereoscopic System (H.E.S.S.) observations towards the supernova remnant (SNR) LMC N132D that were accumulated between December 2004 and March 2016 during a deep survey of the Large Magellanic Cloud, adding 104 hours of observations to the previously published data set to ensure a > 5 sigma detection. To broaden the gamma-ray spectral coverage required for mod… ▽ More We analyzed 252 hours of High Energy Stereoscopic System (H.E.S.S.) observations towards the supernova remnant (SNR) LMC N132D that were accumulated between December 2004 and March 2016 during a deep survey of the Large Magellanic Cloud, adding 104 hours of observations to the previously published data set to ensure a > 5 sigma detection. To broaden the gamma-ray spectral coverage required for modeling the spectral energy distribution, an analysis of Fermi-LAT Pass 8 data was also included. We unambiguously detect N132D at very high energies (VHE) with a significance of 5.7 sigma. We report the results of a detailed analysis of its spectrum and localization based on the extended H.E.S.S. data set. The joint analysis of the extended H.E.S.S and Fermi-LAT data results in a spectral energy distribution in the energy range from 1.7 GeV to 14.8 TeV, which suggests a high luminosity of N132D at GeV and TeV energies. We set a lower limit on a gamma-ray cutoff energy of 8 TeV with a confidence level of 95%. The new gamma-ray spectrum as well as multiwavelength observations of N132D when compared to physical models suggests a hadronic origin of the VHE gamma-ray emission. SNR N132D is a VHE gamma-ray source that shows a spectrum extending to the VHE domain without a spectral cutoff at a few TeV, unlike the younger oxygen-rich SNR Cassiopeia A. The gamma-ray properties of N132D may be affected by an interaction with a nearby molecular cloud that partially lies inside the 95% confidence region of the source position. [Abridged] △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 655, A7 (2021)

TeV emission of Galactic plane sources with HAWC and H.E.S.S

Authors: H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun, P. Brun , et al. (299 additional authors not shown)

Abstract: The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy gamma-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their… ▽ More The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy gamma-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their observation strategy, the size of their field of view and their angular resolution, leading to different analysis approaches. Until now, it has been unclear if the results of observations by both types of instruments are consistent: several of the recently discovered HAWC sources have been followed up by IACTs, resulting in a confirmed detection only in a minority of cases. With this paper, we go further and try to resolve the tensions between previous results by performing a new analysis of the H.E.S.S. Galactic plane survey data, applying an analysis technique comparable between H.E.S.S. and HAWC. Events above 1 TeV are selected for both datasets, the point spread function of H.E.S.S. is broadened to approach that of HAWC, and a similar background estimation method is used. This is the first detailed comparison of the Galactic plane observed by both instruments. H.E.S.S. can confirm the gamma-ray emission of four HAWC sources among seven previously undetected by IACTs, while the three others have measured fluxes below the sensitivity of the H.E.S.S. dataset. Remaining differences in the overall gamma-ray flux can be explained by the systematic uncertainties. Therefore, we confirm a consistent view of the gamma-ray sky between WCD and IACT techniques. △ Less

Submitted 8 September, 2021; v1 submitted 3 July, 2021; originally announced July 2021.

Forensic reconstruction of galaxy colour evolution and population characterisation

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

Abstract: Mapping the evolution of galaxy colours, from blue star-forming to red passive systems, is fundamental to understand the processes involved in galaxy evolution. To this end, we reconstruct the colour evolution of low-redshift galaxies, combining stellar templates with star formation and metallicity histories of galaxies from the Galaxy And Mass Assembly survey and \shark\ semi-analytic model. We u… ▽ More Mapping the evolution of galaxy colours, from blue star-forming to red passive systems, is fundamental to understand the processes involved in galaxy evolution. To this end, we reconstruct the colour evolution of low-redshift galaxies, combining stellar templates with star formation and metallicity histories of galaxies from the Galaxy And Mass Assembly survey and \shark\ semi-analytic model. We use these colour histories to robustly characterise the evolution of red and blue galaxy populations over cosmic time. Using a Gaussian Mixture Model to characterise the colour distribution at any given epoch and stellar mass, we find both observations and simulations strongly favour a model with only two populations (blue and red), with no evidence for a third "green" population. We map the evolution of mean, weight, and scatter of the blue and red populations as a function of both stellar mass and lookback time. Using our simulated galaxy catalogue as a testbed, we find that we can accurately recover galaxies colour histories up to a lookback time of $\sim6$ Gyr. We find that both populations show little change in the mean colour for low-mass galaxies, while the colours at the massive end become significantly redder with time. The stellar mass above which the galaxy population is predominantly red decreases by 0.3 dex in the last 5 Gyrs. We find a good agreement between observations and simulations, with the largest tension being that massive galaxies from \shark\ are too blue (a known issue with many galaxy evolution models). △ Less

Submitted 1 February, 2022; v1 submitted 22 June, 2021; originally announced June 2021.

Comments: 23 pages, 16 figures. Accepted for publication in MNRAS

Evidence of 100 TeV $γがんま$-ray emission from HESS J1702-420: A new PeVatron candidate

Authors: H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun, P. Brun, M. Bryan , et al. (211 additional authors not shown)

Abstract: The identification of PeVatrons, hadronic particle accelerators reaching the knee of the cosmic ray spectrum (few $10^{15}$ eV), is crucial to understand the origin of cosmic rays in the Galaxy. We provide an update on the unidentified source HESS J1702-420, a promising PeVatron candidate. We present new observations of HESS J1702-420 made with the High Energy Stereoscopic System (H.E.S.S.), and p… ▽ More The identification of PeVatrons, hadronic particle accelerators reaching the knee of the cosmic ray spectrum (few $10^{15}$ eV), is crucial to understand the origin of cosmic rays in the Galaxy. We provide an update on the unidentified source HESS J1702-420, a promising PeVatron candidate. We present new observations of HESS J1702-420 made with the High Energy Stereoscopic System (H.E.S.S.), and processed using improved analysis techniques. The analysis configuration was optimized to enhance the collection area at the highest energies. We applied a three-dimensional (3D) likelihood analysis to model the source region and adjust non thermal radiative spectral models to the $γがんま$-ray data. We also analyzed archival data from the Fermi Large Area Telescope (LAT) to constrain the source spectrum at $γがんま$-ray energies >10 GeV. We report the detection of a new source component called HESS J1702-420A, that was separated from the bulk of TeV emission at a $5.4σしぐま$ confidence level. The power law $γがんま$-ray spectrum of HESS J1702-420A extends with an index of $Γがんま=1.53\pm0.19_\text{stat}\pm0.20_\text{sys}$ and without curvature up to the energy band 64-113 TeV, in which it was detected by H.E.S.S. at a $4.0σしぐま$ confidence level. This brings evidence for the source emission up to $100\,\text{TeV}$, which makes HESS J1702-420A a compelling candidate site for the presence of extremely high energy cosmic rays. Remarkably, in a hadronic scenario, the cut-off energy of the proton distribution powering HESS J1702-420A is found to be higher than 0.5 PeV at a 95% confidence level. HESS J1702-420A becomes therefore one of the most solid PeVatron candidates detected so far in H.E.S.S. data, altough a leptonic origin of its emission could not be ruled out either. △ Less

Submitted 14 June, 2021; v1 submitted 11 June, 2021; originally announced June 2021.

Comments: Accepted for publication in the 2. Astrophysical processes section of Astronomy & Astrophysics

Deep Extragalactic VIsible Legacy Survey (DEVILS): Consistent multi-wavelength photometry for the DEVILS regions (COSMOS, XMMLSS & ECDFS)

Authors: L. J. M. Davies, J. E. Thorne, A. S. G. Robotham, S. Bellstedt, S. P. Driver, N. J. Adams, M. Bilicki, R. A. A. Bowler, M. Bravo, L. Cortese, C. Foster, M. W. Grootes, B. Häußler, A. Hashemizadeh, B. W. Holwerda, P. Hurley, M. J. Jarvis, C. Lidman, N. Maddox, M. Meyer, M. Paolillo, S. Phillipps, M. Radovich, M. Siudek, M. Vaccari , et al. (1 additional authors not shown)

Abstract: The Deep Extragalactic VIsible Legacy Survey (DEVILS) is an ongoing high-completeness, deep spectroscopic survey of $\sim$60,000 galaxies to Y$<$21.2 mag, over $\sim$6 deg2 in three well-studied deep extragalactic fields: D10 (COSMOS), D02 (XMM-LSS) and D03 (ECDFS). Numerous DEVILS projects all require consistent, uniformly-derived and state-of-the-art photometric data with which to measure galaxy… ▽ More The Deep Extragalactic VIsible Legacy Survey (DEVILS) is an ongoing high-completeness, deep spectroscopic survey of $\sim$60,000 galaxies to Y$<$21.2 mag, over $\sim$6 deg2 in three well-studied deep extragalactic fields: D10 (COSMOS), D02 (XMM-LSS) and D03 (ECDFS). Numerous DEVILS projects all require consistent, uniformly-derived and state-of-the-art photometric data with which to measure galaxy properties. Existing photometric catalogues in these regions either use varied photometric measurement techniques for different facilities/wavelengths leading to inconsistencies, older imaging data and/or rely on source detection and photometry techniques with known problems. Here we use the ProFound image analysis package and state-of-the-art imaging datasets (including Subaru-HSC, VST-VOICE, VISTA-VIDEO and UltraVISTA-DR4) to derive matched-source photometry in 22 bands from the FUV to 500μみゅーm. This photometry is found to be consistent, or better, in colour-analysis to previous approaches using fixed-size apertures (which are specifically tuned to derive colours), but produces superior total source photometry, essential for the derivation of stellar masses, star-formation rates, star-formation histories, etc. Our photometric catalogue is described in detail and, after internal DEVILS team projects, will be publicly released for use by the broader scientific community. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: 33 pages, 22 figures, Accepted to MNRAS

Search for dark matter annihilation signals from unidentified Fermi-LAT objects with H.E.S.S

Authors: H. E. S. S. Collaboration, H. Abdallah, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun , et al. (205 additional authors not shown)

Abstract: Cosmological $N$-body simulations show that Milky Way-sized galaxies harbor a population of unmerged dark matter subhalos. These subhalos could shine in gamma-rays and be eventually detected in gamma-ray surveys as unidentified sources. We performed a thorough selection among unidentified Fermi-LAT Objects (UFOs) to identify them as possible TeV-scale dark matter subhalo candidates. We search for… ▽ More Cosmological $N$-body simulations show that Milky Way-sized galaxies harbor a population of unmerged dark matter subhalos. These subhalos could shine in gamma-rays and be eventually detected in gamma-ray surveys as unidentified sources. We performed a thorough selection among unidentified Fermi-LAT Objects (UFOs) to identify them as possible TeV-scale dark matter subhalo candidates. We search for very-high-energy (E $\gtrsim$ 100 GeV) gamma-ray emissions using H.E.S.S. observations towards four selected UFOs. Since no significant very-high-energy gamma-ray emission is detected in any dataset of the four observed UFOs nor in the combined UFO dataset, strong constraints are derived on the product of the velocity-weighted annihilation cross section $\langle σしぐまv \rangle$ by the $J$-factor for the dark matter models. The 95% C.L. observed upper limits derived from combined H.E.S.S. observations reach $\langle σしぐまv \rangle J$ values of 3.7$\times$10$^{-5}$ and 8.1$\times$10$^{-6}$ GeV$^2$cm$^{-2}$s$^{-1}$ in the $W^+W^-$ and $τたう^+τたう^-$ channels, respectively, for a 1 TeV dark matter mass. Focusing on thermal WIMPs, the H.E.S.S. constraints restrict the $J$-factors to lie in the range 6.1$\times$10$^{19}$ - 2.0$\times$10$^{21}$ GeV$^2$cm$^{-5}$, and the masses to lie between 0.2 and 6 TeV in the $W^+W^-$ channel. For the $τたう^+τたう^-$ channel, the $J$-factors lie in the range 7.0$\times$10$^{19}$ - 7.1$\times$10$^{20}$ GeV$^2$cm$^{-5}$ and the masses lie between 0.2 and 0.5 TeV. Assuming model-dependent predictions from cosmological N-body simulations on the $J$-factor distribution for Milky Way-sized galaxies, the dark matter models with masses greater than 0.3 TeV for the UFO emissions can be ruled out at high confidence level. △ Less

Submitted 15 June, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: 11 pages, 7 figures, matches accepted version in The Astrophysical Journal

Journal ref: Astrophys.J., 918, 17 (2021)

Search for dark matter annihilation in the dwarf irregular galaxy WLM with H.E.S.S

Authors: H. E. S. S. Collaboration, H. Abdallah, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun , et al. (211 additional authors not shown)

Abstract: We search for an indirect signal of dark matter through very high-energy gamma rays from the Wolf-Lundmark-Melotte (WLM) dwarf irregular galaxy. The pair annihilation of dark matter particles would produce Standard Model particles in the final state such as gamma rays, which might be detected by ground-based Cherenkov telescopes. Dwarf irregular galaxies represent promising targets as they are dar… ▽ More We search for an indirect signal of dark matter through very high-energy gamma rays from the Wolf-Lundmark-Melotte (WLM) dwarf irregular galaxy. The pair annihilation of dark matter particles would produce Standard Model particles in the final state such as gamma rays, which might be detected by ground-based Cherenkov telescopes. Dwarf irregular galaxies represent promising targets as they are dark matter dominated objects with well measured kinematics and small uncertainties on their dark matter distribution profiles. In 2018, the H.E.S.S. five-telescope array observed the dwarf irregular galaxy WLM for 18 hours. We present the first analysis based on data obtained from an imaging atmospheric Cherenkov telescope for this subclass of dwarf galaxy. As we do not observe any significant excess in the direction of WLM, we interpret the result in terms of constraints on the velocity-weighted cross section for dark matter pair annihilation as a function of the dark matter particle mass for various continuum channels as well as the prompt gamma-gamma emission. For the $τたう^+τたう^-$ channel the limits reach a $\langle σしぐまv \rangle$ value of about $4\times 10^{-22}$ cm3s-1 for a dark matter particle mass of 1 TeV. For the prompt gamma-gamma channel, the upper limit reaches a $\langle σしぐまv \rangle$ value of about $5 \times10^{-24}$ cm3s-1 for a mass of 370 GeV. These limits represent an improvement of up to a factor 200 with respect to previous results for the dwarf irregular galaxies for TeV dark matter search. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 13 pages, 5 figures

Journal ref: Phys. Rev. D 103, 102002 (2021)

A journey of exploration to the polar regions of a star: probing the solar poles and the heliosphere from high helio-latitude

Authors: Louise Harra, Vincenzo Andretta, Thierry Appourchaux, Frédéric Baudin, Luis Bellot-Rubio, Aaron C. Birch, Patrick Boumier, Robert H. Cameron, Matts Carlsson, Thierry Corbard, Jackie Davies, Andrew Fazakerley, Silvano Fineschi, Wolfgang Finsterle, Laurent Gizon, Richard Harrison, Donald M. Hassler, John Leibacher, Paulett Liewer, Malcolm MacDonald, Milan Maksimovic, Neil Murphy, Giampiero Naletto, Giuseppina Nigro, Christopher Owen , et al. (7 additional authors not shown)

Abstract: A mission to view the solar poles from high helio-latitudes (above 60$^\circ$) will build on the experience of Solar Orbiter as well as a long heritage of successful solar missions and instrumentation (e.g. SOHO \cite{SOHO}, STEREO \cite{stereo}, Hinode \cite{Hinode}, SDO \cite{SDO}), but will focus for the first time on the solar poles, enabling scientific investigations that cannot be done by an… ▽ More A mission to view the solar poles from high helio-latitudes (above 60$^\circ$) will build on the experience of Solar Orbiter as well as a long heritage of successful solar missions and instrumentation (e.g. SOHO \cite{SOHO}, STEREO \cite{stereo}, Hinode \cite{Hinode}, SDO \cite{SDO}), but will focus for the first time on the solar poles, enabling scientific investigations that cannot be done by any other mission. One of the major mysteries of the Sun is the solar cycle. The activity cycle of the Sun drives the structure and behaviour of the heliosphere and is, of course, the driver of space weather. In addition, solar activity and variability provides fluctuating input into the Earth climate models, and these same physical processes are applicable to stellar systems hosting exoplanets. One of the main obstructions to understanding the solar cycle, and hence all solar activity, is our current lack of understanding of the polar regions. In this White Paper, submitted to the European Space Agency in response to the Voyage 2050 call, we describe a mission concept that aims to address this fundamental issue. In parallel, we recognise that viewing the Sun from above the polar regions enables further scientific advantages, beyond those related to the solar cycle, such as unique and powerful studies of coronal mass ejection processes, from a global perspective, and studies of coronal structure and activity in polar regions. Not only will these provide important scientific advances for fundamental stellar physics research, they will feed into our understanding of impacts on the Earth and other planets' space environment. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 29 pages, 9 figures, White paper for ESA Voyage 2050

Deep Extragalactic VIsible Legacy Survey (DEVILS): Stellar Mass Growth by Morphological Type since $z = 1$

Authors: Abdolhosein Hashemizadeh, Simon P. Driver, Luke J. M. Davies, Aaron S. G. Robotham, Sabine Bellstedt, Rogier A. Windhorst, Malcolm Bremer, Steven Phillipps, Matt Jarvis, Benne W. Holwerda, Claudia del P. Lagos, Soheil Koushan, Malgorzata Siudek, Natasha Maddox, Jessica E. Thorne, Pascal Elahi

Abstract: Using high-resolution Hubble Space Telescope imaging data, we perform a visual morphological classification of $\sim 36,000$ galaxies at $z < 1$ in the DEVILS/COSMOS region. As the main goal of this study, we derive the stellar mass function (SMF) and stellar mass density (SMD) sub-divided by morphological types. We find that visual morphological classification using optical imaging is increasingl… ▽ More Using high-resolution Hubble Space Telescope imaging data, we perform a visual morphological classification of $\sim 36,000$ galaxies at $z < 1$ in the DEVILS/COSMOS region. As the main goal of this study, we derive the stellar mass function (SMF) and stellar mass density (SMD) sub-divided by morphological types. We find that visual morphological classification using optical imaging is increasingly difficult at $z > 1$ as the fraction of irregular galaxies and merger systems (when observed at rest-frame UV/blue wavelengths) dramatically increases. We determine that roughly two-thirds of the total stellar mass of the Universe today was in place by $z \sim 1$. Double-component galaxies dominate the SMD at all epochs and increase in their contribution to the stellar mass budget to the present day. Elliptical galaxies are the second most dominant morphological type and increase their SMD by $\sim 2.5$ times, while by contrast, the pure-disk population significantly decreases by $\sim 85\%$. According to the evolution of both high- and low-mass ends of the SMF, we find that mergers and in-situ evolution in disks are both present at $z < 1$, and conclude that double-component galaxies are predominantly being built by the in-situ evolution in disks (apparent as the growth of the low-mass end with time), while mergers are likely responsible for the growth of ellipticals (apparent as the increase of intermediate/high-mass end). △ Less

Submitted 26 February, 2021; originally announced February 2021.

Comments: 29 pages, 22 figures

GAMA/DEVILS: Constraining the cosmic star-formation history from improved measurements of the 0.3-2.2 micron Extragalactic Background Light

Authors: Soheil Koushan, Simon P. Driver, Sabine Bellstedt, Luke J. Davies, Aaron S. G. Robotham, Claudia del P Lagos, Abdolhosein Hashemizadeh, Danail Obreschkow, Jessica E. Thorne, Malcolm Bremer, B. W. Holwerda, Andrew M. Hopkins, Matt J. Jarvis, Malgorzata Siudek, Rogier A. Windhorst

Abstract: We present a revised measurement of the optical extragalactic background light (EBL), based on the contribution of resolved galaxies to the integrated galaxy light (IGL). The cosmic optical background radiation (COB), encodes the light generated by star-formation, and provides a wealth of information about the cosmic star formation history (CSFH). We combine wide and deep galaxy number counts from… ▽ More We present a revised measurement of the optical extragalactic background light (EBL), based on the contribution of resolved galaxies to the integrated galaxy light (IGL). The cosmic optical background radiation (COB), encodes the light generated by star-formation, and provides a wealth of information about the cosmic star formation history (CSFH). We combine wide and deep galaxy number counts from the Galaxy And Mass Assembly survey (GAMA) and Deep Extragalactic VIsible Legacy Survey (DEVILS), along with the Hubble Space Telescope (HST) archive and other deep survey datasets, in 9 multi-wavelength filters to measure the COB in the range from 0.35 micron to 2.2 micron. We derive the luminosity density in each band independently and show good agreement with recent and complementary estimates of the optical-EBL from very high-energy (VHE) experiments. Our error analysis suggests that the IGL and Gamma-ray measurements are now fully consistent to within ~10%, suggesting little need for any additional source of diffuse light beyond the known galaxy population. We use our revised IGL measurements to constrain the cosmic star-formation history, and place amplitude constraints on a number of recent estimates. As a consistency check, we can now demonstrate convincingly, that the CSFH, stellar mass growth, and the optical-EBL provide a fully consistent picture of galaxy evolution. We conclude that the peak of star-formation rate lies in the range 0.066-0.076 Msol/yr/Mpc^3 at a lookback time of 9.1 to 10.9 Gyrs. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: Accepted for publication in MNRAS

Galaxy And Mass Assembly (GAMA): The inferred mass--metallicity relation from z=0 to 3.5 via forensic SED fitting

Authors: Sabine Bellstedt, Aaron S. G. Robotham, Simon P. Driver, Jessica E. Thorne, Luke J. M. Davies, Benne W. Holwerda, Andrew M. Hopkins, Maritza A. Lara-Lopez, Ángel R. López-Sánchez, Steven Phillipps

Abstract: We analyse the metallicity histories of ~4,500 galaxies from the GAMA survey at z<0.06 modelled by the SED-fitting code ProSpect using an evolving metallicity implementation. These metallicity histories, in combination with the associated star formation histories, allow us to analyse the inferred gas-phase mass--metallicity relation. Furthermore, we extract the mass--metallicity relation at a sequ… ▽ More We analyse the metallicity histories of ~4,500 galaxies from the GAMA survey at z<0.06 modelled by the SED-fitting code ProSpect using an evolving metallicity implementation. These metallicity histories, in combination with the associated star formation histories, allow us to analyse the inferred gas-phase mass--metallicity relation. Furthermore, we extract the mass--metallicity relation at a sequence of epochs in cosmic history, to track the evolving mass--metallicity relation with time. Through comparison with observations of gas-phase metallicity over a large range of redshifts, we show that, remarkably, our forensic SED analysis has produced an evolving mass--metallicity relationship that is consistent with observations at all epochs. We additionally analyse the three dimensional mass--metallicity--SFR space, showing that galaxies occupy a clearly defined plane. This plane is shown to be subtly evolving, displaying an increased tilt with time caused by general enrichment, and also the slowing down of star formation with cosmic time. This evolution is most apparent at lookback times greater than 7 Gyr. The trends in metallicity recovered in this work highlight that the evolving metallicity implementation used within the SED fitting code ProSpect produces reasonable metallicity results over the history of a galaxy. This is expected to provide a significant improvement to the accuracy of the SED fitting outputs. △ Less

Submitted 15 March, 2021; v1 submitted 23 February, 2021; originally announced February 2021.

Comments: 19 pages, 14 figures (inc. 2 animations), accepted for publication in MNRAS