arXiv:1701.02936v1 [quant-ph] 11 Jan 2017 codnl ls system close a Accordingly set rlHamiltonians trol htcnb mlmne nti a sfre ythe by formed group is way Lie this the in of implemented unitaries be can that oepoemr drcin. o finite-dimensional a system For the system quantum enabling “directions.” closed more ones, explore original to linearly the are of which independent generators with effective associated produce evolutions can two one by Hamiltonians induced noncommuuting transformations concatenating properly pc pne yalpsil iercmiain fthe of combinations linear of possible elements all by spanned space aia i algebra Lie namical ayeouini uteuvln oteoegnrtdby uni- generated one resulting of the combination the to linear fields, equivalent a just control is external evolution by tary off and on t ftersligdnmc.Frisac,frtocom- two for instance, For Hamiltonians dynamics. muting complex- resulting the the generators influences of the fields ity of control the noncommutativity with the associated desired Here a time-dependent [2]. realize classical quantum to fields shaping by of system it quantum goal on a noncommu- The transformation steer the to role. is [1] important control ap- an theory Lie-theoretical plays control the tativity quantum in to particular In proach of and/or outcomes operations experiment. the which an influences in performed order are measurements The mechanics. tum ocmuaiiyi n ftekyfaue fquan- of features key the of one is Noncommutativity H hita Arenz, Christian ssi ob ul otolbeif controllable fully be to said is H hd ute ih nterl fmaueetbcato on backaction measurement of role t the conclude syste on to quantum light us open further allows nonaccessible sheds This a measurements sizes. selective system which generic Lindbladians for of pair system universal a the auxilia provide measuring we frequently an Moreover, constructio through with recovered simple space be can larger a dynamics a presenting on measu commutative by frequent Hamiltonians, through dynamics space open-system Hilbert to original a th the Hamiltonians in onto noncommutative dynamics back original the Hamiltonian m The noncommutative thus making be and of purification.” can idea space, This other, Hilbert simple. each extended and with an an noncommutative in opposite: general them the in of think are also can which We other. each with commute h otolblt vrtesse eoetepoeto i cont projection full the even before subspace system this the In over controllability complex. quantum the more a be of can dynamics the dynamics project can measurements frequent narcn ok[.K Burgarth K. [D. work recent a In n hi trtdcmuaos[,3 4]. 3, [1, commutators iterated their and .INTRODUCTION I. 3 H ET coaNraeSproeadIttt Nanoscienze-C Istituto and Superiore Normale Scuola NEST, := L H ( ,2 1, H 1 { H Q H and soitdt h e fcon- of set the to associated ) 4 1 h e feetv evolutions effective of set the , ailBurgarth, Daniel eateto hsc,Wsd nvriy oy 169-8555, Tokyo University, Waseda Physics, of Department 1 2 and Q rc aoaoy rneo nvriy rneo J0854 NJ Princeton University, Princeton Laboratory, Frick H , e L H 1 ( hrceie yacontrol a by characterized 2 H nttt fMteais hsc,adCmue Science, Computer and Physics, Mathematics, of Institute 2 . . . , H ) hc a eswitched be can which , brswt nvriy brswt Y32Z UK 2BZ, SY23 Aberystwyth University, Aberystwyth 2 eeae ytedy- the by generated ntecnrr,by contrary, the On . } .. h elvector real the i.e., , e L ( H idlda purification Lindbladian 1 ) Dtd etme 6 2018) 16, September (Dated: tal. et nldsall includes itroGiovannetti, Vittorio a.Commun. Nat. , n h orsodn i group Lie corresponding the and n omttr feeet of elements combinations of linear commutators scenario, and system closed the with ance akei ysuyn h i algebra to Lie tempted the be studying would by One it tackle problem. unsolved difficult a is eeaos(.. hycno ecs nteGL form), GKLS differently, the it in said cast dynamical be or allowed cannot they longer (e.g., no generators are which super-operators hc r nhscl utemr,ee o h elements the for of even Furthermore, unphysical. are which lsia ed hc ciaete opromalphysi- all perform to them using activate of which ability fields the classical for i.e., facts no- controllability, physical these of the of tion distinguishes view one In systems simply quantum fields. by open them control implement the to with possible playing be would it if clear Lnbain ntefloig.Silas nti case this collection in given control- also by a activated ling Still be can dynamics 8] which following). determining the [7, Gorini- in form of (Lindbladians dy- super-operators (GKLS) of the notion Kossakowski-Lindblad-Sudarshan equipped the i.e., provided are admit generators, is which and namical exception structure processes semigroup partial Markov a A of with subset 6]. the generators [5, by dynamical lost of typically notion is the a and possess structure not do group more operations allowed becomes the controllability, since of complicated notion well associated as the operations, as (realizable) reachable the of terization accessibility. as to referred rtragbaof algebra erator ftednmclLealgebra Lie dynamical the if osbeuiaytasomtoson transformations unitary possible hni oe ooe unu ytm,techarac- the systems, quantum open to comes it When e L 5 eypo ic h oto Hamiltonians control the since poor very s ( 13(04]i a hw htasre of series a that shown was it (2014)] 5173 , L 3 ) ecie n“cesbe pnquantum open “accessible” an describes ioih Nakazato, Hiromichi uiir ytmi o-eetv way. non-selective a in system auxiliary ercvrdb rjcigtesystem the projecting by recovered re hc r hsclyalwd ti ngnrlnot general in is it allowed, physically are which eeal eoe cesbe This accessible. becomes generally m eet.Hr egnrlz hsidea this generalize we Here rements. ysse.W hwta h original the that show We system. ry a hog eiso rqetnon- frequent of series a through hat h oto fqatmsystems. quantum of control the omttv scle “Hamiltonian called is commutative s aitnaso unu system, quantum a of Hamiltonians y omake to n olblt a eahee,although achieved, be can rollability ytmot usaei hc the which in subspace a onto system xeddsaebcmstrivial becomes space extended e d omttv yembedding by commutative ade R -62 ia Italy Pisa, I-56126 NR, Q hsls rpryof property last this , Lindbladians e L L ,US 4, ( := L Japan ) {L ilicuetransformations include will 4 L n auaYuasa Kazuya and 1 ( , L H L swl as well as , e L ili eea produce general in will 2 pn h hl op- whole the spans ) ( . . . , L L ) Q ( oee,a vari- at However, . L } requivalently, or , eeae by generated ) fLindbladians of H en also being L n the and 4 L 2 cally allowed quantum transformations, from the weaker commutative on an extended space by means of an auxil- notion of accessibility, which in this case corresponds to iary system, which, through frequent non-selective mea- have L(L) equal to whole Lie algebra generated by arbi- surements, yields the original noncommutative dynam- trary Lindbladians. In the following, we will refer to this ics. Moreover, we present a universal pair of Lindbladi- as the ‘GKLS algebra’ noting however that it contains ans that generate the full GKLS dynamical Lie algebra many elements which are not of GKLS form. Differently for generic system sizes, the analysis providing us with from the closed quantum system case, it is indeed possi- a short and elementary proof of the generic accessibility ble that a control set L is accessible but not controllable. [7]. Applying hence the Lindbladian purification proce- Still studying the accessibility of a collection of Lind- dure to such a universal set we then show that almost all bladians is a well posed mathematical problem, which open systems become accessible, even though their gen- can also shed light on the controllability issue, with ac- erators are commutative with each other, by performing cessibility being a necessary condition for controllabil- frequent non-selective measurements on a part of the sys- ity. Furthermore accessibility implies that the reachable tem. set has non-zero volume and therefore has physical rele- This article is organized as follows. Along the lines vance: the short time dynamics explores a high dimen- of [12] we begin in Sec. II by reviewing the definition sional space and is therefore of high complexity. of Hamiltonian purification and presenting an explicit It turns out that almost all control sets L are acces- construction for purifying an arbitrary number of Lind- sibile [7]. Analogously to the case of closed systems, bladians and Hamiltonians. In Sec. III we consider the the key ingredient of this result can be identified with accessibility of controlled master equations. Concluding the noncommutativity of the elements of L. But what remarks are given in Sec. IV, and some details on the about models where L includes only mutually commutat- derivation of the projected dynamics and the proof of ing Lindbladians? Is there a way to expand their algebra accessibility are provided in the Appendices. L(L) to cover the full GKLS algebra? For close quantum systems it has been observed that one can substantially change the dimension of the dynamical Lie algebra L(H) II. LINDBLADIAN PURIFICATION AND NON-SELECTIVE ZENO MEASUREMENTS through frequently observing a part of the system [9], or by tampering it with a strong dissipative process that ex- hibits decoherence-free subspaces [10] (the gain being ex- To begin with we first review the definition of Hamil- ponential in some cases). As a matter of fact, on the ba- tonian purification [12]. Suppose that we have n control sis of the quantum Zeno effect [11], starting with a set of Hamiltonians, which are switched on and off to steer a d- H dimensional quantum system Q. Let H = H ,...,H commuting control Hamiltonians , noncommutativity { 1 n} can be enforced through frequently projecting out part of be the set of the control Hamiltonians acting on the Hilbert space of Q, and H˜ = H˜ ,..., H˜ be a cor- the system onto a subspace where accessibility and hence Hd { 1 n} full controllability is achieved. Also it has been observed responding set of Hamiltonians acting on an extended that the projection trick can be reversed: specifically, Hilbert space dE of dimension dE (> d), which includes H H˜ H starting from a set of noncommuutative Hamiltonians H d as a proper subspace. We call a purifying set of H one can construct a new set H˜ formed by commutative if all the elements of H˜ commute with each other, elements on an extended Hilbert space which under pro- H˜ H˜ = H˜ H˜ , i, j 1,...,n , (1) jection reduces to the original one. This mechanism was i j j i ∀ ∈{ } studied in great detail in [12], where, borrowing from the and they are related to those from H through notion of purification of mixed quantum states [13], the term Hamiltonian purification was introduced. H = P H˜ P, j 1,...,n , (2) j j ∀ ∈{ } One may then ask whether a similar procedures can with P being the projection onto . For a generic set H be applied to the algebra of a set L of Lindbladians, d consisting of n linearly independentH Hamiltonians it can namely, if it is possible to enlarge L(L) by means of some be shown [12] that there always exists an H˜ where the projection mechanisms and, on the contrary, if Lind- (min) bladian purification is always achievable. In this arti- minimal dimension dE of the extended Hilbert space is cle we address these issues showing that indeed any set bounded above by d(min) nd. For instance for the case E ≤ L of Lindbladians can be “purified,” i.e., can be made with n = 2 Hamiltonians H1 and H2, Proposition 1 of commutative with each other, by embedding them in Ref. [12] states that a purifying set can be constructed on a larger space (note that the term “pure” was already = with an auxiliary single qubit Hilbert HdE Hd ⊗ HdA used in [14] for Markovian generators in a slightly dif- space dA , the purifications and the projector being ferent way). To this end, we need to employ a different H H˜ = H 11 + H σ , scheme from those for the Hamiltonian purification in- 1 1 ⊗ 2 2 ⊗ x troduced in [12], since the naive application of the lat- H˜2 = H2 112 + H1 σx, (3) ter trivially violates some structural properties of GKLS ⊗ ⊗ generators (more details in the following). Our construc- 112 + σz P = 11 , (4) tion allows one to make Lindbladians and Hamiltonians d ⊗ 2 3 with σ , σ , and 11 the Pauli and the identity opera- identifying the Hamiltonians H˜ and the Lindblad op- x z 2 { j } tors of the auxiliary qubit, respectively. The mapping erators L˜ of the purifying element ˜ = ˜ + ˜ { j,α} Lj Kj Dj H˜j Hj can finally be realized through the quantum as Zeno→ effect [11, 15] by frequently monitoring the extended system via a von Neumann measurement which projects ˜ Hj = nHj j j , the system onto , i.e., ⊗ | ih | j 1,...,n , (10) d  ∈{ } H L˜j,α = √nLj,α j j , ˜ ⊗ | ih | lim (Pe−iHj t/N P )N = e−iHj tP. (5)  N→∞ with j n being an orthonormal basis for . {| i}j=1 Hn The question arises if an analogous construction can be Obviously through such a construction the operators extended to the case of Lindbladians. Specifically con- H˜ and L˜ commute with each other for differ- sider a set L = ,..., of n GKLS generators op- { j} { j,α} {L1 Ln} ent j trivially ensuring the requirement (9). Regard- erating on a target system Q, ing the analog of (5) we focus on non-selective projective measurement [17, 18] operating on the auxiliary system, = + , j 1,...,n , (6) Lj Kj Dj ∈{ } i.e., the completely positive and trace preserving (CPTP) mapping of the form with j and j being the Hamiltonian and dissipator contributions,K D i.e., the super-operators ( )= P ( )P , (11) P ··· X k ··· k ( )= i[H , ], (7) k Kj ··· − j ··· ( )= [2L ( )L† given in terms of a complete set of orthonormal projec- Dj ··· X j,α ··· j,α α † † tion operators Pk corresponding to measurement out- L Lj,α( ) ( )L Lj,α], (8) { } − j,α ··· − ··· j,α comes and satisfying PkPk′ = δkk′ Pk and k Pk = 11. Notice that if we perform (N + 1) of such non-selectiveP L being the Lindblad operators acting on the Hilbert j,α measurements at regular time intervals t/N during the space d of Q. We ask whether if it is possible to asso- evolution driven by a Lindbladian , the system will ciate withH L a purifying set L˜ = ˜ ,..., ˜ formed by 1 n evolve according to the CPTP transformationL GKLS generators possibly acting{ onL an extendedL } system, which are mutually commuting, i.e., Φ(L) := ( eLt/N )N t,N P◦ ◦P 2 N ˜ ˜ ˜ ˜ t t i j = j i, i, j 1,...,n , (9) = id + ( ) + O 2 , (12) L ◦ L L ◦ L ∀ ∈{ } h P◦L◦P N  N i ◦P from which one can recover the original elements via a which in the limit of N converges to projective mapping that should mimics (5) (in the above → ∞ expressions we used the symbol “ ” to indicate the com- (L) (L) (P◦L◦P)t Φt,∞ = lim Φ = e , (13) position of super-operators). ◦ N→∞ t,N ◦P L˜ A natural guess for identifying and the projective where id is the identity map and where we used the idem- mapping would be to simply transporting the purification potent property = = 2 of (11). Equation (13) schemes of Ref. [12] at super-operator level, or equiv- can also be derivedP followingP◦P a pertubativeP approach with alently, to represent the s as operators in Liouville Lj a strong amplitude-damping channel inducing the pro- space [16] and then simply applying to them the Hamilto- jection [19–21]. In our construction Eq. (13) is the nian purication scheme. This simple trick however does formal counterpartP of the Zeno limit (5): it shows that not work because, for instance, mapping as (3) will take alternating the dynamics induced by a GKLS generator positive operators into non-positive one, hence spoiling with induces on the system an evolution which can one fundamental property of GKLS generators. Another beL effectivelyP described in terms of an effective dynami- problem comes from the fact that for Markovian open cal generator described by the projected super-operator systems described by Lindbladians, the quantum Zeno . It should be stressed that this last is not in effect, which as we have seen is responsible for the im- GKLSP◦L◦P form, i.e., it is not a Lindbladian. Indeed it acts as plementation of the mapping H˜ H , does not take j → j a proper Lindbladian only within the subspace specified place: a Markovian system can leak from one subspace by the super-projector , but the map e(P◦L◦P)t itself is specified by the projection operator belonging to a mea- not CPTP (an explicit exampleP of this fact is provided in surement outcome even in the limit of infinitely frequent Appendix A). Still we are going to identify (13) with the projective measurements. In spite of these issues how- mechanism that yields the original Lindbladians j L ever a Lindbladian purification scheme can be obtained expressed in the form (6)–(8) from their purifiedL coun-∈ with the following simple construction: terparts ˜ of with (10). For this purpose we assume the Lj projectors Pk in (11) to be of the form A purifying set L˜ can always be constructed by introduc- ing an auxiliary Hilbert space of dimension n and P = 11 φ φ , (14) Hn k ⊗ | kih k| 4

dA L L where φk k=1 is an orthonormal basis for the auxiliary of ( ) and simple examples, we refer to [7]. Here we fo- Hilbert{| spacei} which is chosen to be mutually unbiased cus instead on studying how the purification mechanism [22] against the orthonormal basis j dA used for the can influence the dimension of L(L). In particular we {| i}j=1 purification (10). Then as shown in Appendix B one shall see how a set of commutative Lindbladians can be can verify that under the transformation (13) a generic turned into a new set of noncommuutative Lindbladians density operator ρQ(0) for the original system, obtained which grant accessibility to the full GKLS algebra via the by taking the trace over the auxiliary Hilbert space dA , projection through frequent measurements. evolves according to H To show this we start by showing that it is possible to identify a set L formed by just a pair of Lindbladians Lj t ρQ(t)= e ρQ(0), (15) whose algebra L(L) spans the full GKLS algebra. We therefore first prove that the pair recovering hence the original dynamics generated by the unpurified Lindbladian . 0 = i adH0 + , (20) Lj L − D|1ih2| = i ad (21) K − |1ih1| III. ACCESSIBILITY with

d−1 We now turn our attention to the question on how fre- H = j j +1 + h.c., (22) quent non-selective measurements can enrich the algebra 0 X | ih | L(L) of a Markovian open quantum system described by j=1 L a collection of controlled generators. Specifically we where shall focus on systems driven by master equations of the form i ad ( )= i[H, ], (23) − H ··· − ··· ∂ ( )=2L( )L† L†L( ) ( )L†L, (24) ρ(t)= (t)ρ(t), (16) DL ··· ··· − ··· − ··· ∂t L does the job, namely, every possible Lindbladian can be where the super-operator (t)= (t)+ is provided by generated by linear combinations and iterated commu- L K D a constant dissipative part represented by Lindblad op- tators of (20) and (21). We only sketch the main steps erators Lα, and by a time-dependent Hamiltonian term here, whereas the details can be found in Appendix C. (t)( )= i[H(t), ] with In the following we also use the notations K ··· − ··· m † † † A,B( )=2B( )A A B( ) ( )A B, (25) H(t)= H0 + uk(t)Hk, (17) D ··· ··· − ··· − ··· X Ad ( )= U( )U †. (26) k=1 U ··· ··· m u (t) being classical control fields that can be op- We first note that terms of the form i ad|jihj| com- { k }k=1 − erated to switch on and off m control Hamiltonians mute with the dissipative part |1ih2| and according to m D Hk k=1. This corresponds to having a control set [23] we can generate every element in i adu(d) with L{ :=} , ,..., consisting of a drift (unmodu- u(d) the Lie algebra of d d hermitian− matrices. Us- {L0 K1 Km} × lated) term ing AdU(d) = exp( i adu(d)) with U(d) being the unitary group, we can get − = + , (18) L0 K0 D AdU |1ih2| AdU † = U|1ih2|U † (27) that includes both the dissipative part and the Hamil- ◦ D ◦ D D tonian contribution 0( )= i[H0, ], and of the set for any U U(d). Now we consider unitaries U that act K ··· − ··· ∈ (j) of Hamiltonian control generators as U j = c k for j = 1, 2 and = 1, 2, 3, 4 . | i Pk∈I k | i I { } We numerically verified that, from U|1ih2|U † , thus cre- k( )= i[Hk, ], k 1,...,m . (19) ated together with i ad for allD Hamiltonians H = K ··· − ··· ∈{ } − H i,j∈I hij i j having support on , all the operators As already mentioned in the introduction, for a closed Pof the form| ih | I quantum system, i.e., without the dissipative part , the algebra L(L) associated with L (i.e., the set of all realD lin- |iihj|,|kihl| + |kihl|,|iihj|, ear combinations and iterated commutators of these ele-  D D i,j,k,l (28)  ∈ I ments, drift term included) will fully characterize the set i( |iihj|,|kihl| |kihl|,|iihj|), of unitary operations that can be implemented through  D − D m shaping the control functions uk(t) k=1. For an open can be generated. Doing the same for different quartets quantum system described by the{ master} equation (16), = i,j,k,l we are able to provide linearly indepen- instead, L(L) only characterizes the accessibility of the dentI operators{ } (28) for all i,j,k,l 1,...,d . Since any system. For a detailed analysis of the general structure Lindbladian can be written in the∈{ Kossakowski} form as a 5 linear combination of those operators, it means that ev- that the noncommutativity distinguishes between quan- ery Lindbladian can be generated through iterated com- tum and classical mechanics, it is tempting to say that mutators and linear combinations of the pair of gener- every quantum system can be made classical by purifying ators 0, in (20) and (21). Given that this specific it to a larger space. pair of{L LindbladiansK} is accessible, it then follows from the standard argument (see, e.g., [9]) that almost all pairs are. This was shown previously in a more abstract way ACKNOWLEDGMENTS by Kurniawan [7]. Now that we have found a pair L = 0, that describes an accessible quantum system in{L arbitraryK} di- We like to thank John Gough for fruitful discussions. mensions, we can make them commutative using a two- This work was supported by the Top Global Univer- sity Project from the Ministry of Education, Culture, dimensional (dA = 2) auxiliary Hilbert space, i.e., we can purify them to Sports, Science and Technology (MEXT), Japan. DB ac- knowledges support from EPSRC grant EP/M01634X/1. KY was supported by the Grant-in-Aid for Scientific Re- ˜ = 2i ad + √2 , (29) L0 − H0⊗|2ih2| D|1ih2|⊗|2ih2| search (C) (No. 26400406) from the Japan Society for the ˜ = 2i ad . (30) Promotion of Science (JSPS) and by the Waseda Uni- K − |1ih1|⊗|1ih1| versity Grant for Special Research Projects (No. 2016K- Obviously on the extended Hilbert space the Lie alge- 215). HN was supported by the Waseda University Grant bra associated with the set L˜ = ˜ , ˜ is just two- for Special Research Project (No. 2016B-173). {L0 K} dimensional, dim L(L˜) = 2, and the system is not ac- cessible. If we perform frequent non-selective projective measurements on the auxiliary system described by the Appendix A: Projected Lindbladians superprojector (11) with P± = 11 , where are defined at the end of Sec. II, the⊗|±ih±| original dynamics|±i is re- covered as (15) and the system becomes accessible. The As an example of the fact that the projected counter- part of a Lindbladian does not generate proper existence of such a specific setup allows us to conclude [9] P◦L◦P L that almost all open quantum systems become accessible quantum dynamics, consider for instance the case where by Zeno measurements. describes a qubit amplitude damping with fixed point L0 0 (this is characterized by a null Hamiltonian term | ih | H = 0 and a unique Lindblad operator Lα = 0 1 ) and where the transformation is the dephasing| mapih | [13] IV. CONCLUSIONS associated with the canonicalP qubit base, i.e.,

We have generalized the work [12] on Hamiltonian ( )= 0 0 ( ) 0 0 + 1 1 ( ) 1 1 . (A1) purification by establishing a new and simple purifica- P ··· | ih | ··· | ih | | ih | ··· | ih | tion scheme for Lindbladians, which is also applicable to Accordingly for an arbitrary density matrix ρ we have Hamiltonians. Given n Lindbladians, they can be made commutative by adding an n-dimensional auxiliary sys- lim (e(P◦L◦P)t )ρ = 0 0 , (A2) tem to extend the Lindblad operators with hermitian pro- t→∞ ◦P | ih | jectors that form an orthonormal basis for the auxiliary while on the contrary space. Through the projection by Zeno measurements for semigroup dynamics the original possibly noncommuta- lim e(P◦L◦P)tρ = 0 0 + (id )ρ tive dynamics can be recovered by frequently measuring t→∞ | ih | −P the auxiliary system in a non-selective way. Moreover, we = 0 0 + 0 0 ρ 1 1 + 1 1 ρ 0 0 , | ih | | ih | | ih | | ih | | ih | have proven that the pair of Lindbladians (20) and (21) (A3) describes an accessible open quantum system for generic system sizes, which tells us that generally a nonacces- which in general is not a valid state. sible open quantum system is turned into an accessible one by frequent non-selective measurements. The model has also potential applications in simulating an arbitrary Markovian open system dynamics [24, 25] by steering it Appendix B: Derivation of the Projected Dynamics through control fields. Clearly, the presented purification scheme also works We start by noticing that given a generic non-selective for observables and density operators, although, except transformation as in (11) and the unitary generator for the partial trace, an operational way that allows us to with HamiltonianP H, the following identity holds K recover the original observables and states is not known to us. Since the noncommutativity is a unique feature of ( )( )= i [H(k), P ( )P ], (B1) P◦K◦P ··· − X k ··· k quantum mechanics, and in fact it was argued in [26, 27] k 6

(k) where H = PkHPk. Similarly, given a dissipator (B9) we obtain characterized by Lindblad operators L we have D α ∂ ρQ(t) = TrA[( ˜j )ρ(t)] ( )( ) ∂t P◦ L ◦P † P◦D◦P ··· ′ ′ ′ ′ (kk ) (kk )† (kk )† (kk ) = j (TrA[ ρ(t)]) + 2 Lj,α TrA[ ρ(t)]Lj,α, = 2L ( )L L L ( )Pk′ L P X T X h α ··· α − α α ··· α α,k,k′ (kk′ )† (kk′ ) (B12) Pk′ ( )L L , − ··· α α i which finally yields the thesis (B2) ∂ ρ (t)= ρ (t) = ρ (t)= eLj tρ (0), (B13) (kk′) ∂t Q j Q Q Q with Lα = PkLαPk′ . Assume next H and Lα as those L ⇒ ˜ associated with the Lindbladian j with (10), and Pk as by noticing that in (14). Since φ is mutuallyL unbiased with respect {| ki} to j the following identity holds, TrA[ ρ(t)] = ρQ(t), TrA[ ρ(t)]=0. (B14) {| i} P T

−iϕjk j φk = e / dA, (B3) h | i p Appendix C: An Accessible Pair of Lindbladians with ϕjk generic phases, and hence Here we show that the pair of Lindbladians , {L0 K} ˜ (k) ˜ in (20) and (21) generates an accessible system. We use Hj = PkHj Pk = Hj φk φk = Hj Pk, (B4) ′ ⊗ | ih | the notations (23)–(26). First of all, we show that = (kk ) i(ϕjk −ϕ ′ ) K L˜ = P L˜ P ′ = e jk L φ φ ′ / d , i ad commutes with the dissipative part of j,α k j,α k j,α k k p A |jihj| |1ih2| ⊗ | ih | − in (20). Using an identity [28] D (B5) L0 (kk′)† (kk′) † † 1 1 L˜ L˜ = L Lj,α φk′ φk′ = L Lj,αPk′ /dA. j,α j,α j,α ⊗ | ih | j,α [ i adH , A]= A−i[H,A] A+i[H,A], (C1) (B6) − D 2D − 2D we have Inserting these into (B1) and (B2) we then obtain 1 [ i ad|jihj|, |1ih2|]= |1ih2|−i|1ih2|δj1+i|1ih2|δj2 ˜ − D 2D ( j )( ) = ( j )( ), (B7) 1 P◦ K ◦P ··· K ◦P ··· . ( ˜ )( ) = ( )( ) |1ih2|+i|1ih2|δj1−i|1ih2|δj2 P◦ Dj ◦P ··· Dj ◦P ··· − 2D +2 L ( )L† , (B8) (C2) X j,αT ··· j,α α For j = 1, 2 it trivially vanishes, while for j = 1, 2 we get 1 (6 1 i 2 1 i 2) = 0, where we have used that is 2 | ∓ | − | ± | D|1ih2| = α 2 . This commutativity implies that we DαA | | DA ˜ † can generate every i adu(d) (see [23]) and thus every ( j )( ) = ( j )( )+2 Lj,α ( )Lj,α, − P◦ L ◦P ··· L ◦P ··· X T ··· Ad Ad † = † for any U U(d). α U ◦ D|1ih2| ◦ U DU|1ih2|U ∈ (B9) (j) Taking unitaries U j = k∈I ck k for j = 1, 2 and where is the super-operator = 1, 2, 3, 4 , we have| i P | i T I { } 1 (1) ∗ (2) (1) (2) ∗ † = (c ) c c (c ) . ( )= 11A TrA[( )] ( ), (B10) DU|1ih2|U X i j k l D|iihj|,|kihl| T ··· dA ··· −P ··· i,j,k,l∈I (C3) with Tr [( )] indicating the partial trace over the aux- A ··· We numerically verified that all the operators of the form iliary system A and 11A being the identity operator on (28) for = 1, 2, 3, 4 can be obtained by linear com- the associated Hilbert space. I { } binations of U|1ih2|U † and i adH with different U and In order to prove (15) let us now focus on the evolu- D − H = i,j∈I hij i j on = 1, 2, 3, 4 . The same argu- (Lj ) P | ih | I { } tion induced by CPTP map Φt,∞ in (13) associated with ment applies to any quartets = i,j,k,l , and all the the jth element of L on a generic density matrix ρ(0) of operators of the form (28) forI all {are available.} Then, I (Lj ) every Lindbladian can be given as a linear combination the joint system Q + A, i.e., ρ(t)=Φt,∞ ρ(0). We are interested in the dynamics of the reduced density matrix of those operators, i.e., of Q, i.e., = [c+ ( + ) L X ijkl D|iihj|,|kihl| D|kihl|,|iihj| i,j,k,l (Lj ) − ρQ(t) = TrA[ρ(t)] = TrA[Φt,∞ ρ(0)]. (B11) + c i( )] (C4) ijkl D|iihj|,|kihl| − D|kihl|,|iihj| ± By taking the first derivative with respect to t and using with some coefficients cijkl. 7

[1] G. Dirr and U. Helmke, Lie Theory for Quantum Control, Cambridge, 2000). GAMM-Mitt. 31, 59 (2008). [14] G. Lindblad, On the Generators of Quantum Dynamical [2] J. S. Glaser, U. Boscain, T. Calarco, P. C. Koch, W. Semigroups, Commun. Math. Phys. 48, 119 (1976). K¨ockenberger, R. Kosloff, I. Kuprov, B. Luy, S. Schirmer, [15] P. Facchi and S. Pascazio, Quantum Zeno Subspaces, T. Schulte-Herbr¨uggen, D. Sugny, and K. F. Wilhelm, Phys. Rev. Lett. 89, 080401 (2002). Training Schr¨odinger’s Cat: Quantum Optimal Control, [16] A. Royer, Wigner Function in Liouville Space: A Canon- Eur. Phys. J. D 69, 279 (2015). ical Formalism, Phys. Rev. A 43, 44 (1991). [3] D. L. Elliott, Bilinear Control Systems: Matrices in Ac- [17] J. Schwinger, The Algebra of Microscopic Measurement, tion (Springer, Heidelberg, 2009). Proc. Natl. Acad. Sci. USA 45, 1542 (1959). [4] D. D’Alessandro, Introduction to Quantum Control and [18] J. Schwinger, Quantum Kinematics and Dynamics Dynamics (Chapman & Hall/CRC, Boca Raton, FL, (Perseus, New York, 1991), p. 26. 2008). [19] P. Zanardi and L. Campos Venuti, Coherent Quantum [5] G. Dirr, U. Helmke, I. Kurniawan, and T. Schulte- Dynamics in Steady-State Manifolds of Strongly Dissi- Herbr¨uggen, Lie-Semigroup Structures for Reachability pative Systems, Phys. Rev. Lett. 113, 240406 (2014). and Control of Open Quantum Systems: Kossakowski- [20] M. Ippoliti, L. Mazza, M. Rizzi, and V. Giovannetti, Per- Lindblad Generators Form Lie Wedge to Markovian turbative Approach to Continuous-Time Quantum Error Channels, Rep. Math. Phys. 64, 93 (2009). Correction, Phys. Rev. A 91, 042322 (2015). [6] C. O’Meara, G. Dirr, and T. Schulte-Herbr¨uggen, Illus- [21] P. Zanardi and L. Campos Venuti, Geometry, Robust- trating the Geometry of Coherently Controlled Unital ness, and Emerging Unitarity in Dissipation-Projected Open Quantum Systems, IEEE Trans. Autom. Control Dynamics, Phys. Rev. A 91, 052324 (2015). 57, 2050 (2012). [22] W. K. Wootters and B. D. Fields, Optimal State- [7] I. Kurniawan, Controllability Aspects of the Lindblad- Determination by Mutually Unbiased Measurements, Kossakowski Master Equation: A Lie-Theoretical Ap- Ann. Phys. (N.Y.) 191, 363 (1989). proach, Ph.D. thesis, Bayerischen Julius-Maximilians- [23] D. Burgarth and K. Yuasa, Quantum System Identifica- Universit¨at W¨urzburg, W¨urzburg, 2009. tion, Phys. Rev. Lett. 108, 080502 (2012). [8] I. Kurniawan, G. Dirr, and U. Helmke, in Proc. 19th [24] D. Bacon, A. M. Childs, I. L. Chuang, J. Kempe, D. W. Int. Symp. Math. Theory Networks Syst. (MTNS 2010), Leung, and X. Zhou, Universal Simulation of Markovian Budapest, Hungary (Budapest, 2010), pp. 2333–2338. Quantum Dynamics, Phys. Rev. A 64, 062302 (2001). [9] D. K. Burgarth, P. Facchi, V. Giovannetti, H. Nakazato, [25] R. Sweke, I. Sinayskiy, D. Bernard, and F. Petruccione, S. Pascazio, and K. Yuasa, Exponential Rise of Dynami- Universal Simulation of Markovian Open Quantum Sys- cal Complexity in Quantum Computing through Projec- tems, Phys. Rev. A 91, 062308 (2015). tions, Nat. Commun. 5, 5173 (2014). [26] D. I. Bondar, R. Cabrera, R. R. Lompay, M. Y. Ivanov, [10] C. Arenz, D. Burgarth, P. Facchi, V. Giovannetti, H. and H. A. Rabitz, Operational Dynamic Modeling Tran- Nakazato, S. Pascazio, and K. Yuasa, Universal Control scending Quantum and Classical Mechanics, Phys. Rev. Induced by Noise, Phys. Rev. A 93, 062308 (2016). Lett. 109, 190403 (2012). [11] P. Facchi and S. Pascazio, Quantum Zeno Dynamics: [27] D. I. Bondar, R. Cabrera, D. V. Zhdanov, and H. A. Ra- Mathematical and Physical Aspects, J. Phys. A: Math. bitz, Wigner phase-space distribution as a wave function, Theor. 41, 493001 (2008). Phys. Rev. A 88, 052108 (2013). [12] D. Orsucci, D. Burgarth, P. Facchi, H. Nakazato, S. Pas- [28] S. Machnes and M. B. Plenio, Surprising Interactions of cazio, K. Yuasa, and V. Giovannetti, Hamiltonian Purifi- Markovian noise and Coherent Driving, arXiv:1408.3056 cation, J. Math. Phys. 56, 122104 (2015). [quant-ph] (2014). [13] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge University Press,