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.

Recent ground-based deep observations of the Universe have discovered large populations of massive quiescent galaxies at $z\sim3\!-\!5$ . With the launch of the James Webb Space Telescope (JWST), the on-board Near-Infrared Spectrograph (NIRSpec) instrument will provide continuous $0.6\!-\!5.3\,\unicode{x03BC}\,\mathrm{m}$ spectroscopic coverage of these galaxies. Here we show that NIRSpec/CLEAR spectroscopy is ideal to probe the completeness of photometrically selected massive quiescent galaxies such as the ones presented by Schreiber et al. (2018b, A&A, 618, A85). Using a subset of the Schreiber et al. (2018b, A&A, 618, A85) sample with deep Keck/MOSFIRE spectroscopy presented by Esdaile J., et al. (2021b, ApJ, 908, L35), we perform a suite of mock JWST/NIRSpec observations to determine optimal observing strategies to efficiently recover the star formation histories (SFHs), element abundances, and kinematics of these massive quiescent galaxies. We find that at $z\sim3$ , medium resolution G235M/FL170LP NIRSpec observations could recover element abundances at an accuracy of ${\sim}15\%$ , which is comparable to local globular clusters. Mimicking ZFOURGE COSMOS photometry, we perform mock spectrophotometric fitting with Prospector to show that the overall shape of the SFHs of our mock galaxies can be recovered well, albeit with a dependency on the number of non-parametric SFH bins. We show that deep high-resolution G235H/FL170LP integral field spectroscopy with a $S/N\sim7$ per spaxel is required to constrain the rotational properties of our sample at $>\!2\sigma$ confidence. Thus, through optimal grism/filter choices, JWST/NIRSpec slit and integral field spectroscopy observations would provide tight constraints to galaxy evolution in the early Universe.

The science cases for incorporating high time resolution capabilities into modern radio telescopes are as numerous as they are compelling. Science targets range from exotic sources such as pulsars, to our Sun, to recently detected possible extragalactic bursts of radio emission, the so-called fast radio bursts (FRBs). Originally conceived purely as an imaging telescope, the initial design of the Murchison Widefield Array (MWA) did not include the ability to access high time and frequency resolution voltage data. However, the flexibility of the MWA’s software correlator allowed an off-the-shelf solution for adding this capability. This paper describes the system that records the 100 μs and 10 kHz resolution voltage data from the MWA. Example science applications, where this capability is critical, are presented, as well as accompanying commissioning results from this mode to demonstrate verification.

We present the most sensitive and detailed view of the neutral hydrogen ( ${\rm H\small I}$ ) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal ${\rm H\small I}$ in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K ( $1.6\,\mathrm{mJy\ beam}^{-1}$ ) $\mathrm{per}\ 0.98\,\mathrm{km\ s}^{-1}$ spectral channel with an angular resolution of $30^{\prime\prime}$ ( ${\sim}10\,\mathrm{pc}$ ). We discuss the calibration scheme and the custom imaging pipeline that utilises a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire ${\sim}25\,\mathrm{deg}^2$ field-of-view. We provide an overview of the data products and characterise several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high-velocity cloud with previous ASKAP+Parkes ${\rm H\small I}$ test observations.

Proposed next-generation networks of gravitational-wave observatories include dedicated kilohertz instruments that target neutron star science, such as the proposed Neutron Star Extreme Matter Observatory, NEMO. The original proposal for NEMO highlighted the need for it to exist in a network of gravitational-wave observatories to ensure detection confidence and sky localisation of sources. We show that NEMO-like observatories have significant utility on their own as coincident electromagnetic observations can provide the detection significance and sky localisation. We show that, with a single NEMO-like detector and expected electromagnetic observatories in the late 2020 s and early 2030 s such as the Vera C. Rubin observatory and SVOM, approximately 40% of all binary neutron star mergers detected with gravitational waves could be confidently identified as coincident multimessenger detections. We show that we expect $2^{+10}_{-1}{yr^{-1}}{}$ coincident observations of gravitational-wave mergers with gamma-ray burst prompt emission, $13^{+23}_{-10}{yr^{-1}}{}$ detections with kilonova observations, and $4^{+18}_{-3}{yr^{-1}}{}$ with broadband afterglows and kilonovae, where the uncertainties are 90% confidence intervals arising from uncertainty in current merger-rate estimates. Combined, this implies a coincident detection rate of $14^{+25}_{-11}{yr^{-1}}{}$ out to $300\,\mathrm{Mpc}$ . These numbers indicate significant science potential for a single kilohertz gravitational-wave detector operating without a global network of other gravitational-wave observatories.

