Source author record

Stefan Gottlöber

Stefan Gottlöber 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

26works
2topics
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

26 published item(s)

preprint2022arXiv

A Machine Learning Approach to Correct for Mass Resolution Effects in Simulated Halo Clustering Statistics

The increase in the observed volume in cosmological surveys imposes various challenges on simulation preparations. Firstly, the volume of the simulations required increases proportionally to the observations. However, large-volume simulations are quickly becoming computationally intractable. Secondly, on-going and future large-volume survey are targeting smaller objects, e.g. emission line galaxies, compared to the earlier focus, i.e. luminous red galaxies. They require the simulations to have higher mass resolutions. In this work we present a machine learning (ML) approach to calibrate the halo catalogue of a low-resolution (LR) simulation by training with a paired high-resolution (HR) simulation with the same background white noise, thus we can build the training data by matching HR haloes to LR haloes in a one-to-one fashion. After training, the calibrated LR halo catalogue reproduces the mass-clustering relation for mass down to $2.5\times 10^{11}~h^{-1}M_\odot$ within $5~{\rm per~cent}$ at scales $k<1~h\,\rm Mpc^{-1}$. We validate the performance of different statistics including halo mass function, power spectrum, two-point correlation function, and bispectrum in both real and redshift space. Our approach generates high-resolution-like halo catalogues ($>200$ particles per halo) from low-resolution catalogues ($>25$ particles per halo) containing corrected halo masses for each object. This allows to bypass the computational burden of a large-volume real high-resolution simulation without much compromise in the mass resolution of the result. The cost of our ML approach ($\sim 1$ CPU-hour) is negligible compared to the cost of a $N$-body simulation (e.g. millions of CPU-hours), The required computing time is cut a factor of 8.

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

Astraeus VII: The environmental-dependent assembly of galaxies in the Epoch of Reionization

Using the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dark matter simUlationS) framework, we explore the impact of environmental density and radiative feedback on the assembly of galaxies and their host halos during the Epoch of Reionization. The ASTRAEUS framework allows us to study the evolution of galaxies with masses ($\rm 10^{8.2}M_\odot < M_{\rm h} < 10^{13}M_\odot$) in wide variety of environment ($-0.5 < {\rm log}(1+δ) < 1.3$ averaged over $(2~{\rm cMpc})^3$). We find that : (i) there exists a mass- and redshift- dependent "characteristic" environment (${\rm log} (1+δ_a(M_{\rm h}, z)) = 0.021\times (M_{\rm h}/M_\odot)^{0.16} + 0.07 z -1.12$, up to $z\sim 10$) at which galaxies are most efficient at accreting dark matter, e.g at a rate of $0.2\%$ of their mass every Myr at $z=5$; (ii) the number of minor and major mergers and their contributions to the dark matter assembly increases with halo mass at all redshifts and is mostly independent of the environment; (iii) at $z=5$ minor mergers contribute slightly more (by up to $\sim 10\%$) to the dark matter assembly while for the stellar assembly, major mergers dominate the contribution from minor mergers for $M_{\rm h}\lesssim 10^{11.5}M_\odot$ galaxies; (iv) radiative feedback quenches star formation more in low-mass galaxies ($M_{\rm h} \lesssim 10^{9.5}M_\odot$) in over-dense environments (${\rm log}(1+δ) > 0.5$); dominated by their major branch, this yields star formation histories biased towards older ages with a slower redshift evolution.

preprint2022arXiv

Numerical simulations of bar formation in the Local Group

More than 50 per cent of present-day massive disc galaxies show a rotating stellar bar. Their formation and dynamics have been widely studied both numerically and observationally. Although numerical simulations in the $Λ$CDM cosmological framework predict the formation of such stellar components, there seems to be a tension between theoretical and observational results. Simulated bars are typically larger in size and have slower pattern speed than observed ones. We study the formation and evolution of barred galaxies, using two $Λ$CDM zoom-in hydrodynamical simulations of the CLUES project that follow the evolution of a cosmological Local Group-like volume. We found that our simulated bars, at $z = 0$, are both shorter and faster rotators than previous ones found in other studies on cosmological simulations alleviating the tension mentioned above. These bars match the short tail-end of the observed bar length distribution. In agreement with previous numerical works, we find that bars form in those systems where the disc self-gravity is dominant over the dark matter halo, making them unstable against bar formation. Our bars developed in the last 3-4 Gyr until they achieve their current length and strength; as bars grow, their lengths increase while their rotation speeds decrease. Despite this slowdown, at redshift $z = 0$ their rotation speeds and size match well the observational data.

