<<

Quantum Transport and Thermodynamics

Giuliano Benenti Center for Nonlinear and Complex Systems, Univ. Insubria, Como, Italy INFN, Milano, Italy Outline 1) Basic thermodynamics of nonequilibrium states (linear response, Onsager relations, efficiency of thermal machines, finite- thermodynamics) 2) Landauer formalism (scattering theory) (energy filtering, scattering theory and the laws of thermodynamics) 3) Rate equations (local detailed balance, examples of thermal machines) 4) Thermodynamic bounds on -to-work conversion (power-efficiency trade-off, thermodynamic uncertainty relations and power-efficiency-fluctuations trade-off) What is quantum in energy conversion? Ex: Traditional versus quantum thermoelectrics

Structures smaller than the Relaxation length (tens of relaxation length (many nanometers at room microns at low ); temperature) of the order of quantum interference effects; the mean free path; inelastic Boltzmann transport theory s c a t t e r i n g ( p h o n o n s ) cannot be applied; efficiency thermalizes the electrons depends on geometry and size [see G. B., G. Casati, K. Saito, R. S. Whitney, Phys. Rep. 694, 1 (2017)] Traditional thermocouple

Quantum thermocouple Ex: Cooling by heating A third (hot) reservoir can help cooling [Pekola & Hekking, PRL 98, 210604 (2007); Mari & Eisert, PRL 108, 120602 (2012), …]

[Cleuren, Rutten, Van den Broeck, PRL 108, 120603 (2012)] I. Basic thermodynamics of nonequilibrium states The Nobel Prize in Chemistry 1968: From the award ceremony speech “Professor Lars Onsager has been awarded this year’s Nobel Prize for Chemistry (1968) for the discovery of the reciprocal relations, named after him, and basic to irreversible thermodynamics… Onsager’s reciprocal relations can be described as a universal natural law…It can be said that Onsager’s reciprocal relations represent a further law making possible a thermodynamic study of irreversible processes…It represents one of the great advances in science during this century. According to Nico Van Kampen Onsager derived his reciprocal relations in a “stroke of genius” Irreversible thermodynamic

Irreversible thermodynamics based on the postulates of equilibrium thermostatics plus the postulate of time- reversal symmetry of physical laws (if time t is replaced by -t and simultaneously applied magnetic field B by -B)

The thermodynamic theory of irreversible processes is based on the Onsager Reciprocity Theorem

Refs.: Thermodynamic and Irreversible processes are driven by thermodynamic forces (or generalized forces or affinities) Xi Fk Fluxes Ji characterize the response of the system to the applied forces production rate given by the sum of the products of each with its associated thermodynamic

S = S(U, V, N1,N2, ...)=S(E0,E1,E2, ...)

dS S dEk = = XkkJJk k dt Ek dt k k F k Linear response Purely resistive systems: fluxes at a given instant depend only on the thermodynamic forces at that instant (memory effects not considered) J = L + L + ... i ijFj ijkFjFk j j,k Fluxes vanish as thermodynamic forces vanish Linear (and purely resistive) processes: J = L i ijFj j Lij Onsager coefficients (first-order kinetic coefficients) depend on intensive quantities (T,P,µ,...) Phenomenological linear Ohm’s, Fourier’s, Fick’s laws Onsager reciprocal relations

Relationship of Onsager theorem to time-reversal symmetry of physical laws Consider delayed correlation moments of fluctuations (without applied magnetic fields)

E (t) E (t) E , E =0, j j j j E (t)E (t + ) = E (t)E (t ) = E (t + )E (t) j k j k j k

Ek(t + ) Ek(t) Ej(t + ) Ej(t) lim Ej(t) = lim Ek(t) 0 0 E E˙ = E˙ E j k j k Assume that fluctuations decay is governed by the same linear dynamical laws as are macroscopic processes E˙ = L k kl Fl l L E = L E kl j Fl jl Fl k

(null) j l Xl Assume that the fluctuation of each thermodynamic force is associated only with the fluctuation of the corresponding extensive variable

E = k h(null) j Fli B jl Onsager relations: L = L (null) jk kj Onsager-Casimir relations

Onsager reciprocal relations reflect at the macroscopic level the time-reversal symmetry of the microscopic dynamics, invariant under the transformation:

With an applied magnetic field one instead obtains Onsager-Casimir relations:

but in principle one could violate the Onsager symmetry: Linear response for coupled (particle and heat) flows

Stochastic baths: ideal gases at fixed temperature and electrochemical potential

µ = µ µ L R Onsager-Casimir relations: T = T T L R (we assume TL >TR,µL <µR)

Onsager reciprocal relations for time-reversal symmetric systems: Positivity of entropy production Entropy production rate given by the sum of the products of each flux with its associated thermodynamic force

S = S(U, V, N1,N2, ...)=S(E0,E1,E2, ...)

dS S dEk = = XkkJJk k dt Ek dt k k F k For thermoelectricity: Linear response:

Positivity of entropy production rate: Onsager and transport coefficients

Note that the positivity of entropy production implies that the (isothermal) electric conductance G>0 and the thermal conductance K>0 Seebeck and Peltier coefficients

Seebeck and Peltier coefficients are related by a Onsager-Casimir reciprocal relation (when time symmetry is not broken, we simply have ) Interpretation of the Peltier coefficient

Entropy current:

entropy transported by the electron flow each electron carries an entropy of advective term in thermal transport (reversible) open-circuit term in thermal transport (by electrons and phonons, irreversible) Entropy production/ heat rate

Joule heating

heat lost by To minimize disappears for time-reversal dissipation large G and symmetric systems small K are needed Linear response?

(exhaust gases)

(room temperature)

[Vining, Nat. Mater. 8, 83 (2009)] Linear response for small temperature and electrochemical potential differences (compared to the average temperature) on the scale of the relaxation length Exhaust pipe: temperature drop over a mm scale: temperature drop of 0.003 K on the relaxation length scale (of 10 nm) Maximum efficiency Within linear response and for steady-state heat to work conversion:

