Pith. sign in

REVIEW 4 major objections 6 minor 1 cited by

Investigating the vertical distribution of the disk as a function of radial action: Results from simulations

T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A simple N-body disk reproduces the observed Milky Way relation between radial action and scale height, $z_0=\sqrt{J_R/a}+b$, matching the thin disk and falling short of the inner thick disk.

desk verdict Useful first simulation test of the J23 radial-action/scale-height relation, but the central 'reproduction' claim is visual, not quantitative. read the letter →

arxiv 2411.12432 v1 pith:AZNK6WHV submitted 2024-11-19 astro-ph.GA

Yunpeng Jia , Chengqun Yang , Yuqin Chen , Cuihua Du , Gang Zhao This is my paper · ORCID
classification astro-ph.GA
keywords galacticdiskheatingradialactionscaleheightN-bodysimulationsevolutionMilkyWaythinthickgiantmolecularclouds
topics Dark Matter
open problems Dark Matter
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper sets out to verify, in N-body simulations, an observed correlation between two ways a galactic disk can be heated: radially, measured by the radial action $J_R$, and vertically, measured by the scale height $z_0$. The authors show that a disk embedded in a fixed dark-matter halo produces the same functional relation $z_0 = \sqrt{J_R/a} + b$ that was measured in Milky Way data. The simulated relation lines up with the Galactic thin disk but sits below the inner thick disk, which they read as evidence that the inner thick disk needs an additional heating source or was born hot. The results also show that radial and vertical heating move together particle by particle, and that adding massive, long-lived particles changes the heating rate but not the underlying functional form. If correct, the thin disk's radial-vertical heating correlation is a natural outcome of disk dynamics rather than the signature of one specific heating agent.

What carries the argument

The central object is the relation $z_0 = \sqrt{J_R/a} + b$, where $J_R$ is the radial action, an adiabatic invariant measuring how eccentric a star's orbit is, and $z_0$ is the vertical scale height of the disk computed in bins of $J_R$. The argument is carried by GADGET-4 N-body simulations of a disk in a static Hernquist halo, with radial actions computed by AGAMA and scale heights estimated by maximum likelihood from the vertical distribution $\rho(z)\propto\exp(-z/z_0)$. Early non-axisymmetric structures heat the disk and leave it on this functional relation; in a second simulation, massive particles placed on spiral-arm-like orbits test whether giant-molecular-cloud-like perturbers alter the relation.

What would settle it

Run the same disk in a live, particle-resolved dark-matter halo of the same mass and concentration, or with a different halo scale length, and check whether $z_0$ versus $J_R$ still follows $z_0=\sqrt{J_R/a}+b$ at 4 Gyr and still matches the thin disk; a clear departure would falsify the claim that the functional form is a robust outcome of the model.

Watch

Extended reading notes

Core claim

Using a static Hernquist dark-matter halo and an exponential stellar disk, the simulations reproduce the functional form $z_0 = \sqrt{J_R/a} + b$ at every snapshot after about 0.4 Gyr, the same form previously reported for the Milky Way thin disk. The simulated curves match the observed thin disk, fall below the observed inner thick disk, and roughly agree with the outer thick disk, which the authors take to mean that the heating processes in the simulations are not enough for the inner thick disk. Mean radial action and scale height rise rapidly in the first roughly one gigayear, when non-axisymmetric irregularities heat the disk, then saturate. Particle tracking shows that a particle that gains radial action oscillates more vertically and migrates outward, while one that loses radial action does the opposite, giving a direct particle-level picture of coupled radial-vertical heating. Adding massive, long-lasting particles representing giant molecular clouds increases the rate at which scale height grows with radial action during the first 2-3 Gyr but leaves the functional form intact, so their effect is mainly vertical.

Load-bearing premise

The paper assumes that one particular initial condition, an exponential disk in a fixed spherical Hernquist halo with the authors' chosen masses and scale lengths, heats like the real Milky Way; if a different halo or disk setup changes the shape or normalization of the $z_0$-$J_R$ curve, the comparison with the observed thin and thick disks would not stand.

