Search

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.

A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for models of galaxy formation. We investigate the impact of different halo merger trees on a semi-analytic model. We find that the z=0 galaxy properties in our model show differences between trees when using a common parameter set. The star formation history of the Universe and the properties of satellite galaxies can show marked differences between trees with different construction methods. Independently calibrating the semi-analytic model for each tree can reduce the discrepancies between the z=0 global galaxy properties, at the cost of increasing the differences in the evolutionary histories of galaxies. Furthermore, the underlying physics implied can vary, resulting in key quantities such as the supernova feedback efficiency differing by factors of 2. Such a change alters the regimes where star formation is primarily suppressed by supernovae. Therefore, halo merger trees extracted from a common halo catalogue using different, but reliable, algorithms can result in a difference in the semi-analytic model. Given the uncertainties in galaxy formation physics, however, these differences may not necessarily be viewed as significant., Comment: 15 pages, 18 figures, accepted for publication in MNRAS

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., Comment: 19 pages, 6 figures, accepted in ApJ

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, Comment: 8 pages, submitted to MN

Taking N-body simulations with volumes and particle densities tuned to match the SDSS DR7 spectroscopic main sample, we assess the ability of current void catalogs (e.g., Sutter et al. 2012b) to distinguish a model of coupled dark matter-dark energy from {\Lambda}CDM cosmology using properties of cosmic voids. Identifying voids with the VIDE toolkit, we find no statistically significant differences in the ellipticities, but find that coupling produces a population of significantly larger voids, possibly explaining the recent result of Tavasoli et al. (2013). In addition, we use the universal density profile of Hamaus et al. (2014) to quantify the relationship between coupling and density profile shape, finding that the coupling produces broader, shallower, undercompensated profiles for large voids by thinning the walls between adjacent medium-scale voids. We find that these differences are potentially measurable with existing void catalogs once effects from survey geometries and peculiar velocities are taken into account., Comment: 5 pages, 4 figures, moderate revision and extended discussion from referee comments, MNRAS letters accepted

We examine the integrated Sachs-Wolfe (ISW) imprint of voids and superclusters on the cosmic microwave background. We first study results from the Jubilee $N$-body simulation. From Jubilee, we obtain the full-sky ISW signal from structures out to redshift $z=1.4$ and a mock luminous red galaxy (LRG) catalogue. We confirm that the expected signal in the concordance \Lambda CDM model is very small and likely to always be much smaller than the anisotropies arising at the last scattering surface. Any current detections of such an imprint must, therefore, predominantly arise from something other than an ISW effect in a \Lambda CDM universe. Using the simulation as a guide, we then look for the signal using a catalogue of voids and superclusters from the Sloan Digital Sky Survey. We find a result that is consistent with the \Lambda CDM model, i.e. a signal consistent with zero., Comment: Now matches version published in MNRAS last October. Conclusions unchanged

We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range of disk scaling relationships. We find that the best fit parameters of a generic double power-law density profile vary in a systematic manner that depends on the stellar-to-halo mass ratio of each galaxy. Thus, the quantity Mstar/Mhalo constrains the inner ($\gamma$) and outer ($\beta$) slopes of dark matter density, and the sharpness of transition between the slopes($\alpha$), reducing the number of free parameters of the model to two. Due to the tight relation between stellar mass and halo mass, either of these quantities is sufficient to describe the dark matter halo profile including the effects of baryons. The concentration of the haloes in the hydrodynamical simulations is consistent with N-body expectations up to Milky Way mass galaxies, at which mass the haloes become twice as concentrated as compared with pure dark matter runs. This mass dependent density profile can be directly applied to rotation curve data of observed galaxies and to semi analytic galaxy formation models as a significant improvement over the commonly used NFW profile., Comment: 10 pages, 4 figures. Accepted for publication in MNRAS