We present a detailed overview of the cosmological surveys that we aim to carry 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 deg2; a wide and deep continuum galaxy and HI intensity mapping (IM) survey over 20 000 deg2 from $z = 0.35$ to 3; and a deep, high-redshift HI IM survey over 100 deg2 from $z = 3$ to 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to $z \sim 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 and near-infrared (NIR) surveys like Euclid, LSST, and WFIRST 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.

The Murchison Widefield Array (MWA) is one of three Square Kilometre Array Precursor telescopes and is located at the Murchison Radio-astronomy Observatory in the Murchison Shire of the mid-west of Western Australia, a location chosen for its extremely low levels of radio frequency interference. The MWA operates at low radio frequencies, 80–300 MHz, with a processed bandwidth of 30.72 MHz for both linear polarisations, and consists of 128 aperture arrays (known as tiles) distributed over a ~3-km diameter area. Novel hybrid hardware/software correlation and a real-time imaging and calibration systems comprise the MWA signal processing backend. In this paper, the as-built MWA is described both at a system and sub-system level, the expected performance of the array is presented, and the science goals of the instrument are summarised.

The Commensal Real-time Australian Square Kilometre Array Pathfinder Fast Transients survey is the first extensive astronomical survey using phased array feeds. Since January 2017, it has been searching for fast radio bursts in fly’s eye mode. Here, we present a calculation of the sensitivity and total exposure of the survey that detected the first 20 of these bursts, using the pulsars B1641-45 and B0833-45 as calibrators. The beamshape, antenna-dependent system noise, and the effects of radio-frequency interference and fluctuations during commissioning are quantified. Effective survey exposures and sensitivities are calculated as a function of the source counts distribution. Statistical ‘stat’ and systematics ‘sys’ effects are treated separately. The implied fast radio burst rate is significantly lower than the 37 sky−1 day−1 calculated using nominal exposures and sensitivities for this same sample by Shannon et al. (2018). At the Euclidean (best-fit) power-law index of −1.5 (−2.2), the rate is $12.7_{-2.2}^{+3.3}$ (sys) ± 3.6 (stat) sky−1 day−1 ( $20.7_{-1.7}^{+2.1}$ (sys) ± 2.8 (stat) sky−1 day−1) above a threshold of 56.6 ± 6.6(sys) Jy ms (40.4 ± 1.2(sys) Jy ms). This strongly suggests that these calculations be performed for other FRB-hunting experiments, allowing meaningful comparisons to be made between them.

The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the centre of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus on the mechanism which transforms the collapse into an explosion. Progress in understanding this mechanism is reviewed with particular attention to its asymmetric character. We highlight a series of successful studies connecting observations of supernova remnants and pulsars properties to the theory of core-collapse using numerical simulations. The encouraging results from first principles models in axisymmetric simulations is tempered by new puzzles in 3D. The diversity of explosion paths and the dependence on the pre-collapse stellar structure is stressed, as well as the need to gain a better understanding of hydrodynamical and MHD instabilities such as standing accretion shock instability and neutrino-driven convection. The shallow water analogy of shock dynamics is presented as a comparative system where buoyancy effects are absent. This dynamical system can be studied numerically and also experimentally with a water fountain. The potential of this complementary research tool for supernova theory is analysed. We also review its potential for public outreach in science museums.

This is an introduction to Bayesian inference with a focus on hierarchical models and hyper-parameters. We write primarily for an audience of Bayesian novices, but we hope to provide useful insights for seasoned veterans as well. Examples are drawn from gravitational-wave astronomy, though we endeavour for the presentation to be understandable to a broader audience. We begin with a review of the fundamentals: likelihoods, priors, and posteriors. Next, we discuss Bayesian evidence, Bayes factors, odds ratios, and model selection. From there, we describe how posteriors are estimated using samplers such as Markov Chain Monte Carlo algorithms and nested sampling. Finally, we generalise the formalism to discuss hyper-parameters and hierarchical models. We include extensive appendices discussing the creation of credible intervals, Gaussian noise, explicit marginalisation, posterior predictive distributions, and selection effects.

We describe an ultra-wide-bandwidth, low-frequency receiver recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band ( ${\sim}60\%$ ), the system temperature is approximately 22 K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver, including its astronomical objectives, as well as the feed, receiver, digitiser, and signal processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified, including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration, and timing stability.

