{"id":"8d3677a9-f2ef-4a2f-88eb-3ec7ebd9356c","arxiv_id":"2505.04975","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Two low-lying 0+ states in 24Mg are identified as credible candidates for a 16O-plus-two-alpha condensate, with small predicted alpha decay widths.","lead":"This theoretical study uses a real-time wave-packet method to look for two alpha particles forming a condensate-like cloud around an oxygen-16 core inside magnesium-24. It identifies two states near the decay threshold as credible candidates and predicts that they may be observable in experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The omitted proton-decay channel undermines the claimed small width and 'firmly establishing' language for the 0+4 state, although the condensate interpretation itself may survive.","rationale":"The paper is a credible methodological advance: REM reduces continuum contamination relative to random basis generation, and ACCC supplies explicit width estimates. The agreement in B(IS0) with Ichikawa et al. is expected because the Hamiltonian and model space are nearly identical; it is consistency, not independent confirmation. I do not see an internal inconsistency in the REM derivation itself. The most load-bearing weakness is the open proton channel for 0+4, which directly affects the width and observability part of the central claim. The reader's frozen-core concern is real, but it is less decisive here: core excitation would modify the model space and could shift states, whereas the proton channel is an explicitly acknowledged open decay channel whose omission undermines the quantitative width prediction. This concern does not warrant rejection; it strengthens the case for a conditional verdict with a coupled-channel calculation as the path to firm establishment.","tokens_in":11118,"tokens_out":9690,"duration_ms":108422,"concrete_test":"Extend the REM+ACCC calculation by coupling the 0+ states to a 23Na+p channel, e.g., by adding a proton-plus-core generator-coordinate basis or by discretizing the proton continuum, and recompute the complex pole energies E_R - i Gamma/2 for 0+3 and 0+4 at the physical coupling constant. If Gamma_p(0+4) is comparable to or larger than the reported Gamma_alpha = 0.17 MeV, or if E_R shifts by more than a few hundred keV, the paper's claim that 0+4 is an observable narrow alpha-decay resonance is not supported. The same calculation would verify the stated 'minor' proton effect for 0+3 at Q_p = 0.14 MeV.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the calculated alpha widths dominate the physical decay widths. Section III.B concedes that the proton-decay channel is outside the model space and that, for the 0+4 state, its influence may be non-negligible. Using the manuscript's own numbers, the 0+4 state lies at -1.26 MeV while the proton threshold is at -5.09 MeV, i.e., about 3.8 MeV above the proton threshold; the 0+3 state has Q_p = 0.14 MeV and is less affected. An open channel with this much Q-value can add a substantial proton width, so the reported Gamma_alpha = 0.17 MeV for 0+4 does not bound the total width. The conclusions 'experimentally observable' and 'firmly establishing' therefore outrun a model space that excludes an open decay channel. The frozen tetrahedral core is a related restriction, but the proton channel is the more immediate threat because the paper itself flags it for 0+4.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies 24Mg in a 16O + 2α cluster model with a frozen tetrahedral 16O core, using the real-time evolution method (REM) to generate basis states from the equations of motion and then solving the Hill-Wheeler equation. The authors examine convergence with evolution time T, reflecting-wall radius R, and initial intrinsic excitation E*, and they use analytical continuation in the coupling constant (ACCC) to extract α-decay widths. They identify the 0+3 and 0+4 states as candidates for 2α condensation around 16O, with strong isoscalar monopole transition strengths B(IS0)=145.0 and 160.7 fm^4, valence radii near 4.1 fm, and α widths ≤0.01 and 0.17 MeV. The paper concludes that these states are firmly established condensate candidates and should be experimentally observable.","tokens_in":11380,"tokens_out":10431,"duration_ms":93105,"significance":"If the central claim holds, this is a valuable step toward a robust identification of core + nα condensate candidates, providing quantitative predictions for transition strengths and widths that can be compared with experiment. The REM approach is a genuine methodological improvement over random basis generation, and the paper is transparent about the model parameters and includes convergence checks over T, R, and E*. The B(IS0) values are in good agreement with Ref. [31], which is reassuring. However, the strength of the conclusion currently exceeds what the model evidence supports, mainly because the proton decay channel is omitted and because the convergence and error estimates for the resonance properties are incomplete. The manuscript is clearly written and the numerical procedure is described in enough detail to be reproduced.","major_comments":[{"comment":"The energies of the unbound states 0+3 and 0+4 continue to drift downward with evolution time up to T = 4000 fm/c, and the text states that the energy 'slowly decreases over time due to coupling with the continuum.' A convergence criterion (for example, a plateau in E(T) over a substantial interval, or an extrapolation in 1/T) is needed before these eigenvalues can be presented as converged resonance energies; without it, the ACCC input depends on the arbitrary cutoff T.","section":"Sec. III.A, Fig. 1"},{"comment":"The reported Γ_α = 0.17 MeV for the 0+4 state does not bound the total decay width because the proton channel is outside the model space; the paper itself concedes that its influence on 0+4 'may be non-negligible.' With the proton threshold at −5.09 MeV and the 0+4 energy at −1.26 MeV, the proton Q-value is about 3.8 MeV, so the proton width may be comparable to or larger than the α width. The conclusions that the states are 'experimentally observable' and 'firmly established' therefore overstate what the model can support; the authors should either estimate the proton width or explicitly restrict the claim to the α-decay channel within the model space.","section":"Sec. III.B, Table I"},{"comment":"No uncertainties are propagated for the ACCC widths or for the adopted energies; the convergence with R and E* is described only qualitatively, and the 0+3 width is quoted as '≤0.01 MeV' without stating whether this is an upper limit from the Padé extrapolation, a numerical resolution, or a statistical error. The choice M = N = 6 is also asserted to be 'large enough' without showing the rank dependence. Because the small widths are load-bearing for the observability claim, the authors should quantify the spread from R, E*, and Padé-rank variations.","section":"Sec. II.D and Sec. III.B, Table I"},{"comment":"The 16O core is frozen in a tetrahedral configuration with a fixed side length of 0.5 fm, and no sensitivity study is provided for this choice. Since the core size and stiffness can affect the B(IS0) values and the barrier that confines the valence alphas, a variation of the side length (or allowing core breathing) would test whether the identification of 0+3 and 0+4 as condensate candidates is robust against this model assumption.","section":"Sec. II.A, Eqs. (5)-(6)"}],"minor_comments":[{"comment":"The denominator in Eq. (18) is written as 2µ, but the reduced mass was denoted by m in Eq. (16); since µ is also used as the ACCC coupling parameter, this is confusing and likely a typo for 2m.","section":"Sec. II.D, Eq. (18)"},{"comment":"The terms 'rebound radius' and 'reflecting wall' are used interchangeably; please choose one consistent term, preferably 'reflection radius.'","section":"Sec. II.C and figure captions"},{"comment":"The label 'α threhsold' should read 'α threshold'; please also check that the red dashed lines indicating Ref. [31] energies are clearly visible in the printed version.","section":"Figures 1-3"},{"comment":"For the 0+1 and 0+2 bound states, the 'width' and 'B(IS0)' columns contain dashes; please state explicitly that these quantities are not applicable for bound states, and note that B(IS0) is given for 0+2 only.","section":"Table I"},{"comment":"The text says sampled wave functions with overlap greater than 0.95 are removed, but the displayed condition requires the normalized squared overlap to be < 0.95 for every pair; please clarify the pruning algorithm, for example by stating that one member of each pair exceeding the threshold is discarded.","section":"Sec. II.C, step 3"},{"comment":"The phrase 'firmly establishing' is stronger than the evidence presented, given the model-space limitations and the ACCC uncertainties; 'supporting the identification of' or 'providing evidence for' would be more appropriate.","section":"Abstract and Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid model study with a useful methodological contribution, but the central claim is currently overstated. The proton-channel caveat is acknowledged in the text, yet the abstract and conclusions still say 'firmly establishing' and 'experimentally observable'; that mismatch should be resolved. The convergence drift of the unbound states with T and the lack of uncertainty estimates for the ACCC widths are additional load-bearing issues that need concrete responses. I would not recommend rejection, as the issues are fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real content here is the method application: real-time evolution generation of basis states for 16O+2α, plus ACCC widths and valence radii. That is new relative to Ichikawa et al., and the convergence checks against the reflection radius and E* are a genuine improvement over the older random-basis approach. The B(IS0) values for 0+3 and 0+4 (145 and 161 fm^4) agree closely with the prior prediction, and the valence radii near 4 fm give a sensible physical picture of α particles trapped outside the core. Credit where due: the paper is transparent about its model limitations and does not hide the proton-channel problem for 0+4.\n\nThe soft spot is not subtle. The authors claim the small α widths make these states experimentally observable, but the model space excludes proton decay, and for 0+4 the proton Q-value is about 3.8 MeV. An open channel with that much energy can add a substantial width, so Γα = 0.17 MeV does not bound the total width. The paper itself concedes this in Section III.B, yet the abstract and conclusion still say \"firmly establishing.\" That language is not supported, and a referee should push them to either estimate the proton width with a larger model or restrict the claim to 0+3, whose proton Q-value is small and whose width is only an upper limit anyway. The frozen tetrahedral core is a related but less immediate issue. Also, the unbound-state energies drift with evolution time, and the ACCC width for 0+3 is reported as an upper limit with no uncertainty estimate, so the quantitative widths should be read as indicative rather than firm.\n\nEven with these caveats, the condensate interpretation for 0+3 survives; the paper's core result, that REM identifies two resonant states with strong monopole strength and small α widths, is a reasonable confirmation of an earlier prediction using a better solver. The overreach is in the conclusions, not in the calculation.\n\nThis paper deserves a serious referee. The method is sound, the claims are falsifiable, and the proton-channel gap can be addressed. I would send it to review, but I would ask the authors to soften the \"firmly establishing\" language and to present the 0+4 width as a partial width conditional on the missing proton channel. For a cluster-physics audience, this is a useful contribution; for experimentalists planning 24Mg searches, the takeaway should be that only 0+3 is a clean candidate right now.","headline":"Solid REM re-examination of 16O+2α candidates that strengthens the case for the 0+3 state but overreaches for 0+4 by ignoring the open proton channel.","tokens_in":11860,"tokens_out":1412,"would_cite":true,"duration_ms":15611,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Two 24Mg states are the strongest core+α condensate candidates yet","keywords":["alpha condensation","core-plus-alpha cluster","real-time evolution method","analytical continuation in coupling constant","isoscalar monopole transition","24Mg","alpha decay width","cluster resonance"],"falsifier":"An experiment or an extended calculation that opens the proton channel would settle the claim: the proton threshold lies at $-5.09$ MeV in the model, only $0.14$ MeV above the $0^+_3$ state, so a proton width comparable to the predicted $\\alpha$ width would rule out the $0^+_3$ assignment; conversely, $\\alpha$ inelastic scattering on $^{24}$Mg that finds no narrow $0^+$ states with monopole strength near $145$-$161$ fm$^4$ would falsify the condensate identification.","tokens_in":10931,"feed_emoji":"⚛️","tokens_out":17012,"duration_ms":153365,"temperature":0.7,"pith_summary":"This paper sets out to settle whether two excited $0^+$ states of $^{24}$Mg are genuine resonances in which two $\\alpha$ particles form a gas-like cloud around an $^{16}$O core, a core+$2\\alpha$ condensate. In earlier random-basis calculations these states were obscured by strong mixing with non-resonant continuum configurations. Using the real-time evolution method to generate the basis and analytic continuation in the coupling constant to extract decay widths, the authors find that the $0^+_3$ state near $-4.95$ MeV and the $0^+_4$ state near $-1.26$ MeV carry very large isoscalar monopole transition strengths ($145.0$ and $160.7$ fm$^4$), valence $\\alpha$ radii near $4$ fm, and small $\\alpha$-decay widths ($\\le 0.01$ and $0.17$ MeV). The point of the paper is that these two states are the most firmly established candidates for $\\alpha$ condensation around a core nucleus and that they should be observable in experiment.","feed_headline":"Two 24Mg states are the strongest core+α condensate candidates yet","feed_subtitle":"New real-time evolution calculation gives them large monopole strength and tiny alpha widths, making them observable.","key_machinery":"The engine is the real-time evolution method (REM): instead of drawing Gaussian cluster centroids at random, it solves the time-dependent variational equations of motion for the $\\alpha$-particle positions and momenta, so the sampled wave packets follow physically relevant trajectories and the Hill-Wheeler basis built from them carries little continuum contamination. The analytic continuation in the coupling constant (ACCC) supplies the resonance parameters: an auxiliary potential is added to artificially bind the state, the bound energies are computed as a function of the coupling, and a Padé approximant in $\\sqrt{\\mu-\\mu_0}$ extrapolates the complex energy to the physical point, giving each resonance energy and width. The isoscalar monopole transition strength $\\mathrm{B(IS0)}$ from the ground state is the observable used to identify spatially extended cluster states.","core_discovery":"Within a microscopic $^{16}$O+$2\\alpha$ cluster model, the paper identifies the $0^+_3$ and $0^+_4$ states of $^{24}$Mg as resonant candidates for $2\\alpha$ condensation around the $^{16}$O core. The real-time evolution method yields a convergent spectrum with far less contamination from non-resonant states than the earlier random-basis calculation, and analytic continuation in the coupling constant gives small $\\alpha$-decay widths for these states. Their isoscalar monopole transition strengths, $\\mathrm{B(IS0)}=145.0$ and $160.7$ fm$^4$, are close to the earlier calculation's values, and their valence $2\\alpha$ distribution radii, about $4.1$ fm, match the expected radius of the Coulomb barrier around $^{16}$O, suggesting two $\\alpha$ particles trapped between the core surface and the barrier. The paper concludes that the $0^+_3$ and $0^+_4$ states, roughly $5$ and $1$ MeV below the $^{16}$O+$2\\alpha$ threshold, are strong candidates whose small $\\alpha$ widths make experimental observation feasible.","pith_inferences":["A natural next step is to include the proton channel, which the model omits; the proton threshold at $-5.09$ MeV lies only $0.14$ MeV above the $0^+_3$ state, so a microscopic calculation including protons could alter the $0^+_4$ width and test the assignment.","The same machinery could be turned on heavier cores and more valence alphas, where random-basis searches have been even more ambiguous; if REM keeps its convergence there, core+$n\\alpha$ condensates may turn out to be a general near-threshold phenomenon.","The predicted valence radius of about $4$ fm, matched to the Coulomb barrier, implies a characteristic oscillatory structure in the monopole transition form factor; extracting that form factor from inelastic $\\alpha$ scattering would test the spatial picture directly.","Because the core is frozen in a tetrahedron, a fully dynamical core could shift the two candidate states; the closeness of the $0^+_4$ state to the proton threshold makes the ordering sensitive to such effects."],"forward_implications":["If the $0^+_3$ and $0^+_4$ assignments are correct, $^{24}$Mg becomes the clearest known case of a core-plus-alpha condensate, and its narrow $0^+$ resonances can be searched for directly.","The small $\\alpha$-decay width of the $0^+_3$ state ($\\le 0.01$ MeV) makes it a particularly clean experimental target, with the $0^+_4$ state ($0.17$ MeV) still narrow enough to separate from the background.","The close agreement of the monopole strengths with the earlier calculation cross-validates both methods, so $\\mathrm{B(IS0)}\\approx 145$ and $161$ fm$^4$ become concrete benchmarks for future experiments.","The same REM-plus-ACCC procedure can be applied to other core-plus-alpha systems, such as $^{16}$O+$3\\alpha$ ($^{28}$Si) or $^{40}$Ca+$3\\alpha$."],"supporting_citations":[{"why":"the earlier random-basis calculation whose $^{16}$O+$2\\alpha$ condensate candidates this paper verifies and sharpens.","marker":"[31]"},{"why":"introduces the real-time evolution method used to generate the model basis with reduced continuum contamination.","marker":"[36]"},{"why":"supplies the analytic continuation in the coupling constant used to extract resonance energies and $\\alpha$-decay widths.","marker":"[37]"},{"why":"establishes the isoscalar monopole transition strength as the signature of spatially extended cluster states.","marker":"[35]"},{"why":"provides the Volkov No. 2 effective interaction used in the Hamiltonian.","marker":"[38]"},{"why":"provides the Brink-Bloch cluster wave function on which the $^{16}$O+$2\\alpha$ model is built.","marker":"[39]"}],"fun_headline_variants":["Real-time method reveals 2α condensate candidates in 24Mg","2α condensate around 16O: two 24Mg states stand out","Small α widths make 24Mg 0+ states prime 2α condensate candidates","Real-time evolution sharpens 2α condensate signal in 24Mg","Two 24Mg states are prime candidates for 2α condensation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The $^{16}$O core is treated as inert, frozen into four $\\alpha$ clusters at the corners of a $0.5$ fm tetrahedron with only its center of mass mobile, so every $^{24}$Mg state is described purely as $^{16}$O+$2\\alpha$ and core excitation, deformation, and proton decay are absent.","fun_headline_variants_meta":{"raw":{"variants":["Real-time method reveals 2α condensate candidates in 24Mg","2α condensate around 16O: two 24Mg states stand out","Small α widths make 24Mg 0+ states prime 2α condensate candidates","Real-time evolution sharpens 2α condensate signal in 24Mg","Two 24Mg states are prime candidates for 2α condensation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000971,"raw_usage":{"total_tokens":4202,"prompt_tokens":1092,"completion_tokens":3110,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":708,"completion_tokens_details":{"reasoning_tokens":3008}},"tokens_in":708,"tokens_out":3110,"duration_ms":22331,"temperature":1.0,"reasoning_tokens":3008,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:16:11.416436+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An experiment or an extended calculation that opens the proton channel would settle the claim: the proton threshold lies at $-5.09$ MeV in the model, only $0.14$ MeV above the $0^+_3$ state, so a proton width comparable to the predicted $\\alpha$ width would rule out the $0^+_3$ assignment; conversely, $\\alpha$ inelastic scattering on $^{24}$Mg that finds no narrow $0^+$ states with monopole strength near $145$-$161$ fm$^4$ would falsify the condensate identification.","supporting_citations":[{"cited_title":"Kokalova, N","cited_arxiv_id":null,"evidence_quote":"the earlier random-basis calculation whose $^{16}$O+$2\\alpha$ condensate candidates this paper verifies and sharpens."},{"cited_title":"Ichikawa, N","cited_arxiv_id":null,"evidence_quote":"introduces the real-time evolution method used to generate the model basis with reduced continuum contamination."},{"cited_title":"Kawabata, H","cited_arxiv_id":null,"evidence_quote":"supplies the analytic continuation in the coupling constant used to extract resonance energies and $\\alpha$-decay widths."},{"cited_title":"Ichikawa, N","cited_arxiv_id":null,"evidence_quote":"establishes the isoscalar monopole transition strength as the signature of spatially extended cluster states."},{"cited_title":"Kanada-En’yo, Negative parity states of 11b and 11c a nd the similarity with 12c, Physical Review C 75, 024302 (2007)","cited_arxiv_id":null,"evidence_quote":"provides the Volkov No. 2 effective interaction used in the Hamiltonian."},{"cited_title":"Yamada, Y","cited_arxiv_id":null,"evidence_quote":"provides the Brink-Bloch cluster wave function on which the $^{16}$O+$2\\alpha$ model is built."}],"review_version":1}