Find the maximum of η over , for fixed (i.e., over the applied voltage ΔV for fixed temperature difference ΔT)

(T T T ) (null) L ⇡ R ⇡ Thermoelectric figure of merit L2 GS2 ZT eh = T det L K Conditions for Carnot efficiency

ZT diverging implies that the Onsager matrix is ill- conditioned, that is, the condition number diverges:

In such case the system is singular (tight coupling limit): J J h e (the ratio Jh/Je is independent of the applied voltage and temperature gradients) Efficiency at maximum power

Output power

Find the maximum of P over , for fixed (over the applied voltage ΔV for fixed ΔT)

Maximum output power 2 T Leh 2 1 2 2 Pmax = h = S G(T ) 4 Lee F 4 Power factor P quadratic function of , with maximum at half of the stopping force:

Efficiency at maximum power ZT ( )= C C max 2 ZT +2 CA 2

ηCA Curzon-Ahlborn upper bound max

(Pmax) Efficiency versus power

⇒ Maximum refrigeration efficiency

Cooling power (heat extracted from the cold reservoir) Coefficient of performance (COP) J (r) = h P

(can be >1) ZT is the figure of merit also for refrigeration ZT is an intrinsic material property? For mesoscopic systems size-dependence for G,K,S can be expected In the diffusive transport regime Ohm’s and Fourier’s scaling laws hold: Local equilibrium Under the assumption of local equilibrium we can write phenomenological equations with ∇T and ∇µ rather than ΔT and Δµ

charge and heat current densities In this case we connect Onsager coefficients to electric and rather than to conductances j j = e , = h V T =0 T j =0 e II. Landauer formalism (scattering theory) Scattering theory

Scattering region connected to N terminals (reservoirs)

Describes elastic scattering (including the effect of a disorder potential), but not electron-electron interactions beyond Hartree approximation and electron-phonon interactions Transmission matrix Probability for an electron with energy E to go from (transverse) mode m of reservoir j to mode n of reservoir i: scattering matrix elements

transmission matrix elements

probabilities Conservation of current and condition of zero current at zero bias from:

From time reversal symmetry of the scatterer Hamiltonian: Landauer approach Electrical current into the scatterer from reservoir i:

Fermi function Energy current into the scatterer from reservoir i:

heat carried by an electron leaving reservoir i Heat current: Kirkhoff’s law of current conservation for (steady state) electrical and energy currents:

Heat current not conserved:

Heat dissipated in the reservoirs: entropy production rate

Heat (not energy) current gauge invariant. The generated power equals heat and so is also gauge invariant. Two-terminal (thermoelectric) power production

Left (L) Right (R) reservoir S reservoir TL , L TR , R P =[(µ µ )/e]J AAAB/nicbVDLSsNAFJ3UV62vqODGTbAIdWFNRFAXQtGNiIsqxhbSECbT23bo5MHMRCixC3/FjQsVt36HO//GSZuFth64l8M59zJ3jh8zKqRpfmuFmdm5+YXiYmlpeWV1TV/fuBdRwgnYJGIRb/pYAKMh2JJKBs2YAw58Bg2/f5H5jQfggkbhnRzE4Aa4G9IOJVgqydO36mdOpRUk3u1+1q/3DsC98sDTy2bVHMGYJlZOyihH3dO/Wu2IJAGEkjAshGOZsXRTzCUlDIalViIgxqSPu+AoGuIAhJuO7h8au0ppG52IqwqlMVJ/b6Q4EGIQ+GoywLInJr1M/M9zEtk5cVMaxomEkIwf6iTMkJGRhWG0KQci2UARTDhVtxqkhzkmUkVWUiFYk1+eJvZh9bRq3hyVa+d5GkW0jXZQBVnoGNXQJaojGxH0iJ7RK3rTnrQX7V37GI8WtHxnE/2B9vkDXKKUlw== R L e (T >T ,µ <µ ) AAAB/3icbVDLSsNAFJ3UV62vqAsXboJFqFBKIoIKIkU3LrqopbGFJoTJdNIOnUzCzEQooRt/xY0LFbf+hjv/xkmbhbYeuJfDOfcyc48fUyKkaX5rhaXlldW14nppY3Nre0ff3XsQUcIRtlFEI971ocCUMGxLIinuxhzD0Ke4449uM7/ziLkgEWvLcYzdEA4YCQiCUkmeflBpe43rtteqOlUnTLzGVdZbJ55eNmvmFMYisXJSBjmanv7l9COUhJhJRKEQPcuMpZtCLgmieFJyEoFjiEZwgHuKMhhi4abTAybGsVL6RhBxVUwaU/X3RgpDIcahryZDKIdi3svE/7xeIoMLNyUsTiRmaPZQkFBDRkaWhtEnHCNJx4pAxIn6q4GGkEMkVWYlFYI1f/IisU9rlzXz/qxcv8nTKIJDcAQqwALnoA7uQBPYAIEJeAav4E170l60d+1jNlrQ8p198Afa5w+yEpS/sha1_base64="C39OhB+IczRcjLNINXH29e9lt8M=">AAAB2HicbZDNSgMxFIXv1L86Vq1rN8EiuCpTN+pOcOOygmML7VAymTttaCYzJHeEMvQFXLhRfDB3vo3pz0KtBwIf5yTk3hMXSloKgi+vtrW9s7tX3/cPGv7h0XGz8WTz0ggMRa5y04+5RSU1hiRJYb8wyLNYYS+e3i3y3jMaK3P9SLMCo4yPtUyl4OSs7qjZCtrBUmwTOmtowVqj5ucwyUWZoSahuLWDTlBQVHFDUiic+8PSYsHFlI9x4FDzDG1ULcecs3PnJCzNjTua2NL9+aLimbWzLHY3M04T+zdbmP9lg5LS66iSuigJtVh9lJaKUc4WO7NEGhSkZg64MNLNysSEGy7INeO7Djp/N96E8LJ90w4eAqjDKZzBBXTgCm7hHroQgoAEXuDNm3iv3vuqqpq37uwEfsn7+Aap5IoMsha1_base64="1MiH6zav+sah2+m9pJXGXl8SvHs=">AAAB9HicbZBPS8MwGMbfzn9zTq1ePHgpDmHCGK0XFUQELx52mGNzg7WUNEu3sKQtSSqMsotfxYsHFT+LN7+N6baDbj6Q8ON5EpL3CRJGpbLtb6Owtr6xuVXcLu2Ud/f2zYPyo4xTgUkHxywWvQBJwmhEOooqRnqJIIgHjHSD8V2ed5+IkDSO2mqSEI+jYURDipHSlm8eVdt+46btt2puzeWp37jO99aZb1bsuj2TtQrOAiqwUNM3v9xBjFNOIoUZkrLv2InyMiQUxYxMS24qSYLwGA1JX2OEOJFeNhtgap1qZ2CFsdArUtbM/X0jQ1zKCQ/0SY7USC5nuflf1k9VeOllNEpSRSI8fyhMmaViK2/DGlBBsGITDQgLqv9q4RESCCvdWUmX4CyPvAqd8/pV3X6woQjHcAJVcOACbuEemtABDFN4gTd4N56NV+Nj3lbBWNR2CH9kfP4ATGOTVQ==sha1_base64="zDW22uEPn9Prthn+nPZd4cAl3CQ=">AAAB/3icbVDLSsNAFJ3UV62vqAsXbgaLUKGUxI0KIkU3LrqopbGFJoTJdNIOnUnCzEQopRt/xY0LFbf+hjv/xkmbhVYP3MvhnHuZuSdIGJXKsr6MwtLyyupacb20sbm1vWPu7t3LOBWYODhmsegGSBJGI+IoqhjpJoIgHjDSCUY3md95IELSOGqrcUI8jgYRDSlGSku+eVBp+42rtt+qulWXp37jMuutE98sWzVrBviX2DkpgxxN3/x0+zFOOYkUZkjKnm0lypsgoShmZFpyU0kShEdoQHqaRogT6U1mB0zhsVb6MIyFrkjBmfpzY4K4lGMe6EmO1FAuepn4n9dLVXjuTWiUpIpEeP5QmDKoYpilAftUEKzYWBOEBdV/hXiIBMJKZ1bSIdiLJ/8lzmntombdWeX6dZ5GERyCI1ABNjgDdXALmsABGEzBE3gBr8aj8Wy8Ge/z0YKR7+yDXzA+vgGw0pS7 L R L R The upper bound to efficiency is given by the Carnot efficiency (expected only at zero power; intuitively, finite currents entail dissipation): TR ⌘C =1

