The halo spin flip refers to the phenomenon that the spin axes of dark matter halos with masses above a certain threshold tend to be preferentially aligned perpendicular to the hosting large-scale filaments, while low-mass halos tend to have their spin axes aligned parallel to such structures. Extensive work has so far been conducted to understand this phenomenon under the assumption of cold dark matter and suggested that its origin should be closely related to the nonlinear evolution of the halo angular momentum in the anisotropic cosmic web. We present, for the first time, a numerical examination of this phenomenon assuming the presence of massive neutrinos, finding a clear and robust dependence of the threshold mass for the spin flip on the total neutrino mass. Our physical explanation is that the presence of more massive neutrinos retard the nonlinear evolution of the cosmic web, which in turn allows the halo spin vectors to better retain their memories of the initial tidal interactions in the nonlinear regime. Our finding implies that the statistical alignment of halo spins with the large-scale structures can be in principle used as a probe of the total neutrino mass.

We perform an analysis of the Cosmic Web as a complex network, which is built on a $\Lambda$CDM cosmological simulation. For each of nodes, which are in this case dark matter halos formed in the simulation, we compute 10 network metrics, which characterize the role and position of a node in the network. The relation of these metrics to topological affiliation of the halo, i.e. to the type of large scale structure, which it belongs to, is then investigated. In particular, the correlation coefficients between network metrics and topology classes are computed. We have applied different machine learning methods to test the predictive power of obtained network metrics and to check if one could use network analysis as a tool for establishing topology of the large scale structure of the Universe. Results of such predictions, combined in the confusion matrix, show that it is not possible to give a good prediction of the topology of Cosmic Web (score is $\approx$ 70 $\%$ in average) based only on coordinates and velocities of nodes (halos), yet network metrics can give a hint about the topological landscape of matter distribution.

This article explores in depth a watershed concept to partition the universe, introduced in \cite{2019MNRAS.489L...1D} and applied to the {\it Cosmicflows-3} observational dataset. We present a series of tests conducted with cosmological dark matter simulations. In particular we are interested in quantifying the evolution with redshift of large scale structures when defined as segmented basins of attraction. This new dynamical definition in the field of measuring standard rulers demonstrates robustness since all basins show a density contrast $\delta$ above one (mean universe density) independently of the simulation spatial resolution or the redshift. Another major finding is that density profiles of the basins show universality in slope. Consequently, there is a unique definition of what is a gravitational watershed at large scale, that can be further used as a probe for cosmology studies.

We study the formation of polar ring galaxies via minor mergers. We used N-body+hydrodynamics simulations to reproduce the dynamics of the observed system AM 2229-735 that is a minor merger whose interaction signals are those of a progenitor for a polar ring galaxy. We used the observational information of the system to get initial conditions for the orbit and numerical realisations of the galaxies to run the simulations. Our simulations reproduce the global characteristics of interaction observed in the system such as arms and a material bridge connecting the galaxies. As a merger remnant, we found a quasi-stable and self gravitating planar tidal stream with dark matter, stars and gas orbiting in a plane approximately perpendicular to the main galaxy disk leading in the future to a polar ring galaxy. We studied the dynamical conditions of the polar structure and found evidence suggesting that this kind of merger remnant can settle down in a disk-like structure with isothermal support, providing inspiring evidence about the process of formation of galactic disks and providing a potentially independent scenario to study the presence of dark matter in this kind of galaxies.

In this work an ensemble of simulated Local Group analogues is used to constrain the properties of the mass assembly history of the Milky Way (MW) and Andromeda (M31) galaxies. These objects have been obtained using the constrained simulation technique, which ensures that simulated LGs live within a large scale environment akin to the observed one. Our results are compared against a standard $\Lambda$ Cold Dark Matter ($\Lambda$CDM) series of simulations which use the same cosmological parameters. This allows us to single out the effects of the constraints on the results. We find that (a) the median constrained merging histories for M31 and MW live above the standard ones at the 1-$\sigma$ level, (b) the median formation time takes place $\approx$ 0.5 Gyr earlier than unconstrained values, while the latest major merger happens on average 1.5 Gyr earlier and (c) the probability for both LG haloes to have experienced their last major merger in the first half of the history of the Universe is $\approx$ 50% higher for the constrained pairs. These results have been estimated to be significant at the 99% confidence level by means of a Kolmogorov-Simirnov test. These results suggest that the particular environment in which the Milky Way and Andromeda form plays a role in shaping their properties, and favours earlier formation and last major merger time values in agreement with other observational and theoretical considerations.

The universe region obscured by the Milky Way is very large and only future blind large HI redshift, and targeted peculiar surveys on the outer borders will determine how much mass is hidden there. Meanwhile, we apply for the first time two independent techniques to the galaxy peculiar velocity catalog $CosmicFlows-3$ in order to explore for the kinematic signature of a specific large-scale structure hidden behind this zone : the Vela supercluster at cz $\sim 18,000$,km s$^{-1}$ . Using the gravitational velocity and density contrast fields, we find excellent agreement when comparing our results to the Vela object as traced in redshift space. The article provides the first kinematic evidence of a major mass concentration (knot of the Cosmic Web) located in the direction behind Vela constellation, pin-pointing that the Zone of Avoidance should be surveyed in detail in the future .

To date at least 10 highly flattened planes of dwarf galaxies are claimed to have been discovered in the Local Universe. The origin of these planes of galaxies remains unknown. One suggestion is that they are related to the large-scale structure of the cosmic web. A recent study found that the normal of a number of these dwarf galaxy planes are very closely aligned with the eigenvector of the shear tensor corresponding to the direction of greatest collapse obtained by reconstructing the full velocity field in the linear regime. Here we extend that work by both considering an additional 5 planes beyond the 5 examined previously and by examining the alignment with respect to the quasi-linear field, a more sophisticated reconstruction, which is a better approximation on smaller (quasi-linear) scales. Our analysis recovers the previous result while not finding a significantly tight alignment with the additional 5 planes. However, the additional 5 planes normals also do not appear to be randomly oriented. We conclude that this could either be due to the normals of the new planes being poorly defined and described; the quasi-linear field at those locations being poorly constrained; or different formation mechanisms for the orientation of planes of dwarf galaxies.

In this work, we explore the idea that substructures like stellar clusters could be formed from the tidal stream produced in galactic minor mergers. We use $N$-body and SPH simulations of satellite galaxies interacting with a larger galaxy. We study the distribution of mass in streams to identify overdensity regions in which a substructure could be formed. We found that without gas, no substructure formed as none of the overdensities shows a definite morphology nor dynamical stability. Including gas we found that several clumps appear and proved to be real long standing physical structures ($t \geq$ 1 Gyr). We analyzed the orbits, ages and masses of these structures, finding its correspondence with the halo subsystems. We conclude that it is possible to form cluster-like structures from the material in tidal streams and found evidence in favour of the presence of dark matter in these systems.

This letter presents a new approach using the cosmic peculiar velocity field to characterize the morphology and size of large scale structures in the local Universe. The algorithm developed uses the three-dimensional peculiar velocity field to compute flow lines, or streamlines. The local Universe is then partitioned into volumes corresponding to gravitational basins, also called watersheds, among the different end-points of the velocity flow lines. This new methodology is first tested on numerical cosmological simulations, used as benchmark for the method, and then applied to the {\it Cosmic-Flows} project observational data in order to to pay particular attention to the nearby superclusters including ours. More extensive tests on both simulated and observational data will be discussed in an accompanying paper.

