We present the Extreme-Horizon (EH) cosmological simulation: EH models galaxy formation with stellar and AGN feedback and uses a very high resolution in the intergalactic and circumgalactic medium. The high resolution in low-density regions results in smaller-size massive galaxies at redshift $z=2$, in better agreement with observations compared to other simulations. This results from the improved modeling of cold gas flows accreting onto galaxies. Besides, the EH simulation forms a population of particularly compact galaxies with stellar masses of $10^{10-11}$\,M$_\sun$ that are reminiscent of observed ultracompact galaxies at $z\simeq2$. These objects form mainly through repeated major mergers of low-mass progenitors, independently of baryonic feedback mechanisms. This formation process can be missed in simulations using a too low resolution in low-density intergalactic regions.

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.

The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for M_star<10^10.5 Msun. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.

Horizon Run 5 (HR5) is a cosmological hydrodynamics simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1 kpc. This enormous dynamic range allows us to simultaneously capture the physics of the cosmic web on very large scales and account for the formation and evolution of dwarf galaxies on much smaller scales. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of 1049x114x114 Mpc^3. The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, reionization, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and active galactic nuclei (AGN) feedback in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OpenMP version of the RAMSES code, specifically targeted for modern many-core many thread parallel architectures. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder to analyse the large outputs of HR5. The simulation successfully reproduces many observations, such as the cosmic star formation history, connectivity of galaxy distribution and stellar mass functions. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one should carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe.

We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light profiles in Hubble Space Telescope images. Our measurements from a large sample of massive clusters presented here provide new constraints on dark matter and cluster astrophysics. We find that DM haloes are on average highly elongated with the mean ellipticity of $0.482\pm 0.028$, and position angles of major axes of DM haloes and their BCGs tend to be aligned well with the mean value of alignment angles of $22.2\pm 3.9$ deg. We find that DM haloes in our sample are on average more elongated than their BCGs with the mean difference of their ellipticities of $0.11\pm 0.03$. In contrast, the Horizon-AGN cosmological hydrodynamical simulation predicts on average similar ellipticities between DM haloes and their central galaxies. While such a difference between the observations and the simulation may well be explained by the difference of their halo mass scales, other possibilities include the bias inherent to strong lensing measurements, limited knowledge of baryon physics, or a limitation of cold dark matter.

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.

Galaxy merger histories correlate strongly with stellar mass, largely regardless of morphology. Thus, at fixed stellar mass, spheroids and discs share similar assembly histories, both in terms of the frequency of mergers and the distribution of their mass ratios. Since mergers are the principal drivers of disc-to-spheroid morphological transformation, and the most massive galaxies typically have the richest merger histories, it is surprising that discs exist at all at the highest stellar masses (e.g. beyond the knee of the mass function). Using Horizon-AGN, a cosmological hydro-dynamical simulation, we show that extremely massive (M*> 10^11.4 MSun) discs are created via two channels. In the primary channel (accounting for ~70% of these systems and ~8% of massive galaxies) the most recent, significant merger (stellar mass ratio > 1:10) between a massive spheroid and a gas-rich satellite `spins up' the spheroid by creating a new rotational stellar component, leaving a massive disc as the remnant. In the secondary channel (accounting for ~30% of these systems and ~3% of massive galaxies), a system maintains a disc throughout its lifetime, due to an anomalously quiet merger history. Not unexpectedly, the fraction of massive discs is larger at higher redshift, due to the Universe being more gas-rich. The morphological mix of galaxies at the highest stellar masses is, therefore, a strong function of the gas fraction of the Universe. Finally, these massive discs have similar black-hole masses and accretion rates to massive spheroids, providing a natural explanation for why a minority of powerful AGN are surprisingly found in disc galaxies.

Improving distance measurements in large imaging surveys is a major challenge to better reveal the distribution of galaxies on a large scale and to link galaxy properties with their environments. Photometric redshifts can be efficiently combined with the cosmic web (CW) extracted from overlapping spectroscopic surveys to improve their accuracy. We apply a similar method using a new generation of photometric redshifts based on a convolution neural network (CNN). The CNN is trained on the SDSS images with the main galaxy sample (SDSS-MGS, $r \leq 17.8$) and the GAMA spectroscopic redshifts up tor $\sim 19.8$. The mapping of the CW is obtained with 680,000 spectroscopic redshifts from the MGS and BOSS surveys. The redshift probability distribution functions (PDF), which are well calibrated (unbiased and narrow, $\leq 120$ Mpc), intercept a few CW structure along the line of sight. Combining these PDFs with the density field distribution provides new photometric redshifts, $z_{web}$, whose accuracy is improved by a factor of two (i.e.,${\sigma} \sim 0.004(1+z)$) for galaxies with $r \leq 17.8$. For half of them, the distance accuracy is better than 10 cMpc. The narrower the original PDF, the larger the boost in accuracy. No gain is observed for original PDFs wider than 0.03. The final $z_{web}$ PDFs also appear well calibrated. The method performs slightly better for passive galaxies than star-forming ones, and for galaxies in massive groups since these populations better trace the underlying large-scale structure. Reducing the spectroscopic sampling by a factor of 8 still improves the photometric redshift accuracy by 25%. Extending the method to galaxies fainter than the MGS limit still improves the redshift estimates for 70% of the galaxies, with a gain in accuracy of 20% at low $z$ where the resolution of the CW is the highest.

Using the Horizon-AGN hydrodynamical simulation and self-organising maps (SOMs), we show how to compress the complex data structure of a cosmological simulation into a 2-d grid which is much easier to analyse. We first verify the tight correlation between the observed 0.3$\!-\!5\mu$m broad-band colours of Horizon-AGN galaxies and their high-resolution spectra. The correlation is found to extend to physical properties such as redshift, stellar mass, and star formation rate (SFR). This direct mapping from colour to physical parameter space is shown to work also after including photometric uncertainties that mimic the COSMOS survey. We then label the SOM grid with a simulated calibration sample and estimate redshift and SFR for COSMOS-like galaxies up to $z\sim3$. In comparison to state-of-the-art techniques based on synthetic templates, our method is comparable in performance but less biased at estimating redshifts, and significantly better at predicting SFRs. In particular our "data-driven" approach, in contrast to model libraries, intrinsically allows for the complexity of galaxy formation and can handle sample biases. We advocate that obtaining the calibration for this method should be one of the goals of next-generation galaxy surveys.

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

The intrinsic correlations of galaxy shapes and orientations across the large-scale structure of the Universe are a known contaminant to weak gravitational lensing. They are known to be dependent on galaxy properties, such as their mass and morphologies. The complex interplay between alignments and the physical processes that drive galaxy evolution remains vastly unexplored. We assess the sensitivity of intrinsic alignments (shapes and angular momenta) to Active Galactic Nuclei -AGN- feedback by comparing galaxy alignment in twin runs of the cosmological hydrodynamical Horizon simulation, which do and do not include AGN feedback respectively. We measure intrinsic alignments in three dimensions and in projection at z=0 and z=1. We find that the projected alignment signal of all galaxies with resolved shapes with respect to the density field in the simulation is robust to AGN feedback, thus giving similar predictions for contamination to weak lensing. The relative alignment of galaxy shapes around galaxy positions is however significantly impacted, especially when considering high-mass ellipsoids. Using a sample of galaxy "twins" across simulations, we determine that AGN changes both the galaxy selection and their actual alignments. Finally, we measure the alignments of angular momenta of galaxies with their nearest filament. Overall, these are more significant in the presence of AGN as a result of the higher abundance of massive pressure-supported galaxies.

High-redshift quasars are believed to reside in highly biased regions of the Universe, where black hole (BH) growth is sustained by an enhanced number of mergers and by being at the intersection of filaments bringing fresh gas. This assumption should be supported by an enhancement of the number counts of galaxies in the field of view of quasars. While the current observations of quasar environments do not lead to a consensus on a possible excess of galaxies, the future missions JWST, WFIRST, and Euclid will provide new insights on quasar environments, and will substantially increase the number of study-cases. We are in a crucial period, where we need to both understand the current observations and predict how upcoming missions will improve our understanding of BH environments. Using the large-scale simulation Horizon-AGN, we find that statistically the most massive BHs reside in environments with the largest number counts of galaxies. However, we find a large variance in galaxy number counts, and some massive BHs do not show enhanced counts in their neighborhood. Interestingly, some massive BHs have a very close galaxy companion but no further enhancement of the galaxy number counts at larger scales, in agreement with recent observations. We find that AGN feedback in the surrounding galaxies is able to decrease their luminosity and stellar mass, and therefore to make them un-observable when using restrictive galaxy selection criteria. Radiation from the quasars can spread over large distances, which could affect the formation history of surrounding galaxies, but a careful analysis of these processes requires radiative transfer simulations.

Recent integral field spectroscopy observations have found that about 11\% of galaxies show star-gas misalignment. The misalignment possibly results from external effects such as gas accretion, interaction with other objects, and other environmental effects, hence providing clues to these effects. We explore the properties of misaligned galaxies using Horizon-AGN, a large-volume cosmological simulation, and compare the result with the result of the Sydney-AAO Multi-object integral field spectrograph (SAMI) Galaxy Survey. Horizon-AGN can match the overall misalignment fraction and reproduces the distribution of misalignment angles found by observations surprisingly closely. The misalignment fraction is found to be highly correlated with galaxy morphology both in observations and in the simulation: early-type galaxies are substantially more frequently misaligned than late-type galaxies. The gas fraction is another important factor associated with misalignment in the sense that misalignment increases with decreasing gas fraction. However, there is a significant discrepancy between the SAMI and Horizon-AGN data in the misalignment fraction for the galaxies in dense (cluster) environments. We discuss possible origins of misalignment and disagreement.

It is known observationally that the major axes of galaxy clusters and their brightest cluster galaxies are roughly aligned with each other. To understand the origin of the alignment, we identify 40 cluster-sized dark matter (DM) haloes with masses higher than $5\times10^{13}~M_{\odot}$ and their central galaxies (CGs) at $z\approx 0$ in the Horizon-AGN cosmological hydrodynamical simulation. We trace the progenitors at 50 different epochs between $0<z<5$. We then fit their shapes and orientations with a triaxial ellipsoid model. While the orientations of both DM haloes and CGs change significantly due to repeated mergers and mass accretions, their relative orientations are well aligned at each epoch even at high redshifts, $z>1$. The alignment becomes tighter with cosmic time; the major axes of the CGs and their host DM haloes at present are aligned on average within $\sim 30^{\circ}$ in the three dimensional space and $\sim 20^{\circ}$ in the projected plane. The orientations of the major axes of DM haloes on average follow one of the eigen-vectors of the surrounding tidal field that corresponds to the {\it slowest collapsing} (or even stretching) mode, and the alignment with the tidal field also becomes tighter. This implies that the orientations of CGs and DM haloes at the present epoch are largely imprinted in the primordial density field of the Universe, whereas strong dynamical interactions such as mergers are important to explain their mutual alignment at each epoch.

Elliptical galaxies today appear aligned with the large-scale structure of the Universe, but it is still an open question when they acquire this alignment. Observational data is currently insufficient to provide constraints on the time evolution of intrinsic alignments, and hence existing models range from assuming that galaxies gain some primordial alignment at formation, to suggesting that they react instantaneously to tidal interactions with the large-scale structure. Using the cosmological hydrodynamical simulation Horizon-AGN, we measure the relative alignments between the major axes of galaxies and eigenvectors of the tidal field as a function of redshift. We focus on constraining the time evolution of the alignment of the main progenitors of massive $z=0$ elliptical galaxies, the main weak lensing contaminant at low redshift. We show that this population, which at $z=0$ has a stellar mass above $10^{10.4}$ M$_\odot$, transitions from having no alignment with the tidal field at $z=3$, to a significant alignment by $z=1$. From $z=0.5$ they preserve their alignment at an approximately constant level until $z=0$. We find a mass-dependence of the alignment signal of elliptical progenitors, whereby ellipticals that are less massive today ($10^{10.4}<M/{\rm M}_\odot<10^{10.7}$) do not become aligned till later redshifts ($z<2$), compared to more massive counterparts. We also present an extended study of progenitor alignments in the parameter space of stellar mass and galaxy dynamics, the impact of shape definition and tidal field smoothing.

Mapping of the large-scale structure through cosmic time has numerous applications in the studies of cosmology and galaxy evolution. At $z > 2$, the structure can be traced by the neutral intergalactic medium (IGM) by way of observing the Ly$\alpha$, forest towards densely-sampled lines-of-sight of bright background sources, such as quasars and star forming galaxies. We investigate the scientific potential of MOSAIC, a planned multi-object spectrograph on the European Extremely Large Telescope (ELT), for the 3D mapping of the IGM at $z \gtrsim 3$. We simulate a survey of $3 \lesssim z \lesssim 4$ galaxies down to a limiting magnitude of $m_{r}\sim 25.5$ mag in an area of 1 degree$^2$ in the sky. Galaxies and their spectra (including the line-of-sight Ly$\alpha$ absorption) are taken from the lightcone extracted from the Horizon-AGN cosmological hydrodynamical simulation. The quality of the reconstruction of the original density field is studied for different spectral resolutions and signal-to-noise ratios of the spectra. We demonstrate that the minimum $S/N$ (per resolution element) of the faintest galaxies that such survey has to reach is $S/N = 4$. We show that a survey with such sensitivity enables a robust extraction of cosmic filaments and the detection of the theoretically-predicted galaxy stellar mass and star-formation rate gradients towards filaments. By simulating the realistic performance of MOSAIC we obtain $S/N(T_{\rm obs}, R, m_{r})$ scaling relations. We estimate that $\lesssim 35~(65)$ nights of observation time are required to carry out the survey with the instrument's high multiplex mode and with the spectral resolution of $R=1000~(2000)$. A survey with a MOSAIC-concept instrument on the ELT is found to enable the mapping of the IGM at $z > 3$ on Mpc scales, and as such will be complementary to and competitive with other planned IGM tomography surveys. [abridged]

Hydrodynamical cosmological simulations have recently made great advances in reproducing galaxy mass assembly over cosmic time - as often quantified from the comparison of their predicted stellar mass functions to observed stellar mass functions from data. In this paper we compare the clustering of galaxies from the hydrodynamical cosmological simulated lightcone Horizon-AGN, to clustering measurements from the VIDEO survey observations. Using mocks built from a VIDEO-like photometry, we first explore the bias introduced into clustering measurements by using stellar masses and redshifts derived from SED-fitting, rather than the intrinsic values. The propagation of redshift and mass statistical and systematic uncertainties in the clustering measurements causes us to underestimate the clustering amplitude. We find then that clustering and halo occupation distribution (HOD) modelling results are qualitatively similar in Horizon-AGN and VIDEO. However at low stellar masses Horizon-AGN underestimates the observed clustering by up to a factor of ~3, reflecting the known excess stellar mass to halo mass ratio for Horizon-AGN low mass haloes, already discussed in previous works. This reinforces the need for stronger regulation of star formation in low mass haloes in the simulation. Finally, the comparison of the stellar mass to halo mass ratio in the simulated catalogue, inferred from angular clustering, to that directly measured from the simulation, validates HOD modelling of clustering as a probe of the galaxy-halo connection.

We investigate the impact of the number of filaments connected to the nodes of the cosmic web on the physical properties of their galaxies using the Sloan Digital Sky Survey. We compare these measurements to the cosmological hydrodynamical simulations Horizon-(no)AGN and Simba. We find that more massive galaxies are more connected, in qualitative agreement with theoretical predictions and measurements in dark matter only simulation. The star formation activity and morphology of observed galaxies both display some dependence on the connectivity of the cosmic web at fixed stellar mass: less star forming and less rotation supported galaxies also tend to have higher connectivity. These results qualitatively hold both for observed and virtual galaxies, and can be understood given that the cosmic web is the main source of fuel for galaxy growth. The simulations show the same trends at fixed halo mass, suggesting that the geometry of filamentary infall impacts galaxy properties beyond the depth of the local potential well. Based on simulations, it is also found that AGN feedback is key in reversing the relationship between stellar mass and connectivity at fixed halo mass. Technically, connectivity is a practical observational proxy for past and present accretion (minor mergers or diffuse infall).

We present the first detection of mass dependent galactic spin alignments with local cosmic filaments with over 2 sigma confidence using IFS kinematics. The 3D network of cosmic filaments is reconstructed on Mpc scales across GAMA fields using the cosmic web extractor DisPerSe. We assign field galaxies from the SAMI survey to their nearest filament segment in 3D and estimate the degree of alignment between SAMI galaxies kinematic spin axis and their nearest filament in projection. Low-mass galaxies align their spin with their nearest filament while higher mass counterparts are more likely to display an orthogonal orientation. The stellar transition mass from the first trend to the second is bracketed between log stellar masses 10.4 and 10.9, with hints of an increase with filament scale. Consistent signals are found in the HorizonAGN cosmological hydrodynamic simulation. This supports a scenario of early angular momentum build-up in vorticity rich quadrants around filaments at low stellar mass followed by progressive flip of spins orthogonal to the cosmic filaments through mergers at high stellar mass. Conversely, we show that dark-matter only simulations post-processed with a semi-analytic model treatment of galaxy formation struggles to reproduce this alignment signal. This suggests that gas physics is key in enhancing the galaxy-filament alignment.

Understanding how rotationally-supported discs transform into dispersion-dominated spheroids is central to our comprehension of galaxy evolution. Morphological transformation is largely merger-driven. While major mergers can efficiently create spheroids, recent work has highlighted the significant role of other processes, like minor mergers, in driving morphological change. Given their rich merger histories, spheroids typically exhibit large fractions of `ex-situ' stellar mass, i.e. mass that is accreted, via mergers, from external objects. This is particularly true for the most massive galaxies, whose stellar masses typically cannot be attained without a large number of mergers. Here, we explore an unusual population of extremely massive (M* > 10^11 MSun) spheroids, in the Horizon-AGN simulation, which exhibit anomalously low ex-situ mass fractions, indicating that they form without recourse to significant merging. These systems form in a single minor-merger event (with typical merger mass ratios of 0.11 - 0.33), with a specific orbital configuration, where the satellite orbit is virtually co-planar with the disc of the massive galaxy. The merger triggers a catastrophic change in morphology, over only a few hundred Myrs, coupled with strong in-situ star formation. While this channel produces a minority (~5 per cent) of such galaxies, our study demonstrates that the formation of at least some of the most massive spheroids need not involve major mergers -- or any significant merging at all -- contrary to what is classically believed.

The origin of the disk and spheroid of galaxies has been a key open question in understanding their morphology. Using the high-resolution cosmological simulation, New Horizon, we explore kinematically decomposed disk and spheroidal components of 144 field galaxies with masses greater than $\rm 10^9\,M_{\odot}$ at $z=0.7$. The origins of stellar particles are classified according to their birthplace (in situ or ex situ) and their orbits at birth. Before disk settling, stars form mainly through chaotic mergers between proto-galaxies and become part of the spheroidal component. When disk settling starts, we find that more massive galaxies begin to form disk stars from earlier epochs; massive galaxies commence to develop their disks at $z\sim1-2$, while low-mass galaxies do after $z\sim1$. The formation of disks is affected by accretion as well, as mergers can trigger gas turbulence or induce misaligned gas infall that prevents galaxies from forming co-rotating disk stars. The importance of accreted stars is greater in more massive galaxies, especially in developing massive spheroids. A significant fraction of the spheroids comes from the disk stars that are perturbed, which becomes more important at lower redshifts. Some ($\sim12.5\%$) of our massive galaxies develop counter-rotating disks from the gas infall misaligned with the existing disk plane, which can last for more than a Gyr until they become the dominant component, and flip the angular momentum of the galaxy in the opposite direction. The final disk-to-total ratio of a galaxy needs to be understood in relation to its stellar mass and accretion history. We quantify the significance of the stars with different origins and provide them as guiding values.

We study the Intra-Halo Stellar Component (IHSC) of Milky Way-mass systems up to galaxy clusters in the Horizon-AGN cosmological hydrodynamical simulation. We identify the IHSC using an improved phase-space galaxy finder algorithm which provides an adaptive, physically motivated and shape-independent definition of this stellar component, that can be applied to halos of arbitrary masses. We explore the IHSC mass fraction-total halo's stellar mass, $f_{M*,IHSC}-M*$, relation and the physical drivers of its scatter. We find that on average the $f_{M*,IHSC}$ increases with $M_{*,tot}$, with the scatter decreasing strongly with mass from 2 dex at $M_{*,tot}\sim10^{11}M_\odot$ to 0.3 dex at group masses. At high masses, $M_{*,tot}>10^{11.5}M_\odot$, $f_{M*,IHSC}$ increases with the number of substructures, and with the mass ratio between the central galaxy and largest satellite, at fixed $M_{*,tot}$. From mid-size groups and systems below $M_{*,tot}<10^{12}M_\odot$, we find that the central galaxy's stellar rotation-to-dispersion velocity ratio, V/{\sigma}, displays the strongest (anti)-correlation with $f_{M*,IHSC}$ at fixed $M_{*,tot}$ of all the galaxy and halo properties explored, transitioning from $f_{M*,IHSC}$<0.1% for high V/{\sigma}, to $f_{M*,IHSC}\sim5$% for low V/{\sigma} galaxies. By studying the $f_{M*,IHSC}$ temporal evolution, we find that, in the former, mergers not always take place, but if they did, they happened early (z>1), while the high $f_{M*,IHSC}$ population displays a much more active merger history. In the case of massive groups and galaxy clusters, $M_{*,tot}>10^{12}M_\odot$, a fraction $f_{M*,IHSC}\sim$10-20% is reached at $z\sim1$ and then they evolve across lines of constant $f_{M*,IHSC}$ modulo some small perturbations. Because of the limited simulation's volume, the latter is only tentative and requires a larger sample of simulated galaxy clusters to confirm.