We investigate the formation, growth, merger history, movement, and destruction of cosmic voids detected via the watershed transform code VIDE in a cosmological N-body dark matter {\Lambda}CDM simulation. By adapting a method used to construct halo merger trees, we are able to trace individual voids back to their initial appearance and record the merging and evolution of their progenitors at high redshift. For the scales of void sizes captured in our simulation, we find that the void formation rate peaks at scale factor 0.3, which coincides with a growth in the void hierarchy and the emergence of dark energy. Voids of all sizes appear at all scale factors, though the median initial void size decreases with time. When voids become detectable they have nearly their present-day volumes. Almost all voids have relatively stable growth rates and suffer only infrequent minor mergers. Dissolution of a void via merging is very rare. Instead, most voids maintain their distinct identity as annexed subvoids of a larger parent. The smallest voids are collapsing at the present epoch, but void destruction ceases after scale factor 0.3. In addition, voids centers tend to move very little, less than 0.01 of their effective radii per ln a, over their lifetimes. Overall, most voids exhibit little radical dynamical evolution; their quiet lives make them pristine probes of cosmological initial conditions and the imprint of dark energy., Comment: 11 pages, 12 figures, MNRAS accepted, minor revisions from referee comments

Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of \subhalos) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced., Comment: Accepted in MNRAS. 15 pages, 10 figures, 1 table

We study the $z=0$ properties of clusters (and large groups) of galaxies within the context of interacting and non-interacting quintessence cosmological models, using a series of adiabatic SPH simulations. Initially, we examine the average properties of groups and clusters, quantifying their differences in LCDM, uncoupled Dark Energy (\ude) and coupled Dark Energy (\cde) cosmologies. In particular, we focus upon radial profiles of the gas density, temperature and pressure, and we also investigate how the standard hydrodynamic equilibrium hypothesis holds in quintessence cosmologies. While we are able to confirm previous results about the distribution of baryons, we also find that the main discrepancy (with differences up to $20\%$) can be seen in cluster pressure profiles. We then switch attention to individual structures, mapping each halo in quintessence cosmology to its \LCDM\ counterpart. We are able to identify a series of small correlations between the coupling in the dark sector and halo spin, triaxiality and virialization ratio. When looking at spin and virialization of dark matter haloes, we find a weak ($5\%$) but systematic deviation in fifth force scenarios from \LCDM., Comment: 12 pages. Accepted for publication in MNRAS

We present the results of a series of adiabatic hydrodynamical simulations of several quintessence models (both with a free and an interacting scalar field) in comparison to a standard \LCDM\ cosmology. For each we use $2\times1024^3$ particles in a $250$\hMpc\ periodic box assuming WMAP7 cosmology. In this work we focus on the properties of haloes in the cosmic web at $z=0$. The web is classified into \emph{voids}, \emph{sheets}, \emph{filaments} and \emph{knots} depending on the eigenvalues of the velocity shear tensor, which are an excellent proxy for the underlying overdensity distribution. We find that the properties of objects classified according to their surrounding environment shows a substantial dependence on the underlying cosmology; for example, while $V_{\rm max}$ shows average deviations of $\approx5$ per cent across the different models when considering the full halo sample, comparing objects classified according to their environment, the size of the deviation can be as large as $20$ per cent. We also find that halo spin parameters are positively correlated to the coupling, whereas halo concentrations show the opposite behaviour. Furthermore, when studying the concentration-mass relation in different environments, we find that in all cosmologies underdense regions have a larger normalization and a shallower slope. While this behaviour is found to characterize all the models, differences in the best-fit relations are enhanced in (coupled) dark energy, thus providing a clearer prediction for this class of models., Comment: 16 pages, 10 figures. Accepted for publication in MNRAS

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness., Comment: 15 pages, 11 figures, 2 tables. 1 figure removed to reduce paper extension, shorter discussion, accepted for publication in MNRAS

