# A Survey of Quantum Lyapunov Control Methods

Authors: Shuang Cong; Fangfang Meng
Journal: The Scientific World Journal
 (2013)

---

## Abstract

The condition of a quantum Lyapunov-based control which can be well used in a closed quantum system is that the method can make the system convergent but not just stable. In the convergence study of the quantum Lyapunov control, two situations are classified: nondegenerate cases and degenerate cases. For these two situations, respectively, in this paper the target state is divided into four categories: the eigenstate, the mixed state which commutes with the internal Hamiltonian, the superposition state, and the mixed state which does not commute with the internal Hamiltonian. For these four categories, the quantum Lyapunov control methods for the closed quantum systems are summarized and analyzed. Particularly, the convergence of the control system to the different target states is reviewed, and how to make the convergence conditions be satisfied is summarized and analyzed.


---

## Body

## 1. Introduction

The theory of quantum mechanics is one of the major discoveries in the history of science in the 20th century. It is a very important issue to study the properties of the quantum mechanical systems and their control. According to whether the system is isolated or not, a quantum mechanical system can be a closed system or an open system. In a closed quantum system, the evolution of the state is unitary. There are mainly two methods to describe the evolution of a closed quantum system’s states. They are the Schrödinger equationi|ψ˙(t)〉=(H0+∑k=1rHkuk(t))|ψ(t)〉 and the quantum Liouville equation iρ˙(t)=[H0+∑k=1rHkuk(t),ρ(t)], in which |ψ(t)〉 is the quantum state vector, ρ(t) is the density operator, H0 is the internal Hamiltonian, and Hk,(k=1,…,r) are the control Hamiltonians. In an open quantum system, the system interacts with the surroundings, thus the loss of the system’s information leads to the non-unitary evolution of the state. The most common method to describe the open system is the Lindblad master equation: iρ˙=[H0+∑k=1rHkuk(t),ρ]+LD(ρ), which is in fact the sum of a closed system and a dissipative term caused by the loss of the information or energy. Obviously, the research of the properties of the closed quantum system and its control is relatively simple. Moreover, there is a more important fact: the research of the closed quantum system is the basis of that of the open quantum system.Quantum control has attracted much attention in recent years and it has been found the potential applications in many fields such as atomic physics [1–4], molecular chemistry [5–9] and quantum information [10, 12]. Up to now, there have been many quantum methods, such as quantum optimal control [13–15], adiabatic control [16–18], the Lyapunov-based control [19–41], and optimal Lyapunov-based quantum control [42]. For the Lyapunov-based quantum control, it is relatively easy to design an analytical but not numerical control law, and the control system based on this control method is at least stable, so it has been a common control method.One of the major concerns of the Lyapunov control is choosing an appropriate Lyapunov function to design the control laws. Ordinarily, the control laws and the control effects are different when the Lyapunov functions are distinct. It’s a good idea to choose the Lyapunov function based on the geometrical and physical meanings. Usually, there are mainly three Lyapunov functions to be selected: the Lyapunov function based on the state distance [19–21, 26–32, 39–41], the state error [22, 23, 37, 38], and the average value of an imaginary mechanical quantity [24, 25, 39, 41]. The so-called imaginary mechanical quantity means that it is a linear Hermitian operator to be designed and may be not a physically meaningful observable quantity such as coordinate and energy. Among these three Lyapunov functions, the Lyapunov-based quantum control methods based on the state distance and state error only need to adjust the scale factors of the control laws. These two Lyapunov control methods are relatively simple and easy to grasp. The Lyapunov-based quantum control based on the average value of an imaginary mechanical quantity contains more adjustable parameters. So it is more flexible and also more complex at the same time.Generally speaking, the Lyapunov-based control method can only ensure that the control system is stable. The probability control in the quantum system requires us to design a control strategy which can make the system convergent, because a stable quantum control method may result in that the control system cannot reach the desired target state. Therefore another major concern of this control strategy for the closed quantum systems is the convergence of the control systems. So far, there have been the following research results on the convergence of the closed quantum systems [19, 23–25, 27–29, 33, 36].(I)
For the Schrödinger equation, the convergence conditions are as follows. (i) The internal Hamiltonian is strongly regular; (ii) All the eigenstates, which are different from the target state, are directly coupled to the target state for the Lyapunov control based on the state distance or the state error, or any two eigenstates are coupled directly for the Lyapunov control based on the average value of an imaginary mechanical quantity.(II)
For the quantum Liouville equation, the convergence conditions are as follows. (i) The internal Hamiltonian is strongly regular, and (ii) the control Hamiltonians are full connected.At first, the Lyapunov control method which could only ensure the convergence to the target eigenstate was studied [19, 23, 25]. Then, the target mixed state, which commutes with the internal Hamiltonian, was studied [33, 36]. Later, the convergence to the target superposition state was solved by means of designing the control laws in the interaction picture of the control system, and the convergence to the target mixed state which does not commute with the internal Hamiltonian was solved by using a unitary transformation.In fact, many actual systems do not satisfy the convergence conditions mentioned above, such as the time domain model of the selective excitation of the stimulated Raman scattering [43], the coupled two spin systems, and one-dimensional oscillators [44]. These systems are called nonideal systems and in the degenerate cases. For the degenerate cases, the convergence of the control systems was solved by introducing a series of implicit function perturbations and choosing an implicit Lyapunov function [37–40]. At first, the convergence to the target eigenstate was only guaranteed [37–40]. Then the convergence to the target mixed state which commutes with the internal Hamiltonian was solved by introducing a series of implicit function perturbations, and the convergence to the target superposition state and the target mixed state which does not commute with the internal Hamiltonian was solved by introducing a series of constant disturbances.The aim of this paper is to summarize and analyze the existing Lyapunov control methods for the nondegenerate and degenerate cases, respectively. Dividing the target state into four categories: the eigenstate, the mixed state which commutes with the internal Hamiltonian, the superposition state, and the mixed state which does not commute with the internal Hamiltonian state, we summarize the design methods of the control laws, analyze the convergence to the target state, and investigate how to make these conditions of the convergence be satisfied.The remainder of this paper is arranged as follows. In Section2, the research results for the nond-egenerate cases are summarized and analyzed. In Section 3, the research results for the degenerate cases are summarized and analyzed. Some concluding remarks are drawn in Section 4.


## 2. Non-Degenerate Cases

The design of the control laws and the analysis of the convergence are very important in the Lyapunov control method. The design of the control laws is based on the Lyapunov stability theorem, which is to design the control laws to make the selected positive semi-definite Lyapunov functionV(t) satisfy V˙(t)≤0. The convergence analysis of this control method is mainly based on the LaSalle invariance principle [44] for the autonomous systems, or the improved Barbalat lemma [45] for the non-autonomous systems.In this Section, we will summarize and analyze the research results on the convergence of the control system in the non-degenerate cases for the target state being an eigenstate, a mixed state which commutes with the internal Hamiltonian, the superposition state, or/and a mixed state which does not commute with the internal Hamiltonian state, respectively.


### 2.1. Target Eigenstate

It is convenient to use the bilinear Schrödinger equation to describe the control systems if the target state is a pure state. Consider theN-level closed quantum system governed by the following bilinear Schrödinger equation:
(1)i|ψ˙(t)〉=(H0+∑k=1rHkuk(t))|ψ(t)〉,
where |ψ(t)〉 is the quantum state vector, H0 is the internal Hamiltonian, Hk,(k=1,…,r) are control Hamiltonians, and uk(t),(k=1,…,r) are scalar and real control laws.In the following sections, in the case that the target state|ψf〉 is an eigenstate, that is, H0|ψf〉=λf|ψf〉, where λf is the eigenvalue of the internal Hamiltonian H0, the research results on the convergence are summarized and analyzed for the Lyapunov control based on the state distance, the state error, and the average value of an imaginary mechanical quantity, respectively.


#### 2.1.1. Lyapunov Control Based on State Distance

Consider the following Lyapunov function based on the state distance:(2)V1(|ψ〉)=12(1-|〈ψ∣ψf〉|2).
The time derivative of the Lyapunov function (2) is
(3)V˙1=-∑k=1ruk(t)·Im[〈ψ∣ψf〉〈ψf|Hk|ψ〉].
The control laws uk(t) which can make V˙1≤0 hold can be designed as
(4)uk(t)=Kkfk(Im(〈ψ∣ψf〉〈ψf|Hk|ψ〉)),(k=1,…,r),
where Kk>0, and yk=fk(xk),(k=1,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.The control laws designed in (4) can only ensure the control system (1) to be stable. One needs to do further study on the convergence of the control system. The control system governed by (1) is an autonomous system, whose convergence can be analyzed based on the LaSalle invariance principle [44]. According to the LaSalle invariance principle, as t→∞, any state trajectory will converge to the largest invariant set contained in the set E in which the states satisfy that the first order derivative of the Lyapunov function equals zero. In fact, the set E contains not only the target state but also other states, thus the system may converge to other states rather than the target state. The main idea to solve this problem is to add restrictions to make the set E as small as possible. Based on the LaSalle invariance principle, the convergence of the control system governed by (1) can be depicted by Theorem 1.Theorem 1 (see [19]).
Consider the control system governed by (1) with control fields uk(t) designed in (4). If (i) The system is strongly regular, that is, ωi′j′≠ωlm,(i′,j′)≠(l,m),i′,j′,l,m∈{1,2,…,N},ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) for any |ϕi〉≠|ψf〉,i∈{1,…,N}, there exists at least a k∈{1,…,r} such that 〈ϕi|Hk|ψf〉≠0, then any state trajectory will converge toward E1={eiθ|ψf〉,θ∈R}.From Theorem1, one can see that in the case that the target state is an eigenstate, if the control system governed by (1) satisfies the conditions (i)-(ii), the control system can converge to the equivalent state of the target eigenstate eiθ|ψf〉 from any initial pure state. These two conditions are relevant to the internal Hamiltonian and the control Hamiltonians, which are system parameters. Once the control system is determined, the Hamiltonians are fixed and cannot be changed by designing the control laws.


#### 2.1.2. Lyapunov Control Based on State Error

Consider the following Lyapunov function based on the state error:(5)V2(|ψ〉)=12〈ψ-ψf∣ψ-ψf〉.In the case of selecting the Lyapunov function based on the state error defined by (5), in order to facilitate to design the control laws based on the Lyapunov stability theorem, the drift item appeared in the first order time derivative of Lyapunov function, which is caused by the internal Hamiltonian, is needed to be eliminated. The existing solution is to add a global phase control item ω into the control system governed by (1). This method will not change the population distribution of the control system. Thus the dynamical equation (1) becomes
(6)i|ψ˙(t)〉=(H0+∑k=1rHkuk(t)+ωI)|ψ(t)〉.After some deduction, one can obtain the time derivative of the Lyapunov function (5) as
(7)V˙2=-(λf+ω)ℑ(〈ψf∣ψ〉)-∑k=1rℑ(〈ψf|Hk|ψ〉)uk(t).
The control laws which can make V˙2≤0 hold can be designed as
(8)ω=-λf+cf0(ℑ(〈ψf∣ψ〉)),(9)uk(t)=Kkfk(ℑ(〈ψf|Hk|ψ〉)),(k=1,…,r),
where Kk>0, and yk=fk(xk),(k=0,…,r) are the monotonic increasing functions through the coordinate origin of the plane xk-yk.Based on the LaSalle invariance principle, the convergence of the control system governed by (6) can be depicted by Theorem 2.Theorem 2 (see [19, 23]).
Consider the control system governed by (6) with control fields uk(t) designed in (9) and ω designed in (8). If (i) ωi′j′≠ωlm,(i′,j′)≠(l,m), i′,j′,l,m∈{1,2,…,N}, ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) for any |ϕi〉≠|ψf〉,i∈{1,…,N}, there exists at least a k∈{1,…,r} such that 〈ϕi|Hk|ψf〉≠0. Then any state trajectory will converge toward E2={eiθ|ψf〉,θ∈R}.From Theorem2, one can see that for the case that the target state is an eigenstate, if the control system governed by (6) satisfies the conditions (i)-(ii), the control system can also converge to the equivalent state of the target eigenstate eiθ|ψf〉 from any initial pure state.


#### 2.1.3. Lyapunov Control Based on Average Value of an Imaginary Mechanical Quantity

Consider the following Lyapunov function based on the average value of an imaginary mechanical quantity:(10)V3(|ψ〉)=〈ψ|P|ψ〉,
where the imaginary mechanical quantity P is a positive definite Hermitian operator.The first order time derivative of the Lyapunov function (10) can be obtained as
(11)V˙3=i〈ψ|[H0,P]|ψ〉+i∑k=1r〈ψ|[Hk,P]|ψ〉uk.Set[H0,P]=0 such that the drift term in the right side of (11) can be eliminated. In order to ensure V˙3≤0, one can design uk(t) as
(12)uk(t)=-Kkfk(i〈ψ|[Hk,P]|ψ〉),(k=1,…,r),
where Kk>0, and yk=fk(xk),(k=1,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.Then based on the LaSalle invariance principle, all the state trajectories of the system will converge to the invariant set contained in the setE in which V˙=0 holds. Denote the state at the time t as |ψ(t)〉=∑l=1Ncl(t)|ϕl〉, where cl(t) is the coefficient corresponding to the lth eigenstate |ϕl〉. For the Schrödinger equation governed by (1), if the control system is strongly regular and any eigenstate is directly coupled to all other eigenstates, that is, for any j≠l, j,l∈{1,2,…,N}, there exists a k∈{1,…,r} such that 〈ϕj|Hk|ϕl〉≠0, then one can deduce that V˙3(t)=0 holds for all t≥t0,t0∈R is equivalent to
(13)(Pl-Pj)cj(t0)cl*(t0)=0,(l,j=1,…,N),
where Pl and Pj are the lth and jth eigenvalues of P, respectively.If the target state|ψf〉 is an eigenstate, and all the eigenvalues of P are designed mutually different, that is, Pl≠Pj for any l≠j,(l,j=1,2…,N), then (13) is equivalent to
(14)cjcl*=0,(l,j=1,…,N).
Equation (14) implies that there is at most one cj≠0,(j=1,…,N); that is, the system will converge to an eigenstate with t→∞. Thus the convergence of the control system (1) can be depicted by Theorem 3.Theorem 3 (see [19, 24, 25]).
Consider the control system governed by (1) with the control fields uk(t) designed in (12). If (i) ωi′j′≠ωlm,(i′,j′)≠(l,m),i′,j′,l,m∈{1,…,N},ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) for any i≠j,i,j∈{1,…,N}, there exsits at least a k∈{1,…,r} such that 〈ϕi|Hk|ϕj〉≠0; (iii) [H0,P]=0; (iv) for any l≠m∈{1,2,…,N}, Pl≠Pm holds, where Pl is the lth eigenvalue of P. Then any state trajectory will converge toward E3={eiθ|ϕi〉,i∈{1,…,N},θ∈R}.From Theorem3, one can see that the control system will converge from any initial pure state to an eigenstate which may not be the target eigenstate. In order to make the system converge to the target eigenstate |ψf〉 from any initial pure state |ψ0〉, as V˙3(t)≤0, one can add a restriction as
(15)V3(|ψf〉)<V3(|ψ0〉)<V3(|ψother〉),
where |ψ0〉 is the initial state, |ψother〉 represents any other state in the set E3 except the target state.In such a way, any state trajectory of the system will converge toeiθ|ψf〉 from any initial pure state |ψ0〉.Next, let us analyze how to make these convergence conditions be satisfied. Conditions (i) and (ii) are only relevant to the internal HamiltonianH0 and the control Hamiltonians Hk,(k=1,…,r) which cannot be changed by designing appropriate control laws. Condition (iii) means that P and H0 have the same eigenstates. In order to make condition (iii) be satisfied, the eigenvalues of P can be designed as
(16)P=∑j=1NPj|ϕj〉〈ϕj|.
Because P should be positive definite, one needs to design Pi>0,(i=1,…,N). The restriction (15) can be satisfied by means of designing an appropriate P. For the restriction (15), Grivopoulos and Bamieh proposed a design principle of P to make V3(|ψf〉)<V3(|ψother〉) hold. This design principle of P can be depicted by Proposition 4.Proposition 4 (see [25]).
With the constraint condition〈ψ∣ψ〉=1, the set of critical points of the Lyapunov function V3(|ψ〉)=〈ψ|P|ψ〉 is given by the normalized eigenvectors of P. The eigenvectors with the largest eigenvalue are the maxima of V3, the eigenvectors with the smallest eigenvalue are the minima and all others are saddle points.According to Proposition4, in order to make V3(|ψf〉)<V3(|ψother〉) hold, Pi>Pf,(i=1,…,N,Pi≠Pf) needs to be designed, where Pf is the eigenvalue of P corresponding to |ψf〉. Then let us consider the whole restriction (15); it is an attraction problem. If the eigenvalues of P except Pf are close together, the attraction region will be very large. For the limiting case Pi=Pj>Pf>0,(i≠j,Pi,Pj≠Pf), the attraction region will be the whole state space. Thus the design principle of P is Pi>Pf,(i=1,…,N,Pi≠Pf) and to make Pi,Pj≠Pf close together.In conclusion, for the target state being an eigenstate, the design principle ofP is as follows: (i)
Pl≠Pm>0 for any l≠m∈{1,2,…,N},(ii)
Pi>Pf,(i=1,…,N,Pi≠Pf), Pi,Pj≠Pf are close together,(iii)
Equation (16): P=∑j=1NPj|ϕj〉〈ϕj|.From the above analyses, we can conclude that the Lyapunov control method based on the imaginary mechanical quantity proposed in [19, 24, 25] can only ensure the convergence to an eigenstate, but cannot guarantee the convergence to the target eigenstate from any initial pure state. However, if there exists a P to make the restriction (15) hold, then any state trajectory of the system will converge to the equivalent state of the target eigenstate eiθ|ψf〉 from any initial pure state.


#### 2.1.4. Relations between Three Lyapunov Functions

