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.

Upcoming wide-field surveys are well-suited to studying the growth of galaxy clusters by tracing galaxy and gas accretion along cosmic filaments. We use hydrodynamic simulations of volumes surrounding 324 clusters from \textsc{The ThreeHundred} project to develop a framework for identifying and characterising these filamentary structures, and associating galaxies with them. We define 3-dimensional reference filament networks reaching $5R_{200}$ based on the underlying gas distribution and quantify their recovery using mock galaxy samples mimicking observations such as those of the WEAVE Wide-Field Cluster Survey. Since massive galaxies trace filaments, they are best recovered by mass-weighting galaxies or imposing a bright limit (e.g. $>L^*$) on their selection. We measure the transverse gas density profile of filaments, derive a characteristic filament radius of $\simeq0.7$--$1~h^{-1}\rm{Mpc}$, and use this to assign galaxies to filaments. For different filament extraction methods we find that at $R>R_{200}$, $\sim15$--$20%$ of galaxies with $M_*>3 \times 10^9 M_{\odot}$ are in filaments, increasing to $\sim60%$ for galaxies more massive than the Milky-Way. The fraction of galaxies in filaments is independent of cluster mass and dynamical state, and is a function of cluster-centric distance, increasing from $\sim13$% at $5R_{200}$ to $\sim21$% at $1.5R_{200}$. As a bridge to the design of observational studies, we measure the purity and completeness of different filament galaxy selection strategies. Encouragingly, the overall 3-dimensional filament networks and $\sim67$% of the galaxies associated with them are recovered from 2-dimensional galaxy positions.

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

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

We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $\Delta z>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...

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

In an effort to better understand the formation of galaxy groups, we examine the kinematics of a large sample of spectroscopically confirmed X-ray galaxy groups in the Cosmic Evolution Survey (COSMOS) with a high sampling of galaxy group members up to $z=1$. We compare our results with predictions from the cosmological hydrodynamical simulation of {\sc Horizon-AGN}. Using a phase-space analysis of dynamics of groups with halo masses of $M_{\mathrm{200c}}\sim 10^{12.6}-10^{14.50}M_\odot$, we show that the brightest group galaxies (BGG) in low mass galaxy groups ($M_{\mathrm{200c}}<2 \times 10^{13} M_\odot$) have larger proper motions relative to the group velocity dispersion than high mass groups. The dispersion in the ratio of the BGG proper velocity to the velocity dispersion of the group, $\sigma_{\mathrm{BGG}}/\sigma_{group}$, is on average $1.48 \pm 0.13$ for low mass groups and $1.01 \pm 0.09$ for high mass groups. A comparative analysis of the {\sc Horizon-AGN} simulation reveals a similar increase in the spread of peculiar velocities of BGGs with decreasing group mass, though consistency in the amplitude, shape, and mode of the BGG peculiar velocity distribution is only achieved for high mass groups. The groups hosting a BGG with a large peculiar velocity are more likely to be offset from the $L_x-\sigma_{v}$ relation; this is probably because the peculiar motion of the BGG is influenced by the accretion of new members.

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

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

We searched for isolated dark matter deprived galaxies within several state-of-the-art hydrodynamical simulations: Illustris, IllustrisTNG, EAGLE, and Horizon-AGN and found a handful of promising objects in all except Horizon-AGN. While our initial goal was to study their properties and evolution, we quickly noticed that all of them were located at the edge of their respective simulation boxes. After carefully investigating these objects using the full particle data, we concluded that they are not merely caused by a problem with the algorithm identifying bound structures. We provide strong evidence that these oddballs were created from regular galaxies that get torn apart due to unphysical processes when crossing the edge of the simulation box. We show that these objects are smoking guns indicating an issue with the implementation of the periodic boundary conditions of the particle data in Illustris, IllustrisTNG, and EAGLE, which was eventually traced down to be a minor bug occurring for a very rare set of conditions.

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

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

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

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

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.

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

Galaxy clusters and groups are thought to accrete material along the preferred direction of cosmic filaments. Yet these structures have proven difficult to detect due to their low contrast with few studies focusing on cluster infall regions. In this work, we detected cosmic filaments around galaxy clusters using photometric redshifts in the range 0.15<z<0.7. We characterised galaxy populations in these structures to study the influence of "pre-processing" by cosmic filaments and galaxy groups on star-formation quenching. The cosmic filament detection was performed using the AMASCFI Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) T0007 cluster sample (Sarron et al. 2018). The filament reconstruction was done with the DISPERSE algorithm in photometric redshift slices. We showed that this reconstruction is reliable for a CFHTLS-like survey at 0.15<z<0.7 using a mock galaxy catalogue. We split our galaxy catalogue in two populations (passive and star-forming) using the LePhare SED fitting algorithm and worked with two redshift bins (0.15<z<=0.4 and 0.4<z<0.7). We showed that the AMASCFI cluster connectivity (i.e. the number of filaments connecting to a cluster) increases with cluster mass M200. Filament galaxies outside R200 are found to be closer to clusters at low redshift, whatever the galaxy type. Passive galaxies in filaments are closer to clusters than star-forming galaxies in the low redshift bin only. The passive fraction of galaxies decreases with increasing clustercentric distance up to d~5 cMpc. Galaxy groups/clusters that are not located at nodes of our reconstruction are mainly found inside cosmic filaments. These results give clues for "pre-processing" in cosmic filaments, that could be due to smaller galaxy groups. This trend could be further explored by applying this method to larger photometric surveys such as HSC-SPP or Euclid.

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.