It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the $\Lambda$CDM model. One type of anisotropy recently detected in the SDSS (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called "lopsided". Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set up when comparing pairs of galaxies in $\Lambda$CDM simulations indicating that such a set up is not inconsistent with $\Lambda$CDM, the origin has yet to be explained. Here we examine the origin of such lopsided setups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at $z=0$ display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards $z=0$. Second, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the $z=0$ signal is driven primarily by satellites that are on first approach, who have yet to experience a "flyby". This suggests that the signal seen in the observations is also dominated by dynamically young accretion events.

It has been shown, both in simulations and observationally, that the tidal field of a large galaxy can torque its satellites such that the major axis of satellite galaxies points towards their hosts. This so-called `shape alignment' has been observed in isolated Milky Way-like galaxies but not in `Local Group'-like pairs. In this study, we investigate the shape alignment of satellite galaxies in galaxy pairs similar to the Local Group identified in the Sloan Digital Sky Survey Data Release 13 (SDSS DR13). By stacking tens of thousands of satellite galaxies around primary galaxy pairs, we find two statistically strong alignment signals. (1) The major axes of satellite galaxies located in the (projected) area between two primaries (the {\it facing} region) tend to be perpendicular to the line connecting the satellite to its host (tangential alignment), while (2) the major axes of satellite galaxies located in regions away from the other host (the {\it away} region) tend to be aligned with the line connecting the satellite to its host (radial alignment). These alignments are confirmed at $\sim5\sigma$ levels. The alignment signal increases with increasing primary brightness, decreasing pair separation, and decreasing satellite distance. The alignment signal is also found to be stronger in filamentary environments. These findings will shed light on understanding the mechanisms of how satellite galaxies are affected by the tidal field in galaxy pairs and will be useful for investigating galaxy intrinsic alignment in the analyses of weak gravitational lensing.

Following Cui et al. 2018 (hereafter Paper I) on the classification of large-scale environments (LSE) at z = 0, we push our analysis to higher redshifts and study the evolution of LSE and the baryon distributions in them. Our aim is to investigate how baryons affect the LSE as a function of redshift. In agreement with Paper I, the baryon models have negligible effect on the LSE over all investigated redshifts. We further validate the conclusion obtained in Paper I that the gas web is an unbiased tracer of total matter -- even better at high redshifts. By separating the gas mainly by temperature, we find that about 40 per cent of gas is in the so-called warm-hot intergalactic medium (WHIM). This fraction of gas mass in the WHIM decreases with redshift, especially from z = 1 (29 per cent) to z = 2.1 (10 per cent). By separating the whole WHIM gas mass into the four large-scale environments (i.e. voids, sheets, filaments, and knots), we find that about half of the WHIM gas is located in filaments. Although the total gas mass in WHIM decreases with redshift, the WHIM mass fractions in the different LSE seem unchanged.

The correlation between the spins of dark matter halos and the large-scale structure (LSS) has been studied in great detail over a large redshift range, while investigations of galaxies are still incomplete. Motivated by this point, we use the state-of-the-art hydrodynamic simulation, Illustris-1, to investigate mainly the spin--LSS correlation of galaxies at redshift of $z=0$. We mainly find that the spins of low-mass, blue, oblate galaxies are preferentially aligned with the slowest collapsing direction ($e_3$) of the large-scale tidal field, while massive, red, prolate galaxy spins tend to be perpendicular to $e_3$. The transition from a parallel to a perpendicular trend occurs at $\sim10^{9.4} M_{\odot}/h$ in the stellar mass, $\sim0.62$ in the g-r color, and $\sim0.4$ in triaxiality. The transition stellar mass decreases with increasing redshifts. The alignment was found to be primarily correlated with the galaxy stellar mass. Our results are consistent with previous studies both in N-body simulations and observations. Our study also fills the vacancy in the study of the galaxy spin--LSS correlation at $z=0$ using hydrodynamical simulations and also provides important insight to understand the formation and evolution of galaxy angular momentum.

HST/ACS observations along the major axis of M33 show that the mean age of its stars decreases with increasing distance from the galaxy center. Such a behavior is consistent with an inside-out growth of the disc. However, in the outermost observed field, at r$\simeq$11.6 kpc, a reversal of this gradient is detected, with old stars found in high percentages beyond this radius. In this work we investigate the origin of such a reversal in stellar age gradient, by using a simulated M33 analogue from the Constrained Local UniversE Simulations (CLUES). The simulated M33 is similar to the observed one in terms of mass, rotation velocity, surface brightness and, similar to what has been reported in observations, shows a stellar age turnaround at large radii. We demonstrate that this reversal is mostly a result of stellar accretion from old satellite galaxies and, to a lesser extent, of stellar migration of in-situ stars. The old accreted stars, with formation times tf < 4 Gyrs, are kinematically hot and can be differentiated from the in-situ stars by their high velocity dispersion and the fact that they do not have rotationally-supported orbits. In the future, obtaining kinematic information of the stars in the outskirt of M33 will help to verify this scenario.

The local Universe provides a unique opportunity for testing cosmology and theories of structure formation. To facilitate this opportunity we present a new method for the reconstruction of the quasi-linear matter density and velocity fields from galaxy peculiar velocities and apply it to the Cosmicflows-2 data. The method consists of constructing an ensemble of cosmological simulations, constrained by the standard cosmological model and the observational data. The quasi-linear density field is the geometric mean and variance of the fully non-linear density fields of the simulations. The main nearby clusters (Virgo, Centaurus, Coma), superclusters (Shapley, Perseus-Pisces) and voids (Dipole Repeller) are robustly reconstructed. Galaxies are born biased with respect to the underlying dark matter distribution. Using our quasi-linear framework we demonstrate that the luminosity-weighted density field derived from the 2M++ redshift compilations is non-linearly biased with respect to the matter density field. The bias diminishes in the linear regime.

The alignment between satellites and central galaxies has been studied in detail both in observational and theoretical works. The widely accepted fact is that the satellites preferentially reside along the major axis of their central galaxy. However, the origin and large-scale environment dependence of this alignment are still unknown. In an attempt to figure out those, we use data constructed from SDSS DR7 to investigate the large-scale environmental dependence of this alignment with emphasis on examining the alignments' dependence on the colour of the central galaxy. We find a very strong large-scale environmental dependence of the satellite-central alignment in groups with blue centrals. Satellites of blue centrals in knots are preferentially located perpendicular to the major axis of the centrals, and the alignment angle decreases with environment namely when going from knots to voids. The alignment angle strongly depend on the ${}^{0.1}(g-r)$ colour of centrals. We suggest that the satellite-central alignment is the result of a competition between satellite accretion within large scale-structure and galaxy evolution inside host haloes. For groups containing red central galaxies, the satellite-central alignment is mainly determined by the evolution effect, while for blue central dominated groups, the effect of large-scale structure plays a more important role, especially in knots. Our results provide an explanation for how the satellite-central alignment forms within different large-scale environments. The perpendicular case in groups and knots with blue centrals may also provide insight into understanding similar polar arrangements such the formation of the Milky Way and Centaurus A's satellite system.

The Milky Way's satellite galaxies and Globular Clusters (GCs) are known to exhibit an anisotropic spatial distribution. We examine in detail this anisotropy by the means of the inertia tensor. We estimate the statistical significance of the results by repeating this analysis for random catalogues which use the radial distribution of the real sample. Our method reproduces the well-known planar structure in the distribution of the satellite galaxies. We show that for GCs several anisotropic structures are observed. The GCs at small distances, $2<R<10$ kpc, show a structure coplanar with the Galactic plane. At smaller and larger distances the whole sample of GCs shows quite weak anisotropy. Nevertheless, at largest distances the orientation of the structure is close to that of the satellite galaxies, i.e. perpendicular to the Galactic plane. We estimate the probability of random realization for this structure of 1.7%. The Bulge-Disk GCs show a clear disk-like structure lying within the galactic disk. The Old Halo GCs show two structures: a well pronounced polar elongated structure at $R<3$ kpc which is perpendicular to the galactic plane, and a less pronounced disk-like structure coplanar with the galactic disk at $6<R<20$ kpc. The Young Halo GCs do not show significant anisotropy.

Using N-body simulations we study the impact of various systematic effects on the bulk flow (BF) and the Cosmic Mach Number (CMN). We consider two types of systematics: those related to survey properties and those induced by observer's location in the Universe. In the former category we model sparse sampling, velocity errors, and survey incompleteness. In the latter, we consider Local Group (LG) analogue observers, placed in a specific location within the Cosmic Web, satisfying various observational criteria. We differentiate such LG observers from Copernican ones, who are at random locations. We report strong systematic effects on the measured BF and CMN induced by sparse sampling, velocity errors and radial incompleteness. For BF most of these effects exceed 10\% for scales $R\leq100h^{-1}$Mpc. For CMN some of these systematics can be catastrophically large ($>50\%$) also on bigger scales. Moreover, we find that the position of the observer in the Cosmic Web significantly affects the locally measured BF (CMN), with effects as large as $\sim20\%$ ($30\%)$ at $R\leq50h^{-1}$Mpc for a LG-like observer as compared to a random one. This effect is comparable to the sample variance. To highlight the importance of these systematics, we additionally study a model of modified gravity (MG) with $\sim15\%$ enhanced growth rate. We found that the systematic effects can mimic the modified gravity signal. The worst-case scenario is realized for a case of a LG-like observer, when the effects induced by local structures are degenerate with the enhanced growth rate fostered by MG. Our results indicate that dedicated constrained simulations and realistic mock galaxy catalogs will be absolutely necessary to fully benefit from the statistical power of the forthcoming peculiar velocity data from surveys such as TAIPAN, WALLABY, Cosmic Flows-4 and SKA.

We revisit the estimation of the combined mass of the Milky Way and Andromeda (M31), which dominate the mass of the Local Group. We make use of an ensemble of 30,190 halo pairs from the Small MultiDark simulation, assuming a $\Lambda$CDM (Cosmological Constant with Cold Dark Matter) cosmology, to investigate the relationship between the bound mass and parameters characterising the orbit of the binary and their local environment with the aid of machine learning methods (artificial neural networks, ANN). Results from the ANN are most successful when information about the velocity shear is provided, which demonstrates the flexibility of machine learning to model physical phenomena and readily incorporate new information as it becomes available. The resulting estimate for the Local Group mass, when shear information is included, is $4.9 \times 10^{12} M_\odot$, with an error of $\pm0.8 \times 10^{12} M_\odot$ from the 68% uncertainty in observables, and a 68% confidence interval of $^{+1.3}_{-1.4} \times 10^{12}M_\odot$ from the intrinsic scatter from the differences between the model and simulation masses. We also consider a recently reported large transverse velocity of M31 relative to the Milky Way, and produce an alternative mass estimate of $3.6\pm0.3\pm1.4 \times 10^{12}M_\odot$. Although different methods predict similar values for the most likely mass of the LG, application of ANN compared to the Timing Argument reduces the scatter in the log mass by over half when tested on samples from the simulation.

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.

The influence of the large scale structure on host halos may be studied by examining the angular infall pattern of subhalos. In particular, since warm and cold dark matter cosmologies predict different abundances and internal properties for halos at the low mass end of the mass function, it is interesting to examine if there are differences in how these low mass halos are accreted. The accretion events are defined as the moment a halo becomes a substructure, namely when it crosses its host's virial radius. We quantify the cosmic web at each point by the shear tensor and examine where, with respect to its eigenvectors, such accretion events occur in cold ($\Lambda$CDM) and warm (1keV sterile neutrino WDM) dark matter cosmological models. We find that the CDM and WDM subhalos are preferentially accreted along the principal axis of the shear tensor corresponding to the direction of weakest collapse. The beaming strength is modulated by the host and subhalo masses and by the redshift at which the accretion event occurs. Although strongest for the most massive hosts and subhalos at high redshift, the preferential infall is found to be always aligned with the axis of weakest collapse, thus we say that it has universal nature. We compare the strength of beaming in the WDM cosmology with the one found in the $\Lambda$CDM scenario. While the main findings remain the same, the accretion in the WDM model for the most massive host halos appears more beamed than in $\Lambda$CDM cosmology across all the redshifts.

The existence of a spatially thin, kinematically coherent Disk of Satellites (DoS) around the Milky Way (MW), is a problem that often garners vivacious debate in the literature or at scientific meetings. One of the most recent incarnations of this wrangle occurred with two papers by Maji et al, who argued that these structures "maybe a misinterpretation of the data". These claims are in stark contrast to previous works. Motivated by this and other recent publications on this problem, we discuss necessary considerations to make, observational effects to consider, and pitfalls to avoid when investigating satellite galaxy planes such as the MW's DoS. In particular, we emphasize that conclusions need to have a statistical basis including a determination of the significance of satellite alignments, observational biases must not be ignored, and measurement errors (e.g. for proper motions) need to be considered. We discuss general problems faced by attempts to determine the dynamical stability of the DoS via orbit integrations of MW satellite galaxies, and demonstrate that to interpret simulations, it is helpful to compare them with a null case of isotropically distributed satellite positions and velocities. Based on these criteria, we find that the conclusions of Maji et al. do not hold up to scrutiny, and that their hydrodynamic cosmological simulation of a single host shows no evidence for a significant kinematic coherence among the simulated satellite galaxies, in contrast to the observed MW system.

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

Galaxy groups and clusters are the main tools used to test cosmological models and to study the environmental effect of galaxy formation. This work provides a catalogue of galaxy groups and clusters, as well as potentially merging systems based on the SDSS main galaxy survey. We identified galaxy groups and clusters using the modified friends-of-friends (FoF) group finder designed specifically for flux-limited galaxy surveys. The FoF group membership is refined by multimodality analysis to find subgroups and by using the group virial radius and escape velocity to expose unbound galaxies. We look for merging systems by comparing distances between group centres with group radii. The analysis results in a catalogue of 88662 galaxy groups with at least two members. Among them are 6873 systems with at least six members which we consider to be more reliable groups. We find 498 group mergers with up to six groups. We performed a brief comparison with some known clusters in the nearby Universe, including the Coma cluster and Abell 1750. The Coma cluster in our catalogue is a merging system with six distinguishable subcomponents. In the case of Abell 1750 we find a clear sign of filamentary infall toward this cluster. Our analysis of mass-to-light ratio (M/L) of galaxy groups reveals that M/L slightly increases with group richness.

Observations of the Milky Way (MW), M31, and their vicinity, known as the Local Group (LG), can provide clues about the sources of reionization. We present a suite of radiative transfer simulations based on initial conditions provided by the Constrained Local UniversE Simulations (CLUES) project that are designed to recreate the Local Universe, including a realistic MW-M31 pair and a nearby Virgo. Our box size (91 Mpc) is large enough to incorporate the relevant sources of ionizing photons for the LG. We employ a range of source models, mimicking the potential effects of radiative feedback for dark matter haloes between $10^{8}-10^{9}$ M$_{\odot}$. Although the LG mostly reionizes in an inside-out fashion, the final 40 per cent of its ionization shows some outside influence. For the LG satellites, we find no evidence that their redshift of reionization is related to the present-day mass of the satellite or the distance from the central galaxy. We find that less than 20 per cent of present-day satellites for MW and M31 have undergone any star formation prior to the end of global reionization. Approximately five per cent of these satellites could be classified as fossils, meaning the majority of star formation occurred at these early times. The more massive satellites have more cumulative star formation prior to the end of global reionization, but the scatter is significant, especially at the low-mass end. Present-day mass and distance from the central galaxy are poor predictors for the presence of ancient stellar populations in satellite galaxies.

The distribution of smaller satellite galaxies around large central galaxies has attracted attention because peculiar spatial and kinematic configurations have been detected in some systems. A particularly striking example of such behavior is seen in the satellite system of the Andromeda galaxy, where around 80\% are on the nearside of that galaxy, facing the Milky Way. Motivated by this departure from anisotropy, we examined the spatial distribution of satellites around pairs of galaxies in the SDSS. By stacking tens of thousands of satellites around galaxy pairs we found that satellites tend to bulge towards the other central galaxy, preferably occupying the space between the pair, rather than being spherically or axis-symmetrically distributed around each host. The bulging is a function of the opening angle examined and is fairly strong -- there are up to $\sim$10\% more satellites in the space between the pair, than expected from uniform. Consequently, it is a statistically very strong signal, being inconsistent with a uniform distribution at the 5$\sigma$ level. The possibility that the observed signal is the result of the overlap of two haloes with extended satellite distributions, is ruled out by testing this hypothesis by performing the same tests on isolated galaxies (and their satellites) artificially placed at similar separations. These findings highlight the unrelaxed and interacting nature of galaxies in pairs.

An open question in cosmology and the theory of structure formation is to what extent does environment affect the properties of galaxies and haloes. The present paper aims at shedding light on this problem. The paper focuses on the analysis of a dark matter only simulation and it addresses the issue of how the environment affects the abundance of haloes, which are are assigned four attributes: their virial mass, an ambient density calculated with an aperture that scales with $R_{vir}$ ($\Delta_M$), a fixed-aperture ($\Delta_R$) ambient density, and a cosmic web classification (i.e. voids, sheets, filaments, and knots, as defined by the V--web algorithm). $\Delta_M$ is the mean density around a halo evaluated within a sphere of a radius of $5$\rvir, where \rvir\ is the virial radius. $\Delta_R$ is the density field Gaussian smoothed with $R=4h^{-1}$Mpc, evaluated at the center of the halo. The main result of the paper is that the difference between haloes in different web elements stems from the difference in their mass functions, and does not depend on their adaptive-aperture ambient density. A dependence on the fixed-aperture ambient density is induced by the cross correlation between the mass of a halo and its fixed-aperture ambient density.

Near field cosmology is practiced by studying the Local Group (LG) and its neighbourhood. The present paper describes a framework for simulating the near field on the computer. Assuming the LCDM model as a prior and applying the Bayesian tools of the Wiener filter (WF) and constrained realizations of Gaussian fields to the Cosmicflows-2 (CF2) survey of peculiar velocities, constrained simulations of our cosmic environment are performed. The aim of these simulations is to reproduce the LG and its local environment. Our main result is that the LG is likely a robust outcome of the LCDM scenario when subjected to the constraint derived from CF2 data, emerging in an environment akin to the observed one. Three levels of criteria are used to define the simulated LGs. At the base level, pairs of halos must obey specific isolation, mass and separation criteria. At the second level the orbital angular momentum and energy are constrained and on the third one the phase of the orbit is constrained. Out of the 300 constrained simulations 146 LGs obey the first set of criteria, 51 the second and 6 the third. The robustness of our LG factory enables the construction of a large ensemble of simulated LGs. Suitable candidates for high resolution hydrodynamical simulations of the LG can be drawn from this ensemble, which can be used to perform comprehensive studies of the formation of the LG

Tidal torque theory suggests that galaxies gain angular momentum in the linear stage of structure formation. Such a theory predicts alignments between the spin of haloes and tidal shear field. However, non-linear evolution and angular momentum acquisition may alter this prediction significantly. In this paper, we use a reconstruction of the cosmic shear field from observed peculiar velocities combined with spin axes extracted from galaxies within $115\, \mathrm{Mpc} $ ($\sim8000 \, {\mathrm {km}}{\mathrm s}^{-1}$) from 2MRS catalog, to test whether or not galaxies appear aligned with principal axes of shear field. Although linear reconstructions of the tidal field have looked at similar issues, this is the first such study to examine galaxy alignments with velocity-shear field. Ellipticals in the 2MRS sample, show a statistically significant alignment with two of the principal axes of the shear field. In general, elliptical galaxies have their short axis aligned with the axis of greatest compression and perpendicular to the axis of slowest compression. Spiral galaxies show no signal. Such an alignment is significantly strengthened when considering only those galaxies that are used in velocity field reconstruction. When examining such a subsample, a weak alignment with the axis of greatest compression emerges for spiral galaxies as well. This result indicates that although velocity field reconstructions still rely on fairly noisy and sparse data, the underlying alignment with shear field is strong enough to be visible even when small numbers of galaxies are considered - especially if those galaxies are used as constraints in the reconstruction.

We present the first model that couples high-resolution simulations of the formation of Local Group galaxies with calculations of the galactic habitable zone (GHZ), a region of space which has sufficient metallicity to form terrestrial planets without being subject to hazardous radiation. These simulations allow us to make substantial progress in mapping out the asymmetric three-dimensional GHZ and its time evolution for the Milky Way (MW) and Triangulum (M33) galaxies, as opposed to works that generally assume an azimuthally symmetric GHZ. Applying typical habitability metrics to MW and M33, we find that while a large number of habitable planets exist as close as a few kiloparsecs from the galactic centre, the probability of individual planetary systems being habitable rises as one approaches the edge of the stellar disc. Tidal streams and satellite galaxies also appear to be fertile grounds for habitable planet formation. In short, we find that both galaxies arrive at similar GHZs by different evolutionary paths, as measured by the first and third quartiles of surviving biospheres. For the Milky Way, this interquartile range begins as a narrow band at large radii, expanding to encompass much of the galaxy at intermediate times before settling at a range of 2-13kpc. In the case of M33, the opposite behaviour occurs - the initial and final interquartile ranges are quite similar, showing gradual evolution. This suggests that galaxy assembly history strongly influences the time evolution of the GHZ, which will affect the relative time lag between biospheres in different galactic locations. We end by noting the caveats involved in such studies and demonstrate that high resolution cosmological simulations will play a vital role in understanding habitability on galactic scales, provided that these simulations accurately resolve chemical evolution.

Intermediate surface brightness (ISB) galaxies are less numerous than their counterparts at high and low surface brightness (HSB and LSB). Investigating ISB characteristics from a sample from the S4G survey, complete down to M_B=-16, we find that they have intermediate stellar, gas and baryonic masses and on average as much gas as stars. They lie on the (baryonic) Tully-Fisher relation between HSBs and LSBs, although they present a higher scatter than the latter. Their stellar to baryonic mass ratios have intermediate values unlike their condensed baryonic fractions. By comparing their environments, as classified by the eigenvalues of the velocity shear tensor of local constrained simulations, ISBs have a 5-10% probability higher (smaller) to be in sheets (filaments) with respect to HSBs and LSBs. Additionally, for galaxies in filaments (with close neighbors), the mass and mu_0 are correlated at 2.5 (2) sigma more than for those in sheets. ISBs live in regions where the divergence of the velocity field is smaller than where HSBs and LSBs live, a result at more than 50% significance. ISBs may exist as an unstable transition state between LSBs and HSBs, the low flow activity environment maximally encouraging their formation. Interaction events altering the central baryon fraction could happen at a lower rate in these less dense environment, whilst in the higher density environments the LSBs are primarily satellite galaxies, whose accretion is sufficiently constrained that it fails to promote them to HSBs.

Recent observational studies have demonstrated that the majority of satellite galaxies tend to orbit their hosts on highly flattened, vast, possibly co-rotating planes. Two nearly parallel planes of satellites have been confirmed around the M31 galaxy and around the Centaurus A galaxy, while the Milky Way also sports a plane of satellites. It has been argued that such an alignment of satellites on vast planes is unexpected in the standard ({\Lambda}CDM) model of cosmology if not even in contradiction to its generic predictions. Guided by {\Lambda}CDM numerical simulations, which suggest that satellites are channeled towards hosts along the axis of the slowest collapse as dictated by the ambient velocity shear tensor, we re-examine the planes of local satellites systems within the framework of the local shear tensor derived from the Cosmicflows-2 dataset. The analysis reveals that the Local Group and Centaurus A reside in a filament stretched by the Virgo cluster and compressed by the expansion of the Local Void. Four out of five thin planes of satellite galaxies are indeed closely aligned with the axis of compression induced by the Local Void. Being the less massive system, the moderate misalignment of the Milky Way's satellite plane can likely be ascribed to its greater susceptibility to tidal torques, as suggested by numerical simulations. The alignment of satellite systems in the local universe with the ambient shear field is thus in general agreement with predictions of the {\Lambda}CDM model.

The cosmic web that characterizes the large-scale structure of the Universe can be quantified by a variety of methods. For example, large redshift surveys can be used in combination with point process algorithms to extract long curvilinear filaments in the galaxy distribution. Alternatively, given a full 3D reconstruction of the velocity field, kinematic techniques can be used to decompose the web into voids, sheets, filaments and knots. In this paper we look at how two such algorithms - the Bisous model and the velocity shear web - compare with each other in the local Universe (within 100 Mpc), finding good agreement. This is both remarkable and comforting, given that the two methods are radically different in ideology and applied to completely independent and different data sets. Unsurprisingly, the methods are in better agreement when applied to unbiased and complete data sets, like cosmological simulations, than when applied to observational samples. We conclude that more observational data is needed to improve on these methods, but that both methods are most likely properly tracing the underlying distribution of matter in the Universe.

The accretion of satellites onto central galaxies along vast cosmic filaments is an apparent outcome of the anisotropic collapse of structure in our Universe. Numerical work (based on gravitational dynamics of N-body simulations) indicates that satellites are beamed towards hosts along preferred directions imprinted by the velocity shear field. Here we use the Sloan Digital Sky Survey to observationally test this claim. We construct 3D filaments and sheets and examine the relative position of satellite galaxies. A statistically significant alignment between satellite galaxy position and filament axis in observations is confirmed. We find a qualitatively compatible alignments by examining satellites and filaments similarly identified in the Millennium simulation, semi-analytical galaxy catalogue. We also examine the dependence of the alignment strength on galaxy properties such as colour, magnitude and (relative) satellite magnitude, finding that the alignment is strongest for the reddest and brightest central and satellite galaxies. Our results confirm the theoretical picture and the role of the cosmic web in satellite accretion. Furthermore our results suggest that filaments identified on larger scales can be reflected in the positions of satellite galaxies that are quite close to their hosts.

Tip of the red giant branch measurements based on Hubble Space Telescope and ground-based imaging have resulted in accurate distances to 29 galaxies in the nearby Centaurus A Group. All but two of the 29 galaxies lie in either of two thin planes roughly parallel with the supergalactic equator. The planes are only slightly tilted from the line-of-sight, leaving little ambiguity regarding the morphology of the structure. The planes have characteristic r.m.s. long axis dimensions of ~300 kpc and short axis dimensions of ~60 kpc, hence axial ratios ~0.2, and are separated in the short axis direction by 303 kpc.

We investigate whether the satellite luminosity function (LF) of primary galaxies identified in the Sloan Digital Sky Survey (SDSS) depends on whether the host galaxy is in a filament or not. Isolated primary galaxies are identified in the SDSS spectroscopic sample while potential satellites (that are up to 4 magnitudes fainter than their hosts) are searched for in the much deeper photometric sample. Filaments are constructed from the galaxy distribution by the "Bisous" process. Isolated primary galaxies are divided into two subsamples: those in filaments and those not in filaments. We examine the stacked mean satellite LF of both the filament and non-filament sample and find that, on average, the satellite LFs of galaxies in filaments is significantly higher than those of galaxies not in filaments. The filamentary environment can increases the abundance of the brightest satellites ($M_\mathrm{sat.} < M_\mathrm{prim.} + 2.0$), by a factor of $\sim 2$ compared with non-filament isolated galaxies. This result is independent of primary galaxy magnitude although the satellite LF of galaxies in the faintest magnitude bin, is too noisy to determine if such a dependence exists. Since our filaments are extracted from a spectroscopic flux-limited sample, we consider the possibility that the difference in satellite LF is due to a redshift, colour or environmental bias, finding these to be insufficient to explain our result. The dependence of the satellite LF on the cosmic web suggests that the filamentary environment may have a strong effect on the efficiency of galaxy formation.

We search for vast planes of satellites (VPoS) in a high resolution simulation of the Local Group performed by the CLUES project, which improves significantly the resolution of former similar studies. We use a simple method for detecting planar configurations of satellites, and validate it on the known plane of M31. We implement a range of prescriptions for modelling the satellite populations, roughly reproducing the variety of recipes used in the literature, and investigate the occurence and properties of planar structures in these populations. The structure of the simulated satellite systems is strongly non-random and contains planes of satellites, predominantly co-rotating, with, in some cases, sizes comparable to the plane observed in M31 by Ibata et al.. However the latter is slightly richer in satellites, slightly thinner and has stronger co-rotation, which makes it stand out as overall more exceptional than the simulated planes, when compared to a random population. Although the simulated planes we find are generally dominated by one real structure, forming its backbone, they are also partly fortuitous and are thus not kinematically coherent structures as a whole. Provided that the simulated and observed planes of satellites are indeed of the same nature, our results suggest that the VPoS of M31 is not a coherent disc and that one third to one half of its satellites must have large proper motions perpendicular to the plane.

We seek to understand the relationship between galaxy properties and their local environment, which calls for a proper formulation of the notion of environment. We analyse the GIMIC suite of cosmological hydrodynamical simulations within the framework of the cosmic web as formulated by Hoffman et al., focusing on properties of simulated DM haloes and luminous galaxies with respect to voids, sheets, filaments and knots - the four elements of the cosmic web. We find that the mass functions of haloes depend on environment, which drives other environmental dependence of galaxy formation. The web shapes the halo mass function, and through the strong dependence of the galaxy properties on the mass of their host haloes, it also shapes the galaxy-(web) environment dependence.

We use high resolution simulations of the formation of the local group post-processed by a radiative transfer code for UV photons, to investigate the reionization of the satellite populations of an isolated Milky Way-M31 galaxy pair in a variety of scenarios. We use an improved version of ATON which includes a simple recipe for radiative feedback. In our baseline models, reionization is initiated by low mass, radiatively regulated haloes at high redshift, until more massive haloes appear, which then dominate and complete the reionization process. We investigate the relation between reionization history and present-day positions of the satellite population. We find that the average reionization redshift (zr) of satellites is higher near galaxy centers (MW and M31). This is due to the inside-out reionization patterns imprinted by massive haloes within the progenitor during the EoR, which end up forming the center of the galaxy. Thanks to incomplete dynamical mixing during galaxy assembly, these early patterns survive down to present day, resulting in a a clear radial gradient in the average satellites reionization redshift, up to the virial radius of MW and M31 and beyond. In the lowest emissivity scenario, the outer satellites are reionized about 180 Myr later than the inner satellites. This delay decreases with increasing source model emissivity, or in the case of external reionization by Virgo or M31, because reionization happens faster overall, and becomes spatially quasi-uniform at the highest emissivity.

We present, for the first time, a Local Universe (LU) characterization using high precision constrained $N$-body simulations based on self-consistent phase-space reconstructions of the large-scale structure in the Two-Micron All-Sky Galaxy Redshift Survey. We analyse whether we live in a special cosmic web environment by estimating cosmic variance from a set of unconstrained $\Lambda$CDM simulations as a function of distance to random observers. By computing volume and mass filling fractions for voids, sheets, filaments and knots, we find that the LU displays a typical scatter of about $1\sigma$ at scales $r>15\,h^{-1}\,{\rm Mpc}$, in agreement with $\Lambda$CDM, converging to a fair unbiased sample when considering spheres of about $60\,h^{-1}\,{\rm Mpc}$ radius. Additionally, we compute the matter density profile of the LU and found a reasonable agreement with the estimates of Karachentsev (2012) only when considering the contribution of dark haloes. This indicates that observational estimates may be biased towards low density values. As a first application of our reconstruction, we investigate the likelihood of different galaxy morphological types to inhabit certain cosmic web environments. In particular, we find that, irrespective of the method used to define the web, either based on the density or the peculiar velocity field, elliptical galaxies show a clear tendency to preferentially reside in clusters as opposed to voids (up to a level of $5.3\sigma$ and $9.8\sigma$ respectively) and conversely for spiral galaxies (up to a level of $5.6\sigma$ and $5.4\sigma$ respectively). These findings are compatible with previous works, albeit at higher confidence levels.

We examine the angular infall pattern of subhaloes onto host haloes in the context of the large-scale structure. We find that this infall pattern is essentially driven by the shear tensor of the ambient velocity field. Dark matter subhaloes are found to be preferentially accreted along the principal axis of the shear tensor which corresponds to the direction of weakest collapse. We examine the dependence of this preferential infall on subhalo mass, host halo mass and redshift. Although strongest for the most massive hosts and the most massive subhaloes at high redshift, the preferential infall of subhaloes is effectively universal in the sense that its always aligned with the axis of weakest collapse of the velocity shear tensor. It is the same shear tensor that dictates the structure of the cosmic web and hence the shear field emerges as the key factor that governs the local anisotropic pattern of structure formation. Since the small (sub-Mpc) scale is strongly correlated with the mid-range (~10 Mpc) scale - a scale accessible by current surveys of peculiar velocities - it follows that findings presented here open a new window into the relation between the observed large scale structure unveiled by current surveys of peculiar velocities and the preferential infall direction of the Local Group. This may shed light on the unexpected alignments of dwarf galaxies seen in the Local Group

The evolution of the large scale distribution of matter in the universe is often characterized by the density field. Here we take a complimentary approach and characterize it using the cosmic velocity field, specifically the deformation of the velocity field. The deformation tensor is decomposed into its symmetric component (known as the "shear tensor") and its anti-symmetric part (the "vorticity"). Using a high resolution cosmological simulation we examine the relative orientations of the shear and the vorticity as a function of spatial scale and redshift. The shear is found to be remarkable stable to the choice of scale, while the vorticity is found to quickly decay with increasing spatial scale or redshift. The vorticity emerges out of the linear regime randomly oriented with respect to the shear eigenvectors. Non-linear evolution drives the vorticity to lie within the plane defined by the eigenvector of the fastest collapse. Within that plane the vorticity first gets aligned with the middle eigenvector and then it moves to be preferentially aligned with the third eigenvector, of slowest collapse. Finally a scale of "non-linearity" to be used when calculating properties of the non-linear deformation tensor at different redshifts is suggested.

We present the results of a set of numerical simulations aimed at studying reionization at galactic scale. We use a high resolution simulation of the formation of the Milky Way-M31 system to simulate the reionization of the local group. The reionization calculation was performed with the post-processing radiative transfer code ATON and the underlying cosmological simulation was performed as part of the CLUES project. We vary the source models to bracket the range of source properties used in the literature. We investigate the structure and propagation of the galatic ionization fronts by a visual examination of our reionization maps. Within the progenitors we find that reionization is patchy, and proceeds locally inside out. The process becomes patchier with decreasing source photon output. It is generally dominated by one major HII region and 1-4 additional isolated smaller bubbles, which eventually overlap. Higher emissivity results in faster and earlier local reionization. In all models, the reionization of the Milky Way and M31 are similar in duration, i.e. between 203 Myr and 22 Myr depending on the source model, placing their zreion between 8.4 and 13.7. In all models except the most extreme, the MW and M31 progenitors reionize internally, ignoring each other, despite being relatively close to each other even during the epoch of reionization. Only in the case of strong supernova feedback suppressing star formation in haloes less massive than 10^9 M_sun, and using our highest emissivity, we find that the MW is reionized by M31.

The large-scale structure of the Universe is characterised by a web-like structure made of voids, sheets, filaments, and knots. The structure of this so-called cosmic web is dictated by the local velocity shear tensor. In particular, the local direction of a filament should be strongly aligned with e3, the eigenvector associated with the smallest eigenvalue of the tensor. That conjecture is tested here on the basis of a cosmological simulation. The cosmic web delineated by the halo distribution is probed by a marked point process with interactions (the Bisous model), detecting filaments directly from the halo distribution (P-web). The detected P-web filaments are found to be strongly aligned with the local e3: the alignment is within 30 degree for 80% of the elements. This indicates that large-scale filaments defined purely from the distribution of haloes carry more than just morphological information, although the Bisous model does not make any prior assumption on the underlying shear tensor. The P-web filaments are also compared to the structure revealed from the velocity shear tensor itself (V-web). In the densest regions, the P- and V-web filaments overlap well (90%), whereas in lower density regions, the P-web filaments preferentially mark sheets in the V-web.

Galaxies properties are known to be affected by their environment. One important question is how their angular momentum reflects the surrounding cosmic web. We use the SDSS to investigate the spin axes of spiral and elliptical galaxies relative to their surrounding filament/sheet orientations. To detect filaments a marked point process with interactions (the "Bisous model") is used. Sheets are found by detecting "flattened" filaments. The minor axes of ellipticals are found to be preferentially perpendicular to hosting filaments. A weak correlation is found with sheets. These findings are consistent with the notion that elliptical galaxies formed via mergers which predominantly occurred along the filaments. The spin axis of spiral galaxies is found to align with the host filament, with no correlation between spiral spin and sheet normal. When examined as a function of distance from the filament axis, a much stronger correlation is found in outer parts, suggesting that the alignment is driven by the laminar infall of gas from sheets to filaments. When compared with numerical simulations our results suggest that the connection between DM halo and galaxy spin is not straightforward. Our results provide an important input to the understanding of how galaxies acquire their angular momentum.

With the ever increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations has simultaneously grown. To this extent we initiated the Halo Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes (Knebe et al. 2011), and one for dark matter subhaloes (Onions et al. 2012). However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also account for baryonic physics (e.g. hydrodynamics, star formation, feedback) object finders should also be capable of taking these additional processes into consideration. Here we report on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the MW, M31, and M33) as well as their satellite populations for a variety of halo finders ranging from phase-space to velocity-space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons, at least for the total, dark, and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.