Editorial extensions

If this is right

  • Because the functional form appears in a disk with no live halo and no giant molecular clouds, the radial-vertical heating correlation is a generic consequence of the disk's own secular evolution, not a fingerprint of any one perturber.
  • The match with the Galactic thin disk supports the idea that the observed thin-disk relation can be produced by internal heating processes of the kind present in the simulations.
  • The shortfall relative to the inner thick disk implies that reproducing the inner thick disk requires either additional heating (for example from mergers) or a population that was born hot.
  • The early rapid rise and later saturation of both mean radial action and scale height imply that most disk heating in this model happens in the first billion years, after which the disk settles.
  • Massive, long-lived particles accelerate vertical heating early on without changing the shape of the $z_0$-$J_R$ relation, so their observable signature is a temporary steepening of the curve.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A testable extension the authors do not run is to vary the halo concentration or replace the static halo with a live one and check whether $z_0 = \sqrt{J_R/a} + b$ still holds; if it does, the functional form is robust and only the parameters $a$ and $b$ shift with total heating.
  • The tracked particles' coupling of radial-action gain with outward migration implies a link between heating and radial migration that the paper does not develop; this could connect the $z_0$-$J_R$ relation to metallicity gradients and age distributions in the disk.
  • If the inner thick disk indeed needs extra heating, a minor-merger run on the same initial conditions should raise the high-$J_R$ end of the curve toward the observed inner thick disk without changing the functional form.
  • One could test whether the same law holds in Milky-Way-mass galaxies that form in full cosmological simulations, a step the paper does not take but that follows naturally from the claim that the correlation emerges from generic disk dynamics.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. This paper uses GADGET-4 N-body simulations of a stellar disk embedded in a static Hernquist dark-matter halo to study the relation between radial action JR and vertical scale height z0. Two setups are considered: a disk without and one with 1000 long-lived massive particles placed on spiral arms as a proxy for giant molecular clouds. The authors report that, after 0.4 Gyr, the simulated JR–z0 relation can be described by the functional form z0 = sqrt(JR/a)+b proposed by Jia et al. (2023, J23), that the relation matches the observed Milky Way thin disk, lies below the inner thick disk, and roughly matches the outer thick disk, implying additional heating is needed for the inner thick disk. They also report rapid early heating, saturation of the relation, and a particle-level correlation between radial heating and vertical heating accompanied by radial migration.

Significance. If quantitatively established, this would be a useful result: it suggests that a simple disk in a static halo generically produces a correlated radial–vertical heating relation with the same functional form as observed, and it would constrain the role of massive perturbers. The paper is transparent in using publicly available tools (GADGET-4, AGAMA, galstep) and shows several example vertical distribution fits. However, because the central functional-form claim is supported only by visual comparison without fitted parameters or residuals, and because the comparison to observations is only qualitative, the significance is currently conditional on additional quantitative analysis.

