BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stefan Gottloeber

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-phApplicationsastro-ph.IMDatabasesgr-qc

Catalog footprint

What is connected

49works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

49 published item(s)

preprint2022arXiv

Reionization time of the Local Group and Local-Group-like halo pairs

Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired halos within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully-coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation and isolation criteria but in other environments, alongside isolated halos within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out given our reionization model. Within this model, halos in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ~625 Myr against 660+/-4 Myr (z~8.25 against 7.87+/-0.02) on average.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxy Three-Point Correlation Functions and Halo/Subhalo Models

We present the measurements of the luminosity-dependent redshift-space three-point correlation functions (3PCFs) for the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample. We compare the 3PCF measurements to the predictions from three different halo and subhalo models. One is the halo occupation distribution (HOD) model and the other two are extensions of the subhalo abundance matching (SHAM) model by allowing the central and satellite galaxies to have different occupation distributions in the host halos and subhalos. Parameters in all the models are chosen to best describe the projected and redshift-space two-point correlation functions (2PCFs) of the same set of galaxies. All three model predictions agree well with the 3PCF measurements for the most luminous galaxy sample, while the HOD model better performs in matching the 3PCFs of fainter samples (with luminosity threshold below $L^*$), which is similar in trend to the case of fitting the 2PCFs. The decomposition of the model 3PCFs into contributions from different types of galaxy triplets shows that on small scales the dependence of the 3PCFs on triangle shape is driven by nonlinear redshift-space distortion (and not by the intrinsic halo shape) while on large scales it reflects the filamentary structure. The decomposition also reveals more detailed differences in the three models, which are related to the radial distribution, the mean occupation function, and the velocity distribution of satellite galaxies inside halos. The results suggest that galaxy 3PCFs can further help constrain the above galaxy-halo relation and test theoretical models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How did the Virgo cluster form?

While the Virgo cluster is the nearest galaxy cluster and therefore the best observed one, little is known about its formation history. In this paper, a set of cosmological simulations that resemble the Local Universe is used to shed the first light on this mystery. The initial conditions for these simulations are constrained with galaxy peculiar velocities of the second catalog of the Cosmicflows project using algorithms developed within the Constrained Local UniversE Simulation project. Boxes of 500 Mpc/h on a side are set to run a series of dark matter only constrained simulations. In each simulation, a unique dark matter halo can be reliably identified as Virgo's counterpart. The properties of these Virgo halos are in agreement at a 10-20% level with the global properties of the observed Virgo cluster. Their zero-velocity masses agree at one-sigma with the observational mass estimate. In all the simulations, the matter falls onto the Virgo objects along a preferential direction that corresponds to the observational filament and the slowest direction of collapse. A study of the mass accretion history of the Virgo candidates reveals the most likely formation history of the Virgo cluster, namely a quiet accretion over the last 7 Gigayears.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modelling galaxy clustering: halo occupation distribution versus subhalo matching

We model the luminosity-dependent projected and redshift-space two-point correlation functions (2PCFs) of the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, using the halo occupation distribution (HOD) model and the subhalo abundance matching (SHAM) model and its extension. All the models are built on the same high-resolution $N$-body simulations. We find that the HOD model generally provides the best performance in reproducing the clustering measurements in both projected and redshift spaces. The SHAM model with the same halo-galaxy relation for central and satellite galaxies (or distinct haloes and subhaloes), when including scatters, has a best-fitting $χ^2/\rm{dof}$ around $2$--$3$. We therefore extend the SHAM model to the subhalo clustering and abundance matching (SCAM) by allowing the central and satellite galaxies to have different galaxy--halo relations. We infer the corresponding halo/subhalo parameters by jointly fitting the galaxy 2PCFs and abundances and consider subhaloes selected based on three properties, the mass $M_{\rm acc}$ at the time of accretion, the maximum circular velocity $V_{\rm acc}$ at the time of accretion, and the peak maximum circular velocity $V_{\rm peak}$ over the history of the subhaloes. The three subhalo models work well for luminous galaxy samples (with luminosity above $L_*$). For low-luminosity samples, the $V_{\rm acc}$ model stands out in reproducing the data, with the $V_{\rm peak}$ model slightly worse, while the $M_{\rm acc}$ model fails to fit the data. We discuss the implications of the modeling results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmicflows Constrained Local UniversE Simulations

