pith. sign in

arxiv: 2506.14338 · v4 · submitted 2025-06-17 · ⚛️ nucl-th · nucl-ex

Phenomenological refinement of p-d elastic scattering descriptions towards the 3NF study in nuclei via the (p,pd) reaction

Yoshiki Chazono , Tokuro Fukui , Futoshi Minato , Yukinobu Watanabe , Kazuyuki Ogata This is my paper

Pith reviewed 2026-05-19 09:43 UTC · model grok-4.3

classification ⚛️ nucl-th nucl-ex
keywords p-d scatteringthree-nucleon forcesLegendre polynomialsphenomenological model(p,pd) reactioneffective interactionsnuclear reactions
0
0 comments X

The pith

The p-d elastic scattering cross section is improved by modeling the residual amplitude with Legendre polynomials whose coefficients have quadratic energy dependence.

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

This paper develops a phenomenological method to refine calculations of proton-deuteron elastic scattering, which is key for studying three-nucleon forces inside nuclei through the (p,pd) reaction. The approach splits the scattering amplitude into a part from two-nucleon effective interactions and a residual part that is represented by a sum of Legendre polynomials. Coefficients for these polynomials are fitted to experimental differential cross-section data at different energies and then approximated by quadratic functions for smooth energy variation. With this, the model reproduces experimental data accurately across a broad energy range. This serves as an initial step for including three-nucleon force effects in nuclear medium descriptions.

Core claim

The p-d elastic amplitude is decomposed into a 2N part using effective interactions and a residual part approximated by a superposition of Legendre polynomials. The coefficients of these polynomials are adjustable parameters determined by fitting to experimental p-d differential cross-section data at various incident energies. These parameters exhibit smooth energy dependence that is captured by quadratic functions, allowing the analytic form to also reproduce the data well.

What carries the argument

Residual part of the p-d amplitude modeled as superposition of Legendre polynomials with energy-dependent coefficients fitted to data.

If this is right

  • The approach improves the quantitativity of p-d scattering cross sections calculated with effective interactions.
  • It enables better theoretical descriptions for the (p,pd) reaction as a probe of 3NFs in nuclei.
  • The quadratic parametrization provides a practical way to use the model at energies where direct data is limited.
  • Results with the analytic energy dependence match experimental observations across wide incident energies.

Where Pith is reading between the lines

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

  • This fitting strategy could be extended to other few-body scattering problems where effective interactions fall short.
  • Applying the refined amplitudes in distorted-wave calculations for (p,pd) on nuclei might reveal medium modifications to 3NFs.
  • Future comparisons with precise data at intermediate energies could test the quadratic approximation's validity.
  • The method offers a bridge between microscopic 2N calculations and phenomenological adjustments for practical nuclear reaction modeling.

Load-bearing premise

The part of the p-d amplitude not captured by 2N effective interactions can be adequately described by a superposition of Legendre polynomials with coefficients that vary smoothly and quadratically with energy.

What would settle it

Measuring p-d differential cross sections at a new energy not used in the fit and finding that the quadratic functions fail to predict the data accurately would falsify the approach.

Figures

Figures reproduced from arXiv: 2506.14338 by Futoshi Minato, Kazuyuki Ogata, Tokuro Fukui, Yoshiki Chazono, Yukinobu Watanabe.

