Superoperator Master Equations and Effective Dynamics

We developed the projection method to derive an analog of the quantum master equation for propagators rather than density matrices themselves. As these propagators are superoperators, we call them superoperator master equations. Furthermore, as the projector maps superoperators to superoperators, we call it a hyperprojector. We gave general perturbative expansions for generators of the weak coupling superoperator master equation and the stroboscopic limit superoperator master equation. After that, we considered a particular example of a hyperprojector that is the infinite time average of unitary dynamics. We call it the averaging hyperprojector. We discussed the properties of this hyperprojector and its physical meaning. Then, we illustrated our general second order superoperator master equations arising in the weak coupling limit and the stroboscopic limit, taking the averaging hyperprojector as an example, which we call effective dynamics. We discussed some properties of these superoperator master equations, in particular, the entropy increase.


Introduction
The Nakajima-Zwanzig projection approach [1,2] is wide-spread for the derivation of master equations in statistical physics.It can be used both to derive Nakajima-Zwanzig integro-differential master equations or time-convolutionless master equations [3].In this work, we are mostly interested in time-convolutionless master equations, as they are widely used in open quantum systems theory [4][5][6][7].Usually, the projectors are applied to density matrix dynamics to obtain the master equation for such projected density matrices.The difference in our approach is that we apply projectors to dynamical maps and derive master equations for the projected dynamical maps.Therefore, we refer to these equations as superoperator master equations (because these dynamical maps are superoperators) and to the projectors as hyperprojectors (because they map superoperators to superoperators).We apply this technique to the hyperprojector of averaging with respect to free dynamics to obtain an effective dynamical map.Thus, such a map can be considered a dynamical analog of the effective Gibbs state introduced in [8].
We have already used a special case of the hyperprojector technique before to derive effective Heisenberg equations for a specific model with a quadratic fermionic Hamiltonian [9].However, we have neither discussed a general approach nor given a systematic perturbative expansion for the time-convolutionless weak coupling master equations.Therefore, in Section 2, we describe a general setup of the hyperprojector technique for general dynamical maps for some linear equations.For this purpose, we first briefly summarize the usual projection approach for deriving the time-convolutionless master equations.Then, we introduce our approach in the same notation to emphasize the similarities and differences.In particular, by choosing identity preserving hyperprojectors, the superoperator master equation becomes homogeneous for arbitrary dynamical maps.In the case of the usual projectors, it is not possible for arbitrary density matrices, which leads to non-homogeneous terms in the usual master equations [10], and there are some approaches [11] that are trying to avoid presence of these terms to obtain better convergence.After that, we gave a systematic perturbative expansion for the generators of time-convolutionless weak coupling master equations.In Section 2.1, we focus on the second order superoperator master equation and consider in more detail the case when the time-dependent equation before projection arises in the interaction picture for some initial time-independent Liouvillian.In Section 2.2, we briefly discuss a superoperator analog of the Nakajima-Zwanzig equation as it is also used in open quantum systems [12].
In Section 3, we discuss a superoperator analog of a master equation that occurs in the stroboscopic limit of repeated non-selective measurements [13].In particular, we discuss that actually there are several such analogs with different operational meanings.In Section 3.1, we focus on the second order master equation for the stroboscopic limit.
In Section 4, we introduce the averaging hyperprojector and summarize its properties.In Section 4.1, we discuss that its physical meaning is the transformation of dynamical maps as a result of the non-selective measurement of increments of the energy between two time moments.In particular, we show that the averaging hyperprojector is a map from quantum channels to quantum channels, i.e., a supermap or superchannel [14,15], which was discussed in the literature from different perspectives [16][17][18][19].We also show that the averaging hyperprojector maps the unitary channel to the entropy non-decreasing channel.In Section 4.2, we show that the averaging hyperprojector can define the effective dynamics of fast observables, which were introduced in [8].
In Section 5, we consider the second order superoperator weak coupling master equation in the special case of the averaging hyperprojector.Using the expansion of the interaction Hamiltonian as a linear combination of eigenoperators of the commutator with the free Hamiltonian, we give more explicit representation of the second order generator.It is not necessarily of the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) form [20,21], but it is represented in the GKSL-like form.As usual, derivations of the Markovian master equations [22,23] and their corrections [24][25][26][27] are based not only on the weak coupling limit but on the Bogolubov-van Hove scaling [28,29]; as well, in Section 5.1 we discuss it for superoperator master equations.We also discuss its connection with algebraic perturbation theory [30][31][32][33][34].
In Section 6, we discuss the second order superoperator stroboscopic limit master equation for the averaging hyperprojector.In particular, we show that the entropy is non-decreasing in such a case.
In Conclusions, we summarize our results and discuss the directions of the further study.

