Showing 1–50 of 124 results for author: Parkinson, D

STag II: Classification of Serendipitous Supernovae Observed by Galaxy Redshift Surveys

Authors: W. Davison, D. Parkinson, S. BenZvi, A. Palmese, J. Aguilar, S. Ahlen, D. Brooks, T. Claybaugh, A. de la Macorra, Arjun Dey, P. Doel, E. Gaztañaga, S. Gontcho A Gontcho, C. Howlett, S. Juneau, T. Kisner, A. Kremin, A. Lambert, M. Landriau, L. Le Guillou, A. Meisner, R. Miquel, J. Moustakas, A. D. Myers, C. Poppett , et al. (10 additional authors not shown)

Abstract: With the number of supernovae observed expected to drastically increase thanks to large-scale surveys like the Dark Energy Spectroscopic Instrument (DESI), it is necessary that the tools we use to classify these objects keep up with this increase. We previously created Supernova Tagging and Classification (STag) to address this problem by employing machine learning techniques alongside logistic re… ▽ More With the number of supernovae observed expected to drastically increase thanks to large-scale surveys like the Dark Energy Spectroscopic Instrument (DESI), it is necessary that the tools we use to classify these objects keep up with this increase. We previously created Supernova Tagging and Classification (STag) to address this problem by employing machine learning techniques alongside logistic regression in order to assign 'tags' to spectra based on spectral features. STag II is a continuation of this work, which now makes use of model supernova spectra combined with real DESI spectra in order to train STag to better deal with realistic data. We also make use of the rlap score as a trustworthiness cut, making for a more robust and accurate supernova classifier than before. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 27 pages, 7 figures. The STag code is available at https://github.com/wdavison909/STag

The MOST Hosts Survey: spectroscopic observation of the host galaxies of ~40,000 transients using DESI

Authors: Maayane T. Soumagnac, Peter Nugent, Robert A. Knop, Anna Y. Q. Ho, William Hohensee, Autumn Awbrey, Alexis Andersen, Greg Aldering, Matan Ventura, Jessica N. Aguilar, Steven Ahlen, Segev Y. Benzvi, David Brooks, Dillon Brout, Todd Claybaugh, Tamara M. Davis, Kyle Dawson, Axel de la Macorra, Arjun Dey, Biprateep Dey, Peter Doel, Kelly A. Douglass, Jaime E. Forero-Romero, Enrique Gaztanaga, Satya Gontcho A Gontcho , et al. (32 additional authors not shown)