major comments (4)
  1. [Section 3, paragraph after Fig. 3] The functional-form claim is not quantitatively tested. The text states that the authors have "chosen not to include the best-fit lines" and that "the perfect agreement between the data points of the thin disk observed by J23 and the GMC simulation results at 0.4 Gyr provides compelling evidence." A two-parameter function z0 = sqrt(JR/a)+b is flexible; asserting it on visual grounds is not a test. Please provide, for a representative set of snapshots (e.g., 0.4, 1, 4, and 9 Gyr for both simulations), the best-fit values of a and b with their uncertainties, the residuals or a goodness-of-fit statistic weighted by the reported scale-height uncertainties (0.01–0.03 kpc), and a comparison between the fitted parameters and the J23 thin-disk parameters.
  2. [Section 4, paragraphs on the inner thick disk] The conclusion that additional heating mechanisms are needed for the inner thick disk depends on the assumed correspondence between a specific simulation snapshot and the Milky Way observations. Since no fitted parameters are given, the reader cannot assess whether the simulation's relation at, say, 4 Gyr is statistically consistent with the J23 thin disk and how far it lies from the inner thick disk. Please state which snapshot(s) are used for this comparison, quote the fitted a and b values, and provide a quantitative measure of the offset (e.g., chi-squared or the difference in z0 over the observed JR range) with uncertainties. Otherwise the claim that the inner thick disk requires extra heating is not supported by the data presented.
  3. [Section 2 and final paragraph of Section 4] The acknowledged absence of an in-depth analysis of how the findings depend on the chosen initial conditions is load-bearing. The static Hernquist halo with M=1e12 Msun and a=47 kpc and the disk parameters are taken from Ruggiero & Lima Neto (2017) without a demonstrated connection to the Milky Way. The abstract and conclusions state a match to the Galactic thin disk and a shortfall for the inner thick disk; both statements are normalized to the simulation curves. If a live halo, a different halo concentration, or a different disk mass changes the normalization or the shape of the JR–z0 relation, these astrophysical conclusions would change. The paper should either present tests over a modest grid of initial conditions or explicitly frame the conclusions as valid only for the adopted ICs, with the Milky Way comparison removed from the abstract.
  4. [Section 3, first paragraph after Fig. 3] The exclusion of all snapshots before 0.4 Gyr requires justification. The authors state that the without-GMC simulation does not show monotonic growth before 0.4 Gyr and that AGAMA failed for the GMC simulation. Because the 0.4 Gyr snapshot is the one used for the "perfect agreement" with the thin disk, and because the disk heats rapidly in the first Gyr (Fig. 4), the claim would be strengthened by showing the early-time behavior or explaining why these snapshots are unphysical. As written, the selection of the starting time makes the comparison to observations partially a matter of choosing a convenient epoch.
minor comments (6)
  1. [Section 2, paragraph on massive particles] The phrase "the massive particles are been placed" is ungrammatical; it should read "were placed."
  2. [Section 2, Eq. (2)] The density profile equation is not typeset clearly; specify the argument of the sech factor as z/(2z0) and add the missing parentheses around the exponential factor.
  3. [Fig. 3] The text states that scale-height uncertainties are 0.01–0.02 kpc (0.03 kpc at high JR), but no error bars are plotted for the simulation points; please add them or include a representative error bar in the figure.
  4. [Fig. 3] The legend in the upper panel repeats "APOGEE+Gaia:thin disk" for two entries, and the lower panel repeats the inner/outer thick-disk labels; simplify the legend so that each observational curve is uniquely labeled.
  5. [Section 4] The notation "√aR(R)R" is confusing; define a_R(R) as the radial acceleration and use parentheses consistently throughout the text.
  6. [Abstract and Section 1] The abstract says "Previous research has established a relationship..."; since the paper builds directly on J23, citing J23 explicitly in the abstract would clarify the lineage of the functional form.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the simulation's JR–z0 relation is an independent N-body output, not a fitted renaming of the J23 input; the J23 citation supplies external APOGEE/Gaia data and the functional form under test, not a load-bearing self-citation chain.

full rationale

The paper's central claim is that N-body simulations of a disk embedded in a static Hernquist halo reproduce the empirical functional form z0 = sqrt(JR/a) + b reported in Jia et al. (2023, J23). That claim is a simulation output, not an input: scale heights are measured in radial-action bins via maximum likelihood for each snapshot, and a and b are free parameters that the paper explicitly does not fit or use to force agreement ('we have chosen not to include the best-fit lines for this functional forms for clarity, the perfect agreement between the data points of the thin disk observed by J23 ... and the GMC simulation results at 0.4 Gyr provides compelling evidence'). No equation in the paper converts the initial conditions into the target relation by construction, and no fitted parameter is renamed as a prediction. The J23 citation is heavily used, and the current authors overlap with J23, but J23 supplies an external, APOGEE/Gaia-based observational comparison and the functional form under test; this is real evidence rather than a self-citation chain. The more serious issues—no reported fitted a and b values, no residuals, and possible dependence on the initial conditions (acknowledged in Section 4: 'we acknowledge the absence of an in-depth analysis of how our findings are shaped by the selection of ICs')—are correctness or validation concerns, not circularity. Therefore no circular step is exhibited.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The central comparison rests on the J23 functional form with two free parameters fitted to each snapshot, plus a set of modeling choices (static halo, fixed ICs, action computation, density profile) that are not varied or directly validated. The massive-particle simulation also assumes a simplified GMC proxy.

