arxiv: 2605.04007 · v2 · submitted 2026-05-05 · 🌌 astro-ph.IM · astro-ph.GA· astro-ph.SR

Recognition: no theorem link

PySME v1.0: improved modelling of stellar spectra for survey-scale applications

Mingjie Jian , Nikolai Piskunov , Jeff Valenti , Ella Xi Wang , Brian Thorsbro , Henrik J\"onsson , Ansgar Wehrhahn

Authors on Pith no claims yet

Pith reviewed 2026-05-12 02:04 UTC · model grok-4.3

classification 🌌 astro-ph.IM astro-ph.GAastro-ph.SR

keywords PySMEstellar spectraspectral synthesisabundance analysisline listsNLTEsurvey-scale modelingequation of state

0 comments

The pith

PySME v1.0 trims large line lists via opacity ratios to enable survey-scale stellar spectrum synthesis.

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

The paper presents PySME v1.0 as an updated Python version of Spectroscopy Made Easy that addresses the computational cost of handling very large atomic line lists in stellar abundance work. It describes a revised selection method that uses opacity ratio and line depth to drop negligible lines, combined with dynamic list construction and wavelength-span limits so that only relevant contributions reach the synthesis core. An updated equation-of-state treatment improves hydrogen-line modeling, especially Balmer features, while metal-line results stay close to the prior version. The interface now supports NLTE departure grids for 17 elements and parameter-dependent derived quantities during fitting. A reader would care because modern surveys observe spectra for hundreds of thousands of stars, and earlier codes become impractical without losing either speed or precision.

Core claim

PySME v1.0 introduces a revised line-selection framework based on opacity ratio and line depth, together with dynamic line list construction and control of the effective wavelength span over which each line contributes to the synthetic spectrum. These changes support parallel preprocessing of weak-line selection and reduce the line list passed to the synthesis core. An updated equation-of-state treatment improves the modelling of hydrogen lines, particularly Balmer features, while maintaining close agreement with previous SME results for metal lines. The Python interface is extended to support parameter-dependent derived quantities updated during optimisation, and non-local thermodynamic eqm

What carries the argument

Revised line-selection framework based on opacity ratio and line depth that reduces the line list passed to the synthesis core while preserving accuracy.

If this is right

High-precision stellar abundance analyses become feasible for the millions of spectra expected from current and future large surveys.
Synthetic spectra for metal lines remain in close agreement with earlier SME calculations.
Balmer line profiles are modeled more accurately without separate post-processing steps.
NLTE effects for 17 elements can be included directly in the fitting process.
Parameter-dependent derived quantities can be updated automatically during optimisation loops.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same line-pruning logic could be ported to other synthesis codes that face growing atomic databases.
Real-time quality checks during survey observations might become possible if the speed gain is large enough.
Users could validate the accuracy claim by re-reducing a public archive of benchmark stars with both versions.

Load-bearing premise

That the opacity-ratio and line-depth selection plus the new equation-of-state treatment preserves synthetic accuracy while substantially improving scalability.

What would settle it

Compare synthetic spectra and derived abundances from PySME v1.0 against the original SME on the same set of high-resolution survey spectra and measure whether abundance differences stay within typical uncertainties while wall-clock time drops by a large factor.

Figures

Figures reproduced from arXiv: 2605.04007 by Ansgar Wehrhahn, Brian Thorsbro, Ella Xi Wang, Henrik J\"onsson, Jeff Valenti, Mingjie Jian, Nikolai Piskunov.

**Figure 1.** Figure 1: Periodic table highlighting the chemical elements for which 1D NLTE departure coefficients are available in PySME. both isolated strong lines and clusters of weak lines are retained whenever their combined contribution may become significant within a local wavelength interval. When q = dc, the threshold is more directly interpretable in terms of the impact on the continuum-normalised spectrum, and qw corre… view at source ↗

**Figure 2.** Figure 2: Synthetic spectrum of the example lines listed in view at source ↗

