Researcher profile

Frank C. van den Bosch

Frank C. van den Bosch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

17 published item(s)

preprint2026arXiv

Universal numerical convergence criteria for subhalo tidal evolution

Dark matter subhalos and satellite galaxies in state-of-the-art cosmological simulations still suffer from the ``overmerging'' problem, where inadequate force and/or mass resolution cause artificially enhanced tidal mass loss and premature disruption. Previous idealized simulations addressing this issue have been restricted to a small subset of the subhalo orbital parameter space, and all assumed subhalos to be isotropic. Here, we present the first extensive simulation suite that quantifies numerical convergence in the tidal evolution of anisotropic subhalos under varying numerical resolutions and orbits. We report a universal force resolution criterion: the subhalo's instantaneous tidal radius must always be resolved by at least 20 cells in adaptive mesh refinement (AMR)-based simulations, or by 20 softening lengths (Plummer equivalent) in tree-based simulations, regardless of refinement details or subhalo physical properties such as concentration or velocity anisotropy. We also report a universal expression for the discreteness-noise-driven scatter in the bound-mass fraction of subhalos that depends only on the subhalo mass resolution at infall and the instantaneous bound mass fraction, agnostic of any further subhalo properties. Such stochastic discreteness noise causes both premature disruption and, notably, spurious survival of poorly mass-resolved subhalos. We demonstrate that as many as 50 percent of all subhalos in state-of-the-art cosmological simulations are likely to be either force and/or mass unresolved. Our findings advocate for adaptive softening or grid refinement based on the instantaneous tidal radius of the subhalo.

preprint2022arXiv

A Comprehensive Perturbative Formalism for Phase-Mixing in Perturbed Disks. I. Phase spirals in an Infinite, Isothermal Slab

Galactic disks are highly responsive systems that often undergo external perturbations and subsequent collisionless equilibration, predominantly via phase-mixing. We use linear perturbation theory to study the response of infinite isothermal slab analogues of disks to perturbations with diverse spatio-temporal characteristics. Without self-gravity of the response, the dominant Fourier modes that get excited in a disk are the bending and breathing modes, which, due to vertical phase-mixing, trigger local phase-space spirals that are one- and two-armed, respectively. We demonstrate how the lateral streaming motion of slab stars causes phase spirals to damp out over time. The ratio of the perturbation timescale ($τ_{\mathrm{P}}$) to the local, vertical oscillation time ($τ_z$) ultimately decides which of the two modes is excited. Faster, more impulsive ($τ_{\mathrm{P}} < τ_z$) and slower, more adiabatic ($τ_{\mathrm{P}} > τ_z$) perturbations excite stronger breathing and bending modes, respectively, although the response to very slow perturbations is exponentially suppressed. For encounters with satellite galaxies, this translates to more distant and more perpendicular encounters triggering stronger bending modes. We compute the direct response of the Milky Way disk to several of its satellite galaxies, and find that recent encounters with all of them excite bending modes in the Solar neighborhood. The encounter with Sagittarius triggers a response that is at least $1-2$ orders of magnitude larger than that due to any other satellite, including the Large Magellanic Cloud. We briefly discuss how ignoring the presence of a dark matter halo and the self-gravity of the response might impact our conclusions.

preprint2022arXiv

Quenching in Cosmic Sheets: Tracing the Impact of Large Scale Structure Collapse on the Evolution of Dwarf Galaxies

