Impact of subhalo dynamical friction heating on the formation of the first structures in the universe

arXiv:2603.19385 · astro-ph.GA, astro-ph.CO · Submitted 2026-03-19 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Today's paper: "Impact of subhalo dynamical friction heating on the formation of the first structures in the universe".

Jocelyn: The study investigates the impact of subhalo dynamical friction heating on structure formation, analyzing various physical processes and statistical distributions within massive halos across cosmic time.

Vera: First, who's behind it and why it matters.

Title and authors: Vera: We’re starting with the title and authors for this paper, "Impact of subhalo dynamical friction heating on the formation of the first structures in the universe." It sounds very technical, focusing on how these smaller clumps affect big structure growth.

Jocelyn: That title suggests they are looking at something specific related to those early structures—the first ones forming in the universe—and connecting it to dynamical friction heating. What kind of observations or simulations are they using as their foundation?

Subrahmanyan: The authors are Zhenyu Wu, Sadegh Khochfar, Muhammad A. Latif, Ben Morton, and Britton Smith from the Institute for Astronomy at the University of Edinburgh. They’re using data from the TNG50 simulation to analyze how dynamical friction acts on gas heating within dark matter halos across a wide range of redshifts.

Vera: That's right, and they're using TNG50 simulation data to calculate the subhalo mass function and derive that dynamical friction heating rate for many halo masses from redshift fifteen all the way to zero. It gives us a lot of context on structure formation history.

Jocelyn: So, if I understand correctly, this paper is taking those complex simulation results—the subhalo mass function and the heating rates—and trying to make sense of how that energy transfer happens across cosmic time?

Subrahmanyan: Precisely; they’re converting gravitational potential energy into thermal energy through dynamical friction from orbiting subhalos within dark matter halos, and their main finding is that this process is an important mechanism for galaxy quenching in massive halos at low redshifts, which aligns with previous studies.

Vera: That connection to quenching in massive halos at lower redshifts seems like a really significant link between the tiny scale dynamics of subhalos and the larger galaxy evolution we observe. It makes the physics feel much more connected than just looking at isolated events.

Jocelyn: It’s powerful when you see how these small-scale interactions feed into the big picture of how galaxies stop forming stars over time, which is what this paper seems to be highlighting in detail.

Subrahmanyan: That's exactly it; it helps build the overall picture by showing that heating isn't just one thing; it reinforces existing theories about how galaxies evolve within their cosmic environment.

The paper's summary: Vera: Moving into the main body of the "Impact of subhalo dynamical friction heating on the formation of the first structures in the universe," we see they’re digging into several complex areas, starting with how they look at statistical scatter in those subhalo mass functions.

Jocelyn: The paper explores non-Poissonian fluctuations in that mass function, which usually means looking for things that aren't just random noise. What kind of statistical deviations are they actually finding there?

Subrahmanyan: They reference prior work showing that in Nbody simulations, the square root scatter of the subhalo mass function becomes super-Poissonian when the mass ratio psi is smaller than five times ten to the three, but then it shifts to a sub-Poissonian nature at higher mass ratios.

Vera: So they’re tracking a transition point where this noise changes its statistical behavior depending on how massive the subhalo is relative to the host halo. That shift from super- to sub-Poissonian scatter seems like a really detailed piece of information.

Jocelyn: And what does that specific transition mean for their overall calculation of the dynamical friction heating rate, since they have to use this mass function as an input?

Subrahmanyan: Because the deviation in this regime is only mild when psi is greater than ten to the negative two, they decide to assume Poisson fluctuations for estimating that heating rate, even though they’ve seen these non-Poissonian effects.

Vera: So they are making a pragmatic choice to use the simpler Poisson model for their main calculation because trying to model the non-Poissonian nature fully is too complex given the current data constraints.

Jocelyn: That makes sense; it shows they’re balancing deep statistical analysis with the need to deliver a clear physical result about how much energy is being injected into these halos.

