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

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.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

The Lyman-$\alpha$ forest is a powerful probe for cosmology, but it is also strongly impacted by galaxy evolution and baryonic processes such as Active Galactic Nuclei (AGN) feedback, which can redistribute mass and energy on large scales. We constrain the signatures of AGN feedback on the 1D power spectrum of the Lyman-$\alpha$ forest using a series of eight hydro-cosmological simulations performed with the Adaptative Mesh Refinement code RAMSES. This series starts from the Horizon-AGN simulation and varies the sub-grid parameters for AGN feeding, feedback and stochasticity. These simulations cover the whole plausible range of feedback and feeding parameters according to the resulting galaxy properties. AGNs globally suppress the Lyman-$\alpha$ power at all scales. On large scales, the energy injection and ionization dominate over the supply of gas mass from AGN-driven galactic winds, thus suppressing power. On small scales, faster cooling of denser gas mitigates the suppression. This effect increases with decreasing redshift. We provide lower and upper limits of this signature at nine redshifts between $z=4.25$ and $z=2.0$, making it possible to account for it at post-processing stage in future work given that running simulations without AGN feedback can save considerable amounts of computing resources. Ignoring AGN feedback in cosmological inference analyses leads to strong biases with 2\% shift on $\sigma_8$ and 1\% shift on $n_s$, which represents twice the standards deviation of the current constraints on $n_s$.

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.

Supermassive black holes (BHs) of millions solar masses and above reside in the center of most local galaxies, but they also power active galactic nuclei and quasars, detected up to z=7. This observational evidence puts strong constraints on the BH growth and the mass of the first BH seeds. The scenario of "direct collapse" is very appealing as it leads to the formation of large mass BH seeds in the range 10^4-10^6 Msun, which eases explaining how quasars at z=6-7 are powered by BHs with masses >10^9 Msun. Direct collapse, however, appears to be rare, as the conditions required by the scenario are that gas is metal-free, the presence of a strong photo-dissociating Lyman-Werner flux, and large inflows of gas at the center of the halo, sustained for 10-100 Myr. We performed several cosmological hydrodynamical simulations that cover a large range of box sizes and resolutions, thus allowing us to understand the impact of several physical processes on the distribution of direct collapse BHs. We identify halos where direct collapse can happen, and derive the number density of BHs. We also investigate the discrepancies between hydrodynamical simulations, direct or post-processed, and semi-analytical studies. We find that for direct collapse to account for BHs in normal galaxies, the critical Lyman-Werner flux required for direct collapse must be much less than predicted by 3D simulations that include detailed chemical models. However, when supernova feedback is relatively weak, enough direct collapse BHs to explain z=6-7 quasars can be obtained for more realistic values of the critical Lyman Werner flux.

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.

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

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.

Using the Ly$\alpha$ mass assignment scheme (LyMAS), we make theoretical predictions for the 3-dimensional 3-point correlation function (3PCF) of the Ly$\alpha$ forest at redshift $z=2.3$. We bootstrap results from the (100 $h^{-1} \mbox{ Mpc}$)$^3$ Horizon hydrodynamic simulation to a (1 $h^{-1}$ Gpc)$^3$ $N$-body simulation, considering both a uniform UV background (UVB) and a fluctuating UVB sourced by quasars with a comoving $n_q \approx 10^{-5}$ $h^3$ Mpc$^{-3}$ placed either in massive halos or randomly. On scales of $10-30$ $h^{-1} \mbox{ Mpc}$, the flux 3PCF displays hierarchical scaling with the square of the 2PCF, but with an unusual value of $Q \equiv \zeta_{123}/(\xi_{12} \xi_{13} + \xi_{12} \xi_{23} + \xi_{13} \xi_{23}) \approx -4.5$ that reflects the low bias of the Ly$\alpha$ forest and the anti-correlation between mass density and transmitted flux. For halo-based quasars and an ionizing photon mean free path of $\lambda = 300$ $h^{-1} \mbox{ Mpc}$ comoving, UVB fluctuations moderately depress the 2PCF and 3PCF, with cancelling effects on $Q$. For $\lambda = 100$ $h^{-1} \mbox{ Mpc}$ or 50 $h^{-1} \mbox{ Mpc}$, UVB fluctuations substantially boost the 2PCF and 3PCF on large scales, shifting the hierarchical ratio to $Q \approx -3$. We scale our simulation results to derive rough estimate of the 3PCF detectability in observational data sets for the redshift range $z=2.1 - 2.6$. At $r = 10$ $h^{-1} \mbox{ Mpc}$ and 20 $h^{-1} \mbox{ Mpc}$, we predict a signal-to-noise (SNR) of $\sim$ 9 and $\sim$ 7, respectively, for both BOSS and eBOSS, and $\sim$ 37 and $\sim$ 25 for DESI. At $r = 40$ $h^{-1} \mbox{ Mpc}$ the predicted SNR is lower by $\sim$ 3$-$5 times. Measuring the flux 3PCF would be a novel test of the conventional paradigm of the Ly$\alpha$ forest and help separate the contributions of UVB fluctuations and density fluctuations to Ly$\alpha$ forest clustering.

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.

