Invariant measure for stochastic Schr\"odinger equations

TL;DR

This paper investigates the invariant measures of quantum trajectories described by stochastic Schrödinger equations, proving uniqueness under certain conditions and exponential convergence, with explicit examples provided.

Contribution

It establishes conditions for the uniqueness and exponential convergence of invariant measures for quantum trajectories governed by stochastic Schrödinger equations.

Findings

Invariant measure is unique under ergodicity and purification conditions.

Quantum trajectories converge exponentially fast to the invariant measure.

Explicit expressions for invariant measures are derived in specific examples.

Abstract

Quantum trajectories are Markov processes that describe the time-evolution of a quantum system undergoing continuous indirect measurement. Mathematically, they are defined as solutions of the so-called "Stochastic Schr\"odinger Equations", which are nonlinear stochastic differential equations driven by Poisson and Wiener processes. This paper is devoted to the study of the invariant measures of quantum trajectories. Particularly, we prove that the invariant measure is unique under an ergodicity condition on the mean time evolution, and a "purification" condition on the generator of the evolution. We further show that quantum trajectories converge in law exponentially fast towards this invariant measure. We illustrate our results with examples where we can derive explicit expressions for the invariant measure.

Equations

$${\mathcal{D}}_k=\{\rho \in {\mathrm{M}}_k(\mathbb{C})\text{mbox}s.t.\rho \ge 0,\mathrm{tr}\rho =1\}.$$

$$\mathrm{d}{\bar{\rho }}_t=\mathcal{L}({\bar{\rho }}_t)\mathrm{d}t,{\bar{\rho }}_0\in {\mathcal{D}}_k,$$

$$\mathcal{L}:\rho\mapsto-\mathrm{i}[H,\rho]+\sum_{i\in I}\big{(}V_{i}\rho V_{i}^{*}-\genfrac{}{}{}{1}{1}{2}\{V_{i}^{*}V_{i},\rho\}\big{)},$$

$$\mathrm{d}{\rho }_t=$$

$$\displaystyle+\quad\sum_{i\in I_{b}}\Big{(}L_{i}\rho_{t-}+\rho_{t-}L_{i}^{*}-\operatorname{tr}\big{(}\rho_{t-}(L_{i}+L_{i}^{*})\big{)}\rho_{t-}\Big{)}\,\mathrm{d}B_{i}(t)$$

$$\displaystyle+\quad\sum_{j\in I_{p}}\Big{(}\frac{C_{j}\rho_{t-}C_{j}^{*}}{\operatorname{tr}(C_{j}\rho_{t-}C_{j}^{*})}-\rho_{t-}\Big{)}\Big{(}\mathrm{d}N_{j}(t)-\operatorname{tr}(C_{j}\rho_{t-}C_{j}^{*})\,\mathrm{d}t\Big{)},$$

$$\mathrm{d}{x}_t$$

$$\displaystyle\quad+\sum_{i\in I_{b}}\big{(}L_{i}-\genfrac{}{}{}{1}{1}{2}v_{i}(t-)\,\mathrm{Id}\big{)}x_{t-}\,\mathrm{d}B_{i}(t)$$

$$\displaystyle\quad+\sum_{j\in I_{p}}\Big{(}\frac{C_{j}}{\sqrt{n_{j}(t-)}}-\mathrm{Id}\Big{)}x_{t-}\,\mathrm{d}N_{j}(t),$$

$$D(x_{t-})=-\big{(}\mathrm{i}H+\frac{1}{2}\sum_{i\in I_{b}}L_{i}^{*}L_{i}+\frac{1}{2}\sum_{j\in I_{p}}C_{j}^{*}C_{j}\big{)}+\frac{1}{2}\sum_{i\in I_{b}}v_{i}(t-)\big{(}L_{i}-\genfrac{}{}{}{1}{1}{4}\,v_{i}(t-)\,\mathrm{Id}\big{)}+\frac{1}{2}\sum_{j\in I_{p}}n_{j}(t-),$$

$$v_i(t-)=⟨x_{t-},(L_i+L_i^{\ast })x_{t-}⟩,n_j(t-)=⟨x_{t-},C_j^{\ast }{C}_j{x}_{t-}⟩=\Vert C_j{x}_{t-}{\Vert }^2.$$

$$\underset{t\to \infty }{\lim }\underset{\hat{y}\in \mathrm{P}{\mathbb{C}}^k}{\inf }\Vert {\rho }_t-{\pi }_{\hat{y}}\Vert =0\text{almost surely.}$$

$$W_1({\mu }_t,{\mu }_{\mathrm{inv}})\le C{\mathrm{e}}^{-\lambda t}$$

$$A\cdot \hat{x}=\widehat{Ax}.$$

$$d(\hat{x},\hat{y})=\sqrt{1-|⟨x,y⟩{|}^2},$$

