{"id":"8fa0d67a-2a26-467a-a20b-6d1c809871f6","arxiv_id":"2508.21768","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A four-stroke anyon-Otto cycle converts anyon exclusion energy into work, and weakly interacting anyons maximize low-temperature work at intermediate statistical angles.","lead":"This paper proposes a four-stroke engine cycle that changes the quantum statistics of particles between boson-like and fermion-like states, extracting work from the energy tied up in anyonic exclusion. It reports that weakly interacting anyons on a lattice produce the most work at intermediate statistics, not at the boson or fermion extremes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Recompute Fig. 4 with Eq. (7)'s W': the intermediate-θ work peak may vanish once the θ-ramp is charged to the work reservoir, and the non-interacting reconciliation does not cover U2>0.","rationale":"The paper's central new result is that in the weakly interacting limit, work output peaks at an intermediate statistical angle rather than at the bosonic or pseudo-fermionic extremes. The authors explicitly define W so that θ-ramping during the thermalization strokes is treated as heat, and they later introduce W' in Eq. (7) to account for the work cost of those ramps. For the interacting case, this correction is not a minor bookkeeping detail: ΔB is computed with the state ρ2 that also determines the expansion-stroke work, and H(θ, U) depends on θ through the density-dependent hopping phase whenever U > 0. Therefore the θ-dependence of W' can differ substantially from that of W. The non-interacting limit cannot validate the interacting curves because H is θ-independent at U = 0, making ΔB vanish identically. This is the most load-bearing assumption because the abstract and Fig. 4 make their advantage claim specifically using W, while the paper's own consistency section suggests W' is the physically appropriate work measure. I am not claiming the result is wrong; I am claiming it is unverified under the more conservative work accounting, and a single exact-diagonalization recomputation would settle it. Other concerns, such as the unverified adiabatic assumption for interacting strokes and the missing experimental protocol, are secondary and addressable, but the work-definition issue directly targets the headline claim.","tokens_in":13891,"tokens_out":4567,"duration_ms":51871,"concrete_test":"Exact-diagonalize the L=8, N=4 anyon-Hubbard model with U1 = 0, U2 = 0.1, J = 1.0, θ2 = 0, TA = TB = 0.1. Reproduce Fig. 4(a) as W(θ1), and also compute W'(θ1) = W − ΔA − ΔB using Eq. (7), with ΔB = Tr[ρ2(H(0, U2) − H(θ1, U2))] and ρ2 the post-expansion state in the adiabatic limit (diagonal in the H(θ1, U2) eigenbasis with initial thermal populations). Check whether W' has an interior maximum at 0 < θ* < π. Repeat for U2 = 0.05 and 0.2 to see whether any non-monotonicity persists. If W' is monotone or peaks at θ = 0 or π, the abstract's advantage claim is not supported under consistent work accounting.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim in the abstract and Fig. 4 uses W = W12 + W21, in which the energy change during the θ-ramp sub-strokes is booked as heat (Eq. 4). The paper's own Eq. (7) defines W' = W − ΔA − ΔB as the work when the θ-ramp is charged to the work reservoir. For the interacting protocol (U1 = 0, θ2 = 0), ΔA = 0, but ΔB = Tr[ρ2(H(0, U2) − H(θ1, U2))] is nonzero for U2 > 0 because the density-dependent hopping term in Eq. (1) is θ-dependent. Since ρ2 is exactly the state whose energy difference −ΔE_A(θ1) produces the intermediate-θ peak, ΔB may compensate or shift the θ-dependence of W'. The non-interacting U = 0 curves are unaffected because H is θ-independent there, so the reconciliation in Fig. 6 does not test the interacting regime. Thus the claimed 'greater work extraction' from interacting anyons may be an artifact of not charging the statistical stroke to the work reservoir. This is a modeling choice, but it is load-bearing for the abstract's main quantitative claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a four-stroke hybrid anyon-Otto (HAO) cycle for the 1D anyon Hubbard model, in which the statistical angle θ is changed during the two thermalization strokes and the energy change associated with that change is booked as heat (Eqs. 4 and 6). For U=0 and low bath temperatures, exact diagonalization shows that the work output W increases monotonically with θ, approaching a maximum in the pseudo-fermionic limit. For weak interactions U2>0 with U1=0, the same exact-diagonalization approach reports that W/NU2 is maximized at an intermediate statistical angle for N ≥ L/2. The authors explain this via a first-order perturbation argument: the interaction-energy cost is smallest at intermediate θ because double-occupancy probability is suppressed there. They also reconcile the apparent second-law violation of the low-temperature inverse-accelerator mode by redefining work as W' = W − ΔA − ΔB (Eq. 7), and report enhanced maximum efficiency for high θ in the non-interacting case.","tokens_in":14272,"tokens_out":19535,"duration_ms":204437,"significance":"If the central claim holds, the paper demonstrates a nontrivial synergy between anyonic statistics and interactions in a quantum thermal machine, which would be a genuinely new thermodynamic effect. The numerical work is transparent and parameter-free: the exact-diagonalization results are direct solutions of the stated model, and the first-order perturbation explanation is supported by the collapse of the data in Fig. 4(b). The paper is also honest about the ambiguous role of the θ-ramp work, explicitly discussing the second-law issue and proposing Eq. (7). The main significance is conditional on resolving whether the intermediate-θ advantage survives under a consistent accounting of the work cost of changing θ.","major_comments":[{"comment":"The central claim that W is maximized at intermediate θ uses W = W12 + W21, in which the energy change during the θ-ramp sub-strokes is booked as heat (Eq. 4). The paper's own reconciliation defines W' = W − ΔA − ΔB (Eq. 7) as the work when the θ-ramps are charged to the work reservoir. For the interacting protocol (U1=0, θ2=0, U2>0), ΔB = Tr[ρ2(H(0,U2) − H(θ1,U2))] is nonzero and θ-dependent because the density-dependent hopping in Eq. (1) depends on θ. Moreover, ΔB is evaluated in exactly the state whose energy difference produces the intermediate peak. The non-interacting reconciliation in Fig. 6 does not cover U2>0. The authors should recompute Fig. 4(a) with W' and report the fate of the intermediate-θ peak; without this, the abstract's 'greater quantum thermodynamic advantage' claim is not established.","section":"§Weak interaction U ≪ J, Fig. 4(a), Eq. (7)"},{"comment":"The interacting work strokes are assumed to be perfectly adiabatic (τ→∞), with the assertion that non-adiabatic excitations 'do not change results qualitatively'. For U1=0, the initial Hamiltonian is non-interacting and gapless at low temperature for the system sizes used, so the adiabatic limit is not guaranteed to be well-defined. Without finite-τ data or a gap analysis along the ramp, the ED curves in Fig. 4 may not represent the claimed quasi-static limit. Since Fig. 4 is the central numerical result, this assumption requires support, e.g., finite-τ calculations or an estimate of the many-body gap along the U-ramp.","section":"§Weak interaction U ≪ J, τ→∞ assumption"}],"minor_comments":[{"comment":"The abstract promises 'we also outline an experimental protocol to realize the HAO cycle', but no such protocol section appears in the text; the final discussion only cites existing experimental realizations [47,57,58].","section":"Abstract vs. main text"},{"comment":"The notation W′12(21) = W12(21) − ΔB(A) is introduced only in the reconciliation section. For clarity, the authors should apply it explicitly to the interacting protocol and state whether Fig. 4 changes.","section":"Eq. (7)"},{"comment":"Typo: 'refgrigerator' should be 'refrigerator'.","section":"Introduction"},{"comment":"The y-axis label 'max' is incomplete; it should be η_max or equivalent.","section":"Fig. 6(d)"},{"comment":"The term 'anyon energy' is used as an analogue of the Pauli energy but is never defined. A formula, e.g., the ground-state energy difference between different θ values at fixed particle number, would make the resource explicit.","section":"General"},{"comment":"The y-axis label W/NU2 and the text 'scaled with U2' may be misread as implying W ∝ U2^2. Please clarify the scaling convention.","section":"Fig. 4(a)"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the work-accounting issue in Fig. 4. In my reading, the companion stress-test's premise that H is θ-independent at U=0 is inconsistent with the paper's own Fig. 3, which shows θ1-dependent W at U=0; however, the broader concern about Eq. (7) is valid and is the basis for my recommendation. The authors should be required to show W' results for U2>0 and to address the adiabaticity assumption. I would not reject on this basis, as the numerical setup is sound and the requested analysis is straightforward."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The paper proposes a four-stroke anyon-Otto cycle on the 1D anyon Hubbard model and reports two results: in the noninteracting limit, low-T work grows monotonically with theta toward the pseudo-fermionic point; with weak interactions and N>=L/2, work peaks at intermediate theta. The cycle construction itself is clean and the numerics (exact diagonalization, no fitted parameters) are internally consistent. The perturbative explanation for the intermediate peak—double-occupancy probability minimized at intermediate theta—is plausible and supported by the density distributions in Fig. 5. That part is genuinely new relative to the earlier anyon-engine papers, which stay noninteracting or few-body.\n\nThe soft spot is the work bookkeeping. The central figure uses W = W12 + W21, where the energy change during the theta-ramp sub-strokes is booked as heat. The paper itself acknowledges this in Eq. (7) and reconciles the noninteracting case with W'. But the reconciliation is only done for U=0, where H is theta-independent, so the correction vanishes there. For the interacting protocol (U1=0, theta2=0), the stroke-2 ramp has Delta_B nonzero because the density-dependent hopping term depends on theta. Since rho2 is the very state whose energy difference drives the intermediate-theta peak, Delta_B may compensate or shift the theta-dependence. The paper does not show Fig. 4 replotted with W'. This is not an error in the numerics—it's a modeling choice about what counts as work—but the abstract's claim that interacting anyons offer \"greater quantum thermodynamic advantage\" depends on that choice. If the theta-ramp is charged to the work reservoir, the advantage may shrink or vanish. The stress-test note is on point. I'd ask for a recomputation.\n\nTwo smaller issues: the abstract says an experimental protocol is outlined, but the text just points to existing anyon Hubbard realizations; there's no actual protocol. That should be either supplied or softened. And the adiabatic assumption for the interacting work strokes is stated without evidence; finite-tau effects are dismissed with a hand-wave. Minor, but worth a sentence.\n\nBottom line: the paper is a solid, clearly written proposal with an interesting new cycle and a plausible mechanism, but the main quantitative claim is currently propped up by a work definition that excludes a cost the authors themselves have identified. It deserves peer review, and the fix is straightforward: compute W' for the interacting case and see if the intermediate-theta peak survives. If it does, the paper is a nice result. If not, the introductory claim needs real revision. I'd send it to a referee with that instruction, and would bring it to a reading group for the thermodynamics discussion.","headline":"Interesting anyon-Otto cycle with clean numerics, but the main 'advantage' figure excludes the cost of changing statistics; that needs fixing before the quantitative claim stands.","tokens_in":14708,"tokens_out":2744,"would_cite":false,"duration_ms":31137,"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":"A four-stroke engine on the 1D anyon Hubbard model converts exclusion statistics into work; with weak interactions and at least half filling, work peaks at an intermediate statistical angle, not the bosonic or pseudo-fermionic endpoints.","keywords":["anyon Hubbard model","quantum Otto cycle","quantum thermal machine","anyonic statistics","low-temperature work extraction","Pauli engine","optical lattice anyons"],"falsifier":"Compute the cycle work W' = W − ΔA − ΔB (Eq. 7) for the interacting case U≪J with N≥L/2, and check whether W' still has an interior maximum in θ. If the interior peak disappears, the claimed anyonic advantage is an artifact of the θ-ramp being classified as heat. Alternatively, measure the on-site double-occupancy probability P(n=2) in an optical-lattice realization of the anyon Hubbard model as a function of θ: the paper's mechanism predicts a dip at intermediate θ, and its absence would falsify the explanation.","tokens_in":13851,"feed_emoji":"⚛️","tokens_out":7048,"duration_ms":71783,"temperature":0.7,"pith_summary":"The paper proposes a four-stroke hybrid anyon-Otto (HAO) cycle for the one-dimensional anyon Hubbard model, in which the statistical parameter θ is tuned during the thermalization strokes, analogous to the Pauli engine's bosonization/fermionization strokes. Its central claim is that in the weakly interacting limit (U ≪ J) at low temperature and with at least half filling, the work output per particle is maximized at an intermediate statistical angle 0 < θ* < π, rather than at the bosonic (θ=0) or pseudo-fermionic (θ=π) extremes. The mechanism is that the ground-state energy cost of switching on the interaction is smallest at intermediate θ, because the density-dependent hopping phase suppresses double occupancy most effectively there. In the non-interacting limit, work grows monotonically with θ and peaks at the pseudo-fermionic endpoint; the paper also reports an \"inverse accelerator\" operating mode whose apparent second-law violation disappears when the work of changing θ is properly accounted for. If the claims hold, interacting anyons would offer a genuine low-temperature thermodynamic advantage over both bosons and pseudo-fermions, and the existing experimental realization of the anyon Hubbard model would make the predictions testable.","feed_headline":"Interacting anyons beat bosons and fermions in a quantum heat engine","feed_subtitle":"A four-stroke anyon engine gets the most low-temperature work from particles at an intermediate statistical angle, not the extremes","key_machinery":"The central object is the anyon Hubbard model, a bosonic lattice model with a density-dependent hopping phase θ that interpolates between bosonic (θ=0) and pseudo-fermionic (θ=π) statistics. The HAO cycle uses two unitary strokes (ramping λ=J or U) and two thermalization strokes during which θ is changed in contact with heat baths. The load-bearing quantity is the anyon energy—the excess low-temperature energy produced by exclusion statistics, analogous to Pauli energy—together with the first-order perturbation correction Eper = ⟨ψ0|Hint|ψ0⟩ to the non-interacting ground-state energy; this correction's θ-dependence explains the intermediate-θ peak. The reconciliation formula ΔB(A)=Tr[ρ2(1)(H","core_discovery":"On the paper's own terms, the core discovery is that interactions and anyonic statistics cooperate to improve the HAO cycle's low-temperature work output. For U=0, W/N increases monotonically with θ1 and is largest in the pseudo-fermionic limit θ1→π, because the anyon energy built up during the anyonization stroke grows with θ. When a weak interaction U2 is switched on during the expansion stroke (with U1=0, J fixed), the change in ground-state energy ΔEG is minimized at intermediate θ for N≥L/2; since W is defined with a negative sign relative to this energy change, W peaks at θ*=θ*(U2) strictly between 0 and π. The paper traces this to a lower probability of double or higher occupancy at i","pith_inferences":["If the θ-ramp work is charged to the work reservoir throughout the interacting regime, the intermediate-θ peak in W may be reduced or shifted; whether a net advantage survives exact work bookkeeping is a question the paper leaves open for the interacting case.","The same double-occupancy suppression mechanism should also shape other figures of merit—efficiency at finite power, coefficient of performance as a refrigerator, and entropy production—so similar intermediate-θ extrema may appear there.","The prediction could be tested directly by measuring the on-site occupation probability P(nj=2) as a function of θ in a realized 1D anyon gas; the dip at intermediate θ is the microscopic signature behind the work peak.","The N≥L/2 threshold suggests the advantage is a filling-dependent correlation effect; at lower fillings, the monotonic trend toward pseudo-fermionic statistics should reappear."],"forward_implications":["At low temperature, the HAO engine can produce finite work even when both baths are cold, because the anyon energy acts as a quantum-statistical fuel.","In the weakly interacting regime at N≥L/2, tuning θ to an intermediate value gives more work per particle than either the bosonic or pseudo-fermionic endpoint.","The apparent inverse accelerator mode is an artifact of classifying θ-ramp energy as heat; with the corrected accounting, the second law is restored and an engine regime with enhanced maximum efficiency appears for θ1≥2.2.","Since the 1D anyon Hubbard model has been realized in optical lattices, the predicted θ-dependence of work and of double-occupancy probabilities is experimentally accessible.","The filling threshold N≥L/2 indicates the effect is tied to interaction-induced density correlations rather than single-particle band structure."],"supporting_citations":[{"why":"Defines the four-stroke quantum Otto cycle and the work/heat accounting convention the HAO cycle modifies.","marker":"[2]"},{"why":"The Pauli engine, whose bosonization/fermionization strokes motivate the HAO cycle and the concept of statistical energy as a fuel.","marker":"[23]"},{"why":"Introduces the 1D anyon Hubbard model and the density-dependent hopping phase that carries the statistical parameter θ.","marker":"[47]"},{"why":"Derives ground-state properties of the anyon Hubbard model that underlie the numerical work-output calculations.","marker":"[49]"},{"why":"Prior analysis of thermodynamics of statistical anyons that this work extends to a many-body lattice model.","marker":"[43]"},{"why":"A prior anyonic Otto engine model, providing a direct baseline for the HAO cycle.","marker":"[44]"},{"why":"Experimental realization of one-dimensional anyons with tunable statistical phase, making the proposed cycle experimentally testable.","marker":"[57]"},{"why":"Observation of anyonization of bosons in a quantum gas, supporting the feasibility of the Pauli-style statistical strokes.","marker":"[58]"}],"fun_headline_variants":["Anyon heat engine: interactions shift peak work to mid-statistics","With weak interactions, anyons at intermediate angles maximize engine work","Anyon-Otto cycle: interactions move work maximum off the fermionic limit","Interacting anyons outdo bosons and pseudo-fermions in low-temperature work"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The central work-output comparison treats the energy cost of changing the statistical parameter θ during the two thermalization strokes as heat, not work; the paper itself redefines work in Eq. (7) to include that cost when reconciling with the second law, and if the same correction is applied to the interacting regime the interior peak may shrink or vanish.","fun_headline_variants_meta":{"raw":{"variants":["Anyon heat engine: interactions shift peak work to mid-statistics","With weak interactions, anyons at intermediate angles maximize engine work","Anyon-Otto cycle: interactions move work maximum off the fermionic limit","Interacting anyons outdo bosons and pseudo-fermions in low-temperature work"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00057,"raw_usage":{"total_tokens":2520,"prompt_tokens":716,"completion_tokens":1804,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":460,"completion_tokens_details":{"reasoning_tokens":1723}},"tokens_in":460,"tokens_out":1804,"duration_ms":14847,"temperature":1.0,"reasoning_tokens":1723,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T13:58:38.656748+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the cycle work W' = W − ΔA − ΔB (Eq. 7) for the interacting case U≪J with N≥L/2, and check whether W' still has an interior maximum in θ. If the interior peak disappears, the claimed anyonic advantage is an artifact of the θ-ramp being classified as heat. Alternatively, measure the on-site double-occupancy probability P(n=2) in an optical-lattice realization of the anyon Hubbard model as a function of θ: the paper's mechanism predicts a dip at intermediate θ, and its absence would falsify the explanation.","supporting_citations":[{"cited_title":"Note that the two sub-strokes can be carried out simultaneously provided that the phase parameter is changed over a time interval much shorter than the thermalization time-scale","cited_arxiv_id":null,"evidence_quote":"Defines the four-stroke quantum Otto cycle and the work/heat accounting convention the HAO cycle modifies."},{"cited_title":"Mukherjee and U","cited_arxiv_id":null,"evidence_quote":"The Pauli engine, whose bosonization/fermionization strokes motivate the HAO cycle and the concept of statistical energy as a fuel."},{"cited_title":"The Anyonic Quantum Carnot Engine","cited_arxiv_id":"2504.20596","evidence_quote":"Derives ground-state properties of the anyon Hubbard model that underlie the numerical work-output calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior analysis of thermodynamics of statistical anyons that this work extends to a many-body lattice model."},{"cited_title":"Baldelli, B","cited_arxiv_id":null,"evidence_quote":"A prior anyonic Otto engine model, providing a direct baseline for the HAO cycle."},{"cited_title":"Zhang, D.-W","cited_arxiv_id":null,"evidence_quote":"Experimental realization of one-dimensional anyons with tunable statistical phase, making the proposed cycle experimentally testable."},{"cited_title":"Bartolomei, M","cited_arxiv_id":null,"evidence_quote":"Observation of anyonization of bosons in a quantum gas, supporting the feasibility of the Pauli-style statistical strokes."}],"review_version":1}