Search

We use deep spectroscopy from the SAMI Galaxy Survey to explore the precision of the fundamental plane of early-type galaxies (FP) as a distance indicator for future single-fibre spectroscopy surveys. We study the optimal trade-off between sample size and signal-to-noise ratio (SNR), and investigate which additional observables can be used to construct hyperplanes with smaller intrinsic scatter than the FP. We add increasing levels of random noise (parametrised as effective exposure time) to the SAMI spectra to study the effect of increasing measurement uncertainties on the FP-and hyperplane-inferred distances. We find that, using direct-fit methods, the values of the FP and hyperplane best-fit coefficients depend on the spectral SNR, and reach asymptotic values for a mean SNR=40 {\AA}$^{-1}$. As additional variables for the FP we consider three stellar-population observables: light-weighted age, stellar mass-to-light ratio and a novel combination of Lick indices (I$_{\rm age}$). For a SNR=45 {\AA}$^{-1}$ (equivalent to 1-hour exposure on a 4-m telescope), all three hyperplanes outperform the FP as distance indicators. Being an empirical spectral index, I$_{\rm age}$ avoids the model-dependent uncertainties and bias underlying age and mass-to-light ratio measurements, yet yields a 10 per cent reduction of the median distance uncertainty compared to the FP. We also find that, as a by-product, the Iage hyperplane removes most of the reported environment bias of the FP. After accounting for the different signal-to-noise ratio, these conclusions also apply to a 50 times larger sample from SDSS-III. However, in this case, only age removes the environment bias., Comment: 24 pages, 18 figures, accepted for publication in MNRAS

Empirical correlations connecting starlight to galaxy dynamics (e.g., the fundamental plane (FP) of elliptical/quiescent galaxies and the Tully--Fisher relation of spiral/star-forming galaxies) provide cosmology-independent distance estimation and are central to local Universe cosmology. In this work, we introduce the mass hyperplane (MH), which is the stellar-to-dynamical mass relation $(M_\star/M_\mathrm{dyn})$ recast as a linear distance indicator. Building on recent FP studies, we show that both star-forming and quiescent galaxies follow the same empirical MH, then use this to measure the peculiar velocities (PVs) for a sample of 2496 galaxies at $z<0.12$ from GAMA. The limiting precision of MH-derived distance/PV estimates is set by the intrinsic scatter in size, which we find to be $\approx$0.1~dex for both quiescent and star-forming galaxies (when modeled independently) and $\approx$0.11~dex when all galaxies are modeled together; showing that the MH is as good as the FP. To empirically validate our framework and distance/PV estimates, we compare the inferred distances to groups as derived using either quiescent or star-forming galaxies. A good agreement is obtained with no discernible bias or offset, having a scatter of $\approx$0.05~dex $\approx$12\% in distance. Further, we compare our PV measurements for the quiescent galaxies to the previous PV measurements of the galaxies in common between GAMA and the Sloan Digital Sky Survey (SDSS), which shows similarly good agreement. Finally, we provide comparisons of PV measurements made with the FP and the MH, then discuss possible improvements in the context of upcoming surveys such as the 4MOST Hemisphere Survey (4HS)., Comment: Accepted: 15th May 2024

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$\alpha$ 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., Comment: 19 pages, 15 figures, 4 tables, accepted for publication in MNRAS

The recent UNCOVER survey with the James Webb Space Telescope (JWST) exploits the nearby cluster Abell 2744 to create the deepest view of our universe to date by leveraging strong gravitational lensing. In this work, we perform photometric fitting of more than 50,000 robustly detected sources out to $z \sim 15$. We show the redshift evolution of stellar ages, star formation rates, and rest-frame colors across the full range of $0.2 \lesssim z \lesssim 15$. The galaxy properties are inferred using the Prospector Bayesian inference framework using informative Prospector-$\beta$ priors on masses and star formation histories to produce joint redshift and stellar population posteriors, and additionally lensing magnification is performed on-the-fly to ensure consistency with the scale-dependent priors. We show that this approach produces excellent photometric redshifts with $\sigma_{\rm NMAD} \sim 0.03$, of a similar quality to the established photometric redshift code EAzY. In line with the open-source scientific objective of the Treasury survey, we publicly release the stellar population catalog with this paper, derived from the photometric catalog adapting aperture sizes based on source profiles. This release includes posterior moments, maximum-likelihood spectra, star-formation histories, and full posterior distributions, offering a rich data set to explore the processes governing galaxy formation and evolution over a parameter space now accessible by JWST., Comment: Corrected typos: Eq.1 should've been (1-kappa)^2, and the lens maps are normalized to D_ds/D_s=1. These errors were only in the writing; no data products or results were affected. The SPS catalogs are accessible via the UNCOVER survey webpage: https://jwst-uncover.github.io/DR2.html#SPSCatalogs, with a copy deposited to Zenodo: https://doi.org/10.5281/zenodo.8401181

