Stabilization of a delayed quantum system:
the photon box casestudy^{1}
Abstract
We study a feedback scheme to stabilize an arbitrary photon number state in a microwave cavity. The quantum nondemolition measurement of the cavity state allows a nondeterministic preparation of Fock states. Here, by the mean of a controlled field injection, we desire to make this preparation process deterministic. The system evolves through a discretetime Markov process and we design the feedback law applying Lyapunov techniques. Also, in our feedback design we take into account an unavoidable pure delay and we compensate it by a stochastic version of a Smith predictor. After illustrating the efficiency of the proposed feedback law through simulations, we provide a rigorous proof of the global stability of the closedloop system based on tools from stochastic stability analysis. A brief study of the Lyapunov exponents of the linearized system around the target state gives a strong indication of the robustness of the method.
1 Introduction
In the aim of achieving a robust processing of quantum information, one of the main tasks is to prepare and to protect various quantum states. Through the last 15 years, the application of quantum feedback paradigms has been investigated by many physicists [21, 19, 5, 10, 16] as a possible solution for this robust preparation. However, most (if not all) of these efforts have remained at a theoretical level and have not been able to be give rise to successful experiments. This is essentially due to the necessity of simulating, in parallel to the system, a quantum filter [1] providing an estimate of the state of the system based on the historic of quantum jumps induced by the measurement process. Indeed, it is, in general, difficult to perform such simulations in real time. In this paper, we consider a prototype of physical systems, the photonbox, where we actually have the time to perform these computations in real time (see [6] for a detailed description of this cavity quantum electrodynamics system).
Taking into account the measurementinduced quantum projection postulate, the most practical measurement protocols in the aim of feedback control are the quantum nondemolition (QND) measurements [2, 18, 20]. These are the measurements which preserve the value of the measured observable. Indeed, by considering a welldesigned QND measurement process where the quantum state to be prepared is an eigenstate of the measurement operator, the measurement process, not only, is not an obstacle for the state preparation but can even help by adding some controllability.
In [4, 9, 8] QND measures are exploited to detect and/or produce highly nonclassical states of light trapped in a superconducting cavity (see [11, chapter 5] for a description of such QED systems and [3] for detailed physical models with QND measures of light using atoms). For such experimental setups, we detail and analyze here a feedback scheme that stabilizes the cavity field towards any photonnumber states (Fock states). Such states are strongly nonclassical since their photon numbers are perfectly defined. The control corresponds to a coherent lightpulse injected inside the cavity between atom passages. The overall structure of the proposed feedback scheme is inspired by [7] using a quantum adaptation of the observer/controller structure widely used for classical systems (see, e.g., [12, chapter 4]). As the measurementinduced quantum jumps and the controlled field injection happen in a discreteintime manner, the observer part of the proposed feedback scheme consists in a discretetime quantum filter. Indeed, the discreteness of the measurement process provides us a first prototype of quantum systems where we, actually, have enough time to perform the quantum filtering and to compute the measurementbased feedback law to be applied as the controller.
From a mathematical modeling point of view, the quantum filter evolves through a discretetime Markov chain. The estimated state is used in a statefeedback, based on a Lyapunov design. Indeed, by considering a natural candidate for the Lyapunov function, we propose a feedback law which ensures the decrease of its expectation over the Markov process. Therefore, the value of the considered Lyapunov function over the Markov chain defines a supermartingale. The convergence analysis of the closedloop system is, therefore, based on some rather classical tools from stochastic stability analysis [13].
One of the particular features of the system considered in this paper corresponds to a nonnegligible delay in the feedback process. In fact, in the experimental setup considered through this paper, we have to take into account a delay of steps between the measurement process and the feedback injection. Indeed, there are, constantly, atoms flying between the photon box (the cavity) to be controlled and the atomdetector (typically ). Therefore, in our feedback design, we do not have access to the measurement results for the last atoms. Through this paper, we propose an adaptation of the quantum filter, based on a stochastic version of the Smith predictor [17], which takes into account this delay by predicting the actual state of the system without having access to the result of last detections.
In the next section, we describe briefly the physical system and the associated quantum MonteCarlo model. In Section 3, we consider the dynamics of the openloop system. We will prove, through theorem 1 that the QND measurement process, without any additional controlled injection, allows a nondeterministic preparation of the Fock states. Indeed, we will see that the associated Markov chain converges, necessarily, towards a Fock state and that the probability of converging towards a fixed Fock state is given by its population over the initial state. Also, through proposition 1, we will show that the linearized openloop system around a fixed Fock state admits strictly negative Lyapunov exponents (see Appendix B for a definition of the Lyapunov exponent).
In Section 4, we propose a Lyapunovbased feedback design allowing to stabilize globally the delayed closedloop system around a desired Fock state. The theorem 2 proves the almost sure convergence of the trajectories of the closedloop system towards the target Fock state. Also, through proposition 2, we will prove that the linearized closedloop system around the target Fock state admits strictly negative Lyapunov exponents.
Finally in Section 5, we propose a brief discussion on the considered quantum filter and by proving a rather general separation principle (theorem 3), we will show a semiglobal robustness with respect to the knowledge of the initial state of the system. Also, through a brief analysis of the linearized systemobserver around the target Fock state and applying the propositions 1 and 2, we show that its largest Lyapunov exponent is also strictly negative (proposition 3).
2 A discretetime Markov process
As illustrated by Figure 1, the system consists in a highQ microwave cavity, a box producing Rydberg atoms, and two lowQ Ramsey cavities, an atom detector and a microwave source. The dynamics model is discret in time and relies on quantum MonteCarlo trajectories (see [11, chapter 4]). Each timestep indexed by the integer corresponds to atom number coming from , submitted then to a first Ramsey pulse in , crossing the cavity and being entangled with it, submitted to a second pulse in and finally being measured in . The state of the cavity is associated to a quantized mode. The control corresponds to a coherent displacement of amplitude that is applied via the microwave source between two atom passages.
In this paper we consider a finite dimensional approximation of this quantized mode and take a truncation to photons. Thus the cavity space is approximated by the Hilbert space . It admits as orthonormal basis. Each basis vector corresponds to a pure state, called Fock state, where the cavity has exactly photons, . In this Fockstates basis the number operator N corresponds to the diagonal matrix
The annihilation operator truncated to photons is denoted by a. It corresponds to the upper diagonal matrix filled with :
The truncated creation operator denoted by is the Hermitian conjugate of a. Notice that we still have , but truncation does not preserve the usual commutation that is only valid when .
Just after the measurement of the atom number , the state of the cavity is described by the density matrix belonging to the following set of welldefined density matrices:
(1) 
The random evolution of this state can be modeled through a discretetime Markov process that will be described bellow (see [6] and the references therein explaining the physical modeling assumptions).
Let us denote by the control at step . Then , the cavity state after measurement of atom is given by
(2) 
where,