Multi-phase filamentary structures around Brightest Cluster Galaxies are likely a key step of AGN-feedback. We observed molecular gas in 3 cool cluster cores: Centaurus, Abell S1101, and RXJ1539.5 and gathered ALMA and MUSE data for 12 other clusters. Those observations show clumpy, massive and long, 3--25 kpc, molecular filaments, preferentially located around the radio bubbles inflated by the AGN (Active Galactic Nucleus). Two objects show nuclear molecular disks. The optical nebula is certainly tracing the warm envelopes of cold molecular filaments. Surprisingly, the radial profile of the H$\alpha$/CO flux ratio is roughly constant for most of the objects, suggesting that (i) between 1.2 to 7 times more cold gas could be present and (ii) local processes must be responsible for the excitation. Projected velocities are between 100--400 km s$^{-1}$, with disturbed kinematics and sometimes coherent gradients. This is likely due to the mixing in projection of several thin unresolved filaments. The velocity fields may be stirred by turbulence induced by bubbles, jets or merger-induced sloshing. Velocity and dispersions are low, below the escape velocity. Cold clouds should eventually fall back and fuel the AGN. We compare the filament's radial extent, r$_{fil}$, with the region where the X-ray gas can become thermally unstable. The filaments are always inside the low-entropy and short cooling time region, where t$_{cool}$/t$_{ff}$<20 (9 of 13 sources). The range t$_{cool}$/t$_{ff}$, 8-23 at r$_{fil}$, is likely due to (i) a more complex gravitational potential affecting the free-fall time (e.g., sloshing, mergers); (ii) the presence of inhomogeneities or uplifted gas in the ICM, affecting the cooling time. For some of the sources, r$_{fil}$ lies where the ratio of the cooling time to the eddy-turnover time, t$_{cool}$/t$_{eddy}$, is approximately unity.

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.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

In the mean field limit, isolated gravitational systems often evolve towards a steady state through a violent relaxation phase. One question is to understand the nature of this relaxation phase, in particular the role of radial instabilities in the establishment/destruction of the steady profile. Here, through a detailed phase-space analysis based both on a spherical Vlasov solver, a shell code and a $N$-body code, we revisit the evolution of collisionless self-gravitating spherical systems with initial power-law density profiles $\rho(r) \propto r^n$, $0 \leq n \leq -1.5$, and Gaussian velocity dispersion. Two sub-classes of models are considered, with initial virial ratios $\eta=0.5$ ("warm") and $\eta=0.1$ ("cool"). Thanks to the numerical techniques used and the high resolution of the simulations, our numerical analyses are able, for the first time, to show the clear separation between two or three well known dynamical phases: (i) the establishment of a spherical quasi-steady state through a violent relaxation phase during which the phase-space density displays a smooth spiral structure presenting a morphology consistent with predictions from self-similar dynamics, (ii) a quasi-steady state phase during which radial instabilities can take place at small scales and destroy the spiral structure but do not change quantitatively the properties of the phase-space distribution at the coarse grained level and (iii) relaxation to non spherical state due to radial orbit instabilities for $n \leq -1$ in the cool case.

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

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.

While various observations measured ellipticities of galaxy clusters and alignments between orientations of the brightest cluster galaxies and their host clusters, there are only a handful of numerical simulations that implement realistic baryon physics to allow direct comparisons with those observations. Here we investigate ellipticities of galaxy clusters and alignments between various components of them and the central galaxies in the state-of-the-art cosmological hydrodynamical simulation Horizon-AGN, which contains dark matter, stellar, and gas components in a large simulation box of $(100 h^{-1}$ Mpc$)^3$ with high spatial resolution ($\sim1$ kpc). We estimate ellipticities of total matter, dark matter, stellar, gas surface mass density distributions, X-ray surface brightness, and the Compton $y$-parameter of the Sunyaev-Zel'dovich effect, as well as alignments between these components and the central galaxies for 120 projected images of galaxy clusters with masses $M_{200}>5\times10^{13}M_{\odot}$. Our results indicate that the distributions of these components are well aligned with the major-axes of the central galaxies, with the root mean square value of differences of their position angles of $\sim 20^\circ$, which vary little from inner to the outer regions. We also estimate alignments of these various components with total matter distributions, and find tighter alignments than those for central galaxies with the root mean square value of $\sim 15^\circ$. We compare our results with previous observations of ellipticities and position angle alignments and find reasonable agreements. The comprehensive analysis presented in this paper provides useful prior information for analyzing stacked lensing signals as well as designing future observations to study ellipticities and alignments of galaxy clusters.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

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 fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

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.

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.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

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.