AAACAHicbVDLSsNAFJ3UV62vqBvBzWAR3FgSEdSFUOzGhYsqjS00IUymk3boZBJmJkIJceOvuHGh4tbPcOffOG2z0NYD93I4515m7gkSRqWyrG+jtLC4tLxSXq2srW9sbpnbO/cyTgUmDo5ZLDoBkoRRThxFFSOdRBAUBYy0g2Fj7LcfiJA05i01SogXoT6nIcVIack391yikN+4tI/dUCCctfy7XLeb3DerVs2aAM4TuyBVUKDpm19uL8ZpRLjCDEnZta1EeRkSimJG8oqbSpIgPER90tWUo4hIL5tckMNDrfRgGAtdXMGJ+nsjQ5GUoyjQkxFSAznrjcX/vG6qwnMvozxJFeF4+lCYMqhiOI4D9qggWLGRJggLqv8K8QDpJJQOraJDsGdPnifOSe2iZt2eVutXRRplsA8OwBGwwRmog2vQBA7A4BE8g1fwZjwZL8a78TEdLRnFzi74A+PzBwSsli0=sha1_base64="C39OhB+IczRcjLNINXH29e9lt8M=">AAAB2HicbZDNSgMxFIXv1L86Vq1rN8EiuCpTN+pOcOOygmML7VAymTttaCYzJHeEMvQFXLhRfDB3vo3pz0KtBwIf5yTk3hMXSloKgi+vtrW9s7tX3/cPGv7h0XGz8WTz0ggMRa5y04+5RSU1hiRJYb8wyLNYYS+e3i3y3jMaK3P9SLMCo4yPtUyl4OSs7qjZCtrBUmwTOmtowVqj5ucwyUWZoSahuLWDTlBQVHFDUiic+8PSYsHFlI9x4FDzDG1ULcecs3PnJCzNjTua2NL9+aLimbWzLHY3M04T+zdbmP9lg5LS66iSuigJtVh9lJaKUc4WO7NEGhSkZg64MNLNysSEGy7INeO7Djp/N96E8LJ90w4eAqjDKZzBBXTgCm7hHroQgoAEXuDNm3iv3vuqqpq37uwEfsn7+Aap5IoMsha1_base64="R99TBDumeJuZVdCYBPgayQaC/qc=">AAAB9XicbZBNS8MwHMb/nW9zTq2eBC/BIXhxtF7UgyDs4sHDlNUN1lHSLN3C0rQkqTBKvfhVvHhQ8at489uYvRx084GEh+dJSP6/MOVMacf5tkorq2vrG+XNylZ1e2fX3qs+qCSThHok4YnshFhRzgT1NNOcdlJJcRxy2g5HjUnffqRSsUS09DilvRgPBIsYwdpEgX3gU42DxpV76kcSk7wV3Bdmuy0Cu+bUnanQsnHnpgZzNQP7y+8nJIup0IRjpbquk+pejqVmhNOi4meKppiM8IB2jRU4pqqXTyco0LFJ+ihKpFlCo2n6+0aOY6XGcWhOxlgP1WI3Cf/rupmOLno5E2mmqSCzh6KMI52gCQ7UZ5ISzcfGYCKZ+SsiQ2xIaAOtYiC4iyMvG++sfll37hwowyEcwQm4cA7XcANN8IDAE7zAG7xbz9ar9TGjVbLm2Pbhj6zPH49ilMI=sha1_base64="+oMm6VNnWuShRDNcJU9+qmq+3bY=">AAACAHicbVBNS8NAEN34WetX1IvgJVgEL5bEi3oQir148FCltYUmhM120i7dbMLuRighXvwrXjyoePVnePPfuG1z0NYHMzzem2F3XpAwKpVtfxsLi0vLK6ultfL6xubWtrmzey/jVBBokZjFohNgCYxyaCmqGHQSATgKGLSDYX3stx9ASBrzphol4EW4z2lICVZa8s19FxT265fOiRsKTLKmf5frdpP7ZsWu2hNY88QpSAUVaPjml9uLSRoBV4RhKbuOnSgvw0JRwiAvu6mEBJMh7kNXU44jkF42uSC3jrTSs8JY6OLKmqi/NzIcSTmKAj0ZYTWQs95Y/M/rpio89zLKk1QBJ9OHwpRZKrbGcVg9KoAoNtIEE0H1Xy0ywDoJpUMr6xCc2ZPnSeu0elG1b+1K7apIo4QO0CE6Rg46QzV0jRqohQh6RM/oFb0ZT8aL8W58TEcXjGJnD/2B8fkDA2yWKQ== TL Scattering theory for two reservoirs