In the Liouville space, the Hilbert-Schmidt distance between two density operatorsρ1 and ρ2 is
(17)dHS(ρ1,ρ2)=tr(ρ1-ρ2)2.The inner product of two operatorsA and B is defined as 〈〈A|B〉〉=tr(A†B), where the operation A† refers to the conjugate transpose of A. Because ρ=|ψ〉〈ψ|, the square of the Hilbert-Schmidt distance between the density operator ρ and the target density operator ρf can be deduced as
(18)dHS2(ρ,ρf)=2(1-|〈ψ∣ψf〉|2).One can conclude from (18) that the Lyapunov function based on the state distance V1(|ψ〉)=(1/2)(1-|〈ψ∣ψf〉|2) and the state error V2(|ψ〉)=(1/2)〈ψ-ψf∣ψ-ψf〉 are equivalent.For the Lyapunov function based on the average value of an imaginary mechanical quantityV3(|ψ〉)=〈ψ|P|ψ〉, if the imaginary mechanical quantity P=(1/2)(I-|ψf〉〈ψf|), V3(|ψ〉) will become V1(|ψ〉) [19]. Therefore V1(|ψ〉) is formally a special case of V3(|ψ〉).In fact,V1(|ψ〉), V2(|ψ〉), and V3(|ψ〉) can be unified in the following quadratic Lyapunov function:
(19)V4(|ψ〉)=〈ψ-αψf|Q|ψ-αψf〉Next, we consider three cases as follows:(i)
α=0,Q=(1/2)(I-|ψf〉〈ψf|), V4(|ψ〉) will reduce to V1(|ψ〉);(ii)
α=1,Q=(1/2)I, V4(|ψ〉) will reduce to V2(|ψ〉);(iii)
α=0, V4(|ψ〉) will reduce to V3(|ψ〉).We can see thatV1(|ψ〉), V2(|ψ〉), and V3(|ψ〉) are special cases of V4(|ψ〉).From the above analyses, we can conclude that all these three Lyapunov control methods can converge to the equivalent state of the target eigen state from any initial pure state. The Lyapunov control methods based on the state distance and the state error have only one adjustable parameter. So these two methods are very easy to grasp and very simple. The Lyapunov control based on the average value of an imaginary mechanical quantity has more adjustable parameters. So it is more flexible and also more complex at the same time. For the target eigenstate,V1(|ψ〉) and V2(|ψ〉) are equivalent, so the Lyapunov control based on the state distance and the state error have similar control effects. Because V1(|ψ〉) is formally a special case of V3(|ψ〉), generally, the control effect of the Lyapunov control based on the average value of an imaginary mechanical quantity is better than that of the state distance. At least, it can get the same control effect as the Lyapunov control of the state distance.


### 2.2. Target Mixed State Which Commutes with the Internal Hamiltonian

The bilinear Schrödinger equation cannot describe the mixed state. Thus for the target mixed state, it needs to use the quantum Liouville equation which can describe the evolution of any state of a closed quantum system. Consider theN-level closed quantum system governed by the following quantum Liouville equation:
(20)iρ˙(t)=[H0+∑k=1rHkuk(t),ρ(t)],
where ρ(t) is the density operator.Consider the Lyapunov function based on the average value of an imaginary mechanical quantity:(21)V5(ρ)=tr(Pρ).
By means of setting [H0,P]=0, the first order time derivative of the Lyapunov function (21) can be deduced as
(22)V˙5=-itr([P,H0]ρ)-i∑k=1rtr([P,Hk]ρ)uk.In order to ensureV˙5≤0,  one can design uk(t) as
(23)uk(t)=iεktr([P,Hk]ρ),(k=1,…,r),
where εk∈R,εk>0.Next, let us analyze the convergence to the target state. For the control system governed by (20) in the non-degenerate case, one can deduce that V˙5(t)=0 holds for all t≥t0,t0∈R is equivalent to
(24)(Hk)jl(Pl-Pj)ρlj(t)=0,j,l=1,…,N,j<l,
where ρlj(t) is the (l;j)th element of the state ρ(t). If the target state ρf commutes with H0, that is, [H0,ρf]=0, and all the eigenvalues of P are designed mutually different, that is, Pl≠Pj for any l≠j,l,j=1,2,…,N, then (24) is equivalent to
(25)ρlj=0,(l,j=1,…,N)
which implies that the system will converge to a state which commutes with H0. Thus based on the LaSalle invariance principle, the convergence of the control system governed by (20) can be depicted by Theorem 5.Theorem 5 (see [37]).
Consider the control system governed by (20) with the control field uk(t) designed in (23). If (i) the internal Hamiltonian is strongly regular, that is, ωi′j′≠ωlm,(i′,j′)≠(l,m),i′,j′,l,m∈{1,…,N},ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) the control Hamiltonians are full coupled; that is, ∀j≠l, for k=1,…,r, there exists at least a (Hk)jl≠0; (iii) [H0,P]=0; (iv) Pl≠Pm,l≠m∈{1,2,…,N}, where Pl is the lth eigenvalue of P. Then any state trajectory will converge toward E4={ρ∣(ρ)ij=0,i,j∈{1,…,N}}, where (ρ)ij is the (l,l)th element of ρ.From Theorem5, one can see that if the control system satisfies the conditions (i)–(iv), the control system will converge from any initial state to a state that commutes with the internal Hamiltonian, which may not be the target state. Next, what we need to do is to make the control system converge to the target state.Denote the state inE4 as ρE4; then [ρE4,H0]=0 holds which implies that ρE4 and H0 have the same eigenstates. Since the evolution of ρ(t) is unitary, ρ(t) for t≥0 are isospectral. So the eigenvalues of ρE4 are a permutation of the eigenvalues of ρ0. Thus the set E4 has countable elements. If the initial state is generic, that is, the eigenvalues of the initial state are mutually different, the set E4 will have N! elements. For the target state ρf which commutes with the internal Hamiltonian, that is, [ρf,H0]=0, in order to make the system converge to the target state ρf which commutes with H0 from any initial state ρ0, Kuang and Cong proposed a restriction as
(26)V5(ρf)<V5(ρ0)<V5(ρ other),
where ρ other represents any other state in the set E4 except the target state ρf.In such a way, if there exists aP to make the restriction (26) hold, any state trajectory of the system will converge to the target state ρf which commutes with H0 from any initial state ρ0.For how to make conditions (iii)-(iv) to be satisfied, please read Section2.1.3. One can deduce the design principle of P such that V5(ρf)<V5(ρ other), the result can be depicted by Proposition 6.Proposition 6.
If(ρf)ii<(ρf)jj,1≤i,j≤N holds, then one can design Pi>Pj; if (ρf)ii=(ρf)jj,1≤i,j≤N holds, then one can design Pi≠Pj; else if (ρf)ii>(ρf)jj,1≤i,j≤N holds, then one can design Pi<Pj, thus V5(ρf)<V5(ρ other) holds, where (ρf)ii is the (i,i)th element of ρf.It is difficult to designP such that (26) holds for any initial state ρ0 and any target state ρf which satisfies [ρf,H0]=0. One possible method is to introduce a series of implicit function perturbations into the control laws, this method will be presented in Section 3.3.For the target stateρf which commutes with the internal Hamiltonian, the design principle of P is as follows: (i)
Pl≠Pm>0 for any l≠m∈{1,2,…,N};(ii)
DesignP according to Proposition 6;(iii)
Equation (16): P=∑j=1NPj|ϕj〉〈ϕj|.In conclusion, if the control system satisfies the conditions (i)–(iv) in Theorem5 and there exists a P to make the restriction (26) hold, any state trajectory of the system will converge to the target mixed state ρf which commutes with the internal Hamiltonian from any initial state ρ0.


### 2.3. Target Superposition State

From Section2.1.3, one can see that for the control system governed by (1), V˙(t)=0 holds for all t≥t0,t0∈R is equivalent to (Pl-Pj)cjcl*=0,(l,j=1,…,N). If the target state is a superposition state as |ψf〉=∑d=k1kocd|ϕd〉,(1≤k1<ko≤N), one can design Pk1=Pk2=⋯=Pko and other eigenvalues of P are mutually different, then (Pl-Pj)cjcl*=0,(l,j=1,…,N) is equivalent to
(27)cjcl*=0,(l,j=1,…,N,l,j≠k1,…,ko).
Equation (27) implies that the system will converge to the set E5=span{|ϕk1〉,|ϕk2〉,…,|ϕko〉,|ϕl〉;l∈{1,…,N},l≠k1,…,ko}, which means that the set E5 contains infinite elements. Therefore this control method cannot ensure the system converge to the target superposition state by adding the restriction defined by (15). But this control method can ensure that the system will converge to the superposition of the eigenstates corresponding to target state. In fact, once one element of P changes, all the populations of the levels will change accordingly. Therefore when there are not very many eigenstates corresponding to target state, the system maybe converge to the target state by regulating the eigenvalues of P. Otherwise, the system maybe cannot converge to the target state. Consider the extreme situation in which the target state is the superposition of all the eigenstates. According to the design principle of P, we should design P=cpI, where cp is a real number, and I is the unit matrix. Obviously, in this special case, the control method proposed in Section 2.1.3 will become invalid.From Section2.2, one can see that, for the control system governed by (20), V˙(t)=0 holds for all t≥t0,t0∈R is equivalent to (Hk)jl(Pl-Pj)ρlj(t)=0,j,l=1,…,N,j<l. Consider the target state which does not commute with H0, which includes the superposition state and the mixed state which does not commute with H0. Without loss of generality, assume (ρf)3N=(ρf)N3≠0; one can design P3=PN; thus the system will converge to E6={ρ|(ρ)ij=0,(i,j)≠(3,N),(i,j)≠(N,3),i,j∈{1,…,N}}. Because the set E6 has infinite elements, this control method also cannot ensure the system will converge to the target state by adding the restriction defined by (26). But when the target state does not have so many nonzero off-diagonal elements, the system maybe converges to the target state by regulating the diagonal elements of P. When there are many nonzero off-diagonal elements in the target state, the degree of freedom of P may be not enough.From the above analyses, we can conclude that by means of using the control methods proposed in Sections2.1 and 2.2, the control system may but cannot be ensured to converge to the target superposition state. This problem can be solved by means of designing the control laws in the interaction picture of the control system.Consider theN-level control system in the interaction picture as
(28)i∂∂tρ(t)=[∑k=1rHk(t)uk(t),ρ(t)].Choose the Lyapunov function based on the average value of an imaginary mechanical quantity defined by (21). The first order time derivative of the Lyapunov function (21) can be obtained as
(29)V˙5=-∑k=1ruk(t)tr(iHk(t)[ρ(t),P]).
To ensure V˙5≤0, one can design uk(t) as
(30)uk(t)=-Kktr(iHk(t)([ρ(t),P])),Kk>0,
where Kk>0 and yk=fk(xk),(k=1,…,r) are monotonic increasing functions through the coordinate origin of plane xk-yk.The control system governed by (28) is a non-autonomous system; thus the LaSalle invariance principle cannot be used to analyze the convergence. One can use the improved Barbalat lemma which can be used for the non-autonomous system. According to the improved Barbalat lemma, the convergence of the control system governed by (28) can be depicted by Theorem 7.Theorem 7.
Consider the control system governed by (28) with the control field uk(t) designed in (30). If (i) the internal Hamiltonian is strongly regular; (ii) the control Hamiltonians are full connected, then any state trajectory will converge to the limit set ℛ1≡{ρs:[ρs,P]=D} at t→∞, where D is a diagonal matrix.For the case thatP is chosen as a diagonal matrix, the limit set is reduced to ℛ1≡{ρs:[ρs,P]=0}. For the case that P is chosen as a nondiagonal matrix, if rankA~(P→)=n2-n holds, then the limit set ℛ1 is regular; namely, ℛ1≡{ρs:[ρs,P]=0}, where A~(P→) is the first n2-n rows of A(P→),A(P→) is the real (n2-1)*(n2-1) matrix corresponding to the Bloch representation of Adp, and Adp is a linear map from Hermitian or anti-Hermitian matrices into su(n). [ρs,P]=0 means that P and H0 have the same eigenstates. For the target state ρf=|ψf〉〈ψf| being a superposition state, in order to make the target state contain in ℛ1≡{ρs:[ρs,P]=0}, P can be designed as
(31)P=P1|ψf〉〈ψf|+∑k=2nPk|ψk〉〈ψk|,
where |ψ1〉=|ψf〉 and 〈ψi|ψj〉=0, for i≠j.Then some deduction shows thatρs=λi|ψi〉〈ψi|(λi=1). In order to make the system converge to the target superposition state ρf from any initial pure state, one can design P such that
(32)V5(ρf)<V5(ρ0)<V5(ρoth),
where ρoth represents any other state in the set ℛ1 except the target state ρf.For the satisfaction ofV5(ρf)<V5(ρoth), one can design P based on Proposition 6. In order to make (32) holds, one must design Pi,(i=1,…,N) such that
(33)0<P1<P1tr(ρfρ0)+∑k=2NPktr(ρkρ0)<Pj(j≠1).However, for any initial pure state the target superposition state, there may not exist aP such that (33) holds.In conclusion, the design principle ofP for the case that the target state ρf=|ψf〉〈ψf| being a superposition state is (31), (33), and Proposition 6. We can also conclude that if the control system satisfies the conditions (i) and (ii) in Theorem 7 and there exists a P to make (33) hold, any state trajectory of the system will converge to the target superposition state ρf from any initial pure state ρ0.Remark 8.
One can also solve the problem of the convergence of the control system governed by (1) to the target superposition state by means of designing the control laws in the interaction picture of the control system.


### 2.4. Target Mixed State Which Does Not Commute with the Internal Hamiltonian

From Sections2.2 and 2.3, one can see that the control method proposed in Section 2.2 cannot guarantee the convergence to the target state which does not commute with the internal Hamiltonian. For the target state ρf being a mixed state which does not commute with the internal Hamiltonian, that is, [ρf,H0]≠0, Cong, Liu and Yang proposed the control system in the interaction picture and used a unitary transformation to solve the convergence of the control system to the target state. The basic idea is to use a unitary transformation ρ^f=Uf†ρfUf to make [ρ^f,H0]=0 hold. Correspondingly, the control system governed by (28) after this unitary transformation will become
(34)iℏρ^˙(t)=[∑k=1rH^k(t)uk(t),ρ^(t)],
where ρ^=Uf†ρUf, H^k(t)=Uf†Hk(t)Uf.Then control laws can be designed according to Section2.3. The research results show that the designed control laws and Theorem 7 are also valid with every physical quantity changing accordingly. In order to make ρ^f contain in ℛ2≡{ρ^s:[ρ^s,P]=0}, P needs to be designed such that [H0,P]=0 holds. [ρ^s,P]=0 and [H0,P]=0 imply that ρ^s, P, and H0 have the same eigenstates. Since the evolution of ρ^(t) is unitary, all the ρ^(t) for t≥0 are isospectral. So the eigenvalues of ρ^s are a permutation of the eigenvalues of ρ^0. Thus the limit set ℛ2≡{ρ^s:[ρ^s,P]=0} has countable elements. For sake of the convergence from any initial mixed state ρ0 to any target mixed state ρf which does not commute with the internal Hamiltonian H0, P needs to be designed such that
(35)V5(ρ^f)<V5(ρ^0)<V5(ρ^oth),
where ρ^oth represents any other state in the set ℛ2≡{ρ^s:[ρ^s,P]=0} except ρ^f.In such a way, if there exists aP to make (35) hold, any state trajectory of the system will converge to the target mixed state ρf which does not commute with the internal Hamiltonian from any initial mixed state ρ0.Next, let us analyze how to make the restriction defined by (35) hold. Some deductions show that if one designs P based on Proposition 6 with changing ρf into ρ^f accordingly, then V5(ρ^f)<V5(ρ^oth) will hold. Assume ρ0=UρfU†. In order to make the whole restriction defined by (35) hold, Pi,(i=1,…,N) needs to be designed such that
(36)(P)ll>(∑i≠ln(P)ii∑j≠kn((ρ^f)jj-(ρ^f)kk)(Uij)2+(P)ll·∑j≠k,j≠ln((ρ^f)jj-(ρ^f)kk)(Ulj)2)×((ρ^f)kk-(ρ^f)ll)-1,
where Uij is the (i,j)th element of U.In conclusion, the design principle ofP for the case that the target state ρf is a mixed state which does not commute with H0 is(i)
Pl≠Pm>0 for any l≠m∈{1,2,…,N};(ii)
Proposition6 with changing ρf into ρ^f;(iii)
Equations (36) and (16).We can conclude from the above analysis that for the case that the target state is a mixed state which does not commute withH0, if the control system is strongly regular and full connected, and one can seek an imaginary mechanical quantity P to make (36) hold, then the control system can converge from any initial mixed state to the target mixed state which does not commute with the internal Hamiltonian.


## 2.1. Target Eigenstate

It is convenient to use the bilinear Schrödinger equation to describe the control systems if the target state is a pure state. Consider theN-level closed quantum system governed by the following bilinear Schrödinger equation:
(1)i|ψ˙(t)〉=(H0+∑k=1rHkuk(t))|ψ(t)〉,
where |ψ(t)〉 is the quantum state vector, H0 is the internal Hamiltonian, Hk,(k=1,…,r) are control Hamiltonians, and uk(t),(k=1,…,r) are scalar and real control laws.In the following sections, in the case that the target state|ψf〉 is an eigenstate, that is, H0|ψf〉=λf|ψf〉, where λf is the eigenvalue of the internal Hamiltonian H0, the research results on the convergence are summarized and analyzed for the Lyapunov control based on the state distance, the state error, and the average value of an imaginary mechanical quantity, respectively.


### 2.1.1. Lyapunov Control Based on State Distance