preprint2022arXiv

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$).

preprint2020arXiv

Galactic ionising photon budget during the Epoch of Reionisation in the Cosmic Dawn II simulation

Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmission inside high-mass haloes is bi-modal -- ionized cells are transparent, while neutral cells absorb the photons their stars produce - and the halo escape fraction $f_{esc}$ reflects the balance of star formation rate ("SFR") between these modes. The latter is increasingly prevalent at higher halo mass, driving down $f_{esc}$ (we provide analytical fits to our results), whereas halo escape luminosity, proportional to $f_{esc} \times$SFR, increases with mass. Haloes with dark matter masses within $6.10^{8} M_\odot < M_h < 3.10^{10} M_\odot$ produce $\sim 80$% of the escaping photons at z=7, when the Universe is 50% ionized, making them the main drivers of cosmic reionization. Less massive haloes, though more numerous, have low SFRs and contribute less than 10% of the photon budget then, despite their high $f_{esc}$. High mass haloes are too few and too opaque, contributing $<10$% despite their high SFRs. The dominant mass range is lower (higher) at higher (lower) redshift, as mass function and reionization advance together (e.g. at z$=8.5$, x$_{\rm HI}=0.9$, $M_h < 5.10^9 M_\odot$ haloes contributed $\sim$80%). Galaxies with UV magnitudes $M_{AB1600}$ between $-12$ and $-19$ dominated reionization between z$=6$ and 8.

preprint2020arXiv

The Past History of Galaxy Clusters told by their present neighbors

Galaxy clusters can play a key role in modern cosmology provided their evolution is properly understood. However, observed clusters give us only a single timeframe of their dynamical state. Therefore, finding present observable data of clusters that are well correlated to their assembly history constitutes an inestimable tool for cosmology. Former studies correlating environmental descriptors of clusters to their formation history are dominated by halo mass - environment relations. This paper presents a mass-free correlation between the present neighbor distribution of cluster-size halos and the latter mass assembly history. From the Big Multidark simulation, we extract two large samples of random halos with masses ranging from Virgo to Coma cluster sizes. Additionally, to find the main environmental culprit for the formation history of the Virgo cluster, we compare the Virgo-size halos to 200 Virgo-like halos extracted from simulations that resemble the local Universe. The number of neighbors at different cluster-centric distances permits discriminating between clusters with different mass accretion histories. Similarly to Virgo-like halos, clusters with numerous neighbors within a distance of about 2 times their virial radius experience a transition at z~1 between an active period of mass accretion, relative to the mean, and a quiet history. On the contrary, clusters with few neighbors share an opposite trend: from passive to active assembly histories. Additionally, clusters with massive companions within about 4 times their virial radius tend to have recent active merging histories. Therefore, the radial distribution of cluster neighbors provides invaluable insights into the past history of these objects.

preprint2019arXiv

Clustering with JWST: Constraining Galaxy Host Halo Masses, Satellite Quenching Efficiencies, and Merger Rates at z=4-10

Galaxy clustering measurements can be used to constrain many aspects of galaxy evolution, including galaxy host halo masses, satellite quenching efficiencies, and merger rates. We simulate JWST galaxy clustering measurements at z$\sim$4$-$10 by utilizing mock galaxy samples produced by an empirical model, the UniverseMachine. We also adopt the survey footprints and typical depths of the planned joint NIRCam and NIRSpec Guaranteed Time Observation program planned for Cycle 1 to generate realistic JWST survey realizations and to model high-redshift galaxy selection completeness. We find that galaxy clustering will be measured with $\gtrsim$5$σ$ significance at z$\sim$4$-$10. Halo mass precisions resulting from Cycle 1 angular clustering measurements will be $\sim$0.2 dex for faint (-18 $\gtrsim$ $\mathit{M}_{\mathrm{UV}}^{ }$ $\gtrsim$ -19) galaxies at z$\sim$4$-$10 as well as $\sim$0.3 dex for bright ($\mathit{M}_{\mathrm{UV}}^{ }$ $\sim$ -20) galaxies at z$\sim$4$-$7. Dedicated spectroscopic follow-up over $\sim$150 arcmin$^2$ would improve these precisions by $\sim$0.1 dex by removing chance projections and low-redshift contaminants. Future JWST observations will therefore provide the first constraints on the stellar-halo mass relation in the epoch of reionization and substantially clarify how this relation evolves at z$>$4. We also find that $\sim$1000 individual satellites will be identifiable at z$\sim$4$-$8 with JWST, enabling strong tests of satellite quenching evolution beyond currently available data (z$\lesssim$2). Finally, we find that JWST observations can measure the evolution of galaxy major merger pair fractions at z$\sim$4$-$8 with $\sim$0.1$-$0.2 dex uncertainties. Such measurements would help determine the relative role of mergers to the build-up of stellar mass into the epoch of reionization.

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.

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