In this empirical work, we aim to quantify the systematic uncertainties in stellar mass $(M_\star)$ estimates made from spectral energy distribution (SED) fitting through stellar population synthesis (SPS), for galaxies in the local Universe, by using the dynamical mass $(M_\text{dyn})$ estimator as an SED-independent check on stellar mass. We first construct a statistical model of the high dimensional space of galaxy properties; size $(R_e)$, velocity dispersion $(\sigma_e)$, surface brightness $(I_e)$, mass-to-light ratio $(M_\star/L)$, rest-frame colour, S\'ersic index $(n)$ and dynamical mass $(M_\text{dyn})$; accounting for selection effects and covariant errors. We disentangle the correlations among galaxy properties and find that the variation in $M_\star/M_\text{dyn}$ is driven by $\sigma_e$, S\'ersic index and colour. We use these parameters to calibrate an SED-independent $M_\star$ estimator, $\hat{M}_\star$. We find the random scatter of the relation $M_\star-\hat{M}_\star$ to be $0.108\text{dex}$ and $0.147\text{dex}$ for quiescent and star-forming galaxies respectively. Finally, we inspect the residuals as a function of SPS parameters (dust, age, metallicity, star formation rate) and spectral indices (H$\alpha$, H$\delta$, $D_n4000)$. For quiescent galaxies, $\sim65\%$ of the scatter can be explained by the uncertainty in SPS parameters, with dust and age being the largest sources of uncertainty. For star-forming galaxies, while age and metallicity are the leading factors, SPS parameters account for only $\sim13\%$ of the scatter. These results leave us with remaining unmodelled scatters of $0.055\text{dex}$ and $0.122\text{dex}$ for quiescent and star-forming galaxies respectively. This can be interpreted as a conservative limit on the precision in $M_\star$ that can be achieved via simple SPS-modelling., Comment: Accepted for publication in the Astrophysical Journal on 14 June 2023

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 $\sigma$ 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 $\sigma$ 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., Comment: 19 pages, 9 figures Accepted for publication in MNRAS

We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass $M_\star > 10^{10.5} {\rm M}_\odot$) 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. We find that at all redshifts where we have large enough samples, $h_4$ anti-correlates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anti-correlation with rotation-to-dispersion ratio, we also find a correlation between $h_4$ and $M_\star$, potentially connected to the assembly history of galaxies. In contrast, after controlling for mass, we find no evidence of independent correlation between $h_4$ and aperture velocity dispersion or galaxy size. These results hold for both star-forming and quiescent galaxies. For quiescent galaxies, $h_4$ also correlates with projected shape, even after controlling for the rotation-to-dispersion ratio. At any given redshift, star-forming galaxies have lower $h_4$ compared to quiescent galaxies, highlighting the link between kinematic structure and star-forming activity., Comment: 26 pages, 15 figures Accepted for publication in MNRAS

Context. Strong lenses are a biased subset of the general population of galaxies. Aims. The goal of this work is to quantify how lens galaxies and lensed sources differ from their parent distribution, namely the strong lensing bias. Methods. We first studied how the strong lensing cross-section varies as a function of lens and source properties. Then, we simulated strong lensing surveys with data similar to that expected for Euclid and measured the strong lensing bias in different scenarios. We focused particularly on two quantities: the stellar population synthesis mismatch parameter, $\alpha_{sps}$, defined as the ratio between the true stellar mass of a galaxy and the stellar mass obtained from photometry, and the central dark matter mass at fixed stellar mass and size. Results. Strong lens galaxies are biased towards larger stellar masses, smaller half-mass radii and larger dark matter masses. The amplitude of the bias depends on the intrinsic scatter in the mass-related parameters of the galaxy population and on the completeness in Einstein radius of the lens sample. For values of the scatter that are consistent with observed scaling relations and a minimum detectable Einstein radius of $0.5''$, the strong lensing bias in $\alpha_{sps}$ is $10\%$, while that in the central dark matter mass is $5\%$. The bias has little dependence on the properties of the source population: samples of galaxy-galaxy lenses and galaxy-quasar lenses that probe the same Einstein radius distribution are biased in a very similar way. Conclusions. Given current uncertainties, strong lensing observations can be used directly to improve our current knowledge of the inner structure of galaxies, without the need to correct for selection effects. Time-delay measurements of $H_0$ from lensed quasars can take advantage of prior information obtained from galaxy-galaxy lenses with similar Einstein radii., Comment: Published on Astronomy & Astrophysics. A two-minute summary video of this paper is available at https://youtu.be/UmS9jRHTmZU

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., Comment: Accepted for publication in the Astrophysical Journal (ApJ)