The increasing precision in the determination of the Hubble parameter has reached a per cent level at which large-scale cosmic flows induced by inhomogeneities of the matter distribution become non-negligible. Here we use large-scale cosmological N-body simulations to study statistical properties of the local Hubble parameter as measured by local observers. We show that the distribution of the local Hubble parameter depends not only on the scale of inhomogeneities, but also on how one defines the positions of observers in the cosmic web and what reference frame is used. Observers located in random dark matter haloes measure on average lower expansion rates than those at random positions in space or in the centres of cosmic voids, and this effect is stronger from the halo rest frames compared to the CMB rest frame. We compare the predictions for the local Hubble parameter with observational constraints based on type Ia supernovae (SNIa) and CMB observations. Due to cosmic variance, for observers located in random haloes we show that the Hubble constant determined from nearby SNIa may differ from that measured from the CMB by 0.8 per cent at 1sigma statistical significance. This scatter is too small to significantly alleviate a recently claimed discrepancy between current measurements assuming a flat LCDM model. However, for observers located in the centres of the largest voids permitted by the standard LCDM model, we find that Hubble constant measurements from SNIa would be biased high by 5 per cent, rendering this tension inexistent in this extreme case., Comment: 9 pages, 6 figures; accepted for publication in MNRAS; updated references

We contend that a single power law halo mass distribution is appropriate for direct matching to the stellar masses of observed Local Group dwarf galaxies, allowing the determination of the slope of the stellar mass-halo mass relation for low mass galaxies. Errors in halo masses are well defined as the Poisson noise of simulated local group realisations, which we determine using constrained local universe simulations (CLUES). For the stellar mass range 10$^7$

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering., Comment: 19 pages, 15 figures

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., Comment: Accepted for publication in ApJ. 14 pages, 4 figures, 1 table

Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results., Comment: 14 pages,11 figures

We present initial results from the Jubilee ISW project, which models the expected \LambdaCDM Integrated Sachs-Wolfe (ISW) effect in the Jubilee simulation. The simulation volume is (6 Gpc/h)^3, allowing power on very large-scales to be incorporated into the calculation. Haloes are resolved down to a mass of 1.5x10^12 M_sun/h, which allows us to derive a catalogue of mock Luminous Red Galaxies (LRGs) for cross-correlation analysis with the ISW signal. We find the ISW effect observed on a projected sky to grow stronger at late times with the evolution of the ISW power spectrum matching expectations from linear theory. Maps of the gravitational lensing effect, including the convergence and deflection fields, are calculated using the same potential as for the ISW. We calculate the redshift dependence of the ISW-LRG cross-correlation signal for a full sky survey with no noise considerations. For l < 30, the signal is strongest for lower redshift bins (z ~ 0.2 to 0.5), whereas for l > 30 the signal is best observed with surveys covering z ~ 0.6-1.0., Comment: 14 pages. 13 Figures. For submission to MNRAS

We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9

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., Comment: 13 Pages, 5 figures, Accepted for publication in Publications of the Astronomical Society of Australia special issue on Warm Dark Matter

While various codes exist to systematically and robustly find haloes and subhaloes in cosmological simulations (Knebe et al., 2011, Onions et al., 2012), this is the first work to introduce and rigorously test codes that find tidal debris (streams and other unbound substructure) in fully cosmological simulations of structure formation. We use one tracking and three non-tracking codes to identify substructure (bound and unbound) in a Milky Way type simulation from the Aquarius suite (Springel et al., 2008) and post-process their output with a common pipeline to determine the properties of these substructures in a uniform way. By using output from a fully cosmological simulation, we also take a step beyond previous studies of tidal debris that have used simple toy models. We find that both tracking and non-tracking codes agree well on the identification of subhaloes and more importantly, the {\em unbound tidal features} associated with them. The distributions of basic properties of the total substructure distribution (mass, velocity dispersion, position) are recovered with a scatter of $\sim20%$. Using the tracking code as our reference, we show that the non-tracking codes identify complex tidal debris with purities of $\sim40%$. Analysing the results of the substructure finders, we find that the general distribution of {\em substructures} differ significantly from the distribution of bound {\em subhaloes}. Most importantly, both bound and unbound {\em substructures} together constitute $\sim18%$ of the host halo mass, which is a factor of $\sim2$ higher than the fraction in self-bound {\em subhaloes}. However, this result is restricted by the remaining challenge to cleanly define when an unbound structure has become part of the host halo. Nevertheless, the more general substructure distribution provides a more complete picture of a halo's accretion history., Comment: 19 pages, 12 figures, accepted for publication in MNRAS

