{"id":"cc288a24-3e78-449f-83da-e4c1d4daea6d","arxiv_id":"2608.04558","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For weakly coupled driven fermionic systems, the steady-state occupation is a weighted sum of Fermi functions shifted by integer multiples of the driving frequency, with weights fixed by the micromotion operator.","lead":"The authors derive a formula for the steady-state electron occupation of a periodically driven quantum system coupled to a bath, generalizing the Fermi-Dirac distribution to driven, or Floquet, systems. It also yields a simple Landauer-type formula for the DC current, useful for interpreting transport in driven nano-devices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Appendix B's matrix Sokhotski–Plemelj step is not justified for non-commuting Γ, so the claimed robustness to arbitrary weakly coupled bath spectra is unsupported; the flat-bath result likely stands.","rationale":"The reader's verdict is CONDITIONAL, and I agree with that verdict. The reader's stated weakest assumptions (non-degeneracy modulo ℏΩ and site-symmetric coupling) are explicit and legitimate, but I see a more direct obstruction in the non-featureless-bath extension. The featureless-bath derivation, Eqs. (51)-(62), is coherent: the identity-proportional self-energy permits simultaneous diagonalization, the Sokhotski-Plemelj selection is valid under the stated non-degeneracy, and the final weighted-sum form has the correct normalization. The time-domain derivation in Appendix A provides additional support. The trouble begins when Γ(ω) is allowed to vary. Appendix B purports to prove the same result by a matrix generalization of Sokhotski-Plemelj, but Eq. (B10) is not the first-order inverse for non-commuting matrices: the trace-denominator term omits the -A^{-1}BA^{-1} correction. Since in the Floquet Keldysh expression the matrix X = ℏΩ - H_S does not commute with the diagonal Γ(ω), the diagonal-pole reduction δ((X - Λ)_{bb}) is not justified. A numerical Keldysh check with a nontrivial Γ(ω) would settle whether Eq. (62) survives; absent that, the abstract's 'broad class of bath spectral functions' claim is stronger than the proof. The appropriate verdict remains CONDITIONAL: retain Eq. (62) for featureless baths, but either prove the matrix SP identity correctly or restrict the claim to the case where the required commutation holds.","tokens_in":21499,"tokens_out":13742,"duration_ms":157155,"concrete_test":"Numerically evaluate the exact Keldysh lesser Green's function (67) for a two-level driven system, e.g. H_S(t)=εσ_z + A cos(Ωt)σ_x, with a Drude spectral function Γ(ω)=Γ_0/(1+ω²τ²), using truncated Floquet matrices (at least 9×9) for small Γ_0 and several τ. Compare the stroboscopic reduced density matrix at multiples of 2π/Ω with Eq. (62). In the flat-bath limit τ→0 the result should match Eq. (62), serving as a baseline; if any O(Γ_0) deviation or τ-dependence appears in the weak-coupling limit, the Appendix B generalization is invalid. A cheaper companion check is to verify Eq. (B10) by directly inverting 2×2 matrices with non-commuting X-Λ and Γ, where the exact first-order inverse is A^{-1} - (i/2)A^{-1}ΓA^{-1}, not A^{-1}[1+(i/2)Tr(ΓA^{-1})]^{-1}.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix B is the load-bearing weakness. The paper needs the matrix Sokhotski–Plemelj identity (B14) to extend Eq. (62) to frequency-dependent baths, but the derivation of (B14) contains an invalid first-order matrix inversion. Eq. (B10) approximates (X - Σ^R)^{-1} by (X - Λ)^{-1} / [1 + (i/2) Tr(Γ(X - Λ)^{-1})]. The exact expansion is (X - Λ)^{-1} - (i/2)(X - Λ)^{-1} Γ (X - Λ)^{-1} + O(Γ²); the trace-only denominator is not the first-order correction unless Γ commutes with X - Λ, which is precisely the case the appendix is meant to go beyond. Equivalently, Eq. (B10) drops the first-order correction to the adjugate matrix, adj(A+B) = adj(A) + [Tr(adj(A)B)I - adj(A)B] + …, and that dropped term is essential for non-commuting B. Consequently Eqs. (B12)-(B14) reduce the inverse to poles controlled by diagonal entries δ((X - Λ)_{bb}) in the basis where Γ is diagonal. For the actual Keldysh integrand (67), X = ℏΩ - H_S is not diagonal in that basis, so the inverse has a different pole structure. The final claim that the distribution is unchanged 'regardless of the specifics of Γ(ω)' is therefore unsupported by the argument as written. This does not affect the featureless-bath derivation leading to Eq. (62), which is internally consistent and cross-checked in Appendix A.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript embeds Floquet theory into the nonequilibrium Green's function formalism and derives an analytical steady-state occupation formula for a weakly driven open fermionic system. For a featureless bath, the retarded and advanced Green's functions factorize as P(...)^-1 P^dagger, and evaluation of the Keldysh equation in the Gamma -> 0+ limit gives the central result rho_F = sum_xi f(H_F + xi hbar Omega) D(P^dagger_xi P_xi), with rho_S(t) = U_S(t) rho_F U_S(t)^dagger. The derivation is cross-checked in Appendix A and reproduces the Bessel-weighted result for a sinusoidally driven resonant level, Eq. (65). The paper further claims that the formula is unchanged for a broad class of frequency-dependent bath spectral functions, using a matrix generalization of the Sokhotski-Plemelj formula in Appendix B, and derives a Floquet generalization of the Landauer formula, Eq. (81), for the DC current.","tokens_in":21756,"tokens_out":7940,"duration_ms":84961,"significance":"If Eq. (62) is correct, it provides a transparent, parameter-free description of Floquet sideband occupations for non-commuting, time-periodic Hamiltonians, going beyond the commuting case treated in Ref. [24]. The result has a clean physical interpretation through the Fourier components of the micromotion operator and reduces exactly to the Fermi-Dirac distribution in the undriven limit. The Landauer-type DC formula Eq. (81) is a testable prediction for two-terminal driven junctions. The featureless-bath derivation is internally consistent, with an independent time-domain derivation in Appendix A and explicit consistency with the known single-level Bessel result. However, the claimed robustness to arbitrary weakly coupled bath spectral functions is not established by the argument as written; only the featureless-bath result is fully supported in the present version.","major_comments":[{"comment":"The matrix inversion step leading to Eq. (B10) is not valid to first order in Gamma for non-commuting Gamma. Writing A = X - Lambda and B = i Gamma/2, the exact expansion is (A+B)^-1 = A^-1 - A^-1 B A^-1 + O(Gamma^2), whereas Eq. (B10) is equivalent to A^-1 / [1 + (i/2) Tr(Gamma A^-1)] = A^-1 - (i/2) A^-1 Tr(Gamma A^-1) + O(Gamma^2). The dropped term -(i/2) A^-1 Gamma A^-1 is of the same order and cannot be neglected unless Gamma commutes with A, which is precisely the case that Appendix B is intended to go beyond. Consequently, Eqs. (B12)-(B14) reduce the pole structure to delta((X-Lambda)_bb) in the eigenbasis of Gamma, but in the Keldysh integrand Eq. (67), X = hbar Omega - H_S does not generally commute with Gamma, so the claimed pole selection is not established. The statement in Section V that the final Floquet distribution is unchanged regardless of the specifics of Gamma(omega) is therefore unsupported by the derivation presented. The featureless-bath derivation of Eq. (62) is not affected by this issue.","section":"Appendix B, Eq. (B10)"},{"comment":"The central diagonal-selection step that converts the Keldysh equation into the closed form Eq. (62) relies on the assumption, stated before Eq. (58), that no two quasienergies differ by an integer multiple of hbar Omega. If this condition fails, the off-diagonal terms in Eq. (56) do not vanish in the Gamma -> 0+ limit, and the distribution acquires inter-sideband corrections. The authors explicitly note that the condition can fail for symmetry-protected degeneracies or a crowded Floquet-Brillouin zone, but the abstract and introduction claim the result for general periodically driven systems. The manuscript should state the precise theorem with the non-degeneracy condition explicitly listed, and should either discuss how Eq. (62) is modified in the degenerate case or clearly delimit the regime of validity.","section":"Section V, Eqs. (57)-(58)"}],"minor_comments":[{"comment":"The phrase \"the assumption that the spectral function universally commutes\" is undefined; since this section has already specialized to featureless baths, clarify what is being assumed and why it holds in that setting.","section":"Section VI, after Eq. (79)"},{"comment":"The sentence following Eq. (75) states that the first term yields the AC current and the second and third terms yield the DC current, but this is only explained after Eq. (82); moving the explanation earlier would improve readability.","section":"Section VI, Eq. (75)"},{"comment":"The notation Tr_t A(omega)B(omega) is nonstandard and initially confusing because the trace is over internal indices only while the Floquet indices are treated differently; the index structure should be defined explicitly at first use.","section":"Section VI, Eq. (71)"},{"comment":"The captions contain garbled control characters (e.g., sequences such as \"/uni00000014/uni00000013/...\") that appear to be rendering artifacts; these should be regenerated so that the captions are readable.","section":"Figures 2 and 3"},{"comment":"The phrase \"the desired form for applying the Sokhotski-Plemelj formula\" should be qualified, because the scalar-like form in Eq. (B10) is only a valid first-order inverse when Gamma commutes with X - Lambda.","section":"Eq. (B10)"}],"recommendation":"major_revision","confidential_remarks":"The main unsupported claim is localized to Appendix B and the corresponding statement in Section V that the Floquet distribution is robust to arbitrary weakly coupled bath spectral functions. If the authors can either repair the matrix Sokhotski-Plemelj argument under the stated assumptions or retract the broad-class claim and present Eq. (62) as specific to the featureless-bath limit, the revision would be acceptable. The flat-bath derivation appears sound and is independently cross-checked in Appendix A."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe one-sentence take: Eq. (62) is real, new, and mostly sound for featureless baths; the claim that it extends unchanged to arbitrary weakly coupled bath spectra is not supported by the argument as written.\n\nWhat's actually new: the non-commuting generalization. Ref. [24] treated diagonally commuting Hamiltonians; Refs. [25,26] rely on Markovian detailed balance. Here they get a closed-form steady state rho_F = sum_xi f(H_F + xi hbar Omega) D(P^dagger_xi P_xi) for general H(t) under the identity-proportional self-energy assumption. The DC Floquet Landauer formula, Eq. (81), is also clean and new. I spot-checked the central chain: Floquet representation, the Eq. (51) diagonalization via identity-proportional Sigma, the Sokhotski-Plemelj selection in Eq. (57), and the time-domain cross-check in Appendix A. All internally consistent. The reduction to the single-level Bessel result and the equilibrium limit are right. That part is worth a careful read.\n\nSoft spots, in proportion. The weakest assumptions—no quasienergy degeneracies mod hbar Omega, and each system level coupled to its own sub-bath—are stated explicitly, so no complaint there. The real problem is Appendix B. Eq. (B10) approximates (X - Sigma^R)^{-1} by (X - Lambda)^{-1} / [1 + (i/2) Tr(Gamma (X - Lambda)^{-1})]. That is not the first-order correction for non-commuting Gamma: the adjugate has a first-order term too, and the trace-only denominator ignores it. The exact expansion has (X - Lambda)^{-1} Gamma (X - Lambda)^{-1} with Gamma sandwiched, not traced. So the pole structure in Eq. (B14) is derived in a basis where Gamma is diagonal, and the claim that the distribution is unchanged regardless of Gamma(omega) is unsupported. The featureless-bath result does not depend on Appendix B, so the main theorem stands. But the abstract and conclusion explicitly advertise robustness to arbitrary bath spectra, which overreaches.\n\nMinor: gauge invariance of rho_F is asserted without proof; that is likely true but needs an argument.\n\nWho is this for? People working on Floquet open systems, quantum transport, and sideband occupations. The flat-bath distribution and the Landauer formula are useful. It deserves a serious referee: the central derivation is careful and the result is likely correct under its stated assumptions, but the bath-robustness claim needs to be either proved properly or softened.\n\nRecommendation: send to peer review, flag Appendix B as the main issue, and ask for a numerical test with an energy-dependent spectral function.","headline":"The flat-bath Floquet Fermi distribution is a genuine and well-derived result; the advertised generalization to arbitrary bath spectra rests on an invalid matrix Sokhotski–Plemelj step, so that part needs serious revision.","tokens_in":22406,"tokens_out":1976,"would_cite":true,"duration_ms":19594,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82C10","81Q15","82C70"],"pacs":["05.30.-d","05.60.Gg"],"model":"deepseek-v4-flash","headline":"This paper derives a closed analytical Floquet Fermi distribution for the steady state of weakly driven, bath-coupled fermionic systems, and a Floquet Landauer formula for the DC current.","keywords":["Floquet Fermi distribution","nonequilibrium Green's functions","micromotion operator","Landauer formula","periodically driven open quantum systems","Keldysh equation","Floquet sidebands"],"falsifier":"Drive a multi-level system with parameters chosen so that two Floquet quasienergies differ by exactly one photon energy $\\hbar\\Omega$ (or are otherwise degenerate in the Floquet–Brillouin zone), then evaluate the Keldysh equation (56) at small but finite $\\Gamma$; Eq. (58) then acquires off-diagonal $\\alpha \\neq \\gamma$ terms that Eq. (62) omits, so a numerically exact Floquet master-equation or Green's-function calculation would show occupations deviating from the weighted-sum formula.","tokens_in":21127,"feed_emoji":"⚛️","tokens_out":7154,"duration_ms":69025,"temperature":0.7,"pith_summary":"The paper derives a closed analytical formula for the steady-state occupation statistics of a weakly driven fermionic quantum system coupled to a heat bath. The result, called the Floquet Fermi distribution, expresses the steady-state occupation as a weighted sum of ordinary Fermi–Dirac functions whose energies are shifted by integer multiples of the driving frequency, with weights fixed by the Fourier components of the micromotion operator. The derivation proceeds through the nonequilibrium Green's function formalism in the Floquet representation and does not require the time-dependent Hamiltonian to commute with itself at different times. If correct, the formula gives a parameter-free description of Floquet sideband occupations and a Floquet generalization of the Landauer formula for DC transport.","feed_headline":"Floquet Fermi distribution fixes steady states of driven systems","feed_subtitle":"Green's-function derivation shows occupations are Fermi functions shifted by drive multiples, with micromotion weights.","key_machinery":"The central object is the Floquet representation of two-time Green's functions, a discrete Fourier transform in the average time and a continuous Fourier transform in the relative time, under which time convolutions become matrix multiplications via the Floquet convolution theorem. Combined with the micromotion operator $P(t)$ from Floquet's theorem, this representation diagonalizes the retarded and advanced Green's functions as $G^R_S = P(\\hbar\\Omega + i\\Gamma/2 - H_F)^{-1}P^\\dagger$ when the bath self-energy is proportional to the identity. Inserting the lesser self-energy $\\Sigma^< = i\\Gamma F$ into the Keldysh equation and taking the $\\Gamma \\to 0^+$ limit with the Sokhotski–Plemelj identity leaves only diagonal sideband terms, producing the weighted sum over $\\xi$. The weights $D(P^\\dagger_\\xi P_\\xi)$ are the squared column norms of the Fourier components of $P$, so they encode how much each Floquet replica contributes to the steady-state occupation.","core_discovery":"The central claim is that in the weak-coupling limit the steady-state density matrix of a periodically driven system in contact with a thermal fermionic bath takes the form $\\rho_F = \\sum_\\xi f(H_F + \\xi\\hbar\\Omega)\\, D(P^\\dagger_\\xi P_\\xi)$, where $f$ is the Fermi–Dirac function, $H_F$ is the Floquet Hamiltonian, $\\Omega$ is the driving frequency, $P_\\xi$ is the $\\xi$-th Fourier component of the micromotion operator $P(t)$, and $D(\\cdot)$ retains only diagonal entries in the eigenbasis of $H_F$. The time-dependent density matrix rotates as $\\rho_S(t) = U_S(t)\\rho_F U_S(t)^\\dagger$, so occupations are constant while the basis rotates when the Hamiltonian does not commute at different times. The formula holds for general non-commuting time-periodic Hamiltonians under weak coupling, extends to a broad class of frequency-dependent bath spectral functions, and reduces to the equilibrium Fermi–Dirac distribution when the drive vanishes. The paper further establishes a Floquet Landauer formula for the DC current, $\\langle J\\rangle_{\\mathrm{DC}} = \\frac{1}{\\hbar}\\frac{\\Gamma_L\\Gamma_R}{\\Gamma_L+\\Gamma_R}\\mathrm{Tr}(\\rho_L-\\rho_R)$, in which the equilibrium Fermi functions are replaced by their Floquet-modified counterparts.","pith_inferences":["Because the derivation takes the strict $\\Gamma \\to 0^+$ limit, the leading finite-broadening corrections should be of order $\\Gamma$ and proportional to the off-diagonal coherent sideband terms that Eq. (62) discards; an explicit first-order-in-$\\Gamma$ expression would give a practical correction for realistic couplings.","The formula suggests a direct experimental probe: spectroscopic reconstruction of the sideband weights $D(P^\\dagger_\\xi P_\\xi)$ in a driven quantum dot or cold-atom system would effectively measure the Fourier components of the micromotion operator.","A bosonic analogue obtained by replacing Fermi–Dirac functions with Bose–Einstein functions is a natural extension, though the fermionic signs and the diagonalization step would need to be re-derived for bosonic baths.","The gauge-invariance argument implies that the physically robust quantities are the sideband weights rather than the Floquet Hamiltonian itself, which may clarify how effective temperatures extracted from driven-state occupancies should be interpreted."],"forward_implications":["For any weakly driven multi-level system whose quasienergies are non-degenerate modulo $\\hbar\\Omega$, the steady-state occupations are determined entirely by $H_F$ and the Fourier coefficients of $P(t)$, with no need to solve a master equation.","When the Hamiltonian commutes at different times, $\\rho_S(t) = \\rho_F$ is time independent; in general the density matrix rotates with $U_S(t)$, so the occupations are constant but the basis rotates.","The result extends earlier diagonal commuting and single-sideband Floquet–Gibbs approximations to all Floquet sidebands, with weights that sum to unity by Parseval's theorem.","In a two-terminal setup, the DC current obeys a Landauer-type formula with Floquet distributions replacing the equilibrium Fermi functions, and increasing driving strength generally suppresses the current.","The same Floquet distribution survives for frequency-dependent bath spectral functions in the weak-coupling limit, not only for the featureless wide-band bath."],"supporting_citations":[{"why":"Supplies the exactly solvable diagonal commuting-Hamiltonian case and the notation that this paper extends to non-commuting H(t).","marker":"[24]"},{"why":"Provides the Floquet representation and convolution theorem used to turn the Keldysh equation into matrix multiplication in Floquet space.","marker":"[29]"},{"why":"Identifies the near-degeneracy problem that motivates the non-degenerate quasienergy assumption behind Eq. (58).","marker":"[35]"},{"why":"Defines periodic thermodynamics whose sideband rate-equation occupations the Floquet distribution Eq. (62) evaluates analytically.","marker":"[47]"},{"why":"Prior Floquet–Gibbs state based on a single-sideband detailed-balance condition, which the all-sideband distribution generalizes.","marker":"[25]"},{"why":"Two-terminal driven-transport setup and current expressions that motivate the Floquet Landauer formula.","marker":"[21]"},{"why":"Establishes the $Γ = i(Σ^R - Σ^A)$ convention and the undriven Landauer/Caroli baseline recovered in the zero-drive limit.","marker":"[44]"},{"why":"Caroli formula for tunneling current that the DC Floquet Landauer formula reproduces when driving corrections vanish.","marker":"[50]"}],"fun_headline_variants":["Floquet Fermi distribution governs driven open quantum systems","Steady states of driven systems obey Floquet Fermi function","Floquet Landauer formula extends transport to driven devices","Driven quantum systems show Floquet Fermi-Dirac occupancy","Analytic Floquet statistics yield Fermi functions with drive shifts"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on the assumption that no two quasienergies of the driven system differ by exactly an integer multiple of the driving quantum $\\hbar\\Omega$, so that off-resonant sideband terms can be dropped as the bath coupling vanishes, together with a symmetric coupling model in which each system level is identically attached to its own copy of the bath.","fun_headline_variants_meta":{"raw":{"variants":["Floquet Fermi distribution governs driven open quantum systems","Steady states of driven systems obey Floquet Fermi function","Floquet Landauer formula extends transport to driven devices","Driven quantum systems show Floquet Fermi-Dirac occupancy","Analytic Floquet statistics yield Fermi functions with drive shifts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000276,"raw_usage":{"total_tokens":1689,"prompt_tokens":1033,"completion_tokens":656,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":649,"completion_tokens_details":{"reasoning_tokens":572}},"tokens_in":649,"tokens_out":656,"duration_ms":7325,"temperature":1.0,"reasoning_tokens":572,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:07:59.060357+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Drive a multi-level system with parameters chosen so that two Floquet quasienergies differ by exactly one photon energy $\\hbar\\Omega$ (or are otherwise degenerate in the Floquet–Brillouin zone), then evaluate the Keldysh equation (56) at small but finite $\\Gamma$; Eq. (58) then acquires off-diagonal $\\alpha \\neq \\gamma$ terms that Eq. (62) omits, so a numerically exact Floquet master-equation or Green's-function calculation would show occupations deviating from the weighted-sum formula.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the exactly solvable diagonal commuting-Hamiltonian case and the notation that this paper extends to non-commuting H(t)."},{"cited_title":"Kohler, T","cited_arxiv_id":null,"evidence_quote":"Identifies the near-degeneracy problem that motivates the non-degenerate quasienergy assumption behind Eq. (58)."},{"cited_title":"Haughian, M","cited_arxiv_id":null,"evidence_quote":"Defines periodic thermodynamics whose sideband rate-equation occupations the Floquet distribution Eq. (62) evaluates analytically."},{"cited_title":"Matsyshyn, J","cited_arxiv_id":null,"evidence_quote":"Prior Floquet–Gibbs state based on a single-sideband detailed-balance condition, which the all-sideband distribution generalizes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Two-terminal driven-transport setup and current expressions that motivate the Floquet Landauer formula."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Caroli formula for tunneling current that the DC Floquet Landauer formula reproduces when driving corrections vanish."}],"review_version":1}