Showing 1–50 of 69 results for author: Bellstedt, S

Wide Area VISTA Extra-galactic Survey (WAVES): Unsupervised star-galaxy separation on the WAVES-Wide photometric input catalogue using UMAP and ${\rm{\scriptsize HDBSCAN}}$

Authors: Todd L. Cook, Behnood Bandi, Sam Philipsborn, Jon Loveday, Sabine Bellstedt, Simon P. Driver, Aaron S. G. Robotham, Maciej Bilicki, Gursharanjit Kaur, Elmo Tempel, Ivan Baldry, Daniel Gruen, Marcella Longhetti, Angela Iovino, Benne W. Holwerda, Ricardo Demarco

Abstract: Star-galaxy separation is a crucial step in creating target catalogues for extragalactic spectroscopic surveys. A classifier biased towards inclusivity risks including spurious stars, wasting fibre hours, while a more conservative classifier might overlook galaxies, compromising completeness and hence survey objectives. To avoid bias introduced by a training set in supervised methods, we employ an… ▽ More Star-galaxy separation is a crucial step in creating target catalogues for extragalactic spectroscopic surveys. A classifier biased towards inclusivity risks including spurious stars, wasting fibre hours, while a more conservative classifier might overlook galaxies, compromising completeness and hence survey objectives. To avoid bias introduced by a training set in supervised methods, we employ an unsupervised machine learning approach. Using photometry from the Wide Area VISTA Extragalactic Survey (WAVES)-Wide catalogue comprising 9-band $u-K_s$ data, we create a feature space with colours, fluxes, and apparent size information extracted by ${\rm P{\scriptsize RO} F{\scriptsize OUND}}$. We apply the non-linear dimensionality reduction method UMAP (Uniform Manifold Approximation and Projection) combined with the classifier ${\rm{\scriptsize HDBSCAN}}$ to classify stars and galaxies. Our method is verified against a baseline colour and morphological method using a truth catalogue from Gaia, SDSS, GAMA, and DESI. We correctly identify 99.72% of galaxies within the AB magnitude limit of $Z = 21.2$, with an F1 score of 0.9970 across the entire ground truth sample, compared to 0.9871 from the baseline method. Our method's higher purity (0.9966) compared to the baseline (0.9780) increases efficiency, identifying 11% fewer galaxy or ambiguous sources, saving approximately 70,000 fibre hours on the 4MOST instrument. We achieve reliable classification statistics for challenging sources including quasars, compact galaxies, and low surface brightness galaxies, retrieving 95.1%, 84.6%, and 99.5% of them respectively. Angular clustering analysis validates our classifications, showing consistency with expected galaxy clustering, regardless of the baseline classification. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS; 17 pages, 11 figures

Unveiling the (in)consistencies among the galaxy stellar mass function, star formation histories, satellite abundances and intracluster light from a semi-empirical perspective

Authors: Hao Fu, Francesco Shankar, Mohammadreza Ayromlou, Ioanna Koutsouridou, Andrea Cattaneo, Caroline Bertemes, Sabine Bellstedt, Ignacio Martín-Navarro, Joel Leja, Viola Allevato, Mariangela Bernardi, Lumen Boco, Paola Dimauro, Carlotta Gruppioni, Andrea Lapi, Nicola Menci, Iván Muñoz Rodríguez, Annagrazia Puglisi, Alba V. Alonso-Tetilla

Abstract: In a hierarchical, dark matter-dominated Universe, stellar mass functions (SMFs), galaxy merger rates, star formation histories (SFHs), satellite abundances, and intracluster light, should all be intimately connected observables. However, the systematics affecting observations still prevent universal and uniform measurements of, for example, the SMF and the SFHs, inevitably preventing theoretical… ▽ More In a hierarchical, dark matter-dominated Universe, stellar mass functions (SMFs), galaxy merger rates, star formation histories (SFHs), satellite abundances, and intracluster light, should all be intimately connected observables. However, the systematics affecting observations still prevent universal and uniform measurements of, for example, the SMF and the SFHs, inevitably preventing theoretical models to compare with multiple data sets robustly and simultaneously. We here present our holistic semi-empirical model DECODE (Discrete statistical sEmi-empiriCal mODEl) that converts via abundance matching dark matter merger trees into galaxy assembly histories, using different SMFs in input and predicting all other observables in output in a fully data-driven and self-consistent fashion with minimal assumptions. We find that: 1) weakly evolving or nearly constant SMFs below the knee ($M_\star \lesssim 10^{11} \, M_\odot$) are the best suited to generate star formation histories aligned with those inferred from MaNGA, SDSS, GAMA, and, more recently, JWST; 2) the evolution of satellites after infall only affects the satellite abundances and star formation histories of massive central galaxies but not their merger histories; 3) the resulting SFR-$M_\star$ relation is lower in normalization by a factor of $\sim 2$ with respect to observations, with a flattening at high masses more pronounced in the presence of mergers; 4) the latest data on intracluster light can be reproduced if mass loss from mergers is included in the models. Our findings are pivotal in acting as pathfinder to test the self-consistency of the high-quality data from, e.g., JWST and Euclid. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 21 pages, 12 figures, MNRAS accepted

The MAGPI Survey: Massive slow rotator population in place by $z \sim 0.3$

Authors: Caro Derkenne, Richard M. McDermid, Francesco D'Eugenio, Caroline Foster, Aman Khalid, Katherine E. Harborne, Jesse van de Sande, Scott M. Croom, Claudia D. P. Lagos, Sabine Bellstedt, J. Trevor Mendel, Marcie Mun, Emily Wisnioski, Ryan S. Bagge, Andrew J. Battisti, Joss Bland-Hawthorn, Anna Ferré-Mateu, Yingjie Peng, Giulia Santucci, Sarah M. Sweet, Sabine Thater, Lucas M. Valenzuela, Bodo Ziegler

