Bedingham and Maroney (New J. Phys. 18, 113050, 2016) proved that implementing a quantum operation — a completely positive trace-preserving map that must act correctly on every input state of an ensemble, not merely on average — carries an operation-specific excess heat ϵQ\epsilon_Q beyond the Landauer term kTln⁡2[S(ρS)−S(ρS′)]kT\ln 2[S(\rho_S)-S(\rho'_S)], whenever the output states do not share a common diagonal basis whose eigenvalues are related to the input diagonal statistics by a left stochastic map, barring exceptional spectral symmetries. Their argument yields nonzero lower bounds on ϵQ\epsilon_Q via trace-norm and spectral-ratio inequalities, but no matching construction, and they explicitly state that demonstrating a protocol which reaches ϵQ\epsilon_Q is an open problem. The task is to determine the exact characterization of the minimal excess heat: whether the derived bounds are tight, which physical protocols attain the minimum for general quantum operations, and which operations possessing eigenvalue-ratio coincidences λi/λj=λk′/λl′\lambda_i/\lambda_j=\lambda'_k/\lambda'_l admit dissipation-free (ϵQ=0\epsilon_Q=0) implementations. The problem remains open: the most direct continuation (Dunlop, Cerisola, Tabanera-Bravo and Anders, Phys. Rev. Res. 6, 033005, 2024) constructs optimisable protocols only for two-input single-qubit operations, states that before their work no prescription for optimal implementation and no upper bound on the energy cost existed, and notes that their own analytic bound is not tight; the general tightness, attainability, and exceptional-case questions remain unresolved.

Background

Landauer's principle states that erasing one bit of information while in contact with an environment at temperature TT dissipates at least kTln⁡2kT\ln 2 of heat, where kk is Boltzmann's constant; reaching this bound is equivalent to performing the operation in a thermodynamically reversible way. In quantum theory the entropy entering the bound is the von Neumann entropy S(ρ)=−Tr[ρlog⁡2ρ]S(\rho)=-\mathrm{Tr}[\rho\log_2\rho] of the density matrix ρ\rho, so an operation taking an average input state ρS\rho_S to an average output state ρS′\rho'_S should satisfy ΔEQ≥kTln⁡2[S(ρS)−S(ρS′)]\Delta E_Q \ge kT\ln 2[S(\rho_S)-S(\rho'_S)], where ΔEQ\Delta E_Q is the heat dumped into the bath. A quantum operation - formally a completely positive trace-preserving map on density matrices, the most general evolution of an open system's state, realizable as unitary interaction with an environment whose part is then discarded - must map each possible input state ρSn\rho^n_S of a known ensemble to its designated output ρS′n\rho'^n_S, not just the average state, and this per-input requirement is what generates the excess cost studied here. Bedingham and Maroney (2016) analyzed an implementation model in which the system SS is coupled to an auxiliary system AA and a heat bath BB, the auxiliary starts in a standard catalyst state ρA\rho_A that the operation must restore on average and the heat bath starts in a canonical thermal state ρB\rho_B at temperature TT, the joint evolution is a unitary, and ΔEQ\Delta E_Q is the resulting energy change of the bath. They proved that for every quantum operation QQ there is a non-negative, implementation-independent quantity ϵQ\epsilon_Q with ΔEQ≥kTln⁡2[S(ρS)−S(ρS′)]+ϵQ\Delta E_Q \ge kT\ln 2[S(\rho_S)-S(\rho'_S)]+\epsilon_Q. Writing the spectral decompositions of the average input and output states as ρS=∑iλi∣ϕi⟩⟨ϕi∣\rho_S=\sum_i\lambda_i|\phi_i\rangle\langle\phi_i| and ρS′=∑kλk′∣ϕk′⟩⟨ϕk′∣\rho'_S=\sum_k\lambda'_k|\phi'_k\rangle\langle\phi'_k|, and the individual inputs and outputs in these bases as ρSn=∑ijμijn∣ϕi⟩⟨ϕj∣\rho^n_S=\sum_{ij}\mu^n_{ij}|\phi_i\rangle\langle\phi_j| and ρS′n=∑klμkl′n∣ϕk′⟩⟨ϕl′∣\rho'^n_S=\sum_{kl}\mu'^n_{kl}|\phi'_k\rangle\langle\phi'_l|, they showed that ϵQ=0\epsilon_Q=0 is possible whenever all output states share a common diagonal basis and the output eigenvalues are related to the input diagonal elements by a common left stochastic map PQ(k∣i)P_Q(k|i) (non-negative coefficients with ∑kPQ(k∣i)=1\sum_k P_Q(k|i)=1 and μkl′n=δkl∑iPQ(k∣i)μiin\mu'^n_{kl}=\delta_{kl}\sum_i P_Q(k|i)\mu^n_{ii}); classical information processing, erasure to a standard state, pure unitary operations, and single-input operations fall in this class, and explicit reversible constructions were given. Otherwise, barring exceptional spectral symmetries — coincidences of eigenvalue ratios λi/λj=λk′/λl′\lambda_i/\lambda_j=\lambda'_k/\lambda'_l for some i≠ji\neq j or k≠lk\neq l — every physical implementation of the operation has ϵQ>0\epsilon_Q>0. The proof rewrites any implementation through complex coefficients q(kl∣ij)q(kl|ij) (with μkl′n=∑ijq(kl∣ij)μijn\mu'^n_{kl}=\sum_{ij}q(kl|ij)\mu^n_{ij}) and yields nonzero lower bounds on ϵQ\epsilon_Q from a trace-norm inequality and from spectral-ratio inequalities of the form ∥ρ′−ρ⋆∥1≥∣λjλk′−λiλl′∣/(λk′+λl′)⋅∣q(kl∣ij)∣\|\rho'-\rho^\star\|_1 \ge |\lambda_j\lambda'_k-\lambda_i\lambda'_l|/(\lambda'_k+\lambda'_l)\cdot|q(kl|ij)|, where ρ′\rho' is the joint final state of system, auxiliary, and bath actually produced by the implementation unitary, ρ⋆=ρS′⊗ρA⊗ρB\rho^\star=\rho'_S\otimes\rho_A\otimes\rho_B is the target decoupled final state (system in the correct average output, auxiliary and bath uncorrelated and back in their initial states ρA\rho_A and ρB\rho_B), and ∥⋅∥1\|\cdot\|_1 is the trace norm (sum of singular values of its argument). The authors explicitly caution that these bounds may not be tight — their results only place a nonzero lower bound on the excess energy cost — and that the mere existence of spectral-ratio coincidences is not sufficient for ϵQ=0\epsilon_Q=0, since an implementation must in addition have q(kl∣ij)=0q(kl|ij)=0 for every index quadruple with λi/λj≠λk′/λl′\lambda_i/\lambda_j\neq\lambda'_k/\lambda'_l. Later work has circumscribed but not resolved the question: Dunlop, Cerisola, Tabanera-Bravo and Anders (2024) developed explicit optimisable protocols, feasibility conditions, and a complementary (explicitly non-tight) bound on dissipated work for the two-input single-qubit case, noting that prior to their work no prescription for thermodynamically optimal multipurpose operations and no upper bound on their energy cost existed; Chiribella, Meng, Renner and Yung (2022) bounded the nonequilibrium resources needed for approximately accurate information processing in a resource-theoretic framework, a complementary setting that does not determine ϵQ\epsilon_Q for exact operations; experimental energy accounting for a driven single-qubit gate (Stevens et al., 2022) and lower bounds on the field energy required by Hamiltonian gates (Stevens and Deffner, 2025) use different resource accounting and leave ϵQ\epsilon_Q untouched.

Problem Statement

In the Bedingham–Maroney implementation model — a system SS coupled to an auxiliary system AA and a heat bath BB, with AA initially in a standard catalyst state ρA\rho_A restored on average and BB initially in a canonical thermal state at temperature TT, joint unitary evolution, and ΔEQ\Delta E_Q the resulting energy change of the bath — determine the minimal excess thermodynamic cost of an arbitrary quantum operation QQ acting on an ensemble of input states {ρSn}\{\rho^n_S\} with probabilities pnp_n and designated outputs ρS′n=Q(ρSn)\rho'^n_S = Q(\rho^n_S). Specifically: (1) For operations failing the reversible condition of Theorem 1 (the outputs do not share a common diagonal basis related to the input diagonal statistics by a common left stochastic map PQ(k∣i)P_Q(k|i)), determine the exact minimum of the excess heat ϵ=ΔEQ−kTln⁡2[S(ρS)−S(ρS′)]\epsilon = \Delta E_Q - kT\ln 2[S(\rho_S)-S(\rho'_S)] over all physical implementations of QQ (the minimum across implementations being the implementation-independent operation-level quantity ϵQ\epsilon_Q of the Background); decide whether the lower bounds derived from the trace-norm inequality ϵQ≥12kT∥ρ′−ρ⋆∥12\epsilon_Q \ge \frac{1}{2}kT\|\rho'-\rho^\star\|_1^2 (with ρ′\rho' the actual joint final state of system, auxiliary, and bath produced by the implementation and ρ⋆=ρS′⊗ρA⊗ρB\rho^{\star}=\rho'_S\otimes\rho_A\otimes\rho_B the decoupled target state) and the spectral-ratio inequalities ∥ρ′−ρ⋆∥1≥∣λjλk′−λiλl′∣/(λk′+λl′)⋅∣q(kl∣ij)∣\|\rho'-\rho^{\star}\|_1 \ge |\lambda_j\lambda'_k-\lambda_i\lambda'_l|/(\lambda'_k+\lambda'_l)\cdot|q(kl|ij)| are tight; and determine whether the minimal excess heat is attained by a physical protocol, constructing attaining protocols wherever they exist. (2) Classify the exceptional spectral-symmetry cases: characterize, in terms of the spectral data {λi}\{\lambda_i\}, {λk′}\{\lambda'_k\}, {μijn}\{\mu^n_{ij}\}, {μkl′n}\{\mu'^n_{kl}\} of the operation, which operations possessing eigenvalue-ratio coincidences λi/λj=λk′/λl′\lambda_i/\lambda_j=\lambda'_k/\lambda'_l (with i≠ji\neq j or k≠lk\neq l) admit implementations with ϵQ=0\epsilon_Q=0, and construct such dissipation-free implementations, noting that the existence of the coincidences alone is not sufficient. The two parts are jointly required and together characterize the minimal excess heat ϵQ\epsilon_Q of an arbitrary quantum operation of arbitrary system dimension acting on arbitrary input ensembles (not only qubits or two-input cases), and the exact boundary between thermodynamically reversible and irreversible quantum operations.