The differences between cold (CDM) and warm (WDM) dark matter in the formation of a group of galaxies is examined by running two identical simulations where in the WDM case the initial power spectrum has been altered to mimic a 1keV dark matter particle. The CDM initial conditions were constrained to reproduce at z = 0 the correct local environment within which a "Local Group" (LG) of galaxies may form. Two significant differences between the two simulations are found. While in the CDM case a group of galaxies that resembles the real LG forms, the WDM run fails to reproduce a viable LG, instead forming a diffuse group which is still expanding at z = 0. This is surprising since, due to the suppression of small scale power in its power spectrum, WDM is naively expected to only affect the collapse of small haloes and not necessarily the dynamics on a scale of a group of galaxies. Furthermore the concentration of baryons in halo center's is greater in CDM than in WDM and the properties of the disks differ.

We use state-of-the-art numerical simulations to explore the observability and the expected physical properties of the progenitors of the local group galaxies at z=6-8, within 1 billion years of the big bang. We find that the most massive progenitors of the Milky Way (MW) and Andromeda (M31) at z=6 and 7 are predicted to have absolute UV continuum magnitudes between -17 and -18, suggesting that their analogues lie close to the detection limits of the deepest near-infrared surveys conducted to date (i.e. HST WFC3/IR UDF12). This in turn confirms that the majority of currently known z=6-8 galaxies are expected to be the seeds of present-day galaxies which are more massive than L* spirals. We also discuss the properties of the local-group progenitors at these early epochs, extending down to absolute magnitudes M_UV = -13. The most massive MW/M31 progenitors at z=7 have stellar masses of about 10^7.5-8 solar masses, stellar metallicities between 3 and 6% of the solar value, and predicted observed UV continuum slopes between -2.4 and -2.5.

