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

Oliver Hahn

Oliver Hahn 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-ph.IMastro-phquant-phastro-ph.HEhep-phMachine Learning

Catalog footprint

What is connected

39works
8topics
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

39 published item(s)

preprint2025arXiv

Accurate $N$-body simulations with local Primordial non-Gaussianities: initial conditions and aliasing

New generation galaxy surveys targeting constraints on local primordial non-Gaussianity (PNG) demand $N$-body simulations that accurately reproduce its effects. In this work, we explore various prescriptions for the initial conditions of simulations with PNG, aiming to optimise accuracy and minimise numerical errors, particularly due to aliasing. We have used $186$ runs that vary the starting redshift, LPT order, and non-Gaussianities ($f^{\rm local}_{\rm NL}$ and $g^{\rm local}_{\rm NL}$). Starting with $3$LPT at a redshift as low as $z_{\rm ini}\simeq 11.5$ reproduces to $<1 \%$ the power spectrum, bispectrum and halo mass function of a high-resolution reference simulation. The aliasing induced by the PNG terms in the power spectrum produces a $ \leq 3 \%$ excess small-scale power at the initial conditions, dropping below $0.1\%$ by $z=0$. State-of-the-art initial condition generators show a sub-percent agreement. We show that initial conditions for simulations with PNG should be established at a lower redshift using higher-order LPT schemes. We also show that removing the PNG aliasing signal is unnecessary for current simulations. The methodology proposed here can accelerate the generation of simulations with PNG while enhancing their accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

A deep learning framework for jointly extracting spectra and source-count distributions in astronomy

Astronomical observations typically provide three-dimensional maps, encoding the distribution of the observed flux in (1) the two angles of the celestial sphere and (2) energy/frequency. An important task regarding such maps is to statistically characterize populations of point sources too dim to be individually detected. As the properties of a single dim source will be poorly constrained, instead one commonly studies the population as a whole, inferring a source-count distribution (SCD) that describes the number density of sources as a function of their brightness. Statistical and machine learning methods for recovering SCDs exist; however, they typically entirely neglect spectral information associated with the energy distribution of the flux. We present a deep learning framework able to jointly reconstruct the spectra of different emission components and the SCD of point-source populations. In a proof-of-concept example, we show that our method accurately extracts even complex-shaped spectra and SCDs from simulated maps.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Renormalization group and UV completion of cosmological perturbations: Gravitational collapse as a critical phenomenon

Cosmological perturbation theory is known to converge poorly for predicting the spherical collapse and void evolution of collisionless matter. Using the exact parametric solution as a testing ground, we develop two asymptotic methods in spherical symmetry that resolve the gravitational evolution to much higher accuracy than Lagrangian perturbation theory (LPT), which is the current gold standard in the literature. One of the methods selects a stable fixed-point solution of the renormalization-group flow equation, thereby predicting already at the leading order the critical exponent of the phase transition to collapsed structures. The other method completes the truncated LPT series far into the UV regime, by adding a non-analytic term that captures the critical nature of the gravitational collapse. We find that the UV method most accurately resolves the evolution of the nonlinear density as well as its one-point probability distribution function. Similarly accurate predictions are achieved with the renormalization-group method, especially when paired with Padé approximants. Further, our results yield new, very accurate, formulas to relate linear and nonlinear density contrasts. Finally, we chart possible ways on how to adapt our methods to the case of cosmological random field initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Analytical growth functions for cosmic structures in a $Λ$CDM Universe

The cosmological fluid equations describe the early gravitational dynamics of cold dark matter (CDM), exposed to a uniform component of dark energy, the cosmological constant $Λ$. Perturbative predictions for the fluid equations typically assume that the impact of $Λ$ on CDM can be encapsulated by a refined growing mode $D$ of linear density fluctuations. Here we solve, to arbitrary high perturbative orders, the nonlinear fluid equations with an {\it Ansatz} for the fluid variables in increasing powers of $D$. We show that $Λ$ begins to populate the solutions starting at the fifth order in this strict $D$-expansion. By applying suitable resummation techniques, we recast these solutions to a standard perturbative series where not $D$, but essentially the initial gravitational potential serves as the bookkeeping parameter within the expansion. Then, by using the refined growth functions at second and third order in standard perturbation theory, we determine the matter power spectrum to one-loop accuracy as well as the leading-order contribution to the matter bispectrum. We find that employing our refined growth functions impacts the total power- and bispectra at a precision that is below one percent at late times. However, for the power spectrum, we find a characteristic scale-dependent suppression that is fairly similar to what is observed in massive neutrino cosmologies. Therefore, we recommend employing our refined growth functions in order to reduce theoretical uncertainties for analysing data in related pipelines.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deterministic Gaussian conversion protocols for non-Gaussian single-mode resources

In the context of quantum technologies over continuous variables, Gaussian states and operations are typically regarded as freely available, as they are relatively easily accessible experimentally. In contrast, the generation of non-Gaussian states, as well as the implementation of non-Gaussian operations, pose significant challenges. This divide has motivated the introduction of resource theories of non-Gaussianity. As for any resource theory, it is of practical relevance to identify free conversion protocols between resources, namely Gaussian conversion protocols between non-Gaussian states. Via systematic numerical investigations, we address the approximate conversion between experimentally relevant single-mode non-Gaussian states via arbitrary deterministic one-to-one mode Gaussian maps. First, we show that cat and binomial states are approximately equivalent for finite energy, while this equivalence was previously known only in the infinite-energy limit. Then we consider the generation of cat states from photon-added and photon-subtracted squeezed states, improving over known schemes by introducing additional squeezing operations. The numerical tools that we develop also allow to devise conversions of trisqueezed into cubic-phase states beyond previously reported performances. Finally, we identify various other conversions which instead are not viable.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantifying Qubit Magic Resource with Gottesman-Kitaev-Preskill Encoding