Abstract: We use the `Middle Ages Galaxy Properties with Integral field spectroscopy' (MAGPI) survey to investigate whether galaxies have evolved in the distribution of their stellar angular momentum in the past 3-4 Gyr, as probed by the observational proxy for spin, $λ_{R}$. We use 2D stellar kinematics to measure $λ_{R}$ along with detailed photometric models to estimate galaxy ellipticity. The combinatio… ▽ More We use the `Middle Ages Galaxy Properties with Integral field spectroscopy' (MAGPI) survey to investigate whether galaxies have evolved in the distribution of their stellar angular momentum in the past 3-4 Gyr, as probed by the observational proxy for spin, $λ_{R}$. We use 2D stellar kinematics to measure $λ_{R}$ along with detailed photometric models to estimate galaxy ellipticity. The combination of these measurements quantifies the kinematic classes of `fast rotators' and the rarer `slow rotators', which show no regular rotation in their line-of-sight velocity fields. We compare 51 MAGPI galaxies with $\log_{10} (M_{\star}/\mathrm{M}_\odot) > 10$ to carefully drawn samples of MaNGA galaxies in the local Universe, selected to represent possible descendants of the MAGPI progenitors. The EAGLE simulations are used to identify possible evolutionary pathways between the two samples, explicitly accounting for progenitor bias in our results and the varied evolutionary pathways a galaxy might take between the two epochs. We find that the occurrence of slow rotating galaxies is unchanged between the MAGPI ($z \sim 0.3$) and MaNGA ($z \sim 0$) samples, suggesting the massive slow rotator population was already in place $\sim 4$ Gyr ago and has not accumulated since. There is a hint of the MAGPI sample having an excess of high $λ_{R}$ galaxies compared to the MaNGA sample, corresponding to more ordered rotation, but statistically the samples are not significantly different. The large-scale stellar kinematics, as quantified through the $λ_{R}$ parameter, of galaxies at $z \sim 0.3$ have already evolved into the diversity of structures seen today in the local Universe. △ Less

Submitted 4 June, 2024; originally announced June 2024.

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

The MAGPI Survey: Using kinematic asymmetries in stars and gas to dissect drivers of galaxy dynamical evolution

Authors: R. S. Bagge, C. Foster, F. D'Eugenio, A. Battisti, S. Bellstedt, C. Derkenne, S. Vaughan, T. Mendel, S. Barsanti, K. E. Harborne, S. M. Croom, J. Bland-Hawthorn, K. Grasha, C. D. P. Lagos, S. M. Sweet, A. Mailvaganam, T. Mukherjee, L. M. Valenzuela, J. van de Sande, E. Wisnioski, T. Zafar

Abstract: We present a study of kinematic asymmetries from the integral field spectroscopic surveys MAGPI and SAMI. By comparing the asymmetries in the ionsied gas and stars, we aim to disentangle the physical processes that contribute to kinematic disturbances. We normalise deviations from circular motion by $S_{05}$, allowing us to study kinematic asymmetries in the stars and gas, regardless of kinematic… ▽ More We present a study of kinematic asymmetries from the integral field spectroscopic surveys MAGPI and SAMI. By comparing the asymmetries in the ionsied gas and stars, we aim to disentangle the physical processes that contribute to kinematic disturbances. We normalise deviations from circular motion by $S_{05}$, allowing us to study kinematic asymmetries in the stars and gas, regardless of kinematic temperature. We find a similar distribution of stellar asymmetries in galaxies where we do and do not detect ionised gas, suggesting that whatever is driving the stellar asymmetries does not always lead to gas removal. In both MAGPI and SAMI, we find an anti-correlation between stellar asymmetry and stellar mass, that is absent in the gas asymmetries. After stellar mass and mean-stellar-age matching distributions, we find that at all stellar masses, MAGPI galaxies display larger stellar asymmetry compared to SAMI galaxies. In both MAGPI and SAMI galaxies, we find that star-forming galaxies with old mean-stellar-ages typically have larger asymmetries in their gas compared to their stars, whereas galaxies with young mean-stellar-ages have larger asymmetries in their stars compared to their gas. We suggest that this results from continuous, clumpy accretion of gas. △ Less

Submitted 18 May, 2024; originally announced May 2024.

Comments: 12 pages, 10 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

The MAGPI Survey: Evolution of radial trends in star formation activity across cosmic time

Authors: Marcie Mun, Emily Wisnioski, Andrew J. Battisti, J. Trevor Mendel, Sara L. Ellison, Edward N. Taylor, Claudia D. P. Lagos, Katherine E. Harborne, Caroline Foster, Scott M. Croom, Sabine Bellstedt, Stefania Barsanti, Anshu Gupta, Lucas M. Valenzuela, Qian-Hui Chen, Kathryn Grasha, Tamal Mukherjee, Hye-Jin Park, Piyush Sharda, Sarah M. Sweet, Rhea-Silvia Remus, Tayyaba Zafar

Abstract: Using adaptive optics with the Multi-Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT), the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey allows us to study the spatially resolved Universe at a crucial time of ~4 Gyr ago ($z$ ~ 0.3) when simulations predict the greatest diversity in evolutionary pathways for galaxies. We investigate the radial tre… ▽ More Using adaptive optics with the Multi-Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT), the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey allows us to study the spatially resolved Universe at a crucial time of ~4 Gyr ago ($z$ ~ 0.3) when simulations predict the greatest diversity in evolutionary pathways for galaxies. We investigate the radial trends in the star formation (SF) activity and luminosity-weighted stellar ages as a function of offset from the star-forming main sequence (SFMS) for a total of 294 galaxies. Using both H$α$ emission and the 4000 Angstrom break (i.e., D4000) as star formation rate (SFR) tracers, we find overall flat radial profiles for galaxies lying on and above the SFMS, suggestive of physical processes that enhance/regulate SF throughout the entire galaxy disc. However, for galaxies lying below the SFMS, we find positive gradients in SF suggestive of inside-out quenching. Placing our results in context with results from other redshift regimes suggests an evolution in radial trends at $z$ ~ 0.3 for SF galaxies above the SFMS, from uniformly enhanced SF at $z$ ~ 1 and $z$ ~ 0.3 to centrally enhanced SF at $z$ ~ 0 (when averaged over a wide range of mass). We also capture higher local SFRs for galaxies below the SFMS compared to that of $z$ ~ 0, which can be explained by a larger population of quenched satellites in the local Universe and/or different treatments of limitations set by the D4000-sSFR relation. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 19 pages, 15 figures, 4 tables, accepted for publication in MNRAS

The Mass Density of MgII Absorbers from the Australian Dark Energy Survey

Authors: Asif Abbas, Christopher W. Churchill, Glenn G. Kacprzak, Christopher Lidman, Susanna Guatelli, Sabine Bellstedt

Abstract: We present an all-southern sky survey for MgII doublet absorbers in 951 z < 4 AGN/quasar spectra from the Australian Dark Energy Survey (OzDES). The spectral resolution ranges from R = 1400-1700 over the wavelengths 3700 A-8800 A. The survey has a 5sigma detection completeness of 50% and above for rest-frame equivalent widths W_r(2796) >= 0.3 A. We studied 656 MgII absorption systems over the reds… ▽ More We present an all-southern sky survey for MgII doublet absorbers in 951 z < 4 AGN/quasar spectra from the Australian Dark Energy Survey (OzDES). The spectral resolution ranges from R = 1400-1700 over the wavelengths 3700 A-8800 A. The survey has a 5sigma detection completeness of 50% and above for rest-frame equivalent widths W_r(2796) >= 0.3 A. We studied 656 MgII absorption systems over the redshift range 0.33 < z < 2.19 with equivalent widths 0.3 < W_r(2796) < 3.45 A. The equivalent width distribution is well fit by an exponential function with W* = 0.76 +/- 0.04 A and the redshift path density exhibits very little evolution. Overall, our findings are consistent with the large, predominantly northern-sky, surveys of MgII absorbers. We developed and implemented a Monte Carlo model informed by a high-resolution MgII survey for determining the MgII mass density, Omega_MgII. We found Omega_MgII ~ 5 x 10^-7 with no evidence of evolution over a ~7 Gyr time span following Cosmic Noon. Incorporating measurements covering 2.0 < z < 6.4 from the literature, we extended our insights into MgII mass density evolution from the end of reionization well past the Cosmic Noon epoch. The presented Monte Carlo model has potential for advancing our knowledge of the evolution of mass densities of metal-ions common to quasar absorption line studies, as it exploits the efficiency of large low-resolution surveys while requiring only small samples from expensive high-resolution surveys. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 27 pages, 15 figures, 3 tables. Accepted on Pi Day 2024

The LOFAR-eFEDS survey: The incidence of radio and X-ray AGN and the disk-jet connection

Authors: Z. Igo, A. Merloni, D. Hoang, J. Buchner, T. Liu, M. Salvato, R. Arcodia, S. Bellstedt, M. Brüggen, J. H. Croston, F. de Gasperin, A. Georgakakis, M. J. Hardcastle, K. Nandra, Q. Ni, T. Pasini, T. Shimwell, J. Wolf

Abstract: Radio jets are present in a diverse sample of AGN. However, the mechanisms of jet powering are not fully understood, and it is yet unclear to what extent they obey mass-invariant scaling relations, similar to those found for the triggering and fuelling of X-ray selected AGN. We study the incidence of eROSITA/eFEDS X-ray and LOFAR radio AGN as a function of several stellar mass normalised AGN power… ▽ More Radio jets are present in a diverse sample of AGN. However, the mechanisms of jet powering are not fully understood, and it is yet unclear to what extent they obey mass-invariant scaling relations, similar to those found for the triggering and fuelling of X-ray selected AGN. We study the incidence of eROSITA/eFEDS X-ray and LOFAR radio AGN as a function of several stellar mass normalised AGN power indicators. A new sample of radio AGN from the LOFAR-eFEDS survey is defined and we publicly release this catalogue, including host galaxy counterparts from the Legacy Survey DR9, LOFAR radio morphologies and host galaxy properties from the complete, spectroscopic (z<0.4) GAMA09 survey. The fraction of GAMA09 galaxies hosting radio, X-ray and both radio and X-ray AGN are calculated as a function of the specific black hole kinetic ($λ_{\rm Jet}$) and radiative ($λ_{\rm Edd}$) power. The incidence of eFEDS X-ray AGN as a function of $λ_{\rm Edd}$ shows the same mass-invariance as found in past studies. Meanwhile, radio AGN, regardless of their morphology, are more likely to be hosted in more massive galaxies, at all $λ_{\rm Jet}$. Across the stellar mass range, the compact radio AGN incidence follows the same power-law distribution, showing that it is not only high mass galaxies that host high power radio AGN and vice versa. On the other hand, the incidence of compact and complex radio AGN is boosted at the highest jet powers, diverging from a simple power-law. Interestingly, this increased incidence cannot be explained by more powerful radio AGN lying in more dense environments which could naturally boost their radio luminosity. Overall, we show that statistical incidence studies are a powerful method to probe disk-jet coupling for different AGN accretion modes, although future work on a more reliable determination of jet power for diverse samples of radio AGN is needed. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 30 pages, 32 figures, accepted for publication in A&A

Journal ref: A&A 686, A43 (2024)

EMU/GAMA: A Technique for Detecting Active Galactic Nuclei in Low Mass Systems

Authors: Jahang Prathap, Andrew M. Hopkins, Aaron S. G. Robotham, Sabine Bellstedt, José Afonso, Ummee T. Ahmed, Maciej Bilicki, Malcolm N. Bremer, Sarah Brough, Michael J. I. Brown, Yjan Gordon, Benne W. Holwerda, Denis Leahy, Ángel R. López-Sánchez, Joshua R. Marvil, Tamal Mukherjee, Isabella Prandoni, Stanislav S. Shabala, Tessa Vernstrom, Tayyaba Zafar

Abstract: We propose a new method for identifying active galactic nuclei (AGN) in low mass ($\rm M_*\leq10^{10}M_\odot$) galaxies. This method relies on spectral energy distribution (SED) fitting to identify galaxies whose radio flux density has an excess over that expected from star formation alone. Combining data in the Galaxy and Mass Assembly (GAMA) G23 region from GAMA, Evolutionary Map of the Universe… ▽ More We propose a new method for identifying active galactic nuclei (AGN) in low mass ($\rm M_*\leq10^{10}M_\odot$) galaxies. This method relies on spectral energy distribution (SED) fitting to identify galaxies whose radio flux density has an excess over that expected from star formation alone. Combining data in the Galaxy and Mass Assembly (GAMA) G23 region from GAMA, Evolutionary Map of the Universe (EMU) early science observations, and Wide-field Infrared Survey Explorer (WISE), we compare this technique with a selection of different AGN diagnostics to explore the similarities and differences in AGN classification. We find that diagnostics based on optical and near-infrared criteria (the standard BPT diagram, the WISE colour criterion, and the mass-excitation, or MEx diagram) tend to favour detection of AGN in high mass, high luminosity systems, while the ``ProSpect'' SED fitting tool can identify AGN efficiently in low mass systems. We investigate an explanation for this result in the context of proportionally lower mass black holes in lower mass galaxies compared to higher mass galaxies and differing proportions of emission from AGN and star formation dominating the light at optical and infrared wavelengths as a function of galaxy stellar mass. We conclude that SED-derived AGN classification is an efficient approach to identify low mass hosts with low radio luminosity AGN. △ Less

Submitted 18 February, 2024; originally announced February 2024.

Comments: 16 pages, 12 figures, 5 tables. Accepted for publication in PASA

The MAGPI Survey: Drivers of kinematic asymmetries in the ionised gas of $z\sim0.3$ star-forming galaxies

Authors: R. S. Bagge, C. Foster, A. Battisti, S. Bellstedt, M. Mun, K. Harborne, S. Barsanti, T. Mendel, S. Brough, S. M. Croom, C. D. P. Lagos, T. Mukherjee, Y. Peng, R-S. Remus, G. Santucci, P. Sharda, S. Thater, J. van de Sande, L. M. Valenzuela E. Wisnioski T. Zafar, B. Ziegler

Abstract: Galaxy gas kinematics are sensitive to the physical processes that contribute to a galaxy's evolution. It is expected that external processes will cause more significant kinematic disturbances in the outer regions, while internal processes will cause more disturbances for the inner regions. Using a subsample of 47 galaxies ($0.27<z<0.36$) from the Middle Ages Galaxy Properties with Integral Field… ▽ More Galaxy gas kinematics are sensitive to the physical processes that contribute to a galaxy's evolution. It is expected that external processes will cause more significant kinematic disturbances in the outer regions, while internal processes will cause more disturbances for the inner regions. Using a subsample of 47 galaxies ($0.27<z<0.36$) from the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, we conduct a study into the source of kinematic disturbances by measuring the asymmetry present in the ionised gas line-of-sight velocity maps at the $0.5R_e$ (inner regions) and $1.5R_e$ (outer regions) elliptical annuli. By comparing the inner and outer kinematic asymmetries, we aim to better understand what physical processes are driving the asymmetries in galaxies. We find the local environment plays a role in kinematic disturbance, in agreement with other integral field spectroscopy studies of the local universe, with most asymmetric systems being in close proximity to a more massive neighbour. We do not find evidence suggesting that hosting an Active Galactic Nucleus (AGN) contributes to asymmetry within the inner regions, with some caveats due to emission line modelling. In contrast to previous studies, we do not find evidence that processes leading to asymmetry also enhance star formation in MAGPI galaxies. Finally, we find a weak anti-correlation between stellar mass and asymmetry (ie. high stellar mass galaxies are less asymmetric). We conclude by discussing possible sources driving the asymmetry in the ionised gas, such as disturbances being present in the colder gas phase (either molecular or atomic) prior to the gas being ionised, and non-axisymmetric features (e.g., a bar) being present in the galactic disk. Our results highlight the complex interplay between ionised gas kinematic disturbances and physical processes involved in galaxy evolution. △ Less

Submitted 28 November, 2023; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: e.g., 20 pages, 19 figures

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 18 June, 2024; v1 submitted 6 July, 2023; originally announced July 2023.

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

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

The MAGPI Survey: Impact of environment on the total internal mass distribution of galaxies in the last 5 Gyr

Authors: Caro Derkenne, Richard M. McDermid, Adriano Poci, J. Trevor Mendel, Francesco D'Eugenio, Seyoung Jeon, Rhea-Silvia Remus, Sabine Bellstedt, Andrew J. Battisti, Joss Bland-Hawthorn, Anna Ferre-Mateu, Caroline Foster, K. E. Harborne, Claudia D. P. Lagos, Yingjie Peng, Piyush Sharda, Gauri Sharma, Sarah Sweet, Kim-Vy H. Tran, Lucas M. Valenzuela, Sam Vaughan, Emily Wisnioski, Sukyoung K. Yi

Abstract: We investigate the impact of environment on the internal mass distribution of galaxies using the Middle Ages Galaxy Properties with Integral field spectroscopy (MAGPI) survey. We use 2D resolved stellar kinematics to construct Jeans dynamical models for galaxies at mean redshift $z \sim 0.3$, corresponding to a lookback time of $3-4$ Gyr. The internal mass distribution for each galaxy is parameter… ▽ More We investigate the impact of environment on the internal mass distribution of galaxies using the Middle Ages Galaxy Properties with Integral field spectroscopy (MAGPI) survey. We use 2D resolved stellar kinematics to construct Jeans dynamical models for galaxies at mean redshift $z \sim 0.3$, corresponding to a lookback time of $3-4$ Gyr. The internal mass distribution for each galaxy is parameterised by the combined mass density slope $γ$ (baryons $+$ dark matter), which is the logarithmic change of density with radius. We use a MAGPI sample of 28 galaxies from low-to-mid density environments and compare to density slopes derived from galaxies in the high density Frontier Fields clusters in the redshift range $0.29 <z < 0.55$, corresponding to a lookback time of $\sim 5$ Gyr. We find a median density slope of $γ= -2.22 \pm 0.05$ for the MAGPI sample, which is significantly steeper than the Frontier Fields median slope ($γ= -2.01 \pm 0.04$), implying the cluster galaxies are less centrally concentrated in their mass distribution than MAGPI galaxies. We also compare to the distribution of density slopes from galaxies in Atlas3D at $z \sim 0$, because the sample probes a similar environmental range as MAGPI. The Atlas3D median total slope is $γ= -2.25 \pm 0.02$, consistent with the MAGPI median. Our results indicate environment plays a role in the internal mass distribution of galaxies, with no evolution of the slope in the last 3-4 Gyr. These results are in agreement with the predictions of cosmological simulations. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 3, July 2023, Pages 3602 - 3626

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

Dynamic Wisp Removal in JWST NIRCam Images

Authors: A. S. G. Robotham, J. C. J. D'Silva, R. A. Windhorst, R. A. Jansen, J. Summers, S. P. Driver, C. N. A. Wilmer, S. Bellstedt

Abstract: The James Webb Space Telescope (JWST) near-infrared camera (NIRCam) has been found to exhibit serious wisp-like structures in four of its eight short-wavelength detectors. The exact structure and strength of these wisps is highly variable with the position and orientation of JWST, so the use of static templates is non-optimal. Here we investigate a dynamic strategy to mitigate these wisps using lo… ▽ More The James Webb Space Telescope (JWST) near-infrared camera (NIRCam) has been found to exhibit serious wisp-like structures in four of its eight short-wavelength detectors. The exact structure and strength of these wisps is highly variable with the position and orientation of JWST, so the use of static templates is non-optimal. Here we investigate a dynamic strategy to mitigate these wisps using long-wavelength reference images. Based on a suite of experiments where we embed a worst-case scenario median-stacked wisp into wisp-free images, we define suitable parameters for our wisp removal strategy. Using this setup we re-process wisp-affected public Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS) data in the North Ecliptic Pole Time Domain Field (NEP-TDF), resulting in significant visual improvement in our detector frames and reduced noise in the final stacked images. △ Less

Submitted 13 July, 2023; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: 20 pages, 13 figures, accepted to PASP

Evolution in the orbital structure of quiescent galaxies from MAGPI, LEGA-C and SAMI surveys: direct evidence for merger-driven growth over the last 7 Gy

Authors: Francesco D'Eugenio, Arjen van der Wel, Joanna M. Piotrowska, Rachel Bezanson, Edward N. Taylor, Jesse van de Sande, William M. Baker, Eric F. Bell, Sabine Bellstedt, Joss Bland-Hawthorn, Asa F. L. Bluck, Sarah Brough, Julia J. Bryant, Matthew Colless, Luca Cortese, Scott M. Croom, Caro Derkenne, Pieter van Dokkum, Deanne Fisher, Caroline Foster, Anna Gallazzi, Anna de Graaff, Brent Groves, Josha van Houdt, Claudia del P. Lagos , et al. (15 additional authors not shown)

Abstract: We present the first study of spatially integrated higher-order stellar kinematics over cosmic time. We use deep rest-frame optical spectroscopy of quiescent galaxies at redshifts z=0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. Conservatively using a redshift-in… ▽ More We present the first study of spatially integrated higher-order stellar kinematics over cosmic time. We use deep rest-frame optical spectroscopy of quiescent galaxies at redshifts z=0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. Conservatively using a redshift-independent cut in stellar mass ($M_\star = 10^{11}\,{\rm M}_\odot$), and matching the stellar-mass distributions of our samples, we find 7 $σ$ evidence of $h_4$ increasing with cosmic time, from a median value of 0.019$\pm$0.002 at z=0.8 to 0.059$\pm$0.004 at z=0.06. Alternatively, we use a physically motivated sample selection, based on the mass distribution of the progenitors of local quiescent galaxies as inferred from numerical simulations; in this case, we find 10 $σ$ evidence. This evolution suggests that, over the last 7 Gyr, there has been a gradual decrease in the rotation-to-dispersion ratio and an increase in the radial anisotropy of the stellar velocity distribution, qualitatively consistent with accretion of gas-poor satellites. These findings demonstrate that massive galaxies continue to accrete mass and increase their dispersion support after becoming quiescent. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: 19 pages, 9 figures Accepted for publication in MNRAS

A New WISE Calibration of Stellar Mass

Authors: T. H. Jarrett, M. E. Cluver, Edward N. Taylor, Sabine Bellstedt, A. S. G Robotham, H. F. M. Yao

Abstract: We derive new empirical scaling relations between WISE mid-infrared galaxy photometry and well-determined stellar masses from SED modeling of a suite of optical-infrared photometry provided by the DR4 Catalogue of the GAMA-KiDS-VIKING survey of the southern G23 field. The mid-infrared source extraction and characterization are drawn from the WISE Extended Source Catalogue (WXSC) and the archival A… ▽ More We derive new empirical scaling relations between WISE mid-infrared galaxy photometry and well-determined stellar masses from SED modeling of a suite of optical-infrared photometry provided by the DR4 Catalogue of the GAMA-KiDS-VIKING survey of the southern G23 field. The mid-infrared source extraction and characterization are drawn from the WISE Extended Source Catalogue (WXSC) and the archival ALLWISE catalog, combining both resolved and compact galaxies in the G23 sample to a redshift of 0.15. Three scaling relations are derived: W1 3.4 micron luminosity versus stellar mass, and WISE W1-W2, W1-W3 colors versus mass-to-light ratio (sensitive to a variety of galaxy types from passive to star-forming). For each galaxy in the sample, we then derive the combined stellar mass from these scaling relations, producing Mstellar estimates with better than $\sim$25-30% accuracy for galaxies with $>$10$^{9}$ Msolar and $<$40 - 50% for lower luminosity dwarf galaxies. We also provide simple prescriptions for rest-frame corrections and estimating stellar masses using only the W1 flux and the W1-W2 color, making stellar masses more accessible to users of the WISE data. Given a redshift or distance, these new scaling relations will enable stellar mass estimates for any galaxy in the sky detected by WISE with high fidelity across a range of mass-to-light. △ Less

Submitted 14 January, 2023; originally announced January 2023.

Comments: Accepted for publication in the Astrophysical Journal (ApJ)

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

Galaxy And Mass Assembly (GAMA): The Dependence of Star Formation on Surface Brightness in Low Redshift Galaxies

Authors: S. Phillipps, S. Bellstedt, M. N. Bremer, R. De Propris, P. A. James, S. Casura, J. Liske, B. W. Holwerda

Abstract: The star formation rate in galaxies is well known to correlate with stellar mass (the `star-forming main sequence'). Here we extend this further to explore any additional dependence on galaxy surface brightness, a proxy for stellar mass surface density. We use a large sample of low redshift ($z \leq 0.08$) galaxies from the GAMA survey which have both SED derived star formation rates and photometr… ▽ More The star formation rate in galaxies is well known to correlate with stellar mass (the `star-forming main sequence'). Here we extend this further to explore any additional dependence on galaxy surface brightness, a proxy for stellar mass surface density. We use a large sample of low redshift ($z \leq 0.08$) galaxies from the GAMA survey which have both SED derived star formation rates and photometric bulge-disc decompositions, the latter providing measures of disc surface brightness and disc masses. Using two samples, one of galaxies fitted by a single component with Sérsic index below 2 and one of the discs from two-component fits, we find that once the overall mass dependence of star formation rate is accounted for, there is no evidence in either sample for a further dependence on stellar surface density. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

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