Cosmological hydrodynamical simulations are rich tools to understand the build-up of stellar mass and angular momentum in galaxies, but require some level of calibration to observations. We compare predictions at $z\sim0$ from the Eagle, Hydrangea, Horizon-AGN, and Magneticum simulations with integral field spectroscopic (IFS) data from the SAMI Galaxy Survey, ATLAS3D, CALIFA and MASSIVE surveys. The main goal of this work is to simultaneously compare structural, dynamical, and stellar population measurements in order to identify key areas of success and tension. We have taken great care to ensure that our simulated measurement methods match the observational methods as closely as possible. We find that the Eagle and Hydrangea simulations reproduce many galaxy relations but with some offsets at high stellar masses. There are moderate mismatches in $R_e$ (+), $\epsilon$ (-), $\sigma_e$ (-), and mean stellar age (+), where a plus sign indicates that quantities are too high on average, and minus sign too low. The Horizon-AGN simulations qualitatively reproduce several galaxy relations, but there are a number of properties where we find a quantitative offset to observations. Massive galaxies are better matched to observations than galaxies at low and intermediate masses. Overall, we find mismatches in $R_e$ (+), $\epsilon$ (-), $\sigma_e$ (-) and $(V/\sigma)_e$ (-). Magneticum matches observations well: this is the only simulation where we find ellipticities typical for disk galaxies, but there are moderate differences in $\sigma_e$ (-), $(V/\sigma)_e$ (-) and mean stellar age (+). Our comparison between simulations and observational data has highlighted several areas for improvement, such as the need for improved modelling resulting in a better vertical disk structure, yet our results demonstrate the vast improvement of cosmological simulations in recent years.

We present analysis of two- and three-point correlation functions of Ly$\alpha$ forest (at $2\le z\le 2.5$) using X-Shooter spectra of three background quasar triplets probing transverse separations of 0.5-1.6 pMpc. We present statistics based on transmitted flux and clouds identified using Voigt profile fitting. We show that the observed two-, three-point correlation functions and reduced three-point correlation (i.e Q) are well reproduced by our simulations. We assign probabilities for realising all the observed correlation properties simultaneously using our simulations. Our simulations suggest an increase in correlation amplitudes and Q with increasing $N_{\rm HI}$. We roughly see this trend in the observations too. We identify a concurrent gap of 17$\mathring{A}$ (i.e 14.2 $h^{-1}$cMpc, one of the longest reported) wide in one of the triplets. Such gap is realised only in 14.2% of our simulated sightlines and most of the time belongs to a void in the matter distribution. In the second triplet, we detect DLAs along all three sightlines (with spatial separations 0.64 to 1.6 pMpc) within a narrow redshift interval (i.e $\Delta z = 0.088$). Detection of a foreground quasar ($\sim$ 1 pMpc from the triplet sightlines) and excess partial Lyman Limit systems around these DLAs suggest that we may be probing a large over-dense region. We also report positive CIV- CIV correlations up to $\sim 500$ $km s^{-1}$ only in the longitudinal direction. Additionally, we conclude a positive CIV-Ly$\alpha$ correlations for higher $N_{\rm HI}$ thresholds up to a scale of $\sim 1000$ $km s^{-1}$ both in transverse and longitudinal directions.

We study the spin alignment of galaxies and halos with respect to filaments and walls of the cosmic web, identified with DisPerSE, using the SIMBA simulation from z=0-2. Massive halos' spins are oriented perpendicularly to their closest filament's axis and walls, while low mass halos tend to have their spins parallel to filaments and in the plane of walls. A similar mass-dependent spin flip is found for galaxies, albeit with a weaker signal particularly at low mass and low-z, suggesting that galaxies' spins retain memory of their larger-scale environment. Low-z star-forming and rotation-dominated galaxies tend to have spins parallel to nearby filaments, while quiescent and dispersion-dominated galaxies show preferentially perpendicular orientation; the star formation trend can be fully explained by the stellar mass correlation, but the morphology trend cannot. There is a strong dependence on HI mass, such that high-HI galaxies tend to have parallel spins while low-HI galaxies are perpendicular, which persists even when matching samples in stellar mass, suggesting that HI content traces anisotropic infall more faithfully than the stellar component. Finally, at fixed stellar mass, the strength of spin alignments correlates with the filament's density, with parallel alignment for galaxies in high density environments. These findings are consistent with conditional tidal torque theory, and highlight a significant correlation between galactic spin and the larger scale tides that are important e.g. for interpreting weak lensing studies. SIMBA allows us to rule out numerical grid locking as the cause of previously-seen low mass alignment.

Cosmic filaments are the channel through which galaxy groups assemble their mass. Cosmic connectivity, namely the number of filaments connected to a given group, is therefore expected to be an important ingredient in shaping group properties. The local connectivity is measured in COSMOS around X-Ray detected groups between redshift 0.5 and 1.2. To this end, large-scale filaments are extracted using the accurate photometric redshifts of the COSMOS2015 catalogue in two-dimensional slices of thickness 120 comoving Mpc centred on the group's redshift. The link between connectivity, group mass and the properties of the brightest group galaxy (BGG) is investigated. The same measurement is carried out on mocks extracted from the lightcone of the hydrodynamical simulation Horizon-AGN in order to control systematics. More massive groups are on average more connected. At fixed group mass in low-mass groups, BGG mass is slightly enhanced at high connectivity, while in high mass groups BGG mass is lower at higher connectivity. Groups with a star-forming BGG have on average a lower connectivity at given mass. From the analysis of the Horizon-AGN simulation, we postulate that different connectivities trace different paths of group mass assembly: at high group mass, groups with higher connectivity are more likely to have grown through a recent major merger, which might be in turn the reason for the quenching of the BGG. Future large-field photometric surveys, such as Euclid and LSST, will be able to confirm and extend these results by probing a wider mass range and a larger variety of environment.

Context. Accurate model predictions including the physics of baryons are required to make the most of the upcoming large cosmological surveys devoted to gravitational lensing. The advent of hydrodynamical cosmological simulations enables such predictions on sufficiently sizeable volumes. Aims. Lensing quantities (deflection, shear, convergence) and their statistics (convergence power spectrum, shear correlation functions, galaxy-galaxy lensing) are computed in the past lightcone built in the Horizon-AGN hydrodynamical cosmological simulation, which implements our best knowledge on baryonic physics at the galaxy scale in order to mimic galaxy populations over cosmic time. Methods. Lensing quantities are generated over a one square degree field of view by performing multiple-lens plane ray-tracing through the lightcone, taking full advantage of the 1 kpc resolution and splitting the line of sight over 500 planes all the way to redshift z~7. Two methods are explored (standard projection of particles with adaptive smoothing, and integration of the acceleration field) to assert a good implementation. The focus is on small scales where baryons matter most. Results. Standard cosmic shear statistics are impacted at the 10% level by the baryonic component for angular scales below a few arcmin. The galaxy-galaxy lensing signal, or galaxy-shear correlation function, is consistent with measurements for the redshift z~0.5 massive galaxy population. At higher redshift z>1, the impact of magnification bias on this correlation is relevant for separations greater than 1 Mpc. Conclusions. This work is pivotal for all current and upcoming weak lensing surveys and represents a first step towards building a full end-to-end generation of lensed mock images from large cosmological hydrodynamical simulations.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Using the ligthcone from the cosmological hydrodynamical simulation Horizon-AGN, we produced a photometric catalogue over $0<z<4$ with apparent magnitudes in COSMOS, DES, LSST-like, and Euclid-like filters at depths comparable to these surveys. The virtual photometry accounts for the complex star formation history and metal enrichment of Horizon-AGN galaxies, and consistently includes magnitude errors, dust attenuation and absorption by inter-galactic medium. The COSMOS-like photometry is fitted in the same configuration as the COSMOS2015 catalogue. We then quantify random and systematic errors of photometric redshifts, stellar masses, and star-formation rates (SFR). Photometric redshifts and redshift errors capture the same dependencies on magnitude and redshift as found in COSMOS2015, excluding the impact of source extraction. COSMOS-like stellar masses are well recovered with a dispersion typically lower than 0.1 dex. The simple star formation histories and metallicities of the templates induce a systematic underestimation of stellar masses at $z<1.5$ by at most 0.12 dex. SFR estimates exhibit a dust-induced bimodality combined with a larger scatter (typically between 0.2 and 0.6 dex). We also use our mock catalogue to predict photometric redshifts and stellar masses in future imaging surveys. We stress that adding Euclid near-infrared photometry to the LSST-like baseline improves redshift accuracy especially at the faint end and decreases the outlier fraction by a factor $\sim$2. It also considerably improves stellar masses, reducing the scatter up to a factor 3. It would therefore be mutually beneficial for LSST and Euclid to work in synergy.

Rotation is ubiquitous in the Universe, and recent kinematic surveys have shown that early type galaxies and globular clusters are no exception. Yet the linear response of spheroidal rotating stellar systems has seldom been studied. This paper takes a step in this direction by considering the behaviour of spherically symmetric systems with differential rotation. Specifically, the stability of several sequences of Plummer spheres is investigated, in which the total angular momentum, as well as the degree and flavour of anisotropy in the velocity space are varied. To that end, the response matrix method is customised to spherical rotating equilibria. The shapes, pattern speeds and growth rates of the systems' unstable modes are computed. Detailed comparisons to appropriate N-body measurements are also presented. The marginal stability boundary is charted in the parameter space of velocity anisotropy and rotation rate. When rotation is introduced, two sequences of growing modes are identified corresponding to radially and tangentially-biased anisotropic spheres respectively. For radially anisotropic spheres, growing modes occur on two intersecting surfaces (in the parameter space of anisotropy and rotation), which correspond to fast and slow modes, depending on the net rotation rate. Generalised, approximate stability criteria are finally presented.

Our statistical understanding of galaxy evolution is fundamentally driven by objects that lie above the surface-brightness limits of current wide-area surveys (mu ~ 23 mag arcsec^-2). While both theory and small, deep surveys have hinted at a rich population of low-surface-brightness galaxies (LSBGs) fainter than these limits, their formation remains poorly understood. We use Horizon-AGN, a cosmological hydrodynamical simulation to study how LSBGs, and in particular the population of ultra-diffuse galaxies (UDGs; mu > 24.5 mag arcsec^-2), form and evolve over time. For M* > 10^8 MSun, LSBGs contribute 47, 7 and 6 per cent of the local number, mass and luminosity densities respectively (~85/11/10 per cent for M* > 10^7 MSun). Today's LSBGs have similar dark-matter fractions and angular momenta to high-surface-brightness galaxies (HSBGs; mu < 23 mag arcsec^-2), but larger effective radii (x2.5 for UDGs) and lower fractions of dense, star-forming gas (more than x6 less in UDGs than HSBGs). LSBGs originate from the same progenitors as HSBGs at z > 2. However, LSBG progenitors form stars more rapidly at early epochs. The higher resultant rate of supernova-energy injection flattens their gas-density profiles, which, in turn, creates shallower stellar profiles that are more susceptible to tidal processes. After z ~ 1, tidal perturbations broaden LSBG stellar distributions and heat their cold gas, creating the diffuse, largely gas-poor LSBGs seen today. In clusters, ram-pressure stripping provides an additional mechanism that assists in gas removal in LSBG progenitors. Our results offer insights into the formation of a galaxy population that is central to a complete understanding of galaxy evolution, and which will be a key topic of research using new and forthcoming deep-wide surveys.

Using the two large cosmological hydrodynamical simulations, Horizon-AGN (H-AGN) and Horizon-noAGN (H-noAGN, no AGN feedback), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the total density profiles (dark matter + stars) at the effective radius of massive early-type galaxies (M*>10^11 Msun). We have studied the dependencies of the mass-weighted density slope gamma'_tot with the effective radius, the galaxy mass and the host halo mass at z~0.3 and found that the inclusion of AGN feedbackalways leads to a much better agreement with observational values and trends. Our analysis suggests also that the inclusion of AGN feedback favours a strong correlation between gamma'_tot and the density slope of the dark matter component while, in the absence of AGN activity, gamma'_tot is rather strongly correlated with the density slope of the stellar component. Finally, we find that gamma'_tot derived from our samples of galaxies increases from z=2 to z=0,in good agreement with the expected observational trend. The derived slopes are slightly lower than in the data when AGN is included because the simulated galaxies tend to be too extended, especially the least massive ones. However, the simulated compact galaxies without AGN feedback have gamma'_tot values that are significantly too high compared to observations.

Using the linearized Boltzmann equation, we investigate how grooves carved in the phase space of a half-mass Mestel disc can trigger the vigorous growth of two-armed spiral eigenmodes. Such grooves result from the collisional dynamics of a disc subject to finite-N shot noise, as swing-amplified noise patterns push stars towards lower-angular momentum orbits at their inner Lindblad radius. Supplementing the linear theory with analytical arguments, we show that the dominant spiral mode is a cavity mode with reflections off the forbidden region around corotation and off the deepest groove. Other subdominant modes are identified as groove modes. We provide evidence that the depletion of near-circular orbits, and not the addition of radial orbits, is the crucial physical ingredient that causes these new eigenmodes. Thus, it is possible for an isolated, linearly stable stellar disc to spontaneously become linearly unstable via the self-induced formation of phase-space grooves through finite-N dynamics. These results may help explain the growth and maintenance of spiral patterns in real disc galaxies.

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.

Identifying galaxies in hydrodynamical simulations is a difficult task, particularly in regions of high density such as galaxy groups and clusters. We present a new scale-free shape-independent algorithm to robustly and accurately identify galaxies in simulation, implemented within the phase-space halo-finder code VELOCIraptor. This is achieved by using the full phase-space dispersion tensor for particle assignment and an iterative adjustment of search parameters, which help us overcome common structure finding problems. We apply our improved method to the Horizon-AGN simulation and compare galaxy stellar masses ($M_*$), star formation rates (SFR) and sizes with the elaborate configuration-space halo finder, HaloMaker. Galaxies living in halos with $> 1$ galaxy are the most affected by the shortcomings of real-space finders, with their mass, SFR, and sizes being $> 2$ times larger (smaller) in the case of host (satellite) galaxies. Thus, our ability to measure minor/major merger rates and disentangle environmental effects in simulations can be generally hindered if the identification of galaxies is not treated carefully. Though large systematic differences are obtained on a one-to-one basis, the overall Galaxy Stellar Mass Function, the Star Formation Rate Function and mass-size relations are not greatly affected. This is due to isolated galaxies being the most abundant population, dominating broad statistics.

Cosmic connectivity and multiplicity, i.e. the number of filaments globally or locally connected to a given cluster is a natural probe of the growth of structure and in particular of the nature of dark energy. It is also a critical ingredient driving the assembly history of galaxies as it controls mass and angular momentum accretion. The connectivity of the cosmic web is investigated here via the persistent skeleton. This tool identifies topologically the set of ridges of the cosmic landscape which allows us to investigate in details how the nodes of the cosmic web are connected together. When applied to Gaussian random fields, it is found that on average the nodes are connected to exactly $\kappa=4$ neighbours in two dimensions and ~6.1 in three dimensions. Investigating spatial dimensions up to d=11, typical departures from a cubic lattice $\kappa=2d$ are shown to scale like the power 7/4 of the dimension. These numbers strongly depend on the height of the peaks: the higher the peak the larger the connectivity. Locally, the number of filaments sticking out from a node is not equal to its number of connections. At small distance, the peak has an ellipsoidal shape so that only two filaments emerge. Further away, filaments are subdivided at bifurcation points. At the resolution of the maps, the local number of filaments is 3 on average in 2D and 4 in 3D, depending on the height of the peak. Predictions from first principles based on peak theory are shown to reproduce well the connectivity and multiplicity of Gaussian random fields and cosmological simulations. As an illustration, this connectivity is quantified in galaxy lensing convergence maps and large dark haloes catalogues. As a function of redshift and scale the mean connectivity decreases in a cosmology dependent way. As a function of halo mass it scales like 10/3 times the log of the mass.

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.