Figure 1
Figure 1. Figure 1: FIG. 1. Incident-energy dependence of [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Differential cross sections of the [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
read the original abstract

The ($p,pd$) reaction is expected to be a powerful tool for probing three-nucleon forces (3NFs) in nuclear medium since it can be essentially regarded as the $p$-$d$ elastic scattering inside nuclei. One of the important points in the theoretical description of the ($p,pd$) reaction is to calculate the $p$-$d$ scattering in a nucleus quantitatively using effective interactions. This work aims to develop a phenomenological approach to improve the quantitativity of the $p$-$d$ scattering cross section in free space calculated with effective interactions. The $p$-$d$ elastic amplitude is decomposed into a 2N part, described using 2N effective interactions, and a residual part, which the 2N part cannot describe. The latter is approximated by a superposition of Legendre polynomials, with coefficients treated as adjustable parameters. These parameters are determined to reproduce experimental $p$-$d$ differential cross-section data at various incident energies. The obtained parameters exhibit smooth energy dependence, which is approximated by quadratic functions. The numerical results with the analytic energy dependence also reproduce the experimental data. The developed approach works well for improving the $p$-$d$ scattering cross section in a wide range of incident energies. This work can be regarded as the first step toward the description of ($p,pd$) reactions taking 3NF effect in nuclear medium into account.

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

3 major / 2 minor

Summary. The manuscript develops a phenomenological approach to refine p-d elastic scattering amplitudes for use in (p,pd) reaction studies. The amplitude is decomposed into a 2N effective-interaction contribution and a residual part; the latter is expanded in Legendre polynomials whose coefficients are fitted directly to experimental differential cross-section data at selected incident energies. These coefficients are then replaced by quadratic functions of energy, and the resulting analytic expressions are shown to reproduce the input data over a range of energies. The work is framed as a first step toward incorporating 3NF effects in nuclear-medium calculations.

Significance. If the quadratic parametrization of the residual coefficients is shown to be robust and transferable beyond the fitted energies, the method could supply a practical, improved phenomenological input for p-d scattering in nuclei, helping to isolate 3NF contributions in (p,pd) analyses. It builds directly on existing effective-interaction frameworks and addresses a concrete calculational need.

major comments (3)
  1. [§3] §3 (Amplitude decomposition): The residual amplitude is defined by subtracting the 2N effective-interaction contribution from the full amplitude, yet the manuscript provides no test or argument that this residual is primarily due to 3NF rather than incomplete treatment of the 2N force, relativistic effects, or other higher-order terms. This distinction is load-bearing for the stated goal of 3NF studies in nuclei.
  2. [§4] §4 (Fitting procedure and quadratic approximation): Coefficients are first adjusted to reproduce cross-section data at discrete energies and then replaced by quadratic fits to those same adjusted values. Consequently, reproduction of the data at the fitted energies follows by construction; the manuscript does not report predictions or comparisons at energies withheld from the fit, leaving the claimed reliability over a 'wide range' untested.
  3. [Results] Results section: No sensitivity study is presented for the choice of quadratic energy dependence (e.g., comparison with linear or cubic forms, or variation of the number of Legendre terms). Without such checks, the smoothness assumption remains an unquantified modeling choice rather than a demonstrated feature.
minor comments (2)
  1. [Abstract] Abstract: The assertion that the approach 'works well' should be tempered to reflect that agreement is achieved through direct fitting followed by interpolation.
  2. [Formalism] Notation: Define explicitly how the residual amplitude is normalized and subtracted in the formalism; the current description leaves the precise relation between the full, 2N, and residual amplitudes ambiguous.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough review and insightful comments on our manuscript. We address each of the major comments below and outline the revisions we intend to implement.

read point-by-point responses

  1. Referee: [§3] §3 (Amplitude decomposition): The residual amplitude is defined by subtracting the 2N effective-interaction contribution from the full amplitude, yet the manuscript provides no test or argument that this residual is primarily due to 3NF rather than incomplete treatment of the 2N force, relativistic effects, or other higher-order terms. This distinction is load-bearing for the stated goal of 3NF studies in nuclei.

    Authors: We agree that the manuscript does not provide a test or argument demonstrating that the residual is primarily due to 3NF effects. The residual is phenomenologically defined as the difference between the full amplitude and the 2N contribution, and may include other contributions. The framing as a step toward 3NF studies in nuclei is motivational, based on the expectation that (p,pd) can probe 3NF. We will revise the discussion in §3 to clarify this point and note that attributing the residual to specific physical mechanisms would require additional theoretical investigations beyond the scope of this work. revision: partial

  2. Referee: [§4] §4 (Fitting procedure and quadratic approximation): Coefficients are first adjusted to reproduce cross-section data at discrete energies and then replaced by quadratic fits to those same adjusted values. Consequently, reproduction of the data at the fitted energies follows by construction; the manuscript does not report predictions or comparisons at energies withheld from the fit, leaving the claimed reliability over a 'wide range' untested.

    Authors: The referee correctly notes that agreement at the fitted energies is expected by construction. The quadratic energy dependence was chosen to capture the smooth variation observed in the fitted coefficients. To better support the reliability over a wide range, we will perform and report an additional comparison at an energy not used in the fitting process, demonstrating the predictive capability of the quadratic parametrization. revision: yes

  3. Referee: [Results] Results section: No sensitivity study is presented for the choice of quadratic energy dependence (e.g., comparison with linear or cubic forms, or variation of the number of Legendre terms). Without such checks, the smoothness assumption remains an unquantified modeling choice rather than a demonstrated feature.

    Authors: We chose the quadratic form and the number of Legendre terms based on the quality of fit and the apparent smoothness in the data. However, we acknowledge the value of a sensitivity study. In the revised manuscript, we will add a subsection or paragraph in the Results section presenting comparisons with alternative functional forms (linear and cubic) and different numbers of Legendre polynomials to assess the robustness of our parametrization. revision: yes

Circularity Check

1 steps flagged

Fitting Legendre coefficients to data then quadratics to those coefficients makes reproduction by construction

specific steps
  1. fitted input called prediction [Abstract]
    "These parameters are determined to reproduce experimental p-d differential cross-section data at various incident energies. The obtained parameters exhibit smooth energy dependence, which is approximated by quadratic functions. The numerical results with the analytic energy dependence also reproduce the experimental data."

    Coefficients are first fitted directly to cross-section data at selected energies; quadratic analytic forms are then fitted to those same coefficient values. The subsequent statement that the quadratic forms reproduce the data is therefore a direct consequence of the fitting chain rather than an independent prediction or validation.

full rationale

The paper decomposes the p-d amplitude into a 2N effective-interaction term plus a residual term. The residual is modeled as a sum of Legendre polynomials whose coefficients are adjusted to fit experimental differential cross sections at discrete energies. These fitted coefficient values are then replaced by quadratic functions of incident energy that are themselves fitted to the discrete fitted values. The resulting analytic expressions are shown to reproduce the same data. Because the quadratic parametrization has no independent derivation and is calibrated directly on the already-fitted coefficients, the reported agreement follows by construction from the two-stage fitting procedure rather than from an external test or first-principles constraint. No other load-bearing steps reduce to self-citation or self-definition.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The approach rests on the decomposition into 2N and residual amplitudes and on the assumption that residual coefficients vary smoothly enough to be captured by low-order polynomials; both are domain assumptions rather than derived results.

free parameters (1)
  • Legendre polynomial coefficients for residual amplitude
    Adjustable parameters introduced to reproduce experimental p-d differential cross sections at multiple incident energies.
axioms (1)
  • domain assumption p-d elastic amplitude decomposes into a 2N effective-interaction part plus a residual part not captured by 2N interactions
    Explicitly stated as the starting point for the phenomenological refinement.

pith-pipeline@v0.9.0 · 5817 in / 1277 out tokens · 38199 ms · 2026-05-19T09:43:47.383443+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

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

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

21 extracted references · 21 canonical work pages

  1. [1]

    Weinberg, Phenomenological Lagrangians, Phys

    S. Weinberg, Phenomenological Lagrangians, Phys. A96, 327 (1979)

    work page 1979

  2. [2]

    Epelbaum, Few-nucleon forces and systems in chiral effective field theory, Prog

    E. Epelbaum, Few-nucleon forces and systems in chiral effective field theory, Prog. Part. Nucl. Phys. 57, 654 (2006)

    work page 2006

  3. [3]

    Machleidt and D

    R. Machleidt and D. R. Entem, Chiral effective field the- ory and nuclear forces, Phys. Rep. 503, 1 (2011)

    work page 2011

  4. [4]

    Weinberg, Three-body interactions among nucleons and pions, Phys

    S. Weinberg, Three-body interactions among nucleons and pions, Phys. Lett. B 295, 114 (1992)

    work page 1992

  5. [5]

    van Kolck, Few-nucleon forces from chiral La- grangians, Phys

    U. van Kolck, Few-nucleon forces from chiral La- grangians, Phys. Rev. C 49, 2932 (1994)

    work page 1994

  6. [6]

    Hanhart, U

    C. Hanhart, U. van Kolck, and G. A. Miller, Chi- ral Three-Nucleon Forces from p-wave Pion Production, Phys. Rev. Lett. 85, 2905 (2000)

    work page 2000

  7. [7]

    Navr´ atil, V

    P. Navr´ atil, V. G. Gueorguiev, J. P. Vary, W. E. Or- mand, and A. Nogga, Structure of A = 10–13 Nuclei with 6 Two- Plus Three-Nucleon Interactions from Chiral Effec- tive Field Theory, Phys. Rev. Lett. 99, 042501 (2007)

    work page 2007

  8. [8]

    Fukui, G

    T. Fukui, G. De Gregorio, and A. Gargano, Uncovering the mechanism of chiral three-nucleon force in driving spin-orbit splitting, Phys. Lett. B 855, 138839 (2024)

    work page 2024

  9. [9]

    Otsuka, T

    T. Otsuka, T. Suzuki, J. D. Holt, A. Schwenk, and Y. Akaishi, Three-Body Forces and the Limit of Oxygen Isotopes, Phys. Rev. Lett. 105, 032501 (2010)

    work page 2010

  10. [10]

    Sammarruca and R

    F. Sammarruca and R. Millerson, Nuclear Forces in the Medium: Insight From the Equation of State, Front. Phys. 7, 213 (2019)

    work page 2019

  11. [11]

    Wakasa, K

    T. Wakasa, K. Ogata, and T. Noro, Proton-induced knockout reactions with polarized and unpolarized beams, Prog. Part. Nucl. Phys. 96, 32 (2017)

    work page 2017

  12. [12]

    Minomo, M

    K. Minomo, M. Kohno, K. Yoshida, and K. Ogata, Prob- ing three-nucleon-force effects via (p, 2p) reactions, Phys. Rev. C 96, 024609 (2017)

    work page 2017

  13. [13]

    Uesaka, On the kinematics of ( p, pX) knockout re- actions in normal and inverse kinematics, Prog

    T. Uesaka, On the kinematics of ( p, pX) knockout re- actions in normal and inverse kinematics, Prog. Theor. Exp. Phys. 2024, 083D01 (2024)

    work page 2024

  14. [14]

    Chazono, K

    Y. Chazono, K. Yoshida, and K. Ogata, Importance of deuteron breakup in the deuteron knockout reaction, Phys. Rev. C 106, 064613 (2022)

    work page 2022

  15. [15]

    Wita la, W

    H. Wita la, W. Gl¨ ockle, D. H¨ uber, J. Golak, and H. Ka- mada, Cross section minima in elastic N d scattering: Possible evidence for three-nucleon force effects, Phys. Rev. Lett. 81, 1183 (1998)

    work page 1998

  16. [16]

    Wita la, J

    H. Wita la, J. Golak, and R. Skibi´ nski, Significance of chi- ral three-nucleon force contact terms for understanding of elastic nucleon-deuteron scattering, Phys. Rev. C 105, 054004 (2022)

    work page 2022

  17. [17]

    M. A. Franey and W. G. Love, Nucleon-nucleon t-matrix interaction for scattering at intermediate energies, Phys. Rev. C 31, 488 (1985)

    work page 1985

  18. [18]

    Levenberg, A method for the solution of certain non- linear problems in least squares, Quart

    K. Levenberg, A method for the solution of certain non- linear problems in least squares, Quart. Appl. Math. 2, 164 (1944)

    work page 1944

  19. [19]

    D. W. Marquardt, An algorithm for least-squares estima- tion of nonlinear parameters, J. Soc. Indust. Appl. Math. 11, 431 (1963)

    work page 1963

  20. [20]

    Ermisch, H

    K. Ermisch, H. R. Amir-Ahmadi, A. M. v. d. Berg, R. Castelijns, B. Davids, A. Deltuva, E. Epelbaum, W. Gl¨ ockle, J. Golak, M. N. Harakeh, M. Hunyadi, M. A. d. Huu, N. Kalantar-Nayestanaki, H. Kamada, M. Kiˇ s, M. Mahjour-Shafiei, A. Nogga, P. U. Sauer, R. Skibi´ nski, H. Wita la, and H. J. W¨ ortche, Systematic investigation of three-nucleon force effect...

    work page 2005

  21. [21]

    Hatanaka, Y

    K. Hatanaka, Y. Shimizu, D. Hirooka, J. Kamiya, Y. Ki- tamura, Y. Maeda, T. Noro, E. Obayashi, K. Sagara, T. Saito, H. Sakai, Y. Sakemi, K. Sekiguchi, A. Tamii, T. Wakasa, T. Yagita, K. Yako, H. P. Yoshida, V. P. La- dygin, H. Kamada, W. Gl¨ ockle, J. Golak, A. Nogga, and H. Wita la, Cross section and complete set of proton spin observables in ⃗ pdelastic...

    work page 2002