Next Article in Journal
On the Existence and Uniqueness of Solutions for Local Fractional Differential Equations
Next Article in Special Issue
Quantum Thermodynamics with Degenerate Eigenstate Coherences
Previous Article in Journal
Simple Harmonic Oscillator Canonical Ensemble Model for Tunneling Radiation of Black Hole
Previous Article in Special Issue
SU(2) Yang–Mills Theory: Waves, Particles, and Quantum Thermodynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Periodic Energy Transport and Entropy Production in Quantum Electronics

by
María Florencia Ludovico
1,
Liliana Arrachea
1,
Michael Moskalets
2 and
David Sánchez
3,*
1
International Center for Advanced Studies, UNSAM, Campus Miguelete, 25 de Mayo y Francia, 1650 Buenos Aires, Argentina
2
Department of Metal and Semiconductor Physics, NTU “Kharkiv Polytechnic Institute”, 61002 Kharkiv, Ukraine
3
Instituto de Física Interdisciplinar y Sistemas Complejos IFISC (UIB-CSIC), E-07122 Palma de Mallorca, Spain
*
Author to whom correspondence should be addressed.
Entropy 2016, 18(11), 419; https://doi.org/10.3390/e18110419
Submission received: 4 October 2016 / Revised: 14 November 2016 / Accepted: 15 November 2016 / Published: 23 November 2016
(This article belongs to the Special Issue Quantum Thermodynamics)

Abstract

:
The problem of time-dependent particle transport in quantum conductors is nowadays a well established topic. In contrast, the way in which energy and heat flow in mesoscopic systems subjected to dynamical drivings is a relatively new subject that cross-fertilize both fundamental developments of quantum thermodynamics and practical applications in nanoelectronics and quantum information. In this short review, we discuss from a thermodynamical perspective recent investigations on nonstationary heat and work generated in quantum systems, emphasizing open questions and unsolved issues.

1. Introduction

The rapid progress in the manipulation of small systems at submicron scales has spurred the interest in the field of quantum thermodynamics. The goals of this discipline is the understanding of the fundamentals of thermodynamics in quantum systems typically driven out of equilibrium by means of external fields. Most of the recent literature focuses on static fields and the resulting stationary transport effects (for a recent review, see [1]). However, there is a growing interest in analyzing the thermodynamic properties of quantum conductors in the presence of time-dependent potentials. In this case, dynamics is the main objective of the theory as fluxes and responses depend explicitly on time [2].
Quantum electronic devices are paradigmatic realizations of open quantum systems. Typically, we can identify in these setups a finite-size piece where a small number of electrons are confined. However, these particles are not isolated. They are in contact with leads, substrates and the electromagnetic environment, which in practice constitute macroscopic thermal and particle reservoirs. We can distinguish two main classes of time-dependent processes in these systems. One of them is related to the transient behavior after a short-time non-equilibrium perturbation, where the main issue is to understand the relaxation and thermalization processes as the system approaches the equilibrium state. The other class is the response to a periodic driving. Both classes of dynamics have received a significant interest for some years now.
From the conceptual point of view, it is quite natural to relate the energy transport in periodically-driven systems to thermodynamic machines, like heat engines and refrigerators, since these typically operate in cycles. There is a lengthy literature on nanoscale quantum thermodynamics based on two-level systems [3,4,5,6,7], atomic [8,9,10,11] and molecular [12,13,14] systems, nanomechanical [15,16,17,18,19,20] systems and other systems represented by harmonic oscillators [21,22,23,24]. The list of works devoted to study cyclic processes in electron systems is shorter but is growing every week. In [25] heat transport induced by adiabatic pumping is analyzed in a quantum dot driven by two AC voltages operating with a phase lag. Two interesting operational modes were identified: (i) the heat pumping against a gradient of temperature and (ii) the exchange of work between the driving forces. Heat pumping by AC driving was also investigated in the non-adiabatic regime [26] and in arrays of quantum dots in the Coulomb blockade regime [27], while in [28] the effect of the exchange of power between two driven quantum capacitors was analyzed. In [29] the exchange of power between electron currents induced by DC bias voltages and adiabatic driving forces was considered and a similar mechanism was later analyzed in the context of a pair of helical edge states coupled to a magnetic island [30]. The effect of AC driving in the efficiency of usual two-terminal thermoelectric setups was addressed in [27,31,32]. An extended thermoelectric framework to describe the work exchanged with the AC driving sources within the adiabatic regime, on equal footing with the exchange of charge and heat between electron reservoirs with thermal and voltage biases was proposed in [33]. Thermal transport in electron systems driven by the coupling to electromagnetic radiation has been studied in [34,35,36,37]. The related problems of heat fluctuations and the definition of effective temperatures to characterize the net heat transport in these systems were addressed in [38,39,40,41,42,43,44,45,46,47,48] and [49,50,51], respectively.
In addition to the interest in the fundamental aspects, it is of paramount importance for potential applications to discriminate which portion of the energy invested to operate these devices is amenable to be used and which one is wasted by dissipation. This distinction is at the heart of thermodynamics and is conventionally addressed in quasi-stationary processes where the system under study is very weakly coupled to the reservoirs. In quantum electronics, however, the generic situation is to have the driven structure strongly coupled to the rest of the circuit, which plays the role of the reservoir. On the other hand, we are typically interested in the generation of currents, which implies non-equilibrium situations.
In this paper, we review recent developments in dynamical quantum thermodynamics applied to nanoelectronic systems. We will discuss in detail the energy transfer in out-of-equilibrium systems coupled to external baths. The analysis is based on a phase coherent mesoscopic sample (a confined system with discrete energy levels) attached to fermionic reservoirs held at a given temperature. The energies of the sample evolve with time due to the coupling with nearby AC gate terminals. Deep inside the reservoirs, electrons relax their excess energy, and the baths can thus be considered in local thermal equilibrium. We will also consider the entropy production in the whole system and will identify the different terms arising in the redistributed energy and heat [52,53]. Importantly, when the energies shift slowly with time, the response is adiabatic, and an exact Joule law can be demonstrated for the time domain [52,53,54]. In addition, we will briefly discuss the nonadiabatic regime [55] in which case the AC frequency is larger or of the order of the inverse dwell time inside the conductor. Our analysis is completely general and does not rely on the particular approach followed to evaluate the relevant dynamical quantities, like the charge and energy fluxes and the rate of entropy production.
Dynamical transport in quantum conductors is a relatively old subject with established theoretical principles [56,57,58,59], which has been investigated with different techniques [60,61,62,63]. It is well understood that a Floquet scattering matrix approach (FSM) [64] can account for photon-assisted phenomena in nanostructures (absolute negative conductance and dynamical localization [62], single-electron emission [54,65], etc.). In the present work, we also include for completeness the details of the calculation of these quantities with the Keldysh formalism using the Floquet-Fourier representation as in [66,67] and discuss the connection to the FSM formulation as in [68].

2. Theoretical Model

We consider that our system is a slowly-driven (adiabatically) quantum conductor (e.g., a quantum dot or an array of quantum dots) coupled to several baths (fermionic reservoirs), see Figure 1 for the two terminal case. In the adiabatic regime, it is possible to propose a thermodynamic function within an extended resonant level model that satisfies both the first and the second laws [52,53]. We do not make any assumption on the coupling strength between the system and the baths.
Note that the problem of defining a system and a bath in the strong coupling regime has been addressed recently [69]. Moreover, the authors of [70] claim that the extended resonant level model is unable to properly describe energy fluctuations despite the fact that the expectation value of the energy is correctly derived. On the other hand, [71] finds that an influence functional approach [6,72] applied to the spin-boson model can yield exact expressions for the energy dissipation in a driven open quantum system strongly coupled to a heat bath.
Here, we take into account the Hamiltonian of the fermionic baths, H res = α H α with H α = k α ε k α c k α c k α , where ε k α is the energy dispersion relation and c k α ( c k α ) creates (destroys) an electron with wavenumber k α , the Hamiltonian of the system, H S ( V ( t ) ) , where the functions V ( t ) = V ( t + τ ) = ( V 1 ( t ) , , V M ( t ) ) describe the AC driving power sources capacitively coupled to the system, τ being the AC period, and the Hamiltonian of the contact region between the system and the baths, H cont = α H c α with H c α = k α w k α c k α d l α + H . c , where w k α are tunneling amplitudes, l α denotes the site of the central part to which the reservoir α is attached, and d l α , and d l α destroys and creates, respectively, an electron in the central system. The three contributions add up to the Hamiltonian of the full system,
H ( t ) = H res + H S ( t ) + H cont ,
which explicitly depends on time due to the applied AC potentials. We stress that the explicit form of H S remains unspecified, because we will below present a thermodynamic discussion that is completely general and independent of the sample details. In fact, our presentation is also valid for interacting electrons provided these carriers are confined to the mesoscopic region. The fermionic baths are metallic leads with good screening properties and electron-electron interactions can be safely neglected there. The latter is an important remark since a gauge-invariant transport theory in the AC regime must include the self-consistent Coulomb interaction [73]. This requirement in the case of the heat current was discussed in [74].
Detailed knowledge of the heat dynamics in quantum electronics is also relevant for time-dependent thermoelectrics, where temperatures may be modulated with time. A calculation for quantum R C circuits shows [75] that thermoelectric charge relaxation resistance significantly differs from its purely electric counterpart. In a different setup, the thermoelectric power factor has been experimentally demonstrated to grow in the photo-induced regime [76]. When the coupling between a dot and its attached leads quickly changes with time, transient currents are predicted to occur in the absence of tunneling [31]. These findings call for a unified description of thermoelectrics and heat for AC-driven conductors. This is accomplished in [33], where the authors generalize the Onsager reciprocity relations for the case of adiabatic drivings. Thermodynamic forces and fluxes can be identified this way [77]. Unfortunately, experimental measurements of thermoelectric responses in the time domain are scarce. An exception is [78], where a frequency-resolved impedance spectroscopy method is used in thermoelectric modules.