The S-band Polarisation All-Sky Survey has observed the entire southern sky using the 64-m Parkes radio telescope at 2.3 GHz with an effective bandwidth of 184 MHz. The surveyed sky area covers all declinations δ ⩽ 0°. To analyse compact sources, the survey data have been re-processed to produce a set of 107 Stokes I maps with 10.75 arcmin resolution and the large scale emission contribution filtered out. In this paper, we use these Stokes I images to create a total intensity southern-sky extragalactic source catalogue at 2.3 GHz. The source catalogue contains 23 389 sources and covers a sky area of 16 600 deg2, excluding the Galactic plane for latitudes |b| < 10°. Approximately, 8% of catalogued sources are resolved. S-band Polarisation All-Sky Survey source positions are typically accurate to within 35 arcsec. At a flux density of 225 mJy, the S-band Polarisation All-Sky Survey source catalogue is more than 95% complete, and ~ 94% of S-band Polarisation All-Sky Survey sources brighter than 500 mJy beam−1 have a counterpart at lower frequencies.

Closure phase is the phase of a closed-loop product of spatial coherences formed by a ${\ge}3$ -element interferometer array. Its invariance to phase corruption attributable to individual array elements acquired during the propagation and the measurement processes, subsequent calibration, and errors therein, makes it a valuable tool in interferometry applications that otherwise require high-accuracy phase calibration. However, its understanding has remained mainly mathematical and limited to the aperture plane (Fourier dual of the image plane). Here, we present a geometrical, image domain view of closure phase, which until now has been lacking. Using the principal triangle in a 3-element interference image formed by a triad of interferometer elements, we show that the properties of closure phase, particularly its invariance to multiplicative element-based corruption factors (even of a large magnitude) and to translation, are intricately related to the conserved properties of the triangle, namely, its shape, orientation, and size, which is referred herein as the ‘shape-orientation-size (SOS) conservation principle’. In the absence of a need for element-based amplitude calibration of the interferometer array (as is typical in optical interferometry), the principal triangle in any 3-element interference image formed from phase-uncalibrated spatial coherences is still a true and uncorrupted representation of the source object’s morphology, except for a possible shift. Based on this knowledge of the triangle SOS conservation principle, we present two geometric methods to measure the closure phase directly from a simple 3-element interference image (without requiring an aperture-plane view): (i) the closure phase is directly measurable from any one of the triangle’s heights, and (ii) the squared closure phase is proportional to the product of the areas enclosed by the triad of array elements and the principal triangle in the aperture and image planes, respectively. We validate the geometric understanding of closure phase in the image plane using observations with the Karl G. Jansky Very Large Array, and the Event Horizon Telescope. These results verify the SOS conservation principle across a wide range of radio interferometric conditions. This geometric insight can be potentially valuable to other interferometric applications, such as optical interferometry. We also generalise these geometric relationships to an N-element interferometer.