We present the results of a search for galaxy clusters and groups in the $\sim2$ square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of $3\times10^{-16}\;ergs\;cm^{-2}\;s^{-1}$ in 0.5--2 keV range, and identify 247 X-ray groups with $M_{200c}=8\times10^{12}-3\times10^{14}\;M_{\odot}$ at a redshift range of $0.08\leq z<1.53$, using the multiband photometric redshift and the master spectroscopic redshift catalogues of the COSMOS. The X-ray centres of groups are determined using high-resolution Chandra imaging. We investigate the relations between the offset of the brightest group galaxies (BGGs) from halo X-ray centre and group properties and compare with predictions from semi-analytic models and hydrodynamical simulations. We find that BGG offset decreases with both increasing halo mass and decreasing redshift with no strong dependence on the X-ray flux and SNR. We show that the BGG offset decreases as a function of increasing magnitude gap with no considerable redshift dependent trend. The stellar mass of BGGs in observations extends over a wider dynamic range compared to model predictions. At $z<0.5$, the central dominant BGGs become more massive than those with large offsets by up to 0.3dex, in agreement with model prediction. The observed and predicted lognormal scatter in the stellar mass of both low- and large-offset BGGs at fixed halo mass is $\sim0.3$dex.

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.

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

We present a multi-wavelength catalog in the Subaru-XMM Deep Field (SXDF) as part of the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH). We include the newly acquired optical data from the Hyper-Suprime Cam Subaru Strategic Program, accompanied by IRAC coverage from the SPLASH survey. All available optical and near-infrared data is homogenized and resampled on a common astrometric reference frame. Source detection is done using a multi-wavelength detection image including the $u$-band to recover the bluest objects. We measure multi-wavelength photometry and compute photometric redshifts as well as stellar masses for $\sim$1.17 million objects over $\sim$4.2 deg$^2$ with $\sim$800,000 objects in the 2.4 deg$^2$ HSC-UltraDeep coverage. Using the available spectroscopic redshifts from various surveys over the range of $0<z<6$, we verify the performance of the photometric redshifts and we find a normalized median absolute deviation of 0.023 and outlier fraction of 3.2%. The SPLASH-SXDF catalog is a valuable, publicly available resource that is perfectly suited for studying galaxies in the early universe and tracing their evolution through cosmic time. The catalog is available for download from https://z.umn.edu/SXDF .

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

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

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

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.

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

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

We measure the stellar mass function (SMF) of galaxies in the COSMOS field up to $z\sim6$. We select them in the near-IR bands of the COSMOS2015 catalogue, which includes ultra-deep photometry from UltraVISTA-DR2, SPLASH, and Subaru/Hyper-SuprimeCam. At $z>2.5$ we use new precise photometric redshifts with error $\sigma_z=0.03(1+z)$ and an outlier fraction of $12\%$, estimated by means of the unique spectroscopic sample of COSMOS. The increased exposure time in the DR2, along with our panchromatic detection strategy, allow us to improve the stellar mass completeness at high $z$ with respect to previous UltraVISTA catalogues. We also identify passive galaxies through a robust colour-colour selection, extending their SMF estimate up to $z=4$. Our work provides a comprehensive view of galaxy stellar mass assembly between $z=0.1$ and 6, for the first time using consistent estimates across the entire redshift range. We fit these measurements with a Schechter function, correcting for Eddington bias. We compare the SMF fit with the halo mass function predicted from $\Lambda$CDM simulations. We find that at $z>3$ both functions decline with a similar slope in the high-mass end. This feature could be explained assuming that the mechanisms that quench star formation in massive haloes become less effective at high redshift; however further work needs to be done to confirm this scenario. Concerning the SMF low-mass end, it shows a progressive steepening as moving towards higher redshifts, with $\alpha$ decreasing from $-1.47_{-0.02}^{+0.02}$ at $z\simeq0.1$ to $-2.11_{-0.13}^{+0.30}$ at $z\simeq5$. This slope depends on the characterisation of the observational uncertainties, which is crucial to properly remove the Eddington bias. We show that there is currently no consensus on the method to quantify such errors: different error models result in different best-fit Schechter parameters. [Abridged]