Superoperator Weak Coupling Master Equation and Hyperprojector Technique
Let us briefly summarize a standard approach to the projection-based derivation of the quantum master equation.Let us consider a linear equation for a density matrix of the general form For the typical open quantum systems setup, for example, it is usually only the Liouvillevon Neumann equation in the interaction picture for a system and a reservoir, but one can consider a more general setup of composite open quantum systems as well [35][36][37][38][39][40].Furthermore, it is assumed that we are not interested in the full dynamics of ρ(t) but take an interest only in the projected dynamics P ρ(t), where P is a projector (idempotent), i.e., P 2 = P.
There are several reasons why we are interested only in the projected dynamics.The first reason is the experimental accessibility.For example, in the open quantum systems setup, it is usually assumed that we have experimental access to open quantum systems only, but we have no experimental access to the heat bath degrees of freedom.We mostly follow this reason in our work.The second reason lies in simplifications of the projected dynamics compared to the full dynamics.For open quantum systems, typically the dynamics of the system only are simpler than the whole dynamics of the system and heat bath (at least after Markovian timescale separation).The third reason is that the given projector absorbs the physical assumptions in its structure.For the standard open quantum systems setup, it is an assumption of the approximately factorized state of the system and heat bath [41], (Paragraph 3.3.1).The fourth reason is that the projected dynamics are interesting for particular quantum technological applications [38,39].Sometimes, even if we are experimentally interested in some degrees in addition to the projected ones, the projected dynamics can be an auxiliary computational tool [42].
Under weak coupling assumptions, we can write a master equation of the form [41], (Subsection 9.2.1) where Q = I − P. We deal with time-convolutionless master equations, but in statistical physics, Nakajima-Zwanzig integro-differential equations could be usually derived.The inhomogeneous term I(t)Qρ(t 0 ) equals zero for the initial condition consistent with the projector, i.e., P ρ(t 0 ) = ρ(t 0 ).For open quantum system theory, the Argyres-Kelley [43] projector P = Tr S ( • ) ⊗ ρ B is typically used, and the consistency of the initial condition with the projector means that the initial state is factorized in this case.However, in general, Equation (2) has inhomogeneous terms [10].Nevertheless, the adapted projection technique [11] can be used to obtain the homogeneous master equations.This technique assumes that the initial state can be approximated as a finite linear combination of the factorized matrices, and a separate homogeneous master equation is obtained for each factorized matrix.
Let us remark that the solution of (1) can always be written in the form where Φ t t 0 is a linear map, which is simply a Cauchy matrix in terms of linear ordinary differential equations theory, or a propagator, satisfying the propagator equation It is a unitary map for all t and t 0 if (1) is the Liouville-von Neumann equation and a completely positive map if (1) is the GKSL equation.The propagator Φ t t 0 satisfies the equation d dt which has the same form as Equation ( 1), but now it is an equation for superoperator Φ t t 0 rather than operator ρ(t), and it has the fixed initial condition while the initial condition ρ(t 0 ) for ρ(t) can be an arbitrary density matrix.Thus, it is natural to make the same procedure for Φ t t 0 , deriving some analog of a master equation for the projected dynamics P(Φ t t 0 ).Here, P is now a projector that maps superoperators to superoperators, so we will call it the "hyperprojector".Due to the initial condition (4), in contrast with master Equation (2), we can avoid inhomogeneity in the equation for P(Φ t t 0 ) for all allowed Φ t t 0 , if one assumes In such a case, we obtain the following superoperator analog of master Equation (2) In this section, we discuss the general setup, trying to emphasize the formal similarities and differences between operator and superoperator master equations.The general physical reasons to consider the projected dynamics are similar to the case of the usual projection superoperators mentioned above.In Section 4, we consider the specific hyperprojector and discuss concrete physical setups that lead to these hyperprojectors arising.Namely, in Section 4.2, we show that they naturally arise when one measures the correlation functions of the fast observables (fast oscillating time-dependent observables in resonance with a free Hamiltonian [8]), which are widely used in spectroscopy [44], (Section 4), [45], (Section 4), making only the projected dynamics experimentally accessible.As we will see in Section 4.1, it also arises as a non-selective measurement of energy (correspondent to free Hamiltonian).
The next proposition gives the perturbative expansion of K(t).
Proposition 1. Asymptotic expansion at fixed t and for K(t) as λ → 0 has the form where K n (t) are defined as where the condition ∑ q j=0 k j = n, k j ⩾ 1 means the sum runs over all compositions of the number n, and The proof follows the proof of Theorem 1 in [40] for a usual time-convolutionless master equation for the case when the initial condition is consistent with the projector.
For clarity, let us give explicit expressions of K n (t) for several first n.
Here, the sum runs over all the permutations σ but with the restriction that inside any P, the terms are time-ordered, i.e., t σ(k) ≥ t σ(k−1) ≥ . . .≥ t σ(j q ) , t σ(j q −1) ≥ t σ(j q −2) ≥ . . .≥ t σ(j q−1 ) , and so on.Many other approaches that are used to calculate the perturbation theory terms of generators in usual master equations can be applied in our superoperator master setup as well.For example, the terms K n (t) can be calculated based on recurrence relations, similarly to [50].Furthermore, the purely algebraic formulae for M k (t) can be obtained in the case when L(t) arises in the iteration picture for dynamics with a time-independent generator in the Schrödinger picture by taking all the integrals with respect to time, similarly to [40].Their generalization to the superoperator master equation case is straightforward, so we will not dwell further on them.

Second Order Superoperator Master Equation
Second order master equations usually play the most important role in the theory of open quantum systems [41], (Sections 3.3-3.4),because the first order one is usually zero in the iteration picture or gives only the terms leading to unitary dynamics [51] and, therefore, has the same form as for closed systems.Therefore, for the superoperator master equation, we concentrate specifically on second order equations.
Taking into account ( 7), ( 10), ( 11) and ( 9), we obtain the following form of the generator in the second order of perturbation theory in λ.
Let us consider the widespread case, where L(t) is Liouvillian in the interaction picture.Namely, we assume L(t) = e −L 0 (t−t 0 ) Le L 0 (t−t 0 ) , where with H 0 = H † 0 as a free Hamiltonian and H I = H † I as an interaction Hamiltonian.Then, taking the integral with respect to s in (15), similarly to [40], (Corollary 1), we have + O(λ 3 ). Here, where P ker[L 0 , • ] is a hyperprojector to the kernel of the map [L 0 , • ], i.e., it is zero on the kernel and the inverse in the usual sense on the orthogonal complement to the kernel.Such a pseudoinverse can be defined by the explicit formula [52], (Equation ( 2)) Thus, within the accuracy of the second order of perturbation theory, the projected propagator Φt t 0 ≡ P(Φ t t 0 ) can be defined as a solution of the Cauchy problem

Superoperator Nakajima-Zwanzig Equation
We are mostly interested in time-convolutionless superoperator master equations, but it is possible to obtain a superoperator analog of the Nakajima-Zwanzig equations as well.Due to condition (5), this analog also has the homogeneous form where G t s is a map from superoperators to superoperators, which is defined as a solution of the Cauchy problem Here, (L(t) • ) is a map from superoperators to superoperators that acts as a left multiplication on superoperator L(t).The derivation is similar to the one from the usual Nakajima-Zwanzig equation [41], (Subsection 9.1.2).By expanding (23) in the Dyson series, we can obtain perturbative expansion of the integral operator in the right-hand side of (22).
Similarly to (14), we can represent the terms of this expansion as time-ordered integrals of (analogs of) Waldenfels cumulants [49].Let us restrict ourselves to the second order Nakajima-Zwanzig superoperator equation.To obtain the second order terms in (22), we can simply use zero order terms in (23), i.e., only G t s ≡ I + O(λ).This leads to the following Cauchy problem for the second order Nakajima-Zwanzig equation for the projected propagator Φt Sometimes [41], (Section 3.3.1),the usual second order time-convolutionless master equations are derived from the usual second order Nakajima-Zwanzig equation by assuming that the kernel of the analog of the integral operator in the right-hand side of ( 24) decays so fast when t − s differs from zero that one can take the density matrix at time t rather than s in this integral operator.However, if we assume a similar approximation here, i.e., that Φs t 0 ≈ Φt t 0 at the right-hand side of ( 24), we nevertheless do not obtain the superoperator master equation with generator (15).The difference is in the fact that Φt t 0 in ( 24) is inside of the action of the hyperprojector P, and this holds after the approximation Φs t 0 ≈ Φt t 0 , but in (20), it is outside of the action of the hyperprojector P. If we take, for example, P = P • P, where P is a usual projection superoperator, then this difference vanishes.However, let us emphasize that, in general, the second order superoperator Nakajima-Zwanzig Equation (24) with "Markovian approximation" Φs t 0 ≈ Φt t 0 is different from the second order time-convolutionless superoperator master equation (20).Therefore, in this aspect, our hyperprojector method and usual projection methods differ.

Stroboscopic Limit Superoperator Master Equation for Hyperprojectors
There are other master equations fully defined by projectors in addition to the weak coupling limit ones.Namely, the master equations that occur in the stroboscopic limit of repeated non-selective measurements [13].Therefore, it is also natural to derive a superoperator analog from the master equations.In a usual setup, it is assumed that non-selective measurement is performed on an ancillary system (probe) at regular short time intervals.Furthermore, the system of interest and the ancillary system interact with each other and unitarily evolve between successive non-selective measurements.
In our setup, we do not need interaction with an ancillary system, but we assume that the measurement procedure can be modeled by a hyperprojector P (we see that it is the case in Section 4.1), and the evolution between measurements can be modeled by the quantum dynamical semigroup Following [13], we consider only the case of a time-independent generator.Then, evolution during small-time N −1 t, N → ∞, but with strong coupling √ Nλ (such a scaling is standard for the stroboscopic limit) is described by the map Φ N −1 t; √ Nλ .If we perform a measurement P, starting before such an evolution and ending after it, then the resulting dynamical map has the form P(Φ N −1 t; √ Nλ ).Now, if we perform such a procedure N times, then the total dynamical map has the form In the limit of an infinite number of such repetitions, it is a superoperator analog of the stroboscopic limit of repeated non-selective measurements [13], and it is described by the following proposition.
Proposition 2. Let Φ t;λ be defined by (26).Then, for sufficiently large N (P(Φ N −1 t; where L eff,str (N) is the effective stroboscopic generator defined by asymptotic expansion as N → ∞.
See Appendix A for the proof.Despite the fact that usually the stroboscopic limit assumes only second order expansion in λ, Proposition 2 allows one to calculate further perturbative correction up to arbitrary order.
Let us remark that there are other analogs of stroboscopic limit setups with a different order of measurements.As we discuss in Section 4.1, a particular hyperprojector P can model the non-selective measurement of increments of "usual" observables between two time moments.For such a hyperprojector, Formula (27) assumes that we measure the differences between each small time-interval of length N −1 t.However, similarly, one can measure the increments between the initial time and successive time moments nN −1 t, where n = 1, . . ., N.Then, instead of (28), for sufficiently large N, we have Similarly, it is possible to consider a setup such that the final time moment for all two-time measurements is fixed, but initial moments are nN −1 t, where n = 0, . . ., N − 1.Then, instead of ( 28), for sufficiently large N, we have One can obtain perturbative expansions for L R,L eff,str (N) similar to the one given by ( 29), but they are more complicated so we do not provide general expressions for them, giving only the first and second terms of the expansions in the next subsection.

Second Order Superoperator Master Equation
Despite the fact that our perturbative expansion is in N − 1 2 , to emphasize the similarity to the weak coupling expansion of superoperator master equations in Section 2.1 we will call the order of perturbation of these equations in terms of order coupling constant λ.In particular, we call expansion (29), which holds the second order terms in λ the second order superoperator master equation.It has the form as N → ∞.
Usually, the stroboscopic limit is formulated in terms of [13] In such a parametrization, the second order generator ( 32) is time-independent, and it takes the form Thus, for Φ str = e tL eff,str (N) we obtain, the following second order superoperator master equation For the setups defined by Formulae ( 30)- (31), similarly to (32), we have and as N → ∞.
The difference between (32), ( 32) and ( 35) is similar to the difference between (15) and the time-convolutionless approximation of (24) discussed at the end of Section 2.2.Furthermore, it vanishes for P = P • P.

Averaging Hyperprojector
One of the first projections used to derive quantum master equations is the Zwanzig projection operator [53], which is nothing else but a particular case of the averaging projector [8] if the energy spectrum is non-degenerate.Therefore, it is natural to consider a hyperprojector analog of such a projector and the superoperator master equations for this specific hyperprojector.
The averaging (with respect to free dynamics e L 0 t ) hyperprojector is defined by the formulae [9], (Definition 2.1) where L 0 is the free Liouvillian defined by Formula (16), where H 0 is a free Hamiltonian, i.e., a Hermitian operator.Let us use the spectral decomposition of H 0 : where ε are (distinct) eigenvalues of H 0 , Π ε are orthogonal projections to the eigenspaces correspondent to these eigenvalues The hyperprojector P can be represented as [9], (Proposition 2.1) Let us summarize some other properties of P as the following proposition.
Proposition 3. The averaging hyperprojector P has the following properties: 1.
Commuting of averaging with free dynamics Commuting of averaging result with free dynamics e L 0 t P(Φ) = P(Φ)e L 0 t .(41)

3.
Identity preserving property where I is the identity superoperator.4.

5.
Unitality preservation P(Φ)I = I for any superoperator Φ such that Φ(I) = I, where I is the identity operator.6.
Coincidence with the hyperprojector to the kernel of the map [L 0 , • ] See Appendix B for the proof.Properties 1 and 2 simplify the transformation from the interaction picture for projected dynamics.Property 3 is important as we assumed (4) to obtain homogeneous superoperator master equations.Property 4 means that P is indeed a hyperprojector.From property 5, it follows that, in particular, for unitary Φ (which is natural to start with), we have unital P(Φ).The last property 6 gives an additional reason to consider such a hyperprojector due to the fact that it occurs anyway in the second order weak coupling superoperator master Equation ( 20) for an arbitrary hyperprojector, but through the definition ( 18)-( 19) of the pseudoinverse [L 0 , • ] (−1) .Therefore, it is natural to consider the particular case when hyperprojectors P and P ker[L 0 , • ] coincide.

Selective and Non-Selective Measurement of Increments
The averaging projector can be interpreted as a quantum operation correspondent to non-selective measurement (see, e.g., [54,55]) of energy (in the sense of the free Hamiltonian).In this subsection, we discuss that it is possible to interpret the averaging hyperprojector as non-selective measurement of increments of energy between two time moments without measuring the final and initial energies themselves.
Usually, multi-time measurements are discussed in terms of linear combinations of multi-time correlation functions (see discussion in [56], (Section 3.5.2)).Therefore, let us interpret the averaging hyperprojector in terms of them.For Markovian dynamics, in particular for the unitary dynamics of a closed system, a multi-time correlation function is defined by the generalized regression formula [41], (Section 3.2.4) for observables Y (k) and X (k) , k = 0, • • • , N, and the time moments Here Φ t s is defined by one-time dynamics ρ(t) = Φ t s ρ(s), t ⩾ s.In general, the multi-time correlations can be non-Markovian and can be defined by process tensors [57,58], but here, we focus mostly on the Markovian case, for which all the multi-time correlations are fully defined by the propagator Φ t s and the initial density matrix.Thus, it is natural to consider the superoperator master equations for the projected propagator only.However, in principle, one can generalize the hyperprojector methods used here as multi-time hyperprojectors to derive multi-time master equations similar to Equations ( 27)-( 28) in [59], where averaging with respect to a classical random variable played the role of such a multi-time hyperprojector.
For simplicity, let us consider the special case of Formula ( 46) for two-time measurements, which takes the form ⟨Y(t)X(s)⟩ = Tr YΦ t s (Xρ(s)), t ⩾ s.
Similarly to Formula (47), let us define Since we interpret hyperprojector P as a result of the non-selective measurement of increments of energy H 0 between times s and t, then it is natural to consider the hyperprojector defining the result of selective measurement as well.Therefore, we define and the correspondent correlation function The hyperprojectors P ω sum up to To verify that correlation functions ⟨Y(t)X(s)⟩ P and ⟨Y(t)X(s)⟩ P ω have operational meaning, let us show that they can be represented as linear combinations of the usual correlation functions (46).Proposition 4. For arbitrary operators X and Y and t ⩾ s, one has

⟨Y(t)X(s)⟩
See Appendix C for the proof.Moreover, analogously to Equations ( 50) and ( 51), we can define multi-time correlations, assuming that the generalized regression Formula (46) changes in a similar way, but only for the two times between which we perform the measurements.Such modified regression formulae fully define the process tensor [57][58][59].Now, let us show that the hyperprojector P ω can be interpreted as a hyperprojector of the (ideal) selective measurement of energy increment (only the increment without measurement of energies themselves).First of all, let us show that it can be used to define a probability of a certain value of the energy increment ω and the posterior state after such a measurement.Proposition 5.If Φ is a completely positive trace non-increasing map, then P ω (Φ) is also a completely positive trace non-increasing map.(P ω is a probabilistic supermap in terms of [14]).In particular, if ρ is a density matrix and Φ is a completely positive trace-preserving map, then 1.
See the proof in Appendix D. Thus, we can interpret Tr(P ω (Φ)ρ) as a probability that some physical quantity equals ω, and the posterior state after such a measurement, if the result of the measurement is ω, is defined by (53).
Finally, the next proposition shows that P ω satisfies several "natural" equalities that could be expected from the energy difference measurement hyperoperator.Proposition 6.For arbitrary superoperator Φ, one has See the proof in Appendix E. Let us discuss the operational meaning of these equalities.It reflects the fact that measuring the energy and energy increment is equivalent to measuring both the initial and final energy.
) means that we started our experiment of measuring energy difference, then just after the start, we measured the energy; after that, the system evolved according to the propagator Φ, and we finished our experiment measuring the energy increment.
) means that we initially measured the energy and then ran the experiment on measuring the energy increment, and so on.
) means that we just measured the energy at the initial and the final moment.Now, let us consider non-selective measurement, i.e., we measure the energy increment but do not look at the result of the measurement.As we obtain state (53) with the probability Tr(P ω (Φ)ρ), then, after such a non-selective measurement, we have the following posterior state ∑ ω Tr(P ω (Φ)ρ) P ω (Φ)ρ Thus, the hyperprojector P can be interpreted as a transformation of the dynamical map Φ between two time moments, which describes the action of non-selective energy increment measurement between these two time moments.
Proposition 5 also leads to some additional properties of hyperprojector P.

1.
If Φ is a completely positive trace-preserving map (quantum channel), then P(Φ) is a completely positive trace-preserving map (quantum channel).(P is a deterministic supermap in terms of [14].) 2.
The usual averaging superoperator [8] can be interpreted as the action of the usual non-selective measurement on a state, and in particular, it transforms density matrices before the measurement to density matrices after the measurement.Similarly, Proposition 7 says that Φ maps channels to channels and unital channels to unital channels.In particular, from [60], (Corollary 7.10), we obtain the following corollary.
Let us remark that S(Φρ) = S(ρ) for Φ = U • U † .Thus, we have also a higher order analog of the H-theorem that can be interpreted as a higher order analog of the second law of thermodynamics from the physical point of view.The usual H-theorem says that entropy of a state increases (or does not change) after transformation by a bistochastic channel.However, here, we see that entropy increment increases (or does not change) if we transform the channel itself via the hyperprojector.It can be interpreted as a second order analog of the second law of thermodynamics if we change channels rather than states.

Correlation Functions of Fast Observables
In work [8], the usual averaging superoperator (55) was used to define the effective Gibbs state for fast observables.Let us show that, similarly, the averaging hyperprojector can be used to define the effective equilibrium correlation functions for these observables.By a fast observable, we mean an observable that explicitly depends on time in the Schrödinger picture in such a way that it becomes time-independent in the interaction picture.Namely, a fast observable X(t) has the following form in the Schrödinger picture [8], (Equation (3)) where X is a constant matrix, where L 0 is defined by (16).Now, let us consider the equilibrium correlation function of two fast observables X(t) and Y(t).By applying the generalized regression Formula ( 46), assuming that the initial state is a Gibbs state ρ β and taking into account the explicit time-dependence of fast observables, we have Similarly to [8], we consider a long timescale such that ωs ≫ 1 for all non-zero Bohr frequencies ω and ∆ωt ≫ 1 for all non-zero differences ∆ω in Bohr frequencies.However, we do not assume that t − s is large with respect to ∆ω −1 or ω −1 .Then, on such a long timescale, we have Thus, on such a long timescale, the correlation functions become approximately stationary, and instead of exact Gibbs state ρ β and dynamical map Φ t s , we can use the effective Gibbs state P (ρ β ) and P(Φ t s ) in regression Formula (47).Such an effective correlation function can be also considered as a multi-time analog of quasi-stationary states [61].
Similarly to [8], we can interpret the entropy increase in (56), in contrast to unitary dynamics, as the information being lost due to our restricted experimental capabilities, which assume that we can measure only "slow" averages of fast observables.
Let us also remark that in (57), X can be time-dependent but "slow" on the long timescale.We can think of X as a "slow envelope".

Weak Coupling Superoperator Master Equation for Averaging Projector
Now, let us consider a second order weak coupling superoperator master equation in the special case of the averaging hyperprojector.Namely, let us consider unitary dynamics Ψ t;λ defined by equation The properties of the averaging hyperprojector simplify the transformation into and from the interaction picture.
Proposition 8.For projected propagator P(Ψ t;λ ), the following Cauchy problem for the superoperator master equation holds where as λ → 0 and fixed t ⩾ t 0 .
See the proof in Appendix F. Let us expand H I as a sum of eigenoperators of the superoperator [H 0 , • ] in the same way as it is usually performed for the Markov master equation derivation [41], (Subsection 3.3.1),and use it to obtain more explicit representation of (60), as we have carried out for the effective (equilibrium) Hamiltonian in [8].We have where the sum runs over all the eigenvalues of [H 0 , • ], which are sometimes called the Bohr frequencies [23], p. 122 and The possibility to choose D ω such that D −ω = D † ω follows from the fact that H 0 is Hermitian.Then, Equation (60) takes the form (see Appendix G for detailed calculation) The first terms of this generator are nothing else but the commutator with the Hamiltonian in the rotation wave approximation H 0 + λD 0 = H RWA .The Hamiltonians obtained by perturbative corrections are usually called effective Hamiltonians [31,32].In our case, we obtain not only Hamiltonian terms, so we call it an effective generator.
Let us remark that despite the fact that any Hermiticity and trace persevering generator has a GKSL-like form [62], the specific form ( 63) is such that Lindblad operators defined by ( 62) are the so-called weak coupling limit type (WCLT) form [63]. Therefore, our derivation shows that the averaging hyperprojector is the essential component for obtaining the weak coupling limit type form of generator.Actually, even in the derivation [22] of usual weak coupling master equations, we initially obtained a non-GKSL Redfield-like equation.However, the Redfield equation has a known physical issue in that it can predict non-positive probabilities in some cases [64].Furthermore, only after using the averaging (but without regarding it as application of averaging hyperprojector explicitly), the WCLT GKSL generator occurs, which has no such issue.This suggests that hyperprojector methods can be essential for obtaining physically consistent equations.
Equation ( 63) has the GKSL-like form but not the GKSL form exactly, if Hence, despite the fact that P(Ψ t;λ ) is completely positive and, due to corollary 1, is entropy non-decreasing from time t 0 to t, P(Ψ t;λ ) is not completely positive divisible and it is not necessarily entropy non-decreasing for any fixed time t.A more detailed description of the local entropy increase is given by the following proposition.
for all density matrices ρ and small enough λ.