3. First Law

Let ρ = e β ( H μ N ) / Z be the density matrix, with Z = Tr e β ( H μ N ) the partition function, N the particle number and β = 1 / ( k B T ) . We start by considering μ L = μ = R = μ and T L = T R = T . Variations δ ρ = ρ ( t + δ t ) ρ ( t ) are due to changes in δ V = V ( t + δ t ) V ( t ) , which we assume to be small. Then, both the particle number and internal energy change according to:
δ N ν = Tr δ ρ N ν , ν = α , S
and:
δ U ν = Tr δ ρ H ν , ν = α , c α , S
respectively. The particle number and energy for the full system remain constant, since the system plus the baths can be considered as a larger, closed system.
δ N = α δ N α + δ N S ,
and
δ U = Tr δ ρ H = α δ U α + δ U c α + δ U S ,
The first law of thermodynamics states that total heat in the system is:
δ Q tot | μ , T = α δ U α + δ U c α μ δ N α + δ U S δ W AC μ δ N S ,
where:
δ W AC = Tr [ ρ H / V ] δ V
is the work developed by the time-dependent driving sources. The corresponding change in the entropy reads:
T δ S | μ , T = δ Q tot | μ , T = δ W AC ,
where we have used in the last equality the conservation of particles and energy expressed by Equations (2) and (3).
We now consider a more general situation, where the reservoirs have slightly different chemical potentials μ α = μ + δ μ α . This implies investing an extra electric work done by DC bias voltages
δ W el = α δ μ α δ N α .
Here, δ μ α is a small departure of the electrochemical potential of lead α with respect to the common chemical potential μ.
Equation (8) generalizes to:
T δ S | μ , T + T δ S | δ μ = δ Q tot
with
T δ S | δ μ = α δ μ α δ N α = δ W el .
Since δ N = δ U = 0 due to particle and energy conservation, the first law implies that the total heat generated in the full system stems from the work done by the AC and the DC sources,
δ Q tot = δ W AC δ W el .

4. Heat Current and Power

The rate of change for the energy corresponding to the different parts of the system (sample S, reservoir α or coupling region c α ) is given by the exact quantum-mechanical expression
J ν E ( t ) = i H , H ν ,
where ν S , α , c α . The time variation of the expected value of the total Hamiltonian H ,
H ˙ = α J α E ( t ) + J c α E ( t ) + J S E ( t ) F · V ˙ ,
is given by the explicit time derivative of H in terms of the generalized force F = H / V . On the other hand, we know that the change in time of the total energy stored in the full setup, containing the central system, contacts and reservoirs, is equal to the total power developed by the external AC sources, H ˙ = P AC ( t ) . The latter reads:
P AC ( t ) = F · V ˙ .
Hence, it follows that the fluxes obey:
α J α E ( t ) + J c α E ( t ) + J S E ( t ) = 0 .
The time variation of the charges present in the different pieces of the system can be similarly inferred. We start from:
I ν C ( t ) = e N ν ˙ = i e H , N ν ,
with ν = α , S . Conservation of the total charge demands that the currents fulfill:
e N ˙ = I S C ( t ) + α I α C ( t ) = 0 .
Noticeably, the coupling region contributes to the energy current balance (Equation (16)) but not to its particle counterpart (Equation (18)). This will have important consequences for the heat current definition, as explained below.
Using Equation (6) it is straightforward to derive the total heat flux,
Q ˙ tot ( t ) = α J α E ( t ) μ α I α C ( t ) e + J c α E ( t ) + J S E ( t ) P AC ( t ) μ I S C ( t ) e .
Now, we define:
P el ( t ) = α δ μ α I α C ( t ) / e ,
which is the electric power developed by the batteries applied to the fermionic baths to keep the electron flow. It follows from Equations (16) and (18) that
Q ˙ tot ( t ) = P AC ( t ) P el ( t ) ,
which states that the total heat flux in the system is produced by the AC and DC driving powers, as expected.
The question now is how to define physically-meaningful time-dependent currents flowing through different pieces of the setup. The work in [52] demonstrates that the conventional definition for the heat flowing in the reservoir α,
Q ˜ ˙ α ( t ) = J α E ( t ) μ α I α C ( t ) e ,
suffers from serious drawbacks. In particular, for a slowly-driven quantum level coupled to a single reservoir the use of Q ˜ ˙ leads to negative heat fluxes at zero temperature. This would imply that heat would be extracted from a reservoir at zero temperature, which makes no sense if we are to interpret this heat as dissipation. In contrast, when one considers the expression
Q ˙ α ( t ) = J α E ( t ) + J c α E ( t ) 2 μ α I α C ( t ) e .
one finds that the heat flux is always positive at zero temperature, i.e., energy enters into the reservoir and is evacuated from the power source, as should be. Moreover, Equation (23) is a sound expression that fulfills four nice properties:
  • It leads at low frequencies to a correct Joule law valid for all times [52].
  • It shows perfect agreement between the Green function approach and the scattering matrix formalism [52].
  • It displays parity symmetry upon reversal of the AC frequency even for interacting quantum conductors [79].
  • It reduces to the conventional definition in the stationary case since the term J c α E ( t ) / 2 vanishes after time averaging [53].
However, [80,81] propose not to use any heat flux definition due to apparent contradictions with the third law of thermodynamics and the onset of thermodynamic inconsistencies. We would like to clarify and highlight that those references consider a very different setup, in which the time dependence of the Hamiltonian also takes place in the contact region (in our case only the central part of the system evolves in time). This might seem as a minor difference, but on the contrary it is a fundamental point, because in that case there is a part of the energy stored in the contact region that should be interpreted as dissipative power and hence it will not contribute to the heat flux.
We emphasize that the extra term J c α E ( t ) / 2 has a pure dynamical origin, since it vanishes when it is averaged over a period. This makes sense from the physical point of view, because in a period the energy can be stored only at the reservoirs. On the other hand, this contribution of the contact region is also unique to the heat current. For example, the particle flux is independent of the charge variation in the coupling region, simply because it does no exist. In stark contrast, the tunneling term does contribute to the energy transfer between the quantum system and the bath and must be taken into account in any thermodynamically consistent calculation of time-dependent heat currents. Analogously, we can define the heat variation in the sample as
Q ˙ S ( t ) = E ˙ S ( t ) μ I S C ( t ) e + α J c α E ( t ) 2 .
such that the total heat satisfies Q ˙ tot ( t ) = α Q ˙ α ( t ) + Q ˙ S ( t ) , see Figure 2.

5. Entropy Production