Deep Investigation of Neutral Gas Origins (DINGO): HI stacking experiments with early science data

Authors: Jonghwan Rhee, Martin Meyer, Attila Popping, Sabine Bellstedt, Simon P. Driver, Aaron S. G. Robotham, Matthew Whiting, Ivan K. Baldry, Sarah Brough, Michael J. I. Brown, John D. Bunton, Richard Dodson, Benne W. Holwerda, Andrew M. Hopkins, Bärbel S. Koribalski, Karen Lee-Waddell, Ángel R. López-Sánchez, Jon Loveday, Elizabeth Mahony, Sambit Roychowdhury, Kristóf Rozgonyi, Lister Staveley-Smith

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

Submitted 20 October, 2022; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 27 pages, 25 figures, 10 tables, accepted for publication in MNRAS

FLASH Pilot Survey: Detections of associated 21 cm HI absorption in GAMA galaxies at 0.42 < z <1.00

Authors: Renzhi Su, Elaine M. Sadler, James R. Allison, Elizabeth K. Mahony, Vanessa A. Moss, Matthew T. Whiting, Hyein Yoon, J. N. H. S. Aditya, Sabine Bellstedt, Aaron S. G. Robotham, Lilian Garratt-Smithson, Minfeng Gu, Baerbel S. Koribalski, Roberto Soria, Simon Weng