Conserved currents:

Heat currents:

First law of thermodynamics: Second law for scattering theory For two terminals:

monotonically decaying function implies For arbitrary number of terminals proof by Nenciu (2007):

The second law implies that when the reservoirs are at the same temperature the system cannot generate electrical power Joule heating: Scattering theory & Nernst’s unattainability principle

Dynamical formulation of the third law of thermodynamics: it is impossible to reach absolute zero temperature in finite time

Extracting heat from reservoir i with rate Ji, we change its temperature:

Heat capacity of a free-electron reservoir: Nernst principle satisfied (in a weak form) How to obtain the best steady-state heat to work conversion? Heat-to-work conversion through energy filtering

Flow of heat from hot to cold but no flow of charge

[see G. B., G. Casati, K. Saito, R. S. Whitney, Phys. Rep. 694, 1 (2017)] Energy filters in a thermocouple geometry What about phonons?

Necessary both: (i) to reduce phonon transport; (ii) to have an efficient working fluid (optimize the electron dynamics) Reducing thermal conductance

[Blanc, Rajabpour, Volz, Fournier, Bourgeois, APL 103, 043109 (2013)] Thermoelectric efficiency (power production) in the Landauer approach Charge current Heat current from reservoirs: 1 JJ = dE(E µ )(E)[f (E) f (E)] (null) h,q.↵ h L R Efficiency: Delta-energy filtering and Carnot efficiency If transmission is possible only inside a tiny energy window around E=E✶ then

Carnot efficiency

Carnot efficiency obtained in the limit of reversible transport (zero entropy production) and zero output power [Mahan and Sofo, PNAS 93, 7436 (1996); Humphrey et al., PRL 89, 116801 (2002)] Example: single-level quantum dot

Dot’s scattering matrix:

The Green’s function is for a non-Hermitian effective Hamiltonian taking into account coupling to the dots operator coupling the single-level dot to reservoirs:

Bekenstein-Pendry bound

There is an purely quantum upper bound on the heat current through a single transverse mode [Bekenstein, PRL 46, 923 (1981); Pendry, JPA 16, 2161 (1983) ]

For a reservoir coupled to another reservoir at T=0 through a -mode constriction which lets particle flow at all energies: Maximum power of a heat engine Since the heat flow must be less than the Bekenstein- Pendry bound and the efficiency smaller than Carnot efficiency also the output power must be bounded

Within scattering theory:

[Whitney, PRL 112, 130601 (2014); PRB 91, 115425 (2015)] Efficiency optimization (at a given power) Find the transmission function that optimizes the heat-engine efficiency for a given output power

[Whitney, PRL 112, 130601 (2014); PRB 91, 115425 (2015)] Trade-off between power and efficiency

Carnot efficiency forbidden 1

Efficiency Maximum max 2 possible power, Pgen increase voltage

power generated, Pgen

Result from (nonlinear) scattering theory

[Whitney, PRL 112, 130601 (2014); PRB 91, 115425 (2015)] Power-efficiency trade-off including phonons

no phonons

weak phonons Efficiency strong phonons 0 Power output, P

[see Whitney, PRB 91, 115425 (2015)] Boxcar transmission in topological insulators

Graphene nanoribbons with heavy adatoms and nanopores [Chang et al., Nanolett., 14, 3779 (2014)] Linear response and Landauer formalism

The Onsager coefficients are obtained from the linear response expansion of the charge and thermal currents

2 Lee = e TI0,Leh = Lhe = eT I1,Lhh = TI2 Wiedemann-Franz law

Phenomenological law: the ratio of the thermal to the electrical conductivity is directly proportional to the temperature, with a universal proportionality factor.

Lorenz number Sommerfeld expansion The Wiedemann-Franz law can be derived for low- temperature non-interacting systems both within kinetic theory or Landauer approach In both cases it is substantiated by Sommerfeld expansion. Within Landauer approach we consider

1 J = dE(E µ )(E)[f (E) f (E)] q. h L R

We assume smooth transmission functions τ(E) in the neighborhood of E=µ: To leading order in kBT/EF with

e2 1 I2 2k2 T G = e2I (µ),K= I 1 B (µ) 0 h T 2 I 3h 0

Neglected I12/I0 with respect to I2, which in turn implies LeeLhh>>(Leh)2 and Wiedemann-Franz law: K 2 k 2 B T G 3 e Wiedemann-Franz law and thermoelectric efficiency

GS2 S2 ZT = T = K L Wiedemann-Franz law derived under the condition LeeLhh>>(Leh)2 and therefore

Wiedemann-Franz law violated in - low-dimensional interacting systems that exhibit non- Fermi liquid behavior - (smll) systems where transmission can show significant energy dependence (Violation of) Wiedemann-Franz law in small systems