This paper combines observational datasets and cosmological simulations to generate realistic numerical replicas of the nearby Universe. These latter are excellent laboratories for studies of the non-linear process of structure formation in our neighborhood. With measurements of radial peculiar velocities in the Local Universe (cosmicflows-2) and a newly developed technique, we produce Constrained Local UniversE Simulations (CLUES). To assess the quality of these constrained simulations, we compare them with random simulations as well as with local observations. The cosmic variance, defined as the mean one-sigma scatter of cell-to-cell comparison between two fields, is significantly smaller for the constrained simulations than for the random simulations. Within the inner part of the box where most of the constraints are, the scatter is smaller by a factor 2 to 3 on a 5 Mpc/h scale with respect to that found for random simulations. This one-sigma scatter obtained when comparing the simulated and the observation-reconstructed velocity fields is only 104 +/- 4 km/s i.e. the linear theory threshold. These two results demonstrate that these simulations are in agreement with each other and with the observations of our neighborhood. For the first time, simulations constrained with observational radial peculiar velocities resemble the Local Universe up to a distance of 150 Mpc/h on a scale of a few tens of megaparsecs. When focusing on the inner part of the box, the resemblance with our cosmic neighborhood extends to a few megaparsecs (< 5 Mpc/h). The simulations provide a proper Large Scale environment for studies of the formation of nearby objects.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mergers and the outside-in formation of dwarf spheroidals

We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few examples where stars form in two populations distinct in age in order to elucidate the origin of these gradients. The spatial distributions of the two components reflect their diverse origin; the old stellar component is assembled through mergers, but the young population forms largely in situ. The older component results from a first episode of star formation that begins early but is quickly shut off by the combined effects of stellar feedback and reionization. The younger component forms when a late accretion event adds gas and reignites star formation. The effect of mergers is to disperse the old stellar population, increasing their radius and decreasing their central density relative to the young population. We argue that dwarf-dwarf mergers offer a plausible scenario for the formation of systems with multiple distinct populations and, more generally, for the origin of age and metallicity gradients in dwarf spheroidals.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Effect of Environment on Milky Way-mass galaxies in a Constrained Simulation of the Local Group

In this letter we present, for the first time, a study of star formation rate, gas fraction and galaxy morphology of a constrained simulation of the Milky Way (MW) and Andromeda (M31) galaxies, compared to other MW-mass galaxies. By combining with unconstrained simulations we cover a sufficient volume to compare these galaxies environmental densities ranging from the field to that of the Local Group (LG). This is particularly relevant as it has been shown that, quite generally, galaxy properties depend intimately upon their environment, most prominently when galaxies in clusters are compared to those in the field. For galaxies in loose groups such as the LG, however, environmental effects have been less clear. We consider the galaxy's environmental density in spheres of 1200 kpc (comoving) and find that whilst environment does not appear to directly affect morphology, there is a positive trend with star formation rates. This enhancement in star formation occurs systematically for galaxies in higher density environments, regardless whether they are part of the LG or in filaments. Our simulations suggest that the richer environment at Mpc-scales may help replenish the star-forming gas, allowing higher specific star formation rates in galaxies such as the MW.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Milky Way and Andromeda galaxies in a constrained hydrodynamical simulation: morphological evolution

We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry out our study. We run two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the impact of this process on the morphologies and star formation rates of the simulated galaxies. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z=0, the same environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final morphologies. The M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms a disc that grows steadily; at z=0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In summary, by comparing the formation histories of the two galaxies, we conclude that the particular merger/accretion history of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of galaxy discs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The distribution of gas in the Local Group from constrained cosmological simulations: the case for Andromeda and the Milky Way galaxies