We report the discovery of a mysterious giant $H_{\alpha}$ blob that is $\sim 8$ kpc away from the main MaNGA target 1-24145, one component of a dry galaxy merger, identified in the first-year SDSS-IV MaNGA data. The size of the $H_{\alpha}$ blob is $\sim$ 3-4 kpc in radius, and the $H_{\alpha}$ distribution is centrally concentrated. However, there is no optical continuum counterpart in deep broadband images reaching $\sim$26.9 mag arcsec$^{-2}$ in surface brightness. We estimate that the masses of ionized and cold gases are $3.3 \times 10^{5}$ $\rm M_{\odot}$ and $< 1.3 \times 10^{9}$ $\rm M_{\odot}$, respectively. The emission-line ratios indicate that the $H_{\alpha}$ blob is photoionized by a combination of massive young stars and AGN. Furthermore, the ionization line ratio decreases from MaNGA 1-24145 to the $H_{\alpha}$ blob, suggesting that the primary ionizing source may come from MaNGA 1-24145, likely a low-activity AGN. Possible explanations of this $H_{\alpha}$ blob include AGN outflow, the gas remnant being tidally or ram-pressure stripped from MaNGA 1-24145, or an extremely low surface brightness (LSB) galaxy. However, the stripping scenario is less favoured according to galaxy merger simulations and the morphology of the $H_{\alpha}$ blob. With the current data, we can not distinguish whether this $H_{\alpha}$ blob is ejected gas due to a past AGN outburst, or a special category of `ultra-diffuse galaxy' (UDG) interacting with MaNGA 1-24145 that further induces the gas inflow to fuel the AGN in MaNGA 1-24145.

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.

The velocity distribution function (VDF) of dark matter (DM) halos in $\Lambda$CDM dissipationless cosmological simulations, which must be non-separable in its radial and tangential components, is still poorly known. We present the first single-parameter, non-separable, anisotropic model for the VDF in $\Lambda$CDM halos, built from an isotropic $q$-Gaussian (Tsallis) VDF of the isotropic set of dimensionless spherical velocity components (after subtraction of streaming motions), normalized by the respective velocity dispersions. We test our VDF on 90 cluster-mass halos of a dissipationless cosmological simulation. Beyond the virial radius, $r_{\rm vir}$, our model VDF adequately reproduces that measured in the simulated halos, but no $q$-Gaussian model can adequately represent the VDF within $r_{\rm vir}$, as the speed distribution function is then flatter-topped than any $q$-Gaussian can allow. Nevertheless, our VDF fits significantly better the simulations than the commonly used Maxwellian (Gaussian) distribution, at virtually all radii within $5\,r_{\rm vir}$. Within 0.4 (1) $r_{\rm vir}$, the non-Gaussianity index $q$ is (roughly) linearly related to the slope of the density profile and also to the velocity anisotropy profile. We provide a parametrization of the modulation of $q$ with radius for both the median fits and the fit of the stacked halo. At radii of a few percent of $r_{\rm vir}$, corresponding to the Solar position in the Milky Way, our best-fit VDF, although fitting better the simulations than the Gaussian one, overproduces significantly the fraction of high velocity objects, indicating that one should not blindly use these $q$-Gaussian fits to make predictions on the direct detection rate of DM particles.

We characterize the non-sphericity of galaxy clusters by the projected axis ratio of spatial distribution of star, dark matter, and X-ray surface brightness (XSB). We select 40 simulated groups and clusters of galaxies with mass larger than 5E13 Msun from the Horizon simulation that fully incorporates the relevant baryon physics, in particular, the AGN feedback. We find that the baryonic physics around the central region of galaxy clusters significantly affects the non-sphericity of dark matter distribution even beyond the central region, approximately up to the half of the virial radius. Therefore it is very difficult to predict the the probability density function (PDF) of the projected axis ratio of XSB from dark-matter only N-body simulations as attempted in previous studies. Indeed we find that the PDF derived from our simulated clusters exhibits much better agreement with that from the observed X-ray clusters. This indicates that our present methodology to estimate the non-sphericity directly from the Horizon simulation is useful and promising. Further improvements in both numerical modeling and observational data will establish the non-sphericity of clusters as a cosmological test complementary to more conventional statistics based on spherically averaged quantities.

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.

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.

Quiescent galaxies with little or no ongoing star formation dominate the galaxy population above $M_{*}\sim 2 \times 10^{10}~M_{\odot}$, where their numbers have increased by a factor of $\sim25$ since $z\sim2$. Once star formation is initially shut down, perhaps during the quasar phase of rapid accretion onto a supermassive black hole, an unknown mechanism must remove or heat subsequently accreted gas from stellar mass loss or mergers that would otherwise cool to form stars. Energy output from a black hole accreting at a low rate has been proposed, but observational evidence for this in the form of expanding hot gas shells is indirect and limited to radio galaxies at the centers of clusters, which are too rare to explain the vast majority of the quiescent population. Here we report bisymmetric emission features co-aligned with strong ionized gas velocity gradients from which we infer the presence of centrally-driven winds in typical quiescent galaxies that host low-luminosity active nuclei. These galaxies are surprisingly common, accounting for as much as $10\%$ of the population at $M_* \sim 2 \times 10^{10}~ M_{\odot}$. In a prototypical example, we calculate that the energy input from the galaxy's low-level active nucleus is capable of driving the observed wind, which contains sufficient mechanical energy to heat ambient, cooler gas (also detected) and thereby suppress star formation.

Modern cosmology predicts that a galaxy overdensity is associated to a large reservoir of the intergalactic gas, which can be traced by the Ly$\alpha$ forest absorption. We have undertaken a systematic study of the relation between Coherently Strong intergalactic Ly$\alpha$ Absorption systems (CoSLAs), which have highest optical depth ($\tau$) in $\tau$ distribution, and mass overdensities on the scales of $\sim$ 10 - 20 $h^{-1}$ comoving Mpc. On such large scales, our cosmological simulations show a strong correlation between the effective optical depth ($\tau_{\rm{eff}}$) of the CoSLAs and the 3-D mass overdensities. In moderate signal-to-noise spectra, however, the profiles of CoSLAs can be confused with high column density absorbers. For $z>2.6$, where the corresponding Ly$\beta$ is redshifted to the optical, we have developed the technique to differentiate between these two alternatives. We have applied this technique to SDSS-III quasar survey at $z = 2.6$ - 3.3, and we present a sample of five CoSLA candidates with $\tau_{\rm{eff}}$ on 15 $h^{-1}$ Mpc greater than $4.5\times$ the mean optical depth. At lower redshifts of $z < 2.6$, where the background quasar density is higher, the overdensity can be traced by intergalactic absorption groups using multiple sight lines. Our overdensity searches fully utilize the current and next generation of Ly$\alpha$ forest surveys which cover a survey volume of $> (1\ h^{-1}$ Gpc)$^3$. In addition, systems traced by CoSLAs will build a uniform sample of the most massive overdensities at $z > 2$ to constrain the models of structure formation, and offer a unique laboratory to study the interactions between galaxy overdensities and the intergalactic medium.

Quasars powered by massive black holes (BHs) with mass estimates above a billion solar masses have been identified at redshift 6 and beyond. The existence of such BHs requires almost continuous growth at the Eddington limit for their whole lifetime, of order of one billion years. In this paper, we explore the possibility that positively skewed scale-dependent non-Gaussian primordial fluctuations may ease the assembly of massive BHs. In particular, they produce more low-mass halos at high redshift, thus altering the production of metals and ultra-violet flux, believed to be important factors in BH formation. Additionally, a higher number of progenitors and of nearly equal-mass halo mergers would boost the mass increase provided by BH-BH mergers and merger-driven accretion. We use a set of two cosmological simulations, with either Gaussian or scale-dependent non-Gaussian primordial fluctuations to perform a proof-of-concept experiment to estimate how BH formation and growth are altered. We estimate the BH number density and the fraction of halos where BHs form, for both simulations and for two popular scenarios of BH formation (remnants of the first generation of stars and direct collapse in the absence of metals and molecular hydrogen). We find that the fractions of halos where BHs form are almost identical, but that non-Gaussian primordial perturbations increase the total number density of BHs for the both BH formation scenarios. We also evolve BHs using merger trees extracted from the simulations and find that non-Gaussianities increase both the BH masses and the number of the most massive BHs.

We perform a detailed comparison of the phase-space density traced by the particle distribution in Gadget simulations to the result obtained with a spherical Vlasov solver using the splitting algorithm. The systems considered are apodized H\'enon spheres with two values of the virial ratio, R ~ 0.1 and 0.5. After checking that spherical symmetry is well preserved by the N-body simulations, visual and quantitative comparisons are performed. In particular we introduce new statistics, correlators and entropic estimators, based on the likelihood of whether N-body simulations actually trace randomly the Vlasov phase-space density. When taking into account the limits of both the N-body and the Vlasov codes, namely collective effects due to the particle shot noise in the first case and diffusion and possible nonlinear instabilities due to finite resolution of the phase-space grid in the second case, we find a spectacular agreement between both methods, even in regions of phase-space where nontrivial physical instabilities develop. However, in the colder case, R=0.1, it was not possible to prove actual numerical convergence of the N-body results after a number of dynamical times, even with N=10$^8$ particles.

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 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 percent-level precision photometric redshifts in the UltraVISTA-DR1 near-infrared survey to investigate the changing relationship be- tween galaxy stellar mass and the dark matter haloes hosting them to $z\sim2$. We achieve this by measuring the clustering properties and abundances of a series of volume-limited galaxy samples selected by stellar mass and star-formation activity. We interpret these results in the framework of a phenomenological halo model and numerical simulations. Our measurements span a uniquely large range in stellar mass and redshift and reach below the characteristic stellar mass to $z\sim2$. Our results are: 1. At fixed redshift and scale, clustering amplitude depends monotonically on sample stellar mass threshold; 2. At fixed angular scale, the projected clustering amplitude decreases with redshift but the co-moving correlation length remains constant; 3. Characteristic halo masses and galaxy bias increase with increasing median stellar mass of the sample; 4. The slope of these relationships is modified in lower mass haloes; 5. Concerning the passive galaxy population, characteristic halo masses are consistent with a simply less-abundant version of the full galaxy sample, but at lower redshifts the fraction of satellite galaxies in the passive population is very different from the full galaxy sample; 6. Finally we find that the ratio between the characteristic halo mass and median stellar mass at each redshift bin reaches a peak at ${\log(M_h/M_\odot)}\sim12.2$ and the position of this peak remains constant out to $z\sim2$. The behaviour of the full and passively evolving galaxy samples can be understood qualitatively by considering the slow evolution of the characteristic stellar mass in the redshift range probed by our survey.

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.

[Abridged] Understanding how the intergalactic medium (IGM) was reionized at z > 6 is one of the big challenges of current high redshift astronomy. It requires modelling the collapse of the first astrophysical objects (Pop III stars, first galaxies) and their interaction with the IGM, while at the same time pushing current observational facilities to their limits. The observational and theoretical progress of the last few years have led to the emergence of a coherent picture in which the budget of hydrogen-ionizing photons is dominated by low-mass star-forming galaxies, with little contribution from Pop III stars and quasars. The reionization history of the Universe therefore critically depends on the number density of low-mass galaxies at high redshift. In this work, we explore how changes in the statistical properties of initial density fluctuations affect the formation of early galaxies. Following Habouzit et al. (2014), we run 5 N-body simulations with Gaussian and (scale-dependent) non-Gaussian initial conditions, all consistent with Planck constraints. By appealing to a galaxy formation model and to a population synthesis code, we compute the far-UV galaxy luminosity function down to M_UV = -14 at redshift 7 < z < 15. We find that models with strong primordial non-Gaussianities on < Mpc scales show a far-UV luminosity function significantly enhanced in low-mass galaxies. We adopt a reionization model calibrated from state-of-the-art hydrodynamical simulations and show that such non-Gaussianities leave a clear imprint on the Universe reionization history and electron Thomson scattering optical depth tau_E. Although current uncertainties in the physics of reionization and on the determination of tau_E still dominate the signatures of non-Gaussianities, our results suggest that tau_E could ultimately be used to constrain the statistical properties of initial density fluctuations.

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.

Primordial non-Gaussianities provide an important test of inflationary models. Although the Planck CMB experiment has produced strong limits on non-Gaussianity on scales of clusters, there is still room for considerable non-Gaussianity on galactic scales. We have tested the effect of local non-Gaussianity on the high redshift galaxy population by running five cosmological N-body simulations down to z=6.5. For these simulations, we adopt the same initial phases, and either Gaussian or scale-dependent non-Gaussian primordial fluctuations, all consistent with the constraints set by Planck on clusters scales. We then assign stellar masses to each halo using the halo - stellar mass empirical relation of Behroozi et al. (2013). Our simulations with non-Gaussian initial conditions produce halo mass functions that show clear departures from those obtained from the analogous simulations with Gaussian initial conditions at z>~10. We observe a >0.3 dex enhancement of the low-end of the halo mass function, which leads to a similar effect on the galaxy stellar mass function, which should be testable with future galaxy surveys at z>10. As cosmic reionization is thought to be driven by dwarf galaxies at high redshift, our findings may have implications for the reionization history of the Universe.

We study of the role of major mergers (mass ratios >1:4) in driving size growth in high-redshift (1<z<2) spheroidal galaxies (SGs) with stellar masses between 10^9.5 MSun and 10^10.7 MSun. This is a largely unexplored mass range at this epoch, containing the progenitors of more massive SGs on which the bulk of the size-evolution literature is based. We visually split our SGs into systems that are relaxed and those that exhibit tidal features indicative of a recent merger. Numerical simulations indicate that, given the depth of our images, only tidal features due to major mergers will be detectable at these epochs (features from minor mergers being too faint), making the disturbed SGs a useful route for estimating major-merger-driven size growth. The disturbed SGs are offset in size from their relaxed counterparts, lying close to the upper envelope of the local size -- mass relation. The mean size ratio of the disturbed SGs to their relaxed counterparts is ~2. Combining this observed size growth with empirical major-merger histories from the literature suggests that the size evolution of a significant fraction (around two-thirds) of SGs in this mass range could be driven by major mergers. If, as is likely, our galaxies are progenitors of more massive M > 10^10.7 MSun SGs at z<1, then major mergers are also likely to play an important role in the size growth of at least some massive SGs in this mass range.

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.

Warped disks are almost ubiquitous among spiral galaxies. Here we revisit and test the `fly-by scenario' of warp formation, in which impulsive encounters between galaxies are responsible for warped disks. Based on N-body simulations, we investigate the morphological and kinematical evolution of the stellar component of disks when galaxies undergo fly-by interactions with adjacent dark matter halos. We find that the so-called `S'-shaped warps can be excited by fly-bys and sustained for even up to a few billion years, and that this scenario provides a cohesive explanation for several key observations. We show that disk warp properties are governed primarily by the following three parameters; (1) the impact parameter, i.e., the minimum distance between two halos, (2) the mass ratio between two halos, and (3) the incident angle of the fly-by perturber. The warp angle is tied up with all three parameters, yet the warp lifetime is particularly sensitive to the incident angle of the perturber. Interestingly, the modeled S-shaped warps are often non-symmetric depending on the incident angle. We speculate that the puzzling U- and L-shaped warps are geometrically superimposed S-types produced by successive fly-bys with different incident angles, including multiple interactions with a satellite on a highly elongated orbit.