Consider a (basic) model of a molecular wire coupled to electrodes:

Transmission: Green’s function: Level broadening functions: Self-energies: Wide band limit: level widths energy independent:

Take

Transmission: Green’s function obtained by inverting

Mott’s formula for thermopower For non-interacting electrons (thermopower vanishes when there is particle-hole symmetry) f 1 I 1 dE(E µ)(E) E S = 1 = f eT I0 eT dE(E) E Consider smooth transmissions

Electron and holes contribute with opposite signs: we want sharp, asymmetric transmission functions to have large S (ex: resonances, Anderson QPT, see Imry and Amir, 2010), violation of WF, possibly large ZT. Metal-insulator 3D Anderson transition

x conductivity critical exponent [G.B., H. Ouerdane, C. Goupil, Comptes Rendus Physique 17, 1072 (2016)] Energy filtering

For good thermoelectric we desire violation of WF law such that:

No dispersion with delta-energy filtering: ZT diverges III. Rate equations Conditions for a

Consider systems weakly coupled to the environment. Non-Markovian effects are neglected as well as the generation of coherences (in the system’s energy eigenbasis). Interaction effects may be included. Rate equations

Transition rate from state a to state b (due to the coupling to reservoir i):

The probability to find the system in state |b⟩ at time t obeys a rate equation: Local detailed balance The rates can be derived from microscopic Hamiltonian v i a F e r m i g o l d e n r u l e , o r b e c o n s i d e r e d phenomenological constants. They obey the local detailed balance principle:

change of entropy in reservoir i when it induces a system’s transition from a to b.

Clausius relation: (Steady-state) currents Probability current at time t for the transition from a to b induced by reservoir i:

Steady-state solution of the rate equations:

Steady-state probability currents:

Kirkhoff’s law: Particle and energy currents

By taking the steady-state probability currents, we obtain the steady-state charge and heat currents:

There is no net flow of particle and energy into the system at equilibrium: Equilibrium (dynamical definition) Connect the system to reservoirs at the same temperature and electrochemical potential

A state at equilibrium must obey the detailed balance:

From the local detailed balance:

We then derive the equilibrium state Output power and the first law of thermodynamics

The power generated at reservoir i depends on the reference electrochemical potential, but the overall power is gauge-invariant: This is a consequence of Kirkoff’s law and of the relation: We obtain the first law of thermodynamics: the rate of heat absorption equals the rate of work production: Second law of thermodynamics Change of entropy in reservoir i:

System in general not in a thermal state: use Shannon entropy rather than Clausius definition:

Total entropy production rate:

Use local detailed balance to prove that Example: single-level quantum dot

Neglect electron spin; charging energy too high to have doble occupancy

Rate equations for the dot’s dynamics: Particle and energy currents

Particle currents: Set Heat currents:

Entropy production at steady state: Power and efficiency

Generated power:

Efficiency of a heat engine:

Coefficient of performance for a refrigerator: Steady-state solution

Use ⇒

Using local detailed balance: Carnot efficiency Reversible engine if the entropy production rate vanishes:

This leads to Carnot efficiency:

The price to pay is that the power is zero:

E E (N N )µ ✏ 0 1 0 1 L = 1

(null) TL TL

E E (N N )µ ✏ + µ ✏ + ✏ (1 T /T ) ✏ = 0 1 0 1 R = 1 = 1 1 R L = 1

(null) TR TR TR TL ⇒ Multilevel interacting quantum dot

Discrete energy levels: ideal to implement energy filtering Study the effects of Coulomb interaction between electrons

[Erdmann, Mazza, Bosisio, G.B., Fazio, Taddei PRB 95, 245432 (2017)] Sequential (single-electron) tunnelling regime single-electron levels of the QC capacitance number of electrons in the dot electrostatic (Coulomb) interaction tunneling rate from level p to reservoir � Weak coupling to the reservoirs: thermal energy , level spacing and charging energy much larger than the coupling energy between the QD and the reservoirs: charge quantized

Electrostatic energy single-electron charging energy Energy conservation

Configuration determined by occupation numbers

Non-equilibrium probability

Energy conservation for tunnelling into or from reservoirs: Kinetic (rate) equations One kinetic (rate) equation for each configuration:

Stationary solution: Steady-state currents Charge current:

Energy current:

Heat current: Quantum limit Energy spacing and charging energy much bigger than Analytical results for equidistant levels:

power factor

(energy filtering) Coulomb interaction may enhance the thermoelectric performance of a QD Compare interacting and non-interacting two-terminal QD with the same energy spacing Thermal conductance suppressed by Coulomb interaction: ZT is greatly increased. For a single level K=0 (charge and heat current proportional). For at least two levels Coulomb blockade prevents a second electron to enter when one is already there (electrostatic energy to be paid). IV. Thermodynamic bounds on heat-to-work conversion Can interactions improve the power-efficiency trade-off?

What is the role of a magnetic field? Is it possible to have Carnot at finite power?

What is the role played by fluctuations? Thermodynamic uncertainty relations Short intermezzo: a reason why interactions might be interesting for thermoelectricity

thermal conductance at zero voltage

If the ratio K’/K diverges, then the Carnot efficiency is achieved Thermodynamic properties of the working fluid

coupled equations: Setting dN=0 in the coupled equations: Thermodynamic cycle

maximum efficiency (over d� at fixed dT):

thermodynamic figure of merit: Analogy with a classical gas

heat capacity at constant p or V

Power-efficiency trade-off: Is it possible to overcome the non-interacting bound?

Noninteracting systems: for P/Pmax<<1,

[Whitney, PRL 112, 130601 (2014); PRB 91, 115425 (2015)] Bound not favorable for power-efficiency trade-off; due to the fact that delta-energy filtering is the only mechanism to achieve Carnot for noninteracting systems For interacting systems it is possible to achieve Carnot without delta-energy filtering Interacting systems, Green-Kubo formula The Green-Kubo formula expresses linear response transport coefficients in terms of dynamic correlation functions of the corresponding current operators, cal- culated at thermodynamic equilibrium