preprint2016arXiv

Dynamical signatures of infall around galaxy clusters: a generalized Jeans equation

We study the internal kinematics of galaxy clusters in the region beyond the sphere of virialization. Galaxies around a virialized cluster are infalling towards the cluster centre with a non-zero mean radial velocity. We develop a new formalism for describing the dynamical state of clusters, by generalizing the standard Jeans formalism with the inclusion of the peculiar infall motions of galaxies and the Hubble expansion as well as the contributions due to background cosmology. Using empirical fits to the radial profiles of density, mean radial velocity and velocity anisotropy of both a stacked cluster-mass halo and two isolated halos of a cosmological dark matter only simulation, we verify that our generalized Jeans equation correctly predicts the radial velocity dispersion out to 4 virial radii. We find that the radial velocity dispersion inferred from the standard Jeans equation is accurate up to 2 virial radii, but overestimated by \approx 20% for the stacked halo and by \approx 40% for the isolated halos, in the range \approx 2-3 virial radii. Our model depends on the logarithmic growth rate of the virial radius (function of halo mass or concentration), which we estimate in 7 different ways, and on the departure from self-similarity of the evolution of the peculiar velocity profile in virial units.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Modeling the clustering and halo occupation distribution of BOSS-CMASS galaxies in the Final Data Release

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $Λ$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$σ$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$σ$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.

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

The distribution of mass components in simulated disc galaxies

Using 22 hydrodynamical simulated galaxies in a LCDM cosmological context we recover not only the observed baryonic Tully-Fisher relation, but also the observed "mass discrepancy--acceleration" relation, which reflects the distribution of the main components of the galaxies throughout their disks. This implies that the simulations, which span the range 52 < V$_{\rm flat}$ < 222 km/s where V$_{\rm flat}$ is the circular velocity at the flat part of the rotation curve, and match galaxy scaling relations, are able to recover the observed relations between the distributions of stars, gas and dark matter over the radial range for which we have observational rotation curve data. Furthermore, we explicitly match the observed baryonic to halo mass relation for the first time with simulated galaxies. We discuss our results in the context of the baryon cycle that is inherent in these simulations, and with regards to the effect of baryonic processes on the distribution of dark matter.

preprint2014arXiv

Ram pressure statistics for bent tail radio galaxies

In this paper we use the MareNostrum Universe Simulation, a large scale, hydrodynamic, non-radiative simulation in combination with a simple abundance matching approach to determine the ram pressure statistics for bent radio sources (BRSs). The abundance matching approach allows us to determine the locations of all galaxies with stellar masses $> 10^{11} MSol$ in the simulation volume. Assuming ram pressure exceeding a critical value causes bent morphology, we compute the ratio of all galaxies exceeding the ram pressure limit (RPEX galaxies) relative to all galaxies in our sample. According to our model 50% of the RPEX galaxies at $z = 0$ are found in clusters with masses larger than $10^{14.5}MSol$ the other half resides in lower mass clusters. Therefore, the appearance of bent tail morphology alone does not put tight constraints on the host cluster mass. In low mass clusters, $M < 10^{14}MSol$, RPEX galaxies are confined to the central 500 kpc whereas in clusters of $> 10^{15}Msol$ they can be found at distances up to 1.5Mpc. Only clusters with masses $> 10^{15}MSol $ are likely to host more than one BRS. Both criteria may prove useful in the search for distant, high mass clusters.

preprint2014arXiv

Supermassive Black Hole Binaries: Environment and Galaxy Host Properties of PTA and eLISA sources

