Source author record

Noam I. Libeskind

Noam I. Libeskind appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2022arXiv

Anisotropic satellite accretion onto the Local Group with HESTIA

How the cosmic web feeds halos, and fuels galaxy formation is an open question with wide implications. This study explores the mass assembly in the Local Group within the context of the local cosmography by employing simulations whose initial conditions have been constrained to reproduce the local environment. The goal of this study is to inspect whether the direction of accretion of satellites on to the Milky Way and Andromeda galaxies, is related to the cosmic web. The analysis considers the three high-resolution simulations available in the HESTIA simulation suite, as well as the derived velocity shear and tidal tensors. We notice two eras in the Local Group accretion history, delimited by an epoch around $z \approx 0.7$. We also find that satellites can travel up to $\sim 4$ Mpc, relative to their parent halo before crossing its viral radius $R_{200}$. Finally, we observe a strong alignment of the infall direction with the axis of slowest collapse $\vec{e_3}$ of both tidal and shear tensors, implying satellites of the Local Group originated from one particular region of the cosmic web and were channeled towards us via the process of accretion.This alignment is dominated by the satellites that enter during the early infall era, i.e $z>0.7$.

preprint2022arXiv

COWS: A filament finder for Hessian cosmic web identifiers

The large scale galaxy and matter distribution is often described by means of the cosmic web made up of voids, sheets, filaments and knots. Many different recipes exist for identifying this cosmic web. Here we focus on a sub-class of cosmic web identifiers, based on the analysis of the Hessian matrix, and proposed a method, called COsmic Web Skeleton (COWS), of separating a set of filaments cells into an ensemble of individual discreet filaments. Specifically, a thinning algorithm is applied to velocity shear tensor based cosmic web (V-web) to identify the spine of the filaments. This results in a set of filaments with well defined end-point and length. It is confirmed that these sit at local density ridges and align with the appropriate direction defined by the underlying velocity field. The radial density profile of these curved cylindrical filaments, as well as the distribution of their lengths is also examined. The robustness of all results is checked against resolution and the V-web threshold. The code for the COWS method has been made publicly available.

preprint2022arXiv

Hermeian haloes: Field haloes that interacted with both the Milky Way and M31

The Local Group is a unique environment in which to study the astrophysics of galaxy formation. The proximity of the Milky Way and M31 enhances the frequency of interactions of the low-mass halo population with more massive dark matter haloes, which increases their concentrations and strips them of gas and other material. Some low-mass haloes pass through the haloes of the Milky Way or M31 and are either ejected into the field or exchanged between the two primary hosts. We use high resolution gas-dynamical simulations to describe a new class of field haloes that passed through the haloes of both the Milky Way and M31 at early times and are almost twice as concentrated as field haloes that do not interact with the primary pair. These 'Hermeian' haloes are distributed anisotropically at larger distances from the Local Group barycentre than the primary haloes and appear to cluster along the line connecting the Milky Way and M31. Hermeian haloes facilitate the exchange of dark matter, gas, and stars between the Milky Way and M31 and can enhance the star formation rate of the gas in the primary haloes during their interactions with them. We also show that some Hermeian haloes can host galaxies that, because they are embedded in haloes that are more concentrated than regular field haloes, are promising targets for indirect dark matter searches beyond the Milky Way virial radius and can produce signals that are competitive with those of some dwarf galaxies. Hermeian galaxies in the Local Group should be detectable by forthcoming wide-field imaging surveys.

preprint2022arXiv

NIHAO-LG: The uniqueness of Local Group dwarf galaxies

Recent observational and theoretical studies of the Local Group (LG) dwarf galaxies have highlighted their unique star-formation history, stellar metallicity, gas content, and kinematics. We investigate the commonality of these features by comparing constrained LG and field central dwarf halo simulations in the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project. Our simulations, performed with NIHAO-like hydrodynamics which track the evolution of the Milky Way (MW) and M31 along with ~100 dwarfs in the LG, reveal the total gas mass and stellar properties (velocity dispersion, evolution history, etc.) of present-day LG dwarfs to be similar to field systems. However, relative to field galaxies, LG dwarfs have more cold gas in their central parts and more metal-rich gas in the halo stemming from interactions with other dwarfs living in a high-density environment like the LG. Interestingly, the direct impact of massive MW/M31 analogues on the metallicity evolution of LG dwarfs is minimal; LG dwarfs accrete high-metallicity gas mostly from other dwarfs at late times. We have also tested for the impact of metal diffusion on the chemical evolution of LG dwarfs, and found that it does not affect the stellar or gaseous content of LG dwarfs. Our simulations suggest that the stellar components of LG dwarfs offer a unique and unbiased local laboratory for galaxy-formation tests and comparisons, especially against the overall dwarf population in the Universe.

preprint2020arXiv

A robust determination of halo environment in the cosmic field

