Glauber-theory calculations of high-energy nuclear scattering observables using variational Monte Carlo wave functions
Pith reviewed 2026-05-16 20:39 UTC · model grok-4.3
The pith
Glauber theory with variational Monte Carlo wave functions accurately describes high-energy scattering for light nuclei.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Glauber-theory calculations that employ variational Monte Carlo wave functions for 6He and 12C reproduce the measured elastic differential cross sections and total reaction cross sections for p+12C, 12C+12C, and 6He+12C collisions at intermediate to high energies. The full multiple-scattering phase-shift function is obtained by Monte Carlo integration, and its cumulant expansion is shown to converge rapidly already at second order.
What carries the argument
The Glauber phase-shift function evaluated by Monte Carlo integration over nucleon-nucleon multiple scatterings, built from variational Monte Carlo nuclear wave functions that include two- and three-body correlations.
Load-bearing premise
The Glauber multiple-scattering approximation itself remains valid for the chosen beam energies and target nuclei.
What would settle it
New high-precision experimental data showing significant deviations from the calculated differential cross sections for any of the three systems at comparable energies would indicate that the description fails.
Figures
read the original abstract
Experiments using intermediate- to high-energy radioactive nuclear beams present numerous findings. Extracting important properties of physical observables relies on a firm theoretical analysis. Though Glauber theory is believed to work well, no convincing calculation has so far been done. We perform ab initio Glauber theory calculations of both elastic differential cross sections and total reaction cross sections for p+12C, 12C+12C, and 6He+12C systems. The wave functions of both 6He and 12C are generated by variational Monte Carlo calculations with spatial and spin-isospin correlations induced by realistic two- and three-nucleon potentials. Glauber's phase-shift function is computed by Monte Carlo integration up to all orders of nucleon-nucleon multiple scatterings. We show an excellent performance of the Glauber description to the selected data on the above systems. We also find that the cumulant expansion of the phase-shift function converges rapidly up to the second order for the above systems. This finding will open up interesting applications for the analysis of high-energy nuclear experiments.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript performs ab initio Glauber-theory calculations of elastic differential cross sections and total reaction cross sections for the p+12C, 12C+12C, and 6He+12C systems at intermediate-to-high energies. Variational Monte Carlo wave functions are generated from realistic two- and three-nucleon potentials; the Glauber phase-shift function is evaluated by direct Monte Carlo integration over all orders of nucleon-nucleon multiple scatterings. The central results are excellent agreement with selected experimental data and rapid convergence of the cumulant expansion of the phase-shift function up to second order.
Significance. If the numerical results hold, the work supplies the first fully ab initio, parameter-free Glauber calculations for these systems using realistic VMC wave functions, thereby strengthening the theoretical foundation for analyzing high-energy radioactive-beam experiments. The explicit Monte Carlo integration to all orders and the reported rapid cumulant convergence constitute concrete technical advances that could simplify future applications.
major comments (2)
- [Abstract] Abstract and results section: the claim of 'excellent performance' of the Glauber description rests exclusively on agreement with the chosen data sets; no independent benchmark against non-eikonal multiple-scattering calculations, exact few-body solutions, or prior calculations with different wave functions is provided to test the validity of the eikonal approximation itself for the selected energies and nuclei.
- [Results] Results section (Monte Carlo integration paragraph): quantitative details on statistical errors, number of samples, and convergence tests for the full phase-shift function (beyond the cumulant expansion) are not reported, making it impossible to assess whether the claimed all-order integration is numerically stable for the 6He+12C system.
minor comments (2)
- [Abstract] The beam energies for each reaction should be stated explicitly in the abstract and introduction rather than referred to only as 'the above systems'.
- [Figures] Figure captions should include the specific experimental references and the kinematic range shown so that the data-theory comparison can be reproduced without consulting external sources.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and have revised the manuscript to incorporate the suggestions where feasible.
read point-by-point responses
-
Referee: [Abstract] Abstract and results section: the claim of 'excellent performance' of the Glauber description rests exclusively on agreement with the chosen data sets; no independent benchmark against non-eikonal multiple-scattering calculations, exact few-body solutions, or prior calculations with different wave functions is provided to test the validity of the eikonal approximation itself for the selected energies and nuclei.
Authors: We appreciate the referee's emphasis on validating the eikonal approximation independently. Our manuscript's central contribution is the first fully ab initio implementation of Glauber theory for these systems using realistic VMC wave functions, with the phase-shift function computed via direct Monte Carlo integration over all orders of multiple scattering. While we agree that additional benchmarks against non-eikonal or exact few-body calculations would be valuable, performing such new calculations lies outside the scope of the present work, which focuses on the integration of ab initio structure inputs with the standard Glauber framework. In the revised manuscript we have added a dedicated paragraph in the introduction and discussion sections referencing prior literature that has tested the eikonal approximation for light nuclei at comparable energies, thereby providing context for the applicability of our approach. revision: partial
-
Referee: [Results] Results section (Monte Carlo integration paragraph): quantitative details on statistical errors, number of samples, and convergence tests for the full phase-shift function (beyond the cumulant expansion) are not reported, making it impossible to assess whether the claimed all-order integration is numerically stable for the 6He+12C system.
Authors: We agree that explicit quantitative information on the Monte Carlo sampling procedure is essential for assessing numerical stability. In the revised manuscript we have expanded the relevant paragraph in the Results section to report the number of Monte Carlo samples used (10^6 configurations per nucleus), the estimated statistical uncertainties on the full phase-shift function (typically 1-3% for 6He+12C), and explicit convergence tests demonstrating that the all-order results stabilize with increasing sample size. These additions directly address the concern for the 6He+12C system. revision: yes
Circularity Check
No circularity: ab initio computation from established potentials and VMC wave functions
full rationale
The derivation begins with realistic two- and three-nucleon potentials and variational Monte Carlo wave functions for 6He and 12C. The Glauber phase-shift function is obtained by direct Monte Carlo integration over these wave functions up to all orders of multiple scattering. No parameters are adjusted to the p+12C, 12C+12C or 6He+12C scattering data under study; the reported agreement with experiment and the rapid convergence of the cumulant expansion to second order are presented as outcomes, not inputs. No self-citation chain or ansatz is invoked to force the central results. The calculation is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Glauber multiple-scattering series is a valid description of high-energy nucleon-nucleus and nucleus-nucleus collisions
- domain assumption Variational Monte Carlo with realistic two- and three-nucleon potentials produces sufficiently accurate ground-state wave functions for 6He and 12C
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Glauber’s phase-shift function is computed by Monte Carlo integration up to all orders of nucleon-nucleon multiple scatterings... cumulant expansion of the phase-shift function converges rapidly up to the second order
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
-
[1]
R. J. Glauber, Lectures in Theoretical Physics , edited by W. E. Brittin and L. G. Dunham (Interscience, New York, 1959), Vol. 1, p. 315
work page 1959
-
[2]
K. C. Mathur, A. N. Tripathi, and S. K. Joshi, Phys. Rev. A 5, 746 (1972)
work page 1972
-
[3]
J. Ryckebusch, D. Debruyne, P. Lava, S. Janssen, B. Van Overmeire, and T. Van Cauteren, Nucl. Phys. A 728, 226 (2003)
work page 2003
-
[4]
D. d’Enterria and C. Loizides, Annu. Rev. Nucl. Part. Sci. 71, 315 (2021). 5
work page 2021
-
[5]
I. Tanihata, H. Hamagaki, O. Hashimoto, Y. Shida, N. Yoshikawa, K. Sugimoto, O. Yamakawa, T. Kobayashi, and N. Takahashi, Phys. Rev. Lett. 55, 2676 (1985)
work page 1985
- [6]
- [7]
-
[8]
I. Tanihata, H. Savajols, and R. Kanungo, Prog. Part. Nucl. Phys. 68, 215 (2013), and references therein
work page 2013
- [9]
- [10]
- [11]
- [12]
-
[13]
W. Horiuchi, Y. Suzuki, P. Capel, and D. Baye, Phys. Rev. C 81, 024606 (2010)
work page 2010
-
[14]
W. Horiuchi, Y. Suzuki, and T. Inakura, Phys. Rev. C 89, 011601(R) (2014)
work page 2014
-
[15]
W. Horiuchi, S. Hatakeyama, S. Ebata, and Y. Suzuki, Phys. Rev. C 93, 044611 (2016)
work page 2016
-
[16]
W. Horiuchi, S. Hatakeyama, S. Ebata, and Y. Suzuki, Phys. Rev. C 96, 024605 (2017)
work page 2017
-
[17]
R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C 51, 38 (1995)
work page 1995
-
[18]
R. B. Wiringa, R. Schiavilla, S. C. Pieper, and J. Carlso n, Phys. Rev. C 89, 024305 (2014)
work page 2014
-
[19]
J. Carlson, S. Gandolfi, F. Pederiva, S. C. Pieper, R. Schiavilla, K. E. Schmidt, and R. B. Wiringa, Rev. Mod. Phys. 87, 1067 (2015)
work page 2015
- [20]
- [21]
-
[22]
L. Andreoli, G. B. King, S. Pastore, M. Piarulli, J. Carl - son, S. Gandolfi, and R. B. Wiringa, Phys. Rev. C 110, 064004 (2024)
work page 2024
- [23]
-
[24]
M. Piarulli, S. Pastore, R. B. Wiringa, S. Brusilow, and R. Lim, Phys. Rev. C 107, 014314 (2023)
work page 2023
-
[25]
K. Kaki, Y. Suzuki, and R. B. Wiringa, Phys. Rev. C 86, 044601 (2012)
work page 2012
- [26]
-
[27]
B. Abu-Ibrahim, W. Horiuchi, A. Kohama, and Y. Suzuki, Phys. Rev. C 77, 034607 (2008); ibid 80, 029903(E) (2009); 81, 019901(E) (2010)
work page 2008
- [28]
-
[29]
R. Kubo, J. Phys. Soc. Jpn. 17, 1100 (1962)
work page 1962
- [30]
- [31]
-
[32]
O. Rubies-Bigorda, S. Ostermann, and S. F. Yelin, Phys. Rev. Res. 5, 013091 (2023)
work page 2023
- [33]
- [34]
-
[35]
J. S. Al-Khalili and J. A. Tostevin, Phys. Rev. Lett. 76, 3903 (1996)
work page 1996
-
[36]
J. S. Al-Khalili and J. A. Tostevin, Phys. Rev. C 54, 1843 (1996)
work page 1996
-
[37]
W. Horiuchi, Y. Suzuki, B. Abu-Ibrahim, and A. Ko- hama, Phys. Rev. C 75, 044607 (2007); ibid 76, 039903(E) (2007)
work page 2007
-
[38]
Yu. M. Shabelski and A. G. Shuvaev, Phys. Rev. C 104, 064607 (2021)
work page 2021
-
[39]
H. O. Meyer, P. Schwandt, H. P. Gubler, W. P. Lee, W. T. H. van Oers, R. Abegg, D. A. Hutcheon, C. A. Miller, P. Helmer, K. P. Jackson, Phys. Rev. C 31, 1569 (1985)
work page 1985
-
[40]
A. Okamoto, T. Yamagata, H. Akimune, M. Fujiwara, K. Fushimi, M. B. Greenfield, K. Hara, K. Y. Hara, H. Hashimoto, R. Hayami et al. , Phys. Rev. C 81, 054604 (2010)
work page 2010
-
[41]
G. W. Hoffmann, M. L. Barlett, D. Ciskowski, G. Pauletta, M. Purcell, L. Ray, J. F. Amann, J. J. Jarmer, K. W. Jones, S. Penttil¨ a et al. , Phys. Rev. C 41, 1651 (1990)
work page 1990
-
[42]
G. S. Blanpied, G. W. Hoffmann, M. L. Barlett, J. A. McGill, S. J. Greene, L. Ray, O. B. Van Dyck, J. Amann, and H. A. Thiessen, Phys. Rev. C 23, 2599 (1991)
work page 1991
-
[43]
H. Palevsky, J. L. Friedes, R. J. Sutter, G. W. Bennett, G. J. Igo, W. D. Simpson, G. C. Phillips, D. M. Corley, N. S. Wall, R. L. Stearns, and B. Gottschalk, Phys. Rev. Lett. 18, 1200 (1967)
work page 1967
-
[44]
G. D. Alkhazov, G. M. Amalsky, S. L. Belostotsky, A. A. Vorobyov, O. A. Domchenkov, Yu. V. Dotsenko, and V. E. Starodubsky, Phys. Lett. B 42 , 121 (1972)
work page 1972
-
[45]
M. Takechi, M. Fukuda, M. Mihara, T. Chinda, T. Mat- sumasa, H. Matsubara, Y. Nakashima, K. Matsuta, T. Minamisono, R. Koyama et al. , Eur. Phys. J. A 25, s01, 217 (2005)
work page 2005
-
[46]
L. Ponnath, T. Aumann, C. A. Bertulani, R. Gernh¨ auser, M. Heil, T. Almusidi, H. Alvarez-Pol, L. Atar, L. Atkins, Y. Ayyad et al. , Phys. Lett. B 855 , 138780 (2024)
work page 2024
- [47]
-
[48]
M. Takechi, S. Suzuki, D. Nishimura, M. Fukuda, T. Ohtsubo, M. Nagashima, T. Suzuki, T. Yamaguchi, A. Ozawa, T. Moriguchi et al , Phys. Rev. C 90, 061305(R) (2014)
work page 2014
- [49]
-
[50]
M. Buenerd, J. Pinston, J. Cole, C. Guet, D. Lebrun, J. M. Loisea, P. Martin, E. Monnard, J. Mougey, H. Nifenecher et al., Phys. Lett. B 102 , 242 (1981)
work page 1981
-
[51]
J. Y. Hostachy, M. Buenerd, J. Chauvin, D. Lebrun, Ph. Martin, B. Bonin, G. Bruge, J. C. Lugol, L. Papineau, P. Roussel, et al. , Phys. Lett. B 184 , 139 (1987)
work page 1987
-
[52]
I. Tanihata, H. Hamagaki, O. Hashimoto, S. Nagamiya, Y. Shida, N. Yoshikawa, O. Yamakawa, K. Sugimoto, T. Kobayashi, D. E. Greiner et al. , Phys. Lett. B 206 , 592 (1985)
work page 1985
-
[53]
J. Arrington, N. Fomin, and A. Schmidt, Annu. Rev. Nucl. Part. Sci. 72, 307 (2022). 6
work page 2022
-
[54]
J. Zenihiro, H. Sakaguchi, T. Murakami, M. Yosoi, Y. Yasuda, S. Terashima, Y. Iwao, H. Takeda, M. Itoh, H. P. Yoshida, and M. Uchida, Phys. Rev. C 82, 044611 (2010)
work page 2010
-
[55]
D. Adhikari, H. Albataineh, D. Androic, K. Aniol, D. S. Armstrong, T. Averett, C. Ayerbe Gayoso, S. Barcus, V. Bellini, B. S. Beminiwattha et al. , Phys. Rev. Lett. 126, 172502 (2021)
work page 2021
-
[56]
G. Giacalone, G. Nijs, and W. van der Schee, Phys. Rev. Lett. 131, 202302 (2023)
work page 2023
- [57]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.