Quantum resource theories are a powerful framework to characterize and quantify relevant quantum phenomena and identify processes that optimize their use for different tasks. Here, we define a resource measure for magic, the sought-after property in most fault-tolerant quantum computers. In contrast to previous literature, our formulation is based on bosonic codes, well-studied tools in continuous-variable quantum computation. Particularly, we use the Gottesman-Kitaev-Preskill code to represent multi-qubit states and consider the resource theory for the Wigner negativity. Our techniques are useful to find resource lower bounds for different applications as state conversion and general unitary synthesis, in which measurements, auxiliary states, and classical feed-forward are allowed. The analytical expression of our magic measure allows us to extend current analysis limited to small dimensions, easily addressing systems of up to 12 qubits.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large-scale dark matter simulations

We review the field of collisionless numerical simulations for the large-scale structure of the Universe. We start by providing the main set of equations solved by these simulations and their connection with General Relativity. We then recap the relevant numerical approaches: discretization of the phase-space distribution (focusing on N-body but including alternatives, e.g., Lagrangian submanifold and Schrödinger-Poisson) and the respective techniques for their time evolution and force calculation (Direct summation, mesh techniques, and hierarchical tree methods). We pay attention to the creation of initial conditions and the connection with Lagrangian Perturbation Theory. We then discuss the possible alternatives in terms of the micro-physical properties of dark matter (e.g., neutralinos, warm dark matter, QCD axions, Bose-Einstein condensates, and primordial black holes), and extensions to account for multiple fluids (baryons and neutrinos), primordial non-Gaussianity and modified gravity. We continue by discussing challenges involved in achieving highly accurate predictions. A key aspect of cosmological simulations is the connection to cosmological observables, we discuss various techniques in this regard: structure finding, galaxy formation and baryonic modelling, the creation of emulators and light-cones, and the role of machine learning. We finalise with a recount of state-of-the-art large-scale simulations and conclude with an outlook for the next decade.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A sea of tides: star formation and the central-satellite dichotomy in a continuum of tidal environments