Accurate cosmology from upcoming weak lensing surveys relies on knowledge of the total matter power spectrum at percent level at scales $k < 10$ $h$/Mpc, for which modelling the impact of baryonic physics is crucial. We compare measurements of the total matter power spectrum from the Horizon cosmological hydrodynamical simulations: a dark matter-only run, one with full baryonic physics, and another lacking Active Galactic Nuclei (AGN) feedback. Baryons cause a suppression of power at $k\simeq 10$ $h/$Mpc of $<15\%$ at $z=0$, and an enhancement of a factor of a few at smaller scales due to the more efficient cooling and star formation. The results are sensitive to the presence of the highest mass haloes in the simulation and the distribution of dark matter is also impacted up to a few percent. The redshift evolution of the effect is non-monotonic throughout $z=0-5$ due to an interplay between AGN feedback and gas pressure, and the growth of structure. We investigate the effectiveness of an analytic `baryonic correction model' in describing our results. We require a different redshift evolution and propose an alternative fitting function with $4$ free parameters that reproduces our results within $5\%$. Compared to other simulations, we find the impact of baryonic processes on the total matter power spectrum to be smaller at $z=0$. Correspondingly, our results suggest that AGN feedback is not strong enough in the simulation. Total matter power spectra from the Horizon simulations are made publicly available at https://www.horizon-simulation.org/catalogues.html

This work investigates the alignment of galactic spins with the cosmic web across cosmic time using the cosmological hydrodynamical simulation Horizon-AGN. The cosmic web structure is extracted via the persistent skeleton as implemented in the DISPERSE algorithm. It is found that the spin of low-mass galaxies is more likely to be aligned with the filaments of the cosmic web and to lie within the plane of the walls while more massive galaxies tend to have a spin perpendicular to the axis of the filaments and to the walls. The mass transition is detected with a significance of 9 sigmas. This galactic alignment is consistent with the alignment of the spin of dark haloes found in pure dark matter simulations and with predictions from (anisotropic) tidal torque theory. However, unlike haloes, the alignment of low-mass galaxies is weak and disappears at low redshifts while the orthogonal spin orientation of massive galaxies is strong and increases with time, probably as a result of mergers. At fixed mass, alignments are correlated with galaxy morphology: the high-redshift alignment is dominated by spiral galaxies while elliptical centrals are mainly responsible for the perpendicular signal. These predictions for spin alignments with respect to cosmic filaments and unprecendently walls are successfully compared with existing observations. The alignment of the shape of galaxies with the different components of the cosmic web is also investigated. A coherent and stronger signal is found in terms of shape at high mass. The two regimes probed in this work induce competing galactic alignment signals for weak lensing, with opposite redshift and luminosity evolution. Understanding the details of these intrinsic alignments will be key to exploit future major cosmic shear surveys like Euclid or LSST.

The classical theory of cluster relaxation is unsatisfactory because it involves the Coulomb logarithm. The Balescu-Lenard (BL) equation provides a rigorous alternative that has no ill-defined parameter. Moreover, the BL equation, unlike classical theory, includes the cluster's self-gravity. A heuristic argument is given that indicates that relaxation does not occur predominantly through two-particle scattering and is enhanced by self-gravity. The BL equation is adapted to a spherical system and used to estimate the flux through the action space of isochrone clusters with different velocity anisotropies. A range of fairly different secular behaviours is found depending on the fraction of radial orbits. Classical theory is also used to compute the corresponding classical fluxes. The BL and classical fluxes are very different because (a) the classical theory materially under-estimates the impact of large-scale collectively amplified fluctuations and (b) only the leading terms in an infinite sum for the BL flux are computed. A complete theory of cluster relaxation likely requires that the sum in the BL equation be decomposed into a sum over a finite number of small wavenumbers complemented by an integral over large wavenumbers analogous to classical theory.

Energetic feedback from active galactic nuclei (AGNs) is often used in simulations to resolve several outstanding issues in galaxy formation, but its impact is still not fully understood. Here we derive new constraints on AGN feedback by comparing observations and simulations of the thermal Sunyaev-Zel'dovich (tSZ) effect. We draw on observational results presented in Spacek et al. (2016, 2017) who used data from the South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) to measure the tSZ signal from >= 10^11 M_Sun and >= 1 Gyr galaxies at z=0.5-1.0 (low-z) and z=1.0-1.5 (high-z). Using the large-scale cosmological hydrodynamical simulations Horizon-AGN and Horizon-NoAGN, which include and omit AGN feedback, we extract simulated tSZ measurements around galaxies equivalent to the observational work. We find that the Horizon-AGN results only differ from the SPT measurements at levels of 0.4 sigma (low-z) and 0.6 sigma (high-z), but differ from the ACT measurements by 3.4 sigma (low-z) and 2.3 sigma (high-z). The Horizon-NoAGN results provide a slightly better fit to the SPT measurements by differing by 0.2 sigma (low-z) and 0.4 sigma (high-z), but a significantly better match to the ACT measurements by differing by only 0.5 sigma (low-z) and 1.4 sigma (high-z). We conclude that, while the lower-mass (<~ 5 x 10^11 M_Sun) SPT results allow for the presence AGN feedback energy, the higher-mass (>~ 5 x 10^11 M_Sun) ACT results show significantly less energy than predicted in the simulation including AGN feedback, while more closely matching the simulation without AGN feedback, indicating that AGN feedback may be milder than often predicted in simulations.

Understanding the processes that trigger morphological transformation is central to understanding how and why the Universe transitions from being disc-dominated at early epochs to having the morphological mix that is observed today. We use Horizon-AGN, a cosmological hydrodynamical simulation, to perform a comprehensive study of the processes that drive morphological change in massive (M > 10^10 MSun) galaxies over cosmic time. We show that (1) essentially all the morphological evolution in galaxies that are spheroids at z=0 is driven by mergers with mass ratios greater than 1:10, (2) major mergers alone cannot produce today's spheroid population -- minor mergers are responsible for a third of all morphological transformation over cosmic time and are its dominant driver after z~1, (3) prograde mergers trigger milder morphological transformation than retrograde mergers -- while both types of events produce similar morphological changes at z>2, the average change due to retrograde mergers is around twice that due to their prograde counterparts at z~0, (4) remnant morphology depends strongly on the gas fraction of a merger, with gas-rich mergers routinely re-growing discs, and (5) at a given stellar mass, discs do not exhibit drastically different merger histories from spheroids -- disc survival in mergers is driven by acquisition of cold gas (via cosmological accretion and gas-rich interactions) and a preponderance of prograde mergers in their merger histories.

Using the Horizon-AGN simulation data, we study the relative role of mergers and environmental effects in shaping the spin of early-type galaxies (ETGs) after $z \simeq 1$. We follow the spin evolution of 10,037 color-selected ETGs more massive than 10$^{10} \rm \, M_{\odot}$ that are divided into four groups: cluster centrals (3%), cluster satellites (33%), group centrals (5%), and field ETGs (59%). We find a strong mass dependence of the slow rotator fraction, $f_{\rm SR}$, and the mean spin of massive ETGs. Although we do not find a clear environmental dependence of $f_{\rm SR}$, a weak trend is seen in the mean value of spin parameter driven by the satellite ETGs as they gradually lose their spin as their environment becomes denser. Galaxy mergers appear to be the main cause of total spin changes in 94% of central ETGs of halos with $M_{vir} > 10^{12.5}\rm M_{\odot}$, but only 22% of satellite and field ETGs. We find that non-merger induced tidal perturbations better correlate with the galaxy spin-down in satellite ETGs than mergers. Given that the majority of ETGs are not central in dense environments, we conclude that non-merger tidal perturbation effects played a key role in the spin evolution of ETGs observed in the local ($z < 1$) universe.

We quantify the quenching impact of the group environment using the spectroscopic survey Galaxy and Mass Assembly (GAMA) to z=0.2. The fraction of red (quiescent) galaxies, whether in groups or isolated, increases with both stellar mass and large-scale (5 Mpc) density. At fixed stellar mass, the red fraction is on average higher for satellites of red centrals than of blue (star-forming) centrals, a galactic conformity effect that increases with density. Most of the signal originates from groups that have the highest stellar mass, reside in the densest environments, and have massive, red only centrals. Assuming a color-dependent halo-to-stellar-mass ratio, whereby red central galaxies inhabit significantly more massive halos than blue ones of the same stellar mass, two regimes emerge more distinctly: at log(Mhalo/Msol) < 13, central quenching is still ongoing, conformity is no longer existent, and satellites and group centrals exhibit the same quenching excess over field galaxies at all mass and density, in agreement with the concept of "group quenching"; at log(Mhalo/Msol) > 13, a cutoff that sets apart massive (log(M*/Msol) > 11), fully quenched group centrals, conformity is meaningless, and satellites undergo significantly more quenching than their counterparts in smaller halos. The latter effect strongly increases with density, giving rise to the density-dependent conformity signal when both regimes are mixed. The star-formation of blue satellites in massive halos is also suppressed compared to blue field galaxies, while blue group centrals and the majority of blue satellites, which reside in low mass halos, show no deviation from the color-stellar mass relation of blue field galaxies.

Understanding the processes that drive the formation of black holes (BHs) is a key topic in observational cosmology. While the observed $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ correlation in bulge-dominated galaxies is thought to be produced by major mergers, the existence of a $M_{\mathrm{BH}}$--$M_{\star}$ relation, across all galaxy morphological types, suggests that BHs may be largely built by secular processes. Recent evidence that bulge-less galaxies, which are unlikely to have had significant mergers, are offset from the $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ relation, but lie on the $M_{\mathrm{BH}}$--$M_{\star}$ relation, has strengthened this hypothesis. Nevertheless, the small size and heterogeneity of current datasets, coupled with the difficulty in measuring precise BH masses, makes it challenging to address this issue using empirical studies alone. Here, we use Horizon-AGN, a cosmological hydrodynamical simulation to probe the role of mergers in BH growth over cosmic time. We show that (1) as suggested by observations, simulated bulge-less galaxies lie offset from the main $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ relation, but on the $M_{\mathrm{BH}}$--$M_{\star}$ relation, (2) the positions of galaxies on the $M_{\mathrm{BH}}$--$M_{\star}$ relation are not affected by their merger histories and (3) only $\sim$35 per cent of the BH mass in today's massive galaxies is directly attributable to merging -- the majority ($\sim$65 per cent) of BH growth, therefore, takes place gradually, via secular processes, over cosmic time.

The anisotropic distribution of satellites around the central galaxy of their host halo is well-documented. However the relative impact of baryons and dark matter in shaping this distribution is still debated. Using the simulation Horizon-AGN, the angular distribution of satellite galaxies with respect to their central counterpart and halo is quantified. Below one Rvir, satellites cluster more strongly in the plane of the central, rather than merely tracing the shape of their host halo. This is due to the increased isotropy of inner haloes acquired through their inside-out assembly in vorticity-rich flows along the cosmic web. While the effect of centrals decreases with distance, halos' triaxiality increases, impacting more and more the satellite's distribution. Effects become comparable just outside one virial radius. Above this scale, the filamentary infall also impacts the satellites distribution, dominating above two virial radii. The central's morphology plays a governing role: the alignment w.r.t. the central plane is four times stronger in haloes hosting stellar discs than in spheroids. But the impact of the galactic plane decreases for lower satellite-to-central mass ratios, suggesting this might not hold for dwarf satellites of the Local group. The orientation of the Milky-Way's satellites traces their cosmic filament, their level of coplanarity is consistent with systems of similar mass and cosmic location in Horizon-AGN. However, the strong impact of galactic planes in massive groups and clusters bounds the likelihood of finding a relaxed region where satellites can be used to infer halo shape. The minor-to-major axis ratios for haloes with log(M0/Msun)>13.5 is underestimated by 10%. This error soars quickly to 30-40% for individual halo measurements.

As endpoints of the hierarchical mass-assembly process, the stellar populations of local early-type galaxies encode the assembly history of galaxies over cosmic time. We use Horizon-AGN, a cosmological hydrodynamical simulation, to study the merger histories of local early-type galaxies and track how the morphological mix of their progenitors evolves over time. We provide a framework for alleviating `progenitor bias' -- the bias that occurs if one uses only early-type galaxies to study the progenitor population. Early-types attain their final morphology at relatively early epochs -- by $z\sim1$, around 60 per cent of today's early-types have had their last significant merger. At all redshifts, the majority of mergers have one late-type progenitor, with late-late mergers dominating at $z>1.5$ and early-early mergers becoming significant only at $z<0.5$. Progenitor bias is severe at all but the lowest redshifts -- e.g. at $z\sim0.6$, less than 50 per cent of the stellar mass in today's early-types is actually in progenitors with early-type morphology, while, at $z\sim2$, studying only early-types misses almost all (80 per cent) of the stellar mass that eventually ends up in local early-type systems. At high redshift, almost all massive late-type galaxies, regardless of their local environment or star-formation rate, are progenitors of local early-type galaxies, as are lower-mass (M$_\star$ $<$ 10$^{10.5}$ M$_{\odot}$) late-types as long as they reside in high density environments. In this new era of large observational surveys (e.g. LSST, JWST), this study provides a framework for studying how today's early-type galaxies have been built up over cosmic time.

We present a derivation of the kinetic equation describing the secular evolution of spatially inhomogeneous systems with long-range interactions, the so-called inhomogeneous Landau equation, by relying on a functional integral formalism. We start from the BBGKY hierarchy derived from the Liouville equation. At the order ${1/N}$, where $N$ is the number of particles, the evolution of the system is characterised by its 1-body distribution function and its 2-body correlation function. Introducing associated auxiliary fields, the evolution of these quantities may be rewritten as a traditional functional integral. By functionally integrating over the 2-body autocorrelation, one obtains a new constraint connecting the 1-body DF and the auxiliary fields. When inverted, this constraint allows us to obtain the closed non-linear kinetic equation satisfied by the 1-body distribution function. This derivation provides an alternative to previous methods, either based on the direct resolution of the truncated BBGKY hierarchy or on the Klimontovich equation. It may turn out to be fruitful to derive more accurate kinetic equations, e.g., accounting for collective effects, or higher order correlation terms.

Large deviation statistics is implemented to predict the statistics of cosmic densities in cylinders applicable to photometric surveys. It yields few percent accurate analytical predictions for the one-point probability distribution function (PDF) of densities in concentric or compensated cylinders; and also captures the density-dependence of their angular clustering (cylinder bias). All predictions are found to be in excellent agreement with the cosmological simulation Horizon Run 4 in the quasi-linear regime where standard perturbation theory normally breaks down. These results are combined with a simple local bias model that relates dark matter and tracer densities in cylinders and validated on simulated halo catalogues. This formalism can be used to probe cosmology with existing and upcoming photometric surveys like DES, Euclid or WFIRST containing billions of galaxies.

The observed massive end of the galaxy stellar mass function is steeper than its predicted dark matter halo counterpart in the standard $\Lambda $CDM paradigm. In this paper, we investigate the impact of active galactic nuclei (AGN) feedback on star formation in massive galaxies. We isolate the impact of AGNs by comparing two simulations from the HORIZON suite, which are identical except that one also includes super massive black holes (SMBH), and related feedback models. This allows us to cross-identify individual galaxies between simulations and quantify the effect of AGN feedback on their properties, including stellar mass and gas outflows. We find that massive galaxies ($ \rm M_{*} \geq 10^{11} M_\odot $) are quenched by AGN feedback to the extent that their stellar masses decrease by up to 80% at $z=0$. SMBHs affect their host halo through a combination of outflows that reduce their baryonic mass, particularly for galaxies in the mass range $ \rm 10^9 M_\odot \leq M_{*} \leq 10^{11} M_\odot $, and a disruption of central gas inflows, which limits in-situ star formation. As a result, net gas inflows onto massive galaxies, $ \rm M_{*} \geq 10^{11} M_\odot $, drop by up to 70%. We measure a redshift evolution in the stellar mass ratio of twin galaxies with and without AGN feedback, with galaxies of a given stellar mass showing stronger signs of quenching earlier on. This evolution is driven by a progressive flattening of the $\rm M_{\rm SMBH}-M_* $ relation with redshift, particularly for galaxies with $\rm M_{*} \leq 10^{10} M_\odot $. $\rm M_{\rm SMBH}/M_*$ ratios decrease over time, as falling average gas densities in galaxies curb SMBH growth.

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.

