{"id":"ccd084b2-b07e-4a5f-b25f-e89a60f65079","arxiv_id":"2509.03778","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Adding a time-relaxing, impact-parameter-dependent contact distance to the DNS-sysu master-equation model brings quasi-elastic and few-nucleon transfer cross sections into agreement with measured heavy-ion data.","lead":"This paper changes how a nuclear reaction model chooses the distance at which nucleons are exchanged, making that distance depend on how glancing the collision is. The update lets the model reproduce quasi-elastic, few-nucleon transfer channels that earlier versions underpredicted, across multiple projectile-target combinations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"QE-channel agreement hinges on an unspecified angular-momentum cutoff more than on the fitted tau_C; the central validation needs the cutoff and sensitivity shown.","rationale":"The reader's weakest_assumption focused on the fitted relaxation time tau_C. While that is a legitimate concern, the QE/grazing regime has t << tau_C, so Rcont is essentially Rclosest and tau_C has little effect on the central QE claim. The more load-bearing issue is the undocumented angular-momentum cutoff applied before comparison with experimental data. This cutoff directly shapes the QE cross sections and, if chosen flexibly, could manufacture the reported agreement. The reader did note this cutoff as an issue but did not identify it as the weakest link. My read therefore partially agrees. The manuscript has independent support: comparisons span multiple systems and energies, and the model builds on the established DNS-sysu framework. However, the missing cutoff specification and the absence of sensitivity tests leave the central claim conditional. The reader's CONDITIONAL verdict remains appropriate, so I recommend no change.","tokens_in":10137,"tokens_out":10072,"duration_ms":104856,"concrete_test":"Request the authors to report the J-cutoff used for each system and recompute Fig. 3 with J_cut varying by ±10 hbar (or equivalent impact-parameter range) and with Rtr varied by ±0.5 fm. If the 1p/2p stripping cross sections change by more than ~50%, the QE improvement is not robust. Also recompute with tau_C = 1e-22 s and 4e-22 s to verify that the QE results are independent of tau_C, as argued.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that the improved DNS-sysu model resolves quasi-elastic (QE) underestimation—is supported by comparisons in Figs. 3–5 after applying \"an angular momentum cutoff to select computational results that fall within the experimentally observable range\" (Results and Discussions). The cutoff value and procedure are not specified. Because QE/grazing cross sections are dominated by high-J partial waves, the choice of cutoff directly controls the magnitude of few-nucleon transfer cross sections. If the cutoff is chosen post hoc to improve agreement, the comparison in Fig. 3 does not test the mechanism. The relaxation time tau_C in Eq. (3) is explicitly fitted (\"can be determined from the analysis of experimental data\"), but in the QE regime t << tau_C, f(t) ~ 1, so Rcont ~ Rclosest and the results are insensitive to tau_C. The actual new input for QE is the semiclassical tail Ptr = exp(-2k[R-Rtr]) applied at Rcont = Rclosest, with Rtr = R_pro+R_tar+2.5 fm. Neither Rtr nor the functional form of Ptr is varied, and the angular-momentum selection is unreported. Thus the \"resolution\" of QE underestimation may be an artifact of the selection procedure rather than a robust consequence of the physical mechanism.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends the DNS-sysu dinuclear-system model to quasi-elastic (QE) and grazing collisions by introducing impact-parameter-dependent nucleon transfer rates. The key formal ingredients are a relaxing contact distance Rcont (Eq. 3), with relaxation time tau_C = 2e-22 s, and a semiclassical transfer-tail factor Ptr = exp(-2k[R-Rtr]) (Eq. 7). The model is compared with isotopic, mass, and charge distributions for 40Ca, 58Ni, 64Ni, 136Xe, and 208Pb projectiles on 208Pb targets, including the previously problematic few-nucleon transfer channels. The authors report that the improved model resolves the long-standing QE underestimation and describes the data reasonably across many observables.","tokens_in":10501,"tokens_out":4942,"duration_ms":52749,"significance":"If the reported mechanism is robust, the paper would make a useful contribution: it extends a widely used DNS framework into the QE/grazing regime and benchmarks it against a broad set of experimental data. Strengths include the extensive data coverage, the use of absolute cross sections, the coupling to GEMINI++ for de-excitation, and the explicit presentation of the master-equation framework. However, the central QE improvement depends on several underdetermined inputs, especially the angular-momentum cutoff used to mimic detector acceptance and the ad hoc parameters in Eqs. (3) and (7). The current manuscript does not yet establish that the improvement is a consequence of the proposed physics rather than of selection/tuning.","major_comments":[{"comment":"The text states: 'we implement an angular momentum cutoff to select computational results that fall within the experimentally observable range', but no cutoff value, selection criterion, or sensitivity is given. Since the QE/grazing cross sections are dominated by high-J partial waves, this cutoff directly controls the magnitude of the few-nucleon transfer cross sections displayed in Fig. 3. Without reporting the cutoff and showing that the conclusions are stable under reasonable variations, the agreement cannot be attributed to the proposed mechanism. Please provide the cutoff value, its experimental justification, and a sensitivity study.","section":"Results and discussions, Fig. 3"},{"comment":"The relaxation time tau_C = 2e-22 s is introduced with the statement that it 'can be determined from the analysis of experimental data', but no fitting procedure, uncertainty, or independent cross-check is given. This is load-bearing because intermediate angular momenta, and thus the DI/QF distributions, depend on tau_C. Moreover, in the QE limit t << tau_C one has f(t) ≈ 1 and Rcont ≈ Rclosest, so the QE improvement is effectively generated by replacing the old contact distance with Rclosest, not by the specific value of tau_C. The authors should separate these two effects and show how tau_C was determined and how sensitive the results are to it.","section":"Eq. (3) and following paragraph"},{"comment":"The improved QE transfer rates hinge on Ptr = exp(-2k[R-Rtr]) with Rtr = Rpro + Rtar + 2.5 fm. The value 2.5 fm and the form of Ptr are presented without justification or sensitivity analysis. Since the central claim is that the model now reproduces one- and two-proton stripping, the results should be tested against variations in Rtr (and in the nucleon separation energies entering k), and the choice should be compared with known sub-barrier transfer systematics.","section":"Eq. (7) and Fig. 2"},{"comment":"For systems without a potential pocket, Rbottom is fixed at a surface separation of approximately 0.7 fm. This is another parameter of the model, and the manuscript does not state whether this value is system-independent or fitted. If it is system-dependent, the predictive content of the comparisons in Figs. 4-5 should be qualified. Please state how Rbottom is assigned for each system and whether the results are sensitive to this offset.","section":"Theoretical framework, after Eq. (3)"}],"minor_comments":[{"comment":"Typo: 'entrace angular momentum' should be 'entrance angular momentum'.","section":"Fig. 1 caption"},{"comment":"Typo: 'The quantitiy W' should be 'The quantity W'.","section":"After Eq. (4)"},{"comment":"The notation 'C1δβ 1 2 = C2δβ 2 2' and 'δβ 1 2 + δβ 2 2 = 2 β2' is hard to parse. Please use explicit superscripts (e.g., C1δβ_1^2 = C2δβ_2^2) and define all symbols.","section":"Eq. (1)"},{"comment":"The line styles 'black dotted' and 'black solid' are described in the caption, but the distinction between the improved DNS-sysu result and the improved DNS-sysu + GEMINI++ result should be clearly visible in the figure and stated in the text for each panel.","section":"Fig. 5"},{"comment":"The abstract says the underestimation of the QE channel is 'especially for the light reaction systems', yet several comparisons involve heavy projectiles such as 136Xe and 208Pb. Please clarify the intended scope.","section":"Abstract and Introduction"}],"recommendation":"major_revision","confidential_remarks":"The central concern is not the phenomenological nature of the model per se, but the missing details of the angular-momentum cutoff. If the cutoff is chosen post hoc, the QE agreement in Fig. 3 may be a selection effect. I would like to see the cutoff value and sensitivity before judging the paper. The authors should also be asked to report how tau_C and Rtr were set for each system, since these are the new inputs that carry the claimed improvement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a real, incremental improvement to the DNS-sysu model that plausibly closes a known gap in few-nucleon transfer, but the central validation is not yet clean because the angular-momentum selection is unreported and the relaxation time is fitted without sensitivity analysis.\n\nWhat's actually new: the relaxation ansatz for the contact distance, Eq. (3), with tau_C = 2e-22 s. The master equation, PES, and transfer rates are standard, and the authors are honest that lambda0 is not derived from first principles. The systematic comparisons across 40Ca, 58Ni, 64Ni, 136Xe, and 208Pb + 208Pb are a real strength; the original-vs-improved curves in Figs. 3 and 4 show a clear difference in the few-nucleon channels, so the model is doing something.\n\nThe soft spots are load-bearing. In the Results section the authors say they applied an angular-momentum cutoff to mimic experimental detection conditions, but no value or procedure is given. QE cross sections are dominated by high partial waves, so this cutoff can control the magnitude of the few-nucleon transfer cross sections. Without the cutoff value, the agreement in Fig. 3 is not reproducible, and the comparison does not cleanly test the mechanism. The fitted tau_C is the second issue: the paper states it can be determined from data but gives no fitting procedure, uncertainty, or cross-validation. The stress-test's strong claim that tau_C is irrelevant in the QE regime is partially off: with t up to ~1e-22 s and tau_C = 2e-22 s, t/tau_C is around 0.5, so Rcont is a mix of Rclosest and Rbottom, not exactly Rclosest. But the deeper point holds: the angular-momentum cutoff is unreported, and the relative roles of Ptr and the relaxation are not disentangled. The authors should show the sensitivity of the QE cross sections to tau_C and to the cutoff, and ideally release the code or at least the cutoff values.\n\nThis is a paper for the DNS/MNT community. It deserves a serious referee, but I would not accept it as is; the missing details are central, not cosmetic.","headline":"A plausible incremental fix for quasi-elastic underestimation in the DNS-sysu model, but the central validation leans on an unreported angular-momentum cutoff and a fitted relaxation time.","tokens_in":10968,"tokens_out":3877,"would_cite":false,"duration_ms":38944,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V35"],"pacs":["25.70.Hi","24.10.-i"],"model":"deepseek-v4-flash","headline":"A relaxing contact distance lets a dinuclear-system model finally reproduce quasi-elastic transfer channels in heavy-ion collisions.","keywords":["quasi-elastic scattering","multinucleon transfer","dinuclear system model","nucleon transfer probability","heavy-ion collisions","grazing collisions","master equation","impact parameter"],"falsifier":"Extract tau_C by fitting the one- and two-proton stripping cross sections in one reaction, say 58Ni + 208Pb at Ec.m. = 256 MeV, then predict the same channels for 40Ca + 208Pb at Ec.m. = 197 MeV with that value unchanged. If no single tau_C matches both light and heavy projectile systems, the relaxation mechanism is not universal and the ansatz fails. Alternatively, an angular-momentum-resolved measurement showing the exchange distance does not stay near the turning point at grazing J would rule out the mechanism.","tokens_in":10033,"feed_emoji":"⚛️","tokens_out":6155,"duration_ms":60407,"temperature":0.7,"pith_summary":"The paper claims that dinuclear-system (DNS) models have been missing the quasi-elastic part of multinucleon transfer because they let nucleon exchange happen at a fixed contact distance deep inside the potential pocket. The authors introduce an impact-parameter-dependent contact distance that starts at the distance of closest approach for grazing trajectories and relaxes exponentially to the pocket bottom on a time scale of 2e-22 seconds. In the improved model, few-nucleon transfer at large angular momentum now receives the large transfer probability that long-range nucleon exchange needs, so the long-standing underestimation of one- and two-proton stripping disappears. The paper benchmarks the model against isotopic, mass, and charge distributions for 40Ca, 58Ni, 64Ni, 136Xe, and 208Pb beams on 208Pb and reports agreement across systems and energies. If right, the same master-equation model can describe quasi-elastic, deep-inelastic, and quasi-fission channels together.","feed_headline":"Relaxing contact distance restores quasi-elastic transfer yields","feed_subtitle":"A single exponential relaxation term brings few-nucleon transfer data back in line with theory.","key_machinery":"The relaxing contact distance of Eq. (3). Rcont is the internuclear separation at which nucleons are exchanged; the paper makes it a function of interaction time t obtained from the deflection function, interpolating between Rclosest (the distance of closest approach at each angular momentum) and Rbottom (the bottom of the potential pocket) with smoothing function f(t) = exp(-t/tau_C) and tau_C = 2e-22 seconds. This converts the transfer probability exp[-2k(Rcont - Rtr)] into an impact-parameter-dependent quantity, giving grazing collisions a large few-nucleon transfer probability while preserving the old behavior for deep-inelastic collisions.","core_discovery":"The central claim is that the mechanism behind the quasi-elastic channel in heavy-ion transfer is the time-dependent distance at which nucleon exchange takes place. Earlier DNS models assumed contact at the bottom of the potential pocket, which works for long-lived deep-inelastic configurations but fails for grazing collisions, where the interaction lasts less than about 10^-22 seconds and the nuclei barely touch. The paper proposes Rcont(t) = Rclosest exp(-t/tau_C) + Rbottom[1 - exp(-t/tau_C)] with tau_C = 2e-22 seconds, so that for large angular momenta the exchange distance stays close to the turning point and the semiclassical transfer tail exp(-2k[Rcont - Rtr]) is active; for violent co","pith_inferences":["The relaxation time tau_C is the main free parameter; a natural next test is to fit it from one dataset and see whether the same value survives across systems and bombarding energies. If it does, the model is predictive rather than interpolative.","One could derive tau_C from nuclear friction or viscosity in the entrance channel; doing so would turn a phenomenological interpolation into a dynamical prediction.","The enhancement mechanism is exponential in Rcont, so it should be sensitive to surface properties such as neutron-skin thickness; reactions with isotopes of differing neutron excess could test whether the contact-distance picture captures the tail of the single-particle density."],"forward_implications":["The DNS-sysu model can now describe quasi-elastic, deep-inelastic, and quasi-fission channels in one master-equation approach rather than treating quasi-elastic scattering by a separate model.","One- and two-proton stripping cross sections in 58Ni + 208Pb at 256 MeV, previously underpredicted, are reproduced in absolute value and slope.","The same model with the relaxation term matches measured isotopic, mass, and charge distributions for 40Ca, 58Ni, 64Ni, 136Xe, and 208Pb on 208Pb targets.","Predictions for few-nucleon transfer, a first step toward producing neutron-rich exotic nuclei, become reliable enough to guide experimental searches."],"supporting_citations":[{"why":"Foundation of the DNS-sysu model: the potential energy surface and three-dimensional master equation this work extends.","marker":"[38]"},{"why":"Supplies the nucleon transfer rate lambda0 proportional to 5 x A_tot^2 x (T/MeV) x 10^16 used in the master equation.","marker":"[39]"},{"why":"Gives the deflection function method from which the interaction time t at each angular momentum is obtained, feeding the relaxation formula.","marker":"[40]"},{"why":"Provides the macroscopic transition rate form W(S',S) proportional to exp[(U(S') - U(S))/(2T)] between adjacent states.","marker":"[48]"},{"why":"Supplies the semiclassical transfer probability exp(-2k[R - Rtr]) for separated nuclei, central to quasi-elastic nucleon exchange.","marker":"[50]"},{"why":"Supports treating nucleon exchange through tails and temporary necks between separated nuclei in grazing collisions.","marker":"[51]"},{"why":"Experimental 58Ni + 208Pb transfer data used as the main benchmark for the quasi-elastic channel improvement.","marker":"[52]"},{"why":"The GEMINI++ code is used to de-excite primary fragments before comparing calculated cross sections with measured isotopic distributions.","marker":"[53]"},{"why":"Experimental 40Ca + 208Pb isotopic distributions used to test the model at multiple center-of-mass energies.","marker":"[55]"},{"why":"Experimental 136Xe + 208Pb mass distributions benchmark many-nucleon transfer and show the relaxation term does not spoil deep-inelastic predictions.","marker":"[58]"}],"fun_headline_variants":["Time-dependent contact restores quasi-elastic transfer","Grazing collisions need dynamical contact distance","DNS model fix recovers few-nucleon transfer yields","Unified description includes quasi-elastic channel","Contact distance relaxation explains quasi-elastic scattering"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"Everything rests on the assumption that the nucleon-exchange distance relaxes exponentially from the distance of closest approach to the pocket bottom with a universal time constant of 2e-22 seconds; this exponential law is an assumed interpolation, not derived from the dynamics.","fun_headline_variants_meta":{"raw":{"variants":["Time-dependent contact restores quasi-elastic transfer","Grazing collisions need dynamical contact distance","DNS model fix recovers few-nucleon transfer yields","Unified description includes quasi-elastic channel","Contact distance relaxation explains quasi-elastic scattering"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000206,"raw_usage":{"total_tokens":1223,"prompt_tokens":722,"completion_tokens":501,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":466,"completion_tokens_details":{"reasoning_tokens":431}},"tokens_in":466,"tokens_out":501,"duration_ms":5953,"temperature":1.0,"reasoning_tokens":431,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T10:40:10.091444+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Extract tau_C by fitting the one- and two-proton stripping cross sections in one reaction, say 58Ni + 208Pb at Ec.m. = 256 MeV, then predict the same channels for 40Ca + 208Pb at Ec.m. = 197 MeV with that value unchanged. If no single tau_C matches both light and heavy projectile systems, the relaxation mechanism is not universal and the ansatz fails. Alternatively, an angular-momentum-resolved measurement showing the exchange distance does not stay near the turning point at grazing J would rule out the mechanism.","supporting_citations":[{"cited_title":"Zhu and J","cited_arxiv_id":null,"evidence_quote":"Foundation of the DNS-sysu model: the potential energy surface and three-dimensional master equation this work extends."},{"cited_title":"Saiko and A","cited_arxiv_id":null,"evidence_quote":"Supplies the nucleon transfer rate lambda0 proportional to 5 x A_tot^2 x (T/MeV) x 10^16 used in the master equation."},{"cited_title":"Wolschin and W","cited_arxiv_id":null,"evidence_quote":"Gives the deflection function method from which the interaction time t at each angular momentum is obtained, feeding the relaxation formula."},{"cited_title":"Moretto and J","cited_arxiv_id":null,"evidence_quote":"Provides the macroscopic transition rate form W(S',S) proportional to exp[(U(S') - U(S))/(2T)] between adjacent states."},{"cited_title":"Zagrebaev, Sub-barrier fusion enhancement due to neutron transfer, Physical Review C 67, 061601 (2003)","cited_arxiv_id":null,"evidence_quote":"Supplies the semiclassical transfer probability exp(-2k[R - Rtr]) for separated nuclei, central to quasi-elastic nucleon exchange."},{"cited_title":"Von Oertzen, H","cited_arxiv_id":null,"evidence_quote":"Supports treating nucleon exchange through tails and temporary necks between separated nuclei in grazing collisions."},{"cited_title":"Corradi, A","cited_arxiv_id":null,"evidence_quote":"Experimental 58Ni + 208Pb transfer data used as the main benchmark for the quasi-elastic channel improvement."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The GEMINI++ code is used to de-excite primary fragments before comparing calculated cross sections with measured isotopic distributions."},{"cited_title":"Szilner, L","cited_arxiv_id":null,"evidence_quote":"Experimental 40Ca + 208Pb isotopic distributions used to test the model at multiple center-of-mass energies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental 136Xe + 208Pb mass distributions benchmark many-nucleon transfer and show the relaxation term does not spoil deep-inelastic predictions."}],"review_version":1}