Further insight can be gained if we examine the instantaneous entropy production, which is presumably a meaningful quantity for low driving frequencies, for which the system appears to be almost stationary and the classical laws of thermodynamics remain valid. Otherwise, for high enough frequencies, the nonequilibrium and quantum effects may be significant, which may hinder the use of the classical thermodynamics. This regime requires further investigation. The total entropy production has been demonstrated to follow an integral fluctuation theorem for a driven system strongly coupled to a heat bath [82]. Here, we stress that the entropy production associated with dissipative processes is always positive, in agreement with the second law. We begin by noticing that the AC forces have, in general, conservative (cons) and nonconservative (nc) components. In fact, in [33] it was shown that F ( t ) = F t + F nc ( t ) , where the first term corresponds to the Born-Oppenheimer component, which is evaluated with the equilibrium density matrix ρ t corresponding to the frozen Hamiltonian at the time t and F nc ( t ) collects the nonequilibrium components. The first term corresponds to an evolution given by a sequence of equilibrium states which are, thus, reversible and do not contribute to the entropy production. Hence P AC ( t ) = P cons ( t ) + P AC nc ( t ) . The electric power of the batteries is purely non-conservative (dissipative). In this way, the total power can be split into dissipative and conservative forces P tot ( t ) = P el ( t ) + P AC ( t ) = P tot diss ( t ) + P cons ( t ) , where P tot diss ( t ) = P el ( t ) + P AC nc ( t ) .
As a result, the total entropy production reads,
S ˙ diss ( t ) = 1 T α Q ˙ α ( t ) + Q ˙ S diss ( t ) = 1 T P tot diss ( t ) .
We have defined Q ˙ S diss ( t ) = Q ˙ S ( t ) P cons ( t ) , i.e., we identify the irreversible heat production at the central system as the result of subtracting the power developed by the conservative AC forces to the total heat instantaneously stored by the electrons that are in direct contact with the AC sources. After averaging over an AC cycle, the conservative part of the heat current vanishes. Furthermore, the heat variation in the sample also vanishes, Q ˙ S diss ¯ = Q ˙ S ¯ = 0 . The remaining terms have nonzero averages, and hence S ˙ diss ( t ) 0 at all times, in agreement with the second law of thermodynamics. This is explicitly shown in [53] within an adiabatic expansion, valid for low driving frequencies. Here, we analyze an example beyond the adiabatic regime and show that S ˙ diss ( t ) 0 even though for some times, the net heat flow at the reservoirs might be negative, α Q ˙ α ( t ) < 0 .

6. Calculation of the Currents within Green’s Function Approach for Non Interacting Systems

We now consider that the central part of the system is composed by a discrete chain of N sites (e.g., an array of quantum dots), to which many time-periodic potentials V i ( t ) are applied. Then, the Hamiltonian of the central piece reads:
H S ( t ) = H S 0 + m V m ( t ) d m d m ,
where
H S 0 = m = 1 N ε 0 m d m d m + m = 1 N 1 t m , m + 1 d m d m + 1 + h . c .
is the stationary part of the Hamiltonian, with ε 0 m the energy of the site m and t m , m + 1 the hopping between two consecutive sites.

6.1. Time Resolved Charge and Energy Currents Entering the Reservoirs

The charge and energy fluxes are defined from the expectation value of the time derivative of the number of particles and the energy of the reservoir, as in Equations (13) and (17). These currents are positive when entering the reservoir α and since [ H r e s , N α ] = [ H S , N α ] = [ H r e s , H α ] = [ H S , H α ] = 0 , they can be written as:
I α ( t ) = i e [ H c α , N α ] = 2 e Re k α w k α G l α , k α < ( t , t ) ,
and
J α E ( t ) = i [ H c α , H α ] = 2 Re k α ε k α w k α G l α , k α < ( t , t ) ,
where G l α , k α < ( t , t ) i c k α ( t ) d l α ( t ) is a Green’s function that involves operators of the reservoir c k α , as well as the central part, d l α . From the Dyson equation and Langreth rules [61], Green’s function above can be expressed as follows
G l α , k α < ( t , t ) = d t 1 w k α * [ G l α , l α R ( t , t 1 ) g k α < ( t 1 t ) + G l α , l α < ( t , t 1 ) g k α A ( t 1 t ) ] ,
with G m , n R ( t , t ) = i θ ( t t ) { d m ( t ) d n ( t ) } and G m , n < ( t , t ) = i { d n ( t ) d m ( t ) } being the retarded and the lesser Green’s functions of the central piece, and
g k α < ( t t ) = i d ε 2 π f α ( ε ) ρ k α ( ε ) e i ε ( t t ) / ,
g k α A ( t t ) = d ε 2 π d ε 2 π ρ k α ( ε ) ε ε i 0 + e i ε ( t t ) / ,
the Green’s functions of the uncoupled reservoirs, with ρ k α ( ε ) = 2 π δ ( ε ε k α ) , and f α ( ε ) = [ e ( ε μ α ) / ( k B T ) + 1 ] 1 the Fermi-Dirac distribution. Then, substituting Equations (30)–(32) in Equation (28), the charge current entering the reservoir α can be written as
I α ( t ) = 2 e h Re d t 1 d ε e i ε ( t 1 t ) / i G l α , l α R ( t , t 1 ) f α ( ε ) Γ α ( ε ) + G l α , l α < ( t , t 1 ) d ε 2 π Γ α ( ε ) ( ε ε i 0 + ) ,
where we have defined the hybridization as Γ α = k α | w k α | 2 ρ k α . We now consider the Dyson equation for the lesser Green’s function [68]
G n , m < ( t , t ) = α d t 1 d t 2 G n , l α R ( t , t 1 ) Σ α < ( t 1 t 2 ) [ G m , l α R ( t , t 2 ) ] * ,
with Σ α < ( ε ) = i f α ( ε ) Γ α ( ε ) , and we introduce the Floquet-Fourier representation [66], which is convenient for periodic fields:
G ^ R ( t , t ) = d ε 2 π G ^ R ( t , ε ) e i ε ( t t ) / ,
with
G ^ R ( t , ε ) = n = G ^ ( n , ε ) e i n ω t ,
where ω = 2 π / τ is the oscillation frequency of the AC parameters V ( t ) , and G ^ ( n , ε ) are the Floquet components. Green’s function G ^ R ( t , ε ) , can be obtained by solving the Dyson equation [66,67]:
G ^ R ( t , ε ) = G ^ 0 ( ε ) + n = G ^ R ( t , ε + n ω ) V ^ ( n ) G ^ 0 ( ε ) e i n ω t ,
with G ^ 0 ( ε ) = [ ε I ^ H ^ S 0 i Γ ^ / 2 ] 1 being the stationary retarded Green’s function and Γ ^ m , n = α = L , R δ m , n δ l α , m Γ α . The elements of the matrix V ^ l , l ( n ) = δ l , l i δ l , i V i ( n ) are the Fourier components of V ^ l , l ( t ) = n V ^ l , l ( n ) e i n ω t .
Replacing Equations (35) and (36) into Equation (33), we get
I α ( t ) = 2 e h Re l e i l ω t d ε i G l α , l α ( l , ε ) f α ( ε ) Γ α ( ε ) + β = L , R n i G l α , l β ( l + n , ε ) G l α , l β * ( n , ε ) f β ( ε ) Γ β ( ε ) d ε 2 π Γ α ( ε ) ( ε ( ε n ω ) i 0 + ) .
It is convenient to express the charge current in terms of differences between the Fermi functions of the reservoirs. To this end, we first explicitly write the real part as the sum of a complex number and its conjugate
I α ( t ) = e h d ε l i e i l ω t f α ( ε ) Γ α ( ε ) G l α , l α ( l , ε ) G l α , l α * ( l , ε ) + β = L , R n f β ( ε ) Γ β ( ε ) G l α , l β ( l + n , ε ) G l α , l β * ( l , ε ) d ε 2 π Γ α ( ε ) ε ( ε n ω ) i 0 + Γ α ( ε + l ω ) ε ( ε n ω ) + i 0 + .
Within the wide band limit, the hybridization can be considered a constant function Γ α ( ε ) Γ α . Taking into account also the following relation
1 ε ε i 0 + = P { 1 ε ε } + i π δ ( ε ε ) ,
we find
I α ( t ) = e h d ε l i e i l ω t Γ α f α ( ε ) G l α , l α ( l , ε ) G l α , l α * ( l , ε ) + i β = L , R n f β ( ε n ω ) Γ β G l α , l β ( l + n , ε n ω ) G l α , l β * ( n , ω n ω ) .
Finally, using the identity for the spectral components [68]
i G l , l ( n , ε ) i G l , l ( n , ε + n ω ) = β = L , R n G l , l β ( n + n , ε n ω ) Γ β G l , l β * ( n , ε n ω ) ,
the charge current reads
I α ( t ) = e h d ε l e i l ω t Γ α i G l α , l α * ( l , ε ) f α ( ε l ω ) f α ( ε ) β = L , R n Γ β G l α , l β ( l + n , ε n ω ) G l α , l β * ( n , ε n ω ) f β ( ε n ω ) f α ( ε ) } .
Following the same procedure, we can express the energy current entering reservoir α in Equation (29) as
J α E ( t ) = 2 e h Re d t 1 d ε e i ε ( t 1 t ) / i G l α , l α R ( t , t 1 ) f α ( ε ) Γ α ( ε ) ε + G l α , l α < ( t , t 1 ) d ε 2 π Γ α ( ε ) ε ( ε ε i 0 + ) ,
where we have used inside the integral that it is possible to make the replacement ε k α δ ( ε ε k α ) ε δ ( ε ε k α ) . Now, using the Floquet representation for Green’s functions and repeating the same steps as for the charge current, we find
J α E ( t ) = e h d ε l e i l ω t Γ α i G l α , l α * ( l , ε ) ( ε l ω ) f α ( ε l ω ) ε f α ( ε ) n , β = L , R n Γ β G l α , l β ( l + n , ε n ω ) G l α , l β * ( n , ε n ω ) ( ε l ω 2 ) f β ( ε n ω ) ε f α ( ε ) .

