We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We study the implications of deformed quantum algebras for the generation of primordial perturbations from slow-roll inflation. Specifically, we assume that the quantum commutator of the inflaton's amplitude and momentum in Fourier space gets modified at energies above some threshold $M_{\star}$. We show that when the commutator is modified to be a function of the momentum only, the problem of solving for the post-inflationary spectrum of fluctuations is formally equivalent to solving a one-dimensional Schr\"odinger equation with a time dependent potential. Depending on the class of modification, we find results either close to or significantly different from nearly scale invariant spectra. For the former case, the power spectrum is characterized by step-like behaviour at some pivot scale, where the magnitude of the jump is $\mathcal{O}(H^{2}/M_{\star}^{2})$. ($H$ is the inflationary Hubble parameter.) We use our calculated power spectra to generate predictions for the cosmic microwave background and baryon acoustic oscillations, hence demonstrating that certain types of deformations are incompatible with current observations.

We use cosmological perturbation theory to study the backreaction effects of a self-consistent and well-defined cosmological averaging on the dynamics and the evolution of the Universe. Working with a perturbed Friedman-Lemaitre-Robertson-Walker Einstein-de Sitter cosmological solution in a comoving volume-preserving gauge, we compute the expressions for the expansion scalar and deceleration parameter to second order, which we use to characterize the backreaction. We find that the fractional shift in the Hubble parameter with respect to the input background cosmological model is Delta~10^{-5}, which leads to an effective energy density of the order of a few times 10^{-5}. In addition, we find that an appropriate measure of the fractional shift in the deceleration parameter Q is very large.

The cosmological backreaction from perturbations is clearly gauge-dependent, and obviously depends on the choice of averaged Hubble rate. We consider two common choices of Hubble rate and advocate the use of comoving volume-preserving gauges. We highlight two examples valid to an appropriate order in perturbation theory, uniform curvature gauge, which is as close to volume-preserving as possible, and a spatially-traceless uniform cold dark matter gauge which preserves the volume to linear order. We demonstrate the strong gauge- and frame-dependences in averaging. In traceless uniform CDM gauge the backreaction exhibits a strong ultra-violet divergence and can be tuned to an arbitrary magnitude with an appropriate choice of smoothing scale. In uniform curvature gauge we find that for a choice of Hubble rate locked to the spatial surface the backreaction vanishes identically, while for a Hubble rate defined from a fluid's expansion scalar the effective energy density at the current epoch in an Einstein-de Sitter universe is Omega_eff~5e-4, slightly bigger than but in broad agreement with previous results in conformal Newtonian gauge.

Non-Gaussianity in the cosmic microwave background and the large-scale structure of galaxies provides an increasingly powerful probe of the universe. I implement an algorithm to generate realisations of fields that possess an arbitrary probability distribution function and an arbitrary power spectrum and demonstrate the code with a number of examples, including the uniform distribution, the Laplace distribution, the $\chi$ and $\chi^2$ distributions, Rayleigh and Maxwell-Boltzmann distributions. The code is available at http://sourceforge.net/projects/nongaussian.

We study the generation of primordial fluctuations in pure de Sitter inflation where the quantum scalar field dynamics are governed by polymer (not Schrodinger) quantization. This quantization scheme is related to, but distinct from, the structures employed in Loop Quantum Gravity; and it modifies standard results above a polymer energy scale $M_{\star}$. We recover the scale invariant Harrison Zel'dovich spectrum for modes that have wavelengths bigger than $M_{\star}^{-1}$ at the start of inflation. The primordial spectrum for modes with initial wavelengths smaller than $M_{\star}^{-1}$ exhibits oscillations superimposed on the standard result. The amplitude of these oscillations is proportional to the ratio of the inflationary Hubble parameter $H$ to the polymer energy scale. For reasonable choices of $M_{\star}$, we find that polymer effects are likely unobservable in CMB angular power spectra due to cosmic variance uncertainty, but future probes of baryon acoustic oscillations may be able to directly constrain the ratio $H/M_{\star}$.

Understanding the non-adiabatic pressure, or relative entropy, perturbation is crucial for studies of early-universe vorticity and Cosmic Microwave Background observations. We calculate the evolution of the linear non-adiabatic pressure perturbation from radiation domination to late times, numerically solving the linear governing equations for a wide range of wavenumbers. Using adiabatic initial conditions consistent with WMAP seven year data, we find nevertheless that the non-adiabatic pressure perturbation is non-zero and grows at early times, peaking around the epoch of matter/radiation equality and decaying in matter domination. At early times or large redshifts (z=10,000) its power spectrum peaks at a comoving wavenumber k~0.2h/Mpc, while at late times (z=500) it peaks at k~0.02 h/Mpc.