Non-zero generalized Drude weights signature of ballistic transport Conservation laws and thermoelectric efficiency Suzuki’s formula (which generalizes Mazur’s inequality) for finite-size Drude weights

Qm relevant (i.e., non-orthogonal to charge and thermal currents), mutually orthogonal conserved quantities

Assuming commutativity of the two limits, Momentum-conserving systems Consider systems with a single relevant constant of motion, notably momentum conservation Ballistic contribution to vanishes since D D D2 =0 ee hh eh

2 S 1 ZT = T when

( < 1) (G.B., G. Casati, J. Wang, PRL 110, 070604 (2013)) For systems with more than a single relevant constant of motion, for instance for integrable systems, due to the Schwarz inequality D D D2 = x 2 x 2 x , x 0 ee hh eh || e|| || h|| e h

1 JiQ1 JiQM xi =(xi1, ..., xiM )= , ..., 2 Q2 Q2 1 M M x , x = x x e h ek hk k=1 Equality arises only in the exceptional case when the two vectors are parallel; in general det L L2,2 , ZT 0 Example: 1D interacting classical gas

Consider a one dimensional gas of elastically colliding particles with unequal masses: m, M

(integrable model) ZT =1(at µ = 0)

ZT depends on the system size Quantum mechanics needed: Relation between density and electrochemical potential Reservoirs modeled as ideal (1D) gases Maxwell-Bolzmann distribution of injection rates

grand partition function density

de Broglie thermal wave length Non-decaying correlation functions Carnot efficiency at the thermodynamic limit

Anomalous S2 ZT = T thermal transport k [R. Luo, G. B., G. Casati, J. Wang, PRL 121, 080602 (2018)] Delta-energy filtering mechanism?

A mechanism for achieving Carnot different from delta-energy filtering is needed Validity of linear response

The agreement with linear response improves with N Non-interacting classical bound (but quantum mechanics needed)

charge current

heat current

Maxwell-Boltzmann distribution

(in 1D) 1.0

0.8

0.6 C h ê h 0.4

0.2

0.0 0.0 0.2 0.4 0.6 0.8 1.0 P Pmax

ê Overcoming the non-interacting bound

Non-interacting bound [by linear response]

AAACA3icbVA9T8MwEHX4LOUrwNjFokJiqhKEBGwVLIxFIrRSG1WOe22tOnZkO4gq6sDCX2FhAMTKn2Dj3+C0GaDlSSc9vXdn370o4Uwbz/t2lpZXVtfWSxvlza3tnV13b/9Oy1RRCKjkUrUiooEzAYFhhkMrUUDiiEMzGl3lfvMelGZS3JpxAmFMBoL1GSXGSl230hlG8iFrR2OcP0EUVqATKTSEk65b9WreFHiR+AWpogKNrvvV6UmaxiAM5UTrtu8lJsyIMoxymJQ7qYaE0BEZQNtSQWLQYTY9YoKPrNLDfalsCYOn6u+JjMRaj+PIdsbEDPW8l4v/ee3U9M/DjIkkNSDo7KN+yrGROE8E95gCavjYEkIVs7tiOiSKUGNzK9sQ/PmTF0lwUruoeTen1fplkUYJVdAhOkY+OkN1dI0aKEAUPaJn9IrenCfnxXl3PmatS04xc4D+wPn8AagdmE4=sha1_base64="C39OhB+IczRcjLNINXH29e9lt8M=">AAAB2HicbZDNSgMxFIXv1L86Vq1rN8EiuCpTN+pOcOOygmML7VAymTttaCYzJHeEMvQFXLhRfDB3vo3pz0KtBwIf5yTk3hMXSloKgi+vtrW9s7tX3/cPGv7h0XGz8WTz0ggMRa5y04+5RSU1hiRJYb8wyLNYYS+e3i3y3jMaK3P9SLMCo4yPtUyl4OSs7qjZCtrBUmwTOmtowVqj5ucwyUWZoSahuLWDTlBQVHFDUiic+8PSYsHFlI9x4FDzDG1ULcecs3PnJCzNjTua2NL9+aLimbWzLHY3M04T+zdbmP9lg5LS66iSuigJtVh9lJaKUc4WO7NEGhSkZg64MNLNysSEGy7INeO7Djp/N96E8LJ90w4eAqjDKZzBBXTgCm7hHroQgoAEXuDNm3iv3vuqqpq37uwEfsn7+Aap5IoMsha1_base64="GwhlYFXBSeBmFoHvBy/FdU0sTd0=">AAAB+HicbZC9TgJBFIXv+ouIirY0E4mJFVls1M7ExhITV0hgQ2aHC0yYndnMzBo3GwobX8XGQo1PYufbOAsUCp5kki/nzN89USK4sb7/7a2tb2xubZd2yruVvf2D6mHl3qhUMwyYEkp3ImpQcImB5VZgJ9FI40hgO5pcF3n7AbXhSt7ZLMEwpiPJh5xR66x+tdYbR+ox70YZKa6gmmg0iZIGw2m/Wvcb/kxkFZoLqMNCrX71qzdQLI1RWiaoMd2mn9gwp9pyJnBa7qUGE8omdIRdh5LGaMJ8NsSUnDhnQIZKuyUtmbm/T+Q0NiaLI7czpnZslrPC/C/rpnZ4EeZcJqlFyeYPDVNBrCJFI2TANTIrMgeUae7+StiYasqs663sSmguj7wKwVnjsuHf+lCCGhzDKTThHK7gBloQAIMneIE3ePeevVfvY97Wmreo7Qj+yPv8ARwBluA=sha1_base64="hsfpqjk/0uYw3XODntxaikWhXUA=">AAACA3icbVA9T8MwEHXKVylfAcYuFhUSU5WyAFsFC2ORCK3URpXjXlurjh3ZDiKKOrDwV1gYALHyJ9j4NzhtBmh50klP793Zdy+MOdPG876d0srq2vpGebOytb2zu+fuH9xpmSgKPpVcqk5INHAmwDfMcOjECkgUcmiHk6vcb9+D0kyKW5PGEERkJNiQUWKs1HervXEoH7JumOL8CaKwAh1LoSGY9t2aV/dmwMukUZAaKtDqu1+9gaRJBMJQTrTuNrzYBBlRhlEO00ov0RATOiEj6FoqSAQ6yGZHTPGxVQZ4KJUtYfBM/T2RkUjrNAptZ0TMWC96ufif103M8DzImIgTA4LOPxomHBuJ80TwgCmghqeWEKqY3RXTMVGEGptbxYbQWDx5mfin9Yu6d+PVmpdFGmVURUfoBDXQGWqia9RCPqLoET2jV/TmPDkvzrvzMW8tOcXMIfoD5/MHpt2YSg== Multiparticle collision dynamics (Kapral model) in 2D

