pith. sign in

arxiv: 2506.16122 · v1 · submitted 2025-06-19 · 🪐 quant-ph

Quadratic Hamiltonian approach to heat transport in fermionic systems

Ilari K. M\"akinen , Ivan M. Khaymovich , Jukka P. Pekola This is my paper

Pith reviewed 2026-05-19 09:31 UTC · model grok-4.3

classification 🪐 quant-ph
keywords quantum heat transportfermionic systemsquadratic HamiltonianPeschel tricksingle-mode heat valvenumerical methodnon-interacting fermionscorrelation functions
0
0 comments X

The pith

Quadratic Hamiltonians combined with the Peschel trick yield an efficient numerical route to heat transport in non-interacting fermionic systems.

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

The paper establishes a numerical technique for quantum heat transport that works by representing the system through a quadratic Hamiltonian and applying the Peschel trick to single-particle correlation functions. This reduces the problem of computing heat currents to operations on manageable matrices rather than the full many-body state. The method is demonstrated on a single-mode heat valve, where numerical outputs are checked against known analytical expressions for several choices of coupling strengths, temperatures, and geometries. A reader cares because the approach offers a scalable way to simulate thermal flows in simple quantum devices while preserving accuracy for the quadratic case.

Core claim

The quadratic Hamiltonian approach combined with the Peschel trick provides an efficient numerical method to describe heat transport in a single mode heat valve, with results that can be systematically compared to analytical formulae across configurations and parameters. The analysis confirms that the reduced correlation matrices capture the essential dynamics of heat flow between fermionic reservoirs.

What carries the argument

The Peschel trick applied to single-particle correlation functions extracted from quadratic fermionic Hamiltonians, which converts heat-transport observables into quantities computable from reduced correlation matrices.

If this is right

  • Numerical results for heat currents in the single-mode heat valve agree with analytical predictions over a range of couplings and temperatures.
  • The method permits direct, systematic comparison of different valve configurations without recomputing the full many-body evolution.
  • It reduces computational cost for heat transport studies while retaining accuracy for any quadratic fermionic system.
  • The same framework applies to other simple geometries once the single-particle spectrum and couplings are specified.

Where Pith is reading between the lines

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

  • If interactions remain absent, the approach could be applied to chains or networks of multiple heat valves to study collective thermal behavior.
  • Time-dependent driving of the quadratic parameters might be added while preserving the correlation-matrix structure for driven transport.
  • Direct comparison with mesoscopic experiments on superconducting or semiconductor quantum dots could test the numerical predictions beyond analytic limits.

Load-bearing premise

The fermionic system must remain quadratic and non-interacting so that single-particle correlations alone determine the heat transport without higher-order effects.

What would settle it

A clear mismatch between the numerical heat currents and the exact analytical formulas for the single-mode heat valve at fixed parameters would show the method fails to capture the transport.

Figures

Figures reproduced from arXiv: 2506.16122 by Ilari K. M\"akinen, Ivan M. Khaymovich, Jukka P. Pekola.

