arxiv: 2605.13374 · v1 · submitted 2026-05-13 · 🌌 astro-ph.SR

Recognition: 2 theorem links

· Lean Theorem

Full non-LTE multi-level radiative transfer II. The case of a 5-level Ca ii atom with broadened excited levels

T. Lagache , M. Sampoorna , F. Paletou

Authors on Pith no claims yet

Pith reviewed 2026-05-14 18:14 UTC · model grok-4.3

classification 🌌 astro-ph.SR

keywords radiative transfernon-LTECa II linespartial frequency redistributionstellar atmospheresspectral line formationmulti-level atomsvelocity-changing collisions

0 comments

The pith

For Ca II H&K and infrared triplet lines in a simplified atmosphere, standard non-LTE with partial redistribution accurately models line formation and makes full non-LTE unnecessary.

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

The paper extends the full non-LTE radiative transfer method to a realistic 5-level Ca II atom that includes natural broadening of excited levels. It solves a coupled set of kinetic equations self-consistently, computing emission and absorption profiles by convolving non-Lorentzian atomic profiles with non-Maxwellian velocity distributions at each step. An efficient iterative scheme based on approximate operator techniques makes the numerical solution practical. Under the conditions examined, the simpler standard NLTE treatment with partial frequency redistribution reproduces the line profiles for the H and K lines and the infrared triplet, even when velocity-changing collisions are included. This result indicates that the added complexity of full non-LTE is not required for these particular lines in the tested setup.

Core claim

Under the conditions studied, for this particular atomic model and for a simplified atmosphere, the standard NLTE with partial redistribution is sufficient to describe the formation of Ca II spectral lines; the more exact treatment of FNLTE is unnecessary in the case of Ca II H & K and infrared triplet lines, even when accounting for velocity-changing collisions.

What carries the argument

An iterative approximate-operator method that self-consistently solves the coupled kinetic equations while convolving non-Lorentzian atomic profiles with non-Maxwellian velocity distributions to obtain emission and absorption profiles for each line.

If this is right

Ca II line profiles in simplified one-dimensional atmospheres can be computed reliably with existing standard NLTE codes.
Velocity-changing collisions do not force the use of full non-LTE for these lines under the tested conditions.
Computational cost for modeling calcium lines in stellar atmospheres can remain at the level of partial-redistribution calculations.
The same numerical framework can be applied to other multi-level atoms once the Ca II case is validated.

Where Pith is reading between the lines

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

Similar conclusions may hold for other resonance lines in low-density stellar atmospheres if their level structure and collision rates are comparable.
Extending the test to three-dimensional, time-dependent atmospheres would be a direct next step to check robustness.
If the result generalizes, existing large-scale spectral synthesis grids for solar and stellar spectra need not incorporate full non-LTE for calcium.

Load-bearing premise

The simplified atmosphere model and chosen physical conditions are representative enough to conclude that full non-LTE is unnecessary for Ca II lines in general.

What would settle it

A direct comparison of observed Ca II H&K or infrared triplet profiles in a real solar or stellar atmosphere against synthetic profiles computed with standard NLTE versus full non-LTE, showing systematic differences only when the full treatment is omitted.

Figures

Figures reproduced from arXiv: 2605.13374 by F. Paletou, M. Sampoorna, T. Lagache.

**Figure 1.** Figure 1: Maximum relative error on atomic densities (see Eq. 36), iteration after iteration with respect to the MALIBU reference solution, calculated with a high degree of precision (8 points per decade, 12 cosine directions µ and 20 azimuths). The red and black lines correspond to the solutions calculated using the LPS25 method, initialised at LTE and CRD respectively. The blue and green lines correspond to the … view at source ↗