We study the gas distribution in the Milky Way and Andromeda using a constrained cosmological simulation of the Local Group (LG) within the context of the CLUES (Constrained Local UniversE Simulations) project. We analyse the properties of gas in the simulated galaxies at $z=0$ for three different phases: `cold', `hot' and HI, and compare our results with observations. The amount of material in the hot halo ($M_{hot}\approx 4-5\times10^{10}\,$M$_{\odot}$), and the cold ($M_{cold}(r\lesssim10\,$kpc$)\approx10^{8}\,$M$_{\odot}$) and HI ($M_{HI}(r\lesssim50\,$kpc$)\approx 3-4\times10^8\,$M$_{\odot}$) components display a reasonable agreement with observations. We also compute the accretion/ejection rates together with the HI (radial and all-sky) covering fractions. The integrated HI accretion rate within $r=50\,$kpc gives $\sim$$0.2-0.3\,$M$_{\odot}\,$yr$^{-1}$, i.e. close to that obtained from high-velocity clouds in the Milky Way. We find that the global accretion rate is dominated by hot material, although ionized gas with $T\lesssim10^5\,$K can contribute significantly too. The $net$ accretion rates of $all$ material at the virial radii are $6-8\,$M$_{\odot}\,$yr$^{-1}$. At $z=0$, we find a significant gas excess between the two galaxies, as compared to any other direction, resulting from the overlap of their gaseous haloes. In our simulation, the gas excess first occurs at $z\sim1$, as a consequence of the kinematical evolution of the LG.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Imprint of Reionization on the Star Formation Histories of Dwarf Galaxies

We explore the impact of cosmic reionization on nearby isolated dwarf galaxies using a compilation of SFHs estimated from deep HST data and a cosmological hydrodynamical simulation of the Local Group. The nearby dwarfs show a wide diversity of star formation histories; from ancient systems that have largely completed their star formation $\sim 10$ Gyr ago to young dwarfs that have formed the majority of their stars in the past $\sim 5$ Gyr to two-component systems characterized by the overlap of comparable numbers of old and young stars. Taken as an ensemble, star formation in nearby dwarfs dips to lower-than-average rates at intermediate times ($4<t$/Gyr $<8$), a feature that we trace in the simulation to the effects of cosmic reionization. Reionization heats the gas and drives it out of the shallow potential wells of low mass halos, affecting especially those below a sharp mass threshold that corresponds to a virial temperature of $\sim 2 \times 10^4 $ $\mathrm{K}$ at $z_{\rm reion}$. The loss of baryons leads to a sharp decline in the star forming activity of early-collapsing systems, which, compounded by feedback from early star formation, empties halos of gas and leaves behind systems where a single old stellar component prevails. In halos below the threshold at $z_{\rm reion}$, reionization heating leads to a delay in the onset of star formation that lasts until the halo grows massive enough to allow some of the remaining gas to cool and form stars. Young stellar components therefore dominate in dwarfs whose halos assemble late and thus form few stars before reionization. Two-component systems may be traced to late mergers of individual examples of the two aforementioned cases. The relative dearth of intermediate-age stars in nearby dwarfs might thus be the clearest signature yet identified of the imprint of cosmic reionization on the star formation history of dwarf galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The universal nature of subhalo accretion

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dark Matter in the Local Universe

We review how dark matter is distributed in our local neighbourhood from an observational and theoretical perspective. We will start by describing first the dark matter halo of our own galaxy and in the Local Group. Then we proceed to describe the dark matter distribution in the more extended area known as the Local Universe. Depending on the nature of dark matter, numerical simulations predict different abundances of substructures in Local Group galaxies, in the number of void regions and in the abundance of low rotational velocity galaxies in the Local Universe. By comparing these predictions with the most recent observations, strong constrains on the physical properties of the dark matter particles can be derived. We devote particular attention to the results from the Constrained Local UniversE Simulations (CLUES) project, a special set of simulations whose initial conditions are constrained by observational data from the Local Universe. The resulting simulations are designed to reproduce the observed structures in the nearby universe. The CLUES provides a numerical laboratory for simulating the Local Group of galaxies and exploring the physics of galaxy formation in an environment designed to follow the observed Local Universe. It has come of age as the numerical analogue of Near-Field Cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

