{"id":"4b9ab5cb-3e00-4e1a-9ff9-4ab5bc0d18f4","arxiv_id":"2411.17462","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Lattice chiral EFT simulations predict that 22Si is bound against two-proton emission and that its protons and neutrons fill the Z=14 and N=8 shells.","lead":"Using lattice simulations with chiral nuclear forces, the authors predict that the proton-rich nucleus 22Si is bound against breaking into 20Mg plus two protons, and that it has closed shells at Z=14 and N=8. Nonspecialists may care because this bears on where the proton dripline sits, whether two-proton radioactivity exists near silicon, and how mirror-nucleus radii can constrain the symmetry energy of neutron-rich matter.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"HO-basis pinhole occupation numbers at tau = 0.2 MeV^-1 are not demonstrated to be ground-state observables; the doubly magic conclusion rests on an unvalidated, single-time measurement.","rationale":"The reader's conditional verdict identifies the HO-basis pinhole method as the weakest assumption; I agree. The strongest claim is composite: 22Si is bound against 2p emission and doubly magic. The binding-energy part (S2p = 0.72(48) MeV) is only 1.5 sigma above zero and depends on the cancellation of systematic errors between 22Si and 20Mg; using the experimental 20Mg mass would reverse the sign. However, the paper's use of the same Hamiltonian for both nuclei is a standard way to cancel correlated errors, and the two interactions agree, so that part is plausible but not definitive. The doubly magic part is weaker because the occupation numbers are not extrapolated in Euclidean time, the method is new and unvalidated, and the short-time result could reflect the trial state rather than the true ground state. Negative occupation numbers in Table S1 confirm the presence of systematic or statistical problems. A tau-convergence test and a benchmark on a closed-shell nucleus would settle whether the shell-closure conclusion holds. Since this is a necessary condition for the paper's headline, the paper should remain CONDITIONAL until such evidence is provided.","tokens_in":14689,"tokens_out":7647,"duration_ms":77338,"concrete_test":"Re-run the HO-basis pinhole occupation calculation for 22Si and for a known doubly magic nucleus such as 16O at several Euclidean times (tau = 0.1, 0.2, 0.3, 0.4, 0.5 MeV^-1) with the same truncation and Hamiltonian. If the occupations shift by more than the statistical errors between tau = 0.2 and tau = 0.4, or if 16O does not reproduce its expected closed-shell occupancies (e.g., 0d5/2 close to 1, 0d3/2 close to 0), then the single-tau results for 22Si are not reliable ground-state observables and the shell-closure claim should be qualified or withdrawn.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that 22Si is doubly magic (Z=14, N=8) rests on the new HO-basis pinhole method. The occupation numbers in Table S1 are computed at a single Euclidean time tau = 0.2 MeV^-1 (Methods), while energies and radii are extrapolated in tau. At this tau, the ground-state energy in Fig. 1 is still about 4 MeV above its extrapolated value, so the evolved state is not yet the ground state. Because the trial wavefunction is a product of HO single-particle states, the short-time HO-pinhole amplitude (S6) naturally returns the trial occupations: 0d5/2 filled and 0d3/2/1s1/2 empty, which is what Table S1 shows. The Supplemental Material concedes 'we did not apply translations to the wavefunctions' and estimates the CoM effect only for the energy (1.3% at tau = 0.3), not for the occupations. There is no benchmark on a known closed-shell nucleus (e.g., 16O or 40Ca), and Table S1 contains unphysical negative occupations (e.g., -0.127(169) for 0d3/2[3/2]). Hence the doubly magic conclusion is not independently established; the pinhole occupations could simply mirror the initial trial state.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports Nuclear Lattice Effective Field Theory (NLEFT) calculations for the proton-rich nucleus 22Si and its neighbors, using a global chiral NN+3N interaction and a second interaction fine-tuned to silicon isotopes. The authors compute ground-state energies, 2+ excitation energies, charge radii, and nucleon spatial distributions. They introduce a harmonic-oscillator-basis variant of the pinhole algorithm to extract shell-model occupation numbers. On this basis they claim that 22Si is more bound than 20Mg, excluding two-proton emission (S2p = 0.72(48) MeV with the global interaction and 1.34(45) MeV with the silicon-tuned interaction), and that 22Si is a doubly magic nucleus with Z=14 and N=8 shell closures.","tokens_in":15010,"tokens_out":7312,"duration_ms":64918,"significance":"If the claims are quantitatively sound, this would be an interesting first extension of NLEFT to a T_z=-3 proton-rich dripline nucleus, with a potentially useful new tool (HO-basis pinhole occupations) for connecting lattice simulations to the shell model. The paper also provides charge radii and a mirror-radius difference that could be tested by future experiments. However, both central claims require stronger support than the manuscript currently provides: the two-proton separation energy is only about 1.5 standard deviations above zero and appears to contradict an indirect experimental mass quoted in the paper, and the HO-basis occupation numbers are presented without a convergence check, center-of-mass correction, or benchmark on a known closed-shell nucleus.","major_comments":[{"comment":"The first central claim, that two-proton emission from 22Si is excluded, is not supported by the quoted numbers. The global interaction gives S2p = 0.72(48) MeV, which is only about 1.5 standard deviations above zero, so the abstract's wording 'excluding the possibility of two-proton emission' is too strong. Moreover, the same section quotes an indirect experimental g.s. energy of -134.51 MeV for 22Si and -134.61 MeV for 20Mg; if these are binding energies, they imply B(22Si) - B(20Mg) = -0.10 MeV, i.e., a negative S2p with the opposite sign from the calculation. The authors should discuss this discrepancy explicitly and replace the exclusion claim with a quantified confidence statement.","section":"Results, binding energies and S2p (Fig. 1)"},{"comment":"The HO-basis pinhole occupation numbers are not demonstrated to be ground-state observables. They are evaluated at a single Euclidean time tau = 0.2 MeV^-1, whereas the energies in Fig. 1 are extrapolated in tau and are still about 4 MeV above the extrapolated ground state at that tau. The supplemental text states that translations were not applied to the wavefunctions and gives only a 1.3% center-of-mass estimate for the ground-state energy, with no analogous check for occupation numbers. Table S1 contains negative occupation numbers, e.g., 0d3/2[3/2] = -0.127(169), and no benchmark on a known closed-shell nucleus such as 16O or 40Ca is provided. Because the trial wavefunction is a product of HO single-particle states with 0d5/2 filled, a short-time projection can partially reproduce the trial occupations. The shell-closure conclusion therefore needs a tau-convergence study and an independent validation before it can be regarded as established.","section":"Supplemental Material, 'The pinhole method in harmonic oscillator basis'; Table S1"},{"comment":"The claim that 22Si is doubly magic is also supported by qualitative indicators: the 2+ energies in N=8 isotones and the spatial compactness of the outermost proton distribution. These indicators are statistically weak; for example, the 2+ energies of 22Si and 20Mg are 2.11(57) MeV and 1.47(32) MeV, whose error bars overlap, and the spatial-distribution argument is not quantified. Even if the binding-energy calculation is correct, the doubly-magic conclusion should be based on a validated occupation-number analysis or another quantitative measure, rather than on these qualitative trends alone.","section":"Results, shell closure evidence from 2+ states and nucleon distributions"}],"minor_comments":[{"comment":"The word 'pertubatively' should be 'perturbatively'.","section":"Methods, paragraph 1"},{"comment":"The caption's explanation that the squares, circles, and stars are plotted as references on the E-axis, with only horizontal errors calibrated by the E-axis, is confusing and should be rewritten to clarify the axis meaning and the uncertainty representation.","section":"Fig. 2 caption"},{"comment":"The procedure for combining the two interaction results into a 3-sigma radius range is not specified; please state how the central value and the confidence interval are obtained from the two calculations.","section":"Results, charge radii"},{"comment":"The main text refers to Ref. [31] for several pieces of information; since [31] is the Supplemental Material, it would be clearer to name the supplement explicitly in the text.","section":"References and supplemental material"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is in scope for the journal and the calculations are potentially valuable. The main risks are overclaiming the S2p result and the doubly-magic conclusion before the HO-pinhole method is validated. I would encourage the editor to seek a revision that addresses the statistical framing and provides the missing convergence/benchmark tests."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead the 22Si NLEFT paper. The short version: it's a credible lattice calculation with a new idea, but the abstract oversells the main result. The two-proton separation energy comes out 0.72(48) MeV with the global interaction and 1.34(45) MeV with the silicon-tuned one. That's consistent and suggestive, but the global value is only 1.5 sigma above zero. Calling that 'excluding the possibility of two-proton emission' is too strong, and it sits awkwardly against the indirect experimental mass that would make S2p negative. The authors mention that tension but don't resolve it. This needs to be softened or argued more carefully.\n\nWhat's genuinely new: the first NLEFT calculation in the proton-rich N=8, T_z=-3 region, a new pinhole variant in the harmonic-oscillator basis, and a set of predictions—charge radius, mirror radius difference, 2+ energies—that are useful for future experiments. The two interactions agree on the qualitative picture, which gives me some confidence in the binding-energy part. The 2+ states track experiment for 18Ne and 20Mg, so the extrapolation scheme looks okay.\n\nThe soft spot is the doubly-magic claim. The occupation numbers in the HO basis are computed at a single Euclidean time tau = 0.2 MeV^-1, while energies are extrapolated from larger tau. At that tau the ground state is not yet fully projected, and the trial state is a product of HO orbitals. Without a convergence check in tau or a benchmark on a known closed-shell nucleus, the occupation numbers could simply be reflecting the trial state. Negative occupations in Table S1 (e.g., -0.127(169) for 0d3/2) underline that this is not a fully equilibrated measurement. The Supplement concedes no translations were applied to the wavefunctions and only estimates the CoM effect on the energy, not on the occupations. So the shell-closure conclusion rests on shakier ground than the binding-energy result.\n\nStill, the paper is serious work. The binding-energy prediction is a real testable statement, and the new method is worth exploring. I'd send it to review, but with a clear request for major revision: soften the abstract, show tau-convergence for the occupations or drop the strong claim, and benchmark the HO pinhole method on a known closed-shell nucleus.\n\nFor the reading group, it's a good discussion piece about how far lattice methods can push toward the dripline. I wouldn't cite it yet for the doubly-magic claim, though the radius predictions might be worth a footnote.","headline":"Credible NLEFT prediction for 22Si binding, but the abstract overstates the exclusion of 2p emission and the doubly-magic claim rests on an unvalidated pinhole method.","tokens_in":15541,"tokens_out":3690,"would_cite":false,"duration_ms":31363,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Lattice simulation finds that the exotic proton-rich nucleus 22Si is bound against two-proton emission and doubly magic with Z=14 and N=8 shell closures.","keywords":["nuclear lattice effective field theory","pinhole method","harmonic oscillator basis","22Si","proton dripline","shell closure","two-proton separation energy","occupation numbers"],"falsifier":"A direct mass measurement of $^{22}$Si that gives a two-proton separation energy $S_{2p} < 0$ would refute the central claim; short of that, applying the new harmonic-oscillator pinhole method to a known closed-shell nucleus such as $^{16}$O and obtaining occupation numbers that disagree with the expected closed-shell occupancies would refute the doubly-magic interpretation.","tokens_in":14441,"feed_emoji":"⚛️","tokens_out":15744,"duration_ms":115768,"temperature":0.7,"pith_summary":"Nuclear Lattice Effective Field Theory with high-fidelity chiral interactions is used to simulate the proton-rich nucleus $^{22}$Si, and the paper argues that $^{22}$Si is more tightly bound than $^{20}$Mg, with two-proton separation energy $S_{2p} = 0.72(48)$ MeV for the global interaction and $S_{2p} = 1.34(45)$ MeV for a silicon-tuned interaction, ruling out two-proton emission. The paper also introduces a pinhole method in a harmonic-oscillator basis that extracts orbital occupation numbers from the Monte Carlo state, and combines these numbers with coordinate-space nucleon distributions and the $2^+$ excitation energies of neighboring $N=8$ isotones to conclude that $^{22}$Si has shell closures at $Z=14$ and $N=8$ and is a doubly magic nucleus at the proton dripline. A sympathetic reader would care because $^{22}$Si sits at the edge of the nuclear landscape with extreme isospin asymmetry, where earlier ab initio calculations disagree about two-proton stability, so a lattice result that ties the dripline position to shell structure tests the chiral forces and many-body correlations in a new regime while offering a shell-model lens on the simulated wavefunction.","feed_headline":"Lattice simulation says 22Si is bound and doubly magic","feed_subtitle":"Two chiral interactions put 22Si above the two-proton decay threshold; occupancies show Z=14 and N=8 magic numbers.","key_machinery":"The central object is the pinhole method in a harmonic-oscillator basis: a new variant of the standard pinhole algorithm in which a complete set of $A$-body harmonic-oscillator basis states is inserted at the middle of the Euclidean-time projection, the occupation configuration is sampled by the Metropolis algorithm along with the auxiliary fields, and the occupation number of any orbital is the phase-weighted average over the sampled configurations. This is what converts the Monte Carlo wavefunction into a shell-model occupation pattern. The supporting machinery is the wavefunction matching of nuclear lattice effective field theory, which maps the chiral Hamiltonian $H_\\chi$ to a simple sign-problem-free Hamiltonian $H_S$ and treats the difference in perturbation theory up to $N^3$LO, plus the coordinate-space pinhole method that yields nucleon distance distributions relative to the center of mass.","core_discovery":"The central claim is that $^{22}$Si, probably the lightest bound nucleus with isospin projection $T_z = -3$, is bound against two-proton decay and is doubly magic. Using the global chiral interaction, the extrapolated ground-state energy of $^{22}$Si is $-134.38(39)$ MeV versus $-133.66(28)$ MeV for $^{20}$Mg, giving $S_{2p} = 0.72(48)$ MeV; the silicon-tuned interaction gives $S_{2p} = 1.34(45)$ MeV, and in both cases $^{22}$Si stays more bound than $^{20}$Mg by enough to exclude two-proton emission. The shell-closure claim rests on three pieces: the computed $2^+$ excitation energy of $^{22}$Si ($2.11(57)$ MeV with the global interaction) lies above those of $^{18}$Ne and $^{20}$Mg, the simulated proton distance distributions show an outer proton in $^{22}$Si that is less extended than in $^{20}$Mg, and the harmonic-oscillator occupation numbers from the new pinhole variant show the neutron $N=8$ shell essentially full, the six valence protons preferentially in $0d_{5/2}$, and only a minor contribution from $1s_{1/2}$, which together support $Z=14$ and $N=8$ closures.","pith_inferences":["The HO-pinhole method is not yet benchmarked on a known closed-shell nucleus; a natural extension is to compute occupation numbers for $^{16}$O or $^{40}$Ca in the same framework and compare with the expected closed-shell occupancies before trusting the $^{22}$Si numbers.","The occupation numbers are evaluated at a single Euclidean time $\\tau = 0.2$ MeV$^{-1}$, while energies and radii are extrapolated in $\\tau$; testing the convergence of the occupancies with $\\tau$ would tell whether the extracted shell structure is asymptotic or still contaminated by excited-state admixtures.","Since $^{22}$O is the well-established doubly magic mirror partner, a direct comparison of the full proton and neutron orbital occupancies of $^{22}$Si and $^{22}$O could isolate the Coulomb and isospin-breaking effects on the shell structure, rather than the radius difference alone.","If the $Z=14$ closure is as strong as implied, the dripline may extend one nucleon further; running the same simulation for $^{23}$Si would give a sharper, testable boundary for the proton dripline and could be a more sensitive falsifier of the predicted shell stabilisation."],"forward_implications":["If $^{22}$Si is truly bound by roughly $S_{2p} \\approx 0.7$–$1.3$ MeV, searches for two-proton radioactivity in $^{22}$Si should find none, and the proton dripline for the $N=8$ isotones is fixed at $^{22}$Si; a direct mass measurement can then discriminate between the two interaction sets.","The predicted charge radius of $^{22}$Si ($3.145$–$3.334$ fm) and the mirror radius difference $\\Delta R_{\\rm ch}(^{22}\\mathrm{Si}{-}^{22}\\mathrm{O}) \\approx 0.36$–$0.37$ fm, if measured, would place a constraint on the symmetry-energy slope $L$ of the nuclear equation of state.","The new HO-basis pinhole method, if reliable, gives lattice simulations direct access to shell-model occupation numbers, so other shell-model observables such as single-particle strengths or transition matrix elements could be extracted from the same wavefunctions.","The similarity between the computed $\\Delta R_{\\rm ch}(^{22}\\mathrm{Si}{-}^{20}\\mathrm{Mg})$ and the measured silicon–magnesium radius differences suggests that the $Z=14$ closure keeps the $^{22}$Si radius normal despite the small separation energy.","Consistency of the two interaction sets on $S_{2p}$ and on the $2^+$ energies indicates the qualitative conclusions are not sensitive to the 3N-force fine-tuning, so the remaining uncertainty is dominated by the simulation parameters and extrapolations rather than the interaction choice."],"supporting_citations":[{"why":"Supplies the high-fidelity chiral nucleon-nucleon and three-nucleon Hamiltonian plus the wavefunction-matching method that maps it to the sign-problem-free simple Hamiltonian used throughout.","marker":"[24]"},{"why":"Provides the alternative three-nucleon interaction fine-tuned to silicon isotope binding energies, used for the sensitivity analysis of the two-proton separation energy and radii.","marker":"[19]"},{"why":"Is the original pinhole algorithm for coordinate-space nucleon distributions, which the paper extends with the harmonic-oscillator basis variant.","marker":"[25]"},{"why":"Establishes the perturbation theory used to evaluate the chiral Hamiltonian on top of the simple Hamiltonian.","marker":"[26]"},{"why":"Gives the indirect experimental ground-state mass of $^{22}$Si ($-134.51$ MeV) that the calculation is compared with.","marker":"[13]"},{"why":"Documents the measured beta-delayed charged-particle decays of $^{22}$Si and the status that two-proton emission had not been ruled out, setting the experimental context for the $S_{2p}$ claim.","marker":"[14]"},{"why":"Are the ab initio calculations that give varying two-proton separation energies for $^{22}$Si, the disagreement this paper aims to resolve.","marker":"[15–17]"}],"fun_headline_variants":["22Si is bound and doubly magic, lattice simulation shows","Lattice: 22Si stable to two-proton emission, doubly magic","Pinhole method unveils 22Si's magic shells","22Si: new lattice data confirms double magic"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The doubly-magic conclusion depends on the unverified premise that the truncated harmonic-oscillator pinhole sampling at a single Euclidean time, with no center-of-mass projection, yields occupation numbers that faithfully represent the shell structure rather than artifacts of the basis and of center-of-mass motion.","fun_headline_variants_meta":{"raw":{"variants":["22Si is bound and doubly magic, lattice simulation shows","Lattice: 22Si stable to two-proton emission, doubly magic","Pinhole method unveils 22Si's magic shells","22Si: new lattice data confirms double magic"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000309,"raw_usage":{"total_tokens":1837,"prompt_tokens":1091,"completion_tokens":746,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":707,"completion_tokens_details":{"reasoning_tokens":676}},"tokens_in":707,"tokens_out":746,"duration_ms":7293,"temperature":1.0,"reasoning_tokens":676,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:06:24.987525+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct mass measurement of $^{22}$Si that gives a two-proton separation energy $S_{2p} < 0$ would refute the central claim; short of that, applying the new harmonic-oscillator pinhole method to a known closed-shell nucleus such as $^{16}$O and obtaining occupation numbers that disagree with the expected closed-shell occupancies would refute the doubly-magic interpretation.","supporting_citations":[{"cited_title":"A proton density bubble in the doubly magic $^{34}$Si nucleus","cited_arxiv_id":"1707.03583","evidence_quote":"Supplies the high-fidelity chiral nucleon-nucleon and three-nucleon Hamiltonian plus the wavefunction-matching method that maps it to the sign-problem-free simple Hamiltonian used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the indirect experimental ground-state mass of $^{22}$Si ($-134.51$ MeV) that the calculation is compared with."},{"cited_title":"Babo, β-delayed charged particle decays of neutron- deficient nuclei 20Mg and 22.23Si, Theses, Universit ´e de Caen Normandie (2016)","cited_arxiv_id":null,"evidence_quote":"Documents the measured beta-delayed charged-particle decays of $^{22}$Si and the status that two-proton emission had not been ruled out, setting the experimental context for the $S_{2p}$ claim."}],"review_version":1}