A number of methods for studying the large-scale cosmic matter distribution exist in the literature. One particularly common method employed to define the cosmic web is to examine the density, velocity or potential field. Such methods are advantageous since a Hessian matrix can be constructed whose eigenvectors (and eigenvalues) indicate the principal directions (and strength) of local collapse or expansion. Technically this is achieved by diagonalizing the Hessian matrix using a fixed finite grid. The resultant large-scale structure quantification is thus inherently limited by the grid's finite resolution. Here, we overcome the obstacle of finite grid resolution by introducing a new method to determine halo environment using an adaptive interpolation which is more robust to resolution than the typical "Nearest Grid Point" (NGP) method. Essentially instead of computing and diagonalizing the Hessian matrix once for the entire grid, we suggest doing so once for each halo or galaxy in question. We examine how the eigenvalues and eigenvector direction's computed using our algorithm and the NGP method converge for different grid resolutions, finding that our new method is convergent faster. Namely changes of resolution have a much smaller effect than in the NGP method. We therefore suggest this method for future use by the community.

preprint2020arXiv

The alignment of satellite systems with cosmic filaments in the SDSS DR12

Galaxies, as well as their satellites, are known to form within the cosmic web: the large, multi-scale distribution of matter in the universe. It is known that the surrounding large scale structure (LSS) can impact and influence the formation of galaxies, e.g. the spin and shape of haloes or galaxies are correlated with the LSS and the correlation depends on halo mass or galaxy morphology. In this work, we use group and filament catalogues constructed from the SDSS DR12 to investigate the correlation between satellite systems and the large scale filaments they are located in. We find that the distribution of satellites is significantly correlated with filaments, namely the major axis of the satellite systems are preferentially aligned with the spine of the closest filament. Stronger alignment signals are found for the cases where the system away from the filament spine, while systems close to the filament spine show significantly weaker alignment. Our results suggest that satellites are accreted along filaments, which agrees with previous works. The case of which away from the filament spine may help us to understand how the filament forms as well as the peculiar satellite distribution in the Local Universe.

preprint2020arXiv

The Hestia project: simulations of the Local Group

We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter only simulations, and selected for higher resolution magneto hydrodynamic simulation using the Arepo code. Baryonic physics follows the Auriga model of galaxy formation. The simulations contain a high resolution region of 3-5 Mpc in radius from the Local Group midpoint embedded in the correct cosmographic landscape. Within this region a simulated Local Group consisting of a Milky Way and Andromeda like galaxy forms, whose description is in excellent agreement with observations. The simulated Local Group galaxies resemble the Milky Way and Andromeda in terms of their halo mass, mass ratio, stellar disc mass, morphology separation, relative velocity, rotation curves, bulge-disc morphology, satellite galaxy stellar mass function, satellite radial distribution and in some cases, the presence of a Magellanic cloud like object. Because these simulations properly model the Local Group in their cosmographic context, they provide a testing ground for questions where environment is thought to play an important role.

preprint2019arXiv

On the Mass Assembly History of the Local Group

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 $Λ$ Cold Dark Matter ($Λ$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-$σ$ 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.

preprint2019arXiv

Simulating the galactic system in interaction AM 2229-735 and the formation of its polar structure

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.

preprint2016arXiv

Constrained Local UniversE Simulations: A Local Group Factory

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

preprint2015arXiv

Filaments from the galaxy distribution and from the velocity field in the local universe

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.

preprint2015arXiv

Galaxies in filaments have more satellites: the influence of the cosmic web on the satellite luminosity function in the SDSS

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.

preprint2015arXiv

Planes of satellite galaxies and the cosmic web

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 (ΛCDM) model of cosmology if not even in contradiction to its generic predictions. Guided by Λ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 ΛCDM model.

preprint2015arXiv

Properties of Galaxies in the Disc Central Surface Brightness Gap

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.

preprint2015arXiv

The alignment of galaxy spin with the shear field in observations

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.

preprint2015arXiv

Two Planes of Satellites in the Centaurus A Group

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.

preprint2014arXiv

Galaxy properties and the cosmic web in simulations

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.

preprint2013arXiv

CLUES to the past: Local Group progenitors amongst high-redshift Lyman Break Galaxies

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.

preprint2013arXiv

Cold versus Warm Dark Matter simulations of a galaxy group

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.

preprint2013arXiv

Galaxies going MAD: The Galaxy-Finder Comparison Project

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.

preprint2013arXiv

Galaxy spin alignment in filaments and sheets: observational evidence

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.

preprint2013arXiv

Size matters: the non-universal density profile of subhaloes in SPH simulations and implications for the Milky Way's dSphs

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.

preprint2012arXiv

A kinematic classification of the cosmic web

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.

preprint2012arXiv

Applying scale-free mass estimators to the Local Group in Constrained Local Universe Simulations

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-γ}$ is directly related to the host's mass profile $M(<r)\propto r^{1-α}$ as $α=γ-2$. The use of this relation for any given $γ$ 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.

preprint2012arXiv

The cosmic web and the orientation of angular momenta

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.

preprint2012arXiv

The velocity shear tensor: tracer of halo alignment

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.

preprint2011arXiv

Too small to succeed? Lighting up massive dark matter subhaloes of the Milky Way

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.