N-body simulations with a cosmic vector for dark energy

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 $Ω_M$ and $σ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Orbital anisotropy in cosmological haloes revisited

The velocity anisotropy of particles inside dark matter (DM) haloes is an important physical quantity, which is required for the accurate modelling of mass profiles of galaxies and clusters of galaxies. It is typically measured using the ratio of the radial-to-tangential velocity dispersions at a given distance from the halo centre. However, this measure is insufficient to describe the dynamics of realistic haloes, which are typically quite elongated. Studying the velocity distribution in massive DM haloes in cosmological simulations, we find that in the inner parts of the haloes the local velocity ellipsoids are strongly aligned with the major axis of the halo, the alignment being stronger for more relaxed haloes. In the outer regions of the haloes, the alignment becomes gradually weaker and the orientation is more random. These two distinct regions of different degree of the alignment coincide with two characteristic regimes of the DM density profile: shallower and steeper than ρr^{-2}. This alignment of the local velocity ellipsoids requires reinterpretation of features found in measurements based on the spherically averaged ratio of the radial-to-tangential velocity dispersions. In particular, we show that the velocity distribution in the central halo regions is highly anisotropic. For cluster-size haloes with mass 10^{14}-10^{15} h^-1 Msun, the velocity anisotropy along the major axis is nearly independent of radius and is equal to β=1-σ^2_{perp}/σ^2_{radial}=0.4, which is significantly larger than the previously estimated spherically averaged velocity anisotropy. The alignment of density and velocity anisotropies, and the radial trends may also have some implications for the mass modelling based on kinematical data of such objects as galaxy clusters or dwarf spheroidals, where the orbital anisotropy is a key element in an unbiased mass inference.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulating Structure Formation of the Local Universe

In this work we present cosmological N-body simulations of the Local Universe with initial conditions constrained by the Two-Micron Redshift Survey (2MRS) within a cubic volume of 180 Mpc/h side-length centred at the Local Group. We use a self-consistent Bayesian based approach to explore the joint parameter space of primordial density fluctuations and peculiar velocity fields, which are compatible with the 2MRS galaxy distribution after cosmic evolution. This method (the KIGEN-code) includes the novel ALPT (Augmented Lagrangian Perturbation Theory) structure formation model which combines second order LPT (2LPT) on large scales with the spherical collapse model on small scales. Furthermore we describe coherent flows with 2LPT and include a dispersion term to model fingers-of-god (fogs). These implementations are crucial to avoid artificial filamentary structures, which appear when using a structure formation model with 2LPT and data with compressed fogs. We assume LCDM cosmology throughout our method. The recovered initial Gaussian fields are used to perform a set of 25 constrained simulations. Statistically this ensemble of simulations is in agreement with a reference set of 25 simulations based on randomly seeded Gaussian fluctuations in terms of matter statistics, power-spectra and mass functions. Considering the entire volume, we obtain correlation coefficients of about 98.3% for the comparison between the simulated density fields and the galaxy density field in log-space with Gaussian smoothing scales of r_S=3.5 Mpc/h (74% for r_S=1.4 Mpc/h). The cross power-spectra show correlations with the galaxy distribution which weaken towards smaller length scales until they vanish at scales of 2.2-3.0 Mpc/h. The simulations we present provide a fully nonlinear density and velocity field with a high level of correlation with the observed galaxy distribution at scales of a few Mpc.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulations of the Local Universe Constrained by Observational Peculiar Velocities