S t r e a m i n g s t e p : free propagation during a time τ

Collision step: random rotations of the velocities of the particles in cells of linear size a with respect to the center of mass :

Momentum is conserved Overcoming the (2D) non-interacting bound Results can be extended to cooling

linear response

numerical data Applications for cold atoms? Power-efficiency trade-off at the verge of phase transitions For heat engines described as Markov processes:

[N. Shiraishi, K. Saito, H. Tasaki, PRL 117, 190601 (2016)]

For a working substance at a critical point:

[M. Campisi, R. Fazio, Nature Comm. 7, 11895 (2016); see also Allahverdyan et al., PRL 111, 050601 (2013)] Results compatible only with diverging amplitude A when approaching the Carnot efficiency Power-efficiency-fluctuations trade-off For classical Markovian dynamics on a finite set of states and overdamped Langevin dynamics, trade-off between power, efficiency, and constancy for steady- state engines:

[P. Pietzonka, U. Seifert, PRL 120, 190602 (2018)] Bound violated in quantum mechanics, e.g. for resonant tunnelling transport (noninteracting system), but not close to Carnot efficiency. The problem for generic quantum systems is still open. [J. Liu. D. Segal, PRE 99, 062141 (2019)] Carnot efficiency at finite power with broken time-reversal symmetry?

J = L (B) + L (B) e ee Fe eh Fh J = L (B) + L (B) h he Fe hh Fh µ = µL µR B applied magnetic field or any parameter breaking time-reversibility T = T T such as the Coriolis force, etc. L R (we assume TL >TR,µL <µR) Constraints from thermodynamics

POSITIVITY OF THE ENTROPY PRODUCTION:

ONSAGER-CASIMIR RELATIONS:

G(B)=G( B) K(B)=K( B) Both maximum efficiency and efficiency at maximum power depend on two parameters

Pmax

The CA limit can be overcome within linear response

When |x| is large the figure of merit y required to get Carnot efficiency becomes small Carnot efficiency could be obtained far from the tight coupling condition [G..B., K. Saito, G. Casati, PRL 106, 230602 (2011)] Output power at maximum efficiency

When time-reversibility is broken, within linear response it is not forbidden from the second law to have simultaneously Carnot efficiency and non-zero power. Terms of higher order in the entropy production, beyond linear response, will generally be non-zero. However, irrespective how close we are to the Carnot efficiency, we could in principle find small enough forces such that the linear theory holds. Reversible part of the currents

The reversible part of the currents does not contribute to entropy production

Possibility of dissipationless transport?

[K.. Brandner, K. Saito, U. Seifert, PRL 110, 070603 (2013)] How to obtain asymmetry in the Seebeck coefficient?

For non-interacting systems, due to the symmetry properties of the scattering matrix

This symmetry does not apply when electron-phonon and electron-electron interactions are taken into account

Let us consider the case of partially coherent transport, with inelastic processes simulated by “conceptual probes” (Buttiker, 1988). Physical model of probe reservoirs Large but finite “reservoirs”

Voltage and temperature probe (mimicking e-e inelastic scattering)

Voltage probe (mimicking e-ph inelastic scattering) Non-interacting three-terminal model

Je,k =0, Ju,k =0 (J = J (µ/e)J ) h,k u,k e,k k k 4 t J =(JeL,JhL,JeP ,JhP ) t J = Ji i µ T µ T F F t = , , P , P i=1 F eT T 2 eT T 2 Three-terminal Onsager matrix Equation connecting fluxes and thermodynamic forces: J = L F

In block-matrix form: J L L = F J L L F Zero-particle and heat current condition through the probe terminal:

1 = L L F F Two-terminal Onsager matrix for partially coherent transport Reduction to 2x2 Onsager matrix when the third terminal is a probe terminal mimicking inelastic scattering

1 J = L , L L L L L . F J L L 1 = 11 12 F1 J L L 2 21 22 F2 First-principle exact calculation within the Landauer-Büttiker approach

Bilinear Hamiltonian Charge and heat current from the left terminal Onsager coefficients from linear response expansion of the currents

Transmission probabilities: Illustrative three-dot example

Asymmetric structure, e.g.. Asymmetric Seebeck coefficient

[K. Saito, G. B., G. Casati, T. Prosen, PRB 84, 201306(R) (2011)] [see also D. Sánchez, L. Serra, PRB 84, 201307(R) (2011)] Asymmetric power generation and refrigeration

When a magnetic field is added, the efficiencies of power generation and refrigeration are no longer equal:

To linear order in the applied magnetic field:

A small magnetic field improves either power generation or refrigeration, and vice versa if we reverse the direction of the field The large-field enhancement of efficiencies is model- dependent, but the small-field asymmetry is generic Transmission windows model

ij(E)= ij(E)=1 The Curzon-Ahlborn limit i j can be overcome (within linear response)

[V.. Balachandran, G. B., G. Casati, PRB 87, 165419 (2013); se also M. Horvat, T. Prosen, G. B., G. Casati, PRE 86, 052102 (2012)] Saturation of bounds from the unitarity of S-matrix

Bounds obtained for non-interacting 3-terminal transport (K. Brandner, K. Saito, U. Seifert, PRL 110, 070603 (2013))