We use dark matter only and full hydrodynamical Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group to study the density profile of subhaloes of the simulated Milky Way and Andromeda galaxies. We show that the Einasto model provides the best description of the subhaloes' density profile, as opposed to the more commonly used NFW profile or any generalisation of it. We further find that the Einasto shape parameter \nEin\ is strongly correlated with the total subhalo mass, pointing towards the notion of a non-universality of the subhaloes' density profile. We observe that the effect of mass loss due to tidal stripping, in both the dark matter only and the hydrodynamical run, is the reduction of the shape parameter \nEin\ between the infall and the present time. Assuming now that the dSphs of our Galaxy follow the Einasto profile and using the maximum and minimum values of \nEin\ from our hydrodynamical simulation as a gauge, we can improve the observational constraints on the \Rmax-\Vmax\ pairs obtained for the brightest satellite galaxies of the Milky Way. When considering only the subhaloes with $-13.2\lesssim M_V\lesssim-8.8$, i.e. the range of luminosity of the classical dwarfs, we find that all our simulated objects are consistent with the observed dSphs if their haloes follow the Einasto model with $1.6\lesssim n_{\rm E} \lesssim5.3$. The numerically motivated Einasto profile for the observed dSphs will alleviate the recently presented "massive failures" problem.

In the standard model of cosmology, structure emerges out of non-rotational flow and the angular momentum of collapsing halos is induced by tidal torques. The growth of halo angular momentum in the linear and quasi-linear phases is associated with a shear, curl-free, flow and it is well described within the linear framework of tidal torque theory (TTT). However, TTT is rendered irrelevant as haloes approach turn around and virialization. At that stage the flow field around halos has non-zero vorticity. Using a cosmological simulation, we have examined the importance of the curl of the velocity field (vorticity) in determining halo spin, finding a strong alignment between the two. We have also examined the alignment of vorticity with the principle axes of the shear tensor, finding that it tends to be perpendicular to the axis along which material is collapsing fastest (e1). This behavior is independent of halo masses and cosmic web environment. Our results agree with previous findings on the tendency of halo spin to be perpendicular to e1, and of the spin of (simulated) halos and (observed) galaxies to be aligned with the large-scale structure. Our results imply that angular momentum growth proceeds in two distinct phases. In the first phase angular momentum emerges out of a shear, curl-free, potential flow, as described by TTT. In the second phase, in which haloes approach virialization, the angular momentum emerges out of a vortical flow and halo spin becomes strongly aligned with the vorticity of the ambient flow field.

