Accretion of stars on massive black holes (MBHs) can feed MBHs and generate tidal disruption events (TDEs). We introduce a new physically motivated model to self-consistently treat TDEs in cosmological simulations, and apply it to the assembly of a galaxy with final mass $3\times 10^{10}\mathrm{\, M}_\odot$ at $z=6$. This galaxy exhibits a TDE rate of $\sim 10^{-5}\mathrm{\, yr}^{-1}$, consistent with local observations but already in place when the Universe was one billion year old. A fraction of the disrupted stars participate in the growth of MBHs, dominating it until the MBH reaches mass $\sim 5 \times 10^5 \mathrm{\, M}_\odot$, but their contribution then becomes negligible compared to gas. TDEs could be a viable mechanism to grow light MBH seeds, but fewer TDEs are expected when the MBH becomes sufficiently massive to reach the luminosity of, and be detected as, an active galactic nucleus. Galaxy mergers bring multiple MBHs in the galaxy, resulting in an enhancement of the global TDE rate in the galaxy by $\sim 1$ order of magnitude during $100\mathrm{\, Myr}$ around mergers. This enhancement is not on the central MBH, but caused by the presence of MBHs in the infalling galaxies. This is the first self consistent study of TDEs in a cosmological environment and highlights that accretion of stars and TDEs are a natural process occurring in a Milky~Way-mass galaxy at early cosmic times.

Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large volume and low resolution which tracks the high-mass (> 1e7 Msun) MBH population, and NewHorizon with a smaller volume but higher resolution that traces the low-mass (< 1e7 Msun) MBH population. While Horizon-AGN can be used to estimate the rate of inspirals for Pulsar Timing Arrays, NewHorizon can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. After the numerical MBH merger at the resolution limit, we post-process MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after the galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the galaxy merger from which the MBH merger originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z~1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies.

We compute the tidal disruption event (TDE) rate around local massive black holes (MBHs) with masses as low as $2.5\times10^4\mathrm{ M}_\odot$, thus probing the dwarf regime for the first time. We select a sample of 37 galaxies for which we have the surface stellar density profile, a dynamical estimate of the mass of the MBH, and 6 of which, including our Milky Way, have a resolved nuclear star cluster (NSC). For the Milky Way, we find a total TDE rate of $\sim 10^{-4}\mathrm{ yr}^{-1}$ when taking the NSC in account, and $\sim 10^{-7} \mathrm{ yr}^{-1}$ otherwise. TDEs are mainly sourced from the NSC for light ($<3\times 10^{10}\mathrm{ M}_\odot$) galaxies, with a rate of few $10^{-5}\mathrm{ yr}^{-1}$, and an enhancement of up to 2 orders of magnitude compared to non-nucleated galaxies. We create a mock population of galaxies using different sets of scaling relations to explore trends with galaxy mass, taking into account the nucleated fraction of galaxies. Overall, we find a rate of few $10^{-5}\mathrm{ yr}^{-1}$ which drops when galaxies are more massive than $10^{11}\mathrm{ M}_\odot$ and contain MBHs swallowing stars whole and resulting in no observed TDE.

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$.

Mergers of galaxies are extremely violent events shaping their evolution. Such events are thought to trigger starbursts and, possibly, black hole accretion. Nonetheless, it is still not clear how to know the fate of a galaxy pair from the data available at a given time, limiting our ability to constrain the exact role of mergers. In this paper we use the lightcone of the Horizon-AGN simulation, for which we know the fate of each pair, to test three selection processes aiming at identifying true merging pairs. We find that the simplest one (selecting objects within two thresholds on projected distance $d$ and redshift difference $\Delta z$) gives similar results than the most complex one (based on a neural network analyzing $d$, $\Delta z$, redshift of the primary, masses/star formation rates/aspect ratio of both galaxies). Our best thresholds are $d_\mathrm{th}\sim100\mathrm{\, kpc}$ and $\Delta z_\mathrm{th} \sim 10^{-3}$, in agreement with recent results.