Figure 1
Figure 1. Figure 1: FIG. 1. The setup of the fermionic single mode heat valve. [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The estimated steady-state heat current through the [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) The total heat current into the cold bath [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The computed anomalous heat current (Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
read the original abstract

We investigate the problem of quantum heat transport, based on the quadratic fermionic systems with help of the Peschel trick of single-particle correlation functions. The efficient numerical method is applied to the particular case of a single mode heat valve and the results are compared to analytical formulae. Comparing several configurations and parameters we perform the systematic analysis of the method to most efficiently and accurately describe the simple quantum heat valve system.

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

2 major / 2 minor

Summary. The manuscript develops a numerical method for quantum heat transport in quadratic (non-interacting) fermionic systems by combining the quadratic Hamiltonian formalism with the Peschel trick, which reduces the dynamics to single-particle correlation matrices. The approach is applied to the specific case of a single-mode heat valve, with results systematically compared to independent analytical formulae across multiple configurations and parameter regimes to evaluate efficiency and accuracy.

Significance. If the numerical implementation is shown to reproduce the analytical benchmarks with controlled errors and to offer clear computational advantages, the work would supply a practical and extensible tool for studying heat currents in quadratic fermionic models, where the Peschel reduction is exact. The explicit restriction to the non-interacting regime and the use of external analytical validation avoid circularity and provide a solid foundation for further applications.

major comments (2)
  1. [Results] Results section (around the discussion of numerical vs. analytical comparisons): the manuscript states that results are compared to analytical formulae and that a systematic analysis is performed, yet no quantitative error metrics (e.g., relative deviations, L2 norms, or convergence plots) or tabulated agreement data are presented. This omission makes it impossible to judge how well the method performs or to substantiate the efficiency claim.
  2. [Method] Method section describing the Peschel reduction and heat-current operator: while the quadratic premise is correctly stated as the domain of applicability, the text does not explicitly show how the heat current expectation value is computed from the single-particle correlation matrix (e.g., via the appropriate trace formula). Adding this short derivation would strengthen the link between the formalism and the reported observables.
minor comments (2)
  1. [Abstract and Results] The abstract mentions 'several configurations and parameters' but the main text would benefit from a clear table or figure caption that enumerates the exact parameter sets used in the comparisons.
  2. [Method] Notation for the single-particle correlation matrix and the Peschel transformation should be introduced once with a consistent symbol and then used uniformly; occasional redefinitions slow readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive summary and constructive major comments. The suggestions have helped us strengthen the manuscript by adding quantitative validation and clarifying the connection between the formalism and observables. We address each point below.

read point-by-point responses

  1. Referee: [Results] Results section (around the discussion of numerical vs. analytical comparisons): the manuscript states that results are compared to analytical formulae and that a systematic analysis is performed, yet no quantitative error metrics (e.g., relative deviations, L2 norms, or convergence plots) or tabulated agreement data are presented. This omission makes it impossible to judge how well the method performs or to substantiate the efficiency claim.

    Authors: We agree that quantitative error metrics are necessary to rigorously assess performance and support the efficiency claim. In the revised manuscript we have added a new subsection in Results containing relative deviation plots, L2-norm errors, and a table of maximum/average deviations across all tested configurations and parameter regimes. These additions directly address the concern and provide the requested substantiation. revision: yes

  2. Referee: [Method] Method section describing the Peschel reduction and heat-current operator: while the quadratic premise is correctly stated as the domain of applicability, the text does not explicitly show how the heat current expectation value is computed from the single-particle correlation matrix (e.g., via the appropriate trace formula). Adding this short derivation would strengthen the link between the formalism and the reported observables.

    Authors: We thank the referee for this observation. We have inserted a short, self-contained derivation in the Methods section that shows how the heat-current expectation value is obtained from the single-particle correlation matrix via the trace formula involving the appropriate current operator. This explicit step clarifies the implementation without altering the quadratic premise or the Peschel reduction. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained within stated quadratic domain

full rationale

The paper restricts its scope to quadratic (non-interacting) fermionic Hamiltonians and applies the Peschel reduction to single-particle correlation matrices, which is an exact rewriting for this class of models. The central claim concerns numerical efficiency for the single-mode heat valve together with direct comparisons against independent analytical expressions across parameter choices. No internal step reduces by construction to a fitted input or self-citation chain; the quadratic premise is the explicitly stated domain of applicability, and the reported comparisons serve as external validation rather than tautological confirmation.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no explicit free parameters, axioms, or invented entities; the approach implicitly assumes quadratic (non-interacting) fermionic Hamiltonians allow exact reduction to single-particle correlations.

pith-pipeline@v0.9.0 · 5592 in / 1083 out tokens · 29024 ms · 2026-05-19T09:31:07.883078+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    A.I.Khinchin, Mathematical foundations of statistical mechanics (Courier Corporation, 1949)

    work page 1949

  2. [2]

    J. M. Deutsch, Quantum statistical mechanics in a closed system, Phys. Rev. A43, 2046 (1991)

    work page 2046

  3. [3]

    Srednicki, Chaos and quantum thermalization, Phys

    M. Srednicki, Chaos and quantum thermalization, Phys. Rev. E 50, 888 (1994)

    work page 1994

  4. [4]

    Srednicki, Thermal fluctuations in quantized chaotic systems, Journal of Physics A: Mathematical and Gen- eral 29, L75 (1996)

    M. Srednicki, Thermal fluctuations in quantized chaotic systems, Journal of Physics A: Mathematical and Gen- eral 29, L75 (1996)

    work page 1996

  5. [5]

    D. N. Page, Average entropy of a subsystem, Phys. Rev. Lett. 71, 1291 (1993)

    work page 1993

  6. [6]

    Łydżba, Y

    P. Łydżba, Y. Zhang, M. Rigol, and L. Vidmar, Single- particle eigenstate thermalization in quantum-chaotic quadratic hamiltonians, Phys. Rev. B 104, 214203 (2021)

    work page 2021

  7. [7]

    Łydżba, M

    P. Łydżba, M. Mierzejewski, M. Rigol, and L. Vid- mar, Generalized thermalization in quantum-chaotic quadratic hamiltonians, Phys. Rev. Lett. 131, 060401 (2023)

    work page 2023

  8. [8]

    Łydżba, R

    P. Łydżba, R. Świętek, M. Mierzejewski, M. Rigol, and L. Vidmar, Normal weak eigenstate thermalization (2024), arXiv:2404.02199 [cond-mat.stat-mech]

    work page arXiv 2024

  9. [9]

    Vidmar, L

    L. Vidmar, L. Hackl, E. Bianchi, and M. Rigol, Entangle- ment entropy of eigenstates of quadratic fermionic hamil- tonians, Phys. Rev. Lett.119, 020601 (2017)

    work page 2017

  10. [10]

    Vidmar and M

    L. Vidmar and M. Rigol, Entanglement entropy of eigen- states of quantum chaotic hamiltonians, Phys. Rev. Lett. 119, 220603 (2017)

    work page 2017

  11. [11]

    Hackl, L

    L. Hackl, L. Vidmar, M. Rigol, and E. Bianchi, Aver- age eigenstate entanglement entropy of the xy chain in a transverse field and its universality for translationally invariant quadratic fermionic models, Phys. Rev. B99, 075123 (2019)

    work page 2019

  12. [12]

    P.Łydżba, M.Rigol,andL.Vidmar,Eigenstateentangle- ment entropy in random quadratic hamiltonians, Phys. Rev. Lett. 125, 180604 (2020)

    work page 2020

  13. [13]

    Łydżba, M

    P. Łydżba, M. Rigol, and L. Vidmar, Entanglement in many-body eigenstates of quantum-chaotic quadratic hamiltonians, Phys. Rev. B103, 104206 (2021)

    work page 2021

  14. [14]

    Zanardi, Quantum entanglement in fermionic lattices, Phys

    P. Zanardi, Quantum entanglement in fermionic lattices, Phys. Rev. A65, 042101 (2002)

    work page 2002

  15. [15]

    Peschel, Calculation of reduced density matrices from correlation functions, Journal of Physics A: Mathemati- cal and General36, L205 (2003)

    I. Peschel, Calculation of reduced density matrices from correlation functions, Journal of Physics A: Mathemati- cal and General36, L205 (2003)

    work page 2003

  16. [16]

    Calabrese, M

    P. Calabrese, M. Mintchev, and E. Vicari, Entanglement entropy of one-dimensional gases, Phys. Rev. Lett.107, 020601 (2011)

    work page 2011

  17. [17]

    Peschel and V

    I. Peschel and V. Eisler, Reduced density matrices and entanglement entropy in free lattice models, Journal of Physics A: Mathematical and Theoretical 42, 504003 (2009)

    work page 2009

  18. [18]

    B. C. Dias, M. Haque, P. Ribeiro, and P. McClarty, Diffu- sive operator spreading for random unitary free fermion circuits (2021)

    work page 2021

  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, Nature Physics 14, 991–995 (2018)

    work page 2018

  20. [20]

    Colpa, Diagonalization of the quadratic boson hamil- tonian, Physica A: Statistical Mechanics and its Appli- cations 93, 327–353 (1978)

    J. Colpa, Diagonalization of the quadratic boson hamil- tonian, Physica A: Statistical Mechanics and its Appli- cations 93, 327–353 (1978). 7

    work page 1978

  21. [21]

    Muller, S

    D. Segal and A. Nitzan, Spin-boson thermal rec- tifier, Physical Review Letters 94, 10.1103/phys- revlett.94.034301 (2005)

    work page doi:10.1103/phys- 2005

  22. [22]

    Senior, A

    J. Senior, A. Gubaydullin, B. Karimi, J. T. Peltonen, J.Ankerhold,andJ.P.Pekola,Heatrectificationviaasu- perconducting artificial atom, Communications Physics 3, 10.1038/s42005-020-0307-5 (2020)

    work page doi:10.1038/s42005-020-0307-5 2020

  23. [23]

    Ojanen and A.-P

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

    work page 2008

  24. [24]

    C.D.Satrya, A.Guthrie, I.K.Mäkinen,andJ.P.Pekola, Electromagnetic simulation and microwave circuit ap- proachofheattransportinsuperconductingqubits,Jour- nal of Physics Communications7, 015005 (2023)

    work page 2023

  25. [25]

    B. K. Agarwalla and D. Segal, Energy current and its statistics in the nonequilibrium spin-boson model: Ma- jorana fermion representation, New Journal of Physics 19, 043030 (2017)

    work page 2017

  26. [26]

    G. E. Topp, T. Brandes, and G. Schaller, Steady-state thermodynamics of non-interacting transport beyond weak coupling, Europhysics Letters110, 67003 (2015). Appendix A: Derivation of the anomalous heat current This appendix elaborates the concept of dividing the total heat current into a normal and an anomalous part, which is used in the main text to exp...

    work page 2015