Peculiar velocities, obtained from direct distance measurements, are data of choice to achieve constrained simulations of the Local Universe reliable down to a scale of a few Megaparsecs. Unlike redshift surveys, peculiar velocities are direct tracers of the underlying gravitational field as they trace both baryonic and dark matter. This paper presents the first attempt to use solely observational peculiar velocities to constrain cosmological simulations of the nearby universe. In order to set up Initial Conditions, a Reverse Zel'dovich Approximation (RZA) is used to displace constraints from their positions at z=0 to their precursors' locations at higher redshifts. An additional new feature replaces original observed radial peculiar velocity vectors by their full 3D reconstructions provided by the Wiener-Filter (WF) estimator. Subsequently, the Constrained Realization of Gaussian fields technique (CR) is applied to build various realizations of the Initial Conditions. The WF/RZA/CR method is first tested on realistic mock catalogs built from a reference simulation similar to the Local Universe. These mocks include errors on peculiar velocities, on data-point positions and a large continuous zone devoid of data in order to mimic galactic extinction. Large scale structures are recovered with a typical accuracy of 5 Megaparsecs/h in position, the best realizations reaching a 2-3 Mpc/h precision, the limit imposed by the RZA linear theory. Then, the method is applied to the first observational radial peculiar velocity catalog of the project Cosmicflows. This paper is a proof of concept that the WF/RZA/CR method can be applied to observational peculiar velocities to successfully build constrained Initial Conditions.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Structure Finding in Cosmological Simulations: The State of Affairs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The clustering of galaxies at z~0.5 in the SDSS-III Data Release 9 BOSS-CMASS sample: a test for the LCDM cosmology

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The kinematics of the Local Group in a cosmological context

Recent observations constrained the tangential velocity of M31 with respect to the Milky Way (MW) to be v_tan<34.4 km/s and the radial velocity to be in the range v_rad=-109+/- 4.4 km/s (van der Marel et al. 2012). In this study we use a large volume high resolution N-body cosmological simulation (Bolshoi) together with three constrained simulations to statistically study this kinematics in the context of the LCDM. The comparison of the ensembles of simulated pairs with the observed LG at the 1-sigma level in the uncertainties has been done with respect to the radial and tangential velocities, the reduced orbital energy (e_tot), angular momentum (l_orb) and the dimensionless spin parameter, lambda. Our main results are: (i) the preferred radial and tangential velocities for pairs in LCDM are v_rad=-80+/-20 km/s, v_tan=50+/-10 km/s, (ii) pairs around that region are 3 to 13 times more common than pairs within the observational values, (iii) 15%to 24% of LG-like pairs in LCDM have energy and angular momentum consistent with observations while (iv) 9% to 13% of pairs in the same sample show similar values in the inferred dimensionless spin parameter. It follows that within current observational uncertainties the quasi-conserved quantities that characterize the orbit of the LG, i.e. e_tot, r_orb and lambda, do not challenge the standard LCDM model, but the model is in tension with regard to the actual values of the radial and tangential velocities. This might hint to a problem of the LCDM model to reproduce the observed LG.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Vector dark energy and high-z massive clusters

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cosmic Structure and Dynamics of the Local Universe

We present a cosmography analysis of the Local Universe based on the recently released Two-Micron All-Sky Redshift Survey (2MRS). Our method is based on a Bayesian Networks Machine Learning algorithm (the Kigen-code) which self-consistently samples the initial density fluctuations compatible with the observed galaxy distribution and a structure formation model given by second order Lagrangian perturbation theory (2LPT). From the initial conditions we obtain an ensemble of reconstructed density and peculiar velocity fields which characterize the local cosmic structure with high accuracy unveiling nonlinear structures like filaments and voids in detail. Coherent redshift space distortions are consistently corrected within 2LPT. From the ensemble of cross-correlations between the reconstructions and the galaxy field and the variance of the recovered density fields we find that our method is extremely accurate up to k ~ 1 h Mpc^-1 and still yields reliable results down to scales of about 3-4 h^-1 Mpc. The motion of the local group we obtain within ~ 80 h^-1 Mpc (v_LG=522+-86 km s^-1, l_LG=291^o +- 16^o, b_LG=34^o+-8^o) is in good agreement with measurements derived from the CMB and from direct observations of peculiar motions and is consistent with the predictions of LambdaCDM.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dwarf Galaxies and the Cosmic Web