The alignment of DM halos and the surrounding large scale structure (LSS) is examined in the context of the cosmic web. Halo spin, shape and the orbital angular momentum of subhaloes is investigated relative to the LSS using the eigenvectors of the velocity shear tensor evaluated on a grid with a scale of 1 Mpc/h, deep within the non-linear regime. Knots, filaments, sheets and voids are associated with regions that are collapsing along 3, 2, 1 or 0 principal directions simultaneously. Each halo is tagged with a web classification (i.e. knot halo, filament halo, etc) according to the nature of the collapse at the halo's position. The full distribution of shear eigenvalues is found to be substantially different from that tagged to haloes, indicating that the observed velocity shear is significantly biased. We find that larger mass haloes live in regions where the shear is more isotropic, namely the expansion or collapse is more spherical. A correlation is found between the halo's shape and the eigenvectors of the shear tensor, with the longest (shortest) axis of the halo's shape being aligned with the slowest (fastest) collapsing eigenvector. This correlation is web independent, suggesting that the velocity shear is a fundamental tracer of the halo alignment. A similar result is found for the alignment of halo spin with the cosmic web. It has been shown that high mass haloes exhibit a spin flip with respect to the LSS: we find the mass at which this spin flip occurs is web dependent and not universal as suggested previously. Although weaker than haloes, subhalo orbits too exhibit an alignment with the LSS, providing a possible insight into the highly correlated co-rotation of the Milky Way's satellite system. The present study suggests that the velocity shear tensor constitutes the natural framework for studying the directional properties of the non-linear LSS and of halos and galaxies.