A key unresolved question is the role that galaxy mergers play in driving stellar mass growth over cosmic time. Recent observational work hints at the possibility that the overall contribution of `major' mergers (mass ratios $\gtrsim$1:4) to cosmic stellar mass growth may be small, because they enhance star formation rates by relatively small amounts at high redshift, when much of today's stellar mass was assembled. However, the heterogeneity and relatively small size of today's datasets, coupled with the difficulty in identifying genuine mergers, makes it challenging to $\textit{empirically}$ quantify the merger contribution to stellar mass growth. Here, we use Horizon-AGN, a cosmological hydrodynamical simulation, to comprehensively quantify the contribution of mergers to the star formation budget over the lifetime of the Universe. We show that: (1) both major and minor mergers enhance star formation to similar amounts, (2) the fraction of star formation directly attributable to merging is small at all redshifts (e.g. $\sim$35 and $\sim$20 per cent at z$\sim$3 and z$\sim$1 respectively) and (3) only $\sim$25 per cent of today's stellar mass is directly attributable to galaxy mergers over cosmic time. Our results suggest that smooth accretion, not merging, is the dominant driver of stellar mass growth over the lifetime of the Universe.

Intrinsic alignments of galaxies are a significant astrophysical systematic affecting cosmological constraints from weak gravitational lensing. Obtaining numerical predictions from hydrodynamical simulations of expected survey volumes is expensive, and a cheaper alternative relies on populating large dark matter-only simulations with accurate models of alignments calibrated on smaller hydrodynamical runs. This requires connecting the shapes and orientations of galaxies to those of dark matter haloes and to the large-scale structure. In this paper, we characterise galaxy-halo alignments in the Horizon-AGN cosmological hydrodynamical simulation. We compare the shapes and orientations of galaxies in the redshift range $0<z<3$ to those of their embedding dark matter haloes, and to the matching haloes of a twin dark-matter only run with identical initial conditions. We find that galaxy ellipticities in general cannot be predicted directly from halo ellipticities. The mean misalignment angle between the minor axis of a galaxy and its embedding halo is a function of halo mass, with residuals arising from the dependence of alignment on galaxy type, but not on environment. Haloes are much more strongly aligned among themselves than galaxies, and they decrease their alignment towards low redshift. Galaxy alignments compete with this effect, as galaxies tend to increase their alignment with haloes towards low redshift. We discuss the implications of these results for current halo models of intrinsic alignments and suggest several avenues for improvement.

The secular thickening of a discrete self-gravitating galactic disc is investigated using the inhomogeneous multi-component Balescu-Lenard equation. The thick WKB limit for the diffusion and drift coefficients is found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields a simple double quadrature for the drift and diffusion coefficients, providing a clear understanding of the positions of maximum vertical orbital diffusion within the disc induced by the effects of a finite number of particles. When applied to a tepid stable tapered disc, the Balescu-Lenard formalism predicts the formation of ridges of resonant orbits towards larger vertical actions, as found in direct numerical simulations, but over-estimates the timescale involved in their appearance. Swing amplication is likely needed to resolve this discrepancy, as demonstrated in the case of razor-thin discs. The joint evolution of a population of giant molecular clouds within the disc may accelerate the secular disc's thickening induced by finite${-N}$ effects, but the observed number of clouds in the Milky Way does not seem to be sufficient to explain its thick disc.

The secular thickening of a self-gravitating stellar galactic disc is investigated using the dressed collisionless Fokker-Planck equation and the inhomogeneous multicomponent Balescu-Lenard equation. The thick WKB limits for the diffusion fluxes are found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields simple quadratures for the drift and diffusion coefficients, providing a clear understanding of the positions of maximum vertical orbital diffusion within the disc, induced by fluctuations either external or due to the finite number of particles. These thick limits also offer a consistent derivation of a thick disc Toomre parameter, which is shown to be exponentially boosted by the ratio of the vertical to radial scale heights. Dressed potential fluctuations within the disc statistically induce a vertical bending of a subset of resonant orbits, triggering the corresponding increase in vertical velocity dispersion. When applied to a tepid stable tapered disc perturbed by shot noise, these two frameworks reproduce qualitatively the formation of ridges of resonant orbits towards larger vertical actions, as found in direct numerical simulations, but overestimates the time-scale involved in their appearance. Swing amplification is likely needed to resolve this discrepancy, as demonstrated in the case of razor-thin discs. Other sources of thickening are also investigated, such as fading sequences of slowing bars, or the joint evolution of a population of giant molecular clouds within the disc.

Using a suite of three large cosmological hydrodynamical simulations, Horizon-AGN, Horizon-noAGN (no AGN feedback) and Horizon-DM (no baryons), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the inner density profiles of massive dark matter haloes and galaxies. Based on direct object-to-object comparisons, we find that the integrated inner mass and density slope differences between objects formed in these three simulations (hereafter, H_AGN, H_noAGN and H_DM) significantly evolve with time. More specifically, at high redshift (z~5), the mean central density profiles of H_AGN and H_noAGN dark matter haloes tend to be much steeper than their H_DM counterparts owing to the rapidly growing baryonic component and ensuing adiabatic contraction. By z~1.5, these mean halo density profiles in H_AGN have flattened, pummelled by powerful AGN activity ("quasar mode"): the integrated inner mass difference gaps with H_noAGN haloes have widened, and those with H_DM haloes have narrowed. Fast forward 9.5 billion years, down to z=0, and the trend reverses: H_AGN halo mean density profiles drift back to a more cusped shape as AGN feedback efficiency dwindles ("radio mode"), and the gaps in integrated central mass difference with H_noAGN and H_DM close and broaden respectively. On the galaxy side, the story differs noticeably. Averaged stellar profile central densities and inner slopes are monotonically reduced by AGN activity as a function of cosmic time, resulting in better agreement with local observations. As both dark matter and stellar inner density profiles respond quite sensitively to the presence of a central AGN, there is hope that future observational determinations of these quantities can be used constrain AGN feedback models.

Non-Gaussianities of dynamical origin are disentangled from primordial ones using the formalism of large deviation statistics with spherical collapse dynamics. This is achieved by relying on accurate analytical predictions for the one-point probability distribution function (PDF) and the two-point clustering of spherically-averaged cosmic densities (sphere bias). Sphere bias extends the idea of halo bias to intermediate density environments and voids as underdense regions. In the presence of primordial non-Gaussianity, sphere bias displays a strong scale dependence relevant for both high and low density regions, which is predicted analytically. The statistics of densities in spheres are built to model primordial non-Gaussianity via an initial skewness with a scale-dependence that depends on the bispectrum of the underlying model. The analytical formulas with the measured nonlinear dark matter variance as input are successfully tested against numerical simulations. For local non-Gaussianity with a range from $f_{\rm NL}=-100$ to $+100$ they are found to agree within 2\% or better for densities $\rho\in[0.5,3]$ in spheres of radius 15 Mpc$/h$ down to $z=0.35$. The validity of the large deviation statistics formalism is thereby established for all observationally relevant local-type departures from perfectly Gaussian initial conditions. The corresponding estimators for the amplitude of the nonlinear variance $\sigma_8$ and primordial skewness $f_{\rm NL}$ are validated using a fiducial joint maximum likelihood experiment. The influence of observational effects and the prospects for a future detection of primordial non-Gaussianity from joint one- and two-point densities-in-spheres statistics are discussed.

General self-consistent expressions for the coefficients of diffusion and dynamical friction in a stable, bound, multicomponent self-gravitating and inhomogeneous system are derived. They account for the detailed dynamics of the colliding particles and their self-consistent dressing by collective gravitational interactions. The associated Fokker-Planck equation is shown to be fully consistent with the corresponding inhomogeneous Balescu-Lenard equation and, in the weak self-gravitating limit, to the inhomogeneous Landau equation. Hence it provides an alternative derivation to both and demonstrates their equivalence. The corresponding stochastic Langevin equations are presented: they can be a practical alternative to numerically solving the inhomogeneous Fokker-Planck and Balescu-Lenard equations. The present formalism allows for a self-consistent description of the secular evolution of different populations covering a spectrum of masses, with a proper accounting of the induced secular mass segregation, which should be of interest to various astrophysical contexts, from galactic centers to protostellar discs.

Context. Upcoming weak lensing surveys such as Euclid will provide an unprecedented opportunity to quantify the geometry and topology of the cosmic web, in particular in the vicinity of lensing clusters. Aims. Understanding the connectivity of the cosmic web with unbiased mass tracers, such as weak lensing, is of prime importance to probe the underlying cosmology, seek dynamical signatures of dark matter, and quantify environmental effects on galaxy formation. Methods. Mock catalogues of galaxy clusters are extracted from the N-body PLUS simulation. For each cluster, the aperture multipolar moments of the convergence are calculated in two annuli (inside and outside the virial radius). By stacking their modulus, a statistical estimator is built to characterise the angular mass distribution around clusters. The moments are compared to predictions from perturbation theory and spherical collapse. Results. The main weakly chromatic excess of multipolar power on large scales is understood as arising from the contraction of the primordial cosmic web driven by the growing potential well of the cluster. Besides this boost, the quadrupole prevails in the cluster (ellipsoidal) core, while at the outskirts, harmonic distortions are spread on small angular modes, and trace the non-linear sharpening of the filamentary structures. Predictions for the signal amplitude as a function of the cluster-centric distance, mass, and redshift are presented. The prospects of measuring this signal are estimated for current and future lensing data sets. Conclusions. The Euclid mission should provide all the necessary information for studying the cosmic evolution of the connectivity of the cosmic web around lensing clusters using multipolar moments and probing unique signatures of, for example, baryons and warm dark matter.

Gravitational lensing allows to quantify the angular distribution of the convergence field around clusters of galaxies to constrain their connectivity to the cosmic web. We describe in this paper the corresponding theory in Lagrangian space where analytical results can be obtained by identifying clusters to peaks in the initial field. We derive the three-point Gaussian statistics of a two-dimensional field and its first and second derivatives. The formalism allows us to study the statistics of the field in a shell around a central peak, in particular its multipolar decomposition. The peak condition is shown to significantly remove power from the dipolar contribution and to modify the monopole and quadrupole. As expected, higher order multipoles are not significantly modified by the constraint. Analytical predictions are successfully checked against measurements in Gaussian random fields. The effect of substructures and radial weighting is shown to be small and does not change the qualitative picture. The non-linear evolution is shown to induce a non-linear bias of all multipoles proportional to the cluster mass.We predict the Gaussian and weakly non-Gaussian statistics of multipolar moments of a two-dimensional field around a peak as a proxy for the azimuthal distribution of the convergence field around a cluster of galaxies. A quantitative estimate of this multipolar decomposition of the convergence field around clusters in numerical simulations of structure formation and in observations will be presented in two forthcoming papers.

Starting from a very accurate model for density-in-cells statistics of dark matter based on large deviation theory, a bias model for the tracer density in spheres is formulated. It adopts a mean bias relation based on a quadratic bias model to relate the log-densities of dark matter to those of mass-weighted dark haloes in real and redshift space. The validity of the parametrised bias model is established using a parametrisation-independent extraction of the bias function. This average bias model is then combined with the dark matter PDF, neglecting any scatter around it: it nevertheless yields an excellent model for densities-in-cells statistics of mass tracers that is parametrised in terms of the underlying dark matter variance and three bias parameters. The procedure is validated on measurements of both the one and two point statistics of subhalo densities in the state-of-the-art Horizon Run 4 simulation showing excellent agreement for measured dark matter variance and bias parameters. Finally, it is demonstrated that this formalism allows for a joint estimation of the nonlinear dark matter variance and the bias parameters using solely the statistics of subhaloes. Having verified that galaxy counts in hydrodynamical simulations sampled on a scale of 10 Mpc/h closely resemble those of subhaloes, this work provides important steps towards making theoretical predictions for density-in-cells statistics applicable to upcoming galaxy surveys like Euclid or WFIRST.

We have measured the probability distribution function (PDF) of cosmic matter density field from a suite of N-body simulations. We propose the generalized normal distribution of version 2 (Nv2) as an alternative fitting formula to the well-known log-normal distribution. We find that Nv2 provides significantly better fit than the log-normal distribution for all smoothing radii (2, 5, 10, 25 [Mpc/h]) that we studied. The improvement is substantial in the underdense regions. The development of non- Gaissianities in the cosmic matter density field is captured by continuous evolution of the skewness and shifts parameters of the Nv2 distribution. We present the redshift evolution of these parameters for aforementioned smoothing radii and various background cosmology models. All the PDFs measured from large and high-resolution N-body simulations that we use in this study can be obtained from a Web site at https://astro.kias.re.kr/jhshin.

Simple parameter-free analytic bias functions for the two-point correlation of densities in spheres at large separation are presented. These bias functions generalize the so-called Kaiser bias to the mildly non-linear regime for arbitrary density contrasts. The derivation is carried out in the context of large deviation statistics while relying on the spherical collapse model. A logarithmic transformation provides a saddle approximation which is valid for the whole range of densities and shown to be accurate against the 30 Gpc cube state-of-the-art Horizon Run 4 simulation. Special configurations of two concentric spheres that allow to identify peaks are employed to obtain the conditional bias and a proxy to BBKS extrema correlation functions. These analytic bias functions should be used jointly with extended perturbation theory to predict two-point clustering statistics as they capture the non-linear regime of structure formation at the percent level down to scales of about 10 Mpc/h at redshift 0. Conversely, the joint statistics also provide us with optimal dark matter two-point correlation estimates which can be applied either universally to all spheres or to a restricted set of biased (over- or underdense) pairs. Based on a simple fiducial survey, this estimator is shown to perform five times better than usual two-point function estimators. Extracting more information from correlations of different types of objects should prove essential in the context of upcoming surveys like Euclid, DESI, PFS or LSST.

The analytical formalism to obtain the probability distribution functions (PDFs) of spherically-averaged cosmic densities and velocity divergences in the mildly non-linear regime is presented. A large-deviation principle is applied to those cosmic fields assuming their most likely dynamics in spheres is set by the spherical collapse model. We validate our analytical results using state-of-the-art dark matter simulations with a phase-space resolved velocity field finding a 2% percent level agreement for a wide range of velocity divergences and densities in the mildly nonlinear regime (~10Mpc/h at redshift zero), usually inaccessible to perturbation theory. From the joint PDF of densities and velocity divergences measured in two concentric spheres, we extract with the same accuracy velocity profiles and conditional velocity PDF subject to a given over/under-density which are of interest to understand the non-linear evolution of velocity flows. Both PDFs are used to build a simple but accurate maximum likelihood estimators for the redshift evolution of the variance of both the density and velocity divergence fields, which have smaller relative errors than their sample variances when non-linearities appear. Given the dependence of the velocity divergence on the growth rate, there is a significant gain in using the full knowledge of both PDFs to derive constraints on the equation of state of dark energy. Thanks to the insensitivity of the velocity divergence to bias, its PDF can be used to obtain unbiased constraints on the growth of structures ($\sigma_8$,f) or it can be combined with the galaxy density PDF to extract bias parameters.

The drift and diffusion coefficients of the inhomogeneous multi-mass degenerate Landau equation are computed to describe the self-induced resonant relaxation of a discrete self-gravitating quasi-Keplerian razor-thin axisymmetric disc orbiting a massive black hole while relying on Gauss' method. For a disc-like configuration in our Galactic centre, secular diffusion induces an adiabatic distortion of orbits. When considering a disc composed of multiple masses similarly distributed, the population of lighter stars will gain eccentricity, driving it closer to the central black hole provided the distribution function increases with angular momentum. The quenching of the diffusion of a test star in the vicinity of the black hole due to the divergence of the relativistic precessions (the "Schwarzschild barrier") is correctly recovered by the kinetic equation. The dual stochastic Langevin formulation yields consistent results and provides a versatile framework in which to incorporate other stochastic processes.

The variations of galaxy stellar masses and colour-types with the distance to projected cosmic filaments are quantified using the precise photometric redshifts of the COSMOS2015 catalogue extracted from COSMOS field (2 deg$^{2}$). Realistic mock catalogues are also extracted from the lightcone of the cosmological hydrodynamical simulation Horizon-AGN. They show that the photometric redshift accuracy of the observed catalogue ($\sigma_z<0.015$ at $M_*>10^{10}{\rm M}_{\odot}$ and $z<0.9$) is sufficient to provide 2D filaments that closely match their projected 3D counterparts. Transverse stellar mass gradients are measured in projected slices of thickness 75 Mpc between $0.5< z <0.9$, showing that the most massive galaxies are statistically closer to their neighbouring filament. At fixed stellar mass, passive galaxies are also found closer to their filament while active star-forming galaxies statistically lie further away. The contributions of nodes and local density are removed from these gradients to highlight the specific role played by the geometry of the filaments. We find that the measured signal does persist after this removal, clearly demonstrating that proximity to a filament is not equivalent to proximity to an over-density. These findings are in agreement with gradients measured both in 2D or 3D in the Horizon-AGN simulation and those observed in the spectroscopic VIPERS survey (which rely on the identification of 3D filaments). They are consistent with a picture in which the influence of the geometry of the large-scale environment drives anisotropic tides which impact the assembly history of galaxies, and hence their observed properties.

We derive the kinetic equation that describes the secular evolution of a large set of particles orbiting a dominant massive object, such as stars bound to a supermassive black hole or a proto-planetary debris disc encircling a star. Because the particles move in a quasi-Keplerian potential, their orbits can be approximated by ellipses whose orientations remain fixed over many dynamical times. The kinetic equation is obtained by simply averaging the BBGKY equations over the fast angle that describes motion along these ellipses. This so-called Balescu-Lenard equation describes self-consistently the long-term evolution of the distribution of quasi-Keplerian orbits around the central object: it models the diffusion and drift of their actions, induced through their mutual resonant interaction. Hence, it is the master equation that describes the secular effects of resonant relaxation. We show how it captures the phenonema of mass segregation and of the relativistic Schwarzschild barrier recently discovered in $N$-body simulations.

We compare the predictions of Horizon-AGN, a hydro-dynamical cosmological simulation that uses an adaptive mesh refinement code, to observational data in the redshift range 0<z<6. We study the reproduction, by the simulation, of quantities that trace the aggregate stellar-mass growth of galaxies over cosmic time: luminosity and stellar-mass functions, the star formation main sequence, rest-frame UV-optical-near infrared colours and the cosmic star-formation history. We show that Horizon-AGN, which is not tuned to reproduce the local Universe, produces good overall agreement with these quantities, from the present day to the epoch when the Universe was 5% of its current age. By comparison to Horizon-noAGN, a twin simulation without AGN feedback, we quantify how feedback from black holes is likely to help shape galaxy stellar-mass growth in the redshift range 0<z<6, particularly in the most massive galaxies. Our results demonstrate that Horizon-AGN successfully captures the evolutionary trends of observed galaxies over the lifetime of the Universe, making it an excellent tool for studying the processes that drive galaxy evolution and making predictions for the next generation of galaxy surveys.

We present the first quantitative detection of large-scale filamentary structure at $z \simeq 0.7$ in the large cosmological volume probed by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We use simulations to show the capability of VIPERS to recover robust topological features in the galaxy distribution, in particular the filamentary network. We then investigate how galaxies with different stellar masses and stellar activities are distributed around the filaments and find a significant segregation, with the most massive or quiescent galaxies being closer to the filament axis than less massive or active galaxies. The signal persists even after down-weighting the contribution of peak regions. Our results suggest that massive and quiescent galaxies assemble their stellar mass through successive mergers during their migration along filaments towards the nodes of the cosmic web. On the other hand, low-mass star-forming galaxies prefer the outer edge of filaments, a vorticity rich region dominated by smooth accretion, as predicted by the recent spin alignment theory. This emphasizes the role of large scale cosmic flows in shaping galaxy properties.

The interplay between cosmic gas accretion on to galaxies and galaxy mergers drives the observed morphological diversity of galaxies. By comparing the state-of-the-art hydrodynamical cosmological simulations Horizon-AGN and Horizon-noAGN, we unambiguously identify the critical role of active galactic nuclei (AGN) in setting up the correct galaxy morphology for the massive end of the population. With AGN feedback, typical kinematic and morpho-metric properties of galaxy populations as well as the galaxy-halo mass relation are in much better agreement with observations. Only AGN feedback allows massive galaxies at the centre of groups and clusters to become ellipticals, while without AGN feedback those galaxies reform discs. It is the merger-enhanced AGN activity that is able to freeze the morphological type of the post-merger remnant by durably quenching its quiescent star formation. Hence morphology is shown to be driven not only by mass but also by the nature of cosmic accretion: at constant galaxy mass, ellipticals are galaxies that are mainly assembled through mergers, while discs are preferentially built from the in situ star formation fed by smooth cosmic gas infall.

In the conference presentation we have reviewed the theory of non-Gaussian geometrical measures for the 3D Cosmic Web of the matter distribution in the Universe and 2D sky data, such as Cosmic Microwave Background (CMB) maps that was developed in a series of our papers. The theory leverages symmetry of isotropic statistics such as Minkowski functionals and extrema counts to develop post- Gaussian expansion of the statistics in orthogonal polynomials of invariant descriptors of the field, its first and second derivatives. The application of the approach to 2D fields defined on a spherical sky was suggested, but never rigorously developed. In this paper we present such development treating effects of the curvature and finiteness of the spherical space $S_2$ exactly, without relying on the flat-sky approximation. We present Minkowski functionals, including Euler characteristic and extrema counts to the first non-Gaussian correction, suitable for weakly non-Gaussian fields on a sphere, of which CMB is the prime example.

We use the Ly-$\alpha$ Mass Association Scheme (LyMAS; Peirani et al. 2014) to predict cross-correlations at $z=2.5$ between dark matter halos and transmitted flux in the Ly-$\alpha$ forest, and compare to cross-correlations measured for quasars and damped Ly-$\alpha$ systems (DLAs) from the Baryon Oscillation Spectroscopic Survey (BOSS) by Font-Ribera et al. (2012, 2013). We calibrate LyMAS using Horizon-AGN hydrodynamical cosmological simulations of a $(100\ h^{-1}\ \mathrm{Mpc})^3$ comoving volume. We apply this calibration to a $(1\ h^{-1}\ \mathrm{Gpc})^3$ simulation realized with $2048^3$ dark matter particles. In the 100 $h^{-1}$ Mpc box, LyMAS reproduces the halo-flux correlations computed from the full hydrodynamic gas distribution very well. In the 1 $h^{-1}$ Gpc box, the amplitude of the large scale cross-correlation tracks the halo bias $b_h$ as expected. We provide empirical fitting functions that describe our numerical results. In the transverse separation bins used for the BOSS analyses, LyMAS cross-correlation predictions follow linear theory accurately down to small scales. Fitting the BOSS measurements requires inclusion of random velocity errors; we find best-fit RMS velocity errors of 399 km s$^{-1}$ and 252 km s$^{-1}$ for quasars and DLAs, respectively. We infer bias-weighted mean halo masses of $M_h/10^{12}\ h^{-1}M_\odot=2.19^{+0.16}_{-0.15}$ and $0.69^{+0.16}_{-0.14}$ for the host halos of quasars and DLAs, with $\sim 0.2$ dex systematic uncertainty associated with redshift evolution, IGM parameters, and selection of data fitting range.

Using the cosmological hydrodynamics simulation Horizon-AGN, we investigate the spatial distribution of satellite galaxies relative to their central counterpart in the redshift range between 0.3 and 0.8. We find that, on average, these satellites tend to be located on the galactic plane of the central object. This effect is detected for central galaxies with a stellar mass larger than 10^10 solar masses and found to be strongest for red passive galaxies, while blue galaxies exhibit a weaker trend. For galaxies with a minor axis parallel to the direction of the nearest filament, we find that the coplanarity is stronger in the vicinity of the central galaxy, and decreases when moving towards the outskirts of the host halo. By contrast, the spatial distribution of satellite galaxies relative to their closest filament follows the opposite trend: their tendency to align with them dominates at large distances from the central galaxy, and fades away in its vicinity. Relying on mock catalogs of galaxies in that redshift range, we show that massive red centrals with a spin perpendicular to their filament also have corotating satellites well aligned with both the galactic plane and the filament. On the other hand, lower-mass blue centrals with a spin parallel to their filament have satellites flowing straight along this filament, and hence orthogonally to their galactic plane. The orbit of these satellites is then progressively bent towards a better alignment with the galactic plane as they penetrate the central region of their host halo. The kinematics previously described are consistent with satellite infall and spin build-up via quasi-polar flows, followed by a re-orientation of the spin of massive red galaxies through mergers.

Intrinsic galaxy shape and angular momentum alignments can arise in cosmological large-scale structure due to tidal interactions or galaxy formation processes. Cosmological hydrodynamical simulations have recently come of age as a tool to study these alignments and their contamination to weak gravitational lensing. We probe the redshift and luminosity evolution of intrinsic alignments in Horizon-AGN between z=0 and z=3 for galaxies with an r-band absolute magnitude of <-20. Alignments transition from being radial at low redshifts and high luminosities, dominated by the contribution of ellipticals, to being tangential at high redshift and low luminosities, where discs dominate the signal. This cannot be explained by the evolution of the fraction of ellipticals and discs alone: intrinsic evolution in the amplitude of alignments is necessary. The alignment amplitude of elliptical galaxies alone is smaller in amplitude by a factor of ~2, but has similar luminosity and redshift evolution as in current observations and in the nonlinear tidal alignment model at projected separations of > 1 Mpc. Alignments of discs are null in projection and consistent with current low redshift observations. The combination of the two populations yields an overall amplitude a factor of ~4 lower than observed alignments of luminous red galaxies with a steeper luminosity dependence. The restriction on accurate galaxy shapes implies that the galaxy population in the simulation is complete only to an r-band absolute magnitude of <-20. Higher resolution simulations will be necessary to avoid extrapolation of the intrinsic alignment predictions to the range of luminosities probed by future surveys.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

We analyse the demographics of black holes (BHs) in the large-volume cosmological hydrodynamical simulation Horizon-AGN. This simulation statistically models how much gas is accreted onto BHs, traces the energy deposited into their environment and, consequently, the back-reaction of the ambient medium on BH growth. The synthetic BHs reproduce a variety of observational constraints such as the redshift evolution of the BH mass density and the mass function. Strong self-regulation via AGN feedback, weak supernova feedback, and unresolved internal processes result in a tight BH-galaxy mass correlation. Starting at z~2, tidal stripping creates a small population of BHs over-massive with respect to the halo. The fraction of galaxies hosting a central BH or an AGN increases with stellar mass. The AGN fraction agrees better with multi-wavelength studies, than single-wavelength ones, unless obscuration is taken into account. The most massive halos present BH multiplicity, with additional BHs gained by ongoing or past mergers. In some cases, both a central and an off-centre AGN shine concurrently, producing a dual AGN. This dual AGN population dwindles with decreasing redshift, as found in observations. Specific accretion rate and Eddington ratio distributions are in good agreement with observational estimates. The BH population is dominated in turn by fast, slow, and very slow accretors, with transitions occurring at z=3 and z=2 respectively.

We present a first principle approach to obtain analytical predictions for spherically-averaged cosmic densities in the mildly non-linear regime that go well beyond what is usually achieved by standard perturbation theory. A large deviation principle allows us to compute the leading-order cumulants of average densities in concentric cells. In this symmetry, the spherical collapse model leads to cumulant generating functions that are robust for finite variances and free of critical points when logarithmic density transformations are implemented. They yield in turn accurate density probability distribution functions (PDFs) from a straightforward saddle-point approximation valid for all density values. Based on this easy-to-implement modification, explicit analytic formulas for the evaluation of the one- and two-cell PDF are provided. The theoretical predictions obtained for the PDFs are accurate to a few percent compared to the numerical integration, regardless of the density under consideration and in excellent agreement with N-body simulations for a wide range of densities. This formalism should prove valuable for accurately probing the quasi-linear scales of low redshift surveys for arbitrary primordial power spectra.

In order to quantify the error budget in the measured probability distribution functions of cell densities, the two-point statistics of cosmic densities in concentric spheres is investigated. Bias functions are introduced as the ratio of their two-point correlation function to the two-point correlation of the underlying dark matter distribution. They describe how cell densities are spatially correlated. They are computed here via the so-called large deviation principle in the quasi-linear regime. Their large-separation limit is presented and successfully compared to simulations for density and density slopes: this regime is shown to be rapidly reached allowing to get sub-percent precision for a wide range of densities and variances. The corresponding asymptotic limit provides an estimate of the cosmic variance of standard concentric cell statistics applied to finite surveys. More generally, no assumption on the separation is required for some specific moments of the two-point statistics, for instance when predicting the generating function of cumulants containing any powers of concentric densities in one location and one power of density at some arbitrary distance from the rest. This exact "one external leg" cumulant generating function is used in particular to probe the rate of convergence of the large-separation approximation.

The recently published analytic probability density function for the mildly non-linear cosmic density field within spherical cells is used to build a simple but accurate maximum likelihood estimate for the redshift evolution of the variance of the density, which, as expected, is shown to have smaller relative error than the sample variance. This estimator provides a competitive probe for the equation of state of dark energy, reaching a few percent accuracy on wp and wa for a Euclid-like survey. The corresponding likelihood function can take into account the configuration of the cells via their relative separations. A code to compute one-cell density probability density functions for arbitrary initial power spectrum, top-hat smoothing and various spherical collapse dynamics is made available online so as to provide straightforward means of testing the effect of alternative dark energy models and initial power-spectra on the low-redshift matter distribution.

We present the COSMOS2015 catalog which contains precise photometric redshifts and stellar masses for more than half a million objects over the 2deg$^{2}$ COSMOS field. Including new $YJHK_{\rm s}$ images from the UltraVISTA-DR2 survey, $Y$-band from Subaru/Hyper-Suprime-Cam and infrared data from the Spitzer Large Area Survey with the Hyper-Suprime-Cam Spitzer legacy program, this near-infrared-selected catalog is highly optimized for the study of galaxy evolution and environments in the early Universe. To maximise catalog completeness for bluer objects and at higher redshifts, objects have been detected on a $\chi^{2}$ sum of the $YJHK_{\rm s}$ and $z^{++}$ images. The catalog contains $\sim 6\times 10^5$ objects in the 1.5 deg$^{2}$ UltraVISTA-DR2 region, and $\sim 1.5\times 10^5$ objects are detected in the "ultra-deep stripes" (0.62 deg$^{2}$) at $K_{\rm s}\leq 24.7$ (3$\sigma$, 3", AB magnitude). Through a comparison with the zCOSMOS-bright spectroscopic redshifts, we measure a photometric redshift precision of $\sigma_{\Delta z/(1+z_s)}$ = 0.007 and a catastrophic failure fraction of $\eta=0.5$%. At $3<z<6$, using the unique database of spectroscopic redshifts in COSMOS, we find $\sigma_{\Delta z/(1+z_s)}$ = 0.021 and $\eta=13.2\% $. The deepest regions reach a 90\% completeness limit of 10$^{10}M_\odot$ to $z=4$. Detailed comparisons of the color distributions, number counts, and clustering show excellent agreement with the literature in the same mass ranges. COSMOS2015 represents a unique, publicly available, valuable resource with which to investigate the evolution of galaxies within their environment back to the earliest stages of the history of the Universe. The COSMOS2015 catalog is distributed via anonymous ftp (ftp://ftp.iap.fr/pub/from_users/hjmcc/COSMOS2015/) and through the usual astronomical archive systems (CDS, ESO Phase 3, IRSA).

The Lagrangian peaks of a 1D cosmological random field representing dark matter are used as a proxy for a catalogue of biased tracers in order to investigate the small-scale exclusion in the two-halo term. The two-point correlation function of peaks of a given height is numerically estimated and analytical approximations that are valid inside the exclusion zone are derived. The resulting power spectrum of these tracers is investigated and shows clear deviations from Poisson noise at low frequencies. On large scales, the convergence of a perturbative bias expansion is discussed. Finally, we go beyond Gaussian statistics for the initial conditions and investigate the subsequent evolution of the two-point clustering of peaks through their Zel'dovich ballistic displacement, to clarify how exclusion effects mix up with scale-dependencies induced by nonlinear gravitational evolution. While the expected large-scale separation limit is recovered, significant deviations are found in the exclusion zone that tends in particular to be reduced at later times. Even though these findings apply to the clustering of one-dimensional tracers, they provide useful insights into halo exclusion and its impact on the two-halo term.

The intrinsic alignments of galaxies are recognised as a contaminant to weak gravitational lensing measurements. In this work, we study the alignment of galaxy shapes and spins at low redshift ($z\sim 0.5$) in Horizon-AGN, an adaptive-mesh-refinement hydrodynamical cosmological simulation box of 100 Mpc/h a side with AGN feedback implementation. We find that spheroidal galaxies in the simulation show a tendency to be aligned radially towards over-densities in the dark matter density field and other spheroidals. This trend is in agreement with observations, but the amplitude of the signal depends strongly on how shapes are measured and how galaxies are selected in the simulation. Disc galaxies show a tendency to be oriented tangentially around spheroidals in three-dimensions. While this signal seems suppressed in projection, this does not guarantee that disc alignments can be safely ignored in future weak lensing surveys. The shape alignments of luminous galaxies in Horizon-AGN are in agreement with observations and other simulation works, but we find less alignment for lower luminosity populations. We also characterize the systematics of galaxy shapes in the simulation and show that they can be safely neglected when measuring the correlation of the density field and galaxy ellipticities.

The secular evolution of an infinitely thin tepid isolated galactic disc made of a finite number of particles is investigated using the inhomogeneous Balescu-Lenard equation expressed in terms of angle-action variables. The matrix method is implemented numerically in order to model the induced gravitational polarization. Special care is taken to account for the amplification of potential fluctuations of mutually resonant orbits and the unwinding of the induced swing amplified transients. Quantitative comparisons with ${N-}$body simulations yield consistent scalings with the number of particles and with the self-gravity of the disc: the fewer particles and the colder the disc, the faster the secular evolution. Secular evolution is driven by resonances, but does not depend on the initial phases of the disc. For a Mestel disc with ${Q \sim 1.5}$, the polarization cloud around each star boosts up its secular effect by a factor of the order of a thousand or more, promoting accordingly the dynamical relevance of self-induced collisional secular evolution. The position and shape of the induced resonant ridge are found to be in very good agreement with the prediction of the Balescu-Lenard equation, which scales with the square of the susceptibility of the disc. In astrophysics, the inhomogeneous Balescu-Lenard equation may describe the secular diffusion of giant molecular clouds in galactic discs, the secular migration and segregation of planetesimals in proto-planetary discs, or even the long-term evolution of population of stars within the Galactic centre. It could be used as a valuable check of the accuracy of ${N-}$body integrators over secular timescales.

The geometry of the cosmic web drives in part the spin acquisition of galaxies. This can be explained in a Lagrangian framework, by identifying the specific long-wavelength correlations within the primordial Gaussian random field which are relevant to spin acquisition. Tidal Torque Theory is revisited in the context of such anisotropic environments, biased by the presence of a filament within a wall. The point process of filament-type saddles represents it most efficiently. The constrained misalignment between the tidal and the inertia tensors in the vicinity of filament-type saddles simply explains the distribution of spin directions. This misalignment implies in particular an azimuthal orientation for the spins of more massive galaxies and a spin alignment with the filament for less massive galaxies. This prediction is found to be in qualitative agreement with measurements in Gaussian random fields and N-body simulations. It relates the transition mass to the geometry of the saddle, and accordingly predicts its measured scaling with the mass of non-linearity. Implications for galaxy formation and weak lensing are briefly discussed, as is the dual theory of spin alignments in walls.

The secular evolution of an infinitely thin tepid isolated galactic disc made of a finite number of particles is described using the inhomogeneous Balescu-Lenard equation. Assuming that only tightly wound transient spirals are present in the disc, a WKB approximation provides a simple and tractable quadrature for the corresponding drift and diffusion coefficients. It provides insight into the physical processes at work during the secular diffusion of a self-gravitating discrete disc and makes quantitative predictions on the initial variations of the distribution function in action space. When applied to the secular evolution of an isolated stationary self-gravitating Mestel disc, this formalism predicts initially the importance of the corotation resonance in the inner regions of the disc leading to a regime involving radial migration and heating. It predicts in particular the formation of a "ridge like" feature in action space, in agreement with simulations, but over-estimates the timescale involved in its appearance. Swing amplification is likely to resolve this discrepancy. In astrophysics, the inhomogeneous Balescu-Lenard equation and its WKB limit may also describe the secular diffusion of giant molecular clouds in galactic discs, the secular migration and segregation of planetesimals in proto-planetary discs, or even the long-term evolution of population of stars within the Galactic center.

The equation describing the secular diffusion of a self-gravitating collisionless system induced by an exterior perturbation is derived while assuming that the timescale corresponding to secular evolution is much larger than that corresponding to the natural frequencies of the system. Its two dimensional formulation for a tepid galactic disc is also derived using the epicyclic approximation. Its WKB limit is found while assuming that only tightly wound transient spirals are sustained by the disc. It yields a simple quadrature for the diffusion coefficients which provides a straightforward understanding of the loci of maximal diffusion within the disc.

The main orbital signatures of the secular evolution of an isolated self-gravitating stellar Mestel disc are recovered using a dressed Fokker-Planck formalism in angle-action variables. The shot-noise-driven formation of narrow ridges of resonant orbits is recovered in the WKB limit of tightly wound transient spirals, for a tepid Toomre-stable tapered disc. The relative effect of the bulge, the halo, the disc temperature and the spectral properties of the shot noise are investigated in turn. For such galactic discs all elements seem to impact the locus and direction of the ridge. For instance, when the halo mass is decreased, we observe a transition between a regime of heating in the inner regions of the disc through the inner Lindblad resonance to a regime of radial migration of quasi-circular orbits via the corotation resonance in the outer part of the disc. The dressed secular formalism captures both the nature of collisionless systems (via their natural frequencies and susceptibility), and their nurture via the structure of the external perturbing power spectrum. Hence it provides the ideal framework in which to study their long term evolution.

Fluctuations in a stellar system's gravitational field cause the orbits of stars to evolve. The resulting evolution of the system can be computed with the orbit-averaged Fokker-Planck equation once the diffusion tensor is known. We present the formalism that enables one to compute the diffusion tensor from a given source of noise in the gravitational field when the system's dynamical response to that noise is included. In the case of a cool stellar disc we are able to reduce the computation of the diffusion tensor to a one-dimensional integral. We implement this formula for a tapered Mestel disc that is exposed to shot noise and find that we are able to explain analytically the principal features of a numerical simulation of such a disc. In particular the formation of narrow ridges of enhanced density in action space is recovered. As the disc's value of Toomre's $Q$ is reduced and the disc becomes more responsive, there is a transition from a regime of heating in the inner regions of the disc through the inner Lindblad resonance to one of radial migration of near-circular orbits via the corotation resonance in the intermediate regions of the disc. The formalism developed here provides the ideal framework in which to study the long-term evolution of all kinds of stellar discs.

Building galaxy merger trees from a state-of-the-art cosmological hydrodynamics simulation, Horizon-AGN, we perform a statistical study of how mergers and smooth accretion drive galaxy morphologic properties above $z > 1$. More specifically, we investigate how stellar densities, effective radii and shape parameters derived from the inertia tensor depend on mergers of different mass ratios. We find strong evidence that smooth accretion tends to flatten small galaxies over cosmic time, leading to the formation of disks. On the other hand, mergers, and not only the major ones, exhibit a propensity to puff up and destroy stellar disks, confirming the origin of elliptical galaxies. We also find that elliptical galaxies are more susceptible to grow in size through mergers than disc galaxies with a size-mass evolution $r \prop M^{1.2}$ instead of $r \prop M^{-0.5} - M^{0.5}$ depending on the merger mass ratio. The gas content drive the size-mass evolution due to merger with a faster size growth for gas-poor galaxies $r \prop M^2$ than for gas-rich galaxies $r \prop M$.

The intrinsic alignment of galaxy shapes (by means of their angular momentum) and their cross-correlation with the surrounding dark matter tidal field are investigated using the 160 000, z=1.2 synthetic galaxies extracted from the high-resolution cosmological hydrodynamical simulation Horizon-AGN. One- and two-point statistics of the spin of the stellar component are measured as a function of mass and colour. For the low-mass galaxies, this spin is locally aligned with the tidal field `filamentary' direction while, for the high-mass galaxies, it is perpendicular to both filaments and walls. The bluest galaxies of our synthetic catalog are more strongly correlated with the surrounding tidal field than the reddest galaxies, and this correlation extends up to 10 Mpc/h comoving distance. We also report a correlation of the projected ellipticities of blue, intermediate mass galaxies on a similar scale at a level of 10^(-4) which could be a concern for cosmic shear measurements. We do not report any measurable intrinsic alignments of the reddest galaxies of our sample. This work is a first step toward the use of very realistic catalog of synthetic galaxies to evaluate the contamination of weak lensing measurement by the intrinsic galactic alignments.