We consider issues related to the conformal mapping between the Einstein and Jordan frames in f(R) cosmology. We consider the impact of the conformal transformation on the gauge of a perturbed system and show that unless the system is written in a restricted set of gauges the mapping could produce an inconsistent result in the target frame. Newtonian gauge lies within the restricted group but synchronous gauge does not. If this is not treated carefully it could in principle contaminate numerical calculations.

Forthcoming datasets from the Planck experiment and others are in a position to probe the CMB non-Gaussianity with higher accuracy than has yet been possible and potentially open a new window into the physics of the very early universe. However, a signal need not necessarily be inflationary in origin, and possible contaminants should be examined in detail. One such is provided by early universe magnetic fields, which can be produced by a variety of models including during an inflationary phase, at phase transitions, or seeded by cosmic defects. Should such fields have been extent in the early universe they provide a natural source of CMB non-Gaussianity. Knowledge of the CMB angular bispectrum requires the complete Fourier-space (or "intrinsic") bispectrum. In this paper I consider in detail the intrinsic bispectra of an early-universe magnetic field for a range of power-law magnetic power spectra.

In this paper we investigate the impact of the coupling between a quintessence field and clustered matter on the average equation of state of the scalar field. We take the NFW profile to be characteristic of bound structures on galactic and cluster scales, and the isothermal distribution to hold for objects on supercluster scales. Solving analytically for the scalar-field profile, we find that the greatest impact on the equation of state comes from the superclusters. Employing numerical case studies, we verify this effect and probe its dependence on the evolutionary state of the supercluster. An estimate across the Hubble volume yields corrections to the homogeneous equation of state of ~3%, increasing with coupling strength.

Magnetic fields appear to be a generic feature of the early universe and are a natural source of secondary CMB non-Gaussianity. In recent years the statistical nature of the stresses of a primordial magnetic field has been well studied. In this paper we confirm and extend these studies at one- and two-point level, and present analytical results for a wide range of power-law spectra. We also consider two non-power law cases of interest: a blue spectrum with an extended damping tail on small scales, which could be generated by the non-linear mixing of density and vorticity; and a red spectrum with a damping tail on large scales. We then briefly consider the CMB impacts that result from such fields. While this paper focuses on the one- and two-point moments, the techniques we employ are designed to ease the analysis of the full bispectra induced by primordial magnetic fields.

Inflation generically produces primordial gravitational waves with a red spectral tilt. In this paper we calculate the backreaction produced by these gravitational waves on the expansion of the universe. We find that in radiation domination the backreaction acts as a relativistic fluid, while in matter domination a small dark energy emerges with an equation of state w=-8/9.

We present a formalism for spatial averaging in cosmology applicable to general spacetimes and coordinates, and allowing the easy incorporation of a wide variety of matter sources. We apply this formalism to a Friedmann-LeMaitre-Robertson-Walker universe perturbed to second-order and present the corrections to the background in an unfixed gauge. We then present the corrections that arise in uniform curvature and conformal Newtonian gauges.

The cosmological backreaction arises when one directly averages the Einstein equations to recover an effective Robertson-Walker cosmology, rather than assuming a background a priori. While usually discussed in the context of dark energy, strictly speaking any cosmological model should be recovered from such a procedure. We apply the Buchert averaging formalism to linear Robertson-Walker universes containing matter, radiation and dark energy and evaluate numerically the discrepancies between the assumed and the averaged behaviour, finding the largest deviations for an Einstein-de Sitter universe, increasing rapidly with Hubble rate to a 0.01% effect for h=0.701. For the LCDM concordance model, the backreaction is of the order of Omega_eff~4x10^-6, with those for dark energy models being within a factor of two or three. The impacts at recombination are of the order of 10^-8 and those in deep radiation domination asymptote to a constant value. While the effective equations of state of the backreactions in Einstein-de Sitter, concordance and quintessence models are generally dust-like, a backreaction with an equation of state w_eff<-1/3 can be found for strongly phantom models.