The environment-dependent bimodality of the distribution of stellar mass ($M_\ast$) and specific star formation rate (sSFR) of galaxies, and its explanation in terms of the central-satellite dichotomy, form a cornerstone of our current understanding of galaxy evolution in the hierarchical structure formation paradigm. We revisit this framework in the IllustrisTNG simulation in the context of the most extreme local tidal anisotropy $α_{\rm peak}$ experienced by each galaxy over cosmic time, which is an excellent proxy for environmental influence. We show that, while sharing a common monotonic $M_\ast$-$v_{\rm peak}$ relation, central, satellite and `splashback' galaxies define a hierarchy of increasing $α_{\rm peak}$. We also find that the sSFR of objects in small haloes unaffected by feedback from an active nucleus typically decreases with increasing $α_{\rm peak}$. Our results support an alternate viewpoint in which a galaxy can be identified by the value of $α_{\rm peak}$; i.e., rather than being placed on the central-satellite dichotomy, a galaxy is better classified by its location in a continuum of tidal environments. This conceptual shift can potentially yield a more robust understanding of galaxy evolution and the galaxy-dark matter connection, e.g., in accurately modelling subtle effects such as sSFR-induced secondary clustering.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accelerated orbital decay of supermassive black hole binaries in merging nuclear star clusters

The coalescence of supermassive black holes (SMBHs) should generate the strongest sources of gravitational waves (GWs) in the Universe. However, the dynamics of their coalescence is the subject of much debate. In this study, we use a suite of $N$-body simulations to follow the merger of two nuclear star clusters (NSCs), each hosting a SMBH in their centre. We find that the presence of distinct star clusters around each SMBH has important consequences for the dynamical evolution of the SMBH binary: (i) The separation between the SMBHs decreases by a few orders of magnitude in the first few Myrs by the combined effects of dynamical friction and a drag force caused by tidally stripped stars. In fact, this is a significant speedup for equal mass ratio binaries, and becomes extreme for unequal mass ratios, e.g. 1:10 or 1:100, which traditional dynamical friction alone would not permit to bind. (ii) The subsequent binary hardening is driven by the gravitational slingshots between the SMBH binary and stars, and also depends on the mass ratio between the SMBHs. Thus, with this additional drag force, we find that all SMBHs in our suite coalesce within a Hubble time. Given that about 50% of Milky Way sized galaxies host NSCs, our results are encouraging for upcoming GW observations with the Laser Interferometer Space Antenna -- LISA -- which will detect SMBH coalescence in the $10^4-10^7$$M_{\rm \odot}$ mass range.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring the Tidal Response of Structure Formation: Anisotropic Separate Universe Simulations using TreePM

We present anisotropic "separate universe" simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 7 h \text{Mpc}^{-1}$. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value at $z=0$, $R_K\approx 0.5$ for $k \geq 3 h\text{Mpc}^{-1}$ up to the scale where we consider our simulations reliable, $k \leq 7 h\text{Mpc}^{-1}$ at $z=0$. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms

At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Excursion set peaks: the role of shear

Recent analytical work on the modelling of dark halo abundances and clustering has demonstrated the advantages of combining the excursion set approach with peaks theory. We extend these ideas and introduce a model of excursion set peaks that incorporates the role of initial tidal effects or shear in determining the gravitational collapse of dark haloes. The model -- in which the critical density threshold for collapse depends on the tidal influences acting on protohaloes -- is well motivated from ellipsoidal collapse arguments and is also simple enough to be analytically tractable. We show that the predictions of this model are in very good agreement with measurements of the halo mass function and traditional scale dependent halo bias in N-body simulations across a wide range of masses and redshift. The presence of shear in the collapse threshold means that halo bias is naturally predicted to be nonlocal, and that protohalo densities at fixed mass are naturally predicted to have Lognormal-like distributions. We present the first direct estimate of Lagrangian nonlocal bias in N-body simulations, finding broad agreement with the model prediction. Finally, the simplicity of the model (which has essentially a single free parameter) opens the door to building efficient and accurate non-universal fitting functions of halo abundances and bias for use in precision cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

General relativistic 'screening' in cosmological simulations

We revisit the issue of interpreting the results of large volume cosmological simulations in the context of large scale general relativistic effects. We look for simple modifications to the nonlinear evolution of the gravitational potential $ψ$ that lead on large scales to the correct, fully relativistic description of density perturbations in the Newtonian gauge. We note that the relativistic constraint equation for $ψ$ can be cast as a diffusion equation, with a diffusion length scale determined by the expansion of the Universe. Exploiting the weak time evolution of $ψ$ in all regimes of interest, this equation can be further accurately approximated as a Helmholtz equation, with an effective relativistic 'screening' scale $\ell$ related to the Hubble radius. We demonstrate that it is thus possible to carry out N-body simulations in the Newtonian gauge by replacing Poisson's equation with this Helmholtz equation, involving a trivial change in the Green's function kernel. Our results also motivate a simple, approximate (but very accurate) gauge transformation - $δ_{\rm N}(\mathbf{k}) \approx δ_{\rm sim}(\mathbf{k})\times (k^2+\ell^{-2})/k^2$ - to convert the density field $δ_{\rm sim}$ of standard collisionless N-body simulations (initialised in the comoving synchronous gauge) into the Newtonian gauge density $δ_{\rm N}$ at arbitrary times. A similar conversion can also be written in terms of particle positions. Our results can be interpreted in terms of a Jeans stability criterion induced by the expansion of the Universe. The appearance of the screening scale $\ell$ in the evolution of $ψ$, in particular, leads to a natural resolution of the 'Jeans swindle' in the presence of super-horizon modes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

RHAPSODY-G simulations II - Baryonic growth and metal enrichment in massive galaxy clusters

We study the evolution of the stellar component and the metallicity of both the intracluster medium and of stars in massive ($M_{\rm vir}\approx 6\times 10^{14}$ M$_{\odot}/h$) simulated galaxy clusters from the Rhapsody-G suite in detail and compare them to observational results. The simulations were performed with the AMR code RAMSES and include the effect of AGN feedback at the sub-grid level. AGN feedback is required to produce realistic galaxy and cluster properties and plays a role in mixing material in the central regions and regulating star formation in the central galaxy. In both our low and high resolution runs with fiducial stellar yields, we find that stellar and ICM metallicities are a factor of two lower than in observations. We find that cool core clusters exhibit steeper metallicity gradients than non-cool core clusters, in qualitative agreement with observations. We verify that the ICM metallicities measured in the simulation can be explained by a simple "regulator" model in which the metallicity is set by a balance of stellar yield and gas accretion. It is plausible that a combination of higher resolution and higher metal yield in AMR simulation would allow the metallicity of simulated clusters to match observed values; however this hypothesis needs to be tested with future simulations. Comparison to recent literature highlights that results concerning the metallicity of clusters and cluster galaxies might depend sensitively on the scheme chosen to solve the hydrodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics

N-body simulations are essential for understanding the formation and evolution of structure in the Universe. However, the discrete nature of these simulations affects their accuracy when modelling collisionless systems. We introduce a new approach to simulate the gravitational evolution of cold collisionless fluids by solving the Vlasov-Poisson equations in terms of adaptively refineable "Lagrangian phase space elements". These geometrical elements are piecewise smooth maps between Lagrangian space and Eulerian phase space and approximate the continuum structure of the distribution function. They allow for dynamical adaptive splitting to accurately follow the evolution even in regions of very strong mixing. We discuss in detail various one-, two- and three-dimensional test problems to demonstrate the performance of our method. Its advantages compared to N-body algorithms are: i) explicit tracking of the fine-grained distribution function, ii) natural representation of caustics, iii) intrinsically smooth gravitational potential fields, thus iv) eliminating the need for any type of ad-hoc force softening. We show the potential of our method by simulating structure formation in a warm dark matter scenario. We discuss how spurious collisionality and large-scale discreteness noise of N-body methods are both strongly suppressed, which eliminates the artificial fragmentation of filaments. Therefore, we argue that our new approach improves on the N-body method when simulating self-gravitating cold and collisionless fluids, and is the first method that allows to explicitly follow the fine-grained evolution in six-dimensional phase space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Four phases of angular-momentum buildup in high-z galaxies: from cosmic-web streams through an extended ring to disc and bulge

We study the angular-momentum (AM) buildup in high-$z$ massive galaxies using high-resolution cosmological simulations. The AM originates in co-planar streams of cold gas and merging galaxies tracing cosmic-web filaments, and it undergoes four phases of evolution. (I) Outside the halo virial radius ($R_{\rm v}\!\sim\!100\,{\rm kpc}$), the elongated streams gain AM by tidal torques with a specific AM (sAM) $\sim\!1.7$ times the dark-matter (DM) spin due to the gas' higher quadrupole moment. This AM is expressed as stream impact parameters, from $\sim\!0.3R_{\rm v}$ to counter rotation. (II) In the outer halo, while the incoming DM mixes with the existing halo of lower sAM to a spin $λ_{\rm dm}\!\sim\!0.04$, the cold streams transport the AM to the inner halo such that their spin in the halo is $\sim\!3λ_{\rm dm}$. (III) Near pericenter, the streams dissipate into an irregular rotating ring extending to $\sim\!0.3R_{\rm v}$ and tilted relative to the inner disc. Torques exerted partly by the disc make the ring gas lose AM, spiral in, and settle into the disc within one orbit. The ring is observable with 30\% probability as a damped Lyman-$α$ absorber. (IV) Within the disc, $<\!0.1R_{\rm v}$, torques associated with violent disc instability drive AM out and baryons into a central bulge, while outflows remove low-spin gas, introducing certain sensitivity to feedback strength. Despite the different AM histories of gas and DM, the disc spin is comparable to the DM-halo spin. Counter rotation can strongly affect disc evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Rhapsody-G simulations: galaxy clusters as baryonic closed boxes and the covariance between hot gas and galaxies

Within a sufficiently large cosmic volume, conservation of baryons implies a simple `closed box' view in which the sum of the baryonic components must equal a constant fraction of the total enclosed mass. We present evidence from Rhapsody-G hydrodynamic simulations of massive galaxy clusters that the closed-box expectation may hold to a surprising degree within the interior, non-linear regions of haloes. At a fixed halo mass, we find a significant anti-correlation between hot gas mass fraction and galaxy mass fraction (cold gas + stars), with a rank correlation coefficient of -0.69 within $R_{500c}$. Because of this anti-correlation, the total baryon mass serves as a low-scatter proxy for total cluster mass. The fractional scatter of total baryon fraction scales approximately as $0.02 (Δ_c/100)^{0.6}$, while the scatter of either gas mass or stellar mass is larger in magnitude and declines more slowly with increasing radius. We discuss potential observational tests using cluster samples selected by optical and hot gas properties; the simulations suggest that joint selection on stellar and hot gas has potential to achieve 5% scatter in total halo mass.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Properties of Cosmic Velocity Fields

