Work, Heat and Internal Energy in Open Quantum Systems: A Comparison of Four Approaches from the Autonomous System Framework

We compare definitions of the internal energy of an open quantum system and strategies to split the internal energy into work and heat contributions as given by four different approaches from autonomous system framework. Our discussion focuses on methods that allow for arbitrary environments (not just heat baths) and driving by a quantum mechanical system. As a simple application we consider an atom as the system of interest and an oscillator field mode as the environment. Three different types of coupling are analyzed. We discuss ambiguities in the definitions and highlight differences that appear if one aims at constructing environments that act as pure heat or work reservoirs. Further, we identify different sources of work (e.g. coherence, correlations, or frequency offset), depending on the underlying framework. Finally, we give arguments to favour the approach based on minimal dissipation.


I. INTRODUCTION
To lay a quantum mechanical foundation of thermodynamics a definition of the terms "internal energy", "work" and "heat" is a prerequisite.The first law of thermodynamics states that the change in internal energy ∆U is composed of work W and heat Q.In thermodynamics, work is defined as the change of internal energy in the absence of entropic changes, while heat is associated to a change in entropy.At first glance it seems straightforward to extend this concept to the quantum regime.However, the devil turns out to be in the quantum mechanical details, and over the years a range of disagreeing definitions have been proposed [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15].Ref. [16] provides a recent overview of the differences and relations between single-shot, external, dynamical map, autonomous and semi-classical approaches.
Broadly speaking, the approaches differ by whether they consider the external driving to be exerted by a classical system or by a quantum mechanical one, whether they include the measurement process explicitly in the description, whether the system is already coupled to an infinite bath or a part of the system is integrated out to form an environment (which does not necessarily have to be a bath in the thermodynamic sense), and whether individual systems or statistical ensembles are considered.We cannot give an exhaustive overview of all approaches here and refer the reader to reviews in Refs.[16][17][18].
Our work focuses on methods in the autonomous framework, in which both, the environment and the coupling between the system and the environment have a rather general form, i.e. we do not restrict the discussion to thermodynamic heat baths, weak coupling or specific initial environmental states.(By the term "thermodynamic heat bath" we understand continuum systems with infinitely many degrees of freedom in a thermal state.)Further, we consider the driving to be exerted by a part of the quantum mechanical system, not by an external classical system.In this regard the methods we discuss here are more general than other approaches.On the other hand, we do not discuss the measurement process, albeit it obviously being necessary for a complete theory.The reason we ignore the measurement process in this article is that we would like to point out some differences between several approaches which have been proposed in the literature.These difference will in general not be removed by a measurement, thus for arguments presented here, the measurement process does not need to be treated.
For a system with a Hamiltonian H characterized by the density matrix ρ the expectation value of the total energy changes according to dU = dTr(Hρ) = Tr(ρdH) + Tr(Hdρ) := dW + dQ, (1) where the last two terms are conventionally identified as dW and dQ [19].In this way work and heat are treated as average quantities and not as distributions.For a closed system, which cannot exchange energy (or anything else) with the surroundings, energy changes can only be due to a time-dependence of the Hamiltonian.Hence the change of the internal energy is equal to work.Usually external driving is imposed by a classical, macroscopic mechanism [8,9] (e.g. a varying magnetic field).
In contrast, the autonomous framework considers one quantum system as the control or environment for another one.Under certain conditions the environment can be seen either as a work or a heat reservoir, but it could also be a hybrid source.Both components can be treated equivalently such that the assignment of primary system and environment can be inverted.Tracing out the degrees of freedom of the environment leads to the reduced dynamics of the system of interest.The coupling of system and environment causes energy transfer and correlations.Now, the problem arises how to partition the total energy into the internal energy of the arXiv:2308.08215v1[quant-ph] 16 Aug 2023 primary system and of the generic environment, how to treat energy stored in the correlations, and how to label energy changes as work and heat fluxes.Roughly speaking, inclusive and exclusive perspectives can be taken in this context [20] as well as perspectives in between the two, i.e. perspectives in which a part of the interaction Hamiltonian is attributed to the energy operator of the system and in which the environment can contribute to heat and work exchange.
We will analyse and contrast four approaches: the attempt based on a local effective measurement basis (LEMBAS) [4], the non-local interaction method [1], the entropy-based formulation derived from spectral decomposition [2,3] and the minimal dissipation method [21].The first approach defines an effective system Hamiltonian which is compatible with the initial Hamiltonian and relates heat to entropy changes [4].The second approach aims at a separation of the Hamiltonian, in which the interaction energy is not accessible from the reduced states [1].The third relies on the spectral decomposition of the system state to relate changes of those eigenvalues which alter the local entropy to heat.And the last approach is based on the minimization of the dissipative part of the equation of motion [21].
To compare these methods we consider the spinoscillator model with different types of interaction as an example.