6.2. Energy Stored in the Contact Regions

The rate of change for the energy in the contact regions J c α E ( t ) = H ˙ c α can be also expressed in terms of Green’s function G l α , k α < ( t , t )
J c α E ( t ) = k α w k α d c k α ( t ) d l α ( t ) d t + h . c = 2 k α Im w k α d G l α , k α < ( t , t ) d t .
Following a similar procedure for the charge and energy currents, we find that the energy stored in the contact region reads
J c α E ( t ) = d ε 2 π ω Γ α l l 2 Im { e i l ω t G l α , l α ( l , ε ) } f α ( ε ) + n β = L , R f β ( ε ) Γ β Re { e i l ω t G l α , l β ( l + n , ε ) G l α , l β * ( n , ε ) } .
It is easy to show that the second term of the above equation vanishes, since:
l > 0 n l Re { e i l ω t G l α , l β ( l + n , ε ) G l α , l β * ( n , ε ) } l Re { e i l ω t G l α , l β ( l + n , ε ) G l α , l β * ( n , ε ) } = l > 0 n l Re { e i l ω t G l α , l β ( l + n , ε ) G l α , l β * ( n , ε ) } Re { e i l ω t G l α , l β * ( l + n , ε ) G l α , l β ( n , ε ) } = 0 .
Hence,
J c α E ( t ) = d ε h f α ( ε ) Γ α l l ω 2 Im { e i l ω t G l α , l α ( l , ε ) } .
Combining this expression with Equation (36), we find the expression:
J c α E ( t ) = d ε 2 π f α ( ε ) Γ α 2 Re { t G l α , l α R ( t , ε ) } ,
which clearly shows that this energy reactance term is of purely an AC nature since it vanishes after time averaging.

6.3. Power Developed by the AC Sources

The power performed by the AC potential can be calculated as P A C ( t ) = H S t , which is
P A C ( t ) = i d V i ( t ) d t d i ( t ) d i ( t ) = i d V i ( t ) d t Im { G i , i < ( t , t ) } .
Now, replacing Equation (34) in the above equation , and using the Floquet representation, we get
P A C ( t ) = i α = L , R l , m , n d ε h n ω f α ( ε ) Γ α Im { V i ( n ) G i , l α ( m + l , ε ) G i , l α * ( l , ε ) e i ω t ( m n ) } .

7. Relation to the Scattering Matrix Formalism

In this section, we show that the definition of the heat current in Equation (23), which includes the contribution J c α E due to the contact α, is fully in agreement with the scattering matrix formalism. From Equations (43), (45) and (48), the heat reads
Q ˙ α ( t ) = l d ε h e i l ω t Γ α i G l α , l α * ( l , ε ) ( ε l ω 2 μ α ) ( f α ( ε ) f α ( ε l ω ) ) n β = L , R ( ε + l ω 2 μ α ) ( f α ( ε ) f β ( ε n ω ) ) Γ β G l α , l β ( l + n , ε n ω ) G l α , l β * ( n , ε n ω )
Within the scattering matrix approach, the heat flux entering the reservoir α reads [75,83]
Q ˙ α S ( t ) = l , n e i l ω t d ε h ( ε + l ω 2 μ α ) β = L , R ( f β ( ε n ) f α ( ε ) ) S α β * ( ε , ε n ) S α , β ( ε l , ε n ) ,
where ε n = ε + n ω and S ( ε n , ε l ) is the Floquet scattering matrix which is related to the Green function via the generalized Fisher-Lee relation [68]
S α β ( ε m , ε n ) = δ α , β δ m , n i Γ α ( ε m ) Γ β ( ε n ) G l α , l β ( m n , ε n ) .
Using this relation, we find that
Q ˙ α S ( t ) = l , n d ε h e i l ω t ( ε + l ω 2 μ α ) β = L , R ( f β ( ε n ω ) f α ( ε ) ) × G l α , l β * ( n , ε n ω ) i δ α β δ l , n Γ α Γ β + Γ α Γ β G l α , l β ( l + n , ε n ω ) ,
Then, after some algebra and by comparing with Equation (52), we have Q ˙ α S ( t ) = Q ˙ α ( t ) .

8. Low Frequency Expansion

For low frequencies of the driving potentials, a solution of the Dyson equation (37) up to O ( ω ) can be obtained by expanding in powers of the AC frequency as
G ^ R ( t , ε ) G ^ 0 ( ε ) + G ^ R ( t , ε ) V ^ ( t ) G ^ 0 ( ε ) + i ε G ^ R ( t , ε ) d V ^ ( t ) d t G ^ 0 ( ε ) .
Now, we can define the frozen Green’s function
G ^ f R ( t , ε ) = G ^ 0 ( ε ) 1 V ^ ( t ) 1 ,
in terms of which the exact solution of Equation (37) up to first order in ω reads
G ^ R ( t , ε ) = G ^ f R ( t , ε ) + i 2 2 t ε G ^ f R ( t , ε ) + A ^ ( t , ε ) ,
with
A ^ ( t , ε ) = ε G ^ f R ( t , ε ) V ^ ( t ) G ^ f R ( t , ε ) G ^ f R ( t , ε ) V ^ ( t ) ε G ^ f R ( t , ε ) .
Note that this expression is valid to lowest order in ω, but can be used for arbitrary values of the amplitude V A C [28,68]. This is specially appealing if one is interested in addressing nonlinear effects [65,84,85]. The validity of the adiabatic approximation considered here was confirmed in [86] by a comparison with an exact calculation in the limit of low level occupations.
We can also expand the Floquet components G ^ ( n , ε ) up to first order in ω as
G ^ ( n , ε ) G ^ ( 0 ) ( n , ε ) + ω G ^ ( 1 ) ( n , ε ) ,
Then, using the representation of Equation (36) in Equation (58), it is possible to identify the first and second terms of the expansion:
G ^ ( 0 ) ( n , ε ) = 0 τ d t τ G ^ f R ( t , ε ) e i n ω t ω G ^ ( 1 ) ( n , ε ) = 0 τ d t τ i 2 2 t ε G ^ f R ( t , ε ) + A ^ ( t , ε ) e i n ω t .

9. Application

We now illustrate our main ideas with a simple but generic example: a single driven level coupled to two fermionic baths kept at the same temperature, but with different chemical potentials, μ L = μ and μ R = μ δ μ . The Hamiltonian of the sample is then
H S = ε d ( t ) d d ,
where the level is periodically driven with a monochromatic potential of frequency ω and amplitude V AC , ε d ( t ) = ε 0 + V ( t ) with V ( t ) = V AC cos ω t . The response of the system will be periodic. Another possibility is to explore heat dynamics in the transient regime, usually in quantum levels whose potential is suddenly changed by an externally-applied voltage pulse. Thermoelectric efficiency can thus be enhanced by a suitable time-dependent potential applied to a quantum dot level far from equilibrium [32,87]. Based on a nonequilibrium Green’s function approach, [88] investigates the transient heat current using a complex absorbing potential, greatly reducing the complex calculation of the heat current through a molecule connected to metallic leads and under a gate voltage pulse. Indeed, the transient heat current of a quantum dot in response to a voltage shift can reveal signatures of electronic interaction despite the inherent Coulomb repulsion [89]. The full counting statistics of the quantum energy flow can also be determined for the transient regime, showing universal scaling laws in the single quantum dot case [90]. When the system consists of two levels coupled to a thermal bath, the transient dynamics is sensitive to the employed theoretical formalism (master equation approach or functional influence method) [71].
Here, we restrict ourselves to mean currents. However, their fluctuations can provide us with useful information. For instance, fluctuations of heat beyond the level dictated by the Callen-Welton fluctuation-dissipation theorem signal the presence of backscattering [41]. Furthermore, the energy noise reveals additional information not contained in the charge current fluctuations [42]. The novel features are related to electron-hole correlations arising in phase-coherent small conductors driven by periodic voltage biases. The mixed noise between correlations of charge and energy currents shows both interference and transport contributions, and is unique as compared with either charge or energy noises [43,44,45,46,47].