Understanding the velocity field is very important for modern cosmology: it gives insights to structure formation in general, and also its properties are crucial ingredients in modelling redshift-space distortions and in interpreting measurements of the kinetic Sunyaev-Zeldovich effect. Unfortunately, characterising the velocity field in cosmological N-body simulations is inherently complicated by two facts: i) The velocity field becomes manifestly multi-valued after shell-crossing and has discontinuities at caustics. This is due to the collisionless nature of dark matter. ii) N-body simulations sample the velocity field only at a set of discrete locations, with poor resolution in low-density regions. In this paper, we discuss how the associated problems can be circumvented by using a phase-space interpolation technique. This method provides extremely accurate estimates of the cosmic velocity fields and its derivatives, which can be properly defined without the need of the arbitrary "coarse-graining" procedure commonly used. We explore in detail the configuration-space properties of the cosmic velocity field on very large scales and in the highly nonlinear regime. In particular, we characterise the divergence and curl of the velocity field, present their one-point statistics, analyse the Fourier-space properties and provide fitting formulae for the velocity divergence bias relative to the non-linear matter power spectrum. We furthermore contrast some of the interesting differences in the velocity fields of warm and cold dark matter models. We anticipate that the high-precision measurements carried out here will help to understand in detail the dynamics of dark matter and the structures it forms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collisionless Dynamics and the Cosmic Web

I review the nature of three-dimensional collapse in the Zeldovich approximation, how it relates to the underlying nature of the three-dimensional Lagrangian manifold and naturally gives rise to a hierarchical structure formation scenario that progresses through collapse from voids to pancakes, filaments and then halos. I then discuss how variations of the Zeldovich approximation (based on the gravitational or the velocity potential) have been used to define classifications of the cosmic large-scale structure into dynamically distinct parts. Finally, I turn to recent efforts to devise new approaches relying on tessellations of the Lagrangian manifold to follow the fine-grained dynamics of the dark matter fluid into the highly non-linear regime and both extract the maximum amount of information from existing simulations as well as devise new simulation techniques for cold collisionless dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new approach to simulating collisionless dark matter fluids

