Manybody effects on the thermodynamics of closed quantum systems
Abstract
Thermodynamics of quantum systems outofequilibrium is very important for the progress of quantum technologies, however, the effects of many body interactions and their interplay with temperature, different drives and dynamical regimes is still largely unknown. Here we present a systematic study of these interplays: we consider a variety of interaction (from noninteracting to strongly correlated) and dynamical (from sudden quench to quasiadiabatic) regimes, and draw some general conclusions in relation to work extraction and entropy production. As treatment of manybody interacting systems is highly challenging, we introduce a simple approximation which includes, for the average quantum work, manybody interactions only via the initial state, while the dynamics is fully noninteracting. We demonstrate that this simple approximation is surprisingly good for estimating both the average quantum work and the related entropy variation, even when manybody correlations are significant.
pacs:
I Introduction
Progress on applications of quantum technologies is linked to acquiring deeper understanding of the outofequilibrium thermodynamic properties of small quantum systems. Small quantum systems operating at finite temperature constitute the hardware for most of these technologies, so that outofequilibrium quantum thermodynamics guides, and ultimately may limit, these technologies Goold et al. [2016], Vinjanampathy and Anders [2015], Millen and Xuereb [2016], Parrondo et al. [2015], LuizzoScorpo et al. [2016]. Quantum thermodynamics has then become a rapidly expanding field, also supported by the advances in experimental techniques. Nonetheless, a remaining challenge is properly accounting for the effects of manybody interactions on quantum thermodynamics properties. Manybody interactions in quantum systems give rise to complex phenomena, such as collective behaviours and quantum phase transitions. From a practical perspective, these effects are challenging to calculate and often require approximations. In this respect, there have been recent works studying outofequilibrium thermodynamics of manybody systems such as quantum harmonic oscillator chains and spin chains Silva [2008], Dorner et al. [2012], Joshi and Campisi [2013], Mascarenhas [2014], Sindona [2014], Fusco et al. [2014], Zhong and Tong [2015], Eisert et al. [2015], Bayat et al. [2016], SolanoCarrillo and Millis [2016], and a proposal for a densityfunctionaltheorybased set of approximations which is in principle applicable to systems of high complexity Herrera et al. [2017, 2018].
Here we present a systematic study of the outof equilibrium thermodynamics of manybody quantum systems subject to a set of qualitatively different driving potentials. For each type of driving potential, we consider dynamical regimes from sudden quench, to finite times, to quasiadiabatic; for each dynamical regime we consider different interaction strengths, from noninteracting to strongly correlated systems. For each driving potential, dynamical regime, and interaction strength, we consider different temperatures: low, intermediate and high temperature. For all cases considered we calculate and discuss the average quantum work extracted and entropy produced in the dynamical trajectory. Our systematic study allows to uncover some important dependencies of work and entropy on the systems’ correlation and dynamical regimes, which cut across the different applied drives.
Afterwards we consider two quite drastic approximations and compare their estimates with the exact results. The first is the completely noninteracting approximation, where manybody interactions are set to zero in all phases of the thermodynamic processes considered. The second approximation assumes knowledge of the initial interacting manybody state, but completely neglects interactions afterwards, during the driven dynamics. Our results show that including interactions just within the initial state provides surprisingly good accuracy. We provide an analytical analysis that explains this accuracy in the sudden quench and adiabatic regimes for a general Hamiltonian and driving potential.
Ii Theory
ii.1 Hubbard model
The onedimensional Hubbard model can depict systems from weakly to strongly correlated and model numerous phases of matter and related phase transitions, such as metallic, antiferromagnetic, Mottinsulator, superconductivity, and FFLO transition França and Capelle [2006], França et al. [2012], de Picoli et al. [2018], Essler et al. [2005]. It is being widely used to study many physical systems, from coupled quantum dots, to molecules, to chains of atoms Coe et al. [2010], Yang et al. [2011], Murmann et al. [2015], Coe et al. [2011], Brown et al. [2019], Nichols et al. [2019]. These are systems of importance, as hardware, to quantum technologies. For small chains, the Hubbard model is numerically exactly solvable, yet still displays nontrivial behaviours, including, for repulsive interactions, the precursor to the metalMott insulator phase transition. Hence it is often the system of choice for exploring approximations to interacting quantum systems Herrera et al. [2018, 2017], Carrascal et al. [2015], Fuks and Maitra [2014].
For a fermionic system of sites, the Hamiltonian of the Hubbard model can be written as ()
(1) 
where is the hopping parameter, () is the creation (annihilation) operator for a fermion with spin ( or ) in site , is the strength of the onsite Coulomb interaction, is the spin , site number operator, , and is the onsite potential.
ii.2 Average quantum work and entropy variation
The average quantum work, much like its classical counter part, is described as the usable energy in a quantum system Vinjanampathy and Anders [2015]. In a closed system at temperature , it can be calculated as Vinjanampathy and Anders [2015]
(2) 
with and the initial (final) system state and Hamiltonian, respectively.
For a given dynamic process, the variation in thermodynamic entropy is defined using the average work and the change in the free energy of the system Batalhão et al. [2015], Goold et al. [2016],
(3) 
where , and the free energy variation is
(4) 
with the partition function at the beginning (end) of the dynamics, . This thermodynamic entropy can be considered a measure of the degree of irreversibility of the system dynamics: in fact it captures an uncompensated heat which would need to be dispersed to the environment for the system to return to thermodynamic equilibrium at the end of the driven process Herrera et al. [2017], Batalhão et al. [2015].
Iii Systems, driving potentials, temperature range and dynamical regimes
We shall calculate the work extracted and entropy produced for different systems, temperatures, driving potentials, and dynamical regimes.
We will consider Hubbard chains of 2, 4, and 6 sites, at halffilling and under open boundary conditions, and explore low (), medium (), and high () temperatures.
For each system size and temperature, we will explore regimes from noninteracting () all the way to strongly correlated (), and dynamics from sudden quench (, the overall driving time) all the way to quasiadiabatic ().
For each parameter combination, we will consider three types of driving potentials ^{1}^{1}1We note that for the 2 site chain the only relevant dynamics becomes the ‘Applied Electric Field’ one., where each potential has a linear time dependency via , with and the timeindependent coefficients for site at time and respectively. With this choice, the character of the dynamics will depend on , while the final Hamiltonian will be independent of it. These driving potentials are:

“Comb”: for each site , at and at , where and .

“Middle Island (MI)”: the inhomogeneity is driven only for the middle two sites of the chain; for where goes from 1 to , and is the chain length. For the middle two sites, , and .

“Applied Electric Field (AEF)”: this potential mimics the application of a potential difference between the extremes of the chain. The sites form a linear slope from to and are described using where , and with .
The and form of the driving potentials for a sixsite Hubbard chain are illustrated in figure 1.
Iv Approximations
iv.1 Average quantum work
The type of approximations we consider for the average quantum work are of the form:
(5) 
Here (initial system) refers to the approximation used to derive the system state at , , and evo is the approximation used for the evolution operator where is the timeordered operator. The final state is then . We note that . In the approximation where and evo are the same, only one acronym shall be written. As (5) indicates, the Hamiltonians and explicitly entering (2) and the evolution Hamiltonian are to be taken in the same approximation: if this is not the case, we found that the mismatch in eigenstates leads to spurious oscillations in the work production (not shown).
We will consider two approximations, as described in Table 1. The first, , corresponds to a completely noninteracting system, the one obtained by setting in the Hubbard Hamiltonian. The second approximation, , uses the exact manybody initial state, but the approximated (noninteracting) Hamiltonian for the evolution of the system, according to the notation previously introduced.
Acronym  Approximation  Hamiltonian  Initial State 

