pith. sign in

arxiv: 2509.12404 · v5 · pith:SHYAXT6Fnew · submitted 2025-09-15 · ⚛️ nucl-th

Thermal effects on stellar neutron capture reactions: a quantum dynamical approach

N. Lightfoot , A. Diaz-Torres , P. Stevenson This is my paper

Pith reviewed 2026-05-25 08:29 UTC · model grok-4.3

classification ⚛️ nucl-th
keywords neutron capturethermal effectsTDCCWPr-processHauser-Feshbach188Osdynamical couplingstellar nucleosynthesis
0
0 comments X

The pith

Thermal effects reduce the neutron capture cross section on 188Os as temperature increases due to dynamical nuclear couplings.

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

This paper applies the time-dependent coupled channels wave-packet method to neutron capture on osmium in a thermal stellar environment. It uses both a many-body nuclear potential and a temperature-dependent initial state to represent the mixed thermal population of nuclear levels. Calculations show the capture cross section falling as temperature rises, with reaction rates also dropping at the highest energies examined. This trend arises because dynamical couplings between the excited states open a dominant pathway involving higher neutron speed. The result differs from statistical Hauser-Feshbach models and carries consequences for the rapid neutron capture process that forms heavy elements.

Core claim

TDCCWP calculations indicate a decrease of the n+188Os capture cross section with increasing temperature, along with a decrease in reaction rates for the highest thermal energies studied, which are contrary to Hauser-Feshbach calculations and important in the rapid neutron capture process. The physical reason for this discrepancy is the key role of the dynamical nuclear coupling between the thermally populated states of the target nucleus, which is neglected in the Hauser-Feshbach approach, but creates a dominant neutron capture pathway with increased neutron speed and thus reduces the neutron capture cross section.

What carries the argument

The time-dependent coupled channels wave-packet (TDCCWP) method with a many-body nuclear potential and a temperature-dependent initial state.

If this is right

  • The n+188Os capture cross section decreases with increasing temperature.
  • Reaction rates decrease for the highest thermal energies studied.
  • This outcome is important for the rapid neutron capture process.
  • The dynamical nuclear coupling between thermally populated states drives the effect and is absent from the Hauser-Feshbach approach.

Where Pith is reading between the lines

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

  • Current statistical models of the r-process may need revision if dynamical couplings prove widespread across other nuclei.
  • The TDCCWP method could be extended to additional target nuclei to map temperature trends in capture rates.
  • Controlled experiments that vary temperature while tracking capture probabilities would provide a direct test.

Load-bearing premise

The many-body nuclear potential together with the chosen temperature-dependent initial state faithfully reproduces the dynamical couplings between thermally populated states without uncontrolled approximations that would reverse the reported temperature trend.

What would settle it

An independent calculation or measurement showing that the n+188Os neutron capture cross section increases rather than decreases with temperature at the studied thermal energies.

Figures

Figures reproduced from arXiv: 2509.12404 by A. Diaz-Torres, N. Lightfoot, P. Stevenson.