In November 2022, the James Webb Space Telescope (JWST) returned deep near-infrared images of Abell 2744 -- a powerful lensing cluster capable of magnifying distant, incipient galaxies beyond it. Together with the existing Hubble Space Telescope (HST) imaging, this publicly available dataset opens a fundamentally new discovery space to understand the remaining mysteries of the formation and evolution of galaxies across cosmic time. In this work, we detect and measure some 60,000 objects across the 49 arcmin$^2$ JWST footprint down to a $5\,\sigma$ limiting magnitude of $\sim$30 mag in 0.32" apertures. Photometry is performed using circular apertures on images matched to the point spread function of the reddest NIRCam broad band, F444W, and cleaned of bright cluster galaxies and the related intra-cluster light. To give an impression of the photometric performance, we measure photometric redshifts and achieve a $\sigma_{\rm NMAD}\approx0.03$ based on known, but relatively small, spectroscopic samples. With this paper, we publicly release our HST and JWST PSF-matched photometric catalog with optimally assigned aperture sizes for easy use, along with single aperture catalogs, photometric redshifts, rest-frame colors, and individual magnification estimates. These catalogs will set the stage for efficient and deep spectroscopic follow-up of some of the first JWST-selected samples in Summer 2023., Comment: 28 pages, 19 figures, resubmitted to ApJS following significant data product improvements. Comments welcome. Catalogs can be accessed at https://jwst-uncover.github.io/DR2.html#PhotometricCatalogs

We have developed a pipeline called \mentari to generate the far-ultraviolet to far-infrared spectral energy distribution (SED) of galaxies from the \dustysage semi-analytic galaxy formation model (SAM). \dustysage incorporates dust-related processes directly on top of the basic ingredients of galaxy formation like gas infall, cooling, star formation, feedback, and mergers. We derive a physically motivated attenuation model from the computed dust properties in \dustysage, so each galaxy has a self-consistent set of attenuation parameters based on the complicated dust physics that occurred across the galaxy's assembly history. Then, we explore several dust emission templates to produce infrared spectra. Our results show that a physically-motivated attenuation model is better for obtaining a consistent multi-wavelength description of galaxy formation and evolution, compared to using a constant attenuation. We compare our predictions with a compilation of observations and find that the fiducial model is in reasonable agreement with: (i) the observed $z=0$ luminosity functions from the far-ultraviolet to far-infrared simultaneously, and hence (ii) the local cosmic SED in the same range, (iii) the rest-frame K-band luminosity function across $0 < z < 3$, and (iv) the rest-frame far-ultraviolet luminosity function across $0 < z < 1$. Our model underproduces the far-ultraviolet emission at $z=2$ and $z=3$, which can be improved by altering the AGN feedback and dust processes in \dustysage. However, this combination thus worses the agreement at $z=0$, which suggests that more detailed treatment of such processes is required., Comment: 20 pages, 14 figures, accepted for publication in MNRAS

In this paper we describe the survey design for the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program, which executed its early imaging component in November 2022. The UNCOVER survey includes ultradeep ($\sim29-30\mathrm{AB}$) imaging of $\sim$45 arcmin$^2$ on and around the well-studied Abell 2744 galaxy cluster at $z=0.308$ and will follow-up ${\sim}500$ galaxies with extremely deep low-resolution spectroscopy with the NIRSpec/PRISM during the summer of 2023. We describe the science goals, survey design, target selection, and planned data releases. We also present and characterize the depths of the first NIRCam imaging mosaic, highlighting previously unparalleled resolved and ultradeep 2-4 micron imaging of known objects in the field. The UNCOVER primary NIRCam mosaic spans 28.8 arcmin$^2$ in seven filters (F115W, F150W, F200W, F277W, F356W, F410M, F444W) and 16.8 arcmin$^2$ in our NIRISS parallel (F115W, F150W, F200W, F356W, and F444W). To maximize early community use of the Treasury data set, we publicly release full reduced mosaics of public JWST imaging including 45 arcmin$^2$ NIRCam and 17 arcmin$^2$ NIRISS mosaics on and around the Abell 2744 cluster, including the Hubble Frontier Field primary and parallel footprints., Comment: 17 pages, 8 figures, 4 tables, submitted to ApJ, comments welcome (v2 with full author list in metadata)