4 4 (P⇥max)= C at x = max 7 3 Bounds for multi-terminal thermoelectricity

Numerical evidence that the power vanishes when the Carnot efficiency is approached

[Brandner and Seifert, NJP 15, 105003 (2013); PRE 91, 012121 (2015)] Bounds with electron-phonon scattering

Efficiency bounded by the non-negativity of the entropy production of the original three-terminal junction. However, the efficiency at maximum power can be enhanced [Yamamoto, Entin-Wohlman, Aharony, Hatano; PRB 94, 121402(R) (2015)] Onsager-Casimir relations

Onsager reciprocal relations reflect at the macroscopic level the time-reversal symmetry of the microscopic dynamics, invariant under the transformation:

With an applied magnetic field one instead obtains Onsager-Casimir relations:

but in principle one could violate the Onsager symmetry: Onsager relations and thermodynamic constraints on heat-to-work conversion For thermoelectricity:

and in principle one could have the Carnot efficiency at finite power: Onsager relations with broken time-reversal symmetry Onsager relations under an applied magnetic field remain valid (for two terminals, Hamiltonian dynamics): 1) for noninteracting systems

2) if the magnetic field is constant [Bonella, Ciccotti, Rondoni, EPL 108, 60004 (2014)] What about for a generic, spatially dependent magnetic field? Symmetry without magnetic field inversion

Analytical result for Landau gauge: Equations of motion invariant under: Numerical results

generic 2D case:

generic 3D case: Theoretical argument: divide the system into small volumes Time-reversal trajectories without reversing the field for No-go theorem for finite power at the Carnot efficiency on purely thermodynamic grounds?

Onsager reciprocal relations much more general than expected so far. Some open problems

Investigate strongly-interacting systems close to electronic phase transitions

Find suitable (numerical) methods to address the strong system-reservoir coupling regime

Further investigate/optimize time-dependent driving

Power-efficiency-fluctuation trade-off for non- Markovian quantum dynamics? Quantum computation and information is a quantum cryptography) up to advanced rapidly developing interdisciplinary field. It topics (like entanglement measures, quan- Benenti is not easy to understand its fundamental tum discord, quantum noise, quantum Casati concepts and central results without facing channels, quantum error correction, quan- numerous technical details. This book tum simulators, and tensor networks). Rossini provides the reader with a useful It can be used as a broad range Strini Giuliano Benenti Giulio Casati guide. In particular, the initial textbook for a course in quantum chapters offer a simple and self- Davide Rossini Giuliano Strini information and computation, contained introduction; no previous both for upper-level undergraduate Principles of Quantum Computation and Information knowledge of quantum mechanics students and for graduate or classical computation is required. students. It contains a large Various important aspects of quan- number of solved exercises, which tum computation and information are an essential complement to the text, as are covered in depth, starting from the foun- they will help the student to become dations (the basic concepts of computational familiar with the subject. The book may also complexity, energy, entropy, and information, be useful as general education for readers quantum superposition and entanglement, who want to know the fundamental elementary quantum gates, the main quan- principles of quantum information and

tum algorithms, quantum teleportation, and computation. A Comprehensive Textbook SATOR

“Thorough introductions to classical computation and irreversibility, and a primer of quantum AREPO theory, lead into the heart of this impressive and substantial book. All the topics – quantum algorithms, quantum error correction, adiabatic quantum computing and decoherence are just TENET a few – are explained carefully and in detail. Particularly attractive are the connections between the conceptual structures and mathematical formalisms, and the different experimental protocols for bringing them to practice. A more wide-ranging, comprehensive, and definitive OPERA text is hard to imagine.” –— Sir Michael Berry, University of Bristol, UK ROTAS “This second edition of the textbook is a timely and very comprehensive update in a rapidly developing field, both in theory as well as in the experimental implementation of quantum information processing. The book provides a solid introduction into the field, a deeper insight in the formal description of quantum information as well as a well laid-out overview on several platforms for quantum simulation and quantum computation. All in all, a well-written and commendable textbook, which will prove very valuable both for the novices and the scholars in the fields of quantum computation and information.” –— Rainer Blatt, Universität Innsbruck and IQOQI Innsbruck, Austria Principles of Quantum Computation

“The book by Benenti, Casati, Rossini and Strini is an excellent introduction to the fascinating and Information field of quantum information, of great benefit for scientists entering the field and a very useful reference for people already working in it. The second edition of the book is considerably A Comprehensive Textbook extended with new chapters, as the one on many-body systems, and necessary updates, most notably on the physical implementations. ” –— Rosario Fazio, Te Abdus Salam International Centre for Teoretical Physics, Trieste, Italy

World Scientifc ISBN 978-981-3237-22-3 www.worldscientific.com World Scientific 10909 hc https://tqsp.lakecomoschool.org/

International School on THERMODYNAMICS OF QUANTUM SYSTEMS AND PROCESSES Como (Italy), August 31 – September 4, 2020 Lectures Alexia Auffèves (Institut Néel, CNRS Grenoble, France) Quantum measurement and thermodynamics Michele Campisi (Università di Firenze, Italy) Fluctuation relations in quantum systems Ronnie Kosloff (Hebrew University of Jerusalem, Israel) Quantum thermal machines Massimo Palma (Università di Palermo, Italy) Non-Markovianity and thermodynamics of open quantum systems Jukka Pekola (Aalto University, Finland) Non equilibrium thermodynamics of quantum circuits Ferdinand Schmidt-Kaler (University of Mainz, Germany) Experimental stochastic and quantum engines: from single particle to many body Robert Whitney (Université Grenoble-Alpes, France) Quantum thermoelectricity

Important dates Deadline for application: May 15, 2020 Acceptance notification: May 25, 2020 Organizers Giuliano Benenti, Università dell’Insubria (IT) Dario Gerace, Università di Pavia (IT) Mauro Paternostro, Queen’s University of Belfast (UK) Jukka Pekola, Aalto University (FI) School website https://tqsp.lakecomoschool.org/