II. METHODS TO DEFINE WORK AND HEAT
A closed quantum system, which is isolated in the thermodynamic sense, exchanges neither matter nor energy with the environment.If we describe its state by a density matrix ρ, and if the time evolution is generated by a Hamiltonian H, it is natural to identify U = Tr(Hρ) as the internal energy.The internal energy can only vary in time if the Hamiltonian is time-dependent which is usually identified as the work rate [7].(Throughout this article we use ℏ = 1.)For a genuine quantum mechanical description and in order to consider the most general processes we work in the quantum autonomous framework [16].We start with a composite system that includes the source of driving, such that the total system dynamics are unitary and generated by the static Hamiltonian H.The total state is given by a density operator ρ on the Hilbert space H and initially prepared in a product state ρ(0) = ρ S (0)⊗ρ E (0), where the index S denotes the system of interest (or primary system) and E an arbitrary environment.In line with this, the Hamiltonian is composed of the bare Hamiltonians H S and H E and an interaction term H SE .The reduced system state is obtained by tracing over the degrees of freedom of the environment, i.e. ρ S = Tr E (ρ).This setup defines a completely-positive trace preserving dynamical map Φ t which maps the initial system state ρ S (0) to with unitary U(t) = e −iHt [22].From hereon, unless indicated otherwise, we drop the dependence of the density matrices on time t to keep the notation compact.
The equation of motion for the system can now be written as where we have introduced The most common approach identifies U S = Tr S (H S ρ S ) as the internal energy of the system.Work W S done on S is then related to changes in energy due to time-dependence of H S and heat Q S is related to the variation of the reduced state [19].Analogously U E = Tr E (H E ρ E ) can be defined for the environment, but in general one cannot retrieve the internal energy of the composite system by adding U E and U S with this choice, because the energy contribution Tr(H SE ρ) is not taken into account.
Alternatively, a completely inclusive view can be taken by assigning the entire interaction energy to the system of interest and defining U S = Tr((H S + H SE ) ρ) [20].This can be inferred if the environment is a heat bath, i.e. it just exchanges heat [7,23].Variations of the internal energy in E are then identified as heat transferred to or from the system Q = − UE = −Tr E (H E ρE ).
With the notion that work contributions only arise due to an explicit time dependence of the Hamiltonian (and H E being time-independent) the whole work done by or performed on the composite system is attributed to S. In the limit of weak coupling Tr((H S +H E ) ρ) ≫ Tr(H SE ρ) we obtain the conventional definition given above in Eq. (1).
The opposite perspective is taken if the system is coupled to a secondary system which can be labeled as work depository, battery, weight or control system.Work is identified as the negative average energy change of this environment [16,24]: Generally, these approaches impose particular requirements on the environment or on the coupling.An ideal heat source should remain in a thermal state at all times, while a perfect work reservoir does not build up correlations with the system such that the reduced dynamics of S are unitary.
The approaches we discuss in the following differ from this idealized view in that they allow the notion of an environment as a hybrid source of work and heat.
Further, the relations above do not take into account that the partition of the total Hamiltonian into system and environment is not unique.The schemes we discuss in the following give distinct motivations for a partition which is applicable to arbitrary setups.We label the effective, rearranged Hamiltonians with superscripts A, B, C and D according to the respective subheadings.