Detecting post-merger features of merger remnants is highly dependent on the depth of observation images. However, it has been poorly discussed how long the post-merger features are visible under different observational conditions. We investigate a merger-feature time useful for understanding the morphological transformation of galaxy mergers via numerical simulations. We use N-body/hydrodynamic simulations, including gas cooling, star formation, and supernova feedback. We run a set of simulations with various initial orbital configurations and with progenitor galaxies having different morphological properties mainly for equal-mass mergers. As reference models, we ran additional simulations for non-equal mass mergers and mergers in a large halo potential. Mock images using the SDSS $r$ band are synthesized to estimate a merger-feature times and compare it between the merger simulations. The mock images suggest that the post-merger features involve a small fraction of stars, and the merger-feature time depends on galaxy interactions. In an isolated environment, the merger-feature time is, on average, $\sim$ 2 times the final coalescence time for a shallow surface bright limit of 25 mag/arcsec^2. For a deeper surface brightness limit of 28 mag/arcsec^2, however, the merger-feature time is a factor of two longer, which is why the detection of post-merger features using shallow surveys has been difficult. Tidal force of a cluster potential is effective in stripping post-merger features out and reduces the merger-feature time.

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.

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

In order to understand the physical mechanisms at work during the formation of massive early-type galaxies, we performed six zoomed hydrodynamical cosmological simulations of halos in the mass range 4.3 10^12 < M_vir < 8.0 10^13 M_sun at z=0, using the Adaptive Mesh Refinement code RAMSES. These simulations explore the role of Active Galactic Nuclei (AGN), through jets powered by the accretion onto supermassive black holes on the formation of massive elliptical galaxies. In the absence of AGN feedback, large amounts of stars accumulate in the central galaxies to form overly massive, blue, compact and rotation-dominated galaxies. Powerful AGN jets transform the central galaxies into red extended and dispersion-dominated galaxies. This morphological transformation of disc galaxies into elliptical galaxies is driven by the efficient quenching of the in situ star formation due to AGN feedback, which transform these galaxies into systems built up by accretion. For galaxies mainly formed by accretion, the proportion of stars deposited farther away from the centre increases, and galaxies have larger sizes. The accretion is also directly responsible for randomising the stellar orbits, increasing the amount of dispersion over rotation of stars as a function of time. Finally, we find that our galaxies simulated with AGN feedback better match the observed scaling laws, such as the size-mass, velocity dispersion-mass, fundamental plane relations, and slope of the total density profiles at z~0, from dynamical and strong lensing constraints.

