{"id":"b05b5cda-5133-4765-9049-3c67033f0da5","arxiv_id":"2608.13308","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For a driven qubit coupled to a bosonic bath, the nonsecular Floquet-Redfield heat current is derived, and an instantaneous-eigenbasis master equation reproduces its low-frequency behavior, whereas the full secular approximation changes peak heights.","lead":"This paper derives and compares several quantum master equations for heat flow out of a small, periodically driven quantum system into a cold bath, using a driven qubit coupled to an ohmic-Drude bath. It shows that a commonly used secular approximation distorts the height of heat-current peaks, while a nonsecular Floquet-Redfield equation and an instantaneous-eigenbasis equation agree at low driving frequencies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Fig. 3 comparison is weakened by the paper's own admission that the instantaneous-eigenbasis rates come from the Floquet basis; at Ad = 0.4 Delta, delta is not small, so the agreement may be an artifact of that rate transfer.","rationale":"The reader's weakest-assumption analysis already identifies the same load-bearing issue, and the paper itself flags it as an inconsistency. The main Floquet-Redfield formulas, Eqs. (28) and (35), appear internally consistent and provide a useful nonsecular steady-state framework; I do not find a soundness problem there that would change the verdict. The comparison with the instantaneous-eigenbasis master equation is the least secure part of the central claim: the dissipative rates are imported from the Floquet basis, so at the parameters of Fig. 3 the apparent agreement cannot be attributed to the instantaneous-eigenbasis master equation alone. The authors are explicit about the procedure, which keeps the issue at the level of an overclaim rather than an error. A targeted numerical test with consistently computed instantaneous rates would settle whether the concern is material. If the test shows a delta-dependent discrepancy, the verdict should remain conditional with a narrowed claim; if the test shows agreement, the concern is resolved. Since the reader already assigned CONDITIONAL, I leave the verdict unchanged.","tokens_in":19623,"tokens_out":15101,"duration_ms":174120,"concrete_test":"Recompute the dashed-blue curves of Fig. 3 with the rates evaluated consistently in the instantaneous eigenbasis: keep the unitary part of Eq. (55) unchanged, but replace the Floquet-basis rates of Eq. (50) by the instantaneous-gap rates 2Re W(+-omega_q(t))/hbar^2 multiplied by |<psi_0(t)|sigma_z|psi_1(t)>|^2 = Delta^2/(Delta^2 + epsilon(t)^2), and propagate to steady state for A_d = 0.4 Delta over omega_d/Delta in [0.05, 1.5]. Compare with the full Floquet-Redfield curves of Fig. 3. If the deviation grows with delta and exceeds the plotted linewidth above omega_d/Delta ~ 0.3, the agreement in Fig. 3 depends on borrowing Floquet rates, and the claim that Eqs. (55)-(57) reproduce the Floquet-Redfield heat current at moderate drive must be downgraded.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section V derives Eq. (55) in the instantaneous eigenbasis |psi_0/1(t)>, but the relaxation rates entering it are taken from Eq. (50), which is a Floquet-basis expression. The paper explicitly states after Eq. (55) that the Floquet basis and the instantaneous eigenbasis coincide only in the adiabatic limit delta = omega_d A_d / Delta^2 << 1, and it calls the procedure an inconsistency. In Fig. 3 the comparison is made at A_d = 0.4 Delta; delta reaches 0.4 at resonance and is 0.08-0.2 across the fractional-peak region. The plotted agreement is therefore not a clean test of the instantaneous-eigenbasis master equation: the dissipative part already contains Floquet-basis information, so the calculation mostly checks that the nonadiabatic unitary corrections are small in the chosen regime, not that the instantaneous-eigenbasis master equation is valid at arbitrary drive strength or frequency. The sentence claiming Eq. (55) is valid in principle for any drive strength/frequency is unsupported by the paper's own caveat, and the central comparison in Fig. 3 does not establish the claimed reproduction of the Floquet-Redfield current at moderate delta.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops weak-coupling (Born, second-order) master-equation descriptions of steady-state heat transport in a periodically driven quantum system coupled to a bosonic bath. Its central object is a non-Markovian, nonsecular Floquet-Redfield equation and the corresponding heat-current formula, Eqs. (28) and (35), derived from the Nakajima-Zwanzig equation without secular or Markov approximations. The authors compare this benchmark with a fully secular Floquet master equation, Eqs. (37)-(39), finding that the secular approximation overestimates the low-frequency fractional peaks and underestimates the resonance peak. They then study a master equation in the instantaneous eigenbasis, Eqs. (55)-(57), together with an analytical weak-drive solution in Appendix A, and claim that the instantaneous-eigenbasis approach reproduces the Floquet-Redfield heat current at low drive frequency and amplitude. The paper also discusses the suppression of even-fractional resonances at zero static bias. The numerical results are for a driven qubit with an ohmic-Drude bath, and comparisons are shown in Figs. 1-4.","tokens_in":19814,"tokens_out":5465,"duration_ms":55155,"significance":"If the claims are established, the main value of the paper is a systematic derivation of a nonsecular Floquet-Redfield heat current with explicit formulas, and a clear diagnosis of secular-approximation errors: full secularization overestimates fractional peaks and roughly halves the resonance peak, while coherence contributions reduce the fractional-peak heights. The analytical weak-drive solution in Appendix A gives a compact explanation of the multiphoton peaks and of the odd-harmonic selection rule at zero bias. The paper is also commendably explicit about the rate-transfer inconsistency in the instantaneous-eigenbasis comparison. However, because of that inconsistency, the comparison in Fig. 3 does not support the broad claim that Eq. (55) is valid for any drive strength and frequency, and the numerical benchmark lacks a convergence study. The significance of the paper therefore depends on a substantial revision of these claims.","major_comments":[{"comment":"The sentence \"Eq. (55) is valid in principle for any drive strength/frequency and only assumes weak system-bath coupling\" is not supported by the manuscript's own analysis. The relaxation rates entering Eq. (55) are taken from Eq. (50), which is a Floquet-basis expression, and the paper itself calls this an inconsistency. At the parameters of Fig. 3 with A_d = 0.4 Delta, the adiabaticity parameter delta = omega_d A_d / Delta^2 reaches 0.4 at resonance and is in the range 0.08-0.2 across the fractional-peak region, so the plotted agreement does not test the instantaneous-eigenbasis master equation at moderate drive. It tests, rather, that the nonadiabatic unitary corrections are small in that regime and that the Floquet-basis rates are transferable there. Because all the compared approaches are solutions of the same Born approximation, the agreement is also a consistency check rather than an independent validation. Please either derive rates consistently in the instantaneous eigenbasis, or explicitly restrict the claim and the Fig. 3 comparison to the adiabatic regime delta << 1 and remove or qualify the \"any drive strength/frequency\" statement. This is load-bearing for the abstract and introductory claim that the instantaneous-eigenbasis equation reproduces the Floquet-Redfield heat current.","section":"Section IV C, Fig. 1 caption"},{"comment":"The numerical benchmark relies on truncation parameters K = 6, L = 20, and N_Mats = 200, but no convergence study is reported. Since the Floquet-Redfield result is the reference against which the secular and instantaneous-eigenbasis approximations are judged, the paper should include a convergence check for these cutoffs at representative frequencies, including the resonance peak and a fractional peak, and state the resulting numerical tolerance. Without this, the quantitative claims such as the factor-of-two underestimation at resonance and the specific peak heights in Figs. 1-3 are not fully established.","section":"Section IV C, Fig. 1 caption"},{"comment":"The claim that the analytical solution \"essentially reproduces\" the numerical solution in the whole frequency range is not quantified. The agreement is shown only graphically, and at A_d = 0.4 Delta the expansion parameter A_d/Delta is not small. Please provide a quantitative measure of the deviation (e.g., relative error as a function of omega_d for each drive amplitude) and state the range of parameters in which the analytical formula is intended to be quantitative rather than illustrative.","section":"Appendix A, Eq. (A4) and Fig. 4"}],"minor_comments":[{"comment":"The text contains typos such as \"time evolutio\" and \"and and\" in the discussion of the bath correlation function; these should be corrected.","section":"Section II"},{"comment":"The caption sentence \"while is the number of terms retained in the sum over the Matsubara frequencies\" is missing the symbol N_Mats and is grammatically incomplete. The phrase \"The truncation of the sums over l in [-L,L] and k' in [-K,K] constants L and K\" is also garbled.","section":"Fig. 1 caption"},{"comment":"There are typos such as \"setady-state\" and \"het current\" (should be \"steady-state\" and \"heat current\", respectively).","section":"Section IV A and IV C"},{"comment":"The caption contains \"master euqation\" (should be \"master equation\"). Also, the abbreviation \"ME\" should be defined on first use in the caption or text.","section":"Fig. 3 caption"},{"comment":"The notation \"full secular approximation\" is used for Eqs. (37)-(39), but the relation between this approximation and the \"secular approximation\" in Eq. (48) could be made more explicit; in particular, Eq. (36) sets coherences to zero before the Markov approximation, while Eq. (48) is a secular Markovian equation for transients. Clarifying this distinction would help the reader.","section":"Section IV B"}],"recommendation":"major_revision","confidential_remarks":"The paper's central Floquet-Redfield derivation appears sound, and the comparison with the full secular approximation is informative. The main issue is the overstatement in Section V: the instantaneous-eigenbasis master equation is not valid for arbitrary drive strength/frequency when its rates are computed in the Floquet basis, as the authors themselves note. I would ask the editor to require the authors to qualify this claim and to add a convergence study for the numerical truncations. The paper fits the journal's scope and, with those revisions, could be publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I'll be direct: the genuinely new pieces here are the explicit nonsecular Floquet-Redfield heat-current formula, Eq. (35), and the analytical weak-drive solution in Appendix A. If you model driven qubit heat transport, those are useful. The paper does a clean job deriving the second-order steady-state formalism and is honest about the one ad hoc step: Eq. (55), the instantaneous-eigenbasis master equation, uses relaxation rates computed in the Floquet basis, which the authors themselves call an inconsistency. That honesty is real and should be credited.\n\nThe main soft spot is Fig. 3. At Ad = 0.4 Delta, the adiabaticity parameter delta = omega_d A_d / Delta^2 is 0.4 at resonance and between 0.08 and 0.2 over the fractional-peak region. That is not the adiabatic limit, so the agreement between Eq. (55) and the Floquet-Redfield result is not as clean as it looks: the dissipative part already knows about the Floquet basis. The comparison mostly demonstrates that nonadiabatic corrections to the unitary dynamics are small in that regime, not that the instantaneous-eigenbasis ME is valid at moderate delta. The text does point out the high-frequency discrepancy and blames delta, but the sentence claiming Eq. (55) is valid in principle for any drive strength/frequency is not supported by the paper's own caveat. That should be fixed.\n\nTwo smaller issues: there is no convergence study for the truncations K=6, L=20, N_Mats=200, and no code or data are provided. And because all three methods are solutions of the same Born approximation, the agreement between them is a consistency check, not independent validation. That is acceptable for a comparison paper, but it deserves a sharper statement.\n\nWho is this for? People doing practical modeling of driven quantum thermal machines, especially anyone deciding between secular and nonsecular master equations. The paper does not overturn prior results, but it quantifies the secular approximation error and gives a benchmark that is worth having. I would send it to review. The revision requests should be modest: soften the 'any drive strength' claim, add a short convergence check, and state more plainly that Fig. 3 is an adiabatic-regime comparison.","headline":"A careful comparison paper with two genuinely new formulas; the instantaneous-eigenbasis benchmark is weaker than it looks because of an admitted rate-transfer inconsistency, but the nonsecular Floquet-Redfield derivation is solid and worth citing.","tokens_in":20441,"tokens_out":3515,"would_cite":true,"duration_ms":34456,"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":"For a periodically driven quantum system weakly coupled to a heat bath, the steady-state heat current is fixed by the nonsecular Floquet-Redfield equations (28) and (35); against this benchmark, the full secular approximation…","keywords":["Floquet-Redfield equation","periodically driven qubit","heat current","spin-boson model","secular approximation","instantaneous eigenbasis","multi-photon resonances","weak coupling"],"falsifier":"Compute the steady-state heat current with the instantaneous-eigenbasis master equation (55)-(57) using relaxation rates evaluated directly in the instantaneous eigenbasis rather than transferred from the Floquet basis, at parameters where $\\delta = \\omega_d A_d / \\Delta^2$ is not small, e.g. $A_d = 0.4\\Delta$ and $\\omega_d = 1.5\\Delta$, giving $\\delta = 0.6$; if the resulting curve departs from the Floquet-Redfield benchmark of Eqs. (28) and (35), the reported agreement is an artifact of the rate transfer the paper itself flags as inconsistent.","tokens_in":2095,"feed_emoji":"⚛️","tokens_out":5628,"duration_ms":121708,"temperature":0.7,"pith_summary":"This paper tries to establish which master-equation schemes correctly describe heat transport in a periodically driven quantum system weakly coupled to a thermal bath. Its benchmark is the Floquet-Redfield equation for the steady state, derived from the Nakajima-Zwanzig equation without Markovian or secular approximations. Against that benchmark, the paper shows that the commonly used full secular approximation distorts the driven-qubit heat current: it inflates the low-frequency multi-photon side peaks and cuts the main resonance peak to about half its true value. A simpler master equation in the instantaneous eigenbasis, solved numerically and analytically, matches the benchmark in the low-frequency, low-amplitude regime. Getting this right matters for designing quantum thermal machines and for interpreting heat measurements on driven superconducting qubits.","feed_headline":"Full secular master equation halves the driven qubit's resonant peak","feed_subtitle":"The simplified instantaneous-eigenbasis equation, by contrast, matches the benchmark at low drive frequency.","key_machinery":"The central object is the Floquet-Redfield tensor $R^{k,k'}_{nmn'm'}$, built from the one-sided Fourier transform $W(\\omega)$ of the bath correlation function and the Fourier components $Q^l_{nm}$ of the coupling operator in the Floquet basis. It couples populations, coherences, and different Fourier harmonics of the steady-state density matrix through Eq. (28); the heat current to the bath is then the weighted sum Eq. (35), with weights set by the transition frequencies $\\omega^l_{nm} = \\varepsilon_n - \\varepsilon_m + l\\omega_d$ and the bath power spectrum. The comparison uses two simplifications of this object: the full secular truncation that keeps only diagonal populations, Eqs. (37)-(39), and the instantaneous-eigenbasis master equation (55)-(57), whose rates are taken from the Floquet calculation but whose unitary evolution follows the instantaneous qubit levels. A second key piece is the weak-drive analytical solution of that equation, Eqs. (A1)-(A4), which reduces the dynamics to a Bloch-type equation with drive-induced harmonics and yields closed-form Lorentzian peaks at $\\omega_d = \\bar{\\omega}_q/l$.","core_discovery":"The paper's central claim is that, for a periodically driven quantum system weakly coupled to a heat bath, the steady-state heat current is given exactly, to second order in the coupling, by the nonsecular Floquet-Redfield equations (28) and (35), with no Markovian or secular approximation needed. Applied to a driven spin-boson qubit with $\\hat{Q}=\\hat{\\sigma}_z$ coupling to an ohmic-Drude bath, this benchmark predicts multi-photon peaks at drive frequencies $\\omega_d = \\Delta/l$ once the qubit frequency is renormalized by the drive, and it predicts that even-fraction peaks are suppressed at zero static bias. The full secular approximation, which discards steady-state coherences and higher Fourier harmonics, overestimates the heights of the low-frequency fractional peaks and underestimates the resonant peak at $\\omega_d = \\Delta$ by roughly a factor of two. The master equation written in the instantaneous eigenbasis, Eqs. (55)-(57), reproduces the Floquet-Redfield heat current at low drive frequency and amplitude, and its weak-drive analytical solution captures the fractional peaks and the resonance. The paper also acknowledges that the rates entering the instantaneous-eigenbasis equation are computed in the Floquet basis, which coincides with the instantaneous eigenbasis only in the adiabatic limit $\\delta = \\omega_d A_d / \\Delta^2 \\ll 1$.","pith_inferences":["Editorial inference: If the Floquet-Redfield benchmark is trusted, studies of driven thermal machines that rely on secular master equations should quantify the nonsecular correction, because the resonant operating point is precisely where the secular value is off by a factor of two.","Editorial inference: The paper's own admitted inconsistency suggests a sharp numerical experiment: compute the relaxation rates in the instantaneous eigenbasis directly and rerun Fig. 3; sustained agreement at moderate $\\delta$ would then be nontrivial, while breakdown would locate the rate transfer as the source of the match.","Editorial inference: The predicted factor-of-two difference at $\\omega_d = \\Delta$ is a measurable signature in superconducting-qubit heat-transport setups; a measured resonance peak near the secular value would indicate that the nonsecular corrections are less relevant than the model claims.","Editorial inference: The even-fraction suppression at $\\epsilon_0 = 0$ could serve as an in-situ symmetry-point diagnostic in experiments, since any nonzero bias immediately activates even-harmonic peaks."],"forward_implications":["The nonsecular Floquet-Redfield equations (28) and (35) provide a steady-state heat-current formula for any periodically driven system at weak coupling, free of Markovian and secular assumptions, so they can serve as a benchmark for all simpler schemes.","For the driven spin-boson qubit, the full secular approximation is quantitatively unreliable for heat currents: it overestimates the fractional side peaks and underestimates the resonance peak by about a factor of two.","The instantaneous-eigenbasis master equation, despite its rate-transfer inconsistency away from the adiabatic limit, reproduces the Floquet-Redfield heat current at low drive frequency and amplitude and captures the magnitude of the resonant maximum.","Heat-current peaks appear at drive frequencies matching fractions of the renormalized qubit frequency, $\\omega_d \\simeq \\bar{\\omega}_q/l$, and at the symmetry point $\\epsilon_0 = 0$ the even-fraction peaks are suppressed by the harmonic content of the drive-induced coupling.","The weak-drive analytical solution (A1)-(A4) reproduces the numerical heat current from the instantaneous-basis equation across the frequency range for not-too-large drive amplitude, giving a closed-form set of Lorentzian resonances."],"supporting_citations":[{"why":"Supplies the projection-operator and Nakajima-Zwanzig formalism, including the second-order kernel and heat-current kernel used to derive Eqs. (28) and (35).","marker":"[14]"},{"why":"Defines the spin-boson model, the bath correlation function, and the modified spectral density used throughout the paper.","marker":"[20]"},{"why":"Provides the Floquet-Liouville supermatrix framework underlying the multi-photon resonance peaks at fractional drive frequencies.","marker":"[21]"},{"why":"Source for the fractional-frequency heat peaks and the odd-harmonic selection rule at the symmetry point.","marker":"[25]"},{"why":"Supplies the analytic expression for $W(\\omega)$ with an ohmic-Drude cutoff used in the numerical evaluation of the heat current.","marker":"[32]"},{"why":"Supplies the density-matrix formulation of quantum transport underlying the exact weak-coupling heat-current expression.","marker":"[38]"},{"why":"Gives the static Redfield+ equation that the steady-state Floquet-Redfield equations reduce to in the static limit.","marker":"[42]"},{"why":"Foundational for the instantaneous-eigenbasis master equation for a driven qubit used in Eqs. (55)-(57).","marker":"[47]"}],"fun_headline_variants":["Floquet-Redfield exact for driven qubit heat current","Secular approximation halves resonant peak in driven qubit","Nonsecular Floquet heat current: benchmark for driven qubits","Instantaneous eigenbasis matches Floquet only in adiabatic limit","Driven qubit heat flow: secular eq misses resonant peak by half"],"cache_read_input_tokens":22528,"weakest_assumption_plain":"The reported agreement of the instantaneous-eigenbasis master equation rests on borrowing relaxation rates computed in the Floquet basis, a step the paper itself calls inconsistent and valid only in the adiabatic limit $\\delta \\ll 1$.","fun_headline_variants_meta":{"raw":{"variants":["Floquet-Redfield exact for driven qubit heat current","Secular approximation halves resonant peak in driven qubit","Nonsecular Floquet heat current: benchmark for driven qubits","Instantaneous eigenbasis matches Floquet only in adiabatic limit","Driven qubit heat flow: secular eq misses resonant peak by half"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000184,"raw_usage":{"total_tokens":1301,"prompt_tokens":908,"completion_tokens":393,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":303}},"tokens_in":524,"tokens_out":393,"duration_ms":3905,"temperature":1.0,"reasoning_tokens":303,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:01:29.244632+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the steady-state heat current with the instantaneous-eigenbasis master equation (55)-(57) using relaxation rates evaluated directly in the instantaneous eigenbasis rather than transferred from the Floquet basis, at parameters where $\\delta = \\omega_d A_d / \\Delta^2$ is not small, e.g. $A_d = 0.4\\Delta$ and $\\omega_d = 1.5\\Delta$, giving $\\delta = 0.6$; if the resulting curve departs from the Floquet-Redfield benchmark of Eqs. (28) and (35), the reported agreement is an artifact of the rate transfer the paper itself flags as inconsistent.","supporting_citations":[{"cited_title":"Mori, Floquet states in open quantum systems, Annu","cited_arxiv_id":null,"evidence_quote":"Supplies the projection-operator and Nakajima-Zwanzig formalism, including the second-order kernel and heat-current kernel used to derive Eqs. (28) and (35)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the spin-boson model, the bath correlation function, and the modified spectral density used throughout the paper."},{"cited_title":"Weiss,Quantum Dissipative Systems(World Scientific, Singapore, 2012, 4th ed.)","cited_arxiv_id":null,"evidence_quote":"Provides the Floquet-Liouville supermatrix framework underlying the multi-photon resonance peaks at fractional drive frequencies."},{"cited_title":"Deppe, M","cited_arxiv_id":null,"evidence_quote":"Source for the fractional-frequency heat peaks and the odd-harmonic selection rule at the symmetry point."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the analytic expression for $W(\\omega)$ with an ohmic-Drude cutoff used in the numerical evaluation of the heat current."},{"cited_title":"Floquet, Sur la th´ eorie des ´ equations diff´ erentielles lin´ eaires, Ann","cited_arxiv_id":null,"evidence_quote":"Gives the static Redfield+ equation that the steady-state Floquet-Redfield equations reduce to in the static limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Foundational for the instantaneous-eigenbasis master equation for a driven qubit used in Eqs. (55)-(57)."}],"review_version":1}