We investigate the star formation properties of a large sample of ~2300 X-ray-selected Type 2 Active Galactic Nuclei (AGNs) host galaxies out to z~3 in the Chandra COSMOS Legacy Survey in order to understand the connection between the star formation and nuclear activity. Making use of the existing multi-wavelength photometric data available in the COSMOS field, we perform a multi-component modeling from far-infrared to near-ultraviolet using a nuclear dust torus model, a stellar population model and a starburst model of the spectral energy distributions (SEDs). Through detailed analysis of SEDs, we derive the stellar masses and the star formation rates (SFRs) of Type 2 AGN host galaxies. The stellar mass of our sample is in the range 9 < log M_{stellar}/M_{\odot} < 12 with uncertainties of ~0.19 dex. We find that Type 2 AGN host galaxies have, on average, similar SFRs compared to the normal star-forming galaxies with similar M_{stellar} and redshift ranges, suggesting no significant evidence for enhancement or quenching of star formation. This could be interpreted in a scenario, where the relative massive galaxies have already experienced substantial growth at higher redshift (z>3), and grow slowly through secular fueling processes hosting moderate-luminosity AGNs.

We examine the behaviour of the infrared-radio correlation (IRRC) over the range $0<z<6$ using new, highly sensitive 3GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2deg$^{2}$ COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4GHz radio luminosity ratio ($q_{\rm TIR}$), decreases with increasing redshift: $q_{\rm TIR}(z)=(2.88\pm0.03)(1+z)^{-0.19\pm0.01}$. Moderate-to-high radiative luminosity AGN do not follow the same $q_{\rm TIR}$$(z)$ trend, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed $q_{\rm TIR}(z)$ trend of the star-forming population. An increasing fractional contribution to the observed 3GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect $K$ corrections in the radio may govern the derived redshift trend of $q_{\rm TIR}$. Lastly, we present a redshift-dependent relation between rest-frame 1.4GHz radio luminosity and star formation rate taking the derived redshift trend into account.

ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample has known far-infrared luminosities and, hence, star formation rates (SFRs), and stellar masses (M$_{\rm *}$) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, M$_{\rm *}$, and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFR$_MS$. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolutionary and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function we estimate the accretion rates of gas required to maintain continuity of the MS evolution ($>100$ \msun yr$^{-1}$ at z $>$ 2.5). Simple power-law dependences are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and AGN activity at early epochs.

(abridged) We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project. The survey covers a 2.6sq.deg. area with a mean rms of ~2.3uJy/b, cataloging 10830 sources (>5sigma). Combining these radio data with optical, near-infrared (UltraVISTA), mid-infrared (Spitzer/IRAC) data, and X-ray data (Chandra), we find counterparts to radio sources for ~93% of the radio sample (in the areas of the COSMOS field not affected by saturated or bright sources in the optical to NIR bands), reaching out to z<6. We further classify the sources as star forming galaxies or AGN based on various criteria, such as X-ray luminosity, observed MIR color, UV-FIR spectral-energy distribution, rest-frame NUV-optical color corrected for dust extinction, and radio-excess relative to that expected from the hosts' star-formation rate. We separate the AGN into sub-samples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, candidates for high-redshift analogues to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these sub-populations down to radio flux levels of ~11uJy at 3GHz (or ~20uJy at 1.4GHz assuming a spectral index of -0.7). We find that the dominant fraction at 1.4GHz flux densities above ~200uJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100uJy the fraction of star-forming galaxies increases to ~60%, followed by the moderate-to-high radiative luminosity AGN (~20%), and low-to-moderate radiative luminosity AGN (~20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the SKA. Our estimates suggest that at the faint flux limits to be reached by the SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.

We make use of the deep Karl G. Jansky Very Large Array (VLA) COSMOS radio observations at 3 GHz to infer radio luminosity functions of star-forming galaxies up to redshifts of z~5 based on approximately 6000 detections with reliable optical counterparts. This is currently the largest radio-selected sample available out to z~5 across an area of 2 square degrees with a sensitivity of rms=2.3 ujy/beam. By fixing the faint and bright end shape of the radio luminosity function to the local values, we find a strong redshift trend that can be fitted with a pure luminosity evolution L~(1+z)^{(3.16 +- 0.2)-(0.32 +- 0.07) z}. We estimate star formation rates (SFRs) from our radio luminosities using an infrared (IR)-radio correlation that is redshift dependent. By integrating the parametric fits of the evolved luminosity function we calculate the cosmic SFR density (SFRD) history since z~5. Our data suggest that the SFRD history peaks between 2<z<3 and that the ultraluminous infrared galaxies (ULIRGs; 100 Msol/yr<SFR<1000 Msol/yr) contribute up to ~25% to the total SFRD in the same redshift range. Hyperluminous infrared galaxies (HyLIRGs; SFR>1000 Msol/yr) contribute an additional <2% in the entire observed redshift range. We find evidence of a potential underestimation of SFRD based on ultraviolet (UV) rest-frame observations of Lyman break galaxies (LBGs) at high redshifts (z>4) on the order of 15-20%, owing to appreciable star formation in highly dust-obscured galaxies, which might remain undetected in such UV observations.