We study ~330 massive (M* > 10^9.5 MSun), newborn spheroidal galaxies (SGs) around the epoch of peak star formation (1<z<3), to explore the high-redshift origin of SGs and gain insight into when and how the old stellar populations that dominate today's Universe formed. The sample is drawn from the HST/WFC3 Early-Release Science programme, which provides deep 10-filter (0.2 - 1.7 micron) HST imaging over a third of the GOODS-South field. We find that the star formation episodes that built the SGs likely peaked in the redshift range 2<z<5 (with a median of z~3) and have decay timescales shorter than ~1.5 Gyr. Starburst timescales and ages show no trend with stellar mass in the range 10^9.5 < M* < 10^10.5 MSun. However, the timescales show increased scatter towards lower values (<0.3 Gyr) for M* > 10^10.5 MSun, and an age trend becomes evident in this mass regime: SGs with M* > 10^11.5 MSun are ~2 Gyrs older than their counterparts with M* < 10^10.5 MSun. Nevertheless, a smooth downsizing trend with galaxy mass is not observed, and the large scatter in starburst ages indicate that SGs are not a particularly coeval population. Around half of the blue SGs appear not to drive their star formation via major mergers, and those that have experienced a recent major merger, show only modest enhancements (~40%) in their specific star formation rates. Our empirical study indicates that processes other than major mergers (e.g. violent disk instability driven by cold streams and/or minor mergers) likely play a dominant role in building SGs, and creating a significant fraction of the old stellar populations that dominate today's Universe.

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.

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.

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.