NI  
exact + NI 
iv.2 Entropy variation
When approximating the entropy, we use , while the free energy is estimated in the same way as . In the NI approximation, this implies that the free energy is constant. For , the exact free energy is then used as this uses the same assumption made for calculating the initial thermal state, that the system Hamiltonian can be exactly – or very accurately – diagonalised.
V Work extraction
The system is considered at equilibrium at time , when the coupling with the thermal bath is switched off. Then the closed system is driven by a timedependent external potential from the initial Hamiltonian to the final Hamiltonian in a time , and the extracted work from this dynamics is calculated according to (2), with .
We stress that with each of the driven dynamics described in section III, the final Hamiltonian is the same for all ’s, so that the latter controls the rate of driving. Therefore, the larger is, the slower the system has evolved and hence more adiabatic the evolution.
For each of the manybody systems, their approximations, temperatures, and driving potentials described in section IV, we will consider the parameter space and . Due to the sheer number of results from all the combination of parameters, we will only explicitly show results for 6 site chains (all 6 site chains’ exact results, and part of them for the approximations), and comment on the rest.
v.1 Exact results
Figure 2 shows the exact average quantum work extracted from a 6 site chain driven via “applied electric field (AEF)” (right column), “comb” (middle column), and “middle island (MI)” (left column) potentials at temperatures of (first row), (second row), and (third row). Each panel shows a wide range of regimes: from noninteracting to strongly correlated systems as increases along the axis; and from sudden quench towards adiabaticity as increases along the axis. A lighter shade of colour corresponds to higher extracted work.
Figure 2 presents a variety of behaviours, with the extracted work varying over a wide range of values, and even from positive to negative. This confirms that the chosen dynamics are a good testbed for understanding work extraction in systems representable via Hubbard chains, and hence a good test bed for related approximations.
At all temperatures, the largest work can be extracted via the AEF dynamics, while work needs to be done on the system in order to perform the MI dynamics. For all driving potentials, increasing temperatures decreases both the range of extracted/provided work and the maximum extractable work (the minimum work to be performed on the system in the case of MI). The maximum applied potential difference at is comparable to the highest temperature, so, as the temperature rises, the systems become less sensitive to the applied field. For all temperatures and dynamics the highest work that can be extracted (the lowest work performed on the system in the case of MI) is reached for large ’s, as the system gets closer to adiabaticity: here the dynamical state better adjusts to the driving force and, compatible with temperature and manybody interactions, the energetics favour lowpotential chain sites.
In the case of AEF, at all temperatures, increasing hampers the transient current dynamics, so the highest work is achieved for zero to weak correlations. The maximum potential step between nearby sites is about so, as increases, the AEF dynamics tends to freeze and lesser and lesser work can be extracted from the system. At the highest temperature, the thermal energy is almost equal to the potential difference between the chain extremes at the end of the dynamics.
At low temperatures, ‘comb’ dynamics presents a work extraction pattern similar to AEF; however, as the temperature increases, maximum work extraction can be achieved for higher manybody interactions, and at high it is achieved only for relatively strong manybody interactions (). This can be understood by realising that in this case the thermal energy is twice the potential barrier between even and odd chain sites, and hence a certain degree of repulsion is necessary to depopulate the highenergy sites completely and maximise work extraction.
Extracted work under MI dynamics is negative, meaning that work must be performed on the systems to achieve the final states. Indeed in this case the drive raises the potential of the central sites and, under the dynamics, the (closed) system cannot decrease its overall energy. As even for the system is not fully adiabatic, a finite amount of repulsion helps to completely deplete the central island, even at low temperature, and more so as temperature increases. The maximum barrier between high and low potential sites is about half of the maximum thermal energy considered.
v.1.1 Adiabatic regime
From figure 2 we see that, as increases and we approach the adiabatic regime, the average work becomes strongly dependent on and weakly dependent on . This behaviour can be easily explained. For a closed system, can be also calculated from the work distribution function, and hence written as Vinjanampathy and Anders [2015]
(6) 
where () is the th (th) eigenvalue of the final (initial) Hamiltonian, and is the joint probability distribution of arriving to the eigenstate of , given the probability of initially being in the eigenstate of .
In this formalism, the effect of the dynamics driven by the external potential (and hence the dependence on ) is contained in the conditional probabilities . In the adiabatic regime , so that we obtain
(7)  
which indeed depends on the interaction coupling but not on . By looking at figure 2 we can then note that the value of at which the different systems enter a nearly adiabatic behaviour depends both on the driving potential and on the temperature, with the ‘comb’ driving potential getting closer to adiabaticity for smaller values of , and increasing temperature seemingly decreasing adiabaticity. The latter can be understood by realising that when temperature increases to the point of dominating by far interparticle interactions, the system behaves as noninteracting, and hence the dependence of any quantity over will be lost (see also discussion in section V.2.1).
v.2 Approximated results
We will firstly approximate the manybody system with a noninteracting system and use this to estimate quantum thermodynamic quantities. Afterwards we will extend this approach to include some memory of the electronelectron interaction through the system initial state.
v.2.1 Noninteracting approximation
For noninteracting (NI) systems, , the average work extracted has no dependence on . This is shown in the upper panels of figure 3, where we plot for ‘comb’ dynamics, with (axis) and (axis), and low (left), intermediate (middle), and high (right) temperatures. The range of variation in decreases with temperature (see colour bar above each panel), as the maximum driving potential becomes comparable to the thermal energy.
The lower panels of figure 3 show the corresponding relative error of the NI approximation, where the darker the blue, the more accurate the approximation is in that regime. This approximation accurately captures the work extraction only in the parameters regions where interactions are weak compared to the other energy scales. These regions include higher values of as the temperature increases, see related discussion in section V.1. However it is worth noting that there is very little work extracted at high temperatures, where the thermal energy is comparable to the driving potential which is then less effective.
As the lowest values of the average extracted work are driven by the freezing of the system dynamics due to strong manybody repulsion (dark areas in figure 2), the NI approximation strongly overestimates the minimum work that can be extracted by a system. To see this compare the work range indicated over each panel for the mid column of figure 2 to the corresponding upper panels of figure 3. At the other end, the value of the maximum average work extracted is quantitatively well captured by the NI approximation, but the corresponding parameter regions are qualitatively wrong, compare the shape of the lightshade areas of the midcolumn panels of figure 2 to the corresponding areas of the upper panels of figure 3.
We find a similar accuracy pattern with increasing temperature for 2 ^{2}^{2}2See Herrera et al. [2017] and Herrera et al. [2018] for 2site chain examples, and especially Herrera et al. [2017] for 2site noninteracting results and 4 site chains.
The ‘comb’ driving potential corresponds to an accuracy of the results consistently inbetween those of the AEF and MI potentials.
For the MI driving potential, the NI approximation works better than for the ‘comb’ potential at all temperatures. As discussed in section V.1, for these dynamics manybody interactions become comparatively dominant at higher values of . As a consequence, for 6 site chains, the MI driving potential results in 10% accuracy (or better) for at low temperatures, for at intermediate temperatures, and for all regimes at high temperatures: here thermal energy dominates so that the average extractedwork range is very narrow and so weakly sensitive to parameter changes.
For the AEF driving potential, the results in the NI approximation are in general worse than with the ‘comb’ driving potential, and comparatively worsen as the temperature increases. At difference with the MI and comb driving potentials, the maximum AEF potential difference between nearby sites becomes at most of , so that even a Coulomb repulsion of just will remain relevant. The NI approximation is then bound to fail, even at high temperatures, where, for 6 site chains, we get an accuracy of 10% for only.
Panels (d) to (f): Relative error for with respect to the exact results for the same parameters as the upper panels. The darker the colour, the more accurate the approximation is in that regime.
v.2.2 Exact initial system with noninteracting evolution operator
To try and improve the estimate of the work extracted, we shall consider to still implement a noninteracting evolution, but this time starting from the exact manybody initial state. The rationale is that a manybody evolution is in general a more challenging part of the calculation (and hence here it is approximated), while an accurate estimate for the initial state would be more readily available. This approximation is referred to as in table 1.
Indeed this simple approximation leads to strikingly improved accuracy. Results are presented in figure 4: in the upper three panels, and its relative error with respect to the exact results in the lower panels. Parameters are the same as in figure 3.
By comparing to the corresponding exact results in panels (b), (e), and (h) of figure 2 we see that including interactions just through the initial state is sufficient to recover the qualitative (and in great part quantitative) behaviour at low and intermediate temperatures. The greatest improvement is seen in the low temperature, where recaptures the correct work to 1020% for most regimes up to .
At high temperature the qualitative behaviour is not recovered, but, as the work extracted varies only slightly at this temperature, quantitatively the approximation remains overall good, as it reproduces rather well the overall work variation range (compare the colour bar limits of figures 2 and 4). In general the ‘exact + NI’ approximation significantly improves for the minimum extracted average work over the NI value and at all temperatures.
A similar pattern is seen in all the other systems considered, i.e. for 2 and 4 site chains, and for the MI and AEF evolutions (see appendix), demonstrating the scalability and versatility of this approximation. The approximation handles weak to medium correlated systems well in all regimes and temperatures, from adiabatic to sudden quench, and from low to high temperatures.
Panels (d) to (f): Relative error for with respect to the exact results for the same parameters as the upper panels. The darker the colour, the more accurate the approximation is in that regime.
v.2.3 Adiabatic behaviour and the ‘exact + NI’ approximation
The most striking improvement of the ‘exact + NI’ approximation with respect to the NI one, is the recovery of the qualitative behaviour of the average quantum work in the quasiadiabatic regime for low to intermediate temperatures. In this regime the exact is becoming independent from – as to be expected in the adiabatic regime – but strongly depends on , see figure 2. The recovery of this behaviour by the ‘exact + NI’ approximation can be explained as follows.
As discussed in section V.1.1, the dependency of over is contained in the conditional probabilities . In the NI approximation these would read which, would lead, as expected, to a result potentially dependent on but completely independent from . However in the ‘exact + NI’ approximation we have , which could lead to dependency on (through ) as well as on ^{3}^{3}3For the ‘exact + NI’ approximation, the results of calculating using (2) or (6) differ for the terms. The term via (2) can be written as , while the corresponding term via (6) is . The terms give instead the same results..
Let us name the th eigenstate of : then , independent of , and , with . Consider
(8) 
In the adiabatic regime , so that (8) becomes
(9) 
which indeed depends on but not on , as observed.
At high temperatures, inclusion of manybody interactions just via the initial state is a tooweak correction to counter the high thermal energy, which becomes even more dominant than in the exact case: at high temperatures for the ‘exact+NI’ approximation a behaviour closer to the fully NI approximation is recovered (compare figure 4 panel (c) to figure 3 panel (c)). We observe this for all driving potentials (see appendix, panel (c) of figures 8 and 9).
v.2.4 Sudden quench and the ‘exact + NI’ approximation
With respect to the NI approximation, the ‘exact + NI’ approximation recovers the qualitative and, in great part, quantitative exact behaviour for the average quantum work also in the sudden quench limit. This can be seen by comparing results for the small parameter region of the upper panels of figure 4 with the corresponding panels in the central column of figure 2. We give an explanation for this below.
In the quasisuddenquench regime , , with the second term a small correction. By using this and the fact that is a constant, we can approximate (2) as
(10)  
(11) 
In the ‘exact + NI’ approximation, (10) and (11) become
(12)  
(13) 
where we have explicitly stressed which terms on the r.h.s. are to be computed exactly and which ones in a specific approximation. In contrast to a fully NI approximation, (13) is now dependent through , which means that this dependence for is recovered. Similarly, the correction for small (but finite) values in (12) contains a dependence through , and will then improve both qualitatively and quantitatively over NI results, as comparison of figures 2, 3, and 4 demonstrates.
Vi Entropy variation
vi.1 Exact results
Let us now examine the exact results for the entropy, from (3). This quantity corresponds to the heat which the system would have to disperse in the environment to return to thermodynamic equilibrium at the end of the dynamics. Apart from the average quantum work, the other key ingredient for is the freeenergy variation, shown by (4). Since our final Hamiltonians are independent of , the free energy only depends on and .
Figure 5, panels (a)(c), shows the variation of free energy as varies at each of the temperatures considered (green line for low, blue for medium, and red for high temperature) and for each of the driving potentials. is weakly dependent on for high temperatures, while, at intermediate and, even more, at low temperatures, it significantly changes with the interaction strength. This confirms that the system behaves more and more like a noninteracting system as the temperature increases.
The exact entropy variation for all driving potentials and temperatures is shown in figure 5, second to fourth rows: left column for MI, middle for ‘comb’, and right for AEF driving potential. increases as the colour shade becomes lighter; note however that the same shade corresponds to different values in different panels, as the overall entropy range significantly changes according to both temperature and type of driving potential.
The temperature affects the entropy production drastically, compare the extent of ranges between the second and last row of figure 5. This can be understood by comparing energy scales. By the end of the dynamics, our driving potentials reach the maximum energy difference of for comb, for MI, and for AEF potential. For the low temperature and the range of parameters explored, both the interaction strength and the driving potentials can be much bigger than the thermal energy, and so they have a large impact on the system evolution. The system can change quite drastically leading to the possibility of large work extraction and large entropy production. However at high temperature, , the interaction strength and driving potential are, at most, comparable to the thermal energy: the system is not as receptive to being driven, it will remain closer to its thermal state, and hence the energy required to be dispersed to return to equilibrium (i.e. the entropy that we are here considering) will be much less.
In general at the lowest temperature the sudden quench with weakmedium strength coupling parameter region corresponds to very high values, while a dramatic reduction in entropy is seen moving towards the adiabatic regime. Given the correspondence between and the heat to be dispersed at the end to recover equilibrium, it stands to reason that a sudden quench would require a larger dissipation of energy to return to an equilibrium state compared to an adiabatic evolution, which remains closer to an equilibrium state at all times.
We note that systems subject to AEF potential show at low a relatively larger entropy production in the stronglycoupled regime and than systems subject to ‘comb’ and MI potentials. As discussed in section V.1.1, the level of adiabaticity reached for the same value of differs with driving potential. Indeed the dynamic induced by AEF at and is less adiabatic than the ones by MI or ‘comb’, as can be observed by comparing panel (c) to panels (a) and (b) of figure 2. This leads to a larger amount of entropy production occurring with AEF even in this stronglycoupled large region.
vi.2 Noninteracting approximation
We extend the approximations for the extracted work to the entropy production. will be calculated from and by setting in the calculation of the free energy ( values in figure 5, upper panels).
Figure 6 presents the results for 6 sites, ‘comb’ driving potential, and increasing temperature (left to right). The upper panels show the noninteracting entropy variation, and the lower ones the relative difference with the exact entropy variation (the darker the purple, the more accurate the approximation is in that region).
For each given temperature, the entropy is just with an added constant, so, much like the noninteracting work, for all driving potentials this approximation is unable to qualitatively describe the exact entropy produced.
Comparing with the NIwork approximation accuracy, the overall quantitative accuracy of the entropy is in general reduced for all three temperatures, which is to be expected since, on top of the work, we are severely approximating the free energy as well. As increases, we are exploring highly correlated systems, including systems subject to a precursor to the Mottinsulator transition, but these Coulomb correlations are completely neglected by the NI approximation.
For the AEF driving potential, the NI approximation works better then for ‘comb’ drive at all temperatures, giving a 10% accuracy (or better) for at low and intermediate temperatures and for almost all regimes at high temperatures. For the MI driving potential, the accuracy of the results is comparable to the ‘comb’ driving potential.
vi.3 Exact initial state with noninteracting evolution
Let us now see how considering the exact initial state affects the estimate of the entropy production. Initially the entropy is calculated from (3) using and the exact free energy variation: in this approximation, we are already assuming that we can diagonalise the initial Hamiltonian to get the exact initial state, we then make the same assumption for , as this operation would have the same calculation costs/difficulties. This leads to the exact free energy variation. However we note that, with the implementation described above, this approximation could lead to the nonphysical occurrence of negative entropy: in fact the two contributions to the entropy have opposite sign, and one of them (the work) has been approximated, so the occurrence of a negative sign cannot apriori be excluded. We then further impose that to correct for it.
Related results are plotted in figure 7 for 6 sites, ‘comb’ driving potential, and increasing temperature (left to right). By comparing the upper panels of figure 7 to the upper panels of figure 6 and to the corresponding ones in the mid column of figure 5, we note a marked improvement in the qualitative behaviour of the approximation. As with the work, we can see an improvement also in the quantitative results. The high temperature accuracy is very much akin to that of the work, and is accurate within 30% for all regimes. As the temperature decreases, however, the quantitative accuracy also decreases: inaccuracy comes into the entropy calculations through the approximation of the average work, so the regimes of greater/lesser accuracy fairly mirror those of the work.
Results for the other two driving potentials confirm these trends and are shown in the appendix.
Similarly to the ‘comb’ potential, the ‘exact+NI’ approximation with the MI and AEF driving potentials recover to a good extent the qualitative behaviour of for low and intermediate temperatures. For high temperatures the qualitative behaviour is recovered only for . Quantitatively, the areas of worse performance are related to the areas of worse performance for the corresponding , however the approximation performs worse for than for MI, and better for AEF. Overall the approximation improves its quantitative performance with temperature, as it reproduces well the limits of the entropy variation range, and especially so at high temperature.
Vii Conclusion
We have presented a comprehensive study of the extracted average quantum work and entropy variation in complex manybody systems subject to a wealth of driving dynamics and dynamical regimes. We have discussed in details the effects of the interplay of the different energy scales governing the systems – driving potentials, manybody interactions, and thermal energy – on and , and compared results as manybody correlations are turned on up to the strongly coupled regime and as the dynamical regime is continuously changed from sudden quench to nearly adiabatic.
For all driving potentials and low to intermediate temperatures, intermediate to strong Coulomb repulsion decreases work extraction by making the system less responsive to the external drive. At weak Coulomb correlations, more work can be extracted as the system becomes adiabatic, while with stronger manybody interactions the work production becomes independent of the overall driving time much sooner. For the same parameter sets, at high temperatures, work extraction is greatly reduced both in absolute values and in variation range.
At low temperatures the entropy variation presents a behaviour somewhat more dependent on the applied driving potential. We observe a general tendency of lower values for weak and strong Coulomb correlations and intermediate to large ’s, but the onset of an adiabaticlike entropy dynamics varies considerably with the driving potential at intermediate coupling strengths. For a zero temperature, openboundary, finite, homogeneous Hubbard chain, intermediate coupling strengths correspond to the transition between metallic and insulating behaviour (precursor to a Mott insulator transition). The behaviour observed at intermediate coupling strengths may be a signal of how the different driving potentials affect this transition.
Similarly to the extracted average quantum work, entropy production decreases with temperature in both absolute value and range of variation, as the system becomes less and less responsive to the applied potential.
The strong effect of the Coulomb interaction on the extracted average quantum work can be also appreciated by comparing exact results to the corresponding noninteracting approximations. Apart from the obvious independence on for all regimes, this approximation strongly overestimates the minimum work that can be extracted from a system, as it cannot reproduce the freezing of the system dynamics due to strong manybody correlations. On the other hand, the value of the maximum average work extracted is well captured, albeit often it is attributed to the wrong parameter regions.
Including manybody interactions exactly for complex system is often an hopeless task. In this paper we have proposed a relatively simple approximation which relies on been able to provide an accurate approximation for the system’s initial state, while interactions are completely neglected in the dynamics.
We found this approximation to behave surprisingly well. Including interactions just through the initial state recovers the qualitative (and in great part quantitative) behaviour at low and intermediate temperatures for all regimes, driving potentials and chain lengths considered. The greatest improvement is seen at low temperatures, where, for example, with the ‘comb’ driving potential the ‘exact+NI’ approximation reproduces the exact work within 1020% up to very strong interactions () for most regimes. At all temperatures, the ‘exact + NI’ approximation reproduces fairly well the range of variation for the value of the average extracted quantum work and in particular it significantly improves the value of its minimum over the NI estimates.
We demonstrate analytically that, independent of system Hamiltonian and driving potential, including interactions in just the initial state is sufficient to recover at low to intermediate temperatures the characteristic dependency on and in the adiabatic region, which was completely missed by the noninteracting approximation. This is important as this is the region in which the most work can be produced, and the dependency on is very strong.
We also demonstrate analytically that this approximation recovers the characteristic dependence of the extracted work in the small parameter region for all Hamiltonians and driving potentials.
Fully noninteracting estimates of the entropy variation strongly underperform, both qualitatively and quantitatively. However at high temperatures the exact entropy variation range becomes very small and its quantitative noninteracting estimate is in the right ballpark due to the decreased influence of manybody interactions: as a result, at high temperatures even the noninteracting approximation gives reasonable quantitative (but not qualitative) results.
When we extend the ‘exact+NI’ approximation to the entropy variation, we find qualitative improvements similar to the work extraction, with the behaviour for low and intermediate temperatures largely recovered. Quantitatively this approximation significantly improves over the noninteracting approximation, albeit not as strikingly as for the average quantum work. Overall the ‘exact+NI’ approximation improves its quantitative performance with temperature, as it captures well the entropy variation range, including at high temperature. As an example, at high temperature and MI driving potential, this approximation would reproduce exact results within 10% for all parameter range, to be contrasted with the performance by the non interacting approximation which gives such an accuracy only for .
Our results show that, even when taking a very crude approximation for the evolution operator, starting from an accurate initial state is sufficient for greatly improving the estimate of thermodynamical quantities such as the average quantum work and the corresponding entropy variation for the wide range of driving potentials tested, all temperatures, and the great part of interaction and dynamical regimes.
Acknowledgements.
AHS thanks EPSRC for financial support. KZ thanks the Schlumberger Foundation for financial support through the Faculty for the Future program. IDA acknowledges support from the Conselho Nacional de Desenvolvimento Cientfico e Tecnologico (CNPq, Grant: PVE Processo: 401414/20140) and from the Royal Society (Grant no. NA140436). We thank M. Herrera and R. M. Serra for useful discussions in the early stages of this work.Appendix A Sixsite ‘exact+NI’ approximation results for MI and AEF driving potentials: extracted average quantum work
To demonstrate how well the the ‘exact+NI’ approximation fares with other driving potentials, we present here the related MI and AEF driving potential results for 6 site chains and all three temperatures.
Figure 8 shows results for the MI driving potential: this particular dynamic is captured very well by the ‘exact+NI’ approximation for all temperatures and regimes. As with all other driving potentials considered, the lowest temperature is where quantitatively the approximation struggles the most, though its accuracy is, in most regions, within 10% of the exact work, otherwise it is 20%. At higher temperatures the accuracy remains always within 10%. We note that this approximation recovers to high accuracy and at all temperatures the overall range of extractable work (compare colourscale range above panels (a) to (c) in figure 8 with the corresponding ranges in panels (a), (d), and (g) of figure 2). This approximation also recovers most of the work qualitative behaviour, at least at low and intermediate temperature. At high temperature, although the trends of the exact and approximate results are qualitatively different for and , the exact amount of work extracted changes very little with and (panel (g) of figure 2), so, as the quantitative range of the approximated work is close to the exact one, the results in figure 8(c) show high quantitative accuracy.
Panels (d) to (f): Relative error for with respect to the exact results for the same parameters as the upper panels. The darker the colour, the more accurate the approximation is in that regime.
Panels (d) to (f): Relative error for with respect to the exact results for the same parameters as the upper panels. The darker the colour, the more accurate the approximation is in that regime.
Results for and the AEF driving potential are shown in figure 9. Once again the approximation accuracy improves with increasing temperature (panels (d) to (f)), while the qualitative trend is well reproduced for low and intermediate temperatures, but missed at high temperatures for intermediatetolarge and values (compare panels (a) to (c) of figure 9 to panels (c), (f), and (i) of figure 2, respectively). For AEF – and at difference with ‘comb’ and MI driving potentials – the Coulomb repulsion for remains a dominant energy scale at any temperature: in fact even at the driving potential difference between nearby sites remains at most of (see figure 1). This results in a quantitatively poor approximation for even at high temperature.
Appendix B Sixsite ‘exact+NI’ approximation results for MI and AEF driving potentials: entropy variation
Estimates for the entropy variation in the ‘exact+NI’ approximation for MI and AEF driving potentials and sixsite chains are shown in figures 10 and 11, respectively.
Similarly to the ‘comb’ potential, the ‘exact+NI’ approximation with the MI and AEF driving potentials recover to a good extent the qualitatively the behaviour of for low and intermediate temperatures. For high temperatures the qualitative behaviour is recovered only for .
Quantitatively, the areas of worse performance are related to the areas of worse performance for the corresponding , however the approximation performs worse for the entropy than for the extracted work for MI, and better for AEF. Overall the approximation improves its quantitative performance with temperature, as it reproduces well the limits of the entropy variation range, and especially so at high temperature.
References
 Goold et al. [2016] J. Goold, M. Huber, L. del Rio, and P. Skrzypczyk. The role of quantum information in thermodynamics  a topical review. J. Phys. A, 49:143001, 2016.
 Vinjanampathy and Anders [2015] S. Vinjanampathy and J. Anders. Quantum thermodynamics. Contem. Phys., 57:1, 2015.
 Millen and Xuereb [2016] J. Millen and A. Xuereb. Perspective on quantum thermodynamics. New J. Phys., 18:011002, 2016.
 Parrondo et al. [2015] J. M. R. Parrondo, J. M. Horowitz, and T. Sagawa. Thermodynamics of information. Nat. Phys., 11:131, 2015.
 LuizzoScorpo et al. [2016] P. LuizzoScorpo, L. A. Correa, R. Schmidt, and G. Adesso. Thermodynamics of quantum feedback cooling. Entropy, 18:48, 2016.
 Silva [2008] A. Silva. Statistics of the work done on a quantum critical system by quenching a control parameter. Phys. Rev. Lett., 101:120603, 2008.
 Dorner et al. [2012] R. Dorner, J. Goold, C. Cormick, M. Paternostro, and V. Vedral. Emergent thermodynamics in a quenched quantum manybody system. Phys. Rev. Lett., 109:160601, 2012.
 Joshi and Campisi [2013] D. Joshi and M. Campisi. Quantum hertz entropy increase in a quenched spin chain. Eur. Phys. J. B., 86:157, 2013.
 Mascarenhas [2014] E. Mascarenhas. Work and quantum phase transitions: Quantum latency. Phys. Rev. E, 89:062103, 2014.
 Sindona [2014] A. Sindona. Statistics of the work distribution for a quenched fermi gas. New J. Phys., 16:045013, 2014.
 Fusco et al. [2014] L. Fusco et al. Assessing the nonequilibrium thermodynamics in a quenched quantum manybody system via single projective measurements. Phys. Rev. X, 4:031029, 2014.
 Zhong and Tong [2015] M. Zhong and P. Tong. Work done and irreversible entropy production in a suddenly quenched quantum spin chain with asymmetrical excitation spectra. Phys. Rev. E, 91:032137, 2015.
 Eisert et al. [2015] J. Eisert, M. Friesdorf, and C. Gogolin. Quantum manybody systems out of equilibrium. Nat. Phys., 11:124, 2015.
 Bayat et al. [2016] A. Bayat et al. Nonequilibrium critical scaling in quantum thermodynamics. Phys. Rev. B, 93:201106(R), 2016.
 SolanoCarrillo and Millis [2016] E. SolanoCarrillo and A. J. Millis. Theory of entropy production in quantum manybody systems. Phys. Rev. B, 93:224305, 2016.
 Herrera et al. [2017] Marcela Herrera, Roberto M. Serra, and Irene D’Amico. DFTinspired methods for quantum thermodynamics. Scientific Reports, 7:4655, 2017.
 Herrera et al. [2018] M. Herrera, K. Zawadzki, and I. D’Amico. Melting a Hubbard dimer: benchmarks of ‘ALDA’ for quantum thermodynamics. ArXiv eprints, March 2018.
 França and Capelle [2006] V. V. França and K. Capelle. Entanglement of strongly interacting lowdimensional fermions in metallic, superfluid, and antiferromagnetic insulating systems. Phys. Rev. A, 74:042325, Oct 2006. doi: 10.1103/PhysRevA.74.042325.
 França et al. [2012] Vivian V. França, Dominik Hörndlein, and Andreas Buchleitner. FuldeFerrellLarkinOvchinnikov critical polarization in onedimensional fermionic optical lattices. Phys. Rev. A, 86:033622, Sep 2012. doi: 10.1103/PhysRevA.86.033622.
 de Picoli et al. [2018] T. de Picoli, I. D’Amico, and V. V. França. Metricspace approach for distinguishing quantum phase transitions in spinimbalanced systems. Brazilian Journal of Physics, 48(5):472–476, Oct 2018. ISSN 16784448. doi: 10.1007/s1353801805926.
 Essler et al. [2005] F.H.L. Essler, H. Frahm, F. Göhmann, A. Klümper, and V.E. Korepin. The OneDimensional Hubbard Model. Cambridge University Press, 2005. ISBN 9781139441582.
 Coe et al. [2010] J. P. Coe, V. V. França, and I. D’Amico. Hubbard model as an approximation to the entanglement in nanostructures. Phys. Rev. A, 81:052321, May 2010. doi: 10.1103/PhysRevA.81.052321.
 Yang et al. [2011] Shuo Yang, Xin Wang, and S. Das Sarma. Generic Hubbard model description of semiconductor quantumdot spin qubits. Phys. Rev. B, 83:161301, Apr 2011. doi: 10.1103/PhysRevB.83.161301.
 Murmann et al. [2015] Simon Murmann, Andrea Bergschneider, Vincent M. Klinkhamer, Gerhard Zürn, Thomas Lompe, and Selim Jochim. Two fermions in a double well: Exploring a fundamental building block of the hubbard model. Phys. Rev. Lett., 114:080402, Feb 2015. doi: 10.1103/PhysRevLett.114.080402.
 Coe et al. [2011] J. P. Coe, V. V. França, and I. D’Amico. Feasibility of approximating spatial and local entanglement in longrange interacting systems using the extended Hubbard model. EPL (Europhysics Letters), 93(1):10001, jan 2011. doi: 10.1209/02955075/93/10001.
 Brown et al. [2019] Peter T. Brown, Debayan Mitra, Elmer GuardadoSanchez, Reza Nourafkan, Alexis Reymbaut, CharlesDavid Hébert, Simon Bergeron, A.M. S. Tremblay, Jure Kokalj, David A. Huse, Peter Schauß, and Waseem S. Bakr. Bad metallic transport in a cold atom fermihubbard system. Science, 363(6425):379–382, 2019. doi: 10.1126/science.aat4134.
 Nichols et al. [2019] Matthew A. Nichols, Lawrence W. Cheuk, Melih Okan, Thomas R. Hartke, Enrique Mendez, T. Senthil, Ehsan Khatami, Hao Zhang, and Martin W. Zwierlein. Spin transport in a mott insulator of ultracold fermions. Science, 363(6425):383–387, 2019. ISSN 00368075. doi: 10.1126/science.aat4387.
 Carrascal et al. [2015] D. J. Carrascal, J. Ferrer, J. C. Smith, and K. Burke. The hubbard dimer: a density functional case study of a manybody problem. J. Phys. Cond. Mat., 27(39):393001, 2015.
 Fuks and Maitra [2014] J. I. Fuks and N. T. Maitra. Challenging adiabatic timedependent density functional theory with a hubbard dimer: the case of timeresolved longrange charge transfer. Phys. Chem. Chem. Phys., 16(28):14504, 2014.
 Batalhão et al. [2015] T. B. Batalhão, A. M. Souza, R. S. Sarthour, I. S. Oliveira, M. Paternostro, E. Lutz, and R. M. Serra. Irreversibility and the arrow of time in a quenched quantum system. Phys. Rev. Lett., 115:190601, Nov 2015. doi: 10.1103/PhysRevLett.115.190601.
 [31] Note1. We note that for the 2 site chain the only relevant dynamics becomes the ‘Applied Electric Field’ one.
 [32] Note2. See Herrera et al. [2017] and Herrera et al. [2018] for 2site chain examples, and especially Herrera et al. [2017] for 2site noninteracting results.
 [33] Note3. For the ‘exact + NI’ approximation, the results of calculating using (2) or (6) differ for the terms. The term via (2) can be written as , while the corresponding term via (6) is . The terms give instead the same results.