Subrahmanyan: They are deferring a systematic investigation of that non-Poissonian nature to future work, which is fair given the current data constraints and the focus on dynamical friction heating itself.

The paper's improvements: Vera: Now, let’s look at what the authors suggest to refine their model, focusing on how they handle the distribution of subhalo Mach numbers, which is crucial for the heating calculation.

Jocelyn: I'm curious if they found that their initial choice for modeling these Mach numbers was inadequate? Did they find a flaw in how they were describing those supersonic subhalos?

Subrahmanyan: They started with a Maxwell-Boltzmann distribution, but that distribution didn't capture the peaks of the higher mass groups accurately, which is what prompted them to change their approach.

Vera: So they moved away from that initial fit and switched to using a Gaussian distribution instead, and they normalized it by Equation (C1) after truncating it at zero mass. That sounds like a solid step toward getting a better fit for the supersonic subhalos.

Jocelyn: And what did this new fitting procedure reveal about the physical properties of these subhalos as you move up in host halo mass? Did anything else change besides just the fit quality?

Subrahmanyan: The analysis showed that both the mean Mach number mu and the standard deviation sigma increase as the host halo mass increases, which is a trend they are tracking.

Vera: And they also noted something about the average Mach number correction factor h I/M i, which is that for a fixed host halo mass, lower redshift ones show a larger scatter in subhalo Mach numbers.

Jocelyn: But within the range of best-fit standard deviation of about zero point two to zero point five, they didn't find that the average Mach number correction factor changed much across those redshifts, did they?

Subrahmanyan: That’s right; within that specific sigma range, the average Mach number correction factor doesn't have a significant change across those redshifts.

Conclusion: Vera: So we’ve covered a lot about the paper "Impact of subhalo dynamical friction heating on the formation of the first structures in the universe," covering everything from statistical scatter to model refinements for Mach numbers and density profiles.

Jocelyn: To wrap up, what is the main message we should be telling our listeners about this research? What’s the biggest implication of this research?

Subrahmanyan: The main implication is that dynamical friction heating from subhalos is an important mechanism for galaxy quenching in massive halos at low redshifts, and it’s consistent with previous studies.

Vera: So, the big takeaway for our audience is that even when we have complex statistical variations in subhalo distribution, the fundamental physics of energy conversion through dynamical friction still points to this as a valid way to quench structures.

Jocelyn: It’s powerful because it grounds this astrophysical concept in concrete calculations derived from simulations and shows us exactly how much energy is being injected into the system.

Subrahmanyan: This work helps build the overall picture by showing that heating isn't just one thing; it reinforces existing theories about how galaxies evolve within their cosmic environment.

Vera: It really solidifies the role of these subhalos as crucial contributors to galaxy evolution across vast stretches of cosmic time.

Jocelyn: I’m just excited to see what other papers come out next that explore these ideas, but it sounds like a lot of foundational work has been laid down here.

Subrahmanyan: Indeed, this paper lays down the groundwork for understanding how subhalo dynamics feed into the larger story of structure formation.

Vera: We'll be ready for whatever comes next on our schedule; thank you all for joining us today.

Zhenyu Wu, Sadegh Khochfar, Muhammad A. Latif, Ben Morton, Britton Smith

Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK · Physics Department, College of Science, United Arab Emirates University

astro-ph.GA, astro-ph.CO

Submitted: 2026-03-19

Updated: 2026-08-15

Comments: 19 pages, 23 figures. Submitted to MNRAS

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 91/100

The gist: The study investigates the impact of subhalo dynamical friction heating on structure formation, analyzing various physical processes and statistical distributions within massive halos across cosmic

Key concepts

Subhalo Dynamical Friction Heating
This process describes how gravitational potential energy from orbiting subhalos within dark matter halos is converted into thermal energy through dynamical friction. This heating affects the gas within the halo.
Subhalo Mass Function Scatter
The study examines statistical scatter in the distribution of subhalo masses. They found that this scatter changes its statistical behavior depending on the mass ratio between a subhalo and its host halo.
Galaxy Quenching
This refers to the process where galaxies stop forming stars. The paper suggests that dynamical friction heating from subhalos is an important mechanism contributing to galaxy quenching in massive halos at low redshifts.

