pith. sign in

arxiv: 2606.30465 · v1 · pith:KA4Z2YUHnew · submitted 2026-06-29 · ❄️ cond-mat.mes-hall

Heat rectification through a quantum two-level system

Pith reviewed 2026-06-30 04:51 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall
keywords heat rectificationspin-boson modelKondo temperaturequantum two-level systemtensor networkOhmic bathsteady-state heat currentinfrared fixed point
0
0 comments X

The pith

Heat rectification through an asymmetrically coupled quantum two-level system follows a universal power law near its infrared fixed point.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies steady-state heat flow in the Ohmic spin-boson model, where a two-level system couples asymmetrically to two thermal baths. It identifies a scaling regime around the Kondo temperature in which the system flows from an ultraviolet regime at higher temperatures to an infrared fixed point at lower temperatures. Perturbation theory applied near that infrared fixed point predicts that the rectification ratio obeys a universal power law, and tensor-network computations of the heat current confirm this prediction across the scaling window. The work shows how dissipation-induced correlations control thermal transport in a quantum impurity setting. A reader would care because the result supplies a concrete, parameter-free signature that links many-body physics to measurable heat rectification.

Core claim

The authors evaluate the steady-state heat current in the Ohmic spin-boson model with tensor networks and benchmark it against analytic limits. They locate a scaling regime that runs from the ultraviolet regime above the Kondo temperature to the infrared regime below it. Perturbation theory around the infrared fixed point then yields a universal power-law form for the rectification ratio, which the numerical data reproduce.

What carries the argument

Perturbation theory around the infrared fixed point of the Ohmic spin-boson model, which produces the universal power-law expression for the rectification ratio.

If this is right

  • The rectification ratio becomes independent of microscopic parameters once the system enters the infrared regime.
  • Tensor-network results match the analytic power-law prediction obtained from infrared perturbation theory.
  • Distinct transport regimes (weak coupling, incoherent tunneling, and strongly correlated scaling) produce qualitatively different rectification behavior.
  • Dissipation-induced many-body effects dominate the heat current once temperatures drop below the Kondo scale.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same infrared scaling might govern rectification in other quantum impurity models with bosonic baths.
  • Temperature-tuned measurements around the Kondo scale could serve as a diagnostic for the fixed-point structure in mesoscopic thermal devices.
  • The power-law form supplies a target signature for experiments that engineer asymmetric couplings in superconducting circuits or quantum dots.

Load-bearing premise

The tensor-network computation accurately captures the steady-state heat current throughout the scaling regime that flows from the ultraviolet to the infrared fixed point around the Kondo temperature.

What would settle it

A plot of rectification ratio versus temperature that deviates from the predicted power law in the regime below the Kondo temperature would falsify the central claim.

Figures

Figures reproduced from arXiv: 2606.30465 by Manuel Houzet, Tsuyoshi Yamamoto.

