We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties and cover the physically plausible range, i.e. $f_\mathrm{esc}$ effectively decreasing, being constant, or increasing with halo mass. We find the 21cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array.

One of the most important parameters in characterizing the Epoch of Reionization, the escape fraction of ionizing photons, $f_\mathrm{esc}$, remains unconstrained both observationally and theoretically. With recent work highlighting the impact of galaxy-scale feedback on the instantaneous value of $f_\mathrm{esc}$, it is important to develop a model in which reionization is self-consistently coupled to galaxy evolution. In this work, we present such a model and explore how physically motivated functional forms of $f_\mathrm{esc}$ affect the evolution of ionized hydrogen within the intergalactic medium. Using the $21$cm power spectrum evolution, we investigate the likelihood of observationally distinguishing between a constant $f_\mathrm{esc}$ and other models that depend upon different forms of galaxy feedback. We find that changing the underlying connection between $f_\mathrm{esc}$ and galaxy feedback drastically alters the large-scale $21$cm power. The upcoming Square Kilometre Array Low Frequency instrument possesses the sensitivity to differentiate between our models at a fixed optical depth, requiring only $200$ hours of integration time focused on redshifts $z = 7.5-8.5$. Generalizing these results to account for a varying optical depth will require multiple $800$ hour observations spanning redshifts $z = 7-10$. This presents an exciting opportunity to observationally constrain one of the most elusive parameters during the Epoch of Reionization.

The exact role of quasars during the Epoch of Reionization remains uncertain. With consensus leaning towards quasars producing a negligible amount of ionizing photons, we pose an alternate question: Can quasars indirectly contribute to reionization by allowing ionizing photons from stars to escape more easily? Using the Semi-Analytic Galaxy Evolution model to evolve a galaxy population through cosmic time, we construct an idealized scenario in which the escape fraction of stellar ionizing photons ($f_\mathrm{esc}$) is boosted following quasar wind events, potentially for several dynamical times. We find that under this scenario, the mean value of $f_\mathrm{esc}$ as a function of galaxy stellar mass peaks for intermediate mass galaxies. This mass dependence will have consequences for the 21cm power spectrum, enhancing power at small scales and suppressing it at large scales. This hints that whilst quasars may not directly contribute to the ionizing photon budget, they could influence reionization indirectly by altering the topology of ionized regions.

As an alternative explanation for the cosmic acceleration, $f(R)$ theories of gravity can predict an almost identical expansion history to standard $\Lambda$CDM, yet make very different predictions for the growth of cosmological structures. Measurements of the cosmic bulk flow provides a method for determining the strength of gravity over the history of structure formation. We use the modified gravity N-body code ECOSMOG to simulate dark matter particles and make predictions for the bulk flow magnitude in both $\Lambda$CDM and $f(R)$ gravity. With the peculiar velocities output by ECOSMOG we determine the bulk flow at depths ranging from $20h^{-1}$Mpc to $50h^{-1}$Mpc, following the redshift and sky distribution of the 2MASS Tully-Fisher survey (2MTF). At each depth, we find that the $\Lambda$CDM and $f_{R0} = 10^{-5}$ simulations produce bulk flow measurements that are consistent with $\Lambda$CDM predictions and the 2MTF survey at a $1\sigma$ level. We also find that adopting an $f(R)$ strength of $f_{R0} = 10^{-3}$ predict a much larger value for the bulk flow, which disagree with $\Lambda$CDM predictions at all depths considered. We conclude that $f_{R0}$ must be constrained to a level no greater than $10^{-4}$ to agree with bulk flow measurements.

We present the first analytical inspiral-merger-ringdown gravitational waveforms from binary black holes (BBHs) with non-precessing spins, that is based on a description of the late-inspiral, merger and ringdown in full general relativity. By matching a post-Newtonian description of the inspiral to a set of numerical-relativity simulations, we obtain a waveform family with a conveniently small number of physical parameters. These waveforms will allow us to detect a larger parameter space of BBH coalescence, including a considerable fraction of precessing binaries in the comparable-mass regime, thus significantly improving the expected detection rates.

Binary black-hole systems with spins aligned or anti-aligned to the orbital angular momentum provide the natural ground to start detailed studies of the influence of strong-field spin effects on gravitational wave observations of coalescing binaries. Furthermore, such systems may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins parallel to the orbital angular momentum but otherwise arbitrary. Our attention is particularly focused on the gravitational-wave emission so as to quantify how much spin effects contribute to the signal-to-noise ratio, to the horizon distances, and to the relative event rates for the representative ranges in masses and detectors. As expected, the signal-to-noise ratio increases with the projection of the total black hole spin in the direction of the orbital momentum. We find that equal-spin binaries with maximum spin aligned with the orbital angular momentum are more than "three times as loud" as the corresponding binaries with anti-aligned spins, thus corresponding to event rates up to 30 times larger. We also consider the waveform mismatch between the different spinning configurations and find that, within our numerical accuracy, binaries with opposite spins S_1=-S_2 cannot be distinguished whereas binaries with spin S_1=S_2 have clearly distinct gravitational-wave emissions. Finally, we derive a simple expression for the energy radiated in gravitational waves and find that the binaries always have efficiencies E_rad/M > 3.6%, which can become as large as E_rad/M = 10% for maximally spinning binaries with spins aligned with the orbital angular momentum.

We provide a compact analytic formula to compute the spin of the black hole produced by the coalescence of two black holes. The expression, which uses an analytic fit of numerical-relativity data and relies on four assumptions, aims at modelling generic initial spin configurations and mass ratios. A comparison with numerical-relativity simulations already shows very accurate agreements with all of the numerical data available to date, but we also suggest a number of ways in which our predictions can be further improved.

Determining the final spin of a black-hole (BH) binary is a question of key importance in astrophysics. Modelling this quantity in general is made difficult by the fact that it depends on the 7-dimensional space of parameters characterizing the two initial black holes. However, in special cases, when symmetries can be exploited, the description can become simpler. For black-hole binaries with unequal masses but with equal spins which are aligned with the orbital angular momentum, we show that the use of recent simulations and basic but exact constraints derived from the extreme mass-ratio limit allow to model this quantity with a simple analytic expression. Despite the simple dependence, the expression models very accurately all of the available estimates, with errors of a couple of percent at most. We also discuss how to use the fit to predict when a Schwarzschild BH is produced by the merger of two spinning BHs, when the total angular momentum of the spacetime ``flips'' sign, or under what conditions the final BH is ``spun-up'' by the merger. Finally, suggest an extension of the fit to include unequal-spin binaries, thus potentially providing a complete description of the final spin from the coalescence of generic black-hole binaries with spins aligned to the orbital angular momentum.