Dwarf galaxies are thought to quench primarily due to environmental processes most typically occurring in galaxy groups and clusters or around single, massive galaxies. However, at earlier epochs, ($5 < z < 2$), the collapse of large scale structure (forming Zel&#39;dovich sheets and subsequently filaments of the cosmic web) can produce volume-filling accretion shocks which elevate large swaths of the intergalactic medium (IGM) in these structures to a hot ($T>10^6$ K) phase. We study the impact of such an event on the evolution of central dwarf galaxies ($5.5 < \log M_* < 8.5$) in the field using a spatially large, high resolution cosmological zoom simulation which covers the cosmic web environment between two protoclusters. We find that the shock-heated sheet acts as an environmental quencher much like clusters and filaments at lower redshift, creating a population of quenched, central dwarf galaxies. Even massive dwarfs which do not quench are affected by the shock, with reductions to their sSFR and gas accretion. This process can potentially explain the presence of isolated quenched dwarf galaxies, and represents an avenue of pre-processing, via which quenched satellites of bound systems quench before infall.

preprint2021arXiv

Dynamical Friction, Buoyancy and Core-Stalling -- I. A Non-perturbative Orbit-based Analysis

We examine the origin of dynamical friction using a non-perturbative, orbit-based approach. Unlike the standard perturbative approach, in which dynamical friction arises from the LBK torque due to pure resonances, this alternative, complementary view nicely illustrates how a massive perturber significantly changes the energies and angular momenta of field particles on near-resonant orbits, with friction arising from an imbalance between particles that gain energy and those that lose energy. We treat dynamical friction in a spherical host system as a restricted three-body problem. This treatment is applicable in the `slow&#39; regime, in which the perturber sinks slowly and the standard perturbative framework fails due to the onset of non-linearities. Hence it is especially suited to investigate the origin of core-stalling: the cessation of dynamical friction in central constant-density cores. We identify three different families of near-co-rotation-resonant orbits that dominate the contribution to dynamical friction. Their relative contribution is governed by the Lagrange points (fixed points in the co-rotating frame). In particular, one of the three families, which we call Pac-Man orbits because of their appearance in the co-rotating frame, is unique to cored density distributions. When the perturber reaches a central core, a bifurcation of the Lagrange points drastically changes the orbital make-up, with Pac-Man orbits becoming dominant. In addition, due to relatively small gradients in the distribution function inside a core, the net torque from these Pac-Man orbits becomes positive (enhancing), thereby effectuating a dynamical buoyancy. We argue that core stalling occurs where this buoyancy is balanced by friction.

preprint2021arXiv

SatGen -- II. Assessing the impact of a disc potential on subhalo populations

The demographics of dark matter substructure depend sensitively on the nature of dark matter. Optimally leveraging this probe requires accurate theoretical predictions regarding the abundance of subhaloes. These predictions are hampered by artificial disruption in numerical simulations, by large halo-to-halo variance, and by the fact that the results depend on the baryonic physics of galaxy formation. In particular, numerical simulations have shown that the formation of a central disc can drastically reduce the abundance of substructure compared to a dark matter-only simulation, which has been attributed to enhanced destruction of substructure due to disc shocking. We examine the impact of discs on substructure using the semi-analytical subhalo model SatGen, which accurately models the tidal evolution of substructure free of the numerical disruption that still hampers $N$-body simulations. Using a sample of 10,000 merger trees of Milky Way-like haloes, we study the demographics of subhaloes that are evolved under a range of composite halo-disc potentials with unprecedented statistical power. We find that the overall subhalo abundance is relatively insensitive to properties of the disc aside from its total mass. For a disc that contains $5\%$ of $M_\mathrm{vir}$, the mean subhalo abundance within $r_\mathrm{vir}$ is suppressed by ${\lesssim}10\%$ relative to the no-disc case, a difference that is dwarfed by halo-to-halo variance. For the same disc mass, the abundance of subhaloes within 50 kpc is reduced by ${\sim}30\%$. We argue that the disc mainly drives excess mass loss for subhaloes with small pericentric radii and that the impact of disc shocking is negligible.

preprint2020arXiv

Dynamical self-friction: how mass loss slows you down

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&#39;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&#39;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.

preprint2020arXiv

Scatter in Sunyaev--Zel&#39;dovich effect scaling relations explained by inter-cluster variance in mass accretion histories

X-ray and microwave cluster scaling relations are immensely valuable for cosmological analysis. However, their power is limited by astrophysical systematics that bias mass estimates and introduce additional scatter. Turbulence injected into the intracluster medium via mass assembly contributes substantially to cluster non-thermal pressure support, a significant source of such uncertainties. We use an analytical model to compute the assembly-driven non-thermal pressure profiles of haloes based on Monte Carlo-generated accretion histories. We introduce a fitting function for the average non-thermal pressure fraction profile, which exhibits minimal dependence on redshift at fixed peak height. Using the model, we predict deviations from self-similarity and the intrinsic scatter in the Sunyaev--Zel&#39;dovich effect observable-mass scaling relation ($Y_\mathrm{SZ}-M$) due solely to inter-cluster variation in mass accretion histories. We study the dependence of $Y_\mathrm{SZ}-M$ on aperture radius, cosmology, redshift, and mass limit. The model predicts $5-9\%$ scatter in $Y_\mathrm{SZ}-M$ at $z=0$, increasing as the aperture used to compute $Y_\mathrm{SZ}$ increases from $R_\mathrm{500c}$ to $5R_\mathrm{500c}$. The predicted scatter lies slightly below that of studies based on non-radiative hydro-simulations, illustrating that assembly history variance is likely responsible for a substantial fraction of scatter in $Y_\mathrm{SZ}-M$. This should be regarded as a lower bound, which will likely increase with the use of an updated gas density model that incorporates a more realistic response to halo assembly. As redshift increases, $Y_\mathrm{SZ}-M$ deviates more from self-similarity and scatter increases. We show that the $Y_\mathrm{SZ}-M$ residuals correlate strongly with the recent halo mass accretion rate, potentially providing an opportunity to infer the latter.

preprint2019arXiv

How to Optimally Constrain Galaxy Assembly Bias: Supplement Projected Correlation Functions with Count-in-cells Statistics

Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy--halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function $w_{\mathrm{p}}(r_{\mathrm{p}})$, the galaxy--galaxy lensing signal $ΔΣ(r_{\mathrm{p}})$, the void probability function $\mathrm{VPF}(r)$, the distributions of counts-in-cylinders $P(N_{\mathrm{CIC}})$, and counts-in-annuli $P(N_{\mathrm{CIA}})$, and the distribution of the ratio of counts in cylinders of different sizes $P(N_2/N_5)$. We find that despite the frequent use of the combination $w_{\mathrm{p}}(r_{\mathrm{p}})+ΔΣ(r_{\mathrm{p}})$ in interpreting galaxy data, the count statistics, $P(N_{\mathrm{CIC}})$ and $P(N_{\mathrm{CIA}})$, are generally more efficient in constraining galaxy assembly bias when combined with $w_{\mathrm{p}}(r_{\mathrm{p}})$. Constraints based upon $w_{\mathrm{p}}(r_{\mathrm{p}})$ and $ΔΣ(r_{\mathrm{p}})$ share common degeneracy directions in the parameter space, while combinations of $w_{\mathrm{p}}(r_{\mathrm{p}})$ with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy--halo connection.

preprint2019arXiv

The Dearth of Difference between Central and Satellite Galaxies III. Environmental Dependence of Mass-Size and Mass-Structure Relations

As demonstrated in Paper I, the quenching properties of central and satellite galaxies are quite similar as long as both stellar mass and halo mass are controlled. Here we extend the analysis to the size and bulge-to-total light ratio (B/T) of galaxies. In general central galaxies have size-stellar mass and B/T-stellar mass relations different from satellites. However, the differences are eliminated when halo mass is controlled. We also study the dependence of size and B/T on halo-centric distance and find a transitional stellar mass (M$_{*,t}$) at given halo mass (M$_h$), which is about one fifth of the mass of the central galaxies in halos of mass M$_h$. The transitional stellar masses for size, B/T and quenched fraction are similar over the whole halo mass range, suggesting a connection between the quenching of star formation and the structural evolution of galaxies. Our analysis further suggests that the classification based on the transitional stellar mass is more fundamental than the central-satellite dichotomy, and provide a more reliable way to understand the environmental effects on galaxy properties. We compare the observational results with the hydro-dynamical simulation, EAGLE and the semi-analytic model, L-GALAXIES. The EAGLE simulation successfully reproduces the similarities of size for centrals and satellites and even M$_{*,t}$, while L-GALAXIES fails to recover the observational results.

preprint2018arXiv

Maturing Satellite Kinematics into a Competitive Probe of the Galaxy-Halo Connection

The kinematics of satellite galaxies moving in a dark matter halo are a direct probe of the underlying gravitational potential. Thus, the phase-space distributions of satellites represent a powerful tool to determine the galaxy-halo connection from observations. By stacking the signal of a large number of satellite galaxies this potential can be unlocked even for haloes hosting a few satellites on average. In this work, we test the impact of various modelling assumptions on constraints derived from analysing satellite phase-space distributions in the non-linear, 1-halo regime. We discuss their potential to explain the discrepancy between average halo masses derived from satellite kinematics and gravitational lensing previously reported. Furthermore, we develop an updated, more robust analysis to extract constraints on the galaxy-halo relation from satellite properties in spectroscopic galaxy surveys such as the SDSS. We test the accuracy of this approach using a large number of realistic mock catalogues. Furthermore, we find that constraints derived from such an analysis are complementary and competitive with respect to the commonly used galaxy clustering and galaxy-galaxy lensing observables.

preprint2018arXiv

Updated Results on the Galaxy-Halo Connection from Satellite Kinematics in SDSS

We present new results on the relationship between central galaxies and dark matter haloes inferred from observations of satellite kinematics in the Sloan Digital Sky Survey (SDSS) DR7. We employ an updated analysis framework that includes detailed mock catalogues to model observational effects in SDSS. Our results constrain the colour-dependent conditional luminosity function (CLF) of dark matter haloes, as well as the radial profile of satellite galaxies. Confirming previous results, we find that red central galaxies live in more massive haloes than blue galaxies at fixed luminosity. Additionally, our results suggest that satellite galaxies have a radial profile less centrally concentrated than dark matter but not as cored as resolved subhaloes in dark matter-only simulations. Compared to previous works using satellite kinematics by More et al., we find much more competitive constraints on the galaxy-halo connection, on par with those derived from a combination of galaxy clustering and galaxy-galaxy lensing. We compare our results on the galaxy-halo connection to other studies using galaxy clustering and group catalogues, showing very good agreement between these different techniques. We discuss future applications of satellite kinematics in the context of constraining cosmology and the relationship between galaxies and dark matter haloes.

preprint2017arXiv

Brightest galaxies as halo centre tracers in SDSS DR7

Determining the positions of halo centres in large-scale structure surveys is crucial for many cosmological studies. A common assumption is that halo centres correspond to the location of their brightest member galaxies. In this paper, we study the dynamics of brightest galaxies with respect to other halo members in the Sloan Digital Sky Survey DR7. Specifically, we look at the line-of-sight velocity and spatial offsets between brightest galaxies and their neighbours. We compare those to detailed mock catalogues, constructed from high-resolution, dark-matter-only $N$-body simulations, in which it is assumed that satellite galaxies trace dark matter subhaloes. This allows us to place constraints on the fraction $f_{\rm BNC}$ of haloes in which the brightest galaxy is not the central. Compared to previous studies we explicitly take into account the unrelaxed state of the host haloes, velocity offsets of halo cores and correlations between $f_{\rm BNC}$ and the satellite occupation. We find that $f_{\rm BNC}$ strongly decreases with the luminosity of the brightest galaxy and increases with the mass of the host halo. Overall, in the halo mass range $10^{13} - 10^{14.5} h^{-1} M_\odot$ we find $f_{\rm BNC} \sim 30\%$, in good agreement with a previous study by Skibba et al. We discuss the implications of these findings for studies inferring the galaxy--halo connection from satellite kinematics, models of the conditional luminosity function and galaxy formation in general.

preprint2017arXiv

The Galaxy Clustering Crisis in Abundance Matching

Galaxy clustering on small scales is significantly under-predicted by sub-halo abundance matching (SHAM) models that populate (sub-)haloes with galaxies based on peak halo mass, $M_{\rm peak}$. SHAM models based on the peak maximum circular velocity, $V_{\rm peak}$, have had much better success. The primary reason $M_{\rm peak}$ based models fail is the relatively low abundance of satellite galaxies produced in these models compared to those based on $V_{\rm peak}$. Despite success in predicting clustering, a simple $V_{\rm peak}$ based SHAM model results in predictions for galaxy growth that are at odds with observations. We evaluate three possible remedies that could &#34;save&#34; mass-based SHAM: (1) SHAM models require a significant population of &#34;orphan&#34; galaxies as a result of artificial disruption/merging of sub-haloes in modern high resolution dark matter simulations; (2) satellites must grow significantly after their accretion; and (3) stellar mass is significantly affected by halo assembly history. No solution is entirely satisfactory. However, regardless of the particulars, we show that popular SHAM models based on $M_{\rm peak}$ cannot be complete physical models as presented. Either $V_{\rm peak}$ truly is a better predictor of stellar mass at $z\sim 0$ and it remains to be seen how the correlation between stellar mass and $V_{\rm peak}$ comes about, or SHAM models are missing vital component(s) that significantly affect galaxy clustering.

preprint2017arXiv

The Immitigable Nature of Assembly Bias: The Impact of Halo Definition on Assembly Bias

Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $Δ$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $Δ$. For conventional halo definitions ($Δ\sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $Δ\sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $Δ$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $Δ\gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass.

preprint2016arXiv

Constraints on Assembly Bias from Galaxy Clustering

We constrain the newly-introduced decorated Halo Occupation Distribution (HOD) model using SDSS DR7 measurements of projected galaxy clustering or r-band luminosity threshold samples. The decorated HOD is a model for the galaxy-halo connection that augments the HOD by allowing for the possibility of galaxy assembly bias: galaxy luminosity may be correlated with dark matter halo properties besides mass, Mvir. We demonstrate that it is not possible to rule out galaxy assembly bias using DR7 measurements of galaxy clustering alone. Moreover, galaxy samples with Mr < -20 and Mr < -20.5 favor strong central galaxy assembly bias. These samples prefer scenarios in which high-concentration are more likely to host a central galaxy relative to low-concentration halos of the same mass. We exclude zero assembly bias with high significance for these samples. Satellite galaxy assembly bias is significant for the faintest sample, Mr < -19. We find no evidence for assembly bias in the Mr < -21 sample. Assembly bias should be accounted for in galaxy clustering analyses or attempts to exploit galaxy clustering to constrain cosmology. In addition to presenting the first constraints on HOD models that accommodate assembly bias, our analysis includes several improvements over previous analyses of these data. Therefore, our inferences supersede previously-published results even in the case of a standard HOD analysis.

preprint2009arXiv

Constraints on the relationship between stellar mass and halo mass at low and high redshift

We use a statistical approach to determine the relationship between the stellar masses of galaxies and the masses of the dark matter halos in which they reside. We obtain a parameterized stellar-to-halo mass (SHM) relation by populating halos and subhalos in an N-body simulation with galaxies and requiring that the observed stellar mass function be reproduced. We find good agreement with constraints from galaxy-galaxy lensing and predictions of semi-analytic models. Using this mapping, and the positions of the halos and subhalos obtained from the simulation, we find that our model predictions for the galaxy two-point correlation function (CF) as a function of stellar mass are in excellent agreement with the observed clustering properties in the SDSS at z=0. We show that the clustering data do not provide additional strong constraints on the SHM function and conclude that our model can therefore predict clustering as a function of stellar mass. We compute the conditional mass function, which yields the average number of galaxies with stellar masses in the range [m, m+dm] that reside in a halo of mass M. We study the redshift dependence of the SHM relation and show that, for low mass halos, the SHM ratio is lower at higher redshift. The derived SHM relation is used to predict the stellar mass dependent galaxy CF and bias at high redshift. Our model predicts that not only are massive galaxies more biased than low mass ones at all redshifts, but the bias increases more rapidly with increasing redshift for massive galaxies than for low mass ones. We present convenient fitting functions for the SHM relation as a function of redshift, the conditional mass function, and the bias as a function of stellar mass and redshift.

preprint2008arXiv

New Constraints on the Efficiencies of Ram-Pressure Stripping and the Tidal Disruption of Satellite Galaxies

Using data from the Sloan Digital Sky Survey (SDSS) it has recently been shown that the red fraction of satellite galaxies increases with stellar mass. Semi-analytical models, however, predict red satellite fractions that are independent of stellar mass, and much higher than observed. It has been argued that this discrepancy owes to the fact that the models assume that satellite galaxies are instantaneously stripped of their hot gas reservoirs at the moment they are accreted into a bigger halo. In this letter we show that the fraction of red satellites can be brought in better agreement with the data by simply decreasing this stripping efficiency. However, this also results in a red fraction of massive centrals that is much too low. This owes to the fact that the massive centrals now accrete satellite galaxies that are bluer and more gas-rich. However, if a significant fraction of low mass satellite galaxies is tidally disrupted before being accreted by their central host galaxy, as suggested by recent studies, the red fractions of both centrals and satellites can be reproduced reasonably well. A problem remains with the red fraction of centrals of intermediate mass, which is likely to reflect an oversimplified treatment of AGN feedback.