9.1. Adiabatic Regime and Linear Response in the Bias Voltage

For small bias voltage, it is possible to combine the procedure described in Section 8 with an expansion in δ μ . To calculate up to second order in ω and δ μ the heat current of Equation (52), the AC power in Equation (51), and the power developed by the batteries in Equation (20), we need to perform an expansion of the Fermi functions entering the integrals as:
f α ( ε + n ω ) f α ( ε ) + ε f α n ω + ε 2 f α ( n ω ) 2 2 ,
with
f L ( ε ) f ( ε )
f R ( ε ) f ( ε ) + ε f δ μ + ε 2 f δ μ 2 2 .
Accordingly, we can compute all of the fluxes entering Equation (25). Then, up to second order in ω and δ μ , the total dissipative power P tot diss reads:
P tot diss ( t ) = L δ μ 2 ( t ) δ μ 2 + L ω 2 ( t ) ω 2 ,
where the transport coefficients are:
L δ μ 2 ( t ) = d ε h ε f Γ L Γ R | G f R ( t , ε ) | 2 L ω 2 ( t ) = 2 ω 2 d ε 2 π ε f Γ 2 | t G f R ( t , ε ) | 2 ,
where Γ = Γ L + Γ R is the total hybridization. Notice that, as is shown is [33] for systems with time-reversal symmetry, the terms that are proportional to ω δ μ do not contribute to the total dissipated power.
Following a similar procedure in Equation (52) and then subtracting Equation (66), we find that the dissipated heat flux of the central piece is
Q ˙ S diss ( t ) = L d ω ( t ) ω + L d ω δ μ ( t ) ω δ μ + L d ω 2 ( t ) ω 2 ,
with
L d ω ( t ) = d ε 2 π ε f ( ε μ ) Γ | G f R ( t , ε ) | 2 V ac sin ( ω t ) , L d ω δ μ ( t ) = d ε 2 π Γ R ε f ( ε μ ) ε | G f R ( t , ε ) | 2 V AC sin ( ω t ) , L d ω 2 ( t ) = 2 ω 2 d ε 2 π Γ 2 ε f ( ε μ ) t ( ε G f R ( t , ε ) t G f R * ( t , ε ) ) .
In Figure 3 (upper panel) we show the total heat generated in the reservoirs as a function of time for zero temperature (red dashed line) and nonzero temperature (blue line). In the latter case, an individual heat flux can become negative but the total dissipated heat must be positive for all temperatures, as shown in Figure 3 (lower panel). Note that one should also take into account the contribution from the sample, which is purely dynamical (Figure 3 (middle panel)) and vanishes for zero temperature.

9.2. Nonadiabatic Regime

We now consider the entropy production for high driving frequencies, where the linear response approach of Section 9.1 is not applicable. The evaluation of the currents implies the evaluation of the exact Dyson equation (37) following the methods of [66,67]. A similar situation was analyzed in [91], where it was found that the rate of change for the nonadiabatic energy agrees with the work done on a molecule subjected to a time-dependent applied field. The work in [92] discusses heat transfer across a triple dot structure with the nonadiabatic driving applied to an edge dot.
In our case, we focus on reservoirs at zero temperature for simplicity. We see in the figure that, unlike the adiabatic case, the heat flux at the reservoirs may attain negative values for short (but finite) time intervals (see Figure 4 (upper panel)). Nevertheless, we emphasize that the total entropy production is always positive at every instant (see Figure 4 (lower panel)). For this to occur the dissipated heat at the sample must be nonvanishing and positive (see Figure 4 (middle panel)).

10. Conclusions

The subject of time-dependent heat in quantum electronic systems presents unique aspects as compared with the better understood macroscopic case. In this work, we have briefly reviewed the recent developments in the field. On the basis of general thermodynamic arguments and conservation laws we have presented exact equations for the time-resolved energy conversion in a quantum system in contact with AC driving sources and reservoirs at which DC bias voltages are applied. These equations lead to the unambiguous definition of the total heat and production in the full system containing the central driven device, the reservoirs and the contacts between the different parts. These results are independent of the particular method used to evaluate the physical quantities involved. In order to get the exact solution of an example we have used, for convenience, a Green function technique as in [52,53]. Of course, the physical picture is beyond the method employed to solve the particular problem.
We have also analyzed the possibility of defining the heat flows not only resolved in time but also in space, in the sense of identifying flows in the reservoirs and through the piece that is in contact with the AC sources. How to properly perform this split in a physically-meaningful way is a non-trivial task. Intuitively, we expect that deep in the reservoirs, we should have an equilibrium system. For adiabatically-driven systems, we have found that this seems to be, indeed the case. Then, by suitably generalizing the definition of the time-dependent heat current. We argue that the proper generalization is given by Equation (23), which implies taking into account the rate of change of the energy stored at the reservoir plus half of the one stored at the contact. Such a definition leads to results compatible with scattering matrix expressions for continuum models, and to the derivation of a Joule law when the reservoirs are at zero temperature.
However, for highly non-equilibrium driving at the central piece and strongly-coupled reservoirs, it is not obvious at all if we should ever recover an equilibrium state at every instant, even deep inside the reservoirs. In fact, we have analyzed the example of a single-level system driven by a high-frequency gate voltage in contact with reservoirs at zero temperature and we have found that the definition of the instantaneous heat current at the reservoirs (Equation (23)) that leads to meaningful physical results in the adiabatic regime may be negative at some instants, meaning a temporary flow of heat exiting the reservoirs. This result is accompanied by a finite value of what we have identified as the heat flow in the central piece. The sum of both contributions leads to a net instantaneous positive value of the total entropy production. Hence, the result is in agreement with the second law, albeit non-intuitive. This is because the natural expectation is to have a dissipative energy flow only into the reservoirs and not in the central device. This is, in fact, the case when we consider the average over one period of the driving. However, instantaneously, and for non-equilibrium situations, the exact equation for the entropy production is given by Equation (25), which contains contributions associated to the reservoirs, contacts and the central device. The second law implies that the sum of all these contributions should be positive but no fundamental law imposes independent conditions on the individual terms. Finally, we have commented on the difficulties of accessing experimentally a reliable test of the time- and spatially-resolved heat transport discussed in the present work.

Acknowledgments

This work was supported by MINECO under Grant No. FIS2014-52564, UBACyT, CONICET and MINCyT, Argentina.

Author Contributions