The WALLABY pilot survey has been conducted using the Australian SKA Pathfinder (ASKAP). The integrated 21-cm HI line spectra are formed in a very different manner compared to usual single-dish spectra Tully-Fisher measurements. It is thus extremely important to ensure that slight differences (e.g. biases due to missing flux) are quantified and understood in order to maximise the use of the large amount of data becoming available soon. This article is based on four fields for which the data are scientifically interesting by themselves. The pilot data discussed here consist of 614 galaxy spectra at a rest wavelength of 21cm. Of these spectra, 472 are of high enough quality to be used to potentially derive distances using the Tully-Fisher relation. We further restrict the sample to the 251 galaxies whose inclination is sufficiently close to edge-on. For these, we derive Tully-Fisher distances using the deprojected WALLABY velocity widths combined with infrared (WISE W1) magnitudes. The resulting Tully-Fisher distances for the Eridanus, Hydra, Norma and NGC 4636 clusters are 21.5, 53.5, 69.4 and 23.0 Mpc respectively, with uncertainties of 5--10\%, which are better or equivalent to the ones obtained in studies using data obtained with giant single dish telescopes. The pilot survey data show the benefits of WALLABY over previous giant single-dish telescope surveys. WALLABY is expected to detect around half a million galaxies with a mean redshift of $z = 0.05 (200 Mpc)$. This study suggests that about 200,000 Tully-Fisher distances might result from the survey., Comment: 21 pages, 12 Figures, 6 Tables, accepted for publication in the MNRAS

We develop a framework for using clustering-based redshift inference (cluster-$z$) to measure the evolving galaxy luminosity function (GLF) and galaxy stellar mass function (GSMF) using WISE W1 ($3.4\mu m$) mid-infrared photometry and positions. We use multiple reference sets from the Galaxy And Mass Assembly (GAMA) survey, Sloan Digital Sky Survey (SDSS) and Baryon Oscillation Spectroscopic Survey (BOSS). Combining the resulting cluster-$z$s allows us to enlarge the study area, and by accounting for the specific properties of each reference set, making best use of each reference set to produce the best overall result. Thus we are able to measure the GLF and GSMF over $\sim 7500\, \mathrm{deg}^2 $ of the Northern Galactic Cap (NGC) up to $z<0.6$. Our method can easily be adapted for new studies with fainter magnitudes, which pose difficulties for the derivation of photo-$z$s. The measurement of the GSMF is currently limited by the models for k-corrections and mass-to-light ratios, rather than more complicated effects tied to the evolution of the differential galaxy bias. With better statistics in future surveys this technique is a strong candidate for studies with new emerging data from, e.g. the Vera C. Rubin Observatory, the Euclid mission or the Nancy Grace Roman Space Telescope., Comment: 17 pages, 16 figures, accepted by MNRAS

We derive single S\'ersic fits and bulge-disk decompositions for 13096 galaxies at redshifts z < 0.08 in the GAMA II equatorial survey regions in the Kilo-Degree Survey (KiDS) g, r and i bands. The surface brightness fitting is performed using the Bayesian two-dimensional profile fitting code ProFit. We fit three models to each galaxy in each band independently with a fully automated Markov-chain Monte Carlo analysis: a single S\'ersic model, a S\'ersic plus exponential and a point source plus exponential. After fitting the galaxies, we perform model selection and flag galaxies for which none of our models are appropriate (mainly mergers/Irregular galaxies). The fit quality is assessed by visual inspections, comparison to previous works, comparison of independent fits of galaxies in the overlap regions between KiDS tiles and bespoke simulations. The latter two are also used for a detailed investigation of systematic error sources. We find that our fit results are robust across various galaxy types and image qualities with minimal biases. Errors given by the MCMC underestimate the true errors typically by factors 2-3. Automated model selection criteria are accurate to > 90 % as calibrated by visual inspection of a subsample of galaxies. We also present g-r component colours and the corresponding colour-magnitude diagram, consistent with previous works despite our increased fit flexibility. Such reliable structural parameters for the components of a diverse sample of galaxies across multiple bands will be integral to various studies of galaxy properties and evolution. All results are integrated into the GAMA database., Comment: 36 pages, 33 figures. Accepted for publication in MNRAS