In the context of count-in-cells statistics, the joint probability distribution of the density in two concentric spherical shells is predicted from first first principle for sigmas of the order of one. The agreement with simulation is found to be excellent. This statistics allows us to deduce the conditional one dimensional probability distribution function of the slope within under dense (resp. overdense) regions, or of the density for positive or negative slopes. The former conditional distribution is likely to be more robust in constraining the cosmological parameters as the underlying dynamics is less evolved in such regions. A fiducial dark energy experiment is implemented on such counts derived from Lambda-CDM simulations.

The kinematic analysis of dark matter and hydrodynamical simulations suggests that the vorticity in large-scale structure is mostly confined to, and predominantly aligned with their filaments, with an excess of probability of 20 per cent to have the angle between vorticity and filaments direction lower than 60 degrees relative to random orientations. The cross sections of these filaments are typically partitioned into four quadrants with opposite vorticity sign, arising from multiple flows, originating from neighbouring walls. The spins of halos embedded within these filaments are consistently aligned with this vorticity for any halo mass, with a stronger alignment for the most massive structures up to an excess of probability of 165 per cent. On large scales, adiabatic/cooling hydrodynamical simulations display the same vorticity in the gas as in the dark matter. The global geometry of the flow within the cosmic web is therefore qualitatively consistent with a spin acquisition for smaller halos induced by this large-scale coherence, as argued in Codis et al. (2012). In effect, secondary anisotropic infall (originating from the vortex-rich filament within which these lower-mass halos form) dominates the angular momentum budget of these halos. The transition mass from alignment to orthogonality is related to the size of a given multi-flow region with a given polarity. This transition may be reconciled with the standard tidal torque theory if the latter is augmented so as to account for the larger scale anisotropic environment of walls and filaments.

We use Horizon-AGN, a hydrodynamical cosmological simulation, to explore the role of mergers in the evolution of massive (M > 10^10 MSun) galaxies around the epoch of peak cosmic star formation (1<z<4). The fraction of massive galaxies in major mergers (mass ratio R<4:1) is around 3%, a factor of ~2.5 lower than minor mergers (4:1<R <10:1) at these epochs, with no trend with redshift. At z~1, around a third of massive galaxies have undergone a major merger, while all such systems have undergone either a major or minor merger. While almost all major mergers at z>3 are 'blue' (i.e. have significant associated star formation), the proportion of 'red' mergers increases rapidly at z<2, with most merging systems at z~1.5 producing remnants that are red in rest-frame UV-optical colours. The star formation enhancement during major mergers is mild (~20-40%) which, together with the low incidence of such events, implies that this process is not a significant driver of early stellar mass growth. Mergers (R < 10:1) host around a quarter of the total star formation budget in this redshift range, with major mergers hosting around two-thirds of this contribution. Notwithstanding their central importance to the standard LCDM paradigm, mergers are minority players in driving star formation at the epochs where the bulk of today's stellar mass was formed.

Using the Horizon-AGN simulation we find a mass dependent spin orientation trend for galaxies: the spin of low-mass, rotation-dominated, blue, star-forming galaxies are preferentially aligned with their closest filament, whereas high-mass, velocity dispersion- supported, red quiescent galaxies tend to possess a spin perpendicular to these filaments. We explore the physical mechanisms driving galactic spin swings and quantify how much mergers and smooth accretion re-orient them relative to their host filaments and impact their shape. In particular, we analyze the effect of dispersion and morphology of galaxies and discuss potential tracers for prospective surveys.

Intrinsic alignments are believed to be a major source of systematics for future generation of weak gravitational lensing surveys like Euclid or LSST. Direct measurements of the alignment of the projected light distribution of galaxies in wide field imaging data seem to agree on a contamination at a level of a few per cent of the shear correlation functions, although the amplitude of the effect depends on the population of galaxies considered. Given this dependency, it is difficult to use dark matter-only simulations as the sole resource to predict and control intrinsic alignments. We report here estimates on the level of intrinsic alignment in the cosmological hydrodynamical simulation Horizon-AGN that could be a major source of systematic errors in weak gravitational lensing measurements. In particular, assuming that the spin of galaxies is a good proxy for their ellipticity, we show how those spins are spatially correlated and how they couple to the tidal field in which they are embedded. We also present theoretical calculations that illustrate and qualitatively explain the observed signals.

Filaments of the cosmic web drive spin acquisition of disc galaxies. The point process of filament-type saddle represent best this environment and can be used to revisit the Tidal Torque Theory in the context of an anisotropic peak (saddle) background split. The constrained misalignment between the tidal tensor and the Hessian of the density field generated in the vicinity of filament saddle points simply explains the corresponding transverse and longitudinal point-reflection symmetric geometry of spin distribution. It predicts in particular an azimuthal orientation of the spins of more massive galaxies and spin alignment with the filament for less massive galaxies. Its scale dependence also allows us to relate the transition mass corresponding to the alignment of dark matter halos spin relative to the direction of their neighboring filament to this geometry, and to predict accordingly it s scaling with the mass of non linearity, as was measured in simulations.

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

The close relationship between mergers and the reorientation of the spin for galaxies and their host dark haloes is investigated using a cosmological hydrodynamical simulation (Horizon-AGN). Through a statistical analysis of merger trees, we show that spin swings are mainly driven by mergers along the filamentary structure of the cosmic web, and that these events account for the preferred perpendicular orientation of massive galaxies with respect to their nearest filament. By contrast, low-mass galaxies (M_s<10^10 M_sun at redshift 1.5) undergoing very few mergers, if at all, tend to possess a spin well aligned with their filament. Haloes follow the same trend as galaxies but display a greater sensitivity to smooth anisotropic accretion. The relative effect of mergers on spin magnitude is qualitatively different for minor and major mergers: mergers (and diffuse accretion) generally increase the magnitude of the angular momentum, but the most massive major mergers also give rise to a population of objects with less spin left. Without mergers secular accretion builds up the spin of galaxies but not that of haloes. It also (re)aligns galaxies with their filament.

In the context of upcoming large-scale structure surveys such as Euclid, it is of prime importance to quantify the effect of peculiar velocities on geometric probes. Hence the formalism to compute in redshift space the geometrical and topological one-point statistics of mildly non-Gaussian 2D and 3D cosmic fields is developed. Leveraging the partial isotropy of the target statistics, the Gram-Charlier expansion of the joint probability distribution of the field and its derivatives is reformulated in terms of the corresponding anisotropic variables. In particular, the cosmic non-linear evolution of the Minkowski functionals, together with the statistics of extrema are investigated in turn for 3D catalogues and 2D slabs. The amplitude of the non-Gaussian redshift distortion correction is estimated for these geometric probes. In 3D, gravitational perturbation theory is implemented in redshift space to predict the cosmic evolution of all relevant Gram-Charlier coefficients. Applications to the estimation of the cosmic parameters sigma(z) and beta=f/b1 from upcoming surveys is discussed. Such statistics are of interest for anisotropic fields beyond cosmology.

[abridged] We describe LyMAS (Ly-alpha Mass Association Scheme), a method of predicting clustering statistics in the Ly-alpha forest on large scales from moderate resolution simulations of the dark matter distribution, with calibration from high-resolution hydrodynamic simulations of smaller volumes. We use the "Horizon MareNostrum" simulation, a 50 Mpc/h comoving volume evolved with the adaptive mesh hydrodynamic code RAMSES, to compute the conditional probability distribution P(F_s|delta_s) of the transmitted flux F_s, smoothed (1-dimensionally) over the spectral resolution scale, on the dark matter density contrast delta_s, smoothed (3-dimensionally) over a similar scale. In this study we adopt the spectral resolution of the SDSS-III BOSS at z=2.5, and we find optimal results for a dark matter smoothing length sigma=0.3 Mpc/h (comoving). In extended form, LyMAS exactly reproduces both the 1-dimensional power spectrum and 1-point flux distribution of the hydro simulation spectra. Applied to the MareNostrum dark matter field, LyMAS accurately predicts the 2-point conditional flux distribution and flux correlation function of the full hydro simulation for transverse sightline separations as small as 1 Mpc/h, including redshift-space distortion effects. It is substantially more accurate than a deterministic density-flux mapping ("Fluctuating Gunn-Peterson Approximation"), often used for large volume simulations of the forest. With the MareNostrum calibration, we apply LyMAS to 1024^3 N-body simulations of a 300 Mpc/h and 1.0 Gpc/h cube to produce large, publicly available catalogs of mock BOSS spectra that probe a large comoving volume. LyMAS will be a powerful tool for interpreting 3-d Ly-alpha forest data, thereby transforming measurements from BOSS and other massive quasar absorption surveys into constraints on dark energy, dark matter, space geometry, and IGM physics.

The joint probability distribution function (PDF) of the density within multiple concentric spherical cells is considered. It is shown how its cumulant generating function can be obtained at tree order in perturbation theory as the Legendre transform of a function directly built in terms of the initial moments. In the context of the upcoming generation of large-scale structure surveys, it is conjectured that this result correctly models such a function for finite values of the variance. Detailed consequences of this assumption are explored. In particular the corresponding one-cell density probability distribution at finite variance is computed for realistic power spectra, taking into account its scale variation. It is found to be in agreement with $\Lambda$-CDM simulations at the few percent level for a wide range of density values and parameters. Related explicit analytic expansions at the low and high density tails are given. The conditional (at fixed density) and marginal probability of the slope -- the density difference between adjacent cells -- and its fluctuations is also computed from the two-cells joint PDF; it also compares very well to simulations, in particular in under-dense regions, with a significant reduced cosmic scatter compared to over-dense regions. It is emphasized that this could prove useful when studying the statistical properties of voids as it can serve as a statistical indicator to test gravity models and/or probe key cosmological parameters.

The polarization transfer coefficients of a relativistic magnetized plasma are derived. These results apply to any momentum distribution function of the particles, isotropic or anisotropic. Particles interact with the radiation either in a non resonant mode when the frequency of the radiation exceeds their characteristic synchrotron emission frequency, or quasi resonantly otherwise. These two classes of particles contribute differently to the polarization transfer coefficients. For a given frequency, this dichotomy corresponds to a regime change in the dependence of the transfer coefficients on the parameters of the particle s population. The derivation of the transfer coefficients involves an exact expression of the conductivity tensor of the relativistic magnetized plasma that has not been used hitherto in this context. Suitable expansions valid at frequencies larger than the cyclotron frequency allow us to analytically perform the summation over all resonances at high harmonics of the relativistic gyrofrequency. The transfer coefficients are represented in the form of two variable integrals that can be conveniently computed for any set of parameters by using Olver s expansion of high order Bessel functions. We particularize our results to a number of distribution functions, isotropic, thermal or powerlaw, with different multipolar anisotropies of low order, or strongly beamed. For isotropic distributions, the Faraday coefficients are expressed in the form of a one variable quadrature over energy, for which we provide the kernels in the high-frequency limit and in the asymptotic low-frequency limit. A similar reduction to a one-variable quadrature over energy is derived at high frequency for a large class of anisotropic distribution functions that may form a basis on which any smoothly anisotropic distribution could be expanded.

An Adaptive Mesh Refinement cosmological resimulation is analyzed in order to test whether filamentary flows of cold gas are responsible for the build-up of angular momentum within a Milky Way like disk at z>=3. A set of algorithms is presented that takes advantage of the high spatial resolution of the simulation (12 pc) to identify: (i) the central gas disk and its plane of orientation; (ii) the complex individual filament trajectories that connect to the disk, and; (iii) the infalling satellites. The results show that two filaments at z>5.5, which later merge to form a single filament at z<4, drive the angular momentum and mass budget of the disk throughout its evolution, whereas luminous satellite mergers make negligible fractional contributions. Combined with the ubiquitous presence of such filaments in all large-scale cosmological simulations that include hydrodynamics, these findings provide strong quantitative evidence that the growth of thin disks in haloes with masses below 10^{12} M_{sun}, which host the vast majority of galaxies, is supported via inflowing streams of cold gas at intermediate and high redshifts.

Supermassive black holes (BH) are powerful sources of energy that are already in place at very early epochs of the Universe (by z=6). Using hydrodynamical simulations of the formation of a massive M_vir=5 10^11 M_sun halo by z=6 (the most massive progenitor of a cluster of M_vir=2 10^15 M_sun at z=0), we evaluate the impact of Active Galactic Nuclei (AGN) on galaxy mass content, BH self-regulation, and gas distribution inside this massive halo. We find that SN feedback has a marginal influence on the stellar structure, and no influence on the mass distribution on large scales. In contrast, AGN feedback alone is able to significantly alter the stellar-bulge mass content by quenching star formation when the BH is self-regulating, and by depleting the cold gas reservoir in the centre of the galaxy. The growth of the BH proceeds first by a rapid Eddington-limited period fed by direct cold filamentary infall. When the energy delivered by the AGN is sufficiently large to unbind the cold gas of the bulge, the accretion of gas onto the BH is maintained both by smooth gas inflow and clump migration through the galactic disc triggered by merger-induced torques. The feedback from the AGN has also a severe consequence on the baryon mass content within the halo, producing large-scale hot superwinds, able to blow away some of the cold filamentary material from the centre and reduce the baryon fraction by more than 30 per cent within the halo's virial radius. Thus in the very young universe, AGN feedback is likely to be a key process, shaping the properties of the most massive galaxies.

Using hydrodynamical zoom simulations in the standard LCDM cosmology, we investigate the evolution of the distribution of baryons (gas and stars) in a local group-type universe. First, with standard star formation and supernova feedback prescriptions, we find that the mean baryonic fraction value estimated at the virial radius of the two main central objects (i.e. the Milky Way and Andromeda) is decreasing over time, and is 10-15% lower than the universal value, 0.166, at z=0. This decrease is mainly due to the fact that the amount of accretion of dissipative gas onto the halo, especially at low redshift, is in general much lower than that of the dissipationless dark matter. Indeed, a significant part of the baryons does not collapse onto the haloes and remains in their outskirts, mainly in the form of warm-hot intergalactic medium (WHIM). Moreover, during the formation of each object, some dark matter and baryons are also expelled through merger events via tidal disruption. In contrast to baryons, expelled dark matter can be more efficiently re-accreted onto the halo, enhancing both the reduction of fb inside Rv, and the increase of the mass of WHIM outside Rv. Varying the efficiency of supernovae feedback at low redshift does not seem to significantly affect these trends. Alternatively, when a significant fraction of the initial gas in the main objects is released at high redshifts by more powerful sources of feedback, such as AGN from intermediate mass black holes in lower mass galaxies, the baryonic fraction at the virial radius can have a lower value (fb~0.12) at low redshift. Hence physical mechanisms able to drive the gas out of the virial radius at high redshifts will have a stronger impact on the deficit of baryons in the mass budget of Milky Way type-galaxies at present times than those that expel the gas in the longer, late phases of galaxy formation.