Supermassive black hole (BH) binaries would comprise the strongest sources of gravitational waves (GW) once they reach <<1 pc separations, for both pulsar timing arrays (PTAs) and space based (SB) detectors. While BH binaries coalescences constitute a natural outcome of the cosmological standard model and galaxy mergers, their dynamical evolution is still poorly understood and therefore their abundances at different stages. We use a dynamical model for the decay of BH binaries coupled with a cosmological simulation and semi-empirical approaches to the occupation of haloes by galaxies and BHs, in order to follow the evolution of the properties distribution of galaxies hosting BH binaries candidates to decay due to GWs emission. Our models allow us to relax simplifying hypothesis about the binaries occupation in galaxies and their mass, as well as redshift evolution. Following previously proposed electromagnetic (EM) signatures of binaries in the subpc regime, that include spectral features and variability, we model possible distributions of such signatures and also set upper limits to their lifespan. We found a bimodal distribution of hosts properties, corresponding to BH binaries suitable to be detected by PTA and the ones detectable only from space missions, as eLISA. Although it has been discussed that the peak of eLISA sources may happen at high z, we show that there must be a population of such sources in the nearby Universe that might show detectable EM signatures, representing an important laboratory for multimessenger astrophysics. We found a weak dependence of galaxy host properties on the binaries occupation, that can be traced back to the BH origin. The combination of the host correlations reported here with the expected EM signal, may be helpful to verify the presence of nearby GW candidates, and to distinguish them from 'regular' intrinsic AGN variability.

preprint2014arXiv

The MUSIC of Galaxy Clusters II: X-ray global properties and scaling relations

We present the X-ray properties and scaling relations of a large sample of clusters extracted from the Marenostrum MUltidark SImulations of galaxy Clusters (MUSIC) dataset. We focus on a sub-sample of 179 clusters at redshift z~0.11, with 3.2e14M_sun/h<M_vir<2e15Msun/h, complete in mass. We employed the X-ray photon simulator PHOX to obtain synthetic Chandra Observations and derive observable-like global properties of the intracluster medium (ICM), as X-ray temperature (T_X) and luminosity (L_X). T_X is found to slightly under-estimate the true mass-weighted temperature, although tracing fairly well the cluster total mass. We also study the effects of T_X on scaling relations with cluster intrinsic properties: total (M_500) and gas (M_g500) mass; integrated Compton parameter (Y_SZ) of the Sunyaev-Zel'dovich (SZ) thermal effect; Y_X=M_g500 T_X. We confirm that Y_X is a very good mass proxy, with a scatter on M_500-Y_X and Y_SZ-Y_X lower than 5%. The study of scaling relations among X-ray, intrinsic and SZ properties indicates that MUSIC clusters reasonably resemble the self-similar prediction, especially for correlations involving T_X. The observational approach also allows for a more direct comparison with real clusters, from which we find deviations mainly due to the physical description of the ICM, affecting T_X and, particularly, L_X.

preprint2014arXiv

The velocity shear and vorticity across redshifts and non-linear scales

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

preprint2013arXiv

The MUSIC of Galaxy Clusters III: Properties, evolution and Y-M scaling relation of protoclusters of galaxies

In this work we study the properties of protoclusters of galaxies by employing the MUSIC set of hydrodynamical simulations, featuring a mass-limited sample of 282 resimulated clusters with available merger trees up to high redshift, and we trace the cluster formation back to $z$ = 1.5, 2.3 and 4. We study the features and redshift evolution of the mass and the spatial distribution for all the cluster progenitors and for the protoclusters, which we define as the most massive progenitors of the clusters identified at $z$ = 0. A natural extension to redshifts larger than 1 is applied to the estimate of the baryon content also in terms of gas and stars budgets: no remarkable variations with redshift are discovered. Furthermore, motivated by the proven potential of Sunyaev-Zel'dovich surveys to blindly search for faint distant objects, we focus on the scaling relation between total object mass and integrated Compton $y$-parameter, and we check for the possibility to extend the mass-observable paradigm to the protocluster regime, far beyond the redshift of 1, to account for the properties of the simulated objects. We find that the slope of this scaling law is steeper than what expected for a self-similarity assumption among these objects, and it increases with redshift mainly for the synthetic clusters where radiative processes, such as radiative cooling, heating processes of the gas due to UV background, star formation and supernovae feedback, are included. We use three different criteria to account for the dynamical state of the protoclusters, and find no significant dependence of the scaling parameters from the level of relaxation. Based on this, we exclude that the dynamical state is the cause of the observed deviations from self-similarity.

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

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

