The merging rate of cosmic structures is computed, relying on the Ansatz that they can be predicted in the initial linear density field from the coalescence of critical points with increasing smoothing scale, used here as a proxy for cosmic time. Beyond the mergers of peaks with saddle points (a proxy for halo mergers), we consider the coalescence and nucleation of all sets of critical points, including wall-saddle to filament-saddle and wall-saddle to minima (a proxy for filament and void mergers respectively), as they impact the geometry of galactic infall, and in particular filament disconnection. Analytical predictions of the one-point statistics are validated against multiscale measurements in 2D and 3D realisations of Gaussian random fields (the corresponding code being available upon request) and compared qualitatively to cosmological $N$-body simulations at early times ($z\geq 10$) and large scales ($\geq 5\, \mathrm{Mpc}/h$). The rate of filament coalescence is compared to the merger rate of haloes and the two-point clustering of these events is computed, along with their cross-correlations with critical points. These correlations are qualitatively consistent with the preservation of the connectivity of dark matter haloes, and the impact of the large scale structures on assembly bias. The destruction rate of haloes and voids as a function of mass and redshift is quantified down to $z=0$ for a $\Lambda$CDM cosmology. The one-point statistics in higher dimensions are also presented, together with consistency relations between critical point and critical event counts.

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$.

Extended filamentary H$\alpha$ emission nebulae are a striking feature of nearby galaxy clusters but the formation mechanism of the filaments, and the processes which shape their morphology remain unclear. We conduct an investigation into the formation, evolution and destruction of dense gas in the center of a simulated, Perseus-like, cluster under the influence of a spin-driven jet. We particularly study the role played by condensation of dense gas from the diffuse intracluster medium, and the impact of direct uplifting of existing dense gas by the jets, in determining the spatial distribution and kinematics of the dense gas. We present a hydrodynamical simulation of an idealised Perseus-like cluster using the adaptive mesh refinement code {\sc ramses}. Our simulation includes a supermassive black hole (SMBH) that self-consistently tracks its spin evolution via its local accretion, and in turn drives a large-scale jet whose direction is based on the black hole's spin evolution. We show that the formation and destruction of dense gas is closely linked to the SMBH's feedback cycle, and that its morphology is highly variable throughout the simulation. While extended filamentary structures readily condense from the hot intra-cluster medium, they are easily shattered into an overly clumpy distribution of gas during their interaction with the jet driven outflows. Condensation occurs predominantly onto infalling gas located 5 - 15 kpc from the center during quiescent phases of the central AGN, when the local ratio of the cooling time to free fall time falls below 20, i.e. when $t_{\rm cool}/t_{\rm ff} < 20$. We find evidence for both condensation and uplifting of dense gas, but caution that purely hydrodynamical simulations struggle to effectively regulate the cluster cooling cycle and produce overly clumpy distributions of dense gas morphologies, compared to observation.

We present a new implementation of the tracer particles algorithm based on a Monte Carlo approach for the Eulerian adaptive mesh refinement code Ramses. The purpose of tracer particles is to keep track of where fluid elements originate in Eulerian mesh codes, so as to follow their Lagrangian trajectories and re-processing history. We provide a comparison to the more commonly used velocity-based tracer particles, and show that the Monte Carlo approach reproduces the gas distribution much more accurately. We present a detailed statistical analysis of the properties of the distribution of tracer particles in the gas and report that it follows a Poisson law. We extend these Monte Carlo gas tracer particles to tracer particles for the stars and black holes, so that they can exchange mass back and forth between themselves. With such a scheme, we can follow the full cycle of baryons, that is, from gas-forming stars to the release of mass back to the surrounding gas multiple times, or accretion of gas onto black holes. The overall impact on computation time is $\sim$ 3% per tracer per initial cell. As a proof of concept, we study an astrophysical science case -- the dual accretion modes of galaxies at high redshifts --, which highlights how the scheme yields information hitherto unavailable. These tracer particles will allow us to study complex astrophysical systems where both efficiency of shock-capturing Godunov schemes and a Lagrangian follow-up of the fluid are required simultaneously.