We investigate the alignment of the spin of dark matter halos relative (i) to the surrounding large-scale filamentary structure, and (ii) to the tidal tensor eigenvectors using the Horizon 4pi dark matter simulation which resolves over 43 million dark matter halos at redshift zero. We detect a clear mass transition: the spin of dark matter halos above a critical mass tends to be perpendicular to the closest filament, and aligned with the intermediate axis of the tidal tensor, whereas the spin of low-mass halos is more likely to be aligned with the closest filament. Furthermore, this critical mass of 5 10^12 is redshift-dependent and scales as (1+z)^-2.5. We propose an interpretation of this signal in terms of large-scale cosmic flows. In this picture, most low-mass halos are formed through the winding of flows embedded in misaligned walls; hence they acquire a spin parallel to the axis of the resulting filaments forming at the intersection of these walls. On the other hand, more massive halos are typically the products of later mergers along such filaments, and thus they acquire a spin perpendicular to this direction when their orbital angular momentum is converted into spin. We show that this scenario is consistent with both the measured excess probabilities of alignment w.r.t. the eigen-directions of the tidal tensor, and halo merger histories. On a more qualitative level, it also seems compatible with 3D visualization of the structure of the cosmic web as traced by "smoothed" dark matter simulations or gas tracer particles. Finally, it provides extra support to the disc forming paradigm presented by Pichon et al (2011) as it extends it by characterizing the geometry of secondary infall at high redshift.

We use cosmological hydrodynamical simulations to show that a significant fraction of the gas in high redshift rare massive halos falls nearly radially to their very centre on extremely short timescales. This process results in the formation of very compact bulges with specific angular momentum a factor 5-30$smaller than the average angular momentum of the baryons in the whole halo. Such low angular momentum originates both from segregation and effective cancellation when the gas flows to the centre of the halo along well defined cold filamentary streams. These filaments penetrate deep inside the halo and connect to the bulge from multiple rapidly changing directions. Structures falling in along the filaments (satellite galaxies) or formed by gravitational instabilities triggered by the inflow (star clusters) further reduce the angular momentum of the gas in the bulge. Finally, the fraction of gas radially falling to the centre appears to increase with the mass of the halo; we argue that this is most likely due to an enhanced cancellation of angular momentum in rarer halos which are fed by more isotropically distributed cold streams. Such an increasingly efficient funnelling of low-angular momentum gas to the centre of very massive halos at high redshift may account for the rapid pace at which the most massive supermassive black holes grow to reach observed masses around $10^9$M$_\odot$ at an epoch when the Universe is barely 1 Gyr old.

The formalism to compute the geometrical and topological one-point statistics of mildly non-Gaussian 2D and 3D cosmological fields is developed. Leveraging the isotropy of the target statistics, the Gram-Charlier expansion is reformulated with rotation invariant variables. This formulation allows us to track the geometrical statistics of the cosmic field to all orders. It then allows us to connect the one point statistics of the critical sets to the growth factor through perturbation theory, which predicts the redshift evolution of higher order cumulants. In particular, the cosmic non-linear evolution of the skeleton's length, together with the statistics of extrema and Euler characteristic are investigated in turn. In 2D, the corresponding differential densities are analytic as a function of the excursion set threshold and the shape parameter. In 3D, the Euler characteristics and the field isosurface area are also analytic to all orders in the expansion. Numerical integrations are performed and simple fits are provided whenever closed form expressions are not available. These statistics are compared to estimates from N-body simulations and are shown to match well the cosmic evolution up to root mean square of the density field of ~0.2. In 3D, gravitational perturbation theory is implemented to predict the cosmic evolution of all the relevant Gram-Charlier coefficients for universes with scale invariant matter distribution. The one point statistics of critical sets could be used to constrain primordial non-Gaussianities and the dark energy equation of state on upcoming cosmic surveys; this is illustrated on idealized experiments.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

State-of-the-art hydrodynamical simulations show that gas inflow through the virial sphere of dark matter halos is focused (i.e. has a preferred inflow direction), consistent (i.e. its orientation is steady in time) and amplified (i.e. the amplitude of its advected specific angular momentum increases with time). We explain this to be a consequence of the dynamics of the cosmic web within the neighbourhood of the halo, which produces steady, angular momentum rich, filamentary inflow of cold gas. On large scales, the dynamics within neighbouring patches drives matter out of the surrounding voids, into walls and filaments before it finally gets accreted onto virialised dark matter halos. As these walls/filaments constitute the boundaries of asymmetric voids, they acquire a net transverse motion, which explains the angular momentum rich nature of the later infall which comes from further away. We conjecture that this large-scale driven consistency explains why cold flows are so efficient at building up high redshift thin discs from the inside out.

The imprint of Baryonic Acoustic Oscillations (BAO) on the matter power spectrum can be constrained using the neutral hydrogen density in the intergalactic medium as a tracer of the matter density. One of the goals of the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey (SDSS-III) is to derive the Hubble expansion rate and the angular scale from the BAO signal in the IGM. To this aim, the Lyman-alpha forest of 10^5 quasars will be observed in the redshift range 2.2<z<3.5 and over 10,000 deg^2. We simulated the BOSS QSO survey to estimate the statistical accuracy on the BAO scale determination provided by such a large scale survey. In particular, we discuss the effect of the poorly constrained estimate of the unabsorbed intrinsic quasar spectrum. The volume of current N-body simulations being too small for such studies, we resorted to Gaussian random field (GRF) simulations. We validated the use of GRFs by comparing the output of GRF simulations with that of the Horizon N-body simulation with the same initial conditions. Realistic mock samples of QSO Lyman-\alpha forest were generated; their power spectrum was computed and fitted to obtain the BAO scale. The rms of the results for 100 different simulations provides an estimate of the statistical error expected from the BOSS survey. We confirm the results from Fisher matrix estimate. In the absence of error on the unabsorbed quasar spectrum, the BOSS quasar survey should measure the BAO scale with an error of the order of 2.3%, or the transverse and radial BAO scales separately with errors of the order of 6.8% and 3.9%, respectively. The significance of the BAO detection is assessed by an average \Delta\chi^2=17 but for individual realizations \Delta\chi^2 ranges from 2 t o 35. The error on the unabsorbed quasar spectrum increases the error on the BAO scale by 10 to 20% and results in a sub percent bias.

In the context of the geometrical analysis of weakly non Gaussian CMB maps, the 2D differential extrema counts as functions of the excursion set threshold is derived from the full moments expansion of the joint probability distribution of an isotropic random field, its gradient and invariants of the Hessian. Analytic expressions for these counts are given to second order in the non Gaussian correction, while a Monte Carlo method to compute them to arbitrary order is presented. Matching count statistics to these estimators is illustrated on fiducial non-Gaussian "Planck" data.

Recent theoretical studies have shown that galaxies at high redshift are fed by cold, dense gas filaments, suggesting angular momentum transport by gas differs from that by dark matter. Revisiting this issue using high-resolution cosmological hydrodynamics simulations with adaptive mesh refinement, we find that at the time of accretion, gas and dark matter do carry a similar amount of specific angular momentum, but that it is systematically higher than that of the dark matter halo as a whole. At high redshift, freshly accreted gas rapidly streams into the central region of the halo, directly depositing this large amount of angular momentum within a sphere of radius r=0.1rvir. In contrast, dark matter particles pass through the central region unscathed, and a fraction of them ends up populating the outer regions of the halo (r/rvir>0.1), redistributing angular momentum in the process. As a result, large-scale motions of the cosmic web have to be considered as the origin of gas angular momentum rather than its virialised dark matter halo host. This generic result holds for halos of all masses at all redshifts, as radiative cooling ensures that a significant fraction of baryons remain trapped at the centre of the halos. Despite this injection of angular momentum enriched gas, we predict an amount for stellar discs which is in fair agreement with observations at z=0. This arises because the total specific angular momentum of the baryons remains close to that of dark matter halos. We propose a new scenario where gas efficiently carries the angular momentum generated by large-scale structure motions deep inside dark matter halos, redistributing it only in the vicinity of the disc.

From a Principal Component Analysis (PCA) of 78 z~3 high quality quasar spectra in the SDSS-DR7, we derive the principal components characterizing the QSO continuum over the full wavelength range available. The shape of the mean continuum, is similar to that measured at low-z (z~1), but the equivalent width of the emission lines are larger at low redshift. We calculate the correlation between fluxes at different wavelengths and find that the emission line fluxes in the red part of the spectrum are correlated with that in the blue part. We construct a projection matrix to predict the continuum in the Lyman-$\alpha$ forest from the red part of the spectrum. We apply this matrix to quasars in the SDSS-DR7 to derive the evolution with redshift of the mean flux in the Lyman-$\alpha$ forest due to the absorption by the intergalactic neutral hydrogen. A change in the evolution of the mean flux is apparent around z~3 in the sense of a steeper decrease of the mean flux at higher redshifts. The same evolution is found when the continuum is estimated from the extrapolation of a power-law continuum fitted in the red part of the quasar spectrum if a correction, derived from simple simulations, is applied. Our findings are consistent with previous determinations using high spectral resolution data. We provide the PCA eigenvectors over the wavelength range 1020-2000 \AA\ and the distribution of their weights that can be used to simulate QSO mock spectra.

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

[Abridged] Cold gas flowing within the "cosmic web" is believed to be an important source of fuel for star formation at high redshift. However, the presence of such filamentary gas has never been observationally confirmed. In this work, we investigate in detail whether such cold gas is detectable using low-ionisation metal absorption lines, such as CII \lambda1334 as this technique has a proven observational record for detecting gaseous structures. Using a large statistical sample of galaxies from the Mare Nostrum N-body+AMR cosmological simulation, we find that the typical covering fraction of the dense, cold gas in 10^12 Msun haloes at z~2.5 is lower than expected (~5%). In addition, the absorption signal by the interstellar medium of the galaxy itself turns out to be so deep and so broad in velocity space that it completely drowns that of the filamentary gas. A detectable signal might be obtained from a cold filament exactly aligned with the line of sight, but this configuration is so unlikely that it would require surveying an overwhelmingly large number of candidate galaxies to tease it out. Finally, the predicted metallicity of the cold gas in filaments is extremely low (\leq 0.001 Zsun). Should this result persist when higher resolution runs are performed, it would significantly increase the difficulty of detecting filamentary gas inflows using metal lines. However, even if we assume that filaments are enriched to Zsun, the absorption signal that we compute is still weak. We are therefore led to conclude that it is extremely difficult to observationally prove or disprove the presence of cold filaments as the favorite accretion mode of galaxies using low-ionisation metal absorption lines. The Ly-alpha emission route looks more promising but due to the resonant nature of the line, radiative transfer simulations are required to fully characterize the observed signal.

We present an analytical calculation of the extreme value statistics for dark matter halos - that is, the probability distribution of the most massive halo within some region of the universe of specified shape and size. Our calculation makes use of the counts-in-cells formalism for the correlation functions, and the halo bias derived from the Sheth-Tormen mass function. We demonstrate the power of the method on spherical regions, comparing the results to measurements in a large cosmological dark matter simulation and achieving good agreement. Particularly good fits are obtained for the most likely value of the maximum mass and for the high-mass tail of the distribution, relevant in constraining cosmologies by observations of most massive clusters.

The recently introduced discrete persistent structure extractor (DisPerSE, Soubie 2010, paper I) is implemented on realistic 3D cosmological simulations and observed redshift catalogues (SDSS); it is found that DisPerSE traces equally well the observed filaments, walls, and voids in both cases. In either setting, filaments are shown to connect onto halos, outskirt walls, which circumvent voids. Indeed this algorithm operates directly on the particles without assuming anything about the distribution, and yields a natural (topologically motivated) self-consistent criterion for selecting the significance level of the identified structures. It is shown that this extraction is possible even for very sparsely sampled point processes, as a function of the persistence ratio. Hence astrophysicists should be in a position to trace and measure precisely the filaments, walls and voids from such samples and assess the confidence of the post-processed sets as a function of this threshold, which can be expressed relative to the expected amplitude of shot noise. In a cosmic framework, this criterion is comparable to friend of friend for the identifications of peaks, while it also identifies the connected filaments and walls, and quantitatively recovers the full set of topological invariants (Betti numbers) {\sl directly from the particles} as a function of the persistence threshold. This criterion is found to be sufficient even if one particle out of two is noise, when the persistence ratio is set to 3-sigma or more. The algorithm is also implemented on the SDSS catalogue and used to locat interesting configurations of the filamentary structure. In this context we carried the identification of an ``optically faint'' cluster at the intersection of filaments through the recent observation of its X-ray counterpart by SUZAKU. The corresponding filament catalogue will be made available online.

The full moments expansion of the joint probability distribution of an isotropic random field, its gradient and invariants of the Hessian is presented in 2 and 3D. It allows for explicit expression for the Euler characteristic in ND and computation of extrema counts as functions of the excursion set threshold and the spectral parameter, as illustrated on model examples.

The correlation between the large-scale distribution of galaxies and their spectroscopic properties at z=1.5 is investigated using the Horizon MareNostrum cosmological run. We have extracted a large sample of 10^5 galaxies from this large hydrodynamical simulation featuring standard galaxy formation physics. Spectral synthesis is applied to these single stellar populations to generate spectra and colours for all galaxies. We use the skeleton as a tracer of the cosmic web and study how our galaxy catalogue depends on the distance to the skeleton. We show that galaxies closer to the skeleton tend to be redder, but that the effect is mostly due to the proximity of large haloes at the nodes of the skeleton, rather than the filaments themselves. This effects translate into a bimodality in the colour distribution of our sample. The origin of this bimodality is investigated and seems to follow from the ram pressure stripping of satellite galaxies within the more massive clusters of the simulation. The virtual catalogues (spectroscopical properties of the MareNostrum galaxies at various redshifts) are available online at http://www.iap.fr/users/pichon/MareNostrum/catalogues

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

We investigate the problem of probing the local spatial structure of the magnetic field of the interstellar medium using multi-frequency polarized maps of the synchrotron emission at radio wavelengths. We focus in this paper on the three-dimensional reconstruction of the largest scales of the magnetic field, relying on the internal depolarization (due to differential Faraday rotation) of the emitting medium as a function of electromagnetic frequency. We argue that multi-band spectroscopy in the radio wavelengths, developed in the context of high-redshift extragalactic HI lines, can be a very useful probe of the 3D magnetic field structure of our Galaxy when combined with a Maximum A Posteriori reconstruction technique. When starting from a fair approximation of the magnetic field, we are able to recover the true one by using a linearized version of the corresponding inverse problem. The spectral analysis of this problem allows us to specify the best sampling strategy in electromagnetic frequency and predicts a spatially anisotropic distribution of posterior errors. The reconstruction method is illustrated for reference fields extracted from realistic magneto-hydrodynamical simulations.

We report on two quantitative, morphological estimators of the filamentary structure of the Cosmic Web, the so-called global and local skeletons. The first, based on a global study of the matter density gradient flow, allows us to study the connectivity between a density peak and its surroundings, with direct relevance to the anisotropic accretion via cold flows on galactic halos. From the second, based on a local constraint equation involving the derivatives of the field, we can derive predictions for powerful statistics, such as the differential length and the relative saddle to extrema counts of the Cosmic web as a function of density threshold (with application to percolation of structures and connectivity), as well as a theoretical framework to study their cosmic evolution through the onset of gravity-induced non-linearities.

Baryonic Acoustic Oscillations (BAO) and their effects on the matter power spectrum can be studied by using the Lyman-alpha absorption signature of the matter density field along quasar (QSO) lines of sight. A measurement sufficiently accurate to provide useful cosmological constraints requires the observation of ~100000 quasars in the redshift range 2.2<z<3.5 over ~8000 deg2. Such a survey is planned by the Baryon Oscillation Spectroscopic Survey (BOSS) project of the Sloan Digital Sky Survey (SDSS-III).In practice, one needs a stellar rejection of more than two orders of magnitude with a selection efficiency for quasars better than 50% up to magnitudes as large as g ~ 22. To obtain an appropriate target list and estimate quasar redshifts, we have developed an Artificial Neural Networks (NN) with a multilayer perceptron architecture. The input variables are photometric measurements, i.e. the object magnitudes and their errors in the five bands (ugriz) of the SDSS photometry. For target selection, we achieve a non-quasar point-like object rejection of 99.6% and 98.5% for a quasar efficiency of, respectively, 50% and 85%. The photometric redshift precision is of the order of 0.1 over the region relevant for BAO studies.

Within the context of upcoming full-sky lensing surveys, the edge-preserving non- linear algorithm Aski is presented. Using the framework of Maximum A Posteriori inversion, it aims at recovering the full-sky convergence map from surveys with masks. It proceeds in two steps: CCD images of crowded galactic fields are deblurred using automated edge-preserving deconvolution; once the reduced shear is estimated, the convergence map is also inverted via an edge- preserving method. For the deblurring, it is found that when the observed field is crowded, this gain can be quite significant for realistic ground-based surveys when both positivity and edge-preserving penalties are imposed during the iterative deconvolution. For the convergence inversion, the quality of the reconstruction is investigated on noisy maps derived from the horizon N-body simulation, with and without Galactic cuts, and quantified using one-point statistics, power spectra, cluster counts, peak patches and the skeleton. It is found that the reconstruction is able to interpolate and extrapolate within the Galactic cuts/non-uniform noise; its sharpness-preserving penalization avoids strong biasing near the clusters of the map; it reconstructs well the shape of the PDF as traced by its skewness and kurtosis; the geometry and topology of the reconstructed map is close to the initial map as traced by the peak patch distribution and the skeleton's differential length; the two-points statistics of the recovered map is consistent with the corresponding smoothed version of the initial map; the distribution of point sources is also consistent with the corresponding smoothing, with a significant improvement when edge preserving prior is applied. The contamination of B-modes when realistic Galactic cuts are present is also investigated. Leakage mainly occurs on large scales.

[Abridged] This paper aims at providing new conservative constraints to the cosmic star-formation history from the empirical modeling of mid- and far-infrared data. We perform a non-parametric inversion of galaxy counts at 15, 24, 70, 160, and 850 microns simultaneously. It is a "blind" search (no redshift information is required) of all possible evolutions of the infrared luminosity function of galaxies, from which the evolution of the star-formation rate density and its uncertainties are derived. The cosmic infrared background (CIRB) measurements are used a posteriori to tighten the range of solutions. The inversion relies only on two hypotheses: (1) the luminosity function remains smooth both in redshift and luminosity, (2) a set of infrared spectral energy distributions (SEDs) of galaxies must be assumed. The range of star-formation histories that we derive is well constrained and consistent with redshift-based measurements from deep surveys. The redshift decompositions of the counts are also recovered successfully. Therefore, multi-wavelength counts and CIRB (both projected observations) alone seem to contain enough information to recover the cosmic star-formation history with quantifiable errors. A peak of the SFRD at z~2 is preferred, although higher redshifts are not excluded. We also find a good consistency between the observed evolution of the stellar mass density and the prediction from our model. Finally, the inability of the inversion to model perfectly and simultaneously all the multi-wavelength infrared counts (especially at 160 microns where an excess is seen around 20 mJ) implies either (i) the existence of a sub-population of colder galaxies, (ii) a larger dispersion of dust temperatures among local galaxies than expected, (iii) or a redshift evolution of the infrared SEDs of galaxies.

We have performed a 70 billion dark-matter particles N-body simulation in a 2 $h^{-1}$ Gpc periodic box, using the concordance, cosmological model as favored by the latest WMAP3 results. We have computed a full-sky convergence map with a resolution of $\Delta \theta \simeq 0.74$ arcmin$^{2}$, spanning 4 orders of magnitude in angular dynamical range. Using various high-order statistics on a realistic cut sky, we have characterized the transition from the linear to the nonlinear regime at $\ell \simeq 1000$ and shown that realistic galactic masking affects high-order moments only below $\ell < 200$. Each domain (Gaussian and non-Gaussian) spans 2 decades in angular scale. This map is therefore an ideal tool for testing map-making algorithms on the sphere. As a first step in addressing the full map reconstruction problem, we have benchmarked in this paper two denoising methods: 1) Wiener filtering applied to the Spherical Harmonics decomposition of the map and 2) a new method, called MRLens, based on the modification of the Maximum Entropy Method on a Wavelet decomposition. While the latter is optimal on large spatial scales, where the signal is Gaussian, MRLens outperforms the Wiener method on small spatial scales, where the signal is highly non-Gaussian. The simulated full-sky convergence map is freely available to the community to help the development of new map-making algorithms dedicated to the next generation of weak-lensing surveys.