Figure 1
Figure 1. Figure 1: FIG. 1. Quantum heat transport through a two-level sys [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Schematic overview of the regimes considered in this [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Linear thermal conductance [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Heat rectification ratio [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Heat rectification ratio [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. MPO representation of the transfer matrix [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Time evolution of the correlation function [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
read the original abstract

We study heat rectification through a quantum two-level system asymmetrically coupled to two thermal baths, as described by the Ohmic spin-boson model. We evaluate the steady-state heat current using a tensor-network approach, which enables us to access the strongly correlated regime, and benchmark the results against analytical formulas in several limiting regimes, including the weak-coupling and incoherent-tunneling regimes. We identify a scaling regime where the studied system flows from an ultraviolet regime, at temperatures larger than the Kondo temperature, to an infrared regime, at temperatures lower than the Kondo temperature. By applying perturbation theory near the infrared fixed point, we find that the rectification ratio follows a universal power law. Our numerical results agree well with this analytical prediction. Our results provide a fundamental understanding of how dissipation-induced many-body physics affects heat transport.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The paper examines heat rectification in the Ohmic spin-boson model with asymmetric bath couplings. Tensor-network methods are used to compute the steady-state heat current in the strongly correlated regime, with benchmarks against analytic limits in weak-coupling and incoherent-tunneling regimes. A scaling regime is identified flowing from UV (T ≫ TK) to IR (T ≪ TK) fixed points around the Kondo temperature. Perturbation theory near the IR fixed point yields a universal power law for the rectification ratio, with reported numerical agreement.

Significance. If the tensor-network results reliably capture the scaling window, the work provides a concrete link between many-body Kondo physics and heat rectification, with potential implications for quantum thermal devices. The combination of an analytic power-law derivation (independent of fitting) and numerical validation is a positive feature.

major comments (1)
  1. [Numerical methods / tensor-network section (as described in abstract)] The central claim of numerical agreement with the analytic power law rests on the tensor-network method accurately extracting the steady-state heat current throughout the UV-to-IR scaling regime around TK. The abstract positions this method as essential for the strongly correlated regime, yet no explicit convergence tests, truncation-error bounds, or comparisons to known limits inside the scaling window are referenced in the provided description; this leaves the validation inconclusive.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on the numerical validation. We address the single major comment below.

read point-by-point responses
  1. Referee: [Numerical methods / tensor-network section (as described in abstract)] The central claim of numerical agreement with the analytic power law rests on the tensor-network method accurately extracting the steady-state heat current throughout the UV-to-IR scaling regime around TK. The abstract positions this method as essential for the strongly correlated regime, yet no explicit convergence tests, truncation-error bounds, or comparisons to known limits inside the scaling window are referenced in the provided description; this leaves the validation inconclusive.

    Authors: We agree that the abstract (and the summary description) does not explicitly reference convergence tests, truncation-error bounds, or direct comparisons inside the scaling window, which leaves the validation of the central claim less conclusive than it could be. The manuscript does benchmark the tensor-network results against analytic formulas in the weak-coupling and incoherent-tunneling regimes, but these lie outside the strongly correlated scaling regime of interest. To strengthen the paper, we will add an appendix (or expanded methods section) that reports explicit convergence tests with respect to bond dimension and time-step, truncation-error estimates, and additional comparisons to known limits or perturbative results inside the UV-to-IR window around TK. This revision will directly address the concern and make the numerical support for the universal power law more robust. revision: yes

Circularity Check

0 steps flagged

Analytical perturbation theory yields independent power-law prediction

full rationale

The paper derives the universal power-law rectification ratio via perturbation theory applied near the infrared fixed point of the Ohmic spin-boson model. This analytical step is performed after identifying the UV-to-IR scaling regime around the Kondo temperature and is benchmarked against separate limiting-case formulas (weak-coupling, incoherent tunneling). The tensor-network numerics are used only for validation in the strongly correlated regime and do not enter the derivation of the power law itself. No self-citations, fitted parameters renamed as predictions, or self-definitional steps are present in the reported chain; the central claim therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the Ohmic spin-boson model with asymmetric bath coupling, the accuracy of the tensor-network method in the strongly correlated regime, and the validity of perturbation theory around the infrared fixed point. No explicit free parameters are identified in the abstract as being fitted to data; the Kondo temperature functions as a crossover scale rather than a fitted constant. No new entities are postulated.

axioms (2)
  • domain assumption The Ohmic spin-boson model with asymmetric coupling to two baths describes the essential physics of the heat-rectifying two-level system.
    Invoked at the outset of the abstract as the model under study.
  • domain assumption The tensor-network method yields the correct steady-state heat current in the strongly correlated regime.
    Required for the numerical results that are then compared to the analytic power law.

pith-pipeline@v0.9.1-grok · 5659 in / 1566 out tokens · 64923 ms · 2026-06-30T04:51:38.728598+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

69 extracted references · 2 canonical work pages

  1. [1]

    = ∆/π(see Eq. (23)). While the leading contributionJ 2 does not induce heat rectification, the next-leading contributionJ (a) 4 can give rise to rectification, since it originates from the anhar- monic term∝˜ρ 4 in Eq. (32b). By contrast, the contri- butionJ (b) 4 does not induce rectification either, because it comes from a harmonic term in the field∝˜ρ∂...

  2. [2]

    Feynman-Vernon path integral First, we show the real-time dynamics of the two-level system using the Feynman-Vernon path-integral formal- ism [44]. The time-evolution of the reduced density ma- trix of the two-level system,ρ(t) = tr B[ϱ(t)], whereϱ(t) is the density matrix of the global system and tr B[O] de- notes tracing over with respect to the degree ...

  3. [3]

    To overcome this, we adopt the tensor network algorithm, the time-evolving matrix product operators (TEMPO) [43]

    TEMPO The memory effects caused by the coupling to the bath prevent us from computing the long-time dynam- ics. To overcome this, we adopt the tensor network algorithm, the time-evolving matrix product operators (TEMPO) [43]. Introducing the augmented density tensor (ADT) Ay0,...,yM = MY m≤n=0 F yn,ym[ρ]y0 ,(A7) whereF yn,ym is the renormalized influence ...

  4. [4]

    Steady-state heat current The steady-state heat current (14) is rewritten as J= αLαR α Z ∞ 0 dt χ(t)[F(t, T L)−F(t, T R)],(A17) whereF(t, T r) is the window function F(t, T r) = Z ∞ 0 dω ω 2e−ω/ωc sin(ωt)nr(ω) =T 3 r Im h ψ(2)(1 + (it+ω −1 c )Tr) i ,(A18) whereψ (2)(z) =d 2ψ(z)/dz 2 is the polygamma func- tion of order 2, and it decays asF(t, T L)−F(t, T ...

  5. [5]

    ContributionJ 2 The commutation [H L 0 , V2] is [H0 L, V2] =−ig 2[α1(:ρ L ˙ρL : + : ˙ρLρL :) + 2√αLαR : ˙ρLρR :],(B6) where we writeρ r = ˜ρr(x= 0) for simplicity and g2 =πc 2TK/α(v/TK)2. Noting that, for the free bosonic Hamiltonian, the different modes are separable, ⟨OLOR⟩0 =⟨O L⟩L ⟨OR⟩R and the odd-order moments are zero,⟨ρ 2n+1 r ⟩r = 0, the steady-s...

  6. [6]

    ContributionJ (a) 4 The steady-state heat currentJ (a) 4 is obtained, in the same manner asJ 2, as J(a) 4 =−8ig 2g4αLαR × Z ∞ −∞ dt n αL ⟨[: ˙ρL(0)ρR(0) :,:ρ 3 L(t)ρR(t) :]⟩0 +α R ⟨[: ˙ρL(0)ρR(0) :,:ρ L(t)ρ3 R(t) :]⟩0 = 24ig2g4αLαR(αL ⟨:ρ 2 L :⟩L +α R ⟨:ρ 2 R :⟩R) × Z ∞ −∞ dt h ˙¯SL(t) ¯SR(t)− ˙SL(t)SR(t) i = 48g2g4αLαR(αL ⟨:ρ 2 L :⟩L +α R ⟨:ρ 2 R :⟩R) × ...

  7. [7]

    In the Fourier space, it reads J(b) 4 =−4g 2g3αLαR × Z ∞ −∞ dω 2π ω n SL(ω)DR(−ω)−D L(−ω)SR(ω) o , (B13) where we used ¯Dr(ω) =D r(−ω)

    ContributionJ (b) 4 The contribution fromV (b) 4 is calculated as J(b) 4 =−2ig 2g3αLαR × Z ∞ −∞ dt n ⟨[: ˙ρL(0)ρR(0) :,;ρ L(t)∂2 xρR(t) :]⟩0 +⟨[: ˙ρL(0)ρR(0) :,:∂ 2 xρL(t)ρR(t) :]⟩0 =−2ig 2g3αL X r αr Z ∞ −∞ dt[ ˙SL(t)Dr(t) + ˙DL(t)Sr(t) − ˙¯SL(t) ¯Dr(t)− ˙¯DL(t) ¯Sr(t)], (B12) withg 3 =−πc 4TK/(2α)(v/TK)4 and the correla- tion functions,D r(t) =⟨∂ 2 xρr(...

  8. [8]

    Spectral function Now, we calculate the spectral functionA r(ω). The bosonic fields in the free bosonic Hamiltonian at the IR fixed point with the lengthLare expressed as ˜ϕr(x) = X n>0 r π knL cos(knx)(akn +a † kn ),(B15a) ˜ρr(x) = X n>0 1 i r kn πL cos(knx)(akn −a † kn ),(B15b) wherek n =nπ/L. Hence, using these fields, the spectral function is Ar(ω) = ...

  9. [9]

    Fluctuations The auto-correlation function of ˜ρr(x= 0) =ρ r at the same time is, using the spectral function (B16), ⟨ρ2 r⟩=S r(t= 0) = Z ∞ −∞ dω 2π Ar(ω) 1−e −ω/Tr = 1 (πv)2 Z ∞ 0 dω ω[2n r(ω) + 1].(B17) Therefore, we obtain the fluctuations ofρ r as ⟨: ˜ρ2 r(x= 0) :⟩= T 3 3v2 .(B18) By plugging the spectral function (B16) and the fluc- tuations (B18) in...

  10. [10]

    J. P. Pekola, Towards quantum thermodynamics in elec- tronic circuits, Nat. Phys.11, 118 (2015)

  11. [11]

    J. P. Pekola and B. Karimi, Colloquium: Quantum heat transport in condensed matter systems, Rev. Mod. Phys. 93, 041001 (2021)

  12. [12]

    J. B. Pendry, Quantum limits to the flow of information and entropy, J. Phys. A: Math. Gen.16, 2161 (1983)

  13. [13]

    Schwab, E

    K. Schwab, E. A. Henriksen, J. M. Worlock, and M. L. Roukes, Measurement of the quantum of thermal con- ductance, Nature404, 974 (2000)

  14. [14]

    Jezouin, F

    S. Jezouin, F. D. Parmentier, A. Anthore, U. Gennser, A. Cavanna, Y. Jin, and F. Pierre, Quantum Limit of Heat Flow Across a Single Electronic Channel, Science 342, 601 (2013)

  15. [15]

    Meschke, W

    M. Meschke, W. Guichard, and J. P. Pekola, Single-mode heat conduction by photons, Nature444, 187 (2006)

  16. [16]

    A. V. Timofeev, M. Helle, M. Meschke, M. M¨ ott¨ onen, and J. P. Pekola, Electronic Refrigeration at the Quantum Limit, Phys. Rev. Lett.102, 200801 (2009)

  17. [17]

    Giazotto, T

    F. Giazotto, T. T. Heikkil¨ a, A. Luukanen, A. M. Savin, and J. P. Pekola, Opportunities for mesoscopics in ther- mometry and refrigeration: Physics and applications, Rev. Mod. Phys.78, 217 (2006)

  18. [18]

    Giazotto and M

    F. Giazotto and M. J. Mart´ ınez-P´ erez, The josephson 14 heat interferometer, Nature492, 401 (2012)

  19. [19]

    Ronzani, B

    A. Ronzani, B. Karimi, J. Senior, Y.-C. Chang, J. T. Peltonen, C. Chen, and J. P. Pekola, Tunable photonic heat transport in a quantum heat valve, Nat. Phys.14, 991 (2018)

  20. [20]

    Maillet, D

    O. Maillet, D. Subero, J. T. Peltonen, D. S. Golubev, and J. P. Pekola, Electric field control of radiative heat transfer in a superconducting circuit, Nat. Commun.11, 4326 (2020)

  21. [21]

    Senior, A

    J. Senior, A. Gubaydullin, B. Karimi, J. T. Peltonen, J. Ankerhold, and J. P. Pekola, Heat rectification via a superconducting artificial atom, Commun. Phys.3, 40 (2020)

  22. [22]

    Gubaydullin, G

    A. Gubaydullin, G. Thomas, D. S. Golubev, D. Lvov, J. T. Peltonen, and J. P. Pekola, Photonic heat transport in three terminal superconducting circuit, Nat. Commun. 13, 1552 (2022)

  23. [23]

    Upadhyay, B

    R. Upadhyay, B. Karimi, D. Subero, C. D. Satrya, J. T. Peltonen, Y.-C. Chang, and J. P. Pekola, Strong-coupling quantum thermodynamics using a superconducting flux qubit (2025), arXiv:2411.10774 [quant-ph]

  24. [24]

    Le Hur, Kondo resonance of a microwave photon, Phys

    K. Le Hur, Kondo resonance of a microwave photon, Phys. Rev. B85, 140506 (2012)

  25. [25]

    Goldstein, M

    M. Goldstein, M. H. Devoret, M. Houzet, and L. I. Glaz- man, Inelastic Microwave Photon Scattering off a Quan- tum Impurity in a Josephson-Junction Array, Phys. Rev. Lett.110, 017002 (2013)

  26. [26]

    Saito and T

    K. Saito and T. Kato, Kondo signature in heat trans- fer via a local two-state system, Phys. Rev. Lett.111, 214301 (2013)

  27. [27]

    Yamamoto and T

    T. Yamamoto and T. Kato, Quantum critical phenomena in heat transport via a two-state system, Phys. Rev. B 98, 245412 (2018)

  28. [28]

    De Filippis, A

    G. De Filippis, A. de Candia, G. Di Bello, C. A. Perroni, L. M. Cangemi, A. Nocera, M. Sassetti, R. Fazio, and V. Cataudella, Signatures of dissipation driven quantum phase transition in rabi model, Phys. Rev. Lett.130, 210404 (2023)

  29. [29]

    Schmid, Diffusion and localization in a dissipative quantum system, Phys

    A. Schmid, Diffusion and localization in a dissipative quantum system, Phys. Rev. Lett.51, 1506 (1983)

  30. [30]

    Houzet, T

    M. Houzet, T. Yamamoto, and L. I. Glazman, Microwave spectroscopy of the schmid transition, Phys. Rev. B109, 155431 (2024)

  31. [31]

    Burshtein and M

    A. Burshtein and M. Goldstein, Inelastic decay from in- tegrability, PRX Quantum5, 020323 (2024)

  32. [32]

    Yamamoto, L

    T. Yamamoto, L. I. Glazman, and M. Houzet, Thermal transport across a josephson junction in a dissipative en- vironment, Phys. Rev. B110, L060512 (2024)

  33. [33]

    Paris, L

    N. Paris, L. Giacomelli, R. Daviet, C. Ciuti, N. Dupuis, and C. Mora, Resilience of the quantum critical line in the schmid transition, Phys. Rev. B111, 064509 (2025)

  34. [34]

    Subero, O

    D. Subero, O. Maillet, D. S. Golubev, G. Thomas, J. T. Peltonen, B. Karimi, M. Mar´ ın-Su´ arez, A. L. Yeyati, R. S´ anchez, S. Park, and J. P. Pekola, Bolometric de- tection of josephson inductance in a highly resistive en- vironment, Nat. Commun.14, 7924 (2023)

  35. [35]

    Nakamura, Y

    Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Coherent control of macroscopic quantum states in a single-cooper- pair box, Nature398, 786 (1999)

  36. [36]

    T. P. Orlando, J. E. Mooij, L. Tian, C. H. van der Wal, L. S. Levitov, S. Lloyd, and J. J. Mazo, Superconducting persistent-current qubit, Phys. Rev. B60, 15398 (1999)

  37. [37]

    J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, Rabi oscillations in a large josephson-junction qubit, Phys. Rev. Lett.89, 117901 (2002)

  38. [38]

    J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. I. Schuster, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, and R. J. Schoelkopf, Charge-insensitive qubit design de- rived from the cooper pair box, Phys. Rev. A76, 042319 (2007)

  39. [39]

    V. E. Manucharyan, J. Koch, L. I. Glazman, and M. H. Devoret, Fluxonium: Single cooper-pair circuit free of charge offsets, Science326, 113 (2009)

  40. [40]

    A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, and W. Zwerger, Dynamics of the dissi- pative two-state system, Rev. Mod. Phys.59, 1 (1987)

  41. [41]

    Weiss,Quantum Dissipative Systems, 4th ed

    U. Weiss,Quantum Dissipative Systems, 4th ed. (WORLD SCIENTIFIC, Singapore, 2012)

  42. [42]

    Magazz` u, P

    L. Magazz` u, P. Forn-D´ ıaz, R. Belyansky, J. L. Orgiazzi, M. A. Yurtalan, M. R. Otto, A. Lupascu, C. M. Wilson, and M. Grifoni, Probing the strongly driven spin-boson model in a superconducting quantum circuit, Nat. Com- mun.9, 1403 (2018)

  43. [43]

    Lepp¨ akangas, J

    J. Lepp¨ akangas, J. Braum¨ uller, M. Hauck, J.-M. Reiner, I. Schwenk, S. Zanker, L. Fritz, A. V. Ustinov, M. Weides, and M. Marthaler, Quantum simulation of the spin-boson model with a microwave circuit, Phys. Rev. A97, 052321 (2018)

  44. [44]

    Nitzan,Chemical Dynamics in Condensed Phases: Relaxation, Transfer and Reactions in Condensed Molec- ular Systems(Oxford University Press, 2006)

    A. Nitzan,Chemical Dynamics in Condensed Phases: Relaxation, Transfer and Reactions in Condensed Molec- ular Systems(Oxford University Press, 2006)

  45. [45]

    Recati, P

    A. Recati, P. O. Fedichev, W. Zwerger, J. von Delft, and P. Zoller, Atomic Quantum Dots Coupled to a Reservoir of a Superfluid Bose-Einstein Condensate, Phys. Rev. Lett.94, 040404 (2005)

  46. [46]

    Yamamoto, M

    T. Yamamoto, M. Kato, T. Kato, and K. Saito, Heat transport via a local two-state system near thermal equi- librium, New J. Phys.20, 093014 (2018)

  47. [47]

    K. A. Velizhanin, M. Thoss, and H. Wang, Meir–wingreen formula for heat transport in a spin-boson nanojunction model, J. Chem. Phys.133, 084503 (2010)

  48. [48]

    Ruokola and T

    T. Ruokola and T. Ojanen, Thermal conductance in a spin-boson model: Cotunneling and low-temperature properties, Phys. Rev. B83, 045417 (2011)

  49. [49]

    Nicolin and D

    L. Nicolin and D. Segal, Quantum fluctuation theorem for heat exchange in the strong coupling regime, Phys. Rev. B84, 161414 (2011)

  50. [50]

    Segal, Heat transfer in the spin-boson model: A com- parative study in the incoherent tunneling regime, Phys

    D. Segal, Heat transfer in the spin-boson model: A com- parative study in the incoherent tunneling regime, Phys. Rev. E90, 012148 (2014)

  51. [51]

    C. Wang, J. Ren, and J. Cao, Nonequilibrium energy transfer at nanoscale: A unified theory from weak to strong coupling, Sci. Rep.5, 11787 (2015)

  52. [52]

    Strathearn, P

    A. Strathearn, P. Kirton, D. Kilda, J. Keeling, and B. W. Lovett, Efficient non-markovian quantum dynamics using time-evolving matrix product operators, Nat. Commun. 9, 3322 (2018)

  53. [53]

    Feynman and F

    R. Feynman and F. Vernon, The theory of a general quan- tum system interacting with a linear dissipative system, Ann. Phys.24, 118 (1963)

  54. [54]

    Popovic, M

    M. Popovic, M. T. Mitchison, A. Strathearn, B. W. Lovett, J. Goold, and P. R. Eastham, Quantum heat statistics with time-evolving matrix product operators, PRX Quantum2, 020338 (2021)

  55. [55]

    A. J. Dunnett and A. W. Chin, Matrix product state simulations of non-equilibrium steady states and tran- sient heat flows in the two-bath spin-boson model at finite temperatures, Entropy23, 10.3390/e23010077 (2021). 15

  56. [56]

    Chen, Heat current in non-markovian open systems, New J

    R. Chen, Heat current in non-markovian open systems, New J. Phys.25, 033035 (2023)

  57. [57]

    Guinea, V

    F. Guinea, V. Hakim, and A. Muramatsu, Bosonization of a two-level system with dissipation, Phys. Rev. B32, 4410 (1985)

  58. [58]

    T. A. Costi, Scaling and universality in the anisotropic kondo model and the dissipative two-state system, Phys. Rev. Lett.80, 1038 (1998)

  59. [59]

    Lesage and H

    F. Lesage and H. Saleur, Perturbation of infra-red fixed points and duality in quantum impurity problems, Nucl. Phys. B546, 585 (1999)

  60. [60]

    Freton, E

    L. Freton, E. Boulat, and H. Saleur, Infrared expansion of entanglement entropy in the interacting resonant level model, Nucl. Phys. B874, 279 (2013)

  61. [61]

    A. C. Hewson,The Kondo Problem to Heavy Fermions (Cambridge University Press, 1993)

  62. [62]

    Lesage and H

    F. Lesage and H. Saleur, Strong-coupling resistivity in the kondo model, Phys. Rev. Lett.82, 4540 (1999)

  63. [63]

    Meir and N

    Y. Meir and N. S. Wingreen, Landauer formula for the current through an interacting electron region, Phys. Rev. Lett.68, 2512 (1992)

  64. [64]

    Ojanen and A.-P

    T. Ojanen and A.-P. Jauho, Mesoscopic photon heat transistor, Phys. Rev. Lett.100, 155902 (2008)

  65. [65]

    We neglected the temperature-dependent term inO 4 in- troduced in Ref. [50]

  66. [66]

    [50, 51] to ensure consistency with our numerical simulations

    We changed the sign ofc 4 relative to the original ex- pression in Refs. [50, 51] to ensure consistency with our numerical simulations

  67. [67]

    Makri and D

    N. Makri and D. E. Makarov, Tensor propagator for itera- tive quantum time evolution of reduced density matrices. I. Theory, J. Chem. Phys.102, 4600 (1995)

  68. [68]

    Makri and D

    N. Makri and D. E. Makarov, Tensor propagator for it- erative quantum time evolution of reduced density ma- trices. II. Numerical methodology, J. Chem. Phys.102, 4611 (1995)

  69. [69]

    Fishman, S

    M. Fishman, S. R. White, and E. M. Stoudenmire, The ITensor Software Library for Tensor Network Calcula- tions, SciPost Phys. Codebases , 4 (2022)