We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.

We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.

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

We determine the radio size distribution of a large sample of 152 SMGs in COSMOS that were detected with ALMA at 1.3 mm. For this purpose, we used the observations taken by the VLA-COSMOS 3 GHz Large Project. One hundred and fifteen of the 152 target SMGs were found to have a 3 GHz counterpart. The median value of the major axis FWHM at 3 GHz is derived to be $4.6\pm0.4$ kpc. The radio sizes show no evolutionary trend with redshift, or difference between different galaxy morphologies. We also derived the spectral indices between 1.4 and 3 GHz, and 3 GHz brightness temperatures for the sources, and the median values were found to be $\alpha=-0.67$ and $T_{\rm B}=12.6\pm2$ K. Three of the target SMGs, which are also detected with the VLBA, show clearly higher brightness temperatures than the typical values. Although the observed radio emission appears to be predominantly powered by star formation and supernova activity, our results provide a strong indication of the presence of an AGN in the VLBA and X-ray-detected SMG AzTEC/C61. The median radio-emitting size we have derived is 1.5-3 times larger than the typical FIR dust-emitting sizes of SMGs, but similar to that of the SMGs' molecular gas component traced through mid-$J$ line emission of CO. The physical conditions of SMGs probably render the diffusion of cosmic-ray electrons inefficient, and hence an unlikely process to lead to the observed extended radio sizes. Instead, our results point towards a scenario where SMGs are driven by galaxy interactions and mergers. Besides triggering vigorous starbursts, galaxy collisions can also pull out the magnetised fluids from the interacting disks, and give rise to a taffy-like synchrotron-emitting bridge. This provides an explanation for the spatially extended radio emission of SMGs, and can also cause a deviation from the well-known IR-radio correlation.

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

We present a structural study of 182 obscured Active Galactic Nuclei (AGNs) at z<=1.5, selected in the COSMOS field from their extreme infrared to X-ray luminosity ratio and their negligible emission at optical wavelengths. We fit optical to far-infrared spectral energy distributions and analyze deep HST imaging to derive the physical and morphological properties of their host galaxies. We find that such galaxies are more compact than normal star-forming sources at similar redshift and stellar mass, and we show that it is not an observational bias related to the emission of the AGN. Based on the distribution of their UVJ colors, we also argue that this increased compactness is not due to the additional contribution of a passive bulge. We thus postulate that a vast majority of obscured AGNs reside in galaxies undergoing dynamical compaction, similar to processes recently invoked to explain the formation of compact star-forming sources at high redshift.

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

We present the clustering properties of a complete sample of 968 radio sources detected at 1.4 GHz by the VLA-COSMOS survey with radio fluxes brighter than 0.15 mJy. 92% have redshift determinations from the Laigle et al. (2016) catalogue. Based on their radio-luminosity, these objects have been divided into two populations of 644 AGN and 247 star-forming galaxies. By fixing the slope of the auto-correlation function to gamma=2, we find r_0=11.7^{+1.0}_{-1.1} Mpc for the clustering length of the whole sample, while r_0=11.2^{+2.5}_{-3.3} Mpc and r_0=7.8^{+1.6}_{-2.1} Mpc (r_0=6.8^{+1.4}_{-1.8} Mpc if we restrict our analysis to z<0.9) are respectively obtained for AGN and star-forming galaxies. These values correspond to minimum masses for dark matter haloes of M_min=10^[13.6^{+0.3}_{-0.6}] M_sun for radio-selected AGN and M_min=10^[13.1^{+0.4}_{-1.6}] M_sun for radio-emitting star-forming galaxies (M_min=10^[12.7^{+0.7}_{-2.2}] M_sun for z<0.9). Comparisons with previous works imply an independence of the clustering properties of the AGN population with respect to both radio luminosity and redshift. We also investigate the relationship between dark and luminous matter in both populations. We obtain <M*>/M_halo<~10^{-2.7} for AGN, and <M*>/M_halo<~10^{-2.4} in the case of star-forming galaxies. Furthermore, if we restrict to z<~0.9 star-forming galaxies, we derive <M*>/M_halo<~10^{-2.1}, result which clearly indicates the cosmic process of stellar build-up as one moves towards the more local universe. Comparisons between the observed space density of radio-selected AGN and that of dark matter haloes shows that about one in two haloes is associated with a black hole in its radio-active phase. This suggests that the radio-active phase is a recurrent phenomenon.

