{'fraction_non_alphanumeric': 0.04793858058275904, 'fraction_numerical': 0.026478055857816468, 'mean_word_length': 4.62687074829932, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 3, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Demonstration of dynamical control of three-level open systems with a superconducting qutrit

We propose a method for the dynamical control in three-level open systems and realize it in the experiment with a superconducting qutrit. Our work demonstrates that in the Markovian environment for a relatively long time (3 us), the systemic populations or coherence can still strictly follow the preset evolution paths. This is the first experiment for precisely controlling the Markovian dynamics of three-level open systems, providing a solid foundation for the future realization of dynamical control in multiple open systems. An instant application of the techniques demonstrated in this experiment is to stabilize the energy of quantum batteries.

Introduction.-The control of quantum systems is always an extremely interesting and essential topic among diverse types of quantum information tasks [1,2]. With the gradual development of experimental techniques, decoherence, mainly resulted from energy relaxation and dephasing, has become a primary barrier lying on the way to experimentally implement the quantum control. To clear this barrier, the researchers choose to take measures to shorten the evolution time and have achieved remarkable successes in the experiments, such as, preparing 20-qubit entanglement [3,4] and simulating quantum walk on 62-qubit processor [5]. Meanwhile, some researchers attempt on constructing qubits that are insensitive to decoherence [5][6][7][8]. So far, the relaxation time T 1 and dephasing time T 2 of superconducting qubit have both been extended to 200 µs [8]. Indeed, the methods of shortening the evolution time and making more perfect qubits naturally reduce the impact of decoherence rather than avoiding it, which still process quantified errors. For realizing the leap of quantum error rates from 10 −3 to 10 −15 [9, 10] and further accomplishing near-term quantum works, it is necessary to explore the control of open quantum systems, tailoring the decoherence as one of the effective elements involved in the quantum control.
The earlier open-system quantum control theory draws supports from a kind of nuclear magnetic resonance technology, called spin-echo [11], in which a π pulse is inserted midway in the evolution time to increase T 2 . However, this inserted π pulse will change the systemic dynamics. Several spin-echo-like techniques, such as bangbang control [12], dynamical decoupling technique [13], and parity-kicks [14], process the same problems with the destruction of dynamics. So these works generally served as techniques for storing coherence of quantum systems [15][16][17][18][19][20]. Researchers alternatively develop an open-system adiabatic theorem [21] to fully control the dynamics. This method utilizes super operators to describe the open systems and can approximately predict the density matrix at each moment. Since the adiabatic theorem needs a long evolution time to guarantee the adiabatic approximation, researchers further try to build a shortcut to adiabaticity (STA) (similar to that in the closed systems [22][23][24][25][26][27][28][29][30][31][32][33][34][35]).
By virtue of super operators, Sarandy et al. [36] have proposed a concept about open-system dynamical invariants to build a STA in open systems. They constructed some 4-dimensional invariants for engineering two-level open systems, nevertheless, with several simplifications, such as, either considering energy relaxation or dephasing. Finding proper invariants will be complicated when someone simultaneously takes energy relaxation and dephasing into account, let alone expansion to three-level open systems with 9-dimensional invariants. Therefore other control methods are eagerly in need of discovery.
Recently, Medina and Semião [37] have directly set the density operator with time-dependent parameters and then submitted them to the Markovian master equation [38], which gives appropriate functions of pulses to control the populations in two-level open systems. The error rates of populations control in Ref. [37] under serious decoherence are almost zero, an exciting result. Some follow-up studies [39,40] have made a lot of additions to the application scenarios, but still remaining in two-level open systems. With the general trend of quantum tasks towards multiplication [41], it is of great significance to explore the dynamics of three (or more)-level open systems.
In this letter, we propose a method for the dynamical control in three-level open systems, and demonstrate it with a superconducting circuit platform [42][43][44][45][46][47]. Two time-dependent parameters are presupposed in the Markovian master equation to singly control the populations or coherence, and then the analytical functions of driving pulses are given to complete desired dynamical control. For proving the correctness of this method, we apply the designed driving pulses to a superconducting Xmon qutrit and measure populations of excited states |1 and |2 , or coherence between ground and excited states. The experimental results are in good agreement with the theoretical design for the dynamical con-trol. This is the first experiment to precisely control the Markovian dynamics of three-level open systems, where the decoherence is effectively dominated. In addition, an interesting freezing phenomenon of populations is observed. That is, in the late stages of evolutions, the populations stabilize at specific values for a relatively long time (>1.8 µs) in the presence of decoherence. An instant application scenario of this freezing phenomenon is the quantum battery (QB) [48][49][50][51][52][53][54][55][56][57][58]. One can steadily lock the energy of QBs for a comparatively long time (>1.2 µs), resisting the relaxation of excitation, and therefore building prototypes of QBs with more stable energy.
Method.-Assume that a three-level system consists of bases |0 , |1 , and |2 . The systemic Hamiltonian in the interaction picture is (h = 1 hereafter) H i (t) = Ω 01 (t)|0 1|+Ω 12 (t)|1 2|+H.c., where Ω 01 (t) and Ω 12 (t) respectively represent the Rabi frequencies of the designed driving pulses for controlling transitions |0 ↔ |1 and |1 ↔ |2 . We assume the density matrix of this three-level system being [basis order {|2 ,|1 ,|0 } and dropping (t) henceforth] in which f 1 , f 2 , h 1 , h 2 , and h 3 are all time-dependent real functions. Here f 1 and f 2 describe the populations of states |1 and |2 , respectively. Additionally, h 1 , h 2 , and h 3 describe the coherence of this three-level system. Equation 1 (5 degrees of freedom) is not the most generalized case for a three-level density operator (8 degrees of freedom), however sufficient as a demonstration. The Markovian master equation [38] here iṡ where Lindblad operators are L k = √ γ k |k k| and L k+2 = √ Γ k |k − 1 k|, with Γ k and γ k being the rates of energy relaxation and dephasing of state |k , respectively (k = 1, 2). By utilizing Eqs. (1, 2), we can derive the Rabi frequencies to control the dynamics of this threelevel system, yielding It is noteworthy that there are three constraint equationṡ Up to now, the populations or coherence of the threelevel system can be effectively controlled by the interaction Hamiltonian H i , i.e., by Rabi frequencies Ω 01 and Ω 12 . A specific example is, an evolution from the ground state |0 to a final state with populations respectively being P 0 (t f ), P 1 (t f ), and P 2 (t f ) in states |0 , |1 , and |2 , can be achieved by adjusting f 1 = f P 1 (t f ) and 27,37], with a = 50/t f determining the gradient of the transformation. By submitting f 1 and f 2 into Eqs. (3, 4), one can obtain Rabi frequencies Ω 01 and Ω 12 to accomplish the desired transformation of populations.
An essential process for this experiment is to precisely measure the coefficients of decoherence, γ k and Γ k . According to the above form of Lindblad operators, one can deduce Γ 1 = 1/T 01 1 , γ 1 = 2/T 01 2 − Γ 1 , Γ 2 = 1/T 12 1 , and 4), preset γ k = Γ k = 0] pulses, respectively. Intuitively, the populations are controlled more precisely in the outer layer, as compared to the inner one. More rigorously, we define a standard deviation to describe the error rate of the population control where P i (t f ) ideal and P i (t f ) exp. are the ideal and experimental values of the population in |i , respectively. In the outer layer of Fig. 3, the populations control achieves an error rate of 1.02%, not small [60] because we intentionally extended the evolution time (3 µs). In contrast, if the control is performed for such a long time by applying closed-system pulses, the error rate reaches 7.49%. This stark difference in error demonstrates the effectiveness of the present control method. We also measure the results of 29 different controls of populations, 6 of them shown in Fig. 4 and all displayed in the Supplemental Material [61]. The error rates of the controls in Fig. 4 are around 1%, which we believe can be further improved with more stable superconducting qutrits [8].
Another exceptional physical phenomenon is that the populations seem to be frozen after 1.8 µs in the outer layers of Fig. 3. This freezing phenomenon of populations does not mean that the driving pulses Ω 01 and Ω 12 have stopped. On the contrary, it is the driving pulses [see the insets of Fig. 2(b)] we applied that caused this freezing phenomenon to occur. Specifically, such a freezing phenomenon arises from the interplay between the dynamics induced by the continuous microwave drives and the two decoherence channels characterized by T 1 and T 2 (both available for |0 ↔ |1 and |1 ↔ |2 transitions), similar to the cases for generation of steady states in most open systems [36,[62][63][64][65][66][67]. But here it differs significantly in that both the energy relaxation and the dephasing are involved in the nonequilibrium dynamical processes and together help freeze the states of three-level systems, as compared to the previous ones which generally consider only one decoherence channel [36,66,67].
Results of the coherence control. Since only two free variables, Ω 01 and Ω 12 in Eqs. (3, 4) can be controlled in the experiment, we can only control two parts of coherence. Here we choose h 3 and h 2 as an example for illuminating the coherence control. Corresponding variations of Eqs. (3, 4) and the forbidden zone of coherence control are shown in the Supplemental Material [61]. Ac- cording to the same preset of parameters and intermediate function f , we show 6 groups of results of the coherence control in Fig. 5 (all the 36 groups of results in the Supplemental Material [61]). The error rate here is larly, the microwaves protect the coherence (lasting 1.2 µs) from the influence of dephasing. If there are no microwaves applied, roughly, after the same 1.2 µs, h 3 (1.2µs)/h 3 (0) ∼ exp(−1.2µs/T 01 2 ) ∼ 0.8 and h 2 (1.2µs)/h 2 (0) ∼ exp(−1.2µs/T 12 2 ) ∼ 0.45, the coherence will be seriously damaged. In contrast, the error rates are all within 2% when we apply microwaves for the coherence control.
Application.-Quantum battery [48][49][50][51][52][53][54][55][56][57][58]. The freezing phenomenon seems helpful to stabilize the energy of the QB for a fairly long time, especially in the energy-storing stage of the QB. For proving this, we design a pulse for a qutrit in the experiment to simulate the charging, storing, discharging processes of the QB (see Fig.  6). A considerably long period, about 1.2 µs, of energy [ǫ = (ω 1 − ω 0 )P 1 ] stability was observed. This stability is hard to achieve by applying a pulse designed without considering decoherence because energy relaxation of the superconducting qutrit will exponentially drop the excitation down. Therefore, the present method may contribute to the construction of more stable QBs.
Conclusion.-We have proposed a method to control the dynamics of three-level open systems and realized it in the experiment with a platform of a superconducting Xmon qutrit [42][43][44][45][46][47]. The populations and coherence could be singly controlled with error rates around 1% under the influence of decoherence for a relatively long time, 3 µs, close to the dephasing time T 12 2 of the qutrit. In some situations where the control is more successful, the error rates can even be less than 0.3%. We believe these error rates can be further reduced with more stable superconducting qutrits [8].
Additionally, an interesting freezing phenomenon of populations was observed. The designed microwave pulses just offset the impact of decoherence and visually freeze the populations. We then applied this phenomenon to make more stable prototypes of QBs, whose energy can hold 1.2 µs with only one charging. Moreover, this freezing phenomenon strongly proves that the Markovian master equation precisely describes the dynamics of three-level open systems, as demonstrated here with a superconducting Xmon qutrit that possesses the specially intrinsic decoherence rates. The present work provides a positive prospect of accurately realizing the dynamical control of three (or more)-level open systems by using the adjustable driving pulses in the Markovian environment.
Note that all the data points of experiments were averaged over 12000 times and the readout errors of the qutrit have been corrected. We thank Ye-Hong Chen and Xin Wang in RIKEN for discussions. This work was supported by the National Natural Science Foundation of China (   The entire control and readout layout of the experiment is shown in Fig. S1. Electronic devices for qutrit readout, Josephson parametric amplifier (JPA) control, and qutrit control are displayed from top to bottom. The readin and XY control signals are the mixtures of lowfrequency signals [from two independent digital-to-analog converter (DAC) channels I and Q] and high-frequency signals (from microwave sources) to achieve nanosecond fast tuning. On the other hand, the qutrit Z control signal is sent directly from the DAC without mixing, whose frequency can also be slowly tuned by the direct current (DC) bias line. Additionally, the output signal of the readout feeder is amplified by the JPA, high electron mobility transistor (HEMT), and room temperature amplifier, and then demodulated by the analog-to-digital converter (ADC). Four cryogenic unidirectional circulators are inserted between the JPA and the 4K HEMT to block the reflection and noise from outside. As for the JPA, it is pumped by an independent microwave signal source and a DC bias, which is converged by the bias tee. At different temperature stages of the dilution refrigerator, each control line is balanced with some attenuators and filters (intuitively shown in Fig. S1) to prevent unwanted noise from affecting the equipment.

II. QUTRIT READOUT AND CALIBRATION
Here we show the results of the qutrit readout in Fig.  S2, where blue, orange, and yellow points represent the results of applying I, X 01 , and X 01 X 12 to the qutrit, i.e., depicting the results of preparing |0 , |1 , |2 , respectively. The readout pulses are 1 µs-long and the repetition is 3000. It can be seen that these three states are clearly separated, though with several error transitions caused by the readout errors. Like that in Refs. [1, 2], the calibration matrix can be defined as (S1) * These two authors contributed equally to this work. † zbyang@fzu.edu.cn ‡ xia-208@163.com § t96034@fzu.edu.cn where F ii ′ (i, i ′ = 0, 1, 2) is the probability of measuring the qutrit in |i when it is prepared in |i ′ . Such that the measured probability is given by P m Equations (3, 4) in the main text are used to control the populations, with designable f k (k = 1, 2) and iterative h i+1 (i = 0, 1, 2). If we alternatively aim to control two parts of the coherence (e.g., h 2 and h 3 ) of systems, we need to alter Eqs. (3, 4) in the main text as in which h k+1 is designable and f k is iterative. Therefore one can design suitable h 2 and h 3 to accomplish desired coherence control.

IV. FEASIBLE AND INFEASIBLE AREAS OF THE COHERENCE CONTROL
For the control of h 2 and h 3 , similarly, there are feasible area (FA) and infeasible area (IFA), which mainly

V. TOMOGRAPHY OF THE COHERENCE CONTROL
The populations of three-level systems can be directly measured by the readout cavity. (In the experiments, only diagonal elements of the density matrix can be directly measured.) While the coherence can be indirectly measured by a simplified tomography, including 4 steps as follows, (i) U 1 = I,