Abstract: We present the results of a search for associated 21 cm HI absorption at redshift 0.42 < z < 1.00 in radio-loud galaxies from three Galaxy And Mass Assembly (GAMA) survey fields. These observations were carried out as part of a pilot survey for the ASKAP First Large Absorption Survey in HI (FLASH). From a sample of 326 radio sources with 855.5 MHz peak flux density above 10 mJy we detected two ass… ▽ More We present the results of a search for associated 21 cm HI absorption at redshift 0.42 < z < 1.00 in radio-loud galaxies from three Galaxy And Mass Assembly (GAMA) survey fields. These observations were carried out as part of a pilot survey for the ASKAP First Large Absorption Survey in HI (FLASH). From a sample of 326 radio sources with 855.5 MHz peak flux density above 10 mJy we detected two associated HI absorption systems, in SDSS J090331+010847 at z= 0.522 and SDSS J113622+004852 at z= 0.563. Both galaxies are massive (stellar mass > 10$^{11}$ M$_{sun}$) and have optical spectra characteristic of luminous red galaxies,though SED fitting implies that SDSS J113622+004852 contains a dust-obscured starburst with SFR ~ 69 M$_{sun}$ yr$^{-1}$. The HI absorption lines have a high optical depth, with $τ_{pk}$ of 1.77 $\pm$ 0.16 for SDSS J090331+010847 (the highest value for any z > 0.1 associated system found to date) and 0.14 $\pm$ 0.01 for SDSS J113622+004852. In the redshift range probed by our ASKAP observations, the detection rate for associated HI absorption lines (with $τ_{pk}$ > 0.1 and at least 3$σ$ significance) is 2.9 (+9.7 -2.6) percent. Although the current sample is small, this rate is consistent with a trend seen in other studies for a lower detection rate of associated 21 cm HI absorption systems at higher redshift. We also searched for OH absorption lines at 0.67 < z < 1.34, but no detection was made in the 145 radio sources searched. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 26 pages, 15 figures, accepted for publication in 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