Recently, we have shown how current cosmological N-body codes already follow the fine grained phase-space information of the dark matter fluid. Using a tetrahedral tesselation of the three-dimensional manifold that describes perfectly cold fluids in six-dimensional phase space, the phase-space distribution function can be followed throughout the simulation. This allows one to project the distribution function into configuration space to obtain highly accurate densities, velocities, and velocity dispersions. Here, we exploit this technique to show first steps on how to devise an improved particle-mesh technique. At its heart, the new method thus relies on a piecewise linear approximation of the phase space distribution function rather than the usual particle discretisation. We use pseudo-particles that approximate the masses of the tetrahedral cells up to quadrupolar order as the locations for cloud-in-cell (CIC) deposit instead of the particle locations themselves as in standard CIC deposit. We demonstrate that this modification already gives much improved stability and more accurate dynamics of the collisionless dark matter fluid at high force and low mass resolution. We demonstrate the validity and advantages of this method with various test problems as well as hot/warm-dark matter simulations which have been known to exhibit artificial fragmentation. This completely unphysical behaviour is much reduced in the new approach. The current limitations of our approach are discussed in detail and future improvements are outlined.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How closely do baryons follow dark matter on large scales?

We investigate the large-scale clustering and gravitational interaction of baryons and dark matter (DM) over cosmic time using a set of collisionless N-body simulations. Both components, baryons and DM, are evolved from distinct primordial density and velocity power spectra as predicted by early-universe physics. We first demonstrate that such two-component simulations require an unconventional match between force and mass resolution (i.e. force softening on at least the mean particle separation scale). Otherwise, the growth on any scale is not correctly recovered because of a spurious coupling between the two species at the smallest scales. With these simulations, we then demonstrate how the primordial differences in the clustering of baryons and DM are progressively diminished over time. In particular, we explicitly show how the BAO signature is damped in the spatial distribution of baryons and imprinted in that of DM. This is a rapid process, yet it is still not fully completed at low redshifts. On large scales, the overall shape of the correlation function of baryons and DM differs by 2% at z = 9 and by 0.2% at z = 0. The differences in the amplitude of the BAO peak are approximately a factor of 5 larger: 10% at z = 9 and 1% at z = 0. These discrepancies are, however, smaller than effects expected to be introduced by galaxy formation physics in both the shape of the power spectrum and in the BAO peak, and are thus unlikely to be detected given the precision of the next generation of galaxy surveys. Hence, our results validate the standard practice of modelling the observed galaxy distribution using predictions for the total mass clustering in the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How to Zoom: Bias, Contamination, and Lagrange Volumes in Multimass Cosmological Simulations

We perform a suite of multimass cosmological zoom simulations of individual dark matter halos and explore how to best select Lagrangian regions for resimulation without contaminating the halo of interest with low-resolution particles. Such contamination can lead to significant errors in the gas distribution of hydrodynamical simulations, as we show. For a fixed Lagrange volume, we find that the chance of contamination increases systematically with the level of zoom. In order to avoid contamination, the Lagrangian volume selected for resimulation must increase monotonically with the resolution difference between parent box and the zoom region. We provide a simple formula for selecting Lagrangian regions (in units of the halo virial volume) as a function of the level of zoom required. We also explore the degree to which a halo's Lagrangian volume correlates with other halo properties (concentration, spin, formation time, shape, etc.) and find no significant correlation. There is a mild correlation between Lagrange volume and environment, such that halos living in the most clustered regions have larger Lagrangian volumes. Nevertheless, selecting halos to be isolated is not the best way to ensure inexpensive zoom simulations. We explain how one can safely choose halos with the smallest Lagrangian volumes, which are the least expensive to resimulate, without biasing one's sample.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mergers and Mass Accretion for Infalling Halos Both End Well Outside Cluster Virial Radii

We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Population III Star Formation In Large Cosmological Simulations I. Halo Temporal and Physical Environment

We present a semi-analytic, computationally inexpensive model to identify halos capable of forming a Population III star in cosmological simulations across a wide range of times and environments. This allows for a much more complete and representative set of Population III star forming halos to be constructed, which will lead to Population III star formation simulations that more accurately reflect the diversity of Population III stars, both in time and halo mass. This model shows that Population III and chemically enriched stars coexist beyond the formation of the first generation of stars in a cosmological simulation until at least z~10, and likely beyond, though Population III stars form at rates that are 4-6 orders of magnitude lower than chemically enriched stars by z=10. A catalog of more than 40,000 candidate Population III forming halos were identified, with formation times temporally ranging from z=30 to z=10, and ranging in mass from 2.3x10^5 M_sun to 1.2x10^10 M_sun. At early times, the environment that Population III stars form in is very similar to that of halos hosting chemically enriched star formation. At later times Population III stars are found to form in low-density regions that are not yet chemically polluted due to a lack of previous star formation in the area. Population III star forming halos become increasingly spatially isolated from one another at later times, and are generally closer to halos hosting chemically enriched star formation than to another halo hosting Population III star formation by z~10.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rhapsody. I. Structural Properties and Formation History From a Statistical Sample of Re-simulated Cluster-size Halos

We present the first results from the Rhapsody cluster re-simulation project: a sample of 96 "zoom-in" simulations of dark matter halos of 10^14.8 +- 0.05 Msun/h, selected from a 1 (Gpc/h)^3 volume. This simulation suite is the first to resolve this many halos with ~5x10^6 particles per halo in the cluster-mass regime, allowing us to statistically characterize the distribution of and correlation between halo properties at fixed mass. We focus on the properties of the main halos and how they are affected by formation history, which we track back to z=12, over five decades in mass. We give particular attention to the impact of the formation history on the density profiles of the halos. We find that the deviations from the Navarro-Frenk-White (NFW) model and the Einasto model depend on formation time. Late-forming halos tend to have considerable deviations from both models, partly due to the presence of massive subhalos, while early-forming halos deviate less but still significantly from the NFW model and are better described by the Einasto model. We find that the halo shapes depend only moderately on formation time. Departure from spherical symmetry impacts the density profiles through the anisotropic distribution of massive subhalos. Further evidence of the impact of subhalos is provided by analyzing the phase-space structure. A detailed analysis of the properties of the subhalo population in Rhapsody is presented in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rhapsody. II. Subhalo Properties and the Impact of Tidal Stripping From a Statistical Sample of Cluster-Size Halos

