{"id":"2c388391-6127-4244-9d5c-9cc033bc50fe","arxiv_id":"2501.01309","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A Wannier-Stark field in the charger increases maximum power of Hubbard-model quantum batteries above a threshold, and can make bosonic batteries beat fermionic ones.","lead":"Researchers show that adding a tilted electric field, known as a Wannier-Stark field, to the charging step can boost the power of small quantum batteries made from ultracold atoms in optical lattices, especially for bosonic atoms. The work suggests a practical knob for improving energy storage in atomic quantum batteries, though the gains are not universal.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Environment-assisted steady-state ergotropy is the load-bearing risk: a uniform-temperature Davies bath should drive the chain to a passive Gibbs state, so the positive steady-state ergotropy in Fig. 6 needs a direct thermalization check.","rationale":"The reader's weakest assumption identifies exactly the same point: the environment-assisted ergotropy claim rests on the validity of the edge-local Davies master equation and on the steady state being non-thermal/non-passive. I agree that this is the most load-bearing concern, because it is the distinctive open-system result advertised in the title and abstract. The closed-system analytical results (Propositions 1-3 and Eq. (7)) appear internally consistent, and the numerical WS-activation claim, while under-documented without code, is not obviously self-contradictory. By contrast, a standard single-temperature Davies bath should thermalize a finite Hubbard chain to the passive Gibbs state, making positive steady-state ergotropy surprising. The paper does not provide a baseline check against the exact thermal state, does not report relaxation timescales, and does not identify any conserved quantity that would block thermalization. The Eq. (17) max/min error strengthens the need for an independent reimplementation, but the thermalization check is the decisive test. If the long-time state turns out to be thermal, the central environment-boosted claim fails; if a non-thermal steady state with ergotropy persists, the authors need to explain the mechanism. Either way, the manuscript's current CONDITIONAL verdict remains appropriate: the concern must be resolved before the open-system claims can be accepted.","tokens_in":16905,"tokens_out":11998,"duration_ms":142215,"concrete_test":"Recompute the open-system evolution for the Fig. 6 fermionic parameters (N=4, n_up=n_down=2, JF=1, UF=0, rF=0, beta=1, eta=0.01) with two initial states: the ground state used in the paper and the exact Gibbs state rho_beta = exp(-beta H_F)/Z. Evolve both to t = 10^5 (or until the population dynamics converge to a fixed point). If the Gibbs state is not stationary under Eqs. (13)-(16), or if the ground-state initial condition relaxes to zero ergotropy at long times, the reported steady-state environment-assisted ergotropy is not a genuine steady-state effect. Also recompute ergotropy with the standard definition Tr(H rho) - min_U Tr(H U rho U^dagger) to exclude the Eq. (17) sign issue.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is the open-system claim in Sec. IV, specifically environment-assisted ergotropy in the fermionic steady state without charging. The dissipator (13)-(16) is presented as a Davies generator built from eigenoperators of the battery Hamiltonian H_x with KMS rates, and the bath operators are local number operators. For such a generator with a single temperature (Fig. 6 uses TE1 = TEL = 1), the Gibbs state rho_beta = exp(-beta H_x)/Z is stationary; for a finite, generic chain with edge density couplings and no driving, the expected steady state is passive, so its ergotropy should be zero. The paper instead reports a positive steady-state ergotropy for fermions, and states that ergotropy increases when both baths have the same temperature. This is not a minor numerical issue: it indicates either (i) the implemented master equation is not the Davies generator described by Eqs. (14)-(16), (ii) the evolution time is too short and the 'steady state' is a slowly decaying transient, or (iii) there is a conserved-charge or decoherence-free structure that prevents full thermalization, which the paper neither identifies nor justifies. Separately, Eq. (17) defines ergotropy with a max over unitaries; the correct definition requires a min. If implemented literally, computed ergotropy would be non-positive, so the positive values in Fig. 6 imply either a typo or an unstated alternative formula. The closed-system power results are less at risk, but the title's 'environment-boosted' claim depends directly on this open-system result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies quantum batteries modelled by Bose- and Fermi-Hubbard chains with Wannier-Stark (WS) fields, in both closed and open dynamics. In the closed case, it reports analytic results for two-site batteries (Proposition 1), a hopping-only charging result for arbitrary size (Proposition 3), and numerical evidence for a critical WS charging strength above which maximum average power increases and bosonic batteries can outperform fermionic ones. In the open case, it claims environment-assisted ergotropy: fermionic batteries reach nonzero steady-state extractable work when the edge sites are coupled to local thermal baths, even without any charging unitary, while bosons show this only transiently. The paper also studies scaling of power with lattice size and particle number, and discharging in the presence of baths.","tokens_in":17212,"tokens_out":23896,"duration_ms":245072,"significance":"If the closed-system results are correct, they provide a concrete and potentially testable mechanism for improving quantum-battery power by adding a WS field to the charging Hamiltonian, and they identify a reversal of the boson/fermion power ordering. The analytic propositions and the single-particle exact solution are useful and appear internally consistent. The open-system claim of environment-assisted ergotropy is more surprising and, if correct, would be a notable thermodynamic phenomenon. However, the environment-assisted result requires a direct reconciliation with Davies-bath thermalization, and the advertised arbitrary-size closed form is currently an ansatz with unspecified coefficients. The numerical scans are extensive, but the manuscript does not include reproducible code or machine-checked proofs.","major_comments":[{"comment":"The ergotropy definition in Eq. (17) uses a maximum over unitaries: E(t)=W_x(t)-max_U Tr[H_x U rho U^dagger]. The standard definition uses a minimum of Tr[H_x U rho U^dagger] (equivalently, a maximum of the extracted work Tr[H_x rho]-Tr[H_x U rho U^dagger]). As written, E(t) is bounded above by W_x(t)-Tr[H_x rho], which is non-positive for nonnegative initial energy, so the positive steady-state values in Figs. 6-8 cannot come from Eq. (17) as stated. Please correct the sign/order of the optimization and state the exact formula actually implemented in the numerics.","section":"Eq. (17), Figs. 6-8"},{"comment":"The positive steady-state ergotropy for fermions at equal bath temperatures needs a thermalization check. For the Davies generator defined by Eqs. (14)-(16) with KMS rates and a single temperature (TE1=TEL=1), the canonical Gibbs state of the battery Hamiltonian H_x is stationary; within fixed particle-number sectors and under number-conserving unitaries, that state has zero ergotropy. The reported saturation of ergotropy therefore indicates either (i) the evolution time is too short and the state is a transient, (ii) the stationary state is non-thermal due to a conserved quantity or decoherence-free subspace, or (iii) the implemented master equation differs from Eqs. (14)-(16). The authors should compare the long-time state with e^{-beta H_x}/Z in each fixed particle-number sector and, if they differ, identify the symmetry or structure responsible.","section":"Sec. IV, Eqs. (13)-(16), Fig. 6"},{"comment":"The abstract promises a closed-form expression for the stored work when the battery is in the ground state of the hopping-only Hubbard model and is charged by onsite interactions and the WS field, irrespective of lattice size. The only arbitrary-size statement in the text is Eq. (5), which is introduced as an ansatz suggested by numerical simulation, with coefficients alpha, beta, and gamma left unspecified. This is not a closed form. Either provide explicit expressions for alpha, beta, and gamma with a derivation, or revise the abstract and the text to state that the arbitrary-size result is a fitted numerical form.","section":"Sec. II, Eq. (5), and Abstract"},{"comment":"Proposition 3 claims that for batteries with only hopping, charged by the WS field, the normalized work is W_x(t)=1-cos(r_c^x t) for any lattice size and particle number. No proof is given. Since this proposition underlies the scaling analysis in Sec. III and the 'arbitrary particle number' aspect of the closed-form claim, a proof or at least a clear derivation is needed. Without it, the proposition is unsupported.","section":"Sec. II, Proposition 3"}],"minor_comments":[{"comment":"The caption of Fig. 1 labels the vertical axis as J_c^F/U_c^F and the text says Case 2 varies J_c^F and r_c^F; these conventions are inconsistent and should be clarified.","section":"Fig. 1 and Sec. II"},{"comment":"References [50] and [62] appear to be the same paper (Konar et al., Phys. Rev. A 106, 022618 (2022)); merge or distinguish them to avoid duplicate citations.","section":"References"},{"comment":"There are several typos and grammatical errors, including 'femionic', 'explicitely', and 'All of them is computed'; the manuscript should be carefully proofread.","section":"Throughout"},{"comment":"The factor 2 on the right-hand side of Eq. (16) is unexplained; if the operators A_i(omega) are defined with an implicit normalization, this should be stated.","section":"Eq. (16)"},{"comment":"The system-bath coupling is called a dephasing model, but the KMS rates in Eq. (15) allow energy exchange and thermalization; a less misleading term would be 'local number-coupling bath' or 'amplitude-damping-type dephasing bath'.","section":"Sec. IV and Eqs. (11)-(12)"},{"comment":"For the fermionic initial state in Proposition 2, the paper uses a mixture of only the two S_z=0 states of the degenerate ground-state manifold; this choice should be justified or explicitly stated as a particular preparation.","section":"Sec. II, Proposition 2"}],"recommendation":"major_revision","confidential_remarks":"The closed-system results appear sound and could be publishable on their own. The open-system environment-assisted ergotropy claim is the main risk: it is currently in tension with standard Davies-bath thermalization, and the ergotropy definition in Eq. (17) is written incorrectly. I would ask the editor to require the authors to perform a direct Gibbs-stationarity check and to correct the ergotropy formula before the paper is reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Closed system, this paper is a decent piece of work. The two-site analytical result (Prop. 1) and the hopping-only closed form (Prop. 3) derive cleanly from the stated Hamiltonians, and the numerical observation that a Wannier-Stark field in the charging Hamiltonian can reverse the boson/fermion power ordering is new and worth reporting. The Eq. (5) ansatz with unspecified α, β, γ is unsatisfying but understandable if this is a numerical study; it doesn't undercut the main closed-system claims.\n\nThe open-system section is where things go wrong. First, Eq. (17) defines ergotropy with a max over unitaries; the correct definition uses a min. As written, the quantity is non-positive, so the positive values in Fig. 6 can't be the same quantity. This is a fixable typo, but it needs to be flagged. More seriously, the 'environment-assisted ergotropy' claim is not physically credible in its current form. The dissipator (13)-(16) is a Davies generator with KMS rates at a single temperature (T_E1 = T_E2 = 1 in Fig. 6). For such a generator, the canonical Gibbs state is stationary, and within the fixed particle-number sector the projected Gibbs state is passive, i.e., has zero ergotropy. The paper reports a positive steady-state ergotropy without ever computing the overlap with or distance to the thermal state. The likely resolutions are that the simulated time (t up to 300 with η=10^-2, so a relaxation scale around 100) is not long enough to reach the true steady state, or that some conserved quantity/decoherence-free structure blocks full thermalization. Either way, the authors need to demonstrate convergence to that state or identify the mechanism explicitly.\n\nNet: the closed-system power results are a modest but honest contribution to the Hubbard-battery subfield. The environment-assisted ergotropy, which is the title's 'environment-boosted' part, is not supported as it stands. This paper deserves peer review, but with the strong expectation of major revision: fix Eq. (17), check the steady state against the Gibbs state (or explain the non-thermal steady state), and provide either longer-time convergence data or the simulation code. I would not cite the open-system claim as is.\n\nFor a reading group, it's a decent case study in what can go wrong with 'steady-state' claims in weak-coupling master equations. I'd say maybe bring it, mostly for the open-system discussion.","headline":"Solid closed-system numerics and a clean two-site result, but the environment-assisted ergotropy claim is not credible: ergotropy is defined with max instead of min, and the alleged steady state is never checked against the thermal passive state.","tokens_in":17739,"tokens_out":6669,"would_cite":false,"duration_ms":66984,"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":"Adding a Wannier-Stark static field to the charging step raises the maximum average power of Bose- and Fermi-Hubbard quantum batteries once the field exceeds a critical strength, and can make bosonic batteries beat fermionic ones.","keywords":["quantum battery","Wannier-Stark field","Bose-Hubbard model","Fermi-Hubbard model","ergotropy","maximum average power","open quantum systems","ultracold atoms"],"falsifier":"Compute the steady state of the Lindblad equation (13)-(16) for a small fermionic chain (say $N=2$ or $4$) by exact diagonalization of the Liouvillian and evaluate its ergotropy; if the steady state is passive (zero ergotropy), environment-assisted ergotropy is absent. A second check is to replace the edge-only number coupling with a full coupling to a thermal reservoir and see whether the asymptotic ergotropy vanishes toward the Gibbs value.","tokens_in":16695,"feed_emoji":"🔋","tokens_out":8194,"duration_ms":70731,"temperature":0.7,"pith_summary":"This paper studies quantum batteries built from ultracold atoms in an optical lattice, modeled by the Bose- and Fermi-Hubbard Hamiltonians, and asks whether a static Wannier-Stark (WS) field—a linear tilt of the site energies—can help rather than hurt. The paper's central claim is that when the WS field is put into the charging Hamiltonian, the maximum average power rises once the field exceeds a critical strength, and bosonic batteries, which without the field lag behind fermionic ones at moderate interactions, overtake them. It also derives a closed-form expression for the stored work in two-site and effectively infinite lattices, and reports a new open-system effect: fermionic batteries with edge sites coupled to thermal baths reach a steady state with nonzero extractable work (ergotropy) even with no charging, while bosonic batteries show only transient ergotropy unless the WS field is present. If the claims hold, static fields and boundary thermal baths become resources for improving energy storage and extraction in cold-atom quantum batteries.","feed_headline":"Wannier-Stark field unlocks more power in Hubbard quantum batteries","feed_subtitle":"Above a critical field strength, maximum average power rises and bosonic batteries overtake fermionic ones.","key_machinery":"The central object is the Wannier-Stark field, the linear site-energy tilt $-r\\sum_i i\\,\\hat n_i$ added to the Hubbard Hamiltonian, used in the charging step together with the onsite interaction $U^c_x$. The load-bearing identity is the periodic work formula for a two-site battery, $W_x(t)=1-\\cos(r^c_x t)\\cos(U^c_x t)$ (with the fermionic and bosonic expressions coinciding), whose fitted many-site generalization $W^N_x(t)=\\alpha+\\beta\\cos(r^c_x t)+\\gamma\\cos(r^c_x t)\\cos(U^c_x t)$ makes the threshold and activation visible; the maximum average power is $P^{\\max}_x=\\max_t W_x(t)/t$. In the open-system part, the machinery is a local dephasing GKSL master equation in which bosonic baths couple to the edge-site number operators with KMS (detailed-balance) transition rates, and extractable work is measured by ergotropy $\\mathcal{E}(t)=\\mathrm{Tr}[H_x(\\hat\\rho-\\hat\\rho_{\\rm passive})]$.","core_discovery":"On the paper's own terms, the discovery is that the Wannier-Stark potential is a constructive control knob for Hubbard-model quantum batteries. For a battery initialized in the ground state of the Bose- or Fermi-Hubbard model with hopping and onsite interactions, charging with a Hamiltonian that includes both an onsite interaction $U^c_x$ and a WS field $r^c_x$ yields maximum average power $P^{\\max}_x$ that first dips and then rises with $|r^c_x|$; above a threshold $|r^{\\prime c}_x|$ the power is higher than without the field. For two sites with a hopping-only ground state the stored work is exactly $W_x(t)=1-\\cos(r^c_x t)\\cos(U^c_x t)$ for both statistics, and for larger lattices the numerics support $W^N_x(t)=\\alpha+\\beta\\cos(r^c_x t)+\\gamma\\cos(r^c_x t)\\cos(U^c_x t)$. The field reverses the usual ordering $P^{\\max}_F\\ge P^{\\max}_B$: for sufficient $r^c_x/U^c_x$, $\\Delta P^{\\max}_{F-B}=P^{\\max}_F-P^{\\max}_B$ becomes negative, an effect the paper calls activation of power, and the effect survives finite-temperature initial states. With edge baths, the paper claims 'environment-assisted ergotropy': fermions reach a steady state with nonzero ergotropy without any unitary charger, bosons reach it only transiently when the WS field is present, and under active charging both species show nonmonotonic steady-state ergotropy versus $U^c_x$ and $r^c_x$.","pith_inferences":["A consequence the authors leave implicit: the positive steady-state ergotropy in the fermionic battery requires the edge-bath dephasing model's steady state to be non-thermal; checking whether a full thermal bath coupling preserves the effect would settle how physical it is.","The two-site work formula suggests a resonance test: if $r^c_x$ and $U^c_x$ are commensurate, the work is periodic with the least common multiple of their periods, so a single lattice experiment could map the predicted revival structure.","The activation threshold likely tracks the tilt strength at which the Wannier-Stark ladder localizes single-particle eigenstates; if so, the boson-over-fermion inversion should persist for weak harmonic confinement or weak disorder, which could be probed in existing optical-lattice setups.","One could extend the environment-assisted idea to discharging: if edge baths can store work without a charger, the same steady state should be dischargeable by reversing the bath temperature bias, which the paper does not address."],"forward_implications":["A cold-atom battery designer can use the WS charging field as a switch: below threshold it hurts, above threshold it helps, and the threshold can be located from the ratio $r^c_x/U^c_x$.","Bosonic batteries, which are usually the weaker choice at moderate onsite interactions, become the better choice once the charging WS field exceeds the critical strength, even at finite temperature.","The WS field changes the scaling behavior of maximum average power: without it power falls with lattice size, with it power grows sublinearly with $N$, so larger optical lattices become beneficial.","In the open setting, attaching hot edge baths can store extractable work in fermionic batteries with no charging step, and moderate WS and onsite charging strengths can raise steady-state ergotropy after a threshold.","The closed-form work formula gives explicit periodic charging times and resonances, so the maximum-work and maximum-power times can be separated and tuned."],"supporting_citations":[{"why":"defines ergotropy, the extractable-work measure central to the environment-assisted charging claims.","marker":"[1]"},{"why":"sets the quantum-battery framework and the maximum-average-power figure of merit used throughout.","marker":"[2]"},{"why":"provides the earlier ultracold-atom Hubbard battery result that bosonic power falls below fermionic, the baseline the WS-field activation reverses.","marker":"[62]"},{"why":"supplies the Wannier-Stark ladder spectrum of the tilted Hubbard model that the charging protocol exploits.","marker":"[76]"},{"why":"gives the GKSL master-equation formalism used for the open-system evolution with edge baths.","marker":"[77]"},{"why":"provides the open-quantum-systems background for the Lindblad dissipator and its rates.","marker":"[78]"}],"fun_headline_variants":["Wannier-Stark field boosts power in Hubbard quantum batteries above critical strength","Critical Wannier-Stark field flips bosonic-fermionic battery power order","Hubbard battery power rises with Wannier-Stark field past a threshold","Wannier-Stark field enhances quantum batteries, bosons benefit most","Environment heat baths enable work extraction from fermionic batteries"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The environment-assisted ergotropy claim rests on the assumption that the edge-coupled dephasing baths drive the chain to a steady state with positive ergotropy, rather than to the thermal Gibbs state at the bath temperature, which has zero ergotropy.","fun_headline_variants_meta":{"raw":{"variants":["Wannier-Stark field boosts power in Hubbard quantum batteries above critical strength","Critical Wannier-Stark field flips bosonic-fermionic battery power order","Hubbard battery power rises with Wannier-Stark field past a threshold","Wannier-Stark field enhances quantum batteries, bosons benefit most","Environment heat baths enable work extraction from fermionic batteries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000471,"raw_usage":{"total_tokens":2443,"prompt_tokens":1144,"completion_tokens":1299,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":760,"completion_tokens_details":{"reasoning_tokens":1201}},"tokens_in":760,"tokens_out":1299,"duration_ms":11530,"temperature":1.0,"reasoning_tokens":1201,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:31:34.585377+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the steady state of the Lindblad equation (13)-(16) for a small fermionic chain (say $N=2$ or $4$) by exact diagonalization of the Liouvillian and evaluate its ergotropy; if the steady state is passive (zero ergotropy), environment-assisted ergotropy is absent. A second check is to replace the edge-only number coupling with a full coupling to a thermal reservoir and see whether the asymptotic ergotropy vanishes toward the Gibbs value.","supporting_citations":[{"cited_title":"To achieve this goal, the initial state is pre- pared in the ground state of the battery Hamiltonian, HF = −JF P ⟨ij⟩ ˆc† iσˆcjσ + h.c","cited_arxiv_id":null,"evidence_quote":"defines ergotropy, the extractable-work measure central to the environment-assisted charging claims."},{"cited_title":"Therefore, it is crucial to exam- ine the charging Hamiltonian, which governs the unitary evo- lution of the system","cited_arxiv_id":null,"evidence_quote":"sets the quantum-battery framework and the maximum-average-power figure of merit used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the earlier ultracold-atom Hubbard battery result that bosonic power falls below fermionic, the baseline the WS-field activation reverses."},{"cited_title":"Carisch, A","cited_arxiv_id":null,"evidence_quote":"supplies the Wannier-Stark ladder spectrum of the tilted Hubbard model that the charging protocol exploits."},{"cited_title":"Scarlatella, A","cited_arxiv_id":null,"evidence_quote":"gives the GKSL master-equation formalism used for the open-system evolution with edge baths."}],"review_version":1}