It has been postulated that the delivery of radiotherapy at ultra-high dose rates ("FLASH") reduces normal tissue toxicities by depleting them of oxygen. The fraction of normal tissue and cancer cells surviving radiotherapy depends on dose and oxygen levels in an exponential manner and even a very small fraction of tissue at low oxygen levels can determine radiotherapy response. The effect of FLASH on radiation-induced normal and tumour tissue cell killing was studied by simulating oxygen diffusion, metabolism, and radiolytic oxygen depletion over domains with simulated capillary architectures. Two architectural models were used: 1.) randomly distributed capillaries and 2.) capillaries forming a regular square lattice array. The resulting oxygen partial pressure distribution histograms were used to simulate normal and tumour tissue cell survival using the linear quadratic model of cell survival, modified to incorporate oxygen-enhancement ratio (OER) effects. Tumour cell survival was found to be increased by FLASH as compared to conventional radiotherapy, with a 0-1 order of magnitude increase for expected levels of tumour hypoxia, depending on the relative magnitudes of radiolytic oxygen depletion and tissue oxygen metabolism. Interestingly, for the random capillary model, the impact of FLASH on well-oxygenated (normal) tissues was found to be much greater, with an estimated increase in cell survival by up to 10 orders of magnitude, even though reductions in mean tissue partial pressure were modest, less than 7 mmHg for the parameter values studied. The presence of very small nearly hypoxic regions in otherwise well-perfused normal tissues with high mean oxygen levels resulted in a greater proportional sparing of normal tissue than tumour cells during FLASH irradiation, possibly explaining empirical normal tissue sparing and iso-tumour control results., Comment: Published version

There is evidence that the population of cells that experience fluctuating oxygen levels are more radioresistant than chronically hypoxic ones and hence, this population may determine radiotherapy (RT) response, in particular for hypofractionated RT, where reoxygenation may not be as prominent. A considerable effort has been devoted to examining the impact of hypoxia on hypofractionated RT; however, much less attention has been paid to cyclic hypoxia specifically and the role its kinetics may play in determining the efficacy of these treatments. Here, a simple mathematical model of cyclic hypoxia and fractionation effects was worked out to quantify this. Cancer clonogen survival fraction was estimated using the linear quadratic model, modified to account for oxygen enhancement effects and inter-fraction tissue oxygen kinetics. The resulting survival fraction formula was used to derive an expression for the iso-survival biologically effective dose, $\bedeff$. These were computed for some common extra-cranial hypofractionated radiotherapy regimens. Using relevant literature parameter values, inter-fraction fluctuations in oxygenation were found to result in an added 1-2 logs of clonogen survival fraction in going from five fractions to one for the same nominal biologically effective dose. $\bedeff$'s for most ultra-hypofractionated (five or fewer fractions) regimens in a given tumour site are similar in magnitude, suggesting iso-efficacy for common fractionation schedules. Although significant, the loss of cell-killing with increasing hypofractionation is not nearly as large as previous estimates based on the assumption of complete reoxygenation between fractions. Most ultra-hypofractionated regimens currently in place offer sufficiently high doses to counter this loss of cell killing, although care should be taken in implementing single-fraction regimens., Comment: Published version

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., Comment: Accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS. GAMA Data Release 4 is available at: http://www.gama-survey.org/dr4/

Understanding the connection between nuclear activity and galaxy environment remains critical in constraining models of galaxy evolution. By exploiting extensive catalogued data from the Galaxy and Mass Assembly (GAMA) survey, we identify a representative sample of 205 quasars at 0.1 < z < 0.35 and establish a comparison sample of galaxies, closely matched to the quasar sample in terms of both stellar mass and redshift. On scales <1 Mpc, the galaxy number counts and group membership of quasars appear entirely consistent with those of the matched galaxy sample. Despite this, we find that quasars are ~1.5 times more likely to be classified as the group center, indicating a potential link between quasar activity and cold gas flows or galaxy interactions associated with rich group environments. On scales of ~a few Mpc, the clustering strength of both samples are statistically consistent and beyond 10 Mpc we find no evidence that quasars trace large scale structures any more than the galaxy control sample. Both populations are found to prefer intermediate-density sheets and filaments to either very high- or very low- density environments. This weak dependence of quasar activity on galaxy environment supports a paradigm in which quasars represent a phase in the lifetime of all massive galaxies and in which secular processes and a group-centric location are the dominant trigger of quasars at low redshift., Comment: 26 pages, 20 figures, accepted to ApJ