We characterise the physical nature of a 1.1 mm-selected, flux-limited, and interferometrically followed up sample of SMGs in COSMOS. We used the MAGPHYS code to fit the multiwavelength (UV-radio) SEDs of 16 of the target SMGs. We also constructed the pure radio SEDs of our SMGs using three different radio bands (325 MHz, 1.4 GHz, and 3 GHz). Moreover, since two SMGs in our sample, AzTEC1 and AzTEC3, benefit from previous CO line observations, we studied their properties in more detail. We found that 63% of our target SMGs lie above the galaxy main-sequence by more than a factor of 3, and hence are starbursts. The 3 GHz radio sizes we have previously measured for the target SMGs were compared with the present stellar mass estimates, and we found that the z>3 SMGs are fairly consistent with the mass-size relationship of z~2 compact, quiescent galaxies (cQGs). The median IR-radio correlation parameter is found to be q=2.27, which is lower than measured locally (median q=2.64). AzTEC1 is found to have a sub-Eddington SFR surface density (by a factor of 2.6), while AzTEC3 appears to be an Eddington-limited starburster. The gas reservoir in these two high-z SMGs would be exhausted in only ~86 and 19 Myr at the current SFR, respectively. A comparison of the MAGPHYS-based properties of our SMGs with those of equally bright ALESS SMGs suggests that the two populations share fairly similar physical characteristics, including the q parameter. A hint of negative correlation is found between the 3 GHz size and the level of starburstiness, and hence cosmic-ray electrons in more compact starbursts might be more susceptible to free-free absorption. Some of the derived low and high q values (compared to the local median) could be the result of a specific merger/post-starburst phase of galaxy evolution. Overall, our results support the scenario where z>3 SMGs evolve into today's giant ellipticals.

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

ALMA Cycle 2 observations of the long wavelength dust emission in 145 star-forming galaxies are used to probe the evolution of star-forming ISM. We also develop the physical basis and empirical calibration (with 72 low-z and z ~ 2 galaxies) for using the dust continuum as a quantitative probe of interstellar medium (ISM) masses. The galaxies with highest star formation rates (SFRs) at <z> = 2.2 and 4.4 have gas masses up to 100 times that of the Milky Way and gas mass fractions reaching 50 to 80%, i.e. gas masses 1 - 4 times their stellar masses. We find a single high-z star formation law: SFR = 35 M_ mol^0.89 x (1+z)_{z=2}^0.95 x (sSFR)_{MS}^0.23 \msun yr^-1 -- an approximately linear dependence on the ISM mass and an increased star formation efficiency per unit gas mass at higher redshift. Galaxies above the Main Sequence (MS) have larger gas masses but are converting their ISM into stars on a timescale only slightly shorter than those on the MS -- thus these 'starbursts' are largely the result of having greatly increased gas masses rather than and increased efficiency for converting gas to stars. At z $> 1$, the entire population of star-forming galaxies has $\sim$ 2 - 5 times shorter gas depletion times than low-z galaxies. These shorter depletion times indicate a different mode of star formation in the early universe -- most likely dynamically driven by compressive, high-dispersion gas motions -- a natural consequence of the high gas accretion rates.

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

We study the environment of 23 submillimetre galaxies (SMGs) drawn from the JCMT/AzTEC 1.1mm S/N-limited sample in the COSMOS field, as well as 4 COSMOS SMGs at z_spec>4.5, and 1 at z_spec=2.49, yielding a sample of 28 SMGs. We search for overdensities using the COSMOS photometric redshifts based on over 30 UV-NIR photometric bands, reaching an accuracy of sigma(Delta z/(1+z))=0.0067 (0.0155) at z<3.5 (>3.5). To identify overdensities we apply the Voronoi tessellation analysis, and estimate the overdensity estimator delta_g as a function of distance from the SMG and/or overdensity center. We test and validate our approach via simulations, X-ray detected groups, and spectroscopic verifications using VUDS and zCOSMOS catalogues showing that even with photometric redshifts in COSMOS we can efficiently retrieve overdensities out to z~5. Our results yield that 11/23 (48%) JCMT/AzTEC 1.1mm SMGs occupy overdense environments. Considering the entire JCMT/AzTEC 1.1mm S/N>4 sample, and accounting for the expected fraction of spurious detections, yields that 35-61% of the SMGs in the S/N-limited sample occupy overdense environments. We perform an X-ray stacking analysis in the 0.5-2keV band using a 32" aperture and our SMG positions, and find statistically significant detections. For our z<2 [z>2] subsample we find an average flux of (4.0+/-0.8)x10^{-16} [(1.3+/-0.5)x10^{-16}] erg/s/cm^2 and a corresponding total mass of M_200 = 2.8x10^{13} [2x10^{13}] MSol. Our results suggest a higher occurrence of SMGs occupying overdense environments at z>3, than at z<3. This may be understood if highly star forming galaxies can only be formed in the highest peaks of the density field tracing the most massive dark matter haloes at early cosmic epochs, while at later times cosmic structure may have matured sufficiently that more modest overdensities correspond to sufficiently massive haloes to form SMGs.

