Pith. sign in

REVIEW 3 major objections 6 minor 16 references

A method of numerical calculation of the effect of short-range correlations for a wide range of nuclei

T0 review · 3 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A numerical method now computes short-range-correlation effects on proton densities and form factors for any closed-shell nucleus from helium-4 to lead-208.

desk verdict A useful numerical extension of the cluster expansion that overreaches in its abstract: the claimed Bethe–Goldstone validation is never actually shown. read the letter →

arxiv 2506.05512 v1 pith:YCG2S4GA submitted 2025-06-05 nucl-th

Tianyang Ma , Shalom Shlomo This is my paper
classification nucl-th
keywords short-rangecorrelationsprotondensityformfactorclusterexpansionJastrowcorrelationBethe-Goldstoneequationclosed-shellnucleiGaussianfunction
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper introduces a numerical method for computing how short-range nucleon–nucleon correlations change the proton density, root-mean-square radius, and elastic form factor of closed-shell nuclei. Earlier first-order cluster-expansion treatments of Jastrow correlations stalled at light nuclei because the angular integrals become intractable once orbitals with high angular momentum are occupied. The new method integrates the angular parts analytically, using Legendre expansions and the spherical-harmonic addition theorem, leaving only two-dimensional radial integrals that are easy to evaluate numerically for any spherical-symmetric mean-field potential. It also adds a Gaussian attractive term to the correlation function $g(r)$, not just the repulsive term. Applied to helium-4 through lead-208, the method with adjusted parameters yields SRC-modified proton densities similar to those obtained from the Bethe-Goldstone equation, the standard two-body scattering equation of nuclear many-body theory.

What carries the argument

The central object is the Jastrow correlation factor $f(r)$ and the associated density change $g(r) = f(r)^2 - 1$, taken here as the Gaussian $g(r) = -G e^{-\beta^2 r^2}$ and extended to a two-Gaussian form with an attractive part. The argument is carried by the first-order cluster expansion of the correlated one-body density, which splits the SRC correction into four integrals, $A$, $B$, $C$, and $D$, over the shell-model density matrix and total density. Three technical identities make those integrals tractable: the expansion of $\rho_0(r,r')$ into Legendre polynomials with radial coefficients $k_l(r,r')$; the product identity $P_l(x)P_m(x) = \sum_n a_n x^n$, which turns angular integrals into sums of exponentials divided by powers of the Gaussian width; and the spherical-harmonic addition theorem, which reduces the three-point term $D$ to a single angular integration. Together they convert a six-dimensional problem into two-dimensional radial quadrature.

What would settle it

Compute the term $D$ in Eq. (4d) by direct numerical integration in three dimensions, for example with Monte Carlo sampling, using the same Wood-Saxon wave functions and $g(r)$ used in the paper for $^{40}$Ca, and compare it with the closed-form angular reduction; any difference beyond integration error would show the reduction is wrong. Alternatively, measure the elastic electron scattering form factor of $^{208}$Pb at high momentum transfer and check whether the predicted SRC-modified diffraction minima move in the direction and magnitude the method gives.

Watch

Extended reading notes

Core claim

