REVIEW 3 major objections 3 minor 50 references
The d-electron contribution to the stopping power of transition metals
T0 review · 3 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The low-energy stopping of protons in six transition metals is reproduced by the atomic d-subshell momentum distribution, not by a free-electron gas.
desk verdict A genuinely new d-electron stopping model for transition metals; the atomic-to-solid momentum distribution assumption is the key thing a referee should test. 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 object is the inhomogeneous momentum distribution $f_{nl}(p)=(2\pi)^{3/2}|\Phi_{nl}(p)|^2$: the squared Fourier transform of the atomic d-subshell wavefunction, normalized to the number $N_d$ of bound d electrons. It replaces the Fermi-sphere step function of the free-electron-gas model, so d-electrons of all momenta contribute, and it is what carries the low-energy d-electron response. The stopping cross section is then computed by folding this distribution with the transport cross section $\sigma_{\mathrm{tr}}(v_r)$ (the momentum-transfer cross section) generated by a velocity-dependent screened potential that satisfies the cusp condition; the Fourier transforms are evaluated analytically after expanding the wavefunctions in exponential-type basis functions. The model becomes a full stopping theory by adding free-electron-gas and inner-shell contributions, covering 0.1 keV to 100 MeV.
What would settle it
Measure proton stopping in a single crystal of Ni or Au along a well-characterized off-channelling direction at velocities below about 0.3 atomic units with a few-percent uncertainty and compare with the model curve: a deviation larger than the stated few-percent agreement, or the appearance of a sharp slope break that the model excludes, would falsify the free-atom d-electron description. A less costly test is to recompute $f_{3d}(p)$ using a band-structure d-projected wavefunction and check whether the low-velocity total stopping shifts by more than the model's current scatter.
Extended reading notes
Core claim
The paper's central discovery is that the stopping of protons in the later transition metals is not governed by how the d-electrons sit in the Fermi sphere but by their full atomic momentum profile. Concretely, the d-subshell stopping cross section is obtained by replacing the Fermi step function with $f_{nl}(p) = (2\pi)^{3/2}|\Phi_{nl}(p)|^2$ in the transport-cross-section integral, with $\Phi_{nl}(p)$ the Fourier transform of the atomic d wavefunction normalized to the number $N_d$ of d electrons that remain bound. The integer $N_{\mathrm{FEG}}$ of electrons in the free-electron gas is taken as the integer closest to the value implied by the measured plasmon frequency, so the only solid-state input is that frequency. Adding the resulting d-curve to a free-electron-gas curve and to inner-shell curves yields total stopping cross sections that match recent low-energy experiments closely for Ni, Pd, and Pt, lie between the conflicting data sets for Cu, Ag, and Au, and agree with available TDDFT results. The authors interpret the low-energy 'break of linearity' seen in group 11 measurements as a soft nonlinearity generated by the inhomogeneous momentum distribution, with the d-curve peaking near the mean d-electron velocity $v_d$.
Load-bearing premise
The load-bearing premise is that the d electrons that stay bound in the solid respond like free-atom d electrons, with the number of such electrons fixed by rounding the conduction-electron count from the measured plasmon frequency to an integer; if solid-state band effects or that electron-count assignment are wrong, the low-energy d contribution changes.
Editorial extensions
If this is right
- For Ni, Pd, and Pt the calculated total stopping cross sections fall within a few percent of recent low-energy measurements, so the model provides a quantitative account of the group 10 low-energy response without a slope break.
- For Cu, Ag, and Au the calculated curves sit between the two conflicting families of low-energy data, clarifying that the historical 'break of linearity' is not reproduced as a sharp kink but as a smooth d-electron-driven rise.
- The d-subshell contribution to the stopping cross section is largest when the projectile speed is close to the mean d-electron speed $v_d$; this places the d-curve maximum at a distinct impact energy for each metal (about 350 keV/amu for Ni, about 100 keV for Pt, with the other metals in between).
- Combining the d-model with free-electron-gas and inner-shell models yields a fully theoretical stopping-power curve over 0.1 keV to 100 MeV, with the only electron-count parameter fixed by the measured plasmon frequency.
Reading between the lines
- Beyond the paper: the same inhomogeneous-momentum construction could be tested on other d-block metals (for example, group 9 or group 12), where low-energy stopping data are sparser; the only inputs needed are a d-subshell wavefunction and a measured or computed plasmon frequency.
- Beyond the paper: replacing the free-atom d wavefunctions in Eq. (7) with solid-state d-projected wavefunctions would isolate the size of band-structure and band-gap effects the authors identify as a possible source of low-velocity overestimation.
- Beyond the paper: a systematic comparison of d-only stopping curves (not total stopping) between this model and TDDFT across all six metals would test whether the atomic momentum profile captures the same physics as explicit time-dependent electron dynamics.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a non-perturbative model for the d-electron contribution to the electronic stopping cross section of the late transition metals Ni, Pd, Pt, Cu, Ag, and Au. The stopping integral of Eq. (1) is evaluated with an atomic "inhomogeneous" momentum distribution f_nl(p) built from Hartree-Fock or relativistic wave functions (Eq. (7)), and this d contribution is combined with a free-electron-gas description of the valence electrons and a shellwise local approximation for inner shells. The authors compare their low-velocity results with experimental data from the IAEA database and with available TDDFT calculations, and present extended-energy curves from 0.1 keV to 100 MeV. The central claim is that the d-electron stopping response can be captured by atomic momentum distributions rather than by a free-electron gas or full TDDFT, and that this yields good agreement with experiment and with TDDFT where available.
Significance. If validated, the model is attractive because it is simple and analytic, it requires no fitting to stopping data, and the electron-count inputs (N_FEG, N_d) are inferred from experimental plasmon frequencies. The paper's strengths include a transparent derivation, extensive comparison with a large experimental database, and direct comparison with TDDFT results for several targets. The main obstacle is that the central assumption — that free-atom d-orbital momentum distributions represent the solid-state d-electron response — is not tested against any solid-state calculation, and the sensitivity of the results to the N_FEG/N_d partition is not quantified. These issues are addressable, but they are load-bearing for the claimed agreement at low velocities.
major comments (3)
- [Sec. II, Eq. (7); Sec. III A] The central assumption is that the d-electron response in the solid is described by the free-atom momentum distribution f_nl(p). In the low-velocity limit, Eqs. (5)-(6) weight the distribution at p ~ v, so the low-energy stopping power is controlled by the low-momentum tail of f_nd(p). In a transition metal this tail is affected by band structure, s-d hybridization, and d-band occupation, none of which is captured by an isolated-atom orbital. The authors' caveat at the end of Sec. III A that a band-gap effect may cause overestimation at the lowest velocities is a symptom of this atomic-versus-solid gap, but it does not test the shape of f_nd(p) itself. I request a quantitative comparison with a DFT-based d-projected momentum density (or equivalent solid-state calculation) for at least one Group-10 and one Group-11 target, with a discussion of how hybridization changes the low-p weight and how that would affect the computed stopping curves.
- [Table I; Sec. III] The partition of outer electrons into N_FEG and N_d is based on integer rounding of the experimental plasmon frequency through the free-electron relation, together with the free-atom configuration; it is not validated by a solid-state occupation calculation. Because the d contribution is proportional to N_d through the normalization in Eq. (7), the low-velocity results depend directly on this assignment. The authors state that N_FEG and N_d 'may be the subject of discussion' and compare with an N_FEG=10 variant for Pt in Figs. 6 and 12, but this is not a systematic sensitivity analysis for the present model. I would like to see either a quantitative sensitivity study over the admissible integer choices of N_FEG/N_d for all six targets, or an independent estimate of the solid d-occupation from DFT projected densities of states.
- [Eqs. (6), (A11)-(A13)] As printed, the integration limits in I'(v_r), g5, and g7 are reversed. After exchanging the order of integration in Eqs. (3)-(4), the p-integral should run from |v_r - v| to v_r + v. The printed limits |v_r + v| to |v_r - v| decay over the integration interval and would give the negative of the intended value. This is likely typographical, but it appears in the central derivation and in the analytical expressions of Appendix A, so it must be corrected and the numerical implementation must be confirmed to use the correct limits.
minor comments (3)
- [Eq. (9)] The Heaviside function is written as Θ(p - p_F), but the text says the step function lies 'within the Fermi sphere'; this should be Θ(p_F - p). In addition, the numerical constant (2π)^{3/2} on the right-hand side is inconsistent with the Fermi-sphere volume factor: with p_F = (3π^2 n_e)^{1/3}, one has ∫ d^3p = 4π^3 n_e, not (2π)^{3/2} n_e. Please clarify the normalization convention and reconcile Eq. (9) with Eq. (8).
- [Table I caption] The phrase 'the dump value γ_p^exp' appears to be a typo for 'damping value'; the text around Table I should be corrected.
- [Sec. III A, Figs. 2-7] The comparison with TDDFT for Au is complicated by the fact that the cited TDDFT results are channelling calculations for Au<100>; the authors note this, but the figure caption could make the distinction clearer to avoid the impression of a direct off-channelling comparison.
Circularity Check
No significant circularity: the d-electron stopping is computed from atomic wave functions and independently inferred electron counts, and then compared with external experimental and TDDFT benchmarks.
full rationale
The paper's derivation chain is a forward calculation, not a fit to the target data. Equation (7) defines the d-electron momentum distribution from Hartree-Fock or relativistic atomic wave functions, and Eqs. (5)-(6) then compute the stopping cross section by integrating that distribution against a transport cross section. No stopping-power measurement is used to set any parameter in this chain. The electron counts NFEG and Nd are inferred from experimental plasmon frequencies quoted from Werner et al. [40], an independent measurement, with NFEG taken as the nearest integer and Nd as the complement; this affects the magnitude of the d contribution but is not tuned to the stopping data. The velocity-dependent screened potential is taken from the authors' prior work [4], but it is a fixed input model that is itself being tested by the present comparisons rather than a result derived from the target claim. The validation step, where the computed curves are compared with IAEA data and with TDDFT, is external to the derivation. Even the acknowledged low-velocity caveat about the missing band gap is a stated limitation of the atomic approximation, not evidence that the prediction is equivalent to an input. The only parts of the paper that involve self-citations are the use of the authors' earlier screened-potential and shell-model ingredients and the comparison with their own earlier FEG-only calculation; these are normal building blocks and do not make the central claim reduce to its inputs. The central claim has independent content: it predicts quantitative stopping curves from independent atomic-structure and plasmon-frequency inputs, and those predictions are then checked against experiment and TDDFT.
Assumptions & free parameters
free parameters (2)
- NFEG (integer number of electrons in the free-electron gas) =
Ni=3, Cu=3, Pd=7, Pt=7, Ag=3, Au=7
- Nd (number of d-bound electrons) =
Ni=7, Cu=8, Pd=3, Pt=3, Ag=8, Au=4
assumptions (4)
- standard math The standard kinetic-theory formula, Eq. (1), with transport cross-section and phase shifts, gives the stopping cross-section from a momentum distribution.
- domain assumption The velocity-dependent screened potential of Ref. [4] verifies the cusp condition and correctly describes the projectile-target interaction.
- domain assumption Free-atom Hartree-Fock or relativistic atomic wave functions represent the bound d-electrons in the solid target.
- domain assumption Independent-shell approximation: total electronic stopping is the sum of d-electron, free-electron-gas, and inner-shell contributions.
Cite this review
Pith. "Pith review of The d-electron contribution to the stopping power of transition metals." pith.science (2026). https://pith.science/paper/OJ7ZBE4Y
@misc{pith2026241112810,
author = {Pith},
title = {Pith review of: The d-electron contribution to the stopping power of transition metals},
year = {2026},
howpublished = {\url{https://pith.science/paper/OJ7ZBE4Y}},
note = {Machine review of arXiv:2411.12810}
}
abstract
We present a new non-perturbative model to describe the stopping power by ionization of the $d$-electrons of transition metals. These metals are characterized by the filling of the d-subshell and the promotion of part of the electrons to the conduction band. The contribution of d-electrons at low-impact energies has been noted experimentally in the past as a break of the linear dependence of the stopping power with the ion velocity. In this contribution, we describe the response of these electrons considering the atomic "inhomogeneous" momentum distribution. We focus on the transition metals of Groups 10 and 11 in the periodic table: Ni, Pd, Pt, Cu, Ag, and Au. Results describe the low energy-stopping power, with good agreement with the experimental data and available TDDFT results. By combining the present non-perturbative model for the $d$-subshell contribution with other approaches for the valence electrons and for the inner shells, we provide a coherent theoretical method capable of describing the stopping power of these transition metals from the very low to the high energy region.
Figures
Figures from the paper (7 more)
Reference graph
Works this paper leans on
-
[1]
G. J. Leigh, ed., Principles of Chemical Nomenclature A Guide to IUPAC Recommendations (RSC Publishing, 2011)
work page 2011
-
[2]
Electron Binding Energies of the Ele- ments,
G. P. Williams, “Electron Binding Energies of the Ele- ments,” (CRC Press, 2014) Chap. 10, pp. 200–205, 95th ed
work page 2014
-
[3]
C. C. Montanari and J. E. Miraglia, Nuclear Instruments and Methods in Physics Research Section B: Beam Inter- actions with Materials and Atoms 460, 27 (2019)
work page 2019
-
[4]
Only two sets of data are available at low velocities,
-
[5]
J. E. Vald´ es, P. Vargas, and V. A. Esaulov, Radiation Effects and Defects in Solids 171, 60 (2016)
work page 2016
-
[6]
C. C. Montanari and J. E. Miraglia, Phys. Rev. A 96, 012707 (2017)
work page 2017
-
[7]
E. D. Cantero, G. H. Lantschner, J. C. Eckardt, and N. R. Arista, Phys. Rev. A 80, 032904 (2009)
work page 2009
-
[8]
S. N. Markin, D. Primetzhofer, M. Spitz, and P. Bauer, Phys. Rev. B 80, 205105 (2009)
work page 2009
Show all 50 references
-
[9]
S. N. Markin, D. Primetzhofer, S. Prusa, M. Brunmayr, G. Kowarik, F. Aumayr, and P. Bauer, Phys. Rev. B 78, 195122 (2008)
2008
-
[10]
J. E. Vald´ es, J. C. Eckardt, G. H. Lantschner, and N. R. Arista, Phys. Rev. A 49, 1083 (1994)
1994
-
[11]
Vald´ es, G
J. Vald´ es, G. Tamayo, G. Lantschner, J. Eckardt, and N. Arista, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 73, 313 (1993)
1993
-
[12]
D. Roth, C. Celedon, D. Goebl, E. Sanchez, B. Bruckner, R. Steinberger, J. Guimpel, N. Arista, and P. Bauer, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 437, 1 (2018)
2018
-
[13]
Our results de- scribe these values quite well; the total stopping power is slightly larger than the experimental data, with differ- ences that are less than 5%
and [17], which agree pretty well among them and show an almost perfect linear dependency of the stop- ping power with the impact velocity. Our results de- scribe these values quite well; the total stopping power is slightly larger than the experimental data, with differ- ence...
2022
-
[14]
Goebl, D
D. Goebl, D. Roth, and P. Bauer, Phys. Rev. A 87, 062903 (2013)
2013
-
[15]
Celed´ on, E
C. Celed´ on, E. A. S´ anchez, M. S. Moreno, N. R. Arista, J. D. Uribe, M. Mery, J. E. Vald´ es, and P. Vargas, Phys. Rev. A 88, 012903 (2013)
2013
-
[16]
C. E. Celed´ on, E. A. S´ anchez, L. Salazar Alarc´ on, J. Guimpel, A. Cort´ es, P. Vargas, and N. R. Arista, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 360, 103 (2015)
2015
-
[17]
S. P. M¨ oller, A. Csete, T. Ichioka, H. Knudsen, U. I. Uggerhøj, and H. H. Andersen, Phys. Rev. Lett. 88, 193201 (2002)
2002
-
[18]
Bruckner, D
B. Bruckner, D. Roth, D. Goebl, P. Bauer, and D. Primetzhofer, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Ma- terials and Atoms 423, 82 (2018)
2018
-
[19]
Vald´ es, C
J. Vald´ es, C. Agurto, F. Ortega, P. Vargas, R. Labb´ e, and N. Arista, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Ma- terials and Atoms 164-165, 268 (2000)
2000
-
[20]
Primetzhofer, Phys
D. Primetzhofer, Phys. Rev. B 86, 094102 (2012)
2012
-
[21]
C. D. Denton, I. Abril, J. C. Moreno-Mar ´ ın, S. Heredia- Avalos, and R. Garcia-Molina, physica status solidi (b) 245, 1498 (2008)
2008
-
[22]
Abril, P
I. Abril, P. de Vera, and R. Garcia-Molina, Journal of Physics: Conference Series 2326, 012016 (2022)
2022
-
[23]
M. A. Zeb, J. Kohanoff, D. S´ anchez-Portal, A. Arnau, J. I. Juaristi, and E. Artacho, Phys. Rev. Lett. 108, 225504 (2012)
2012
-
[24]
E. E. Quashie, B. C. Saha, and A. A. Correa, Phys. Rev. B 94, 155403 (2016)
2016
-
[25]
E. E. Quashie and A. A. Correa, Phys. Rev. B 98, 235122 (2018)
2018
-
[26]
S.-M. Li, F. Mao, X.-D. Zhao, W.-Q. Jin, W.-Q. Zuo, B.-S. Li, F. Wang, and F.-S. Zhang, Phys. Rev. B 106, 014103 (2022)
2022
-
[27]
C.-K. Li, X. Guo, J.-M. Xue, and F.-S. Zhang, Phys. Rev. A 107, 052814 (2023)
2023
-
[28]
X.-D. Zhao, F. Mao, and H. Deng, Phys. Rev. A 109, 032807 (2024)
2024
-
[29]
Runge and E
E. Runge and E. K. U. Gross, Phys. Rev. Lett. 52, 997 (1984)
1984
-
[30]
Electronic Stop- ping Power of Matter for Ions Graphs, Data, Comments and Programs,
IAEA, H. Paul, and C. C. Montanari, “Electronic Stop- ping Power of Matter for Ions Graphs, Data, Comments and Programs,” (1928-2024)
1928
-
[31]
Montanari, P
C. Montanari, P. Dimitriou, L. Marian, A. Mendez, J. Peralta, and F. Bivort-Haiek, Nuclear Instruments and Methods in Physics Research Section B: Beam Inter- actions with Materials and Atoms 551, 165336 (2024)
2024
-
[32]
N. D. Mermin, Phys. Rev. B 1, 2362 (1970)
1970
-
[33]
C. C. Montanari and J. E. Miraglia, inAdvances in Quan- tum Chemistry: Theory of Heavy Ion Collision Physics in Hadron Therapy, Vol. 2, edited by D. Belkic (Elsevier, New York, 2013) Chap. 7, pp. 165–201
2013
-
[34]
J. P. Peralta, M. Fiori, A. M. P. Mendez, and C. C. Montanari, Phys. Rev. A 105, 062814 (2022)
2022
-
[35]
Nagy and A
I. Nagy and A. Bergara, Nuclear Instruments and Meth- ods in Physics Research Section B: Beam Interactions with Materials and Atoms 115, 58 (1996), atomic Colli- sions in Solids
1996
-
[36]
N. Wang, I. Nagy, and P. M. Echenique, Phys. Rev. B 58, 2357 (1998)
1998
-
[37]
Bunge, J
C. Bunge, J. Barrientos, and A. Bunge, Atomic Data and Nuclear Data Tables 53, 113 (1993)
1993
-
[38]
Klapisch, J
M. Klapisch, J. L. Schwob, B. S. Fraenkel, and J. Oreg, Journal of the Optical Society of America 67, 148 (1977). 10
1977
-
[39]
Bar-Shalom, M
A. Bar-Shalom, M. Klapisch, and J. Oreg, Journal of Quantitative Spectroscopy and Radiative Transfer 71, 169 (2001)
2001
-
[40]
A. M. P. Mendez, C. C. Montanari, and D. M. Mitnik, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 460, 114 (2019)
2019
-
[41]
M. R. Flannery and I. Levy, H., The Journal of Chemical Physics 50, 2938 (1969)
1969
-
[42]
W. S. M. Werner, K. Glantschnig, and C. Ambrosch- Draxl, Journal of Physical and Chemical Reference Data 38, 1013 (2009)
2009
-
[43]
Z. Lin, L. V. Zhigilei, and V. Celli, Phys. Rev. B 77, 075133 (2008)
2008
-
[44]
Nomura and S
A. Nomura and S. Kiyono, Journal of Physics D: Applied Physics 8, 1551 (1975)
1975
-
[45]
Morita, H
K. Morita, H. Akimune, and T. Suita, Journal of the Physical Society of Japan 25, 1525 (1968)
1968
-
[46]
Mart ´ ınez-Tamayo, J
G. Mart ´ ınez-Tamayo, J. C. Eckardt, G. H. Lantschner, and N. R. Arista, Phys. Rev. A 54, 3131 (1996)
1996
-
[47]
Semrad, P
D. Semrad, P. Mertens, and P. Bauer, Nuclear Instru- ments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 15, 86 (1986)
1986
-
[48]
M. V. Moro, P. Bauer, and D. Primetzhofer, Phys. Rev. A 102, 022808 (2020)
2020
-
[49]
Eppacher and D
C. Eppacher and D. Semrad, Nuclear Instruments and Methods in Physics Research Section B: Beam Interac- tions with Materials and Atoms 69, 33 (1992)
1992
-
[50]
F. F. Selau, H. Trombini, G. G. Marmitt, A. M. De An- drade, J. Morais, P. L. Grande, I. Alencar, M. Vos, and R. Heller, Physical Review A 102, 032812 (2020). Appendix A: Momentum distribution calculations In the present model, we propose an inhomoge- neous momentum distributi...
2020
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.