A ‘pulsar timing array’ (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of ‘global’ phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 ms pulsars is being observed at three radio-frequency bands, 50 cm (~700 MHz), 20 cm (~1400 MHz), and 10 cm (~3100 MHz), with observations at intervals of two to three weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters, and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For 10 of the 20 pulsars, rms timing residuals are less than 1 μs for the best band after fitting for pulse frequency and its first time derivative. Significant ‘red’ timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array and a PTA based on the Square Kilometre Array. We also present an ‘extended PPTA’ data set that combines PPTA data with earlier Parkes timing data for these pulsars.

We present VELOCIraptor, a massively parallel galaxy/(sub)halo finder that is also capable of robustly identifying tidally disrupted objects and separate stellar halos from galaxies. The code is written in C++11, use the Message Passing Interface (MPI) and OpenMP Application Programming Interface (API) for parallelisation, and includes python tools to read/manipulate the data products produced. We demonstrate the power of the VELOCIraptor (sub)halo finder, showing how it can identify subhalos deep within the host that have negligible density contrasts to their parent halo. We find a subhalo mass-radial distance dependence: large subhalos with mass ratios of ≳10−2 are more common in the central regions than smaller subhalos, a result of dynamical friction and low tidal mass loss rates. This dependence is completely absent in (sub)halo finders in common use, which generally search for substructure in configuration space, yet is present in codes that track particles belonging to halos as they fall into other halos, such as hbt+. VELOCIraptor largely reproduces the dependence seen without tracking, finding a similar radial dependence to hbt+ in well-resolved halos from our limited resolution fiducial simulation.

The Galactic bulge of the Milky Way is made up of stars with a broad range of metallicity, –3.0 < [Fe/H] < 1 dex. The mean of the metallicity distribution function decreases as a function of height z from the plane and, more weakly, with galactic radius RGC. The most metal-rich stars in the inner Galaxy are concentrated to the plane and the more metal-poor stars are found predominantly further from the plane, with an overall vertical gradient in the mean of the metallicity distribution function of about − 0.45 dex kpc−1. This vertical gradient is believed to reflect the changing contribution with height of different populations in the innermost region of the Galaxy. The more metal-rich stars of the bulge are part of the boxy/peanut structure and comprise stars in orbits which trace out the underlying X-shape. There is still a lack of consensus on the origin of the metal-poor stars ([Fe/H] < −0.5) in the region of the bulge. Some studies attribute the more metal-poor stars of the bulge to the thick disk and stellar halo that are present in the inner region, and other studies propose that the metal-poor stars are a distinct ‘old spheroid’ bulge population. Understanding the origin of the populations that make up the metallicity distribution function of the bulge, and identifying if there is a unique bulge population which has formed separately from the disk and halo, has important consequences for identifying the relevant processes in the formation and evolution of the Milky Way.

GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination + 25° at frequencies between 72 and 231 MHz, made with the MWA using a drift scan method that makes efficient use of the MWA’s very large field-of-view. We present the observation details, imaging strategies, and theoretical sensitivity for GLEAM. The survey ran for two years, the first year using 40-kHz frequency resolution and 0.5-s time resolution; the second year using 10-kHz frequency resolution and 2 s time resolution. The resulting image resolution and sensitivity depends on observing frequency, sky pointing, and image weighting scheme. At 154 MHz, the image resolution is approximately 2.5 × 2.2/cos (δ + 26.7°) arcmin with sensitivity to structures up to ~ 10° in angular size. We provide tables to calculate the expected thermal noise for GLEAM mosaics depending on pointing and frequency and discuss limitations to achieving theoretical noise in Stokes I images. We discuss challenges, and their solutions, that arise for GLEAM including ionospheric effects on source positions and linearly polarised emission, and the instrumental polarisation effects inherent to the MWA’s primary beam.

Epoch of Reionisation (EoR) data analysis requires unprecedented levels of accuracy in radio interferometer pipelines. We have developed an imaging power spectrum analysis to meet these requirements and generate robust 21 cm EoR measurements. In this work, we build a signal path framework to mathematically describe each step in the analysis, from data reduction in the Fast Holographic Deconvolution (FHD) package to power spectrum generation in the εppsilon package. In particular, we focus on the distinguishing characteristics of FHD/εppsilon: highly accurate spectral calibration, extensive data verification products, and end-to-end error propagation. We present our key data analysis products in detail to facilitate understanding of the prominent systematics in image-based power spectrum analyses. As a verification to our analysis, we also highlight a full-pipeline analysis simulation to demonstrate signal preservation and lack of signal loss. This careful treatment ensures that the FHD/εppsilon power spectrum pipeline can reduce radio interferometric data to produce credible 21 cm EoR measurements.

In recent years, the discovery of massive quasars at $z\sim7$ has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on the discussions held at the Monash Prato Centre from November 20 to 24, 2017, during the workshop ‘Titans of the Early Universe: The Origin of the First Supermassive Black Holes’.

We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain ( ${\sim}2.8\,\mbox{K Jy}^{-1}$ ) low-system temperature ( ${\sim}18\,\mbox{K at }20\,\mbox{cm}$ ) radio array that currently operates at 580–1 670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar $\mbox{J}0737{-}3039\mbox{A}$ , pulse profiles from 34 millisecond pulsars (MSPs) from a single 2.5-h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR $\mbox{J}0540{-}6919$ , and nulling identified in the slow pulsar PSR J0633–2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright MSPs confirm that MeerKAT delivers exceptional timing. PSR $\mbox{J}2241{-}5236$ exhibits a jitter limit of $<4\,\mbox{ns h}^{-1}$ whilst timing of PSR $\mbox{J}1909{-}3744$ over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1 000 pulsars per day and the future deployment of S-band (1 750–3 500 MHz) receivers will further enhance its capabilities.