A new approach for the classification of the cosmic web is presented. In extension of the previous work of Hahn et al. (2007) and Forero-Romero et al. (2009) the new algorithm is based on the analysis of the velocity shear tensor rather than the gravitational tidal tensor. The procedure consists of the construction of the the shear tensor at each (grid) point in space and the evaluation of its three eigenvectors. A given point is classified to be either a void, sheet, filament or a knot according to the number of eigenvalues above a certain threshold, 0, 1, 2, or 3 respectively. The threshold is treated as a free parameter that defines the web. The algorithm has been applied to a dark matter only, high resolution simulation of a box of side-length 64$h^{-1}$Mpc and N = $1024^3$ particles with the framework of the WMAP5/LCDM model. The resulting velocity based cosmic web resolves structures down to <0.1$h^{-1}$Mpc scales, as opposed to the ~1$h^{-1}$Mpc scale of the tidal based web. The under-dense regions are made of extended voids bisected by planar sheets, whose density is also below the mean. The over-dense regions are vastly dominated by the linear filaments and knots. The resolution achieved by the velocity based cosmic web provides a platform for studying the formation of halos and galaxies within the framework of the cosmic web.

We use the recently proposed scale-free mass estimators to determine the masses of the Milky Way (MW) and Andromeda (M31) galaxy in a dark matter only Constrained Local UniversE Simulation (CLUES). While these mass estimators work rather well for isolated spherical host systems, we examine here their applicability to a simulated binary system with a unique satellite population similar to the observed satellites of MW and M31. We confirm that the scale-free estimators work also very well in our simulated Local Group galaxies with the right number of satellites which follow the observed radial distribution. In the isotropic case and under the assumption that the satellites are tracking the total gravitating mass, the power-law index of the radial satellite distribution $N(<r)\propto r^{3-\gamma}$ is directly related to the host's mass profile $M(<r)\propto r^{1-\alpha}$ as $\alpha=\gamma-2$. The use of this relation for any given $\gamma$ leads to highly accurate mass estimations which is a crucial point for observer, since they do not know a priori the mass profile of the MW and M31 haloes. We discuss possible bias in the mass estimators and conclude that the scale-free mass estimators can be satisfactorily applied to the real MW and M31 system.