free parameters (3)
  • a (in z0 = sqrt(JR/a) + b) = not reported (free per snapshot)
    The functional form has two free parameters fitted to each snapshot; no best-fit values are given, so the 'agreement' with the functional form cannot be quantitatively assessed.
  • b (in z0 = sqrt(JR/a) + b) = not reported (free per snapshot)
    Same as a; the intercept parameter is adjusted per snapshot.
  • J-hat_R (pseudo-isothermal scale for radial action distribution) = matches mean radial action of each snapshot
    Used to characterize the disk temperature; the paper verifies the best-fit matches the mean, so it is not independent.
assumptions (6)
  • standard math Epicycle approximation and isothermal-sheet relation sigma_z^2 = 8*pi*G*rho0*z0^2 are used to connect the observed z0 vs JR relation to velocity dispersions (Section 1).
    The paper invokes these to argue that a constant value of a and b implies a linear relation between sigma_R and sigma_z.
  • domain assumption A static Hernquist dark matter halo with M=1e12 Msun and a=47 kpc adequately represents the Galaxy's halo for heating purposes (Section 2).
    A live halo is not simulated; the authors choose a fixed potential to avoid particle noise, but this removes dynamical friction and halo shot noise that could affect vertical heating.
  • domain assumption The disk's own non-axisymmetric irregularities (spiral-like features) are the sole radial and vertical heating source in the without-GMC simulation (Section 4).
    The paper attributes early heating to these irregularities without quantifying their amplitude or comparing with other heating agents.
  • ad hoc to paper Massive, long-lived particles of 1e6 Msun placed on circular orbits along spiral arms are a valid proxy for giant molecular clouds (Section 2).
    The paper acknowledges GMCs are short-lived, so this proxy is a simplification; results from the GMC simulation are interpreted as qualitative.
  • domain assumption The vertical density profile rho(z) proportional to exp(-z/z0) is appropriate for measuring scale heights in all radial action bins and at all snapshots (Section 3).
    The paper verifies this visually for one example snapshot (4 Gyr) but claims consistency across all snapshots.
  • domain assumption Radial actions computed with AGAMA from simulation snapshots are reliable, despite the potential being non-axisymmetric and time-varying during the early epoch (Section 3).
    Actions are exact adiabatic invariants only for slowly varying potentials; during the first 0.4-1 Gyr the disk has strong irregularities, and AGAMA's axisymmetric approximation may be poor. The paper does not test this.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Investigating the vertical distribution of the disk as a function of radial action: Results from simulations." pith.science (2026). https://pith.science/paper/AZNK6WHV

@misc{pith2026241112432,
  author       = {Pith},
  title        = {Pith review of: Investigating the vertical distribution of the disk as a function of radial action: Results from simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AZNK6WHV}},
  note         = {Machine review of arXiv:2411.12432}
}
read the original abstract