, , with operators and ( constant parameters). For any we set .

where the unitary displacement operator is given by In openloop, , (identity operator) and . Notice that .

is a random variable taking the value when the atom is detected in (resp. when the atom is detected in ) with probability
(3) 
The control elaborated at step , , is subject to a delay of steps, being the number of flying atoms between the cavity and the detector .
We will assume through out the paper that the parameters , are chosen in order to have , invertible and such that the spectrum of and are not degenerate. This implies that the nonlinear operators and are well defined for all and that and belongs also to the state space defined by (1). Notice that and commute, are diagonal in the Fock basis and satisfy . The Kraus map associated to this Markov process is given by:
(4) 
It corresponds to the expectation value of knowing and :
(5) 
3 Open loop dynamics
3.1 Simulations
We consider in this section the following dynamics
(6) 
obtained from (2) when . Figure 2 corresponds to 100 realizations of this Markov process with photons, and . For each realization, is initialized to the the same coherent state with as mean photon number. We observe that either tends to or . Since the ensemble average curve is almost constant, the proportion of trajectories for which tends to is given approximatively by .
3.2 Global convergence analysis
The following theorem underlies the observations made for simulations of Figure 2
Theorem 1.
Proof.
Let us prove that is a martingale. Set . We have
Since commutes with and and we have .
Considering the following function:
where Notice that is convexe, on and satisfies
(7) 
The function is increasing and convex and is a martingale. Thus is submartingale. Since
we have
Then (7), together with
yields to
Thus we recover that is a submartingale, . We have also shown that implies that either or (assumption and invertible is used here).
We apply now the invariance theorem established by Kushner [13] (recalled in the Appendix A) for the Markov process and the submartingale . This theorem implies that the Markov process converges in probability to the largest invariant subset of
But the set is invariant. It remains thus to characterized the largest invariant subset denoted by and included in .
Take . Invariance means that and belong to (the fact that and are invertible ensures that probabilities to jump with or are strictly positive for any ). Consequently . This means that . By CauchySchwartz inequality,
with equality if, and only if, and are colinear. being nondegenerate, is necessarily a projector over an eigenstate of , i.e., for some . Since , and thus is reduced to . Therefore the only possibilities for the limit set are or and
The convergence in probability together with the fact that is a positive bounded () random process implies the convergence in expectation. Indeed
where for the last inequality, we have applied the convergence in probability of towards . As the above inequality is valid for any , we have
Furthermore, by the first part of the Theorem, we know that is a bounded martingale and therefore by the Doob’s first martingale convergence theorem (see the Theorem 4 of the Appendix A), converges almost surely towards a random variable . This implies that converges almost surely towards the random variable . We apply now the dominated convergence theorem
This implies that vanishes almost surely and therefore
As we can repeat this same analysis for any choice of , converges almost surely to the set of of Fock states
which ends the proof of the second part.
We have shown that the probability measure associated to the random variable converges to the probability measure
where denotes the Dirac distribution at and is the probability of convergence towards . In particular, we have
But is a martingale and . Thus
which ends the proof of the third and last part.
∎
3.3 Local convergence rate
According to theorem 1, the limit set of the Markov process (6) is the discrete set of Fock states . We investigate here the local convergence rate around one of these Fock states denoted by for some .
Since , we can develop the dynamics (6) around the fixed point . We write with small, Hermitian and with zero trace. Keeping only the first order terms in (6), we have
Thus the linearized Markov process around the fixed point reads
(8) 
where the random matrices are given by :