The massive galaxies in the young universe, ten billion years ago, formed stars at surprising intensities. Although this is commonly attributed to violent mergers, the properties of many of these galaxies are incompatible with such events, showing gas-rich, clumpy, extended rotating disks not dominated by spheroids (Genzel et al. 2006, 2008). Cosmological simulations and clustering theory are used to explore how these galaxies acquired their gas. Here we report that they are stream-fed galaxies, formed from steady, narrow, cold gas streams that penetrate the shock-heated media of massive dark matter haloes (Dekel & Birnboim 2006; Keres et al. 2005). A comparison with the observed abundance of star-forming galaxies implies that most of the input gas must rapidly convert to stars. One-third of the stream mass is in gas clumps leading to mergers of mass ratio greater than 1:10, and the rest is in smoother flows. With a merger duy cycle of 0.1, three-quarters of the galaxies forming stars at a given rate are fed by smooth streams. The rarer, submillimetre galaxies that form stars even more intensely are largely merger-induced starbursts. Unlike destructive mergers, the streams are likely to keep the rotating disk configuration intact, although turbulent and broken into giant star-forming clumps that merge into a central spheroid (Noguchi 1999; Genzel et al. 2008, Elmegreen, Bournaud & Elmegreen 2008, Dekel, Sari & Ceverino 2009). This stream-driven scenario for the formation of disks and spheroids is an alternative to the merger picture.

A method to compute the full hierarchy of the critical subsets of a density field is presented. It is based on a watershed technique and uses a probability propagation scheme to improve the quality of the segmentation by circumventing the discreteness of the sampling. It can be applied within spaces of arbitrary dimensions and geometry. This recursive segmentation of space yields, for a $d$-dimensional space, a $d-1$ succession of $n$-dimensional subspaces that fully characterize the topology of the density field. The final 1D manifold of the hierarchy is the fully connected network of the primary critical lines of the field : the skeleton. It corresponds to the subset of lines linking maxima to saddle points, and provides a definition of the filaments that compose the cosmic web as a precise physical object, which makes it possible to compute any of its properties such as its length, curvature, connectivity etc... When the skeleton extraction is applied to initial conditions of cosmological N-body simulations and their present day non linear counterparts, it is shown that the time evolution of the cosmic web, as traced by the skeleton, is well accounted for by the Zel'dovich approximation. Comparing this skeleton to the initial skeleton undergoing the Zel'dovich mapping shows that two effects are competing during the formation of the cosmic web: a general dilation of the larger filaments that is captured by a simple deformation of the skeleton of the initial conditions on the one hand, and the shrinking, fusion and disappearance of the more numerous smaller filaments on the other hand. Other applications of the N dimensional skeleton and its peak patch hierarchy are discussed.

The local theory of the critical lines of 2D and 3D Gaussian fields that underline the cosmic structures is presented. In the context of cosmological matter distribution the subset of critical lines of the 3D density field serves to delineate the skeleton of the observed filamentary structure at large scales. A stiff approximation used to quantitatively describe the filamentary skeleton shows that the flux of the skeleton lines is related to the average Gaussian curvature of the 1D (2D) sections of the field, much in the same way as the density of the peaks. The distribution of the length of the critical lines with threshold is analyzed in detail, while the extended descriptors of the skeleton - its curvature and its singular points, are introduced and briefly described. Theoretical predictions are compared to measurements of the skeleton in realizations of Gaussian random fields in 2D and 3D. It is found that the stiff approximation predicts accurately the shape of the differential length, allows for analytical insight, and explicit closed form solutions. Finally, it provides a simple classification of the singular points of the critical lines: i) critical points; ii) bifurcation points; iii) slopping plateaux.

abridged] A method to rapidly estimate the Fourier power spectrum of a point distribution is presented. This method relies on a Taylor expansion of the trigonometric functions. It yields the Fourier modes from a number of FFTs, which is controlled by the order N of the expansion and by the dimension D of the system. In three dimensions, for the practical value N=3, the number of FFTs required is 20. We apply the method to the measurement of the power spectrum of a periodic point distribution that is a local Poisson realization of an underlying stationary field. We derive explicit analytic expression for the spectrum, which allows us to quantify--and correct for--the biases induced by discreteness and by the truncation of the Taylor expansion, and to bound the unknown effects of aliasing of the power spectrum. We show that these aliasing effects decrease rapidly with the order N. The only remaining significant source of errors is reduced to the unavoidable cosmic/sample variance due to the finite size of the sample. The analytical calculations are successfully checked against a cosmological N-body experiment. We also consider the initial conditions of this simulation, which correspond to a perturbed grid. This allows us to test a case where the local Poisson assumption is incorrect. Even in that extreme situation, the third-order Fourier-Taylor estimator behaves well. We also show how to reach arbitrarily large dynamic range in Fourier space (i.e., high wavenumber), while keeping statistical errors in control, by appropriately "folding" the particle distribution.

This paper reviews and analyses methods used to identify neighbours in 6D space and estimate the corresponding phase-space density. It compares SPH methods to 6D Delaunay tessellation on statical and dynamical realisation of single halo profiles, paying attention to the unknown scaling, S_G, used to relate the spatial dimensions to the velocity dimensions. The methods with local adaptive metric provide the best phase-space estimators. They make use of a Shannon entropy criterion combined with a binary tree partitioning and with SPH interpolation using 10-40 neighbours. Local scaling implemented by such methods, which enforces local isotropy of the distribution function, can vary by about one order of magnitude in different regions within the system. It presents a bimodal distribution, in which one component is dominated by the main part of the halo and the other one is dominated by the substructures. While potentially better than SPH techniques, since it yields an optimal estimate of the local softening volume (and the local number of neighbours required to perform the interpolation), the Delaunay tessellation in fact poorly estimates the phase-space distribution function. Indeed, it requires, the choice of a global scaling S_G. We propose two methods to estimate S_G that yield a good global compromise. However, the Delaunay interpolation still remains quite sensitive to local anisotropies in the distribution. We also compare 6D phase-space density estimation with the proxy, Q=rho/sigma^3, where rho is the local density and 3 sigma^2 is the local 3D velocity dispersion. We show that Q only corresponds to a rough approximation of the true phase-space density, and is not able to capture all the details of the distribution in phase-space, ignoring, in particular, filamentation and tidal streams.

We present a statistical study of the effects induced by substructures on the deflection potential of dark matter halos in the strong lensing regime. This investigation is based on the pertubative solution around the Einstein radius (Alard 2007) in which all the information on the deflection potential is specified by only a pair of one-dimensional functions on this ring. Using direct comparison with ray-tracing solutions, we found that the iso-contours of lensed images predicted by the pertubative solution is reproduced with a mean error on their radial extension of less than 1% - in units of the Einstein radius, for reasonable substructure masses. It demonstrates the efficiency of the approximation to track possible signatures of substructures. We have evaluated these two fields and studied their properties for different lens configurations modelled either through massive dark matter halos from a cosmological N-body simulation, or via toy models of Monte Carlo distribution of substructures embedded in a triaxial Hernquist potential. As expected, the angular power spectra of these two fields tend to have larger values for larger harmonic numbers when substructures are accounted for and they can be approximated by power-laws, whose values are fitted as a function of the profile and the distribution of the substructures.

(abridged) The physics of diffuse gas accretion and the properties of the cold and hot modes of accretion onto proto-galaxies between z=2 and z=5.4 is investigated using the large cosmological simulation performed with the RAMSES code on the MareNostrum supercomputing facility. Galactic winds, chemical enrichment, UV background heating and radiative cooling are taken into account in this very high resolution simulation. Using accretion--weighted temperature histograms, we have perfomed two different measurements of the thermal state of the gas accreted towards the central galaxy. The first measurement, performed using accretion--weighted histograms on a spherical surface of radius 0.2 Rvir centred on the densest gas structure near the halo centre of mass, is a good indicator of the presence of an accretion shock in the vicinity of the galactic disc. We define the hot shock mass, Mshock, as the typical halo mass separating cold dominated from hot dominated accretion in the vicinity of the galaxy. The second measurement is performed by radially averaging histograms between 0.2 Rvir and Rvir, in order to detect radially extended structures such as gas filaments: this is a good proxy for detecting cold streams feeding the central galaxy. We define Mstream as the transition mass separating cold dominated from hot dominated accretion in the outer halo, marking the disappearance of these cold streams. We find a hot shock transition mass of Mshock=10^{11.6} Msun (dark matter), with no significant evolution with redshift. Conversely, we find that Mstream increases sharply with z. Our measurements are in agreement with the analytical predictions of Birnboim & Dekel (2003) and Dekel & Birnboim (2006), if we correct their model by assuming low metallicity (< 10^{-3} Zsun) for the filaments, correspondingly to our measurements.

This technical paper describes a software package that was designed to produce initial conditions for large cosmological simulations in the context of the Horizon collaboration. These tools generalize E. Bertschinger's Grafic1 software to distributed parallel architectures and offer a flexible alternative to the Grafic2 software for ``zoom'' initial conditions, at the price of large cumulated cpu and memory usage. The codes have been validated up to resolutions of 4096^3 and were used to generate the initial conditions of large hydrodynamical and dark matter simulations. They also provide means to generate constrained realisations for the purpose of generating initial conditions compatible with, e.g. the local group, or the SDSS catalog.

The onset of stochasticity is measured in $\Lambda$CDM cosmological simulations using a set of classical observables. It is quantified as the local derivative of the logarithm of the dispersion of a given observable (within a set of different simulations differing weakly through their initial realization), with respect to the cosmic growth factor. In an Eulerian framework, it is shown here that chaos appears at small scales, where dynamic is non-linear, while it vanishes at larger scales, allowing the computation of a critical transition scale corresponding to ~ 3.5 Mpc/h. This picture is confirmed by Lagrangian measurements which show that the distribution of substructures within clusters is partially sensitive to initial conditions, with a critical mass upper bound scaling roughly like the perturbation's amplitude to the power 0.15. The corresponding characteristic mass, $M_{\rm crit}=2 10^{13} M_{\odot}$, is roughly of the order of the critical mass of non linearities at z=1 and accounts for the decoupling induced by the dark energy triggered acceleration. The sensitivity to detailed initial conditions spills to some of the overall physical properties of the host halo (spin and velocity dispersion tensor orientation) while other "global" properties are quite robust and show no chaos (mass, spin parameter, connexity and center of mass position). This apparent discrepancy may reflect the fact that quantities which are integrals over particles rapidly average out details of difference in orbits, while the other observables are more sensitive to the detailed environment of forming halos and reflect the non-linear scale coupling characterizing the environments of halos.

We investigate how well the 3D density field of neutral hydrogen in the Intergalactic Medium (IGM) can be reconstructed using the Lyman-alpha absorptions observed along lines of sight to quasars separated by arcmin distances in projection on the sky. We use cosmological hydrodynamical simulations to compare the topologies of different fields: dark matter, gas and neutral hydrogen optical depth and to investigate how well the topology of the IGM can be recovered from the Wiener interpolation method implemented by Pichon et al. (2001). The global statistical and topological properties of the recovered field are analyzed quantitatively through the power-spectrum, the probability distribution function (PDF), the Euler characteristics, its associated critical point counts and the filling factor of underdense regions. The local geometrical properties of the field are analysed using the local skeleton by defining the concept of inter-skeleton distance. At scales larger than ~1.4 <d_LOS>, where <d_LOS> is the mean separation between lines of sight, the reconstruction accurately recovers the topological features of the large scale density distribution of the gas, in particular the filamentary structures. At scales larger than the intrinsic smoothing length of the inversion procedure, the power spectrum of the recovered HI density field matches well that of the original one and the low order moments of the PDF are well recovered as well as the shape of the Euler characteristic. The integral errors on the PDF and the critical point counts are indeed small, less than 20% for <d_LOS>~2.5 arcmin. The small deviations between the reconstruction and the exact solution mainly reflect departures from the log-normal behaviour that are ascribed to highly non-linear objects in overdense regions.

The skeleton formalism aims at extracting and quantifying the filamentary structure of the universe is generalized to 3D density fields; a numerical method for computating a local approximation of the skeleton is presented and validated here on Gaussian random fields. This method manages to trace well the filamentary structure in 3D fields such as given by numerical simulations of the dark matter distribution on large scales and is insensitive to monotonic biasing. Two of its characteristics, namely its length and differential length, are analyzed for Gaussian random fields. Its differential length per unit normalized density contrast scales like the PDF of the underlying density contrast times the total length times a quadratic Edgeworth correction involving the square of the spectral parameter. The total length scales like the inverse square smoothing length, with a scaling factor given by 0.21 (5.28+ n) where n is the power index of the underlying field. This dependency implies that the total length can be used to constrain the shape of the underlying power spectrum, hence the cosmology. Possible applications of the skeleton to galaxy formation and cosmology are discussed. As an illustration, the orientation of the spin of dark halos and the orientation of the flow near the skeleton is computed for dark matter simulations. The flow is laminar along the filaments, while spins of dark halos within 500 kpc of the skeleton are preferentially orthogonal to the direction of the flow at a level of 25%.

Systems whose potential energies consists of pieces that scale as r^-2 together with pieces that scale as r^2, show no violent relaxation to Virial equilibrium but may pulsate at considerable amplitude for ever. Despite this pulsation these systems form lattices when the non-pulsational `energy' is low, and these disintegrate as that energy is increased. The `specific heats' show the expected halving as the `solid' is gradually replaced by the `fluid' of independent particles. The forms of the lattices are described here for N ~ 20 and they become hexagonal close packed for large N. In the larger N limit, a shell structure is formed. Their large N behaviour is analogous to a gamma=5/3 polytropic fluid with a quasi-gravity such that every element of fluid attracts every other in proportion to their separation. For such a fluid, we study the `rotating pulsating equilibria' and their relaxation back to uniform but pulsating rotation. We also compare the rotating pulsating fluid to its discrete counter part, and study the rate at which the rotating crystal redistributes angular momentum and mixes as a function of extra heat content.

