{"id":"820d0f58-6038-43a0-b375-ceb485766104","arxiv_id":"2608.08875","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In the first microscopic simulations of the pre-saddle stage of 235U(n,f), the compound nucleus thermalizes rapidly in both the ground state and isomeric wells, while octupole mass asymmetry oscillates periodically in the isomeric well when its initial value is large enough.","lead":"Using time-dependent density functional theory, this paper simulates the early stages of neutron-induced fission of uranium-235, when the compound nucleus sits in its ground state or first isomeric well before crossing the fission barrier. It finds the motion in both wells is strongly dissipative and that the nuclear shape asymmetry oscillates almost harmonically in the isomeric well, suggesting a new observable window into a stage of fission that is hard to measure directly.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Triaxiality restriction is the load-bearing soft spot: the central Q30 oscillation and dissipation claims could be altered by omitted triaxial degrees of freedom, and the paper itself flags triaxiality as important near the first barrier.","rationale":"The reader's weakest assumption matches my independent reading: the axial-symmetry restriction is the most load-bearing vulnerability in the paper's central argument. The paper is otherwise internally consistent: the dissipation and Q30 oscillation statements are direct simulation outputs, not fitted results, and the authors do not overclaim beyond their mean-field framework, explicitly calling for future beyond-mean-field studies. However, because the central qualitative difference between the two wells (chaotic Q30 in the ground-state well versus persistent harmonic Q30 in the isomeric well) is established only under axial symmetry, and because the paper itself acknowledges triaxiality matters in the first-barrier region where some trajectories originate, the conclusions are appropriately treated as conditional rather than established. The proposed 3D rerun is a concrete, feasible check that would settle whether the omission changes the central picture. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":9459,"tokens_out":5565,"duration_ms":58908,"concrete_test":"Rerun the eight isomer-well and three ground-state-well trajectories with the same energy density functional and initial densities in full 3D (no axial symmetry imposed), using the authors' LISE code or an equivalent TDDFT solver, for at least 6000 fm/c. Compare Eflow(t), Q20(t), Q30(t), and, if present, the oscillation period and critical threshold. If the large-amplitude Q30 oscillation damps within ~2000 fm/c, becomes aperiodic, or its period shifts by more than ~10% from 960±70 fm/c, the axial-symmetry assumption is material to the central claim. Also vary initial Q30 at fixed excitation energy to check that the claimed threshold is not a proxy for initial energy.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claims—strong dissipative damping in both wells and persistent large-amplitude Q30 oscillations in the isomeric well—are direct outputs of TDDFT trajectories restricted to axially symmetric shapes. The paper explicitly concedes, in the passage after Fig. 1, that 'This investigation is restricted to axially symmetric shapes, although it is known triaxiality will play an important role in fission in the region surrounding the first barrier [52,53].' Two of the eight isomer-well trajectories start near the first barrier, precisely the region where the omitted triaxial degree of freedom is known to matter. In a full triaxial TDDFT treatment, the Q30 mode can couple to gamma-type quadrupole vibrations and additional single-particle channels; this could damp the ~960 fm/c oscillation, convert it to quasiperiodic motion, or shift the claimed critical Q30 threshold. Since the paper's proposed experimental observable (the asymmetric/symmetric fission ratio as a probe of Q30 fluctuations) rests on the persistence of this oscillation, the axial restriction is load-bearing. The absence of a precisely quantified threshold and the small number of trajectories make the oscillation claim harder to validate, but the more fundamental vulnerability is the omitted triaxial degree of freedom.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper applies time-dependent density functional theory (TDDFT) with superfluid local density approximation to simulate the dynamics of 236U inside the ground-state and first isomeric wells of the potential energy surface, as a first step toward a microscopic description of the pre-saddle, compound-nucleus stage of neutron-induced fission. Ten trajectories are prepared at various deformations (three in the ground-state well, eight in the isomeric well), with initial conditions marked in Fig. 1. The key reported results are: (i) the collective flow energy, defined in the caption of Fig. 2, decays rapidly from about 0.5 MeV to about 0.1 MeV in both wells, which the authors interpret as strong dissipation and thermalization; (ii) the quadrupole moment Q20 quickly relaxes to a well-specific central value and then shows seemingly chaotic fluctuations; (iii) in the isomeric well, trajectories with a sufficiently large initial octupole moment Q30 undergo large-amplitude, almost harmonic oscillations with period T = 960 ± 70 fm/c, while ground-state-well trajectories show only chaotic Q30 fluctuations; and (iv) a small number of neutrons are emitted in both wells, with a transient burst followed by roughly linear emission. The authors propose that the ratio of asymmetric to symmetric fission serves as an experimental probe of the strength of Q30 fluctuations in the isomeric well.","tokens_in":9596,"tokens_out":3421,"duration_ms":38113,"significance":"If the central results hold, this is a useful first microscopic exploration of a stage of fission that has so far been treated only through schematic or statistical models. The reported dissipation timescale and the persistent, large-amplitude Q30 oscillation in the isomeric well would provide concrete benchmarks for more complete theories and suggest a measurable observable (the asymmetric/symmetric fission ratio) tied to pre-saddle dynamics. The study uses no fitted parameters to produce the central observations: the flow-energy decay and the Q30 oscillations are direct outputs of the TDDFT equations of motion, which is a strength. The paper also connects its findings to earlier saddle-to-scission calculations, giving a coherent picture of dissipative, non-Markovian dynamics throughout the fission process. However, the significance is tempered by the acknowledged axial-symmetry restriction, the small number of trajectories, and the qualitative basis of the 'harmonic' and 'critical threshold' claims. These issues currently limit the strength of the conclusions beyond a proof-of-principle demonstration.","major_comments":[{"comment":"The manuscript explicitly states, after Fig. 1, that the investigation is restricted to axially symmetric shapes although triaxiality is known to play an important role in the region surrounding the first barrier [52,53]. Two of the eight isomeric-well trajectories start near that barrier (the dashed-dotted red lines in Fig. 3). Since the central claims of this work concern the dissipation rate, the stability of the isomeric-well trajectories, and the period and persistence of the Q30 oscillations, the omission of triaxial degrees of freedom is load-bearing. In a full triaxial TDDFT treatment, the Q30 mode can couple to gamma-type quadrupole vibrations and additional single-particle channels, which could damp the oscillation, shift its period, or change the critical threshold. The authors should either include triaxial dynamics for at least a representative subset of trajectories or provide a quantitative argument, based on e.g. static triaxial PES curvature or coupling matrix elements, that the axial restriction does not materially alter the reported period or the dissipative decay of E_flow. As it stands, the paper's own caveat undermines the robustness of the main claims.","section":"Section 2, near Fig. 1"},{"comment":"The claim that Q30 oscillates 'in almost harmonic motion' with period T = 960 ± 70 fm/c is based on visual inspection of Fig. 3(b), and the 'critical value' of the initial Q30 separating oscillatory from chaotic behavior is never defined quantitatively. This is a central result because the proposed experimental probe (the asymmetric/symmetric fission ratio as a measure of Q30 fluctuations) depends on the persistence and regularity of this oscillation. The authors should provide a quantitative characterization, for example by computing an autocorrelation function, a power spectrum, or by fitting the time series to a damped harmonic oscillator with extracted frequency, amplitude, and damping time. They should also state how the period and its uncertainty were obtained (e.g., number of zero crossings, number of trajectories used), and define the critical initial Q30 value with an uncertainty or a range. Without this, the distinction between 'harmonic' and 'chaotic' remains subjective and cannot be independently verified.","section":"Section 3, Fig. 3(b) and surrounding text"},{"comment":"The interpretation of the decay of E_flow as 'strong dissipation' and 'rapid thermalization' needs a more careful statement of what is meant by dissipation in a unitary mean-field evolution. Equation defining E_flow (Fig. 2 caption) measures collective flow kinetic energy, and its decay reflects dephasing and Landau damping within the single-determinant TDDFT, not true entropy production or approach to a thermal ensemble as in the compound-nucleus picture. The paper does note the mean-field limitation and cites non-Markovian behavior, but the abstract and summary state that the system 'thermalizes.' I recommend either softening the language or adding a direct comparison to the single-particle occupation redistribution σ_1(t) (Fig. 4) to make explicit that the dissipation is a mean-field effect and to estimate what fraction of the initial flow energy is transferred to single-particle degrees of freedom. This is particularly important because the paper simultaneously argues that Bohr's compound-nucleus hypothesis requires beyond-mean-field complexity; the reader should not infer that the TDDFT trajectories themselves thermalize in the statistical-mechanical sense.","section":"Section 3, Fig. 2 and surrounding text"},{"comment":"The quantitative support for the central claims is thin: only three trajectories in the ground-state well and eight in the isomeric well, with the period quoted as 960 ± 70 fm/c. No error bars are shown on the E_flow curves in Fig. 2 beyond a single shaded standard-deviation band, and the statement that trajectories 'become indistinguishable, at least visually' is not backed by a quantitative metric (e.g., a distance in Q20-Q30 space or a correlation measure). Given that the phenomenon of large-amplitude octupole oscillations rests on a small number of trajectories, the paper would be substantially strengthened by reporting, for each trajectory, the fitted oscillation parameters (period, amplitude, damping time) and a measure of the variance across trajectories. If the spread is too large to allow a robust period extraction, that should be stated explicitly. This is not a request for excessive statistics, but for the minimum quantitative detail needed to evaluate the strength of the claimed effect.","section":"Section 3, Figs. 3 and 4"}],"minor_comments":[{"comment":"Reference [8] contains a typo: 'xparimental' should be 'experimental'. Reference [32] contains 'dimentions' in the title; it should be 'dimensions'.","section":"References"},{"comment":"There are several grammatical and typographical slips, including 'it's' used as a possessive pronoun (e.g., 'it's magnitude', 'it's evolution') and inconsistent comma use. A careful proofread is recommended.","section":"Throughout"},{"comment":"The definitions of Q20 and Q30 are given without factors that are often used to make them dimensionless; this is fine, but the unit labels in Fig. 1 ('b' and 'b^(3/2)') should be introduced in the text or figure caption for clarity.","section":"Section 2, Eq. (3)"},{"comment":"The caption states 'The color bar denotes the difference in energy of a point on the PES and the ground state,' but in the printed figure the color bar is not visible or legible in panels (a)-(d). Please ensure the color scale is reproducible and, if possible, add contour labels or a separate legend.","section":"Section 2, Fig. 1 caption"},{"comment":"The description of the neutron emission criterion ('17 fm or more beyond the center of mass along the fission axis') would benefit from a brief justification, since the result is stated in the abstract as a notable finding and the interpretation of the early-time jump as removal of constraints is only a hypothesis.","section":"Section 3, Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a fresh and potentially interesting application of TDDFT to a largely unexplored stage of fission, and the observed dissipation and Q30 oscillations are direct outputs of the model with no fitted parameters. The main concern is that the axial-symmetry restriction, which the authors themselves flag as important in the barrier region, directly affects the two most novel claims (dissipation timescale and persistent octupole oscillation). Two of eight isomer-well trajectories start near the first barrier, so the effect is not merely a distant caveat. In addition, the 'harmonic' and 'critical threshold' characterizations are too qualitative for the claims to be fully assessed. I would like to see the authors either add triaxial benchmarks, provide a quantitative oscillation analysis, or explicitly reframe the paper as a proof-of-principle study whose conclusions are conditional on the axial symmetry assumption. The self-citation pattern is not excessive for this group, but the paper relies heavily on the authors' earlier work for the method; that is appropriate given the novelty of the approach."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first microscopic picture of what happens in the ground-state and isomer wells before a fissioning nucleus reaches the outer saddle. The headline observation—large-amplitude, roughly harmonic Q30 oscillations in the isomer well when the initial octupole deformation is large enough—is new and, if it holds up, gives a mechanism for mass asymmetry that is generated pre-saddle and persists. The dissipation claim (flow energy dropping from about 0.5 to about 0.1 MeV) is expected but has not been shown directly before.\n\nThe paper does a few things right. The trajectories are direct outputs of TDDFT, not fitted to reproduce the claims. The period is consistent across the eight isomer-well runs. The authors also honestly flag the main limitation: everything is axially symmetric, even though triaxiality is known to matter near the first barrier. That is exactly where two of the eight trajectories start.\n\nThe soft spots are real but not fatal. The 'harmonic' label comes from visually inspecting Fig. 3; no quantitative measure of oscillation quality is given. The critical Q30 threshold is not defined. With eleven trajectories total, the scatter is large enough that I would not bet on the 960 fm/c period surviving a triaxial calculation without change. And no code or data are released, which makes it hard to check the flow-energy decay and neutron counts.\n\nStill, the central qualitative picture is coherent: mean-field dynamics in both wells is strongly dissipative, and the isomeric well is special because it can sustain a large oscillating octupole moment. The proposed experimental connection—asymmetric-to-symmetric yield ratio as a probe of Q30 fluctuations in the isomer well—is speculative but worth stating, and the authors label it as such.\n\nI would send this out for review. It is exploratory, but it opens a stage of fission that has not been studied microscopically, and the authors are candid about what would need to be redone (triaxiality, beyond-mean-field, more trajectories). The referee should ask for a quantitative definition of the oscillation threshold, a test of sensitivity to the initial conditions, and ideally a triaxial run or two. I would not require code release, but it would help.","headline":"First microscopic look at pre-saddle fission well dynamics, with a notable harmonic octupole oscillation in the isomer well; the axial-symmetry restriction is the main thing to test before the claims become robust.","tokens_in":10193,"tokens_out":2383,"would_cite":true,"duration_ms":23841,"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":"The paper claims that 236U formed by neutron capture dissipates its collective energy to ~0.1 MeV in both pre-saddle wells, and that in the isomeric well the octupole moment oscillates with a 960 fm/c period, generating persistent mass…","keywords":["time-dependent density functional theory","nuclear fission","compound nucleus","collective dissipation","isomeric well","octupole oscillation","mass asymmetry","235U neutron-induced fission"],"falsifier":"Evolve the same isomeric-well initial configurations with a triaxial TDDFT solver that allows all three deformation degrees of freedom; if the octupole moment no longer oscillates with period 960±70 fm/c beyond the critical initial value, or the flow-energy decay rate changes by more than about a factor of two, the paper's central picture fails.","tokens_in":9184,"feed_emoji":"⚛️","tokens_out":7651,"duration_ms":74576,"temperature":0.7,"pith_summary":"This paper uses time-dependent density functional theory to follow the compound nucleus 236U during the two earliest stages of neutron-induced fission: the ground-state well and the first isomeric well, which together dominate the pre-saddle lifetime. It claims that collective motion in both wells is strongly dissipative, with the collective flow energy dropping from below 0.5 MeV to about 0.1 MeV, so the freshly formed system thermalizes quickly and loses memory of its initial deformation. It also claims a qualitative difference between the wells: in the isomeric well, if the initial octupole moment exceeds a threshold, it oscillates harmonically with period 960±70 fm/c, generating a mass asymmetry that persists; in the ground-state well the octupole signal is chaotic. If correct, this is the first microscopic demonstration that dissipation precedes the outer fission barrier and that mass asymmetry can be born before the saddle. The result would connect microscopic dynamics to the compound-nucleus picture and give an experimental handle via the ratio of asymmetric to symmetric fission yields.","feed_headline":"Fission thermalizes before the outer barrier","feed_subtitle":"In 236U, both pre-saddle wells dissipate to ~0.1 MeV and octupole shape oscillates with 960 fm/c period.","key_machinery":"The machinery is time-dependent superfluid density functional theory (TDDFT), evolving a single generalized Slater determinant in a $30^{2}$×60 fm³ box, restricted to axially symmetric shapes. The key observables are the quadrupole moment Q20 and octupole moment Q30 (defined as spatial moments of the total nucleon density), the collective flow energy E_flow = Σ_{q=n,p} ∫ d³r ħ² j_q²/(2m), and σ₁(t) = Σ_k |n_k(t) − n_k(0)|, which tracks the redistribution of single-particle occupation numbers. These tools let the authors see whether the cold initial deformation, a quantum mean-field state formed when a low-energy neutron is absorbed, is quickly forgotten or develops persistent collective patterns.","core_discovery":"The central discovery is that the two pre-saddle wells are not quiet waiting rooms: in both the ground-state and first-isomeric wells of 236U, the collective flow energy—initially below 500 keV—decays to about 100 keV on a timescale of a few thousand fm/c, showing that the mean-field dynamics dissipates the energy of the collective motion. At the same time, the quadrupole moment Q20 relaxes to a common central value in each well, and the single-particle occupation numbers keep redistributing approximately linearly in time, a non-Markovian signature. The notable result is in the isomeric well: when the initial octupole moment Q30 is large enough, it oscillates with a large amplitude and period T = 960 ± 70 fm/c, i.e., mass asymmetry is generated, persists, and behaves almost harmonically; below a critical value of Q30, the motion is instead mostly chaotic. The paper interprets this as evidence that the compound nucleus thermalizes during the early stages, that tunneling through the barriers would be hindered by dissipation, and that the ratio of asymmetric to symmetric fission yields can probe the Q30 fluctuations in the isomeric well.","pith_inferences":["A triaxial TDDFT calculation in the isomeric well is the natural next test: if the ~960 fm/c octupole oscillation disappears when non-axial shapes are allowed, the harmonic persistence would be a consequence of the axial restriction rather than a robust property of 236U.","The predicted threshold in Q30 implies that the ratio of asymmetric to symmetric fission yields should depend on how the system enters the isomeric well, for instance on the excitation energy of the compound nucleus; existing yield systematics could be examined for such a correlation.","The pre-saddle neutron emission seen in the wells could be included in statistical models of fission as an additional pre-scission neutron source, affecting predicted neutron multiplicities.","The fast thermalization time suggests that entrance-channel memory is lost long before the outer saddle, strengthening the statistical treatment of fission fragment distributions but also predicting that detailed entrance-channel effects are washed out early."],"forward_implications":["The compound-nucleus picture acquires microscopic support: the freshly formed 236U loses its initial deformation and thermalizes in both pre-saddle wells.","Dissipation in the wells implies that quantum tunneling through the fission barriers is hindered, lengthening the pre-saddle lifetime and making spontaneous fission slower.","The ratio of asymmetric to symmetric fission yields becomes a possible experimental probe of octupole fluctuations in the isomeric well, a stage that is otherwise inaccessible to direct measurement.","The roughly linear growth of the occupation-redistribution measure σ₁(t) in both wells extends the non-Markovian, dissipative character of fission from saddle-to-scission back to the earliest stages.","Low-probability neutron emission from the ground-state and isomeric wells contributes a pre-scission neutron component, though the paper estimates the numbers are conservative."],"supporting_citations":[{"why":"Defines the TDDFT framework for fission and establishes the saddle-to-scission stage that this paper extends backward.","marker":"[1]"},{"why":"Prior demonstration of strong dissipation in the saddle-to-scission dynamics of 240Pu, the baseline the wells are compared against.","marker":"[14]"},{"why":"The LISE solver used to evolve the time-dependent superfluid density functional equations in three dimensions.","marker":"[32]"},{"why":"Shows Lyapunov exponents in TDDFT are negligible, supporting the interpretation that the observed fluctuations are not classical chaos.","marker":"[33]"},{"why":"Caldeira-Leggett result invoked to argue that the measured dissipation would hinder tunneling through the fission barriers.","marker":"[17]"},{"why":"Provides the quantum extension of the collision term that describes how TDDFT captures dissipation beyond Landau damping.","marker":"[24]"},{"why":"One of the references establishing that triaxiality plays an important role near the first barrier, the main caveat to the axial-symmetry restriction.","marker":"[52]"},{"why":"Another reference on the importance of triaxiality near the first barrier and on numerical accuracy of mean-field calculations.","marker":"[53]"}],"fun_headline_variants":["Pre-saddle wells dissipate energy in 236U","Isomeric well oscillates with 960 fm/c period","236U fission wells dissipate and oscillate","Both wells thermalize before the outer barrier","Neutron emission from early wells of 236U"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire study assumes axially symmetric nuclear shapes, even though the paper itself notes that triaxiality is expected to matter near the first fission barrier, so the reported dissipation rates and the isomeric-well octupole oscillations could change materially when triaxial shapes are allowed.","fun_headline_variants_meta":{"raw":{"variants":["Pre-saddle wells dissipate energy in 236U","Isomeric well oscillates with 960 fm/c period","236U fission wells dissipate and oscillate","Both wells thermalize before the outer barrier","Neutron emission from early wells of 236U"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.002039,"raw_usage":{"total_tokens":7968,"prompt_tokens":995,"completion_tokens":6973,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":611,"completion_tokens_details":{"reasoning_tokens":6897}},"tokens_in":611,"tokens_out":6973,"duration_ms":51526,"temperature":1.0,"reasoning_tokens":6897,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:21:18.289131+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evolve the same isomeric-well initial configurations with a triaxial TDDFT solver that allows all three deformation degrees of freedom; if the octupole moment no longer oscillates with period 960±70 fm/c beyond the critical initial value, or the flow-energy decay rate changes by more than about a factor of two, the paper's central picture fails.","supporting_citations":[{"cited_title":"Fission dynamics of 240Pufrom saddle to scission and beyond,","cited_arxiv_id":null,"evidence_quote":"Prior demonstration of strong dissipation in the saddle-to-scission dynamics of 240Pu, the baseline the wells are compared against."},{"cited_title":"The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimentions,","cited_arxiv_id":null,"evidence_quote":"The LISE solver used to evolve the time-dependent superfluid density functional equations in three dimensions."},{"cited_title":"Sen- sitivity of time-dependent density functional theory to initial conditions,","cited_arxiv_id":null,"evidence_quote":"Shows Lyapunov exponents in TDDFT are negligible, supporting the interpretation that the observed fluctuations are not classical chaos."},{"cited_title":"Quantum tunnelling in a dissipative system,","cited_arxiv_id":null,"evidence_quote":"Caldeira-Leggett result invoked to argue that the measured dissipation would hinder tunneling through the fission barriers."},{"cited_title":"Pure quantum extension of the semiclassical Boltzmann-Uehling-Uhlenbeck equation,","cited_arxiv_id":null,"evidence_quote":"Provides the quantum extension of the collision term that describes how TDDFT captures dissipation beyond Landau damping."},{"cited_title":"Equilibrium configurations of rotating charged or gravitating liquid masses with surface tension. ii,","cited_arxiv_id":null,"evidence_quote":"One of the references establishing that triaxiality plays an important role near the first barrier, the main caveat to the axial-symmetry restriction."},{"cited_title":"Numerical accuracy of mean-field calculations in coordinate space,","cited_arxiv_id":null,"evidence_quote":"Another reference on the importance of triaxiality near the first barrier and on numerical accuracy of mean-field calculations."}],"review_version":1}