(Abridged) Our aim is to investigate whether the presence of baryons can have any significant influence on the properties of the local Hubble flow which has proved to be "cold". We use two cosmological zoom simulations in the standard LCDM cosmology with the same set of initial conditions to study the formation of a local group-like system within a sphere of ~7 Mpc/h. The first one is a pure dark matter simulation (runDM) while a complete treatment of the physics of baryons is introduced in the second one (runB). We found that galaxies identified in runB and their corresponding dark matter haloes in runDM have very similar spatial distributions and dynamical properties on large scales. Then, when analyzing the velocity field and the deviation from a pure Hubble flow in both simulations, namely when computing the dispersion of peculiar velocities of galaxies \sigma*(R) and those of their corresponding dark matter haloes in runDM, we found no particular differences for distances R=1 to 8 Mpc from the local group mass center. The results indicate that the "true" \sigma*(R) values can be estimated from the pure dark matter simulation with a mean error of 3 km/s when dark matter haloes are selected with maximum circular velocities of Vc$\ge$30 km/s, corresponding to a population of dark matter haloes in runB that host galaxies. By investigating the properties of the Hubble flow at distances R~0.7 to 3 Mpc, we also found that the estimation of the total mass enclosed at the radius of the zero-velocity surface R0, using the spherical infall model adapted to LCDM, can be underestimated by at least 50%.

The star formation history of nearby early-type galaxies is investigated via numerical modelling. Idealized hydrodynamical N-body simulations with a star formation prescription are used to study the minor merger process between a giant galaxy (host) and a less massive spiral galaxy (satellite) with reasonable assumptions for the ages and metallicities of the merger progenitors. We find that the evolution of the star formation rate is extended over several dynamical times and shows peaks which correspond to pericentre passages of the satellite. The newly formed stars are mainly located in the central part of the satellite remnant while the older stars of the initial disk are deposited at larger radii in shell-like structures. After the final plunge of the satellite, star formation in the central part of the remnant can continue for several Gyrs depending on the star formation efficiency. Although the mass fraction in new stars is small, we find that the half-mass radius differs from the half-light radius in the V and H bands. Moreover synthetic 2D images in J, H, NUV, Hb and V bands, using the characteristic filters of the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST), reveal that residual star formation induced by gas-rich minor mergers can be clearly observed during and after the final plunge, especially in the NUV band, for interacting systems at (z<0.023) over moderate numbers of orbits (~2 orbits correspond to typical exposure times of ~3600 sec). This suggests that WFC3 has the potential to resolve these substructures, characterize plausible past merger episodes, and give clues to the formation of early-type galaxies.

At intermediate redshifts, many galaxies seem to be perturbed or suffering from an interaction. Considering that disk galaxies may have formed and evolved through minor mergers or through major mergers, it is important to understand the mechanisms at play during each type of merger in order to be able to establish the outcome of such an event. In some cases, only the use of both morphological and kinematical information can disentangle the actual configuration of an encounter at intermediate redshift. In this work, we present the morphological and kinematical analysis of a system at z=0.74 in order to understand its configuration, interacting stage and evolution. Using the integral field spectrograph GIRAFFE, long-slit spectroscopy by FORS2 and direct optical images from the HST-ACS and ISAAC near-infrared images, we disentangle the morphology of this system, its star-formation history and its extended kinematics in order to propose a possible configuration for the system. Numerical simulations are used to test different interacting scenarii. We identify this system as a face-on disk galaxy with a very bright bar in interaction with a smaller companion with a mass ratio of 3:1. The relevance of kinematical information and the constraints it imposes on the interpretation of the observations of distant galaxies is particularly strengthened in this case. This object is amongst the best example on how one may misinterpret morphology in the absence of kinematical information.

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.

(Context) In a Universe dominated by dark matter, halos are continuously accreting mass (violently or not) and such mechanism affects their dynamical state. (Aims) The evolution of dark matter halos in phase-space, and using the phase-space density indicator Q=rho/sigma^3 as a tracer, is discussed. (Methods) We have performed cosmological N-body simulations from which we have carried a detailed study of the evolution of ~35 dark halos in the interval 0<z<10. (Results)The follow up of individual halos indicates two distinct evolutionary phases. First, an early and fast decrease of Q associated to virialization after the gravitational collapse takes place. The nice agreement between simulated data and theoretical expectations based on the spherical collapse model support such a conjecture. The late and long period where a slow decrease of the phase-space density occurs is related to accretion and merger episodes. The study of some merger events in the phase-space (radial velocity versus radial distance) reveals the formation of structures quite similar to caustics generated in secondary infall models of halo formation. After mixing in phase-space, halos in quasi-equilibrium have flat-topped velocity distributions (negative kurtosis) with respect to Gaussians. The effect is more noticiable for captured satellites and/or substructures than for the host halo.

