{"id":"f08c1337-d25f-4e8d-9e1a-cc07a3945168","arxiv_id":"2606.10773","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"δNEGF reformulates NEGF using quantum fluctuations of field operators to represent correlations, enabling stable linear-time simulations of GW and T-matrix dynamics up to 10^4 basis states.","lead":"This paper introduces δNEGF, a new fluctuation-based version of nonequilibrium Green functions that represents two-particle correlations via operator fluctuations to achieve stable, time-linear scaling for large systems. If it works, it could let researchers simulate correlated dynamics in systems with thousands of states, such as 2D materials and nanostructures, that were previously inaccessible.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Whether the stochastic low-rank decomposition of δĜ fluctuations preserves positivity and dynamical correlations without uncontrolled errors at Nb∼10^4 and strong coupling","rationale":"The reader's weakest_assumption directly isolates the same point. Because the full manuscript is now available, the concrete_test above can be executed on the existing benchmark data and code; if it passes, the central claim is internally consistent on this axis. No other internal inconsistency (e.g., in the time-linear scaling derivation or the mapping to GW/T-matrix self-energies) appears load-bearing from the abstract and stated claims.","tokens_in":1877,"tokens_out":397,"duration_ms":13446,"concrete_test":"Take the 2D Hubbard lattice benchmark (already performed at small Nb) and recompute the same trajectory at fixed stochastic rank r=32 while increasing Nb from 100 to 400; monitor the smallest eigenvalue of the one-body reduced density matrix at each time step. If any eigenvalue drops below −10^{-8} or the deviation from the exact/HF-GKBA reference grows faster than O(1/√r), the positivity/stability guarantee does not survive the stochastic approximation at larger basis sizes.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on representing two-particle correlations via fluctuations δĜ of field-operator products, which is asserted to guarantee positivity of reduced density matrices and thus stable propagation, while the stochastic low-rank decomposition reduces storage and cost to enable Nb∼10^4. This combination is what allows extension beyond the G1-G2 scheme's Nb∼10^2 limit. The load-bearing step is whether the stochastic sampling of the δĜ ensemble, when truncated at low rank, still reproduces the exact positivity properties and the essential dynamical content of the two-particle correlations for inhomogeneous systems; any violation or truncation error here would undermine both stability and accuracy claims even if the formal δNEGF equations are correct.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript introduces δNEGF, a quantum-fluctuation formulation of nonequilibrium Green functions that represents two-particle correlations via fluctuations δĜ of field-operator products. This is claimed to guarantee stable dynamics through positivity preservation of reduced density matrices, avoid explicit two-particle Green function storage, and reduce propagation to an ensemble of Hartree-Fock-like trajectories. Combined with stochastic low-rank decomposition of the correlation functions, the approach retains time-linear scaling and extends dynamical GW and particle-particle/particle-hole T-matrix simulations to Nb∼10^4. Benchmarks against exact and HF-GKBA results on lattice systems are reported to show stable correlated dynamics at strong coupling, with demonstrations on diffusion in 2D Hubbard lattices and ultrafast relaxation in graphene nanoribbon heterostructures.","tokens_in":2014,"tokens_out":508,"duration_ms":21136,"significance":"If the stochastic low-rank approximation to δĜ preserves positivity and essential dynamical correlations without uncontrolled errors, the method would represent a substantial advance by overcoming the Nb∼10^2 limit of the G1-G2 scheme and enabling scalable NEGF simulations for large inhomogeneous systems. The time-linear scaling and extension to strong coupling are notable strengths.","major_comments":[{"comment":"The central stability claim rests on the δĜ representation guaranteeing positivity of reduced density matrices even after stochastic low-rank truncation; however, the manuscript provides no explicit derivation or numerical test quantifying positivity violations or truncation errors as a function of rank and Nb in the benchmark comparisons to exact results.","section":"Method and benchmarks"},{"comment":"The extension to Nb∼10^4 for inhomogeneous systems with long-range interactions is load-bearing for the main result, yet the abstract and reported demonstrations lack quantitative error analysis (e.g., deviation from exact or HF-GKBA data at strong coupling) that would confirm the low-rank decomposition captures the essential two-particle correlations without significant uncontrolled approximations.","section":"Results on lattice systems and applications"}],"minor_comments":[{"comment":"Notation for δĜ and the stochastic ensemble should be defined more explicitly at first use to improve readability for readers unfamiliar with fluctuation-based formulations.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The manuscript appears to fit the scope of cond-mat.str-el, but the absence of detailed error metrics in the provided text raises questions about whether the benchmarks sufficiently support the scalability claims for the targeted journal audience."},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful reading, positive assessment of the work's significance, and constructive major comments. We address each point below and will implement revisions to strengthen the presentation of the method and results.","responses":[{"response":"We agree that an explicit treatment of positivity under truncation is needed. The exact δĜ dynamics preserve positivity of the reduced density matrices by construction, as the fluctuation operators are defined directly from the full many-body state (see Sec. II). For the stochastic low-rank case we will add a short derivation in the methods section showing the conditions for approximate preservation, together with numerical tests in the benchmarks (new panel or appendix) that report the lowest eigenvalue of the one-particle density matrix versus rank for the lattice systems, directly compared to exact results. This will quantify any violations as a function of rank and Nb.","revision_made":"yes","referee_comment":"[Method and benchmarks] The central stability claim rests on the δĜ representation guaranteeing positivity of reduced density matrices even after stochastic low-rank truncation; however, the manuscript provides no explicit derivation or numerical test quantifying positivity violations or truncation errors as a function of rank and Nb in the benchmark comparisons to exact results."},{"response":"The referee correctly notes the value of quantitative error metrics. While the manuscript already shows qualitative stability and agreement with exact/HF-GKBA data for small systems, we will revise the abstract to include a brief statement on observed error magnitudes and add a table or figure panel in the results section reporting relative deviations in observables (density, energy, or correlation functions) versus rank and interaction strength for the benchmark lattices. For the Nb∼10^4 demonstrations we will include rank-convergence checks where feasible.","revision_made":"yes","referee_comment":"[Results on lattice systems and applications] The extension to Nb∼10^4 for inhomogeneous systems with long-range interactions is load-bearing for the main result, yet the abstract and reported demonstrations lack quantitative error analysis (e.g., deviation from exact or HF-GKBA data at strong coupling) that would confirm the low-rank decomposition captures the essential two-particle correlations without significant uncontrolled approximations."}],"tokens_in":1517,"tokens_out":465,"duration_ms":23726,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The new element is a fluctuation formulation (δNEGF) that encodes two-particle correlations through δĜ operator fluctuations instead of propagating the full two-particle Green function. Paired with stochastic low-rank decomposition this keeps linear time scaling and reaches basis sizes around 10^4, well beyond the G1-G2 limit of roughly 100. They benchmark against exact results and HF-GKBA on lattice models, report stable runs at strong coupling, and show example calculations for diffusion on 2D Hubbard lattices plus relaxation in graphene nanoribbon heterostructures with long-range interactions.\n\nThe positivity argument for stability is plausible because the fluctuation representation is constructed to respect reduced-density-matrix positivity. The demonstrations on inhomogeneous systems are the clearest practical gain.\n\nThe soft spot is exactly the stochastic low-rank step. The abstract asserts that the sampled ensemble plus truncation still captures the essential dynamical content without uncontrolled errors, but no error bounds, convergence tests with rank, or checks on positivity violation at large Nb and strong coupling are visible here. If the truncation introduces systematic bias or occasional negativity, both the stability claim and the accuracy at Nb~10^4 would be compromised even if the formal equations are correct. The benchmarks are mentioned but not quantified in the abstract, so it is hard to judge how far the method has actually been pushed.\n\nThis is aimed at groups already running real-time NEGF or GKBA calculations who hit the basis-size wall on inhomogeneous correlated systems. It is worth sending to referees because the scaling target is concrete and the formulation is distinct from prior work; a serious review would need to see the derivations, the rank-convergence data, and larger-system error metrics.","headline":"δNEGF scales NEGF to Nb~10^4 via fluctuations and stochastic low-rank but the truncation's effect on positivity and accuracy at strong coupling remains the key untested piece.","tokens_in":2512,"tokens_out":400,"would_cite":false,"duration_ms":13128,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"A fluctuation-based reformulation of nonequilibrium Green functions enables stable correlated simulations up to basis sizes of order 10,000 while retaining linear time scaling.","keywords":["nonequilibrium Green functions","correlated dynamics","stochastic decomposition","two-particle correlations","Hubbard model","graphene nanoribbons","large-scale simulations","positivity preservation"],"falsifier":"A calculation on a lattice system with roughly 1000 basis states at strong coupling where δNEGF results diverge from exact or HF-GKBA benchmarks by more than a few percent in observables such as double occupancy or current.","tokens_in":2755,"feed_emoji":"⚛️","tokens_out":714,"duration_ms":12502,"temperature":0.7,"pith_summary":"The paper presents δNEGF, a quantum-fluctuation version of nonequilibrium Green functions that encodes two-particle correlations as fluctuations δĜ of field-operator products rather than storing the full two-particle function. This change preserves positivity of reduced density matrices and converts the propagation into an ensemble of Hartree-Fock-like trajectories. When paired with a stochastic low-rank decomposition, the scheme keeps the favorable linear scaling in the number of time steps and extends dynamical GW and T-matrix approximations to systems with Nb approximately 10,000. Benchmarks on lattice models confirm stable dynamics even at strong coupling, and the method is applied to diffusion on two-dimensional Hubbard lattices and ultrafast relaxation in graphene nanoribbon heterostructures with long-range interactions.","feed_headline":"Fluctuation reformulation scales NEGF to 10,000 basis states","feed_subtitle":"δNEGF encodes correlations as operator fluctuations to keep linear time scaling and stable dynamics at strong coupling.","key_machinery":"The quantum-fluctuation formulation that encodes two-particle correlations as fluctuations δĜ of field-operator products, together with stochastic low-rank decomposition.","core_discovery":"δNEGF represents dynamical two-particle correlations through fluctuations of field-operator products δĜ, which guarantees stable dynamics by preserving positivity of the reduced density matrices, avoids explicit storage of the two-particle Green function, and reduces propagation to a finite ensemble of Hartree-Fock-like trajectories; combined with stochastic low-rank decomposition this retains time-linear scaling and extends GW and T-matrix simulations to Nb∼10^4.","pith_inferences":["The positivity-preserving property may allow δNEGF to serve as a stable platform for embedding other self-energy approximations beyond GW and T-matrix.","The reduction to an ensemble of single-particle trajectories suggests possible hybrid schemes that combine δNEGF with classical or semiclassical sampling methods for even larger systems.","Extension to three-dimensional or disordered geometries could be tested by applying the same stochastic decomposition to systems with broken translational symmetry."],"forward_implications":["Dynamical GW and particle-particle or particle-hole T-matrix simulations become feasible for basis sizes up to order 10,000.","Time propagation remains linear in the number of time steps.","Stable correlated dynamics are obtained even at strong coupling.","Large-scale simulations of diffusion in two-dimensional Hubbard lattices and ultrafast relaxation in graphene nanoribbon heterostructures are demonstrated."],"fun_headline_variants":["δNEGF uses operator fluctuations for stable large system dynamics","Correlations encoded as δĜ fluctuations enable 10000 basis NEGF","Ensemble of Hartree Fock trajectories scales NEGF to ten thousand states","Stochastic low rank method extends T matrix NEGF to Nb 10^4"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The representation of two-particle correlations through fluctuations of field-operator products together with the stochastic low-rank decomposition accurately captures the essential dynamical correlations and preserves positivity without significant uncontrolled errors.","fun_headline_variants_meta":{"raw":{"variants":["δNEGF uses operator fluctuations for stable large system dynamics","Correlations encoded as δĜ fluctuations enable 10000 basis NEGF","Ensemble of Hartree Fock trajectories scales NEGF to ten thousand states","Stochastic low rank method extends T matrix NEGF to Nb 10^4"]},"model":"grok-4.3","cost_usd":0.005552,"raw_usage":{"total_tokens":2714,"prompt_tokens":771,"num_sources_used":0,"completion_tokens":76,"cost_in_usd_ticks":55524500,"prompt_tokens_details":{"text_tokens":771,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1867,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":771,"tokens_out":76,"duration_ms":14518,"temperature":1.0,"reasoning_tokens":1867,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-27T11:41:47.944511+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A calculation on a lattice system with roughly 1000 basis states at strong coupling where δNEGF results diverge from exact or HF-GKBA benchmarks by more than a few percent in observables such as double occupancy or current.","supporting_citations":[],"review_version":1}