A. Local Effective Measurement Basis
We start with the local effective measurement basis (LEMBAS) approach introduced in Ref. [4].The term "local" indicates that we are only interested in the properties of one component (here either the system S or the environment E) of the bipartite system.The choice of the measurement basis depends on the context, thus work and heat will be basis dependent.Since in this article we are interested in defining energies in S, the eigenbasis of the bare Hamiltonian H S is a natural choice.
The chosen basis is used as reference for the evaluation of work or heat [4].The effect of the actual measurement on the dynamics is not considered in this scheme.The effective operator is where we have split the correction Eq. ( 6) into , and H ′ S,a is the part that commutes with H S , while H ′ S,b is the part that does not commute.In other words, H ′ S,a is the part of H ′ S which is diagonal in the energybasis of H S .The internal energy of the system is determined by U S = Tr(H A S ρ S ).Heat flux is then identified as the contribution to the change in internal energy that has an effect on the local von Neumann entropy of the system state S S = −Tr S (ρ S ln ρ S ), while work flux is the contribution that does not affect S S .

QS (t) = −iTr([H
Note, that the work flux is not simply defined as Tr S ( ḢA S (t)ρ S ), but related to unitary dynamics.Further, there can be a change in entropy without a corresponding heat flux [5].
The evaluation of the heat flux still requires knowledge of the entire system state.It vanishes if the dynamics does not induce bipartite correlations, i.e. χ SE (t) = 0, or H A S commutes with H SE [5].
Since commutation relations clearly govern the energy fluxes it can be useful to break the interaction into H SE = k S k ⊗ E k where S k and E k are operators that act on the respective subsystem, such that ) If both, the system and the environment operators commute with the free Hamiltonians neither heat nor work flows.
In the opposite case ([S k , H S ] ̸ = 0, [E k , H E ] ̸ = 0) the subsystems exchange work and heat with each other.Then we have The fluxes are, in general, unequal [25].
It is apparent that the initial state of the composite system determines if there are work fluxes.At least one of the subsystems has to feature non-diagonal elements with respect to the energy eigenbasis.

B. Non-local Interaction
To determine how two systems exchange work and heat with each other, for a start the partition of the interaction energy could be considered.In general, correlations χ SE are built up between the two components and thus part of the total energy will neither be accessible through the system nor the environment state alone.In Ref. [1] the internal energy of each subsystem is identified as the part which is accessible by arbitrary local measurements.
Instead of looking for a compatible correction to the system Hamiltonian in this approach the idea is to rearrange the Hamiltonian into local parts that are associated with system and environment respectively and a third component that is only available through the composite system state.This requires Tr(H B SE (ρ S ⊗ I E )) = Tr(H B SE (I S ⊗ρ E )) = 0 for the effective interaction Hamiltonian.Defined like this, the interaction energy is neither accessible from the system nor the environment by itself.It is given by and accordingly where α S and α E are free parameters that add up to one.In this framework the internal energy of the entire composite system is made up of the contributions from the subsystems and, additionally, the "binding energy" U χ = Tr(H B SE (t) χ SE ).With these choices the work flux can be expressed as The last equality is obtained by treating the environment analogously to the system.Thus the work fluxes always add up to zero.They clearly depend on the choice of the free parameter and are thus not invariant under gauge transformations [13].The fact that Hamiltonians H and H + g(t)I generate the same dynamics but lead to different conclusions for work was pointed out in Ref. [26] and has caused some debate about the physical significance of the non-local approach.Different setups can imply a specific gauge [1,13], but a clear prescription for generic situations is not available.One possibility would be to remove contributions proportional to the identity which clearly do not influence the dynamics by the effective Hamiltonian.In the non-local approach this would amount to choosing α S = 0.But this, implicitly, identifies α E = 1, such that Tr(H SE (ρ S ⊗ ρ E ))I is the last term in H B S , Eq. ( 15).We see this definition, which combines properties of system and environment, as a drawback of the non-local scheme.
In contrast, the change in heat does not depend on the free parameters and it has the same structure as Eq. ( 9): but with a different effective Hamiltonian.The heat fluxes add up to the change in binding energy ( QS + QE = − Uχ ).
As in the LEMBAS approach, correlations are necessary but not sufficient for a finite heat flux.Or conversely, the environment can also be considered a work source if the commutator in Eq. ( 17) vanishes.

C. Decomposition
Recently, a purely entropic motivation for a split of internal energy changes into heat (related to changes in local von Neumann entropy S S ) and work (not related to changes in S S ) was proposed in Ref. [2] and Ref. [3].The authors base their approach on the spectral decomposition of the system state and the observation that the entropy varies if the eigenvalues r k change: ṠS (t) = − k ṙk (t) ln r k (t), but remains constant if just the eigenstates and not the eigenvalues are modified by the dynamics.Thus, instead of identifying heat flux as Tr S (H S ρS ) only the change in eigenvalues is taken into account: The remainder due to the change of eigenstates is attributed to work, thus adding a term to the conventional definition Eq. ( 2), i.e.

ẆS (t) = Tr
The second term is labeled "environment-induced dissipative work" caused by counter-diabatic dynamics [2].Internal energy is conventionally identified with U S = Tr S (H S ρ S ).
These assignments have the advantage that they solely rely on knowledge of the system state and the bare Hamiltonian.However, it seems impractical to decompose the density matrix at each point in time, especially for systems beyond two-level-systems.
The decomposition can be exploited to define a trajectory-based Lindblad-like master equation ρS = where and d is the system's dimension [2].With these identifications the heat flux, Eq. ( 19), can also be expressed as Clearly, the effective Hamiltonian K S obtained in this way, the rates and the Lindblad operators depend on the initial state of S. We do not label the effective Hamiltonian H C S , because it is not associated with the internal energy in this approach.However, it would seem natural to use K S as the operator for an inclusive identification of internal energy.Especially, since H S can be time-dependent in the original works [2,3] and thus, it would be consequent to include the time-dependence, induced by the embedding in an environment, as well.That would be in line with other approaches to define work, heat and internal energy from a master equation perspective [7,21], as discussed in the next section II D.

D. Minimal dissipation
In Ref. [21] an approach based on the master equation description for the dynamics of an open system is introduced.With the time-convolutionless projection operator technique it is possible to construct an exact, time local master equation for arbitrary coupling strengths and environmental temperatures [27,28].
The propagation of an initial system state to some later time t ≥ 0 can be expressed via the dynamical map For varying t we have a family of completely positive and trace-preserving maps [22].The generator of the dynamics is now ρS (t) = L t [ρ S ] = Φt Φ −1 t [ρ S (t)] if the time-dependence of the map is sufficiently smooth and the inverse exists, which is typically the case [29].From the requirement of Hermiticity and trace preservation follows that the generator can be rewritten as ] and a dissipator with generalized Lindblad structure with time dependent rate functions {γ k (t)} and operators L k (t).This decomposition is not unique, i.e. for each time t the transformation to Lindblad operators and Hamiltonian with arbitrary scalar functions {α k (t)} and β(t) leaves the generator invariant [21].However, an unambiguous choice can be made for D t if we demand a minimal dissipator which corresponds to traceless Lindblad operators [30].
Once the minimal dissipator is determined the corresponding K S defines H D S .Note, that the minimization does not fix β(t) but, because it is irrelevant for the dynamics, it can be set to zero.
The change of internal energy, work and heat is now given by δQ The advantage of this identification is that it just depends on the system state ρ S and could, in general, be applied without knowledge about the environment.In addition, it is based on an extremum principle.However, the useful correspondence between minimal dissipator and traceless Lindblad operators is only established for finite Hilbert spaces.This limitation already causes problems if a single field mode is considered as the system.Even though this can be circumvented by appropriate truncation of the Hilbert space if the dynamics are restricted to a finite subspace, it remains a minor short-coming of an otherwise very convincing approach.