**Figure 3.** Figure 3: Left and middle panels: growth of dc , da and λv for two example Fe lines with metallicity in solar parameters. Right panel: dc vs da for the example Fe lines. reaching approximately 75% around Teff ∼ 3000 K (and even higher for log g ∼ 0), where nearly all lines must be included and synthesis becomes significantly slower. A similar trend is observed at lower metallicities, although the temperature thresh… view at source ↗

**Figure 4.** Figure 4: Flux contribution function for the stronger (6003.011 Å) and weaker (6004.045 Å) Fe i lines, plotted as a function of column mass density ρx. Lower values of ρx correspond to higher layers in the stellar atmosphere, while higher values represent deeper layers. 5.4. Comparison between the old and new SMElib The differences between spectra synthesised with PySME v0.4 (using SMElib v6.0) and v1.0 (using SMEli… view at source ↗

**Figure 5.** Figure 5: Central depths of optical spectral lines for the Sun (top) and Arcturus (bottom). Coloured points mark molecular features, and grey points show all other atomic lines view at source ↗

**Figure 6.** Figure 6: Difference between synthetic spectra computed with negligible lines removed (fs) and those including all lines (fall), shown for solar parameters at infinite resolution (top) and resolution 50,000 (bottom). Horizontal dashed lines indicate the dw thresholds, and fall is plotted in grey as a reference. Amarsi, A. M., Barklem, P. S., Collet, R., Grevesse, N., & Asplund, M. 2019, A&A, 624, A111 Amarsi, A. M.,… view at source ↗

**Figure 7.** Figure 7: Same as view at source ↗

**Figure 8.** Figure 8: The ratio of non-negligible lines across the Kiel diagram for MARCS grid points. Grid points marked with orange circles indicate locations where more than 75% of TiO lines are non-negligible, while orange crosses denote points where all TiO lines are negligible. Blanco-Cuaresma, S., Soubiran, C., Jofré, P., & Heiter, U. 2014, A&A, 566, A98 Buder, S., Sharma, S., Kos, J., et al. 2021, MNRAS, 506, 150 Calisk… view at source ↗

**Figure 9.** Figure 9: Comparison between the observed solar spectrum in the Hα region and synthetic spectra generated using PySME v0.4, PySME v1.0, and Turbospectrum view at source ↗

**Figure 10.** Figure 10: Wall-time comparison of line selection and post-selection synthesis in PySME for synthesis windows of different widths starting at 5880 Å. Left: serial selection time for the internal (default), ALMAX, and dc-based workflows. Right: post-selection synthesis time for PySME v0.4 as well as PySME v1.0 using ALMAX-selected (internal) and dc-selected line lists. Lind, K., Melendez, J., Asplund, M., Collet, R.,… view at source ↗

read the original abstract

Stellar abundance analysis relies on flexible, high-performance spectral synthesis. To meet these needs, we present PySME v1.0, an updated Python implementation of Spectroscopy Made Easy (SME) designed for precise and survey-scale modelling of stellar spectra.A central challenge in SME based synthesis is the efficient treatment of very large line lists, including both the preselection of negligible lines and the subsequent formal synthesis. PySME v1.0 introduces a revised line-selection framework based on opacity ratio and line depth, together with dynamic line list construction and control of the effective wavelength span over which each line contributes to the synthetic spectrum. These workflows support parallel preprocessing of weak-line selection and reduce the line list passed to the synthesis core, thereby improving scalability while preserving synthetic accuracy. PySME v1.0 also incorporates an updated equation-of-state treatment that improves the modelling of hydrogen lines, particularly Balmer features, while maintaining close agreement with previous SME results for metal lines. The Python interface has further been extended to support parameter-dependent derived quantities updated during optimisation, and PySME provides non-local thermodynamic equilibrium (NLTE) departure-coefficient grids for 17 elements. Together, these developments establish PySME v1.0 as a robust and efficient framework for high-precision stellar abundance analyses in large spectroscopic surveys.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

PySME v1.0 updates SME with opacity-based line selection, dynamic lists, wavelength control, and better hydrogen EOS handling to support survey-scale work, but the abstract supplies no benchmarks or residuals to show the accuracy and speed claims actually hold.

read the letter

PySME v1.0 is a practical update to the Spectroscopy Made Easy code, now in Python, aimed at stellar abundance work on large surveys. The concrete changes are a revised line-selection approach that uses opacity ratios and line depths, dynamic construction of the active line list, and explicit control over the wavelength span each line contributes to the spectrum. They also revised the equation-of-state treatment to improve Balmer-line modeling while keeping metal-line results close to earlier SME versions. The Python layer now carries parameter-dependent derived quantities through optimization and ships NLTE departure-coefficient grids for 17 elements. These steps target the real bottleneck of very large line lists in survey pipelines and add some flexibility for fitting workflows. The description of the implementation is straightforward and stays within the existing SME framework, which should make adoption easier for current users. The engineering choices around parallel preprocessing and reduced line lists passed to the synthesis core look like sensible responses to scale problems. The main gap is the complete absence of numbers. The abstract asserts that synthetic accuracy is preserved and scalability improves, yet it contains no timing comparisons, no spectrum residuals against the old code or against observations, no scaling curves with line-list size, and no abundance tests. Without those checks it is impossible to judge whether the new selection rules introduce systematic offsets or whether the speed gain is large enough to matter in practice. The stress-test note correctly flags this as the load-bearing claim that remains untested from the given text. This paper is aimed at people who already run SME or similar synthesis codes inside survey pipelines and need to process thousands of spectra without hitting memory or time walls. A reader building or maintaining such a pipeline would find the line-handling details and the NLTE additions directly relevant. It deserves a serious referee. Tool updates to established codes can be worth publishing when they solve documented bottlenecks, provided the full manuscript includes the missing validation data and confirms code availability. I would send it to review rather than desk-reject, with the expectation that referees will require the benchmarks before acceptance.

Referee Report

1 major / 0 minor

Summary. The manuscript presents PySME v1.0, an updated Python implementation of Spectroscopy Made Easy (SME) for stellar spectral synthesis. It introduces a revised line-selection framework based on opacity ratio and line depth, with dynamic line list construction and wavelength-span control; an updated equation-of-state treatment that improves hydrogen-line (particularly Balmer) modeling; an extended Python interface supporting parameter-dependent derived quantities; and NLTE departure-coefficient grids for 17 elements. The central claim is that these changes improve scalability for survey-scale applications while preserving synthetic accuracy and maintaining close agreement with prior SME results for metal lines.

Significance. If the asserted scalability gains and accuracy preservation are validated, the work would provide a useful practical update to an established spectral synthesis tool, directly addressing the computational demands of large line lists in modern spectroscopic surveys.

major comments (1)

Abstract: the central claim that the revised line-selection framework (opacity ratio + line depth, dynamic construction, wavelength-span control) together with the updated EOS treatment 'improve scalability while preserving synthetic accuracy' and yield 'close agreement with previous SME results for metal lines' is asserted without any quantitative benchmarks, timing deltas, residual statistics, line-by-line comparisons, or validation spectra. Because the headline conclusion that PySME v1.0 constitutes 'a robust and efficient framework for high-precision stellar abundance analyses in large spectroscopic surveys' depends directly on these two improvements holding simultaneously, the absence of supporting data renders the claim untestable from the manuscript as presented.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive review and recommendation. We address the single major comment below.

read point-by-point responses

Referee: [—] Abstract: the central claim that the revised line-selection framework (opacity ratio + line depth, dynamic construction, wavelength-span control) together with the updated EOS treatment 'improve scalability while preserving synthetic accuracy' and yield 'close agreement with previous SME results for metal lines' is asserted without any quantitative benchmarks, timing deltas, residual statistics, line-by-line comparisons, or validation spectra. Because the headline conclusion that PySME v1.0 constitutes 'a robust and efficient framework for high-precision stellar abundance analyses in large spectroscopic surveys' depends directly on these two improvements holding simultaneously, the absence of supporting data renders the claim untestable from the manuscript as presented.

Authors: We agree that the abstract, as a concise summary, would be strengthened by direct reference to quantitative validation results. The manuscript body contains dedicated validation material, including line-list size reductions, timing comparisons, residual statistics between PySME and prior SME versions for metal lines, and example synthetic spectra. In the revised manuscript we will update the abstract to incorporate brief, specific quantitative statements (e.g., typical line-list compression factors and agreement metrics) drawn from those sections, thereby making the central claims directly testable from the abstract itself. revision: yes

Circularity Check

0 steps flagged

No circularity: abstract describes concrete software implementation changes without derivation or self-referential reduction

full rationale

The abstract presents PySME v1.0 as an updated Python implementation of SME, introducing a revised line-selection framework (opacity ratio and line depth, dynamic construction, wavelength-span control) and an updated equation-of-state treatment. These are described as supporting parallel preprocessing, reducing the line list for synthesis, improving scalability, preserving synthetic accuracy, and maintaining agreement with prior SME results for metal lines. No equations, fitted parameters, predictions, or derivation chain are provided that reduce to inputs by construction. The text contains no self-citations used as load-bearing justification for uniqueness or ansatz. This is a standard software-update description rather than a theoretical derivation, so the central claim does not collapse into circularity. The work is self-contained against external benchmarks in the sense that its assertions concern observable code behavior.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The contribution is an implementation update rather than new physics; no free parameters, new axioms, or invented entities are introduced beyond standard stellar atmosphere and radiative transfer assumptions.

axioms (1)

domain assumption Standard assumptions of local thermodynamic equilibrium and radiative transfer in stellar atmospheres hold for the modeled spectra.
The tool relies on established physics for line formation; the abstract does not introduce or question these.

pith-pipeline@v0.9.0 · 5541 in / 1227 out tokens · 63580 ms · 2026-05-12T02:04:20.681582+00:00 · methodology

Review history (3 revisions) →

discussion (0)

Reference graph

Works this paper leans on

47 extracted references · 47 canonical work pages

[1]

1992, Sol

Achmad, L. 1992, Sol. Phys., 138, 411

work page 1992
[2]

Amarsi, A. M. & Asplund, M. 2017, MNRAS, 464, 264

work page 2017
[3]

M., Barklem, P

Amarsi, A. M., Barklem, P. S., Asplund, M., Collet, R., & Zatsarinny, O. 2018, A&A, 616, A89 Article number, page 8 of 12 M. Jian et al.:PySMEv1.0 Fig. 5.Central depths of optical spectral lines for the Sun (top) and Arcturus (bottom). Coloured points mark molecular features, and grey points show all other atomic lines. Fig. 6.Difference between synthetic...

work page 2018
[4]

M., Barklem, P

Amarsi, A. M., Barklem, P. S., Collet, R., Grevesse, N., & Asplund, M. 2019, A&A, 624, A111

work page 2019
[5]

M., Li, W., Grevesse, N., & Jurewicz, A

Amarsi, A. M., Li, W., Grevesse, N., & Jurewicz, A. J. G. 2025, A&A, 703, A35

work page 2025
[6]

M., Liljegren, S., & Nissen, P

Amarsi, A. M., Liljegren, S., & Nissen, P. E. 2022, A&A, 668, A68

work page 2022
[7]

M., Lind, K., Osorio, Y ., et al

Amarsi, A. M., Lind, K., Osorio, Y ., et al. 2020, A&A, 642, A62

work page 2020
[8]

A., Lind, K., & Bergemann, M

Bautista, M. A., Lind, K., & Bergemann, M. 2017, A&A, 606, A127

work page 2017
[9]

J., Eitner, P., et al

Bergemann, M., Gallagher, A. J., Eitner, P., et al. 2019, A&A, 631, A80 Bijavara Seshashayana, S., Jönsson, H., D’Orazi, V ., et al. 2025, A&A, 704, A220

work page 2019
[10]

2014, A&A, 566, A98 Article number, page 9 of 12 A&A proofs:manuscript no

Blanco-Cuaresma, S., Soubiran, C., Jofré, P., & Heiter, U. 2014, A&A, 566, A98 Article number, page 9 of 12 A&A proofs:manuscript no. aa59986-26 Fig. 7.Same as Figure 6, but for the stellar parameters of Arcturus. Fig. 8.The ratio of non-negligible lines across the Kiel diagram for MARCS grid points. Grid points marked with orange circles indicate locatio...

work page 2014
[11]

2021, MNRAS, 506, 150

Buder, S., Sharma, S., Kos, J., et al. 2021, MNRAS, 506, 150

work page 2021
[12]

M., Racca, M., et al

Caliskan, S., Amarsi, A. M., Racca, M., et al. 2025, A&A, 696, A210

work page 2025
[13]

2012, Research in Astronomy and Astrophysics, 12, 1197

Cui, X.-Q., Zhao, Y .-H., Chu, Y .-Q., et al. 2012, Research in Astronomy and Astrophysics, 12, 1197

work page 2012
[14]

C., Abrams, D

Dalton, G., Trager, S. C., Abrams, D. C., et al. 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, V ol. 8446, Ground- based and Airborne Instrumentation for Astronomy IV , ed. I. S. McLean, S. K. Ramsay, & H. Takami, 84460P de Jong, R. S., Agertz, O., Berbel, A. A., et al. 2019, The Messenger, 175, 3 De Silva, G. M., F...

work page 2012
[15]

J., Liu, C., et al

Deng, L.-C., Newberg, H. J., Liu, C., et al. 2012, Research in Astronomy and Astrophysics, 12, 735

work page 2012
[16]

J., Bergemann, M., Collet, R., et al

Gallagher, A. J., Bergemann, M., Collet, R., et al. 2020, A&A, 634, A55

work page 2020
[17]

M., Magg, E., Plez, B., et al

Gerber, J. M., Magg, E., Plez, B., et al. 2023, A&A, 669, A43

work page 2023
[18]

2012, The Messenger, 147, 25

Gilmore, G., Randich, S., Asplund, M., Binney, J., et al. 2012, The Messenger, 147, 25

work page 2012
[19]

Gray, D. F. 2008, The Observation and Analysis of Stellar Photospheres

work page 2008
[20]

A., Sheminova, V

Gurtovenko, E. A., Sheminova, V . A., & Sarychev, A. P. 1991, Sol. Phys., 136, 239

work page 1991
[21]

2008, A&A, 486, 951

Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

work page 2008
[22]

& Lanz, T

Hubeny, I. & Lanz, T. 2011, Synspec: General Spectrum Synthesis Program, Astrophysics Source Code Library, record ascl:1109.022

work page 2011
[23]

C., Dalton, G

Jin, S., Trager, S. C., Dalton, G. B., et al. 2024, MNRAS, 530, 2688 Jofré, P., Heiter, U., & Soubiran, C. 2019, ARA&A, 57, 571

work page 2024
[24]

Kurucz, R. L. 1993b, VizieR Online Data Catalog: Model Atmospheres (Ku- rucz, 1979), VizieR On-line Data Catalog: VI/39. Originally published in: 1979ApJS...40....1K

work page 1979
[25]

S., & Belyaev, A

Lind, K., Asplund, M., Barklem, P. S., & Belyaev, A. K. 2011, A&A, 528, A103

work page 2011
[26]

2013, A&A, 554, A96

Lind, K., Melendez, J., Asplund, M., Collet, R., & Magic, Z. 2013, A&A, 554, A96

work page 2013
[27]

1986, A&A, 163, 135 Article number, page 10 of 12 M

Magain, P. 1986, A&A, 163, 135 Article number, page 10 of 12 M. Jian et al.:PySMEv1.0 Fig. 9.Comparison between the observed solar spectrum in the Hαregion and syn- thetic spectra generated using PySME v0.4, PySME v1.0, and Turbospectrum. Fig. 10.Wall-time comparison of line selection and post-selection synthesis inPySMEfor synthesis windows of different ...

work page 1986
[28]

R., Schiavon, R

Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., Allende Prieto, C., et al. 2017, AJ, 154, 94

work page 2017
[29]

Mallinson, J. W. E., Lind, K., Amarsi, A. M., & Youakim, K. 2024, A&A, 687, A5

work page 2024
[30]

& Lind, K

Nordlander, T. & Lind, K. 2017, A&A, 607, A75

work page 2017
[31]

S., Lind, K., et al

Osorio, Y ., Barklem, P. S., Lind, K., et al. 2015, A&A, 579, A53

work page 2015
[32]

S., Allende Prieto, C., & Zatsarinny, O

Osorio, Y ., Lind, K., Barklem, P. S., Allende Prieto, C., & Zatsarinny, O. 2019, A&A, 623, A103

work page 2019
[33]

& Valenti, J

Piskunov, N. & Valenti, J. A. 2017, A&A, 597, A16

work page 2017
[34]

E., Kupka, F., Ryabchikova, T

Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 112, 525

work page 1995
[35]

2012, Turbospectrum: Code for spectral synthesis, Astrophysics Source Code Library, record ascl:1205.004

Plez, B. 2012, Turbospectrum: Code for spectral synthesis, Astrophysics Source Code Library, record ascl:1205.004

work page 2012
[36]

& Gilmore, G

Randich, S. & Gilmore, G. 2013, The Messenger, 154, 47

work page 2013
[37]

M., Lind, K., et al

Reggiani, H., Amarsi, A. M., Lind, K., et al. 2019, A&A, 627, A177

work page 2019
[38]

2024, A&A, 681, A107

Royer, P., Merle, T., Dsilva, K., et al. 2024, A&A, 681, A107

work page 2024
[39]

L., et al

Ryabchikova, T., Piskunov, N., Kurucz, R. L., et al. 2015, Phys. Scr, 90, 054005

work page 2015
[40]

2012, MOOG: LTE line analysis and spectrum synthesis, Astrophysics Source Code Library, record ascl:1202.009

Sneden, C., Bean, J., Ivans, I., Lucatello, S., & Sobeck, J. 2012, MOOG: LTE line analysis and spectrum synthesis, Astrophysics Source Code Library, record ascl:1202.009

work page 2012
[41]

2011, Science, 333, 176

Tolstoy, E. 2011, Science, 333, 176

work page 2011
[42]

Valenti, J. A. & Piskunov, N. 1996, A&AS, 118, 595

work page 1996
[43]

X., Nordlander, T., Asplund, M., et al

Wang, E. X., Nordlander, T., Asplund, M., et al. 2021, MNRAS, 500, 2159

work page 2021
[44]

2023, A&A, 671, A171

Wehrhahn, A., Piskunov, N., & Ryabchikova, T. 2023, A&A, 671, A171

work page 2023
[45]

J., Abruzzo, M

Wheeler, A. J., Abruzzo, M. W., Casey, A. R., & Ness, M. K. 2023, AJ, 165, 68

work page 2023
[46]

Zatsarinny, O., Bartschat, K., Fernandez-Menchero, L., & Tayal, S. S. 2019, Phys. Rev. A, 99, 023430

work page 2019
[47]

M., Aguirre Børsen-Koch, V ., et al

Zhou, Y ., Amarsi, A. M., Aguirre Børsen-Koch, V ., et al. 2023, A&A, 677, A98 Article number, page 11 of 12 A&A proofs:manuscript no. aa59986-26 Appendix A: NLTE Grid References Table A.1 lists the references for the model atoms used to compute the NLTE departure coefficient grids implemented in PySME. Table A.1.NLTE model atom references. Element Refere...

work page 2023