Previous research has established a relationship between radial action and scale height in Galactic disks, unveiling a correlation between radial and vertical heating. This finding poses a challenge to our existing comprehension of heating theories and consequently encodes crucial insights into the formation and heating history of Galactic disks. In this study, we perform N-body simulations with the aim of verifying the existence of this correlation between radial action and scale height, thereby enhancing our comprehension of the heating history of Galactic disks. We find that the relationship between radial action and scale height in our simulations can be described by the same functional form observed in previous work. Furthermore, the relationships derived from our simulations align well with those of the Galactic thin disk. However, they do not coincide with the inner thick disk but exhibit a rough correspondence with the outer thick disk, suggesting the possibility that additional heating mechanisms may be required to explain the inner thick disk. We also find that the mean radial action and scale height undergo rapid increases during the initial stages of the simulation, yet remain relatively unchanged as the disk evolves further. By tracing example particles, we uncover a correlation between radial and vertical heating in our simulation: as a particle in the disk gains or loses radial action, its vertical motion tends to oscillate on a more or less extended orbit, accompanied by a tendency to migrate outward or inward, respectively. The massive, long-lasting particles in our simulation contribute to disk heating by solely enhancing the rate of increase in scale height with radial action, while maintaining the functional form that describes the relationship between these two variables.

Figures

Figures reproduced from arXiv: 2411.12432 by the authors.