Terminology

Summary

The study investigates the impact of subhalo dynamical friction heating on structure formation, analyzing various physical processes and statistical distributions within massive halos across cosmic time.

Host Halo Selection and Subhalo Mass Function (SHMF) Analysis:

For host halo selection, the authors target masses above 10 5.5 M/h, which is noted as roughly the threshold mass for Pop III star formation. Specific selections were made, yielding 14 host halos with M vir in [10 5.5, 10 6.7], M/h at z=12, and 10 host halos with M vir in [10 5.5, 10 6.2], M/h at z=15. The analysis compares the SHMF for individual host halos (shown as grey histograms) against the mean SHMF (blue histogram). In terms of the mass ratio psi, the Pop2Prime SHMF appears to lie between the high-redshift TNG result and the z=0 TNG result in the higher mass-ratio regime (psi about 10-1). Conversely, In the lower- psi regime, the slope of the TNG z=0 SHMF shows better agreement. While fitting with Equation 4 is suggested for minihalos, the authors caution that due to the limited number of host halos, the data lack the statistical power required to provide a precise fit.

Non-Poissonian Fluctuations and Statistical Assumptions:

The analysis addresses non-Poisson fluctuations in SHMF. It references prior work stating that in Nbody simulations, "the sqrt scatter of SHMF becomes super-Poissonian, i.e. sigma > sigma Poisson = N, when the mass ratio psi is smaller than 5 times 10-3. Furthermore, JvBIII further find[s] that at the higher mass ratio end, or when the average number of subhalos becomes small (h N (psi) i < 2), the scatter of SHMF has a sub-Poissonian nature, i.e. sigma < sigma Poisson. A transition in scatter behavior is observed: Figure 6, we observe that the scatter of the SHMF transitions from super-Poissonian to sub-Poissonian as psi increases, with h N (psi) i about 2 approximately marking the boundary. Given that dynamical friction heating is predominantly contributed by subhalos with psi > 10-2, and the deviation in this regime is only mild, the authors conclude: We therefore assume Poisson fluctuations for estimating the impact of SHMF scatter on the dynamical-friction heating rate," deferring a systematic investigation of non-Poissonian nature to future work.

Mach Number Distribution Fitting:

The study compares different fitting procedures for subhalo Mach numbers. While a Maxwell-Boltzmann distribution was initially used, it was found that the Maxwell-Boltzmann distribution doesn’t capture the peaks of the higher mass groups accurately. Consequently, the authors present fitting results using a Gaussian distribution, which is truncated at M=0 and normalized by Equation (C1). Analysis of this fit shows that both mu and sigma increase with host halo mass. Regarding the average Mach number correction factor (h I/M i), it is noted that for a fixed host halo mass, the lower redshift ones exhibit a larger scatter of subhalo Mach numbers. However, within the best-fit sigma range (about 0.2 - 0.5), the average Mach number correction h I/M i doesn’t have a significant change. Overall, the authors conclude that different fitting procedures will not qualitatively change our results in the paper.

Cooling and Heating Profiles:

The investigation also compares cooling and heating processes based on halo gas density profiles. The default profile used is the core profile (alpha = 0). When comparing cooling and heating using a cuspy NFW profile (alpha = 1), the authors find that "The cuspy profile leads to a higher cooling rate by a factor of a few, but doesn’

Improvements for AI systems

The following improvements are designed to enhance an AI system's ability to model complex astrophysical phenomena, moving beyond simple data correlation toward dynamic, physically constrained predictive modeling, as demonstrated by the findings in this paper.