The paper's central claim is that the first-order cluster expansion of a Jastrow-correlated many-body wave function can be evaluated numerically for any doubly closed-shell nucleus, independent of the maximum orbital angular momentum $l$, and that including the attractive part of the N-N interaction in $g(r)$ gives SRC corrections consistent with Bethe-Goldstone results. The reduction writes the shell-model one-body density matrix as $\rho_0(r,r') = \sum_l k_l(r,r') P_l(\cos\theta)$, expands products of Legendre polynomials into powers of $x$, and evaluates $\int_{-1}^{1} P_l(x)P_m(x)e^{cx}\,dx$ in closed form. The three-point term $D$ is handled with the spherical-harmonic addition theorem and an orthogonality step that removes all but one angular mode. What remains is a set of radial double integrals, so the same code applies to $^4$He, $^{16}$O, $^{28}$Si, $^{32}$S, $^{40}$Ca, $^{60}$Ni, $^{90}$Zr, $^{140}$Ce, and $^{208}$Pb. With the two-Gaussian correlation function $g_{\rm sum}(r) = K e^{-\beta_1^2 r^2} - (G+K)e^{-\beta^2 r^2}$, the computed proton RMS radii shift by less than about $1.5\%$, with repulsion alone increasing the radius and the full repulsive-plus-attractive correlation decreasing it.

Load-bearing premise

The whole calculation depends on the step in which the angle between the position vectors $\mathbf{r}$ and $\mathbf{r}_2$ is treated as the sum of the angles between $\mathbf{r}$ and $\mathbf{r}_1$ and between $\mathbf{r}_1$ and $\mathbf{r}_2$; if that planar-geometry replacement is invalid for general three-dimensional configurations, every computed density and form factor inherits the error.

Editorial extensions

If this is right

  • Proton density, RMS radius, and elastic form factor can be computed with SRC included for any doubly magic nucleus with a spherical-symmetric mean field, including $^{208}$Pb.
  • With the paper's parameters, full SRC (repulsion plus attraction) lowers the proton RMS radius by about $0.18\%$ for $^{208}$Pb up to about $1.0\%$ for $^4$He, while repulsion alone raises it.
  • The method reproduces earlier harmonic-oscillator results for light nuclei when the same correlation parameters are used, confirming consistency with prior cluster-expansion calculations.
  • Because $g(r)$ is written as a sum of Gaussians, any short-range correlation that can be expanded in Gaussians, including realistic N-N interaction defect functions, can be fed into the same formalism.
  • Form factors computed from the SRC-modified densities can be compared directly with elastic electron scattering data within the first Born approximation.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same angular-integration reduction should transfer to momentum distributions and two-body densities, where short-range correlation effects are considerably larger than in the one-body density.
  • Treating the two Gaussian widths and strengths as fit parameters against realistic N-N interactions, rather than against final radii, would give the method predictive power for nuclei outside the closed-shell set.
  • Extending the calculation to neutron densities would yield SRC-modified neutron-skin thicknesses and Coulomb energy shifts, which are measurable and would provide independent tests of the correlation model.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper presents a numerical method for computing the first-order cluster-expansion correction to proton densities and form factors due to short-range correlations (SRC), using a Jastrow-type correlation function. The method is designed for arbitrary spherical single-particle potentials, and the authors apply it to nine closed-shell nuclei from 4He to 208Pb using a Woods-Saxon mean field. The correlation function is a sum of two Gaussians, one repulsive and one attractive, with four free parameters (G, beta, K, beta1) that are stated to be adjusted to reproduce Bethe-Goldstone results. The authors report RMS radii and density/form-factor plots for the nine nuclei and claim in the conclusion that the full SRC model yields results similar to the Bethe-Goldstone equation.

Significance. If the method is correct and properly validated, it would fill a useful niche: it extends the Jastrow cluster-expansion treatment of SRC to heavy nuclei with arbitrary mean-field potentials, including the attractive part of the N-N interaction, which earlier works using harmonic-oscillator or light-nucleus approximations did not cover. The consistency check against Ref. [5] for Ca40 is a positive internal validation, and the multipole reduction of the angular integrals, while obscurely presented, appears mathematically sound in its essential steps. However, the paper's central claim—that the parameters are tuned to reproduce Bethe-Goldstone results—is not demonstrated anywhere in the manuscript, and the printed correlation functions in Eqs. (27)-(28) are unphysical as written. The numerical results in Table 1 are therefore not yet predictive; they are a demonstration with an arbitrary parameter set. The potential significance is real, but the current manuscript does not establish it.

major comments (3)
  1. [Section 3 and Conclusions] The central claim that the SRC parameters are adjusted to reproduce Bethe-Goldstone results is not supported by any data in the paper. Section 3 presents densities, form factors, and RMS radii for nine nuclei, but it never shows a Bethe-Goldstone density, radius, or form factor for any of them, and it does not describe any fitting procedure for G, beta, K, and beta1. Table 1 compares only shell-model and SRC models among themselves. The concluding sentence that the full model 'can yield a result similar to the one provided by Bethe-Goldstone Equation' therefore rests on no evidence presented in the manuscript. This is the load-bearing validation for the parameter choice and for the claimed usefulness of the method; it must be added before the central claim can be assessed.
  2. [Eqs. (27)-(28) and Section 2] The correlation functions printed in Eqs. (27) and (28), g_attraction(r) = e^{(2.4r)^2} - e^{-(4.0r)^2} and g_sum(r) = e^{(2.4r)^2} - 2e^{-(4.0r)^2}, contain positive exponents and diverge as r goes to infinity. That violates requirement (ii) stated in Section 2 that g_attraction must vanish rapidly at large r, and it would make the integrals in Eq. (18) divergent. Since Eq. (16) defines g_attraction with a negative exponent, the printed expressions are likely sign typos, but as written they describe an unphysical model. The authors must correct these equations and confirm that all reported results were obtained with e^{-(2.4r)^2}, not e^{+(2.4r)^2}.
  3. [Appendix A, Eq. (4d*)] The derivation of term D uses the notation P_{l3}(cos(theta+theta2)) as if the angle between r and r2 were the sum of the angles between r and r1 and r1 and r2. That is not generally true for three arbitrary vectors, and the text's explanation of the subsequent reduction is too terse. The underlying multipole reduction is, in fact, valid: after expanding in spherical harmonics and integrating over the orientation of r1, only l3 = l and m = 0 survive, leading to Eq. (A8). However, the manuscript's presentation obscures this and should be rewritten to use the spherical-harmonic addition theorem explicitly, both to remove the apparent planar-geometry assumption and to make the orthogonality step clear.
minor comments (6)
  1. [Eq. (12) and throughout] The word 'polynormal' should be 'polynomial' in Eq. (12) and in Appendix A.
  2. [Eq. (14) and Eq. (18)] The typesetting of these multi-line equations is extremely hard to read; the sums over k and n are interleaved ambiguously with the integrands. The authors should reformat the equations with proper brackets, alignment, and explicit ranges so that a reader can verify the terms.
  3. [Eqs. (22)-(25)] The potential parameters are given as V0 = 50(1+0.72(N-Z)/A) MeV and then written again in Eqs. (23) and (25); this is redundant and the notation should be made consistent. Also, 'Wood-Saxon' should be 'Woods-Saxon'.
  4. [Table 1] The column header 'Nucleon' should be 'Nucleus'; the table would also benefit from an explicit statement of which parameter set was used and where the Bethe-Goldstone comparison is (or is not) shown.
  5. [References] The conclusion cites Ref. [10] (S. Shlomo, 'Nuclear Coulomb Energies') as the source of Bethe-Goldstone results, but that reference does not appear to provide the densities or form factors needed for the claimed comparison; a direct reference to Bethe-Goldstone calculations (e.g., Ref. [11]) should be used, and the actual comparison should be shown.
  6. [Figure captions] The phrase 'British Flag shaped dots' is unclear; please use standard marker names (e.g., 'asterisks', 'crosses', 'diamonds').

Circularity Check

1 steps flagged · score 6.0 of 10

The four SRC parameters are adjusted to reproduce Bethe-Goldstone results, so the concluding agreement with Bethe-Goldstone is a restatement of that fit rather than an independent prediction.

  1. fitted input called prediction [Abstract; Eq. (13); Eq. (16); Section 3 Results; Table 1; Section 4 Conclusions]
    "In the short-range correlations, we have included the effects of repulsive and attractive parts of the N-N interaction and adjusted the SRC parameters to reproduce the results obtained by solving the Bethe-Goldstone Equation. ... with a specific set of parameters that include both the repulsive and attractive part of the N-N interaction can yield a result similar to the one provided by Bethe-Goldstone Equation[10]."

    Eq. (13) introduces g(r) with free parameters G and beta, and Eq. (16) introduces g_attraction with free parameters K and beta1; all four are adjustable. The abstract states these parameters are adjusted to reproduce Bethe-Goldstone results, and the conclusion then offers agreement with the Bethe-Goldstone Equation as the outcome of the model. Thus the 'similar' result is guaranteed up to fitting error by the calibration procedure; it is not an independent numerical prediction. No Bethe-Goldstone density, RMS radius, or form factor is shown in Section 3 or Table 1, so the claimed agreement cannot be checked against an external curve. The Appendix B check against Ref. [5] is a genuine external benchmark, but the full-repulsion-plus-attraction model is only compared with its own fit target.

full rationale

The paper's core technical content, the numerical evaluation of the first-order cluster expansion for spherical shell-model wave functions, is not circular: Eq. (14) and Eq. (18) are obtained by substituting the chosen Gaussian g(r) into the cluster-expansion density formula adopted from Ref. [5], and Appendix B reproduces the density and form factor of Ref. [5] as an external benchmark. The terse angular reduction in Appendix A is an algebraic integration step, not a circular use of the target result, so no circularity is established there. The single clear circular step is the calibration of the four correlation parameters (G, beta, K, beta1) to Bethe-Goldstone results and the subsequent presentation of agreement with Bethe-Goldstone as a conclusion of the work. Because the agreement is the direct output of the fit, the central validation claim reduces by construction to the fit input, giving partial circularity. The paper is transparent about the adjustment, which is why the overall score is 6 rather than higher: the method itself retains independent content, but the headline Bethe-Goldstone consistency is not an independent test.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central claim rests on four adjustable parameters (G, beta, K, beta1) that are fitted to Bethe-Goldstone results, plus the validity of the first-order cluster expansion and a specific questionable angular reduction for the D term. No new physical entities are introduced.

free parameters (4)
  • G = 1
    Strength of the repulsive SRC correlation factor in Eq. (13); adjusted to reproduce Bethe-Goldstone results (Section 3).
  • beta = 4.0 fm^-1
    Range parameter of the repulsive Gaussian in Eq. (13); adjusted to reproduce Bethe-Goldstone results (Section 3).
  • K = 1
    Strength of the attractive SRC component in Eq. (16); adjusted together with beta1 to reproduce Bethe-Goldstone results (Section 3).
  • beta1 = 2.4 fm^-1
    Range parameter of the attractive Gaussian in Eq. (16); adjusted to reproduce Bethe-Goldstone results (Section 3).
assumptions (4)
  • domain assumption The cluster expansion truncated at first order in g is a valid approximation for all nuclei considered.
    Used to derive Eq. (3) from the Jastrow wave function; no convergence check is provided for A>40.
  • ad hoc to paper The correlation factor g(r) is a sum of two Gaussians, Eqs. (13) and (16).
    The Gaussian form is chosen for analytic tractability; it is fitted to Bethe-Goldstone results, not derived from a realistic N-N interaction.
  • ad hoc to paper The angular integration in term D (Eq. 4d) can be reduced using orthogonality after integrating over d r1, yielding Eq. (A8).
    The addition theorem is applied with the angle between r and r2 replaced by theta+theta2, which is not generally valid; this step is not rigorously justified in Appendix A.
  • domain assumption Woods-Saxon parameters from global systematics (Eqs. 22-25) provide adequate single-particle wave functions.
    Standard parameterization; no comparison with experimental single-particle energies or radii is provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A method of numerical calculation of the effect of short-range correlations for a wide range of nuclei." pith.science (2026). https://pith.science/paper/YCG2S4GA

@misc{pith2026250605512,
  author       = {Pith},
  title        = {Pith review of: A method of numerical calculation of the effect of short-range correlations for a wide range of nuclei},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YCG2S4GA}},
  note         = {Machine review of arXiv:2506.05512}
}
abstract

We introduce a method to calculate the effect of short-range correlation (SRC) on the proton density distribution numerically, up to the first order of the cluster expansion, that can be used for wide range of closed shell nuclei and determine the effect on proton density, root-mean square (RMS) radius, and form factor of many different closed-shell nuclei: $^4$He, $^{16}$O, $^{28}$Si, $^{32}$S, $^{40}$Ca, $^{60}$Ni, $^{90}$Zr, $^{140}$Ce and $^{208}$Pb. In the short-range correlations, we have included the effects of repulsive and attractive parts of the N-N interaction and adjusted the SRC parameters to reproduce the results obtained by solving the Bethe-Goldstone Equation.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

  1. [5]

    The parameters used in [5] is like this: 23 V= ℏ2α4 2m r2 wher α=0.55fm−1, and the correlation factor they used is β=1.4fm−1, g(r)=-𝑒−(1.4𝑟)2

    to check the consistency. The parameters used in [5] is like this: 23 V= ℏ2α4 2m r2 wher α=0.55fm−1, and the correlation factor they used is β=1.4fm−1, g(r)=-𝑒−(1.4𝑟)2 . Although [5] did not give RMS radius results, it provided the proton density and form factors derived from such parameters as the Fig.(9) and Fig.(8b) of [5]. And the proton density and f...

    work page

  2. [1]

    Hence, we need to consider the effects that mean field methods omit

    Introduction Traditionally, mean-field methods, such as the Shell Model, have been used to describe the nuclear properties, but since the 1980s experiments like elastic scattering experiments that measure charge distributions [1] and knockout experiments that measure spectral functions [2] have shown that mean-field methods have difficulty in describing t...

    work page

  3. [2]

    Formalism 5 Under the independent particle approximation, the shell model many-body wave function is given by the Slater determinant of the occupied single-particle wave functions 𝜓𝑖(𝒓𝒋), ψsm = N √A! det(ψ1(𝐫𝟏) … ψA(𝐫𝐀)), (1) which assumes that all the nucleons are only subject to a mean field. To account for the short-range correlation, we use the Jastro...

    work page

  4. [3]

    Results We show that in Appendix B that, the method can calculate the nucleon densities for light nuclei with results consistent with preexisting works see also [5], but unlike those works, our method can be applied to heavier nuclei with any spherical- symmetric shell model potentials. When applied to heavier nuclei, the effect of short-range correlation...

    work page

  5. [4]

    Conclusions In section 2, we derived a method that can be used to calculate the effect of Short- range Correlation (SRC) on proton density and form factor, from light to heavy nuclei, numerically for any spherical symmetric nuclei. In section 3 we have applied that method to several nuclei that under shell model should have closed shells: He 4 , O 16 , Si...

    work page

  6. [6]

    Is the Shell-Model Concept Relevant for the Nuclear Interior?

    J. M. Cavedon et al, “Is the Shell-Model Concept Relevant for the Nuclear Interior?” Phys. Rev. Lett. 49, 978 (1982)

    work page 1982

  7. [7]

    Relative 3𝑠 Spectroscopic Strength in Pb 206 and Pb 208 Studied with the (𝑒,𝑒′𝑝) Knockout Reaction

    E. N. M. Quint et al, “Relative 3𝑠 Spectroscopic Strength in Pb 206 and Pb 208 Studied with the (𝑒,𝑒′𝑝) Knockout Reaction.” Phys. Rev. Lett. 57, 186 (1986)

    work page 1986

  8. [8]

    Many-Body Problem with Strong Forces

    R. Jastrow, “Many-Body Problem with Strong Forces.” Phys. Rev. 98, 1479 (1955)

    work page 1955