Deep near-infrared photometric surveys are efficient in identifying high-redshift galaxies, however they can be prone to systematic errors in photometric redshift. This is particularly salient when there is limited sampling of key spectral features of a galaxy's spectral energy distribution (SED), such as for quiescent galaxies where the expected age-sensitive Balmer/4000 A break enter the $K$-band at $z>4$. With single filter sampling of this spectral feature, degeneracies between SED models and redshift emerge. A potential solution to this comes from splitting the $K$-band into multiple filters. We use simulations to show an optimal solution is to add two medium-band filters, $K_\mathrm{blue}$ ($\lambda_\mathrm{cen}$=2.06 $\mu$m, $\Delta\lambda$=0.25 $\mu$m) and $K_\mathrm{red}$ ($\lambda_\mathrm{cen}$=2.31 $\mu$m, $\Delta\lambda$=0.27 $\mu$m), that are complementary to the existing $K_\mathrm{s}$ filter. We test the impact of the $K$-band filters with simulated catalogues comprised of galaxies with varying ages and signal-to-noise. The results suggest that the $K$-band filters do improve photometric redshift constraints on $z>4$ quiescent galaxies, increasing precision and reducing outliers by up to 90$\%$. We find that the impact from the $K$-band filters depends on the signal-to-noise, the redshift and the SED of the galaxy. The filters we designed were built and used to conduct a pilot of the FLAMINGOS-2 Extra-galactic Near-Infrared $K$-band Split (FENIKS) survey. While no new $z>4$ quiescent galaxies are identified in the limited area pilot, the $K_\mathrm{blue}$ and $K_\mathrm{red}$ filters indicate strong Balmer/4000 A breaks in existing candidates. Additionally we identify galaxies with strong nebular emission lines, for which the $K$-band filters increase photometric redshift precision and in some cases indicate extreme star-formation., Comment: 27 pages, 18 figures. Accepted for publication in the Astronomical Journal

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 $\alpha \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., Comment: 21 pages, 13 figures, accepted by MNRAS

Dose painting of hypoxic tumour sub-volumes using positron-emission tomography (PET) has been shown to improve tumour control in silico in several sites. Pancreatic cancer presents a more stringent challenge, given its proximity to critical organs-at-risk (OARs) and anatomic motion. A radiobiological model was developed to estimate clonogen survival fraction (SF), using 18F-fluoroazomycin arabinoside PET (FAZA PET) images from ten patients with pancreatic cancer to quantify oxygen enhancement effects. For each patient, four simulated five-fraction stereotactic body radiotherapy (SBRT) plans were generated: 1) a standard SBRT plan aiming to cover the planning target volume with 40 Gy, 2) dose painting plans delivering escalated doses to FAZA-avid hypoxic sub-volumes, 3) dose painting plans with simulated spacer separating the duodenum and pancreatic head, and 4), plans with integrated boosts to geometric contractions of the tumour (GTV). All plans saturated at least one OAR dose limit. SF was calculated for each plan and sensitivity of SF to simulated hypoxia quantification errors was evaluated. Dose painting resulted in a 55% reduction in SF as compared to standard SBRT; 78% with spacer. Integrated boosts to hypoxia-blind geometric contractions resulted in a 41% reduction in SF. The reduction in SF for dose-painting plans persisted for all hypoxia quantification parameters studied, including registration and rigid motion errors that resulted in shifts and rotations of the GTV and hypoxic sub-volumes by as much as 1 cm and 10 degrees. Although proximity to OARs ultimately limited dose escalation, with estimated SFs (~10^-5) well above levels required to completely ablate a ~10 cm^3 tumour, dose painting robustly reduced clonogen survival when accounting for expected treatment and imaging uncertainties and thus, may improve local response and associated morbidity.