We discuss the properties of subhalos in cluster-size halos, using a high-resolution statistical sample: the Rhapsody simulations introduced in Wu et al. (2012). We demonstrate that the criteria applied to select subhalos have significant impact on the inferred properties of the sample, including the scatter in the number of subhalos, the correlation between the subhalo number and formation time, and the shape of subhalos' spatial distribution and velocity structure. We find that the number of subhalos, when selected using the peak maximum circular velocity in their histories (a property expected to be closely related to the galaxy luminosity), is uncorrelated with the formation time of the main halo. This is in contrast to the previously reported correlation from studies where subhalos are selected by the current maximum circular velocity; we show that this difference is a result of the tidal stripping of the subhalos. We also find that the dominance of the main halo and the subhalo mass fraction are strongly correlated with halo concentration and formation history. These correlations are important to take into account when interpreting results from cluster samples selected with different criteria. Our sample also includes a fossil cluster, which is presented separately and placed in the context of the rest of the sample.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The locations of halo formation and the peaks formalism

We investigate the problem of predicting the halo mass function from the properties of the Lagrangian density field. We focus on a perturbation spectrum with a small-scale cut-off (as in warm dark matter cosmologies). This cut-off results in a strong suppression of low mass objects, providing additional leverage to rigorously test which perturbations collapse and to what mass. We find that all haloes are consistent with forming near peaks of the initial density field, with a strong correlation between proto-halo density and ellipticity. We demonstrate that, while standard excursion set theory with correlated steps completely fails to reproduce the mass function, the inclusion of the peaks constraint leads to the correct number of haloes but significantly underpredicts the masses of low-mass objects (with the predicted halo mass function at low masses behaving like dn/dln m ~ m^{2/3}). This prediction is very robust and cannot be easily altered within the framework of a single collapse barrier. The nature of collapse in the presence of a small-scale cut-off thus reveals that excursion set calculations require a more detailed understanding of the collapse-time of a general ellipsoidal perturbation to predict the ultimate collapsed mass of a peak -- a problem that has been hidden in the large abundance of small-scale structure in CDM. We demonstrate how this problem can be resolved within the excursion set framework.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Warm DM halo mass function below the cut-off scale

Warm Dark Matter (WDM) cosmologies are a viable alternative to the Cold Dark Matter (CDM) scenario. Unfortunately, an accurate scrutiny of the WDM predictions with N-body simulations has proven difficult due to numerical artefacts. Here, we report on cosmological simulations that, for the first time, are devoid of those problems, and thus, are able to accurately resolve the WDM halo mass function well below the cut-off. We discover a complex picture, with perturbations at different evolutionary stages populating different ranges in the halo mass function. On the smallest mass scales we can resolve, identified objects are typically centres of filaments that are starting to collapse. On intermediate mass scales, objects typically correspond to fluctuations that have collapsed and are in the process of relaxation, whereas the high mass end is dominated by objects similar to haloes identified in CDM simulations. We then explicitly show how the formation of low-mass haloes is suppressed, which translates into a strong cut-off in the halo mass function. This disfavours some analytic formulations that predict a halo mass function that would extend well below the free streaming mass. We argue for a more detailed exploration of the formation of the smallest structures expected to form in a given cosmology, which, we foresee, will advance our overall understanding of structure formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Virial scaling of galaxies in clusters: bright to faint is cool to hot

By combining galaxy tracers from high-resolution N-body and hydrodynamical simulations, we present a consistent picture of the behaviour of galaxy velocities in massive clusters. In haloes above ~ 10^14 Msun, the brightest satellite galaxies are slightly cooler compared to the dark matter, while fainter satellites are hotter. Within the virial radius of a cluster, the mean velocity dispersion based on the 100 brightest galaxies is a factor of 1.065 +/- 0.005 (stat) +/- 0.027 (sys) higher than that of the dark matter (corresponding to a ~10-15 per cent bias in the dynamical mass estimate) while that based on only the five brightest galaxies is 0.868 +/- 0.039 (stat) +/- 0.035 (sys). These trends are approximately independent of redshift. The velocity structure is sensitive to the modelling of galaxies in clusters, indicative of the complex interplay of tidal stripping, dynamical friction, and merging. Velocity dispersions derived from instantaneous subhalo properties are larger than those employing either peak subhalo properties or hydrodynamical galaxy tracers. The latter two methods are consistent, implying that stacked spectroscopic analysis of cluster samples should, after correction for projection, show a trend towards slightly higher velocities when fainter galaxies are included, with an unbiased measure of dark matter velocity dispersion coming from approximately 30 galaxies per cluster. We show evidence that the velocity distribution function of bright galaxies near the cluster centre has a low-velocity tail due to strong dynamical friction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Novel Approach to Visualizing Dark Matter Simulations