Abreg: By combining HST/UDF imagery with kinematics from VLT/GIRAFFE we derive a physical model of distant galaxy J033245.11-274724.0 in a way similar to what can be done in the nearby Universe. Here we study the properties of a distant compact LIRGs galaxy. Given the photometric and spectro photometric accuracies, we can decompose the galaxy in sub components and correct them for reddening. The galaxy is dominated by a dust enshrouded disk revealed by UDF imagery. The disk radius is half that of the Milky Way and the galaxy have a SFR=20Mo/yr. Morphology and kinematics show that gas and stars together spiral inwards rapidly to feed the disk and the central regions. A combined system of a bar and two non rotating spiral arms regulates the material accretion, induces large sigma, with sigma larger than 100 km/s and redistributes the angular momentum (AM). The detailed physical properties resemble to the expectations from modeling a merger of two equal mass, gaseous rich galaxies, 0.5 Gyr after the merger. In its later evolution, this galaxy could become a late type galaxy which falls on the T-F relation, with an AM mostly induced by the orbital AM of the merger.

We demonstrate that the large scatter in the ultra-violet (UV) colours of intermediate-mass early-type galaxies in the local Universe and the inferred low-level recent star formation in these objects can be reproduced by minor mergers in the standard LCDM cosmology. Numerical simulations of mergers with mass ratios less than or equal to 1:4, with reasonable assumptions for the ages, metallicities and dust properties of the merger progenitors, produce good agreement to the observed UV colours of the early-type population, if the infalling satellites are assumed to have (cold) gas fractions greater than 20%. Early-types that satisfy (NUV-r) < 3.8 are likely to have experienced mergers with mass ratios between 1:4 and 1:6 within the last ~1.5 Gyrs, while those that satisfy 3.8<(NUV-r)<5.5 are consistent with either recent mergers with mass ratios < 1:6 or mergers with higher mass ratios that occurred more than ~1.5 Gyrs in the past. We demonstrate that the early-type colour-magnitude relations and colour distributions in both the UV and optical spectral ranges are consistent with the expected frequency of minor merging activity in the standard LCDM cosmology at low redshift. We present a strong plausibility argument for minor mergers to be the principal mechanism behind the large UV scatter and associated low-level recent star formation observed in early-type galaxies in the nearby Universe.

(Aims) We present a physical model of the formation of J033239.72-275154.7, a galaxy observed at z=0.41 and characterized by a big young bar of size 6 kpc. The study of this system is particularly interesting for understanding the connection between mergers and bars as well as the properties and fate of this system as it relates to disk galaxy formation. (Methods) We compare the morphological and kinematic properties of J033239.72-275154.7, the latter obtained by the GIRAFFE spectrograph, to those derived from the merger of two spiral galaxies described by idealized N-body simulations including a star formation prescription. (Results) We found that the general morphological shape and most of the dynamical properties of the object can be well reproduced by a model in which the satellite is initially put in a retrograde orbit and the mass ratio of the system is 1:3. In such a scenario, a bar forms in the host galaxy after the first passage of the satellite where an important fraction of available gas is consumed in an induced burst. In its later evolution, however, we find that J033239.72-275154.7, whose major progenitor was an Sab galaxy, will probably become a S0 galaxy. This is mainly due to the violent relaxation and the angular momentum loss experienced by the host galaxy during the merger process, which is caused by the adopted orbital parameters. This result suggests that the building of the Hubble sequence is significantly influenced by the last major collision. In the present case, the merger leads to a severe damage of the disk of the progenitor, leading to an evolution towards a more bulge dominated galaxy.

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.

(context) Different cosmological data are consistent with an accelerated expansion produced by an exotic matter-energy component, dubbed "dark-energy''. A cosmological constant is a possibility since it satisfies most of the observational constraints. (aims) In this work, the consequences of such a component in the dynamics of groups of galaxies is investigated, aiming to detect possible effects in scales of the order of few Mpc. (methods) The Lema\^itre-Tolman model was modified by the inclusion of the cosmological constant term and, from the numerical solution of the equations of motion, a velocity-distance relation was obtained. This relation depends on two parameters: the central core mass and the Hubble parameter. The non-linear fit of such a relation to available data permitted to obtain masses for five nearby groups of galaxies and for the Virgo cluster as well as estimates of the Hubble constant. (results) The analysis of the present results indicates that the velocity-distance relation derived from the modified Lema\^itre-Tolman model as well as that derived from the "canonical'' model give equally acceptable fits to the existent data. As a consequence, any robust conclusion on the effects of the cosmological constant in the dynamics of groups could be established. The mean value of the Hubble parameter derived from the present study of local flows is H_0 = 65\pm 7 km/s/Mpc.

The relativistic accretion rate of dark matter by a black hole is revisited. Under the assumption that the phase space density indicator, $Q=\rho_{\infty}/\sigma^3_{\infty}$, remains constant during the inflow, the derived accretion rate can be higher up to five orders of magnitude than the classical accretion formula, valid for non-relativistic and non-interacting particles, when typical dark halo conditions are considered. For these typical conditions, the critical point of the flow is located at distances of about 30-150 times the horizon radius. Application of our results to black hole seeds hosted by halos issued from cosmological simulations indicate that dark matter contributes to no more than ~10% of the total accreted mass, confirming that the bolometric quasar luminosity is related to the baryonic accretion history of the black hole.

(Context) Central active galactic nuclei (AGN) are supposed to play a key role in the evolution of their host galaxies. In particular, the dynamical and physical properties of the gas core must be affected by the injected energy. (Aims) Our aim is to study the effects of an AGN on the dark matter profile and on the central stellar light distribution in massive early type galaxies. (Methods) By performing self-consistent N-body simulations, we assume in our analysis that periodic bipolar outbursts from a central AGN can induce harmonic oscillatory motions on both sides of the gas core. (Results) Using realistic AGN properties, we find that the motions of the gas core, driven by such feedback processes, can flatten the dark matter and/or stellar profiles after 4-5 Gyr. These results are consistent with recent observational studies that suggest that most giant elliptical galaxies have cores or are "missing light" in their inner part. Since stars behave as a "collisionless" fluid similar to dark matter, the density profile both of stars and dark matter should be affected in a similar way, leading to an effective reduction in the central brightness.

