REVIEW 5 minor 1 cited by
The PGPUC horizontal branch evolutionary tracks
T0 review · 0 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read The PGPUC database now holds a fine grid of over 19,000 horizontal-branch evolutionary tracks, plus derived zero-, middle-, and terminal-age loci, enabling precise interpolation and comparison with observations.
desk verdict A solid, useful database paper: 19,000+ PGPUC HB tracks with derived loci, public download, and an honest BaSTI comparison; the main limits are closed code and inherited older input physics, neither of which sinks the resource. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing quantity is the helium core mass at the tip of the red giant branch, $M_{\rm cHe}$: for a given envelope composition, practically the entire HB evolution depends on this one number, and the paper tabulates it for every combination of composition and progenitor mass. The machinery that produces those core masses is the PGPUC stellar evolution code, a code for low-mass stars whose input physics was deliberately kept unchanged — most importantly the semiconvection treatment at the edge of the convective helium core (following Sweigart & Demarque 1972) and the suppression of 'breathing pulses' (following Sweigart 1971), together with the adopted opacities, equation of state, and nuclear reaction rates. The third ingredient is the progenitor mass $M_{\rm MS}$, which fixes $M_{\rm cHe}$ for a given composition and thereby encodes the age of the stellar population in the shape and position of the HB.
What would settle it
Recompute a slice of the grid with a modern stellar-evolution code using updated opacities, equation of state, and a different semiconvection treatment, and check whether the zero-age horizontal branch luminosities and helium core masses shift by more than the grid's own interpolation precision; alternatively, count red HB stars and RR Lyrae variables in ultra-metal-poor systems, where the grid predicts that only stars near core-helium exhaustion can appear.
Extended reading notes
Core claim
The paper's discovery, stated on its own terms, is that one internally consistent evolutionary grid can now carry low-mass stars from the main sequence through the red giant tip and across the entire HB, with the HB section dense enough for precise interpolation: 19,441 tracks over 11 metallicities, 7 helium abundances, two $\alpha$-enhancements, and three progenitor masses, with HB masses stepped by $0.01\,M_\odot$ down to $0.55\,M_\odot$ and by $0.005\,M_\odot$ below that, down to a minimum mass that sits only about $0.0002\,M_\odot$ above the helium core mass. Using these tracks, the authors derive zero-age (ZAHB), middle-age (MAHB), and terminal-age (TAHB) loci for every composition and quantify their dependence on the helium core mass $M_{\rm cHe}$ inherited from the RGB tip: $M_{\rm cHe}$ decreases with increasing helium abundance and progenitor mass and varies weakly with metallicity, the hottest reachable HB model moves to lower masses as $Z$ and $Y$ rise, and higher progenitor masses yield redder, fainter HB models. A central systematic result is that at the lowest metallicities the ZAHB never becomes cool enough to enter the RR Lyrae instability strip, implying that extremely metal-poor red HB stars and RR Lyrae variables can only exist close to the point of core-helium exhaustion and should be very rare — a conclusion the paper notes is consistent with the ultra-metal-poor RR Lyrae stars observed so far.
Load-bearing premise
Every one of the 19,000-plus tracks inherits the code's deliberately retained older input physics — most importantly the semiconvection treatment at the helium-core edge and the suppression of breathing pulses, together with the adopted opacities and equation of state — so any systematic error in those prescriptions is baked uniformly into the whole grid.
Editorial extensions
If this is right
- Any composition inside the grid can be interpolated to the full HB evolution — from zero age to terminal age — and converted into roughly 50 photometric filter systems, so observed star clusters can be matched to models directly in color–magnitude space.
- HB masses read off the loci without choosing the right progenitor mass can be systematically wrong, with effective-temperature differences reaching about 6,000 K at the hot end, so the progenitor (hence age) must be matched to the population being studied.
- At the lowest metallicities the grid predicts that red HB stars and RR Lyrae variables should be very rare, and any that exist should be close to exhausting the helium in their cores — a prediction the paper finds consistent with the known ultra-metal-poor RR Lyrae stars.
- The luminosity gap between the zero-age and middle-age loci depends strongly on helium abundance and metallicity, so distance, helium, and age estimates based on HB fitting will be biased if the wrong composition is adopted.
Reading between the lines
- An inter-code comparison at fixed composition could separate the grid's interpolation smoothness from genuine physical accuracy, since the authors deliberately kept older physics for consistency with the existing database; if modern codes shift ZAHB luminosities by more than the mass stepping does, interpolated values could be precisely wrong rather than approximately right.
- The strong sensitivity of track temperatures to progenitor mass at the hot end suggests extreme HB stars could serve as an age diagnostic for old populations, independent of main-sequence turnoff fitting — a direction the paper does not develop.
- The rarity prediction becomes a statistical test: as samples of ultra-metal-poor RR Lyrae grow, their positions along the HB should cluster at the terminal-age rather than the zero-age locus, and a significant population of zero-age stars would challenge the grid's semiconvection-dependent core masses.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a new grid of 19,441 horizontal-branch (HB) evolutionary tracks computed with the PGPUC stellar evolution code. The grid spans 11 metallicities, 7 initial helium abundances, 2 alpha-element enhancement values ([α/Fe]=0.0 and 0.3), and 3 progenitor masses (0.7, 0.8, and 0.9 M_sun). The tracks are made publicly available through the PGPUC Online database, and the authors use them to construct zero-age, middle-age, and terminal-age HB loci. The paper describes the input physics, the grid structure, and a comparison with BaSTI models, and it discusses how HB morphology depends on composition and progenitor mass. The scientific content is primarily a database resource rather than a new physical result.
Significance. If the grid is reliable, it is a valuable public resource for HB studies, particularly because it systematically includes progenitor mass as a parameter and extends to very low metallicities and high helium abundances. The grid's mass resolution (0.005-0.01 M_sun) is fine, and the integration into an online interpolation tool makes it a practical tool for comparing theory with observations. The authors are transparent about the older input physics (Sweigart & Demarque 1972 semiconvection, breathing-pulse suppression, I99/F05 opacities, FreeEOS) and justify retaining it for consistency with the existing PGPUC database. Appendix A honestly addresses the known faintness of PGPUC RGB tips and corrects a misattribution about electron screening. The comparison with BaSTI, while not identical in input physics or alpha-enhancement, gives a sanity check that the loci are broadly consistent. No load-bearing technical errors are evident; the main weaknesses are presentation and verification details.
minor comments (5)
- [Abstract / Sect. 3] The claim of 'precise interpolation' should be substantiated: the grid spacing is described, but no interpolation test or error estimate is provided, and the interpolation method used by the PGPUC Online database is not described or referenced.
- [Sect. 4 / Fig. 3] The BaSTI comparison uses [α/Fe]=0.4 while the PGPUC grid uses [α/Fe]=0.3, and the Z,Y values shown in the figure panels (e.g., Z=0.00002, Y=0.247) are not PGPUC grid nodes; the text should state explicitly that the PGPUC loci were obtained by interpolation and should mention the alpha-enhancement mismatch when interpreting the differences.
- [Sect. 3] The total of 19,441 tracks is not derived from the stated grid parameters; a short formula or a table giving the number of tracks per (Z, Y, [α/Fe], M_MS) combination would allow readers to verify the count.
- [Sect. 2.2 / 3] The database URLs should be accompanied by a persistent identifier (e.g., a DOI) and a version number so that the specific track set described in this paper can be cited and retrieved unambiguously.
- [Sects. 2.2 and 4] The statement in Sect. 2.2 that higher M_MS progenitors 'tend to appear redder and more luminous' and the statement in Sect. 4 that higher M_MS 'results in a lower luminosity at a given Teff' should be reconciled by specifying the relevant mass range, to avoid apparent contradiction.
Circularity Check
No circularity found: the HB tracks and loci are direct computational products of the PGPUC code, not fitted to observations, and the BaSTI comparison provides an external benchmark.
full rationale
The paper's central claim is a database-resource claim: it presents 19,441 newly computed HB evolutionary tracks covering specified ranges of Z, Y, [alpha/Fe], and M_MS, plus derived ZAHB, MAHB, and TAHB loci. These quantities are direct outputs of stellar evolution calculations with the PGPUC code; they are not fitted to any subset of data and then re-predicted. The MAHB and TAHB loci are defined as time-percentage loci from the tracks (following Catelan et al. 2009), so they are constructions from the tracks themselves rather than independent predictions that reduce to inputs. The paper explicitly compares the PGPUC ZAHB and TAHB loci with BaSTI models (Fig. 3) and reports small luminosity differences, which is an external, independent benchmark rather than a self-referential validation. Self-citations to Valcarce et al. (2012) for the code description and to Catelan et al. (2009) for locus definitions are not load-bearing in a circular sense: the code is a tool whose outputs are the content of the paper, and the locus definitions are standard descriptive conventions. The retention of older input physics is disclosed as a consistency choice, not hidden behind a citation. No equation or fitted parameter is renamed as a prediction, and no uniqueness claim is imported from the authors' prior work. The paper is self-contained in its derivation chain: the tracks are computed, not inferred from the conclusions they are used to illustrate.
Assumptions & free parameters
free parameters (1)
- RGB mass loss efficiency (eta) =
1.0 (default), 0.0 (alternative)
assumptions (5)
- domain assumption Standard stellar physics inputs: OPAL and Ferguson opacities, Cassisi conductive opacities, Haft neutrino losses, NACRE reaction rates with Kunz 12C(alpha,gamma)16O.
- domain assumption FreeEOS equation of state (Irwin 2007), option EOS4, is accurate for HB conditions.
- domain assumption Semiconvection treatment at the outer edge of the convective core (Sweigart & Demarque 1972) and breathing pulse suppression (Sweigart 1971) correctly model HB core evolution.
- domain assumption Electron screening of nuclear reactions is described by DeWitt et al. (1973) and Graboske et al. (1973), not by Salpeter (1954).
- domain assumption Mass loss is turned off during the HB phase.
Cite this review
Pith. "Pith review of The PGPUC horizontal branch evolutionary tracks." pith.science (2026). https://pith.science/paper/QPNZ4KVX
@misc{pith2026250616562,
author = {Pith},
title = {Pith review of: The PGPUC horizontal branch evolutionary tracks},
year = {2026},
howpublished = {\url{https://pith.science/paper/QPNZ4KVX}},
note = {Machine review of arXiv:2506.16562}
}
read the original abstract
The horizontal branch (HB) phase of stellar evolution plays a critical role in understanding the life cycle of stars, particularly in the context of low-mass stars. However, generating theoretical evolutionary tracks for HB stars is computationally intensive, especially when attempting to match a wide range of observational data. This work presents an extensive grid of HB evolutionary tracks computed using the PGPUC code, covering a broad range of chemical compositions, progenitor masses, and alpha-element distributions. The aim is to provide a robust tool for interpreting HB stellar populations and advancing our understanding of their diverse properties. The evolutionary tracks are made publicly available through the PGPUC Online database for easy access and interpolation. We computed over 19,000 HB evolutionary tracks encompassing a wide range in terms of mass, metallicity, helium abundance, and alpha-element enhancement, along with different progenitor masses. Using these tracks, we calculated zero-age, middle-age, and terminal-age loci. The PGPUC database now includes a fine grid of HB evolutionary tracks, allowing for precise interpolation. Key findings include the dependence of HB morphology on progenitor mass, metallicity, and helium abundance.
Figures
Forward citations
Cited by 1 Pith paper
-
Neighbors, Not Kin: Kinematic evidence for tidal tails from NGC 7492 along the Sagittarius stream
GIRAFFE radial velocities and metallicities confirm NGC 7492 tidal tails extending ≥1.8 deg, kinematically and chemically distinct from the overlapping Sagittarius stream.
Reference graph
Works this paper leans on
- [1]
-
[2]
Bono , G., Caputo , F., Cassisi , S., Castellani , V., & Marconi , M. 1997 a , , 489, 822
work page 1997
-
[3]
Bono , G., Caputo , F., Cassisi , S., Castellani , V., & Marconi , M. 1997 b , , 479, 279
work page 1997
-
[4]
Caputo , A. & Raffelt , G. 2024, in 1st General Meeting and 1st Training School of the COST Action COSMIC WISPers, Vol. 1, 41
work page 2024
- [5]
-
[6]
Carenza , P., Giannotti , M., Isern , J., Mirizzi , A., & Straniero , O. 2025, , 1117, 1
work page 2025
-
[7]
& VandenBerg , D
Casagrande , L. & VandenBerg , D. A. 2018, , 475, 5023
2018
- [8]
Show all 74 references
-
[9]
1996, , 459, 298
Cassisi , S., Castellani , V., & Tornambe , A. 1996, , 459, 298
1996
-
[10]
Y., Pietrinferni , A., Catelan , M., & Salaris , M
Cassisi , S., Potekhin , A. Y., Pietrinferni , A., Catelan , M., & Salaris , M. 2007, , 661, 1094
2007
-
[11]
Cassisi , S., Salaris , M., & Irwin , A. W. 2003, , 588, 862
2003
-
[12]
1994, , 282, 771
Castellani , M., Castellani , V., Pulone , L., & Tornambe , A. 1994, , 282, 771
1994
-
[13]
2009, , 320, 261
Catelan , M. 2009, , 320, 261
2009
-
[14]
V., Valcarce , A
Catelan , M., Grundahl , F., Sweigart , A. V., Valcarce , A. A. R., & Cort \'e s , C. 2009, , 695, L97
2009
-
[15]
& Smith , H
Catelan , M. & Smith , H. A. 2015, Pulsating Stars (Wiley-VCH, Weinheim)
2015
-
[16]
2019, , 632, A105
Chen , Y., Girardi , L., Fu , X., et al. 2019, , 632, A105
2019
-
[17]
& Mengel , J
Demarque , P. & Mengel , J. G. 1972, , 171, 583
1972
-
[18]
E., Graboske , H
DeWitt , H. E., Graboske , H. C., & Cooper , M. S. 1973, , 181, 439
1973
-
[19]
2025, , 694, A158
D'Orazi , V., Braga , V., Bono , G., et al. 2025, , 694, A158
2025
-
[20]
1992, , 80, 701
Dorman , B. 1992, , 80, 701
1992
-
[21]
2007, , 134, 376
Dotter , A., Chaboyer , B., Jevremovi \'c , D., et al. 2007, , 134, 376
2007
-
[22]
2008, , 178, 89
Dotter , A., Chaboyer , B., Jevremovi \'c , D., et al. 2008, , 178, 89
2008
-
[23]
Eggleton , P. P. 1968, , 140, 387
1968
-
[24]
1966, , 144, 978
Faulkner , J. 1966, , 144, 978
1966
-
[25]
& Iben , Jr., I
Faulkner , J. & Iben , Jr., I. 1966, , 144, 995
1966
-
[26]
W., Alexander , D
Ferguson , J. W., Alexander , D. R., Allard , F., et al. 2005, , 623, 585
2005
-
[27]
C., Dewitt , H
Graboske , H. C., Dewitt , H. E., Grossman , A. S., & Cooper , M. S. 1973, , 181, 457
1973
-
[28]
& Sauval , A
Grevesse , N. & Sauval , A. J. 1998, , 85, 161
1998
-
[29]
Gross , P. G. 1973, , 164, 65
1973
-
[30]
1994, , 425, 222
Haft , M., Raffelt , G., & Weiss , A. 1994, , 425, 222
1994
-
[31]
& Schwarzschild , M
H \"a rm , R. & Schwarzschild , M. 1966, , 145, 496
1966
-
[32]
L., Pietrinferni , A., Cassisi , S., et al
Hidalgo , S. L., Pietrinferni , A., Cassisi , S., et al. 2018, , 856, 125
2018
-
[33]
1964, , 140, 1631
Iben , Jr., I. 1964, , 140, 1631
1964
-
[34]
& Faulkner , J
Iben , Jr., I. & Faulkner , J. 1968, , 153, 101
1968
-
[35]
& Rood , R
Iben , Jr., I. & Rood , R. T. 1970, , 161, 587
1970
-
[36]
Iglesias , C. A. & Rogers , F. J. 1996, , 464, 943
1996
-
[37]
Irwin , A. W. 2007, http://freeeos.sourceforge.net
2007
-
[38]
Krishna Swamy , K. S. 1966, , 145, 174
1966
-
[39]
2002, , 567, 643
Kunz , R., Fey , M., Jaeger , M., et al. 2002, , 567, 643
2002
-
[40]
2023, , 677, A86
Latour , M., H \"a mmerich , S., Dorsch , M., et al. 2023, , 677, A86
2023
-
[41]
2004, , 612, 168
Pietrinferni , A., Cassisi , S., Salaris , M., & Castelli , F. 2004, , 612, 168
2004
-
[42]
2013, , 558, A46
Pietrinferni , A., Cassisi , S., Salaris , M., & Hidalgo , S. 2013, , 558, A46
2013
-
[43]
2021, , 908, 102
Pietrinferni , A., Hidalgo , S., Cassisi , S., et al. 2021, , 908, 102
2021
-
[44]
1977, in Saas-Fee Advanced Course 7: Advanced Stages in Stellar Evolution, ed
Renzini , A. 1977, in Saas-Fee Advanced Course 7: Advanced Stages in Stellar Evolution, ed. P. Bouvier & A. Maeder , 151
1977
-
[45]
Rood , R. T. 1973, , 184, 815
1973
-
[46]
& Cassisi , S
Salaris , M. & Cassisi , S. 2005, Evolution of Stars and Stellar Populations
2005
-
[47]
Salpeter , E. E. 1954, Australian Journal of Physics, 7, 373
1954
-
[48]
Schr \"o der , K. P. & Cuntz , M. 2005, , 630, L73
2005
-
[49]
& H \"a rm , R
Schwarzschild , M. & H \"a rm , R. 1965, , 142, 855
1965
-
[50]
2017, , 606, A33
Serenelli , A., Weiss , A., Cassisi , S., Salaris , M., & Pietrinferni , A. 2017, , 606, A33
2017
-
[51]
2020, , 644, A166
Straniero , O., Pallanca , C., Dalessandro , E., et al. 2020, , 644, A166
2020
-
[52]
Sweigart , A. V. 1971, , 168, 79
1971
-
[53]
Sweigart , A. V. 1973, , 24, 459
1973
-
[54]
Sweigart , A. V. 1987, , 65, 95
1987
-
[55]
Sweigart , A. V. 1994, , 426, 612
1994
-
[56]
Sweigart , A. V. 1997, , 474, L23
1997
-
[57]
Sweigart , A. V. & Catelan , M. 1998, , 501, L63
1998
-
[58]
Sweigart , A. V. & Demarque , P. 1972, , 20, 445
1972
-
[59]
Sweigart , A. V. & Gross , P. G. 1974, , 190, 101
1974
-
[60]
Sweigart , A. V. & Gross , P. G. 1976, , 32, 367
1976
-
[61]
Sweigart , A. V. & Gross , P. G. 1978, , 36, 405
1978
-
[62]
Valcarce , A. A. R. 2011, Study of the Helium Enrichment in Globular Clusters (Ph.D. Thesis, Pontificia Universidad Católica de Chile)
2011
-
[63]
Valcarce , A. A. R., Catelan , M., Alonso-Garc \' a , J., Contreras Ramos , R., & Alves , S. 2016, , 589, A126
2016
-
[64]
Valcarce , A. A. R., Catelan , M., Alonso-Garc \' a , J., Cort \'e s , C., & De Medeiros , J. R. 2014, , 782, 85
2014
-
[65]
Valcarce , A. A. R., Catelan , M., & Sweigart , A. V. 2012, , 547, A5
2012
-
[66]
VandenBerg , D. A. 1992, , 391, 685
1992
-
[67]
VandenBerg , D. A. 2024, , 527, 6888
2024
-
[68]
A., Bergbusch , P
VandenBerg , D. A., Bergbusch , P. A., & Dowler , P. D. 2006, , 162, 375
2006
-
[69]
A., Swenson , F
VandenBerg , D. A., Swenson , F. J., Rogers , F. J., Iglesias , C. A., & Alexander , D. R. 2000, , 532, 430
2000
-
[70]
B., et al
Viaux , N., Catelan , M., Stetson , P. B., et al. 2013 a , , 111, 231301
2013
-
[71]
B., et al
Viaux , N., Catelan , M., Stetson , P. B., et al. 2013 b , , 558, A12
2013
-
[72]
2000, , 538, 289
Zoccali , M., Cassisi , S., Bono , G., et al. 2000, , 538, 289
2000
-
[73]
, " * write output.state after.block = add.period write newline
ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sent...
-
[74]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 gl...
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.