In the last decades cosmological N-body dark matter simulations have enabled ab initio studies of the formation of structure in the Universe. Gravity amplified small density fluctuations generated shortly after the Big Bang, leading to the formation of galaxies in the cosmic web. These calculations have led to a growing demand for methods to analyze time-dependent particle based simulations. Rendering methods for such N-body simulation data usually employ some kind of splatting approach via point based rendering primitives and approximate the spatial distributions of physical quantities using kernel interpolation techniques, common in SPH (Smoothed Particle Hydrodynamics)-codes. This paper proposes three GPU-assisted rendering approaches, based on a new, more accurate method to compute the physical densities of dark matter simulation data. It uses full phase-space information to generate a tetrahedral tessellation of the computational domain, with mesh vertices defined by the simulation's dark matter particle positions. Over time the mesh is deformed by gravitational forces, causing the tetrahedral cells to warp and overlap. The new methods are well suited to visualize the cosmic web. In particular they preserve caustics, regions of high density that emerge, when several streams of dark matter particles share the same location in space, indicating the formation of structures like sheets, filaments and halos. We demonstrate the superior image quality of the new approaches in a comparison with three standard rendering techniques for N-body simulation data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Co-planar streams, pancakes, and angular-momentum exchange in high-z disc galaxies

We study the feeding of massive galaxies at high redshift through streams from the cosmic web using the Mare Nostrum hydro-cosmological simulation. Our statistical sample consists of 350 dark-matter haloes of ~10^12 Msun at z = 2.5. We find that ~70% of the influx into the virial radius Rv is in narrow streams covering 10% of the virial shell. On average 64% of the stream influx is in one stream, and 95% is in three dominant streams. The streams that feed a massive halo tend to lie in a plane that extends from half to a few Rv, hereafter "the stream plane" (SP). The streams are typically embedded in a thin sheet of low-entropy gas, a Zel'dovich pancake, which carries ~20% of the influx into Rv. The filaments-in-a-plane configuration about the massive haloes at the nodes of the cosmic web differs from the large- scale structure of the web where the filaments mark the intersections of slanted sheets. The stream plane is only weakly aligned with the angular momentum (AM) near Rv, consistent with the fact that typically 80% of the AM is carried by one dominant stream. The galactic disc plane shows a weak tendency to be perpendicular to the large-scale SP, consistent with tidal-torque theory. Most interesting, the direction of the disc AM is only weakly correlated with the AM direction at Rv. This indicates a significant AM exchange at the interphase between streams and disc in the greater environment of the disc inside an "AM sphere" of radius ~0.3Rv . The required large torques are expected based on the perturbed morphology and kinematics within this interaction sphere. This AM exchange may or may not require a major modification of the standard disc modeling based on AM conservation, depending on the extent to which the amplitude of the disc AM is affected, which is yet to be studied.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Halo-to-Halo Similarity and Scatter in the Velocity Distribution of Dark Matter

We examine the Velocity Distribution Function (VDF) in dark matter halos from Milky Way to cluster mass scales. We identify an empirical model for the VDF with a wider peak and a steeper tail than a Maxwell--Boltzmann distribution, and discuss physical explanations. We quantify sources of scatter in the VDF of cosmological halos and their implication for direct detection of dark matter. Given modern simulations and observations, we find that the most significant uncertainty in the VDF of the Milky Way arises from the unknown radial position of the solar system relative to the dark matter halo scale radius.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tracing the Dark Matter Sheet in Phase Space

The primordial velocity dispersion of dark matter is small compared to the velocities attained during structure formation. The initial density distribution is close to uniform and it occupies an initial sheet in phase space that is single valued in velocity space. Because of gravitational forces this three dimensional manifold evolves in phase space without ever tearing, conserving phase-space volume and preserving the connectivity of nearby points. N-body simulations already follow the motion of this sheet in phase space. This fact can be used to extract full fine-grained phase-space-structure information from existing cosmological N-body simulations. Particles are considered as the vertices of an unstructured three dimensional mesh, moving in six dimensional phase-space. On this mesh, mass density and momentum are uniquely defined. We show how to obtain the space density of the fluid, detect caustics, and count the number of streams as well as their individual contributions to any point in configuration-space. We calculate the bulk velocity, local velocity dispersions, and densities from the sheet - all without averaging over control volumes. This gives a wealth of new information about dark matter fluid flow which had previously been thought of as inaccessible to N-body simulations. We outline how this mapping may be used to create new accurate collisionless fluid simulation codes that may be able to overcome the sparse sampling and unphysical two-body effects that plague current N-body techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multi-scale initial conditions for cosmological simulations

We discuss a new algorithm to generate multi-scale initial conditions with multiple levels of refinements for cosmological "zoom-in" simulations. The method uses an adaptive convolution of Gaussian white noise with a real space transfer function kernel together with an adaptive multi-grid Poisson solver to generate displacements and velocities following first (1LPT) or second order Lagrangian perturbation theory (2LPT). The new algorithm achieves RMS relative errors of order 10^(-4) for displacements and velocities in the refinement region and thus improves in terms of errors by about two orders of magnitude over previous approaches. In addition, errors are localized at coarse-fine boundaries and do not suffer from Fourier-space induced interference ringing. An optional hybrid multi-grid and Fast Fourier Transform (FFT) based scheme is introduced which has identical Fourier space behaviour as traditional approaches. Using a suite of re-simulations of a galaxy cluster halo our real space based approach is found to reproduce correlation functions, density profiles, key halo properties and subhalo abundances with per cent level accuracy. Finally, we generalize our approach for two-component baryon and dark-matter simulations and demonstrate that the power spectrum evolution is in excellent agreement with linear perturbation theory. For initial baryon density fields, it is suggested to use the local Lagrangian approximation in order to generate a density field for mesh based codes that is consistent with Lagrangian perturbation theory instead of the current practice of using the Eulerian linearly scaled densities.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Large-Scale Orientations of Disk Galaxies