We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project, a radio continuum survey performed at 10 cm wavelength. The survey covers the full 2 square degree COSMOS field with mean $rms\sim2.3$ $\mu$Jy/beam, cataloging 10,899 source components above $5\times rms$. By combining these radio data with UltraVISTA, optical, near-infrared, and Spitzer/IRAC mid-infrared data, as well as X-ray data from the Chandra Legacy, and Chandra COSMOS surveys, we gain insight into the emission mechanisms within our radio sources out to redshifts of $z\sim5$. From these emission characteristics we classify our souces as star forming galaxies or AGN. Using their multi-wavelength properties we further separate the AGN into sub-samples dominated by radiatively efficient and inefficient AGN, often referred to as high- and low-excitation emission line AGN. We compare our method with other results based on fitting of the sources' spectral energy distributions using both galaxy and AGN spectral models, and those based on the infrared-radio correlation. We study the fractional contributions of these sub-populations down to radio flux levels of $\sim$10 $\mu$Jy. We find that at 3 GHz flux densities above $\sim$400 $\mu$Jy quiescent, red galaxies, consistent with the low-excitation radio AGN class constitute the dominant fraction. Below densities of $\sim$200 $\mu$Jy star-forming galaxies begin to constitute the largest fraction, followed by the low-excitation, and X-ray- and IR-identified high-excitation radio AGN.

We measure the H{\alpha} and [OIII] emission line properties as well as specific star-formation rates (sSFR) of spectroscopically confirmed 3<z<6 galaxies in COSMOS from their observed colors vs. redshift evolution. Our model describes consistently the ensemble of galaxies including intrinsic properties (age, metallicity, star-formation history), dust-attenuation, and optical emission lines. We forward-model the measured H{\alpha} equivalent-widths (EW) to obtain the sSFR out to z~6 without stellar mass fitting. We find a strongly increasing rest-frame H{\alpha} EW that is flattening off above z~2.5 with average EWs of 300-600A at z~6. The sSFR is increasing proportional to (1+z)^2.4 at z<2.2 and (1+z)^1.5 at higher redshifts, indicative of a fast mass build-up in high-z galaxies within e-folding times of 100-200Myr at z~6. The redshift evolution at z>3 cannot be fully explained in a picture of cold accretion driven growth. We find a progressively increasing [OIII]{\lambda}5007/H{\beta} ratio out to z~6, consistent with the ratios in local galaxies selected by increasing H{\alpha} EW (i.e., sSFR). This demonstrates the potential of using "local high-z analogs" to investigate the spectroscopic properties and relations of galaxies in the re-ionization epoch.

