Researcher profile

Pascal Elahi

Pascal Elahi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
12works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2021arXiv

Performance Enhancement of Tree-based Friends-of-friend Galaxy-finder for High-resolution Simulations of Galaxy Formation

Cosmological simulations are useful tools for studying the evolution of galaxies, and it is critical to accurately identify galaxies and their halos from raw simulation data. The friends-of-friend (FoF) algorithm has been widely adopted for this purpose because of its simplicity and expandability to higher dimensions. However, it is cost-inefficient when applied to high-resolution simulations because standard FoF implementation leads to too many distance calculations in dense regions. We confirm this through our exercise of applying the 6-dimensional (6D) FoF galaxy finder code, VELOCIraptor (Elahi et al. 2019), on the NewHorizon simulation (Dubois et al. 2021). The high particle resolution of NewHorizon ($M_{\rm star} \sim 10^4 M_{\odot}$) allows a large central number density ($10^{6}\,{\rm kpc}^{-3}$) for typical galaxies, resulting in a few days to weeks of galaxy searches for just one snapshot. Even worse, we observed a significant decrease in the FoF performance in the high-dimensional 6D searches: "the curse of dimensionality" problem. To overcome these issues, we have developed several implementations that can be readily applied to any tree-based FoF code. They include limiting visits to tree nodes, re-ordering the list of particles for searching neighbor particles, and altering the tree structure. Compared to the run with the original code, the new run with these implementations results in the identical galaxy detection with the ideal performance, $O(N \log{N})$, $N$ being the number of particles in a galaxy -- with a speed gain of a factor of 2700 in 3D or 12 in 6D FoF search.

preprint2015arXiv

Major Mergers Going Notts: Challenges for Modern Halo Finders

Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z>1). This raises significant uncertainties for theoretical models that require major mergers for, e.g., galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.

preprint2014arXiv

A thousand shadows of Andromeda: rotating planes of satellites in the Millennium-II cosmological simulation

In a recent contribution, Bahl \& Baumgardt investigated the incidence of planar alignments of satellite galaxies in the Millennium-II simulation, and concluded that vast thin planes of dwarf galaxies, similar to that observed in the Andromeda galaxy (M31), occur frequently by chance in $Λ$-Cold Dark Matter cosmology. However, their analysis did not capture the essential fact that the observed alignment is simultaneously radially extended, yet thin, and kinematically unusual. With the caveat that the Millennium-II simulation may not have sufficient mass resolution to identify confidently simulacra of low-luminosity dwarf galaxies, we re-examine that simulation for planar structures, using the same method as employed by Ibata et al. (2013) on the real M31 satellites. We find that 0.04\% of host galaxies display satellite alignments that are at least as extreme as the observations, when we consider their extent, thickness and number of members rotating in the same sense. We further investigate the angular momentum properties of the co-planar satellites, and find that the median of the specific angular momentum derived from the line of sight velocities in the real M31 structure ($1.3\times10^4$ km/s kpc) is very high compared to systems drawn from the simulations. This analysis confirms that it is highly unlikely that the observed structure around the Andromeda galaxy is due to a chance occurrence. Interestingly, the few extreme systems that are similar to M31 arise from the accretion of a massive sub-halo with its own spatially-concentrated entourage of orphan satellites.

preprint2014arXiv

Accretion in action: phase space coherence of stellar debris and globular clusters in Andromeda's South-West Cloud

A central tenet of the current cosmological paradigm is that galaxies grow over time through the accretion of smaller systems. Here, we present new kinematic measurements near the centre of one of the densest pronounced substructures, the South-West Cloud, in the outer halo of our nearest giant neighbour, the Andromeda galaxy. These observations reveal that the kinematic properties of this region of the South-West Cloud are consistent with those of PA-8, a globular cluster previously shown to be co-spatial with the stellar substructure. In this sense the situation is reminiscent of the handful of globular clusters that sit near the heart of the Sagittarius dwarf galaxy, a system that is currently being accreted into the Milky Way, confirming that accretion deposits not only stars but also globular clusters into the halos of large galaxies.

preprint2014arXiv

Gravitational lensing in WDM cosmologies: The cross section for giant arcs

The nature of the dark sector of the Universe remains one of the outstanding problems in modern cosmology, with the search for new observational probes guiding the development of the next generation of observational facilities. Clues come from tension between the predictions from ΛCDM and observations of gravitationally lensed galaxies. Previous studies showed that galaxy clusters in the ΛCDM are not strong enough to reproduce the observed number of lensed arcs. This work aims to constrain the warm dark matter cosmologies by means of the lensing efficiency of galaxy clusters drawn from these alternative models. The lensing characteristics of two samples of simulated clusters in the warm dark matter (ΛWDM) and cold dark matter (ΛCDM) cosmologies have been studied. The results show that even though the CDM clusters are more centrally concentrated and contain more substructures, the WDM clusters have slightly higher lensing efficiency than their CDM counterparts. The key difference is that WDM clusters have more extended and more massive subhaloes than CDM analogues. These massive substructures significantly stretch the critical lines and caustics and hence they boost the lensing efficiency of the host halo. Despite the increase in the lensing efficiency due to the contribution of massive substructures in the WDM clusters, this is not enough to resolve the arc statistics problem.

preprint2014arXiv

Subhaloes gone Notts: Subhaloes as tracers of the dark matter halo shape

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

preprint2014arXiv

Subhaloes gone Notts: the clustering properties of subhaloes

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

preprint2014arXiv

The life and death of cosmic voids

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

preprint2013arXiv

Galaxies going MAD: The Galaxy-Finder Comparison Project

With the ever increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations has simultaneously grown. To this extent we initiated the Halo Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes (Knebe et al. 2011), and one for dark matter subhaloes (Onions et al. 2012). However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also account for baryonic physics (e.g. hydrodynamics, star formation, feedback) object finders should also be capable of taking these additional processes into consideration. Here we report on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the MW, M31, and M33) as well as their satellite populations for a variety of halo finders ranging from phase-space to velocity-space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons, at least for the total, dark, and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.

preprint2013arXiv

Structure Finding in Cosmological Simulations: The State of Affairs

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

preprint2012arXiv

Subhaloes gone Notts: Spin across subhaloes and finders

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