We use a 380 h-1 pc resolution hydrodynamic AMR simulation of a cosmic filament to investigate the orientations of a sample of ~100 well-resolved galactic disks spanning two orders of magnitude in both stellar and halo mass. We find: (i) At z=0, there is an almost perfect alignment at a median angle of 18 deg, in the inner dark matter halo regions where the disks reside, between the spin vector of the gaseous and stellar galactic disks and that of their inner host haloes. The alignment between galaxy spin and spin of the entire host halo is however significantly weaker, ranging from a median of ~46 deg at z=1 to ~50 deg at z=0. (ii) The most massive galaxy disks have spins preferentially aligned so as to point along their host filaments. (iii) The spin of disks in lower-mass haloes shows, at redshifts above z~0.5 and in regions of low environmental density, a clear signature of alignment with the intermediate principal axis of the large-scale tidal field. This behavior is consistent with predictions of linear tidal torque theory. This alignment decreases with increasing environmental density, and vanishes in the highest density regions. Non-linear effects in the high density environments are plausibly responsible for establishing this density-alignment correlation. We expect that our numerical results provide important insights for both understanding intrinsic alignment in weak lensing from the astrophysical perspective and formation and evolution processes of galactic disks in a cosmological context.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Halo mass function and scale-dependent bias from N-body simulations with non-Gaussian initial conditions

We perform a series of high-resolution N-body simulations of cosmological structure formation starting from Gaussian and non-Gaussian initial conditions. We adopt the best-fitting cosmological parameters of WMAP (3rd- and 5th-year) and we consider non-Gaussianity of the local type parameterised by 8 different values of the non-linearity parameter F_NL. Building upon previous work based on the Gaussian case, we show that, expressed in terms of suitable variables, the mass function of friends-of-friends haloes is approximately universal (i.e. independent of redshift, cosmology, and matter transfer function) to good precision (nearly 10 per cent) also in non-Gaussian scenarios. We provide fitting formulae for the high-mass end (M>10^13 M_sol/h) of the universal mass function in terms of F_NL, and we also present a non-universal fit in terms of both F_NL and z to be used for applications requiring higher accuracy. In the Gaussian case, we extend our fit to a wider range of halo masses (M>2.4 x 10^10 M_sol/h) and we also provide a consistent fit of the linear halo bias. We show that, for realistic values of F_NL, the matter power-spectrum in non-Gaussian cosmologies departs from the Gaussian one by up to 2 per cent on the scales where the baryonic- oscillation features are imprinted on the 2-point statistics. We confirm the strong k-dependence of the halo bias on large scales (k<0.05 h Mpc^-1) which was already detected in previous studies. However, we find that commonly used parameterisations based on the peak-background split do not provide an accurate description of our simulations which present extra dependencies on the wavenumber, the non-linearity parameter and, possibly, the clustering strength. We provide an accurate fit of the simulation data that can be used as a benchmark for future determinations of F_NL with galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lagrangian Statistics of Dark Halos in a LCDM Cosmology

New statistical properties of dark matter halos in Lagrangian space are presented. Tracing back the dark matter particles constituting bound halos resolved in a series of N-body simulations, we measure quantitatively the correlations of the proto-halo's inertia tensors with the local tidal tensors and investigate how the correlation strength depends on the proto-halo's sphericity, local density and filtering scale. It is shown that the majority of the proto-halos exhibit strong correlations between the two tensors provided that the tidal field is smoothed on the proto-halo's mass scale. The correlation strength is found to increase as the proto-halo's sphericity increases, as the proto-halo's mass increases, and as the local density becomes close to the critical value, delta_{ec}. It is also found that those peculiar proto-halos which exhibit exceptionally weak correlations between the two tensors tend to acquire higher specific angular momentum in Eulerian space, which is consistent with the linear tidal torque theory. In the light of our results, it is intriguing to speculate a hypothesis that the low surface brightness galaxies observed at present epoch correspond to the peculiar proto-halos with extreme low-sphericity whose inertia tensors are weakly correlated with the local tidal tensors.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Tidal Effects and the Environment Dependence of Halo Assembly

We explore a possible origin for the puzzling anti-correlation between the formation epoch of galactic dark-matter haloes and their environment density. This correlation has been revealed from cosmological N-body simulations and is in conflict with the Extended Press-Schechter model of halo clustering. Using similar simulations, we first quantify the straightforward association of an early formation epoch with a reduced mass growth rate at late times. We then find that a primary driver of suppressed growth, by accretion and mergers, is tidal effects dominated by a neighbouring massive halo. The tidal effects range from a slowdown of the assembly of haloes due to the shear along the large-scale filaments that feed the massive halo to actual mass loss in haloes that pass through the massive halo. Using the restricted three-body problem, we show that haloes are prone to tidal mass loss within 1.5 virial radii of a larger halo. Our results suggest that the dependence of formation epoch on environment density is a secondary effect induced by the enhanced density of haloes in filaments near massive haloes where the tides are strong. Our measures of assembly rate are particularly correlated with the tidal field at high redshifts z~1.

0 citations0 reviewsNo code linkedOpen access