Consider the following Lyapunov function based on the state distance:(2)V1(|ψ〉)=12(1-|〈ψ∣ψf〉|2).
The time derivative of the Lyapunov function (2) is
(3)V˙1=-∑k=1ruk(t)·Im[〈ψ∣ψf〉〈ψf|Hk|ψ〉].
The control laws uk(t) which can make V˙1≤0 hold can be designed as
(4)uk(t)=Kkfk(Im(〈ψ∣ψf〉〈ψf|Hk|ψ〉)),(k=1,…,r),
where Kk>0, and yk=fk(xk),(k=1,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.The control laws designed in (4) can only ensure the control system (1) to be stable. One needs to do further study on the convergence of the control system. The control system governed by (1) is an autonomous system, whose convergence can be analyzed based on the LaSalle invariance principle [44]. According to the LaSalle invariance principle, as t→∞, any state trajectory will converge to the largest invariant set contained in the set E in which the states satisfy that the first order derivative of the Lyapunov function equals zero. In fact, the set E contains not only the target state but also other states, thus the system may converge to other states rather than the target state. The main idea to solve this problem is to add restrictions to make the set E as small as possible. Based on the LaSalle invariance principle, the convergence of the control system governed by (1) can be depicted by Theorem 1.Theorem 1 (see [19]).
Consider the control system governed by (1) with control fields uk(t) designed in (4). If (i) The system is strongly regular, that is, ωi′j′≠ωlm,(i′,j′)≠(l,m),i′,j′,l,m∈{1,2,…,N},ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) for any |ϕi〉≠|ψf〉,i∈{1,…,N}, there exists at least a k∈{1,…,r} such that 〈ϕi|Hk|ψf〉≠0, then any state trajectory will converge toward E1={eiθ|ψf〉,θ∈R}.From Theorem1, one can see that in the case that the target state is an eigenstate, if the control system governed by (1) satisfies the conditions (i)-(ii), the control system can converge to the equivalent state of the target eigenstate eiθ|ψf〉 from any initial pure state. These two conditions are relevant to the internal Hamiltonian and the control Hamiltonians, which are system parameters. Once the control system is determined, the Hamiltonians are fixed and cannot be changed by designing the control laws.


### 2.1.2. Lyapunov Control Based on State Error

Consider the following Lyapunov function based on the state error:(5)V2(|ψ〉)=12〈ψ-ψf∣ψ-ψf〉.In the case of selecting the Lyapunov function based on the state error defined by (5), in order to facilitate to design the control laws based on the Lyapunov stability theorem, the drift item appeared in the first order time derivative of Lyapunov function, which is caused by the internal Hamiltonian, is needed to be eliminated. The existing solution is to add a global phase control item ω into the control system governed by (1). This method will not change the population distribution of the control system. Thus the dynamical equation (1) becomes
(6)i|ψ˙(t)〉=(H0+∑k=1rHkuk(t)+ωI)|ψ(t)〉.After some deduction, one can obtain the time derivative of the Lyapunov function (5) as
(7)V˙2=-(λf+ω)ℑ(〈ψf∣ψ〉)-∑k=1rℑ(〈ψf|Hk|ψ〉)uk(t).
The control laws which can make V˙2≤0 hold can be designed as
(8)ω=-λf+cf0(ℑ(〈ψf∣ψ〉)),(9)uk(t)=Kkfk(ℑ(〈ψf|Hk|ψ〉)),(k=1,…,r),
where Kk>0, and yk=fk(xk),(k=0,…,r) are the monotonic increasing functions through the coordinate origin of the plane xk-yk.Based on the LaSalle invariance principle, the convergence of the control system governed by (6) can be depicted by Theorem 2.Theorem 2 (see [19, 23]).
Consider the control system governed by (6) with control fields uk(t) designed in (9) and ω designed in (8). If (i) ωi′j′≠ωlm,(i′,j′)≠(l,m), i′,j′,l,m∈{1,2,…,N}, ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) for any |ϕi〉≠|ψf〉,i∈{1,…,N}, there exists at least a k∈{1,…,r} such that 〈ϕi|Hk|ψf〉≠0. Then any state trajectory will converge toward E2={eiθ|ψf〉,θ∈R}.From Theorem2, one can see that for the case that the target state is an eigenstate, if the control system governed by (6) satisfies the conditions (i)-(ii), the control system can also converge to the equivalent state of the target eigenstate eiθ|ψf〉 from any initial pure state.


### 2.1.3. Lyapunov Control Based on Average Value of an Imaginary Mechanical Quantity

Consider the following Lyapunov function based on the average value of an imaginary mechanical quantity:(10)V3(|ψ〉)=〈ψ|P|ψ〉,
where the imaginary mechanical quantity P is a positive definite Hermitian operator.The first order time derivative of the Lyapunov function (10) can be obtained as
(11)V˙3=i〈ψ|[H0,P]|ψ〉+i∑k=1r〈ψ|[Hk,P]|ψ〉uk.Set[H0,P]=0 such that the drift term in the right side of (11) can be eliminated. In order to ensure V˙3≤0, one can design uk(t) as
(12)uk(t)=-Kkfk(i〈ψ|[Hk,P]|ψ〉),(k=1,…,r),
where Kk>0, and yk=fk(xk),(k=1,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.Then based on the LaSalle invariance principle, all the state trajectories of the system will converge to the invariant set contained in the setE in which V˙=0 holds. Denote the state at the time t as |ψ(t)〉=∑l=1Ncl(t)|ϕl〉, where cl(t) is the coefficient corresponding to the lth eigenstate |ϕl〉. For the Schrödinger equation governed by (1), if the control system is strongly regular and any eigenstate is directly coupled to all other eigenstates, that is, for any j≠l, j,l∈{1,2,…,N}, there exists a k∈{1,…,r} such that 〈ϕj|Hk|ϕl〉≠0, then one can deduce that V˙3(t)=0 holds for all t≥t0,t0∈R is equivalent to
(13)(Pl-Pj)cj(t0)cl*(t0)=0,(l,j=1,…,N),
where Pl and Pj are the lth and jth eigenvalues of P, respectively.If the target state|ψf〉 is an eigenstate, and all the eigenvalues of P are designed mutually different, that is, Pl≠Pj for any l≠j,(l,j=1,2…,N), then (13) is equivalent to
(14)cjcl*=0,(l,j=1,…,N).
Equation (14) implies that there is at most one cj≠0,(j=1,…,N); that is, the system will converge to an eigenstate with t→∞. Thus the convergence of the control system (1) can be depicted by Theorem 3.Theorem 3 (see [19, 24, 25]).
Consider the control system governed by (1) with the control fields uk(t) designed in (12). If (i) ωi′j′≠ωlm,(i′,j′)≠(l,m),i′,j′,l,m∈{1,…,N},ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) for any i≠j,i,j∈{1,…,N}, there exsits at least a k∈{1,…,r} such that 〈ϕi|Hk|ϕj〉≠0; (iii) [H0,P]=0; (iv) for any l≠m∈{1,2,…,N}, Pl≠Pm holds, where Pl is the lth eigenvalue of P. Then any state trajectory will converge toward E3={eiθ|ϕi〉,i∈{1,…,N},θ∈R}.From Theorem3, one can see that the control system will converge from any initial pure state to an eigenstate which may not be the target eigenstate. In order to make the system converge to the target eigenstate |ψf〉 from any initial pure state |ψ0〉, as V˙3(t)≤0, one can add a restriction as
(15)V3(|ψf〉)<V3(|ψ0〉)<V3(|ψother〉),
where |ψ0〉 is the initial state, |ψother〉 represents any other state in the set E3 except the target state.In such a way, any state trajectory of the system will converge toeiθ|ψf〉 from any initial pure state |ψ0〉.Next, let us analyze how to make these convergence conditions be satisfied. Conditions (i) and (ii) are only relevant to the internal HamiltonianH0 and the control Hamiltonians Hk,(k=1,…,r) which cannot be changed by designing appropriate control laws. Condition (iii) means that P and H0 have the same eigenstates. In order to make condition (iii) be satisfied, the eigenvalues of P can be designed as
(16)P=∑j=1NPj|ϕj〉〈ϕj|.
Because P should be positive definite, one needs to design Pi>0,(i=1,…,N). The restriction (15) can be satisfied by means of designing an appropriate P. For the restriction (15), Grivopoulos and Bamieh proposed a design principle of P to make V3(|ψf〉)<V3(|ψother〉) hold. This design principle of P can be depicted by Proposition 4.Proposition 4 (see [25]).
With the constraint condition〈ψ∣ψ〉=1, the set of critical points of the Lyapunov function V3(|ψ〉)=〈ψ|P|ψ〉 is given by the normalized eigenvectors of P. The eigenvectors with the largest eigenvalue are the maxima of V3, the eigenvectors with the smallest eigenvalue are the minima and all others are saddle points.According to Proposition4, in order to make V3(|ψf〉)<V3(|ψother〉) hold, Pi>Pf,(i=1,…,N,Pi≠Pf) needs to be designed, where Pf is the eigenvalue of P corresponding to |ψf〉. Then let us consider the whole restriction (15); it is an attraction problem. If the eigenvalues of P except Pf are close together, the attraction region will be very large. For the limiting case Pi=Pj>Pf>0,(i≠j,Pi,Pj≠Pf), the attraction region will be the whole state space. Thus the design principle of P is Pi>Pf,(i=1,…,N,Pi≠Pf) and to make Pi,Pj≠Pf close together.In conclusion, for the target state being an eigenstate, the design principle ofP is as follows: (i)
Pl≠Pm>0 for any l≠m∈{1,2,…,N},(ii)
Pi>Pf,(i=1,…,N,Pi≠Pf), Pi,Pj≠Pf are close together,(iii)
Equation (16): P=∑j=1NPj|ϕj〉〈ϕj|.From the above analyses, we can conclude that the Lyapunov control method based on the imaginary mechanical quantity proposed in [19, 24, 25] can only ensure the convergence to an eigenstate, but cannot guarantee the convergence to the target eigenstate from any initial pure state. However, if there exists a P to make the restriction (15) hold, then any state trajectory of the system will converge to the equivalent state of the target eigenstate eiθ|ψf〉 from any initial pure state.


### 2.1.4. Relations between Three Lyapunov Functions

In the Liouville space, the Hilbert-Schmidt distance between two density operatorsρ1 and ρ2 is
(17)dHS(ρ1,ρ2)=tr(ρ1-ρ2)2.The inner product of two operatorsA and B is defined as 〈〈A|B〉〉=tr(A†B), where the operation A† refers to the conjugate transpose of A. Because ρ=|ψ〉〈ψ|, the square of the Hilbert-Schmidt distance between the density operator ρ and the target density operator ρf can be deduced as
(18)dHS2(ρ,ρf)=2(1-|〈ψ∣ψf〉|2).One can conclude from (18) that the Lyapunov function based on the state distance V1(|ψ〉)=(1/2)(1-|〈ψ∣ψf〉|2) and the state error V2(|ψ〉)=(1/2)〈ψ-ψf∣ψ-ψf〉 are equivalent.For the Lyapunov function based on the average value of an imaginary mechanical quantityV3(|ψ〉)=〈ψ|P|ψ〉, if the imaginary mechanical quantity P=(1/2)(I-|ψf〉〈ψf|), V3(|ψ〉) will become V1(|ψ〉) [19]. Therefore V1(|ψ〉) is formally a special case of V3(|ψ〉).In fact,V1(|ψ〉), V2(|ψ〉), and V3(|ψ〉) can be unified in the following quadratic Lyapunov function:
(19)V4(|ψ〉)=〈ψ-αψf|Q|ψ-αψf〉Next, we consider three cases as follows:(i)
α=0,Q=(1/2)(I-|ψf〉〈ψf|), V4(|ψ〉) will reduce to V1(|ψ〉);(ii)
α=1,Q=(1/2)I, V4(|ψ〉) will reduce to V2(|ψ〉);(iii)
α=0, V4(|ψ〉) will reduce to V3(|ψ〉).We can see thatV1(|ψ〉), V2(|ψ〉), and V3(|ψ〉) are special cases of V4(|ψ〉).From the above analyses, we can conclude that all these three Lyapunov control methods can converge to the equivalent state of the target eigen state from any initial pure state. The Lyapunov control methods based on the state distance and the state error have only one adjustable parameter. So these two methods are very easy to grasp and very simple. The Lyapunov control based on the average value of an imaginary mechanical quantity has more adjustable parameters. So it is more flexible and also more complex at the same time. For the target eigenstate,V1(|ψ〉) and V2(|ψ〉) are equivalent, so the Lyapunov control based on the state distance and the state error have similar control effects. Because V1(|ψ〉) is formally a special case of V3(|ψ〉), generally, the control effect of the Lyapunov control based on the average value of an imaginary mechanical quantity is better than that of the state distance. At least, it can get the same control effect as the Lyapunov control of the state distance.


## 2.1.1. Lyapunov Control Based on State Distance

Consider the following Lyapunov function based on the state distance:(2)V1(|ψ〉)=12(1-|〈ψ∣ψf〉|2).
The time derivative of the Lyapunov function (2) is
(3)V˙1=-∑k=1ruk(t)·Im[〈ψ∣ψf〉〈ψf|Hk|ψ〉].
The control laws uk(t) which can make V˙1≤0 hold can be designed as
(4)uk(t)=Kkfk(Im(〈ψ∣ψf〉〈ψf|Hk|ψ〉)),(k=1,…,r),
where Kk>0, and yk=fk(xk),(k=1,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.The control laws designed in (4) can only ensure the control system (1) to be stable. One needs to do further study on the convergence of the control system. The control system governed by (1) is an autonomous system, whose convergence can be analyzed based on the LaSalle invariance principle [44]. According to the LaSalle invariance principle, as t→∞, any state trajectory will converge to the largest invariant set contained in the set E in which the states satisfy that the first order derivative of the Lyapunov function equals zero. In fact, the set E contains not only the target state but also other states, thus the system may converge to other states rather than the target state. The main idea to solve this problem is to add restrictions to make the set E as small as possible. Based on the LaSalle invariance principle, the convergence of the control system governed by (1) can be depicted by Theorem 1.Theorem 1 (see [19]).
Consider the control system governed by (1) with control fields uk(t) designed in (4). If (i) The system is strongly regular, that is, ωi′j′≠ωlm,(i′,j′)≠(l,m),i′,j′,l,m∈{1,2,…,N},ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) for any |ϕi〉≠|ψf〉,i∈{1,…,N}, there exists at least a k∈{1,…,r} such that 〈ϕi|Hk|ψf〉≠0, then any state trajectory will converge toward E1={eiθ|ψf〉,θ∈R}.From Theorem1, one can see that in the case that the target state is an eigenstate, if the control system governed by (1) satisfies the conditions (i)-(ii), the control system can converge to the equivalent state of the target eigenstate eiθ|ψf〉 from any initial pure state. These two conditions are relevant to the internal Hamiltonian and the control Hamiltonians, which are system parameters. Once the control system is determined, the Hamiltonians are fixed and cannot be changed by designing the control laws.


## 2.1.2. Lyapunov Control Based on State Error

Consider the following Lyapunov function based on the state error:(5)V2(|ψ〉)=12〈ψ-ψf∣ψ-ψf〉.In the case of selecting the Lyapunov function based on the state error defined by (5), in order to facilitate to design the control laws based on the Lyapunov stability theorem, the drift item appeared in the first order time derivative of Lyapunov function, which is caused by the internal Hamiltonian, is needed to be eliminated. The existing solution is to add a global phase control item ω into the control system governed by (1). This method will not change the population distribution of the control system. Thus the dynamical equation (1) becomes
(6)i|ψ˙(t)〉=(H0+∑k=1rHkuk(t)+ωI)|ψ(t)〉.After some deduction, one can obtain the time derivative of the Lyapunov function (5) as
(7)V˙2=-(λf+ω)ℑ(〈ψf∣ψ〉)-∑k=1rℑ(〈ψf|Hk|ψ〉)uk(t).
The control laws which can make V˙2≤0 hold can be designed as
(8)ω=-λf+cf0(ℑ(〈ψf∣ψ〉)),(9)uk(t)=Kkfk(ℑ(〈ψf|Hk|ψ〉)),(k=1,…,r),
where Kk>0, and yk=fk(xk),(k=0,…,r) are the monotonic increasing functions through the coordinate origin of the plane xk-yk.Based on the LaSalle invariance principle, the convergence of the control system governed by (6) can be depicted by Theorem 2.Theorem 2 (see [19, 23]).
Consider the control system governed by (6) with control fields uk(t) designed in (9) and ω designed in (8). If (i) ωi′j′≠ωlm,(i′,j′)≠(l,m), i′,j′,l,m∈{1,2,…,N}, ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) for any |ϕi〉≠|ψf〉,i∈{1,…,N}, there exists at least a k∈{1,…,r} such that 〈ϕi|Hk|ψf〉≠0. Then any state trajectory will converge toward E2={eiθ|ψf〉,θ∈R}.From Theorem2, one can see that for the case that the target state is an eigenstate, if the control system governed by (6) satisfies the conditions (i)-(ii), the control system can also converge to the equivalent state of the target eigenstate eiθ|ψf〉 from any initial pure state.


## 2.1.3. Lyapunov Control Based on Average Value of an Imaginary Mechanical Quantity