All the authors conceived the theoretical model and performed the calculations. They analyzed the results and wrote the paper. All the authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Benenti, G.; Casati, G.; Saito, K.; Whitney, R.S. Fundamental aspects of steady-state conversion of heat to work at the nanoscale. arXiv, 2016; arXiv:1608.05595. [Google Scholar]
  2. Kosloff, R. Quantum thermodynamics: A dynamical viewpoint. Entropy 2013, 15, 2100–2128. [Google Scholar] [CrossRef]
  3. Campisi, M.; Fazio, R.J. Dissipation, Correlation and Lags in Heat Engines. Phys. A Math. Theor. 2016, 49, 345002. [Google Scholar] [CrossRef]
  4. Gasparinetti, S.; Solinas, P.; Pugnetti, S.; Fazio, R.; Pekola, J.P. Environment-governed dynamics in driven quantum systems. Phys. Rev. Lett. 2013, 110, 150403. [Google Scholar] [CrossRef] [PubMed]
  5. Segal, D. Heat transfer in the spin-boson model: A comparative study in the incoherent tunneling regime. Phys. Rev. E 2014, 90, 012148. [Google Scholar] [CrossRef] [PubMed]
  6. Carrega, M.; Solinas, P.; Sassetti, M.; Weiss, U. Energy Exchange in Driven Open Quantum Systems at Strong Coupling. Phys. Rev. Lett. 2016, 116, 240403. [Google Scholar] [CrossRef] [PubMed]
  7. Schmidt, R.; Carusela, M.F.; Pekola, J.P.; Suomela, S.; Ankerhold, J. Work and heat for two-level systems in dissipative environments: Strong driving and non-Markovian dynamics. Phys. Rev. B 2015, 91, 224303. [Google Scholar] [CrossRef]
  8. Ponomarev, A.V.; Denisov, S.; Hänggi, P. AC-Driven Atomic Quantum Motor. Phys. Rev. Lett. 2009, 102, 230601. [Google Scholar] [CrossRef] [PubMed]
  9. Abah, O.; Roßnagel, J.; Jacob, G.; Deffner, S.; Schmidt- Kaler, F.; Singer, K.; Lutz, E. Single-ion heat engine at maximum power. Phys. Rev. Lett. 2012, 109, 203006. [Google Scholar] [CrossRef] [PubMed]
  10. Brantut, J.-P.; Grenier, C.; Meineke, J.; Stadler, D.; Krinner, S.; Kollath, C.; Esslinger, T.; Georges, A. A thermoelectric heat engine with ultracold atoms. Science 2013, 342, 713–715. [Google Scholar] [CrossRef] [PubMed]
  11. Uzdin, R.; Levy, A.; Kosloff, R. Equivalence of quan- tum heat machines, and quantum-thermodynamic signatures. Phys. Rev. X 2015, 5, 031044. [Google Scholar]
  12. Blickle, V.; Bechinger, C. Realization of a micrometer-sized stochastic heat engine. Nat. Phys. 2011, 8, 143–146. [Google Scholar] [CrossRef]
  13. Martínez, I.A.; Roldán, E.; Dinis, L.; Petrov, D.; Rica, R.A. Adiabatic Processes Realized with a Trapped Brownian Particle. Phys. Rev. Lett. 2015, 114, 120601. [Google Scholar] [CrossRef] [PubMed]
  14. Scully, M.O. Quantum photocell: Using quantum co-herence to reduce radiative recombination and increase efficiency. Phys. Rev. Lett. 2010, 104, 207701. [Google Scholar] [CrossRef] [PubMed]
  15. Steeneken, P.G.; Le Phan, K.; Goossens, M.J.; Koops, G.E.J.; Brom, G.J.A.M.; van der Avoort, C.; van Beek, J.T.M. Piezoresistive heat engine and refrigerator. Nat. Phys. 2010, 7, 354–359. [Google Scholar] [CrossRef]
  16. Scully, M.O.; Zubairy, M.S.; Agarwal, G.S.; Walther, H. Extracting work from a single heat bath via vanishing quantum coherence. Science 2003, 299, 862–864. [Google Scholar] [CrossRef] [PubMed]
  17. Chamon, C.; Mucciolo, E.R.; Arrachea, L.; Capaz, R.B. Heat Pumping in Nanomechanical Systems. Phys. Rev. Lett. 2011, 106, 135504. [Google Scholar] [CrossRef] [PubMed]
  18. Arrachea, L.; Mucciolo, E.R.; Chamon, C.; Capaz, R.B. Microscopic model of a phononic refrigerator. Phys. Rev. B 2012, 86, 125424. [Google Scholar] [CrossRef]
  19. Steele, G.A.; Hẗtel, A.K.; Witkamp, B.; Poot, M.; Meerwaldt, H.B.; Kouwenhoven, L.P.; van der Zant, H.S.J. Strong coupling between single-electron tunneling and nanomechanical motion. Science 2009, 325, 1103–1107. [Google Scholar] [CrossRef] [PubMed]
  20. Lotze, C.; Corso, M.; Franke, K.J.; von Oppen, F.; Pascual, J.I. Driving a macroscopic oscillator with the stochastic motion of a hydrogen molecule. Science 2012, 338, 779–782. [Google Scholar] [CrossRef] [PubMed]
  21. Brandner, K.; Saito, K.; Seifert, U. Thermodynamics of micro- and nano-systems driven by periodic temperature variations. Phys. Rev. X 2015, 5, 031019. [Google Scholar] [CrossRef]
  22. Horowitz, J.M. Quantum-trajectory approach to the stochastic thermodynamics of a forced harmonic oscillator. Phys. Rev. E 2012, 85, 031110. [Google Scholar] [CrossRef] [PubMed]
  23. Martinez, E.A.; Paz, J.P. Dynamics and thermodynamics of linear quantum open systems. Phys. Rev. Lett. 2013, 110, 130406. [Google Scholar] [CrossRef] [PubMed]
  24. Freitas, N.; Paz, J.P. Fundamental limits for cooling of linear quantum refrigerators. arXiv, 2016; arXiv:1607.04234. [Google Scholar]
  25. Arrachea, L.; Moskalets, M.; Martin-Moreno, L. Heat production and energy balance in nanoscale engines driven by time-dependent fields. Phys. Rev. B 2007, 75, 245420. [Google Scholar] [CrossRef]
  26. Rey, M.; Strass, M.; Kohler, S.; Hänggi, P.; Sols, F. Nonadiabatic electron heat pump. Phys. Rev. B 2007, 76, 08533. [Google Scholar] [CrossRef]
  27. Juergens, S.; Haupt, F.; Moskalets, M.; Splettstoesser, J. Thermoelectric performance of a driven double quantum dot. Phys. Rev. B 2013, 87, 245423. [Google Scholar] [CrossRef]
  28. Moskalets, M.; Büttiker, M. Heat production and current noise for single- and double-cavity quantum capacitors. Phys. Rev. B 2009, 80, 081302. [Google Scholar] [CrossRef]
  29. Bustos-Marun, R.; Refael, G.; von Oppen, F. Adiabatic Quantum Motors. Phys. Rev. Lett. 2013, 111, 060802. [Google Scholar] [CrossRef] [PubMed]
  30. Arrachea, L.; von Oppen, F. Nanomagnet coupled to quantum spin Hall edge: An adiabatic quantum motor. Physica E 2015, 74, 596–602. [Google Scholar] [CrossRef]
  31. Daré, A.M.; Lombardo, P. Time-dependent thermoelectric transport for nanoscale thermal machines. Phys. Rev. B 2016, 93, 035303. [Google Scholar] [CrossRef]
  32. Zhou, H.; Thingna, J.; Hänggi, P.; Wang, J.-S.; Li, B. Boosting thermoelectric efficiency using time-dependent control. Sci. Rep. 2015, 5, 14870. [Google Scholar] [CrossRef] [PubMed]
  33. Ludovico, M.F.; Battista, F.; von Oppen, F.; Arrachea, L. Adiabatic response and quantum thermoelectrics for AC-driven quantum systems. Phys. Rev. B 2016, 93, 075136. [Google Scholar] [CrossRef]
  34. Bergenfeldt, C.; Samuelsson, P.; Sothmann, B.; Flindt, C.; Büttiker, M. Hybrid microwave-cavity heat engine. Phys. Rev. Lett. 2014, 112, 076803. [Google Scholar] [CrossRef] [PubMed]
  35. Biele, R.; D’Agosta, R.; Rubio, A. Time-dependent thermal transport theory. Phys. Rev. Lett. 2015, 115, 056801. [Google Scholar] [CrossRef] [PubMed]
  36. Alicki, R.; Gelbwaser-Klimovsky, D.; Szczygielski, K. Solar cell as self-oscillating heat engine. J. Phys. A Math. Theor. 2016, 49, 015002. [Google Scholar] [CrossRef]
  37. Alicki, R. Thermoelectric generators as self-oscillating heat engines. J. Phys. A Math. Theor. 2016, 49, 085001. [Google Scholar] [CrossRef]
  38. Kindermann, M.; Pilgram, S. Statistics of heat transfer in mesoscopic circuits. Phys. Rev. B 2004, 69, 155334. [Google Scholar] [CrossRef]
  39. Averin, D.V.; Pekola, J.P. Violation of the Fluctuation-Dissipation Theorem in Time-Dependent Mesoscopic Heat Transport. Phys. Rev. Lett. 2010, 104, 220601. [Google Scholar] [CrossRef] [PubMed]
  40. Sergi, D. Energy transport and fluctuations in small conductors. Phys. Rev. B 2011, 83, 033401. [Google Scholar] [CrossRef]
  41. Moskalets, M. Floquet scattering matrix theory of heat fluctuations in dynamical quantum conductors. Phys. Rev. Lett. 2014, 112, 206801. [Google Scholar] [CrossRef]
  42. Battista, F.; Haupt, F.; Splettstoesser, J. Energy and power fluctuations in AC-driven coherent conductors. Phys. Rev. B 2014, 90, 085418. [Google Scholar] [CrossRef]
  43. Sánchez, R.; Sothmann, B.; Jordan, A.N.; Büttiker, M. Correlations of heat and charge currents in quantum-dot thermoelectric engines. New J. Phys. 2013, 15, 125001. [Google Scholar] [CrossRef]
  44. Battista, F.; Haupt, F.; Splettstoesser, J. Correlations between charge and energy current in AC-driven coherent conductors. J. Phys. Conf. Ser. 2014, 568, 052008. [Google Scholar] [CrossRef]
  45. Crépieux, A.; Michelini, F. Mixed, charge and heat noises in thermoeleectric nanosystems. J. Phys. Condens. Matter 2015, 27, 015302. [Google Scholar] [CrossRef] [PubMed]
  46. Crépieux, A.; Michelini, F. Heat-charge mixed noise and thermoelectric efficiency fluctuations. J. Stat. Mech. arXiv, 2016; arXiv:1510.01056. [Google Scholar]
  47. Eyméoud, P.; Crépieux, A. Mixed electrical-heat noise spectrum in a quantum dot. arXiv, 2016; arXiv:1608.02414. [Google Scholar]
  48. Campisi, M.; Pekola, J.; Fazio, R. Nonequilibrium fluctuations in quantum heat engines: Theory, example, and possible. New J. Phys. 2015, 17, 035012. [Google Scholar] [CrossRef]
  49. Caso, A.; Arrachea, L.; Lozano, G.S. Local temperatures and heat flow in quantum driven systems. Phys. Rev. B 2011, 83, 165419. [Google Scholar] [CrossRef]
  50. Arrachea, L.; Cugliandolo, L.F. Study of a fluctuation-dissipation relation of a dissipative driven mesoscopic system. EPL Europhys. Lett. 2005, 70, 642–648. [Google Scholar] [CrossRef]
  51. Caso, A.; Arrachea, L.; Lozano, G.S. Local and effective temperatures of quantum driven systems. Phys. Rev. B 2010, 81, 041301. [Google Scholar] [CrossRef]
  52. Ludovico, M.F.; Lim, J.S.; Moskalets, M.; Arrachea, L.; Sánchez, D. Dynamical energy transfer in AC-driven quantum systems. Phys. Rev. B 2014, 89, 161306. [Google Scholar] [CrossRef]
  53. Ludovico, M.F.; Moskalets, M.; Sánchez, D.; Arrachea, L. Dynamics of energy transport and entropy production in AC-driven quantum electron systems. Phys. Rev. B 2016, 94, 035436. [Google Scholar] [CrossRef]
  54. Moskalets, M.; Haack, G. Heat and charge transport measurements to access single-electron quantum characteristics. arXiv, 2016; arXiv:1609.04544. [Google Scholar]
  55. Ludovico, M.F.; Lim, J.S.; Moskalets, M.; Arrachea, L.; Sánchez, D. Time resolved heat exchange in driven quantum systems. J. Phys. Conf. Ser. 2014, 568, 052017. [Google Scholar] [CrossRef]
  56. Büttiker, M.; Prêtre, A.; Thomas, H. Mesoscopic capacitors. Phys. Lett. A 1993, 180, 364–369. [Google Scholar] [CrossRef]
  57. Brouwer, P.W. Scattering approach to parametric pumping. Phys. Rev. B 1998, 58, 10135. [Google Scholar] [CrossRef]
  58. Avron, J.E.; Elgart, A.; Graf, G.M.; Sadum, L. Optimal Quantum Pumps. Phys. Rev. B 2001, 87, 236601. [Google Scholar] [CrossRef] [PubMed]
  59. Vavilov, M.G.; Ambegoakar, V.; Aleiner, I.L. Charge pumping and photovoltaic effect in open quantum dots. Phys. Rev. B 2001, 63, 195313. [Google Scholar] [CrossRef]
  60. Pastawski, H.M. Classical and quantum transport from generalized Landauer-Bẗtiker equations II. Time-dependent resonant tunneling. Phys. Rev. B 1992, 63, 4053–4070. [Google Scholar] [CrossRef]
  61. Jauho, A.-P.; Wingreen, N.S.; Meir, Y. Time-dependent transport in interacting and noninteracting resonant-tunneling systems. Phys. Rev. B 1994, 50, 5528. [Google Scholar] [CrossRef] [Green Version]
  62. Platero, G.; Aguado, R. Photon-assisted transport in semiconductor nanostructures. Phys. Rep. 2004, 395, 1–157. [Google Scholar] [CrossRef]
  63. Kohler, S.; Lehmann, J.; Hänggi, P. Driven quantum transport on the nanoscale. Phys. Rep. 2005, 406, 379–443. [Google Scholar] [CrossRef]
  64. Moskalets, M. Scattering Matrix Approach to Non-Stationary Quantum Transport; Imperial College Press: London, UK, 2011. [Google Scholar]
  65. Moskalets, M.; Samuelsson, P.; Büttiker, M. Quantized Dynamics of a Coherent Capacitor. Phys. Rev. B 2008, 100, 086601. [Google Scholar] [CrossRef] [PubMed]
  66. Arrachea, L. Green-function approach to transport phenomena in quantum pumps. Phys. Rev. B 2005, 72, 125349. [Google Scholar] [CrossRef]
  67. Arrachea, L. Exact Green’s function renormalization approach to spectral properties of open quantum systems driven by harmonically time-dependent fields. Phys. Rev. B 2007, 75, 035319. [Google Scholar] [CrossRef]
  68. Arrachea, L.; Moskalets, M. Relation between scattering-matrix and Keldysh formalisms for quantum transport driven by time-periodic fields. Phys. Rev. B 2006, 74, 245322. [Google Scholar] [CrossRef]
  69. Bruch, A.; Thomas, M.; Kusminskiy, S.V.; von Oppen, F.; Nitzan, A. Quantum thermodynamics of the driven resonant level model. Phys. Rev. B 2016, 93, 115318. [Google Scholar] [CrossRef]
  70. Ochoa, M.A.; Bruch, A.; Nitzan, A. Energy distribution and local fluctuations in strongly coupled open quantum systems: The extended resonant level model. Phys. Rev. B 2016, 94, 035420. [Google Scholar] [CrossRef]
  71. Carrega, M.; Solinas, P.; Braggio, A.; Sassetti, M. Computation of transient dynamics of energy power for a dissipative two state system. J. Stat. Mech. 2016, 054014. [Google Scholar] [CrossRef]
  72. Carrega, M.; Solinas, P.; Braggio, A.; Sassetti, M.; Weiss, U. Functional integral approach to time-dependent heat exchange in open quantum systems: General method and applications. New J. Phys. 2015, 17, 045030. [Google Scholar] [CrossRef]
  73. Büttiker, M.; Thomas, H.; Prêtre, A. Current partitioning in multiprobe conductors in the presence of slowly oscillating external potentials. Z. Phys. B 1994, 94, 133–137. [Google Scholar] [CrossRef]
  74. Chen, J.; ShangGuan, M.; Wang, J. A gauge invariant theory for time dependent heat current. New J. Phys. 2015, 17, 053034. [Google Scholar] [CrossRef]
  75. Lim, J.S.; López, R.; Sánchez, D. Dynamic thermoelectric and heat transport in mesoscopic capacitors. Phys. Rev. B 2013, 88, 201304. [Google Scholar] [CrossRef]
  76. Lv, Y.; Chen, J.; Zhen, R.-K.; Song, J.; Zhang, T.; Li, X.; Shi, X.; Chen, L. Photo-induced enhancement of the power factor of Cu2S thermoelectric films. Sci. Rep. 2015, 5, 16291. [Google Scholar] [CrossRef] [PubMed]
  77. Proesmans, K.; Cleuren, B.; Van den Broeck, C. Linear stochastic thermodynamics for periodically driven systems. J. Stat. Mech. 2016, 023202. [Google Scholar] [CrossRef]
  78. García-Cañadas, J.; Min, G. Thermal dynamics of thermoelectric phenomena from frequency resolved methods. AIP Adv. 2016, 6, 035008. [Google Scholar] [CrossRef]
  79. Rosselló, G.; Lim, J.S.; López, R. Time-dependent heat flow in interacting quantum conductors. Phys. Rev. B 2015, 92, 115402. [Google Scholar] [CrossRef]
  80. Esposito, M.; Ochoa, M.A.; Galperin, M. Quantum thermodynamics: A nonequilibrium Green’s function approach. Phys. Rev. Lett. 2015, 114, 080602. [Google Scholar] [CrossRef] [PubMed]
  81. Esposito, M.; Ochoa, M.A.; Galperin, M. Nature of heat in strongly coupled open quantum systems. Phys. Rev. B 2015, 92, 235440. [Google Scholar] [CrossRef]
  82. Seifert, U. First and second law of thermodynamics at strong coupling. Phys. Rev. Lett. 2016, 116, 020601. [Google Scholar] [CrossRef] [PubMed]
  83. Battista, F.; Moskalets, M.; Albert, M.; Samuelsson, P. Quantum heat fluctuations of single-particle sources. Phys. Rev. Lett. 2013, 110, 126602. [Google Scholar] [CrossRef] [PubMed]
  84. Kashuba, O.; Schoeller, H.; Splettstoesser, J. Nonlinear adiabatic response of interacting quantum dots. Europhys. Lett. 2012, 98, 57003. [Google Scholar] [CrossRef]
  85. Alomar, M.I.; Lim, J.S.; Sánchez, D. Coulomb-blockade effect in nonlinear mesoscopic capacitors. Phys. Rev B 2016, 94, 165425. [Google Scholar] [CrossRef]
  86. Nocera, A.; Perroni, C.A.; Ramaglia, V.M.; Cataudella, V. Charge and heat transport in soft nanosystems in the presence of time-dependent perturbations. Beilstein J. Nanotechnol. 2016, 7, 439–464. [Google Scholar] [CrossRef] [PubMed]
  87. Crépieux, A.; Simkovic, F.; Cambon, B.; Michelini, F. Enhanced thermopower under a time-dependent gate voltage. Phys. Rev. B 2011, 83, 153417. [Google Scholar] [CrossRef]
  88. Yu, Z.; Zhang, L.; Xing, Y.; Wang, J. Investigation of transient heat current from first principles using complex absorbing potential. Phys. Rev. B 2014, 90, 115428. [Google Scholar] [CrossRef]
  89. Schulenborg, J.; Saptsov, R.B.; Haupt, F.; Splettstoesser, J.; Wegewijs, M.R. Fermion-parity duality and energy relaxation in interacting open systems. Phys. Rev. B 2016, 93, 081411. [Google Scholar] [CrossRef]
  90. Yu, Z.; Tang, G.M.; Wang, J. Full-counting statistics of transient energy current in mesoscopic systems. Phys. Rev. B 2016, 93, 195419. [Google Scholar] [CrossRef]
  91. Mandal, A.; Hunt, K.L.C. Non-adiabatic current densities, transitions, and power absorbed by a molecule in a time-dependent electromagnetic field. J. Chem. Phys. 2015, 143, 034102. [Google Scholar] [CrossRef] [PubMed]
  92. Gallego-Marcos, F.; Platero, G. Coherent long-range thermoelectrics in nonadiabatic driven quantum systems. arXiv, 2016; arXiv:1607.06000. [Google Scholar]