**Figure 2.** Figure 2: Schematic view of the atomic model associated with the H & K lines (purple lines) and the infrared triplet (red lines) of Ca ii. The atomic levels are numbered in increasing order of energy. This figure is a slightly modified version of the one shown in [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗

**Figure 3.** Figure 3: Velocity distributions functions, normalised to the Maxwellian, associated with the fraction of atoms in levels 5 (top) and 2 (bottom). Their variation is shown for different optical depths τ = 0, 1, 103 , 105 , 107 . The Maxwellian distribution is shown as a dashed line. Clearly, the VDFs calculated in FNLTE are not Maxwellians. able to provide this new quantities [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗

**Figure 4.** Figure 4: Frequency variation of the emergent intensity at µ = 1, normalised to the Wien function BW , for the H & K lines and the infrared triplet of Ca ii. The coloured lines are the results obtained with the FNLTE formalism, using the MALIBU method. These results are compared with those of the XRD (circles). viates only slightly from the Maxwellian distribution, as we observe in bottom panel of [PITH_FULL_IMA… view at source ↗

**Figure 5.** Figure 5: Normalised emergent intensity at µ = 1 of the K line (left) and velocity distributions at τ = 0 of the fraction of atoms in level 5 (right) shown for different values of the (total) collision rate QE and its velocity-changing component QV . The amount of elastic collision increases from the top to the bottom with γ5,coh = 0.95, 0.80, 0.50, 0.10. The dashed blue (red) line represents the emerging intensity … view at source ↗

read the original abstract

The so-called full non-local thermodynamic equilibrium (FNLTE) radiative transfer problem allows us to take into account not only deviations of the radiation field from the Planckian but also deviations of the densities and velocity distributions of massive particles from Maxwell-Boltzmann statistics. This article discusses the extension of this formalism to physically realistic multi-level atoms, including natural broadening of the excited levels. In practice, we must solve self-consistently a coupled set of kinetic equations and determine, for each line, an emission and absorption profile by convolving a non-Lorentzian atomic profile with a non-Maxwellian velocity distribution at each iteration. To solve this numerically challenging problem, we have developed a new efficient iterative method based on well-known approximate operator techniques. After validating our numerical strategy, we present the results obtained for the H & K lines and the infrared triplet of the Ca II. Under the conditions studied, for this particular atomic model and for a simplified atmosphere, we find that the standard NLTE with partial redistribution is sufficient to describe the formation of Ca II spectral lines. The more exact treatment of FNLTE is unnecessary in the case of Ca II H & K, and infrared triplet lines, even when accounting for velocity-changing collisions.

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

The paper extends FNLTE to a 5-level Ca II atom with a new iterative solver and finds standard NLTE plus partial redistribution works fine for H&K and IR triplet lines in their simplified test case.

read the letter

The core result is that, under the conditions they ran, adding the full non-LTE treatment with non-Maxwellian velocity distributions and velocity-changing collisions does not change the Ca II line profiles enough to matter compared with ordinary NLTE plus partial redistribution. They limit the claim to their particular 5-level model and simplified atmosphere, which keeps the conclusion honest rather than overstated. The new piece is the extension of the FNLTE equations to multi-level atoms that include natural broadening, plus the efficient iterative scheme built on approximate operator methods to handle the repeated profile convolutions. That numerical strategy is the part worth looking at if you need to solve similar coupled kinetic and radiative problems. They report that the solver was validated before the Ca II runs, which is the right order. The main limitation is the atmosphere itself. A simplified static model without depth-dependent velocity gradients or realistic stratification leaves open whether small non-Maxwellian effects could grow in more structured cases, such as moving atmospheres or multi-dimensional simulations. The paper does not claim broader applicability, so the gap is clear rather than hidden. This is useful reading for people who maintain or extend NLTE radiative transfer codes for stellar spectra, especially those already using Ca II diagnostics. It is not a field-changing result, but the method and the specific numerical check are worth having in the literature. I would send it to peer review so the validation details and the exact conditions can be examined by specialists who work on these codes daily.

Referee Report

1 major / 2 minor

Summary. The manuscript extends the full non-LTE (FNLTE) radiative transfer formalism to multi-level atoms with natural broadening of excited levels. It develops an efficient iterative numerical method based on approximate operator techniques to solve the coupled kinetic equations self-consistently, computing emission and absorption profiles via convolution of non-Lorentzian atomic profiles with non-Maxwellian velocity distributions. After validation of the strategy, the method is applied to a 5-level Ca II atom in a simplified atmosphere, yielding the result that standard NLTE with partial redistribution suffices to describe the formation of the Ca II H&K and infrared triplet lines, even when velocity-changing collisions are included.

Significance. If the central result holds, the work is significant because it shows that the more computationally intensive FNLTE treatment is unnecessary for these important Ca II lines under the tested conditions, enabling simpler and faster modeling in stellar atmosphere codes while preserving accuracy. The validated iterative solver for multi-level FNLTE with broadened levels and non-Maxwellian distributions provides a reusable numerical framework that can be extended to other atoms. The finding is grounded in explicit numerical validation and direct comparison to standard NLTE-PRD, strengthening its utility for practical applications in solar and stellar spectroscopy.

major comments (1)

[Application to Ca II and conclusions] The conclusion that standard NLTE with partial redistribution is sufficient (and FNLTE unnecessary) for Ca II H&K and IR triplet lines rests on results obtained exclusively in a simplified atmosphere. The manuscript does not test whether the same holds under realistic stratification, velocity gradients, or multi-dimensional effects, where non-Maxwellian deviations could produce observable profile differences; this limits the load-bearing strength of the general claim for Ca II line formation.

minor comments (2)

[Discussion] The abstract and conclusions already note the simplified atmosphere, but a short dedicated paragraph in the discussion section quantifying the range of temperatures, densities, and velocity fields explored would help readers assess representativeness.
[Numerical method] Clarify in the methods how the convolution of the non-Lorentzian atomic profile with the non-Maxwellian velocity distribution is discretized and updated at each iteration to ensure numerical stability for the 5-level model.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. We address the single major comment below and have incorporated a partial revision to clarify the scope of our conclusions.

read point-by-point responses

Referee: [Application to Ca II and conclusions] The conclusion that standard NLTE with partial redistribution is sufficient (and FNLTE unnecessary) for Ca II H&K and IR triplet lines rests on results obtained exclusively in a simplified atmosphere. The manuscript does not test whether the same holds under realistic stratification, velocity gradients, or multi-dimensional effects, where non-Maxwellian deviations could produce observable profile differences; this limits the load-bearing strength of the general claim for Ca II line formation.

Authors: We agree that the numerical results are obtained in a simplified one-dimensional static atmosphere and that this restricts the generality of the claim. The manuscript already qualifies the finding explicitly (abstract: 'under the conditions studied, for this particular atomic model and for a simplified atmosphere'; similar statements appear in the introduction and conclusions). The central purpose of the paper is the development and validation of the new iterative solver for multi-level FNLTE with natural broadening and non-Maxwellian velocity distributions; the Ca II application serves as a controlled demonstration that isolates the effect of velocity-changing collisions. We do not assert that FNLTE is unnecessary under all conditions. To address the referee's concern, we will add a dedicated paragraph in the discussion section that (i) reiterates the limitations of the current atmospheric model, (ii) notes that non-Maxwellian deviations could become more important in the presence of strong velocity gradients or multi-dimensional flows, and (iii) indicates that the validated method is now available for such extensions. This revision makes the scope of the conclusion unambiguous without altering the reported numerical results. revision: partial

Circularity Check

0 steps flagged

No circularity: numerical results from independent iterative solver

full rationale

The paper extends the FNLTE formalism to a 5-level Ca II atom by solving the coupled kinetic and radiative transfer equations via a new iterative scheme built on established approximate operator techniques. The central claim—that standard NLTE with partial redistribution suffices—is obtained directly from the computed line profiles under the stated simplified atmosphere and collision model; it does not reduce to any fitted parameter, self-definition, or load-bearing self-citation. Validation consists of numerical convergence checks on the same equations, not a renaming or tautological prediction. No step in the derivation chain collapses to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract, the work rests on standard radiative transfer and kinetic equations extended to non-Maxwellian velocities and broadened profiles; no explicit free parameters or invented entities are mentioned.

axioms (1)

standard math Radiative transfer equation coupled with statistical equilibrium equations for level populations and velocity distributions
Invoked as the foundation for the FNLTE problem throughout the abstract.

pith-pipeline@v0.9.0 · 5535 in / 1331 out tokens · 42367 ms · 2026-05-14T18:14:30.598885+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We must solve self-consistently a coupled set of kinetic equations and determine, for each line, an emission and absorption profile by convolving a non-Lorentzian atomic profile with a non-Maxwellian velocity distribution at each iteration.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Under the conditions studied, for this particular atomic model and for a simplified atmosphere, we find that the standard NLTE with partial redistribution is sufficient

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

39 extracted references

[1]

1987, JQSRT, 38, 427

Atanackovi c , O., Borsenberger, J., Oxenius, J., & Simonneau, E. 1987, JQSRT, 38, 427

1987
[2]

Avrett, E. H. 1968, in Resonance Lines in Astrophysics, eds. R. G. Athay, J. Mathis, A. Skumanich (Boulder: NCAR), 27

1968
[3]

1997 a , , 328, 706

Bommier , V. 1997 a , , 328, 706

1997
[4]

1997 b , , 328, 726

Bommier , V. 1997 b , , 328, 726

1997
[5]

1986, JQSRT, 35, 303

Borsenberger, J., Oxenius, J., & Simonneau, E. 1986, JQSRT, 35, 303

1986
[6]

1987, JQSRT, 37, 331

Borsenberger, J., Oxenius, J., & Simonneau, E. 1987, JQSRT, 37, 331

1987
[7]

1987, JQSRT, 37, 397

Huben \'y . 1987, JQSRT, 37, 397

1987
[8]

1982, JQSRT, 27, 593

Huben \'y , I. 1982, JQSRT, 27, 593

1982
[9]

& Cooper, J

Huben \'y , I. & Cooper, J. 1986, ApJ, 305, 852

1986
[10]

& Lites, B

Huben \'y , I. & Lites, B. 1995, ApJ, 455, 376

1995
[11]

& Mihalas, D

Huben \'y , I. & Mihalas, D. 2014, Theory of stellar atmospheres (Princeton University Press)

2014
[12]

& Oxenius, J

Huben \'y , I. & Oxenius, J. 1987, JQSRT, 37, 65

1987
[13]

1983, JQSRT, 29, 477

Huben \'y , I., Oxenius, J., & Simonneau, E. 1983, JQSRT, 29, 477

1983
[14]

Hummer, D. G. 1962, MNRAS, 125, 21

1962
[15]

Jefferies, J. T. 1968, Spectral line formation (A Blaisdell Book in the Pure and Applied Sciences, Waltham, Mass.: Blaisdell, 1968)

1968
[16]

2025a, , 699, A198

Lagache , T., Paletou , F., & Sampoorna , M. 2025a, , 699, A198
[17]

2025b, RAS Techniques and Instruments, Volume 4, id.rzaf023

Lagache, T., Sampoorna, M., & Paletou, F. 2025b, RAS Techniques and Instruments, Volume 4, id.rzaf023
[18]

Atomic Spectra Database, 2026, https://physics.nist.gov/asd

NIST. Atomic Spectra Database, 2026, https://physics.nist.gov/asd

2026
[19]

L., Auer , L

Olson , G. L., Auer , L. H., & Buchler , J. R. 1986, JQSRT, 35, 431

1986
[20]

1972, ApJ, vol

Omont, A., Smith, E., & Cooper, J. 1972, ApJ, vol. 175, p. 185, 175, 185

1972
[21]

1965, JQSRT, 5, 771

Oxenius, J. 1965, JQSRT, 5, 771

1965
[22]

1986, Kinetic theory of particles and photons (Berlin: Springer)

Oxenius, J. 1986, Kinetic theory of particles and photons (Berlin: Springer)

1986
[23]

& Simonneau, E

Oxenius, J. & Simonneau, E. 1994, Annals of Physics, 234, 60

1994
[24]

1995, A&A, 302, 587

Paletou, F. 1995, A&A, 302, 587

1995
[25]

& Auer, L

Paletou, F. & Auer, L. H. 1995, A&A, 297, 771

1995
[26]

& Lagache, T

Paletou, F. & Lagache, T. 2025, https://hal.science/hal-05051883

2025
[27]

& Peymirat, C

Paletou, F. & Peymirat, C. 2021, A&A, 649, A165

2021
[28]

2020, A&A, 633, A111

Paletou, F., Peymirat, C., Anterrieu, E., & B \"o hm, T. 2020, A&A, 633, A111

2020
[29]

2023, A&A, 671, A93

Paletou, F., Sampoorna, M., & Peymirat, C. 2023, A&A, 671, A93

2023
[30]

Rybicki, G. B. & Hummer, D. G. 1991, A&A, 245, 171

1991
[31]

Sampoorna, M., Nagendra, K., & Stenflo, J. O. 2013, , 770, 92

2013
[32]

2024, A&A, 690, A213

Sampoorna, M., Paletou, F., Bommier, V., & Lagache, T. 2024, A&A, 690, A213

2024
[33]

Shine, R. A. & Linsky, J. L. 1974, Solar Physics, 39, 49

1974
[34]

1989, A&A, 213, 360

Uitenbroek, H. 1989, A&A, 213, 360

1989
[35]

2001, ApJ, 557, 389

Uitenbroek, H. 2001, ApJ, 557, 389

2001
[36]

1931, Annalen der Physik, 401, 23

Weisskopf, V. 1931, Annalen der Physik, 401, 23

1931
[37]

1933, Zeitschrift f \"u r Physik, 85, 451

Weisskopf, V. 1933, Zeitschrift f \"u r Physik, 85, 451

1933
[38]

1931, MNRAS, Vol

Woolley, R. 1931, MNRAS, Vol. 91, p. 977, 91, 977

1931
[39]

1938, MNRAS, Vol

Woolley, R. 1938, MNRAS, Vol. 98, p. 624, 98, 624

1938