Consider the following Lyapunov function based on the average value of an imaginary mechanical quantity:(10)V3(|ψ〉)=〈ψ|P|ψ〉,
where the imaginary mechanical quantity P is a positive definite Hermitian operator.The first order time derivative of the Lyapunov function (10) can be obtained as
(11)V˙3=i〈ψ|[H0,P]|ψ〉+i∑k=1r〈ψ|[Hk,P]|ψ〉uk.Set[H0,P]=0 such that the drift term in the right side of (11) can be eliminated. In order to ensure V˙3≤0, one can design uk(t) as
(12)uk(t)=-Kkfk(i〈ψ|[Hk,P]|ψ〉),(k=1,…,r),
where Kk>0, and yk=fk(xk),(k=1,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.Then based on the LaSalle invariance principle, all the state trajectories of the system will converge to the invariant set contained in the setE in which V˙=0 holds. Denote the state at the time t as |ψ(t)〉=∑l=1Ncl(t)|ϕl〉, where cl(t) is the coefficient corresponding to the lth eigenstate |ϕl〉. For the Schrödinger equation governed by (1), if the control system is strongly regular and any eigenstate is directly coupled to all other eigenstates, that is, for any j≠l, j,l∈{1,2,…,N}, there exists a k∈{1,…,r} such that 〈ϕj|Hk|ϕl〉≠0, then one can deduce that V˙3(t)=0 holds for all t≥t0,t0∈R is equivalent to
(13)(Pl-Pj)cj(t0)cl*(t0)=0,(l,j=1,…,N),
where Pl and Pj are the lth and jth eigenvalues of P, respectively.If the target state|ψf〉 is an eigenstate, and all the eigenvalues of P are designed mutually different, that is, Pl≠Pj for any l≠j,(l,j=1,2…,N), then (13) is equivalent to
(14)cjcl*=0,(l,j=1,…,N).
Equation (14) implies that there is at most one cj≠0,(j=1,…,N); that is, the system will converge to an eigenstate with t→∞. Thus the convergence of the control system (1) can be depicted by Theorem 3.Theorem 3 (see [19, 24, 25]).
Consider the control system governed by (1) with the control fields uk(t) designed in (12). If (i) ωi′j′≠ωlm,(i′,j′)≠(l,m),i′,j′,l,m∈{1,…,N},ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) for any i≠j,i,j∈{1,…,N}, there exsits at least a k∈{1,…,r} such that 〈ϕi|Hk|ϕj〉≠0; (iii) [H0,P]=0; (iv) for any l≠m∈{1,2,…,N}, Pl≠Pm holds, where Pl is the lth eigenvalue of P. Then any state trajectory will converge toward E3={eiθ|ϕi〉,i∈{1,…,N},θ∈R}.From Theorem3, one can see that the control system will converge from any initial pure state to an eigenstate which may not be the target eigenstate. In order to make the system converge to the target eigenstate |ψf〉 from any initial pure state |ψ0〉, as V˙3(t)≤0, one can add a restriction as
(15)V3(|ψf〉)<V3(|ψ0〉)<V3(|ψother〉),
where |ψ0〉 is the initial state, |ψother〉 represents any other state in the set E3 except the target state.In such a way, any state trajectory of the system will converge toeiθ|ψf〉 from any initial pure state |ψ0〉.Next, let us analyze how to make these convergence conditions be satisfied. Conditions (i) and (ii) are only relevant to the internal HamiltonianH0 and the control Hamiltonians Hk,(k=1,…,r) which cannot be changed by designing appropriate control laws. Condition (iii) means that P and H0 have the same eigenstates. In order to make condition (iii) be satisfied, the eigenvalues of P can be designed as
(16)P=∑j=1NPj|ϕj〉〈ϕj|.
Because P should be positive definite, one needs to design Pi>0,(i=1,…,N). The restriction (15) can be satisfied by means of designing an appropriate P. For the restriction (15), Grivopoulos and Bamieh proposed a design principle of P to make V3(|ψf〉)<V3(|ψother〉) hold. This design principle of P can be depicted by Proposition 4.Proposition 4 (see [25]).
With the constraint condition〈ψ∣ψ〉=1, the set of critical points of the Lyapunov function V3(|ψ〉)=〈ψ|P|ψ〉 is given by the normalized eigenvectors of P. The eigenvectors with the largest eigenvalue are the maxima of V3, the eigenvectors with the smallest eigenvalue are the minima and all others are saddle points.According to Proposition4, in order to make V3(|ψf〉)<V3(|ψother〉) hold, Pi>Pf,(i=1,…,N,Pi≠Pf) needs to be designed, where Pf is the eigenvalue of P corresponding to |ψf〉. Then let us consider the whole restriction (15); it is an attraction problem. If the eigenvalues of P except Pf are close together, the attraction region will be very large. For the limiting case Pi=Pj>Pf>0,(i≠j,Pi,Pj≠Pf), the attraction region will be the whole state space. Thus the design principle of P is Pi>Pf,(i=1,…,N,Pi≠Pf) and to make Pi,Pj≠Pf close together.In conclusion, for the target state being an eigenstate, the design principle ofP is as follows: (i)
Pl≠Pm>0 for any l≠m∈{1,2,…,N},(ii)
Pi>Pf,(i=1,…,N,Pi≠Pf), Pi,Pj≠Pf are close together,(iii)
Equation (16): P=∑j=1NPj|ϕj〉〈ϕj|.From the above analyses, we can conclude that the Lyapunov control method based on the imaginary mechanical quantity proposed in [19, 24, 25] can only ensure the convergence to an eigenstate, but cannot guarantee the convergence to the target eigenstate from any initial pure state. However, if there exists a P to make the restriction (15) hold, then any state trajectory of the system will converge to the equivalent state of the target eigenstate eiθ|ψf〉 from any initial pure state.


## 2.1.4. Relations between Three Lyapunov Functions

In the Liouville space, the Hilbert-Schmidt distance between two density operatorsρ1 and ρ2 is
(17)dHS(ρ1,ρ2)=tr(ρ1-ρ2)2.The inner product of two operatorsA and B is defined as 〈〈A|B〉〉=tr(A†B), where the operation A† refers to the conjugate transpose of A. Because ρ=|ψ〉〈ψ|, the square of the Hilbert-Schmidt distance between the density operator ρ and the target density operator ρf can be deduced as
(18)dHS2(ρ,ρf)=2(1-|〈ψ∣ψf〉|2).One can conclude from (18) that the Lyapunov function based on the state distance V1(|ψ〉)=(1/2)(1-|〈ψ∣ψf〉|2) and the state error V2(|ψ〉)=(1/2)〈ψ-ψf∣ψ-ψf〉 are equivalent.For the Lyapunov function based on the average value of an imaginary mechanical quantityV3(|ψ〉)=〈ψ|P|ψ〉, if the imaginary mechanical quantity P=(1/2)(I-|ψf〉〈ψf|), V3(|ψ〉) will become V1(|ψ〉) [19]. Therefore V1(|ψ〉) is formally a special case of V3(|ψ〉).In fact,V1(|ψ〉), V2(|ψ〉), and V3(|ψ〉) can be unified in the following quadratic Lyapunov function:
(19)V4(|ψ〉)=〈ψ-αψf|Q|ψ-αψf〉Next, we consider three cases as follows:(i)
α=0,Q=(1/2)(I-|ψf〉〈ψf|), V4(|ψ〉) will reduce to V1(|ψ〉);(ii)
α=1,Q=(1/2)I, V4(|ψ〉) will reduce to V2(|ψ〉);(iii)
α=0, V4(|ψ〉) will reduce to V3(|ψ〉).We can see thatV1(|ψ〉), V2(|ψ〉), and V3(|ψ〉) are special cases of V4(|ψ〉).From the above analyses, we can conclude that all these three Lyapunov control methods can converge to the equivalent state of the target eigen state from any initial pure state. The Lyapunov control methods based on the state distance and the state error have only one adjustable parameter. So these two methods are very easy to grasp and very simple. The Lyapunov control based on the average value of an imaginary mechanical quantity has more adjustable parameters. So it is more flexible and also more complex at the same time. For the target eigenstate,V1(|ψ〉) and V2(|ψ〉) are equivalent, so the Lyapunov control based on the state distance and the state error have similar control effects. Because V1(|ψ〉) is formally a special case of V3(|ψ〉), generally, the control effect of the Lyapunov control based on the average value of an imaginary mechanical quantity is better than that of the state distance. At least, it can get the same control effect as the Lyapunov control of the state distance.


## 2.2. Target Mixed State Which Commutes with the Internal Hamiltonian

The bilinear Schrödinger equation cannot describe the mixed state. Thus for the target mixed state, it needs to use the quantum Liouville equation which can describe the evolution of any state of a closed quantum system. Consider theN-level closed quantum system governed by the following quantum Liouville equation:
(20)iρ˙(t)=[H0+∑k=1rHkuk(t),ρ(t)],
where ρ(t) is the density operator.Consider the Lyapunov function based on the average value of an imaginary mechanical quantity:(21)V5(ρ)=tr(Pρ).
By means of setting [H0,P]=0, the first order time derivative of the Lyapunov function (21) can be deduced as
(22)V˙5=-itr([P,H0]ρ)-i∑k=1rtr([P,Hk]ρ)uk.In order to ensureV˙5≤0,  one can design uk(t) as
(23)uk(t)=iεktr([P,Hk]ρ),(k=1,…,r),
where εk∈R,εk>0.Next, let us analyze the convergence to the target state. For the control system governed by (20) in the non-degenerate case, one can deduce that V˙5(t)=0 holds for all t≥t0,t0∈R is equivalent to
(24)(Hk)jl(Pl-Pj)ρlj(t)=0,j,l=1,…,N,j<l,
where ρlj(t) is the (l;j)th element of the state ρ(t). If the target state ρf commutes with H0, that is, [H0,ρf]=0, and all the eigenvalues of P are designed mutually different, that is, Pl≠Pj for any l≠j,l,j=1,2,…,N, then (24) is equivalent to
(25)ρlj=0,(l,j=1,…,N)
which implies that the system will converge to a state which commutes with H0. Thus based on the LaSalle invariance principle, the convergence of the control system governed by (20) can be depicted by Theorem 5.Theorem 5 (see [37]).
Consider the control system governed by (20) with the control field uk(t) designed in (23). If (i) the internal Hamiltonian is strongly regular, that is, ωi′j′≠ωlm,(i′,j′)≠(l,m),i′,j′,l,m∈{1,…,N},ωlm=λl-λm, where λl is the lth eigenvalue of H0 corresponding to the eigenstate |ϕl〉; (ii) the control Hamiltonians are full coupled; that is, ∀j≠l, for k=1,…,r, there exists at least a (Hk)jl≠0; (iii) [H0,P]=0; (iv) Pl≠Pm,l≠m∈{1,2,…,N}, where Pl is the lth eigenvalue of P. Then any state trajectory will converge toward E4={ρ∣(ρ)ij=0,i,j∈{1,…,N}}, where (ρ)ij is the (l,l)th element of ρ.From Theorem5, one can see that if the control system satisfies the conditions (i)–(iv), the control system will converge from any initial state to a state that commutes with the internal Hamiltonian, which may not be the target state. Next, what we need to do is to make the control system converge to the target state.Denote the state inE4 as ρE4; then [ρE4,H0]=0 holds which implies that ρE4 and H0 have the same eigenstates. Since the evolution of ρ(t) is unitary, ρ(t) for t≥0 are isospectral. So the eigenvalues of ρE4 are a permutation of the eigenvalues of ρ0. Thus the set E4 has countable elements. If the initial state is generic, that is, the eigenvalues of the initial state are mutually different, the set E4 will have N! elements. For the target state ρf which commutes with the internal Hamiltonian, that is, [ρf,H0]=0, in order to make the system converge to the target state ρf which commutes with H0 from any initial state ρ0, Kuang and Cong proposed a restriction as
(26)V5(ρf)<V5(ρ0)<V5(ρ other),
where ρ other represents any other state in the set E4 except the target state ρf.In such a way, if there exists aP to make the restriction (26) hold, any state trajectory of the system will converge to the target state ρf which commutes with H0 from any initial state ρ0.For how to make conditions (iii)-(iv) to be satisfied, please read Section2.1.3. One can deduce the design principle of P such that V5(ρf)<V5(ρ other), the result can be depicted by Proposition 6.Proposition 6.
If(ρf)ii<(ρf)jj,1≤i,j≤N holds, then one can design Pi>Pj; if (ρf)ii=(ρf)jj,1≤i,j≤N holds, then one can design Pi≠Pj; else if (ρf)ii>(ρf)jj,1≤i,j≤N holds, then one can design Pi<Pj, thus V5(ρf)<V5(ρ other) holds, where (ρf)ii is the (i,i)th element of ρf.It is difficult to designP such that (26) holds for any initial state ρ0 and any target state ρf which satisfies [ρf,H0]=0. One possible method is to introduce a series of implicit function perturbations into the control laws, this method will be presented in Section 3.3.For the target stateρf which commutes with the internal Hamiltonian, the design principle of P is as follows: (i)
Pl≠Pm>0 for any l≠m∈{1,2,…,N};(ii)
DesignP according to Proposition 6;(iii)
Equation (16): P=∑j=1NPj|ϕj〉〈ϕj|.In conclusion, if the control system satisfies the conditions (i)–(iv) in Theorem5 and there exists a P to make the restriction (26) hold, any state trajectory of the system will converge to the target mixed state ρf which commutes with the internal Hamiltonian from any initial state ρ0.


## 2.3. Target Superposition State

From Section2.1.3, one can see that for the control system governed by (1), V˙(t)=0 holds for all t≥t0,t0∈R is equivalent to (Pl-Pj)cjcl*=0,(l,j=1,…,N). If the target state is a superposition state as |ψf〉=∑d=k1kocd|ϕd〉,(1≤k1<ko≤N), one can design Pk1=Pk2=⋯=Pko and other eigenvalues of P are mutually different, then (Pl-Pj)cjcl*=0,(l,j=1,…,N) is equivalent to
(27)cjcl*=0,(l,j=1,…,N,l,j≠k1,…,ko).
Equation (27) implies that the system will converge to the set E5=span{|ϕk1〉,|ϕk2〉,…,|ϕko〉,|ϕl〉;l∈{1,…,N},l≠k1,…,ko}, which means that the set E5 contains infinite elements. Therefore this control method cannot ensure the system converge to the target superposition state by adding the restriction defined by (15). But this control method can ensure that the system will converge to the superposition of the eigenstates corresponding to target state. In fact, once one element of P changes, all the populations of the levels will change accordingly. Therefore when there are not very many eigenstates corresponding to target state, the system maybe converge to the target state by regulating the eigenvalues of P. Otherwise, the system maybe cannot converge to the target state. Consider the extreme situation in which the target state is the superposition of all the eigenstates. According to the design principle of P, we should design P=cpI, where cp is a real number, and I is the unit matrix. Obviously, in this special case, the control method proposed in Section 2.1.3 will become invalid.From Section2.2, one can see that, for the control system governed by (20), V˙(t)=0 holds for all t≥t0,t0∈R is equivalent to (Hk)jl(Pl-Pj)ρlj(t)=0,j,l=1,…,N,j<l. Consider the target state which does not commute with H0, which includes the superposition state and the mixed state which does not commute with H0. Without loss of generality, assume (ρf)3N=(ρf)N3≠0; one can design P3=PN; thus the system will converge to E6={ρ|(ρ)ij=0,(i,j)≠(3,N),(i,j)≠(N,3),i,j∈{1,…,N}}. Because the set E6 has infinite elements, this control method also cannot ensure the system will converge to the target state by adding the restriction defined by (26). But when the target state does not have so many nonzero off-diagonal elements, the system maybe converges to the target state by regulating the diagonal elements of P. When there are many nonzero off-diagonal elements in the target state, the degree of freedom of P may be not enough.From the above analyses, we can conclude that by means of using the control methods proposed in Sections2.1 and 2.2, the control system may but cannot be ensured to converge to the target superposition state. This problem can be solved by means of designing the control laws in the interaction picture of the control system.Consider theN-level control system in the interaction picture as
(28)i∂∂tρ(t)=[∑k=1rHk(t)uk(t),ρ(t)].Choose the Lyapunov function based on the average value of an imaginary mechanical quantity defined by (21). The first order time derivative of the Lyapunov function (21) can be obtained as
(29)V˙5=-∑k=1ruk(t)tr(iHk(t)[ρ(t),P]).
To ensure V˙5≤0, one can design uk(t) as
(30)uk(t)=-Kktr(iHk(t)([ρ(t),P])),Kk>0,
where Kk>0 and yk=fk(xk),(k=1,…,r) are monotonic increasing functions through the coordinate origin of plane xk-yk.The control system governed by (28) is a non-autonomous system; thus the LaSalle invariance principle cannot be used to analyze the convergence. One can use the improved Barbalat lemma which can be used for the non-autonomous system. According to the improved Barbalat lemma, the convergence of the control system governed by (28) can be depicted by Theorem 7.Theorem 7.
Consider the control system governed by (28) with the control field uk(t) designed in (30). If (i) the internal Hamiltonian is strongly regular; (ii) the control Hamiltonians are full connected, then any state trajectory will converge to the limit set ℛ1≡{ρs:[ρs,P]=D} at t→∞, where D is a diagonal matrix.For the case thatP is chosen as a diagonal matrix, the limit set is reduced to ℛ1≡{ρs:[ρs,P]=0}. For the case that P is chosen as a nondiagonal matrix, if rankA~(P→)=n2-n holds, then the limit set ℛ1 is regular; namely, ℛ1≡{ρs:[ρs,P]=0}, where A~(P→) is the first n2-n rows of A(P→),A(P→) is the real (n2-1)*(n2-1) matrix corresponding to the Bloch representation of Adp, and Adp is a linear map from Hermitian or anti-Hermitian matrices into su(n). [ρs,P]=0 means that P and H0 have the same eigenstates. For the target state ρf=|ψf〉〈ψf| being a superposition state, in order to make the target state contain in ℛ1≡{ρs:[ρs,P]=0}, P can be designed as
(31)P=P1|ψf〉〈ψf|+∑k=2nPk|ψk〉〈ψk|,
where |ψ1〉=|ψf〉 and 〈ψi|ψj〉=0, for i≠j.Then some deduction shows thatρs=λi|ψi〉〈ψi|(λi=1). In order to make the system converge to the target superposition state ρf from any initial pure state, one can design P such that
(32)V5(ρf)<V5(ρ0)<V5(ρoth),
where ρoth represents any other state in the set ℛ1 except the target state ρf.For the satisfaction ofV5(ρf)<V5(ρoth), one can design P based on Proposition 6. In order to make (32) holds, one must design Pi,(i=1,…,N) such that
(33)0<P1<P1tr(ρfρ0)+∑k=2NPktr(ρkρ0)<Pj(j≠1).However, for any initial pure state the target superposition state, there may not exist aP such that (33) holds.In conclusion, the design principle ofP for the case that the target state ρf=|ψf〉〈ψf| being a superposition state is (31), (33), and Proposition 6. We can also conclude that if the control system satisfies the conditions (i) and (ii) in Theorem 7 and there exists a P to make (33) hold, any state trajectory of the system will converge to the target superposition state ρf from any initial pure state ρ0.Remark 8.
One can also solve the problem of the convergence of the control system governed by (1) to the target superposition state by means of designing the control laws in the interaction picture of the control system.


## 2.4. Target Mixed State Which Does Not Commute with the Internal Hamiltonian

