{"id":"c0415d76-4f25-4e80-986b-3cc99c4816b8","arxiv_id":"2411.10846","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A stochastic mean-field plus statistical decay calculation reproduces measured transuranium yields in 238U+248Cm collisions and predicts that elements with Z=102-105 are produced below the microbarn level.","lead":"This paper calculates how often a beam of uranium-238 colliding with a curium-248 target produces heavy, neutron-rich elements like californium, einsteinium, fermium, and mendelevium. The authors find their model matches measured yields to within an order of magnitude and predict that elements 102-105 would be produced at rates below one microbarn.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted sub-microbarn ceiling for Z=102-105 rests on curvature parameters α and β extracted by hand-picked time windows, with α borrowed from 250U+236Cm; no sensitivity analysis is provided, so the central predictive claim is not yet robust.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing point: the curvature parameters α and β, which control the quantal diffusion and hence the width and correlation of the primary fragment distribution, are extracted through manually selected time windows and α is imported from a neighboring reaction. Because the Z=102-105 predictions are extrapolations into the tail of this distribution, they are directly sensitive to these choices. The paper provides no uncertainty quantification and no sensitivity scan, so the central predictive claim is conditional on the unverified representativeness of the surrogate system and the chosen windows. An additional, more minor issue is that Eq. (12) contains an apparent extraneous factor of (2ℓ+1) that would make Psec not a probability and would double-count angular momentum when combined with Eq. (13); this should be corrected in a revision, but the primary concern remains the missing sensitivity analysis. The reader's CONDITIONAL verdict is therefore appropriate: the qualitative picture is plausible, but the quantitative statements, especially the sub-microbarn bound for Z=102-105, need to be demonstrated robust under variations of the extracted parameters.","tokens_in":17191,"tokens_out":10330,"duration_ms":111113,"concrete_test":"Recompute the SMF+GEMINI++ secondary cross-sections for Z=102-105 using the extreme curvature parameters from Table I: (α, β) = (0.108, 0.010) and (0.177, 0.004) for all four geometries, and also with the time windows tA, tB shifted by ±50 fm/c. If the integrated cross-sections for Z=102-105 remain below 1 μb in all tested cases, the prediction is robust; if any case exceeds 1 μb or changes by more than an order of magnitude relative to the published values, the central predictive claim requires revision and the stated sub-microbarn ceiling is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The production cross-sections are computed from the correlated Gaussian in Eqs. (2)-(3), whose variances and covariance evolve via Eqs. (4)-(6). These equations require the derivatives of the drift coefficients, which are not obtained directly from TDHF but from the parabolic potential energy surface (Eq. 10) through the reduced curvature parameters α and β (Eqs. 14-15). β is fitted to manually selected time intervals tA and tB on the 238U+248Cm drift path (Figs. 3-4), while α is taken from the neighboring system 250U+236Cm because the actual system has equal initial charge asymmetries (Figs. 5-6). Table I shows that the extracted parameters vary substantially across the four orientations: α ranges from 0.108 to 0.177, and β from 0.004 to 0.010. The yields for Z=102-105 are tail values of this correlated Gaussian, so even modest changes in the variance or correlation can shift these tail cross-sections by orders of magnitude. The paper claims no adjustable parameters, but the manual time-window selection is effectively a tunable input, and the transfer of α from a surrogate system is an unvalidated assumption. Without a sensitivity scan, the statement that Z=102-105 cross-sections remain below the microbarn level is not supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies the stochastic mean-field (SMF) quantal diffusion approach, combined with TDHF mean-field evolution and GEMINI++ statistical de-excitation, to the 238U + 248Cm reaction at Ec.m. = 898.7 MeV. It computes primary and secondary isotope production cross-sections for Z >= 98 (Cf through Db), compares the secondary yields for Z = 98-101 with the measured data of Kratz et al., and predicts that neutron-rich isotopes up to Z = 101 are produced with sizable cross-sections while Z = 102-105 remain below the microbarn level. The central claim is that the SMF approach 'explains' the available experimental results and is a predictive, essentially parameter-free microscopic tool for multi-nucleon transfer.","tokens_in":17415,"tokens_out":6666,"duration_ms":64546,"significance":"If the central claim held, this would be a valuable validation of the SMF quantal diffusion method for actinide multi-nucleon transfer, with concrete predictions for unexplored transuranium isotopes. The calculation is not fitted to the benchmark data, so the comparison is a genuine out-of-sample test; the inclusion of neutron-proton correlations and the detailed four-orientation TDHF tables are strengths. However, the agreement with the measured cross-sections is only order-of-magnitude (factor 3-10, with shifted peak mass numbers), and the predicted sub-microbarn ceiling for Z = 102-105 rests on curvature parameters extracted from manually chosen time windows, with no sensitivity analysis. These issues do not invalidate the method but do undermine the overstrong 'explains' claim and the robustness of the prediction.","major_comments":[{"comment":"The comparison in Table IV shows only order-of-magnitude agreement, not the 'explain' stated in the Abstract and Sec. IV: for Z=98 the secondary peak is at A=250 with 0.60 mb vs measured A=251 with 1.87 mb; for Z=100 the secondary peak is 0.02 mb vs 0.002 mb, a factor of 10. Please rephrase the claim to 'reproduced within a factor of roughly 3-10' and discuss possible sources of the peak shift (e.g., GEMINI++ decay treatment, equal-weight orientation averaging).","section":"Sec. III, Table IV"},{"comment":"The reduced curvature parameters alpha and beta, which control the drift and hence the correlated Gaussian variances in Eqs. (4)-(6), are extracted from manually selected time windows tA and tB, and alpha is borrowed from the surrogate system 250U + 236Cm. The Z=102-105 predictions are tail values of this correlated Gaussian and are highly sensitive to these parameters, yet no sensitivity analysis is provided with respect to the window choice, the surrogate alpha, or the effective temperature T*. Without such an analysis, the sub-microbarn ceiling conclusion is not robust.","section":"Sec. II.C, Eqs. (14)-(15), Table I"},{"comment":"As printed, Eq. (12) includes an explicit (2*ell+1) factor in the definition of Psec_ell, and Eq. (13) multiplies Psec_ell by (2*ell+1) again, double-counting the angular momentum weight. This makes the normalization of the secondary cross-section inconsistent; please clarify the intended definitions, since the absolute cross-section values in Table IV depend on this factor.","section":"Sec. II.A, Eqs. (12)-(13)"},{"comment":"The claim that the SMF theory 'does not contain any adjustable parameters' is not supported: the extraction of alpha and beta involves manual selection of the time intervals tA and tB (Table I), and the Einstein relations in Eqs. (11) introduce an effective temperature T* whose value is not given. At minimum, the time-window choices should be identified as tunable inputs, or a sensitivity study should demonstrate that the results are insensitive to them.","section":"Sec. IV and Abstract"}],"minor_comments":[{"comment":"The notation alternates between 'C_l' and 'C_ell'; please use a single symbol, and ensure the subscript is typeset consistently.","section":"Eqs. (2)-(3)"},{"comment":"The ell=360 rows for XX and XY list A2f values (257.1 and 284.7) that are inconsistent with mass conservation; these are likely typos.","section":"Table II"},{"comment":"The caption contains a garbled string '/uni00000014/...' that appears to be a font or encoding artifact; please replace it with the proper caption text.","section":"Fig. 10 caption"},{"comment":"The phrase 'closed to the neutron and proton drip lines' should read 'close to the neutron and proton drip lines'.","section":"Sec. I"},{"comment":"The statement that 'The SMF results are consistent with the experimental values in terms of magnitude and peak points' is contradicted by Table IV, where the secondary peak mass numbers differ from the measured ones for Z=98-101.","section":"Sec. III"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for a nuclear physics journal. The gap between the claimed 'explains' and the actual factor-of-3-10 agreement, together with the unquantified sensitivity of the Z=102-105 predictions to the manually selected curvature parameters, is the main reason for major revision. I would be willing to see a revised version with a sensitivity analysis and tempered claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plain take: this is a solid application of the authors' own SMF quantal diffusion machinery to a new actinide+actinide system, 238U+248Cm, with concrete cross-section maps for Z=98-101 and a first quantitative statement that Z=102-105 fall below a microbarn. That is genuinely useful for experiment planning. But the paper overclaims when it says the calculations 'explain' the data: Table IV shows order-of-magnitude agreement, with factors of 3-10 between calculated and measured secondary peaks. The central prediction for Z=102-105 is conditional, not robust, because the curvature parameters α and β that set the drift and variances are extracted by hand-picked time windows, and α is borrowed from 250U+236Cm. No sensitivity scan is given.\n\nWhat is new: first SMF calculation for this system; predictions for No, Lr, Rf, and Db are new. What is done well: the four collision orientations are treated separately; TDHF drift paths and diffusion coefficients are well documented; σNZ correlations are included; GEMINI++ de-excitation is standard; the qualitative trend—cross sections fall steeply with Z and shift to higher masses—is reproduced. The comparison to the measured data is at least in the right ballpark.\n\nSoft spots, in order:\n1. The 'no adjustable parameters' claim is too strong. The reduced curvature parameters are fit to the model's own drift paths over manually chosen windows (tA, tB), and α is transferred from a neighboring system because the actual system has equal initial charge asymmetries. Table I shows these parameters vary substantially with orientation (α 0.108-0.177, β 0.004-0.010). The yields for Z=102-105 are tail values of a correlated Gaussian; even modest changes can shift them by orders of magnitude. The paper needs a sensitivity scan over α, β, and time windows before the sub-microbarn ceiling can be taken as quantitative.\n2. 'Explain' is an overstatement. In Table IV, the secondary peak for Z=98 is 0.60 mb at A=250 versus measured 1.87 mb at A=251; for Z=101, 0.001 mb versus 0.0002 mb. That is order-of-magnitude reproduction, which is worth reporting, but not explanation.\n3. There are numerical typos and corrupted text: Table II XX ℓ=360 A_f2=257.1 looks off, Table III XY ℓ=360 A_f2=284.7 is impossible next to neighboring ~248, Eq. (15) has a sign issue with tA-tB, and the Figure 10 caption contains stray '/uni...' code. These are addressable but need fixing.\n\nWho it is for: MNT experimentalists planning transuranium production, and theorists benchmarking SMF. It deserves a serious referee, but should come back with a rebenchmarked abstract, a clear fitted-versus-predicted statement, and sensitivity tests on the tail predictions.","headline":"A useful but overclaimed application of the authors' SMF quantal diffusion framework to 238U+248Cm; the qualitative picture and new Z=102-105 predictions are worth having, but the 'explains the data' statement and the sub-microbarn ceiling need to be hedged and rebenchmarked.","tokens_in":18079,"tokens_out":3071,"would_cite":true,"duration_ms":34124,"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":"The paper claims that the stochastic mean-field approach combined with GEMINI++ explains the measured multi-nucleon transfer cross-sections in $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ at $E_\\mathrm{c.m.}=898.7$ MeV and predicts sizable…","keywords":["multi-nucleon transfer reactions","stochastic mean-field theory","quantal diffusion","time-dependent Hartree-Fock","transuranium isotopes","GEMINI++","actinide collisions","neutron-rich nuclei"],"falsifier":"Measure the secondary cross-sections for $Z=102$--$105$ isotopes in $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ at $E_\\mathrm{c.m.}=898.7$ MeV; finding any of them above 1 microbarn, or finding the $Z=98$--$101$ yields more than roughly an order of magnitude away from the predicted millibarn-to-microbarn values, would show the calculation's central prediction is wrong.","tokens_in":16844,"feed_emoji":"⚛️","tokens_out":8402,"duration_ms":75347,"temperature":0.7,"pith_summary":"This paper argues that the stochastic mean-field (SMF) approach, a quantal extension of time-dependent Hartree-Fock that includes neutron and proton fluctuations, can describe multi-nucleon transfer in the heaviest actinide collisions. Using $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ at $E_\\mathrm{c.m.}=898.7$ MeV, the authors compute primary fragment cross-sections from quantal diffusion coefficients and then pass the hot fragments through the GEMINI++ statistical decay code to obtain secondary (surviving) yields. They find that the calculation reproduces the measured cross-sections for transuranium products from Cf ($Z=98$) through Md ($Z=101$), and that the same approach predicts these elements are producible with sizable cross-sections while $Z=102$--$105$ (No, Lr, Rf, Db) remain below one microbarn. The significance, if correct, is a parameter-free microscopic route to choosing reactions that can extend the nuclear chart to new neutron-rich heavy isotopes.","feed_headline":"At 898.7 MeV, 238U+248Cm makes neutron-rich isotopes up to Z=101","feed_subtitle":"A stochastic mean-field calculation matches measured transfer yields and sets Z=102-105 below the microbarn level.","key_machinery":"The machinery is the quantal diffusion description of multi-nucleon transfer. Diffusion coefficients for neutrons and protons are computed directly from the occupied single-particle wave functions of the TDHF evolution, including quantal shell effects and Pauli blocking; the drift coefficients are linked through Einstein relations to a parabolic potential energy surface in the $(N,Z)$ plane, with an isoscalar curvature parameter $\\beta$ and an isovector curvature parameter $\\alpha$. Since the isovector curvature for the $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ system cannot be obtained directly, it is estimated from the neighbouring $^{250}\\mathrm{U}+^{236}\\mathrm{Cm}$ reaction, and both parameters are averaged over selected time windows. These inputs feed a correlated Gaussian (Fokker-Planck) probability distribution for primary fragments, whose variances satisfy coupled differential equations, and the resulting hot fragments are decayed statistically with GEMINI++. The nonzero neutron-proton covariance in the distribution is what aligns the yield ellipses along the valley of stability, a feature that standard mean-field treatments miss.","core_discovery":"The central claim is that a stochastic mean-field description, with no adjustable parameters beyond the Skyrme energy functional, reproduces the measured isotope cross-sections for the $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ reaction and can therefore be used to predict unreachable regions. For the four initial geometries (tip-tip XX, tip-side XY, side-tip YX, side-side YY), the drift paths and quantal diffusion coefficients are generated from TDHF trajectories; a correlated Gaussian probability distribution for neutron and proton transfers is built from the resulting variances, including a nonzero neutron-proton covariance, and the de-excitation of primary fragments is followed with GEMINI++. The paper finds primary peaks at $A=251$ (Cf, 10.4 mb), $A=253$ (Es, 6.20 mb), $A=256$ (Fm, 4.40 mb), and $A=259$ (Md, 3.20 mb), with secondary peaks shifted down by neutron evaporation, and reports agreement with experimental data in magnitude and location. For $Z=102$--$105$, where no data exist, the calculation predicts cross-sections below 1 microbarn, so those elements are not expected to be produced at observable rates in this reaction.","pith_inferences":["One could extend the same machinery to other actinide combinations, such as $^{238}\\mathrm{U}+^{238}\\mathrm{U}$ or $^{248}\\mathrm{Cm}+^{250}\\mathrm{Cf}$, to see whether the method identifies a system with higher predicted yields for $Z=102$--$105$.","The strong sensitivity to the borrowed $\\alpha$ parameter suggests that a systematic uncertainty band, obtained by varying the time window and the surrogate system, would show whether the 'below microbarn' conclusion is robust.","Because the diffusion coefficients come from the TDHF mean field, the predicted yields can also serve as a probe of the isovector part of the Skyrme energy functional used in the calculation."],"forward_implications":["If the calculation is correct, the same SMF+GEMINI++ procedure can rank candidate projectile-target combinations for producing new neutron-rich heavy isotopes without tuning reaction parameters.","The computed cross-section ladder (millibarn for $Z=93$--$98$, microbarn for $Z=99$--$102$, nanobarn for $Z=103$--$110$) sets a quantitative expectation of where multi-nucleon transfer becomes impractical in these actinide systems.","The secondary peaks being shifted several nucleons below the primary peaks means the observable nuclei are not the primary transfer products; neutron evaporation and fission compete strongly, so survival probabilities must be included in any yield estimate.","The prediction that $Z=102$--$105$ stays below 1 microbarn at this energy implies that simply repeating the same reaction with higher statistics would not open the region; different systems or energies would be needed."],"supporting_citations":[{"why":"Supplies the quantal diffusion description of multinucleon transfers from which the primary cross-section formula and diffusion coefficients are taken.","marker":"[34]"},{"why":"Provides the experimental isotope cross-sections for $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ that the SMF+GEMINI++ results are compared against.","marker":"[43]"},{"why":"Gives the earlier $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ actinide production data that motivates the system choice and anchors the comparison.","marker":"[44]"},{"why":"The GEMINI++ statistical de-excitation code that converts the primary fragment yields into secondary cross-sections.","marker":"[56]"},{"why":"Introduces the inverse quasifission mechanism and the XX/XY/YX/YY geometry notation used to define the four collision orientations.","marker":"[61]"},{"why":"The three-dimensional TDHF code used to generate mean-field trajectories, fragment drifts, and total kinetic energies.","marker":"[62,63]"},{"why":"The SLy4d Skyrme energy density functional that fixes the mean-field Hamiltonian and hence all transport coefficients.","marker":"[64]"}],"fun_headline_variants":["Parameter-free SMF predicts neutron-rich isotopes up to Z=101 in U+Cm","Stochastic model reproduces U+Cm transfer yields, predicts Z=101","Neutron-rich isotopes up to Z=101 from U+Cm at 898.7 MeV","SMF with no adjustable params nails U+Cm yields, predicts Z=101","Microscopic model predicts new isotopes to Z=101 in U+Cm reaction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction rests on taking the isovector curvature parameter $\\alpha$ from a neighbouring reaction and on the time windows chosen for averaging, so a mismatch in that surrogate potential would shift the computed yields.","fun_headline_variants_meta":{"raw":{"variants":["Parameter-free SMF predicts neutron-rich isotopes up to Z=101 in U+Cm","Stochastic model reproduces U+Cm transfer yields, predicts Z=101","Neutron-rich isotopes up to Z=101 from U+Cm at 898.7 MeV","SMF with no adjustable params nails U+Cm yields, predicts Z=101","Microscopic model predicts new isotopes to Z=101 in U+Cm reaction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001611,"raw_usage":{"total_tokens":6521,"prompt_tokens":1157,"completion_tokens":5364,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":773,"completion_tokens_details":{"reasoning_tokens":5254}},"tokens_in":773,"tokens_out":5364,"duration_ms":38872,"temperature":1.0,"reasoning_tokens":5254,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:14:30.998780+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the secondary cross-sections for $Z=102$--$105$ isotopes in $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ at $E_\\mathrm{c.m.}=898.7$ MeV; finding any of them above 1 microbarn, or finding the $Z=98$--$101$ yields more than roughly an order of magnitude away from the predicted millibarn-to-microbarn values, would show the calculation's central prediction is wrong.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the quantal diffusion description of multinucleon transfers from which the primary cross-section formula and diffusion coefficients are taken."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental isotope cross-sections for $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ that the SMF+GEMINI++ results are compared against."},{"cited_title":"Schädel, W","cited_arxiv_id":null,"evidence_quote":"Gives the earlier $^{238}\\mathrm{U}+^{248}\\mathrm{Cm}$ actinide production data that motivates the system choice and anchors the comparison."},{"cited_title":"Charity, GEMINI: A code to simulate the decay of a com- pound nucleus by a series of binary decays , Tech","cited_arxiv_id":null,"evidence_quote":"The GEMINI++ statistical de-excitation code that converts the primary fragment yields into secondary cross-sections."},{"cited_title":"Kedziora and Cédric Simenel, New inverse quasifis- sion mechanism to produce neutron-rich transfermium nuclei, Phys","cited_arxiv_id":null,"evidence_quote":"Introduces the inverse quasifission mechanism and the XX/XY/YX/YY geometry notation used to define the four collision orientations."}],"review_version":1}