Abstract: We present the MOST Hosts survey (Multi-Object Spectroscopy of Transient Hosts). The survey is planned to run throughout the five years of operation of the Dark Energy Spectroscopic Instrument (DESI) and will generate a spectroscopic catalog of the hosts of most transients observed to date, in particular all the supernovae observed by most public, untargeted, wide-field, optical surveys (PTF/iPTF,… ▽ More We present the MOST Hosts survey (Multi-Object Spectroscopy of Transient Hosts). The survey is planned to run throughout the five years of operation of the Dark Energy Spectroscopic Instrument (DESI) and will generate a spectroscopic catalog of the hosts of most transients observed to date, in particular all the supernovae observed by most public, untargeted, wide-field, optical surveys (PTF/iPTF, SDSS II, ZTF, DECAT, DESIRT). Scientific questions for which the MOST Hosts survey will be useful include Type Ia supernova cosmology, fundamental plane and peculiar velocity measurements, and the understanding of the correlations between transients and their host galaxy properties. Here, we present the first release of the MOST Hosts survey: 21,931 hosts of 20,235 transients. These numbers represent 36% of the final MOST Hosts sample, consisting of 60,212 potential host galaxies of 38,603 transients (a transient can be assigned multiple potential hosts). Of these galaxies, 40% do not appear in the DESI primary target list and therefore require a specific program like MOST Hosts. Of all the transients in the MOST Hosts list, only 26.7% have existing classifications, and so the survey will provide redshifts (and luminosities) for nearly 30,000 transients. A preliminary Hubble diagram and a transient luminosity-duration diagram are shown as examples of future potential uses of the MOST Hosts survey. The survey will also provide a training sample of spectroscopically observed transients for photometry-only classifiers, as we enter an era when most newly observed transients will lack spectroscopic classification. The MOST Hosts DESI survey data will be released through the Wiserep platform on a rolling cadence and updated to match the DESI releases. Dates of future releases and updates are available through the https://mosthosts.desi.lbl.gov website. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to ApJS

Planetary Complexity Revealed by the Joint Differential Entropy of Eigencolours

Authors: Gary Segal, David Parkinson, Stuart Bartlett

Abstract: We propose a measure, the joint differential entropy of eigencolours, for determining the spatial complexity of exoplanets using only spatially unresolved light curve data. The measure can be used to search for habitable planets, based on the premise of a potential association between life and exoplanet complexity. We present an analysis using disk-integrated light curves from Earth, developed in… ▽ More We propose a measure, the joint differential entropy of eigencolours, for determining the spatial complexity of exoplanets using only spatially unresolved light curve data. The measure can be used to search for habitable planets, based on the premise of a potential association between life and exoplanet complexity. We present an analysis using disk-integrated light curves from Earth, developed in previous studies, as a proxy for exoplanet data. We show that this quantity is distinct from previous measures of exoplanet complexity due to its sensitivity to spatial information that is masked by features with large mutual information between wavelengths, such as cloud cover. The measure has a natural upper limit and appears to avoid a strong bias toward specific planetary features. This makes it a candidate for being used as a generalisable measure of exoplanet habitability, since it is agnostic regarding the form that life could take. △ Less

Submitted 15 November, 2023; v1 submitted 12 November, 2023; originally announced November 2023.

Comments: 10 pages, 3 figures, 2 tables

Journal ref: AJ 167 114 (2024)

The Galaxy Number Density Profile of Haloes

Authors: Fei Qin, David Parkinson, Adam R. H. Stevens, Cullan Howlett

Abstract: More precise measurements of galaxy clustering will be provided by the next generation of galaxy surveys such as DESI, WALLABY and SKA. To utilize this information to improve our understanding of the Universe, we need to accurately model the distribution of galaxies in their host dark matter halos. In this work we present a new galaxy number density profile of haloes, which makes predictions for t… ▽ More More precise measurements of galaxy clustering will be provided by the next generation of galaxy surveys such as DESI, WALLABY and SKA. To utilize this information to improve our understanding of the Universe, we need to accurately model the distribution of galaxies in their host dark matter halos. In this work we present a new galaxy number density profile of haloes, which makes predictions for the positions of galaxies in the host halo, different to the widely adopted Navarro-Frenk-White (NFW) profile, since galaxies tend to be found more in the outskirts of halos (nearer the virial radius) than an NFW profile. The parameterised galaxy number density profile model of haloes is fit and tested using the DARKSAGE semi-analytic model of galaxy formation. We find that our galaxy number density profile model of haloes can accurately reproduce the halo occupation distribution and galaxy two-point correlation function of the DARKSAGE simulation. We also derive the analytic expressions for the circular velocity and gravitational potential energy for this profile model. We use the SDSS DR10 galaxy group catalogue to validate this galaxy number density profile model of haloes. Compared to the NFW profile, we find that our model more accurately predicts the positions of galaxies in their host halo and the galaxy two-point correlation function. △ Less

Submitted 25 September, 2023; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 13 pages. 10 figures. Appear on ApJ

NANCY: Next-generation All-sky Near-infrared Community surveY

Authors: Jiwon Jesse Han, Arjun Dey, Adrian M. Price-Whelan, Joan Najita, Edward F. Schlafly, Andrew Saydjari, Risa H. Wechsler, Ana Bonaca, David J Schlegel, Charlie Conroy, Anand Raichoor, Alex Drlica-Wagner, Juna A. Kollmeier, Sergey E. Koposov, Gurtina Besla, Hans-Walter Rix, Alyssa Goodman, Douglas Finkbeiner, Abhijeet Anand, Matthew Ashby, Benedict Bahr-Kalus, Rachel Beaton, Jayashree Behera, Eric F. Bell, Eric C Bellm , et al. (184 additional authors not shown)

Abstract: The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GAL… ▽ More The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

The Early Data Release of the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (240 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes… ▽ More The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes good-quality spectral information from 466,447 objects targeted as part of the Milky Way Survey, 428,758 as part of the Bright Galaxy Survey, 227,318 as part of the Luminous Red Galaxy sample, 437,664 as part of the Emission Line Galaxy sample, and 76,079 as part of the Quasar sample. In addition, the release includes spectral information from 137,148 objects that expand the scope beyond the primary samples as part of a series of secondary programs. Here, we describe the spectral data, data quality, data products, Large-Scale Structure science catalogs, access to the data, and references that provide relevant background to using these spectra. △ Less

Submitted 15 June, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 43 pages, 7 figures, 17 tables, submitted to AJ, DESI EDR references added

Validation of the Scientific Program for the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (239 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of… ▽ More The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar (MWS), bright galaxy (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the five-year program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a `One-Percent survey' conducted at the conclusion of Survey Validation covering 140 deg$^2$ using the final target selection algorithms with exposures of a depth typical of the main survey. The Survey Validation indicates that DESI will be able to complete the full 14,000 deg$^2$ program with spectroscopically-confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval $z<1.1$, 0.39% over the redshift interval $1.1<z<1.9$, and 0.46% over the redshift interval $1.9<z<3.5$. △ Less

Submitted 12 January, 2024; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 42 pages, 18 figures, accepted by AJ

Maximum Likelihood based Phase-Retrieval using Fresnel Propagation Forward Models with Optional Constraints

Authors: K. Aditya Mohan, Jean-Baptiste Forien, Venkatesh Sridhar, Jefferson A. Cuadra, Dilworth Parkinson

Abstract: X-ray phase-contrast tomography (XPCT) is widely used for high contrast 3D imaging using either synchrotron or laboratory microfocus X-ray sources. XPCT enables an order of magnitude improvement in image contrast of the reconstructed material interfaces with low X-ray absorption contrast. The dominant approaches to 3D reconstruction using XPCT relies on the use of phase-retrieval algorithms that m… ▽ More X-ray phase-contrast tomography (XPCT) is widely used for high contrast 3D imaging using either synchrotron or laboratory microfocus X-ray sources. XPCT enables an order of magnitude improvement in image contrast of the reconstructed material interfaces with low X-ray absorption contrast. The dominant approaches to 3D reconstruction using XPCT relies on the use of phase-retrieval algorithms that make one or more limiting approximations for the experimental configuration and material properties. Since many experimental scenarios violate such approximations, the resulting reconstructions contain blur, artifacts, or other quantitative inaccuracies. Our solution to this problem is to formulate new iterative non-linear phase-retrieval (NLPR) algorithms that avoid such limiting approximations. Compared to the widely used state-of-the-art approaches, we show that our proposed algorithms result in sharp and quantitatively accurate reconstruction with reduced artifacts. Unlike existing NLPR algorithms, our approaches avoid the laborious manual tuning of regularization hyper-parameters while still achieving the stated goals. As an alternative to regularization, we propose explicit constraints on the material properties to constrain the solution space and solve the phase-retrieval problem. These constraints are easily user-configurable since they follow directly from the imaged object's dimensions and material properties. △ Less

Submitted 2 October, 2023; v1 submitted 29 April, 2023; originally announced May 2023.

Comments: This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

tomoCAM: Fast Model-based Iterative Reconstruction via GPU Acceleration and Non-Uniform Fast Fourier Transforms

Authors: Dinesh Kumar, Dilworth Y. Parkinson, Jeffrey J. Donatelli

Abstract: X-Ray based computed tomography (CT) is a well-established technique for determining the three-dimensional structure of an object from its two-dimensional projections. In the past few decades, there have been significant advancements in the brightness and detector technology of tomography instruments at synchrotron sources. These advancements have led to the emergence of new observations and disco… ▽ More X-Ray based computed tomography (CT) is a well-established technique for determining the three-dimensional structure of an object from its two-dimensional projections. In the past few decades, there have been significant advancements in the brightness and detector technology of tomography instruments at synchrotron sources. These advancements have led to the emergence of new observations and discoveries, with improved capabilities such as faster frame rates, larger fields of view, higher resolution, and higher dimensionality. These advancements have enabled the material science community to expand the scope of tomographic measurements towards increasingly in-situ and in-operando measurements. In these new experiments, samples can be rapidly evolving, have complex geometries, and restrictions on the field of view, limiting the number of projections that can be collected. In such cases, standard filtered back-projections (FBP) for the reconstructions often result in poor-quality reconstructions. Iterative reconstruction algorithms, such as model-based iterative reconstructions (MBIR), have demonstrated considerable success in producing high-quality reconstructions under such restrictions, but typically require high-performance computing resources with hundreds of compute nodes to solve the problem in a reasonable time. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Measurement of the matter-radiation equality scale using the extended Baryon Oscillation Spectroscopic Survey Quasar Sample

Authors: Benedict Bahr-Kalus, David Parkinson, Eva-Maria Mueller

Abstract: The position of the peak of the matter power spectrum, the so-called turnover scale, is set by the horizon size at the epoch of matter-radiation equality. It can easily be predicted in terms of the physics of the Universe in the relativistic era, and so can be used as a standard ruler, independent of other features present in the matter power spectrum, such as baryon acoustic oscillations (BAO). W… ▽ More The position of the peak of the matter power spectrum, the so-called turnover scale, is set by the horizon size at the epoch of matter-radiation equality. It can easily be predicted in terms of the physics of the Universe in the relativistic era, and so can be used as a standard ruler, independent of other features present in the matter power spectrum, such as baryon acoustic oscillations (BAO). We use the distribution of quasars measured by the extended Baryon Oscillation Spectroscopic Survey (eBOSS) to determine the turnover scale in a model-independent fashion statistically. We avoid modelling the BAO by down-weighting affected scales in the covariance matrix using the mode deprojection technique. We measure the wavenumber of the peak to be $k_\mathrm{TO} = \left( 17.6^{+1.9}_{-1.8} \right) \times 10^{-3}h/\mathrm{Mpc}$, corresponding to a dilation scale of $ D_\mathrm{V}(z_\mathrm{eff} = 1.48) = \left({36.2^{+4.1}_{-4.4}}\right)r_\mathrm{H}$. This is not competitive with current BAO distance measures in terms of determining the expansion history but does provide a useful cross-check. We combine this measurement with low-redshift distance measurements from type-Ia supernova data from Pantheon and BAO data from eBOSS to make a sound-horizon free estimate of the Hubble-Lemaître parameter and find it to be $H_0=\left({74.7\pm 9.6}\right) \ \mathrm{km/s/Mpc}$ with Pantheon, and $H_0=\left({72.9^{+10.0}_{-8.6}}\right) \ \mathrm{km/s/Mpc}$ with eBOSS BAO. We make predictions for the measurement of the turnover scale by the Dark Energy Spectroscopic Instrument (DESI) survey, the Maunakea Spectroscopic Explorer (MSE) and MegaMapper, which will make more precise and accurate distance determinations. △ Less

Submitted 5 October, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: 14 pages, 13 figures, includes erratum past equation (24) published by MNRAS (DOI: 10.1093/mnras/stad3000)

Reconstructing the cosmological density and velocity fields from redshifted galaxy distributions using V-net

Authors: Fei Qin, David Parkinson, Sungwook E. Hong, Cristiano G. Sabiu

Abstract: The distribution of matter that is measured through galaxy redshift and peculiar velocity surveys can be harnessed to learn about the physics of dark matter, dark energy, and the nature of gravity. To improve our understanding of the matter of the Universe, we can reconstruct the full density and velocity fields from the galaxies that act as tracer particles. In this paper, we use the simulated ha… ▽ More The distribution of matter that is measured through galaxy redshift and peculiar velocity surveys can be harnessed to learn about the physics of dark matter, dark energy, and the nature of gravity. To improve our understanding of the matter of the Universe, we can reconstruct the full density and velocity fields from the galaxies that act as tracer particles. In this paper, we use the simulated halos as proxies for the galaxies. We use a convolutional neural network, a V-net, trained on numerical simulations of structure formation to reconstruct the density and velocity fields. We find that, with detailed tuning of the loss function, the V-net could produce better fits to the density field in the high-density and low-density regions, and improved predictions for the probability distribution of the amplitudes of the velocities. However, the weights will reduce the precision of the estimated $β$ parameter. We also find that the redshift-space distortions of the halo catalogue do not significantly contaminate the reconstructed real-space density and velocity field. We estimate the velocity field $β$ parameter by comparing the peculiar velocities of halo catalogues to the reconstructed velocity fields, and find the estimated $β$ values agree with the fiducial value at the 68\% confidence level. △ Less

Submitted 19 June, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

Comments: 25 pages, 18 figures, 2 tables. Published JCAP

Journal ref: JCAP06(2023)062

Constraining Cosmic Inflation with Observations: Prospects for 2030

Authors: Benedict Bahr-Kalus, David Parkinson, Richard Easther

Abstract: The ability to test and constrain theories of cosmic inflation will advance substantially over the next decade. Key data sources include cosmic microwave background (CMB) measurements and observations of the distribution of matter at low-redshift from optical, near-infrared, and 21cm intensity surveys. A positive detection of a CMB B-mode consistent with a primordial stochastic gravitational wave… ▽ More The ability to test and constrain theories of cosmic inflation will advance substantially over the next decade. Key data sources include cosmic microwave background (CMB) measurements and observations of the distribution of matter at low-redshift from optical, near-infrared, and 21cm intensity surveys. A positive detection of a CMB B-mode consistent with a primordial stochastic gravitational wave background (SGWB) is widely viewed as a smoking gun for an inflationary phase. Still, a null result does not exclude inflation. However, in a significant class of inflationary scenarios, a low SGWB amplitude is correlated with a more significant running, $α_s$, in the primordial density perturbations than is seen with the simplest inflationary potentials. With this motivation, we forecast the precision with which the spectral index $n_{\rm{s}}$ and $α_{\rm{s}}$ can be constrained by currently envisaged observations, including CMB (Simons Observatory, CMB-S4 and LiteBIRD), optical/near infra-red (DESI and SPHEREx), and 21cm intensity mapping (Tianlai and CHIME) surveys. We identify optimal combinations of datasets for constraining the running and show that they may yield additional and informative constraints on the overall inflationary parameter space if the SGWB remains undetected. △ Less

Submitted 10 January, 2023; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: 12 pages, 8 figures. Version matches the one accepted by MNRAS

Detecting complex sources in large surveys using an apparent complexity measure

Authors: David Parkinson, Gary Segal

Abstract: Large area astronomical surveys will almost certainly contain new objects of a type that have never been seen before. The detection of 'unknown unknowns' by an algorithm is a difficult problem to solve, as unusual things are often easier for a human to spot than a machine. We use the concept of apparent complexity, previously applied to detect multi-component radio sources, to scan the radio conti… ▽ More Large area astronomical surveys will almost certainly contain new objects of a type that have never been seen before. The detection of 'unknown unknowns' by an algorithm is a difficult problem to solve, as unusual things are often easier for a human to spot than a machine. We use the concept of apparent complexity, previously applied to detect multi-component radio sources, to scan the radio continuum Evolutionary Map of the Universe (EMU) Pilot Survey data for complex and interesting objects in a fully automated and blind manner. Here we describe how the complexity is defined and measured, how we applied it to the Pilot Survey data, and how we calibrated the completeness and purity of these interesting objects using a crowd-sourced 'zoo'. The results are also compared to unexpected and unusual sources already detected in the EMU Pilot Survey, including Odd Radio Circles, that were found by human inspection. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 6 pages, 4 figures. Prepared for the proceedings of the International Astronomical Union Symposium 368 "Machine Learning in Astronomy: Possibilities and Pitfalls". arXiv admin note: text overlap with arXiv:2206.14677

Panel Discussion: Practical Problem Solving for Machine Learning

Authors: Guillermo Cabrera, Sungwook E. Hong, Lilianne Nakazono, David Parkinson, Yuan-Sen Ting

Abstract: Machine Learning is a powerful tool for astrophysicists, which has already had significant uptake in the community. But there remain some barriers to entry, relating to proper understanding, the difficulty of interpretability, and the lack of cohesive training. In this discussion session we addressed some of these questions, and suggest how the field may move forward. Machine Learning is a powerful tool for astrophysicists, which has already had significant uptake in the community. But there remain some barriers to entry, relating to proper understanding, the difficulty of interpretability, and the lack of cohesive training. In this discussion session we addressed some of these questions, and suggest how the field may move forward. △ Less

Submitted 30 November, 2022; originally announced November 2022.

Comments: 6 pages. Prepared for the proceedings of the International Astronomical Union Symposium 368 "Machine Learning in Astronomy: Possibilities and Pitfalls"

HI HOD -- I. The Halo Occupation Distribution of HI Galaxies

Authors: Fei Qin, Cullan Howlett, Adam R. H. Stevens, David Parkinson

Abstract: The next generation of galaxy surveys will provide more precise measurements of galaxy clustering than have previously been possible. The 21-cm radio signals that are emitted from neutral atomic hydrogen (HI) gas will be detected by large-area radio surveys such as WALLABY and the SKA, and deliver galaxy positions and velocities that can be used to measure galaxy clustering statistics. But, to har… ▽ More The next generation of galaxy surveys will provide more precise measurements of galaxy clustering than have previously been possible. The 21-cm radio signals that are emitted from neutral atomic hydrogen (HI) gas will be detected by large-area radio surveys such as WALLABY and the SKA, and deliver galaxy positions and velocities that can be used to measure galaxy clustering statistics. But, to harness this information to improve our cosmological understanding, and learn about the physics of dark matter and dark energy, we need to accurately model the manner in which galaxies detected in HI trace the underlying matter distribution of the Universe. For this purpose, we develop a new HI-based Halo Occupation Distribution (HOD) model, which makes predictions for the number of galaxies present in dark matter halos conditional on their HI mass. The parameterised HOD model is fit and validated using the Dark Sage semi-analytic model, where we show that the HOD parameters can be modelled by simple linear and quadratic functions of HI mass. However, we also find that the clustering predicted by the HOD depends sensitively on the radial distributions of the HI galaxies within their host dark matter halos, which does not follow the NFW profile in the Dark Sage simulation. As such, this work enables -- for the first time -- a simple prescription for placing galaxies of different HI mass within dark matter halos in a way that is able to reproduce the HI mass-dependent galaxy clustering and HI mass function simultaneously and without requiring knowledge of the optical properties of the galaxies. Further efforts are required to demonstrate that this model can be used to produce large ensembles of mock galaxy catalogues for upcoming surveys. △ Less

Submitted 1 September, 2022; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: 18 pages, 13 figures, 1 table. Published in ApJ

The DESI Survey Validation: Results from Visual Inspection of Bright Galaxies, Luminous Red Galaxies, and Emission Line Galaxies

Authors: Ting-Wen Lan, R. Tojeiro, E. Armengaud, J. Xavier Prochaska, T. M. Davis, David M. Alexander, A. Raichoor, Rongpu Zhou, Christophe Yeche, C. Balland, S. BenZvi, A. Berti, R. Canning, A. Carr, H. Chittenden, S. Cole, M. -C. Cousinou, K. Dawson, Biprateep Dey, K. Douglass, A. Edge, S. Escoffier, A. Glanville, S. Gontcho A Gontcho, J. Guy , et al. (57 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) Survey has obtained a set of spectroscopic measurements of galaxies to validate the final survey design and target selections. To assist in these tasks, we visually inspect (VI) DESI spectra of approximately 2,500 bright galaxies, 3,500 luminous red galaxies (LRGs), and 10,000 emission line galaxies (ELGs), to obtain robust redshift identifications.… ▽ More The Dark Energy Spectroscopic Instrument (DESI) Survey has obtained a set of spectroscopic measurements of galaxies to validate the final survey design and target selections. To assist in these tasks, we visually inspect (VI) DESI spectra of approximately 2,500 bright galaxies, 3,500 luminous red galaxies (LRGs), and 10,000 emission line galaxies (ELGs), to obtain robust redshift identifications. We then utilize the VI redshift information to characterize the performance of the DESI operation. Based on the VI catalogs, our results show that the final survey design yields samples of bright galaxies, LRGs, and ELGs with purity greater than $99\%$. Moreover, we demonstrate that the precision of the redshift measurements is approximately 10 km/s for bright galaxies and ELGs and approximately 40 km/s for LRGs. The average redshift accuracy is within 10 km/s for the three types of galaxies. The VI process also helps improve the quality of the DESI data by identifying spurious spectral features introduced by the pipeline. Finally, we show examples of unexpected real astronomical objects, such as Ly$α$ emitters and strong lensing candidates, identified by VI. These results demonstrate the importance and utility of visually inspecting data from incoming and upcoming surveys, especially during their early operation phases. △ Less

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

Comments: 22 pages, 14 figures, 3 tables, one of a suite of 8 papers detailing targeting for DESI. ApJ accepted version with minor textual updates

Identifying anomalous radio sources in the EMU Pilot Survey using a complexity-based approach

Authors: Gary Segal, David Parkinson, Ray Norris, Andrew M. Hopkins, Heinz Andernach, Emma L. Alexander, Ettore Carretti, Bärbel S. Koribalski, Letjatji S. Legodi, Sarah Leslie, Yan Luo, Jonathon C. S. Pierce, Hongming Tang, Eleni Vardoulaki, Tessa Vernstrom

Abstract: The Evolutionary Map of the Universe (EMU) large-area radio continuum survey will detect tens of millions of radio galaxies, giving an opportunity for the detection of previously unknown classes of objects. To maximise the scientific value and make new discoveries, the analysis of this data will need to go beyond simple visual inspection. We propose the coarse-grained complexity, a simple scalar q… ▽ More The Evolutionary Map of the Universe (EMU) large-area radio continuum survey will detect tens of millions of radio galaxies, giving an opportunity for the detection of previously unknown classes of objects. To maximise the scientific value and make new discoveries, the analysis of this data will need to go beyond simple visual inspection. We propose the coarse-grained complexity, a simple scalar quantity relating to the minimum description length of an image, that can be used to identify unusual structures. The complexity can be computed without reference to the broader sample or existing catalogue data, making the computation efficient on new surveys at very large scales (such as the full EMU survey). We apply our coarse-grained complexity measure to data from the EMU Pilot Survey to detect and confirm anomalous objects in this data set and produce an anomaly catalogue. Rather than work with existing catalogue data using a specific source detection algorithm, we perform a blind scan of the area, computing the complexity using a sliding square aperture. The effectiveness of the complexity measure for identifying anomalous objects is evaluated using crowd-sourced labels generated via the Zooniverse.org platform. We find that the complexity scan identifies unusual sources, such as odd radio circles, by partitioning on complexity. We achieve partitions where 5\% of the data is estimated to be 86\% complete, and 0.5\% is estimated to be 94\% pure, with respect to anomalies and use this to produce an anomaly catalogue. △ Less

Submitted 3 March, 2023; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: 19 pages, 14 figures. The EMU Pilot Survey anomaly catalogues are available on direct request from the first author

Journal ref: MNRAS 521, 1429-1447 (2023)

Nested sampling for physical scientists

Authors: Greg Ashton, Noam Bernstein, Johannes Buchner, Xi Chen, Gábor Csányi, Andrew Fowlie, Farhan Feroz, Matthew Griffiths, Will Handley, Michael Habeck, Edward Higson, Michael Hobson, Anthony Lasenby, David Parkinson, Livia B. Pártay, Matthew Pitkin, Doris Schneider, Joshua S. Speagle, Leah South, John Veitch, Philipp Wacker, David J. Wales, David Yallup

Abstract: We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the applic… ▽ More We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: 20 pages + supplementary information, 5 figures. preprint version; published version at https://www.nature.com/articles/s43586-022-00121-x

Journal ref: Nature Reviews Methods Primers volume 2, Article number: 39 (2022)

Venus as an Exoplanet: I. An Initial Exploration of the 3-D Energy Balance for a CO$_{2}$ Exoplanetary Atmosphere Around an M-Dwarf Star

Authors: Christopher D. Parkinson, Stephen W. Bougher, Franklin P. Mills, Amanda Brecht, James Li, Yuk L. Yung, Renyu Hu, Guillaume Gronoff

Abstract: The chemical evolution of an exoplanetary Venus-like atmosphere is dependent upon the ultraviolet to near ultraviolet (FUV-NUV) radiation ratio from the parent star, the balance between CO$_{2}$ photolysis and recombination via reactions that depend on the water abundance, and various catalytic chemical cycles. In this study, we use a three-dimensional (3-D) model to simulate conditions for a Venu… ▽ More The chemical evolution of an exoplanetary Venus-like atmosphere is dependent upon the ultraviolet to near ultraviolet (FUV-NUV) radiation ratio from the parent star, the balance between CO$_{2}$ photolysis and recombination via reactions that depend on the water abundance, and various catalytic chemical cycles. In this study, we use a three-dimensional (3-D) model to simulate conditions for a Venus-like exoplanet orbiting the M-dwarf type star GJ 436 by varying the star/planet distance and considering the resultant effects on heating/cooling and dynamics. The simulation includes the middle and upper atmosphere (<40 mbar). Overall, these model comparisons reveal that the impact of extreme ultraviolet to ultraviolet (EUV-UV) heating on the energy balance shows both radiative and dynamical processes are responsible for driving significant variations in zonal winds and global temperature profiles at < 10$^{-5}$ mbar. More specifically, CO$_{2}$ 15-$μ$m cooling balances EUV/UV and Near InfraRed (NIR) heating at altitudes below 10$^{-7}$ mbar pressure with a strong maximum balance for pressures at ~10$^{-5}$ mbar, thus explaining the invariance of the temperature distribution at altitudes below 10$^{-5}$mbar pressure for all cases. Our model comparisons also show that moderate changes in NIR heating result in relatively small changes in neutral temperature in the upper atmosphere, and virtually no change in the middle atmosphere. However, with larger changes in the NIR heating profile, much greater changes in neutral temperature occur in the entire upper and middle atmosphere studied. △ Less

Submitted 30 May, 2024; v1 submitted 22 May, 2022; originally announced May 2022.

Comments: 44 pages, 11 Figures, 3 tables

A measurement of the Integrated Sachs-Wolfe Effect with the Rapid ASKAP Continuum Survey

Authors: Benedict Bahr-Kalus, David Parkinson, Jacobo Asorey, Stefano Camera, Catherine Hale, Fei Qin

Abstract: The evolution of the gravitational potentials on large scales due to the accelerated expansion of the Universe is an important and independent probe of dark energy, known as the integrated Sachs-Wolfe (ISW) effect. We measure this ISW effect through cross-correlating the cosmic microwave background maps from the Planck satellite with a radio continuum galaxy distribution map from the recent Rapid… ▽ More The evolution of the gravitational potentials on large scales due to the accelerated expansion of the Universe is an important and independent probe of dark energy, known as the integrated Sachs-Wolfe (ISW) effect. We measure this ISW effect through cross-correlating the cosmic microwave background maps from the Planck satellite with a radio continuum galaxy distribution map from the recent Rapid ASKAP Continuum Survey (RACS). We detect a positive cross-correlation at $\sim 2.8\,σ$ relative to the null hypothesis of no correlation. We parameterise the strength of the ISW effect through an amplitude parameter and find the constraints to be $A_{\mathrm{ISW}} = 0.94^{+0.42}_{-0.41}$, which is consistent with the prediction of an accelerating universe within the current concordance cosmological model, $Λ$CDM. The credible interval on this parameter is independent of the different bias models and redshift distributions that were considered when marginalising over the nuisance parameters. We also detect a power excess in the galaxy auto-correlation angular power spectrum on large scales ($\ell \leq 40$), and investigate possible systematic causes. △ Less

Submitted 28 July, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Comments: 18 pages, 18 figures. Angular power spectra and covariance matrices can be found at https://github.com/racs-cosmology/isw

A workflow for segmenting soil and plant X-ray CT images with deep learning in Googles Colaboratory

Authors: Devin A. Rippner, Pranav Raja, J. Mason Earles, Alexander Buchko, Mina Momayyezi, Fiona Duong, Dilworth Parkinson, Elizabeth Forrestel, Ken Shackel, Jeffrey Neyhart, Andrew J. McElrone

Abstract: X-ray micro-computed tomography (X-ray microCT) has enabled the characterization of the properties and processes that take place in plants and soils at the micron scale. Despite the widespread use of this advanced technique, major limitations in both hardware and software limit the speed and accuracy of image processing and data analysis. Recent advances in machine learning, specifically the appli… ▽ More X-ray micro-computed tomography (X-ray microCT) has enabled the characterization of the properties and processes that take place in plants and soils at the micron scale. Despite the widespread use of this advanced technique, major limitations in both hardware and software limit the speed and accuracy of image processing and data analysis. Recent advances in machine learning, specifically the application of convolutional neural networks to image analysis, have enabled rapid and accurate segmentation of image data. Yet, challenges remain in applying convolutional neural networks to the analysis of environmentally and agriculturally relevant images. Specifically, there is a disconnect between the computer scientists and engineers, who build these AI/ML tools, and the potential end users in agricultural research, who may be unsure of how to apply these tools in their work. Additionally, the computing resources required for training and applying deep learning models are unique, more common to computer gaming systems or graphics design work, than to traditional computational systems. To navigate these challenges, we developed a modular workflow for applying convolutional neural networks to X-ray microCT images, using low-cost resources in Googles Colaboratory web application. Here we present the results of the workflow, illustrating how parameters can be optimized to achieve best results using example scans from walnut leaves, almond flower buds, and a soil aggregate. We expect that this framework will accelerate the adoption and use of emerging deep learning techniques within the plant and soil sciences. △ Less

Submitted 21 July, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 58 pages, 9 figures, 2 Tables

Running primordial perturbations: Inflationary Dynamics and Observational Constraints

Authors: Richard Easther, Benedict Bahr-Kalus, David Parkinson

Abstract: Inflationary cosmology proposes that the early Universe undergoes accelerated expansion, driven, in simple scenarios, by a single scalar field, or inflaton. The form of the inflaton potential determines the initial spectra of density perturbations and gravitational waves. We show that constraints on the duration of inflation together with the BICEP3/Keck bounds on the gravitational wave background… ▽ More Inflationary cosmology proposes that the early Universe undergoes accelerated expansion, driven, in simple scenarios, by a single scalar field, or inflaton. The form of the inflaton potential determines the initial spectra of density perturbations and gravitational waves. We show that constraints on the duration of inflation together with the BICEP3/Keck bounds on the gravitational wave background imply that higher derivatives of the potential are nontrivial with a confidence of 99%. Such terms contribute to the scale-dependence, or running, of the density perturbation spectrum. We clarify the ``universality classes'' of inflation in this limit showing that a very small gravitational wave background can be correlated with a larger running. If pending experiments do not observe a gravitational wave background the running will be at the threshold of detectability if inflation is well-described at third-order in the slow roll expansion. △ Less

Submitted 15 September, 2022; v1 submitted 20 December, 2021; originally announced December 2021.

Comments: 5 pages; 5 figures; as published in PRD -- change of title, minor clarifications

Journal ref: Physical Review D 106 (6), L061301 (2022)

Cosmic Flow Measurement and Mock Sampling Algorithm of Cosmicflows-4 Tully-Fisher Catalogue

Authors: Fei Qin, David Parkinson, Cullan Howlett, Khaled Said

Abstract: Measurements of cosmic flows enable us to test whether cosmological models can accurately describe the evolution of the density field in the nearby Universe. In this paper, we measure the low-order kinematic moments of the cosmic flow field, namely bulk flow and shear moments, using the Cosmicflows-4 Tully-Fisher catalogue (CF4TF). To make accurate cosmological inferences with the CF4TF sample, it… ▽ More Measurements of cosmic flows enable us to test whether cosmological models can accurately describe the evolution of the density field in the nearby Universe. In this paper, we measure the low-order kinematic moments of the cosmic flow field, namely bulk flow and shear moments, using the Cosmicflows-4 Tully-Fisher catalogue (CF4TF). To make accurate cosmological inferences with the CF4TF sample, it is important to make realistic mock catalogues. We present the mock sampling algorithm of CF4TF. These mock can accurately realize the survey geometry and luminosity selection function, enabling researchers to explore how these systematics affect the measurements. These mocks can also be further used to estimate the covariance matrix and errors of power spectrum and two-point correlation function in future work. In this paper, we use the mocks to test the cosmic flow estimator and find that the measurements are unbiased. The measured bulk flow in the local Universe is 376 $\pm$ 23 (error) $\pm$ 183 (cosmic variance) km s$^{-1}$ at depth $d_{\text{MLE}}=35$ Mpc $h^{-1}$, to the Galactic direction of $(l,b)=(298\pm 3^{\circ}, -6\pm 3^{\circ})$. Both the measured bulk and shear moments are consistent with the concordance $Λ$ Cold Dark Matter cosmological model predictions. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 15 pages, 13 figures, 2 tables. Appear in ApJ

STag: Supernova Tagging and Classification

Authors: William Davison, David Parkinson, Brad E. Tucker

Abstract: Supernovae classes have been defined phenomenologically, based on spectral features and time series data, since the specific details of the physics of the different explosions remain unrevealed. However, the number of these classes is increasing as objects with new features are observed, and the next generation of large-surveys will only bring more variety to our attention. We apply the machine le… ▽ More Supernovae classes have been defined phenomenologically, based on spectral features and time series data, since the specific details of the physics of the different explosions remain unrevealed. However, the number of these classes is increasing as objects with new features are observed, and the next generation of large-surveys will only bring more variety to our attention. We apply the machine learning technique of multi-label classification to the spectra of supernovae. By measuring the probabilities of specific features or `tags' in the supernova spectra, we can compress the information from a specific object down to that suitable for a human or database scan, without the need to directly assign to a reductive `class'. We use logistic regression to assign tag probabilities, and then a feed-forward neural network to filter the objects into the standard set of classes, based solely on the tag probabilities. We present STag, a software package that can compute these tag probabilities and make spectral classifications. △ Less

Submitted 12 November, 2021; v1 submitted 23 August, 2021; originally announced August 2021.

Comments: 20 pages, 9 figures. Pages 13-20 are long tables. The code can be found at https://github.com/wdavison909/STag

Journal ref: The Astrophysical Journal (2022) 925, 186

The Evolutionary Map of the Universe Pilot Survey

Authors: Ray P. Norris, Joshua Marvil, J. D. Collier, Anna D. Kapinska, Andrew N. O'Brien, L. Rudnick, Heinz Andernach, Jacobo Asorey, Michael J. I. Brown, Marcus Bruggen, Evan Crawford, Jayanne English, Syed Faisal ur Rahman, Miroslav D. Filipovic, Yjan Gordon, Gulay Gurkan, Catherine Hale, Andrew M. Hopkins, Minh T. Huynh, Kim HyeongHan, M. James Jee, Baerbel S. Koribalski, Emil Lenc, Kieran Luken, David Parkinson , et al. (23 additional authors not shown)

Abstract: We present the data and initial results from the first Pilot Survey of the Evolutionary Map of the Universe (EMU), observed at 944 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The survey covers 270 \sqdeg of an area covered by the Dark Energy Survey, reaching a depth of 25--30 \ujybm\ rms at a spatial resolution of $\sim$ 11--18 arcsec, resulting in a catalogue of… ▽ More We present the data and initial results from the first Pilot Survey of the Evolutionary Map of the Universe (EMU), observed at 944 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The survey covers 270 \sqdeg of an area covered by the Dark Energy Survey, reaching a depth of 25--30 \ujybm\ rms at a spatial resolution of $\sim$ 11--18 arcsec, resulting in a catalogue of $\sim$ 220,000 sources, of which $\sim$ 180,000 are single-component sources. Here we present the catalogue of single-component sources, together with (where available) optical and infrared cross-identifications, classifications, and redshifts. This survey explores a new region of parameter space compared to previous surveys. Specifically, the EMU Pilot Survey has a high density of sources, and also a high sensitivity to low surface-brightness emission. These properties result in the detection of types of sources that were rarely seen in or absent from previous surveys. We present some of these new results here. △ Less

Submitted 1 August, 2021; originally announced August 2021.

Comments: Accepted by PASA

Future Radio Continuum Cosmology Clustering Surveys

Authors: Jacobo Asorey, David Parkinson

Abstract: The use of continuum emission radio galaxies as cosmological tracers of the large-scale structure will soon move into a new phase. Upcoming surveys from the Australian Square Kilometre Array Pathfinder (ASKAP), MeerKAT, and the Square Kilometre Array project (SKA) will survey the entire available sky down to an ~100uJy flux limit, increasing the number of detected extra-galactic radio sources by s… ▽ More The use of continuum emission radio galaxies as cosmological tracers of the large-scale structure will soon move into a new phase. Upcoming surveys from the Australian Square Kilometre Array Pathfinder (ASKAP), MeerKAT, and the Square Kilometre Array project (SKA) will survey the entire available sky down to an ~100uJy flux limit, increasing the number of detected extra-galactic radio sources by several orders of magnitude. External data and machine learning algorithms will also enable some low resolution radial selection (photometric redshift binning) of the sample, increasing the cosmological utility of the sample observed. In this paper, we discuss the flux limit required to detect enough galaxies to decrease the shot noise term in the error to be 10% of the total. We show how future surveys of this type will be limited by available technology. The confusion generated by the intrinsic sizes of galaxies may have the consequence that surveys of this type eventually reach a hard flux limit of ~100nJy, as is predicted by the current modelling of AGN sizes by simulations such as the Tiered Radio Extragalactic Continuum Simulation (T-RECS). Finally, when considering the multi-tracer approach, where galaxies are split by type to measure some bias ratio, we find that there are not enough AGN present to achieve a reasonable level of shot noise for this kind of measurement. △ Less

Submitted 4 June, 2021; originally announced June 2021.

Comments: 10 pages, 8 figures, submitted to the journal

KiDS-1000 Cosmology: constraints beyond flat $Λ$CDM

Authors: Tilman Tröster, Marika Asgari, Chris Blake, Matteo Cataneo, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Chieh-An Lin, Ariel G. Sánchez, Angus H. Wright, Maciej Bilicki, Benjamin Bose, Martin Crocce, Andrej Dvornik, Thomas Erben, Benjamin Giblin, Karl Glazebrook, Henk Hoekstra, Shahab Joudaki, Arun Kannawadi, Fabian Köhlinger, Konrad Kuijken, Chris Lidman, Lucas Lombriser, Alexander Mead , et al. (4 additional authors not shown)

Abstract: We present constraints on extensions to the flat $Λ$CDM cosmological model by varying the spatial curvature $Ω_K$, the sum of the neutrino masses $\sum m_ν$, the dark energy equation of state parameter $w$, and the Hu-Sawicki $f(R)$ gravity $f_{R0}$ parameter. With the combined $3\times2$pt measurements of cosmic shear from the Kilo-Degree Survey (KiDS-1000), galaxy clustering from the Baryon Osci… ▽ More We present constraints on extensions to the flat $Λ$CDM cosmological model by varying the spatial curvature $Ω_K$, the sum of the neutrino masses $\sum m_ν$, the dark energy equation of state parameter $w$, and the Hu-Sawicki $f(R)$ gravity $f_{R0}$ parameter. With the combined $3\times2$pt measurements of cosmic shear from the Kilo-Degree Survey (KiDS-1000), galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing from the overlap between KiDS-1000, BOSS, and the spectroscopic 2-degree Field Lensing Survey (2dFLenS), we find results that are fully consistent with a flat $Λ$CDM model with $Ω_K=0.011^{+0.054}_{-0.057}$, $\sum m_ν<1.76$ eV (95% CL), and $w=-0.99^{+0.11}_{-0.13}$. The $f_{R0}$ parameter is unconstrained in our fully non-linear $f(R)$ cosmic shear analysis. Considering three different model selection criteria, we find no clear preference for either the fiducial flat $Λ$CDM model or any of the considered extensions. Besides extensions to the flat $Λ$CDM parameter space, we also explore restrictions to common subsets of the flat $Λ$CDM parameter space by fixing the amplitude of the primordial power spectrum to the Planck best-fit value, as well as adding external data from supernovae and lensing of the CMB. Neither the beyond-$Λ$CDM models nor the imposed restrictions explored in this analysis are able to resolve the $\sim 3σ$ tension in $S_8$ between the $3\times2$pt constraints and Planck, with the exception of $w$CDM, where the $S_8$ tension is resolved. The tension in the $w$CDM case persists, however, when considering the joint $S_8$-$w$ parameter space. The joint flat $Λ$CDM CMB lensing and $3\times2$pt analysis is found to yield tight constraints on $Ω_{\rm m}=0.307^{+0.008}_{-0.013}$, $σ_8=0.769^{+0.022}_{-0.010}$, and $S_8=0.779^{+0.013}_{-0.013}$. △ Less

Submitted 2 March, 2021; v1 submitted 30 October, 2020; originally announced October 2020.

Comments: 15 pages, 9 figures, 3 tables, accepted in A&A. This paper concludes the KiDS-1000 series of papers: Heymans, Tröster et al. (arXiv:2007.15632), Asgari et al. (arXiv:2007.15633), Hildebrandt et al. (arXiv:2007.15635), Joachimi et al. (arXiv:2007.01844), and Giblin et al. (arXiv:2007.01845). Data products, likelihoods, and posteriors can be found http://kids.strw.leidenuniv.nl/DR4/lensing.php

Journal ref: A&A 649, A88 (2021)

Cosmology with the Einstein Telescope: No Slip Gravity Model and Redshift Specifications

Authors: Ayan Mitra, Jurgen Mifsud, David F. Mota, David Parkinson

Abstract: The Einstein Telescope and other third generation interferometric detectors of gravitational waves are projected to be operational post $2030$. The cosmological signatures of gravitational waves would undoubtedly shed light on any departure from the current gravitational framework. We here confront a specific modified gravity model, the No Slip Gravity model, with forecast observations of gravitat… ▽ More The Einstein Telescope and other third generation interferometric detectors of gravitational waves are projected to be operational post $2030$. The cosmological signatures of gravitational waves would undoubtedly shed light on any departure from the current gravitational framework. We here confront a specific modified gravity model, the No Slip Gravity model, with forecast observations of gravitational waves. We compare the predicted constraints on the dark energy equation of state parameters $w_0^{}-w_a^{}$, between the modified gravity model and that of Einstein gravity. We show that the No Slip Gravity model mimics closely the constraints from the standard gravitational theory, and that the cosmological constraints are very similar. The use of spectroscopic redshifts, especially in the low--redshift regime, lead to significant improvements in the inferred parameter constraints. We test how well such a prospective gravitational wave dataset would function at testing such models, and find that there are significant degeneracies between the modified gravity model parameters, and the cosmological parameters that determine the distance, due to the gravitational wave dimming effect of the modified theory. △ Less

Submitted 30 September, 2020; originally announced October 2020.

Comments: 15 pages, 11 figures

Constrained Non-Linear Phase Retrieval for Single Distance X-ray Phase Contrast Tomography

Authors: K. Aditya Mohan, Dilworth Y. Parkinson, Jefferson A. Cuadra

Abstract: X-ray phase contrast tomography (XPCT) is widely used for 3D imaging of objects with weak contrast in X-ray absorption index but strong contrast in refractive index decrement. To reconstruct an object imaged using XPCT, phase retrieval algorithms are first used to estimate the X-ray phase projections, which is the 2D projection of the refractive index decrement, at each view. Phase retrieval is fo… ▽ More X-ray phase contrast tomography (XPCT) is widely used for 3D imaging of objects with weak contrast in X-ray absorption index but strong contrast in refractive index decrement. To reconstruct an object imaged using XPCT, phase retrieval algorithms are first used to estimate the X-ray phase projections, which is the 2D projection of the refractive index decrement, at each view. Phase retrieval is followed by refractive index decrement reconstruction from the phase projections using an algorithm such as filtered back projection (FBP). In practice, phase retrieval is most commonly solved by approximating it as a linear inverse problem. However, this linear approximation often results in artifacts and blurring when the conditions for the approximation are violated. In this paper, we formulate phase retrieval as a non-linear inverse problem, where we solve for the transmission function, which is the negative exponential of the projections, from XPCT measurements. We use a constraint to enforce proportionality between phase and absorption projections. We do not use constraints such as large Fresnel number, slowly varying phase, or Born/Rytov approximations. Our approach also does not require any regularization parameter tuning since there is no explicit sparsity enforcing regularization function. We validate the performance of our non-linear phase retrieval (NLPR) method using both simulated and real synchrotron datasets. We compare NLPR with a popular linear phase retrieval (LPR) approach and show that NLPR achieves sharper reconstructions with higher quantitative accuracy. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: 2020 IS&T International Symposium on Electronic Imaging

Journal ref: Electronic Imaging, no. 14, pp. 146-1-146-8, 2020

KiDS-1000 Cosmology: Multi-probe weak gravitational lensing and spectroscopic galaxy clustering constraints

Authors: Catherine Heymans, Tilman Tröster, Marika Asgari, Chris Blake, Hendrik Hildebrandt, Benjamin Joachimi, Konrad Kuijken, Chieh-An Lin, Ariel G. Sánchez, Jan Luca van den Busch, Angus H. Wright, Alexandra Amon, Maciej Bilicki, Jelte de Jong, Martin Crocce, Andrej Dvornik, Thomas Erben, Maria Cristina Fortuna, Fedor Getman, Benjamin Giblin, Karl Glazebrook, Henk Hoekstra, Shahab Joudaki, Arun Kannawadi, Fabian Köhlinger , et al. (9 additional authors not shown)

Abstract: We present a joint cosmological analysis of weak gravitational lensing observations from the Kilo-Degree Survey (KiDS-1000), with redshift-space galaxy clustering observations from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing observations from the overlap between KiDS-1000, BOSS and the spectroscopic 2-degree Field Lensing Survey (2dFLenS). This combination of larg… ▽ More We present a joint cosmological analysis of weak gravitational lensing observations from the Kilo-Degree Survey (KiDS-1000), with redshift-space galaxy clustering observations from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing observations from the overlap between KiDS-1000, BOSS and the spectroscopic 2-degree Field Lensing Survey (2dFLenS). This combination of large-scale structure probes breaks the degeneracies between cosmological parameters for individual observables, resulting in a constraint on the structure growth parameter $S_8=σ_8 \sqrt{Ω_{\rm m}/0.3} = 0.766^{+0.020}_{-0.014}$, that has the same overall precision as that reported by the full-sky cosmic microwave background observations from Planck. The recovered $S_8$ amplitude is low, however, by $8.3 \pm 2.6$ % relative to Planck. This result builds from a series of KiDS-1000 analyses where we validate our methodology with variable depth mock galaxy surveys, our lensing calibration with image simulations and null-tests, and our optical-to-near-infrared redshift calibration with multi-band mock catalogues and a spectroscopic-photometric clustering analysis. The systematic uncertainties identified by these analyses are folded through as nuisance parameters in our cosmological analysis. Inspecting the offset between the marginalised posterior distributions, we find that the $S_8$-difference with Planck is driven by a tension in the matter fluctuation amplitude parameter, $σ_8$. We quantify the level of agreement between the CMB and our large-scale structure constraints using a series of different metrics, finding differences with a significance ranging between $\sim\! 3\,σ$, when considering the offset in $S_{8}$, and $\sim\! 2\,σ$, when considering the full multi-dimensional parameter space. △ Less

Submitted 8 December, 2020; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 25 pages, 11 figures, 5 tables, A&A accepted, including a new appendix on Intrinsic Alignments. The KiDS-1000 data products are available for download at http://kids.strw.leidenuniv.nl/DR4/lensing.php. This data release includes open source software, the shear-photo-z catalogue, the cosmic shear and 3x2pt data vectors and covariances, and posteriors in the form of Multinest chains

Journal ref: A&A 646, A140 (2021)

Discovery of a Radio Relic in the Massive Merging Cluster SPT-CL 2023-5535 from the ASKAP-EMU PILOT SURVEY

Authors: Kim HyeongHan, M. James Jee, Lawrence Rudnick, David Parkinson, Kyle Finner, Mijin Yoon, Wonki Lee, Gianfranco Brunetti, Marcus Brüggen, Jordan D. Collier, Andrew M. Hopkins, Michał J. Michałowski, Ray P. Norris, Chris Riseley

Abstract: The ASKAP-EMU survey is a deep wide-field radio continuum survey designed to cover the entire southern sky and a significant fraction of the northern sky up to $+30^{\circ}$. Here, we report a discovery of a radio relic in the merging cluster SPT-CL 2023-5535 at z=0.23 from the ASKAP-EMU pilot 300 sq. deg survey (800-1088 MHz). The deep high-resolution data reveal a $\sim2$ Mpc-scale radio halo el… ▽ More The ASKAP-EMU survey is a deep wide-field radio continuum survey designed to cover the entire southern sky and a significant fraction of the northern sky up to $+30^{\circ}$. Here, we report a discovery of a radio relic in the merging cluster SPT-CL 2023-5535 at z=0.23 from the ASKAP-EMU pilot 300 sq. deg survey (800-1088 MHz). The deep high-resolution data reveal a $\sim2$ Mpc-scale radio halo elongated in the east-west direction, coincident with the intracluster gas. The radio relic is located at the western edge of this radio halo stretched $\sim0.5$ Mpc in the north-south orientation. The integrated spectral index of the radio relic within the narrow bandwidth is $α^{\scriptstyle \rm 1088~MHz}_{\scriptstyle \rm 800~MHz}=-0.76 \pm 0.06$. Our weak-lensing analysis shows that the system is massive ($M_{200}=1.04\pm0.36\times 10^{15} M_{\odot}$) and composed of at least three subclusters. We suggest a scenario, wherein the radio features arise from the collision between the eastern and middle subclusters. Our discovery illustrates the effectiveness of the ASKAP-EMU survey in detecting diffuse emissions in galaxy clusters and when completed, the survey will greatly increase the number of merging cluster detections with diffuse radio emissions. △ Less

Submitted 17 July, 2020; v1 submitted 16 July, 2020; originally announced July 2020.

Comments: Accepted to ApJ

HIR4: Cosmological signatures imprinted on the cross correlation between 21-cm map and galaxy clustering

Authors: Feng Shi, Yong-Seon Song, Jacobo Asorey, David Parkinson, Kyungjin Ahn, Jian Yao, Le Zhang, Shifan Zuo

Abstract: We explore the cosmological multitracer synergies between an emission line galaxy distribution from the Dark Energy Spectroscopic Instrument and a Tianlai Project 21-cm intensity map. We use simulated maps generated from a particle simulation in the light-cone volume (Horizon Run 4), sky-trimmed and including the effects of foreground contamination, its removal and instrument noise. We first valid… ▽ More We explore the cosmological multitracer synergies between an emission line galaxy distribution from the Dark Energy Spectroscopic Instrument and a Tianlai Project 21-cm intensity map. We use simulated maps generated from a particle simulation in the light-cone volume (Horizon Run 4), sky-trimmed and including the effects of foreground contamination, its removal and instrument noise. We first validate how the foreground residual affects the recovered 21-cm signal by putting different levels of foreground contamination into the 21-cm maps. We find that the contamination cannot be ignored in the angular autocorrelation power spectra of HI even when it is small, but has no influence on the accuracy of the angular cross-power spectra between HI and galaxies. In the foreground-cleaned map case, as information is lost in the cleaning procedure, there is also a bias in the cross-correlation power spectrum. However, we found that the bias from the cross-correlation power spectrum is scale-independent, which is easily parameterized as part of the model, while the offset in the HI autocorrelation power spectrum is non-linear. In particular, we tested that the cross-correlation power also benefits from the cancellation of the bias in the power spectrum measurement that is induced by the instrument noise, which changes the shape of the autocorrelation power spectra but leaves the cross-correlation power unaffected. We then modelled the angular cross-correlation power spectra to fit the baryon acoustic oscillation feature in the broad-band shape of the angular cross-correlation power spectrum, including contamination from the residual foreground and the effect of instrument noise. We forecast a constraint on the angular diameter distance $D_\mathrm{A}$ for the Tianlai Pathfinder redshift $0.775<z<1.03$, giving a distance measurement with a precision of 2.7\% at that redshift. △ Less

Submitted 16 November, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: Published in MNRAS. 14 pages, 10 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 499, Issue 4, December 2020, Pages 4613-4625

OzDES multi-object fibre spectroscopy for the Dark Energy Survey: Results and second data release

Authors: C. Lidman, B. E. Tucker, T. M. Davis, S. A. Uddin, J. Asorey, K. Bolejko, D. Brout, J. Calcino, D. Carollo, A. Carr, M. Childress, J. K. Hoormann, R. J. Foley, L. Galbany, K. Glazebrook, S. R. Hinton, R. Kessler, A. G. Kim, A. King, A. Kremin, K. Kuehn, D. Lagattuta, G. F. Lewis, E. Macaulay, U. Malik , et al. (79 additional authors not shown)

Abstract: We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the… ▽ More We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30,000 sources, some as faint as $r_{\mathrm AB}=24$ mag, and 375,000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to $z\sim1.2$ and the masses of hundreds of black holes out to $z\sim4$. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs. △ Less

Submitted 31 May, 2020; originally announced June 2020.

Comments: 20 pages, 16 figures. Accepted for publication by MNRAS. Data release available at https://datacentral.org.au

Report number: Report-no: FERMILAB-PUB-20-029-AE; DES 2019-0517

Testing gravity using galaxy-galaxy lensing and clustering amplitudes in KiDS-1000, BOSS and 2dFLenS

Authors: Chris Blake, Alexandra Amon, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Thomas Erben, Benjamin Giblin, Karl Glazebrook, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Shahab Joudaki, Arun Kannawadi, Konrad Kuijken, Chris Lidman, David Parkinson, HuanYuan Shan, Tilman Tröster, Jan Luca van den Busch, Christian Wolf, Angus H. Wright

Abstract: The physics of gravity on cosmological scales affects both the rate of assembly of large-scale structure, and the gravitational lensing of background light through this cosmic web. By comparing the amplitude of these different observational signatures, we can construct tests that can distinguish general relativity from its potential modifications. We used the latest weak gravitational lensing data… ▽ More The physics of gravity on cosmological scales affects both the rate of assembly of large-scale structure, and the gravitational lensing of background light through this cosmic web. By comparing the amplitude of these different observational signatures, we can construct tests that can distinguish general relativity from its potential modifications. We used the latest weak gravitational lensing dataset from the Kilo-Degree Survey, KiDS-1000, in conjunction with overlapping galaxy spectroscopic redshift surveys BOSS and 2dFLenS, to perform the most precise existing amplitude-ratio test. We measured the associated E_G statistic with 15-20% errors, in five dz = 0.1 tomographic redshift bins in the range 0.2 < z < 0.7, on projected scales up to 100 Mpc/h. The scale-independence and redshift-dependence of these measurements are consistent with the theoretical expectation of general relativity in a Universe with matter density Omega_m = 0.27 +/- 0.04. We demonstrate that our results are robust against different analysis choices, including schemes for correcting the effects of source photometric redshift errors, and compare the performance of angular and projected galaxy-galaxy lensing statistics. △ Less

Submitted 27 August, 2020; v1 submitted 28 May, 2020; originally announced May 2020.

Comments: 23 pages, 14 figures, version accepted for publication by A&A

Journal ref: A&A 642, A158 (2020)

Testing the violation of the equivalence principle in the electromagnetic sector and its consequences in $f(T)$ gravity

Authors: Jackson Levi Said, Jurgen Mifsud, David Parkinson, Emmanuel N. Saridakis, Joseph Sultana, Kristian Zarb Adami

Abstract: A violation of the distance-duality relation is directly linked with a temporal variation of the electromagnetic fine-structure constant. We consider a number of well-studied $f(T)$ gravity models and we revise the theoretical prediction of their corresponding induced violation of the distance-duality relationship. We further extract constraints on the involved model parameters through fine-struct… ▽ More A violation of the distance-duality relation is directly linked with a temporal variation of the electromagnetic fine-structure constant. We consider a number of well-studied $f(T)$ gravity models and we revise the theoretical prediction of their corresponding induced violation of the distance-duality relationship. We further extract constraints on the involved model parameters through fine-structure constant variation data, alongside with supernovae data, and Hubble parameter measurements. Moreover, we constrain the evolution of the effective $f(T)$ gravitational constant. Finally, we compare with revised constraints on the phenomenological parametrisations of the violation of the equivalence principle in the electromagnetic sector. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Sparse Matrix-Based HPC Tomography

Authors: Stefano Marchesini, Anuradha Trivedi, Pablo Enfedaque, Talita Perciano, Dilworth Parkinson

Abstract: Tomographic imaging has benefited from advances in X-ray sources, detectors and optics to enable novel observations in science, engineering and medicine. These advances have come with a dramatic increase of input data in the form of faster frame rates, larger fields of view or higher resolution, so high performance solutions are currently widely used for analysis. Tomographic instruments can vary… ▽ More Tomographic imaging has benefited from advances in X-ray sources, detectors and optics to enable novel observations in science, engineering and medicine. These advances have come with a dramatic increase of input data in the form of faster frame rates, larger fields of view or higher resolution, so high performance solutions are currently widely used for analysis. Tomographic instruments can vary significantly from one to another, including the hardware employed for reconstruction: from single CPU workstations to large scale hybrid CPU/GPU supercomputers. Flexibility on the software interfaces and reconstruction engines are also highly valued to allow for easy development and prototyping. This paper presents a novel software framework for tomographic analysis that tackles all aforementioned requirements. The proposed solution capitalizes on the increased performance of sparse matrix-vector multiplication and exploits multi-CPU and GPU reconstruction over MPI. The solution is implemented in Python and relies on CuPy for fast GPU operators and CUDA kernel integration, and on SciPy for CPU sparse matrix computation. As opposed to previous tomography solutions that are tailor-made for specific use cases or hardware, the proposed software is designed to provide flexible, portable and high-performance operators that can be used for continuous integration at different production environments, but also for prototyping new experimental settings or for algorithmic development. The experimental results demonstrate how our implementation can even outperform state-of-the-art software packages used at advanced X-ray sources worldwide. △ Less

Submitted 22 April, 2020; v1 submitted 27 March, 2020; originally announced March 2020.

Atmospheric Escape Processes and Planetary Atmospheric Evolution

Authors: Guillaume Gronoff, Phil Arras, Suleiman M Baraka, Jared M Bell, Gaël Cessateur, Ofer Cohen, Shannon M. Curry, Jeremy J Drake, Meredith K Elrod, Justin T. Erwin, Katherine Garcia-Sage, Cecilia Garraffo, Alex Glocer, Nicholas Gray Heavens, Kylie Lovato, Romain Maggiolo, Christopher D. Parkinson, Cyril L. Simon Wedlund, Daniel R Weimer, William B. Moore

Abstract: The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to und… ▽ More The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to understand the broad range of processes that lead to atmospheric escape and evolution under a wide range of conditions if we are to assess the habitability of worlds around other stars. One problem encountered between the planetary and the astrophysics communities is a lack of common language for describing escape processes. Each community has customary approximations that may be questioned by the other, such as the hypothesis of H-dominated thermosphere for astrophysicists, or the Sun-like nature of the stars for planetary scientists. Since exoplanets are becoming one of the main targets for the detection of life, a common set of definitions and hypotheses are required. We review the different escape mechanisms proposed for the evolution of planetary and exoplanetary atmospheres. We propose a common definition for the different escape mechanisms, and we show the important parameters to take into account when evaluating the escape at a planet in time. We show that the paradigm of the magnetic field as an atmospheric shield should be changed and that recent work on the history of Xenon in Earth's atmosphere gives an elegant explanation to its enrichment in heavier isotopes: the so-called Xenon paradox. △ Less

Submitted 6 March, 2020; originally announced March 2020.

HIR4: cosmology from a simulated neutral hydrogen full sky using Horizon Run 4

Authors: Jacobo Asorey, David Parkinson, Feng Shi, Yong-Seon Song, Kyungjin Ahn, Juhan Kim, Jian Yao, Le Zhang, Shifan Zuo

Abstract: The distribution of cosmological neutral hydrogen will provide a new window into the large-scale structure of the Universe with the next generation of radio telescopes and surveys. The observation of this material, through 21cm line emission, will be confused by foreground emission in the same frequencies. Even after these foregrounds are removed, the reconstructed map may not exactly match the or… ▽ More The distribution of cosmological neutral hydrogen will provide a new window into the large-scale structure of the Universe with the next generation of radio telescopes and surveys. The observation of this material, through 21cm line emission, will be confused by foreground emission in the same frequencies. Even after these foregrounds are removed, the reconstructed map may not exactly match the original cosmological signal, which will introduce systematic errors and offset into the measured correlations. In this paper, we simulate future surveys of neutral hydrogen using the Horizon Run 4 (HR4) cosmological N-body simulation. We generate HI intensity maps from the HR4 halo catalogue, and combine with foreground radio emission maps from the Global Sky Model, to create accurate simulations over the entire sky. We simulate the HI sky for the frequency range 700-800 MHz, matching the sensitivity of the Tianlai pathfinder. We test the accuracy of the fastICA, PCA and log-polynomial fitting foreground removal methods to recover the input cosmological angular power spectrum and measure the parameters. We show the effect of survey noise levels and beam sizes on the recovered the cosmological constraints. We find that while the reconstruction removes power from the cosmological 21cm distribution on large-scales, we can correct for this and recover the input parameters in the noise-free case. However, the effect of noise and beam size of the Tianlai pathfinder prevents accurate recovery of the cosmological parameters when using only intensity mapping information. △ Less

Submitted 27 August, 2020; v1 submitted 3 January, 2020; originally announced January 2020.

Comments: Published in MNRAS. 21 pages, 20 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 495, Issue 2, June 2020, pp.1788-1806

The SPTpol Extended Cluster Survey

Authors: L. E. Bleem, S. Bocquet, B. Stalder, M. D. Gladders, P. A. R. Ade, S. W. Allen, A. J. Anderson, J. Annis, M. L. N. Ashby, J. E. Austermann, S. Avila, J. S. Avva, M. Bayliss, J. A. Beall, K. Bechtol, A. N. Bender, B. A. Benson, E. Bertin, F. Bianchini, C. Blake, M. Brodwin, D. Brooks, E. Buckley-Geer, D. L. Burke, J. E. Carlstrom , et al. (113 additional authors not shown)

Abstract: We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 c… ▽ More We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $ξ\ge 5$ and an additional 204 systems at $4<ξ<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $ξ>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($λ-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $σ$ level---with the relations intersecting at $λ=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses. △ Less

Submitted 13 December, 2019; v1 submitted 9 October, 2019; originally announced October 2019.

Comments: 49 pages, 14 figures, 10 tables. Minor changes to match accepted version in ApJS

Developing a unified pipeline for large-scale structure data analysis with angular power spectra -- II. A case study for magnification bias and radio continuum surveys

Authors: Konstantinos Tanidis, Stefano Camera, David Parkinson

Abstract: Following on our purpose of developing a unified pipeline for large-scale structure data analysis with angular (i.e. harmonic-space) power spectra, we now include the weak lensing effect of magnification bias on galaxy clustering in a publicly available, modular parameter estimation code. We thus forecast constraints on the parameters of the concordance cosmological model, dark energy, and modifie… ▽ More Following on our purpose of developing a unified pipeline for large-scale structure data analysis with angular (i.e. harmonic-space) power spectra, we now include the weak lensing effect of magnification bias on galaxy clustering in a publicly available, modular parameter estimation code. We thus forecast constraints on the parameters of the concordance cosmological model, dark energy, and modified gravity theories from galaxy clustering tomographic angular power spectra. We find that a correct modelling of magnification is crucial in order not to bias the estimation of cosmological parameters, especially in the case of deep galaxy surveys. Our case study adopts specifications of the Evolutionary Map of the Universe (EMU), which is a full-sky, deep radio-continuum survey, and is expected to probe the Universe up to redshift $z\sim6$. We assume the Limber approximation, and include magnification bias on top of density fluctuations and redshift-space distortions. By restricting our analysis to the regime where the Limber approximation holds true, we significantly minimise the computational time needed, compared to that of the exact calculation. We also show that there is a trend for more biased parameter estimates from neglecting magnification when the redshift bins are very wide. We conclude that this result implies a strong dependence on the lensing contribution, which is an integrated effect and becomes dominant when wide redshift bins are considered. Finally, we note that instead of being considered a contaminant, magnification bias encodes important cosmological information, and its inclusion leads to an alleviation of the degeneracy between the galaxy bias and the amplitude normalisation of the matter fluctuations. △ Less

Submitted 24 September, 2020; v1 submitted 23 September, 2019; originally announced September 2019.

Comments: 16 pages, 6 figures, 12 tables. Version matching publication at journal

Journal ref: MNRAS 491, 4869 (2020)

Research and Development for HI Intensity Mapping

Authors: Zeeshan Ahmed, David Alonso, Mustafa A. Amin, Réza Ansari, Evan J. Arena, Kevin Bandura, Adam Beardsley, Philip Bull, Emanuele Castorina, Tzu-Ching Chang, Romeel Davé, Joshua S. Dillon, Alexander van Engelen, Aaron Ewall-Wice, Simone Ferraro, Simon Foreman, Josef Frisch, Daniel Green, Gilbert Holder, Daniel Jacobs, Dionysios Karagiannis, Alexander A. Kaurov, Lloyd Knox, Emily Kuhn, Adrian Liu , et al. (29 additional authors not shown)

Abstract: Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and an… ▽ More Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift. △ Less

Submitted 29 July, 2019; originally announced July 2019.

Comments: 10 pages + references, 2 figures, 1 table; APC white paper submitted to the Astro2020 decadal survey. arXiv admin note: text overlap with arXiv:1810.09572

DASH: Deep Learning for the Automated Spectral Classification of Supernovae and their Hosts

Authors: Daniel Muthukrishna, David Parkinson, Brad Tucker

Abstract: We present DASH (Deep Automated Supernova and Host classifier), a novel software package that automates the classification of the type, age, redshift, and host galaxy of supernova spectra. DASH makes use of a new approach that does not rely on iterative template matching techniques like all previous software, but instead classifies based on the learned features of each supernova's type and age. It… ▽ More We present DASH (Deep Automated Supernova and Host classifier), a novel software package that automates the classification of the type, age, redshift, and host galaxy of supernova spectra. DASH makes use of a new approach that does not rely on iterative template matching techniques like all previous software, but instead classifies based on the learned features of each supernova's type and age. It has achieved this by employing a deep convolutional neural network to train a matching algorithm. This approach has enabled DASH to be orders of magnitude faster than previous tools, being able to accurately classify hundreds or thousands of objects within seconds. We have tested its performance on four years of data from the Australian Dark Energy Survey (OzDES). The deep learning models were developed using TensorFlow, and were trained using over 4000 supernova spectra taken from the CfA Supernova Program and the Berkeley SN Ia Program as used in SNID (Supernova Identification software, Blondin & Tonry 2007). Unlike template matching methods, the trained models are independent of the number of spectra in the training data, which allows for DASH's unprecedented speed. We have developed both a graphical interface for easy visual classification and analysis of supernovae, and a Python library for the autonomous and quick classification of several supernova spectra. The speed, accuracy, user-friendliness, and versatility of DASH presents an advancement to existing spectral classification tools. We have made the code publicly available on GitHub and PyPI (pip install astrodash) to allow for further contributions and development. The package documentation is available at https://astrodash.readthedocs.io. △ Less

Submitted 7 October, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: 15 pages plus appendices, 8 figures, 2 tables. Accepted for publication in the Astrophysical Journal

Journal ref: ApJ 885, 85 (2019)

4MOST: Project overview and information for the First Call for Proposals

Authors: R. S. de Jong, O. Agertz, A. Agudo Berbel, J. Aird, D. A. Alexander, A. Amarsi, F. Anders, R. Andrae, B. Ansarinejad, W. Ansorge, P. Antilogus, H. Anwand-Heerwart, A. Arentsen, A. Arnadottir, M. Asplund, M. Auger, N. Azais, D. Baade, G. Baker, S. Baker, E. Balbinot, I. K. Baldry, M. Banerji, S. Barden, P. Barklem , et al. (313 additional authors not shown)

Abstract: We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolut… ▽ More We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = λ/Δλ\sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations. △ Less

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 3

Cosmology with Phase 1 of the Square Kilometre Array; Red Book 2018: Technical specifications and performance forecasts

Authors: Square Kilometre Array Cosmology Science Working Group, David J. Bacon, Richard A. Battye, Philip Bull, Stefano Camera, Pedro G. Ferreira, Ian Harrison, David Parkinson, Alkistis Pourtsidou, Mario G. Santos, Laura Wolz, Filipe Abdalla, Yashar Akrami, David Alonso, Sambatra Andrianomena, Mario Ballardini, Jose Luis Bernal, Daniele Bertacca, Carlos A. P. Bengaly, Anna Bonaldi, Camille Bonvin, Michael L. Brown, Emma Chapman, Song Chen, Xuelei Chen , et al. (22 additional authors not shown)

Abstract: We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z… ▽ More We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1. △ Less

Submitted 6 November, 2018; originally announced November 2018.

Comments: Red Book 2018 of the Square Kilometre Array Cosmology Science Working Group; 35 pages, 27 figures; To be submitted to PASA

Journal ref: Publ. Astron. Soc. Austral. 37 (2020) e007

Probing $Λ$CDM cosmology with the Evolutionary Map of the Universe survey

Authors: José Luis Bernal, Alvise Raccanelli, Ely D. Kovetz, David Parkinson, Ray P. Norris, George Danforth, Courtney Schmitt

Abstract: The Evolutionary Map of the Universe (EMU) is an all-sky survey in radio-continuum which uses the Australian SKA Pathfinder (ASKAP). Using galaxy angular power spectrum and the integrated Sachs-Wolfe effect, we study the potential of EMU to constrain models beyond $Λ$CDM (i.e., local primordial non-Gaussianity, dynamical dark energy, spatial curvature and deviations from general relativity), for d… ▽ More The Evolutionary Map of the Universe (EMU) is an all-sky survey in radio-continuum which uses the Australian SKA Pathfinder (ASKAP). Using galaxy angular power spectrum and the integrated Sachs-Wolfe effect, we study the potential of EMU to constrain models beyond $Λ$CDM (i.e., local primordial non-Gaussianity, dynamical dark energy, spatial curvature and deviations from general relativity), for different design sensitivities. We also include a multi-tracer analysis, distinguishing between star-forming galaxies and galaxies with an active galactic nucleus, to further improve EMU's potential. We find that EMU could measure the dark energy equation of state parameters around 35\% more precisely than existing constraints, and that the constraints on $f_{\rm NL}$ and modified gravity parameters will improve up to a factor $\sim2$ with respect to Planck and redshift space distortions measurements. With this work we demonstrate the promising potential of EMU to contribute to our understanding of the Universe. △ Less

Submitted 18 February, 2019; v1 submitted 15 October, 2018; originally announced October 2018.

Comments: 15 pages (29 with references and appendices), 6 figures and 10 tables. Matches the published version. Minimal changes from previous version

Report number: CERN-TH-2018-223

Journal ref: JCAP02(2019)030

Fundamental Physics with the Square Kilometre Array

Authors: A. Weltman, P. Bull, S. Camera, K. Kelley, H. Padmanabhan, J. Pritchard, A. Raccanelli, S. Riemer-Sørensen, L. Shao, S. Andrianomena, E. Athanassoula, D. Bacon, R. Barkana, G. Bertone, C. Bonvin, A. Bosma, M. Brüggen, C. Burigana, C. Bœhm, F. Calore, J. A. R. Cembranos, C. Clarkson, R. M. T. Connors, Á. de la Cruz-Dombriz, P. K. S. Dunsby , et al. (28 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA a… ▽ More The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime. △ Less

Submitted 22 October, 2019; v1 submitted 5 October, 2018; originally announced October 2018.

Comments: 69 pages, PASA, in press

Journal ref: Publ. Astron. Soc. Aust. 37 (2020) e002

Impact of bias and redshift-space modelling for the halo power spectrum: Testing the effective field theory of large-scale structure

Authors: Lucia Fonseca de la Bella, Donough Regan, David Seery, David Parkinson

Abstract: We study the impact of different bias and redshift-space models on the halo power spectrum, quantifying their effect by comparing the fit to a subset of realizations taken from the WizCOLA suite. These provide simulated power spectrum measurements between $k_{\rm min}$ = 0.03 h/Mpc and $k_{\rm max}$ = 0.29 h/Mpc, constructed using the comoving Lagrangian acceleration method. For the bias prescript… ▽ More We study the impact of different bias and redshift-space models on the halo power spectrum, quantifying their effect by comparing the fit to a subset of realizations taken from the WizCOLA suite. These provide simulated power spectrum measurements between $k_{\rm min}$ = 0.03 h/Mpc and $k_{\rm max}$ = 0.29 h/Mpc, constructed using the comoving Lagrangian acceleration method. For the bias prescription we include (i) simple linear bias; (ii) the McDonald & Roy model and (iii) its coevolution variant introduced by Saito et al.; and (iv) a very general model including all terms up to one-loop and corrections from advection. For the redshift-space modelling we include the Kaiser formula with exponential damping and the power spectrum provided by (i) tree-level perturbation theory and (ii) the Halofit prescription; (iii) one-loop perturbation theory, also with exponential damping; and (iv) an effective field theory description, also at one-loop, with damping represented by the EFT subtractions. We quantify the improvement from each layer of modelling by measuring the typical improvement in chi-square when fitting to a member of the simulation suite. We attempt to detect overfitting by testing for compatibility between the best-fit power spectrum per realization and the best-fit over the entire WizCOLA suite. For both bias and the redshift-space map we find that increasingly permissive models yield improvements in chi-square but with diminishing returns. The most permissive models show modest evidence for overfitting. Accounting for model complexity using the Bayesian Information Criterion, we argue that standard perturbation theory up to one-loop, or a related model such as that of Taruya, Nishimichi & Saito, coupled to the coevolution bias model, is likely to provide a good compromise for near-future galaxy surveys operating with comparable $k_{\rm max}$. △ Less

Submitted 28 May, 2020; v1 submitted 31 May, 2018; originally announced May 2018.

Comments: 34 pages, 8 figures. Appendices give further details of the one-loop SPT calculations

Journal ref: JCAP07(2020)011

Identifying complex sources in large astronomical data using a coarse-grained complexity measure

Authors: Gary Segal, David Parkinson, Ray P. Norris, Jesse Swan

Abstract: The volume of data that will be produced by the next generation of astrophysical instruments represents a significant opportunity for making unplanned and unexpected discoveries. Conversely, finding unexpected objects or phenomena within such large volumes of data presents a challenge that may best be solved using computational and statistical approaches. We present the application of a coarse-gra… ▽ More The volume of data that will be produced by the next generation of astrophysical instruments represents a significant opportunity for making unplanned and unexpected discoveries. Conversely, finding unexpected objects or phenomena within such large volumes of data presents a challenge that may best be solved using computational and statistical approaches. We present the application of a coarse-grained complexity measure for identifying interesting observations in large astronomical data sets. This measure, which has been termed apparent complexity, has been shown to model human intuition and perceptions of complexity. Apparent complexity is computationally efficient to derive and can be used to segment and identify interesting observations in very large data sets based on their morphological complexity. We show, using data from the Australia Telescope Large Area Survey, that apparent complexity can be combined with clustering methods to provide an automated process for distinguishing between images of galaxies which have been classified as having simple and complex morphologies. The approach generalizes well when applied to new data after being calibrated on a smaller data set, where it performs better than tested classification methods using pixel data. This generalizability positions apparent complexity as a suitable machine-learning feature for identifying complex observations with unanticipated features. △ Less

Submitted 9 September, 2019; v1 submitted 27 May, 2018; originally announced May 2018.

Comments: 15 pages, 8 figures. New version to match that appearing in the journal. Appearing in PASP Focus on Machine Intelligence in Astronomy and Astrophysics

Journal ref: PASP, 131, 108007 (2019)

KiDS+2dFLenS+GAMA: Testing the cosmological model with the $E_{\rm G}$ statistic

Authors: A. Amon, C. Blake, C. Heymans, C. D. Leonard, M. Asgari, M. Bilicki, A. Choi, T. Erben, K. Glazebrook, J. Harnois-Deraps, H. Hildebrandt, H. Hoekstra, B. Joachimi, S. Joudaki, K. Kuijken, C. Lidman, J. Loveday, D. Parkinson, E. A. Valentijn, C. Wolf

Abstract: We present a new measurement of $E_{\rm G}$, which combines measurements of weak gravitational lensing, real-space galaxy clustering and redshift space distortions. This statistic was proposed as a consistency test of General Relativity (GR) that is insensitive to linear, deterministic galaxy bias and the matter clustering amplitude. We combine deep imaging data from KiDS with overlapping spectros… ▽ More We present a new measurement of $E_{\rm G}$, which combines measurements of weak gravitational lensing, real-space galaxy clustering and redshift space distortions. This statistic was proposed as a consistency test of General Relativity (GR) that is insensitive to linear, deterministic galaxy bias and the matter clustering amplitude. We combine deep imaging data from KiDS with overlapping spectroscopy from 2dFLenS, BOSS DR12 and GAMA and find $E_{\rm G}(\overline{z}=0.267)=0.43 \pm 0.13$ (GAMA), $E_{\rm G}(\overline{z}=0.305)=0.27 \pm 0.08$ (LOWZ+2dFLOZ) and $E_{\rm G}(\overline{z}=0.554)=0.26 \pm 0.07$ (CMASS+2dFHIZ). We demonstrate that the existing tension in the value of the matter density parameter hinders the robustness of this statistic as solely a test of GR. We find that our $E_{\rm G}$ measurements, as well as existing ones in the literature, favour a lower matter density cosmology than the Cosmic Microwave Background. For a flat $Λ$CDM Universe and assuming GR, we find $Ω_{\rm m}(z=0)=0.25\pm0.03$. With this paper we publicly release the 2dFLenS dataset at: \url{http://2dflens.swin.edu.au}. △ Less

Submitted 25 October, 2018; v1 submitted 29 November, 2017; originally announced November 2017.

Comments: 17 pages, 8 figures

Journal ref: 2018, MNRAS, 479, 3422

Photometric redshifts for the Kilo-Degree Survey. Machine-learning analysis with artificial neural networks

Authors: M. Bilicki, H. Hoekstra, M. J. I. Brown, V. Amaro, C. Blake, S. Cavuoti, J. T. A. de Jong, C. Georgiou, H. Hildebrandt, C. Wolf, A. Amon, M. Brescia, S. Brough, M. V. Costa-Duarte, T. Erben, K. Glazebrook, A. Grado, C. Heymans, T. Jarrett, S. Joudaki, K. Kuijken, G. Longo, N. Napolitano, D. Parkinson, C. Vellucci , et al. (2 additional authors not shown)

Abstract: We present a machine-learning photometric redshift analysis of the Kilo-Degree Survey Data Release 3, using two neural-network based techniques: ANNz2 and MLPQNA. Despite limited coverage of spectroscopic training sets, these ML codes provide photo-zs of quality comparable to, if not better than, those from the BPZ code, at least up to zphot<0.9 and r<23.5. At the bright end of r<20, where very co… ▽ More We present a machine-learning photometric redshift analysis of the Kilo-Degree Survey Data Release 3, using two neural-network based techniques: ANNz2 and MLPQNA. Despite limited coverage of spectroscopic training sets, these ML codes provide photo-zs of quality comparable to, if not better than, those from the BPZ code, at least up to zphot<0.9 and r<23.5. At the bright end of r<20, where very complete spectroscopic data overlapping with KiDS are available, the performance of the ML photo-zs clearly surpasses that of BPZ, currently the primary photo-z method for KiDS. Using the Galaxy And Mass Assembly (GAMA) spectroscopic survey as calibration, we furthermore study how photo-zs improve for bright sources when photometric parameters additional to magnitudes are included in the photo-z derivation, as well as when VIKING and WISE infrared bands are added. While the fiducial four-band ugri setup gives a photo-z bias $δz=-2e-4$ and scatter $σ_z<0.022$ at mean z = 0.23, combining magnitudes, colours, and galaxy sizes reduces the scatter by ~7% and the bias by an order of magnitude. Once the ugri and IR magnitudes are joined into 12-band photometry spanning up to 12 $μ$, the scatter decreases by more than 10% over the fiducial case. Finally, using the 12 bands together with optical colours and linear sizes gives $δz<4e-5$ and $σ_z<0.019$. This paper also serves as a reference for two public photo-z catalogues accompanying KiDS DR3, both obtained using the ANNz2 code. The first one, of general purpose, includes all the 39 million KiDS sources with four-band ugri measurements in DR3. The second dataset, optimized for low-redshift studies such as galaxy-galaxy lensing, is limited to r<20, and provides photo-zs of much better quality than in the full-depth case thanks to incorporating optical magnitudes, colours, and sizes in the GAMA-calibrated photo-z derivation. △ Less

Submitted 11 May, 2018; v1 submitted 13 September, 2017; originally announced September 2017.

Comments: A&A, in press. Data available from the KiDS website http://kids.strw.leidenuniv.nl/DR3/ml-photoz.php#annz2

Journal ref: A&A 616, A69 (2018)