We have entered a new era where integral-field spectroscopic surveys of galaxies are sufficiently large to adequately sample large-scale structure over a cosmologically significant volume. This was the primary design goal of the SAMI Galaxy Survey. Here, in Data Release 3 (DR3), we release data for the full sample of 3068 unique galaxies observed. This includes the SAMI cluster sample of 888 unique galaxies for the first time. For each galaxy, there are two primary spectral cubes covering the blue (370-570nm) and red (630-740nm) optical wavelength ranges at spectral resolving power of R=1808 and 4304 respectively. For each primary cube, we also provide three spatially binned spectral cubes and a set of standardized aperture spectra. For each galaxy, we include complete 2D maps from parameterized fitting to the emission-line and absorption-line spectral data. These maps provide information on the gas ionization and kinematics, stellar kinematics and populations, and more. All data are available online through Australian Astronomical Optics (AAO) Data Central., Comment: Accepted for publication in MNRAS. 27 pages, 21 figures. Data available at https://datacentral.org.au/ . See also http://sami-survey.org/

We analyse the first measurements from precision weak lensing (PWL): a new methodology for measuring individual galaxy-galaxy weak lensing through velocity information. Our goal is to understand the observed shear distribution from PWL, which is broader than can be explained by the statistical measurement errors. We identify two possible sources of scatter to explain the observed distribution: a shape noise term associated with the underlying assumption of circular stable rotation, and an astrophysical signal consistent with a log-normal dispersion around the stellar-to-halo mass relation (SHMR). We have modelled the observed distribution as the combination of these two factors and quantified their most likely values given our data. For the current sample, we measure an effective shape noise of $\sigma_\gamma = 0.024 \pm 0.007$, highlighting the low noise impact of the method and positioning PWL as $\sim 10$ times more precise than conventional weak lensing. We also measure an average dispersion in shears of $\xi_\gamma = 0.53^{+0.26}_{-0.28}$\,dex over the range of $8.5 < \log M_\star < 11$. This measurement is higher than expected, which is suggestive of a relatively high dispersion in halo mass and/or profile.

We present catalogues of stellar masses, star formation rates, and ancillary stellar population parameters for galaxies spanning $02.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., Comment: Accepted for publication in MNRAS. 28 pages, 20 figures

We present the velocity dispersion measurements of four massive $\sim10^{11}M_\odot$ quiescent galaxies at $3.2 < z < 3.7$ based on deep H and K$-$band spectra using the Keck/MOSFIRE near-infrared spectrograph. We find high velocity dispersions of order $\sigma_e\sim250$ km/s based on strong Balmer absorption lines and combine these with size measurements based on HST/WFC3 F160W imaging to infer dynamical masses. The velocity dispersion are broadly consistent with the high stellar masses and small sizes. Together with evidence for quiescent stellar populations, the spectra confirm the existence of a population of massive galaxies that formed rapidly and quenched in the early universe $z>4$. Investigating the evolution at constant velocity dispersion between $z\sim3.5$ and $z\sim2$, we find a large increase in effective radius $0.35\pm0.12$ dex and in dynamical-to-stellar mass ratio $<$log(M$_{\textrm{dyn}}$/M*)$>$ of 0.33$\pm0.08$ dex, with low expected contribution from dark matter. The dynamical masses for our $z\sim3.5$ sample are consistent with the stellar masses for a Chabrier initial mass function (IMF), with the ratio $<$log(M$_{\textrm{dyn}}$/M$^*_{\textrm{Ch}})>$ = -0.13$\pm$0.10 dex suggesting an IMF lighter than Salpeter may be common for massive quiescent galaxies at $z>3$. This is surprising in light of the Salpeter or heavier IMFs found for high velocity dispersion galaxies at $z\sim2$ and cores of present-day ellipticals, which these galaxies are thought to evolve into. Future imaging and spectroscopic observations with resolved kinematics using the upcoming James Webb Space Telescope could rule out potential systematics from rotation, and confirm these results., Comment: 11 pages, 3 figures. Accepted to ApJ Letters

We present an end-to-end methodology to measure the effects of weak lensing on individual galaxy-galaxy systems exploiting their kinematic information. Using this methodology, we have measured a shear signal from the velocity fields of 18 weakly-lensed galaxies. We selected a sample of systems based only on the properties of the sources, requiring them to be bright (apparent $i$-band magnitude $ < 17.4$) and in the nearby Universe ($z < 0.15$). We have observed the velocity fields of the sources with WiFeS, an optical IFU on a 2.3m telescope, and fitted them using a simple circular motion model with an external shear. We have measured an average shear of $\langle \gamma \rangle = 0.020 \pm 0.008$ compared to a predicted $\langle \gamma_{pred} \rangle = 0.005$ obtained using median stellar-to-halo relationships from the literature. While still a statistical approach, our results suggest that this new weak lensing methodology can overcome some of the limitations of traditional stacking-based techniques. We describe in detail all the steps of the methodology and make publicly available all the velocity maps for the weakly-lensed sources used in this study.