ProFuse: Physical Multi-Band Structural Decomposition of Galaxies and the Mass-Size-Age Plane

Authors: A. S. G. Robotham, S. Bellstedt, S. P. Driver

Abstract: We present the new ProFuse R package, a simultaneous spectral (ultraviolet to far infrared) and spatial structural decomposition tool that produces physical models of galaxies and their components. This combines the functionality of the recently released ProFound (for automatic source extraction), ProFit (for extended source profiling) and ProSpect (for stellar population modelling) software packa… ▽ More We present the new ProFuse R package, a simultaneous spectral (ultraviolet to far infrared) and spatial structural decomposition tool that produces physical models of galaxies and their components. This combines the functionality of the recently released ProFound (for automatic source extraction), ProFit (for extended source profiling) and ProSpect (for stellar population modelling) software packages. The key novelty of ProFuse is that it generates images using a self-consistent model for the star formation and metallicity history of the bulge and disk separately, and uses target images across a range of wavelengths to define the model likelihood and optimise our physical galaxy reconstruction. The first part of the paper explores the ProFuse approach in detail, and compares results to published structural and stellar population properties. The latter part of the paper applies ProFuse to 6,664 z < 0.06 GAMA galaxies. Using re-processed ugriZYJHKs imaging we extract structural and stellar population properties for bulges and disks in parallel. As well as producing true stellar mass based mass-size relationships, we further extend this correlation to explore the third dimensions of age and gas phase metallicity. The disks in particular demonstrate strong co-dependency between mass-size-age in a well defined plane, where at a given disk stellar mass younger disks tend to be larger. These findings are in broad agreement with work at higher redshift suggesting disks that formed earlier are physically smaller. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: Accepted for publication in MNRAS. 29 Pages, 37 Figures, 2 Tables

An empirical measurement of the Halo Mass Function from the combination of GAMA DR4, SDSS DR12, and REFLEX II data

Authors: Simon P. Driver, Aaron S. G. Robotham, Danail Obreschkow, John A. Peacock, Ivan K. Baldry, Sabine Bellstedt, Joss Bland-Hawthorn, Sarah Brough, Michelle Cluver, Benne W. Holwerda, Andrew Hopkins, Claudia Lagos, Jochen Liske, Jon Loveday, Steven Phillipps, Edward N. Taylor