We present the first results from the JUropa huBbLE volumE (Jubilee) project, based a large N-body, dark matter-only cosmological simulation with a volume of $V=(6 h^{-1}\mathrm{Gpc})^3$, containing 6000$^3$ particles, performed within the concordance $\Lambda$CDM cosmological model. The simulation volume is sufficient to probe extremely large length scales in the universe, whilst at the same time the particle count is high enough so that dark matter haloes down to $1.5\times10^{12} h^{-1}\mathrm{M}_\odot$ can be resolved. At $z = 0$ we identify over 400 million haloes. The cluster mass function is derived using three different halofinders and compared to fitting functions in the literature. The distribution of clusters of maximal mass across redshifts agrees well with predicted masses of extreme objects, and we explicitly confirm that the Poisson distribution is very good at describing the distribution of rare clusters. The Poisson distribution also matches well the level to which cosmic variance can be expected to affect number counts of high mass clusters. We find that objects like the Bullet cluster exist in the far-tail of the distribution of mergers in terms of relative collisional speed. We also derive the number counts of voids in the simulation box for $z = 0$, $0.5$ and $1$., Comment: Version 2. 12 pages, 9 figures. Accepted by MNRAS

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general., Comment: 28 pages containing 13 figures & 4 tables + 9 pages appendix containing another 4 tables + 4 pages of references, accepted for publication in MNRAS

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., Comment: 8 pages, 3 figures, Submitted to ApJ Letters

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter., Comment: 10 pages, 14 figures and 3 tables. Accepted for publication in MNRAS