2.
If sin ω(t for all density matrices ρ and small enough λ. It follows from the fact that in the first case of this proposition generator, (63) has the GKSL form and P(Ψ t;λ ) is completely positive divisible and L eff (t, t 0 ; λ)I = 0 so it is unital, then by [60], (Corollary 7.10), we obtain (64).Similarly, in the second case of this proposition generator (63) has the opposite GKSL form and (P(Ψ t;λ )) −1 is completely positive divisible and unital, so we have (65).This proposition describes some extreme situations, when the entropy increment at the time moment t has a definite sign for all the initial states.In general, it cannot be the case that sin ω(t − t 0 ) ⩾ 0 or sin ω(t − t 0 ) ⩽ 0 for all positive Bohr frequencies ω and the entropy increment at the time moment t can have different signs depending on state.However, the total entropy increment is non-negative by Corollary 1.
To illustrate (63), let us consider the simplest case of a two-level system then, (63) takes the form

Bogolubov-Van Hove Limit
The mathematically strict derivation of Markovian master equations is based on the Bogolubov-van Hove scaling [22,23] and from the physical point of view the Bogolubovvan Hove limit identifies the timescale consistent with the projector.So let us discuss the Bogolubov-van Hove scaling analog in our case as well.We obtain the following proposition.
Proposition 10.Let L eff,BvH (λ) be defined by Using expansion ( 61)- (62), we obtain a more explicit form of the generator in the Bogolubov-van Hove limit: Thus, after the Bogolubov-van Hove scaling in the limit λ → 0, the generator becomes a generator of fully unitary dynamics without any dissipator-like terms, which is not obvious from the abstract form (68). Let us also remark that there is a connection of Formula (69) with the algebraic perturbation theory [30,32,52].The algebraic perturbation theory for given L 0 and L I finds perturbatively such U (λ) and L eff,sec (λ) that where U (λ) = U(λ) • U † (λ) and U(λ) is a unitary transformation and L eff,sec contains only the "secular" terms, which can be formalized [52] as the condition L eff,sec = P(L eff,sec ).
Comparing with Formula (13) from [52] for the effective Hamiltonian, we obtain However, let us remark that after the Bogolubov-van Hove scaling, the unitary transformation U(λ), which can also be found in [52], (Section IV.B) seems to play no role.This is because L eff,BvH (λ) in ( 69) has higher order corrections O(λ), so with the same accuracy, U(λ) = I + O(λ) and the unitary transformation can be neglected.It is the result of the Bogolubov-van Hove scaling, but U(λ) plays the role for the higher order expansions of L eff,BvH (λ).
Therefore, despite the fact that the RWA (or secular approximation) and perturbative corrections to it are usually not formulated in terms of hyperprojectors and the Bogolubovvan Hove limit, these widespread physical assumptions can be compactly formulated via the hyperprojector method.For deriving systematic perturbative corrections, such a compact formulation can be important.
Let us also remark that from the point of view of interpretation discussed in Section 4.2, another analog of the Bogolubov-van Hove scaling can be considered, where both initial and final time are taken at the long timescale.This means that both t and t 0 are subjected to scaling.Although λ −2 t and λ −2 t 0 are large as λ → 0, it is possible to consider such sufficiently close times that λ −2 (t − t 0 ) has a finite limit as λ → 0. Furthermore, no simplification occurs in (63) in such a case.However, the case when λ −2 (t − t 0 ) is also large reduces to Proposition 10.
Let us also remark that in this work, we have focused on the finite-dimensional case.However, in the infinite-dimensional case, the sum in the dissipator-like term in (63) becomes the integral, which, similarly to [65], can give non-trivial terms in the limit due to the Fermi golden rule; therefore, it can be a non-trivial dissipator on the timescale identified by the Bogolubov-van Hove scaling.Owing to the connection with the algebraic perturbation theory mentioned in Section 5.1, it is also natural to ask if effective generators of unitary dynamics for the infinite-dimensional case always generate the unitary dynamics.This discussion allows us to hypothesize that possibly it is not the case.Some works interpreting resonances [66][67][68] in the infinite-dimensional case as spontaneous time-symmetry breaking [69] are also indirect evidence in favor of this hypothesis.

Stroboscopic Limit Superoperator Master Equation for Averaging Hyperprojector
Now, let us consider the stroboscopic limit superoperator master equation in the special case of the averaging hyperprojector.In an abstract form (32), it is not overly simplified in such a special case; therefore, we turn to a more explicit form.Using expansion (61) and (62), generator (32) takes the form Let us remark that we hold the same notation H 0 , H and for expansion (61) and (62) to emphasize the similarity between the calculations.However, we do not have to understand H 0 as some Hamiltonian, but we could consider it as any other quantity.For the quantum thermodynamics applications, it can be volume or voltage or some other quantity, the increment of which we repeatedly measure.
Moreover, we can think about it not as a repeated measurement but as a constraint on dynamics, which indirectly describes the influence of the environment on an actually open system.For example, in non-quantum physics, we can just assume the ideal mechanical constraints or fixed voltage in a socket without describing a microscopic model of a surface or a rope that sets the ideal mechanical constraints and without a detailed model of a power plant that supplies the fixed voltage to the socket.In a certain sense, we also simply "project" unitary dynamics to constrained dynamics.Moreover, in our quantum case, the consistency with P for evolution between arbitrary close time moments can be considered as a constraints-based description of the environment's influence on the system.Generator (71) has the GKSL form.If we apply our projected dynamics to an initial density matrix, we can say that we have derived a quantum master equation without an explicit description of the reservoir, but only in terms of constraints on dynamics, which can be useful in the case where we do not have a detailed model of the reservoir and have access only to the system.To some extent, our approach is close to the discussion of the decoherence without a reservoir in [65].
As ( 71) is a generator of a unital semigroup [60], (Corollary 7.10), then we have the following proposition.
Proposition 11.For large enough N and arbitrary density matrix ρ, we have d dt S e L eff,str (N)t ρ ⩾ 0.

Conclusions
We have introduced superoperator master equations for the cases of the weak coupling limit and the stroboscopic limit in Sections 2 and 3, respectively.We think that these equations can be interesting in a broad context, but in this work, we have focused on the particular case of the averaging hyperprojector.In Section 4.1, we have shown that it has physical meaning as a non-selective measurement of energy increments between two times without measuring the initial and final energies themselves.In Section 4.2, we have shown that it arises naturally if we measure the correlation functions of fast observables.It is important because in widely used spectroscopy setups [44], (Section 4), [45], (Section 4), we have no direct access to the density matrix but only to correlation functions.In particular, the multi-time correlations have now attracted great attention in open quantum systems theory [57,58,70].Under the conditions of Corollary 1, we have also shown that such a hyperprojector leads to an increase in entropy increments, which can be regarded as the second order analog of the second law of thermodynamics.In Sections 5 and 6, we apply the general approaches developed in Sections 2 and 3 to the specific averaging hyperprojector.In Section 5, we have obtained the generator that has the WCLT GKSL-like form, which suggests that the averaging hyperprojector is essential for obtaining such a form.In Section 5.1, we have shown that in the Bogolubov-van Hove limit, this generator coincides with the second order effective generator from the algebraic perturbation theory, which suggests that RWA (or secular approximation) and its perturbative corrections can be absorbed by hyperprojector methods.In Section 6, we have derived the superoperator master equation arising due to the repeated measurement of increments of the observable that defines the averaging hyperprojector.It has the GKSL form and leads to dynamics with monotone entropy.
As a possible direction for further development, we should mention the consideration of analogs of other limits leading to Markovian master equations for open quantum systems, e.g., singular coupling limit [41], (Section 3.3.3),a singular and weak coupling limit combination leading to a unified master equation [71], low density limit [72][73][74] and so on.It can be interesting to find hyperprojectors that lead to the generators satisfying quantum detailed balance [75] rather than simply the unital, as in our work.