III. APPLICATION TO EXAMPLES
To compare these different approaches we study the implications for a simple model and three types of interaction.A two-level system (TLS) or atom is considered as the system of interest and a mode of the radiation field as its environment.This is a standard example in quantum optics and it has been considered in various variations in quantum thermodynamics.It is also frequently encountered as an analogue to a steam engine with a working medium (gas, TLS) coupled to a work depository (piston, field mode) [16,31].One advantage of this system, and probably the reason for its frequent use as a model, is that it can be treated analytically.Further, semi-classical analogues are available, that can be exploited for a comparison with semi-classical results The free Hamiltonian is H 0 = H S + H E with H S = ω S σ † and H E = σ + ω E a † a where ω S and ω E denote the frequencies and σ and a the lowering operators of atom and field, respectively.
Atom and field are initially uncorrelated, i.e. they are in a product state ρ(0 and the field is in a thermal state, i.e. ρ E (0) = (1 − e −βω E )e −βH E where β = 1/k B T .
For the minimal dissipation approach we have to determine the master equation for the reduced density matrix and identify the effective Hamiltonian and traceless Lindblad operators with corresponding rates.The procedure for this setup is detailed in Appendix A.

A. Jaynes-Cummings
At first we consider the interaction Hamiltonian given by the Jaynes-Cummings model: with interaction parameter g.H SE does not commute with the free Hamiltonians of atom and field: We obtain the correction Hamiltonian S,b which does not commute with H S (see Eqs. ( 6) and ( 8)).Consequently, the effective atomic Hamiltonian equals the bare Hamiltonian in A: H A S = H S .In this way U S corresponds to the conventional identification, Eq. ( 1).The work and heat fluxes are given by Eq. ( 11) and Eq. ( 12).Only on resonance we obtain zero net work flux ẆS + ẆE = 0 and heat flux QS + QE = 0, because then [H SE , H S + H E ] = 0.
In B the correction H ′ S contributes to the effective local Hamiltonian, provided correlations are built up.In the explored parameter range this results only in a small deviation for the integrated fluxes compared to LEMBAS (compare the lines labelled A and B in Fig. 1, where α S = 0).For large interaction strengths g and considerable coherence in both, atom and field state, the difference between the quantities increases, because it arises due to H SE and the development of correlations.
The identification of work and heat in the decomposition scheme (lines labeled C in Fig. 1) is qualitatively similar to the previous two approaches.The change of the average internal energy is identical to the LEMBAS picture for this non-commuting interaction.However, the work estimate is in general lower, and the heat estimate higher.By design the decomposition of the state yields a heat flux that is closely related to the local change of entropy (compare the upper panel of Fig. 4 with δQ S in 1: the extrema coincide).
All three approaches agree that there is no work done by or on the system if the two-level system (and the field) is initially in a diagonal state with respect to the energybasis, i.e. p eg = 0.
The average energy fluxes determined with the minimal dissipation framework are rather different (see lines labeled D in Fig. 1).Especially, the behaviour of internal energy and work can be more "erratic", i.e. in comparison to the methods A, B and C we observe peaked features for certain setup parameters.As we approach resonance these features become more prominent.The peaks are not related to the change in entropy or the build-up of mutual information (see Fig. 4) but closely linked to the dynamical map that generates the dynamics.The peaks appear if a singularity is approached, i.e. a point in time where Lρ = ρ has no unique solution or equivalently, the map Φ t is not invertible (see the Appendix A for details).Whenever the determinant of F t (the matrix representation of Φ t ) approaches 0 a peak forms, which does not occur in the other frameworks.A possible interpretation is that a high energetic cost has to be paid to prepare states close to a crossing, where the information about the initial state would be lost.Or conversely a large amount of energy could be extracted by a transformation of the form U ⊗ I on system and environment [32].
Contrary to methods A, B and C, in the minimal dissipation approach a finite energy change associated with work is only obtained if there is a frequency offset of the bare Hamiltonians, i.e. ω E − ω S ̸ = 0.This bias leads to a shift of the populations in the field mode away from a thermal state but remaining diagonal in the eigenenergy basis.Thus one could associate local temperatures to each energy level.The bias offers a resource of work which can be consumed (provided) by S and increase (decreases) the energy stored in the correlations [33].Initial coherence with respect to the energy eigenbasis of the TLS or the field is not necessary, since the timedependence of the effective Hamiltonian is determined by the initial state of the environment and the interaction Hamiltonian (see Appendix A).But coherences enhance the exchange of work as well and can be identified as a resource.
The heat defined by minimal dissipation is always smaller than in the other approaches and thus serves as a lower bound.In the evolution without coherences β∆S S and δQ S share the same extrema and on resonance they coincide.Comparison of the integrated energy fluxes for Jaynes-Cummings coupling defined by the methods presented in the previous section II: A LEMBAS, B non-local (where αS = 0 is used), C decomposition, D minimal dissipation for ωE = 0.9ωS, kBT = 1.0ωS, g = 0.1ωS, pe = 0.25, peg = 0.1i.

B. Displaced spin oscillator
Next, we consider another kind of coupling: which is known as the displaced spin oscillator model.It has already been regarded as an analogue to a steam engine where the oscillator represents the piston that compresses and expands the spin as the "working gas" [31].
The displacement a + a † of the oscillator modulates the splitting of the spin's energy levels [34].During the interaction the occupations remain constant and mutual information is build up and destroyed again at the field frequency ω E even without initial coherence in ρ S and on resonance.Likewise, the local entropies vary as depicted in Fig. 4. For p eg = 0 the local von Neumann entropy of the system remains constant ∆S S = 0 for p eg = 0 since the .This interaction is partially commuting since [σ z , H S ] = 0 and [a + a † , H E ] ̸ = 0 [5].The correction H ′ S for the atomic Hamiltonian is time-dependent and compatible with the bare system Hamiltonian H S such that with A Thus there is no heat flux, because the effective Hamiltonian commutes with the interaction even if correlations are built up between atom and field (see Eq. ( 9)).Still, the local entropy S S oscillates (see Fig. 4).So the entire entropy change is due to entropy production Note, that the entropy production rate ΣS is temporarily negative.This violates the second law of thermodynamics, if it is understood as the positivity of the entropy production rate at all times [7].The violation can be related to non-Markovianity and memory effects in the dynamics [21].
This is also true for B, where H B S only differs by a scalar (if α S = 0, they coincide).Now, coherence in the initial state is not a prerequisite to obtain a timedependent correction to H S and in fact it has no influence on it.Because of the time-dependence both approaches define a work flux and varying internal energy.We observe that ω E determines the oscillation frequency of ∆U S , δW S and the amplitude is given by −g 2 ω E and (α S − 1)g 2 ω E respectively (see Fig. 2).
If the instantaneous basis of the system state is employed to infer the thermodynamic energies in C, the results suggest different interpretations.As internal energy is the expectation value of the bare atomic Hamiltonian H S and the occupations remain constant during the interaction, we have ∆U S = 0.Although there is no net energy change, there is still an oscillatory exchange of work and heat QS = − ẆS .
For D the internal energy does not change and there is no work flux because we do not obtain a time-dependent (non-trivial) correction to H S (see Appendix A).So in this formulation, it is not sufficient to have coherences or a frequency-offset for a work-like variation of U S .
Although the dynamics is described with a dissipator with non-vanishing rate, it does not contribute to heat since the relevant Lindblad operator σ z commutes with H S such that: So the change of local system entropy S S is clearly not related to energetic changes and solely due to entropy production.

C. Dispersive spin oscillator
For the sake of completeness, we also consider the dispersive spin-oscillator, an interaction that is fully commuting with the free Hamiltonian: This coupling leaves the occupations of the atomic levels unchanged and only influences the S S but not S E .
Since σ z and a † a commute with H S and H E the corrections H ′ S (t) = gσ z ⟨a † a⟩, H ′ E = g⟨σ z ⟩a † a commute with the free Hamiltonian and this results in vanishing energy fluxes in the LEMBAS setting A.
By construction, the non-local approach B does not show any energy fluxes, either.For the heat flux this is obvious because it also depends on the commutator [H B S , H SE ] which vanishes, and because the additional scalar in Eq. ( 14) has no effect in the commutator.Correlations between the subsystems are necessary but not sufficient to yield a heat flux.The work flux is also zero, because the dynamics only varies the off-diagonal elements in the reduced system states.
The analysis of the instantaneous eigenbasis C leads to a different result: heat and work are again oscillating and they add up to zero.I.e. each influx of heat is compensated by an outflux of work and vice-versa.This interpretation suggests, that this setup is not suitable to model a heat or work reservoir.
In contrast, the minimal dissipation approach yields an effective Hamiltonian H D S (t) that varies in time and leads to a periodic variation of internal energy and work.As in the displaced interaction we also obtain a dissipator with a single time-dependent rate and operator σ z (see Eq. ( 36) and Appendix A).Thus there is no heat transfer due to these incoherent dynamics.Since the entropy is still changing, we can infer that the entropy production rate is directly related to this part of the dynamics.

IV. CONCLUSION
We have compared four different approaches to identify work, heat and internal energy in bipartite quantum systems.All four can be applied to a system coupled to an arbitrary environment, which does not need to be at equilibrium with a coupling that does not need to be weak.All methods provide a motivation to find a partition of internal energy into a component that locally has a work effect and one that locally has a heat effect.Further, they agree that heat is associated with a change in von Neumann entropy, but still reach different partitions.In infinitely large systems, at weak coupling and close to thermodynamic equilibrium, the von Neumann entropy of the subsystem coincides with the thermodynamic entropy [34].Accordingly all identifications of heat (and work) should transition to their thermodynamic counterparts in this limit.
In approaches A (LEMBAS), B (non-local) and D (minimal dissipation) work is related to energy exchanges due to changing parameters of the system Hamiltonian.Nonetheless the effective, time-dependent Hamiltonian that arises due to the interaction is defined differently.The LEMBAS framework (A) gives an operational view on work, by presuming that an observer knows the bare system Hamiltonian and will perform measurements with respect to its eigenbasis.So only the internal energy and work accessible in this way is identified.We note that B gives an ambiguous definition, which is to some extent problematic.The local effective Hamiltonians of S and E are inferred in tandem such that at least one of them includes a term proportional to (I) and varying in time, an arguably unphysical situation.Criticism might also be brought forward against C (the decomposition scheme) since the treatment of the time-dependence in the coherent part of the dynamics is inconsistent.
Method D offers a unique allocation [35] in terms of system variables and it is widely applicable.Apart from exact time-convolutionless master equations the approach can be used for approximate time-local master equations as well.The drawback is the limitation of the correspondence between minimal dissipator and traceless Lindblad operators to finite Hilbert spaces [36].
If we consider the ability of the methods to identify when an environment acts as a classical driver, A is especially transparent.The variation of the internal energy is solely due to work if the interaction Hamiltonian commutes with the effective system Hamiltonian or if the composite system remains in a factorized state.We observe this for the partially and fully commuting interaction discussed in the spin-oscillator model (see Sec. III B and Sec.III C).But for the latter we also have ∆U S = 0.In D the heat contribution vanishes if the rates in the dissipator vanish or the Lindblad operators L K commute with the effective system Hamiltonian.This is equally found for partially and fully commuting couplings, but now ∆U S ̸ = 0 for the dispersive interaction and it vanishes for the other.
Quite differently, the diagonalization scheme C does not render any of the interaction types suitable for just driving.To obtain a qualitative comparison of how well an environment performs, a measure for the quality of a work source as suggested in Ref. [34] could be used.
Further, we conclude that the different schemes identify different sources.In A and B coherences are essential, while D also ascertains that the off-set in eigenenergies of system and environment can yield an energy change in the form of work.
In addition, the different notions of heat also lead to different interpretations of entropy production if we consider Σ S = ∆S S − βδQ S .Thus, statements about the second law of thermodynamics depend strongly on the underlying definition of thermodynamic quantities.A detailed discussion would require a particular definition of effective system temperatures applicable in nonequilibrium situations but this is beyond the scope of this investigation.Further investigation is required to determine which of the formulations is suited best for a discussion of accessibility in measurements and which part of the allocated work can actually be used.

FIG. 2 .
FIG.2.Integrated energy fluxes in the displaced spin oscillator model determined for the same parameter set as in Figure1..

FIG. 3 .
FIG. 3. Integrated energy fluxes obtained in the dispersive spin oscillator model.The same parameter set and labels as in Figure 1 are used.