with probability

with probability .
The following proposition shows that the convergence of the linearized dynamics is exponential (a crucial robustness indication).
Proposition 1.
Consider the linear Markov chain (8) of state belonging to the set of Hermitian matrices with zero trace. Then the largest Lyapunov exponent is given by ()
and is strictly negative: .
Proof.
Set for any . Since , we exclude here the case because . Since and are diagonal matrices, we have
(9) 
where (resp. ) with probability (resp. ) and where and .
Denote by the Lyapunov exponent of (9) for . By the law of large numbers, we know that converges almost surely towards
Thus, we have
The function
increases strictly from to when goes from to and decreases strictly from to when goes from to . Since , . Denote by for :
Since , we have and is strictly negative. ∎
4 Feedback stabilization with delays
4.1 Feedback scheme and closedloop simulations
Through out this section we assume that we have access at each step to the cavity state . The goal is to design a causal feedback law that stabilizes globally the Markov chain (2) towards a goal Fock state with photon(s), . To be consistent with truncation to photons, has to be far from (typically with in the simulations below).
The feedback is based on the fact that, in openloop when , is a martingale. When , [15] proves global almost sure convergence of the following feedback law
(10) 
for any when are small enough. This feedback law ensures that is a submartingale.
When , we cannot set since will depend on and the feedback law is not causal. In [6], this feedback law is made causal by replacing by its expectation value (average prediction) knowing and the past controls , …, :
where the Kraus map is defined by (4).
We will thus consider here the following causal feedback based on an average compensation of the delay
(11) 
with
The closedloop system, i.e. Markov chain (2) with the causal feedback (11) is still a Markov chain but with as state at step . More precisely, denote by this state where stands for the control delayed steps. Then the state form of the closedloop dynamics reads
(12) 
where the control law defined by (11) corresponds to a static state feedback since
(13) 
Notice that .
Simulations displayed on Figures 3 and 4 correspond to 100 realizations of the above closedloop systems with and . The goal state contains photons and , and are those used for the openloop simulations of Figure 2. Each realization starts with the same coherent state and . The feedback parameters appearing in (11) are as follows:
This simulations illustrate the influence of the delay on the average convergence speed: the longer the delay is the slower convergence speed becomes.
Remark 1.
The choice of the feedback law whenever might seem complicated for realtime simulation issues. However, this choice is only technical. Actually, any nonzero constant feedback law will seems to achieve the task here (see for instance the simulations of [6]). However, the convergence proof for such simplified control scheme is more complicated and not considered in this paper.