Figure 1. Sketch of the system under consideration. A quantum conductor (described by the Hamiltonian H S ), is coupled to two reservoirs ( H L and H R ) kept at the same temperature T, but with different chemical potentials μ L and μ R . The conductor is also driven out of equilibrium by the application of AC local power sources, which are all collected in the vector V ( t ) . The Hamiltonians representing the left and right contact regions are H c L and H c R , respectively.
Figure 1. Sketch of the system under consideration. A quantum conductor (described by the Hamiltonian H S ), is coupled to two reservoirs ( H L and H R ) kept at the same temperature T, but with different chemical potentials μ L and μ R . The conductor is also driven out of equilibrium by the application of AC local power sources, which are all collected in the vector V ( t ) . The Hamiltonians representing the left and right contact regions are H c L and H c R , respectively.
Entropy 18 00419 g001
Figure 2. The total heat flux in the system is produced by the AC and DC driving sources, which is summarized by the relation Q ˙ tot ( t ) = P AC ( t ) P el ( t ) . The definition of the heat currents entering the reservoirs, as well as the heat variation in the central part, takes into account half of the energy stored in the contact regions. The heat flux of the central piece contributes purely dynamically, since it is zero when averaged over one period of the AC driving potentials. In the stationary state, charge and energy can be stored only at the reservoirs.
Figure 2. The total heat flux in the system is produced by the AC and DC driving sources, which is summarized by the relation Q ˙ tot ( t ) = P AC ( t ) P el ( t ) . The definition of the heat currents entering the reservoirs, as well as the heat variation in the central part, takes into account half of the energy stored in the contact regions. The heat flux of the central piece contributes purely dynamically, since it is zero when averaged over one period of the AC driving potentials. In the stationary state, charge and energy can be stored only at the reservoirs.
Entropy 18 00419 g002
Figure 3. Adiabatic regime of a single driven quantum level coupled to two fermionic baths. The different components of the total heat production Q tot diss ( t ) = P tot diss ( t ) as a function of time. Dashed lines correspond to reservoirs at T = 0 , while solid lines are for k B T = 0.02 . Energies are expressed in units of Γ L . The energy of the level evolves in time with V A C = 8 and ω = 1 × 10 3 . Parameters: μ = 1 , δ μ = 0.004 , ε 0 = 0 and hybridization widths Γ L = 1 and Γ R = 0.5 .
Figure 3. Adiabatic regime of a single driven quantum level coupled to two fermionic baths. The different components of the total heat production Q tot diss ( t ) = P tot diss ( t ) as a function of time. Dashed lines correspond to reservoirs at T = 0 , while solid lines are for k B T = 0.02 . Energies are expressed in units of Γ L . The energy of the level evolves in time with V A C = 8 and ω = 1 × 10 3 . Parameters: μ = 1 , δ μ = 0.004 , ε 0 = 0 and hybridization widths Γ L = 1 and Γ R = 0.5 .
Entropy 18 00419 g003
Figure 4. Nonadiabatic regime of a single driven quantum level coupled to two fermionic baths at T = 0 . The different components of the total heat production Q tot diss ( t ) = P tot diss ( t ) as a function of time. The energy of the level evolves in time with V A C = 0.7 and ω = 0.3 . Parameters: μ = 0.2 , δ μ = 0 , ε 0 = 0 , and hybridization widths Γ L = Γ R = 0.5 . Energies are expressed in units of Γ = Γ L + Γ R . The upper panel shows that the heat flux at the reservoirs may instantaneously attain negative values. The dissipative heat flux at the driven dot Q S diss ( t ) is highly fluctuating but the sum of the two contributions satisfies Q tot diss ( t ) 0 , consistent with the second law.
Figure 4. Nonadiabatic regime of a single driven quantum level coupled to two fermionic baths at T = 0 . The different components of the total heat production Q tot diss ( t ) = P tot diss ( t ) as a function of time. The energy of the level evolves in time with V A C = 0.7 and ω = 0.3 . Parameters: μ = 0.2 , δ μ = 0 , ε 0 = 0 , and hybridization widths Γ L = Γ R = 0.5 . Energies are expressed in units of Γ = Γ L + Γ R . The upper panel shows that the heat flux at the reservoirs may instantaneously attain negative values. The dissipative heat flux at the driven dot Q S diss ( t ) is highly fluctuating but the sum of the two contributions satisfies Q tot diss ( t ) 0 , consistent with the second law.
Entropy 18 00419 g004

Share and Cite

MDPI and ACS Style

Ludovico, M.F.; Arrachea, L.; Moskalets, M.; Sánchez, D. Periodic Energy Transport and Entropy Production in Quantum Electronics. Entropy 2016, 18, 419. https://doi.org/10.3390/e18110419

AMA Style

Ludovico MF, Arrachea L, Moskalets M, Sánchez D. Periodic Energy Transport and Entropy Production in Quantum Electronics. Entropy. 2016; 18(11):419. https://doi.org/10.3390/e18110419

Chicago/Turabian Style

Ludovico, María Florencia, Liliana Arrachea, Michael Moskalets, and David Sánchez. 2016. "Periodic Energy Transport and Entropy Production in Quantum Electronics" Entropy 18, no. 11: 419. https://doi.org/10.3390/e18110419

APA Style

Ludovico, M. F., Arrachea, L., Moskalets, M., & Sánchez, D. (2016). Periodic Energy Transport and Entropy Production in Quantum Electronics. Entropy, 18(11), 419. https://doi.org/10.3390/e18110419

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop