REVIEW 3 major objections 6 minor 52 references
np spin correlations in the deuteron ground state
T0 review · 3 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read The deuteron's proton and neutron carry near-maximal spin entanglement in the M=0 channel, with negativity close to the two-qubit bound.
desk verdict Solid fixed-projection numbers, but the zero-field entanglement claim is phase-dependent and rests on an arbitrary coherent superposition. 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 machinery is the Clebsch-Gordan decoupling of the coupled deuteron state into the uncoupled two-nucleon spin basis, followed by a partial trace over the radial, orbital, and isospin degrees of freedom. That yields the 4x4 reduced spin density matrix $\hat{\rho}_{\rm spin}$, which is then scored by two standard bipartite measures: mutual information $I(A:B)$ (Eq. B3) and entanglement negativity $N(\hat{\rho}_{\rm AB})$ (Eq. B4). For the superposition case, the additional input is Eq. (7), an equal-probability coherent superposition ansatz over M=0,±1 with free phase parameters $\zeta_{\pm 1}$; the phases become the control knobs that govern how much entanglement survives.
What would settle it
Measure the spin state of an unpolarized, field-free deuteron ensemble by full tomography of the proton-neutron spin correlations: if the observed mutual information and negativity match the incoherent-mixture averages over M rather than the coherent-equal-superposition values of Table II, the equal-prior coherent ansatz is falsified. A second check is to vary the L=2 admixture, for instance in a dineutron or a different two-nucleon channel, and test whether M=0 negativity follows the monotone curve in the right panel of Fig. 1.
Extended reading notes
Core claim
On the paper's own terms, the deuteron ground state, written as an S-D admixture $|\Psi^{(M)}\rangle = \sum_L \alpha_L |R_L\rangle |LS;JM\rangle |T M_T\rangle$, yields after partial tracing a spin density matrix $\hat{\rho}_{\rm spin}^{(M)}$ whose entanglement is largest for M=0 and essentially vanishing for M=±1. For the six potentials listed in Table I, the M=0 mutual information lies between 1.1829 and 1.2671 and the negativity between 0.4580 and 0.4784, the latter approaching the 0.5 ceiling for two-qubit pure states. Extending to the field-free case by Eq. (7), an equal coherent superposition $\frac{1}{\sqrt{3}}\sum_M e^{i\zeta_M} |\Psi^{(M)}\rangle$, the neutron-proton spin state remains strongly entangled; the entanglement is maximal when the phases of the M=+1 and M=-1 components are complementary and minimal when they are equal, as quantified in Table II and Figs. 2-3.
Load-bearing premise
The result depends on Eq. (7): that in the absence of a magnetic field the deuteron is an equal coherent superposition of the three degenerate M projections with arbitrary phases; if the real field-free deuteron is instead an incoherent mixture of those projections, the superposition-dependent entanglement numbers in Table II and Figs. 2-3 do not follow.
Editorial extensions
If this is right
- Because the M=0 spin state sits close to the two-qubit entanglement ceiling across all six potential models, the deuteron offers a concrete small-system benchmark in which nuclear-structure predictions and entanglement measures can be compared directly.
- In the superposition case, the entanglement is controlled by a single phase parameter $\theta_{+1}$ after the redundant combination is removed, so measuring spin correlations in a field-free deuteron would probe the relative coherence of the degenerate angular-momentum projections.
- The M=±1 channels contribute essentially no negativity, so spin entanglement in the deuteron is concentrated in the M=0 (spin-singlet-like) component; experiments that separate M components should therefore see the entanglement signal only in that channel.
- The same reduced-density-matrix construction can be applied to other two-body nuclear systems where the S-D admixture differs, giving an entanglement-versus-$|\alpha_2|^2$ curve that could serve as a model discriminator.
Reading between the lines
- If the field-free state is instead described as a thermal ensemble over M, the strong phase-dependent entanglement of Table II would collapse to a much weaker average; the paper's central quantitative claim for the superposition case therefore hinges on a preparation assumption that a measurement could distinguish.
- The near-maximal M=0 negativity implies a testable correlation witness: a polarized deuteron beam prepared in M=0 should exhibit proton-neutron spin correlations that exceed the classical bound, which could be probed in scattering or breakup experiments.
- The same Clebsch-Gordan decoupling plus partial-trace construction could be applied to the dineutron continuum or to heavier two-cluster systems, mapping out how entanglement tracks the orbital admixture rather than the binding energy.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives the reduced 4×4 spin density matrix for the proton-neutron spins in the deuteron ground state, first for a fixed total-angular-momentum projection M and then for an equal coherent superposition of M=0, ±1 (Eq. (7)). Using D-state probabilities |α_2|^2 taken from six potential models, it reports the mutual information and entanglement negativity of the reduced spin state. For fixed M it finds that the M=0 component is nearly maximally entangled (negativity ≈0.46–0.48), while M=±1 have negligible negativity. For the field-free superposition, it reports a strong phase dependence and, for complementary phases, values identical to the M=0 case; the abstract nevertheless states unconditionally that strong entanglement persists when all spin projections are superposed.
Significance. If the field-free superposition claim were justified, the paper would provide a clean quantitative link between nuclear structure and quantum information measures in the simplest nucleus. The fixed-M calculation is standard and internally consistent, and the numerical values for six independent potentials are plausible; importantly, no parameter is fitted to the entanglement measures themselves, and the M=0 result is robust across all six potentials. The limitation is that the superposition result is phase-dependent and rests on an unjustified pure-state preparation, so the headline field-free claim is not a property of the deuteron ground state as such.
major comments (3)
- [II, Eq. (7)] The 'principle of equal prior probabilities' cannot justify the coherent superposition |Ψ> = (1/√3) Σ_M e^{iζ_M} |Ψ(M)>. For degenerate M levels in the absence of a magnetic field, equal a priori probabilities describe a statistical mixture ρ_mix = (1/3) Σ_M |Ψ(M)><Ψ(M)|, not a pure state with arbitrary phases. The phases ζ_{±1} are free state-preparation parameters, and their values are not selected by any property of the deuteron. Since Table II and Fig. 2 show that I(A:B) and N vary by roughly a factor of 5 and 3 between the complementary-phase and equal-phase branches, the abstract's unconditional statement that 'strong entanglement still exists even when all spin states are superposed' is phase-selected and requires either a physical preparation mechanism or an explicit conditionalization.
- [III, Table II] For every potential, the ζ_{+1} ≠ ζ_{-1} columns of Table II reproduce the M=0 columns of Table I to four decimals. Thus the 'maximum correlation occurs when the phases are complementary' branch is not a genuinely mixed-M superposition at the level of the spin density matrix; it reduces to the already reported fixed-M=0 state. The branch that actually represents an equal-weight superposition (ζ_{+1}=ζ_{-1}) gives I(A:B) ≈ 0.26–0.28 and N ≈ 0.15; for the Woods-Saxon row these are 0.2781 and 0.1468, compared with 1.1921 and 0.4604 for M=0. The authors should compute and report the entanglement of the rotationally invariant incoherent mixture ρ_mix, or explicitly restrict the field-free claim to a stated phase preparation.
- [IV and Abstract] The summary statement that in the degenerate field-free case 'the nucleons are close to maximally entangled in their spins' is not supported for the equal-phase branch, where the negativity is roughly 0.15, far below the two-qubit maximum of 0.5. The abstract and conclusions should be reworded so that 'strong entanglement' refers either to the fixed M=0 projection or to a superposition with complementary phases, not to the generic superposition of all projections.
minor comments (6)
- [I] The Introduction contains several typos ('The is becuse of', 'experimently', 'resulant') and would benefit from a careful proofread.
- [III] The text repeats 'we shall later study the case' and uses 'momentum projections' where 'angular momentum projections' is meant.
- [IV] The Summary calls the M=0 component a 'spin singlet state'; for S=1, M_S=0 the correct term is a triplet Bell state, not a spin singlet.
- [II, Eq. (7)] The notation ζ_M ∈ (0,2π) should be [0,2π), and the text should state explicitly that a global phase has been fixed so only the relative phase ζ_{+1}−ζ_{-1} is physical.
- [III] The sentence about the 'absence of closed contours' and 'both phases equally favored' is unclear; the authors should state directly that the correlation appears to depend only on the relative phase ζ_{+1}−ζ_{-1}.
- [Table I] The caption says 'These calculation are for fixed value of projection quantum number M'; this should be corrected grammatically, and the method used to extract |α_2|^2 from the cited potential models should be described or referenced precisely.
Circularity Check
Superposition claim partly re-labels the M=0 result; the fixed-M analysis remains independent and self-contained.
-
renaming known result
[Abstract; Eq. (7); Table I vs Table II]
"Our findings show that the spins are most entangled when the total projection is zero, and that strong entanglement still exists even when all spin states are superposed. ... TABLE I: ... These calculation are for fixed value of projection quantum number M. ... TABLE II: ... the entanglement quantifiers for different potential models is calculated for the case considering all possible projection of J as equally probable. ζ±1 = 0,π representing the sign of coefficient in Eq.(8) to ±1."
For the ζ+1≠ζ−1 branch, every entry in Table II exactly matches the M=0 column of Table I (e.g., NLO: I=1.2671, N=0.4784 in both tables), not merely approximately. Since Table I is for a fixed projection M and Table II is for an equal coherent superposition of all M, the numerical identity means the complementary-phase superposition has the same reduced spin density matrix as the M=0 fixed-projection state. The abstract's unconditional 'strong entanglement still exists even when all spin states are superposed' therefore reports the already-computed M=0 result under the label 'superposition'; it is not an independent consequence of superposing the three projections.
full rationale
The fixed-M calculation (Table I) is self-contained: it uses |α2|2 values from six independent potential models, applies the standard J=1, S=1, T=0 deuteron wavefunction, traces out orbital and isospin degrees of freedom, and computes mutual information and negativity. No parameter is fitted to the entanglement measures, and no load-bearing self-citation supplies the result. The near-maximal M=0 entanglement follows from the triplet Bell-like structure and is not circular. The superposition section, however, contains one circular/renaming step: the maximal 'all-M superposed' column (ζ+1≠ζ−1) is numerically identical to the M=0 fixed-projection column, so the abstract's broad claim about strong entanglement under superposition inherits the M=0 input by construction. Because the paper explicitly discloses the phase dependence in Table II and Figs. 2–3, this is a partial overstatement and a re-labeling rather than a hidden fitting of parameters; the independent fixed-M content keeps the overall circularity moderate.
Assumptions & free parameters
free parameters (2)
- D-state probability |alpha_2|^2 =
0.0362 to 0.0702 (six potential models); scanned 0 to 1 in Fig. 1
- Relative phases zeta_+1, zeta_-1 (or theta_+1) =
scanned over [0, 2pi); Table II uses 0 and pi
assumptions (4)
- domain assumption The deuteron ground state wavefunction has the form of Eq. (1): pure state with S=1, T=0, J=1, L=0,2, with probabilities |alpha_L|^2 from potential models.
- standard math Partial tracing over spatial, orbital, and isospin degrees yields a spin density matrix via angular momentum orthogonality (Eqs. 5-6).
- ad hoc to paper In the absence of a magnetic field, the deuteron is in an equal coherent superposition of all M projections (Eq. 7).
- domain assumption Nucleon spins can be treated as subsystems of distinguishable particles.
Cite this review
Pith. "Pith review of np spin correlations in the deuteron ground state." pith.science (2026). https://pith.science/paper/6CEDS33Z
@misc{pith2026250616621,
author = {Pith},
title = {Pith review of: np spin correlations in the deuteron ground state},
year = {2026},
howpublished = {\url{https://pith.science/paper/6CEDS33Z}},
note = {Machine review of arXiv:2506.16621}
}
read the original abstract
The deuteron is the simplest atomic nucleus made of two particles - a proton and a neutron. In this work, we study how their spins are quantum entangled with each other. We study two cases: when the deuteron is in a fixed projection of total angular momentum, and when it exists in a superposition of all projections. Our findings show that the spins are most entangled when the total projection is zero, and that strong entanglement still exists even when all spin states are superposed.
Figures
Reference graph
Works this paper leans on
-
[27]
A. T. Kruppa, J. Kovács, P. Salamon, and O. Legeza, Journal of Physics G: Nuclear and Particle Physics48, 025107 (2021)
work page 2021
-
[1]
Amico, R
L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. Phys.80, 517 (2008)
2008
-
[2]
A. K. Ekert, Phys. Rev. Lett.67, 661 (1991)
1991
-
[3]
R. Jozsa and N. Linden, Proceedings of the Royal So- ciety of London. Series A. Mathematical, Physical and Engineering Sciences459, 2011 (2003)
work page 2003
-
[4]
Vidal, Phys
G. Vidal, Phys. Rev. Lett.91, 147902 (2003)
2003
-
[5]
F. W. Strauch, American Journal of Physics84, 495 (2016), https://pubs.aip.org/aapt/ajp/article- pdf/84/7/495/13122295/495_1_online.pdf
work page 2016
-
[6]
S. Szalay, M. Pfeffer, V. Murg, G. Barcza, F. Ver- straete, R. Schneider, and O. Legeza, International Journal of Quantum Chemistry115, 1342 (2015), https://onlinelibrary.wiley.com/doi/pdf/10.1002/qua.24898
-
[7]
J. S. Dehesa, T. Koga, R. J. Yáñez, A. R. Plastino, and R. O. Esquivel, Journal of Physics B: Atomic, Molecular and Optical Physics45, 015504 (2011)
work page 2011
Show all 52 references
-
[8]
Qvarfort, S
S. Qvarfort, S. Bose, and A. Serafini, New Journal of Physics22, 093062 (2020)
2020
-
[9]
Eisert, M
J. Eisert, M. Cramer, and M. B. Plenio, Rev. Mod. Phys. 82, 277 (2010)
2010
-
[10]
Laflorencie, Physics Reports646, 1 (2016), quantum entanglement in condensed matter systems
N. Laflorencie, Physics Reports646, 1 (2016), quantum entanglement in condensed matter systems
2016
-
[11]
Horodecki, P
R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, Rev. Mod. Phys.81, 865 (2009)
2009
-
[12]
Horodecki, Physics Letters A232, 333 (1997)
P. Horodecki, Physics Letters A232, 333 (1997)
1997
-
[13]
Peres, Phys
A. Peres, Phys. Rev. Lett.77, 1413 (1996)
1996
-
[14]
Vedral, M
V. Vedral, M. B. Plenio, M. A. Rippin, and P. L. Knight, Phys. Rev. Lett.78, 2275 (1997)
1997
-
[15]
S. A. Hill and W. K. Wootters, Phys. Rev. Lett.78, 5022 (1997)
1997
-
[16]
Schliemann, J
J. Schliemann, J. I. Cirac, M. Kuś, M. Lewenstein, and D. Loss, Phys. Rev. A64, 022303 (2001)
2001
-
[17]
Eckert, J
K. Eckert, J. Schliemann, D. Bruß, and M. Lewenstein, Annals of Physics299, 88 (2002)
2002
-
[18]
Debarba, R
T. Debarba, R. O. Vianna, and F. Iemini, Phys. Rev. A 95, 022325 (2017)
2017
-
[19]
Gigena and R
N. Gigena and R. Rossignoli, Phys. Rev. A94, 042315 (2016)
2016
-
[20]
Gigena and R
N. Gigena and R. Rossignoli, Phys. Rev. A95, 062320 (2017)
2017
-
[21]
Shapourian and S
H. Shapourian and S. Ryu, Phys. Rev. A99, 022310 (2019)
2019
-
[22]
Benatti, R
F. Benatti, R. Floreanini, and U. Marzolino, Phys. Rev. A89, 032326 (2014)
2014
-
[23]
Friis, A
N. Friis, A. R. Lee, and D. E. Bruschi, Phys. Rev. A87, 022338 (2013)
2013
-
[24]
Legeza, L
O. Legeza, L. Veis, A. Poves, and J. Dukelsky, Phys. Rev. C92, 051303 (2015)
2015
-
[25]
Di Tullio, R
M. Di Tullio, R. Rossignoli, M. Cerezo, and N. Gigena, Phys. Rev. A100, 062104 (2019)
2019
-
[26]
J. Faba, V. Martín, and L. Robledo, Phys. Rev. A104, 032428 (2021)
2021
-
[28]
Robin, M
C. Robin, M. J. Savage, and N. Pillet, Phys. Rev. C103, 034325 (2021)
2021
-
[29]
A. T. Kruppa, J. Kovács, P. Salamon, O. Legeza, and G. Zaránd, Phys. Rev. C106, 024303 (2022)
2022
-
[30]
Lacroix, A
D. Lacroix, A. B. Balantekin, M. J. Cervia, A. V. Patward- han, and P. Siwach, Phys. Rev. D106, 123006 (2022)
2022
-
[31]
Pazy, Phys
E. Pazy, Phys. Rev. C107, 054308 (2023)
2023
-
[32]
C. Gu, Z. H. Sun, G. Hagen, and T. Papenbrock, Phys. Rev. C108, 054309 (2023)
2023
-
[33]
Bulgac, M
A. Bulgac, M. Kafker, and I. Abdurrahman, Phys. Rev. C107, 044318 (2023)
2023
-
[34]
C. W. Johnson and O. C. Gorton, Journal of Physics G: Nuclear and Particle Physics50, 045110 (2023)
2023
-
[35]
Bai, Phys
D. Bai, Phys. Rev. C109, 034001 (2024)
2024
- [36]
-
[37]
J. M. Blatt and V. F. Weisskopf,Theoretical nuclear physics(Springer, New York, 1952)
1952
-
[38]
Levchuk and A
M. Levchuk and A. L’vov, Nuclear Physics A674, 449 (2000)
2000
-
[39]
Fujiwara, T
Y. Fujiwara, T. Fujita, M. Kohno, C. Nakamoto, and Y. Suzuki, Phys. Rev. C65, 014002 (2001)
2001
-
[40]
Gross and A
F. Gross and A. Stadler, Phys. Rev. C82, 034004 (2010)
2010
-
[41]
Veerasamy and W
S. Veerasamy and W. N. Polyzou, Phys. Rev. C84, 034003 (2011)
2011
-
[42]
Epelbaum, W
E. Epelbaum, W. Glöckle, and U.-G. Meißner, Nuclear Physics A671, 295 (2000)
2000
-
[43]
Arenhövel and H
H. Arenhövel and H. G. Miller, Zeitschrift für Physik266, 13 (1974)
1974
-
[44]
R. B. Wiringa, R. A. Smith, and T. L. Ainsworth, Phys. Rev. C29, 1207 (1984)
1984
-
[45]
Reid, Roderick V., Annals of Physics50, 411 (1968)
J. Reid, Roderick V., Annals of Physics50, 411 (1968)
1968
-
[46]
Rezaei and A
B. Rezaei and A. Dashtimoghadam, Journal of Theoretical and Applied Physics8, 203 (2014)
2014
-
[47]
Hamada and I
T. Hamada and I. Johnston, Nuclear Physics34, 382 (1962)
1962
-
[48]
C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. Wootters, Phys. Rev. A54, 3824 (1996). 6
1996
-
[49]
Życzkowski, P
K. Życzkowski, P. Horodecki, A. Sanpera, and M. Lewen- stein, Phys. Rev. A58, 883 (1998)
1998
-
[50]
J. Lee, M. S. Kim, Y. J. Park, and S. Lee, Journal of Modern Optics47, 2151 (2000), https://doi.org/10.1080/09500340008235138
2000 doi
-
[51]
T. A. Palomaki, J. D. Teufel, R. W. Simmonds, and K. W. Lehnert, Science342, 710 (2013), https://www.science.org/doi/pdf/10.1126/science.1244563
2013 doi
-
[52]
Vidal and R
G. Vidal and R. F. Werner, Phys. Rev. A65, 032314 (2002). 7 APPENDICES npspin correlations in the deuteron ground state Appendix A: Spin density matrices When the system assumes a fixed projectionM, the reduced density matrix encompassingnp spin correlations is given by ˆρ(M) ...
2002
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.