Disentangling the dark matter halo from the stellar halo

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

preprint2011arXiv

Large scale environmental bias of the QSO line of sight proximity effect

We analyse the proximity zone of the intergalactic matter around high-redshift quasars in a cosmological environment. In a box of 64 h-1 Mpc base length we employ dark matter only simulations. For estimating the hydrogen temperature and density distribution we use the effective equation of state. Hydrogen is assumed to be in photoionisation equilibrium with a model background flux which is fit to recent observations of the mean optical depth and transmission flux statistics. At redshifts z = 3, 4, and 4.8, we select model quasar positions at the centre of the 20 most massive halos and 100 less massive halos identified in the simulation. From each assumed quasar position we cast 100 random lines of sight for two box length including the changes in the ionisation fractions by the QSO flux field and derive mock Lyα spectra. The proximity effect describes the dependence of the mean normalised optical depth ξ = τeff, QSO/τeff, Lyα as a function of the ratio of the ionisation rate by the QSO and the background field, ω = ΓQSO/ΓUVB, i.e. the profile ξ = (1 + ω/a)-0.5, where a strength parameter a is introduced. The strength parameter measures the deviation from the theoretical background model and is used to quantify any influence of the environmental density field. We reproduce an unbiased measurement of the proximity effect which is not affected by the host halo mass. The scatter between different lines of sight and different quasar host positions increases with decreasing redshift. Around the host halos, we find only a slight average overdensity in the proximity zone at comoving radii of 1 < rc < 10h-1 Mpc. However, a clear power-law correlation of the strength parameter with the average overdensity in rc is found, showing an overestimation of the ionising background in overdense regions and an underestimation in underdense regions.

preprint2011arXiv

The overdensity and masses of the friends-of-friends halos and universality of the halo mass function

The friends-of-friends algorithm (hereafter, FOF) is a percolation algorithm which is routinely used to identify dark matter halos from N-body simulations. We use results from percolation theory to show that the boundary of FOF halos does not correspond to a single density threshold but to a range of densities close to a critical value that depends upon the linking length parameter, b. We show that for the commonly used choice of b = 0.2, this critical density is equal to 81.62 times the mean matter density. Consequently, halos identified by the FOF algorithm enclose an average overdensity which depends on their density profile (concentration) and therefore changes with halo mass contrary to the popular belief that the average overdensity is ~180. We derive an analytical expression for the overdensity as a function of the linking length parameter b and the concentration of the halo. Results of tests carried out using simulated and actual FOF halos identified in cosmological simulations show excellent agreement with our analytical prediction. We also find that the mass of the halo that the FOF algorithm selects crucially depends upon mass resolution. We find a percolation theory motivated formula that is able to accurately correct for the dependence on number of particles for the mock realizations of spherical and triaxial Navarro-Frenk-White halos. However, we show that this correction breaks down when applied to the real cosmological FOF halos due to presence of substructures. Given that abundance of substructure depends on redshift and cosmology, we expect that the resolution effects due to substructure on the FOF mass and halo mass function will also depend on redshift and cosmology and will be difficult to correct for in general. Finally, we discuss the implications of our results for the universality of the mass function.

preprint2009arXiv

Constrained simulations of the local universe: II. The nature of the local Hubble flow

Using a suite of N-body simulations in different Cold Dark Matter (CDM) scenarios, with cosmological constant (\LCDM) and without (OCDM, SCDM), we study the Hubble flow (\sigh) in Local Volumes (LV) around Local Group (LG) like objects found in these simulations, and compare the numerical results with the most recent observations. We show that \LCDM and OCDM models exhibit the same behavior of \sigh. Hence, we demonstrate that the observed coldness of the Hubble flow is not likely to be a manifestation of the dark energy, contrary to previous claims. The coldness does not constitute a problem by itself but it poses a problem to the standard \LCDM model only if the mean density within the Local Volume is greater than twice the mean matter cosmic density. The lack of blueshifted galaxies in the LV, outside of the LG can be considered as another manifestation of the coldness of the flow. Finally, we show that the main dynamical parameter that affects the coldness of the flow is the relative isolation of the LG, and the absence of nearby Milky Way like objects within a distance of about $3\mpc$.