A survey of quantum Lyapunov control methods.

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.

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 equation and the quantum Liouville equation , in which |ψ(t)〉 is the quantum state vector, ρ(t) is the density operator, H 0 is the internal Hamiltonian, and H , (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: , 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].

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.

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 Section 2, 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 function V(t) satisfy . 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 the N-level closed quantum system governed by the following bilinear Schrödinger equation: where |ψ(t)〉 is the quantum state vector, H 0 is the internal Hamiltonian, H , (k = 1,…, r) are control Hamiltonians, and u (t), (k = 1,…, r) are scalar and real control laws.

In the following sections, in the case that the target state |ψ 〉 is an eigenstate, that is, H 0|ψ 〉 = λ |ψ 〉, where λ is the eigenvalue of the internal Hamiltonian H 0, 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: The time derivative of the Lyapunov function (2) is The control laws u (t) which can make hold can be designed as where K > 0, and y = f (x ), (k = 1,…, r) are monotonic increasing functions through the coordinate origin of the plane x − y .

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 u (t) designed in (4). If (i) The system is strongly regular, that is, ω ≠ ω , (i′, j′)≠(l, m), i′, j′, l, m ∈ {1,2,…, N}, ω = λ − λ , where λ is the lth eigenvalue of H 0 corresponding to the eigenstate |ϕ l〉; (ii) for any |ϕ 〉 ≠ |ψ 〉, i ∈ {1,…, N}, there exists at least a k ∈ {1,…, r} such that 〈ϕ |H |ψ 〉 ≠ 0, then any state trajectory will converge toward E 1 = {e |ψ 〉, θ ∈ R}.

From Theorem 1, 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 e |ψ 〉 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:

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

After some deduction, one can obtain the time derivative of the Lyapunov function (5) as The control laws which can make hold can be designed as where K > 0, and y = f (x ), (k = 0,…, r) are the monotonic increasing functions through the coordinate origin of the plane x − y .

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 u (t) designed in (9) and ω designed in (8). If (i) ω ≠ ω , (i′, j′)≠(l, m), i′, j′, l, m ∈ {1,2,…, N}, ω = λ − λ , where λ l is the lth eigenvalue of H 0 corresponding to the eigenstate |ϕ l〉; (ii) for any |ϕ 〉 ≠ |ψ 〉, i ∈ {1,…, N}, there exists at least a k ∈ {1,…, r} such that 〈ϕ |H |ψ 〉 ≠ 0. Then any state trajectory will converge toward E 2 = {e |ψ 〉, θ ∈ R}.

From Theorem 2, 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 e |ψ 〉 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: 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

Set [H 0, P] = 0 such that the drift term in the right side of (11) can be eliminated. In order to ensure , one can design u (t) as where K > 0, and y = f (x ), (k = 1,…, r) are monotonic increasing functions through the coordinate origin of the plane x − y .

