We extend our MCMC sampler of 3D EoR simulations, 21CMMC, to perform parameter estimation directly on light-cones of the cosmic 21cm signal. This brings theoretical analysis one step closer to matching the expected 21-cm signal from next generation interferometers like HERA and the SKA. Using the light-cone version of 21CMMC, we quantify biases in the recovered astrophysical parameters obtained from the 21cm power spectrum when using the co-eval approximation to fit a mock 3D light-cone observation. While ignoring the light-cone effect does not bias the parameters under most assumptions, it can still underestimate their uncertainties. However, significant biases (∼few – 10 σ) are possible if all of the following conditions are met: (i) foreground removal is very efficient, allowing large physical scales (k ∼ 0.1 Mpc−1) to be used in the analysis; (ii) theoretical modelling is accurate to ∼10 per cent in the power spectrum amplitude; and (iii) the 21cm signal evolves rapidly (i.e. the epochs of reionisation and heating overlap significantly

This review presents briefly two aspects of Galactic interstellar medium science that seem relevant for studying EoR. First, we give some statistical properties of the Galactic foreground emission in the diffuse regions of the sky. The properties of the emission observed in projection on the plane of the sky are then related to how matter is organised along the line of sight. The diffuse atomic gas is multi-phase, with dense filamentary structures occupying only about 1% of the volume but contributing to about 50% of the emission. The second part of the review presents aspect of structure formation in the Galactic interstellar medium that could be relevant for the subgrid physics used to model the formation of the first stars.

We study the inhomogeneity of the reionization process by comparing the reionization times of z = 0 galaxies as a function of their mass. For this purpose, we combine the results of the CODA-I AMR radiative hydrodynamics simulation of the Reionization with the halo merger trees of a pure dark matter tree-code z = 0 simulation evolved from the same set of initial conditions. We find that galaxies with M(z = 0) > 1011M⊙ are reionized earlier than the whole Universe, with e.g. MW-like haloes reionized between 100 and 300 million years before the diffuse IGM. Lighter galaxies reionized as late as the global volume, probably from external radiation.

The latest generation of low frequency radio interferometers, e.g. LOFAR, MWA, PAPER, has been pushing down the detection limits on the hydrogen signal from the Epoch of Reionisation. However, due to the challenges posed by foregrounds and instrumental systematics the signal has eluded detection thus far. To overcome these challenges we require a detailed understanding of the calibration of these relatively new telescopes. This led to a renewed interest in redundant calibration. Classical calibration schemes depend on sky models based on limited knowledge of the low frequency sky. Redundant calibration, however, allows us to escape our ignorance as it is sky model independent. We will review the field of redundant calibration, and present work we have undertaken to understand the limitations of this calibration method.

The infrared-radio correlation (IRRC) offers a way to assess star formation from radio emission. Multiple studies found the IRRC to decrease with increasing redshift. This may in part be due to the lack of knowledge about the possible radio spectral energy distributions (SEDs) of star-forming galaxies. We constrain the radio SED of a complete sample of highly star-forming galaxies (SFR > 100 M⊙/ yr) based on the VLA-COSMOS 1.4 GHz Joint and 3 GHz Large Project catalogs. We reduce archival GMRT 325 MHz and 610 MHz observations, broadening the rest-frame frequency range to 0.3-15 GHz. Employing survival analysis and fitting a double power law SED, we find that the slope steepens from a spectral index of α1 = 0.51±0.04 below 4.5 GHz to α2 = 0.98±0.07 above 4.5 GHz. Our results suggest that the use of a K-correction assuming a single power-law radio SED for star forming galaxies is likely not the root cause of the IRRC trend.

The paper presents a proof of concept method of background rejection based on convolutional neural networks (CNN). The method was tested on simulated data and achieved very high accuracy (100%). What is more, method based on CNN is very fast and could be easily applied to wide field surveys. Since early stage results suggest method is very accurate and robust, it could be helpful in creating very low-latency pipelines for EM Follow-up purposes, which will be needed in LIGO-Virgo O3 EM Follow-up.

Characterizing the diffuse Galactic synchrotron emission (DGSE) at arcminute angular scales is needed to remove this foregrounds in cosmological 21-cm measurements. Here, we present the angular power spectrum (Cℓ) measurement of the diffuse Galactic synchrotron emission using two fields observed by the TIFR GMRT Sky Survey (TGSS). We apply 2D Tapered Gridded Estimator (TGE) to estimate the Cℓ from the visibilities. We find that the residual data after subtracting the point sources is likely dominated by the diffuse Galactic synchrotron radiation across the angular multipole range 240 ≤ ℓ ≲ 500. We fit a power law to the measured Cℓ over this ℓ range. We find that the slopes in both fields are consistent with earlier measurements. For the second field, however, we interpret the measured Cℓ as an upper limit for the DGSE as there is an indication of a significant residual point source contribution.

The EOR1 field at (RA,DEC)=(4hrs, −30.0°) is one of the main targets of the MWA EOR experiment. It is notable for being in one of the coldest regions of the southern radio sky, as well as for containing the bright radio galaxy Fornax A. We report an early demonstration that the distance of this field from the Galactic Centre may make it a prime field for EOR observations.

Blazar OJ287 exhibits large thermal flares at least twice every 12 years. The times of these flares have been predicted successfully using the model of a quasi-Keplerian eccentric black hole binary where the secondary impacts the accretion disk of the primary, creating the thermal flares. New measurements of the historical light curve have been combined with the observations of the 2015 November/December flare to identify the impact record since year 1886, and to constrain the orbit of the binary. The orbital solution shows that the binary period, now 12.062 years, is decreasing at the rate of 36 days per century. This corresponds to an energy loss to gravitational waves that is 6.5 ± 4 % less than the rate predicted by the standard quadrupolar gravitational wave (GW) emission. We show that the difference is due to higher order gravitational radiation reaction terms that include the dominant order tail contributions.

Direct detection of the Epoch of Reionization via the redshifted 21-cm line will have unprecedented implications on the study of structure formation in the early Universe. To fulfill this promise current and future 21-cm experiments will need to detect the weak 21-cm signal over foregrounds several order of magnitude greater. This requires accurate modeling of the galactic and extragalactic emission and of its contaminants due to instrument chromaticity, ionosphere and imperfect calibration. To solve for this complex modeling, we propose a new method based on Gaussian Process Regression (GPR) which is able to cleanly separate the cosmological signal from most of the foregrounds contaminants. We also propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere. Using this method, chromatic effects causing the so-called “wedge” are effectively eliminated (i.e. deconvolved) in the cylindrical (k⊥, k∥) power spectrum.

Cosmic reionization is thought to be initiated by the first generation of stars and blackholes. We review recent progress in theoretical studies of early structure formation. Cosmic structure formation is driven by gravitational instability of primeval density fluctuations left over from Big Bang. At early epochs, there are baryonic streaming motions with significant relative velocity with respect to dark matter. The formation of primordial gas clouds is typically delayed by the streaming motions, but then physical conditions for the so-called direct collapse blackhole formation are realized in proto-galactic halos. We present a promising model in which intermediate mass blackholes are formed as early as z = 30.

Here I discuss progress in both the theory and practice of data analysis for the Hydrogen Epoch of Reionization Array (HERA), focusing on techniques to calibrate the instrumental response and preserve the spectral smoothness that is essential to separating the cosmological 21 cm signal from foregrounds that are five orders of magnitude brighter. I explain how mis-calibration can create ruinous spectral structure and how we take advantage of HERA’s highly-redundant configuration for calibration. This proceeding draws from a talk I gave on October 3, 2017. Slides for it and all my talks are available at joshdillon.net.

We conducted a systematic search for periodically varying quasars, which are predicted manifestations of sub-pc supermassive black hole binaries (SMBHBs), in the Pan-STARRS1 Medium Deep Survey (PS1 MDS). Since the normal variability of quasars can also mimic periodicity over a small number of cycles, we have extended the temporal baseline by monitoring the candidates with the Discovery Channel Telescope and the Las Cumbres Observatory telescopes. We have also adopted a more rigorous method to evaluate the significance of the periodic candidates, by considering in the light curves a “red noise” background modeled as the Damped Random Walk process. While none of the candidates can be resolved by the current pulsar timing arrays (PTAs) as individual gravitational wave sources, the Large Synoptic Survey Telescope is capable of finding more periodic candidates, some of which are likely to be detected by the PTA experiment with the Square Kilometre Array.

The matched filtering technique is an efficient method to detect H ii bubbles and absorption regions in radio interferometric observations of the redshifted 21-cm signal from the epoch of reionization and the Cosmic Dawn. Here, we present an implementation of this technique to the upcoming observations such as the SKA1-low for a blind search of absorption regions at the Cosmic Dawn. The pipeline explores four dimensional parameter space on the simulated mock visibilities using a MCMC algorithm. The framework is able to efficiently determine the positions and sizes of the absorption/H ii regions in the field of view.

The identification of the electromagnetic(EM) counterpart of gravitational wave(GW) trigger in the sky localization is a very difficult task because of the large uncertainty. Two complementary approaches are used in order to search for EM counterpart of GW signal with a typical large sky localization uncertainty: wide-field tilling search on high probability GW region, e.g. Gravitational Wave Inaf Team(GRAWITA) project or pointed search of selected galaxies in high probability GW region, e.g. Distance Less Than 40 Mpc survey(DLT40) project.

The DESGW program is a collaboration between members of the Dark Energy Survey, the wider astronomical community, and the LIGO-Virgo Collaboration to search for optical counterparts of gravitational wave events, such as those expected from binary neutron star mergers or neutron star-black hole mergers. While binary black hole (BBH) events are not expected to produce an electromagnetic (EM) signature, emission is certainly not impossible. The DESGW program has performed follow-up observations of four BBH events detected by LIGO in order to search for any possible EM counterpart. Failure to find such counterparts is still relevant in that it produces limits on optical emission from such events. This is a review of follow-up results from O1 BBH events and a discussion of the status of ongoing uniform re-analysis of all BBH events that DESGW has followed up to date.

The MWA EoR is one of a small handful of experiments designed to detect the statistical signal from the Epoch of Reionisation. Each of these experiments has reached a level of maturity, where the challenges, in particular of foreground removal, are being more fully understood. Over the past decade, the MWA EoR Collaboration has developed expertise and an understanding of the elements of the telescope array, the end-to-end pipelines, ionospheric conditions, and and the foreground emissions. Sufficient data has been collected to detect the theoretically predicted EoR signal. Limits have been published regularly, however we still several orders of magnitude from a possible detection. This paper outlines recent progress and indicates directions for future efforts.