From Sections2.2 and 2.3, one can see that the control method proposed in Section 2.2 cannot guarantee the convergence to the target state which does not commute with the internal Hamiltonian. For the target state ρf being a mixed state which does not commute with the internal Hamiltonian, that is, [ρf,H0]≠0, Cong, Liu and Yang proposed the control system in the interaction picture and used a unitary transformation to solve the convergence of the control system to the target state. The basic idea is to use a unitary transformation ρ^f=Uf†ρfUf to make [ρ^f,H0]=0 hold. Correspondingly, the control system governed by (28) after this unitary transformation will become
(34)iℏρ^˙(t)=[∑k=1rH^k(t)uk(t),ρ^(t)],
where ρ^=Uf†ρUf, H^k(t)=Uf†Hk(t)Uf.Then control laws can be designed according to Section2.3. The research results show that the designed control laws and Theorem 7 are also valid with every physical quantity changing accordingly. In order to make ρ^f contain in ℛ2≡{ρ^s:[ρ^s,P]=0}, P needs to be designed such that [H0,P]=0 holds. [ρ^s,P]=0 and [H0,P]=0 imply that ρ^s, P, and H0 have the same eigenstates. Since the evolution of ρ^(t) is unitary, all the ρ^(t) for t≥0 are isospectral. So the eigenvalues of ρ^s are a permutation of the eigenvalues of ρ^0. Thus the limit set ℛ2≡{ρ^s:[ρ^s,P]=0} has countable elements. For sake of the convergence from any initial mixed state ρ0 to any target mixed state ρf which does not commute with the internal Hamiltonian H0, P needs to be designed such that
(35)V5(ρ^f)<V5(ρ^0)<V5(ρ^oth),
where ρ^oth represents any other state in the set ℛ2≡{ρ^s:[ρ^s,P]=0} except ρ^f.In such a way, if there exists aP to make (35) hold, any state trajectory of the system will converge to the target mixed state ρf which does not commute with the internal Hamiltonian from any initial mixed state ρ0.Next, let us analyze how to make the restriction defined by (35) hold. Some deductions show that if one designs P based on Proposition 6 with changing ρf into ρ^f accordingly, then V5(ρ^f)<V5(ρ^oth) will hold. Assume ρ0=UρfU†. In order to make the whole restriction defined by (35) hold, Pi,(i=1,…,N) needs to be designed such that
(36)(P)ll>(∑i≠ln(P)ii∑j≠kn((ρ^f)jj-(ρ^f)kk)(Uij)2+(P)ll·∑j≠k,j≠ln((ρ^f)jj-(ρ^f)kk)(Ulj)2)×((ρ^f)kk-(ρ^f)ll)-1,
where Uij is the (i,j)th element of U.In conclusion, the design principle ofP for the case that the target state ρf is a mixed state which does not commute with H0 is(i)
Pl≠Pm>0 for any l≠m∈{1,2,…,N};(ii)
Proposition6 with changing ρf into ρ^f;(iii)
Equations (36) and (16).We can conclude from the above analysis that for the case that the target state is a mixed state which does not commute withH0, if the control system is strongly regular and full connected, and one can seek an imaginary mechanical quantity P to make (36) hold, then the control system can converge from any initial mixed state to the target mixed state which does not commute with the internal Hamiltonian.


## 3. Degenerate Cases

The convergence conditions (i) and (ii) proposed in the Theorem1 through Theorem 7 are relevant to the internal Hamiltonian H0 and the control Hamiltonians Hk,(k=1,…,r). They are system parameters which cannot be changed. And in practice, many actual systems do not satisfy these convergence conditions. These systems are called in the degenerate cases. In order to solve the convergence of the control systems in degenerate cases, the existing method is to introduce a series of implicit function perturbations into the control laws and choose a Lyapunov function which is an implicit function [37–40].In this section, we also divide the target state into four categories as: target eigenstate, the target mixed state which commute with the internal Hamiltonian, target superposition state, and the target mixed state which does not commute with the internal Hamiltonian. For these four target state categories, respectively, research results for the degenerate cases are summarized and analyzed.


### 3.1. Target Eigenstate

In this Section, the convergence of the control system based on the Lyapunov control method in the degenerate cases to any target eigenstate from any initial pure state will be summarized and analyzed.


#### 3.1.1. Implicit Lyapunov Control Based on State Distance

In order to solve the convergence of the control systems in the degenerate cases, several researchers introduced a series of implicit function perturbations into the control laws [37–40]. After a series of perturbations γk(t) introducing into the control laws, the dynamical equation (1) becomes
(37)i|ψ˙(t)〉=(H0+∑k=1rHk(γk(t)+vk(t)))|ψ(t)〉,
where γk(t)+vk(t)=uk(t),(k=1,…,r) are the total control laws.Without loss of generality, assume|ψf〉=|ϕg〉,1≤g≤N. In order to solve the problem of convergence for the degenerate cases, the perturbations γk(t) were introduced into the control laws. The basic idea is as follows: Denote the system with the internal Hamiltonian H0, the control Hamiltonians Hk,(k=1,…,r), and the control laws uk(t)=γk(t)+vk(t) as system 1, and the system with the internal Hamiltonian H02=H0+∑k=1rHkγk(t), the control Hamiltonians Hk,(k=1,…,r), and the control laws vk(t) as system 2. All these two systems can be depicted by (37). Denote the eigenvalues and eigenstates of H02 as λ1,γ1(t),…,γr(t),λ2,γ1(t),…,γr(t),…,λN,γ1(t),…,γr(t) and |ϕ1,γ1(t),…,γr(t)〉,|ϕ2,γ1(t),…,γr(t)〉,…,|ϕN,γ1(t),…,γr(t)〉, respectively, which are the functions of the perturbations γk(t). Assume |ψf,γ1(t),…,γr(t)〉=|ϕg,γ1(t),…,γr(t)〉. If one can design the perturbations γk(t) such that (i) ωl,m,γ1,…,γr≠ωi,j,γ1,…,γr,(l,m)≠(i,j),i,j,l,m∈{1,2,…,N}, where ωl,m,γ1,…,γr=λl,γ1,…,γr-λm,γ1,…,γr holds; (ii) for any |ϕγ1,…,γr〉≠|ψf,γ1,…,γr〉, there exists at least a k∈{1,…,r} satisfying 〈ϕγ1,…,γr|Hk|ψf,γ1,…,γr〉≠0, and select the specific Lyapunov function based on the state distance as
(38)V6(t)=12(1-|〈ψ∣ψf,γ1(t),…,γr(t)〉|2).
Then according to Section 2.1.1, system 2 will converge to |ψf,γ1(t),…,γr(t)〉. And when system 2 converge to |ψf,γ1(t),…,γr(t)〉, if the perturbations γk(t) at |ψf,γ1(t),…,γr(t)〉 are designed to equal zero, then system 2 will become system 1, and |ψf,γ1(t),…,γr(t)〉 will become |ψf〉. Then the convergence of system 1 to |ψf〉 will be ensured. In fact, the evolution of system 1 can be viewed as a composite of two evolution processes. One is system 2 converges to |ψf,γ1(t),…,γr(t)〉 from the initial state |ψ0〉, another one is γk(t) converges to 0. In order to make the introduced perturbations take effect to make system 1 in the non-degenerate case converge to |ψf〉, the speed of γk(t),(k=1,…,r) converging to 0 must be slower than the speed at which system 2 converges toward |ψf,γ1(t),…,γr(t)〉. For convenience, the control system in the following section means system 1.The existing design method ofγk(t) is to design it to be a monotonically increasing functional of V(t), that is,
(39)γk(t)=θk(V(t))=θk(12(1-|〈ψ∣ψf,γ1(t),…,γr(t)〉|2)),(k=1,…,r),
where the function θk(·) satisfies θk(0)=0, θk(s)>0 and θk′(s)>0 for every s>0 and s is the independent variable of the function θk(·).The right side of (39) contains |ψf,γ1(t),…,γr(t)〉 which is a functional of the perturbations γk(t),(k=1,…,r). One can see that the relation between γk and the time t is defined by r equations and cannot be expressed by an explicit expression, so γk is the implicit function of the time t. From (38) and (39), one can also see that the Lyapunov function V6 is the implicit function of the time t. The existence of γk(t) can be depicted by Lemma 9.Lemma 9 (see [37, 40]).
Letθk∈C∞(R+;[0,γk*]),k=1,…,r,(γk*∈R,γk*>0) satisfy θk(0)=0, θk(s)>0, θk′(s)>0 for every s>0, ∥θk′∥∞<1/C*, and C*=1+max{∥(∂ψf,γ1(t),…,γr(t)/∂γk)|(γ10,…,γr0)∥;γk0∈[0,γk*]}, where γk0=γk(0). Then for every state |ψ〉∈S2N-1={x∈CN;∥x∥=1}, there exists a unique γ1,γ2,…, and γr with γk∈C∞(γk∈[0,γk*]),(k=1,…,r) satisfying γk(t)=θk((1/2)(1-|〈ψ∣ψf,γ1(t),…,γr(t)〉|2)),(k=1,…,r),γk(|ψf〉)=0.By (38) and Lemma 9, one can see that for every state |ψ〉∈S2N-1={x∈CN;∥x∥=1}; there also exists a unique V.Assume(40)∂|ψf,γ1(t),…,γr(t)〉∂γ1=∂|ψf,γ1(t),…,γr(t)〉∂γ2=⋯=∂|ψf,γ1(t),…,γr(t)〉∂γr.
The first order time derivative of the Lyapunov function (38) is
(41)V˙6=-∑k=1rvk(t)ℑ(〈ψ∣ψf,γ1,…,γr〉〈ψf,γ1,…,γr|Hk|ψ〉)(1+ℜ(〈ψ∣ψf,γ1,…,γr〉〈(∂|ψf,γ1,…,γr〉/∂γk)∣ψ〉)∑j=1rθj′).
In order to make V˙6(t)≤0, one can design γk(t) such that ∥θj′∥∞<1/(rC*) to make 1/(1+ℜ(〈ψ∣ϕγ1,…,γr〉〈(∂ϕγ1,…,γr/∂γk)∣ψ〉)∑j=1rθj′)>0 hold and design vk(t) as
(42)vk(t)=Kkfk(ℑ(ei∠〈ψ∣ψf,γ1,…,γr〉〈ψf,γ1,…,γr|Hk|ψ〉)),
where Kk>0, and vk(t),(k=1,…,r) are monotonic functions through the coordinate origin and in the first quadrant and the third quadrant.Based on the LaSalle invariance principle, the convergence of the control system governed by (37) can be depicted by Theorem 10.Theorem 10 (see [37, 40]).
Consider the control system (37) with control fields uk(t)=γk(t)+vk(t),γk(t) designed in (39), Lemma 9 and ∥θj′∥∞<1/(rC*), and vk(t) designed in (42). Assume that the target state |ψf〉 is an eigenstate of H0. If the control system satisfies (i) ωl,m,γ1,…,γr≠ωi,j,γ1,…,γr,(l,m)≠(i,j), i,j,l,m∈{1,2,…,N}, ωl,m,γ1,…,γr=λl,γ1,…,γr-λm,γ1,…,γr, where λl,γ1,…,γr is an eigenvalue of (H0+∑k=1rHkγk), (ii) For any |ϕγ1,…,γr〉≠|ψf,γ1,…,γr〉, there exists at least a k∈{1,…,r} satisfying 〈ϕγ1,…,γr|Hk|ψf,γ1,…,γr〉≠0; then the largest invariant set is S2N-1∩E7, where E7={|ψ〉=eiθ|ψf〉,θ∈R}. And the control system will converge toward |ψf〉eiθ,(θ∈R).Conditions (i) and (ii) in Theorem10 are associated with H0, Hk,(k=1,…,r) and γk(t). By designing appropriate γk(t), these two conditions can be satisfied in most cases. one can see that if one designs appropriate control laws uk(t)=γk(t)+vk(t),(k=1,…,r) to make the conditions (i) and (ii) hold, the control system depicted by (37) can converge to the equivalent state of the target eigenstate |ψf〉 from any initial pure state.


#### 3.1.2. Implicit Lyapunov Control Based on the State Error

The basic idea of this method is similar to that of the state distance. In the Lyapunov control based on the state error, a global phase control itemω is added into the control system to facilitate the design of the control laws. Thus the control system can be depicted by
(43)i|ψ˙(t)〉=(H0+∑k=1rHk(γk(t)+vk(t))+ωI)|ψ(t)〉,
where uk(t)=γk(t)+vk(t),(k=1,…,r) and ω is the control laws which need to design.Consider the Lyapunov function based on the state error as(44)V7(t)=12〈ψ-ψf,γ1,…,γr∣ψ-ψf,γ1,…,γr〉.
Thus the implicit function perturbation γk(t) needs to be designed as
(45)γk(t)=θk(V(t))=θk(12〈ψ-ψf,γ1,…,γr∣ψ-ψf,γ1,…,γr〉),
where functions θk(·) satisfy θk(0)=0, θk(s)>0, θk′(s)>0 for every s>0.The existence ofγk(t) can be depicted by Lemma 11.Lemma 11 (see [38, 39]).
Letθk∈C∞(R+;[0,γk*]),k=1,…,r,(γk*∈R,γk*>0) satisfy θk(0)=0,θk(s)>0,θk′(s)>0 for every s>0,∥θk′∥∞<1/rC*, and C*=1+max{∥(∂|ψf,γ1,…,γr〉/∂γk)|(γ10,…,γr0)∥;γk0∈[0,γk*]}. Then for every state |ψ〉∈S2N-1={x∈CN;∥x∥=1}, there exists a unique γ1,γ2,…, and γr with γk∈C∞(γk∈[0,γk*]),(k=1,…,r) satisfying γk(|ψ〉)=θk((1/2)〈ψ-ψf,γ1(|ψ〉),…,γr(|ψ〉)∣ψ-ψf,γ1(|ψ〉),…,γr(|ψ〉)〉), γk(|ψf〉)=0.Assume (40) holds, the first order time derivative of the Lyapunov function (44) is
(46)V˙7=-(1(1+ℜ(〈(∂|ψf,γ1,…,γr〉/∂γk)∣ψ〉)∑j=1rθj′))·(∑k=1rℑ(〈ψf,γ1,…,γr|Hk|ψ〉)vk(t)(λγ1,…,γr+ω)ℑ(ψf,γ1,…,γr∣ψ)+∑k=1rℑ(〈ψf,γ1,…,γr|Hk|ψ〉)vk(t)).
In order to make V˙7(t)≤0, ω and vk(t) need to be designed as
(47)ω=-λγ1,…,γr+cf0(ℑ(〈ψf,γ1,…,γr∣ψ〉)),(48)vk(t)=Kkfk(ℑ(〈ψf,γ1,…,γr|Hk|ψ〉)),(k=1,…,r),
where Kk>0, c>0, and yk=fk(xk),(k=0,1,…,r) are monotonic functions through the coordinate origin and in the first quadrant and the third quadrant.Based on the LaSalle invariance principle, the convergence of the control system governed by (44) can be depicted by Theorem 12.Theorem 12 (see [38, 39]).
Consider the control system (43) with control fields γk(t) designed in (45) and Lemma 11, vk(t) designed in (48), and ω designed in (47). Assume that the target state |ψf〉 is an eigenstate of H0. If the control system satisfies (i) ωl,m,γ1,…,γr≠ωi,j,γ1,…,γr,(l,m)≠(i,j), i,j,l,m∈{1,…,N}, ωl,m,γ1,…,γr=λl,γ1,…,γr-λm,γ1,…,γr,λl,γ1,…,γr is the lth eigenvalue of (H0+∑k=1rHk(γk+ηk)) corresponding to the eigenstate |ϕl,γ1,…,γr〉; (ii) There is at least a k∈{1,…,r}, 〈ϕγ1(t),…,γr(t)|Hk|ψf,γ1(t),…,γr(t)〉≠0, for |ϕγ1(t),…,γr(t)〉≠|ψf,γ1(t),…,γr(t)〉; then the largest invariant set is S2N-1∩E8, where E8={|ψ〉=eiθ|ψf〉,θ∈R}. And the control system will converge toward the equivalent state of the target state: |ψf〉eiθ,(θ∈R).By designing appropriate perturbationsγk(t), conditions (i) and (ii) in Theorem 12 can be satisfied in most cases. From Theorem 12, one can see that if one design appropriate control laws uk(t)=γk(t)+vk(t),(k=1,…,r) and ω to make the conditions (i) and (ii) in Theorem 12 hold, the control system depicted by (43) can converge to the equivalent state of the target eigenstate |ψf〉 from any initial pure state.


#### 3.1.3. Implicit Lyapunov Control Based on Average Value of an Imaginary Mechanical Quantity

