Motivated by the prospect of the wealth of data arising from the inauguration of the era of gravitational wave detection by ground-based interferometers the DES collaboration, in partnership with members of the LIGO collaboration and members of the astronomical community at large, have established a research program to search for their optical counterparts and to explore their use as cosmological probes. In this talk we present the status of our program and discuss prospects for establishing this new probe as part of the portfolio of the Dark Energy research program in the future, in particular for the next generation survey, LSST.

The evolution of neutral hydrogen (HI) across redshifts is a powerful probe of cosmology, large scale structure in the universe and the intergalactic medium. Using a data-driven halo model to describe the distribution of HI in the post-reionization universe (z ∼ 5 to 0), we obtain the best-fitting parameters from a rich sample of observational data: low redshift 21-cm emission line studies, intermediate redshift intensity mapping experiments, and higher redshift Damped Lyman Alpha (DLA) observations. Our model describes the abundance and clustering of neutral hydrogen across redshifts 0 - 5, and is useful for investigating different aspects of galaxy evolution and for comparison with hydrodynamical simulations. The framework can be applied for forecasting future observations with neutral hydrogen, and extended to the case of intensity mapping with molecular and other line transitions at intermediate redshifts.

Observations of HI 21cm transition line is a promising probe into the Dark Ages and Epoch-of-Reionization. Detection of this redshifted 21cm signal is one of the key science goal for several upcoming low-frequency radio telescopes like HERA, SKA and DARE. Other global signal experiments include EDGES, LEDA, BIGHORNS, SCI-HI, SARAS. One of the major challenges for the detection of this signal is the accuracy of the foreground source removal. Several novel techniques have been explored already to remove bright foregrounds from both interferometric as well as total power experiments. Here, we present preliminary results from our investigation on application of ANN to detect 21cm global signal amidst bright galactic foreground. Following the formalism of representing the global 21cm signal by ’tanh’ model, this study finds that the global 21cm signal parameters can be accurately determined even in the presence of bright foregrounds represented by 3rd order log-polynomial or higher.

Compact object mergers are promising candidates for the progenitor system of short gamma-ray bursts (GRBs). Using gravitational wave (GW) triggers to identify a merger, any electromagnetic (EM) counterparts from the jet can be used to constrain the dynamics and structure of short GRB jets. GW triggered searches could reveal a hidden population of optical transients associated with the short-lived jets from the merger object. If the population of merger-jets is dominated by low-Lorentz-factors, then a GW triggered search will reveal the on-axis orphan afterglows from these failed GRBs. By considering the EM counterparts from a jet, with or without the prompt GRB, the jet structure and dynamics can be constrained. By modelling the afterglow of various jet structures with viewing angle, we provide observable predictions for the on- and off- axis EM jet counterparts. The predictions provide an indication for the various features expected from the proposed jet structure models.

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.