References

  1. D. J. Bedingham and O. J. E. Maroney, "The thermodynamic cost of quantum operations", New Journal of Physics 18, 113050 (2016). DOI: 10.1088/1367-2630/18/11/113050; arXiv:1604.03749. https://doi.org/10.1088/1367-2630/18/11/113050
  2. Joe Dunlop, Federico Cerisola, Jorge Tabanera-Bravo, and Janet Anders, "Thermodynamically optimal protocols for dual-purpose qubit operations", Physical Review Research 6, 033005 (2024). DOI: 10.1103/PhysRevResearch.6.033005; arXiv:2306.09088. https://doi.org/10.1103/PhysRevResearch.6.033005
  3. Giulio Chiribella, Fei Meng, Renato Renner, and Man-Hong Yung, "The nonequilibrium cost of accurate information processing", Nature Communications 13, 7155 (2022). DOI: 10.1038/s41467-022-34541-w. https://doi.org/10.1038/s41467-022-34541-w
  4. J. Stevens, D. Szombati, M. Maffei, C. Elouard, R. Assouly, N. Cottet, R. Dassonneville, Q. Ficheux, S. Zeppetzauer, A. Bienfait, A. N. Jordan, A. Auffèves, and B. Huard, "Energetics of a Single Qubit Gate", Physical Review Letters 129, 110601 (2022). DOI: 10.1103/PhysRevLett.129.110601; arXiv:2109.09648. https://doi.org/10.1103/PhysRevLett.129.110601
  5. Josey Stevens and Sebastian Deffner, "Hamiltonian quantum gates — energetic advantage from entangleability", Quantum Science and Technology 10, 04LT03 (2025). DOI: 10.1088/2058-9565/ae0daf; arXiv:2507.01758. https://doi.org/10.1088/2058-9565/ae0daf
  6. M. H. Mohammady, A. Auffèves, and J. Anders, "Energetic footprints of irreversibility in the quantum regime", Communications Physics 3, 89 (2020). DOI: 10.1038/s42005-020-0356-9. https://doi.org/10.1038/s42005-020-0356-9
  7. Kais Abdelkhalek, Yoshifumi Nakata, and David Reeb, "Fundamental energy cost for quantum measurement", arXiv:1609.06981 (2016). https://arxiv.org/abs/1609.06981
  8. Ji-Yong Sung and Jae-Ho Cheong, "Quantum thermodynamics with information and coherence", Quantum Information Processing 25, 200 (2026). DOI: 10.1007/s11128-026-05206-0. https://doi.org/10.1007/s11128-026-05206-0