Magnetic fields have been observed in galaxies, clusters of galaxies and probably in superclusters. While mechanisms exist to generate these in the late universe, it is possible that magnetic fields have existed since very early times. This thesis is concerned with methods to predict the form of such imprints. We review in detail a standard, linearised cosmology before introducing an electromagnetic field. We then consider the intrinsic statistics of the magnetic stresses in two ways, analytically and via static realisations. We construct the power spectra, some of which we present for the first time. At the one- and three-point level we find significant intrinsic non-Gaussianities. Finally we turn to the observable impacts a primordial magnetic field. Assuming coherence, the statistics of the source can be mapped onto the CMB in a simple manner. We demonstrate that our approach is valid by reproducing the signals for Gaussian power law fields on the microwave sky. [ABRIDGED]

We reformulate the averaged Einstein equations in a form suitable for use with Newtonian gauge linear perturbation theory and track the size of the modifications to standard Robertson-Walker evolution on the largest scales as a function of redshift for both Einstein de-Sitter and Lambda CDM cosmologies. In both cases the effective energy density arising from linear perturbations is of the order of 10^-5 the matter density, as would be expected, with an effective equation of state w ~ -1/19. Employing a modified Halofit code to extend our results to quasilinear scales, we find that, while larger, the deviations from Robertson-Walker behaviour remain of the order of 10^-5.

We report the first observation of a transient relativistic jet from the canonical black hole candidate, Cygnus X-1, obtained with the Multi-Element Radio-Linked Interferometer Network (MERLIN). The jet was observed in only one of six epochs of MERLIN imaging of the source during a phase of repeated X-ray spectral transitions in 2004 Jan--Feb, and this epoch corresponded to the softest 1.5-12 keV X-ray spectrum. With only a single epoch revealing the jet, we cannot formally constrain its velocity. Nevertheless, several lines of reasoning suggest that the jet was probably launched 0.5-4.0 days before this brightening, corresponding to projected velocities of 0.2c < v_app < 1.6c, and an intrinsic velocity of > 0.3c. We also report the occurrence of a major radio flare from Cyg X-1, reaching a flux density of ~120 mJy at 15 GHz, and yet not associated with any resolvable radio emission, despite a concerted effort with MERLIN. We discuss the resolved jet in terms of the recently proposed 'unified model' for the disc-jet coupling in black hole X-ray binaries, and tentatively identify the 'jet line' for Cyg X-1. The source is consistent with the model in the sense that a steady jet appears to persist initially when the X-ray spectrum starts softening, and that once the spectral softening is complete the core radio emission is suppressed and transient ejecta / shock observed. However, there are some anomalies, and Cyg X-1 clearly does not behave like a normal black hole transient in progressing to the canonical soft / thermal state once the ejection event has happened.

Magnetic fields in the early universe could have played an important role in sourcing cosmological perturbations. While not the dominant source, even a small contribution might be traceable through its intrinsic non-Gaussianity. Here we calculate analytically the one, two and three point statistics of the magnetic stress energy resulting from tangled Gaussian fields, and confirm these with numerical realizations of the fields. We find significant non-Gaussianity, and importantly predict higher order moments that will appear between the scalar, vector and tensor parts of the stress energy (e.g. scalar-tensor-tensor moments). Such higher order cross correlations are a generic feature of non-linear theories and could prove to be an important probe of the early universe.

A series of observations of microquasar GRS1915+105 were made by the MERLIN telescope array in the spring of 2003 at 18 cm. The images show polarization in the jet components, with two epochs displaying Faraday rotation of about 2500 rad/m2. This effect is a demonstration of the similarities between microquasars and extragalactic jets.

We present the international collaboration MINE (Multi-lambda Integral NEtwork) aimed at conducting multi-wavelength observations of X-ray binaries and microquasars simultaneously with the INTEGRAL gamma-ray satellite. We will focus on the 2003 March-April campaign of observations of the peculiar microquasar GRS 1915+105 gathering radio, IR and X-ray data. The source was observed 3 times in the plateau state, before and after a major radio and X-ray flare. It showed strong steady optically thick radio emission corresponding to powerful compact jets resolved in the radio images, bright near-infrared emission, a strong QPO at 2.5 Hz in the X-rays and a power law dominated spectrum without cutoff in the 3-300 keV range. We compare the different observations, their multi-wavelength light curves, including JEM-X, ISGRI and SPI, and the parameters deduced from fitting the spectra obtained with these instruments on board INTEGRAL.

We present the international collaboration MINE (Multi-lambda INTEGRAL NEtwork) aimed at conducting multi-wavelength observations of microquasars simultaneously with the INTEGRAL satellite. The first results on GRS 1915+105 are encouraging and those to come should help us to understand the physics of the accretion and ejection phenomena around a compact object.