We use a 64$h^{-1}$Mpc dark matter (DM) only cosmological simulation to examine the large scale orientation of haloes and substructures with respect the cosmic web. A web classification scheme based on the velocity shear tensor is used to assign to each halo in the simulation a web type: knot, filament, sheet or void. Using $\sim10^6$ haloes that span ~3 orders of magnitude in mass the orientation of the halo's spin and the orbital angular momentum of subhaloes with respect to the eigenvectors of the shear tensor is examined. We find that the orbital angular momentum of subhaloes tends to align with the intermediate eigenvector of the velocity shear tensor for all haloes in knots, filaments and sheets. This result indicates that the kinematics of substructures located deep within the virialized regions of a halo is determined by its infall which in turn is determined by the large scale velocity shear, a surprising result given the virilaized nature of haloes. The non-random nature of subhalo accretion is thus imprinted on the angular momentum measured at z = 0. We also find that haloes' spin axis is aligned with the third eigenvector of the velocity shear tensor in filaments and sheets: the halo spin axis points along filaments and lies in the plane of cosmic sheets.

We present a novel approach of identifying the Milky Way (MW) and Andromeda (M31) progenitors that could be visible as LAEs at z~6: we couple a snapshot from the Constrained Local UniversE Simulations (CLUES) project, that successfully reproduces the MW and M31 galaxies situated in their correct environment, to a Lyman Alpha Emitter (LAE) model. Exploring intergalactic medium (IGM) ionization states ranging from an almost neutral to a fully ionized one, we find that including (excluding) the effects of clustered sources the first local group progenitor appears as a LAE for a neutral hydrogen fraction of 0.4 (0.1). This number increases to 5 progenitors each of the MW and M31 being visible as LAEs for a completely ionized IGM; the contribution from clustered sources is crucial in making many of the progenitors visible in the Lyman Alpha, for all the ionization states considered. We also put constrains on the physical properties of such progenitors. Finally, we find that the number density of these LAEs are higher than that of general field LAEs (observed in cosmological volumes) by about two (one) orders of magnitude for a neutral hydrogen fraction of 10^-5 (0.4). Detections of such high LAE number densities at z ~ 6 would be a clear signature of an over-dense region that could evolve and resemble the local group volume at z=0.

Using Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group in a cosmological context we investigate the recently highlighted problem that the majority of the most massive dark subhaloes of the Milky Way are too dense to host any of its bright satellites. In particular, we examine the influence of baryonic processes and find that they leave a twofold effect on the relation between the peak of the rotation curve and its position (Vmax and Rmax). Satellites with a large baryon fraction experience adiabatic contraction thus decreasing Rmax while leaving Vmax more or less unchanged. Subhaloes with smaller baryon fractions undergo a decrease in Vmax possibly due to outflows of material. Furthermore, the situation of finding subhaloes in simulations that lie outside the confidence interval for possible hosts of the bright MW dwarf spheroidals, appears to be far more prominent in cosmologies with a high sigma8 normalisation and depends on the mass of the host. We conclude that the problem cannot be simply solved by including baryonic processes and hence demands further investigations.

The outer haloes of the Milky Way and Andromeda galaxies contain as much important information on their assembly and formation history as the properties of the discs resident in their centres. In this paper we have used the Constrained Local UniversE Simulation project to disentangle the stellar and DM component of three galaxies that resemble the MW, M31 and M33 using both DM only and DM + gas-dynamical simulations. Stars that are accreted in substructures and then stripped follow a completely different radial distribution than the stripped DM: the stellar halo is much more centrally concentrated than DM. In order to understand how the same physical process can lead to different z=0 radial profiles, we examined the potential at accretion of each stripped particle. We found that star particles sit at systematically higher potentials than DM, making them harder to strip. We then searched for a threshold in the potential of accreted particles \phi_[th], above which DM particles behave as star particles. We found such a threshold at >16 \phi_{subhalo}, where \phi_{subhalo} is the potential at a subhaloes edge at the time of accretion. Thus a rule as simple as selecting particles according to their potential at accretion is able to reproduce the effect that the complicated physics of star formation has on the stellar distribution. This niversal result reproduces the stellar halo to an accuracy of within ~2%. Studies which make use of DM particles as a proxy for stars will undoubtedly miscalculate their proper radial distribution and structure unless particles are selected according to their potential at accretion. Furthermore, we have examined the time it takes to strip a given star or DM particle after accretion. We find that, owing to their higher binding energies, stars take longer to be stripped than DM. The stripped DM halo is thus considerably older than the stripped stellar halo.

Using a dark matter only Constrained Local UniversE Simulation (CLUES) we examine the existence of subhaloes that change their affiliation from one of the two prominent hosts in the Local Group (i.e. the Milky Way and the Andromeda galaxy) to the other, and call these objects "renegade subhaloes". In light of recent claims that the two Magellanic Clouds (MCs) may have originated from another region (or even the outskirts) of the Local Group or that they have been spawned by a major merger in the past of the Andromeda galaxy, we investigate the nature of such events. However, we cannot confirm that renegade subhaloes enter as deep into the potential well of their present host nor that they share the most simplest properties with the MCs, namely mass and relative velocity. Our simulation rather suggests that these renegade subhaloes appear to be flying past one host before being pulled into the other. A merger is not required to trigger such an event, it is rather the distinct environment of our simulated Local Group facilitating such behavior. Since just a small fraction of the full z=0 subhalo population are renegades, our study indicates that it will be intrinsically difficult to distinguish them despite clear differences in their velocity, radial distribution, shape and spin parameter distributions.