Abstract: We construct the halo mass function (HMF) from the GAMA galaxy group catalogue over the mass range 10^12.7M_sol to 10^15.5M_sol, and find good agreement with the expectation from LambdaCDM. In comparison to previous studies, this result extends the mass range over which the HMF has now been measured over by an order of magnitude. We combine the GAMA DR4 HMF with similar data from the SDSS DR12 and… ▽ More We construct the halo mass function (HMF) from the GAMA galaxy group catalogue over the mass range 10^12.7M_sol to 10^15.5M_sol, and find good agreement with the expectation from LambdaCDM. In comparison to previous studies, this result extends the mass range over which the HMF has now been measured over by an order of magnitude. We combine the GAMA DR4 HMF with similar data from the SDSS DR12 and REFLEX II surveys, and fit a four-parameter Murray-Robotham-Power (MRP) function, valid at z~0.1, yielding: a density normalisation of: log10 (phi Mpc^3)=-3.96[+0.55,-0.82], a high mass turn-over of: log10(M/M_sol)=14.13[+0.43,-0.40], a low mass power law slope of: alpha=-1.68[+0.21,-0.24] , and a high mass softening parameter of: beta= 0.63[+0.25,-0.11]. If we fold in the constraint on Omega_M from Planck 2018 Cosmology, we are able to reduce these uncertainties further, but this relies on the assumption that the power-law trend can be extrapolated from 10^12.7M_sol to zero mass. Throughout, we highlight the effort needed to improve on our HMF measurement: improved halo mass estimates that do not rely on calibration to simulations; reduced halo mass uncertainties needed to mitigate the strong Eddington Bias that arises from the steepness of the HMF low mass slope; and deeper wider area spectroscopic surveys. To our halo mass limit of 10^12.7 M_sol, we are directly resolving (`seeing') 41+/-5 per cent of the total mass density, i.e. Omega_[M>12.7]=0.128+/-0.016, opening the door for the direct construction of 3D dark matter mass maps at Mpc resolution. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS

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/

HST WFC3/Grism Observations of the Candidate Ultra-High-Redshift Radio Galaxy GLEAM J0917-0012

Authors: N. Seymour, G. Drouart, G. Noirot, J. W. Broderick, R. J. Turner, S. S. Shabala, D. K. Stern, S. Bellstedt, S. Driver, L. Davies, C. A. De Breuck, J. Afonso, J. D. R. Vernet, T. J. Galvin

Abstract: We present Hubble Space Telescope Wide Field Camera 3 photometric and grism observations of the candidate ultra-high-redshift (z>7) radio galaxy, GLEAM J0917-0012. This radio source was selected due to the curvature in its 70-230 MHz, low-frequency Murchison Widefield Array radio spectrum and its faintness in K-band. Follow-up spectroscopic observations of this source with the VLA and ALMA were in… ▽ More We present Hubble Space Telescope Wide Field Camera 3 photometric and grism observations of the candidate ultra-high-redshift (z>7) radio galaxy, GLEAM J0917-0012. This radio source was selected due to the curvature in its 70-230 MHz, low-frequency Murchison Widefield Array radio spectrum and its faintness in K-band. Follow-up spectroscopic observations of this source with the VLA and ALMA were inconclusive as to its redshift. Our F105W and F0986M imaging observations detect the host of GLEAM J0917-0012 and a companion galaxy, ~one arcsec away. The G102 grism observations reveal a single weak line in each of the spectra of the host and the companion. To help identify these lines we utilised several photometric redshift techniques including template fitting to the grism spectra, fitting the UV-to-radio photometry with galaxy templates plus a synchrotron model, fitting of the UV-to-near-infrared photometry with EAZY, and fitting the radio data alone with RAiSERed. For the host of GLEAM J0917-0012 we find a line at 1.12 micron and the UV-to-radio spectral energy distribution fitting favours solutions at z~2 or z~8. While this fitting shows a weak preference for the lower redshift solution, the models from the higher redshift solution are more consistent with the strength of the spectral line. The redshift constraint by RAiSERed of z>6.5 also supports the interpretation that this line could be Lyman-alpha at z=8.21; however EAZY favours the z~2 solution. We discuss the implications of both solutions. For the companion galaxy we find a line at 0.98 micron and the spectral energy distribution fitting favours solutions at z<3 implying that the line could be the [OII]3727 doublet at z=1.63 (although the EAZY solution is z~2.6+/-0.5). Further observations are still required to unambiguously determine the redshift of this intriguing candidate ultra-high-redshift radio galaxy (abridged). △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 18 pages, 7 figures, 6 tables, accepted for publication in PASA

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.

Galaxy and Mass Assembly (GAMA): Mid-infrared properties as tracers of galaxy environment

Authors: U. Sureshkumar, A. Durkalec, A. Pollo, M. Bilicki, M. E. Cluver, S. Bellstedt, D. J. Farrow, J. Loveday, E. N. Taylor, J. Bland-Hawthorn

Abstract: We investigate how different mid-infrared (mid-IR) properties of galaxies trace the environment in which the galaxies are located. For this purpose, we first study the dependence of galaxy clustering on the absolute magnitude at 3.4 $μ$m and redshift. Then, we look into the environmental dependence of mid-IR luminosities and the galaxy properties derived from these luminosities. We also explore ho… ▽ More We investigate how different mid-infrared (mid-IR) properties of galaxies trace the environment in which the galaxies are located. For this purpose, we first study the dependence of galaxy clustering on the absolute magnitude at 3.4 $μ$m and redshift. Then, we look into the environmental dependence of mid-IR luminosities and the galaxy properties derived from these luminosities. We also explore how various infrared galaxy luminosity selections influence the galaxy clustering measurements. We use a set of W1 (3.4 $μ$m) absolute magnitude ($M_\text{W1}$) selected samples from the Galaxy and Mass Assembly (GAMA) survey matched with mid-IR properties from the Wide-field Infrared Survey Explorer (WISE) in the redshift range $0.07 \leq z < 0.43$. We compute the galaxy two-point correlation function (2pCF) and compare the clustering lengths between subsamples binned in $M_\text{W1}$ and in redshift. We also measure the marked correlation functions (MCFs) using the luminosities in the WISE W1 to W4 (3.4 to 22 $μ$m) bands as marks. Additionally, we compare the measurements of MCFs with different estimates of stellar mass and star formation rate used as marks. Finally, we check how different selections applied to the sample affect the clustering measurements. We show strong clustering dependence on the W1 absolute magnitude: galaxies brighter in the W1 band are more strongly clustered than their fainter counterparts. We also observe a lack of significant redshift dependence of clustering in the redshift range $0.07 \leq z < 0.43$. We show that although W1 and W2 bands are direct indicators of stellar mass, a galaxy sample selected based on W1 or W2 bands does not perfectly show the clustering behaviour of a stellar mass selected sample. Similar is the case with the proxy relation between W3 and W4 bands and star formation rate. △ Less

Submitted 4 January, 2023; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 18 pages, 12 figures, matches the version published in A&A

Journal ref: A&A 669, A27 (2023)

The eROSITA Final Equatorial Depth Survey (eFEDS): X-ray emission around star-forming and quiescent galaxies at $0.05<z<0.3$

Authors: Johan Comparat, Nhut Truong, Andrea Merloni, Annalisa Pillepich, Gabriele Ponti, Simon Driver, Sabine Bellstedt, Joe Liske, James Aird, Marcus Brüggen, Esra Bulbul, Luke Davies, Justo Antonio González Villalba, Antonis Georgakakis, Frank Haberl, Teng Liu, Chandreyee Maitra, Kirpal Nandra, Paola Popesso, Peter Predehl, Aaron Robotham, Mara Salvato, Jessica E. Thorne, Yi Zhang

Abstract: We aim at characterizing the hot phase of the Circum-Galactic Medium in a large sample of galaxies. We stack X-ray events from the SRG/eROSITA eFEDS survey around central galaxies in the GAMA 9hr field to construct radially projected soft X-ray luminosity profiles as a function of their stellar mass and specific star formation rate. We consider samples of quiescent (star-forming) galaxies in the s… ▽ More We aim at characterizing the hot phase of the Circum-Galactic Medium in a large sample of galaxies. We stack X-ray events from the SRG/eROSITA eFEDS survey around central galaxies in the GAMA 9hr field to construct radially projected soft X-ray luminosity profiles as a function of their stellar mass and specific star formation rate. We consider samples of quiescent (star-forming) galaxies in the stellar mass range $2\times 10^{10}$ -- $10^{12}$ M$_\odot$ ($3\times 10^9$ -- $6\times10^{11}$ M$_\odot$). For quiescent galaxies, the X-ray profiles are clearly extended throughout the available mass range; however, the measured profile is likely biased high due to projection effects, as these galaxies tend to live in dense and hot environments. For the most massive star forming samples ($\geq10^{11}$ M$_\odot$), there is a hint of detection of extended emission. For star-forming galaxies with $< 10^{11}$ M$_\odot$ the X-ray stacked profiles are compatible with unresolved sources and consistent with the expected emission from faint Active Galactic Nuclei and X-ray binaries. We measure for the first time the mean relation between average X-ray luminosity and stellar mass separately for quiescent and star-forming galaxies. High-mass ($\geq 10^{11}$ M$_\odot$) star-forming or quiescent galaxies follow the expected scaling of virialized hot haloes, while lower mass star-forming galaxies show a less prominent luminosity and a weaker dependence on stellar mass, consistent with empirical models of the weak AGN population. When comparing our results with state-of-the art numerical simulations, we find an overall consistency on large ($>80$ kpc) scales at masses $\geq 10^{11}$ M$_\odot$, but disagreement on the small scales, where brighter than observed compact cores are predicted. Simulations also do not predict the clear differentiation that we observe between quiescent and star-forming galaxies. △ Less

Submitted 17 August, 2022; v1 submitted 13 January, 2022; originally announced January 2022.

Comments: 23 pages, 11 figures, 4 tables, accepted in A&A

Journal ref: A&A 666, A156 (2022)

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

Unsupervised Domain Adaptation for Constraining Star Formation Histories

Authors: Sankalp Gilda, Antoine de Mathelin, Sabine Bellstedt, Guillaume Richard

Abstract: The prevalent paradigm of machine learning today is to use past observations to predict future ones. What if, however, we are interested in knowing the past given the present? This situation is indeed one that astronomers must contend with often. To understand the formation of our universe, we must derive the time evolution of the visible mass content of galaxies. However, to observe a complete st… ▽ More The prevalent paradigm of machine learning today is to use past observations to predict future ones. What if, however, we are interested in knowing the past given the present? This situation is indeed one that astronomers must contend with often. To understand the formation of our universe, we must derive the time evolution of the visible mass content of galaxies. However, to observe a complete star life, one would need to wait for one billion years! To overcome this difficulty, astrophysicists leverage supercomputers and evolve simulated models of galaxies till the current age of the universe, thus establishing a mapping between observed radiation and star formation histories (SFHs). Such ground-truth SFHs are lacking for actual galaxy observations, where they are usually inferred -- with often poor confidence -- from spectral energy distributions (SEDs) using Bayesian fitting methods. In this investigation, we discuss the ability of unsupervised domain adaptation to derive accurate SFHs for galaxies with simulated data as a necessary first step in developing a technique that can ultimately be applied to observational data. △ Less

Submitted 26 August, 2022; v1 submitted 28 December, 2021; originally announced December 2021.

Comments: Accepted for oral presentation at the 1st Annual AAAI Workshop on AI to Accelerate Science and Engineering (AI2ASE). Journal article to follow

Journal ref: Astronomy 2024, 3(3), 189-207

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

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

Galaxy And Mass Assembly (GAMA): $\mathbf{z \sim 0}$ Galaxy Luminosity Function down to $\mathbf{L \sim 10^{6}~L_\odot}$ via Clustering Based Redshift Inference

Authors: Geray S. Karademir, Edward N. Taylor, Chris Blake, Ivan K. Baldry, Sabine Bellstedt, Maciej Bilicki, Michael J. I. Brown, Michelle E. Cluver, Simon P. Driver, Hendrik Hildebrandt, Benne W. Holwerda, Andrew M. Hopkins, Jonathan Loveday, Steven Phillipps, Angus H. Wright

Abstract: In this study we present a new experimental design using clustering-based redshift inference to measure the evolving galaxy luminosity function (GLF) spanning 5.5 decades from $L \sim 10^{11.5}$ to $ 10^6 ~ \mathrm{L}_\odot$. We use data from the Galaxy And Mass Assembly (GAMA) survey and the Kilo-Degree Survey (KiDS). We derive redshift distributions in bins of apparent magnitude to the limits of… ▽ More In this study we present a new experimental design using clustering-based redshift inference to measure the evolving galaxy luminosity function (GLF) spanning 5.5 decades from $L \sim 10^{11.5}$ to $ 10^6 ~ \mathrm{L}_\odot$. We use data from the Galaxy And Mass Assembly (GAMA) survey and the Kilo-Degree Survey (KiDS). We derive redshift distributions in bins of apparent magnitude to the limits of the GAMA-KiDS photometric catalogue: $m_r \lesssim 23$; more than a decade in luminosity beyond the limits of the GAMA spectroscopic redshift sample via clustering-based redshift inference. This technique uses spatial cross-correlation statistics for a reference set with known redshifts (in our case, the main GAMA sample) to derive the redshift distribution for the target ensemble. For the calibration of the redshift distribution we use a simple parametrisation with an adaptive normalisation factor over the interval $0.005 < z < 0.48$ to derive the clustering redshift results. We find that the GLF has a relatively constant power-law slope $α\approx -1.2$ for $-17 \lesssim M_r \lesssim -13$, and then appears to steepen sharply for $-13 \lesssim M_r \lesssim -10$. This upturn appears to be where Globular Clusters (GCs) take over to dominate the source counts as a function of luminosity. Thus we have mapped the GLF across the full range of the $z \sim 0$ field galaxy population from the most luminous galaxies down to the GC scale. △ Less

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

Comments: 21 pages, 13 figures, accepted by MNRAS

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

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

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

Remnant Radio Galaxies Discovered in a Multi-frequency Survey

Authors: B. Quici, N. Hurley-Walker, N. Seymour, R. J. Turner, S. S. Shabala, M. Huynh, H. Andernach, A. D. Kapińska, J. D. Collier, M. Johnston-Hollitt, S. V. White, I. Prandoni, T. J. Galvin, T. Franzen, C. H. Ishwara-Chandra, S. Bellstedt, S. J. Tingay, B. M. Gaensler, A. O'Brien, J. Rogers, K. Chow, S. Driver, A. Robotham

Abstract: The remnant phase of a radio galaxy begins when the jets launched from an active galactic nucleus are switched off. To study the fraction of radio galaxies in a remnant phase, we take advantage of a $8.31$\,deg$^2$ sub-region of the GAMA~23~field which comprises of surveys covering the frequency range 0.1--9\,GHz. We present a sample of 104 radio galaxies compiled from observations conducted by th… ▽ More The remnant phase of a radio galaxy begins when the jets launched from an active galactic nucleus are switched off. To study the fraction of radio galaxies in a remnant phase, we take advantage of a $8.31$\,deg$^2$ sub-region of the GAMA~23~field which comprises of surveys covering the frequency range 0.1--9\,GHz. We present a sample of 104 radio galaxies compiled from observations conducted by the Murchison Wide-field Array (216\,MHz), the Australia Square Kilometer Array Pathfinder (887\,MHz), and the Australia Telescope Compact Array (5.5\,GHz). We adopt an `absent radio core' criterion to identify 10 radio galaxies showing no evidence for an active nucleus. We classify these as new candidate remnant radio galaxies. Seven of these objects still display compact emitting regions within the lobes at 5.5\,GHz; at this frequency the emission is short-lived, implying a recent jet switch-off. On the other hand, only three show evidence of aged lobe plasma by the presence of an ultra-steep spectrum ($α<-1.2$) and a diffuse, low surface-brightness radio morphology. The predominant fraction of young remnants is consistent with a rapid fading during the remnant phase. Within our sample of radio galaxies, our observations constrain the remnant fraction to $4\%\lesssim f_{\mathrm{rem}} \lesssim 10\%$; the lower limit comes from the limiting case in which all remnant candidates with hotspots are simply active radio galaxies with faint, undetected radio cores. Finally, we model the synchrotron spectrum arising from a hotspot to show they can persist for 5--10\,Myr at 5.5\,GHz after the jets switch off -- radio emission arising from such hotspots can therefore be expected in an appreciable fraction of genuine remnants. △ Less

Submitted 24 January, 2021; originally announced January 2021.

Deep Extragalactic VIsible Legacy Survey (DEVILS): SED Fitting in the D10-COSMOS Field and the Evolution of the Stellar Mass Function and SFR-$M_\star$ relation

Authors: Jessica E. Thorne, Aaron. S. G. Robotham, Luke J. M. Davies, Sabine Bellstedt, Simon P. Driver, Matias Bravo, Malcolm N. Bremer, Benne W. Holwerda, Andrew M. Hopkins, Claudia del P. Lagos, Steven Phillipps, Malgorzata Siudek, Edward N. Taylor, Angus H. Wright

Abstract: We present catalogues of stellar masses, star formation rates, and ancillary stellar population parameters for galaxies spanning $0<z<9$ from the Deep Extragalactic VIsible Legacy Survey (DEVILS). DEVILS is a deep spectroscopic redshift survey with very high completeness, covering several premier deep fields including COSMOS (D10). Our stellar mass and star formation rate estimates are self-consis… ▽ More We present catalogues of stellar masses, star formation rates, and ancillary stellar population parameters for galaxies spanning $0<z<9$ from the Deep Extragalactic VIsible Legacy Survey (DEVILS). DEVILS is a deep spectroscopic redshift survey with very high completeness, covering several premier deep fields including COSMOS (D10). Our stellar mass and star formation rate estimates are self-consistently derived using the spectral energy distribution (SED) modelling code ProSpect, using well-motivated parameterisations for dust attenuation, star formation histories, and metallicity evolution. We show how these improvements, and especially our physically motivated assumptions about metallicity evolution, have an appreciable systematic effect on the inferred stellar masses, at the level of $\sim$\,0.2 dex. To illustrate the scientific value of these data, we map the evolving galaxy stellar mass function (SMF) and the SFR-$M_\star$ relation for $0<z<4.25$. In agreement with past studies, we find that most of the evolution in the SMF is driven by the characteristic density parameter, with little evolution in the characteristic mass and low-mass slopes. Where the SFR-$M_\star$ relation is indistinguishable from a power-law at $z>2.6$, we see evidence of a bend in the relation at low redshifts ($z<0.45$). This suggests evolution in both the normalisation and shape of the SFR-$M_\star$ relation since cosmic noon. It is significant that we only clearly see this bend when combining our new DEVILS measurements with consistently derived values for lower redshift galaxies from the Galaxy And Mass Assembly (GAMA) survey: this shows the power of having consistent treatment for galaxies at all redshifts. △ Less

Submitted 4 May, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS. 28 pages, 20 figures

The MAGPI Survey -- science goals, design, observing strategy, early results and theoretical framework

Authors: C. Foster, J. T. Mendel, C. D. P. Lagos, E. Wisnioski, T. Yuan, F. D'Eugenio, T. M. Barone, K. E. Harborne, S. P. Vaughan, F. Schulze, R. -S. Remus, A. Gupta, F. Collacchioni, D. J. Khim, P. Taylor, R. Bassett, S. M. Croom, R. M. McDermid, A. Poci, A. J. Battisti, J. Bland-Hawthorn, S. Bellstedt, M. Colless, L. J. M. Davies, C. Derkenne , et al. (18 additional authors not shown)

Abstract: We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on ESO/VLT. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive… ▽ More We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on ESO/VLT. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive optics aided MUSE observations of fields selected from the GAMA survey, providing a wealth of publicly available ancillary multi-wavelength data. With these data, MAGPI will map the kinematic and chemical properties of stars and ionised gas for a sample of 60 massive (> 7 x 10^10 M_Sun) central galaxies at 0.25 < z <0.35 in a representative range of environments (isolated, groups and clusters). The spatial resolution delivered by MUSE with Ground Layer Adaptive Optics (GLAO, 0.6-0.8 arcsec FWHM) will facilitate a direct comparison with Integral Field Spectroscopy surveys of the nearby Universe, such as SAMI and MaNGA, and at higher redshifts using adaptive optics, e.g. SINS. In addition to the primary (central) galaxy sample, MAGPI will deliver resolved and unresolved spectra for as many as 150 satellite galaxies at 0.25 < z <0.35, as well as hundreds of emission-line sources at z < 6. This paper outlines the science goals, survey design, and observing strategy of MAGPI. We also present a first look at the MAGPI data, and the theoretical framework to which MAGPI data will be compared using the current generation of cosmological hydrodynamical simulations including EAGLE, Magneticum, HORIZON-AGN, and Illustris-TNG. Our results show that cosmological hydrodynamical simulations make discrepant predictions in the spatially resolved properties of galaxies at z ~ 0.3. MAGPI observations will place new constraints and allow for tangible improvements in galaxy formation theory. △ Less

Submitted 14 June, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 31 pages, 13 figures, PASA accepted

K-CLASH: spatially-resolving star-forming galaxies in field and cluster environments at $z \approx 0.2$-$0.6$

Authors: Alfred L. Tiley, Sam P. Vaughan, John P. Stott, Roger L. Davies, Laura J. Prichard, Andrew Bunker, Martin Bureau, Michele Cappellari, Matt Jarvis, Aaron Robotham, Luca Cortese, Sabine Bellstedt, Behzad Ansarinejad

Abstract: We present the KMOS-CLASH (K-CLASH) survey, a K-band Multi-Object Spectrograph (KMOS) survey of the spatially-resolved gas properties and kinematics of 191 (predominantly blue) H$α$-detected galaxies at $0.2 \lesssim z \lesssim 0.6$ in field and cluster environments. K-CLASH targets galaxies in four Cluster Lensing And Supernova survey with Hubble (CLASH) fields in the KMOS $IZ$-band, over $7'$ ra… ▽ More We present the KMOS-CLASH (K-CLASH) survey, a K-band Multi-Object Spectrograph (KMOS) survey of the spatially-resolved gas properties and kinematics of 191 (predominantly blue) H$α$-detected galaxies at $0.2 \lesssim z \lesssim 0.6$ in field and cluster environments. K-CLASH targets galaxies in four Cluster Lensing And Supernova survey with Hubble (CLASH) fields in the KMOS $IZ$-band, over $7'$ radius ($\approx2$-$3$ Mpc) fields-of-view. K-CLASH aims to study the transition of star-forming galaxies from turbulent, highly star-forming disc-like and peculiar systems at $z\approx1$-$3$, to the comparatively quiescent, ordered late-type galaxies at $z\approx0$, and to examine the role of clusters in the build-up of the red sequence since $z\approx1$. In this paper, we describe the K-CLASH survey, present the sample, and provide an overview of the K-CLASH galaxy properties. We demonstrate that our sample comprises star-forming galaxies typical of their stellar masses and epochs, residing both in field and cluster environments. We conclude K-CLASH provides an ideal sample to bridge the gap between existing large integral-field spectroscopy surveys at higher and lower redshifts. We find that star-forming K-CLASH cluster galaxies at intermediate redshifts have systematically lower stellar masses than their star-forming counterparts in the field, hinting at possible "downsizing" scenarios of galaxy growth in clusters at these epochs. We measure no difference between the star-formation rates of H$α$-detected, star-forming galaxies in either environment after accounting for stellar mass, suggesting that cluster quenching occurs very rapidly during the epochs probed by K-CLASH, or that star-forming K-CLASH galaxies in clusters have only recently arrived there, with insufficient time elapsed for quenching to have occured. △ Less

Submitted 27 May, 2020; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS, 31 pages, 15 figures, 2 tables; updated to correct typos in Section 5.2

Galaxy And Mass Assembly (GAMA): A forensic SED reconstruction of the cosmic star formation history and metallicity evolution by galaxy type

Authors: Sabine Bellstedt, Aaron S. G. Robotham, Simon P. Driver, Jessica E. Thorne, Luke J. M. Davies, Claudia del P. Lagos, Adam R. H. Stevens, Edward N. Taylor, Ivan K. Baldry, Amanda J. Moffett, Andrew M. Hopkins, Steven Phillipps

Abstract: We apply the spectral energy distribution (SED) fitting code ProSpect to multiwavelength imaging for $\sim$7,000 galaxies from the GAMA survey at $z<0.06$, in order to extract their star formation histories. We combine a parametric description of the star formation history with a closed-box evolution of metallicity where the present-day gas-phase metallicity of the galaxy is a free parameter. We s… ▽ More We apply the spectral energy distribution (SED) fitting code ProSpect to multiwavelength imaging for $\sim$7,000 galaxies from the GAMA survey at $z<0.06$, in order to extract their star formation histories. We combine a parametric description of the star formation history with a closed-box evolution of metallicity where the present-day gas-phase metallicity of the galaxy is a free parameter. We show with this approach that we are able to recover the observationally determined cosmic star formation history (CSFH), an indication that stars are being formed in the correct epoch of the Universe, on average, for the manner in which we are conducting SED fitting. We also show the contribution to the CSFH of galaxies of different present-day visual morphologies, and stellar masses. Our analysis suggests that half of the mass in present-day elliptical galaxies was in place 11 Gyr ago. In other morphological types, the stellar mass formed later, up to 6 Gyr ago for present-day irregular galaxies. Similarly, the most massive galaxies in our sample were shown to have formed half their stellar mass by 11 Gyr ago, whereas the least massive galaxies reached this stage as late as 4 Gyr ago (the well-known effect of "galaxy downsizing"). Finally, our metallicity approach allows us to follow the average evolution in gas-phase metallicity for populations of galaxies, and extract the evolution of the cosmic metal mass density in stars and in gas, producing results in broad agreement with independent, higher-redshift observations of metal densities in the Universe. △ Less

Submitted 7 October, 2020; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: 24 pages, 16 figures, accepted for publication by MNRAS

Journal ref: MNRAS, 2020, 498, 5581

Galaxy And Mass Assembly (GAMA): Assimilation of KiDS into the GAMA database

Authors: Sabine Bellstedt, Simon P. Driver, Aaron S. G. Robotham, Luke J. M. Davies, Cameron R. J. Bogue, Robin H. W. Cook, Abdolhosein Hashemizadeh, Soheil Koushan, Edward N. Taylor, Jessica E. Thorne, Ryan J. Turner, Angus H. Wright

Abstract: The Galaxy And Mass Assembly Survey (GAMA) covers five fields with highly complete spectroscopic coverage ($>95$ per cent) to intermediate depths ($r<19.8$ or $i < 19.0$ mag), and collectively spans 250 square degrees of Equatorial or Southern sky. Four of the GAMA fields (G09, G12, G15 and G23) reside in the ESO VST KiDS and ESO VISTA VIKING survey footprints, which combined with our GALEX, WISE… ▽ More The Galaxy And Mass Assembly Survey (GAMA) covers five fields with highly complete spectroscopic coverage ($>95$ per cent) to intermediate depths ($r<19.8$ or $i < 19.0$ mag), and collectively spans 250 square degrees of Equatorial or Southern sky. Four of the GAMA fields (G09, G12, G15 and G23) reside in the ESO VST KiDS and ESO VISTA VIKING survey footprints, which combined with our GALEX, WISE and Herschel data provide deep uniform imaging in the $FUV\,NUV\,ugriZYJHK_s\,W1\,W2\,W3\,W4\,P100\,P160\,S250\,S350\,S500$ bands. Following the release of KiDS DR4, we describe the process by which we ingest the KiDS data into GAMA (replacing the SDSS data previously used for G09, G12 and G15), and redefine our core optical and near-IR catalogues to provide a complete and homogeneous dataset. The source extraction and analysis is based on the new ProFound image analysis package, providing matched-segment photometry across all bands. The data are classified into stars, galaxies, artefacts, and ambiguous objects, and objects are linked to the GAMA spectroscopic target catalogue. Additionally, a new technique is employed utilising ProFound to extract photometry in the unresolved MIR-FIR regime. The catalogues including the full FUV-FIR photometry are described and will be fully available as part of GAMA DR4. They are intended for both standalone science, selection for targeted follow-up with 4MOST, as well as an accompaniment to the upcoming and ongoing radio arrays now studying the GAMA $23^h$ field. △ Less

Submitted 22 May, 2020; originally announced May 2020.

Comments: 23 pages, 18 figures, accepted for publication in MNRAS

The assembly history of the nearest S0 galaxy NGC 3115 from its kinematics out to six half-light radii

Authors: Arianna Dolfi, Duncan A. Forbes, Warrick J. Couch, Anna Ferre-Mateu, Sabine Bellstedt, Kenji Bekki, Jonathan Diaz, Aaron J. Romanowsky, Jean P. Brodie

Abstract: Using new and archival data, we study the kinematic properties of the nearest field S0 galaxy, NGC 3115, out to $\sim6.5$ half-light radii ($R_\mathrm{e}$) from its stars (integrated starlight), globular clusters (GCs) and planetary nebulae (PNe). We find evidence of three kinematic regions with an inner transition at $\sim0.2\ R_\mathrm{e}$ from a dispersion-dominated bulge ($V_\mathrm{rot}/σ<1$)… ▽ More Using new and archival data, we study the kinematic properties of the nearest field S0 galaxy, NGC 3115, out to $\sim6.5$ half-light radii ($R_\mathrm{e}$) from its stars (integrated starlight), globular clusters (GCs) and planetary nebulae (PNe). We find evidence of three kinematic regions with an inner transition at $\sim0.2\ R_\mathrm{e}$ from a dispersion-dominated bulge ($V_\mathrm{rot}/σ<1$) to a fast-rotating disk ($V_\mathrm{rot}/σ>1$), and then an additional transition from the disk to a slowly rotating spheroid at $\sim2-2.5\, R_\mathrm{e}$, as traced by the red GCs and PNe (and possibly by the blue GCs beyond $\sim5\, R_\mathrm{e}$). From comparison with simulations, we propose an assembly history in which the original progenitor spiral galaxy undergoes a gas-rich minor merger that results in the embedded kinematically cold disk that we see today in NGC 3115. At a later stage, dwarf galaxies, in mini mergers (mass-ratio $<$ 1:10), were accreted building-up the outer slowly rotating spheroid, with the central disk kinematics largely unaltered. Additionally, we report new spectroscopic observations of a sample of ultra-compact dwarfs (UCDs) around NGC 3115 with the Keck/KCWI instrument. We find that five UCDs are inconsistent with the general rotation field of the GCs, suggesting an \textit{ex-situ} origin for these objects, i.e. perhaps the remnants of tidally stripped dwarfs. A further seven UCDs follow the GC rotation pattern, suggesting an \textit{in-situ} origin and, possibly a GC-like nature. △ Less

Submitted 24 May, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: 22 pages (including 3 pages of Appendix material), 14 figures, published in MNRAS

Journal ref: MNRAS 495, 1321-1339 (2020)