Cosmological N-body simulations were performed to study the evolution of the phase-space density Q = rho/sigma^3 of dark matter halos. No significant differences in the scale relations Q ~ sigma^(-2.1) or Q ~ M^(-0.82) are seen for "cold" or "warm" dark matter models. The follow up of individual halos from z = 10 up to the present time indicate the existence of two main evolutionary phases: an early and fast one (10 > z > 6.5), in which Q decreases on the average by a factor of 40 as a consequence of the randomization of bulk motions and a late and long one (6.5 > z > 0), in which Q decreases by a factor of 20 because of mixing induced by merger events. The study of these halos has also evidenced that rapid and positive variations of the velocity dispersion, induced by merger episods, are related to a fast decrease of the phase density Q.

Possible dark matter candidates are reviewed as well as indirect search methods based on annihilation or decay channels of these particles. Neutralino is presently the best particle candidate and its annihilation produces high energy neutrinos, antiprotons, positrons and gamma-rays. To date, only upper limits on neutrino fluxes from the center of the Earth or the Sun, were established by different experiments. Antiprotons detected by the BESS collaboration, if issued from the follow up hadronization of the annihilation process, exclude neutralino masses higher than 100 GeV. The EGRET gamma-ray residual emission seen at high galactic latitudes above 1 GeV could be explained by neutralino annihilations if: i) the dark matter profile is "cored" and ii) the neutralino mass is < 50 GeV. Sterile neutrinos in the keV mass range are a possible candidate to constitute warm dark matter. These particles may provide an adequate free streaming mass able to solve "some" difficulties present in the cold dark matter scenario at small scales and could also explain the natal kick of pulsars. MeV particles, dubbed "light" dark matter, proposed to explain the extended 511 keV line emission from the galactic center will also be discussed.

The spherical infall model first developed by Lema\^{i}tre and Tolman was modified in order to include the effects of a dark energy term. The resulting velocity-distance relation was evaluated numerically. This equation, when fitted to actual data, permits the simultaneous evaluation of the central mass and of the Hubble parameter. Application of this relation to the Local Group, when the dark energy is modeled by a cosmological constant, yields a total mass for the M31-Milky Way pair of (2.5 +/- 0.7) x 10^12 M\_sun, a Hubble parameter H\_0 = 74 +/- 4 km s^-1 Mpc^-1 and a 1-D velocity dispersion for the flow of about 39 km s^-1. The zero-velocity and the marginally bound surfaces of the Local Group are at about 1.0 and 2.3 Mpc respectively from the center of mass. A similar analysis for the Virgo cluster yields a mass of (1.10 +/- 0.12) x 10^15 M\_sun and H\_0 = 65 +/- 9 km s^-1 Mpc^-1. The zero-velocity is located at a distance of 8.6 +/- 0.8 Mpc from the center of the cluster. The predicted peculiar velocity of the Local Group towards Virgo is about 190 kms^-1, in agreement with other estimates. Slightly lower masses are derived if the dark energy is represented by a fluid with an equation of state P = w\epsilon with w = -2/3.

Possible indirect detection of neutralino, through its gamma-ray annihilation product, by the forthcoming GLAST satellite from our galactic halo, M31, M87 and the dwarf galaxies Draco and Sagittarius is studied. Gamma-ray fluxes are evaluated for the two representative energy thresholds, 0.1 GeV and 1.0 GeV, at which the spatial resolution of GLAST varies considerably. Apart from dwarfs which are described either by a modified Plummer profile or by a tidally-truncated King profiles, fluxes are compared for halos with central cusps and cores. It is demonstrated that substructures, irrespective of their profiles, enhance the gamma-ray emission only marginally. The expected gamma-ray intensity above 1 GeV at high galactic latitudes is consistent with the residual emission derived from EGRET data if the density profile has a central core and the neutralino mass is less than 50 GeV, whereas for a central cusp only a substantial enhancement would explain the observations. From M31, the flux can be detected above 0.1 GeV and 1.0 GeV by GLAST only if the neutralino mass is below 300 GeV and if the density profile has a central cusp, case in which a significant boost in the gamma-ray emission is produced by the central black hole. For Sagittarius, the flux above 0.1 GeV is detectable by GLAST provided the neutralino mass is below 50 GeV. From M87 and Draco the fluxes are always below the sensitivity limit of GLAST.

We present new results on the angular momentum evolution of dark matter halos. Halos, from N-body simulations, are classified according to their mass growth histories into two categories: the accretion category contains halos whose mass has varied continuously and smoothly, while the merger category contains halos which have undergone sudden and significant mass variations (greater than 1/3 of their initial mass per event). We observe for the growth of angular momentum that it continues well into the nonlinear regime, that it is correlated with mass growth and that it occurs in two phases: a fast followed by a slow phase. The spin parameter decreases with time for the accretion catalog whereas slightly increases for the merger catalog, but does not change when two catalogs are considered together, and furthermore it does not depend on halo mass or cosmological model. We also develop a simple model for the formation of a disk galaxy similar to Milky Way, and conclude that our own halo must have captured satellites in order to acquire the required angular momentum and to achieve most of the disk around z ~ 1.6. Angular momentum distribution indicates that at z ~ 1.6 only 22% of the halos have angular momentum of magnitude comparable to that of disk galaxies in the mass range of $(0.1-5)*10^{11}$ of solar mass, clearly insufficient to explain the present observed abundance of these objects.