The strikingly anisotropic large-scale distribution of matter made of an extended network of voids delimited by sheets, themselves segmented by filaments, within which matter flows towards compact nodes where they intersect, imprints its geometry on the dynamics of cosmic flows, ultimately shaping the distribution of galaxies and the redshift evolution of their properties. The (filament-type) saddle points of this cosmic web provide a local frame in which to quantify the induced physical and morphological evolution of galaxies on large scales. The properties of virtual galaxies within the Horizon-AGN simulation are stacked in such a frame. The iso-contours of the galactic number density, mass, specific star formation rate (sSFR), kinematics and age are clearly aligned with the filament axis with steep gradients perpendicular to the filaments. A comparison to a simulation without feedback from active galactic nuclei (AGN) illustrates its impact on quenching star formation of centrals away from the saddles. The redshift evolution of the properties of galaxies and their age distribution are consistent with the geometry of the bulk flow within that frame. They compare well with expectations from constrained Gaussian random fields and the scaling with the mass of non-linearity, modulo the redshift dependent impact of feedback processes. Physical properties such as sSFR and kinematics seem not to depend only on mean halo mass and density: the residuals trace the geometry of the saddle, which could point to other environment-sensitive physical processes, such as spin advection, and AGN feedback at high mass.

The role of the cosmic web in shaping galaxy properties is investigated in the GAMA spectroscopic survey in the redshift range $0.03 \leq z \leq 0.25$. The stellar mass, $u - r$ dust corrected colour and specific star formation rate (sSFR) of galaxies are analysed as a function of their distances to the 3D cosmic web features, such as nodes, filaments and walls, as reconstructed by DisPerSE. Significant mass and type/colour gradients are found for the whole population, with more massive and/or passive galaxies being located closer to the filament and wall than their less massive and/or star-forming counterparts. Mass segregation persists among the star-forming population alone. The red fraction of galaxies increases when closing in on nodes, and on filaments regardless of the distance to nodes. Similarly, the star-forming population reddens (or lowers its sSFR) at fixed mass when closing in on filament, implying that some quenching takes place. Comparable trends are also found in the state-of-the-art hydrodynamical simulation Horizon-AGN. These results suggest that on top of stellar mass and large-scale density, the traceless component of the tides from the anisotropic large-scale environment also shapes galactic properties. An extension of excursion theory accounting for filamentary tides provides a qualitative explanation in terms of anisotropic assembly bias: at a given mass, the accretion rate varies with the orientation and distance to filaments. It also explains the absence of type/colour gradients in the data on smaller, non-linear scales.

The mass, accretion rate and formation time of dark matter haloes near proto-filaments (identified as saddle points of the potential) are analytically predicted using a conditional version of the excursion set approach in its so-called "upcrossing" approximation. The model predicts that at fixed mass, mass accretion rate and formation time vary with orientation and distance from the saddle, demonstrating that assembly bias is indeed influenced by the tides imposed by the cosmic web. Starved, early forming haloes of smaller mass lie preferentially along the main axis of filaments, while more massive and younger haloes are found closer to the nodes. Distinct gradients for distinct tracers such as typical mass and accretion rate occur because the saddle condition is anisotropic, and because the statistics of these observables depend on both the conditional means and their covariances. The theory is extended to other critical points of the potential field. The response of the mass function to variations of the matter density field (the so-called large scale bias) is computed, and its trend with accretion rate is shown to invert along the filament. The signature of this model should correspond at low redshift to an excess of reddened galactic hosts at fixed mass along preferred directions, as recently reported in spectroscopic and photometric surveys and in hydrodynamical simulations. The anisotropy of the cosmic web emerges therefore as a significant ingredient to describe jointly the dynamics and physics of galaxies, e.g. in the context of intrinsic alignments or morphological diversity.