We use the stellar-mass-selected catalog from the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH) in the COSMOS field to study the environments of galaxies via galaxy density and clustering analyses up to $z \sim 2.5$. The clustering strength of quiescent galaxies exceeds that of star-forming galaxies, implying that quiescent galaxies are preferentially located in more massive halos. When using local density measurement, we find a clear positive quiescent fraction--density relation at $z < 1$, consistent with earlier results. However, the quiescent fraction--density relation reverses its trend at intermediate redshifts ($1 < z < 1.5$) with marginal significance (<1.8$\sigma$) and is found to be scale dependent (1.6$\sigma$). The lower fraction of quiescent galaxies seen in large-scale dense environments, if confirmed to be true, may be associated with the fact that the star formation can be more easily sustained via cold stream accretion in `large-scale' high-density regions, preventing galaxies from permanent quenching. Finally, at $z > 1.5$, the quiescent fraction depends little on the local density, even though clustering shows that quiescent galaxies are in more massive halos. We argue that at high redshift the typical halo size falls below $10^{13}$ $M_{\odot}$, where intrinsically the local density measurements are so varied that they do not trace the halo mass. Our results thus suggest that in the high-redshift Universe, halo mass may be the key in quenching the star formation in galaxies, rather than the conventionally measured galaxy density.

We present the catalog of optical and infrared counterparts of the Chandra COSMOS-Legacy Survey, a 4.6 Ms Chandra program on the 2.2 square degrees of the COSMOS field, combination of 56 new overlapping observations obtained in Cycle 14 with the previous C-COSMOS survey. In this Paper we report the i, K, and 3.6 micron identifications of the 2273 X-ray point sources detected in the new Cycle 14 observations. We use the likelihood ratio technique to derive the association of optical/infrared (IR) counterparts for 97% of the X-ray sources. We also update the information for the 1743 sources detected in C-COSMOS, using new K and 3.6 micron information not available when the C-COSMOS analysis was performed. The final catalog contains 4016 X-ray sources, 97% of which have an optical/IR counterpart and a photometric redshift, while 54% of the sources have a spectroscopic redshift. The full catalog, including spectroscopic and photometric redshifts and optical and X-ray properties described here in detail, is available online. We study several X-ray to optical (X/O) properties: with our large statistics we put better constraints on the X/O flux ratio locus, finding a shift towards faint optical magnitudes in both soft and hard X-ray band. We confirm the existence of a correlation between X/O and the the 2-10 keV luminosity for Type 2 sources. We extend to low luminosities the analysis of the correlation between the fraction of obscured AGN and the hard band luminosity, finding a different behavior between the optically and X-ray classified obscured fraction.

We measure a relation between the depth of four prominent rest-UV absorption complexes and metallicity for local galaxies and verify it up to z~3. We then apply this relation to a sample of 224 galaxies at 3.5 < z < 6.0 (<z> = 4.8) in COSMOS, for which unique UV spectra from DEIMOS and accurate stellar masses from SPLASH are available. The average galaxy population at z~5 and log(M/Msun) > 9 is characterized by 0.3-0.4 dex (in units of 12+log(O/H)) lower metallicities than at z~2, but comparable to z~3.5. We find galaxies with weak/no Ly-alpha emission to have metallicities comparable to z~2 galaxies and therefore may represent an evolved sub-population of z~5 galaxies. We find a correlation between metallicity and dust in good agreement with local galaxies and an inverse trend between metallicity and star-formation rate (SFR) consistent with observations at z~2. The relation between stellar mass and metallicity (MZ relation) is similar to z~3.5, however, there are indications of it being slightly shallower, in particular for the young, Ly-alpha emitting galaxies. We show that, within a "bathtub" approach, a shallower MZ relation is expected in the case of a fast (exponential) build-up of stellar mass with an e-folding time of 100-200 Myr. Due to this fast evolution, the process of dust production and metal enrichment as a function of mass could be more stochastic in the first billion years of galaxy formation compared to later times.

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

We have analysed a sample of 574 Spitzer 4.5 micron-selected galaxies with [4.5]<23 and Ks_auto>24 (AB) over the UltraVISTA ultra-deep COSMOS field. Our aim is to investigate whether these mid-IR bright, near-IR faint sources contribute significantly to the overall population of massive galaxies at redshifts z>=3. By performing a spectral energy distribution (SED) analysis using up to 30 photometric bands, we have determined that the redshift distribution of our sample peaks at redshifts z~2.5-3.0, and ~32% of the galaxies lie at z>=3. We have studied the contribution of these sources to the galaxy stellar mass function (GSMF) at high redshifts. We found that the [4.5]<23, Ks_auto>24 galaxies produce a negligible change to the GSMF previously determined for Ks_auto<24 sources at 3=<z<4, but their contribution is more important at 4=<z<5, accounting for >~50% of the galaxies with stellar masses Mst>~6 x 10^10 Msun. We also constrained the GSMF at the highest-mass end (Mst>~2 x 10^11 Msun) at z>=5. From their presence at 5=<z<6, and virtual absence at higher redshifts, we can pinpoint quite precisely the moment of appearance of the first most massive galaxies as taking place in the ~0.2 Gyr of elapsed time between z~6 and z~5. Alternatively, if very massive galaxies existed earlier in cosmic time, they should have been significantly dust-obscured to lie beyond the detection limits of current, large-area, deep near-IR surveys.

Current constraints on dark matter density profiles from weak lensing are typically limited to radial scales greater than 50-100 kpc. In this paper, we explore the possibility of probing the very inner regions of galaxy/halo density profiles by measuring stacked weak lensing on scales of only a few tens of kpc. Our forecasts focus on scales smaller than the equality radius (Req) where the stellar component and the dark matter component contribute equally to the lensing signal. We compute the evolution of Req as a function of lens stellar mass and redshift and show that Req=7-34 kpc for galaxies with the stellar mass of 10^{9.5}-10^{11.5} solar masses. Unbiased shear measurements will be challenging on these scales. We introduce a simple metric to quantify how many source galaxies overlap with their neighbours and for which shear measurements will be challenging. Rejecting source galaxies with close-by companions results in about a 20 per cent decrease in the overall source density. Despite this decrease, we show that Euclid and WFIRST will be able to constrain galaxy/halo density profiles at Req with signal-to-noise ratio >20 for the stellar mass of >10^{10} solar masses. Weak lensing measurements at Req, in combination with stellar kinematics on smaller scales, will be a powerful means by which to constrain both the inner slope of the dark matter density profile as well as the mass and redshift dependence of the stellar initial mass function.

We study the physical properties of a sample of 6 SMGs in the COSMOS field, spectroscopically confirmed to lie at z>4. We use new GMRT 325 MHz and 3 GHz JVLA data to probe the rest-frame 1.4 GHz emission at z=4, and to estimate the sizes of the star-forming (SF) regions of these sources, resp. Combining our size estimates with those available in the literature for AzTEC1 and AzTEC3 we infer a median radio-emitting size for our z>4 SMGs of (0.63"+/-0.12")x(0.35"+/-0.05") or 4.1x2.3 kpc^2 (major times minor axis; assuming z=4.5) or lower if we take the two marginally resolved SMGs as unresolved. This is consistent with the sizes of SF regions in lower-redshift SMGs, and local normal galaxies, yet higher than the sizes of SF regions of local ULIRGs. Our SMG sample consists of a fair mix of compact and more clumpy systems with multiple, perhaps merging, components. With an average formation time of ~280 Myr, derived through modeling of the UV-IR SEDs, the studied SMGs are young systems. The average stellar mass, dust temperature, and IR luminosity we derive are M*~1.4x10^11 M_sun, T_dust~43 K, and L_IR~1.3x10^13L_sun, resp. The average L_IR is up to an order of magnitude higher than for SMGs at lower redshifts. Our SMGs follow the correlation between dust temperature and IR luminosity as derived for Herschel-selected 0.1<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=1.95+/-0.26 for our sample, compared to q~2.6 for IR luminous galaxies at z<2). In summary, we find that the physical properties derived for our z>4 SMGs put them at the high end of the L_IR-T_dust distribution of SMGs, and that our SMGs form a morphologically heterogeneous sample. Thus, further in-depth analyses of large, statistical samples of high-redshift SMGs are needed to fully understand their role in galaxy formation and evolution.

We present the cluster mass-richness scaling relation calibrated by a weak lensing analysis of >18000 galaxy cluster candidates in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). Detected using the 3D-Matched-Filter cluster-finder of Milkeraitis et al., these cluster candidates span a wide range of masses, from the small group scale up to $\sim10^{15} M_{\odot}$, and redshifts 0.2 $\lesssim z\lesssim$ 0.9. The total significance of the stacked shear measurement amounts to 54$\sigma$. We compare cluster masses determined using weak lensing shear and magnification, finding the measurements in individual richness bins to yield 1$\sigma$ compatibility, but with magnification estimates biased low. This first direct mass comparison yields important insights for improving the systematics handling of future lensing magnification work. In addition, we confirm analyses that suggest cluster miscentring has an important effect on the observed 3D-MF halo profiles, and we quantify this by fitting for projected cluster centroid offsets, which are typically $\sim$ 0.4 arcmin. We bin the cluster candidates as a function of redshift, finding similar cluster masses and richness across the full range up to $z \sim$ 0.9. We measure the 3D-MF mass-richness scaling relation $M_{200} = M_0 (N_{200} / 20)^\beta$. We find a normalization $M_0 \sim (2.7^{+0.5}_{-0.4}) \times 10^{13} M_{\odot}$, and a logarithmic slope of $\beta \sim 1.4 \pm 0.1$, both of which are in 1$\sigma$ agreement with results from the magnification analysis. We find no evidence for a redshift-dependence of the normalization. The CFHTLenS 3D-MF cluster catalogue is now available at cfhtlens.org.

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.

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.

Gravitational lensing magnification is measured with a significance of 9.7 sigma on a large sample of galaxy clusters in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). This survey covers ~154 deg^2 and contains over 18,000 cluster candidates at redshifts 0.2 <= z <= 0.9, detected using the 3D-Matched Filter cluster-finder of Milkeraitis et al. (2010). We fit composite-NFW models to the ensemble, accounting for cluster miscentering, source-lens redshift overlap, as well as nearby structure (the 2-halo term), and recover mass estimates of the cluster dark matter halos in range of ~10^13 M_sun to 2*10^14 M_sun. Cluster richness is measured for the entire sample, and we bin the clusters according to both richness and redshift. A mass-richness relation M_200 = M_0 (N_200 / 20)^beta is fit to the measurements. For two different cluster miscentering models we find consistent results for the normalization and slope, M_0 = (2.3 +/- 0.2)*10^13 M_sun, beta = 1.4 +/- 0.1 and M_0 = (2.2 +/- 0.2)*10^13 M_sun, beta = 1.5 +/- 0.1. We find that accounting for the full redshift distribution of lenses and sources is important, since any overlap can have an impact on mass estimates inferred from flux magnification.