-
Examining the Effects of Dark Matter Spikes on Eccentric Intermediate Mass Ratio Inspirals Using $N$-body Simulations
Authors:
Diptajyoti Mukherjee,
A. Miguel Holgado,
Go Ogiya,
Hy Trac
Abstract:
Recent studies have postulated that the presence of dark matter (DM) spikes around IMBHs could lead to observable dephasing effects in gravitational wave (GW) signals emitted by Intermediate Mass Ratio Inspirals (IMRIs). While prior investigations primarily relied on non-self-consistent analytic methods to estimate the influence of DM spikes on eccentric IMRIs, our work introduces the first self-c…
▽ More
Recent studies have postulated that the presence of dark matter (DM) spikes around IMBHs could lead to observable dephasing effects in gravitational wave (GW) signals emitted by Intermediate Mass Ratio Inspirals (IMRIs). While prior investigations primarily relied on non-self-consistent analytic methods to estimate the influence of DM spikes on eccentric IMRIs, our work introduces the first self-consistent treatment of this phenomenon through $N$-body simulations. Contrary to previous studies, which suggested that dynamical friction (DF), a cumulative effect of two-body encounters, is the primary mechanism responsible for energy dissipation, we reveal that the slingshot mechanism, a three-body effect, plays a more significant role in driving the binary system's energy loss and consequent orbital shrinkage, similar to stellar loss cone scattering in Massive Black Hole (MBH) binaries. Furthermore, our work extends its analysis to include rotation in DM spikes, a factor often overlooked in previous studies. We find that binaries that counter-rotate with respect to the spike particles merge faster, while binaries that co-rotate merge slower, in opposition to the expectation from DF theory. While our models are idealistic, they offer findings that pave the way for a more comprehensive understanding of the complex interactions between DM spikes, IMRIs, GW emission, and the ability to constrain DM microphysics. Our work systematically includes Post-Newtonian (PN) effects until 2.5 order and our results are accurate and robust.
△ Less
Submitted 6 February, 2024; v1 submitted 4 December, 2023;
originally announced December 2023.
-
Formation of dense filaments induced by runaway supermassive black holes
Authors:
Go Ogiya,
Daisuke Nagai
Abstract:
A narrow linear object extending $\sim 60 \,{\rm kpc}$ from the centre of a galaxy at redshift $z \sim 1$ has recently been discovered and interpreted as shocked gas filament forming stars. The host galaxy presents an irregular morphology, implying recent merger events. Supposing that each of the progenitor galaxies has a central supermassive black hole (SMBH) and the SMBHs are accumulated at the…
▽ More
A narrow linear object extending $\sim 60 \,{\rm kpc}$ from the centre of a galaxy at redshift $z \sim 1$ has recently been discovered and interpreted as shocked gas filament forming stars. The host galaxy presents an irregular morphology, implying recent merger events. Supposing that each of the progenitor galaxies has a central supermassive black hole (SMBH) and the SMBHs are accumulated at the centre of the merger remnant, a fraction of them can be ejected from the galaxy with a high velocity due to interactions between SMBHs. When such a runaway SMBH (RSMBH) passes through the circumgalactic medium (CGM), converging flows are induced along the RSMBH path, and star formation could eventually be ignited. We show that the CGM temperature prior to the RSMBH perturbation should be below the peak temperature in the cooling function to trigger filament formation. While the gas is temporarily heated due to compression, the cooling efficiency increases, and gas accumulation becomes allowed along the path. When the CGM density is sufficiently high, the gas can cool down and develop a dense filament by $z = 1$. The mass and velocity of the RSMBH determine the scale of filament formation. Hydrodynamical simulations validate the analytical expectations. Therefore, we conclude that the perturbation by RSMBHs is a viable channel to form the observed linear object. Using the analytic model validated by simulations, we show that the CGM around the linear object to be warm ($T < 2 \times 10^5 \, K$) and dense ($n > 2 \times 10^{-5} (T/2 \times 10^5 \, K)^{-1} \, {\rm cm^{-3}}$).
△ Less
Submitted 17 November, 2023; v1 submitted 16 September, 2023;
originally announced September 2023.
-
Searching for clues of past binary supermassive black hole mergers in nuclear star clusters
Authors:
Alessandra Mastrobuono-Battisti,
Go Ogiya,
Oliver Hahn,
Mathias Schultheis
Abstract:
Galaxy mergers are common processes in the Universe. As a large fraction of galaxies hosts at their centres a central supermassive black hole (SMBH), mergers can lead to the formation of a supermassive black hole binary (SMBHB). The formation of such a binary is more efficient when the SMBHs are embedded in a nuclear star cluster (NSC). NSCs are dense and massive stellar clusters present in the ma…
▽ More
Galaxy mergers are common processes in the Universe. As a large fraction of galaxies hosts at their centres a central supermassive black hole (SMBH), mergers can lead to the formation of a supermassive black hole binary (SMBHB). The formation of such a binary is more efficient when the SMBHs are embedded in a nuclear star cluster (NSC). NSCs are dense and massive stellar clusters present in the majority of the observed galaxies. Their central densities can reach up to $10^7\,M_\odot/{\rm pc}^3$ and their masses can be as large as a few $10^7\,M_\odot$. The direct detection of an SMBHB is observationally challenging. In this work, we illustrate how the large scale structural and dynamical properties of an NSC can help to identify nucleated galaxies that recently went through a merger that possibly led to the formation of a central SMBHB. Our models show that the merger can imprint signatures on the shape, density profile, rotation and velocity structure of the NSC. The strength of the signatures depends on the mass ratio between the SMBHs and on the orbital initial conditions of the merger. In addition, the number of hypervelocity stars produced in the mergers is linked to the SMBHB properties. The merger can also contribute to the formation of the nuclear stellar disc of the galaxy.
△ Less
Submitted 22 March, 2023;
originally announced March 2023.
-
The effect of stellar encounters on the dark matter annihilation signal from prompt cusps
Authors:
Jens Stücker,
Go Ogiya,
Simon D. M. White,
Raul E. Angulo
Abstract:
Prompt cusps are the densest quasi-equilibrium dark matter objects; one forms at the instant of collapse within every isolated peak of the initial cosmological density field. They have power-law density profiles, $ρ\propto r^{-1.5}$ with central phase-space density set by the primordial velocity dispersion of the dark matter. At late times they account for $\sim 1\%$ of the dark matter mass but fo…
▽ More
Prompt cusps are the densest quasi-equilibrium dark matter objects; one forms at the instant of collapse within every isolated peak of the initial cosmological density field. They have power-law density profiles, $ρ\propto r^{-1.5}$ with central phase-space density set by the primordial velocity dispersion of the dark matter. At late times they account for $\sim 1\%$ of the dark matter mass but for $>90\%$ of its annihilation luminosity in all but the densest regions, where they are tidally disrupted. Here we demonstrate that individual stellar encounters, rather than the mean galactic tide, are the dominant disruptors of prompt cusps within galaxies. Their cumulative effect is fully (though stochastically) characterised by an impulsive shock strength $B_* = 2πG\intρ_*({\bf x}(t))\, \mathrm{d}t$ where $ρ_*$, the total mass density in stars, is integrated over a cusp's entire post-formation trajectory. Stellar encounters and mean tides have only a small effect on the halo annihilation luminosity seen by distant observers, but this is not true for the Galactic halo because of the Sun's position. For a 100 GeV WIMP, Earth-mass prompt cusps are predicted, and stellar encounters suppress their mean annihilation luminosity by a factor of two already at 20 kpc, so that their annihilation emission is predicted to appear almost uniform over the sky. The Galactic Center $γ$-ray Excess is thus unaffected by cusps. If it is indeed dark matter annihilation radiation, then prompt cusps in the outer Galactic halo and beyond must account for 20-80% of the observed isotropic $γ$-ray background in the 1 to 10 GeV range.
△ Less
Submitted 25 April, 2023; v1 submitted 11 January, 2023;
originally announced January 2023.
-
Testing the galaxy collision induced formation scenario for the trail of dark matter deficient galaxies with the susceptibility of globular clusters to the tidal force
Authors:
Go Ogiya,
Frank C. van den Bosch,
Andreas Burkert,
Xi Kang
Abstract:
It has been suggested that a trail of diffuse galaxies, including two dark matter deficient galaxies (DMDGs), in the vicinity of NGC1052 formed because of a high-speed collision between two gas-rich dwarf galaxies, one bound to NGC1052 and the other one on an unbound orbit. The collision compresses the gas reservoirs of the colliding galaxies, which in turn triggers a burst of star formation. In c…
▽ More
It has been suggested that a trail of diffuse galaxies, including two dark matter deficient galaxies (DMDGs), in the vicinity of NGC1052 formed because of a high-speed collision between two gas-rich dwarf galaxies, one bound to NGC1052 and the other one on an unbound orbit. The collision compresses the gas reservoirs of the colliding galaxies, which in turn triggers a burst of star formation. In contrast, the dark matter and pre-existing stars in the progenitor galaxies pass through it. Since the high pressures in the compressed gas are conducive to the formation of massive globular clusters (GCs), this scenario can explain the formation of DMDGs with large populations of massive GCs, consistent with the observations of NGC1052-DF2 (DF2) and NGC1052-DF4. A potential difficulty with this `mini bullet cluster' scenario is that the observed spatial distributions of GCs in DMDGs are extended. GCs experience dynamical friction causing their orbits to decay with time. Consequently, their distribution at formation should have been even more extended than that observed at present. Using a semi-analytic model, we show that the observed positions and velocities of the GCs in DF2 imply that they must have formed at a radial distance of 5-10kpc from the center of DF2. However, as we demonstrate, the scenario is difficult to reconcile with the fact that the strong tidal forces from NGC1052 strip the extendedly distributed GCs from DF2, requiring 33-59 massive GCs to form at the collision to explain observations.
△ Less
Submitted 10 November, 2022;
originally announced November 2022.
-
Shedding light on low-mass subhalo survival and annihilation luminosity with numerical simulations
Authors:
Alejandra Aguirre-Santaella,
Miguel A. Sánchez-Conde,
Go Ogiya,
Jens Stücker,
Raul E. Angulo
Abstract:
In this work, we carry out a suite of specially-designed numerical simulations to shed further light on dark matter (DM) subhalo survival at mass scales relevant for gamma-ray DM searches, a topic subject to intense debate nowadays. Specifically, we have developed and employed an improved version of DASH, a GPU $N$-body code, to study the evolution of low-mass subhaloes inside a Milky Way-like hal…
▽ More
In this work, we carry out a suite of specially-designed numerical simulations to shed further light on dark matter (DM) subhalo survival at mass scales relevant for gamma-ray DM searches, a topic subject to intense debate nowadays. Specifically, we have developed and employed an improved version of DASH, a GPU $N$-body code, to study the evolution of low-mass subhaloes inside a Milky Way-like halo with unprecedented accuracy, reaching solar-mass and sub-parsec resolution in our simulations. We simulate subhaloes with varying mass, concentration, and orbital properties, and consider the effect of the gravitational potential of the Milky Way galaxy itself. More specifically, we analyze the evolution of both the bound mass fraction and annihilation luminosity of subhaloes, finding that most subhaloes survive until present time, even though in some cases they lose more than 99% of their mass at accretion. Baryons in the host induce a much more severe mass loss, especially when the subhalo orbit is more parallel to the galactic disk. Many of these subhaloes cross the solar galactocentric radius, thus making it easier to detect their annihilation fluxes from Earth. We find subhaloes orbiting a DM-only halo with a pericentre in the solar vicinity to lose 70-90% of their initial annihilation luminosity at redshift zero, which increases up to 99% when baryons are also included in the host. We find a strong relation between subhalo's mass loss and the effective tidal field at pericentre. Indeed, much of the dependence on concentration, orbital parameters, host potential and baryonic components can be explained through this single parameter. In addition to shedding light on the survival of low-mass galactic subhaloes, our results can provide detailed predictions that will aid current and future quests for the nature of DM.
△ Less
Submitted 18 July, 2022;
originally announced July 2022.
-
Tidal Stripping in the Adiabatic Limit
Authors:
Jens Stücker,
Go Ogiya,
Raul E. Angulo,
Alejandra Aguirre-Santaella,
Miguel A. Sánchez-Conde
Abstract:
We present a model for the remnants of haloes that have gone through an adiabatic tidal stripping process. We show that this model exactly reproduces the remnant of an NFW halo that is exposed to a slowly increasing isotropic tidal field and approximately for an anisotropic tidal field. The model can be used to predict the asymptotic mass loss limit for orbiting subhaloes, solely as a function of…
▽ More
We present a model for the remnants of haloes that have gone through an adiabatic tidal stripping process. We show that this model exactly reproduces the remnant of an NFW halo that is exposed to a slowly increasing isotropic tidal field and approximately for an anisotropic tidal field. The model can be used to predict the asymptotic mass loss limit for orbiting subhaloes, solely as a function of the initial structure of the subhalo and the value of the tidal field at pericentre. Predictions can easily be made for differently concentrated host-haloes with and without baryonic components, which differ most notably in their relation between pericentre radius and tidal field. The model correctly predicts several empirically measured relations such as the `tidal track' and the `orbital frequency relation' that was reported by Errani & Navarro (2021) for the case of an isothermal sphere. Further, we propose applications of the `structure-tide' degeneracy which implies that increasing the concentration of a subhalo has exactly the same impact on tidal stripping as reducing the amplitude of the tidal field. Beyond this, we find that simple relations hold for the bound mass, truncation radius, WIMP annihilation luminosity and tidal ratio of tidally stripped NFW haloes in relation to quantities measured at the radius of maximum circular velocity. Finally, we note that NFW haloes cannot be completely disrupted when exposed adiabatically to tidal fields of arbitrary magnitudes. We provide an open-source implementation of our model and suggest that it can be used to improve predictions of dark matter annihilation.
△ Less
Submitted 17 March, 2023; v1 submitted 1 July, 2022;
originally announced July 2022.
-
Evolution of Massive Black Hole Binaries in Collisionally Relaxed Nuclear Star Clusters -- Impact of Mass Segregation
Authors:
Diptajyoti Mukherjee,
Qirong Zhu,
Go Ogiya,
Carl L. Rodriguez,
Hy Trac
Abstract:
Massive Black Hole (MBH) binaries are considered to be one of the most important sources of Gravitational Waves (GW) that can be detected by GW detectors like LISA. However, there are a lot of uncertainties in the dynamics of MBH binaries in the stages leading up to the GW-emission phase. It has been recently suggested that Nuclear Star Clusters (NSCs) could provide a viable route to overcome the…
▽ More
Massive Black Hole (MBH) binaries are considered to be one of the most important sources of Gravitational Waves (GW) that can be detected by GW detectors like LISA. However, there are a lot of uncertainties in the dynamics of MBH binaries in the stages leading up to the GW-emission phase. It has been recently suggested that Nuclear Star Clusters (NSCs) could provide a viable route to overcome the final parsec problem for MBH binaries at the center of galaxies. NSCs are collisional systems where the dynamics would be altered by the presence of a mass spectrum. In this study, we use a suite of high-resolution $N$-body simulations with over 1 million particles to understand how collisional relaxation under the presence of a mass spectrum of NSC particles affects the dynamics of the MBH binary under the merger of two NSCs. We consider MBH binaries with different mass ratios and additional non-relaxed models. We find that mass-segregation driven by collisional relaxation can lead to accelerated hardening in lower mass ratio binaries but has the opposite effect in higher mass ratio binaries. Crucially, the relaxed models also demonstrate much lower eccentricities at binary formation and negligible growth during hardening stages leading to longer merger timescales. The results are robust and highlight the importance of collisional relaxation on changing the dynamics of the binary. Our models are state-of-the-art, use zero softening, and high enough particle numbers to model NSCs realistically.
△ Less
Submitted 16 November, 2022; v1 submitted 24 May, 2022;
originally announced May 2022.
-
Dark matter cores in massive high-$z$ galaxies formed by baryonic clumps
Authors:
Go Ogiya,
Daisuke Nagai
Abstract:
The rotation curves of some star forming massive galaxies at redshift two decline over the radial range of a few times the effective radius, indicating a significant deficit of dark matter (DM) mass in the galaxy centre. The DM mass deficit is interpreted as the existence of a DM density core rather than the cuspy structure predicted by the standard cosmological model. A recent study proposed that…
▽ More
The rotation curves of some star forming massive galaxies at redshift two decline over the radial range of a few times the effective radius, indicating a significant deficit of dark matter (DM) mass in the galaxy centre. The DM mass deficit is interpreted as the existence of a DM density core rather than the cuspy structure predicted by the standard cosmological model. A recent study proposed that a galaxy merger, in which the smaller satellite galaxy is significantly compacted by dissipative contraction of the galactic gas, can heat the centre of the host galaxy and help make a large DM core. By using an $N$-body simulation, we find that a large amount of DM mass is imported to the centre by the merging satellite, making this scenario an unlikely solution for the DM mass deficit. In this work, we consider giant baryonic clumps in high redshift galaxies as alternative heating source for creating the baryon dominated galaxies with a DM core. Due to dynamical friction, the orbit of clumps decays in a few Gyr and the baryons condensate at the galactic centre. As a back-reaction, the halo centre is heated up and the density cusp is flattened out. The combination of the baryon condensation and core formation makes the galaxy baryon dominated in the central 2-5 kpc, comparable to the effective radius of the observed galaxies. Thus, the dynamical heating by giant baryonic clumps is a viable mechanism for explaining the observed dearth of DM in high redshift galaxies.
△ Less
Submitted 9 May, 2022;
originally announced May 2022.
-
On the Tidal Formation of Dark Matter Deficient Galaxies
Authors:
Go Ogiya,
Frank C. van den Bosch,
Andreas Burkert
Abstract:
Previous studies have shown that dark matter deficient galaxies (DMDG) such as NGC1052-DF2 (hereafter DF2) can result from tidal stripping. An important question, though, is whether such a stripping scenario can explain DF2's large specific frequency of globular clusters (GCs). After all, tidal stripping and shocking preferentially remove matter from the outskirts. We examine this using idealized,…
▽ More
Previous studies have shown that dark matter deficient galaxies (DMDG) such as NGC1052-DF2 (hereafter DF2) can result from tidal stripping. An important question, though, is whether such a stripping scenario can explain DF2's large specific frequency of globular clusters (GCs). After all, tidal stripping and shocking preferentially remove matter from the outskirts. We examine this using idealized, high-resolution simulations of a regular dark matter dominated galaxy that is accreted onto a massive halo. As long as the initial (pre-infall) dark matter halo of the satellite is cored, which is consistent with predictions of cosmological, hydrodynamical simulations, the tidal remnant can be made to resemble DF2 in all its properties, including its GC population. The required orbit has a peri-centre at the 8.3 percentile of the distribution for subhaloes at infall, and thus is not particularly extreme. On this orbit the satellite loses 98.5 (30) percent of its original dark matter (stellar) mass, and thus evolves into a DMDG. The fraction of GCs that is stripped off depends on the initial radial distribution. If, at infall, the median projected radius of the GC population is roughly two times that of the stars, consistent with observations of isolated galaxies, only $\sim 20$ percent of the GCs are stripped off. This is less than for the stars, which is due to dynamical friction counteracting the tidal stirring. We predict that, if indeed DF2 was crafted by strong tides, its stellar outskirts should have a very shallow metallicity gradient.
△ Less
Submitted 10 December, 2021; v1 submitted 23 November, 2021;
originally announced November 2021.
-
Evolution of subhalo orbits in a smoothly-growing host halo potential
Authors:
Go Ogiya,
James E. Taylor,
Michael J. Hudson
Abstract:
The orbital parameters of dark matter (DM) subhaloes play an essential role in determining their mass-loss rates and overall spatial distribution within a host halo. Haloes in cosmological simulations grow by a combination of relatively smooth accretion and more violent mergers, and both processes will modify subhalo orbits. To isolate the impact of the smooth growth of the host halo from other re…
▽ More
The orbital parameters of dark matter (DM) subhaloes play an essential role in determining their mass-loss rates and overall spatial distribution within a host halo. Haloes in cosmological simulations grow by a combination of relatively smooth accretion and more violent mergers, and both processes will modify subhalo orbits. To isolate the impact of the smooth growth of the host halo from other relevant mechanisms, we study subhalo orbital evolution using numerical calculations in which subhaloes are modelled as massless particles orbiting in a time-varying spherical potential. We find that the radial action of the subhalo orbit decreases over the first few orbits, indicating that the response to the growth of the host halo is not adiabatic during this phase. The subhalo orbits can shrink by a factor of $\sim$1.5 in this phase. Subsequently, the radial action is well conserved and orbital contraction slows down. We propose a model accurately describing the orbital evolution. Given these results, we consider the spatial distribution of the population of subhaloes identified in high-resolution cosmological simulations. We find that it is consistent with this population having been accreted at z < 3, indicating that any subhaloes accreted earlier are unresolved in the simulations. We also discuss tidal stripping as a formation scenario for NGC1052-DF2, an ultra diffuse galaxy significantly lacking DM, and find that its expected DM mass could be consistent with observational constraints if its progenitor was accreted early enough, z > 1.5, although it should still be a relatively rare object.
△ Less
Submitted 4 February, 2021;
originally announced February 2021.
-
Dynamical self-friction: how mass loss slows you down
Authors:
Tim B. Miller,
Frank C. van den Bosch,
Sheridan B. Green,
Go Ogiya
Abstract:
We investigate dynamical self-friction, the process by which material that is stripped from a subhalo torques its remaining bound remnant, which causes it to lose orbital angular momentum. By running idealized simulations of a subhalo orbiting within an analytical host halo potential, we isolate the effect of self-friction from traditional dynamical friction due to the host halo. While at some poi…
▽ More
We investigate dynamical self-friction, the process by which material that is stripped from a subhalo torques its remaining bound remnant, which causes it to lose orbital angular momentum. By running idealized simulations of a subhalo orbiting within an analytical host halo potential, we isolate the effect of self-friction from traditional dynamical friction due to the host halo. While at some points in a subhalo's orbit the torque of the stripped material can boost the orbital angular momentum of the remnant, the net effect over the long term is orbital decay regardless of the initial orbital parameters or subhalo mass. In order to quantify the strength of self-friction, we run a suite of simulations spanning typical host-to-subhalo mass ratios and orbital parameters. We find that the time-scale for self-friction, defined as the exponential decay time of the subhalo's orbital angular momentum, scales with mass ratio and orbital circularity similar to standard dynamical friction. The decay time due to self-friction is roughly an order of magnitude longer, suggesting that self-friction only contributes at the 10 percent level. However, along more radial orbits, self-friction can occasionally dominate over dynamical friction close to pericentric passage, where mass stripping is intense. This is also the epoch at which the self-friction torque undergoes large and rapid changes in both magnitude and direction, indicating that self-friction is an important process to consider when modeling pericentric passages of subhaloes and their associated satellite galaxies.
△ Less
Submitted 17 January, 2020;
originally announced January 2020.
-
Accelerated orbital decay of supermassive black hole binaries in merging nuclear star clusters
Authors:
Go Ogiya,
Oliver Hahn,
Chiara M. F. Mingarelli,
Marta Volonteri
Abstract:
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 sta…
▽ More
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.
△ Less
Submitted 11 February, 2020; v1 submitted 26 November, 2019;
originally announced November 2019.
-
DASH: a library of dynamical subhalo evolution
Authors:
Go Ogiya,
Frank C. van den Bosch,
Oliver Hahn,
Sheridan B. Green,
Tim B. Miller,
Andreas Burkert
Abstract:
The abundance and demographics of dark matter substructure is important for many areas in astrophysics and cosmological $N$-body simulations have been the primary tool used to investigate them. However, it has recently become clear that the simulations are subject to numerical artefacts, which hampers a proper treatment of the tidal evolution of subhaloes. Unfortunately, no analytical models that…
▽ More
The abundance and demographics of dark matter substructure is important for many areas in astrophysics and cosmological $N$-body simulations have been the primary tool used to investigate them. However, it has recently become clear that the simulations are subject to numerical artefacts, which hampers a proper treatment of the tidal evolution of subhaloes. Unfortunately, no analytical models that accurately describe subhalo evolution exist either. We therefore present a library of idealized, high resolution $N$-body simulations of the tidal evolution of individual subhaloes that can be used to calibrate semi-analytical models and to complement cosmological simulations. The simulations focus on minor mergers, i.e., the mass of the subhalo is much smaller than that of the host halo, such that the impact of dynamical friction is negligible. This setup allows the adoption of a fixed analytical potential for modelling the host halo. The dynamical evolution of subhaloes is followed with $N$-body computations. In the library, four parameters, two of which characterize the subhalo orbit with respect to the host halo, and the two concentrations of the host- and subhalo, are varied over the ranges encountered in cosmological simulations. We show several representative examples from the library that illustrate the evolution of the subhalo mass and velocity dispersion profiles. Additionally, we make publicly available a pre-trained non-parametric model of the subhalo mass evolution based on random forest regression. This model is able to interpolate the simulation data at the 0.1\,dex level and provides efficient access to the data for further use in modelling.
△ Less
Submitted 24 January, 2019;
originally announced January 2019.
-
Tidal stripping as a possible origin of the ultra diffuse galaxy lacking dark matter
Authors:
Go Ogiya
Abstract:
Recent observations revealed a mysterious ultra diffuse galaxy, NGC1052-DF2, in the group of a large elliptical galaxy, NGC1052. Compared to expectations from abundance matching models, the dark matter mass contained in NGC1052-DF2 is smaller by a factor of $\sim 400$. We utilize controlled $N$-body simulations of the tidal interaction between NGC1052 and a smaller satellite galaxy, that we suppos…
▽ More
Recent observations revealed a mysterious ultra diffuse galaxy, NGC1052-DF2, in the group of a large elliptical galaxy, NGC1052. Compared to expectations from abundance matching models, the dark matter mass contained in NGC1052-DF2 is smaller by a factor of $\sim 400$. We utilize controlled $N$-body simulations of the tidal interaction between NGC1052 and a smaller satellite galaxy, that we suppose as the progenitor of NGC1052-DF2, to test if tidal stripping can explain dark-matter deficiency at such levels. We find that when assuming a tightly bound and quite radial orbit and cored density structure for the dark halo of the satellite, our simulations reproduce well both the mass profile and the effective radius inferred from the observations of NGC1052-DF2. Orbital parameters are in the tail, but still consistent with measurements of their distributions from cosmological simulations. Such strongly dark-matter deficient galaxies, in our scenario, are thus expected to be relatively rare in groups and clusters, and not present in the field.
△ Less
Submitted 24 July, 2018; v1 submitted 17 April, 2018;
originally announced April 2018.
-
Physical and numerical stability and instability of AGN bubbles in a hot intracluster medium
Authors:
Go Ogiya,
Pawel Biernacki,
Oliver Hahn,
Romain Teyssier
Abstract:
While feedback from Active Galactic Nuclei (AGN) is an important heating source in the centre of galaxy clusters, it is still unclear how the feedback energy is injected into the intracluster medium (ICM) and what role different numerical approaches play. Here, we compare four hydrodynamical schemes in idealized simulations of a rising bubble inflated by AGN feedback in a hot stratified ICM: (trad…
▽ More
While feedback from Active Galactic Nuclei (AGN) is an important heating source in the centre of galaxy clusters, it is still unclear how the feedback energy is injected into the intracluster medium (ICM) and what role different numerical approaches play. Here, we compare four hydrodynamical schemes in idealized simulations of a rising bubble inflated by AGN feedback in a hot stratified ICM: (traditional) smoothed particle hydrodynamics (TSPH), a pressure flavour of SPH (PSPH), a meshless finite mass (MFM) scheme, as well as an Eulerian code with adaptive mesh refinement. In the absence of magnetic fields, the bubble is Kelvin-Helmholtz unstable on short enough time scales to dissolve it fully in the ICM, which is captured by MFM and RAMSES simulations, while in the TSPH simulation the bubble survives. When the ICM is turbulent, mixing of the bubble with the ICM is accelerated. This occurs if the numerical scheme can capture the instabilities well. The differences in the evolution of the bubble has a surprisingly small influence on the thermal structure of the ICM. However, in the simulations with MFM and RAMSES the bubble disruption leads to turbulent stirring of the ICM which is suppressed in SPH. In the latter the thermal energy remains trapped in the bubble and is transported to large radii. We discuss if the choice of hydrodynamical schemes can lead to systematic differences in the outcomes of cosmological simulations.
△ Less
Submitted 6 February, 2018;
originally announced February 2018.
-
Dark Matter Substructure in Numerical Simulations: A Tale of Discreteness Noise, Runaway Instabilities, and Artificial Disruption
Authors:
Frank C. van den Bosch,
Go Ogiya
Abstract:
To gain understanding of the complicated, non-linear and numerical processes associated with the tidal evolution of dark matter subhaloes in numerical simulation, we perform a large suite of idealized simulations that follow individual N-body subhaloes in a fixed, analytical host halo potential. By varying both physical and numerical parameters, we investigate under what conditions the subhaloes u…
▽ More
To gain understanding of the complicated, non-linear and numerical processes associated with the tidal evolution of dark matter subhaloes in numerical simulation, we perform a large suite of idealized simulations that follow individual N-body subhaloes in a fixed, analytical host halo potential. By varying both physical and numerical parameters, we investigate under what conditions the subhaloes undergo disruption. We confirm the conclusions from our more analytical assessment in van den Bosch et al. that most disruption is numerical in origin; as long as a subhalo is resolved with sufficient mass and force resolution, a bound remnant survives. This implies that state-of-the-art cosmological simulations still suffer from significant overmerging. We demonstrate that this is mainly due to inadequate force-softening, which causes excessive mass loss and artificial tidal disruption. In addition, we show that subhaloes in N-body simulations are susceptible to a runaway instability triggered by the amplification of discreteness noise in the presence of a tidal field. These two processes conspire to put serious limitations on the reliability of dark matter substructure in state-of-the-art cosmological simulations. We present two criteria that can be used to assess whether individual subhaloes in cosmological simulations are reliable or not, and advocate that subhaloes that satisfy either of these two criteria be discarded from further analysis. We discuss the potential implications of this work for several areas in astrophysics.
△ Less
Submitted 16 January, 2018;
originally announced January 2018.
-
Disruption of Dark Matter Substructure: Fact or Fiction?
Authors:
Frank C. van den Bosch,
Go Ogiya,
Oliver Hahn,
Andreas Burkert
Abstract:
Accurately predicting the demographics of dark matter (DM) substructure is of paramount importance for many fields of astrophysics, including gravitational lensing, galaxy evolution, halo occupation modeling, and constraining the nature of dark matter. Because of its strongly non-linear nature, DM substructure is typically modeled using N-body simulations, which reveal that large fractions of DM s…
▽ More
Accurately predicting the demographics of dark matter (DM) substructure is of paramount importance for many fields of astrophysics, including gravitational lensing, galaxy evolution, halo occupation modeling, and constraining the nature of dark matter. Because of its strongly non-linear nature, DM substructure is typically modeled using N-body simulations, which reveal that large fractions of DM subhaloes undergo complete disruption. In this paper we use both analytical estimates and idealized numerical simulations to investigate whether this disruption is mainly physical, due to tidal heating and stripping, or numerical (i.e., artificial). We show that, contrary to naive expectation, subhaloes that experience a tidal shock $ΔE$ that exceeds the subhalo's binding energy, $|E_{\rm b}|$, do not undergo disruption, even when $ΔE/|E_{\rm b}|$ is as large as 100. Along the same line, and contrary to existing claims in the literature, instantaneously stripping matter from the outskirts of a DM subhalo also does not result in its complete disruption, even when the instantaneous remnant has positive binding energy. In addition, we show that tidal heating due to high-speed (impulsive) encounters with other subhaloes (`harassment'), is negligible compared to the tidal effects due to the host halo. Hence, we conclude that, in the absence of baryonic processes, the complete, physical disruption of CDM substructure is extremely rare, and that most disruption in numerical simulations therefore must be artificial. We discuss various processes that have been associated with numerical overmerging, and conclude that inadequate force-softening is the most likely culprit.
△ Less
Submitted 14 November, 2017;
originally announced November 2017.
-
What sets the central structure of dark matter haloes?
Authors:
Go Ogiya,
Oliver Hahn
Abstract:
Dark matter (DM) haloes forming near the thermal cut-off scale of the density perturbations are unique, since they are the smallest objects and form through monolithic gravitational collapse, while larger haloes contrastingly have experienced mergers. While standard cold dark matter (CDM) simulations readily produce haloes that follow the universal Navarro-Frenk-White (NFW) density profile with an…
▽ More
Dark matter (DM) haloes forming near the thermal cut-off scale of the density perturbations are unique, since they are the smallest objects and form through monolithic gravitational collapse, while larger haloes contrastingly have experienced mergers. While standard cold dark matter (CDM) simulations readily produce haloes that follow the universal Navarro-Frenk-White (NFW) density profile with an inner slope, $ρ\propto r^{-α}$, with $α=1$, recent simulations have found that when the free-streaming cut-off expected for the CDM model is resolved, the resulting haloes follow nearly power-law density profiles of $α\sim1.5$. In this paper, we study the formation of density cusps in haloes using idealized $N$-body simulations of the collapse of proto-haloes. When the proto-halo profile is initially cored due to particle free-streaming at high redshift, we universally find $\sim r^{-1.5}$ profiles irrespective of the proto-halo profile slope outside the core and large-scale non-spherical perturbations. Quite in contrast, when the proto-halo has a power-law profile, then we obtain profiles compatible with the NFW shape when the density slope of the proto-halo patch is shallower than a critical value, $α_{\rm ini} \sim 0.3$, while the final slope can be steeper for $α_{\rm ini}\ga 0.3$. We further demonstrate that the $r^{-1.5}$ profiles are sensitive to small scale noise, which gradually drives them towards an inner slope of $-1$, where they become resilient to such perturbations. We demonstrate that the $r^{-1.5}$ solutions are in hydrostatic equilibrium, largely consistent with a simple analytic model, and provide arguments that angular momentum appears to determine the inner slope.
△ Less
Submitted 9 October, 2017; v1 submitted 24 July, 2017;
originally announced July 2017.
-
Universal dark halo scaling relation for the dwarf spheroidal satellites
Authors:
Kohei Hayashi,
Tomoaki Ishiyama,
Go Ogiya,
Masashi Chiba,
Shigeki Inoue,
Masao Mori
Abstract:
Motivated by a recently found interesting property of the dark halo surface density within a radius, $r_{\rm max}$, giving the maximum circular velocity, $V_{\rm max}$, we investigate it for dark halos of the Milky Way's and Andromeda's dwarf satellites based on cosmological simulations. We select and analyze the simulated subhalos associated with Milky Way-sized dark halos and find that the value…
▽ More
Motivated by a recently found interesting property of the dark halo surface density within a radius, $r_{\rm max}$, giving the maximum circular velocity, $V_{\rm max}$, we investigate it for dark halos of the Milky Way's and Andromeda's dwarf satellites based on cosmological simulations. We select and analyze the simulated subhalos associated with Milky Way-sized dark halos and find that the values of their surface densities, $Σ_{V_{\rm max}}$, are in good agreement with those for the observed dwarf spheroidal satellites even without employing any fitting procedures. This implies that this surface density would not be largely affected by any baryonic feedbacks and thus universal. Moreover, all subhalos on the small scales of dwarf satellites are expected to obey the relation $Σ_{V_{\rm max}}\propto V_{\rm max}$, irrespective of differences in their orbital evolutions, host halo properties, and observed redshifts. Therefore, we find that the universal scaling relation for dark halos on dwarf galaxy mass scales surely exists and provides us important clues to understanding fundamental properties of dark halos. We also investigate orbital and dynamical evolutions of subhalos to understand the origin of this universal dark halo relation and find that most of subhalos evolve generally along the $r_{\rm max}\propto V_{\rm max}$ sequence, even though these subhalos have undergone different histories of mass assembly and tidal stripping. This sequence, therefore, should be the key feature to understand the nature of the universality of $Σ_{V_{\rm max}}$.
△ Less
Submitted 7 March, 2017;
originally announced March 2017.
-
Dynamical Evolution of Primordial Dark Matter Haloes through Mergers
Authors:
Go Ogiya,
Daisuke Nagai,
Tomoaki Ishiyama
Abstract:
Primordial dark matter (DM) haloes are the smallest gravitationally bound DM structures from which the first stars, black holes, and galaxies form and grow in the early universe. However, their structures are sensitive to the free streaming scale of DM, which in turn depends on the nature of DM particles. In this work, we test the hypothesis that the slope of the central cusps in primordial DM hal…
▽ More
Primordial dark matter (DM) haloes are the smallest gravitationally bound DM structures from which the first stars, black holes, and galaxies form and grow in the early universe. However, their structures are sensitive to the free streaming scale of DM, which in turn depends on the nature of DM particles. In this work, we test the hypothesis that the slope of the central cusps in primordial DM haloes near the free streaming scale depends on the nature of merging process. By combining and analysing data from a cosmological simulation with the cutoff in the small-scale matter power spectrum as well as a suite of controlled, high-resolution simulations of binary mergers, we find that (1) the primordial DM haloes form preferentially through major mergers in radial orbits; (2) their central DM density profile is more susceptible to a merging process compared to that of galaxy and cluster-size DM haloes; (3) consecutive major mergers drive the central density slope to approach the {\it universal} form characterized by the NFW profile, which is shown to be robust to the impacts of mergers and serves an attractor solution for the density structure of DM haloes. Our work highlights the importance of dynamical processes on the structure formation during the Dark Ages.
△ Less
Submitted 24 June, 2016; v1 submitted 11 April, 2016;
originally announced April 2016.
-
Dynamical friction and scratches of orbiting satellite galaxies on host systems
Authors:
Go Ogiya,
Andreas Burkert
Abstract:
We study the dynamical response of extended systems, hosts, to smaller systems, satellites, orbiting around the hosts using extremely high-resolution N-body simulations with up to one billion particles. This situation corresponds to minor mergers which are ubiquitous in the scenario of hierarchical structure formation in the universe. According to Chandrasekhar (1943), satellites create density wa…
▽ More
We study the dynamical response of extended systems, hosts, to smaller systems, satellites, orbiting around the hosts using extremely high-resolution N-body simulations with up to one billion particles. This situation corresponds to minor mergers which are ubiquitous in the scenario of hierarchical structure formation in the universe. According to Chandrasekhar (1943), satellites create density wakes along the orbit and the wakes cause a deceleration force on satellites, i.e. dynamical friction. This study proposes an analytical model to predict the dynamical response of hosts as reflected in their density distribution and finds not only traditional wakes but also mirror images of over- and underdensities centered on the host. Our controlled N-body simulations with high resolutions verify the predictions of the analytical model. We apply our analytical model to the expected dynamical response of nearby interacting galaxy pairs, the Milky Way - Large Magellanic Cloud system and the M31 - M33 system.
△ Less
Submitted 8 January, 2016; v1 submitted 27 October, 2015;
originally announced October 2015.
-
Re-examining the Too-Big-To-Fail Problem for Dark Matter Haloes with Central Density Cores
Authors:
Go Ogiya,
Andreas Burkert
Abstract:
Recent studies found the densities of dark matter (DM) subhaloes which surround nearby dwarf spheroidal galaxies (dSphs) to be significantly lower than those of the most massive subhaloes expected around Milky Way sized galaxies in cosmological simulations, the so called "too-big-to-fail" (TBTF) problem. A caveat of previous work has been that dark substructures were assumed to contain steep densi…
▽ More
Recent studies found the densities of dark matter (DM) subhaloes which surround nearby dwarf spheroidal galaxies (dSphs) to be significantly lower than those of the most massive subhaloes expected around Milky Way sized galaxies in cosmological simulations, the so called "too-big-to-fail" (TBTF) problem. A caveat of previous work has been that dark substructures were assumed to contain steep density cusps in the center of DM haloes even though the central density structure of DM haloes is still under debate. In this study, we re-examine the TBTF problem for models of DM density structure with cores or shallowed cusps. Our analysis demonstrates that the TBTF problem is alleviated as the logarithmic slope of the central cusp becomes shallower. We find that the TBTF problem is avoided if the central cusps of DM haloes surrounding dSphs are shallower than $r^{-0.6}$.
△ Less
Submitted 28 October, 2014; v1 submitted 27 August, 2014;
originally announced August 2014.
-
The connection between the cusp-to-core transformation and observational universalities of DM halos
Authors:
Go Ogiya,
Masao Mori,
Tomoaki Ishiyama,
Andreas Burkert
Abstract:
Observations have revealed interesting universal properties of dark matter (DM) halos especially around low-mass galaxies. Strigari et al. (2008) showed that DM halos have common enclosed masses within 300pc (Strigari relation). Kormendy & Freeman (2004) reported DM halos having almost identical central surface densities (the $μ_{\rm 0D}$ relation). In addition, there exists a core-cusp problem,…
▽ More
Observations have revealed interesting universal properties of dark matter (DM) halos especially around low-mass galaxies. Strigari et al. (2008) showed that DM halos have common enclosed masses within 300pc (Strigari relation). Kormendy & Freeman (2004) reported DM halos having almost identical central surface densities (the $μ_{\rm 0D}$ relation). In addition, there exists a core-cusp problem, a discrepancy of the central density distribution between simulated halos and observations. We investigate whether a scenario where cuspy halos transform into cores by some dynamical processes can also explain their universal structural properties. It is shown that a cusp-to-core transformation model naturally reproduces the $μ_{\rm 0D}$ relation and that Strigari relation follows from the $μ_{\rm 0D}$ relation for dwarf galaxies. We also show that the central densities of cored dark halos provide valuable information about their formation redshifts.
△ Less
Submitted 4 February, 2014; v1 submitted 6 September, 2013;
originally announced September 2013.
-
The Core-Cusp Problem in Cold Dark Matter Halos and Supernova Feedback: Effects of Oscillation
Authors:
Go Ogiya,
Masao Mori
Abstract:
This study investigates the dynamical response of dark matter (DM) halos to recurrent starbursts in forming less-massive galaxies to solve the core-cusp problem. The gas, which is heated by supernova feedback after a starburst, expands and the star formation then terminates. This expanding gas loses energy by radiative cooling and then falls back toward the galactic center. Subsequently, the starb…
▽ More
This study investigates the dynamical response of dark matter (DM) halos to recurrent starbursts in forming less-massive galaxies to solve the core-cusp problem. The gas, which is heated by supernova feedback after a starburst, expands and the star formation then terminates. This expanding gas loses energy by radiative cooling and then falls back toward the galactic center. Subsequently, the starburst is enhanced again. This cycle of expansion and contraction of the interstellar gas leads to a repetitive change in the gravitational potential of the gas. The resonance between DM particles and the density wave excited by the oscillating potential plays a key role in understanding the physical mechanism of the cusp-core transition of DM halos. DM halos effectively gain kinetic energy from the baryon potential through the energy transfer driven by the resonance between the particles and density waves. We determine that the critical condition for the cusp-core transition is such that the oscillation period of the gas potential is approximately the same as the local dynamical time of DM halos. We present the resultant core radius of a DM halo after the cusp-core transition induced by the resonance by using the conventional mass density profile predicted by the cold dark matter models. Moreover, we verify the analytical model by using $N$-body simulations, and the results validate the resonance model.
△ Less
Submitted 28 October, 2014; v1 submitted 23 June, 2012;
originally announced June 2012.
-
Astrophysical Particle Simulations on Heterogeneous CPU-GPU Systems
Authors:
Naohito Nakasato,
Go Ogiya,
Yohei Miki,
Masao Mori,
Ken'ichi Nomoto
Abstract:
A heterogeneous CPU-GPU node is getting popular in HPC clusters. We need to rethink algorithms and optimization techniques for such system depending on the relative performance of CPU vs. GPU. In this paper, we report a performance optimized particle simulation code "OTOO", that is based on the octree method, for heterogenous systems. Main applications of OTOO are astrophysical simulations such as…
▽ More
A heterogeneous CPU-GPU node is getting popular in HPC clusters. We need to rethink algorithms and optimization techniques for such system depending on the relative performance of CPU vs. GPU. In this paper, we report a performance optimized particle simulation code "OTOO", that is based on the octree method, for heterogenous systems. Main applications of OTOO are astrophysical simulations such as N-body models and the evolution of a violent merger of stars. We propose optimal task split between CPU and GPU where GPU is only used to compute the calculation of the particle force. Also, we describe optimization techniques such as control of the force accuracy, vectorized tree walk, and work partitioning among multiple GPUs. We used OTOO for modeling a merger of two white dwarf stars and found that OTOO is powerful and practical to simulate the fate of the process.
△ Less
Submitted 6 June, 2012;
originally announced June 2012.
-
The core-cusp problem in cold dark matter halos and supernova feedback: Effects of Mass Loss
Authors:
Go Ogiya,
Masao Mori
Abstract:
The core-cusp problem remains as one of the unsolved discrepancies between observations and theories predicted by the standard paradigm of cold dark matter (CDM) cosmology. To solve this problem, we perform N-body simulations to study the nonlinear response of CDM halos to the variance of the gravitational potential induced by gas removal from galaxy centers. In this study, we focus on the timesca…
▽ More
The core-cusp problem remains as one of the unsolved discrepancies between observations and theories predicted by the standard paradigm of cold dark matter (CDM) cosmology. To solve this problem, we perform N-body simulations to study the nonlinear response of CDM halos to the variance of the gravitational potential induced by gas removal from galaxy centers. In this study, we focus on the timescale of the gas ejection, which is strongly correlated with stellar activities, and demonstrate that it is one of the key factors in determining the dynamical response of CDM halos. The results of simulations show that the power-low index of the mass-density profile of the dark matter halo correlated with the timescale of the mass loss, and it is flatter when the mass loss occurs over a short time than when it occurs over a long time. However, it is still larger than typical observational values; in other words, the central cusp remains for any mass loss model in the simulations. Moreover, for the slow mass-loss case, the final density profile of the dark matter halo recovers the universal density profiles predicted by the CDM cosmology. Therefore, mass loss driven by stellar feedback may not be an effective mechanism to flatten the central cusp.
△ Less
Submitted 14 June, 2011;
originally announced June 2011.