Improvement: Implement a probabilistic energy balance module that treats Dynamical Friction (DF) heating not as a static input, but as a function of the evolving Subhalo Mass Function (SHMF). The AI must be able to calculate the instantaneous ratio of Energy DF to Energy Cool across redshifts (z=0 to z=15).

What the Improved AI System Can Do: This system can predict when a halo will transition from being cooling-dominated (allowing star formation) to heating-dominated (quenching). It moves beyond simple quenching/not quenching classifications and provides a quantitative measure of the energy balance threshold for galaxy quenching, distinguishing between gravitational heating and DF heating.

Improvement: Integrate a Monte Carlo sampling engine using the inverse transform method (as described in Section 2.3) to quantify the variance in DF heating rates (sigma/sigma P). The AI must propagate this inherent Poissonian and non-Poissonian scatter through its predictive models.

What the Improved AI System Can Do: Instead of providing a single best-fit prediction for gas heating, the system outputs a probability distribution (e.g., 16–84% or 2 sigma range). This allows researchers to understand that even if the average DF heating is negligible, there is a non-zero probability of rare, high- DF subhalos causing rapid quenching. This significantly improves the reliability of its predictions in low-resolution simulations.

Improvement: Develop a conditional logic module that links the local Lyman-Werner (LW) background intensity (J 21) to the critical molecular hydrogen fraction (f crit). The AI must dynamically adjust f crit based on both environmental input and internal dynamics.

What the Improved AI System Can Do: This system can accurately predict the conditions required for Pop III star formation (or failure thereof) in minihalos. It determines if a halo's internal H 2 fraction will remain above f crit (allowing collapse) or be suppressed by external radiation, providing a robust mechanism to identify potential Direct Collapse Black Hole (DCBH) seeds versus Pop III star formation pathways.

Improvement: Incorporate the nonlinear drag force correction factor (chi) and the Mach number dependence (h I/M) into the core DF calculation (Ostriker's model). The AI must dynamically adjust h I/M based on the host halo's virial temperature (T vir) and its subhalo velocity distribution.

What the Improved AI System Can Do: The system can provide a more accurate, physically grounded estimate of DF heating that accounts for how the gas interacts with the perturber (e.g., detached bow shock formation), rather than relying on linear approximations. This ensures that predictions are physically consistent with high-velocity, supersonic subhalo encounters in X-ray clusters or early universe halos.

Improvement: Utilize an automated optimization routine to compare results derived from different gas density profiles (e.g., core alpha=0 vs. cusp alpha=1) and the Shifting Subhalo Mass Function (SHMF) against observational constraints.

What the Improved AI System Can Do: The system can autonomously determine which physical model—a highly concentrated NFW profile or a shallower core profile—is most consistent with observed cooling rates at a given redshift, providing an objective assessment of the uncertainty introduced by theoretical assumptions about gas density profiles.

Abstract

We present a model for gas heating, driven by dynamical friction from orbiting subhalos within dark matter halos. Using data from the TNG50 simulation, we derive the subhalo mass function and calculate the dynamical friction heating rate for a wide range of halo masses and redshifts from z = 15 to 0. Our results show that, by converting gravitational potential energy into thermal energy, dynamical friction is an important mechanism for galaxy quenching in massive halos at low redshifts, consistent with previous studies. Additionally, we find that in the early universe at z about 15, heating rates can be comparable to the molecular hydrogen cooling rates in metal-free minihalos. This can suppress gas cooling and fragmentation and does increase the critical molecular fraction for Pop III star formation by up to one order of magnitude, thereby making Pop III star formation more difficult. In combination with the Lyman-Werner background, the dynamical friction heating mechanism favors the formation of direct-collapse black hole (DCBH) seeds in atomic cooling halos, even when the average H 2 fraction is about 10-5 during the minihalo progenitor phase. Dynamical friction heating at a fixed host halo mass can vary by two orders of magnitude due to the scatter in the number of subhalos. To capture dynamical friction heating in simulations, it is necessary to resolve subhalos with a subhalo to host halo mass ratio psi 0.05.

Sources

Related papers