In a serie of three papers, the dynamical interplay between environments and dark matter haloes is investigated, while focussing on the dynamical flows through their virial sphere. Our method relies on both cosmological simulations, to constrain the environments, and an extension to the classical matrix method to derive the response of the halo (see Pichon & Aubert (2006), paper I). The current paper focuses on the statistical characterisation of the environments surrounding haloes, using a set of large scale simulations. Our description relies on a `fluid' halocentric representation where the interactions between the halo and its environment are investigated in terms of a time dependent external tidal field and a source term characterizing the infall. The method is applied to 15000 haloes, with masses between 5 x 10^12 Ms and 10^14 Ms evolving between z = 1 and z = 0. The net accretion at the virial radius is found to decrease with time, resulting from both an absolute decrease of infall and from a growing contribution of outflows. Infall is found to be mainly radial and occurring at velocities ~ 0.75 V200. Outflows are also detected through the virial sphere and occur at lower velocities ~ 0.6 V200 on more circular orbits. The external tidal field is found to be strongly quadrupolar and mostly stationnary, possibly reflecting the distribution of matter in the halo's near environment. The coherence time of the small scale fluctuations of the potential hints a possible anisotropic distribution of accreted satellites. The flux density of mass on the virial sphere appears to be more clustered than the potential while the shape of its angular power spectrum seems stationnary.

We present the transverse flux correlation function of the Lyman-alpha forest in quasar absorption spectra at z~2.1 from VLT-FORS and VLT-UVES observations of a total of 32 pairs of quasars; 26 pairs with separations in the range 0.6<theta<4 arcmin and 6 pairs with 4<theta<10 arcmin. Correlation is detected at the 3sigma level up to separations of the order of ~4 arcmin (or ~4.4 h^{-1} Mpc comoving at z=2.1 for Omega_{m}=0.3 and Omega_{Lambda}=0.7). We have, furthermore, measured the longitudinal correlation function at a somewhat higher mean redshift (z=2.39) from 20 lines of sight observed with high spectral resolution and high signal-to-noise ratio with VLT-UVES. We compare the observed transverse and longitudinal correlation functions to that obtained from numerical simulations and illustrate the effect of spectral resolution, thermal broadening and peculiar motions. The shape and correlation length of the correlation functions are in good agreement with those expected from absorption by the filamentary and sheet-like structures in the photoionized warm intergalactic medium predicted in CDM-like models for structures formation. Using a sample of 139 C_{iv} systems detected along the lines of sight toward the pairs of quasars we also investigate the transverse correlation of metals on the same scales. The observed transverse correlation function of intervening C_{iv} absorption systems is consistent with that of a randomly distributed population of absorbers. This is likely due to the small number of pairs with separation less than 2 arcmin. We detect, however, a significant overdensity of systems in the sightlines towards the quartet Q0103-294A&B, Q0102-2931 and Q0102-293 which extends over the redshift range 1.5<z<2.2 and an angular scale larger than 10 arcmin.

The length of the three-dimensional filaments observed in the fourth public data-release of the SDSS is measured using the {\em local skeleton} method. It consists in defining the set of points where the gradient of the smoothed density field is extremal along its isocontours, with some additional constraints on local curvature to probe actual ridges in the galaxy distribution. A good fit to the mean filament length per unit volume, $\cal{L}$, in the SDSS survey is found to be ${\cal{L}}=(52500\pm6500) (L/{\rm Mpc})^{-1.75\pm0.06}\rm{Mpc}/(100 \rm{Mpc})^{3}$ for $8.2 \leq L \leq 16.4$ Mpc, where $L$ is the smoothing length in Mpc. This result, which deviates only slightly, as expected, from the trivial behavior ${\cal{L}} \propto L^{-2}$, is in excellent agreement with a $\Lambda$CDM cosmology, as long as the matter density parameter remains in the range $0.25 < \Omega_{\rm matter} < 0.4$ at one sigma confidence level, considering the universe is flat. These measurements, which are in fact dominated by linear dynamics, are not significantly sensitive to observational biases such as redshift distortion, edge effects, incompleteness, and biasing between the galaxy distribution and the dark matter distribution. Hence it is argued that the local skeleton is a rather promising and discriminating tool for the analysis of filamentary structures in three-dimensional galaxy surveys.

We investigate statistically the dynamical consequences of cosmological fluxes of matter and related moments on progenitors of today's dark matter haloes. Their dynamics is described via canonical perturbation theory which accounts for two types of perturbations: the tidal field corresponding to fly-bys and accretion of dark matter through the halo's outer boundary. he dynamical equations are solved linearly, order by order, projecting on a biorthogonal basis to consistently satisfy the field equation. Since our solution of the Boltzmann Poisson equations is explicit, it allows statistical predictions for the ensemble distribution of the inner dynamical features of haloes. The secular evolution of open galactic haloes is investigated: we derive the kinetic equation which governs the quasi-linear evolution of dark matter profile induced by infall and its corresponding gravitational correlations. This yields a Fokker Planck-like equation for the angle-averaged underlying distribution function. We show how these extensions to the classical theory could be used to (i) observationally constrain the statistical nature of the infall (ii) predict the observed distribution and correlations of substructures in upcoming surveys, (iii) predict the past evolution of the observed distribution of clumps, and finally (iv) weight the relative importance of the intrinsic (via the unperturbed distribution function) and external (tidal and/ or infall) influence of the environment in determining the fate of galaxies.

We study the clustering properties of metals in the intergalactic medium (IGM) as traced by 619 CIV, 81 SiIV, N >= 10^12 cm^-2 and 316 MgII, and 82 FeII N >= 10^11.5 cm^-2 absorption components in 19 high signal-to-noise (60-100 per pixel), high resolution (R = 45000) quasar spectra. Over the redshift range probed (1.5-3.0), CIV and SiIV trace each other closely and their line-of-sight correlation functions exhibit a steep decline at large separations and a flatter profile below ~ 150 km s^-1, with a large overall bias. These features do not depend on column depth. Carrying out a detailed SPH simulation (2 X 320^3, 57 Mpc^3 comoving), we show that this behavior can not be reproduced by models in which the IGM metallicity is constant or a local function of density. However, the CIV correlation function is consistent with a model in which metals are confined within bubbles with a typical radius Rs = 2 comoving Mpc about sources of mass >= Ms = 10^12 solar masses at z=3. Our lower redshift (0.5-2) measurements of the MgII and FeII correlation functions also uncover a steep decline at large separations and a flatter profile at small separations, but the clustering is even higher, and the turn-over is shifted to ~ 75 km s^-1. Again these features do not change with column depth. We describe an analytical bubble model for these species, which come from regions that are too compact to be simulated numerically, deriving best-fit values of R_s ~ 2.4 Mpc and M_s ~ 10^12 solar masses. Equally good fits to all four species are found in a similarly biased high-redshift enrichment model in which metals are placed within 2.4 comoving Mpc of 3 x 10^9 solar mass sources at z = 7.5.

We derive the transverse flux correlation function in the Lyman-alpha forest at z~2.1 from VLT-FORS observations of a total of 32 pairs of quasars. The shape and correlation length of the transverse correlation function are in good agreement with the paradigm of intergalactic medium predicted in CDM-like models for structures formation. Using a sample of 139 Civ systems detected along the lines of sight toward the pairs of quasars we investigate the transverse correlation of metals on the same scales. We find that the correlation function is consistent with that of a randomly distributed population of Civ systems. However, we detect an important overdensity of systems in front of a quartet.

We introduce STECKMAP (STEllar Content and Kinematics via Maximum A Posteriori), a method to recover the kinematical properties of a galaxy simultaneously with its stellar content from integrated light spectra. It is an extension of STECMAP (astro-ph/0505209) to the general case where the velocity distribution of the underlying stars is also unknown. %and can be used as is for the analysis of large sets of data. The reconstructions of the stellar age distribution, the age-metallicity relation, and the Line-Of-Sight Velocity Distribution (LOSVD) are all non-parametric, i.e. no specific shape is assumed. The only a propri we use are positivity and the requirement that the solution is smooth enough. The smoothness parameter can be set by GCV according to the level of noise in the data in order to avoid overinterpretation. We use single stellar populations (SSP) from PEGASE-HR (R=10000, lambda lambda = 4000-6800 Angstrom, Le Borgne et al. 2004) to test the method through realistic simulations. Non-Gaussianities in LOSVDs are reliably recovered with SNR as low as 20 per 0.2 Angstrom pixel. It turns out that the recovery of the stellar content is not degraded by the simultaneous recovery of the kinematic distribution, so that the resolution in age and error estimates given in Ocvirk et al. 2005 remain appropriate when used with STECKMAP. We also explore the case of age-dependent kinematics (i.e. when each stellar component has its own LOSVD). We separate the bulge and disk components of an idealized simplified spiral galaxy in integrated light from high quality pseudo data (SNR=100 per pixel, R=10000), and constrain the kinematics (mean projected velocity, projected velocity dispersion) and age of both components.

This paper describes STECMAP (STEllar Content via Maximum A Posteriori), a flexible, non-parametric inversion method for the interpretation of the integrated light spectra of galaxies, based on synthetic spectra of single stellar populations (SSPs). We focus on the recovery of a galaxy's star formation history and stellar age-metallicity relation. We use the high resolution SSPs produced by PEGASE-HR to quantify the informational content of the wavelength range 4000 - 6800 Angstroms. A detailed investigation of the properties of the corresponding simplified linear problem is performed using singular value decomposition. It turns out to be a powerful tool for explaining and predicting the behaviour of the inversion. We provide means of quantifying the fundamental limitations of the problem considering the intrinsic properties of the SSPs in the spectral range of interest, as well as the noise in these models and in the data. We performed a systematic simulation campaign and found that, when the time elapsed between two bursts of star formation is larger than 0.8 dex, the properties of each episode can be constrained with a precision of 0.04 dex in age and 0.02 dex in metallicity from high quality data (R=10 000, signal-to-noise ratio SNR=100 per pixel), not taking model errors into account. The described methods and error estimates will be useful in the design and in the analysis of extragalactic spectroscopic surveys.

One way to extract more information from the integrated light of galaxies is to improve the spectral resolution at which observations and analysis are carried out. The population synthesis code currently providing the highest spectral resolution is Pegase-HR, which was made available by D. Le Borgne et al. in 2004. Based on an empirical stellar library, it provides synthetic spectra between 4000 and 6800 A at lambda/d(lambda)=10000 for any star formation history, with or without chemical evolution. Such a resolution is particularly useful for the study of low mass galaxies, massive star clusters, or other galaxy regions with low internal velocity dispersions. After a summary of the main features of Pegase-HR and comparisons with other population synthesis codes, this paper focuses on the inversion of optical galaxy spectra. We explore the limits of what information can or can not be recovered, based on theoretical principles and extensive simulations. First applications to extragalactic objects are shown.

We measure the anisotropy of dark matter flows on small scales in the near environment of haloes using a large set of simulations. We rely on two different approaches to quantify the anisotropy of the cosmic infall: we measure the flows at the haloes' virial radius while describing the infalling matter via fluxes through a spherical shell; we measure the spatial and kinematical distributions of satellites and substructures around haloes detected by the subclump finder ADAPTAHOP first described in the appendix. The two methods are found to be in agreement both qualitatively and quantitatively via one and two points statistics.The infall takes place preferentially in the plane perpendicular to the direction defined by the halo's spin. We computed the excess of equatorial accretion both through rings and via a harmonic expansion of the infall. The level of anisotropy of infalling matter is found to be ~15 %. The substructures have their spin orthogonal to their velocity vector in the halo's rest frame at a level of about 5%, suggestive of an image of a flow along filamentary structures which provides an explanation for the measured anisotropy. We conclude that a halo does not see its environment as an isotropic perturbation, investigate how the anisotropy is propagated inwards using perturbation theory, and discuss briefly implications for weak lensing, warps and the thickness of galactic disks.

We use high spectral resolution (R=45000) and high signal-to-noise ratio (S/N \~ 35-70 per pixel) spectra of 19 high-redshift (2.1<z_em<3.2) quasars to investigate the metal content of the low-density intergalactic medium using pixel-by-pixel procedures. This high quality homogeneous survey gives the possibility to statistically search for metals at HI optical depths smaller than unity. We find that the gas is enriched in carbon and oxygen for neutral hydrogen optical depths tau(HI)>1. Our observations strongly suggest that the CIV/HI ratio decreases with decreasing tau(HI) with log[tau(CIV)]=1.3xlog[tau(HI)]-3.2. We do not detect CIV absorption statistically associated with gas of tau(HI)<1. However, we observe that a small fraction of the low density gas is associated with strong metal lines as a probable consequence of the IGM enrichment being highly inhomogeneous. We detect the presence of OVI down to tau(HI)~0.2 with log[tau(OVI)/tau(HI)] -2.0. We show that OVI absorption in the lowest density gas is located within ~300 km/s from strong HI lines. This suggests that this OVI phase may be part of winds flowing away from overdense regions. This effect is more important at the largest redshifts (z>2.4). Therefore, at the limit of present surveys, the presence of metals in the underdense regions of the IGM is still to be demonstrated.

We present inversion techniques aiming at recovering the stellar content (age, metallicity, extinction) and kinematics of a stellar population from its absorption line spectrum. These techniques use new synthetic high resolution spectra produced by PEGASE_HR and minimization algorithms. Application to SDSS data and to a simulated disc-bulge system are presented as examples.

We measure the spatial clustering of metals in the intergalactic medium from z = 1.7 to 3.0, as traced by 643 CIV and 104 SiIV N >= 10^12 cm^-2 absorption systems in 19 high signal-to-noise (40-80) and high resolution (R = 45000) quasar spectra. The number densities and two-point correlation functions of both these species are largely constant with redshift, suggesting the bulk of metal ejection occurred at z >= 3. However, at z <= 1.9 some additional signature appears in the CIV correlation function at 500 km/s, associated with four strong and peculiar systems. At all redshifts, the CIV and SiIV correlation functions exhibit a steep rise at large separations and a flatter profile at small separations, with an elbow occurring at ~150 km/s. We show that these properties are consistent with metals confined within bubbles with a typical radius R_s about sources of mass >= M_s, and use numerical simulations to derive best-fit values of R_s ~ 2 comoving Mpc and M_s ~ 5x10^11 solar masses at z = 3. This does not exclude that metals could have been produced at higher redshifts in smaller, but equally rare, objects. At the level of detection of this survey, IGM enrichment is likely to be incomplete and inhomogeneous, with a filling factor ~ 10%.

We study the clustering of matter in the intergalactic medium from the Lyman-alpha forests seen in the spectra of pairs or groups of z ~ 2 quasars observed with FORS2 and UVES at the VLT-UT2 Kueyen ESO telescope. The sample consists of five pairs with separations 0.6, 1, 2.1, 2.6 and 4.4 arcmin and a group of four quasars with separations from 2 up to 10 arcmin. This unprecedented data set allows us to measure the transverse flux correlation function for a range of angular scales. Correlations are clearly detectable at separations smaller than 3 arcmin. The shape and correlation length of the transverse correlation function on these scales is in good agreement with those expected from absorption by the photoionized warm intergalactic medium associated with the filamentary and sheet-like structures predicted in CDM-like models for structure formation. At larger separation no significant correlation is detected. Assuming that the absorbing structures are randomly orientated with respect to the line of sight, the comparison of transverse and longitudinal correlation lengths constrains the cosmological parameters (as a modified version of the Alcock & Paczy\'nski test). The present sample is too small to get significant constraints. Using N-body simulations, we investigate the possibility to constrain Ol from future larger samples of QSO pairs with similar separations. The observation of a sample of 30 pairs at 2, 4.5 and 7.5 arcmin should constrain the value of Ol at +- 15 % (2 sigma level). We further use the observed spectra of the group of four quasars, to search for underdense regions in the intergalactic medium. We find a quasi-spherical structure of reduced absorption with radius 12.5 h^{-1} Mpc which we identify as an underdense region.

Results are presented for our search for warm-hot gas towards the quasar Q 0329-385. We identify ten OVI systems of which two are within 5000 km/s of z_em and a third one should be of intrinsic origin. The seven remaining systems have HI column densities 10^13.7<N(HI)<10^15.6 cm^-2. At least ~1/3 of the individual OVI sub-systems have temperatures T<1x10^5 K and cannot originate in collisionally ionized gas. Photoionization by a hard UV background field reproduces well the ionic ratios for metallicities in the range 10^-2.5 - 10^-0.5 solar, with possibly sub-solar N/C relative abundance. For [O/C]=0, the sizes inferred for the OVI clouds are in some cases larger than the maximum extent implied by the Hubble flow. This constraint is fulfilled assuming a moderate overabundance of oxygen relative to carbon. For a soft UV ionizing spectrum, an overabundance of O/C is required, [O/C]~0.0-1.3. For a hard(soft) UV spectrum and [O/C]=0(1), the OVI regions have overdensities in 10-40.

A method for inverting the statistical star counts equation, including proper motions, is presented; in order to break the degeneracy in that equation it uses the supplementary constraints required by dynamical consistency. The inversion gives access to both the kinematics and the luminosity function of each population in three r\'egimes: the singular ellipsoid, the constant ratio Schwarzschild ellipsoid plane parallel models and the epicyclic model. This more realistic model is taylored to account for local neighbourhood density and velocity distribution. The first model is fully investigated both analytically and via means of a non-parametric inversion technique, while the second model is shown to be formally its equivalent. The effect of noise and incompleteness in apparent magnitude is investigated. The third model is investigated via a 5D+2D non-parametric inversion technique where positivity of the underlying luminosity function is explicitely accounted for. It is argued that its future application to data such as the Tycho catalogue (and in the upcoming satellite GAIA) could lead -- provided the vertical potential, and/or the asymmetric drift or w_0 are known -- to a non-parametric determination of the local neighbourhood luminosity function without any reference to stellar evolution tracks. It should also yield the proportion of stars for each kinematic component and a kinematic diagnostic to split the thin disk from the thick disk or the halo.

Very recently a new inversion method has been developped to analyze the intergalactic medium seen in absorption in quasar spectra (the so-called Lyman-alpha forest). This method is applied to recover the temperature of the gas and the underlying density field. Using constraints from the Lyman-beta forest, it is possible to recover this field up to over-densities delta=10. By inverting the HI and CIV absorptions together it has been shown that the CIV/HI ratio varies through the profile of strong lines, beeing larger in the wings. The method can be applied to reconstruct the 3D density field from multiple lines of sight and is shown to give good results up to mean separations of 3 arcmin. Results from a survey of QSO pairs performed with HST/STIS and VLT/UVES-FORS are summarized.

We study the physical properties of the Lyman-alpha forest by applying the inversion method described by Pichon et al. (2001) to the high resolution and high S/N ratio spectrum of the z_em=2.40 quasar HE 1122-1628 obtained during Science Verification of UVES at the VLT. We compare the column densities obtained with the new fitting procedure with those derived using standard Voigt profile methods. The agreement is good and gives confidence in the new description of the Lyman-alpha forest as a continuous field as derived from our method. We show that the observed number density of lines with logN>13 and 14 is, respectively, 50 and 250 per unit redshift at z~2. We study the physical state of the gas, neglecting peculiar velocities, assuming a relation between the overdensity and the temperature. T=Tbar * rho^(2beta). There is an intrinsic degeneracy between the parameters beta and Tbar. We demonstrate that, at a fixed beta, the temperature at mean density, Tbar, can be uniquely extracted however. While applying the method to HE 1122-1628, we conclude that for 0.2<beta<0.3, 6000<Tbar<15000 K at z~2. We investigate the small scale structure of strong absorption lines using the information derived from the Lyman-beta, Lyman-gamma and Civ profiles. Introducing the Lyman-beta line in the fit allows us to reconstruct the density field up to rho~10 instead of 5 for the Lyman-alpha line only. There may be small velocity shifts ~10km.s^{-1} between the peaks in the Civ and Hi density profiles. Although the statistics is small, it seems that Civ/Hi and n_HI are anti-correlated. This could be a consequence of the high sensitivity of the Civ/Hi ratio to temperature. The presence of associated Ovi absorption, with similar profile, confirms that the gas is photo-ionized and at a temperature of T~10^5 K.

The abundance of high-redshift galaxy clusters depends sensitively on the matter density $\OmM$ and, to a lesser extent, on the cosmological constant $\Lambda$. Measurements of this abundance therefore constrain these fundamental cosmological parameters, and in a manner independent and complementary to other methods, such as observations of the cosmic microwave background and distance measurements. Cluster abundance is best measured by the X-ray temperature function, as opposed to luminosity, because temperature and mass are tightly correlated, as demonstrated by numerical simulations. Taking advantage of the sensitivity of XMM-Newton, our Guaranteed Time program aims at measuring the temperature of the highest redshift (z>0.4) SHARC clusters, with the ultimate goal of constraining both $\OmM$ and $\Lambda$.

We discuss the implementation of Bayesian inversion methods in order to recover the properties of the intergalactic medium from observations of the neutral hydrogen Lyman-$\alpha$ absorptions observed in the spectra of high redshift quasars (the so-called Lyman-$\alpha$ forest). We use two complementary schemes (i) a constrained Gaussian random field linear approach and (ii) a more general non-linear explicit Bayesian deconvolution method which offers in particular the possibility to constrain the parameters of the equation of state for the gas. While relying on prior assumption for the two-point correlation functions, we show how to recover, at least qualitatively, the 3D topology of the large scale structures in redshift space by inverting a suitable network of adjacent, low resolution lines of sight. We also discuss the inversion of single lines of sight observed at high spectral resolution. Our preliminary investigations suggest that the explicit Bayesian method can be used to derive quantitative information on the physical state of the gas. Redshift distortion is considered by simultaneous constrained reconstruction of the velocity and the density field in real space while assuming statistical correlation between the two fields. The method seems to work well in the strong prior regime where peculiar velocities are assumed to be the most likely realization in the density field. Finally, we investigate the effect of line of sight separation and number of lines of sight. Our analyses suggest that multiple low resolution lines of sight could be used to improve most likely velocity reconstruction on a high resolution line of sight.

We report a new analysis of stellar dynamics in the Galactic centre, based on improved sky and LOS velocities for >100 stars within a few arcsec of SgrA*. Overall the motions do not deviate strongly from isotropy. For 32 stars with all 3 components determined the absolute, LOS and sky velocities agree well, as for a spherical cluster. The projected radial and tangential motions of all 104 pm stars are also consistent with overall isotropy. However, the projected velocities of the young, early type stars have a strong radial dependence. Most of the bright HeI stars 1-10" from SgrA* are on tangential orbits. This anisotropy of the HeI stars and most of the brighter IRS16 complex members is largely caused by a CW and counter-rotating, coherent rotation pattern. The overall rotation of the young star cluster probably is a remnant of the angular momentum in the cloud these stars formed from. The fainter, fast stars within \~1" of SgrA* appear to be largely moving on radial or very elliptical orbits. We have not detected nonlinear motion for any of them. Most of the SgrA* cluster members also are on CW orbits. Spectroscopy shows them to be early type stars. We propose that the SgrA* cluster stars are those members of the early type cluster with small angular momentum which can plunge to the vicinity of SgrA*. Our anisotropy-independent estimate of the Sun-GC distance is 7.8-8.2 kpc (+/- 0.9 kpc). We include velocity anisotropy in estimating the central mass distribution. We confirm previous conclusions that a compact central mass concentration is present and dominates the potential from 0.01-1 pc. The central mass ranges from 2.6-3.3E6 M_sun. (abridged)

A fundamental hypothesis for the interpretation of the measured large-scale line-of-sight peculiar velocities of galaxies is that the large-scale cosmic flows are irrotational. In order to assess the validity of this assumption, we estimate, within the frame of the gravitational instability scenario, the amount of vorticity generated after the first shell crossings in large-scale caustics. In the Zel'dovich approximation the first emerging singularities form sheet like structures. Here we compute the expectation profile of an initial overdensity under the constraint that it goes through its first shell crossing at the present time. We find that this profile corresponds to rather oblate structures in Lagrangian space. Assuming the Zel'dovich approximation is still adequate not only at the first stages of the evolution but also slightly after the first shell crossing, we calculate the size and shape of those caustics and their vorticity content as a function of time and for different cosmologies. The average vorticity created in these caustics is small: of the order of one (in units of the Hubble constant). To illustrate this point we compute the contribution of such caustics to the probability distribution function of the filtered vorticity at large scales. We find that this contribution that this yields a negligible contribution at the 10 to 15 $h^{-1}$Mpc scales. It becomes significant only at the scales of 3 to 4 $h^{-1}$Mpc, that is, slightly above the galaxy cluster scales.

This paper is the fourth of a series dealing with the cluster of galaxies ABCG 85. Using our two extensive photometric and spectroscopic catalogues (with 4232 and 551 galaxies respectively), we discuss here three topics derived from optical data. First, we present the properties of emission line versus non-emission line galaxies, showing that their spatial distributions somewhat differ; emission line galaxies tend to be more concentrated in the south region where groups appear to be falling onto the main cluster, in agreement with the hypothesis (presented in our previous paper) that this infall may create a shock which can heat the X-ray emitting gas and also enhance star formation in galaxies. Then, we analyze the luminosity function in the R band, which shows the presence of a dip similar to that observed in other clusters at comparable absolute magnitudes; this result is interpreted as due to comparable distributions of spirals, ellipticals and dwarfs in these various clusters. Finally, we present the dynamical analysis of the cluster using parametric and non-parametric methods and compare the dynamical mass profiles obtained from the X-ray and optical data.

A general inversion technique for the recovery of the underlying distribution function for observed galactic disks is presented and illustrated. Under the assumption that these disks are axi-symmetric and thin, the proposed method yields the unique distribution compatible with all the observables available. The derivation may be carried out from the measurement of the azimuthal velocity distribution arising from positioning the slit of a spectrograph along the major axis of the galaxy. More generally, it may account for the simultaneous measurements of velocity distributions corresponding to slits presenting arbitrary orientations with respect to the major axis. The approach is non-parametric, i.e. it does not rely on a particular algebraic model for the distribution function. Special care is taken to account for the fraction of counter-rotating stars which strongly affects the stability of the disk. An optimisation algorithm is devised -- generalising the work of Skilling & Bryan (1984) -- to carry this truly two-dimensional ill-conditioned inversion efficiently. The performances of the overall inversion technique with respect to the noise level and truncation in the data set is investigated with simulated data. Reliable results are obtained up to a mean signal to noise ratio of~5 and when measurements are available up to $4 R_{e}$. A discussion of the residual biases involved in non parametric inversions is presented. Prospects of application to observed galaxies and other inversion problems are discussed.

Hipparcos data provide the first, volume limited and absolute magnitude limited homogeneous tracer of stellar density and velocity distributions in the solar neighbourhood. The density of A-type stars more luminous than $M_v=2.5$ can be accurately mapped within a sphere of 125 pc radius, while proper motions in galactic latitude provide the vertical velocity distribution near the galactic plane. The potential well across the galactic plane is traced practically hypothesis-free and model-free. The local dynamical density comes out as $\rho_{0}=0.076 \pm0.015~M_{\sun}~{pc}^{-3}$ a value well below all previous determinations leaving no room for any disk shaped component of dark matter.

In a recent Letter, Cowsik, Ratnam and Bhattacharjee (1996a) have built a dynamically self-consistent spatial distribution of particles of galactic dark matter. They have come up with the rather unorthodox conclusion that the mean velocity dispersion of dark matter particles ``should be $600 kms$ or larger''. Their letter triggered immediate comments (Evans 1996, Gates et al. 1996). Here we find, as did Cowsik et al. (1996a), that models of the dark matter halo can be made consistent with velocity dispersion much larger than that expected from a simple application of the virial theorem in the solar neighbourhood. But in contrast to their conclusions, we show that using their model, we also obtain solutions with smaller velocity dispersion down to $\sim 270 kms$. These more orthodox dispersions arise because of less constraining boundary conditions for the central density but do not rely on indirect or model dependent measurements of the large scale distance behaviour of the rotation curve.

The method originally developed by Kalnajs for the numerical linear stability analysis of round galactic disks is implemented in the regimes of non-analytic transformations between position space and angle-action space, and of vanishing growth rates. This allows effectively any physically plausible disk to be studied, rather than only those having analytic transformations into angle-action space which have formed the primary focus of attention to date. The transformations are constructed numerically using orbit integrations in real space, and the projections of orbit radial actions on a given potential density basis are Fourier transformed to obtain a dispersion relation in matrix form. To verify the implementation, the fastest $m=2$ growth rates of the isochrone and the Kuzmin-Toomre disks are recovered, and the weaker $m=2$ modes are computed. The evolution of those growth rates as a function of the halo mass is also calculated, and some $m=1$ modes are derived as illustration.

A general method is presented for constructing distribution functions for flat systems whose surface density and Toomre's Q number profile is given. The purpose of these functions is to provide plausible galactic models and assess their critical stability with respect to global non axi-symmetric modes. The derivation may be carried out for an azimuthal velocity distribution (or a given specific energy distribution) which may either be observed or chosen to match a specified temperature profile. Distribution functions describing stable models with realistic velocity distributions for power law disks, the Isochrone and the Kuzmin disks are provided. Specially simple inversion formulae are also given for finding distribution functions for flat systems whose surface densities are known.

Complete sequences of new analytic solutions of Einstein's equations which describe thin super massive disks are constructed. These solutions are derived geometrically. The identification of points across two symmetrical cuts through a vacuum solution of Einstein's equations defines the gradient discontinuity from which the properties of the disk can be deduced. The subset of possible cuts which lead to physical solutions is presented. At large distances, all these disks become Newtonian, but in their central regions they exhibit relativistic features such as velocities close that of light, and large redshifts. Sections with zero extrinsic curvature yield cold disks. Curved sections may induce disks which are stable against radial instability. The general counter rotating flat disk with planar pressure tensor is found. Owing to gravomagnetic forces, there is no systematic method of constructing vacuum stationary fields for which the non-diagonal component of the metric is a free parameter. However, all static vacuum solutions may be extended to fully stationary fields via simple algebraic transformations. Such disks can generate a great variety of different metrics including Kerr's metric with any ratio of a to m. A simple inversion formula is given which yields all distribution functions compatible with the characteristics of the flow, providing formally a complete description of the stellar dynamics of flattened relativistic disks. It is illustrated for the Kerr disk.