Figure 1
Figure 1. Figure 1: FIG. 1. The nuclear interaction potential between a neutron [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. All of the different potentials in the model for the [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Temperature-dependent capture cross sections for a [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The same as in Fig. 4, however, the range of incident [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. A comparison between the capture cross sections [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Reaction rates with different types of cross sections [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
read the original abstract

The neutron capture process plays a vital role in creating the heavy elements in the universe. Astrophysical environments involved in these processes are characterized by two distinct reaction mechanisms: the slow and rapid neutron capture processes. In this work, the slow neutron capture process is described with the time-dependent coupled channels wave-packet (TDCCWP) method that uses both a many-body nuclear potential and an initial temperature-dependent state to account for the thermal environment. To evaluate the role of a mixed and entangled initial state in the temperature-dependent neutron capture cross section, TDCCWP calculations are compared with those from the coupled-channels density matrix (CCDM) method based on the Lindblad equation. The importance of including temperature in the initial wave-function of the TDCCWP approach is compared to a thermalisation of the reaction rate using a Hauser-Feshbach style approach. TDCCWP calculations indicate a decrease of the n+$^{188}$Os capture cross section with increasing temperature, along with a decrease in reaction rates for the highest thermal energies studied, which are contrary to Hauser-Feshbach calculations and important in the rapid neutron capture process. The physical reason for this discrepancy is the key role of the dynamical nuclear coupling between the thermally populated states of the target nucleus, which is neglected in the Hauser-Feshbach approach, but creates a dominant neutron capture pathway with increased neutron speed and thus reduces the neutron capture cross section.

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.

Referee Report

2 major / 0 minor

Summary. The paper introduces the time-dependent coupled-channels wave-packet (TDCCWP) method, which incorporates a many-body nuclear potential and a temperature-dependent initial state, to compute thermal effects on neutron capture. For the n + 188Os reaction it reports that the capture cross section decreases with rising temperature and that reaction rates fall at the highest thermal energies examined; this trend is opposite to Hauser-Feshbach statistical-model results and is ascribed to dynamical couplings among thermally populated target states that open a faster-neutron capture pathway.

Significance. If the reported temperature dependence survives detailed numerical validation, the result would be significant for r-process nucleosynthesis modeling because it indicates that quantum-dynamical couplings neglected in standard statistical treatments can reverse the expected thermal trend in neutron-capture rates.

major comments (2)
  1. Abstract and Results: the central numerical claim—a decrease of the n+188Os capture cross section with temperature—is stated without any tabulated cross sections, error bars, convergence tests, or benchmark values, so it is impossible to judge whether the TDCCWP propagation actually supports the reported trend.
  2. Methods (TDCCWP initial-state construction): the temperature-dependent initial state is asserted to capture the mixed/entangled thermal ensemble whose dynamical couplings produce the faster-neutron pathway; however, the manuscript supplies no explicit checks (normalization, phase stability, or comparison with the thermal density matrix) that would rule out preparation artifacts capable of reversing the temperature dependence.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed review and constructive feedback on our manuscript. The comments correctly identify areas where additional numerical documentation would strengthen the presentation of the TDCCWP results. We address each major comment below and will incorporate the requested material in a revised version. The central finding—that dynamical couplings in the thermal ensemble reverse the temperature trend relative to Hauser-Feshbach—remains supported by the TDCCWP versus CCDM comparison, but we agree that explicit tables and checks are needed for full transparency.

read point-by-point responses

  1. Referee: Abstract and Results: the central numerical claim—a decrease of the n+188Os capture cross section with temperature—is stated without any tabulated cross sections, error bars, convergence tests, or benchmark values, so it is impossible to judge whether the TDCCWP propagation actually supports the reported trend.

    Authors: We agree that tabulated values, uncertainties, and convergence information are essential for assessing the numerical reliability of the reported trend. In the revised manuscript we will add a table listing capture cross sections at T = 0, 0.25, 0.5, 0.75, and 1.0 MeV (with statistical uncertainties from the wave-packet ensemble), together with convergence tests versus channel basis size, spatial grid spacing, and propagation time step. Low-temperature benchmarks against known experimental data will also be included to confirm that the TDCCWP implementation reproduces established results before the thermal trend is examined. revision: yes

  2. Referee: Methods (TDCCWP initial-state construction): the temperature-dependent initial state is asserted to capture the mixed/entangled thermal ensemble whose dynamical couplings produce the faster-neutron pathway; however, the manuscript supplies no explicit checks (normalization, phase stability, or comparison with the thermal density matrix) that would rule out preparation artifacts capable of reversing the temperature dependence.

    Authors: The initial state is constructed as a coherent superposition of target eigenstates with Boltzmann weights, which by design encodes the mixed thermal ensemble; its norm is preserved by unitary propagation. The consistency of the temperature trend between TDCCWP and the independent CCDM (Lindblad) calculation already provides evidence that the result is not an artifact of state preparation. Nevertheless, we will add an explicit subsection in the Methods that reports (i) the initial-state norm at t=0 and after propagation, (ii) phase stability of the dominant components, and (iii) direct comparison of expectation values of the initial density operator with the corresponding thermal density matrix, thereby ruling out preparation-induced reversals of the temperature dependence. revision: partial

Circularity Check

0 steps flagged

No significant circularity; temperature trend follows from explicit TDCCWP propagation

full rationale

The derivation computes neutron capture cross sections via direct time-dependent propagation of an initial temperature-dependent wave packet under a many-body potential, then extracts the reaction probability from the final state. This is compared to independent CCDM and Hauser-Feshbach results without any step that renames a fitted quantity as a prediction or reduces the central claim to a self-citation. The reported decrease with temperature is an output of the dynamical evolution rather than an input by construction, rendering the chain self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the TDCCWP propagation and the Lindblad description for the comparison method; both are standard domain tools whose accuracy for this thermal neutron-capture regime is assumed rather than re-derived.

axioms (2)
  • domain assumption The time-dependent coupled-channels wave-packet method with the chosen many-body potential accurately captures dynamical couplings between thermally populated nuclear states.
    Invoked as the basis for the TDCCWP results and the claimed physical mechanism.
  • domain assumption The Lindblad master equation provides a faithful open-system description for comparing entangled versus mixed initial states in the CCDM calculations.
    Used to benchmark the role of the temperature-dependent initial state.

pith-pipeline@v0.9.0 · 5786 in / 1568 out tokens · 42045 ms · 2026-05-25T08:29:32.684497+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

35 extracted references · 35 canonical work pages

  1. [1]

    Initialize the wave-packet att= 0 on the radial grid

    work page

  2. [2]

    Time steps can be ob- served to understand the fluctuations in popula- tions of the internal energy levels of the 188Os tar- get

    Propagate the wave-packet along the grid using the time-evolution operator. Time steps can be ob- served to understand the fluctuations in popula- tions of the internal energy levels of the 188Os tar- get

    work page

  3. [3]

    The above process is carried out using a modified Chebyshev propagation scheme [15–17] to apply the ex- ponential time-evolution operator to the wave-packet

    Allow the reflected wave-packet to move back to its original position and extract the transmission coef- ficients from the model to generate capture cross sections and reaction rates. The above process is carried out using a modified Chebyshev propagation scheme [15–17] to apply the ex- ponential time-evolution operator to the wave-packet. Throughout the ...

    work page

  4. [4]

    The WS potential is unable to fit the HF potential at small radii because of shell effects in the Hartree-Fock po- tential

    This potential is a Woods-Saxon potential and can be fitted for the parameters,V 0,a 0, andR 0, using a New- tonian least squares regression, where the Woods-Saxon potential is defined as: VN = V0 1 + exp(r−R0 a0 ) .(1) HereV 0 is the strength of the potential,a 0 is the value for the diffuseness of the potential andR 0 is the radius of the nucleus define...

    work page

  5. [5]

    TheV N(r) defines the interaction potential of the system, and the W(r) is the absorption potential

    and the second is the centrifugal term. TheV N(r) defines the interaction potential of the system, and the W(r) is the absorption potential. The complete coupled channels Hamiltonian can be described by a matrix as shown below: ˆH|ψ⟩= ˆH0 +V 0−0 V0−2 V2−0 ˆH0 +V 2−2 +ϵ 2 |ψ0⟩ |ψ2⟩ .(6) Here and for many other even mass nuclei, 0 + is the ground state and ...

    work page

  6. [6]

    Pollak and P

    E. Pollak and P. Talkner, Reaction rate theory: What it was, where is it today, and where is it going?, Chaos: An Interdisciplinary Journal of Nonlinear Science15(2005)

    work page 2005

  7. [7]

    Langanke and M

    K. Langanke and M. Wiescher, Nuclear reactions and stellar processes, Reports on Progress in Physics64, 1657 (2001)

    work page 2001

  8. [8]

    Wiescher, C

    M. Wiescher, C. Bertulani, C. Brune, R. Deboer, A. Diaz-Torres, L. Gasques, K. Langanke, P. Navr´ atil, W. Nazarewicz, J. Oko lowicz,et al., Quantum physics of stars, Reviews of Modern Physics97, 025003 (2025)

    work page 2025

  9. [9]

    Bracci, G

    L. Bracci, G. Fiorentini, V. Melezhik, G. Mezzorani, and P. Quarati, Atomic effects in the determination of nuclear cross sections of astrophysical interest, Nuclear Physics A513, 316 (1990)

    work page 1990

  10. [10]

    Segawa, T

    M. Segawa, T. Masaki, Y. Nagai, Y. Temma, T. Shima, K. Mishima, M. Igashira, S. Goriely, A. Koning, and S. Hilaire, Neutron capture cross sections of os 186, os 187, and os 189 for the re-os chronology, Physical Re- view C76, 022802 (2007)

    work page 2007

  11. [11]

    Kazakov, S

    L. Kazakov, S. Savvidis, C. Eleftheriadis, A. Policarpo, Y. Popov, H. Leeb, C. Pretel, P. K¨ ohler, S. Raman, V. Vlachoudis,et al.,The Re Os clock revisited, Tech. Rep. (2000)

    work page 2000

  12. [12]

    Hagino, N

    K. Hagino, N. Rowley, and A. Kruppa, A program for coupled-channel calculations with all order couplings for heavy-ion fusion reactions, Computer Physics Communi- cations123, 143 (1999)

    work page 1999

  13. [13]

    Bustreo, G

    C. Bustreo, G. Casini, G. Zollino, T. Bolzonella, and R. Piovan, Fresco, a simplified code for cost analysis of fusion power plants, Fusion Engineering and Design88, 3141 (2013)

    work page 2013

  14. [14]

    Lee and A

    I. Lee and A. Diaz-Torres, Coherence dynamics in low- energy nuclear fusion, Physics Letters B827, 136970 (2022)

    work page 2022

  15. [15]

    Vockerodt and A

    T. Vockerodt and A. Diaz-Torres, Describing heavy-ion fusion with quantum coupled-channels wave-packet dy- namics, Physical Review C100, 034606 (2019)

    work page 2019

  16. [16]

    Boselli and A

    M. Boselli and A. Diaz-Torres, Quantifying low-energy fusion dynamics of weakly bound nuclei from a time- dependent quantum perspective, Physical Review C92, 044610 (2015)

    work page 2015

  17. [17]

    Diaz-Torres, L

    A. Diaz-Torres, L. Gasques, and N. Antonenko, Cluster effects on low-energy carbon burning, Physics Letters B 849, 138476 (2024)

    work page 2024

  18. [18]

    Close, P

    G. Close, P. Stevenson, and A. Diaz-Torres, Quantum dy- namical microscopic approach to stellar carbon burning, Physics Letters B , 139881 (2025). 9

    work page 2025

  19. [19]

    Winters and R

    R. Winters and R. Macklin, Maxwellian-averaged neu- tron capture cross sections for tc-99 and mo-95-98, As- trophysical Journal, Part 1 (ISSN 0004-637X), vol. 313, Feb. 15, 1987, p. 808-812.313, 808 (1987)

    work page 1987

  20. [20]

    Chen and H

    R. Chen and H. Guo, The chebyshev propagator for quantum systems, Computer physics communications 119, 19 (1999)

    work page 1999

  21. [21]

    V. A. Mandelshtam, T. P. Grozdanov, and H. S. Taylor, Bound states and resonances of the hydroperoxyl radical ho2: An accurate quantum mechanical calculation using filter diagonalization, The Journal of chemical physics 103, 10074 (1995)

    work page 1995

  22. [22]

    Diaz-Torres and M

    A. Diaz-Torres and M. Wiescher, Characterizing the as- trophysical s factor for c 12+ c 12 fusion with wave-packet dynamics, Physical Review C97, 055802 (2018)

    work page 2018

  23. [23]

    J. A. Maruhn, P.-G. Reinhard, P. Stevenson, and A. S. Umar, The tdhf code sky3d, Computer Physics Commu- nications185, 2195 (2014)

    work page 2014

  24. [24]

    Seideman and W

    T. Seideman and W. H. Miller, Quantum mechanical re- action probabilities via a discrete variable representation- absorbing boundary condition green’s function, The Journal of chemical physics97, 2499 (1992)

    work page 1992

  25. [25]

    Potnis, Nuclear energy levels of os 188, Physical Re- view104, 722 (1956)

    V. Potnis, Nuclear energy levels of os 188, Physical Re- view104, 722 (1956)

    work page 1956

  26. [26]

    I. Lee, G. Gosselin, and A. Diaz-Torres, Thermal and atomic effects on coupled-channels heavy-ion fusion, Physical Review C107, 054609 (2023)

    work page 2023

  27. [27]

    D. H. Boal, Approach to chemical equilibrium in thermal models, Physical Review C29, 967 (1984)

    work page 1984

  28. [28]

    Thompson, J

    I. Thompson, J. E. Escher, and G. Arbanas, Coupled- channel models of direct-semidirect capture via giant- dipole resonances, Nuclear Data Sheets118, 292 (2014)

    work page 2014

  29. [29]

    I. J. Thompson, Coupled reaction channels calculations in nuclear physics, Computer Physics Reports7, 167 (1988)

    work page 1988

  30. [30]

    Hodgson, The neutron optical potential, Reports on Progress in Physics47, 613 (1984)

    P. Hodgson, The neutron optical potential, Reports on Progress in Physics47, 613 (1984)

    work page 1984

  31. [31]

    B. R. Martin and G. Shaw,Nuclear and particle physics: an introduction(John Wiley & Sons, 2019)

    work page 2019

  32. [32]

    J. N. Avila, M. Lugaro, T. R. Ireland, F. Gyngard, E. Zinner, S. Cristallo, P. Holden, J. Buntain, S. Amari, and A. Karakas, Tungsten isotopic compositions in star- dust sic grains from the murchison meteorite: constraints on the s-process in the hf–ta–w–re–os region, The Astro- physical Journal744, 49 (2011)

    work page 2011

  33. [33]

    Humayun and A

    M. Humayun and A. D. Brandon, s-process implications from osmium isotope anomalies in chondrites, The As- trophysical Journal664, L59 (2007)

    work page 2007

  34. [34]

    Grundl, V

    J. Grundl, V. Spiegel, C. Eisenhauer, H. Heaton, D. Gilliam, and J. Bigelow, A californium-252 fission spectrum irradiation facility for neutron reaction rate measurements, Nuclear Technology32, 315 (1977)

    work page 1977

  35. [35]

    Vertebny, M

    V. Vertebny, M. Vlasov, A. Dadakina, R. Zatserkovsky, A. Kirilyuk, M. Pasechnik, and N. Trofimova, Slow neu- tron total cross-sections and resonance levels of osmium- 187, 188 and 189, The Ukranian Journal of Physics14, 1971 (1969)

    work page 1971