{"id":"7e1487ba-2c13-42e6-90e9-8ea39c560b47","arxiv_id":"2607.07325","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":2,"one_line_summary":"A bulk-driven nonlinear KMP model is shown to exhibit controllable nonlinear energy diffusion and programmable time-crystal phases via packing fields.","lead":"This paper introduces a bulk-driven generalization of the KMP heat conduction model where biased microscopic energy collisions induce a macroscopic current. It demonstrates that this mechanism can be used to control nonlinear energy transport and induce time-crystalline phases with traveling energy condensates.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Local equilibrium approximation is the weakest link, but it is standard, well-motivated, and numerically validated; no more significant concern identified.","rationale":"The reader's identification of the LE approximation as the weakest assumption is correct and well-calibrated. This is the standard closure for KMP-type hydrodynamics, justified by time-scale separation in the diffusive scaling limit, and validated numerically in the paper. The derivation from the exact master equation (Eq. 6) through the balance equation (Eq. 11) to the constitutive relation (Eq. 23) is clean and internally consistent. The Einstein relation σ = 2ρ²D emerges naturally as a consistency check. The time-crystal results are an application of an established framework to the new model, with the β-independence of λ_c providing a nontrivial falsifiable prediction that is confirmed. The paper is appropriately scoped as a proof-of-concept study. No adjustment to the ACCEPT verdict is warranted.","tokens_in":16305,"tokens_out":2515,"duration_ms":89934,"concrete_test":"To further stress-test the LE approximation beyond the constant-field regime, measure the two-point local energy correlation function ⟨ϵ_i ϵ_{i+1}⟩_st under constant driving (E ≠ 0, homogeneous profile) and compare with the LE prediction ⟨ϵ_i⟩⟨ϵ_{i+1}⟩ = ρ₀². If the ratio ⟨ϵ_i ϵ_{i+1}⟩/ρ₀² deviates from 1 by more than a few percent at L = 200, subleading corrections to the LE measure may be non-negligible at the system sizes simulated, potentially affecting the quantitative agreement in Fig. 1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader correctly identifies the local equilibrium (LE) approximation (Eq. 16) as the weakest assumption. The derivation of the constitutive relation (Eq. 23) and transport coefficients (Eqs. 24–25) depends on replacing the true nonequilibrium stationary measure with a locally Gibbsian product of exponentials. This is the standard hydrodynamic closure for KMP-type models and has been validated for the symmetric nonlinear KMP variant in prior work [55]. In the present driven case, the bias is weak (χ = E/L), so deviations from the equilibrium product measure are O(L⁻¹) and subleading corrections do not affect the leading-order hydrodynamics — this argument is sound. The paper provides direct numerical validation: Fig. 1 shows excellent agreement between the predicted mobility σ(ρ) = (β+3)/3 · ρ^(β+2) and simulations across multiple β and E values. The Einstein relation σ = 2ρ²D provides an additional internal consistency check. The time-crystal results build on the same group's prior framework [64–67]; the critical coupling λ_c = 2πm/ρ₀ (Eq. 31) is β-independent, a nontrivial prediction confirmed by the finite-size scaling in Fig. 2. I do not find a more significant concern than the LE approximation, which is already flagged by the reader and is appropriately scoped as a proof-of-concept study.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"The paper introduces a bulk-driven nonlinear generalization of the KMP model of stochastic energy diffusion. The microscopic dynamics is modified in two ways: (i) the collision kernel depends on the total pair energy via a function π(ν), and (ii) the redistribution parameter α is drawn from a biased distribution f_χ(α|ν) derived from a local detailed balance condition, inducing a net energy current. Starting from the exact master equation, the authors derive the hydrodynamic constitutive relation j = −D(ρ)∂_xρ + σ(ρ)E using a local equilibrium (LE) approximation, obtaining closed-form expressions for D(ρ) and σ(ρ) in terms of π(ν). For power-law kernels π(ν) ∝ ν^β, these yield D(ρ) = (β+3)/6 · ρ^β and σ(ρ) = (β+3)/3 · ρ^{β+2}, satisfying an Einstein relation σ = 2ρ²D. Kinetic Monte Carlo simulations validate the mobility prediction across multiple β and E values. As a proof of concept, the authors apply a packing field coupled to energy-field fluctuations and demonstrate a second-order phase transition to a time-crystal phase with traveling energy condensates, consistent with the predicted critical coupling λ_c = 2πm/ρ₀.","tokens_in":16997,"tokens_out":1683,"duration_ms":157442,"significance":"The paper makes a solid contribution to the nonequilibrium statistical mechanics of transport models. The derivation of the transport coefficients (Eqs. 24–25) is parameter-free given the microscopic kernel π(ν), and the Einstein relation σ = 2ρ²D provides an internal consistency check. The numerical validation in Fig. 1 is quantitative, testing the predicted mobility σ(ρ) = (β+3)/3 · ρ^{β+2} across three β values and three field strengths. The critical coupling λ_c = 2πm/ρ₀ (Eq. 31) is a nontrivial, β-independent prediction confirmed by finite-size scaling in Fig. 2. The programmable time-crystal phases (Fig. 4) serve as an effective demonstration of the control afforded by the bulk-driving mechanism. The framework is clearly situated within the established KMP/hydrodynamic literature.","major_comments":[{"comment":"§III, Eqs. (19)–(22): The derivation of the gradient and drift contributions to the current relies on expanding the LE measure to O(L⁻¹) (Eq. 18) and computing two integrals. The O(L⁰) term is stated to vanish by symmetry. However, the manuscript does not explicitly show or verify that the O(L⁻¹) corrections to the LE measure itself (which arise from the weak field χ = E/L and are distinct from the spatial gradient correction already in Eq. 18) do not contribute additional terms at the same order. The text argues that subleading corrections are O(L⁻¹) and do not affect leading hydrodynamics, but the drift term being computed is itself O(L⁻¹), so this argument requires more care. A brief justification that field-induced corrections to the product measure are absent at the order retained would strengthen the derivation.","section":null},{"comment":"§IV.B, Eq. (31): The critical coupling λ_c^{(m)} = 2πm/ρ₀ is stated to be β-independent, which is presented as a nontrivial prediction confirmed by the inset of Fig. 2. However, the inset shows ⟨|z₁|⟩ vs. λ/λ_c^{(1)} for different β, and since λ_c itself is β-independent, plotting against the reduced coupling λ/λ_c by construction removes the β-dependence from the horizontal axis. The β-dependence visible in the inset reflects differences in the order parameter growth above criticality, not in λ_c itself. The claim that this 'confirms' the β-independence of λ_c would be more convincing with an explicit measurement of the critical point as a function of β (e.g., via a Binder-type crossing analysis), rather than relying on the reduced-coupling plot. This is a moderate concern because the β-independence follows analytically from Eq. (31) and the transport coefficients (28), so the claim is,","section":null},{"comment":"§IV.B, Fig. 2 (main panel): The finite-size scaling shows curves for L = 100–1600 converging toward a step-like transition, but the data near the critical point (λ/λ_c ≈ 1) do not show a clear crossing or data-collapse analysis that would quantitatively locate the critical coupling. Given that the paper claims to 'confirm the predicted critical threshold λ_c^{(1)}' (Eq. 31), a more rigorous finite-size scaling analysis (e.g., plotting L^{1/ν}⟨|z₁|⟩ vs. (λ−λ_c)L^{1/(νβ_c)} or a Binder cumulant crossing) would substantially strengthen this claim. As it stands, the visual agreement between the apparent transition region and the predicted λ_c is suggestive but not quantitatively decisive.","section":null}],"minor_comments":[{"comment":"§III, Eq. (16): The LE approximation is identified as the weakest assumption. The paper acknowledges this is a proof-of-concept study, but a brief discussion of the conditions under which LE is expected to hold for the driven case (beyond the O(L⁻¹) argument) would help the reader.","section":null},{"comment":"§IV.A: Only the mobility σ(ρ) is validated numerically; the diffusivity D(ρ) is stated to have been validated in prior boundary-driving work [55]. A brief note on the expected accuracy of D(ρ) in the driven case, or a cross-check via the Einstein relation using the measured σ(ρ), would improve transparency.","section":null},{"comment":"§IV.B, Eq. (31): The β-independence of λ_c is a striking result. A sentence explaining physically why the β-dependence of D and σ cancels exactly in the ratio would help the reader understand the mechanism.","section":null},{"comment":"Fig. 4: The raster plots are effective but the color scale and axis labels could be more clearly annotated. In particular, panel (c) shows a dynamic superposition of packing fields, but the connection between the modulation in panel (d) and the observed pattern is not immediately obvious from a casual reading of the text.","section":null},{"comment":"§V: The connection to the Doob transform and rare-event physics is mentioned as a future direction. This is an interesting idea but is currently undeveloped; a brief sentence on the specific scenario envisioned would help readers unfamiliar with that connection.","section":null},{"comment":"The paper builds substantially on the authors' prior work on packing fields and time crystals [64–67]. While the novelty of the time-crystal results per se is limited (the mechanism is transferred from prior work on other models), the application to the nonlinear KMP framework with bulk driving is a legitimate extension. The authors should ensure that the incremental nature of this contribution relative to [64–67] is clearly stated.","section":null},{"comment":"Typo: §II, the sentence containing 'the normalization factor Ω(ϵ_n) ≡ L⁻¹ Σ π(ν_{l,n})' could benefit from clarifying that the L⁻¹ is part of the definition of Ω (making it an average), not the sum itself.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is a competent extension of the authors' established program on KMP-type models and packing-field time crystals. The hydrodynamic derivation is the most novel element; the time-crystal application is largely a transfer of prior results to a new model class. The LE approximation concern is standard for this type of work and is adequately scoped as proof-of-concept. The main quantitative gap is the lack of a rigorous finite-size scaling analysis for the critical coupling claim, which is why I recommend minor revision rather than acceptance. The self-citation pattern is heavy but reflects a genuine research program; I do not see evidence of inappropriate exclusion of competing work."},"author_rebuttal":null,"desk_editor":{"model":"glm-5.2","letter":"Bottom line: this paper derives the hydrodynamic transport coefficients for a bulk-driven nonlinear KMP model, with no free parameters, and validates them in simulation. The time-crystal part is a proof of concept that works but leans on the authors' prior framework rather than breaking new theoretical ground there. It's a clean, competent paper worth a serious referee. The genuinely new result is the bulk-driving mechanism itself. The authors bias the local energy redistribution rule (Eq. 5) via local detailed balance, derive the constitutive relation (Eq. 23) and transport coefficients (Eqs. 24-25) from the master equation using a local equilibrium closure, and get closed-form expressions for power-law kernels. The mobility sigma(rho) = (beta+3)/3 * rho^(beta+2) is parameter-free given the microscopic kernel, and Fig. 1 shows quantitative agreement across multiple beta and E values. The Einstein relation sigma = 2*rho^2*D provides an internal consistency check. This is real, falsifiable work. The time-crystal section applies the packing-field machinery from the same group's prior papers [64-67] to this new model class. The critical coupling lambda_c = 2*pi*m/rho_0 (Eq. 31) is derived from the transport coefficients, not fitted, and its beta-independence is a nontrivial prediction confirmed by finite-size scaling in Fig. 2. That said, the time-crystal results are largely a demonstration that the existing framework transfers to the new model, not a new mechanism. The soft spot is the local equilibrium approximation (Eq. 16). The reader and stress-test both flag it, and I agree it is the weakest link but not a serious problem. The bias is weak (chi = E/L), so deviations from the equilibrium product measure are O(L^-1) and subleading. The numerical agreement in Fig. 1 validates the closure empirically. I would not demand a rigorous proof here; the paper is appropriately scoped as a proof-of-concept study. One minor point: the paper is self-citing heavily in the time-crystal section, but the cited results are externally published and the transport coefficient derivation is independent, so this does not concern me. This paper is for researchers in nonequilibrium statistical mechanics working on driven diffusive systems, fluctuation theory, or time crystals. It deserves a serious referee who can check the derivation in Section III carefully and assess whether the time-crystal application adds enough beyond the prior framework to stand on its own. I lean toward accept for review.","headline":"Solid derivation of bulk-driven nonlinear KMP hydrodynamics with parameter-free transport coefficients; time-crystal application is a proof of concept on prior framework.","tokens_in":17032,"tokens_out":597,"would_cite":true,"duration_ms":39452,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Biased energy collisions produce programmable time-crystal phases","keywords":[],"falsifier":"Measure the mobility σ(ρ₀) = ⟨j⟩_st / E in kinetic Monte Carlo simulations for various densities and nonlinearity exponents β. If the power-law scaling σ(ρ) ∝ ρ^(β+2) fails, or if the Einstein relation σ = 2ρ²D is violated, the local equilibrium approximation underlying the hydrodynamic derivation does not hold. Separately, if the packing order parameter |z₁| does not show a continuous onset above λ_c = 2π/ρ₀ with finite-size scaling consistent with a second-order transition, the time-crystal mechanism fails for this model class.","tokens_in":16432,"feed_emoji":"🔮","tokens_out":3236,"duration_ms":125923,"temperature":0.7,"pith_summary":"The authors generalize the Kipnis-Marchioro-Presutti (KMP) model of stochastic energy diffusion by replacing its uniform energy redistribution rule with one that is (a) biased by an external field via local detailed balance, and (b) controlled by an energy-dependent collision kernel. Starting from the exact master equation, they derive the hydrodynamic constitutive relation j = -D(ρ)∂_xρ + σ(ρ)E using a local equilibrium approximation, obtaining closed-form expressions for the diffusivity D(ρ) and mobility σ(ρ) as integrals over the collision kernel. For power-law kernels π(ν) ∝ ν^β, both coefficients reduce to simple power laws satisfying an Einstein relation σ = 2ρ²D. Kinetic Monte Carlo simulations confirm these predictions quantitatively. As a proof of concept, the authors then apply a packing field — a spatially structured feedback field that pushes energy lagging behind the instantaneous center of mass and restrains energy moving ahead — to show that this bulk driving mechanism can induce a continuous phase transition to a time-crystal phase. Above a critical coupling λ_c = 2πm/ρ₀, the homogeneous state becomes unstable and traveling energy condensates emerge, breaking continuous time-translation symmetry. Multiple packing fields can be superposed and modulated in time to create programmable spatiotemporal patterns. The central object is the bulk-driven nonlinear KMP model, and the central mechanism is the biased redistribution rule that converts an external field into a bulk drift while preserving local energy conservation.","feed_headline":"Biased energy collisions produce programmable time-crystal phases","feed_subtitle":"Nonlinear transport coefficients derived from microscopic rules, with packing fields inducing traveling energy condensates.","key_machinery":"Bulk-driven nonlinear KMP model: a 1D lattice where nearest-neighbor energy collisions conserve total pair energy but are biased by an external field χ through an asymmetric redistribution distribution f_χ(α|ν) derived from local detailed balance. The hydrodynamic limit yields a constitutive relation j = -D(ρ)∂_xρ + σ(ρ)E with D(ρ) and σ(ρ) given by integrals over the microscopic collision kernel π(ν). A packing field E_x^(m)[ρ] = |z_m|sin(φ_m - 2πmx), coupled to the mth Fourier mode of the energy density, provides nonlinear feedback that amplifies density packing and drives the time-crystal transition.","core_discovery":"The paper shows that biasing the stochastic energy redistribution in a KMP-type lattice model — while preserving local energy conservation — produces a bulk drift term in the hydrodynamic equation, with transport coefficients that are explicitly computable nonlinear functions of the local energy density. Furthermore, coupling this bulk drive to a packing field that selectively amplifies density fluctuations around the energy field's center of mass triggers a second-order phase transition to a time-crystal phase: a state with traveling energy condensates that break continuous time-translation symmetry. The critical coupling for this transition is λ_c^(m) = 2πm/ρ₀, independent of the nonlinear","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Biasing stochastic energy transport produces time-crystal phases","Packing fields trigger traveling energy condensates in nonlinear diffusion","Bulk-driven KMP model yields time-crystalline order with traveling condensates","Nonlinear energy diffusion with biased collisions produces time crystals","Microscopic bias rules yield macroscopic time-crystal transport"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The derivation of the hydrodynamic constitutive relation assumes that the system rapidly relaxes to a locally Gibbsian measure — a product of exponential distributions controlled by instantaneous local averages — which is standard for weakly interacting systems but is here assumed rather than rigorously proven for the strongly nonlinear, driven regime.","fun_headline_variants_meta":{"raw":{"variants":["Biasing stochastic energy transport produces time-crystal phases","Packing fields trigger traveling energy condensates in nonlinear diffusion","Bulk-driven KMP model yields time-crystalline order with traveling condensates","Nonlinear energy diffusion with biased collisions produces time crystals","Microscopic bias rules yield macroscopic time-crystal transport","Packing fields drive phase transition to traveling energy condensates","Traveling energy condensates emerge from biased nonlinear diffusion","Critical coupling for time-crystal transition in driven energy diffusion"]},"model":"glm-5.2","effort":"low","cost_usd":0.0,"raw_usage":{"total_tokens":1718,"prompt_tokens":495,"completion_tokens":1223,"prompt_tokens_details":null},"tokens_in":495,"tokens_out":1223,"duration_ms":17668,"temperature":1.0,"reasoning_tokens":1190,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T14:07:21.856471+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"Measure the mobility σ(ρ₀) = ⟨j⟩_st / E in kinetic Monte Carlo simulations for various densities and nonlinearity exponents β. If the power-law scaling σ(ρ) ∝ ρ^(β+2) fails, or if the Einstein relation σ = 2ρ²D is violated, the local equilibrium approximation underlying the hydrodynamic derivation does not hold. Separately, if the packing order parameter |z₁| does not show a continuous onset above λ_c = 2π/ρ₀ with finite-size scaling consistent with a second-order transition, the time-crystal mechanism fails for this model class.","supporting_citations":[],"review_version":1}