Show all 16 references
  1. [9]

    Cluster Development Method in the Quantum Mechanics of Many Particle System

    Fumiaki Iwamoto, Masami Yamada, “Cluster Development Method in the Quantum Mechanics of Many Particle System.” Progress of Theoretical Physics, 17, 543(1957)

    1957

  2. [10]

    Jastrow correlations

    M. Gaudin, J. Gillespie, and G. Ripka, “Jastrow correlations.” Nucl. Phys. A, 176, 237 (1971)

    1971

  3. [11]

    Stringari,”Effects of short-range correlations on nuclear mass and momentum distributions.” Phys

    Bohigas and S. Stringari,”Effects of short-range correlations on nuclear mass and momentum distributions.” Phys. Lett., 95, 9 (1980)

    1980

  4. [12]

    Geminal wave functions with Jastrow correlation: A first application to atoms

    J. Chem, “Geminal wave functions with Jastrow correlation: A first application to atoms.” Phys. 119, 6500 (2003). 25

    2003

  5. [13]

    Information on Short-Range Nucleon- Nucleon Correlations and Woods-Saxon Nuclear Potential from High-Energy Elastic Electron Scattering

    C. Ciofi, Degli Atti, and N. M. Kabachnik, “Information on Short-Range Nucleon- Nucleon Correlations and Woods-Saxon Nuclear Potential from High-Energy Elastic Electron Scattering.” Phys. Rev. C 1, 809 (1970)

    1970

  6. [14]

    Particles and Nuclei: An Introduction to the Physical Concepts

    B. Povh, K. Rith, C. Scholz, and F . Zetsche, “Particles and Nuclei: An Introduction to the Physical Concepts.” Berlin: Springer-Verlag. 73. ISBN 978-3-540-43823-6

  7. [15]

    Nuclear Coulomb Energies

    S.Shlomo, “Nuclear Coulomb Energies.” Rep. Prog. Phys. 41, 957 (1978)

    1978

  8. [16]

    Solution of the Bethe-Goldstone equation

    G. Dahll, E. Østgaard, B. Brandow, “Solution of the Bethe-Goldstone equation.” Nucl. Phys. A, 124, 481 (1969). Table list: Table 1: The theoretical proton root mean square (RMS) radius obtained by our method of calculating short range correlation (SRC) effects with the paramet...

    1969

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.