For the implicit Lyapunov control based on the average value of an imaginary mechanical quantity, one can consider the control system depicted by (37). In this control method, the introduced implicit function perturbations γk(t) mainly have two tasks. One task is to solve the convergent problem of the control system in the degenerate cases. The basic idea is similar to that of the state distance. Another one is to choose a simpler restriction V(|ψf〉)<V(|ψother〉), which can be satisfied for any initial state and any target state by designing the imaginary mechanical quantity. In order to ensure the system converge to the target state by adding V(|ψf〉)<V(|ψother〉), we can design all the perturbations γk(t)=0 hold for k=1,…,r only at |ψf〉, that is, (1) γk(|ψf〉)=0,(k=1,…,r), and (2) for |ψ〉≠|ψf〉, there exists at least one k such that γk(|ψ〉)≠0. For sake of completing these two tasks, we can design γk(t) as a monotonically increasing functional of V(t):
(49)γk(|ψ〉)=Ck·θk(V(|ψ〉)-V(|ψf〉)),
where Ck≥0, and for k=1,…,r, there exists at least a Ck>0, and the function θk(·) satisfies θk(0)=0, θk(s)>0 and θk′(s)>0 for every s>0, s is the independent variable of the function θk(·).The specific Lyapunov function based on the average value of an imaginary mechanical quantity can be selected as:(50)V8(|ψ〉)=〈ψ|Pγ1,…,γr|ψ〉,
where Pγ1,…,γr=f(γ1(t),…,γr(t)) is a functional of γk(t) and positive definite.The existence ofγk(t) can be depicted by Lemma 13.Lemma 13.
IfCk=0, γk(|ψ〉)=0. Else if Ck>0,θk∈C∞(R+;[0,γk*]),k=1,…,r(γk* is a positive constant) satisfy θk(0)=0,θk(s)>0 and θk′(s)>0 for every s>0, and |θk′|<1/(2CkC*), C*=1+C, C=max{∥∂Pγ1,…,γr/∂γk∥∞,γk∈[0,γk*]}, then for every |ψ〉∈S2N-1, there is a unique γk∈C∞(γk∈[0,γk*]) satisfying γk(|ψ〉)=Ck·θk(〈ψ|Pγ1,…,γr|ψ〉-〈ψf|Pγ1,…,γr|ψf〉).
For the sake of simplicity,γk(t) can be designed zero for some k, and other γk(t) equal, denoted by γ(t); that is, set
(51)γk(t)=γ(t)=θ(〈ψ|Pγ|ψ〉-〈ψf|Pγ|ψf〉),k=k1,…,km;γk(t)=0,k≠k1,…,km(1≤k1,…,km≤r),
where θ(·)=θk1(·)=⋯=θkm(·) and Pγ are functionals of γ(t).
The time derivative of the selected Lyapunov function is(52)V˙8=∑k=1rivk(t)〈ψ|[Hk,Pγ]|ψ〉·(1+θ′(〈ψf|(∂Pγ/∂γ)|ψf〉))(1-θ′(〈ψ|(∂Pγ/∂γ)|ψ〉-〈ψf|(∂Pγ/∂γ)|ψf〉)).
In order to ensure V˙8(t)≤0,vk(t) can be designed as
(53)vk(t)=-Kkfk(i〈ψ|[Hk,Pγ]|ψ〉),(k=1,…,r),
where Kk is a constant and Kk>0, and yk=fk(xk),(k=1,2,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.Based on LaSalle’s invariance principle, the convergence of the control system governed by (37) can be depicted by Theorem 14.Theorem 14.
Consider the control system governed by (37) with control fields uk(t)=γk(t)+vk(t),(k=1,…,r), where γk(t) is defined by Lemma 13 and (51) and vk(t) is defined by (53). If the control system satisfies (i) ωl,m,γ≠ωi,j,γ,(l,m)≠(i,j), i,j,l,m∈{1,2,…,N}, ωl,m,γ=λl,γ-λm,γ, where λl,γ is the lth eigenvalue of H0+∑n=k1kmHnγ(t) corresponding to the eigenstate |ϕl,γ〉; (ii) for any i≠j,i,j∈{1,2,…,N}, there exists at least one k such that (H^k)lm≠0, where (H^k)lm is the (l,m)th element of H^k=U1†HkU1, U1=(|ϕ1,γ〉,…,|ϕN,γ〉); (iii) [Pγ,H0+∑n=k1kmHnγ(t)]=0; (iv) (P^γ)ll≠(P^γ)mm,l≠m, where (P^γ)ll is the (l,l)th element of U1HPγU1, then any trajectory will converge toward E9={|ψt0〉|eiθl|ϕl,γ(|ψt0〉)〉;θl∈R,l∈{1,…,N}}.From Theorem14, one can see that if the target state |ψf〉 is an eigenstate, |ψf〉 is contained in E9. In order to make the system converge to the target eigenstate from any initial pure state |ψ0〉, on the one hand, Pγ needs to be designed to make
(54)V8(|ψf〉)<V8(|ψother〉)
hold, where |ψother〉 represents any other state in the set E9 except the target state. On the other hand, because ∂γ/∂V>0,V˙≤0,γ≥0 holds, when vk(t)=0,γ(t)=γ-≠0 holds for some time, we can design γ=γ--α,(0<α≪γ-) to make the state trajectory evolve but not stay in E9 until |ψf〉eiθl is reached.Next, let us analyze how to realize the convergence conditions in Theorem14 and restriction defined by (54). By designing appropriate γk(t), conditions (i) and (ii) in Theorem 14 can be satisfied in most cases. In order to make condition (iii) be satisfied, the eigenvalues of Pγ can be designed constant, denoted by P1,P2,…,PN, and Pγ can be designed as
(55)Pγ=∑j=1NPj|ϕj,γ〉〈ϕj,γ|.To make condition (iv) satisfied,Pl≠Pj(∀l≠j;1≤l,j≤N) can be designed. The results on how to make (54) hold can be depicted by Theorem 15.Theorem 15.
If one designsPi>Pf,(i=1,…,N,Pi≠Pf), then V8(|ψf〉)<V8(|ψother〉) holds, where Pf is the eigenvalue of Pγ(|ψf〉) corresponding to |ψf〉.We can conclude from the above analyses that by using the implicit Lyapunov control based on the imaginary mechanical quantity, if one designs appropriate control lawsuk(t)=γk(t)+vk(t),(k=1,…,r) to make the conditions (i)–(iv) in Theorem 15 and (54) hold, the control system governed by (37) in the degenerate cases can converge from any initial pure state to the target eigenstate. The design principle of the imaginary mechanical quantity Pγ is Theorem 15 and (55).


### 3.2. Target Superposition State

The method proposed in Section3.1 cannot guarantee the control system governed by (37) or (43) in the degenerate cases converge from any initial pure state to the target superposition state. In order to solve this problem, one can introduce a series of constant disturbances ηk into the control laws. Thus the mechanical equation (37) becomes
(56)i|ψ˙(t)〉=(H0+∑k=1rHk(ηk+γk(t)+vk(t)))|ψ(t)〉.
And the mechanical equation (43) in Section 3.1.2 becomes
(57)i|ψ˙(t)〉=(H0+∑k=1rHk(γk(t)+vk(t)+ηk)+ωI)|ψ(t)〉,
where ηk∈R.The basic idea of solving the convergence to the target superposition state is to designηk to make the target state |ψf〉 be an eigenstate of H0′=H0+∑k=1rHkηk·H0′ can be viewed as the new internal Hamiltonian of the control system. If the number of the control Hamiltonians r is large enough, by designing appropriate ηk, (H0+∑k=1rHkηk)|ψf〉=λf′|ψf〉 can be satisfied in most cases, where λf′ is the eigenvalue of H0+∑k=1rHkηk corresponding to |ψf〉. Then one can design the control laws and analyze the convergence according to the method for the target eigen state cases. Research results show that every conclusion in Section 3.1 also holds with changing H0 into H0′.


### 3.3. Target Mixed State Which Commutes with the Internal Hamiltonian

Consider theN-level closed quantum system governed by the following quantum Liouville equation:
(58)iρ˙(t)=[H0+∑k=1rHk(γk(t)+vk(t)),ρ(t)],
where γk(t)+vk(t)=uk(t) are the total control laws.The design ideas are similar to those of Section3.1. Consider the Lyapunov function based on the imaginary mechanical quantity as
(59)V9(ρ)=tr(Pγ1,…,γrρ),γk(t) can be designed as
(60)γk(t)=γ(t)=θ(V9(ρ)-V9(ρf)),k=k1,…,km;γk(t)=0,k≠k1,…,km(1≤k1,…,km≤r).The existence ofγ(t) can be depicted by Lemma 16.Lemma 16.
Ifθ∈C∞(R+;[0,γ*]),k=1,…,r (γ* is a positive constant) satisfy θ(0)=0,θ(s)>0 and θ′(s)>0 for every s>0, and |θ′|<1/(2C*), C*=1+C, C=max{∥∂Pγ/∂γ∥m1,γ∈[0,γ*]}, then for every ρ, there is a unique γ∈C∞(γ∈[0,γ*]) satisfying γ(ρ)=θ(tr(Pγρ)-tr(Pγρf)).According to the Lyapunov stabilty Theorem,vk(t) can be designed as
(61)vk(t)=Kkfk(itr([Pγ,Hk]ρ)),(k=1,…,r),
where Kk is a constant and Kk>0, and yk=fk(xk),(k=1,2,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.Based on the LaSalle invariance principle, the convergence of the control system can be depicted by Theorem17.Theorem 17.
Consider the control system depicted by (58) with control laws uk(t)=γk(t)+vk(t), where γk(t) is defined by Lemma 16 and (60), and vk(t) is defined by (61). If the control system satisfies (i) ωl,m,γ≠ωi,j,γ,(l,m)≠(i,j),i,j,l,m∈{1,2,…,N}, ωl,m,γ=λl,γ-λm,γ, where λl,γ is the lth eigenvalue of (H0+∑n=k1kmHnγ(t)) corresponding to the eigenstate |ϕl,γ〉; (ii) ∀j≠l, for k=1,…,r, there exists at least a (H^k)jl≠0, where (H^k)jl is the (j,l)th element of H^k=U2HHkU2 with U2=(|ϕ1,γ〉,…,|ϕN,γ〉); (iii) [Pγ,H0+∑n=k1kmHnγ(t)]=0,1≤k1,…,km≤r; (iv) for any l≠j,(1≤l,j≤N),(P^γ)ll≠(P^γ)jj holds, where (P^γ)ll is the (l,l)th element of P^γ=U2HPγU2, then the control system will converge toward E10={ρt0|(U2†ρt0U2)ij=0,γ=γ(ρt0),t0∈R}.For the target stateρf which commutes with the internal Hamiltonian; that is, [ρf,H0]=0, ρf is contained in E10. In order to make the system converge to the target state ρf which commutes with H0 from any initial stat; on the one hand, Pγ needs to be designed to make
(62)V9(ρf)<V9(ρother)
hold, where ρother represents any other state in the set E10 except the target state. On the other hand, γ was designed as γ=γ--α,(0<α≪γ-) when vk(t)=0,γ(t)=γ-≠0 holds for some time to make the state trajectory evolve but not stay in E10 until ρf is reached.For the satisfaction of conditions (i)–(iv) in Theorem17, one can follow that of Section 3.1.3. In order to make (62) hold, one can design Pγ based on Proposition 6. The total design principle of Pγ is (55) and Proposition 6.We can conclude from the above analyses that for the target state which commutes with the internal Hamiltonian, by using the implicit Lyapunov control method based on the imaginary mechanical quantity, if one designs appropriate control lawsuk(t)=γk(t)+vk(t),(k=1,…,r) to make the conditions (i)–(iv) in Theorem 17 and (62) hold, the control system governed by (58) in the degenerate cases can converge from any initial state to the target state which commutes with the internal Hamiltonian, which contains the target eigenstate and the target mixed state which commutes with the internal Hamiltonian.


### 3.4. Target Mixed State Which Does Not Commute with Internal Hamiltonian

For the target mixed state which does not commute with the internal Hamiltonian, the design idea is similar to Section3.2. The difference is to design ηk to make [ρf,H0′]=0. If the number of the control Haimltonians r is large enough, by designing appropriate ηk, [ρf,H0′]=0 can be satisfied in most cases. Then one can design the control laws and analyze the convergence according to the method mentioned in Section 3.3. Research results show that every conclusion in Section 3.3 also holds with changing H0 into H0′.


## 3.1. Target Eigenstate

In this Section, the convergence of the control system based on the Lyapunov control method in the degenerate cases to any target eigenstate from any initial pure state will be summarized and analyzed.


### 3.1.1. Implicit Lyapunov Control Based on State Distance

In order to solve the convergence of the control systems in the degenerate cases, several researchers introduced a series of implicit function perturbations into the control laws [37–40]. After a series of perturbations γk(t) introducing into the control laws, the dynamical equation (1) becomes
(37)i|ψ˙(t)〉=(H0+∑k=1rHk(γk(t)+vk(t)))|ψ(t)〉,
where γk(t)+vk(t)=uk(t),(k=1,…,r) are the total control laws.Without loss of generality, assume|ψf〉=|ϕg〉,1≤g≤N. In order to solve the problem of convergence for the degenerate cases, the perturbations γk(t) were introduced into the control laws. The basic idea is as follows: Denote the system with the internal Hamiltonian H0, the control Hamiltonians Hk,(k=1,…,r), and the control laws uk(t)=γk(t)+vk(t) as system 1, and the system with the internal Hamiltonian H02=H0+∑k=1rHkγk(t), the control Hamiltonians Hk,(k=1,…,r), and the control laws vk(t) as system 2. All these two systems can be depicted by (37). Denote the eigenvalues and eigenstates of H02 as λ1,γ1(t),…,γr(t),λ2,γ1(t),…,γr(t),…,λN,γ1(t),…,γr(t) and |ϕ1,γ1(t),…,γr(t)〉,|ϕ2,γ1(t),…,γr(t)〉,…,|ϕN,γ1(t),…,γr(t)〉, respectively, which are the functions of the perturbations γk(t). Assume |ψf,γ1(t),…,γr(t)〉=|ϕg,γ1(t),…,γr(t)〉. If one can design the perturbations γk(t) such that (i) ωl,m,γ1,…,γr≠ωi,j,γ1,…,γr,(l,m)≠(i,j),i,j,l,m∈{1,2,…,N}, where ωl,m,γ1,…,γr=λl,γ1,…,γr-λm,γ1,…,γr holds; (ii) for any |ϕγ1,…,γr〉≠|ψf,γ1,…,γr〉, there exists at least a k∈{1,…,r} satisfying 〈ϕγ1,…,γr|Hk|ψf,γ1,…,γr〉≠0, and select the specific Lyapunov function based on the state distance as
(38)V6(t)=12(1-|〈ψ∣ψf,γ1(t),…,γr(t)〉|2).
Then according to Section 2.1.1, system 2 will converge to |ψf,γ1(t),…,γr(t)〉. And when system 2 converge to |ψf,γ1(t),…,γr(t)〉, if the perturbations γk(t) at |ψf,γ1(t),…,γr(t)〉 are designed to equal zero, then system 2 will become system 1, and |ψf,γ1(t),…,γr(t)〉 will become |ψf〉. Then the convergence of system 1 to |ψf〉 will be ensured. In fact, the evolution of system 1 can be viewed as a composite of two evolution processes. One is system 2 converges to |ψf,γ1(t),…,γr(t)〉 from the initial state |ψ0〉, another one is γk(t) converges to 0. In order to make the introduced perturbations take effect to make system 1 in the non-degenerate case converge to |ψf〉, the speed of γk(t),(k=1,…,r) converging to 0 must be slower than the speed at which system 2 converges toward |ψf,γ1(t),…,γr(t)〉. For convenience, the control system in the following section means system 1.The existing design method ofγk(t) is to design it to be a monotonically increasing functional of V(t), that is,
(39)γk(t)=θk(V(t))=θk(12(1-|〈ψ∣ψf,γ1(t),…,γr(t)〉|2)),(k=1,…,r),
where the function θk(·) satisfies θk(0)=0, θk(s)>0 and θk′(s)>0 for every s>0 and s is the independent variable of the function θk(·).The right side of (39) contains |ψf,γ1(t),…,γr(t)〉 which is a functional of the perturbations γk(t),(k=1,…,r). One can see that the relation between γk and the time t is defined by r equations and cannot be expressed by an explicit expression, so γk is the implicit function of the time t. From (38) and (39), one can also see that the Lyapunov function V6 is the implicit function of the time t. The existence of γk(t) can be depicted by Lemma 9.Lemma 9 (see [37, 40]).
Letθk∈C∞(R+;[0,γk*]),k=1,…,r,(γk*∈R,γk*>0) satisfy θk(0)=0, θk(s)>0, θk′(s)>0 for every s>0, ∥θk′∥∞<1/C*, and C*=1+max{∥(∂ψf,γ1(t),…,γr(t)/∂γk)|(γ10,…,γr0)∥;γk0∈[0,γk*]}, where γk0=γk(0). Then for every state |ψ〉∈S2N-1={x∈CN;∥x∥=1}, there exists a unique γ1,γ2,…, and γr with γk∈C∞(γk∈[0,γk*]),(k=1,…,r) satisfying γk(t)=θk((1/2)(1-|〈ψ∣ψf,γ1(t),…,γr(t)〉|2)),(k=1,…,r),γk(|ψf〉)=0.By (38) and Lemma 9, one can see that for every state |ψ〉∈S2N-1={x∈CN;∥x∥=1}; there also exists a unique V.Assume(40)∂|ψf,γ1(t),…,γr(t)〉∂γ1=∂|ψf,γ1(t),…,γr(t)〉∂γ2=⋯=∂|ψf,γ1(t),…,γr(t)〉∂γr.
The first order time derivative of the Lyapunov function (38) is
(41)V˙6=-∑k=1rvk(t)ℑ(〈ψ∣ψf,γ1,…,γr〉〈ψf,γ1,…,γr|Hk|ψ〉)(1+ℜ(〈ψ∣ψf,γ1,…,γr〉〈(∂|ψf,γ1,…,γr〉/∂γk)∣ψ〉)∑j=1rθj′).
In order to make V˙6(t)≤0, one can design γk(t) such that ∥θj′∥∞<1/(rC*) to make 1/(1+ℜ(〈ψ∣ϕγ1,…,γr〉〈(∂ϕγ1,…,γr/∂γk)∣ψ〉)∑j=1rθj′)>0 hold and design vk(t) as
(42)vk(t)=Kkfk(ℑ(ei∠〈ψ∣ψf,γ1,…,γr〉〈ψf,γ1,…,γr|Hk|ψ〉)),
where Kk>0, and vk(t),(k=1,…,r) are monotonic functions through the coordinate origin and in the first quadrant and the third quadrant.Based on the LaSalle invariance principle, the convergence of the control system governed by (37) can be depicted by Theorem 10.Theorem 10 (see [37, 40]).
Consider the control system (37) with control fields uk(t)=γk(t)+vk(t),γk(t) designed in (39), Lemma 9 and ∥θj′∥∞<1/(rC*), and vk(t) designed in (42). Assume that the target state |ψf〉 is an eigenstate of H0. If the control system satisfies (i) ωl,m,γ1,…,γr≠ωi,j,γ1,…,γr,(l,m)≠(i,j), i,j,l,m∈{1,2,…,N}, ωl,m,γ1,…,γr=λl,γ1,…,γr-λm,γ1,…,γr, where λl,γ1,…,γr is an eigenvalue of (H0+∑k=1rHkγk), (ii) For any |ϕγ1,…,γr〉≠|ψf,γ1,…,γr〉, there exists at least a k∈{1,…,r} satisfying 〈ϕγ1,…,γr|Hk|ψf,γ1,…,γr〉≠0; then the largest invariant set is S2N-1∩E7, where E7={|ψ〉=eiθ|ψf〉,θ∈R}. And the control system will converge toward |ψf〉eiθ,(θ∈R).Conditions (i) and (ii) in Theorem10 are associated with H0, Hk,(k=1,…,r) and γk(t). By designing appropriate γk(t), these two conditions can be satisfied in most cases. one can see that if one designs appropriate control laws uk(t)=γk(t)+vk(t),(k=1,…,r) to make the conditions (i) and (ii) hold, the control system depicted by (37) can converge to the equivalent state of the target eigenstate |ψf〉 from any initial pure state.


### 3.1.2. Implicit Lyapunov Control Based on the State Error

The basic idea of this method is similar to that of the state distance. In the Lyapunov control based on the state error, a global phase control itemω is added into the control system to facilitate the design of the control laws. Thus the control system can be depicted by
(43)i|ψ˙(t)〉=(H0+∑k=1rHk(γk(t)+vk(t))+ωI)|ψ(t)〉,
where uk(t)=γk(t)+vk(t),(k=1,…,r) and ω is the control laws which need to design.Consider the Lyapunov function based on the state error as(44)V7(t)=12〈ψ-ψf,γ1,…,γr∣ψ-ψf,γ1,…,γr〉.
Thus the implicit function perturbation γk(t) needs to be designed as
(45)γk(t)=θk(V(t))=θk(12〈ψ-ψf,γ1,…,γr∣ψ-ψf,γ1,…,γr〉),
where functions θk(·) satisfy θk(0)=0, θk(s)>0, θk′(s)>0 for every s>0.The existence ofγk(t) can be depicted by Lemma 11.Lemma 11 (see [38, 39]).
Letθk∈C∞(R+;[0,γk*]),k=1,…,r,(γk*∈R,γk*>0) satisfy θk(0)=0,θk(s)>0,θk′(s)>0 for every s>0,∥θk′∥∞<1/rC*, and C*=1+max{∥(∂|ψf,γ1,…,γr〉/∂γk)|(γ10,…,γr0)∥;γk0∈[0,γk*]}. Then for every state |ψ〉∈S2N-1={x∈CN;∥x∥=1}, there exists a unique γ1,γ2,…, and γr with γk∈C∞(γk∈[0,γk*]),(k=1,…,r) satisfying γk(|ψ〉)=θk((1/2)〈ψ-ψf,γ1(|ψ〉),…,γr(|ψ〉)∣ψ-ψf,γ1(|ψ〉),…,γr(|ψ〉)〉), γk(|ψf〉)=0.Assume (40) holds, the first order time derivative of the Lyapunov function (44) is
(46)V˙7=-(1(1+ℜ(〈(∂|ψf,γ1,…,γr〉/∂γk)∣ψ〉)∑j=1rθj′))·(∑k=1rℑ(〈ψf,γ1,…,γr|Hk|ψ〉)vk(t)(λγ1,…,γr+ω)ℑ(ψf,γ1,…,γr∣ψ)+∑k=1rℑ(〈ψf,γ1,…,γr|Hk|ψ〉)vk(t)).
In order to make V˙7(t)≤0, ω and vk(t) need to be designed as
(47)ω=-λγ1,…,γr+cf0(ℑ(〈ψf,γ1,…,γr∣ψ〉)),(48)vk(t)=Kkfk(ℑ(〈ψf,γ1,…,γr|Hk|ψ〉)),(k=1,…,r),
where Kk>0, c>0, and yk=fk(xk),(k=0,1,…,r) are monotonic functions through the coordinate origin and in the first quadrant and the third quadrant.Based on the LaSalle invariance principle, the convergence of the control system governed by (44) can be depicted by Theorem 12.Theorem 12 (see [38, 39]).
Consider the control system (43) with control fields γk(t) designed in (45) and Lemma 11, vk(t) designed in (48), and ω designed in (47). Assume that the target state |ψf〉 is an eigenstate of H0. If the control system satisfies (i) ωl,m,γ1,…,γr≠ωi,j,γ1,…,γr,(l,m)≠(i,j), i,j,l,m∈{1,…,N}, ωl,m,γ1,…,γr=λl,γ1,…,γr-λm,γ1,…,γr,λl,γ1,…,γr is the lth eigenvalue of (H0+∑k=1rHk(γk+ηk)) corresponding to the eigenstate |ϕl,γ1,…,γr〉; (ii) There is at least a k∈{1,…,r}, 〈ϕγ1(t),…,γr(t)|Hk|ψf,γ1(t),…,γr(t)〉≠0, for |ϕγ1(t),…,γr(t)〉≠|ψf,γ1(t),…,γr(t)〉; then the largest invariant set is S2N-1∩E8, where E8={|ψ〉=eiθ|ψf〉,θ∈R}. And the control system will converge toward the equivalent state of the target state: |ψf〉eiθ,(θ∈R).By designing appropriate perturbationsγk(t), conditions (i) and (ii) in Theorem 12 can be satisfied in most cases. From Theorem 12, one can see that if one design appropriate control laws uk(t)=γk(t)+vk(t),(k=1,…,r) and ω to make the conditions (i) and (ii) in Theorem 12 hold, the control system depicted by (43) can converge to the equivalent state of the target eigenstate |ψf〉 from any initial pure state.


### 3.1.3. Implicit Lyapunov Control Based on Average Value of an Imaginary Mechanical Quantity

For the implicit Lyapunov control based on the average value of an imaginary mechanical quantity, one can consider the control system depicted by (37). In this control method, the introduced implicit function perturbations γk(t) mainly have two tasks. One task is to solve the convergent problem of the control system in the degenerate cases. The basic idea is similar to that of the state distance. Another one is to choose a simpler restriction V(|ψf〉)<V(|ψother〉), which can be satisfied for any initial state and any target state by designing the imaginary mechanical quantity. In order to ensure the system converge to the target state by adding V(|ψf〉)<V(|ψother〉), we can design all the perturbations γk(t)=0 hold for k=1,…,r only at |ψf〉, that is, (1) γk(|ψf〉)=0,(k=1,…,r), and (2) for |ψ〉≠|ψf〉, there exists at least one k such that γk(|ψ〉)≠0. For sake of completing these two tasks, we can design γk(t) as a monotonically increasing functional of V(t):
(49)γk(|ψ〉)=Ck·θk(V(|ψ〉)-V(|ψf〉)),
where Ck≥0, and for k=1,…,r, there exists at least a Ck>0, and the function θk(·) satisfies θk(0)=0, θk(s)>0 and θk′(s)>0 for every s>0, s is the independent variable of the function θk(·).The specific Lyapunov function based on the average value of an imaginary mechanical quantity can be selected as:(50)V8(|ψ〉)=〈ψ|Pγ1,…,γr|ψ〉,
where Pγ1,…,γr=f(γ1(t),…,γr(t)) is a functional of γk(t) and positive definite.The existence ofγk(t) can be depicted by Lemma 13.Lemma 13.
IfCk=0, γk(|ψ〉)=0. Else if Ck>0,θk∈C∞(R+;[0,γk*]),k=1,…,r(γk* is a positive constant) satisfy θk(0)=0,θk(s)>0 and θk′(s)>0 for every s>0, and |θk′|<1/(2CkC*), C*=1+C, C=max{∥∂Pγ1,…,γr/∂γk∥∞,γk∈[0,γk*]}, then for every |ψ〉∈S2N-1, there is a unique γk∈C∞(γk∈[0,γk*]) satisfying γk(|ψ〉)=Ck·θk(〈ψ|Pγ1,…,γr|ψ〉-〈ψf|Pγ1,…,γr|ψf〉).
For the sake of simplicity,γk(t) can be designed zero for some k, and other γk(t) equal, denoted by γ(t); that is, set
(51)γk(t)=γ(t)=θ(〈ψ|Pγ|ψ〉-〈ψf|Pγ|ψf〉),k=k1,…,km;γk(t)=0,k≠k1,…,km(1≤k1,…,km≤r),
where θ(·)=θk1(·)=⋯=θkm(·) and Pγ are functionals of γ(t).
The time derivative of the selected Lyapunov function is(52)V˙8=∑k=1rivk(t)〈ψ|[Hk,Pγ]|ψ〉·(1+θ′(〈ψf|(∂Pγ/∂γ)|ψf〉))(1-θ′(〈ψ|(∂Pγ/∂γ)|ψ〉-〈ψf|(∂Pγ/∂γ)|ψf〉)).
In order to ensure V˙8(t)≤0,vk(t) can be designed as
(53)vk(t)=-Kkfk(i〈ψ|[Hk,Pγ]|ψ〉),(k=1,…,r),
where Kk is a constant and Kk>0, and yk=fk(xk),(k=1,2,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.Based on LaSalle’s invariance principle, the convergence of the control system governed by (37) can be depicted by Theorem 14.Theorem 14.
Consider the control system governed by (37) with control fields uk(t)=γk(t)+vk(t),(k=1,…,r), where γk(t) is defined by Lemma 13 and (51) and vk(t) is defined by (53). If the control system satisfies (i) ωl,m,γ≠ωi,j,γ,(l,m)≠(i,j), i,j,l,m∈{1,2,…,N}, ωl,m,γ=λl,γ-λm,γ, where λl,γ is the lth eigenvalue of H0+∑n=k1kmHnγ(t) corresponding to the eigenstate |ϕl,γ〉; (ii) for any i≠j,i,j∈{1,2,…,N}, there exists at least one k such that (H^k)lm≠0, where (H^k)lm is the (l,m)th element of H^k=U1†HkU1, U1=(|ϕ1,γ〉,…,|ϕN,γ〉); (iii) [Pγ,H0+∑n=k1kmHnγ(t)]=0; (iv) (P^γ)ll≠(P^γ)mm,l≠m, where (P^γ)ll is the (l,l)th element of U1HPγU1, then any trajectory will converge toward E9={|ψt0〉|eiθl|ϕl,γ(|ψt0〉)〉;θl∈R,l∈{1,…,N}}.From Theorem14, one can see that if the target state |ψf〉 is an eigenstate, |ψf〉 is contained in E9. In order to make the system converge to the target eigenstate from any initial pure state |ψ0〉, on the one hand, Pγ needs to be designed to make
(54)V8(|ψf〉)<V8(|ψother〉)
hold, where |ψother〉 represents any other state in the set E9 except the target state. On the other hand, because ∂γ/∂V>0,V˙≤0,γ≥0 holds, when vk(t)=0,γ(t)=γ-≠0 holds for some time, we can design γ=γ--α,(0<α≪γ-) to make the state trajectory evolve but not stay in E9 until |ψf〉eiθl is reached.Next, let us analyze how to realize the convergence conditions in Theorem14 and restriction defined by (54). By designing appropriate γk(t), conditions (i) and (ii) in Theorem 14 can be satisfied in most cases. In order to make condition (iii) be satisfied, the eigenvalues of Pγ can be designed constant, denoted by P1,P2,…,PN, and Pγ can be designed as
(55)Pγ=∑j=1NPj|ϕj,γ〉〈ϕj,γ|.To make condition (iv) satisfied,Pl≠Pj(∀l≠j;1≤l,j≤N) can be designed. The results on how to make (54) hold can be depicted by Theorem 15.Theorem 15.
If one designsPi>Pf,(i=1,…,N,Pi≠Pf), then V8(|ψf〉)<V8(|ψother〉) holds, where Pf is the eigenvalue of Pγ(|ψf〉) corresponding to |ψf〉.We can conclude from the above analyses that by using the implicit Lyapunov control based on the imaginary mechanical quantity, if one designs appropriate control lawsuk(t)=γk(t)+vk(t),(k=1,…,r) to make the conditions (i)–(iv) in Theorem 15 and (54) hold, the control system governed by (37) in the degenerate cases can converge from any initial pure state to the target eigenstate. The design principle of the imaginary mechanical quantity Pγ is Theorem 15 and (55).


## 3.1.1. Implicit Lyapunov Control Based on State Distance

In order to solve the convergence of the control systems in the degenerate cases, several researchers introduced a series of implicit function perturbations into the control laws [37–40]. After a series of perturbations γk(t) introducing into the control laws, the dynamical equation (1) becomes
(37)i|ψ˙(t)〉=(H0+∑k=1rHk(γk(t)+vk(t)))|ψ(t)〉,
where γk(t)+vk(t)=uk(t),(k=1,…,r) are the total control laws.Without loss of generality, assume|ψf〉=|ϕg〉,1≤g≤N. In order to solve the problem of convergence for the degenerate cases, the perturbations γk(t) were introduced into the control laws. The basic idea is as follows: Denote the system with the internal Hamiltonian H0, the control Hamiltonians Hk,(k=1,…,r), and the control laws uk(t)=γk(t)+vk(t) as system 1, and the system with the internal Hamiltonian H02=H0+∑k=1rHkγk(t), the control Hamiltonians Hk,(k=1,…,r), and the control laws vk(t) as system 2. All these two systems can be depicted by (37). Denote the eigenvalues and eigenstates of H02 as λ1,γ1(t),…,γr(t),λ2,γ1(t),…,γr(t),…,λN,γ1(t),…,γr(t) and |ϕ1,γ1(t),…,γr(t)〉,|ϕ2,γ1(t),…,γr(t)〉,…,|ϕN,γ1(t),…,γr(t)〉, respectively, which are the functions of the perturbations γk(t). Assume |ψf,γ1(t),…,γr(t)〉=|ϕg,γ1(t),…,γr(t)〉. If one can design the perturbations γk(t) such that (i) ωl,m,γ1,…,γr≠ωi,j,γ1,…,γr,(l,m)≠(i,j),i,j,l,m∈{1,2,…,N}, where ωl,m,γ1,…,γr=λl,γ1,…,γr-λm,γ1,…,γr holds; (ii) for any |ϕγ1,…,γr〉≠|ψf,γ1,…,γr〉, there exists at least a k∈{1,…,r} satisfying 〈ϕγ1,…,γr|Hk|ψf,γ1,…,γr〉≠0, and select the specific Lyapunov function based on the state distance as
(38)V6(t)=12(1-|〈ψ∣ψf,γ1(t),…,γr(t)〉|2).
Then according to Section 2.1.1, system 2 will converge to |ψf,γ1(t),…,γr(t)〉. And when system 2 converge to |ψf,γ1(t),…,γr(t)〉, if the perturbations γk(t) at |ψf,γ1(t),…,γr(t)〉 are designed to equal zero, then system 2 will become system 1, and |ψf,γ1(t),…,γr(t)〉 will become |ψf〉. Then the convergence of system 1 to |ψf〉 will be ensured. In fact, the evolution of system 1 can be viewed as a composite of two evolution processes. One is system 2 converges to |ψf,γ1(t),…,γr(t)〉 from the initial state |ψ0〉, another one is γk(t) converges to 0. In order to make the introduced perturbations take effect to make system 1 in the non-degenerate case converge to |ψf〉, the speed of γk(t),(k=1,…,r) converging to 0 must be slower than the speed at which system 2 converges toward |ψf,γ1(t),…,γr(t)〉. For convenience, the control system in the following section means system 1.The existing design method ofγk(t) is to design it to be a monotonically increasing functional of V(t), that is,
(39)γk(t)=θk(V(t))=θk(12(1-|〈ψ∣ψf,γ1(t),…,γr(t)〉|2)),(k=1,…,r),
where the function θk(·) satisfies θk(0)=0, θk(s)>0 and θk′(s)>0 for every s>0 and s is the independent variable of the function θk(·).The right side of (39) contains |ψf,γ1(t),…,γr(t)〉 which is a functional of the perturbations γk(t),(k=1,…,r). One can see that the relation between γk and the time t is defined by r equations and cannot be expressed by an explicit expression, so γk is the implicit function of the time t. From (38) and (39), one can also see that the Lyapunov function V6 is the implicit function of the time t. The existence of γk(t) can be depicted by Lemma 9.Lemma 9 (see [37, 40]).
Letθk∈C∞(R+;[0,γk*]),k=1,…,r,(γk*∈R,γk*>0) satisfy θk(0)=0, θk(s)>0, θk′(s)>0 for every s>0, ∥θk′∥∞<1/C*, and C*=1+max{∥(∂ψf,γ1(t),…,γr(t)/∂γk)|(γ10,…,γr0)∥;γk0∈[0,γk*]}, where γk0=γk(0). Then for every state |ψ〉∈S2N-1={x∈CN;∥x∥=1}, there exists a unique γ1,γ2,…, and γr with γk∈C∞(γk∈[0,γk*]),(k=1,…,r) satisfying γk(t)=θk((1/2)(1-|〈ψ∣ψf,γ1(t),…,γr(t)〉|2)),(k=1,…,r),γk(|ψf〉)=0.By (38) and Lemma 9, one can see that for every state |ψ〉∈S2N-1={x∈CN;∥x∥=1}; there also exists a unique V.Assume(40)∂|ψf,γ1(t),…,γr(t)〉∂γ1=∂|ψf,γ1(t),…,γr(t)〉∂γ2=⋯=∂|ψf,γ1(t),…,γr(t)〉∂γr.
The first order time derivative of the Lyapunov function (38) is
(41)V˙6=-∑k=1rvk(t)ℑ(〈ψ∣ψf,γ1,…,γr〉〈ψf,γ1,…,γr|Hk|ψ〉)(1+ℜ(〈ψ∣ψf,γ1,…,γr〉〈(∂|ψf,γ1,…,γr〉/∂γk)∣ψ〉)∑j=1rθj′).
In order to make V˙6(t)≤0, one can design γk(t) such that ∥θj′∥∞<1/(rC*) to make 1/(1+ℜ(〈ψ∣ϕγ1,…,γr〉〈(∂ϕγ1,…,γr/∂γk)∣ψ〉)∑j=1rθj′)>0 hold and design vk(t) as
(42)vk(t)=Kkfk(ℑ(ei∠〈ψ∣ψf,γ1,…,γr〉〈ψf,γ1,…,γr|Hk|ψ〉)),
where Kk>0, and vk(t),(k=1,…,r) are monotonic functions through the coordinate origin and in the first quadrant and the third quadrant.Based on the LaSalle invariance principle, the convergence of the control system governed by (37) can be depicted by Theorem 10.Theorem 10 (see [37, 40]).
Consider the control system (37) with control fields uk(t)=γk(t)+vk(t),γk(t) designed in (39), Lemma 9 and ∥θj′∥∞<1/(rC*), and vk(t) designed in (42). Assume that the target state |ψf〉 is an eigenstate of H0. If the control system satisfies (i) ωl,m,γ1,…,γr≠ωi,j,γ1,…,γr,(l,m)≠(i,j), i,j,l,m∈{1,2,…,N}, ωl,m,γ1,…,γr=λl,γ1,…,γr-λm,γ1,…,γr, where λl,γ1,…,γr is an eigenvalue of (H0+∑k=1rHkγk), (ii) For any |ϕγ1,…,γr〉≠|ψf,γ1,…,γr〉, there exists at least a k∈{1,…,r} satisfying 〈ϕγ1,…,γr|Hk|ψf,γ1,…,γr〉≠0; then the largest invariant set is S2N-1∩E7, where E7={|ψ〉=eiθ|ψf〉,θ∈R}. And the control system will converge toward |ψf〉eiθ,(θ∈R).Conditions (i) and (ii) in Theorem10 are associated with H0, Hk,(k=1,…,r) and γk(t). By designing appropriate γk(t), these two conditions can be satisfied in most cases. one can see that if one designs appropriate control laws uk(t)=γk(t)+vk(t),(k=1,…,r) to make the conditions (i) and (ii) hold, the control system depicted by (37) can converge to the equivalent state of the target eigenstate |ψf〉 from any initial pure state.


## 3.1.2. Implicit Lyapunov Control Based on the State Error

The basic idea of this method is similar to that of the state distance. In the Lyapunov control based on the state error, a global phase control itemω is added into the control system to facilitate the design of the control laws. Thus the control system can be depicted by
(43)i|ψ˙(t)〉=(H0+∑k=1rHk(γk(t)+vk(t))+ωI)|ψ(t)〉,
where uk(t)=γk(t)+vk(t),(k=1,…,r) and ω is the control laws which need to design.Consider the Lyapunov function based on the state error as(44)V7(t)=12〈ψ-ψf,γ1,…,γr∣ψ-ψf,γ1,…,γr〉.
Thus the implicit function perturbation γk(t) needs to be designed as
(45)γk(t)=θk(V(t))=θk(12〈ψ-ψf,γ1,…,γr∣ψ-ψf,γ1,…,γr〉),
where functions θk(·) satisfy θk(0)=0, θk(s)>0, θk′(s)>0 for every s>0.The existence ofγk(t) can be depicted by Lemma 11.Lemma 11 (see [38, 39]).
Letθk∈C∞(R+;[0,γk*]),k=1,…,r,(γk*∈R,γk*>0) satisfy θk(0)=0,θk(s)>0,θk′(s)>0 for every s>0,∥θk′∥∞<1/rC*, and C*=1+max{∥(∂|ψf,γ1,…,γr〉/∂γk)|(γ10,…,γr0)∥;γk0∈[0,γk*]}. Then for every state |ψ〉∈S2N-1={x∈CN;∥x∥=1}, there exists a unique γ1,γ2,…, and γr with γk∈C∞(γk∈[0,γk*]),(k=1,…,r) satisfying γk(|ψ〉)=θk((1/2)〈ψ-ψf,γ1(|ψ〉),…,γr(|ψ〉)∣ψ-ψf,γ1(|ψ〉),…,γr(|ψ〉)〉), γk(|ψf〉)=0.Assume (40) holds, the first order time derivative of the Lyapunov function (44) is
(46)V˙7=-(1(1+ℜ(〈(∂|ψf,γ1,…,γr〉/∂γk)∣ψ〉)∑j=1rθj′))·(∑k=1rℑ(〈ψf,γ1,…,γr|Hk|ψ〉)vk(t)(λγ1,…,γr+ω)ℑ(ψf,γ1,…,γr∣ψ)+∑k=1rℑ(〈ψf,γ1,…,γr|Hk|ψ〉)vk(t)).
In order to make V˙7(t)≤0, ω and vk(t) need to be designed as
(47)ω=-λγ1,…,γr+cf0(ℑ(〈ψf,γ1,…,γr∣ψ〉)),(48)vk(t)=Kkfk(ℑ(〈ψf,γ1,…,γr|Hk|ψ〉)),(k=1,…,r),
where Kk>0, c>0, and yk=fk(xk),(k=0,1,…,r) are monotonic functions through the coordinate origin and in the first quadrant and the third quadrant.Based on the LaSalle invariance principle, the convergence of the control system governed by (44) can be depicted by Theorem 12.Theorem 12 (see [38, 39]).
Consider the control system (43) with control fields γk(t) designed in (45) and Lemma 11, vk(t) designed in (48), and ω designed in (47). Assume that the target state |ψf〉 is an eigenstate of H0. If the control system satisfies (i) ωl,m,γ1,…,γr≠ωi,j,γ1,…,γr,(l,m)≠(i,j), i,j,l,m∈{1,…,N}, ωl,m,γ1,…,γr=λl,γ1,…,γr-λm,γ1,…,γr,λl,γ1,…,γr is the lth eigenvalue of (H0+∑k=1rHk(γk+ηk)) corresponding to the eigenstate |ϕl,γ1,…,γr〉; (ii) There is at least a k∈{1,…,r}, 〈ϕγ1(t),…,γr(t)|Hk|ψf,γ1(t),…,γr(t)〉≠0, for |ϕγ1(t),…,γr(t)〉≠|ψf,γ1(t),…,γr(t)〉; then the largest invariant set is S2N-1∩E8, where E8={|ψ〉=eiθ|ψf〉,θ∈R}. And the control system will converge toward the equivalent state of the target state: |ψf〉eiθ,(θ∈R).By designing appropriate perturbationsγk(t), conditions (i) and (ii) in Theorem 12 can be satisfied in most cases. From Theorem 12, one can see that if one design appropriate control laws uk(t)=γk(t)+vk(t),(k=1,…,r) and ω to make the conditions (i) and (ii) in Theorem 12 hold, the control system depicted by (43) can converge to the equivalent state of the target eigenstate |ψf〉 from any initial pure state.


## 3.1.3. Implicit Lyapunov Control Based on Average Value of an Imaginary Mechanical Quantity

For the implicit Lyapunov control based on the average value of an imaginary mechanical quantity, one can consider the control system depicted by (37). In this control method, the introduced implicit function perturbations γk(t) mainly have two tasks. One task is to solve the convergent problem of the control system in the degenerate cases. The basic idea is similar to that of the state distance. Another one is to choose a simpler restriction V(|ψf〉)<V(|ψother〉), which can be satisfied for any initial state and any target state by designing the imaginary mechanical quantity. In order to ensure the system converge to the target state by adding V(|ψf〉)<V(|ψother〉), we can design all the perturbations γk(t)=0 hold for k=1,…,r only at |ψf〉, that is, (1) γk(|ψf〉)=0,(k=1,…,r), and (2) for |ψ〉≠|ψf〉, there exists at least one k such that γk(|ψ〉)≠0. For sake of completing these two tasks, we can design γk(t) as a monotonically increasing functional of V(t):
(49)γk(|ψ〉)=Ck·θk(V(|ψ〉)-V(|ψf〉)),
where Ck≥0, and for k=1,…,r, there exists at least a Ck>0, and the function θk(·) satisfies θk(0)=0, θk(s)>0 and θk′(s)>0 for every s>0, s is the independent variable of the function θk(·).The specific Lyapunov function based on the average value of an imaginary mechanical quantity can be selected as:(50)V8(|ψ〉)=〈ψ|Pγ1,…,γr|ψ〉,
where Pγ1,…,γr=f(γ1(t),…,γr(t)) is a functional of γk(t) and positive definite.The existence ofγk(t) can be depicted by Lemma 13.Lemma 13.
IfCk=0, γk(|ψ〉)=0. Else if Ck>0,θk∈C∞(R+;[0,γk*]),k=1,…,r(γk* is a positive constant) satisfy θk(0)=0,θk(s)>0 and θk′(s)>0 for every s>0, and |θk′|<1/(2CkC*), C*=1+C, C=max{∥∂Pγ1,…,γr/∂γk∥∞,γk∈[0,γk*]}, then for every |ψ〉∈S2N-1, there is a unique γk∈C∞(γk∈[0,γk*]) satisfying γk(|ψ〉)=Ck·θk(〈ψ|Pγ1,…,γr|ψ〉-〈ψf|Pγ1,…,γr|ψf〉).
For the sake of simplicity,γk(t) can be designed zero for some k, and other γk(t) equal, denoted by γ(t); that is, set
(51)γk(t)=γ(t)=θ(〈ψ|Pγ|ψ〉-〈ψf|Pγ|ψf〉),k=k1,…,km;γk(t)=0,k≠k1,…,km(1≤k1,…,km≤r),
where θ(·)=θk1(·)=⋯=θkm(·) and Pγ are functionals of γ(t).
The time derivative of the selected Lyapunov function is(52)V˙8=∑k=1rivk(t)〈ψ|[Hk,Pγ]|ψ〉·(1+θ′(〈ψf|(∂Pγ/∂γ)|ψf〉))(1-θ′(〈ψ|(∂Pγ/∂γ)|ψ〉-〈ψf|(∂Pγ/∂γ)|ψf〉)).
In order to ensure V˙8(t)≤0,vk(t) can be designed as
(53)vk(t)=-Kkfk(i〈ψ|[Hk,Pγ]|ψ〉),(k=1,…,r),
where Kk is a constant and Kk>0, and yk=fk(xk),(k=1,2,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.Based on LaSalle’s invariance principle, the convergence of the control system governed by (37) can be depicted by Theorem 14.Theorem 14.
Consider the control system governed by (37) with control fields uk(t)=γk(t)+vk(t),(k=1,…,r), where γk(t) is defined by Lemma 13 and (51) and vk(t) is defined by (53). If the control system satisfies (i) ωl,m,γ≠ωi,j,γ,(l,m)≠(i,j), i,j,l,m∈{1,2,…,N}, ωl,m,γ=λl,γ-λm,γ, where λl,γ is the lth eigenvalue of H0+∑n=k1kmHnγ(t) corresponding to the eigenstate |ϕl,γ〉; (ii) for any i≠j,i,j∈{1,2,…,N}, there exists at least one k such that (H^k)lm≠0, where (H^k)lm is the (l,m)th element of H^k=U1†HkU1, U1=(|ϕ1,γ〉,…,|ϕN,γ〉); (iii) [Pγ,H0+∑n=k1kmHnγ(t)]=0; (iv) (P^γ)ll≠(P^γ)mm,l≠m, where (P^γ)ll is the (l,l)th element of U1HPγU1, then any trajectory will converge toward E9={|ψt0〉|eiθl|ϕl,γ(|ψt0〉)〉;θl∈R,l∈{1,…,N}}.From Theorem14, one can see that if the target state |ψf〉 is an eigenstate, |ψf〉 is contained in E9. In order to make the system converge to the target eigenstate from any initial pure state |ψ0〉, on the one hand, Pγ needs to be designed to make
(54)V8(|ψf〉)<V8(|ψother〉)
hold, where |ψother〉 represents any other state in the set E9 except the target state. On the other hand, because ∂γ/∂V>0,V˙≤0,γ≥0 holds, when vk(t)=0,γ(t)=γ-≠0 holds for some time, we can design γ=γ--α,(0<α≪γ-) to make the state trajectory evolve but not stay in E9 until |ψf〉eiθl is reached.Next, let us analyze how to realize the convergence conditions in Theorem14 and restriction defined by (54). By designing appropriate γk(t), conditions (i) and (ii) in Theorem 14 can be satisfied in most cases. In order to make condition (iii) be satisfied, the eigenvalues of Pγ can be designed constant, denoted by P1,P2,…,PN, and Pγ can be designed as
(55)Pγ=∑j=1NPj|ϕj,γ〉〈ϕj,γ|.To make condition (iv) satisfied,Pl≠Pj(∀l≠j;1≤l,j≤N) can be designed. The results on how to make (54) hold can be depicted by Theorem 15.Theorem 15.
If one designsPi>Pf,(i=1,…,N,Pi≠Pf), then V8(|ψf〉)<V8(|ψother〉) holds, where Pf is the eigenvalue of Pγ(|ψf〉) corresponding to |ψf〉.We can conclude from the above analyses that by using the implicit Lyapunov control based on the imaginary mechanical quantity, if one designs appropriate control lawsuk(t)=γk(t)+vk(t),(k=1,…,r) to make the conditions (i)–(iv) in Theorem 15 and (54) hold, the control system governed by (37) in the degenerate cases can converge from any initial pure state to the target eigenstate. The design principle of the imaginary mechanical quantity Pγ is Theorem 15 and (55).


## 3.2. Target Superposition State

The method proposed in Section3.1 cannot guarantee the control system governed by (37) or (43) in the degenerate cases converge from any initial pure state to the target superposition state. In order to solve this problem, one can introduce a series of constant disturbances ηk into the control laws. Thus the mechanical equation (37) becomes
(56)i|ψ˙(t)〉=(H0+∑k=1rHk(ηk+γk(t)+vk(t)))|ψ(t)〉.
And the mechanical equation (43) in Section 3.1.2 becomes
(57)i|ψ˙(t)〉=(H0+∑k=1rHk(γk(t)+vk(t)+ηk)+ωI)|ψ(t)〉,
where ηk∈R.The basic idea of solving the convergence to the target superposition state is to designηk to make the target state |ψf〉 be an eigenstate of H0′=H0+∑k=1rHkηk·H0′ can be viewed as the new internal Hamiltonian of the control system. If the number of the control Hamiltonians r is large enough, by designing appropriate ηk, (H0+∑k=1rHkηk)|ψf〉=λf′|ψf〉 can be satisfied in most cases, where λf′ is the eigenvalue of H0+∑k=1rHkηk corresponding to |ψf〉. Then one can design the control laws and analyze the convergence according to the method for the target eigen state cases. Research results show that every conclusion in Section 3.1 also holds with changing H0 into H0′.


## 3.3. Target Mixed State Which Commutes with the Internal Hamiltonian

Consider theN-level closed quantum system governed by the following quantum Liouville equation:
(58)iρ˙(t)=[H0+∑k=1rHk(γk(t)+vk(t)),ρ(t)],
where γk(t)+vk(t)=uk(t) are the total control laws.The design ideas are similar to those of Section3.1. Consider the Lyapunov function based on the imaginary mechanical quantity as
(59)V9(ρ)=tr(Pγ1,…,γrρ),γk(t) can be designed as
(60)γk(t)=γ(t)=θ(V9(ρ)-V9(ρf)),k=k1,…,km;γk(t)=0,k≠k1,…,km(1≤k1,…,km≤r).The existence ofγ(t) can be depicted by Lemma 16.Lemma 16.
Ifθ∈C∞(R+;[0,γ*]),k=1,…,r (γ* is a positive constant) satisfy θ(0)=0,θ(s)>0 and θ′(s)>0 for every s>0, and |θ′|<1/(2C*), C*=1+C, C=max{∥∂Pγ/∂γ∥m1,γ∈[0,γ*]}, then for every ρ, there is a unique γ∈C∞(γ∈[0,γ*]) satisfying γ(ρ)=θ(tr(Pγρ)-tr(Pγρf)).According to the Lyapunov stabilty Theorem,vk(t) can be designed as
(61)vk(t)=Kkfk(itr([Pγ,Hk]ρ)),(k=1,…,r),
where Kk is a constant and Kk>0, and yk=fk(xk),(k=1,2,…,r) are monotonic increasing functions through the coordinate origin of the plane xk-yk.Based on the LaSalle invariance principle, the convergence of the control system can be depicted by Theorem17.Theorem 17.
Consider the control system depicted by (58) with control laws uk(t)=γk(t)+vk(t), where γk(t) is defined by Lemma 16 and (60), and vk(t) is defined by (61). If the control system satisfies (i) ωl,m,γ≠ωi,j,γ,(l,m)≠(i,j),i,j,l,m∈{1,2,…,N}, ωl,m,γ=λl,γ-λm,γ, where λl,γ is the lth eigenvalue of (H0+∑n=k1kmHnγ(t)) corresponding to the eigenstate |ϕl,γ〉; (ii) ∀j≠l, for k=1,…,r, there exists at least a (H^k)jl≠0, where (H^k)jl is the (j,l)th element of H^k=U2HHkU2 with U2=(|ϕ1,γ〉,…,|ϕN,γ〉); (iii) [Pγ,H0+∑n=k1kmHnγ(t)]=0,1≤k1,…,km≤r; (iv) for any l≠j,(1≤l,j≤N),(P^γ)ll≠(P^γ)jj holds, where (P^γ)ll is the (l,l)th element of P^γ=U2HPγU2, then the control system will converge toward E10={ρt0|(U2†ρt0U2)ij=0,γ=γ(ρt0),t0∈R}.For the target stateρf which commutes with the internal Hamiltonian; that is, [ρf,H0]=0, ρf is contained in E10. In order to make the system converge to the target state ρf which commutes with H0 from any initial stat; on the one hand, Pγ needs to be designed to make
(62)V9(ρf)<V9(ρother)
hold, where ρother represents any other state in the set E10 except the target state. On the other hand, γ was designed as γ=γ--α,(0<α≪γ-) when vk(t)=0,γ(t)=γ-≠0 holds for some time to make the state trajectory evolve but not stay in E10 until ρf is reached.For the satisfaction of conditions (i)–(iv) in Theorem17, one can follow that of Section 3.1.3. In order to make (62) hold, one can design Pγ based on Proposition 6. The total design principle of Pγ is (55) and Proposition 6.We can conclude from the above analyses that for the target state which commutes with the internal Hamiltonian, by using the implicit Lyapunov control method based on the imaginary mechanical quantity, if one designs appropriate control lawsuk(t)=γk(t)+vk(t),(k=1,…,r) to make the conditions (i)–(iv) in Theorem 17 and (62) hold, the control system governed by (58) in the degenerate cases can converge from any initial state to the target state which commutes with the internal Hamiltonian, which contains the target eigenstate and the target mixed state which commutes with the internal Hamiltonian.


## 3.4. Target Mixed State Which Does Not Commute with Internal Hamiltonian

For the target mixed state which does not commute with the internal Hamiltonian, the design idea is similar to Section3.2. The difference is to design ηk to make [ρf,H0′]=0. If the number of the control Haimltonians r is large enough, by designing appropriate ηk, [ρf,H0′]=0 can be satisfied in most cases. Then one can design the control laws and analyze the convergence according to the method mentioned in Section 3.3. Research results show that every conclusion in Section 3.3 also holds with changing H0 into H0′.


## 4. Conclusion

In this paper, for the non-degenerate and degenerate cases, the existing quantum Lyapunov control based on the state distance, state error, and average value of an imaginary mechanical quantity for the control systems have been summarized and analyzed. For the target state being the eigenstate, the mixed state which commutes with the internal Hamiltonian, the superposition state, and the mixed state which does not commute with the internal Hamiltonian, respectively, the design methods of the control laws have been summarized; the convergence to the target state has been summarized and analyzed. Research results show that the Lyapunov-based quantum control method can make the control system converge from any initial state to the target state in both non-degenerate and degenerate cases. After ten years of development, the quantum control theory based on the Lyapunov stability theorem has been established.



---

*Source: 967529-2013-05-27.xml*