(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0.

We study the differences and similarities in the luminosities of bound, infalling and the so-called backsplash (Gill et al. 2005) galaxies of the Milky Way and M31 using a hydrodynamical simulation performed within the Constrained Local UniversE Simulation (CLUES) project. The simulation models the formation of the Local Group within a self-consistent cosmological framework. We find that even though backsplash galaxies passed through the virial radius of their host halo and hence may have lost a (significant) fraction of their mass, their stellar populations are hardly affected. This leaves us with comparable luminosity functions for infalling and backsplash galaxies and hence little hope to decipher their past (and different) formation and evolutionary histories by luminosity measurements alone. Nevertheless, due to the tidal stripping of dark matter we find that the mass-to-light ratios have changed when comparing the various populations against each other: they are highest for the infalling galaxies and lowest for the bound satellites with the backsplash galaxies in-between.

Using a high resolution DM simulation of the Local Group (LG), conducted within the framework of the Constrained Local UniversE Simulation (CLUES) project, we investigate the nature of how satellites of the MW and M31 are accreted. Satellites of these galaxies are accreted anisotropically, entering the virial radius of their hosts from specific "spots" with respect to the large scale structure. Furthermore, material which is tidally stripped from these satellites is also, at z=0, distributed anisotropically and is characterized by an ellipsoidal sub-volume embedded in the halo. The angular pattern created by the locus of satellite infall points and the projected stripped DM is investigated within a coordinate system determined by the location of the LG companion and the simulated Virgo cluster across concentric shells ranging from 0.1 to 5 r_vir. Remarkably, the principal axis of the ellipsoidal sub-volume shows a coherent alignment extending from well within the halo to a few r_vir. A spherical harmonics transform applied to the angular distributions confirms the visual impression: namely, the angular distributions of both the satellites entry points and stripped DM is dominated by the l=2 quadrupole term, whose major principal axis is aligned across the shells considered. It follows that the structure of the outer (r>0.5 r_vir) regions of the main halos is closely related to the cosmic web within which it is embedded. Given the hypothesis that a large fraction of the stellar halo of the Milky Way has been accreted from satellites, the present results can be applied to the stellar halo of the MW and M31. We predict that the remnants of tidally stripped satellites should be embedded in streams of material composed of dark matter and stars. The present results can therefore shed light on the existence of satellites embedded within larger streams of matter, such as the Segue 2 satellite.

We use two simulations performed within the Constrained Local UniversE Simulation (CLUES) project to study both the shape and radial alignment of (the dark matter component of) subhaloes; one of the simulations is a dark matter only model while the other run includes all the relevant gas physics and star formation recipes. We find that the involvement of gas physics does not have a statistically significant effect on either property -- at least not for the most massive subhaloes considered in this study. However, we observe in both simulations including and excluding gasdynamics a (pronounced) evolution of the dark matter shapes of subhaloes as well as of the radial alignment signal since infall time. Further, this evolution is different when positioned in the central and outer regions of the host halo today; while subhaloes tend to become more aspherical in the central 50% of their host's virial radius, the radial alignment weakens in the central regime while strengthening in the outer parts. We confirm that this is due to tidal torquing and the fact that subhaloes at pericentre move too fast for the alignment signal to respond.

We examine the properties of satellites found in high resolution simulations of the local group. We use constrained simulations designed to reproduce the main dynamical features that characterize the local neighborhood, i.e. within tens of Mpc around the Local Group (LG). Specifically, a LG-like object is found located within the 'correct' dynamical environment and consisting of three main objects which are associated with the Milky Way, M31 and M33. By running two simulations of this LG from identical initial conditions - one with and one without baryons modeled hydrodynamically - we can quantify the effect of gas physics on the $z=0$ population of subhaloes in an environment similar to our own. We find that above a certain mass cut, $M_{\rm sub} > 2\times10^{8}h^{-1} M_{\odot}$ subhaloes in hydrodynamic simulations are more radially concentrated than those in simulations with out gas. This is caused by the collapse of baryons into stars that typically sit in the central regions of subhaloes, making them denser. The increased central density of such a subhalo, results in less mass loss due to tidal stripping than the same subhalo simulated with only dark matter. The increased mass in hydrodynamic subhaloes with respect to dark matter ones, causes dynamical friction to be more effective, dragging the subhalo towards the centre of the host. This results in these subhaloes being effectively more radially concentrated then their dark matter counterparts.

The highly flattened distribution of satellite galaxies in the Milky Way presents a number of puzzles. Firstly, its polar alignment stands out from the planar alignments commonly found in other galaxies. Secondly, recent proper motion measurements reveal that the orbital angular momentum of at least 3, and possibly as many as 8, of the Milky Ways satellites point (within 30 degrees) along the axis of their flattened configuration, suggesting some form of coherent motion. In this paper we use a high resolution cosmological simulation to investigate whether this pattern conflicts with the expectations of the cold dark matter model of structure formation. We find that this seemingly unlikely set up occurs often: approximately 35% of the time we find systems in which the angular momentum of 3 individual satellites point along, or close to, the short axis of the satellite distribution. In addition, in 30% of the systems we find that the net angular momentum of the 6 best aligned satellites lies within 35 degrees of the short axis of the satellite distribution, as observed for the Milky Way.

Available proper motion measurements of Milky Way (MW) satellite galaxies are used to calculate their orbital poles and projected uncertainties. These are compared to a set of recent cold dark-matter (CDM) simulations, tailored specifically to solve the MW satellite problem. We show that the CDM satellite orbital poles are fully consistent with being drawn from a random distribution, while the MW satellite orbital poles indicate that the disc-of-satellites of the Milky Way is rotationally supported. Furthermore, the bootstrapping analysis of the spatial distribution of theoretical CDM satellites also shows that they are consistent with being randomly drawn. The theoretical CDM satellite population thus shows a significantly different orbital and spatial distribution than the MW satellites, most probably indicating that the majority of the latter are of tidal origin rather than being DM dominated sub-structures. A statistic is presented that can be used to test a possible correlation of satellite galaxy orbits with their spatial distribution.

We investigate the properties of satellite galaxies formed in N-body/SPH simulations of galaxy formation in the $\Lambda$CDM cosmology. The simulations include the main physical effects thought to be important in galaxy formation and, in several cases, produce realistic spiral discs. In total, a sample of 9 galaxies of luminosity comparable to the Milky Way was obtained. At magnitudes brighter than the resolution limit, $M_V=-12$, the luminosity function of the satellite galaxies in the simulations is in excellent agreement with data for the Local Group. The radial number density profile of the model satellites, as well as their gas fractions also match observations very well. In agreement with previous N-body studies, we find that the satellites tend to be distributed in highly flattened configurations whose major axis is aligned with the major axis of the (generally triaxial) dark halo. In 2 out of 3 systems with sufficiently large satellite populations, the satellite system is nearly perpendicular to the plane of the galactic disc, a configuration analogous to that observed in the Milk Way. The discs themselves are perpendicular to the minor axis of their host halos in the inner parts, and the correlation between the orientation of the galaxy and the shape of the halo persists even out to the virial radius. However, in one case the disc's minor axis ends up, at the virial radius, perpendicular to the minor axis of the halo. The angular momenta of the galaxies and their host halo tend to be well aligned.

Galactic bulges are known to harbour central black holes whose mass is tightly correlated with the stellar mass and velocity dispersion of the bulge. In a hierarchical universe, mergers of subgalactic units are accompanied by the amalgamation of bulges and the likely coalescence of galactocentric black holes. In these mergers, the beaming of gravitational radiation during the plunge phase of the black hole collision can impart a linear momentum kick or ``gravitational recoil'' to the remnant. If large enough, this kick will eject the remnant from the galaxy and populate intergalactic space with wandering black holes. Using a semi-analytic model of galaxy formation, we investigate the effect of black hole ejections on the scatter in the relation between black hole and bulge mass. We find that although not the dominant source of the measured scatter, they do make a significant contribution and may be used to set a constraint, v_kick<500 km/s, on the typical kick velocity, in agreement with values found from general relativistic calculations. Even for the more modest kick velocities implied by these calculations, we find that a substantial number of central black holes are ejected from the progenitors of present day galaxies, giving rise to a population of wandering intrahalo and intergalactic black holes whose distribution we investigate in high-resolution N-body simulations of Milk-Way mass halos. We find that intergalactic black holes make up only ~2-3% of the total galactic black hole mass but, within a halo, wandering black holes can contribute up to about half of the total black hole mass orbiting the central galaxy. Intrahalo black holes offer a natural explanation for the compact X-ray sources often seen near the centres of galaxies and for the hyperluminous non-central X-ray source in M82.

The 11 known satellite galaxies within 250 kpc of the Milky Way lie close to a great circle on the sky. We use high resolution N-body simulations of galactic dark matter halos to test if this remarkable property can be understood within the context of the cold dark matter cosmology. We construct halo merger trees from the simulations and use a semianalytic model to follow the formation of satellite galaxies. We find that in all 6 of our simulations, the 11 brightest satellites are indeed distributed along thin, disk-like structures analogous to that traced by the Milky Way's satellites. This is in sharp contrast to the overall distributions of dark matter in the halo and of subhalos within it which, although triaxial, are not highly aspherical. We find that the spatial distribution of satellites is significantly different from that of the most massive subhalos but is similar to that of the subset of subhalos that had the most massive progenitors at earlier times. The elongated disk-like structure delineated by the satellites has its long axis aligned with the major axis of the dark matter halo. We interpret our results as reflecting the preferential infall of satellites along the spines of a few filaments of the cosmic web.