{"id":"01b6fc84-5e02-4218-92d1-10386868937e","arxiv_id":"1908.06488","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In finite-time driven Hubbard chains, crossing the precursor to the Mott insulator transition reverses the sign of the work-distribution skewness and makes entropy production exceed work fluctuations.","lead":"This paper simulates small chains of interacting electrons driven by a time-varying electric field and computes the probability distribution of quantum work. It finds that crossing a many-body phase transition changes the shape of this distribution and makes energy dissipation dominate over fluctuations, which matters for designing tiny quantum devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed reversal of inequality (5) compares closed-system simulations to an open-system bound from ref. 32; the bound's hypotheses are not shown to hold, so the headline conclusion is not supported by the data.","rationale":"The paper's central claim is that entropy production dominates work fluctuations after crossing the pM-QPT, even for slow driving, in contrast to the classical work fluctuation-dissipation relation and to recent predictions for slowly-driven open quantum systems. The data supporting this claim come from exact diagonalization of a closed Hubbard chain, with the entropy production defined by Eq. (4) as a relative entropy to the thermal state of the final Hamiltonian. Inequality (5) is explicitly presented as a result for open quantum systems from ref. 32. A closed unitary evolution does not automatically satisfy the assumptions of an open-system fluctuation-dissipation bound, so observing a ratio above 1 is not by itself evidence of a reversal of that bound. This is the most load-bearing gap because the abstract and conclusion lean on the contrast with open-system predictions. The reader's identified weakest assumption, the imported pM-QPT precursor identification, is a real secondary concern, but even if the precursor were firmly established, the inequality-comparison issue would remain. The proposed test, repeating the protocol under an explicit weak thermal coupling and computing both quantities under the hypotheses of ref. 32, would settle whether the effect is a genuine reversal or simply an out-of-domain comparison. Since the reader already assigned a CONDITIONAL verdict and my concern reinforces that condition rather than changing it, the verdict should remain unchanged.","tokens_in":9200,"tokens_out":7449,"duration_ms":86989,"concrete_test":"Run the same L=4,6,8 Hubbard protocol with the system weakly coupled to a thermal bath at inverse temperature beta, using a Lindblad master equation with local thermal dissipators (e.g., dephasing or amplitude damping) for the same U and tau ranges as in Fig. 4. Compute the entropy production from the full system-plus-bath state and the work variance from the two-point measurement distribution, and check whether 2<Sigma>/(beta^2 Var(W)) exceeds 1 in the slow-driving regime where the assumptions of ref. 32 are satisfied. If the ratio remains at or below 1 whenever the bound applies, the claimed reversal is an artifact of comparing closed-system relative entropy to an open-system bound rather than a genuine violation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that after the pM-QPT the ratio 2<Sigma>/(beta^2 <(W-<W>)^2>) exceeds 1, reversing inequality (5). However, the simulations are closed-system unitary dynamics: footnote 50 states that the driving time tau is much smaller than decoherence or relaxation times, and Eq. (4) defines <Sigma> as the relative entropy between the final state and the equilibrium state of Hf. Inequality (5), by the paper's own description, is a result for slow quantum processes in open systems, imported from ref. 32. The simulated Hubbard chain has no bath and no dissipative channel, so the conditions under which (5) was derived are not satisfied. Observing a ratio larger than 1 in Fig. 4 (top right) is therefore not a 'reversing' of (5); it is an out-of-domain comparison. The ratio itself is a legitimate numerical fact, but the conclusion that entropy production dominates work fluctuations 'in contrast to ... recent predictions for slowly-driven open quantum systems' requires either an open-system simulation under the hypotheses of ref. 32 or a proof that (5) also holds for closed unitary driving. Without one of these, the paper's central novelty is unsupported. The pM-QPT precursor identification is also imported from refs. 39-41, but even if the precursor is accepted, the applicability of inequality (5) remains the load-bearing gap.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the statistics of quantum work and related irreversibility measures in finite half-filled Hubbard chains (L=4, 6, 8) driven by a time-dependent linear potential over finite times, from sudden quench to nearly adiabatic regimes. Work distributions are computed with the two-point measurement protocol using exact diagonalization and QuTiP. The authors analyze the first three central moments, the relative-entropy entropy production ⟨Σ⟩=S(ρτ||ρτ_eq), and trace distances to the final equilibrium and adiabatic states. Their central observations are that the skewness of the work distribution changes sign across the precursor to the metal–Mott insulator transition (pM-QPT) and that in the Mott-insulating regime the ratio 2⟨Σ⟩/(β^2 ⟨(W−⟨W⟩)^2⟩) exceeds 1, which they interpret as entropy production dominating work fluctuations and as a reversal of a fluctuation–dissipation inequality recently proposed for slow quantum processes in open systems. They propose the skewness as a witness of irreversibility and discuss implications for small-scale quantum technologies.","tokens_in":9470,"tokens_out":8714,"duration_ms":89157,"significance":"If the central claims were fully established, the paper would provide a concrete, finite-size example in which strong correlations and a phase-transition precursor alter the fluctuation–dissipation balance in finite-time quantum thermodynamics, with potential relevance for the design of driven nanoscale devices. The numerical work is carefully executed and transparent: the moments, entropy production, and trace distances are computed independently from the same exact-diagonalization data, without parameter fitting, and the code is made available through a GitHub repository. However, the headline comparison to inequality (5) is made with closed-system unitary dynamics, whereas the inequality is imported from a theory of slow processes in open quantum systems. This out-of-domain comparison leaves the main conclusion unsupported. The paper is therefore potentially significant, but it requires substantial additional argument before its central claim can be accepted.","major_comments":[{"comment":"The central quantitative claim—that after the pM-QPT the ratio 2⟨Σ⟩/(β^2⟨(W−⟨W⟩)^2⟩) exceeds 1 and thereby reverses inequality (5)—is not supported by the simulations as presented. The simulations are closed-system unitary dynamics (footnote 50), while inequality (5) is stated to be a result for slow quantum processes in open systems, imported from ref. 32. The numerical ratio in the upper right panel of Fig. 4 is a valid fact for the unitary protocol, but because the hypotheses under which (5) was derived are not satisfied, the comparison is out of domain. The conclusion that entropy production dominates work fluctuations “in contrast to recent predictions for slowly-driven open quantum systems” therefore requires either an open-system simulation that meets the assumptions of ref. 32 or a proof that inequality (5) holds for closed unitary driving of the kind studied here.","section":"Entropy production and work fluctuation-dissipation relation; Eq. (5), Fig. 4"},{"comment":"The interpretation of the skewness sign change and of the fluctuation–dissipation reversal as signatures of crossing the precursor to the metal–Mott insulator transition rests on the assumption that the finite Hubbard chains studied here do undergo such a precursor, as suggested by refs. 39–41. The paper does not provide a direct characterization of the transition: no order parameter, no energy-gap or fidelity-susceptibility calculation, and no finite-size scaling collapse. The extrema of the computed quantities in Fig. 4 shift with L, but this alone does not establish a critical-point estimate or the existence of a quasi-QPT. Without such evidence, the statements that these are signatures of crossing a quantum phase transition lose their quantitative meaning.","section":"Introduction, paragraph beginning 'Finite Hubbard chains may undergo a precursor...'; Fig. 4"}],"minor_comments":[{"comment":"The abstract uses the phrase “quasi quantum phase transition” while the body consistently uses “precursor to the metal–Mott insulator transition (pM-QPT)”; please unify the terminology and define pM-QPT at first use.","section":"Abstract and Introduction"},{"comment":"The color scale in the skewness heatmaps is asymmetric and makes the white zero line difficult to locate for small τ; a symmetric diverging colormap would make the sign-change claim easier to verify.","section":"Fig. 3"},{"comment":"The two right-hand panels use different horizontal ranges (U/J up to 20 in the upper panel and up to 10 in the lower panel); aligning the ranges would make the comparison between the ratio and the trace distance to the adiabatic state more transparent.","section":"Fig. 4"},{"comment":"The supplemental material is hosted on a GitHub repository without a versioned DOI; for reproducibility, consider depositing a snapshot in a permanent archive.","section":"Footnote 45"},{"comment":"Ref. 32 is cited as arXiv:1905.07328v1; if the work has been published in the interim, the final journal reference should be used, and the hypotheses of the published version should be checked against the present closed-system setup.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper's numerical core is competently executed and the manuscript is transparent about its methods, but the headline comparison to inequality (5) is out of domain for closed-system unitary dynamics. The issues are fixable in revision—by adding an open-system calculation under the hypotheses of ref. 32 or by directly addressing whether inequality (5) extends to the closed setting—but as written the central claim overreaches. I would also ask the authors to verify the current published status of ref. 32 and to assess its hypotheses accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does a clean numerical job on work statistics in driven Hubbard chains and finds a real-looking sign change in the skewness across the interaction-driven precursor. The GitHub-hosted data and code are a plus. But the headline claim that the Mott phase 'reverses' inequality (5) is not supported: the simulations are closed unitary (footnote 50), and inequality (5) is imported from ref. 32 for slow open quantum processes. The ratio 2<Sigma>/(beta^2 var(W)) > 1 is a legitimate numerical fact, but calling it a reversal is an interpretive leap.\n\nWhat is actually new: applying the third moment of the work distribution to finite-time driven Hubbard chains, and showing that the skewness changes sign across the pM-QPT for tau > 0.5/J while staying small for sudden quenches. The scaling of trace distance and entropy production extrema with L is also suggestive, though not a finite-size collapse. The main observations are internally consistent and the paper is honest about its closed-system assumption.\n\nThe soft spots are two. First, the central quantitative claim about entropy production dominating work fluctuations 'in contrast to ... slowly-driven open quantum systems' requires either an open-system simulation satisfying the hypotheses of ref. 32, or a proof that inequality (5) holds for closed unitary driving. Without that, the comparison is out-of-domain. The paper could easily be fixed by rephrasing the claim as a closed-system feature that contrasts with the open-system prediction, rather than a reversal. Second, the pM-QPT identification is imported from refs. 39–41; the paper does not compute an order parameter or scaling collapse. That is a minor gap for a Letter, and the authors do call it a 'precursor,' so I would not make it the main objection.\n\nWho this is for: people working on quantum thermodynamics of small correlated systems, especially those interested in work distribution signatures. The paper deserves a serious referee, but the referee should insist on the inequality issue being addressed before publication. My inclination is to recommend major revision, not rejection.","headline":"Solid numerics and a genuinely interesting skewness signature, but the claimed 'reversal' of the open-system fluctuation-dissipation inequality (5) rests on comparing closed unitary dynamics to a bound that has not been shown to apply.","tokens_in":9998,"tokens_out":1853,"would_cite":false,"duration_ms":21197,"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":"Driving a Hubbard chain across the precursor to the metal–Mott-insulator transition flips the work distribution's skewness and makes entropy production dominate work fluctuations even at slow driving.","keywords":["quantum work distribution","Hubbard model","Mott insulator transition","entropy production","fluctuation-dissipation relation","skewness","finite-time driving","irreversibility"],"falsifier":"Perform an exact-diagonalization or tensor-network finite-size scaling study of the skewness extrema and of $2\\langle \\Sigma\\rangle/(\\beta^2\\langle (W-\\langle W\\rangle)^2\\rangle)$ for chains $L=4,6,8,10,12$: if the extrema do not move systematically toward $U=0$ as $L$ grows, or if the reversal of the fluctuation–dissipation inequality disappears once finite-size artifacts are removed, the attribution of these signatures to crossing a quantum phase transition is falsified.","tokens_in":8994,"feed_emoji":"⚛️","tokens_out":8300,"duration_ms":79068,"temperature":0.7,"pith_summary":"This paper aims to show that crossing the precursor to the metal–Mott-insulator quantum phase transition (pM-QPT) inside a finite-time driving protocol leaves a strong thermodynamic fingerprint in the quantum work distribution of a finite Hubbard chain. For half-filled chains driven by a time-dependent electric potential, from sudden quench to near-adiabatic ramps at finite temperature, the third moment (skewness) of the work distribution changes sign across the transition, a signature that survives as the chain grows from four to eight sites. The central quantitative claim is that in the Mott-insulating phase the irreversible entropy production dominates over work fluctuations, so the slow-driving relation $\\langle \\Sigma \\rangle \\leq \\beta^2 \\langle (W - \\langle W \\rangle)^2\\rangle/2$ is reversed even close to adiabaticity. A sympathetic reader would care because it says that correlations alone can make a driven small quantum system dissipate more than its work fluctuations indicate, which is a concrete constraint for designing quantum thermodynamic devices.","feed_headline":"Mott transition flips work asymmetry, reverses dissipation bound","feed_subtitle":"Driven Hubbard chains: entropy production beats work fluctuations even near adiabaticity, reversing the usual bound.","key_machinery":"The load-bearing object is the two-point-measurement quantum work distribution $P(W)=\\sum_{n,m} p_n^0 p_{m|n}^\\tau \\delta[W-(\\varepsilon_m^\\tau-\\varepsilon_n^0)]$ for the driven Hubbard Hamiltonian with hopping $J$, on-site repulsion $U$, and a linear potential that ramps over time $\\tau$. Its first three central moments carry the argument: the mean and variance give the work statistics, while the third moment (skewness) is shown to change sign across the pM-QPT and to track irreversibility. The quantitative reversal is carried by the dimensionless ratio $2\\langle \\Sigma\\rangle / (\\beta^2 \\overline{W}_2)$, where $\\langle \\Sigma\\rangle=S(\\rho_\\tau\\|\\rho_\\tau^{\\mathrm{eq}})$ is the relative-entropy entropy production and $\\overline{W}_2$ is the work variance; the trace distance $D_{\\mathrm{Tr}}(\\rho_\\tau,\\rho_\\tau^{\\mathrm{eq}})$ is used to show the final state remains far from equilibrium. All results come from exact diagonalization of finite chains, so the machinery is fully microscopic.","core_discovery":"The paper's central claim, stated in its conclusion, is that in the Mott-insulating phase entropy production dominates work fluctuations even for slow processes, in contrast with the classical work fluctuation–dissipation relation and with recent predictions for slowly driven open quantum systems. Concretely, for a half-filled Hubbard chain of length $L=4, 6, 8$ prepared in equilibrium at $\\beta^{-1}=2.5J$ and driven by a linear potential ramp over times $0.2/J \\leq \\tau \\leq 10/J$, exact diagonalization shows that the ratio $2\\langle \\Sigma \\rangle / (\\beta^2 \\langle (W-\\langle W\\rangle)^2\\rangle)$ exceeds one after the pM-QPT. The same calculations show the work distribution becomes nearly independent of driving rate in the insulating phase, its skewness changes sign across the transition for finite-time dynamics, and the trace distance from the equilibrium state peaks at the transition. These features are presented as signatures of the precursor to the metal–Mott-insulator quantum phase transition in finite chains.","pith_inferences":["A direct experimental translation would be to ramp cold atoms in an optical lattice across the Mott transition and measure the work distribution via projective number or energy measurements; the prediction is a skewness sign flip and a ratio $2\\langle \\Sigma\\rangle/(\\beta^2\\overline{W}_2)>1$ in the insulating phase.","If the reversal is generic for interaction-driven transitions, quantum thermal machines using strongly correlated media will carry an unavoidable dissipation that standard work-variance estimates miss; protocols could be designed to avoid crossing the transition.","The drift of the extrema toward $U=0$ as $L$ grows hints that, in the thermodynamic limit, the skewness flip and the inequality reversal coincide with the Mott transition point itself; larger-scale numerical methods could test this."],"forward_implications":["The skewness of the work distribution can act as a finite-time witness of the pM-QPT, showing a clear sign change even for chains short enough for exact diagonalization.","Sudden-quench protocols hide the dynamical part of the QPT signature; the effect only appears for finite-time driving with $\\tau J\\gtrsim 0.5$.","In the Mott phase, slow near-adiabatic driving does not restore the classical fluctuation–dissipation balance: dissipation remains larger than work fluctuations.","After crossing the transition, the driven system remains far from equilibrium, as measured by trace distance, even when the work distribution looks close to linear response.","For weak interactions, increasing the chain length makes $P(W)$ more classical (bell-shaped), while the Mott phase restores strong quantum features regardless of size."],"supporting_citations":[{"why":"defines the quantum work distribution used throughout as a two-point measurement over initial and final eigenstates.","marker":"9"},{"why":"provide the classical work fluctuation–dissipation relation $\\langle\\Sigma\\rangle=\\beta^2\\overline{W}_2/2$ that the paper's result contrasts with.","marker":"30,31"},{"why":"states the inequality $\\langle\\Sigma\\rangle\\leq\\beta^2\\overline{W}_2/2$ for slow open quantum systems that the paper finds reversed in the Mott phase.","marker":"32"},{"why":"establishes that finite half-filled Hubbard chains exhibit a precursor to the metal–Mott-insulator transition, the phenomenon the paper studies.","marker":"39–41"},{"why":"motivates the driven Hubbard model by reporting promising thermodynamics in a driven Hubbard dimer.","marker":"43"},{"why":"supplies the numerical time-propagation routines used for the exact-diagonalization results.","marker":"44"}],"fun_headline_variants":["Mott phase reverses entropy-work balance","Irreversibility beats work fluctuations at Mott transition","Mott transition flips work fluctuation-dissipation bound","Entropy production dominates work fluctuations in Mott phase"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's central attribution rests on treating the finite chains as undergoing a genuine precursor to the metal–Mott-insulator quantum phase transition at the interaction strengths studied, based on earlier work on short Hubbard chains; the paper itself computes no order parameter, critical exponent, or finite-size collapse, so if that precursor is only a generic finite-size crossover, the claimed QPT signatures would lose their special meaning.","fun_headline_variants_meta":{"raw":{"variants":["Mott phase reverses entropy-work balance","Irreversibility beats work fluctuations at Mott transition","Mott transition flips work fluctuation-dissipation bound","Entropy production dominates work fluctuations in Mott phase"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001396,"raw_usage":{"total_tokens":5662,"prompt_tokens":979,"completion_tokens":4683,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":4623}},"tokens_in":595,"tokens_out":4683,"duration_ms":33684,"temperature":1.0,"reasoning_tokens":4623,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:44:00.913700+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform an exact-diagonalization or tensor-network finite-size scaling study of the skewness extrema and of $2\\langle \\Sigma\\rangle/(\\beta^2\\langle (W-\\langle W\\rangle)^2\\rangle)$ for chains $L=4,6,8,10,12$: if the extrema do not move systematically toward $U=0$ as $L$ grows, or if the reversal of the fluctuation–dissipation inequality disappears once finite-size artifacts are removed, the attribution of these signatures to crossing a quantum phase transition is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"defines the quantum work distribution used throughout as a two-point measurement over initial and final eigenstates."}],"review_version":1}