In this paper we investigate how the halo mass function evolves with redshift, based on a suite of very large (with N_p = 3072^3 - 6000^3 particles) cosmological N-body simulations. Our halo catalogue data spans a redshift range of z = 0-30, allowing us to probe the mass function from the dark ages to the present. We utilise both the Friends-of-Friends (FOF) and Spherical Overdensity (SO) halofinding methods to directly compare the mass function derived using these commonly used halo definitions. The mass function from SO haloes exhibits a clear evolution with redshift, especially during the recent era of dark energy dominance (z < 1). We provide a redshift-parameterised fit for the SO mass function valid for the entire redshift range to within ~20% as well as a scheme to calculate the mass function for haloes with arbitrary overdensities. The FOF mass function displays a weaker evolution with redshift. We provide a `universal' fit for the FOF mass function, fitted to data across the entire redshift range simultaneously, and observe redshift evolution in our data versus this fit. The relative evolution of the mass functions derived via the two methods is compared and we find that the mass functions most closely match at z=0. The disparity at z=0 between the FOF and SO mass functions resides in their high mass tails where the collapsed fraction of mass in SO haloes is ~80% of that in FOF haloes. This difference grows with redshift so that, by z>20, the SO algorithm finds a ~50-80% lower collapsed fraction in high mass haloes than does the FOF algorithm, due in part to the significant over-linking effects known to affect the FOF method., Comment: v4, 16 pages, 16 colour figures. Changed to match MNRAS print version. NOTE: v1 of this paper has a typo in the fitting function. Please ensure you use the latest version

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., Comment: 14 Pages, 10 Figures; accepted MNRAS 2012

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., Comment: 14 interesting pages, 17 beautiful figures, and 2 informative tables accepted for publication in MNRAS (matches published version)

We examine the properties and evolution of a simulated polar disc galaxy. This galaxy is comprised of two orthogonal discs, one of which contains old stars (old stellar disc), and the other, containing both younger stars and the cold gas (polar disc) of the galaxy. By exploring the shape of the inner region of the dark matter halo, we are able to confirm that the halo shape is a oblate ellipsoid flattened in the direction of the polar disc. We also note that there is a twist in the shape profile, where the innermost 3 kpc of the halo flattens in the direction perpendicular to the old disc, and then aligns with the polar disc out until the virial radius. This result is then compared to the halo shape inferred from the circular velocities of the two discs. We also use the temporal information of the simulation to track the system's evolution, and identify the processes which give rise to this unusual galaxy type. We confirm the proposal that the polar disc galaxy is the result of the last major merger, where the angular moment of the interaction is orthogonal to the angle of the infalling gas. This merger is followed by the resumption of coherent gas infall. We emphasise that the disc is rapidly restored after the major merger and that after this event the galaxy begins to tilt. A significant proportion of the infalling gas comes from filaments. This infalling gas from the filament gives the gas its angular momentum, and, in the case of the polar disc galaxy, the direction of the gas filament does not change before or after the last major merger., Comment: Accepted for publication in MNRAS; 14 pages; 14 figures

We present the results of a series of cosmological $N$-body simulations of a Vector Dark Energy (VDE) model, performed using a suitably modified version of the publicly available \texttt{GADGET}-2 code. The setups of our simulations were calibrated pursuing a twofold aim: 1) to analyze the large scale distribution of massive objects and 2) to determine the properties of halo structure in this different ramework.We observe that structure formation is enhanced in VDE, since the mass function at high redshift is boosted up to a factor of ten with respect to \LCDM, possibly alleviating tensions with the observations of massive clusters at high redshifts and early reionization epoch. Significant differences can also be found for the value of the growth factor, that in VDE shows a completely different behaviour, and in the distribution of voids, which in this cosmology are on average smaller and less abundant. We further studied the structure of dark matter haloes more massive than $5\times10^{13}$\hMsun, finding that no substantial difference emerges when comparing spin parameter, shape, triaxiality and profiles of structures evolved under different cosmological pictures. Nevertheless, minor differences can be found in the concentration-mass relation and the two point correlation function; both showing different amplitudes and steeper slopes.Using an additional series of simulations of a \LCDM\ scenario with the same $\Omega_M$ and $\sigma_8$ used in the VDE cosmology, we have been able to establish whether the modifications induced in the new cosmological picture were due to the particular nature of the dynamical dark energy or a straightforward consequence of the cosmological parameters., Comment: Accepted for publication in MNRAS. 18 pages, 19 figures

We study the non-linear gravitational collapse of dark matter into halos through numerical N-body simulations of Lemaitre-Tolman-Bondi void models. We extend the halo mass function formalism to these models in a consistent way. This extension not only compares well with the simulated data at all times and radii, but it also gives interesting clues about the impact of the background shear on the growth of perturbations. Our results give hints about the possibility of constraining the background shear via cluster number counts, which could then give rise to strong constraints on general inhomogeneous models, of any scale., Comment: 5 pages, 3 figures, Accepted in Physics of the Dark Universe, preprint IFT-UAM/CSIC-12-31

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., Comment: 11 pages, 5 figs, MNRAS accepted version, figures added, section 4 expanded

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(

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important., Comment: 16 pages, 7 figures, 2 tables, Accepted for MNRAS

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., Comment: 5 pages, 2 figures, MNRAS Letters in Accepted

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., Comment: 8 pages, 4 Figures, MNRAS Accepted 2012 June 19. Received 2012 May 10; in original form 2011 August 25

We investigate how long-range scalar interactions affect the properties of dark matter haloes. For doing so we employ the ReBEL model which implements an additional interaction between dark matter particles. On the phenomenological level this is equivalent to a modification of gravity. We analyse the differences between five ReBEL models and $\Lambda$CDM using a series of high resolution cosmological simulations. Emphasis is placed on investigating how halo properties change in the presence of a fifth force. We report that the density profile of ReBEL haloes is well described by the NFW profile but with mean concentrations from $5\%$ to a few times higher than the standard $\Lambda$CDM value. We also find a slight increase of the halo spin for haloes more massive than $5\times10^{11}\M_{\odot}$, reflecting a higher rotational support of those haloes due to scalar forces. In addition, the dark matter haloes in our models are more spherical than their counterparts in $\Lambda$CDM. The ReBEL haloes are also more virialised, with a large difference from $\Lambda$CDM for strong fifth forces and a much smaller change for weak scalar interactions., Comment: Published in JCAP (Issue 10, article id. 012, pp. 2013),32 pages (JCAP format), 13 figures, major revision version, some conclusions changed

Using a statistical sample of dark matter haloes drawn from a suite of cosmological N-body simulations of the Cold Dark Matter (CDM) model, we quantify the impact of a simulated halo's mass accretion and merging history on two commonly used measures of its dynamical state, the virial ratio eta and the centre of mass offset Delta r. Quantifying this relationship is important because the degree to which a halo is dynamically equilibrated will influence the reliability with which we can measure characteristic equilibrium properties of the structure and kinematics of a population of haloes. We begin by verifying that a halo's formation redshift zform correlates with its virial mass Mvir and we show that the fraction of its recently accreted mass and the likelihood of it having experienced a recent major merger increases with increasing Mvir and decreasing zform. We then show that both eta and Delta r increase with increasing Mvir and decreasing zform, which implies that massive recently formed haloes are more likely to be dynamically unrelaxed than their less massive and older counterparts. Our analysis shows that both eta and Delta r are good indicators of a halo's dynamical state, showing strong positive correlations with recent mass accretion and merging activity, but we argue that Delta r provides a more robust and better defined measure of dynamical state for use in cosmological N-body simulations at z~0. We find that Delta r < 0.04 is sufficient to pick out dynamically relaxed haloes at z=0. Finally, we assess our results in the context of previous studies, and consider their observational implications., Comment: Accepted for publication in MNRAS. 13 pages, 15 figures

The detection of extremely massive clusters at $z>1$ such as SPT-CL J0546-5345, SPT-CL J2106-5844, and XMMU J2235.3-2557 has been considered by some authors as a challenge to the standard \LCDM$\;$cosmology. In fact,assuming Gaussian initial conditions, the theoretical expectation of detecting such objects is as low as $\leq 1%$. In this \textit{Letter} we discuss the probability of the existence of such objects in the light of the Vector Dark Energy (VDE) paradigm, showing by means of a series of $N$-body simulations that chances of detection are substantially enhanced in this non-standard framework., Comment: Accepted for publications in MNRAS, 6 pages, 3 figures. Replaced 1st version due to problems with the first submitted .tex file

We perform collisionless N-body simulations to investigate whether binary mergers between rotationally-supported dwarfs can lead to the formation of dwarf spheroidal galaxies (dSphs). Our simulation campaign is based on a hybrid approach combining cosmological simulations and controlled numerical experiments. We select merger events from a Constrained Local UniversE (CLUES) simulation of the Local Group (LG) and record the properties of the interacting dwarf-sized halos. This information is subsequently used to seed controlled experiments of binary encounters between dwarf galaxies consisting of exponential stellar disks embedded in cosmologically-motivated dark matter halos. These simulations are designed to reproduce eight cosmological merger events, with initial masses of the interacting systems in the range ~ (5-60) x 10^7 Mo, occurring quite early in the history of the LG, more than 10 Gyr ago. We compute the properties of the merger remnants as a distant observer would and demonstrate that at least three of the simulated encounters produce systems with kinematic and structural properties akin to those of the classic dSphs in the LG. Tracing the history of the remnants in the cosmological simulation to z=0, we find that two dSph-like objects remain isolated at distances larger than 800 kpc from either the Milky Way or M31. These systems constitute plausible counterparts of the remote dSphs Cetus and Tucana which reside in the LG outskirts, far from the tidal influence of the primary galaxies. We conclude that merging of rotationally-supported dwarfs represents a viable mechanism for the formation of dSphs in the LG and similar environments., Comment: 6 pages, 2 figures, LaTeX (uses emulateapj.cls). Replaced to match the version accepted for publication in ApJL. Conclusions remain unchanged

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., Comment: 6 pages, 3 figures. Accepted as MNRAS Letter 2011 July 21. Received 2011 July 20; in original form 2011 June 6

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., Comment: 13 pages, 6 figures, 1 table. Accepted by MNRAS 2011 July 20. Received 2011 June 24; in original form 2011 April 21

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., Comment: 5 pages, 5 figures, 1 table, accepted for publication as MNRAS Letter

Semi-analytic models (SAMs) are currently one of the primary tools with which to model statistically significant ensembles of galaxies. The underlying physical prescriptions inherent to each SAM are, in many cases, different from one another. Several SAMs have been applied to the dark matter merger trees extracted from the Millennium Run, including those associated with the well-known Munich and Durham lineages. We compare the predicted luminosity distributions of galaxy groups using four publicly available SAMs (De Lucia et al. 2006; Bower et al. 2006; Bertone et al. 2007; Font et al. 2008), in order to explore a galactic environment in which the models have not been explored to the same degree as they have in the field or in rich clusters. We identify a characteristic "wiggle" in the group galaxy luminosity function generated using the De Lucia et al. (2006) SAM, that is not present in the Durham-based models, consistent to some degree with observations. However, a comparison between conditional luminosity functions of groups between the models and observations of Yang et al. (2007) suggest that neither model is a particularly good match. The luminosity function wiggle is interpreted as the result of the two-mode AGN feedback implementation used in the Munich models, which itself results in flattened magnitude gap distribution. An associated analysis of the magnitude gap distribution between first- and second-ranked group galaxies shows that while the Durham models yield distributions with approximately equal luminosity first- and second-ranked galaxies, in agreement with observations, the De Lucia et al. models favours the scenario in which the second-ranked galaxy is approximately one magnitude fainter than the primary,especially when the dynamic range of the mock data is limited to 3 magnitudes., Comment: 16 pages, 12 figures, 4 tables, accepted MNRAS; Revision:removed bold paragraphs in Sections 3.1 and 6

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve)., Comment: 27 interesting pages, 20 beautiful figures, and 4 informative tables accepted for publication in MNRAS. The high-resolution version of the paper as well as all the test cases and analysis can be found at the web site http://popia.ft.uam.es/HaloesGoingMAD

Relaxed, massive galactic objects have been identified at redshifts z = 4;5; and 6 in hydrodynamical simulations run in a large cosmological volume. This allowed us to analyze the assembly patterns of the high mass end of the galaxy distribution at these high zs, by focusing on their structural and dynamical properties. Our simulations indicate that massive objects at high redshift already follow certain scaling relations. These relations define virial planes at the halo scale, whereas at the galactic scale they define intrinsic dynamical planes that are, however, tilted relative to the virial plane. Therefore, we predict that massive galaxies must lie on fundamental planes from their formation. We briefly discuss the physical origin of the tilt in terms the physical processes underlying massive galaxy formation at high z, in the context of a two-phase galaxy formation scenario. Specifically, we have found that it lies on the different behavior of the gravitationally heated gas as compared with cold gas previously involved in caustic formation, and the mass dependence of the energy available to heat the gas., Comment: 12 pages, 1 figure. Accepted to Astrophysical Journal Letters

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., Comment: 9 pages, 10 figures, 1 table, accepted for publication in MNRAS

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., Comment: 14 Pages, 8 Figures, accpeted by MNRAS

We investigate the effect of long-range scalar interactions in dark matter (DM) models of cosmic structure formation with a particular focus on the formation times of haloes. Utilising $N$-body simulations with $512^3$ DM particles we show that in our models dark matter haloes form substantially earlier: tracing objects up to redshift $z\sim6$ we find that the formation time, as characterised by the redshift $z_{1/2}$ at which the halo has assembled half of its final mass, is gradually shifted from $z_{1/2}\approx 1.83$ in the fiducial LCDM model to $z_{1/2}\approx 2.54$ in the most extreme self-interaction model. This is accompanied by a shift of the redshift that marks the transition between merger and steady accretion epochs from $z_{*}\approx 4.32$ in the \lcdm\ halos to $z_{*}\approx 6.39$ in our strongest interaction model. In other words, the scalar-interacting model employed in this work produces more structures at high redshifts, prolonging at the same time the steady accretion phases. These effects taken together can help the LCDM model to account for a high redshift reionisation as indicated by the WMAP data and can alleviate issues related to the survival of the thin-disk dominated galaxies at low redshifts., Comment: 5 pages, 3 figures, 1 table. Final version to be printed in MNRAS Letters. 4 new references added