We use a cosmological simulation of the formation of the Local Group of Galaxies to identify a mechanism that enables the removal of baryons from low-mass halos without appealing to feedback or reionization. As the Local Group forms, matter bound to it develops a network of filaments and pancakes. This moving web of gas and dark matter drifts and sweeps a large volume, overtaking many halos in the process. The dark matter content of these halos is unaffected but their gas can be efficiently removed by ram-pressure. The loss of gas is especially pronounced in low-mass halos due to their lower binding energy and has a dramatic effect on the star formation history of affected systems. This "cosmic web stripping" may help to explain the scarcity of dwarf galaxies compared with the numerous low-mass halos expected in ΛCDM and the large diversity of star formation histories and morphologies characteristic of faint galaxies. Although our results are based on a single high-resolution simulation, it is likely that the hydrodynamical interaction of dwarf galaxies with the cosmic web is a crucial ingredient so far missing from galaxy formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The MUSIC of Galaxy Clusters I: Baryon properties and Scaling Relations of the thermal Sunyaev-Zel'dovich Effect

We introduce the Marenostrum-MultiDark SImulations of galaxy Clusters (MUSIC) Dataset, one of the largest sample of hydrodynamically simulated galaxy clusters with more than 500 clusters and 2000 groups. The objects have been selected from two large N-body simulations and have been resimulated at high resolution using SPH together with relevant physical processes (cooling, UV photoionization, star formation and different feedback processes). We focus on the analysis of the baryon content (gas and star) of clusters in the MUSIC dataset both as a function of aperture radius and redshift. The results from our simulations are compared with the most recent observational estimates of the gas fraction in galaxy clusters at different overdensity radii. When the effects of cooling and stellar feedbacks are included, the MUSIC clusters show a good agreement with the most recent observed gas fractions quoted in the literature. A clear dependence of the gas fractions with the total cluster mass is also evident. The impact of the aperture radius choice, when comparing integrated quantities at different redshifts, is tested: the standard definition of radius at a fixed overdensity with respect to critical density is compared with a definition based on the redshift dependent overdensity with respect to background density. We also present a detailed analysis of the scaling relations of the thermal SZ (Sunyaev Zel'dovich) Effect derived from MUSIC clusters. The integrated SZ brightness, Y, is related to the cluster total mass, M, as well as, the M-Y counterpart, more suitable for observational applications. Both laws are consistent with predictions from the self-similar model, showing a very low scatter. The effects of the gas fraction on the Y-M scaling and the presence of a possible redshift dependence on the Y-M scaling relation are also explored.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Towards a panchromatic picture of galaxy evolution during the reionization epoch

There are thousands of confirmed detections of star forming galaxies at high redshift (z > 4). These observations rely primarily on the detection of the spectral Lyman Break and the Lyman-alpha emission line. Theoretical modelling of these sources helps to interpret the observations in the framework of the standard cosmological paradigm. We present results from the High-z MareNostrum Project, aimed at constructing a panchromatic picture of the high redshift galaxy evolution that will improve our understanding of young star forming galaxies. Our simulation successfully reproduces the observational constraints from Lyman Break Galaxies and Lyman-alpha emitters at 5 < z < 7 . Based on this model we make predictions on the expected Far Infrared (FIR) emission that should be observed for LAEs. These predictions will help to settle down the question on the dust content of massive high-z galaxies, an issue that will be feasible to probe observationally with the Atacama Large Millimetre Array (ALMA).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Haloes gone MAD: The Halo-Finder Comparison Project

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Renegade Subhaloes in the Local Group

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The MultiDark Database: Release of the Bolshoi and MultiDark Cosmological Simulations

We present the online MultiDark Database -- a Virtual Observatory-oriented, relational database for hosting various cosmological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation (www.multidark.org). The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250/h Mpc simulation box, 1/h kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000/h Mpc simulation box, 7/h kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z=0-15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merging trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constrained Local UniversE Simulations (CLUES)

The local universe is the best known part of our universe. Within the CLUES project (http://clues-project.org - Constrained Local UniversE Simulations) we perform numerical simulations of the evolution of the local universe. For these simulations we construct initial conditions based on observational data of the galaxy distribution in the local universe. Here we review the technique of these constrained simulations. In the second part we summarize our predictions of a possible Warm Dark Matter cosmology for the observed local distribution of galaxies and the local spectrum of mini-voids as well as a study of the satellite dynamics in a simulated Local Group.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dwarf Galaxies in Voids: Dark Matter Halos and Gas Cooling

Galaxy surveys have shown that luminous galaxies are mainly distributed in large filaments and galaxy clusters. The remaining large volumes are virtually devoid of luminous galaxies. This is in concordance with the formation of the large-scale structure in Universe as derived from cosmological simulations. However, the numerical results indicate that cosmological voids are abundantly populated with dark matter haloes which may in principle host dwarf galaxies. Observational efforts have in contrast revealed, that voids are apparently devoid of dwarf galaxies. We investigate the formation of dwarf galaxies in voids by hydrodynamical cosmological simulations. Due to the cosmic ultra-violet background radiation low-mass haloes show generally are reduced baryon fraction. We determine the characteristic mass below which dwarf galaxies are baryon deficient. We show that the circular velocity below which the accretion of baryons is suppressed is approximately 40 km/s. The suppressed baryon accretion is caused by the photo-heating due to the UV-background. We set up a spherical halo model and show that the effective equation of state of the gas in the periphery of dwarf galaxies determines the characteristic mass. This implies that any process which heats the gas around dwarf galaxies increases the characteristic mass and thus reduces the number of observable dwarf galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reionization of the Local Group of Galaxies

We present the first detailed structure formation and radiative transfer simulations of the reionization history of our cosmic neighbourhood. To this end, we follow the formation of the Local Group of galaxies and nearby clusters by means of constrained simulations, which use the available observational constraints to construct a representation of those structures which reproduces their actual positions and properties at the present time. We find that the reionization history of the Local Group is strongly dependent on the assumed photon production efficiencies of the ionizing sources, which are still poorly constrained. If sources are relatively efficient, i.e. the process is 'photon-rich', the Local Group is primarily ionized externally by the nearby clusters. Alternatively, if the sources are inefficient, i.e. reionization is 'photon-poor' the Local Group evolves largely isolated and reionizes itself. The mode of reionization, external vs. internal, has important implications for the evolution of our neighbourhood, in terms of e.g. its satellite galaxy populations and primordial stellar populations. This therefore provides an important avenue for understanding the young universe by detailed studies of our nearby structures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Strong lensing in the MareNostrum Universe: biases in the cluster lens population

Strong lensing is one of the most direct probes of the mass distribution in the inner regions of galaxy clusters. It can be used to constrain the density profiles and to measure the mass of the lenses. Moreover, the abundance of strong lensing events can be used to constrain the structure formation and the cosmological parameters through the so-called "arc-statistics" approach. However, several issues related to the usage of strong lensing clusters in cosmological applications are still controversial, leading to the suspect that several biases may affect this very peculiar class of objects. With this study we aim at better understanding the properties of galaxy clusters which can potentially act as strong lenses. We do so by investigating the properties of a large sample of galaxy clusters extracted from the N-body/hydrodynamical simulation MareNostrum Universe. We explore the correlation between the cross section for lensing and many properties of clusters, like the mass, the three-dimensional and projected shapes, their concentrations, the X-ray luminosity and the dynamical activity. We find that the probability of strong alignments between the major axes of the lenses and the line of sight is a growing function of the lensing cross section. In projection, the strong lenses appear rounder within R200, but we find that their cores tend to be more elliptical as the lensing cross section increases. We also find that the cluster concentrations estimated from the projected density profiles tend to be biased high. The X-ray luminosity of strong lensing clusters is higher than that of normal lenses of similar mass and redshift. This is particular significant for the least massive lenses. Finally, we find that the strongest lenses generally exhibit an excess of kinetic energy within the virial radius, indicating that they are more dynamically active than usual clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The grouping, merging and survival of subhaloes in the simulated Local Group

We use a simulation performed within the Constrained Local UniversE Simulation (CLUES) project to study a realistic Local Group-like object. We employ this group as a numerical laboratory for studying the evolution of the population of its subhaloes from the point of view of the effects it may have on the origin of different types of dwarf galaxies. We focus on the processes of tidal stripping of the satellites, their interaction, merging and grouping before infall. The tidal stripping manifests itself in the transition between the phase of mass accretion and mass loss seen in most subhaloes, which occurs at the moment of infall on to the host halo, and the change of the shape of their mass function with redshift. Although the satellites often form groups, they are loosely bound within them and do not interact with each other. The infall of a large group could however explain the observed peculiar distribution of the Local Group satellites, but only if it occurred recently. Mergers between prospective subhaloes are significant only during an early stage of evolution, i.e. more than 7 Gyr ago, when they are still outside the host haloes. Such events could thus contribute to the formation of more distant early type Milky Way companions. Once the subhaloes enter the host halo the mergers become very rare.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The impact of baryonic physics on the shape and radial alignment of substructures in cosmological dark matter haloes

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The luminosities of backsplash galaxies in constrained simulations of the Local Group

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The preferred direction of infalling satellite galaxies in the Local Group

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Constrained simulations of the Local Group: on the radial distribution of substructures

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The mass and anisotropy profiles of galaxy clusters from the projected phase space density: testing the method on simulated data

We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase space density which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters we analyze mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard LCDM model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 percent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Is WMAP3 normalization compatible with the X-Ray cluster abundance?

We present the mass and X-ray temperature functions derived from a sample of more than 15,000 galaxy clusters of the MareNostrum Universe cosmological SPH simulations. In these simulations, we follow structure formation in a cubic volume of 500/h Mpc on a side assuming cosmological parameters consistent with either the first or third year WMAP data and gaussian initial conditions. We compare our numerical predictions with the most recent observational estimates of the cluster X-ray temperature functions and find that the low normalization cosmological model inferred from the 3 year WMAP data results is barely compatible with the present epoch X-ray cluster abundances. We can only reconcile the simulations with the observational data if we assume a normalization of the Mass-Temperature relation which is a factor of 2.5--3 smaller than our non-radiative simulations predict. This deviation seems to be too large to be accounted by the effects of star formation or cooling in the ICM, not taken into account in these simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Supersonic Motions of Galaxies in Clusters

We study motions of galaxies in galaxy clusters formed in the concordance LCDM cosmology. We use high-resolution cosmological simulations that follow dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing motions of galaxies and the properties of intracluster gas in the sample of eight simulated clusters at z=0, we study velocity dispersion profiles of the dark matter, gas, and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ~1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move faster than the dark matter particles in clusters. The magnitude of the velocity bias, b_v~1.1, is, however, smaller than the bias estimated for subhalos in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ~20-30% of the gas sound speed. These random motions provide about 10% of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Optimizing Higher-Order Lagrangian Perturbation Theory for Standard CDM and BSI models

We investigate the performance of Lagrangian perturbation theory up to the second order for two scenarios of cosmological large-scale structure formation, SCDM (standard cold dark matter) and BSI (broken scale invariance). The latter model we study as a representative of COBE-normalized CDM models which fit the small-scale power of galaxy surveys. In this context we optimize the performance of the Lagrangian perturbation schemes by smoothing the small-scale fluctuations in the initial data. The results of the so obtained Lagrangian mappings are computed for a set of COBE-normalized SCDM and BSI initial data of different sizes and at different times. We compare these results against those obtained with a numerical PM-code. We find an excellent performance of the optimized Lagrangian schemes down to scales close to the correlation length. This is explained by the counterintuitive fact that nonlinearities in the model can produce more small-scale power, if initially such power is removed. The optimization scheme can be expressed in a way which is independent of the type of fluctuation spectrum and of the size of the simulations.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Confrontation of a Double Inflationary Cosmological Model with Observations

CDM models with non-scale-free step-like spectra of adiabatic perturbations produced in a realistic double inflationary model are compared with recent observational data. The model contains two additional free parameters relatively to the standard CDM model with the flat ($n=1$) initial spectrum. Results of the COBE experiment are used for the determination of a free overall spectrum normalization. Then predictions for the galaxy biasing parameter, the variance for "counts in cells", the galaxy angular correlation function, bulk flow peculiar velocities and the Mach number test are obtained. Also considered are conditions for galaxy and quasar formation. Observational data strongly restricts allowed values for the two remaining model parameters. However, a non-empty region for them satisfying all considered tests is found.

0 citations0 reviewsNo code linkedOpen access