Quantum dissipative adaptation

Dissipative adaptation is a general thermodynamic mechanism that explains self-organization in a broad class of driven classical many-body systems. It establishes how the most likely (adapted) states of a system subjected to a given drive tend to be those following trajectories of highest work absorption, followed by dissipated heat to the reservoir. Here, we extend the dissipative adaptation phenomenon to the quantum realm. We employ a fully-quantized exactly solvable model, where the source of work on a three-level system is a single-photon pulse added to a zero-temperature infinite environment, a scenario that cannot be treated by the classical framework. We find a set of equalities relating adaptation likelihood, absorbed work, heat dissipation and variation of the informational entropy of the environment. Our proof of principle provides the starting point towards a quantum thermodynamics of driven self-organization.


II. INTRODUCTION
When a physical system is simultaneously subjected to both predictable and random energy exchanges, what dictates the likelihood of a given state to be found? From a classical thermodynamic perspective, where energy exchanges are classified as work or heat, the concept of dissipative adaptation has recently been put forward by J. England as the expected answer, at least in the context of driven self-assembly [1]. Qualitatively speaking, dissipative adaptation establishes that, given a certain drive (an external work source), the most adapted (most likely and lasting) states of a fluctuating system tend to be those with a history of exceptionally high work absorption followed by heat dissipation to the environment. Because heat dissipation is an irreversible process, the higher the dissipation, the less likely the reverse trajectory is. In the long run, the system may appear to us as self-organized in this drive-dependent state of highest energy-consuming history. Dissipative adaptation is the recent theoretical development of a long search for the emergence of order from disorder, as inspired by life-like behaviour [2,3]. Examples revealing this general mechanism of * Electronic address: valente.daniel@gmail.com energy-consuming irreversible self-organization span diverse systems, environments, lengths and timescales, as shown both theoretically [4][5][6] and experimentally [7][8][9][10][11][12].
The dissipative adaptation phenomenon has been originally formulated in terms of fluctuation theorems. Fluctuation theorems are equalities relating out-ofequilibrium processes with thermal-equilibrium variables, giving evidence that the fluctuations present in many realizations of a process can provide useful knowledge when summed up. The pioneering example is the so called Jarzynski equality [13], exp(−βW abs ) = exp(−β∆F ), where β = 1/(k B T ) is the inverse temperature, W abs is the work absorbed by the system as described by a time-dependent Hamiltonian, ∆F is the variation in the Helmholtz free energy and the brackets is the ensemble average over realizations of the process, initially departing from thermal equilibrium. The Jarzynski equality has been lately derived by Crooks from what he called a microscopically reversible condition [14], where the forward p i→j (t) and backward p * j→i (τ −t) probabilities for the system to follow paths linking states i and j are related with the amount of heat stochastically dissipated to the environment, Q diss (a functional of the phase-space trajectory and of the driving protocol performing work on the system). p * j→i (τ − t) is computed with the reversed time-dependent protocol. βQ diss relates the statistical irreversibility with the thermodynamic entropy production. Let us call E ij = E j − E i the energy difference between the specific final and initial states. Energy conservation during each stochastic realization, gives the hint on how the work source fuels the dissipative adaptation. The higher the absorbed work, the more heat can be released, hence the more irreversible the path can become. To emphasize this adaptation-work relation, England rearranges Eq. (1) as [1] p i→j (t) p i→k (t) = e −βE kj p * j→i (τ − t) p * k→i (τ − t) e −βW abs ik e −βW abs ij , where the angle brackets denote a weighted average over all microtrajectories with fixed start (i) and end arXiv:2111.08605v1 [quant-ph] 16 Nov 2021 (j, k) points (the coarse graining over microtrajectories describes experimentally accessible states of out-ofequilibrium self-organizing classical many-body systems; see Perunov et al. [15] for further details). Equation (3) establishes the classical theoretical framework behind  dissipative adaptation by showing how a given final state  can be statistically privileged by work consumption. Here, we extend the concept of dissipative adaptation to the quantum realm. Our main goal is to test the robustness of the dissipative adaptation concept beyond its original theoretical framework discussed above. From a technical viewpoint, at vanishing temperatures (β → ∞), where quantum fluctuations and correlations usually prevail, Eqs. (1) and (3) are ill defined. We employ a systemplus-reservoir approach to derive the exact equations of motion of a three-level lambda (Λ) system driven by a single-photon pulse added to a zero-temperature environment. We find that the most adapted (self-organized) quantum state of the lambda system is indeed that with a history of highest work absorption followed by maximal heat dissipation, thus characterizing a dissipative adaptation. As a consequence of our work, we establish the starting point of a quantum thermodynamics of driven self-organization, so far unexplored to the best of our knowledge. We hope that the notion of a quantum dissipative adaptation may also provide fresh insights to quantum biology [16], not only because adaptation and self-organization are concepts inspired in life-like behaviour, but also because our results may find applications to discussions on quantum signatures in photosynthesis [17][18][19][20].

III. RESULTS
Self-organized quantum state. To achieve our main goal, we look for the most elementary scenario where the quantum state of a certain physical system irreversibly self-organizes, as induced by the absorption of energy from an external drive, the excess of which is dissipated to the environment. We choose a three-level system in Λ configuration, labeled as |a , |b and |e (being |e the most excited state, respectively with transition frequencies ω a,b ) (see Fig.1). To keep the model exactly solvable for both the system and the environment, we assume that the drive source is a single-photon pulse added to the vacuum state of an infinite environment at zero temperature, T = 0. The environment induces spontaneous emission rates, Γ a and Γ b . Our main results are that showing that adaptation likelihood at long times, p a→b (t → ∞), is linearly proportional to the average absorbed work W abs a from a single-photon pulse of arbi-trary shape, resonant with ω a . We call attention to the fact that this result is not immediately expected, since (i) the excitation probability is minimized (p a→e (t) 1, at all times) when the work absorption is maximized (p a→b (∞) → 1), so the process cannot be thought of as a simple absorption-plus-emission picture, and (ii) the final self-organized state is not restrained to be the ground state, which amounts to say that W abs a does not depend on ω b , so the work is not related to (ω a − ω b ), in particular. The second line in Eqs.(4) characterizes the irreversibility of the process. The absorbed work is partly dissipated to the environment in the form of heat, Q diss a , satisfying energy conservation H S here is the Hamiltonian of the system. Finally, we find an exact expression for the informational entropy of the environment at long times, S E (∞), as a function of the average dissipated heat, Q diss a . The entropy analysis here provides us with both a clearer physical picture of the process and an additional signature of the dissipative adaptation.
|ai |bi Figure 1: Quantum dissipative adaptation in a driven self-organized quantum state. a Three-level system in lambda configuration is described by the density operator ρS(t) at time t, starting at t = 0. The energy eigenstates are |a , |b and |e (blue horizontal bars). The transition frequencies are ω a,b (full gray arrows). The environmentinduced spontaneous emission rates are Γ a,b (dashed gray arrows). The initial state is a mixture between |a and |b , b |b b| (smaller black dots). b A singlephoton pulse |1a (the work source) drives the lambda system dynamics, inducing the time-dependent transition probability p a→b (t) from |a to |b (the full black horizontal arrow represents the forward dynamics from time t = 0 to t → ∞). The backward probability (with a time-reversed pulse), p * b→a (t) = 0, vanishes at zero temperature (the dashed gray horizontal arrow represents the prohibited time-reversed transition). c The driven lambda system undergoes an ideal irreversible self-organizing dynamics (p a→b (∞) = 1 and p * b→a (t) = 0, in panel b), so that the asymptotic state is pure, ρS(∞) = |b b| (larger black dot), conditioned to maximizing the work absorbed and the heat dissipated in the |a to |b transition.
Let us suppose that our generic three-level system in Λ configuration, with lowest-energy states |a and |b and excited state |e , is initially in a non-driven steady state, in contact with the environment at temperature T = 0. In this case, at precisely T = 0, the steady state is not uniquely defined, even for non-degenerated cases. We can choose its initial quantum state to be a mixture of the lowest-energy states, as described by the density operator where p (0) a,b depend on the preparation scheme. To give a concrete example, in the preparation by means of a spontaneous emission starting at |e , one has p (0) a,b = Γ a,b /(Γ a + Γ b ). Now we look for the most elementary out-of-equilibrium stochastic process that drives the system from this (generally mixed) initial state into a final (ideally pure) target state, let us say into state |b , In order to guarantee that the process is irreversible, in the light of Crooks condition, we should also apply the time-reversed drive on the system departing from state ρ S (∞) = |b b| and find that ρ * We employ a system-plus-reservoir approach, where we assume a global time-independent Hamiltonian of the system and its environment, As we show in what follows, we find a self-organized quantum state in the well-known dipolar model of lightmatter interaction in the rotating-wave approximation [21,22], Here, σ k = |k e| (for k = a, b) and H.c. is the Hermitian conjugate. We consider a continuum of frequencies, ω → dω ω , with density of modes ω . Modes {a ω } and {b ω } are the quantized field modes respectively interacting with the transitions |a to |e and |b to |e . The continuum of frequencies allows us to employ a Wigner-Weisskopf approximation to obtain dissipation rates Γ a = 2πg 2 a ωa and Γ where δ ab = (ω a − ω b ) is the energy difference between states |b and |a . It is worth emphasizing that our main results in this paper are independent of δ ab . To keep the model exactly solvable for both the system and the environment, we choose the drive as provided by a propagating pulse containing a single photon. We choose the photon to initially populate only the continuum of modes {a ω }, so the vacuum state of {b ω } allows us to avoid depleting our target state |b . The initial state of the field is ω is the vacuum state of all the field modes. The initial state of the Λ system is the mixed ρ S (0) given in Eq. (6). The global quantum state is given by where U = exp(−iHt/ ) for the (time-independent) global Hamiltonian H [Eq. (8)]. The quantum states of the system and the environment are obtained by the par- We obtain analytical expressions for the probabilities p k (t) = k|ρ S (t)|k , for k = a, b, e. Equation (6) allows us to write the reduced state of the system as ρ where n = a, b, and p n→k (t) ≡ k|tr E U |n, 1 a n, 1 a |U † |k .
Equations (13) and (14), which follow directly from the mixed-state structure of the initial state, unravel a close analogy with the notation used in Eq.(1) for the microscopically reversible condition. Equation (13) has no relation with the degree of Markovianity in the dynamics of the reduced state ρ S (t), though. On the contrary, Valente et al. [23] characterize non-Markovianity in a quite similar scenario. The photon initially at modes {a ω } does not interact with the three-level system initially at |b [as derived from Eqs.(6), (9) and (12)], so we find that p b→a (t) = p b→e (t) = 0 and p b→b (t) = 1. Because the transition probabilities vanish identically, regardless of the initial pulse shape, the backwards probability also vanishes, We have defined the reverse drive protocol here as the mirrored shape of the initial pulse (further discussed below). Put simply, we revert only the drive, not the entire universe. The final global state of the system plus the environment is obviously reversible in our model since it is governed by the global unitary operator U in Eq. (12).
Equation (15) explains the second line in Eqs. (4). We now compute where |ξ(t) ≡ U |a, 1 a . Since our H conserves the total number of excitations, we can restrict our model to the one-excitation subspace, and similarly with p a→a (t). Without loss of generality, we have defined a one-dimensional real-space representation for the amplitudes, φ n (z, t) ≡ ω φ n ω (t) exp(ik ω z), which characterizes the pulse shape. We have also employed a linear dispersion relation, ω = ck ω , and approximated the density of modes by a constant, ω ≈ . Note that, in Eq.(9), we have implicitly assumed the threelevel system to be positioned at z S = 0; otherwise, a phase term such as exp(ik ω z S ) should have been included in the sum over modes.
To go a step further, as to obtain explicit expressions for the transition probabilities, we solve the Schrödinger equation for |ξ(t) (see Methods). Our intention with keeping modes {b ω } initially in the vacuum state, as we did in Eq. (11), was to minimize excitations promoting the unwanted backward (|b → |a ) transitions. Now that we have defined the amplitudes φ n (z, t), we see that Eq.(11) implies φ b (z, 0) = 0, which we combine with (34) and substitute in (17). After a change of variables, we find that Although the main results in this paper do not depend on the initial pulse shape in modes {a ω } (i.e., on the choice of |φ a (z, 0)|), it is worth working out an explicit example.
To that end, we now set φ a (z, 0) to have an exponential envelope profile of linewidth ∆ and a central frequency ω L (see Eqs. (36) and (37) in the Methods), as typical in spontaneous emission. We are particularly interested in the resonant condition ω L = ω a (see the Methods for the more general solution). We find, in the monochromatic limit, ∆ Γ a + Γ b , and at long times t ∆ −1 , that When Γ a = Γ b , we have that p a→b (∞) = 1. Equation (19) reveals the ideal driven self-organization we have been looking for. As stated before, the self-organization in our model does not depend on δ ab . As a final step, we combine Eq. (19) and p b→b (t) = 1 to see that, in the ideal (monochromatic and resonant) regime, that is, ρ S (∞) = |b b| and ρ * S (∞ − t) = |b b| for all times t. Hence, in the t → ∞ limit, ρ * S (∞ − t) = ρ S (0), as was to be shown.
Energetics of the self-organization. We now need to verify that the self-organized quantum state we have found can indeed be classified as dissipative adaptation. We shall find that ideal self-organization costs maximal work absorption, followed by maximal dissipation. This is not an obvious relation because in the ideal selforganization (which takes place in the monochromatic limit, as we have shown above), the excitation probability is minimized (rather than maximized), For instance, in the case of an incoming pulse with exponential profile (as used in Eq. (19)), |ψ(t) Extremely low excitation probabilities may leave the false impression that no energy is absorbed neither dissipated at all. To address this issue, we need to resolve energy transfer into work and heat in our model. As mentioned earlier, work and heat can be regarded as the predictable versus the random energy exchanges, respectively. We restrict ourselves once again to the resonant case, ω L = ω a , keeping in mind however an arbitrary pulse shape. The average energy of the Λ system driven by this resonant photon is given by The resonant condition here avoids dynamic Stark shifts that would otherwise bring extra (dispersive) energetic contributions from time-dependent frequencies, as shown by Valente et al. [24]. To be more precise, dispersive (refractive) energetic contributions depend on the average interaction energy [24], which vanishes at resonance (ω L = ω a implies that H I (t) ≡ tr[ρ(t)H I ] = 0; see details in the Methods). That justifies why we have used only H S in Eq. (22). As inspired by Eq.(3), we are interested in the work performed on the system during the dynamical transition starting from |a and arriving at |b . Therefore, we now consider ρ S (0) = |a a|, for which p e (t) = |ψ(t)| 2 . From the Schrödinger equation (see Eq.(32) in the Methods), we find that where Re [•] stands for the real part. Equation (23) clearly shows us that the excited state of the Λ system is governed by a (predictable) drive-dependent term, related to φ a (−ct, 0), and a (random) spontaneous emission term, proportional to Γ a + Γ b . This motivates us to define the total average absorbed work W abs a and the total average absorbed heat Q abs a in the dynamics starting at ρ S (0) = |a a| as where (25) and Our definition of work, Eq. (25), can also be written as W abs a = ∞ 0 (∂ t d(t)) E in (t) dt (in the Heisenberg picture and in the rotating-wave approximation), revealing a more clear link to our classical notion of work. Here, d(t) = U † d U is the dipole operator, where d = k d ek σ k +H.c., the driving field operator at the system's position is E in (t) = ω i ω (a ω + b ω ) exp (−iωt) + H.c., the coupling is g a = d ea ωa / , and the average is calculated at the initial state |a, 1 a . The Heisenberg picture also explains why work can be finite even though the average driving field is precisely zero at the single-photon state (i.e., 1 a |E in (t)|1 a = 0). This shows how the dissipative adaptation can be robust to a phase-incoherent work source, in contrast to the classical forces of welldefined phases used in [4] (in our model, a semiclassical driving field with a well-defined phase would correspond to an initial coherent, or Glauber, state |α = ω |α ω , where a ω |α ω = α ω |α ω ). In Eq.(26), Γ a,b ∝ g 2 a,b reveals that the heat is related with the variance of the interaction energy, explaining why the vanishing average interaction energy does not hinder energy exchange in the form of heat and confirming its stochastic nature. The energy conservation in Eq.(5) results, of course, from defining Q abs a ≡ − Q diss a in Eq. (24). Finally, H S (0) = 0 and H S (∞) = p a→b (∞) δ ab . We remind that, during the dynamical transition from |a to |b , we have p e (t) = |ψ(t)| 2 , allowing us to establish an exact adaptation-energy relation between Eqs. (18) and (26). We finally find our quantum adaptation-work relation, valid for a photon of arbitrary pulse shape and resonant with ω a , as well as for arbitrary δ ab , as stated in Eqs. (4).
In the case of the exponential pulse used in Eq. (19), for instance, we find that W abs a = ω a 4Γ a /(Γ a + Γ b + ∆), maximized in the monochromatic regime (∆ Γ a + Γ b ). If in addition Γ a = Γ b , we get W abs a = 2 ω a , which is twice the initial average energy contained in the singlephoton pulse. This counterintuitive result reinforces the notion that work is the amount of energy transferred during a process, rather than the average energy stored in a system at a given time (after all, unitary evolutions with time-independent Hamiltonians conserve energy; here, ∂ t H(t) = 0 implies that H S (t) + H I (t) + H E (t) = p b (0) δ ab + ω L , so at resonance and for p b (0) = 0 we find H S (t) ≤ ω L ). Equation (27) is our key signature of a quantum dissipative adaptation. The seeming low-excitation issue (|ψ(t)| 2 1, in the monochromatic limit) has finally been clarified, given that the time integral of |ψ(t)| 2 (in Eq. (18)) is not only finite, but also linearly proportional to the work required for the quantum dissipative adaptation to take place.
Entropy in the self-organization. In the classical formulation of dissipative adaptation, entropy production plays a key role in establishing the connection between energy transfer and statistical irreversibility, requiring no detailed knowledge on the state of the environment. Here, irreversibility is readily characterized by the asymmetry between the forward, p a→b (t), and the backward, p * b→a (t), processes. Nevertheless, we have the advantage that we keep the full description of the quantum state of the system plus the environment, ρ(t) (at the expense of a greater degree of generality in our global Hamiltonian H). With ρ E (t) = tr S [ρ(t)] at hands (see Methods), we now seek to describe what happens to the environment during and after the drive interacts with the three-level system. We calculate the exact expression for the von Neumann entropy, See Eqs.(45) and (46) in the Methods for the analytic expression. We have found that our classical intuition, namely, that better adaptation produces more entropy in the environment, can be found by an appropriate distinction between the classical and the quantum contributions to S E . The idea behind this distinction is the following. Let us first suppose that the system is initially at |a . A highly monochromatic incoming photon will fully induce transition from state |a to |b . Hence, the outgoing photon will be detected at modes {b ω }, as well. Now, by considering an initially mixed state of the three-level system as given by (1/2)|a a|+(1/2)|b b|, a highly monochromatic incoming photon will have 1/2 probability of leaving at {a ω } (in the cases where it encounters the system at |b ) and 1/2 of leaving at {b ω } (in the cases where it encounters the system at |a ), so the final global state would be mixed and separable, |1 b and |1 free a have been defined in Eqs.(43) and (44) below (see Methods). This is what we call the classical contribution to the final mixed state of the field: the initial mixture of the system is fully transferred to the environment.
Let us now consider again that the system is initially at |a (i.e., p (0) a = 1). However, let us assume that the linewidth of the pulse is of the order of the dissipation rates of the three-level system. In that case, the final state of the global system becomes entangled, N k being the average number of photons at modes k (see Eq.(42) in the Methods). Therefore, the quantum state of the environment in this case, ρ E (∞) = tr S [|ξ(∞) ξ(∞)|], is also mixed between modes {a ω } and {b ω }. However, the mixture in this case arises from a sustained system-environment quantum entanglement rather than from the statistical mixture in the initial state of the system. To unravel these two entropy contributions, we define the classical contribution to the environment entropy as  This monotonic behaviour shows that the more organized is the three-level system, the higher is the classical contribution to the entropy of the environment at long times t → ∞. Besides having derived S c E (∞) as a function of p a→b (∞), we have also expressed p a→b (∞) as a function of the average dissipated heat in the |a to |b transition, Q diss a , with the help of Eq. (27). This establishes the function S E (∞) vs. Q diss a , as mentioned earlier. Additionally, we have derived S c E (∞) as a function of Q diss a . This is even more meaningful, because we have found that the function S c E (∞) vs. Q diss a is monotonic. To provide an example illustrating this monotonicity, let us take the degenerate case (δ ab = 0) with equal dissipation rates (Γ a = Γ b ). In that case, we have that p a→b (∞) = Q diss a /(2 ω a ), showing that S c E (∞) vs. Q diss a is monotonic, since S c E (∞) vs. p a→b (∞) is monotonic as well (Eq.(50) in the Methods provides the more general relation). The monotonicity here strengthens the signature of dissipative adaptation. Namely, maximal adaptation (irreversible self-organization) not only costs maximal work absorption (as we have shown in the energetics analysis), but also maximizes the dissipated heat which, in turn, maximizes the classical contribution to the environment entropy. , increases monotonically as a function of the transition probability from |a to |b , p a→b (∞) (solid black), providing an additional signature of the dissipative adaptation. The environment entropy, SE(∞), presents a non-monotonic behaviour as a function of p a→b (∞) (dashed blue), due to a sustained systemenvironment quantum entanglement that can contribute to the entropy production of the environment. In both curves, we refer to the asymptotic state t → ∞ of the dynamics. We use Γa = Γ b as the environment-induced spontaneous emission rates and p b |b b| (the curves do not qualitatively depend on these choices). The dotted gray line is ln (2). Note that both S c E (∞) as a function of p a→b (∞) and SE(∞) as a function of p a→b (∞) do not explicitly depend on the photon pulse shape, though p a→b (t) does; e.g., attaining p a→b (∞) → 1 requires an extremely monochromatic (long) photon, whereas the p a→b (∞) → 0 limit is attained in the extremely broadband (short) pulse regime.

IV. DISCUSSION
Our results establish the quantum dissipative adaptation underlying the driven self-organization of a quantum state, going beyond the classical formulation. We have explored an elementary fully-quantized model, exactly solvable for both the system and the environment, where the irreversibly self-organized quantum state of a three-level system is due to the work absorption from a single-photon pulse, part of which is dissipated to the environment (as shown in Eqs. (4) and (5)). The irreversibility of this self-organization became clear from the asymmetry in the transition probabilities in the forward and backward processes. Finally, with the purpose of providing an additional signature of dissipative adapta-tion, we have analytically investigated the environment entropy. We have found that the classical contribution to the environment entropy variation is a monotonic function of p a→b (∞) (as illustrated in Fig.2), showing that maximizing adaptation not only requires maximal work absorption, but also leads to maximal increase in the classical contribution to the environment entropy, due to maximal dissipated heat.
We remind that the meaning of this increase in the environment entropy is that of a statistical mixing between the field modes {a ω } and {b ω }, as we have analysed in Eq. (29). Our model's dynamics do not lead to the thermalization of the environment, which always remains out of equilibrium. If we wish to take the analogy further, so as to give a thermodynamic meaning to the entropy increase, we will have to add extra ingredients. In practice, there could be for instance a slight increase in the local temperature of the environment in the vicinity of the Λ system. To describe that kind of effect, the model should consider a finite heat capacity environment (in contrast to our infinite-size environment) and some auxiliary light-matter coupling mechanism that could effectively create an interaction between the many frequency modes of the light field. Such features provide means for the environment state to eventually approximate a Gibbs state at some new finite temperature T > T = 0 (producing even more entropy than that we have calculated from our model), in the spirit of Timpanaro et al. [25].
Before concluding, we would like to point out how our model can be significant to quantum many-body systems. Our intention is to indicate how broadly applicable the concept of quantum dissipative adaptation may become. Firstly, we notice that the Λ structure of energy levels can arise from the quantization of collective degrees of freedom describing many interacting atoms and electrons, as happens in the so-called artificial atoms (e.g., electron-hole pairs in semiconductor quantum dots [26] and the quantized magnetic flux in superconducting rings [27]). Lodahl et al. [26] and Gu et al. [27] also discuss how these artificial atoms can be driven by single-photon pulses propagating in one-dimensional waveguides, building the closest possible scenario to that in our model.
We can also envision the quantum dissipative adaptation in ensembles of non-interacting atoms or spins. That is more easily seen when we realize that our idea of a driven self-organized quantum state is notably similar to the dynamics induced by stimulated Raman adiabatic passage (STIRAP) [28], a versatile and robust technique that has been performed, e.g., in ultracold gases, in doped crystals and in nitrogen-vacancy centers. STI-RAP consists of an efficient population transfer between two discrete quantum states of an ensemble of emitters (usually the lowest levels of Λ systems) by coherently coupling them with two radiation fields (well-controlled classical pulses) through an unpopulated intermediate state. The connection we have in mind between the driven self-organization provided by STIRAP and dissipative adaptation becomes more evident in light of the recent proposal for using STIRAP as a tool for spectral hole burning (SHB) in inhomogeneously broadened systems [29]. The reason is that the mechanism behind standard SHB [30] turns out to be precisely that of classical dissipative adaptation, as discussed by Kedia et al. [31]. Namely, those dipoles that get excited by the resonant drive (whose frequency can be swept on demand) can become irreversibly trapped in dark states (at sufficiently low temperatures). The difference concerning the newly proposed STIRAP-based SHB [29], or the single-photon pulse that we have studied here, is the quantum coherent nature of the process. Understanding how STIRAP-based SHB depends on temperature seems a valuable opportunity for widening the concept of quantum dissipative adaptation (similarly, Ropp et al. [10] show a groundbreaking experiment on temperaturedependent dissipative self-organization in optical space).
Optomechanical nanoresonators [32] also hold the promise of displaying some kind of quantum dissipative adaptation in the lines of our results. We glean this notion from Kedia et al. [31], where signatures of dissipative adaptation are shown in disordered networks of classical bistable springs. In the mechanical nanoresonators by Yeo et al. [32], bistability can arise from a strain-mediated coupling between the centre of mass of an oscillating nanowire and the quantum state of a single semiconductor quantum dot embedded therein. This kind of coupling between an optically controllable microscopic degree of freedom (within the quantum dot) and a mesoscopic degree of freedom (the nanoresonator centre of mass) opens appealing perspectives for studying dissipative adaptation at a quantum-classical boundary of a many-body system. As a last word on resonators, it seems relevant investigating whether the vibrationassisted exciton transport found in photosynthetic lightharvesting antennae [17] could be related with a quantum dissipative adaptation.
To conclude, the quantum dissipative adaptation we have found can be regarded as a proof of principle, in need of generalizations towards many directions. As a first example, it would be worth investigating dissipative adaptation in larger Hilbert spaces (i.e., in a vast energy landscape, in the words of [1,4]). Multistability makes self-organization in classical many-body systems a fascinating problem [1,4,31]. Curiously, the multistability from complex classical systems is reminiscent of our zerotemperature model. Among the huge amount of available "non-organized" stable states that our Λ system can occupy, only a single (exceptional) state is populated when the work source is optimal. In other words, although we find an infinite number of (pure or mixed) combinations of states |a and |b that are equally stationary in the nondriven zero-temperature environment, the system always ends up in the rare state |a (irrespective of the infinitely many possible initial states) once it has been driven by the suitable pulse (which maximizes the absorbed work in the |a to |b transition).
It also remains to be investigated how other sources of work and finite temperatures affect dissipative adaptation in quantum systems. In our Λ system, for instance, we expect that the forward (|a to |b ) transition will keep being favored with respect to the backward (|b to |a ) transition at finite temperatures, as long as the work source is strong enough. Also, the time dependence of this asymmetry should be transient or stationary, depending on the work source type (pulsed or continuous). In quantum thermodynamics, the influence of the environment temperature on entropy production at the quantum regime, on the one hand, has been the subject of many recent studies [25,[33][34][35][36]. The focus on modelindependent aspects is justified not only due to the interest in generality, but also due to the closest possible analogies with classical fluctuation theorems and with the Landauer erasure principle (that establishes how information erasure is connected to thermodynamic entropy production due to heat dissipation). On the other hand, Manzano et al. [37], for instance, provide an alternative formalism for quantum fluctuation theorems that go beyond thermal-equilibrium states of the environment. It may be the case that the quantum fluctuation theorems developed in the papers above (and in the references therein) present fruitful methods in the search for a general quantum theory of dissipative adaptation. In summary, we provide the starting point towards a quantum thermodynamics of driven self-organization.
Entropy of the environment. The quantum state of the environment in our model is obtained from the global initial state ρ(0) = (p (0) a |a a| + p (0) b |b b|) ⊗ |1 a 1 a |. We obtain that ρ E (t) = p (0) a |ψ(t)| 2 |0 0| + p (0) a N b (t)|1 b 1 b | + p (0) a N a (t)|1 a 1 a | + p We have defined N k (t) ≡ ω |φ k ω (t)| 2 , with and To explicitly calculate the von Neumann entropy, we need to find the eigenvalues λ j of ρ E . The exact diagonalization of ρ E gives us four non-zero eigenvalues, namely, λ 1 = p (0) a |ψ(t)| 2 , λ 2 = p (0) a N b (t), and Terms N a and N b can be interpreted as the average number of photons in each continuum of modes. Their mathematical expressions, however, are invariably connected with the three-level system probabilities, namely, N k (t) = p a→k (t).
We calculate the overlap between the states representing the free propagation and the reemitted photon at mode a in the real-space representation, where φ a,free (z, t) ≡ φ a (z − ct, 0). With the aid of (32) and (33), and using integration by parts, we find that valid at long times t → ∞. Eq.(49) provides the core connection between the overlap and the transition probability at long times that we needed. Additionally, N a (∞) = 1 − N b (∞) = 1 − p a→b (∞) and λ 1 (∞) = 0. Since Eqs. (27) and (24) link p a→b (∞) to Q diss a , we have thus analytically established the function S E (∞) vs. Q diss a , valid at long times t → ∞, under the resonant condition ω L = ω a , for any arbitrary pulse shape. Finally, the classical contribution to the entropy is S c E ≡ S E − p for which we also have the analytic solution.

VI. DATA AVAILABILITY
Data sharing not applicable to this article as no datasets were generated or analysed during the current study.

VII. CODE AVAILABILITY
The code used to produce the figure in this article is available from the corresponding author upon request.
VIII. AUTHOR CONTRIBUTIONS D. V., F. B. and T. W. contributed to all aspects of the research, with the leading input of D. V.