{"id":"892a41df-e374-4170-bb98-26c802de7923","arxiv_id":"2604.01556","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"A generalized TEMPO algorithm builds the Feynman-Vernon influence functional as a matrix product operator for off-diagonal, non-Hermitian system-bath couplings, and shows the secular approximation can fail for sub-ohmic baths.","lead":"This paper extends TEMPO, a numerical method for simulating a small quantum system interacting with a surrounding bosonic bath, to off-diagonal couplings where the system operator is non-Hermitian. The new version lets researchers study a broader class of quantum impurity models, and shows that the commonly used secular (rotating-wave) approximation can fail badly for structured baths.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No exact benchmark for the JC spin-boson results: the fixed n=20 exponential-sum hybridization fit is never tested against the analytically solvable single-excitation sector, so δt/χ convergence cannot detect a systematic MPO-IF error.","rationale":"The paper's formal derivation is careful, and the two validation benchmarks (JC single mode, noninteracting boson) give real support. The code is publicly available, which is a plus. However, the central physical conclusion rests on the JC spin-boson computation, and the only testable exact benchmark for that model—the single-excitation analytic solution—is not used. The convergence studies vary δt and χ but not the fixed exponential-sum order n; hence they cannot rule out a systematic error in the MPO-IF. This is not an accusation of error; it is an identification of the least-supported step in the chain from derivation to physical claim. The concern is addressable by a straightforward numerical experiment, so the conditional verdict is appropriate rather than rejection.","tokens_in":21675,"tokens_out":16924,"duration_ms":176314,"concrete_test":"Using the published Julia code, recompute Fig. 7 for α=0.01, 0.04, 0.08, 0.12 with Prony fit order n=20 and n=40 (and, if feasible, n=10) at fixed δt=0.0125 and converged χ, and compare all runs against the exact single-excitation-sector solution of the JC spin-boson model (or a fine-discretization ED of the equivalent Lee model). If the n=20 and n=40 results differ by more than a few percent in ⟨σ_z⟩(t), or if the n=40 result deviates from the exact solution, then the n=20 assumption in Appendix C is the limiting error and the secular-approximation claim is not yet supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing premise is that the discrete off-diagonal Feynman-Vernon IF is faithfully captured by the MPO-IF built with a fixed n=20 exponential-sum approximation of the hybridization function (Appendix C, Eq. C1) and XTRG with m=7. The δt and χ convergence studies in Figs. 8-9 validate only Trotter error and bond-truncation error; the n=20 fitting error is independent of both and is never quantified. The one place this matters most is the central physical application: the JC spin-boson model in the single-excitation sector, which the authors state is analytically solvable (Sec. IV.C). Instead of comparing against that solution, Fig. 7 compares the extended TEMPO data against standard TEMPO, which is itself not fully converged at α=0.12. If n=20 gives an inaccurate representation of the sub-ohmic s=0.5 hybridization at the longest times or strongest couplings studied, the observed ⟨σ_z⟩ revivals and the conclusion that the secular approximation is very poor would not be established. The code is available, so this is checkable. A secondary technical point is that the equal-time ordering corrections of Appendix B are not shown explicitly in the generic site tensor C2 for arbitrary A; but the tested cases avoid the general issue because for σ± and a the correction reduces to A†A plus a constant, so I do not make this the primary concern.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the time-evolving matrix product operator (TEMPO) method to bosonic quantum impurity problems with off-diagonal system-bath coupling of the form H_hyb = sum_k (V_k A b_k^\\dagger + H.c.), where A may be non-Hermitian. The authors derive the Feynman-Vernon influence functional on the Keldysh contour for this general coupling, discretize it, and represent it as an MPO built with the exponential thermal tensor renormalization group (XTRG) algorithm using an exponential-sum approximation of the hybridization function. They validate the method against exact diagonalization in two solvable cases: the Jaynes-Cummings model with a single bosonic mode and a noninteracting bosonic impurity coupled to a sub-ohmic bath. As a physical application, they compute the real-time spin dynamics for a Jaynes-Cummings-type spin-boson model with a sub-ohmic bath and compare it with the standard spin-boson model, concluding that the secular (rotating-wave) approximation can be very poor even at weak coupling. The paper also claims generality over earlier TEMPO variants and suggests a fermionic extension.","tokens_in":21993,"tokens_out":15713,"duration_ms":144556,"significance":"If the central claims hold, this is a useful and fairly general extension of TEMPO: it unifies previously separate treatments of diagonal/commutative and diagonal/non-commutative couplings, and it removes the need to discretize the bath or truncate its local Hilbert space. Two strengths deserve emphasis: (i) the method is validated against external exact-diagonalization benchmarks in two distinct settings, including a continuous sub-ohmic bath; and (ii) the authors provide a public code repository, which makes the numerical results reproducible and the remaining concerns checkable. The potential payoff for bosonic dynamical mean-field theory and for simulating JC-type impurity models gives the work significance beyond a purely incremental algorithmic step. However, as detailed below, some written parts of the derivation and the advertised scope contain inconsistencies that need to be resolved before the paper is publication-ready.","major_comments":[{"comment":"The equal-time ordering correction derived in Appendix B (Eqs. B7-B10) replaces the naive term e^{-A†_+ Δ^{++}_{j,j} A_+} by an exponent containing both A†A and A A† with different coefficients. The generic site tensor in Eq. (C2), however, contains only a single local operator Δ0 A†A. For arbitrary non-Hermitian A at finite β, A A† is not proportional to A†A, so the MPO-IF as written does not implement the correction. The examples avoid the issue (β=∞, or low-T with the AA† coefficient exponentially small), but the claim in Sec. III.A that the method applies to any linearly-coupled noninteracting bath is then too strong. Please specify how the two ordering terms are encoded in the MPO, or state the restrictions.","section":"Appendix B and Appendix C, Eq. (C2)"},{"comment":"The factor (e^{-βω}-1)^{-1} appears in Δ^{+-}_{j,k} and in the equal-time A A† term. The thermal occupation of the bath is n_B(ω)=(e^{βω}-1)^{-1}, which vanishes at zero temperature. The printed factor instead tends to -1 as β→∞, implying a spurious nonzero contribution from ⟨b†b⟩=0 in the vacuum state. If the formulas as written are used to construct the IF, results at β=∞ would be incorrect; presumably the code uses the correct sign, but the appendix must be corrected or the convention explained.","section":"Appendix B, Eqs. (B5), (B8)"},{"comment":"The abstract states that the paper \"study[ies] the imaginary-time evolution of a bosonic impurity with nonzero on-site interaction that is coupled to a sub-ohmic bath.\" No such imaginary-time simulation appears in Section IV, which contains only real-time results (JC model, free bosonic impurity, JC spin-boson model). Either add the promised imaginary-time study or remove the sentence from the abstract.","section":"Abstract vs. Sec. IV"},{"comment":"The JC spin-boson model is stated to be analytically solvable in the single-excitation sector (β=∞, ρ_S=|e⟩⟨e|). Fig. 7 compares only against standard TEMPO, which itself is not fully converged at α=0.12. A direct comparison with the exact single-excitation solution would provide a stronger benchmark for the new MPO-IF and would directly validate the physical conclusion that the secular approximation is very poor. The existing ED benchmarks for the same spectral function in Sec. IV.B reduce but do not eliminate this concern.","section":"Sec. IV.C, Fig. 7"}],"minor_comments":[{"comment":"The phrase \"communication relation\" should be \"commutation relation\".","section":"Sec. III.D"},{"comment":"The exponential-sum approximation with n=20 is asserted as sufficient without a fit-error plot. A small panel showing the maximum relative error of the fitted hybridization function over the relevant time/frequency range would make the fixed n=20 choice transparent and would directly address the n-convergence question.","section":"Appendix C, Eq. (C1)"},{"comment":"The JC Hamiltonian uses coefficients 2Ω and 2λ without explanation. A sentence connecting this convention to the more standard notation λ and Ω would help readers.","section":"Sec. IV.A, Eq. (16)"},{"comment":"The symbol ̃Tr for the modified trace of the process tensor is used before being explicitly defined. Please define it at first use.","section":"Sec. II, Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The core algorithmic idea is sound and the numerical validations give reasonable confidence in the method, but the manuscript's written appendices contain a sign/thermal-factor issue and an incomplete description of the equal-time MPO site tensor. These are fixable but are load-bearing for the advertised generality. The missing imaginary-time results in the abstract and the absence of the analytic single-excitation benchmark are additional concerns that should be addressed before acceptance. The code availability is a strong positive and makes the requested checks feasible."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Good paper. The extension of TEMPO to off-diagonal system-bath couplings is real, and the construction is more general than anything else in the literature they cite: it works for non-Hermitian A, reduces to the diagonal cases, and is based only on the Trotter decomposition of the Feynman-Vernon IF. The derivation in Appendix A is a standard Gaussian cumulant expansion, done carefully, and the equal-time ordering corrections in Appendix B look right. The two benchmarks — JC model with a single mode and a noninteracting bosonic mode coupled to a sub-ohmic bath — are convincing. The code is on GitHub, so the method is checkable.\n\nThe soft spots are proportionate. The n=20 exponential-sum fit to the hybridization (Eq. C1) is a free parameter that is never tested. The δt and χ convergence studies in Figs. 8-9 validate Trotter error and bond truncation, but not the fitting error. For the main physical application, the JC spin-boson model, the paper itself states the single-excitation sector is analytically solvable. Instead of comparing to that solution, Fig. 7 compares to standard TEMPO, which at α=0.12 is not fully converged even at the largest χ and smallest δt they used. That means the observed revivals and the claim that the secular approximation is very poor for structural baths are not pinned down by a controlled benchmark. This is not a fatal flaw: the method is validated on other solvable cases, and the weak-coupling (α=0.01) result is likely robust since both codes converge there. But it is an addressable omission, and the authors should close it before publication. A second, minor point: the generic site tensor in Eq. C2 does not explicitly show the equal-time ordering corrections for arbitrary A; the tested cases avoid the issue because for σ± and a the correction reduces to A†A plus a constant. Worth a sentence.\n\nThe reader's conditional verdict is about right. The stress-test note is fair: the n=20 fit is an uncontrolled approximation in the current manuscript. But it does not undercut the method's correctness, given the two clean ED benchmarks and the code release.\n\nThis paper deserves a serious referee. It is a real advance for the TEMPO community and a plausible building block for BDMFT impurity solvers. The referee should ask for the single-excitation analytic benchmark and an n-convergence test, not for a rewrite from scratch.","headline":"Solid extension of TEMPO to off-diagonal couplings with sound benchmarks; the main gap is an untested n=20 hybridization fit and a missing analytic benchmark for the flagship application.","tokens_in":22497,"tokens_out":2933,"would_cite":true,"duration_ms":27443,"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":"This paper extends the TEMPO tensor-network method to off-diagonal system-bath couplings, where the system operator may be non-Hermitian, by building the Feynman-Vernon influence functional as a matrix product operator.","keywords":["time-evolving matrix product operators","off-diagonal system-bath coupling","Feynman-Vernon influence functional","process tensor","spin-boson model","Jaynes-Cummings coupling","sub-ohmic bath","secular approximation"],"falsifier":"Apply the algorithm to a bath with two widely separated narrow spectral peaks, so the hybridization function clearly needs more than the n=20 exponential terms, and compare the predicted spin dynamics against exact diagonalization; if the results fail to converge as bond dimension and XTRG steps increase, the exponential-sum representation is the cause.","tokens_in":21561,"feed_emoji":"⚛️","tokens_out":5040,"duration_ms":45156,"temperature":0.7,"pith_summary":"The paper sets out to generalize TEMPO, a tensor-network path-integral method for quantum impurity problems, from diagonal system-bath couplings to off-diagonal couplings of the form sum_k (V_k A b_k^dagger + H.c.) with A possibly non-Hermitian. The key move is to treat the discretized Feynman-Vernon influence functional as the thermal state of an effective one-dimensional quantum Hamiltonian and build it as a matrix product operator, so the construction depends only on the Trotter decomposition and the hybridization function, not on the details of the system operator. The authors validate the method on two exactly solvable models and then use it to show that the rotating-wave (secular) approximation can badly fail for a spin coupled to a sub-ohmic bath, even at weak coupling. A sympathetic reader would care because this supplies a unified route to a broad class of bosonic impurity problems, including the type central to bosonic dynamical mean-field theory.","feed_headline":"New TEMPO extension solves off-diagonal system-bath couplings","feed_subtitle":"Same influence functional handles non-Hermitian couplings; shows the rotating-wave approximation can break down.","key_machinery":"The central object is the MPO-IF: a matrix product operator representation of the discretized Feynman-Vernon influence functional I[A†,A] = exp(−∫dτ∫dτ′ A†(τ)Δ(τ,τ′)A(τ′)). The construction interprets the discretized exponent as the Hamiltonian of an effective one-dimensional system at inverse temperature 1, builds e^{−H_eff/2^m} with a short-step exponential MPO algorithm, and squares it m times via exponential thermal tensor-network renormalization. A crucial ingredient is the exponential-sum approximation of the hybridization function (n=20 terms in this work), whose parameters enter the MPO site tensor and keep the bond dimension compact. Appendix B supplies the equal-time operator-order","core_discovery":"The central claim is that the Feynman-Vernon influence functional for off-diagonal system-bath coupling, after QuAPI discretization on the Keldysh contour, is the exponential of an effective long-range quantum Hamiltonian, and therefore can be represented as a matrix product operator using exponential thermal tensor-network techniques. This MPO acts as the process tensor of the impurity, and contracting it with system propagators yields all multi-time correlations. Because the construction rests only on the Trotter decomposition of the influence functional and not on the specifics of the coupling operator, it automatically contains vanilla TEMPO and its diagonal and non-commutative extension","pith_inferences":["Editorial inference: The equal-time operator-ordering correction in Appendix B is likely essential for any non-commuting system operator, and its necessity could be tested by benchmarking against pseudo-mode or chain-mapping methods on the same model.","Editorial inference: The method's exponential-sum approximation of the hybridization function suggests that spectral functions with multiple sharp peaks or very low-frequency structure may require more than 20 terms, so the claimed n=20 sufficiency should be checked case by case.","Editorial inference: Because the process tensor encodes the bath's influence on the system, the same MPO-IF could be recycled to study different system Hamiltonians or different initial system states without recomputation, a practical advantage the paper notes but does not demonstrate numerically.","Editorial inference: The demonstrated failure of the secular approximation for sub-ohmic baths suggests that previous spin-boson studies using the rotating-wave approximation in structured or sub-ohmic environments may need to be revisited, though this is an inference beyond the paper's own claims."],"forward_implications":["Any bosonic impurity with linear coupling to a noninteracting bath, Hermitian or not, can in principle be solved by the same influence-functional MPO, removing the diagonal-coupling restriction of earlier TEMPO variants.","The unified construction provides a practical impurity solver for bosonic dynamical mean-field theory in the normal phase with a scalar hybridization function.","The numerical comparison indicates that the rotating-wave approximation can be quantitatively unreliable for sub-ohmic baths even in the weak-coupling regime, so conclusions drawn from RWA spin-boson models should be rechecked.","The same strategy suggests a direct fermionic extension, since the discretized influence functional for fermions has the same exponential structure, a direction not yet explored in the paper.","The explicit preservation of time-translational invariance opens the door to infinite-time and steady-state versions of the method for bosonic problems."],"fun_headline_variants":["TEMPO now handles off-diagonal system-bath coupling","Off-diagonal coupling tackled with time-evolving MPO","Extending TEMPO to non-Hermitian couplings","TEMPO exposes rotating-wave approximation failures","Unified TEMPO for off-diagonal system-bath coupling"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"Everything rests on the assumption that a bath's effect can be captured by about twenty exponential decay terms and that the resulting tensor network stays compact; if a real bath needs more terms or generates more entanglement than the examples shown, the method's efficiency and accuracy are not guaranteed.","fun_headline_variants_meta":{"raw":{"variants":["TEMPO now handles off-diagonal system-bath coupling","Off-diagonal coupling tackled with time-evolving MPO","Extending TEMPO to non-Hermitian couplings","TEMPO exposes rotating-wave approximation failures","Unified TEMPO for off-diagonal system-bath coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000196,"raw_usage":{"total_tokens":1214,"prompt_tokens":779,"completion_tokens":435,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":523,"completion_tokens_details":{"reasoning_tokens":359}},"tokens_in":523,"tokens_out":435,"duration_ms":4180,"temperature":1.0,"reasoning_tokens":359,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T16:53:25.425658+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Apply the algorithm to a bath with two widely separated narrow spectral peaks, so the hybridization function clearly needs more than the n=20 exponential terms, and compare the predicted spin dynamics against exact diagonalization; if the results fail to converge as bond dimension and XTRG steps increase, the exponential-sum representation is the cause.","supporting_citations":[],"review_version":1}