$$\mathrm{d}S_{t}^{s}=\big{(}K+{\textstyle{\frac{\#I_{p}}{2}}}\,{\mathrm{Id}}\big{)}S_{t-}^{s}\,\mathrm{d}t+\sum_{i\in I_{b}}L_{i}S_{t-}^{s}\,\mathrm{d}W_{i}(t)+\sum_{j\in I_{p}}(C_{j}-\mathrm{Id})S_{t-}^{s}\,\mathrm{d}N_{j}(t),\qquad S_{s}^{s}=\mathrm{Id}$$

$$K=-\mathrm{i}H-\frac{1}{2}(\sum _{i\in I_b}{L}_i^{\ast }{L}_i+\sum _{j\in I_p}{C}_j^{\ast }{C}_j).$$

$$S_t^s{S}_s^r=S_t^r.$$

$$Z_t^{\rho }=\mathrm{tr}(S_t^{\ast }{S}_t\rho ),$$

$${\rho }_t=\frac{S_t\rho S_t^{\ast }}{\mathrm{tr}(S_t\rho S_t^{\ast })}$$

$$\mathrm{d}Z_{t}^{\rho}=Z_{t-}^{\rho}\Big{(}\sum_{i\in I_{b}}\operatorname{tr}\big{(}(L_{i}+L_{i}^{*})\rho_{t-}\big{)}\mathrm{d}W_{i}(t)+\sum_{j\in I_{p}}\big{(}\operatorname{tr}(C_{j}^{*}C_{j}\rho_{t-})-1\big{)}\big{(}\mathrm{d}N_{j}(t)-\mathrm{d}t\big{)}\Big{)},\quad Z_{0}^{\rho}=1.$$

$$\mathrm{d}{\mathbb{P}}_t^{\rho }=Z_t^{\rho }\mathrm{d}\mathbb{P}{|}_{{\mathcal{F}}_t}.$$

$$B_{i}^{\rho}(t)=W_{i}(t)-\int_{0}^{t}\operatorname{tr}\big{(}(L_{i}+L_{i}^{*})\rho_{s-}\big{)}\,\mathrm{d}s.$$

$$\mathrm{d}{\rho }_t$$

$$\displaystyle\quad+\sum_{i\in I_{b}}\Big{(}L_{i}\rho_{t-}+\rho_{t-}L_{i}^{*}-\operatorname{tr}\big{(}\rho_{t-}(L_{i}+L_{i}^{*})\big{)}\rho_{t-}\Big{)}\mathrm{d}B_{i}^{\rho}(t)$$

$$\displaystyle\quad+\sum_{j\in I_{p}}\Big{(}\frac{C_{j}\rho_{t-}C_{j}^{*}}{\operatorname{tr}(C_{j}\rho_{t-}C_{j}^{*})}-\rho_{t-}\Big{)}\Big{(}\mathrm{d}N_{j}(t)-\operatorname{tr}(C_{j}\rho_{t-}C_{j}^{*})\,\mathrm{d}t\Big{)}.$$

$$\mathbb{E}(S_t\rho S_t^{\ast })=\mathbb{E}({\rho }_t{Z}_t^{\rho })={\mathbb{E}}^{\rho }({\rho }_t)={\mathrm{e}}^{t\mathcal{L}}(\rho ).$$

$${\mathbb{Q}}_{\mu }(E\times A)=\int {\mathbb{P}}^{{\pi }_{\hat{x}}}(E)\text{mathbbm}{1}_{\hat{x}\in A}\mathrm{d}\mu (\hat{x}).$$

$$\mathbb{Q}_{\mu}\big{(}\{S_{t}x=0\}\big{)}=0$$

$${\hat{x}}_t=S_t\cdot x$$

$${\rho }_{\mu }:={\mathbb{E}}_{\mu }({\pi }_{\hat{x}}).$$

$${\mathbb{Q}}_{\mu }(E\times \mathrm{P}{\mathbb{C}}^k)={\mathbb{P}}^{{\rho }_{\mu }}(E).$$

$${\mathbb{Q}}_{\mu }(E\times \mathrm{P}{\mathbb{C}}^k)=\int {\mathbb{P}}^{{\pi }_{\hat{x}}}(E)\mathrm{d}\mu (x)=\int {\int }_E\mathrm{tr}(S_t^{\ast }{S}_t{\pi }_{\hat{x}})\mathrm{d}\mathbb{P}\mathrm{d}\mu (x).$$

Full text

1. 1 Introduction
2. 2 Construction of the model
3. 3 Invariant measure and exponential convergence in Wasserstein distance
4. 4 Set of invariant measures under (Pur)
5. 5 (Pur) is not necessary for purification
6. 6 Examples

1. Introduction

Under a Markov approximation, the evolution of an open quantum system $\mathcal{S}$ in interaction with an environment $\mathcal{E}$ is described by the Gorini–Kossakowski–Sudarshan–Lindblad Master (GKSL) equation [23, 30]. More precisely, assuming that the system is described by the Hilbert space $\mathbb{C}^{k}$, the set of its states is defined as the set $\mathcal{D}_{k}$ of density matrices, i.e. positive semidefinite matrices with trace one:

[TABLE]

The evolution $t\in\mathbb{R}_{+}\mapsto\bar{\rho}_{t}\in\mathcal{D}_{k}$ of states of the system is then determined by the GKSL equation (also called quantum master equation):

[TABLE]

where $\mathcal{L}$ is a bounded linear operator on $\mathrm{M}_{k}(\mathbb{C})$ of the form

[TABLE]

with $I$ a finite set, $H\in\mathrm{M}_{k}(\mathbb{C})$ self adjoint, and $V_{i}\in\mathrm{M}_{k}(\mathbb{C})$ for each $i\in I$ ($[\cdot,\cdot]$ and $\{\cdot,\cdot\}$ are respectively the commutator and anticommutator). Such an $\mathcal{L}$ is called a Lindblad operator.

Since $\mathcal{L}$ is linear, $t\mapsto\bar{\rho}_{t}$ is given by $\bar{\rho}_{t}=\mathrm{e}^{t\mathcal{L}}(\bar{\rho}_{0})$. The flow is therefore a semigroup $(\mathrm{e}^{t\mathcal{L}})_{t}$, which consists of completely positive, trace-preserving maps (see [38]). In particular, $\mathcal{L}$ is the generator of a semigroup of contractions, thus $\operatorname{spec}\mathcal{L}\subset\{\lambda\in\mathbb{C}\ \mbox{ s.t.}\ \operatorname{Re}\lambda\leq 0\}$. Since $\mathrm{e}^{t\mathcal{L}}$ is trace preserving, $0\in\operatorname{spec}\mathcal{L}$. The following assumption is equivalent to the simplicity of the eigenvalue [math] [38, Proposition 7.6]:

**($\mathcal{L}$-erg)****: **

There exists a unique non zero minimal orthogonal projection $\pi$ such that $\mathcal{L}(\pi\mathrm{M}_{k}(\mathbb{C})\pi)\subset\pi\mathrm{M}_{k}(\mathbb{C})\pi$.

Assumption ($\mathcal{L}$-erg) implies directly that there exists a unique $\rho_{\mathrm{inv}}\in\mathcal{D}_{k}$ such that $\mathcal{L}\rho_{\mathrm{inv}}=0$. Moreover, one can show that ($\mathcal{L}$-erg) implies the existence of $\lambda>0$ such that for any $\rho\in\mathcal{D}_{k}$, $\mathrm{e}^{t\mathcal{L}}(\rho)=\rho_{\mathrm{inv}}+O(\mathrm{e}^{-\lambda t})$ (see [38, Proposition 7.5]).

The above framework generalizes that of continuous-time Markov semigroups on a finite number of sites: density matrices $\rho$ over $\mathbb{C}^{k}$ generalize probability distributions over $k$ classical states, while Lindbladians $\mathcal{L}$ generalize generators of Markov jump processes. In Section 6.4, we show how a classical finite state Markov jump process can be encoded in the present formalism.

The family $(\bar{\rho}_{t})_{t}$ describes the reduced evolution of the system $\mathcal{S}$ when coupled to an environment $\mathcal{E}$ in a conservative manner. This evolution can be derived by considering the full Hamiltonian of $\mathcal{S}+\mathcal{E}$ in relevant limiting regimes, e.g. the weak coupling or fast repeated interactions regimes, and tracing out the environment degrees of freedom (see [17, 18] and [1] respectively). It can also be described by a stochastic unravelling, i.e. a stochastic process $(\rho_{t})_{t}$ with values in $\mathcal{D}_{k}$ such that the expectation $\overline{\rho}_{t}$ of $\rho_{t}$ satisfies (1.1); this method was developed in [4, 5, 6]. One possible choice of a stochastic unravelling is described by the following stochastic differential equation (SDE), called a stochastic master equation:

[TABLE]

where

• •

$I=I_{b}\cup I_{p}$ is a partition of $I$ such that $L_{i}=V_{i}$ for $i\in I_{b}$ and $C_{j}=V_{j}$ for $j\in I_{p}$,

• •

each $B_{i}$ is a Brownian motion,

• •

each $N_{j}$ is a Poisson process of intensity $t\mapsto\int_{0}^{t}\operatorname{tr}(C_{j}\rho_{s-}C_{j}^{*})\mathrm{d}s$.

Remark 1*.*

The processes \big{(}B_{j}(t)\big{)}_{t} and $\left(N_{j}(t)-\int_{0}^{t}\operatorname{tr}(C_{j}\rho_{s-}C_{j}^{*})\mathrm{d}s\right)_{t}$ are actually martingales. Then assuming that (1.3) accepts a solution, it is easy to check that for any $t\geq 0$, the expectation of $\rho_{t}$ is equal to $\bar{\rho}_{t}$ whenever $\rho_{0}=\bar{\rho}_{0}$.

Proper definitions of these Poisson processes and proofs of existence and uniqueness of the solution to (1.3) can be found in [5, 6, 33, 34, 35]. A solution $(\rho_{t})_{t}$ of Equation (1.3) is called a quantum trajectory.

Equations of the form (1.3) are used to model experiments in quantum optics (photo-detection, heterodyne or homodyne interferometry), particularly for measurement and control (see [15, 24, 37]). They were also introduced as stochastic collapse models (see [19, 22]) and as numerical tools to compute $\overline{\rho}_{t}$ (see [16]). Here we are interested in the fact that they model the evolution of the system $\mathcal{S}$ when continuous measurements are done on the environment $\mathcal{E}$. This can be shown starting from quantum stochastic differential equations using quantum filtering [3, 10, 13, 21, 25]. An approach using the notion of a priori and a posteriori states has been also developed using “classical” stochastic calculus (see the reference book by Barchielli and Gregoratti [5], and references therein). Continuous-time limits of discrete-time models can also be considered, see [33, 34, 35].

Equation (1.3) has the property that if $\rho_{0}$ is an extreme point of $\mathcal{D}_{k}$, then ${\rho}_{t}$ is almost surely an extreme point of $\mathcal{D}_{k}$ for any $t\in\mathbb{R}_{+}$. Since we will extensively use this property, let us make it explicit. The extreme points of $\mathcal{D}_{k}$ are the rank-one orthogonal projectors of $\mathbb{C}^{k}$; for any $x\in\mathbb{C}^{k}\setminus\{0\}$, let $\hat{x}$ be its equivalence class in $\mathrm{P}\mathbb{C}^{k}$, the projective space of $\mathbb{C}^{k}$. For $\hat{x}\in\mathrm{P}\mathbb{C}^{k}$, let $\pi_{\hat{x}}$ be the orthogonal projector onto $\mathbb{C}x$. Then $\hat{x}\in\mathrm{P}\mathbb{C}^{k}\mapsto\pi_{\hat{x}}$ is a bijective map from $\mathrm{P}\mathbb{C}^{k}$ to the set of extreme points. Assume now that $\rho_{0}=\pi_{\hat{x}_{0}}$ for some $\hat{x}_{0}\in\mathrm{P}\mathbb{C}^{k}$. Then it is easy to check that $\rho_{t}=\pi_{\hat{x}_{t}}$ almost surely for any $t\in\mathbb{R}_{+}$, with $t\mapsto x_{t}$ the unique solution to the following SDE, called a stochastic Schrödinger equation:

[TABLE]

for $x_{0}\in\hat{x}_{0}$ of norm one, where the operator $D(x_{t-})$ is defined as

[TABLE]

with

[TABLE]

The brackets $\langle\cdot,\cdot\rangle$ denote the scalar product in $\mathbb{C}^{k}$. Without possible confusion, a solution $(x_{t})_{t}$ will be also called a quantum trajectory. Remark that $\|x_{0}\|=1$ implies $\|x_{t}\|=1$ almost surely for any $t\in\mathbb{R}_{+}$; remark also that the numerical computation of $\rho_{t}$ involves only multiplications of matrices with vectors and not multiplications of matrices (this is the motivation for the use of quantum trajectories as numerical tools mentioned above).

In the physics literature, extreme points of $\mathcal{D}_{k}$ are called pure states. In particular, the preceding paragraph shows that the evolution dictated by Eq. (1.3) preserves pure states. It actually has also the property that quantum trajectories (solution of (1.3)) tend to “purify”. This has been formalized by Maassen and Kümmerer in [31] for discrete-time quantum trajectories, and extended to the continuous-time case by Barchielli and Paganoni in [7]. Purification is related to the following assumption (here $A\propto B$ means there exists $\lambda\in\mathbb{C}$ such that $A=\lambda B$ or $\lambda A=B$. Particularly we allow for $\lambda=0$).

**(Pur)****: **

Any non zero orthogonal projector $\pi$ such that for all $i\in I_{b}$, $\pi(L_{i}+L_{i}^{*})\pi\propto\pi$ and for all $j\in I_{p}$, $\pi C_{j}^{*}C_{j}\pi\propto\pi$ has rank one.

As shown in [7], (Pur) implies that for any $\rho_{0}\in\mathcal{D}_{k}$

[TABLE]

The main goal of this article is to show how the exponential convergence of the solution $(\overline{\rho}_{t})_{t}$ of Eq. (1.1), induced by ($\mathcal{L}$-erg), translates for its stochastic unravelling $(\rho_{t})_{t}$ solution of Eq. (1.3). We prove uniqueness of the invariant measure for continuous-time quantum trajectories assuming both ($\mathcal{L}$-erg) and (Pur)****. From (1.5), under these assumptions, the invariant measure will be concentrated on pure states, so we only need to prove uniqueness of the invariant measure for $(\hat{x}_{t})_{t}$ equivalence class of $(x_{t})_{t}$ solution of (1.4) (since $\mathrm{P}\mathbb{C}^{k}$ is compact and the involved process is Feller, the existence of an invariant measure is obvious). The difficulty of this proof lies in the failure of usual techniques like $\varphi$-irreducibility. Note that this question has already been partially addressed in the literature: essentially, only diffusive equations have been considered, i.e. equations for which Eq. (1.3) or (1.4) contain no jump term (in our notation, $I_{p}=\emptyset$). The results of [7] were, to our knowledge, the most advanced ones so far. In that article, algebraic conditions on the vector fields describing the stochastic differential equation are imposed to obtain the uniqueness of the invariant measure. This allows the authors to apply directly standard results from the analysis of stochastic differential equations. Unfortunately their assumptions are hard to check for a given family of matrices $(L_{i})_{i\in I_{b}}$.

The main result of the present paper is the following theorem.

Theorem 1.1**.**

Assume that (Pur) and ($\mathcal{L}$-erg) hold. Then the Markov process $(\hat{x}_{t})_{t}$ has a unique invariant probability measure $\mu_{\mathrm{inv}}$, and there exist $C>0$ and $\lambda>0$ such that for any initial distribution $\mu$ of $\hat{x}_{0}$ over $\mathrm{P}\mathbb{C}^{k}$, for all $t\geq 0$, the distribution $\mu_{t}$ of $\hat{x}_{t}$ satisfies

[TABLE]

where $W_{1}$ is the Wasserstein distance of order $1$.

This theorem is more general than previous similar results in different ways. First, we consider stochastic Schrödinger equations involving both Poisson and Wiener processes. Second, our assumptions are standard for quantum trajectories and are easy to check for a given family of operators \big{(}H,(L_{i})_{i\in I_{b}},(C_{j})_{j\in I_{p}}\big{)}. Last, we prove an exponential convergence towards the invariant measure. As a byproduct, we also provide a simple proof of the purification expressed in Eq. (1.5) (see Proposition 2.5). To complete the picture, assuming only (Pur), we show that ($\mathcal{L}$-erg) is necessary. We also provide a complete characterization of the set of invariant measures of $(\hat{x}_{t})$ whenever ($\mathcal{L}$-erg) does not hold (see Proposition 4.2). Arguments in Sections 3 and 4 are adaptations of [11], where similar results for discrete-time quantum trajectories are considered.

The paper is structured as follows. In Section 2, we give a precise description of the model of quantum trajectories with a proper definition of the underlying probability space. In particular, we introduce a new martingale which is central to our proofs. In Section 3, we prove Theorem 1.1. In Section 4 we derive the full set of invariant measures assuming only (Pur). In Section 5 we show that (Pur) is not necessary even if ($\mathcal{L}$-erg) holds. In Section 6, we provide some examples of explicit invariant measures. In Section 6.4 we provide an encoding of any classical finite state Markov jump process into a stochastic master equation.

2. Construction of the model

2.1. Construction of quantum trajectories

In this section we fix the notations and introduce the probability space we use to study $(\hat{x}_{t})_{t}$. First, for an element $x\neq 0$ of $\mathbb{C}^{k}$, and for an operator $A$ with $Ax\neq 0$ we denote

[TABLE]

We consider the following distance on $\mathrm{P}\mathbb{C}^{k}$:

[TABLE]

for all $\hat{x},\hat{y}\in\mathrm{P}\mathbb{C}^{k}$, where $x$ and $y$ are norm-one representatives of $\hat{x}$ and $\hat{y}$ respectively. We equip $\mathrm{P}\mathbb{C}^{k}$ with the associated Borel $\sigma$-algebra denoted by $\mathcal{B}$.

Now we introduce a stochastic process with values in $\mathrm{M}_{k}(\mathbb{C})$. Let \big{(}\Omega,(\mathcal{F}_{t})_{t},\mathbb{P}\big{)} be a filtered probability space with standard brownian motions $W_{i}$ for $i\in I_{b}$, and standard Poisson processes $N_{j}$ for $j\in I_{p}$, such that the full family \big{(}W_{i},N_{j};\,i\in I_{b},j\in I_{p}\big{)} is independent. The filtration $(\mathcal{F}_{t})_{t}$ is assumed to satisfy the standard conditions, and we denote $\mathcal{F}_{\infty}$ by $\mathcal{F}$ and the processes \big{(}W_{i}(t)\big{)}_{t} and \big{(}N_{j}(t)-t\big{)}_{t} are $(\mathcal{F}_{t})_{t}$-martingales under $\mathbb{P}$. We denote by $\mathbb{E}$ the expectation with respect to $\mathbb{P}$.

On \big{(}\Omega,(\mathcal{F}_{t})_{t},\mathbb{P}\big{)}, for $s\in\mathbb{R}_{+}$, let $(S_{t}^{s})_{t\in[s,\infty)}$ be the solution to the following SDE:

[TABLE]

($\#I_{p}$** is the cardinal of $I_{p}$), where**

[TABLE]

Since standard Cauchy–Lipschitz conditions are fulfilled, the SDE defining $(S_{t}^{s})_{t}$ has indeed a unique (strong) solution. We denote $S_{t}:=S_{t}^{0}$. Note that for $s$ fixed the process $(S_{t}^{s})_{t}$ is independent of $\mathcal{F}_{s}$, and we have that for all $0\leq r\leq s\leq t$

[TABLE]

In addition, for any $\rho\in\mathcal{D}_{k}$, let $(Z_{t}^{\rho})_{t}$ be the positive real-valued process defined by

[TABLE]

and let $(\rho_{t})_{t}$ be the $\mathcal{D}_{k}$-valued process defined by

[TABLE]

if $Z_{t}^{\rho}\neq 0$, taking an arbitrarily fixed value whenever $Z_{t}^{\rho}=0$ (this value will always appear with probability zero in the sequel).

The following results on the properties of $(Z_{t}^{\rho})_{t}$ were proven in ****[6]****. We give short proofs adapted to our restricted setting where the Hilbert space is finite-dimensional, and $I=I_{b}\cup I_{p}$ is a finite set.

Lemma 2.1**.**

For any $\rho\in\mathcal{D}_{k}$, the stochastic process $(Z_{t}^{\rho})_{t}$ is the unique solution of the SDE

[TABLE]

Moreover, $(Z_{t}^{\rho})_{t}$ is a nonnegative martingale under $\mathbb{P}$.

Proof.

The fact that $(Z_{t}^{\rho})_{t}$ verifies the given SDE is a direct application of the Itô formula. Since $(\rho_{t})_{t}$ takes its values in the compact space $\mathcal{D}_{k}$, that SDE verifies standard Cauchy–Lipschitz conditions, ensuring the uniqueness of the solution. Since the processes \big{(}W_{i}(t)\big{)}_{t} and \big{(}N_{j}(t)-t\big{)}_{t} are $\mathbb{P}$-martingales, it follows that $(Z_{t}^{\rho})_{t}$ is a $\mathbb{P}$-local martingale. Since $\operatorname{tr}(C_{j}^{*}C_{j}\rho)\geq 0$ for any $j\in I_{p}$ and $\rho\in\mathcal{D}_{k}$, and $(\rho_{t})_{t}$ takes value in the compact space $\mathcal{D}_{k}$, it follows from [27, Theorem 12] that $(Z_{t}^{\rho})_{t\in[0,T]}$ is a $\mathbb{P}$-nonnegative martingale for all $T$. ∎

For any $\rho\in\mathcal{D}_{k}$, we define a probability $\mathbb{P}^{\rho}_{t}$ on $(\Omega,\mathcal{F}_{t})$:

[TABLE]

Since $(Z_{t}^{\rho})_{t}$ is a $\mathbb{P}$-martingale from Lemma 2.1, the family $(\mathbb{P}^{\rho}_{t})_{t}$ is consistent, that is $\mathbb{P}^{\rho}_{t}(E)=\mathbb{P}^{\rho}_{s}(E)$ for $t\geq s$ and $E\in\mathcal{F}_{s}$. Kolmogorov’s extension theorem defines a unique probability on $(\Omega,\mathcal{F}_{\infty})$, which we denote by $\mathbb{P}^{\rho}$. We will denote by $\mathbb{E}^{\rho}$ the expectation with respect to $\mathbb{P}^{\rho}$.

The following proposition makes explicit the relationship between $\mathbb{P}$ and $\mathbb{P}^{\rho}$. It follows from a direct application of Girsanov’s change of measure Theorem (see ****[26, Theorems III.3.24 and III.5.19]****). For all $i\in I_{b}$ and $t\in\mathbb{R}_{+}$, let

[TABLE]

Proposition 2.2**.**

Let $\rho\in\mathcal{D}_{k}$. Then, with respect to $\mathbb{P}^{\rho}$, the processes $\{B_{i}^{\rho}\}_{i\in I_{b}}$ are independent Wiener processes and the processes $\{N_{j}\}_{j\in I_{p}}$ are point processes of respective stochastic intensity $\{t\mapsto\operatorname{tr}(C_{j}^{*}C_{j}\rho_{t-})\}_{j\in I_{p}}$.

The process $(\rho_{t})_{t}$ considered under $\mathbb{P}^{\rho}$ models the evolution of a Markov open quantum system subject to indirect measurements. We refer the reader to ****[5, 14, 15]**** and references therein for a more detailed discussion of this interpretation.

From Itô calculus, $(\rho_{t})_{t}$ is solution of the SDE

[TABLE]

Proposition 2.2 then implies that with respect to $\mathbb{P}^{\rho}$, the process $(\rho_{t})_{t}$ is indeed the unique solution of (1.3) with $\rho_{0}=\rho$. Similarly, if $\rho_{0}=\pi_{\hat{x}}$ for some $\hat{x}\in\mathrm{P}\mathbb{C}^{k}$, then with respect to $\mathbb{P}^{\pi_{\hat{x}}}$, the process \big{(}\frac{S_{t}{x}}{\|S_{t}{x}\|}\big{)}_{t} is the solution of (1.4) with $x$ any norm one representative of $\hat{x}$.

Remark also that for any $\rho\in\mathcal{D}_{k}$, using (2.3), one has from Remark 1

[TABLE]

Our strategy of proof is based on the study of the joint distribution of $S_{t}$ and a random initial state $\hat{x}$. To this end, we consider the product space $\Omega\times\mathrm{P}\mathbb{C}^{k}$ equipped with the filtration $(\mathcal{F}_{t}\otimes\mathcal{B})_{t}$ and the full $\sigma$-algebra $\mathcal{F}\otimes\mathcal{B}$. For any probability measure $\mu$ on $\mathrm{P}\mathbb{C}^{k}$, and for all $E\in\mathcal{F}$ and $A\in\mathcal{B}$, let

[TABLE]

We will denote by $\mathbb{E}_{\mu}$ the expectation with respect to $\mathbb{Q}_{\mu}$. Note that $\mathrm{d}\mathbb{P}^{\pi_{\hat{x}}}_{t}=\|S_{t}x\|^{2}\,\mathrm{d}\mathbb{P}$ for any $\hat{x}\in\mathrm{P}\mathbb{C}^{k}$, so that $\mathbb{P}^{\pi_{\hat{x}}}(\{S_{t}x=0\})=0$ for all $x\in\hat{x}$. Therefore

[TABLE]

and there exists a process $(\hat{x}_{t})_{t}$ for which

[TABLE]

holds almost surely. It has the same distribution as the image by the map $x\mapsto\hat{x}$ of the solution $(x_{t})_{t}$ to (1.4) with $x_{0}\in\hat{x}$, $\|x_{0}\|=1$.

The following proposition shows that the laws of any $\mathcal{F}$-measurable random variables are given by a marginal of $\mathbb{Q}_{\mu}$. For a probability measure $\mu$ on $\mathrm{P}\mathbb{C}^{k}$, we define

[TABLE]

Proposition 2.3**.**

Let $\mu$ be a probability measure on $\mathrm{P}\mathbb{C}^{k}$, then $\rho_{\mu}\in\mathcal{D}_{k}$ and for any $E\in\mathcal{F}$,

[TABLE]

Proof.

The fact that $\rho_{\mu}\in\mathcal{D}_{k}$ follows from the positivity and linearity of the expectation. Concerning the second part, let $t\geq 0$ and $E\in\mathcal{F}_{t}$, then

[TABLE]

Fubini’s Theorem implies

[TABLE]

The uniqueness of the extended measure in Kolmogorov’s extension Theorem yields the proposition. ∎

Remark 2*.*

Any $\mathcal{F}$-measurable random variables $X$ can be extended canonically to a $\mathcal{F}\otimes\mathcal{B}$-measurable random variables setting $X(\omega,\hat{x})=X(\omega)$. Proposition 2.3 then implies that the distribution of a $\mathcal{F}$-measurable random variable under $\mathbb{Q}_{\mu}$ depends on $\mu$ only through $\rho_{\mu}$. The central idea of our proof is that assumption (Pur) will allow us to find a $\mathcal{F}$-measurable process approximating $(\hat{x}_{t})_{t}$. The $\mathcal{F}$-measurability of the process will then imply that it inherits some ergodicity properties from assumption ($\mathcal{L}$-erg)****.

Remark 3*.*

If $\mu_{\mathrm{inv}}$ is an invariant measure for the Markov chain $(\hat{x}_{t})_{t}$, then with the above notation, $\rho_{\mu_{\mathrm{inv}}}$ is an invariant state for $(\mathrm{e}^{t\mathcal{L}})_{t}$. In particular, if ($\mathcal{L}$-erg) holds then $\rho_{\mu_{\mathrm{inv}}}=\rho_{\mathrm{inv}}$. This follows from the identities

[TABLE]

where the second identity uses (2.5).

2.2. Key martingale

The following process is the key to construct a $\mathcal{F}$-measurable process approximating $(\hat{x}_{t})_{t}$. For any $t\geq 0$, let

[TABLE]

whenever $\operatorname{tr}(S_{t}^{*}S_{t})\neq 0$, and give $M_{t}$ a fixed arbitrary value whenever $\operatorname{tr}(S_{t}^{*}S_{t})=0$. Since, by definition, for any $\rho\in\mathcal{D}_{k}$, \mathbb{P}^{\rho}\big{(}\{\operatorname{tr}S_{t}^{*}S_{t}=0\}\big{)}=0, the arbitrary definition of $M_{t}$ on this set of vanishing probability is irrelevant. It turns out that with respect to $\mathbb{P}^{\mathrm{Id}/k}$, $(M_{t})_{t}$ is a martingale. For convenience we write $\mathbb{P}^{\mathrm{ch}}=\mathbb{P}^{\mathrm{Id}/k}$ and similarly for any other $\rho$-dependent object, whenever $\rho=\mathrm{Id}/k$.

Theorem 2.4**.**

With respect to $\mathbb{P}^{\mathrm{ch}}$, the stochastic process $(M_{t})_{t}$ is a bounded martingale. Therefore, it converges $\mathbb{P}^{\mathrm{ch}}$-almost surely and in $\mathrm{L}^{1}$ to a random variable $M_{\infty}$. Moreover, for any $\rho\in\mathcal{D}_{k}$,

[TABLE]

and $(M_{t})_{t}$ converges almost surely and in $\mathrm{L}^{1}$ to $M_{\infty}$ with respect to $\mathbb{P}^{\rho}$.

Proof.

Expressing $(S_{t})_{t}$ in terms of $B_{i}^{\mathrm{ch}}$ for $i\in I_{b}$, we have that

[TABLE]

Recall that the distributions of the $B_{i}^{\mathrm{ch}}$ and $N_{j}$ under $\mathbb{P}^{\mathrm{ch}}$ are given by Proposition 2.2.

Since $\operatorname{tr}(S_{t}^{*}S_{t})$ is $\mathbb{P}^{\mathrm{ch}}$-almost surely non zero, we can define $R_{t}$ by $R_{t}=S_{t}/\sqrt{\operatorname{tr}(S_{t}^{*}S_{t})}$ almost surely for $\mathbb{P}^{\mathrm{ch}}$, and therefore for $\mathbb{P}^{\rho}$ and $\mathbb{Q}_{\mu}$. The Itô formula implies

[TABLE]

Hence, with respect to $\mathbb{P}^{\mathrm{ch}}$, $(M_{t})_{t}$ is a local martingale. By definition, it is positive-semidefinite, and is also bounded since $\operatorname{tr}(M_{t})=1$ almost surely. Thus $(M_{t})_{t}$ is a martingale and standard theorems of convergence for martingales imply the convergence almost surely and in $\mathrm{L}^{1}$.

By direct computation, we get $\mathrm{d}\mathbb{P}^{\rho}|_{\mathcal{F}_{t}}=k\operatorname{tr}(\rho M_{t})\,\mathrm{d}\mathbb{P}^{\mathrm{ch}}|_{\mathcal{F}_{t}}$. The $\mathrm{L}^{1}$ convergence of $(M_{t})_{t}$ with respect to $\mathbb{P}^{\mathrm{ch}}$ then implies $\mathrm{d}\mathbb{P}^{\rho}=k\operatorname{tr}(\rho M_{\infty})\,\mathrm{d}\mathbb{P}^{\mathrm{ch}}$. Finally the inequality $\operatorname{tr}(AB)\leq\|A\|\operatorname{tr}(B)$ for any two positive semidefinite matrices implies $\mathbb{P}^{\rho}\leq k\,\mathbb{P}^{\mathrm{ch}}$, which yields the $\mathrm{L}^{1}$ and almost sure convergence with respect to $\mathbb{P}^{\rho}$. ∎

Now we are in the position to show that under the assumption (Pur) the limit $M_{\infty}$ is a rank-one projector. To this end let us introduce the polar decomposition of $(S_{t})_{t}$: there exists a process $(U_{t})_{t}$ with values in the set of $k\times k$ unitary matrices such that for all $t\geq 0$

[TABLE]

Proposition 2.5**.**

Assume that (Pur) holds. Then, for any $\rho\in\mathcal{D}_{k}$, $\mathbb{P}^{\rho}$-almost surely, the random variable $M_{\infty}$ is a rank-one orthogonal projector on $\mathbb{C}^{k}$.

Proof.

First, since $\mathbb{P}^{\rho}$ is absolutely continuous with respect to $\mathbb{P}^{\mathrm{ch}}$, proving the result with $\rho=\mathrm{Id}/k$ is sufficient. To achieve this, remark that the $\mathbb{P}^{\mathrm{ch}}$-almost sure convergence of $(M_{t})_{t}$ and the $\mathbb{P}^{\mathrm{ch}}$-almost sure bound $\sup_{t\geq 0}\|M_{t}\|\leq 1$ imply the convergence of $\mathbb{E}^{\mathrm{ch}}(M_{t}^{2})$. Now recall that $R_{t}=S_{t}/\sqrt{\operatorname{tr}(S_{t}^{*}S_{t})}$. The Itô isometry implies

[TABLE]

Therefore, the convergence of $\mathbb{E}^{\mathrm{ch}}(M_{t}^{2})$ to $\mathbb{E}^{\mathrm{ch}}(M_{\infty}^{2})$ implies that

[TABLE]

for all $i\in I_{b}$ and

[TABLE]

for all $j\in I_{p}$. Since the integrands are nonnegative, their inferior limits at infinity are [math]. Hence there exists an unbounded increasing sequence $(t_{n})_{n}$ such that for any $i\in I_{b}$,

[TABLE]

and for any $j\in I_{p}$,

[TABLE]

Since convergence in $\mathrm{L}^{1}$ implies the almost sure convergence of a subsequence, there exists an unbounded increasing sequence, which we denote also by $(t_{n})_{n}$, such that $\mathbb{P}^{\mathrm{ch}}$-almost surely,

[TABLE]

and

[TABLE]

for all $i\in I_{b}$ and $j\in I_{p}$.

Now and for the rest of this paragraph, fix a realization (i.e. an element of $\Omega$) such that $(M_{t_{n}})_{n}$ converges to $M_{\infty}$. The polar decomposition of $R_{t}$ is $R_{t}=U_{t}\sqrt{M_{t}}$. Since the set of $k\times k$ unitary matrices is compact, there exists a subsequence $(s_{n})_{n}$ of $(t_{n})_{n}$ such that $(U_{s_{n}})_{n}$ converges to $U_{\infty}$. We therefore have

[TABLE]

and

[TABLE]

for all $i\in I_{b}$ and $j\in I_{p}$. Denoting $P_{\infty}$ the orthogonal projector onto the range of $M_{\infty}$, it follows that there exist real numbers $(\alpha_{i})_{i\in I_{b}}$ and $(\beta_{j})_{j\in I_{p}}$ such that

[TABLE]

and

[TABLE]

Assumption (Pur) implies that the orthogonal projector $U_{\infty}P_{\infty}U_{\infty}^{*}$ has rank one, thus so does $P_{\infty}$. Since $\operatorname{tr}(M_{\infty})=1$, $M_{\infty}$ is a rank one orthogonal projector.

Since $(M_{t_{n}})_{n}$ converges $\mathbb{P}^{\mathrm{ch}}$-almost surely, the above paragraph and the absolute continuity of $\mathbb{P}^{\rho}$ with respect to $\mathbb{P}^{\mathrm{ch}}$ show that $M_{\infty}$ is $\mathbb{P}^{\rho}$-almost surely a rank one orthogonal projector. ∎

3. Invariant measure and exponential convergence in Wasserstein distance

This section is devoted to the main result of the paper, which concerns the exponential convergence to the invariant measure for the Markov process $(\hat{x}_{t})_{t}$. We first show a convergence result for ${\mathcal{F}}$-measurable random variables. The following theorem is a transcription of ****[38, Proposition 7.5]****.

Theorem 3.1**.**

Assume that ($\mathcal{L}$-erg) holds. Then there exist two constants $C>0$ and $\lambda>0$ such that for any $\rho\in\mathcal{D}_{k}$ and any $t\geq 0$,

[TABLE]

Our next proposition requires the introduction of a shift semigroup. From now on we assume that \big{(}\Omega,(\mathcal{F}_{t})_{t},\mathbb{P}\big{)} is a canonical realization of the processes $W_{i}$ and $N_{j}$, in particular $\Omega$ is a subset of $(\mathbb{R}^{I_{b}\cup I_{p}})^{\mathbb{R}_{+}}$. We can then define for every $t\geq 0$ the map $\theta^{t}$ on $\Omega$ by

[TABLE]

From the previous theorem we deduce the following proposition for $\mathcal{F}$-measurable random variables.

Proposition 3.2**.**

Assume ($\mathcal{L}$-erg) holds. Then there exist two constants $C>0$ and $\lambda>0$ such that for any $\mathcal{F}$-measurable, essentially bounded function $f:\Omega\mapsto\mathbb{C}$ with essential bound $\|f\|_{\infty}$, any $t\geq 0$ and any $\rho\in\mathcal{D}_{k}$,

[TABLE]

Proof.

Recall that by definition $\mathbb{P}$ is the law of processes with independent increments. It follows that if $g$ is $\mathcal{F}_{t}$-measurable and $h$ is $\mathcal{F}$-measurable, $\mathbb{E}(h\circ\theta^{t}g)=\mathbb{E}(h\circ\theta^{t})\mathbb{E}(g)$. Then, by definition of $\mathbb{P}^{\rho}$,

[TABLE]

Since $Z_{t+s}^{\rho}=\operatorname{tr}(S_{t+s}^{t\,*}S_{t+s}^{t}S_{t}\rho S_{t}^{*})$ where $S_{t}^{*}\rho S_{t}$ is $\mathcal{F}_{t}$-measurable and $S_{t+s}^{t\,*}S_{t+s}^{t}=S_{s}^{*}S_{s}\circ\theta^{t}$ by (2.2)

[TABLE]

Then relation (2.5), the $\theta$-invariance of $\mathbb{P}$, and the definition of the measures $\mathbb{P}^{\rho}$ yield

[TABLE]

with $\bar{\rho}_{t}=\mathrm{e}^{t\mathcal{L}}(\rho)$. It follows from Theorem 2.4 that

[TABLE]

For any matrix $A$, denoting $\|A\|_{1}$ its trace norm, \big{|}\operatorname{tr}(M_{\infty}A)\big{|}\leq\|A\|_{1}. Therefore,

[TABLE]

Theorem 3.1 then yields the proposition. ∎

The main strategy to show Theorem 1.1 is to construct a $\mathcal{F}$-measurable process $(\hat{y}_{t})_{t}$ approximating the process $(\hat{x}_{t})_{t}$. Let $(\hat{z}_{t})_{t}$ be the maximum likelihood process:

[TABLE]

where $x$ is a norm one representative of $\hat{x}$. If the largest eigenvalue of $S_{t}^{*}S_{t}$ is not simple, the choice of $\hat{z}_{t}$ may not be unique. However we can always choose an appropriate $\hat{z}_{t}$ in an $(\mathcal{F}_{t})_{t}$-adapted way. If (Pur) holds, Proposition 2.5 ensures that the definition of $\hat{z}_{t}$ is almost surely unambiguous for large enough $t$: it is the equivalence class of eigenvectors of $M_{t}$ corresponding to its largest eigenvalue.

Let now $(\hat{y}_{t})_{t}$ be the evolution of this maximum likelihood estimate:

[TABLE]

We shall also use the notation $\hat{z}_{t+s}^{s}:=\hat{z}_{t}\circ\theta^{s}$ and $\hat{y}_{t+s}^{s}:=\hat{y}_{t}\circ\theta^{s}$, that is, processes defined in the same fashion but substituting $S_{t+s}^{s}$ for $S_{t}$. It is worth noticing that these processes are all $\mathcal{F}$-measurable.

Our proof that $(\hat{y}_{t})_{t}$ is an exponentially good approximation of $(\hat{x}_{t})_{t}$ relies in part on the use of the exterior product of $\mathbb{C}^{k}$. We recall briefly the relevant definitions: for $x_{1},x_{2}\in\mathbb{C}^{k}$ we denote by $x_{1}\wedge x_{2}$ the alternating bilinear form

[TABLE]

Then, the set of all $x_{1}\wedge x_{2}$ is a generating family for the set $\wedge^{2}\mathbb{C}^{k}$ of alternating bilinear forms on $\mathbb{C}^{k}$. We equip it with a complex inner product by

[TABLE]

and denote by $\|x_{1}\wedge x_{2}\|$ the associated norm (there should be no confusion with the norm on vectors). It is immediate to verify that our metric $d(\cdot,\cdot)$ on $\mathrm{P}\mathbb{C}^{k}$ satisfies

[TABLE]

For $A\in\mathrm{M}_{k}(\mathbb{C})$, we write $\wedge^{2}A$ for the operator on $\wedge^{2}\mathbb{C}^{k}$ defined by

[TABLE]

It follows that $\wedge^{2}(AB)=\wedge^{2}A\,\wedge^{2}B$, so that $\|\wedge^{2}\!(AB)\|\leq\|\wedge^{2}A\|\|\wedge^{2}B\|$. There exists a useful relationship between the operator norm on $\wedge^{2}\mathrm{M}_{k}(\mathbb{C})$ and singular values of matrices. From e.g.  Chapter XVI of ****[32]****,

[TABLE]

where $a_{1}(A)\geq a_{2}(A)$ are the two first singular values of $A$, i.e.  the square roots of eigenvalues of $A^{*}A$. We recall that the operator norm is defined such that $\|A\|:=a_{1}(A)$.

The exponential decrease of $d(\hat{x}_{t},\hat{y}_{t})$ is derived from the exponential decay of the following function:

[TABLE]

Lemma 3.3**.**

Assume that (Pur) holds. Then there exist two constants $C>0$ and $\lambda>0$ such that for all $t\geq 0$

[TABLE]

Proof.

First, we show that $f$ converges to zero as $t$ grows to $\infty$. To this end recall that $R_{t}=S_{t}/\sqrt{kZ_{t}^{\mathrm{ch}}}$, so that

[TABLE]

Furthermore, since $R_{t}^{*}R_{t}=M_{t}$, we have from Theorem 2.4 and Proposition 2.5 that

[TABLE]

Indeed, since $a_{1}(R_{t})$ and $a_{2}(R_{t})$ are the largest two eigenvalues of $\sqrt{M_{t}}$, the fact that it converges to a rank-one projector implies that $a_{1}(R_{t})$ converges to $1$ and $a_{2}(R_{t})$ to zero. The inequality $\|S_{t}\|^{2}\leq\operatorname{tr}(S_{t}^{*}S_{t})$ implies $\|\wedge^{2}R_{t}\|\leq 1$ almost surely. Then Lebesgue’s dominated convergence theorem yields $\lim_{t\to\infty}f(t)=0$.

Second, we show $f$ is submultiplicative. By the semi-group property, $S_{t+s}=S_{t+s}^{s}S_{s}$ for all $t,s\geq 0$. Using that the norm is submultiplicative, for any $t,s\geq 0$,

[TABLE]

Since $\mathbb{P}$ has independent increments, $\|\wedge^{2}S_{t+s}^{s}\|=\|\wedge^{2}S_{t}\|\circ\theta^{s}$ is $\mathcal{F}_{s}$-independent and $\|\wedge^{2}S_{s}\|$ is $\mathcal{F}_{s}$-measurable,

[TABLE]

The measure $\mathbb{P}$ being shift-invariant,

[TABLE]

which yields that $f$ is submultiplicative.

Since $f$ is measurable, submultiplicative and $0\leq f(t)\leq k$ for all $t$, Fekete’s subadditive lemma ensures that there exists $\lambda\in(-\infty,\infty]$ such that

[TABLE]

Since $f$ converges towards [math], this $\lambda$ belongs to $(0,\infty]$. This yields the lemma. ∎

Proposition 3.4**.**

Assume that (Pur) holds. Then there exist two constants $C>0$ and $\lambda>0$ such that for any $s,t\in\mathbb{R}_{+}$ and for any probability measure $\mu$ on $(\mathrm{P}\mathbb{C}^{k},\mathcal{B})$,

[TABLE]

Proof.

Recall that $\mathbb{E}_{\mu}$ is the expectation with respect to $\mathbb{Q}_{\mu}$. Using the Markov property, we have

[TABLE]

with $\mu_{s}$ the distribution of $\hat{x}_{s}$ conditioned on $\hat{x}_{0}\sim\mu$. Then it is sufficient to prove the proposition for $s=0$. For any $t\geq 0$, using the fact that $\|S_{t}z_{t}\|=\|S_{t}\|$ for $z_{t}$ a norm one representative of $\hat{z}_{t}$,

[TABLE]

Using this inequality and the fact that $\mathrm{d}\mathbb{Q}_{\mu}{}|_{\mathcal{F}_{t}\otimes\mathcal{B}}=\|S_{t}x_{0}\|^{2}\,\mathrm{d}\mathbb{P}\,\mathrm{d}\mu(\hat{x}_{0})$,

[TABLE]

Finally Lemma 3.3 yields the proposition. ∎

We turn to the proof of our main theorem, Theorem 1.1. The speed of convergence is expressed in terms of the Wasserstein distance $W_{1}$. Let us recall the definition of this distance for compact metric spaces: for $X$ a compact metric space equipped with its Borel $\sigma$-algebra, the Wasserstein distance of order $1$ between two probability measures $\sigma$ and $\tau$ on $X$ can be defined using the Kantorovich–Rubinstein duality Theorem as

[TABLE]

where $\mathrm{Lip}_{1}(X)=\{f:X\rightarrow\mathbb{R}\ \mathrm{s.t.}\ |f(x)-f(y)|\leq d(x,y)\}$ is the set of Lipschitz continuous functions with constant one, and $d$ is the metric on $X$. Here we use this for $X=\mathrm{P}\mathbb{C}^{k}$ and $d$ defined in (2.1) (see also (3.4)).

We recall our main theorem before proving it.

Theorem 3.5**.**

Assume that (Pur) and ($\mathcal{L}$-erg) hold. Then the Markov process $(\hat{x}_{t})_{t}$ has a unique invariant probability measure $\mu_{\mathrm{inv}}$, and there exist $C>0$ and $\lambda>0$ such that for any initial distribution $\mu$ of $\hat{x}_{0}$ over $\mathrm{P}\mathbb{C}^{k}$, for all $t\geq 0$, the distribution $\mu_{t}$ of $\hat{x}_{t}$ satisfies

[TABLE]

where $W_{1}$ is the Wasserstein distance of order $1$.

Proof.

Let $f\in\mathrm{Lip}_{1}(\mathrm{P}\mathbb{C}^{k})$. From the definition of Wasserstein distance, we can restrict ourselves to functions $f$ that vanish at some point. Remark that since $\sup_{\hat{x},\hat{y}\in\mathrm{P}\mathbb{C}^{k}}d(\hat{x},\hat{y})=1$, restricting to this set of functions implies $\|f\|_{\infty}\leq 1$. Let $\mu_{\mathrm{inv}}$ be an invariant probability measure for $(\hat{x}_{t})_{t}$. We will prove the exponential convergence of $(\mu_{t})_{t}$ towards $\mu_{\mathrm{inv}}$ for any initial $\mu_{0}$, and that will imply that $(\hat{x}_{t})_{t}$ accepts a unique invariant probability measure. Let $t\geq 0$, and recall that $\hat{y}_{t}^{t/2}=\hat{y}_{t/2}\circ\theta^{t/2}$. We have

[TABLE]

The two terms on the right hand side of line (3.9) are bounded using Proposition 3.4. Using Proposition 2.3, the difference on line (3.10) satisfies

[TABLE]

Then bounding the right hand side using Proposition 3.2, it follows there exist $C>0$ and $\lambda>0$ such that

[TABLE]

Adapting the two constants yields the theorem. ∎

4. Set of invariant measures under (Pur)

The results and proofs of this section are a direct translation of ****[11, Appendix B]****. We reproduce the proofs for the reader’s convenience.

Whenever ($\mathcal{L}$-erg) does not hold, $\dim\ker\mathcal{L}>1$ and the semigroup $(\mathrm{e}^{t\mathcal{L}})_{t}$ accepts more than one fixed point in $\mathcal{D}_{k}$. The convex set of invariant states can be explicitly classified given the matrices $(L_{i})_{i\in I_{b}}$ and $(C_{j})_{j\in I_{b}}$. Following ****[9, Theorem 7]**** (alternatively see Theorem 7.2 and Proposition 7.6 in ****[38]****, and ****[36]****), there exists a decomposition

[TABLE]

with the following properties:

1. (1)

The range of any invariant states is a subspace of $V=\mathbb{C}^{n_{1}}\oplus\cdots\oplus\mathbb{C}^{n_{d}}\oplus\{0\}$; 2. (2)

The restriction of the operators $L_{i}$ and $C_{j}$ to $\mathbb{C}^{n_{1}}\oplus\cdots\oplus\mathbb{C}^{n_{d}}$ are block-diagonal, with

[TABLE] 3. (3)

For each $\ell=1,\ldots,d$ there is a decomposition $\mathbb{C}^{n_{\ell}}=\mathbb{C}^{k_{\ell}}\otimes\mathbb{C}^{m_{\ell}},\,n_{\ell}=k_{\ell}\times m_{\ell}$, a unitary matrix $U_{\ell}$ on $\mathbb{C}^{n_{\ell}}$ and matrices $\{\hat{L}_{\ell,i}\}_{i\in I_{b}}$ and $\{\hat{C}_{\ell,j}\}_{j\in I_{p}}$ on $\mathbb{C}^{k_{\ell}}$ such that

[TABLE] 4. (4)

There exists a positive definite matrix $\rho_{\ell}$ on $\mathbb{C}^{k_{\ell}}$ such that

[TABLE]

is a fixed point of $(\mathrm{e}^{t\mathcal{L}})_{t}$.

Then, the set of fixed points for $(\mathrm{e}^{t\mathcal{L}})$ is

[TABLE]

The decomposition simplifies under the purification assumption.

Proposition 4.1**.**

Assume that (Pur) holds. Then there exists a set $\{\rho_{\ell}\}_{\ell=1}^{d}$ of positive definite matrices and an integer $D$ such that the set of fixed points of $(\mathrm{e}^{t\mathcal{L}})_{t}$ is

[TABLE]

Proof.

The statement follows from the discussion preceding the proposition if we show that (Pur) implies $m_{1}=\ldots=m_{d}=1$. Assume that one of the $m_{\ell}$, e.g. $m_{1}$, is greater than $1$. Let $x$ be a norm one vector in $\mathbb{C}^{k_{1}}$. Then $\pi=U_{1}(\pi_{\hat{x}}\otimes I_{\mathbb{C}^{m_{1}}})U_{1}^{*}\oplus 0\oplus\cdots\oplus 0$ is an orthogonal projection of rank $m_{1}>1$, and

[TABLE]

and this contradicts (Pur)****. ∎

It is clear from Proposition 4.1 that to each extremal fixed point $0\oplus\cdots\oplus\rho_{\ell}\oplus\cdots\oplus 0$ corresponds a unique invariant measure $\mu_{\ell}$ supported on its range $\mathop{\rm ran}\nolimits\rho_{\ell}$. The converse is the subject of the next proposition.

Proposition 4.2**.**

Assume (Pur) holds. Then any invariant probability measure of $(\hat{x}_{t})_{t}$ is a convex combination of the measures $\mu_{\ell}$, $\ell=1,\ldots,d$.

Proof.

Let $\mu$ be an invariant probability measure for $(\hat{x}_{t})_{t}$ and $f$ be a continuous function on $\mathrm{P}\mathbb{C}^{k}$. Proposition 3.4 implies that

[TABLE]

Since $(\hat{y}_{t})_{t}$ is $\mathcal{F}$-measurable, Proposition 2.3 implies \int f\,\mathrm{d}\mu=\lim_{t\to\infty}\mathbb{E}^{\rho_{\mu}}\big{(}f(\hat{y}_{t})\big{)}, and by Remark 3, $\rho_{\mu}\in\mathcal{D}_{k}$ is a fixed point of $(\mathrm{e}^{t\mathcal{L}})_{t}$. Proposition 4.1 ensures that there exist nonnegative numbers $p_{1},\ldots,p_{d}$ summing up to one such that $\rho_{\mu}=p_{1}\,\rho_{1}\oplus\cdots\oplus p_{d}\,\rho_{d}\oplus 0_{M_{D}(\mathbb{C})}$. From the definition of $\mathbb{P}^{\rho_{\mu}}$,

[TABLE]

with the abuse of notation $\rho_{\ell}:=0\oplus\cdots\oplus\rho_{\ell}\oplus\cdots\oplus 0$, so that

[TABLE]

The same argument gives \int f\,\mathrm{d}\mu_{\ell}=\lim_{t\to\infty}\mathbb{E}^{\rho_{\ell}}\big{(}f(\hat{y}_{t})\big{)}, and we have $\mu=p_{1}\,\mu_{1}+\ldots+p_{d}\,\mu_{d}$. ∎

5. (Pur) is not necessary for purification

As shown by the following example, the condition (Pur) is sufficient but not necessary for (1.5) to hold.

Let $k=3$ and fix an orthonormal basis $\{e_{1},e_{2},e_{3}\}$ of $\mathbb{C}^{3}$. Let $I_{b}=\{0,1\}$, $I_{p}=\{2\}$, $u=(e_{1}+e_{2}+e_{3})/\sqrt{3}$ and $v=(e_{1}+e_{3})/\sqrt{2}$. Let

[TABLE]

Proposition 5.1**.**

Let $\mathcal{L}$ be the Lindblad operator given by (1.2) with $H$, $V_{1}$, $V_{2}$, $V_{3}$ defined in (5.1). Then ($\mathcal{L}$-erg) holds and the unique invariant state $\rho_{\mathrm{inv}}$ is positive definite.

Proof.

Using [38, Proposition 7.6], it is sufficient to prove that if $\pi$ is a non null orthogonal projector such that $(\mathrm{Id}-\pi)L_{0}\pi=(\mathrm{Id}-\pi)L_{1}\pi=(\mathrm{Id}-\pi)C_{2}\pi=0$, then $\pi=\mathrm{Id}$. Assume $\operatorname{rank}\pi<3$. Since $\pi\in M_{3}(\mathbb{C})$, there exist $\hat{x}\in\mathrm{P}\mathbb{C}^{3}$ such that either $\pi=\pi_{\hat{x}}$ or $\pi=\mathrm{Id}-\pi_{\hat{x}}$. If the first alternative holds, $(\mathrm{Id}-\pi)L_{0}\pi=(\mathrm{Id}-\pi)L_{1}\pi=(\mathrm{Id}-\pi)C_{2}\pi=0$ implies $\hat{x}$ is the equivalence class of a common eigenvector of $L_{0}$, $L_{1}$ and $C_{2}$. If the second alternative holds, $\hat{x}$ is the equivalence class of a common eigenvector of $L_{0}^{*}$, $L_{1}^{*}$ and $C_{2}^{*}$. The only common eigenvectors of $L_{0}$ and $C_{2}$ or $L_{0}^{*}$ and $C_{2}^{*}$ are elements of $\mathbb{C}(e_{2}-e_{3})$. Since $L_{1}$ is self adjoint, and this eigenspace is not an eigenspace of $L_{1}$, the proposition holds. ∎

In the orthonormal basis $\{e_{1},e_{2},e_{3}\}$,

[TABLE]

Taking $\pi$ the orthogonal projector onto the subspace spanned by $\{e_{2},e_{3}\}$ it follows that (Pur) does not hold. Yet we have the following proposition.

Proposition 5.2**.**

Consider the family of processes $(\rho_{t})_{t}$ defined by (1.3) with $H$, $L_{0}$, $L_{1}$, $C_{2}$ defined in (5.1). Then for any $\rho\in\mathcal{D}_{k}$,

[TABLE]

Proof.

Proposition 5.1 implies that $\rho_{\mathrm{inv}}$, the unique element of $\mathcal{D}_{k}$ invariant by $(\mathrm{e}^{t\mathcal{L}})_{t}$ is positive definite. Then, $\operatorname{tr}(C_{2}^{*}C_{2}\rho_{\mathrm{inv}})>0$. The results of [29] thus ensure that for any $\rho\in\mathcal{D}_{k}$,

[TABLE]

Let $T=\inf\{t\geq 0:N_{2}(t)\geq 1\}$. Then $\mathbb{P}^{\rho}(T<\infty)=1$ and from the definition of $C_{2}$,

[TABLE]

Hence $\inf_{\hat{y}\in\mathrm{P}\mathbb{C}^{k}}\|\rho_{t}-\pi_{\hat{y}}\|=0$ for any $t\geq T$ and $\mathbb{P}^{\rho}(T<\infty)=1$ yield the proposition. ∎

Corollary 5.3**.**

Consider the process $(\hat{x}_{t})_{t}$ defined by (1.4) with $H$, $L_{0}$, $L_{1}$, $C_{2}$ defined in (5.1). Then $(\hat{x}_{t})_{t}$ accepts a unique invariant probability measure $\mu_{\mathrm{inv}}$ and there exist $C>0$ and $\lambda>0$ such that for any initial distribution $\mu$ of $\hat{x}_{0}$ over $\mathrm{P}\mathbb{C}^{3}$, for all $t\geq 0$, the distribution $\mu_{t}$ of $\hat{x}_{t}$ satisfies

[TABLE]

Proof.

It is a direct adaptation of our proof of Theorem 1.1. Indeed Theorem 1.1 holds if one substitutes the conclusion of Proposition 2.5 for (Pur)****. Taking $\rho_{0}=\mathrm{Id}/3$ in the latter proposition yields $\rho_{t}=\frac{S_{t}S_{t}^{*}}{\operatorname{tr}(S_{t}S_{t}^{*})}$ and $M_{t}=\frac{S_{t}^{*}S_{t}}{\operatorname{tr}(S_{t}S_{t}^{*})}$. Therefore, $\rho_{t}$ and $M_{t}$ are unitarily equivalent. Following the arguments and notation of the proof of Proposition 5.2 we see that $\mathbb{P}^{\mathrm{ch}}$-almost surely, $M_{T}$ has rank one, and so does any $M_{t}$ for $t\geq T$. Hence the conclusion of Proposition 2.5 holds and the corollary is proven. ∎

Following the proofs of the discrete-time results of ****[11**]****, we can prove that the implication in Proposition 2.5 is an equivalence if (Pur) is replaced by

---

**(NSC-Pur)****: **

Any non zero orthogonal projector $\pi$ that satisfies $\pi S_{t}^{*}S_{t}\pi\propto\pi$ $\mathbb{P}$-almost surely for any $t\geq 0$ has rank one.

Alas, in practice, such a condition is hard to check.

6. Examples

In the following examples $k=2$. We recall the definition of the Pauli matrices:

[TABLE]

A standard orthonormal basis of $M_{2}(\mathbb{C})$ equipped with the Hilbert–Schmidt inner product is

[TABLE]

In the basis of Pauli matrices one can write in a unique way any projection $\pi_{\hat{x}}$ as

[TABLE]

where

[TABLE]

We denote in particular by $\mathcal{X}_{t},\mathcal{Y}_{t},\mathcal{Z}_{t}$ respectively the coordinates associated with $\pi_{\hat{x}_{t}}$.

6.1. Unitarily perturbed non demolition diffusive measurement

Our first example consists of a $\genfrac{}{}{}{1}{1}{2}$-spin (or qbit) in a magnetic field oriented along the $y$-axis and subject to indirect non demolition measurement along the $z$-axis. It is a typical quantum optics experimental situation (see for example ****[20]****). In terms of the parameters defining the related quantum trajectories, we get $H=\sigma_{y}$, $I_{b}=\{0\}$, $I_{p}=\emptyset$ and $L_{0}=\sqrt{\gamma}\,\sigma_{z}$ with $\gamma>0$. Then $(\pi_{\hat{x}_{t}})_{t}$ conditioned on $\hat{x}_{0}$ is the solution of

[TABLE]

For this quantum trajectory it is immediate to verify (Pur), and solving $\mathcal{L}(\rho)=0$ shows that $\rho_{\mathrm{inv}}=\genfrac{}{}{}{1}{1}{2}\mathrm{Id}$ is the unique invariant state, so that ($\mathcal{L}$-erg) holds. Hence, by Theorem 1.1 $(\hat{x}_{t})_{t}$ has a unique invariant measure. In the following we derive an explicit expression for this invariant measure.

The next Lemma allows us to restrict the state space.

Lemma 6.1**.**

If $\mu(\mathcal{Y}_{0}=0)=1$ then $\mathbb{Q}_{\mu}(\mathcal{Y}_{t}=0)=1$ for all $t$ in $\mathbb{R}$.

Proof.

From equation (6.1), $(\mathcal{Y}_{t})_{t}$ is the solution of $\mathrm{d}\mathcal{Y}_{t}=2\mathcal{Y}_{t}(\gamma\,\mathrm{d}t-\sqrt{\gamma}\,\mathcal{Z}_{t}\,\mathrm{d}B_{t})$. It is therefore a Doléans-Dade exponential:

[TABLE]

and the conclusion follows. ∎

Now we prove that the invariant measure admits a rotational symmetry.

Lemma 6.2**.**

Assume that the distribution $\mu$ of $\hat{x}_{0}$ is invariant with respect to the mapping $\hat{x}\mapsto\sigma_{y}\cdot\hat{x}$. Then $\mu_{t}$ is invariant with respect to the same mapping.

Proof.

Since $\sigma_{y}$ is unitary and self-adjoint, we have $\pi_{\sigma_{y}\cdot\hat{x}}=\sigma_{y}\pi_{\hat{x}}\sigma_{y}$. Since $\sigma_{y}\sigma_{z}=-\sigma_{z}\sigma_{y}$, it follows from (6.1) that

[TABLE]

Then it follows from $\sigma_{y}\cdot\hat{x}_{0}\sim\hat{x}_{0}$ and $(B_{t})_{t}\sim(-B_{t})_{t}$ that $(\hat{x}_{t})_{t}$ and $(\sigma_{y}\cdot\hat{x}_{t})_{t}$ are both weak solutions to the same SDE with the same initial condition. Since this SDE has a unique solution, they have the same distributions. ∎

Proposition 6.3**.**

Let $(\hat{x}_{t})_{t}$ be the process defined by (6.1). Then its unique invariant measure is the normalized image measure by

[TABLE]

of the measure $\tau(\theta)\,\mathrm{d}\theta$ on $(-\pi,\pi]$ with

[TABLE]

for $\theta\in[0,\pi]$ and $\tau(\theta)=\tau(\theta+\pi)$ for $\theta\in(-\pi,0]$.

Proof.

The convergence results in Theorem 1.1 and Lemma 6.1 imply that the invariant measure $\mu_{\mathrm{inv}}$ is the image by $\iota$ of a probability measure $\uptau$ on $(-\pi,\pi]$. Let $(\theta_{t})_{t}$ be the solution of

[TABLE]

with initial condition $\theta_{0}$. Remark that $(\theta_{t})_{t}$ is $2\pi$-periodic with respect to its initial condition, namely, $(\theta_{t}+2\pi)_{t}$ is solution of (6.3) with initial condition $\theta_{0}+2\pi$. Now, using the Itô formula,

[TABLE]

for $(\pi_{\hat{x}_{t}})_{t}$ solution of (6.1) with initial condition \hat{x}_{0}=\genfrac{}{}{}{1}{1}{2}\big{(}\mathrm{Id}+\sin\theta_{0}\,\sigma_{x}+\cos\theta_{0}\,\sigma_{z}\big{)}. Hence \big{(}\iota(\theta_{t})\big{)}_{t} has the same distribution as $(\hat{x}_{t})_{t}$. Therefore $\uptau$ is an invariant measure for the diffusion defined by (6.3); in addition, Theorem 1.1 shows that this invariant measure is unique, and Lemma 6.2 shows that it is $\pi$-periodic. Following standard methods (see [28]), one shows that the restriction of $\uptau$ to $[0,\pi)$ has a density of the form $\tau(\theta)=C_{1}\tau_{1}(\theta)+C_{2}\tau_{2}(\theta)$ with $C_{1},C_{2}\in\mathbb{R}$ and

[TABLE]

Now, straightforward analysis shows that $\int_{0}^{\pi}\tau_{1}(\theta)\,\mathrm{d}\theta<\infty$ while $\int_{0}^{\pi}\tau_{2}(\theta)\,\mathrm{d}\theta=\infty$. Therefore, $\tau$ is proportional to $\tau_{1}$ and the result follows.

∎

Remark 4*.*

For $\gamma\to\infty$, the invariant measure $\mu_{\mathrm{inv}}$ in Proposition 6.3 is a Dirac measure at [math] and $\pi$. To describe the scaling for $\gamma$ large we embed $\tau(\theta)$ into $L^{1}(\mathbb{R})$ by defining it to be zero outside the region $(-\pi,\pi]$. Then on the positive half line, in the $L^{1}$ norm,

[TABLE]

Hence, for large $\gamma$, the stationary probability distribution has two peaks of width (of order) $1/\gamma$ located $1/\gamma$ radians clockwise from the limit points [math] and $\pi$. Furthermore the probability to find the particle around the limit points is exponentially suppressed.

The strong noise limit, $\gamma\to\infty$, was recently studied in various models [2, 8, 12]. This is the first model that allows for an explicit calculation of the shape of the stationary probability measure. The density of the invariant probability distribution is plotted in Figure 1 for three values of $\gamma$, and for $\theta\in[0,\pi]$.

6.2. Thermal qubit, diffusive case

The following second example corresponds to the evolution of a qubit interacting weakly with the electromagnetic field at a fixed temperature. The emission and absorption of photons by the qubit are stimulated by a resonant coherent field (laser). In the limit of a strong stimulating laser, the measurement of emitted photons results in a diffusive signal whose drift depends on the instantaneous average value of the raising and lowering operators of the qubit (see ****[37, §4.4]**** for a more detailed physical derivation). We obtain an analytically solvable model if we assume that the unitary rotation of the qubit is compensated for and thus frozen. In terms of the parameters defining the related quantum trajectories, we get $H=0$, $I=I_{b}=\{0,1\}$, $L_{0}=\sqrt{a}\,\sigma_{+}$ and $L_{1}=\sqrt{b}\,\sigma_{-}$ with $a,b\in\mathbb{R}_{+}\setminus\{0\}$ and $\sigma_{\pm}=\genfrac{}{}{}{1}{1}{2}(\sigma_{x}\pm\mathrm{i}\sigma_{y})$, so that $\sigma_{+}=\begin{pmatrix}0&1\\ 0&0\end{pmatrix}$ and $\sigma_{-}=\begin{pmatrix}0&0\\ 1&0\end{pmatrix}$.

The stochastic master equation satisfied by $\pi_{\hat{x}_{t}}$ is

[TABLE]

Again it is immediate to verify (Pur), and solving for $\mathcal{L}(\rho)=0$ shows that ($\mathcal{L}$-erg) holds.

Lemma 6.4**.**

If $\mu(\mathcal{Y}_{0}=0)=1$ then $\mathbb{Q}_{\mu}(\mathcal{Y}_{t}=0)=1$ for all $t$ in $\mathbb{R}$.

Proof.

From (6.4), $(\mathcal{X}_{t})_{t}$ and $(\mathcal{Y}_{t})_{t}$ satisfy

[TABLE]

Therefore, if one defines

[TABLE]

then one has $\mathcal{Y}_{t}=\mathcal{Y}_{0}\,\mathrm{e}^{-\genfrac{}{}{}{1}{1}{2}(a+b)t}M_{t}$, and this proves Lemma 6.4. ∎

Proposition 6.5**.**

Let $(\hat{x}_{t})_{t}$ be the process defined by (6.4). Then its unique invariant measure is the normalized image measure by $\iota$ (defined by (6.2)) of the measure $\tau(\theta)\,\mathrm{d}\theta$ on $(-\pi,\pi]$ with

[TABLE]

with $z=\frac{a-b}{a+b}$ and $\varsigma=\frac{a+b}{2\sqrt{ab}}.$

Proof.

As in the proof of Proposition 6.3, Theorem 1.1 and Lemma 6.4 imply that the invariant measure $\mu_{\mathrm{inv}}$ is the image by $\iota$ of a probability measure $\uptau$ on $(-\pi,\pi]$. Let $(\theta_{t})_{t}$ be the solution of

[TABLE]

The Itô formula implies once again

[TABLE]

for $(\pi_{\hat{x}_{t}})_{t}$ solution of (6.4) with initial condition \hat{x}_{0}=\genfrac{}{}{}{1}{1}{2}\big{(}\mathrm{Id}+\sin\theta_{0}\,\sigma_{x}+\cos\theta_{0}\,\sigma_{z}\big{)}. Hence \big{(}\iota(\theta_{t})\big{)} has the same distribution as $(\hat{x}_{t})$. As in the proof of Proposition 6.3, standard techniques show that the unique invariant distribution for (6.5) has density proportional to the function $\tau$ above. ∎

The density of the invariant probability distribution for three values of the pair $(a,b)$ is plotted in Figure 2.

6.3. Thermal qubit, jump case

Our third example is the second one where the stimulating coherent field has relatively small amplitude and is filtered out. Then, the signal is composed only of the photons absorbed or emitted by the qubit. The resulting trajectory involves only jumps related to these events. The parameters defining the model are then, $H=0$, $I_{b}=\emptyset$ and $I_{p}=\{0,1\}$, $C_{0}=\sqrt{a}\,\sigma_{+}$ and $C_{1}=\sqrt{b}\,\sigma_{-}$ with $a,b\in\mathbb{R}_{+}\setminus\{0\}$.

The process $(\pi_{\hat{x}_{t}})_{t}$ is solution of

[TABLE]

where $N_{0}$ and $N_{1}$ are Poisson processes of stochastic intensities

[TABLE]

Assumptions (Pur) and ($\mathcal{L}$-erg) hold as in Example 6.2.

Proposition 6.6**.**

Let $\{e_{1},e_{2}\}$ denote the canonical basis of $\mathbb{C}^{2}$. The invariant measure for Equation (6.6) is

[TABLE]

Proof.

It is enough to check from (6.6) that, if ${\hat{x}_{0}}$ is either $\hat{e}_{1}$ or $\hat{e}_{2}$ then $(\pi_{\hat{x}_{t}})_{t}$ is a jump process on $(\pi_{\hat{e}_{1}},\pi_{\hat{e}_{2}})$ with intensity $b$ for the jumps from $\pi_{\hat{e}_{1}}$ to $\pi_{\hat{e}_{2}}$ and intensity $a$ for the reverse jumps. ∎

6.4. Finite state space Markov process embedding

In this last example we show how we can recover all the usual continuous-time Markov chains using special quantum trajectories.

Let $\{e_{\ell}\}_{\ell=1}^{k}$ be an orthonormal basis of $\mathbb{C}^{k}$, and $(X_{t})_{t}$ a $\{e_{1},\ldots,e_{k}\}$-valued Markov process with generator $Q$ (we recall that $Q$ is a $k\times k$ real matrix such that $\mathbb{E}(\langle v,X_{t}\rangle|X_{0})=\langle v,e^{tQ}X_{0}\rangle$ for any vector $v\in\mathbb{C}^{k}$). Let $H$ be diagonal in the basis $\{e_{\ell}\}_{\ell=1}^{k}$, let $I_{b}=\emptyset$ and $I_{p}=\{(i,j);i\neq j\mbox{ in }1,\ldots,k\}$ and for any $(i,j)\in I_{p}$ let $C_{i,j}=\sqrt{Q_{i,j}}\,e_{j}e_{i}^{*}$.

Proposition 6.7**.**

Let $(\hat{x}_{t})_{t}$ be the quantum trajectory defined by Equation (1.4) and the above parameters. Then assumption (Pur)** holds. In addition,**

1. (i)

Let T=\inf\big{\{}t\geq 0:\hat{x}_{t}\in\{\hat{e}_{1},\ldots,\hat{e}_{k}\}\big{\}}. If for all $i$ there exists $j$ with $Q_{i,j}>0$ then for any probability measure $\mu$ over $\mathrm{P}\mathbb{C}^{k}$, $\mathbb{P}_{\mu}(T<\infty)=1$. 2. (ii)

Conditionally on $\hat{x}_{0}\in\{\hat{e}_{\ell}\}_{\ell=1}^{k}$, the process $(\hat{x}_{t})_{t}$ has the same distribution as the image by $x\mapsto\hat{x}$ of $(X_{t})_{t}$. 3. (iii)

The assumption ($\mathcal{L}$-erg) holds if and only if $(X_{t})_{t}$ accepts a unique invariant measure. In that case, the unique invariant measure $\nu_{\mathrm{inv}}$ for $(\hat{x}_{t})_{t}$ is the image by $x\mapsto\hat{x}$ of the unique invariant measure for $(X_{t})_{t}$.

Proof.

Note first that any $C_{i,j}^{*}C_{i,j}=Q_{i,j}\,e_{i}e_{i}^{*}$ , so that (Pur)** holds trivially.**

To prove (i), let $T_{1}=\inf\{t>0;\ \exists(i,j)\in I_{p}\mbox{ such that }N_{i,j}(t)>0\}$. Remark that because $\operatorname{tr}(C_{i,j}\pi_{\hat{x}_{s-}}C_{i,j}^{*})=Q_{i,j}|\langle e_{i},x_{s-}\rangle|^{2}$, the sum $\sum_{i,j}N_{i,j}$ of independent Poisson processes has intensity

[TABLE]

where $Q_{i}=\sum_{j}Q_{i,j}$ is positive by assumption, so that $T_{1}$ is almost surely finite. Now consider the almost surely unique $(i,j)$ in $I_{p}$ such that $N_{i,j}(T_{1})>0$; necessarily $\operatorname{tr}(C_{i,j}\pi_{\hat{x}_{t-}}C_{i,j}^{*})>0$, and then $\frac{C_{i,j}\pi_{\hat{x}_{t-}}C_{i,j}^{*}}{\operatorname{tr}(C_{i,j}\pi_{\hat{x}_{t-}}C_{i,j}^{*})}=\pi_{\hat{e}_{j}}$, so that $T\leq T_{1}$. This proves (i).

Now, to prove (ii), remark that equation (1.3) can be rewritten in the form

[TABLE]

Let $T_{1}$ be defined as above; then for $t<T_{1}$ the process $(\pi_{\hat{x}_{t}})_{t}$ satisfies

[TABLE]

Starting with an initial condition $\hat{x}_{0}\in\{\hat{e}_{\ell}\}_{\ell=1}^{k}$, one proves easily that the integrand is zero, which means that $\pi_{\hat{x}_{t}}=\pi_{\hat{x}_{0}}$ for $t<T_{1}$. This shows in addition that for $t<T_{1}$, the intensity of $N_{i,j}$ is

[TABLE]

Therefore, conditionally on $x_{0}=e_{i}$, $T_{1}=\inf\{t>0;\ \exists j\neq i\mbox{ such that }N_{i,j}(t)>0\}$ and there exists an almost surely unique $j$ such that $N_{i,j}(T_{1})>0$. One then has

[TABLE]

This shows that for $t\in[0,T_{1}]$ the process $(\hat{x}_{t})_{t}$ has the same distribution as the process of equivalence classes of $X_{t}$. This extends to all $t$ by the Markov property of the Poisson processes. This proves (ii).

Points (i) and (ii) show that for $t>T_{1}$, the process $(\hat{x}_{t})_{t}$ has the same distribution as $(X_{t})_{t}$ with initial condition $X_{T_{1}}$ satisfying $\hat{X}_{T_{1}}=\hat{x}_{T_{1}}$. Therefore any invariant measure for $(\hat{x}_{t})_{t}$ is the image by $x\mapsto\hat{x}$ of an invariant measure for $(X_{t})_{t}$. Theorem 1.1 and Section 4 show that $(\hat{x}_{t})_{t}$ admits at least one invariant measure, and that the invariant measure is unique if and only if ($\mathcal{L}$-erg) holds. This implies that $(X_{t})_{t}$ has a unique invariant measure if and only if ($\mathcal{L}$-erg) holds. ∎

Acknowledgments

The research of T.B., Y.P. and C.P. has been supported by the ANR project StoQ ANR-14-CE25-0003-01. The research of T.B. has been supported by ANR-11-LABX-0040-CIMI within the program ANR-11-IDEX-0002-02. The research of M.F. was supported in part by funding from the Simons Foundation and the Centre de Recherches Mathématiques, through the Simons-CRM scholar-in-residence program. Y.P. acknowledges the support of ANR project NonStops ANR-17-CE40-0006, and of the Cantab Capital Institute for the Mathematics of Information at the University of Cambridge.