The dynamics of black hole seeds in high redshift galaxies is key to understand their ability to grow via accretion and to pair in close binaries during galactic mergers. To properly follow the dynamics of black holes we develop a physically motivated model to capture unresolved dynamical friction from stars, dark matter and gas. We first validate the model and then we use it to investigate the dynamics of seed black holes born at $z\sim9$ in dwarf proto-galaxies. We perform a suite of zoom cosmological simulations with spatial resolution as high as 10 pc and with a stellar and dark matter mass resolution of $2\times10^3 \, \textrm{M}_\odot$ and $2\times10^5 \, \textrm{M}_\odot$ respectively. We first explore the dynamics of a seed black hole in the galaxy where it is born and show that it is highly erratic if the seed mass is less than $10^5\, \textrm{M}_\odot$. The dynamics is dominated by the stellar component, whose distribution is irregular and patchy, thus inducing stochasticity in the orbits: the black hole may be anywhere in the proto-galaxy. When this dwarf merges into a larger galaxy, it is paramount to simulate the process with very high spatial and mass resolution in order to correctly account for the stripping of the stellar envelope of the satellite black hole. The outcome of the encounter could be either a tight binary or, at least temporary, a wandering black hole, leading to multiple black holes in a galaxy, each inherited from a different merger.

The rate of tidal disruption events (TDEs) depends sensitively on the stellar properties of the central galactic regions. Simulations show that galaxy mergers cause gas inflows that trigger nuclear starbursts that increase the central stellar density. Motivated by these numerical results, and by the observed over-representation of post-starburst galaxies among TDE hosts, we study the evolution of the TDE rate in high-resolution hydrodynamical simulations of a galaxy merger, in which we capture the evolution of the stellar density around the massive black holes (BHs) in the merging galaxies and in the final remnant. We post-process the simulation to obtain the TDE rate. At the second pericentre, a nuclear starburst enhances the stellar density around the BH in the least massive galaxy, leading to an enhancement of the TDE rate around the secondary BH, although the magnitude and the duration of the increase depend on the stochasticity of star formation on very small scales. The central stellar density around the primary BH remains instead fairly constant, and so is its TDE rate, until the two BHs bind and are surrounded by the same stellar environment. After the formation of the binary, the stellar density decreases, and so does the TDE rate.

We study the dynamical evolution of supermassive black holes, in the late stage of galaxy mergers, from kpc to pc scales. In particular, we capture the formation of the binary, a necessary step before the final coalescence, and trace back the main processes causing the decay of the orbit. We use hydrodynamical simulations of galaxy mergers with different resolutions, from $20\,\rm pc$ down to $1\,\rm pc$, in order to study the effects of the resolution on our results, remove numerical effects, and assess that resolving the influence radius of the orbiting black hole is a minimum condition to fully capture the formation of the binary. Our simulations include the relevant physical processes, namely star formation, supernova feedback, accretion onto the black holes and the ensuing feedback. We find that, in these mergers, dynamical friction from the smooth stellar component of the nucleus is the main process that drives black holes from kpc to pc scales. Gas does not play a crucial role and even clumps do not induce scattering or perturb the orbits. We compare the time needed for the formation of the binary to analytical predictions and suggest how to apply such analytical formalism to obtain estimates of binary formation times in lower resolution simulations.

It is shown that the excellent Murchison Radio-astronomy Observatory site allows the Murchison Widefield Array to employ a simple RFI blanking scheme and still calibrate visibilities and form images in the FM radio band. The techniques described are running autonomously in our calibration and imaging software, which is currently being used to process an FM-band survey of the entire southern sky.

The Murchison Widefield Array (MWA) is a next-generation radio telescope currently under construction in the remote Western Australia Outback. Raw data will be generated continuously at 5GiB/s, grouped into 8s cadences. This high throughput motivates the development of on-site, real time processing and reduction in preference to archiving, transport and off-line processing. Each batch of 8s data must be completely reduced before the next batch arrives. Maintaining real time operation will require a sustained performance of around 2.5TFLOP/s (including convolutions, FFTs, interpolations and matrix multiplications). We describe a scalable heterogeneous computing pipeline implementation, exploiting both the high computing density and FLOP-per-Watt ratio of modern GPUs. The architecture is highly parallel within and across nodes, with all major processing elements performed by GPUs. Necessary scatter-gather operations along the pipeline are loosely synchronized between the nodes hosting the GPUs. The MWA will be a frontier scientific instrument and a pathfinder for planned peta- and exascale facilities.

The MWA is a next-generation radio interferometer under construction in remote Western Australia. The data rate from the correlator makes storing the raw data infeasible, so the data must be processed in real-time. The processing task is of order ~10 TFLOPS. The remote location of the MWA limits the power that can be allocated to computing. We describe the design and implementation of elements of the MWA real-time data processing system which leverage the computing abilities of modern graphics processing units (GPUs). The matrix algebra and texture mapping capabilities of GPUs are well suited to the majority of tasks involved in real-time calibration and imaging. Considerable performance advantages over a conventional CPU-based reference implementation are obtained.