We present a detailed study of emission-line systems in the GAMA G23 region, making use of $\textit{WISE}$ photometry that includes carefully measured resolved sources. After applying several cuts to the initial catalogue of $\sim$41,000 galaxies, we extract a sample of 9,809 galaxies. We then compare the spectral diagnostic (BPT) classification of 1154 emission-line galaxies (38$\%$ resolved in W1) to their location in the $\textit{WISE}$ colour-colour diagram, leading to the creation of a new zone for mid-infrared "warm" galaxies located 2$\sigma$ above the star-forming sequence, below the standard $\textit{WISE}$ AGN region. We find that the BPT and $\textit{WISE}$ diagrams agree on the classification for 85$\%$ and 8$\%$ of the galaxies as non-AGN (star forming = SF) and AGN, respectively, and disagree on $\sim$7$\%$ of the entire classified sample. 39$\%$ of the AGN (all types) are broad-line systems for which the [\ion{N}{ii}] and [H$\alpha$] fluxes can barely be disentangled, giving in most cases spurious [\ion{N}{ii}]/[H$\alpha$] flux ratios. However, several optical AGN appear to be completely consistent with SF in $\textit{WISE}$. We argue that these could be low power AGN, or systems whose hosts dominate the IR emission. Alternatively, given the sometimes high [\ion{O}{iii}] luminosity in these galaxies, the emission lines may be generated by shocks coming from super-winds associated with SF rather than the AGN activity. Based on our findings, we have created a new diagnostic: [W1-W2] vs [\ion{N}{ii}]/[H$\alpha$], which has the virtue of separating SF from AGN and high-excitation sources. It classifies 3$\sim$5 times more galaxies than the classic BPT, Comment: 43 pages, 32 figures, 4 tables

We simultaneously present constraints on the stellar-to-halo mass relation for central and satellite galaxies through a weak lensing analysis of spectroscopically classified galaxies. Using overlapping data from the fourth data release of the Kilo-Degree Survey (KiDS), and the Galaxy And Mass Assembly survey (GAMA), we find that satellite galaxies are hosted by halo masses that are $0.53 \pm 0.39$ dex (68\% confidence, $3\sigma$ detection) smaller than those of central galaxies of the same stellar mass (for a stellar mass of $\log(M_{\star}/M_{\odot}) = 10.6$). This is consistent with galaxy formation models, whereby infalling satellite galaxies are preferentially stripped of their dark matter. We find consistent results with similar uncertainties when comparing constraints from a standard azimuthally averaged galaxy-galaxy lensing analysis and a two-dimensional likelihood analysis of the full shear field. As the latter approach is somewhat biased due to the lens incompleteness and as it does not provide any improvement to the precision when applied to actual data, we conclude that stacked tangential shear measurements are best-suited for studies of the galaxy-halo connection., Comment: 12 pages, 10 figures, accepted for publication in A&A

We use KiDS weak lensing data to measure variations in mean halo mass as a function of several key galaxy properties (namely: stellar colour, specific star formation rate, Sersic index, and effective radius) for a volume-limited sample of GAMA galaxies in a narrow stellar mass range ($M_* \sim 2$--$5 \times 10^{10}$ Msol). This mass range is particularly interesting, inasmuch as it is where bimodalities in galaxy properties are most pronounced, and near to the break in both the galaxy stellar mass function and the stellar-to-halo mass relation (SHMR). In this narrow mass range, we find that both size and Sersic index are better predictors of halo mass than either colour or SSFR, with the data showing a slight preference for Sersic index. In other words, we find that mean halo mass is more tightly correlated with galaxy structure than either past star formation history or current star formation rate. Our results lead to an approximate lower bound on the dispersion in halo masses among $\log M_* \approx {10.5}$ galaxies: we find that the dispersion is $\gtrsim 0.3$ dex. This would imply either that offsets from the mean SHMR are closely coupled to size/structure, or that the dispersion in the SHMR is larger than past results have suggested. Our results thus provide new empirical constraints on the relationship between stellar and halo mass assembly at this particularly interesting mass range., Comment: 16 pages, 7 figures; submitted to MNRAS; updated to reflect accepted version

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., Comment: 24 pages, 16 figures, accepted for publication by MNRAS

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., Comment: 23 pages, 18 figures, accepted for publication in MNRAS