Figure 1
Figure 1. Snapshots of the disk from the face-on view. The upper left panel specifically depicts the distribution of massive, long-lasting particles within the IC (indicated by red dots), arranged in distinct spiral-like patterns. The snapshots presented in the other panels are derived from the “without-GMC” simulation (a simulation that does not include massive, long-lasting particles). for the GMC simulations in [PITH_FULL… view at source ↗
Figure 2
Figure 2. Vertical distributions of disk particles in GMC simulations at 4 Gyr. The solid lines in each panel denote the best-fit density profiles, ρ(z) ∝ exp(−z/z0). The best-fit scale heights, z0, are listed in the legends. 4. Discussion and conclusions By simulating a disk in a fixed spherical dark halo potential (the without-GMC simulation), we have successfully replicated the functional form that has been reported in J23… view at source ↗
Figure 3
Figure 3. Relationships between radial action and scale height for GMC and without-GMC simulations. The lines with dots show the results of the GMC simulation, while the triangles show the results of the without-GMC simulation. To more clearly demonstrate that the relationships between radial action and scale height in the GMC and without-GMC simulations are nearly identical after 4 Gyr, the upper panel exclusively displays r… view at source ↗
Figures from the paper (1 more)
Figure 5
Figure 5. Figure 5: Tracing a population of disk particles with radial action within the range of 40 to 60 kpc km s−1 in the without-GMC simulation from 0.5 Gyr to 2 Gyr. A group of particles with JR < 20 kpc km s−1 indi￾cates a loss of radial action, whereas a group of particles with JR …

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Vertical Structure and Dynamics of a Galactic Disk

    astro-ph.GA 2025-07 conditional novelty 1.0 of 10

    A review of a multi-component disk plus halo model, arguing that gas and dark matter vertically confine the stellar disk, producing steeper-than-sech^2 profiles and flaring.

Reference graph

Works this paper leans on

30 extracted references · 15 canonical work pages · cited by 1 Pith paper

  1. [1]

    2016, MNRAS, 462, 1697

    Aumer, M., Binney, J., & Schönrich, R. 2016, MNRAS, 462, 1697

    2016

  2. [2]

    & Woltjer, L

    Barbanis, B. & Woltjer, L. 1967, ApJ, 150, 461

    1967

  3. [3]

    C., & Creze, M

    Bienayme, O., Robin, A. C., & Creze, M. 1987, A&A, 180, 94

    work page 1987

  4. [4]

    & Tremaine, S

    Binney, J. & Tremaine, S. 2008, Galactic Dynamics: Second Edition (Princeton University Press)

    2008

  5. [5]

    G., & Martig, M

    Bournaud, F., Elmegreen, B. G., & Martig, M. 2009, ApJ, 707, L1

    2009

  6. [6]

    B., Gibson, B

    Brook, C. B., Gibson, B. K., Martel, H., & Kawata, D. 2005, ApJ, 630, 298

    work page 2005

  7. [7]

    B., Kawata, D., Gibson, B

    Brook, C. B., Kawata, D., Gibson, B. K., & Freeman, K. C. 2004, ApJ, 612, 894

    work page 2004

  8. [8]

    Carlberg, R. G. & Sellwood, J. A. 1985, ApJ, 292, 79

    work page 1985

Show all 30 references
  1. [9]

    2011, A&A, 530, A138

    Casagrande, L., Schönrich, R., Asplund, M., et al. 2011, A&A, 530, A138

    2011

  2. [10]

    N., Frankel, N., Zari, E., Xiang, M., & Rix, H.-W

    Garzon, D. N., Frankel, N., Zari, E., Xiang, M., & Rix, H.-W. 2024, ApJ, 972, 155

    2024

  3. [11]

    Grand, R. J. J., Springel, V ., Gómez, F. A., et al. 2016, MNRAS, 459, 199 Hänninen, J. & Flynn, C. 2002, MNRAS, 337, 731

    2016

  4. [12]

    1990, ApJ, 356, 359

    Hernquist, L. 1990, ApJ, 356, 359

    1990

  5. [13]

    2009, A&A, 501, 941

    Holmberg, J., Nordström, B., & Andersen, J. 2009, A&A, 501, 941

    2009

  6. [14]

    Hou, L. G. & Han, J. L. 2014, A&A, 569, A125

    2014

  7. [15]

    L., Brook, C

    House, E. L., Brook, C. B., Gibson, B. K., et al. 2011, MNRAS, 415, 2652

    2011

  8. [16]

    2023, A&A, 669, A107

    Jia, Y ., Chen, Y ., Du, C., & Zhao, G. 2023, A&A, 669, A107

    2023

  9. [17]

    & Davies, M

    Kokaia, G. & Davies, M. B. 2019, MNRAS, 489, 5165

    2019

  10. [18]

    Lacey, C. G. 1984, MNRAS, 208, 687

    1984

  11. [19]

    D., Fall, S

    Ludlow, A. D., Fall, S. M., Schaye, J., & Obreschkow, D. 2021, MNRAS, 508, 5114

    2021

  12. [20]

    & Kalnajs, A

    Lynden-Bell, D. & Kalnajs, A. J. 1972, MNRAS, 157, 1

    1972

  13. [21]

    T., Bovy, J., Leung, H

    Mackereth, J. T., Bovy, J., Leung, H. W., et al. 2019, MNRAS, 489, 176 Nordström, B., Mayor, M., Andersen, J., et al. 2004, A&A, 418, 989

    2019

  14. [22]

    J., Hernquist, L., & Fullagar, D

    Quinn, P. J., Hernquist, L., & Fullagar, D. P. 1993, ApJ, 403, 74

    1993

  15. [23]

    & Lima Neto, G

    Ruggiero, R. & Lima Neto, G. B. 2017, MNRAS, 468, 4107

    2017

  16. [24]

    Saha, K., Tseng, Y .-H., & Taam, R. E. 2010, ApJ, 721, 1878

    2010

  17. [25]

    Sanders, J. L. & Das, P. 2018, MNRAS, 481, 4093

    2018

  18. [26]

    & Schwarzschild, M

    Spitzer, Lyman, J. & Schwarzschild, M. 1953, ApJ, 118, 106

    1953

  19. [27]

    2021, MNRAS, 506, 2871

    Springel, V ., Pakmor, R., Zier, O., & Reinecke, M. 2021, MNRAS, 506, 2871

    2021

  20. [28]

    & Rix, H.-W

    Ting, Y .-S. & Rix, H.-W. 2019, The Astrophysical Journal, 878, 21

    2019

  21. [29]

    2019, MNRAS, 482, 1525 Villalobos, Á

    Vasiliev, E. 2019, MNRAS, 482, 1525 Villalobos, Á. & Helmi, A. 2008, MNRAS, 391, 1806

    2019

  22. [30]

    T., Su, G.-M., et al

    Yang, H.-Y ., Chiang, B. T., Su, G.-M., et al. 2024, MNRAS, 530, 129 Article number, page 6 of 6

    2024

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.