Then based on the LaSalle invariance principle, all the state trajectories of the system will converge to the invariant set contained in the set E in which holds. Denote the state at the time t as |ψ(t)〉 = ∑ c (t)|ϕ 〉, where c (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 〈ϕ |H |ϕ 〉 ≠ 0, then one can deduce that holds for all t ≥ t 0, t 0 ∈ R is equivalent to where P and P are the lth and jth eigenvalues of P, respectively.

If the target state |ψ 〉 is an eigenstate, and all the eigenvalues of P are designed mutually different, that is, P ≠ P for any l ≠ j, (l, j = 1,2 …, N), then (13) is equivalent to Equation (14) implies that there is at most one c ≠ 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 u (t) designed in (12). If (i) ω ≠ ω , (i′, j′)≠(l, m), i′, j′, l, m ∈ {1,…, N}, ω = λ − λ , where λ is the lth eigenvalue of H 0 corresponding to the eigenstate |ϕ l〉; (ii) for any i ≠ j, i, j ∈ {1,…, N}, there exsits at least a k ∈ {1,…, r} such that 〈ϕ |H |ϕ 〉 ≠ 0; (iii) [H 0, P] = 0; (iv) for any l ≠ m ∈ {1,2,…, N}, P ≠ P holds, where P is the lth eigenvalue of P. Then any state trajectory will converge toward E 3 = {e |ϕ 〉, i ∈ {1,…, N}, θ ∈ R}.

From Theorem 3, 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 |ψ 〉 from any initial pure state |ψ 0〉, as , one can add a restriction as where |ψ 0〉 is the initial state, |ψ other 〉 represents any other state in the set E 3 except the target state.

In such a way, any state trajectory of the system will converge to e |ψ 〉 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 Hamiltonian H 0 and the control Hamiltonians H , (k = 1,…, r) which cannot be changed by designing appropriate control laws. Condition (iii) means that P and H 0 have the same eigenstates. In order to make condition (iii) be satisfied, the eigenvalues of P can be designed as Because P should be positive definite, one needs to design P > 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 V 3(|ψ 〉) < V 3(|ψ 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 V 3(|ψ〉) = 〈ψ|P|ψ〉 is given by the normalized eigenvectors of P. The eigenvectors with the largest eigenvalue are the maxima of V 3, the eigenvectors with the smallest eigenvalue are the minima and all others are saddle points.

According to Proposition 4, in order to make V 3(|ψ 〉) < V 3(|ψ other 〉) hold, P > P , (i = 1,…, N, P ≠ P ) needs to be designed, where P is the eigenvalue of P corresponding to |ψ 〉. Then let us consider the whole restriction (15); it is an attraction problem. If the eigenvalues of P except P are close together, the attraction region will be very large. For the limiting case P = P > P > 0, (i ≠ j, P , P ≠ P ), the attraction region will be the whole state space. Thus the design principle of P is P > P , (i = 1,…, N, P ≠ P ) and to make P , P ≠ P close together.

In conclusion, for the target state being an eigenstate, the design principle of P is as follows:

P ≠ P > 0 for any l ≠ m ∈ {1,2,…, N},

P > P , (i = 1,…, N, P ≠ P ), P , P ≠ P are close together,

Equation (16): P = ∑ P |ϕ 〉〈ϕ |.

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 e |ψ 〉 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

The inner product of two operators A 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 ρ can be deduced as

One can conclude from (18) that the Lyapunov function based on the state distance V 1(|ψ〉) = (1/2)(1 − |〈ψ∣ψ 〉|2) and the state error V 2(|ψ〉) = (1/2)〈ψ − ψ | ψ − ψ 〉 are equivalent.

For the Lyapunov function based on the average value of an imaginary mechanical quantity V 3(|ψ〉) = 〈ψ|P|ψ〉, if the imaginary mechanical quantity P = (1/2)(I − |ψ 〉〈ψ |), V 3(|ψ〉) will become V 1(|ψ〉) [19]. Therefore V 1(|ψ〉) is formally a special case of V 3(|ψ〉).

In fact, V 1(|ψ〉), V 2(|ψ〉), and V 3(|ψ〉) can be unified in the following quadratic Lyapunov function:

Next, we consider three cases as follows:

α = 0, Q = (1/2)(I − |ψ 〉〈ψ |), V 4(|ψ〉) will reduce to V 1(|ψ〉);

α = 1, Q = (1/2)I, V 4(|ψ〉) will reduce to V 2(|ψ〉);

α = 0, V 4(|ψ〉) will reduce to V 3(|ψ〉).

We can see that V 1(|ψ〉), V 2(|ψ〉), and V 3(|ψ〉) are special cases of V 4(|ψ〉).

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, V 1(|ψ〉) and V 2(|ψ〉) are equivalent, so the Lyapunov control based on the state distance and the state error have similar control effects. Because V 1(|ψ〉) is formally a special case of V 3(|ψ〉), 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 the N-level closed quantum system governed by the following quantum Liouville equation: where ρ(t) is the density operator.

Consider the Lyapunov function based on the average value of an imaginary mechanical quantity: By means of setting [H 0, P] = 0, the first order time derivative of the Lyapunov function (21) can be deduced as

In order to ensure , one can design u (t) as where ε ∈ R, ε > 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 holds for all t ≥ t 0, t 0 ∈ R is equivalent to where ρ (t) is the (l; j)th element of the state ρ(t). If the target state ρ commutes with H 0, that is, [H 0, ρ ] = 0, and all the eigenvalues of P are designed mutually different, that is, P ≠ P for any l ≠ j, l, j = 1,2,…, N, then (24) is equivalent to which implies that the system will converge to a state which commutes with H 0. 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 u (t) designed in (23). If (i) the internal Hamiltonian is strongly regular, that is, ω ≠ ω , (i′, j′)≠(l, m), i′, j′, l, m ∈ {1,…, N}, ω = λ − λ , where λ l is the lth eigenvalue of H 0 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 (H ) ≠ 0; (iii) [H 0, P] = 0; (iv) P ≠ P , l ≠ m ∈ {1,2,…, N}, where P is the lth eigenvalue of P. Then any state trajectory will converge toward E 4 = {ρ | (ρ) = 0, i, j ∈ {1,…, N}}, where (ρ) is the (l, l)th element of ρ.

From Theorem 5, 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 in E 4 as ρ ; then [ρ , H 0] = 0 holds which implies that ρ and H 0 have the same eigenstates. Since the evolution of ρ(t) is unitary, ρ(t) for t ≥ 0 are isospectral. So the eigenvalues of ρ are a permutation of the eigenvalues of ρ 0. Thus the set E 4 has countable elements. If the initial state is generic, that is, the eigenvalues of the initial state are mutually different, the set E 4 will have N! elements. For the target state ρ which commutes with the internal Hamiltonian, that is, [ρ , H 0] = 0, in order to make the system converge to the target state ρ which commutes with H 0 from any initial state ρ 0, Kuang and Cong proposed a restriction as where ρ  other represents any other state in the set E 4 except the target state ρ .

In such a way, if there exists a P to make the restriction (26) hold, any state trajectory of the system will converge to the target state ρ which commutes with H 0 from any initial state ρ 0.

For how to make conditions (iii)-(iv) to be satisfied, please read Section 2.1.3. One can deduce the design principle of P such that V 5(ρ ) < V 5(ρ  other ), the result can be depicted by Proposition 6.

Proposition 6

If (ρ ) < (ρ ), 1 ≤ i, j ≤ N holds, then one can design P > P ; if (ρ ) = (ρ ), 1 ≤ i, j ≤ N holds, then one can design P ≠ P ; else if (ρ ) > (ρ ), 1 ≤ i, j ≤ N holds, then one can design P < P , thus V 5(ρ ) < V 5(ρ  other ) holds, where (ρ ) is the (i, i)th element of ρ .

It is difficult to design P such that (26) holds for any initial state ρ 0 and any target state ρ which satisfies [ρ , H 0] = 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 ρ which commutes with the internal Hamiltonian, the design principle of P is as follows:

P ≠ P > 0 for any l ≠ m ∈ {1,2,…, N};

Design P according to Proposition 6;

Equation (16): P = ∑ P |ϕ 〉〈ϕ |.

In conclusion, if the control system satisfies the conditions (i)–(iv) in Theorem 5 and there exists a P to make the restriction (26) hold, any state trajectory of the system will converge to the target mixed state ρ which commutes with the internal Hamiltonian from any initial state ρ 0.

2.3. Target Superposition State

From Section 2.1.3, one can see that for the control system governed by (1), holds for all t ≥ t 0, t 0 ∈ R is equivalent to (P − P )c c * = 0, (l, j = 1,…, N). If the target state is a superposition state as |ψ 〉 = ∑ c |ϕ 〉, (1 ≤ k1 < ko ≤ N), one can design P = P = ⋯ = P and other eigenvalues of P are mutually different, then (P − P )c c * = 0, (l, j = 1,…, N) is equivalent to Equation (27) implies that the system will converge to the set E 5 = span⁡{|ϕ 〉, |ϕ 〉,…, |ϕ 〉, |ϕ 〉; l ∈ {1,…, N}, l ≠ k1,…, ko}, which means that the set E 5 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 = c I, where c 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 Section 2.2, one can see that, for the control system governed by (20), holds for all t ≥ t 0, t 0 ∈ R is equivalent to (H )(P − P )ρ (t) = 0, j, l = 1,…, N, j < l. Consider the target state which does not commute with H 0, which includes the superposition state and the mixed state which does not commute with H 0. Without loss of generality, assume (ρ )3 = (ρ ) ≠ 0; one can design P 3 = P ; thus the system will converge to E 6 = {ρ|(ρ) = 0, (i, j) ≠ (3, N), (i, j) ≠ (N, 3), i, j ∈ {1,…, N}}. Because the set E 6 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 Sections 2.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 the N-level control system in the interaction picture as

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 To ensure , one can design u (t) as where K > 0 and y = f (x ), (k = 1,…, r) are monotonic increasing functions through the coordinate origin of plane x − y .

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 u (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 ≡ {ρ : [ρ , P] = D} at t → ∞, where D is a diagonal matrix.

For the case that P is chosen as a diagonal matrix, the limit set is reduced to ℛ 1 ≡ {ρ : [ρ , P] = 0}. For the case that P is chosen as a nondiagonal matrix, if holds, then the limit set ℛ 1 is regular; namely, ℛ 1 ≡ {ρ : [ρ , P] = 0}, where is the first n 2 − n rows of , is the real (n 2 − 1)∗(n 2 − 1) matrix corresponding to the Bloch representation of Ad , and Ad is a linear map from Hermitian or anti-Hermitian matrices into su(n). [ρ , P] = 0 means that P and H 0 have the same eigenstates. For the target state ρ = |ψ 〉〈ψ | being a superposition state, in order to make the target state contain in ℛ 1 ≡ {ρ : [ρ , P] = 0}, P can be designed as where |ψ 1〉 = |ψ 〉 and 〈ψ |ψ 〉 = 0, for i ≠ j.

Then some deduction shows that ρ = λ |ψ 〉〈ψ |(λ = 1). In order to make the system converge to the target superposition state ρ from any initial pure state, one can design P such that where ρ oth represents any other state in the set ℛ 1 except the target state ρ .

For the satisfaction of V 5(ρ ) < V 5(ρ ), one can design P based on Proposition 6. In order to make (32) holds, one must design P , (i = 1,…, N) such that

However, for any initial pure state the target superposition state, there may not exist a P such that (33) holds.

In conclusion, the design principle of P for the case that the target state ρ = |ψ 〉〈ψ | 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 ρ 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 Sections 2.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 ρ being a mixed state which does not commute with the internal Hamiltonian, that is, [ρ , H 0] ≠ 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 to make hold. Correspondingly, the control system governed by (28) after this unitary transformation will become where , .

Then control laws can be designed according to Section 2.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 contain in , P needs to be designed such that [H 0, P] = 0 holds. and [H 0, P] = 0 imply that , P, and H 0 have the same eigenstates. Since the evolution of is unitary, all the for t ≥ 0 are isospectral. So the eigenvalues of are a permutation of the eigenvalues of . Thus the limit set has countable elements. For sake of the convergence from any initial mixed state ρ 0 to any target mixed state ρ which does not commute with the internal Hamiltonian H 0, P needs to be designed such that where represents any other state in the set except .

In such a way, if there exists a P to make (35) hold, any state trajectory of the system will converge to the target mixed state ρ 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 ρ into accordingly, then will hold. Assume ρ 0 = Uρ U †. In order to make the whole restriction defined by (35) hold, P , (i = 1,…, N) needs to be designed such that where U is the (i, j)th element of U.

In conclusion, the design principle of P for the case that the target state ρ is a mixed state which does not commute with H 0 is

P ≠ P > 0 for any l ≠ m ∈ {1,2,…, N};

Proposition 6 with changing ρ into ;

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 with H 0, 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 Theorem 1 through Theorem 7 are relevant to the internal Hamiltonian H 0 and the control Hamiltonians H , (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 γ (t) introducing into the control laws, the dynamical equation (1) becomes where γ (t) + v (t) = u (t), (k = 1,…, r) are the total control laws.

Without loss of generality, assume |ψ 〉 = |ϕ 〉, 1 ≤ g ≤ N. In order to solve the problem of convergence for the degenerate cases, the perturbations γ (t) were introduced into the control laws. The basic idea is as follows: Denote the system with the internal Hamiltonian H 0, the control Hamiltonians H , (k = 1,…, r), and the control laws u (t) = γ (t) + v (t) as system 1, and the system with the internal Hamiltonian H 02 = H 0 + ∑ H γ (t), the control Hamiltonians H , (k = 1,…, r), and the control laws v (t) as system 2. All these two systems can be depicted by (37). Denote the eigenvalues and eigenstates of H 02 as λ 1,, λ 2,,…, λ and |ϕ 1,〉, |ϕ 2,〉,…, |ϕ 〉, respectively, which are the functions of the perturbations γ (t). Assume |ψ 〉 = |ϕ 〉. If one can design the perturbations γ (t) such that (i) ω ≠ ω , (l, m)≠(i, j), i, j, l, m ∈ {1,2,…, N}, where ω = λ − λ holds; (ii) for any |ϕ 〉 ≠ |ψ 〉, there exists at least a k ∈ {1,…, r} satisfying 〈ϕ |H |ψ 〉 ≠ 0, and select the specific Lyapunov function based on the state distance as Then according to Section 2.1.1, system 2 will converge to |ψ 〉. And when system 2 converge to |ψ 〉, if the perturbations γ (t) at |ψ 〉 are designed to equal zero, then system 2 will become system 1, and |ψ 〉 will become |ψ 〉. Then the convergence of system 1 to |ψ 〉 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 |ψ 〉 from the initial state |ψ 0〉, another one is γ (t) converges to 0. In order to make the introduced perturbations take effect to make system 1 in the non-degenerate case converge to |ψ 〉, the speed of γ (t), (k = 1,…, r) converging to 0 must be slower than the speed at which system 2 converges toward |ψ 〉. For convenience, the control system in the following section means system 1.

The existing design method of γ (t) is to design it to be a monotonically increasing functional of V(t), that is, where the function θ (·) satisfies θ (0) = 0, θ (s) > 0 and θ ′(s) > 0 for every s > 0 and s is the independent variable of the function θ (·).

The right side of (39) contains |ψ 〉 which is a functional of the perturbations γ (t), (k = 1,…, r). One can see that the relation between γ and the time t is defined by r equations and cannot be expressed by an explicit expression, so γ is the implicit function of the time t. From (38) and (39), one can also see that the Lyapunov function V 6 is the implicit function of the time t. The existence of γ (t) can be depicted by Lemma 9.

Lemma 9 (see [37, 40])

Let θ ∈ C (R +; [0, γ *]), k = 1,…, r, (γ * ∈ R, γ * > 0) satisfy θ (0) = 0, θ (s) > 0, θ ′(s) > 0 for every s > 0, ||θ ′|| < 1/C*, and C* = 1 + max⁡{||(∂ψ /∂γ )|(||; γ ∈ [0, γ *]}, where γ = γ (0). Then for every state |ψ〉 ∈ S 2 = {x ∈ C ; ||x|| = 1}, there exists a unique γ 1, γ 2,…, and γ with γ ∈ C (γ ∈ [0, γ *]), (k = 1,…, r) satisfying γ (t) = θ ((1/2)(1 − |〈ψ∣ψ 〉|2)), (k = 1,…, r), γ (|ψ 〉) = 0.

By (38) and Lemma 9, one can see that for every state |ψ〉 ∈ S 2 = {x ∈ C ; ||x|| = 1}; there also exists a unique V.

Assume The first order time derivative of the Lyapunov function (38) is In order to make , one can design γ (t) such that ||θ ′|| < 1/(rC*) to make 1/(1 + ℜ(〈ψ | ϕ 〉〈(∂ϕ /∂γ ) | ψ〉)∑ θ ′) > 0 hold and design v (t) as where K > 0, and v (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 u (t) = γ (t) + v (t), γ (t) designed in (39), Lemma 9 and ||θ ′|| < 1/(rC*), and v (t) designed in (42). Assume that the target state |ψ 〉 is an eigenstate of H0. If the control system satisfies (i) ω ≠ ω , (l, m)≠(i, j), i, j, l, m ∈ {1,2,…, N}, ω = λ − λ , where λ is an eigenvalue of (H 0 + ∑ H γ ), (ii) For any |ϕ 〉 ≠ |ψ 〉, there exists at least a k ∈ {1,…, r} satisfying 〈ϕ |H |ψ 〉 ≠ 0; then the largest invariant set is S 2∩E 7, where E 7 = {|ψ〉 = e |ψ 〉, θ ∈ R}. And the control system will converge toward |ψ 〉e , (θ ∈ R).

Conditions (i) and (ii) in Theorem 10 are associated with H 0, H , (k = 1,…, r) and γ (t). By designing appropriate γ (t), these two conditions can be satisfied in most cases. one can see that if one designs appropriate control laws u (t) = γ (t) + v (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 |ψ 〉 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 where u (t) = γ (t) + v (t), (k = 1,…, r) and ω is the control laws which need to design.

Consider the Lyapunov function based on the state error as Thus the implicit function perturbation γ (t) needs to be designed as where functions θ (·) satisfy θ (0) = 0, θ (s) > 0, θ ′(s) > 0 for every s > 0.

The existence of γ (t) can be depicted by Lemma 11.

Lemma 11 (see [38, 39])

Let θ ∈ C (R +; [0, γ *]), k = 1,…, r, (γ * ∈ R, γ * > 0) satisfy θ (0) = 0, θ (s) > 0, θ ′(s) > 0 for every s > 0, ||θ ′|| < 1/rC*, and C* = 1+max⁡{||(∂|ψ 〉/∂γ )|(||; γ ∈ [0, γ *]}. Then for every state |ψ〉 ∈ S 2 = {x ∈ C ; ||x|| = 1}, there exists a unique γ 1, γ 2,…, and γ with γ ∈ C (γ ∈ [0, γ *]),(k = 1,…, r) satisfying γ (|ψ〉) = θ ((1/2)〈ψ − ψ | ψ − ψ 〉), γ (|ψ 〉) = 0.

Assume (40) holds, the first order time derivative of the Lyapunov function (44) is In order to make , ω and v (t) need to be designed as where K > 0, c > 0, and y = f (x ), (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 γ (t) designed in (45) and Lemma 11, v (t) designed in (48), and ω designed in (47). Assume that the target state |ψ 〉 is an eigenstate of H 0. If the control system satisfies (i) ω ≠ ω , (l, m)≠(i, j), i, j, l, m ∈ {1,…, N}, ω = λ − λ , λ is the lth eigenvalue of (H 0 + ∑ H (γ + η )) corresponding to the eigenstate |ϕ l,〉; (ii) There is at least a k ∈ {1,…, r}, 〈ϕ |H |ψ 〉 ≠ 0, for |ϕ 〉 ≠ |ψ 〉; then the largest invariant set is S 2∩E 8, where E 8 = {|ψ〉 = e |ψ 〉, θ ∈ R}. And the control system will converge toward the equivalent state of the target state: |ψ 〉e , (θ ∈ R).

By designing appropriate perturbations γ (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 u (t) = γ (t) + v (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 |ψ 〉 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 γ (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(|ψ 〉) < 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(|ψ 〉) < V(|ψ other 〉), we can design all the perturbations γ (t) = 0 hold for k = 1,…, r only at |ψ 〉, that is, (1) γ (|ψ 〉) = 0, (k = 1,…, r), and (2) for |ψ〉 ≠ |ψ 〉, there exists at least one k such that γ (|ψ〉) ≠ 0. For sake of completing these two tasks, we can design γ (t) as a monotonically increasing functional of V(t): where C ≥ 0, and for k = 1,…, r, there exists at least a C > 0, and the function θ (·) satisfies θ (0) = 0, θ (s) > 0 and θ ′(s) > 0 for every s > 0, s is the independent variable of the function θ (·).

The specific Lyapunov function based on the average value of an imaginary mechanical quantity can be selected as: where P = f(γ 1(t),…, γ (t)) is a functional of γ (t) and positive definite.

The existence of γ (t) can be depicted by Lemma 13.

Lemma 13

If C = 0, γ (|ψ〉) = 0. Else if C > 0, θ ∈ 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*), C* = 1 + C, C = max⁡{||∂P /∂γ ||, γ ∈ [0, γ *]}, then for every |ψ〉 ∈ S 2, there is a unique γ ∈ C (γ ∈ [0, γ *]) satisfying γ (|ψ〉) = C · θ (〈ψ|P |ψ〉 − 〈ψ |P |ψ 〉).

For the sake of simplicity, γ (t) can be designed zero for some k, and other γ (t) equal, denoted by γ(t); that is, set where θ(·) = θ (·) = ⋯ = θ (·) and P are functionals of γ(t).

The time derivative of the selected Lyapunov function is In order to ensure , v (t) can be designed as where K is a constant and K > 0, and y = f (x ), (k = 1,2,…, r) are monotonic increasing functions through the coordinate origin of the plane x − y .

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 u (t) = γ (t) + v (t), (k = 1,…, r), where γ (t) is defined by Lemma 13 and (51) and v (t) is defined by (53). If the control system satisfies (i) ω ≠ ω , (l, m)≠(i, j), i, j, l, m ∈ {1,2,…, N}, ω = λ − λ , where λ is the lth eigenvalue of H 0 + ∑ H γ(t) corresponding to the eigenstate |ϕ 〉; (ii) for any i ≠ j, i, j ∈ {1,2,…, N}, there exists at least one k such that , where is the (l, m)th element of , U 1 = (|ϕ 1,〉,…, |ϕ 〉); (iii) [P , H 0 + ∑ H γ(t)] = 0; (iv) , where is the (l, l)th element of U 1 P U 1, then any trajectory will converge toward E 9 = {|ψ 〉|e |ϕ 〉; θ ∈ R, l ∈ {1,…, N}}.

From Theorem 14, one can see that if the target state |ψ 〉 is an eigenstate, |ψ 〉 is contained in E 9. 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 hold, where |ψ other 〉 represents any other state in the set E 9 except the target state. On the other hand, because holds, when holds for some time, we can design to make the state trajectory evolve but not stay in E 9 until |ψ 〉e is reached.

Next, let us analyze how to realize the convergence conditions in Theorem 14 and restriction defined by (54). By designing appropriate γ (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 P 1, P 2,…, P , and P can be designed as

To make condition (iv) satisfied, P ≠ P (∀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 designs P > P , (i = 1,…, N, P ≠ P ), then V 8(|ψ 〉) < V 8(|ψ other 〉) holds, where P is the eigenvalue of P corresponding to |ψ 〉.

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 laws u (t) = γ (t) + v (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 Section 3.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 η into the control laws. Thus the mechanical equation (37) becomes And the mechanical equation (43) in Section 3.1.2 becomes where η ∈ R.

The basic idea of solving the convergence to the target superposition state is to design η to make the target state |ψ 〉 be an eigenstate of H 0′ = H 0 + ∑ H η · H 0′ 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 η , (H 0 + ∑ H η )|ψ 〉 = λ ′|ψ 〉 can be satisfied in most cases, where λ ′ is the eigenvalue of H 0 + ∑ H η corresponding to |ψ 〉. 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 H 0 into H 0′.

3.3. Target Mixed State Which Commutes with the Internal Hamiltonian

Consider the N-level closed quantum system governed by the following quantum Liouville equation: where γ (t) + v (t) = u (t) are the total control laws.

The design ideas are similar to those of Section 3.1. Consider the Lyapunov function based on the imaginary mechanical quantity as γ (t) can be designed as

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 /∂γ||, γ ∈ [0, γ*]}, then for every ρ, there is a unique γ ∈ C (γ ∈ [0, γ*]) satisfying γ(ρ) = θ(tr⁡(P ρ) − tr⁡(P ρ )).

According to the Lyapunov stabilty Theorem, v (t) can be designed as where K is a constant and K > 0, and y = f (x ), (k = 1,2,…, r) are monotonic increasing functions through the coordinate origin of the plane x − y .

Based on the LaSalle invariance principle, the convergence of the control system can be depicted by Theorem 17.

Theorem 17

Consider the control system depicted by (58) with control laws u (t) = γ (t) + v (t), where γ (t) is defined by Lemma 16 and (60), and v (t) is defined by (61). If the control system satisfies (i) ω ≠ ω , (l, m)≠(i, j), i, j, l, m ∈ {1,2,…, N}, ω = λ − λ , where λ is the lth eigenvalue of (H 0 + ∑ H γ(t)) corresponding to the eigenstate |ϕ l,〉; (ii) ∀j ≠ l, for k = 1,…, r, there exists at least a , where is the (j, l)th element of with U 2 = (|ϕ 1,〉,…, |ϕ 〉); (iii) [P , H 0 + ∑ H γ(t)] = 0,1 ≤ k 1,…, k ≤ r; (iv) for any l ≠ j, (1 ≤ l, j ≤ N), holds, where is the (l, l)th element of , then the control system will converge toward E 10 = {ρ |(U 2 † ρ U 2) = 0, γ = γ(ρ ), t 0 ∈ R}.

For the target state ρ which commutes with the internal Hamiltonian; that is, [ρ , H 0] = 0, ρ is contained in E 10. In order to make the system converge to the target state ρ which commutes with H 0 from any initial stat; on the one hand, P needs to be designed to make hold, where ρ other represents any other state in the set E 10 except the target state. On the other hand, γ was designed as when holds for some time to make the state trajectory evolve but not stay in E 10 until ρ is reached.

For the satisfaction of conditions (i)–(iv) in Theorem 17, 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 laws u (t) = γ (t) + v (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 Section 3.2. The difference is to design η to make [ρ , H 0′] = 0. If the number of the control Haimltonians r is large enough, by designing appropriate η , [ρ , H 0′] = 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 H 0 into H 0′.

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.