Invariant Measure for Quantum Trajectories

TL;DR

This paper investigates the long-term behavior of quantum systems under repeated measurements, establishing conditions for unique invariant measures and demonstrating geometric convergence using novel techniques.

Contribution

It introduces new conditions for the existence and uniqueness of invariant measures for quantum trajectory Markov chains and develops innovative methods for analyzing their convergence.

Findings

Unique invariant measure exists under certain conditions.

Proven geometric convergence in Wasserstein metric.

Developed new techniques beyond standard random matrix theory.

Abstract

We study a class of Markov chains that model the evolution of a quantum system subject to repeated measurements. Each Markov chain in this class is defined by a measure on the space of matrices. It is then given by a random product of correlated matrices taken from the support of the defining measure. We give natural conditions on this support that imply that the Markov chain admits a unique invariant probability measure. We moreover prove the geometric convergence towards this invariant measure in the Wasserstein metric. Standard techniques from the theory of products of random matrices cannot be applied under our assumptions, and new techniques are developed, such as maximum likelihood-type estimations.

Equations

$${\int }_{M_k(\mathbb{C})}{v}^{\ast }v\mathrm{d}\mu (v)={\mathrm{Id}}_{{\mathbb{C}}^k}.$$

$${\hat{x}}_{n+1}=V_n\cdot {\hat{x}}_n,$$

$$\Pi (\hat{x},S)={\int }_{M_k(\mathbb{C})}{1}_S(v\cdot \hat{x})\Vert vx{\Vert }^2\mathrm{d}\mu (v),$$

$${\hat{x}}_n=V_n\dots V_1\cdot {\hat{x}}_0$$

$$\log \Vert S_{n}x\Vert ,$$

$${\hat{x}}_n=S_n\cdot \hat{x}.$$

$${\rho }_n=\frac{v_n\dots v_1{\rho }_0{v}_1^{\ast }\dots v_n^{\ast }}{\mathrm{tr}(v_n\dots v_1{\rho }_0{v}_1^{\ast }\dots v_n^{\ast })}$$

$$\mathrm{tr}(v{\pi }_{{\hat{x}}_0}{v}^{\ast })=\Vert v{x}_0{\Vert }^2,\frac{v{\pi }_{{\hat{x}}_0}{v}^{\ast }}{\mathrm{tr}(v{\pi }_{{\hat{x}}_0}{v}^{\ast })}={\pi }_{v\cdot {\hat{x}}_0},$$

$$\varphi (\rho )={\int }_{M_k(\mathbb{C})}v\rho v^{\ast }\mathrm{d}\mu (v).$$

$$W_1\left(\frac{1}{m}\sum _{r=0}^{m-1}\nu {\Pi }^{mn+r},{\nu }_{\mathrm{inv}}\right)\le C{\lambda }^n,$$

$$d(\hat{x},\hat{y})=\big{(}1-|\langle x,y\rangle|^{2}\big{)}^{\frac{1}{2}},$$

$${\Pi }_s(\hat{x},S)\propto {\int }_{M_k(\mathbb{C})}{1}_S(v\cdot \hat{x})\Vert vx{\Vert }^s\mathrm{d}\mu (v).$$

$$\text{mbox}\text{hyperlink}pur(\mathbf{Pur})⟺T_{\mu }\text{mbox}iscontracting$$

$${\mathcal{J}}_n=\mathcal{B}\otimes {\mathcal{O}}_n,\mathcal{J}=\mathcal{B}\otimes \mathcal{O}.$$

$$V_i(\omega )=v_i\text{mbox}for\omega =(v_1,v_2,\dots ),$$

$$W_n=V_n\dots V_1.$$

$${\mathbb{P}}_{\nu }(S\times O_n):={\int }_{S\times O_n}\Vert W_n(\omega )x{\Vert }^2\mathrm{d}\nu (\hat{x})\mathrm{d}{\mu }^{\otimes n}(\omega ).$$

$${\mathbb{P}}_{\nu }(W_n{x}_0=0)=0.$$

$${\hat{x}}_n:=W_n\cdot \hat{x}$$

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

$${\mathbb{P}}^{\rho}(O_{n}):=\int_{O_{n}}{\operatorname{tr}}\big{(}W_{n}(\omega)\rho\hskip 0.80002ptW_{n}^{*}(\omega)\big{)}\mathrm{d}\mu^{\otimes n}(\omega).$$

$${\mathbb{P}}_{\nu }(S\times A)={\int }_S{\mathbb{P}}^{{\pi }_{\hat{x}}}(A)\mathrm{d}\nu (\hat{x}).$$

$${\mathbb{P}}_{\nu }(\mathrm{P}({\mathbb{C}}^k)\times O_n)$$

$$\displaystyle=\int_{O_{n}}{\operatorname{tr}}\big{(}W_{n}^{*}(\omega)W_{n}(\omega)\int_{{\mathrm{P}}({\mathbb{C}}^{k})}\pi_{\hat{x}}\,{\rm d}\nu(\hat{x})\big{)}\,{\rm d}\mu^{\otimes n}(\omega)$$

$$\displaystyle=\int_{O_{n}}{\operatorname{tr}}\big{(}W_{n}^{*}(\omega)W_{n}(\omega)\rho_{\nu}\big{)}\,{\rm d}\mu^{\otimes n}(\omega).$$

$${\mathbb{E}}_{\nu }({\pi }_{{\hat{x}}_1})={\mathbb{E}}_{\nu }({\pi }_{{\hat{x}}_0})={\rho }_{\nu }.$$

$${\mathbb{E}}_{\nu }({\pi }_{{\hat{x}}_1})=$$

$$=$$

$${\mathbb{P}}^{\rho }\ll {\mathbb{P}}^{\mathrm{ch}}.$$

$$M_n:=\frac{W_n^{\ast }{W}_n}{\mathrm{tr}(W_n^{\ast }{W}_n)}\text{mbox}if\mathrm{tr}(W_n^{\ast }{W}_n)\ne 0$$

$$W_n=U_n\sqrt{\mathrm{tr}(W_n^{\ast }{W}_n)}{M}_n^{\frac{1}{2}}.$$

Full text

1. 1 Introduction

2. 2 Uniqueness of the invariant measure

3. 3 Convergence

4. 3.1 Convergence for $\mathcal{O}$-measurable random variables

5. 3.2 Convergence to a ${\mathcal{O}}$-measurable process

6. 3.3 Convergence in Wasserstein metric

7. 4 Lyapunov exponents

8. A Equivalence of (Pur) and contractivity

9. B Set of invariant measures under assumption (Pur)

10. C Products of special unitary matrices

1. Introduction

We consider a complex vector space $\mathbb{C}^{k}$ and its projective space ${\mathrm{P}}(\mathbb{C}^{k})$ equipped with its Borel $\sigma$-algebra $\mathcal{B}$. For a non zero vector $x\in\mathbb{C}^{k}$, we denote $\hat{x}$ the corresponding equivalence class of $x$ in ${\mathrm{P}}(\mathbb{C}^{k})$. For a linear map $v\in M_{k}(\mathbb{C})$ we denote $v\cdot\hat{x}$ the element of the projective space represented by $v\,x$ whenever $v\,x\neq 0$. We equip $M_{k}(\mathbb{C})$ with its Borel $\sigma$-algebra and let $\mu$ be a measure on $M_{k}(\mathbb{C})$ with a finite second moment, $\int_{M_{k}(\mathbb{C})}\|v\|^{2}\,{\rm d}\mu(v)<\infty$, that satisfies the stochasticity condition

[TABLE]

In this article we are interested in particular Markov chains $(\hat{x}_{n})$ on ${\mathrm{P}}(\mathbb{C}^{k})$, defined by

[TABLE]

where $V_{n}$ is a $M_{k}(\mathbb{C})$-valued random variable with a probability density $||vx_{n}||^{2}/||x_{n}||^{2}\mathrm{d}\mu(v).$ More precisely, such a Markov chain is associated with the transition kernel given for a set $S\in\mathcal{B}$ and $\hat{x}\in{\mathrm{P}}(\mathbb{C}^{k})$ by

[TABLE]

where $x$ is an arbitrary normalized vector representative of $\hat{x}$. Note that the normalization condition (1) imposed on $\mu$ is equivalent to the conservation of probability, \Pi\big{(}\hat{x},{\mathrm{P}}(\mathbb{C}^{k})\big{)}=1. We are interested in the large-time distribution of $(\hat{x}_{n})$.

Note that $\hat{x}_{n}$ can be written as

[TABLE]

so that the study of $\hat{x}_{n}$ can be formulated in terms of random products of matrices. Markov chains associated to random products of matrices were studied in a more general setting where the weight appearing in the transition kernel (2) is proportional to $\|vx\|^{s}$ for some $s\geq 0$, instead of $\|vx\|^{2}$. The classical case of products of independent, identically distributed random matrices pioneered by Kesten, Furstenberg and Guivarc’h corresponds to $s=0$. In that case, for i.i.d. invertible random matrices $Y_{1},Y_{2},\ldots$, denoting $S_{n}=Y_{n}\ldots Y_{1}$, one is usually interested in the asymptotic properties of

[TABLE]

for any $x\neq 0$. In particular, a law of large numbers, a central limit theorem and a large deviation principle have been obtained for this quantity, under contractivity and strong irreducibility assumptions [8, 11, 14]. Such results are closely linked to the uniqueness of the invariant measure of the Markov chain

[TABLE]

These results were generalized to the case $s>0$ in [10]. Our framework corresponds to the case $s=2$; in this case, and with the additional assumption (1), we provide a new method to study this Markov chain, and use it to derive the above results without assuming invertibility of the matrices, and with an optimal irreducibility assumption. We compare our approach with respect to that of [10] at the end of this section.

The method that we employ is motivated by an interpretation of this process as statistics of a quantum system being repeatedly indirectly measured. Let us expand on this as we introduce more notation and terminology. The set of states of a quantum system described by a finite dimensional Hilbert space ${\mathbb{C}}^{k}$ is the set of density matrices $\mathcal{D}_{k}:=\{\rho\in M_{k}({\mathbb{C}})\ |\ \rho\geq 0,\ {\operatorname{tr}}\,\rho=1\}$. This set is convex and the set of its extreme points is called the set of pure states. This latter set is in one to one correspondence with the projective space ${\mathrm{P}}({\mathbb{C}}^{k})$ by the bijection ${\mathrm{P}}({\mathbb{C}}^{k})\ni\hat{x}\mapsto\pi_{\hat{x}}\in\mathcal{D}_{k}$ with $\pi_{\hat{x}}$ the orthogonal projector on the corresponding ray in ${\mathbb{C}}^{k}$.

The time evolution of the system conditioned on a measurement outcome is encoded in a matrix $v$ that updates the state of the system. The support of $\mu$ is endowed with the meaning of the possible updates, and the system is updated according to $v$ with a probability density ${\operatorname{tr}}(v\rho v^{*})\,{\rm d}\mu(v)$. Given $v$, a state $\rho$ is mapped to a state $v\rho v^{*}/{\operatorname{tr}}(v\rho v^{*})$. Iterating this procedure defines a random sequence $(\rho_{n})$ in $\mathcal{D}_{k}$ called a quantum trajectory: after $n$ measurements with resulting matrices $v_{1},\dots,v_{n}$ the state of the system becomes

[TABLE]

where $(v_{1},\ldots,v_{n})$ has probability density ${\operatorname{tr}}(v_{n}\dots v_{1}\rho_{0}v_{1}^{*}\ldots v_{n}^{*})\,{\rm d}\mu^{\otimes n}(v_{1},\ldots,v_{n})$. In other words, the process Eq. 3 describes an evolution of a repeatedly measured quantum system.

A key result in the theory of quantum trajectories is a purification theorem obtained by Kümmerer and Maassen [15] showing that quantum trajectories $(\rho_{n})$ defined on $\mathcal{D}_{k}$ almost surely approach the set of pure states, the extreme points of $\mathcal{D}_{k}$ if and only if the following purification condition is satisfied:

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

Any orthogonal projector $\pi$ such that for any $n\in{\mathbb{N}}$, $\pi v_{1}^{*}\ldots v_{n}^{*}v_{n}\ldots v_{1}\pi\propto\pi$ for $\mu^{\otimes n}$-almost all $(v_{1},\ldots,v_{n})$, is of rank one

(we write $X\propto Y$ for $X,Y$ two operators if there exists $\lambda\in{\mathbb{C}}$ such that $X=\lambda Y$).

Under this assumption, the long-time behavior of the Markov chain is essentially dictated by its form on the set of pure states, i.e. for $\rho_{0}=\pi_{\hat{x}_{0}}$. It is an immediate observation that

[TABLE]

for all $v\in M_{k}(\mathbb{C})$. This way our Markov chain $(\hat{x}_{n})$ corresponds to the quantum trajectory $(\rho_{n})$ described above when $\rho_{0}$ is a pure state $\pi_{\hat{x}_{0}}$.

Although ideas underlying our method are based on the connection of $(\hat{x}_{n})$, with this physical problem, we will not explicitly use it in the paper. The notion of quantum trajectory originates in quantum optics [6], and Haroche’s Nobel prize winning experiment [9] is arguably the most prominent example of a system described by the above formalism. The reader interested in the involved mathematical structures might consult for example the review book [13] or the pioneering article [15].

We will show that under condition (Pur), the set of all invariant measures of the Markov chain (3) can be completely classified, depending on the operator $\phi$ on $\mathcal{D}_{k}$ describing the average evolution:

[TABLE]

The map $\phi$ on $\mathcal{D}_{k}$ is completely positive and trace-preserving.111Complete positivity is stronger than positivity; namely $\phi$ is completely positive if $\phi\otimes{\rm Id}_{M_{n}({\mathbb{C}})}$ is positive for all $n\in{\mathbb{N}}$. Such a map is often called a quantum channel (see e.g.  [20]). It has in particular the property of mapping states to states. Brouwer’s fixed point Theorem shows that there exists an invariant state, i.e. $\rho\in\mathcal{D}_{k}$ such that $\phi(\rho)=\rho$. A necessary and sufficient algebraic condition for uniqueness of this invariant state is (see e.g. [7, 20, 5])

**($\phi$-Erg)****: **

There exists a unique minimal non trivial $\operatorname{supp}\mu$-invariant subspace $E$ of ${\mathbb{C}}^{k}$.

If ($\phi$-Erg) holds with $E={\mathbb{C}}^{k}$, then $\phi$ is said irreducible. We chose the name ($\phi$-Erg) to avoid confusion with the notion of irreducibility for Markov chains. We moreover emphasize that we call this assumption ($\phi$-Erg) because it relies only on $\phi$ and not on the different operators $v$ in the support of $\mu$: an equivalent statement of ($\phi$-Erg) is that there exists a unique minimal nonzero orthogonal projector $\pi$ such that $\phi(\pi)\leq\lambda\pi$ for some $\lambda\geq 0$ (see e.g. [18]).

We now state the main result of the paper:

Theorem 1.1**.**

Assume that $\mu$ satisfies assumptions (Pur) and ($\phi$-Erg). Then, the transition kernel $\Pi$ has a unique invariant probability measure $\nu_{\mathrm{inv}}$ and there exists $m\in\{1,\ldots,k\}$, $C>0$ and $0<\lambda<1$ such that for any probability measure $\nu$ over \big{(}{\mathrm{P}}({\mathbb{C}}^{k}),\mathcal{B}\big{)},

[TABLE]

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

The Wasserstein metric is constructed with respect to a natural metric on the complex projective space. This metric is defined, for $\hat{x},\,\hat{y}$ in ${\mathrm{P}}({\mathbb{C}}^{k})$, by

[TABLE]

where $x,\,y$ are unit length representative vectors and $\langle\,\cdot\,,\,\cdot\,\rangle$ is the canonical hermitian inner product on $\mathbb{C}^{k}$.

Let us now compare our results to those of the article [10] of Guivarc’h and Le Page. They consider a probability distribution $\mu$ with support in $GL_{k}(\mathbb{C})$, without requiring the normalization condition (1), and study the transition kernel on ${\mathrm{P}}({\mathbb{C}}^{k})$ given, for $S\in\mathcal{B}$, by

[TABLE]

In the case $s=2$, Theorem A of [10] implies the conclusions of Theorem 1.1 under two assumptions:

• •

strong irreducibility, in the sense that there is no non-trivial finite union of proper subspaces of ${\mathbb{C}}^{k}$ left invariant by all $v\in\mathrm{supp}\,\mu$,

• •

contractivity, in the sense that there exists a sequence $(a_{n})$ in $T_{\mu}$, the smallest closed sub-semigroup of $GL_{k}({\mathbb{C}})$ containing $\operatorname{supp}\mu$, such that $\lim_{n\to\infty}a_{n}/\|a_{n}\|$ exists and is of rank one.

It is, however, a simple exercise to prove that strong irreducibility of $\mu$ implies ($\phi$-Erg) with $E={\mathbb{C}}^{k}$. In addition, if we assume $\operatorname{supp}\mu\subset GL_{k}(\mathbb{C})$ and $\operatorname{supp}\mu$ is strongly irreducible, the equivalence

[TABLE]

holds (see Appendix A). Our results therefore offer a strong refinement of [10] in the restricted framework of $s=2$ with $\int v^{*}v\,{\rm d}\mu(v)={\rm Id}_{{\mathbb{C}}^{k}}$.

The article is structured as follows. Section 2 is devoted to the first part of Theorem 1.1, that is the uniqueness of the invariant measure. In Section 3 we show the geometric convergence towards the invariant measure with respect to the $1$-Wasserstein metric. In Section 4 we discuss the Lyapunov exponents of the process and relate them to the convergence between the Markov chain and an estimate of the chain used in our proofs.

Notation

In all of the following, for $x\in{\mathbb{C}}^{k}\setminus\{0\}$, $\hat{x}$ is its equivalence class in ${\mathrm{P}}({\mathbb{C}}^{k})$ and, for $\hat{x}$ in ${\mathrm{P}}({\mathbb{C}}^{k})$, $x$ is an arbitrary norm one vector representative of $\hat{x}$. If e.g. ${\mathbb{P}}_{\nu}$ (resp. $\mathbb{P}^{\rho}$) is a probability measure (depending on some a priori object $\nu$ (resp. $\rho$)) then ${\mathbb{E}}_{\nu}$ (resp. ${\mathbb{E}}^{\rho}$) is the expectation with respect to ${\mathbb{P}}_{\nu}$ (resp. ${\mathbb{P}}^{\rho}$). The set ${\mathbb{N}}$ represents the set of positive integers $\{1,2,\ldots\}$.

2. Uniqueness of the invariant measure

This section concerns essentially the first part of Theorem 1.1. More precisely, under ($\phi$-Erg) and (Pur) we show that the Markov chain has a unique invariant measure. We note that an invariant measure always exists since ${\mathrm{P}}(\mathbb{C}^{k})$ is compact.

We now proceed to introduce some additional notation. We consider the space of infinite sequences $\Omega:=M_{k}({\mathbb{C}})^{\mathbb{N}}$, write $\omega=(v_{1},v_{2},\dots)$ for any such infinite sequence, and denote by $\pi_{n}$ the canonical projection on the first $n$ components, $\pi_{n}(\omega)=(v_{1},\ldots,v_{n})$. Let $\mathcal{M}$ be the Borel $\sigma$-algebra on $M_{k}({\mathbb{C}})$. For $n\in{\mathbb{N}}$, let $\mathcal{O}_{n}$ be the $\sigma$-algebra on $\Omega$ generated by the $n$-cylinder sets, i.e. $\mathcal{O}_{n}=\pi_{n}^{-1}(\mathcal{M}^{\otimes n})$. We equip the space $\Omega$ with the smallest $\sigma$-algebra $\mathcal{O}$ containing $\mathcal{O}_{n}$ for all $n\in{\mathbb{N}}$. We let $\mathcal{B}$ be the Borel $\sigma$-algebra on ${\mathrm{P}}(\mathbb{C}^{k})$, and denote

[TABLE]

This makes \big{(}{\mathrm{P}}(\mathbb{C}^{k})\times\Omega,\mathcal{J}\big{)} a measurable space. With a small abuse of notation we denote the sub-$\sigma$-algebra $\{\emptyset,{\mathrm{P}}(\mathbb{C}^{k})\}\times\mathcal{O}$ by $\mathcal{O}$ and equivalently identify any $\mathcal{O}$-measurable function $f$ with a $\mathcal{J}$-measurable function $f$ satisfying $f(\hat{x},\omega)=f(\omega)$.

For $i\in\mathbb{N}$, we consider the random variables $V_{i}:\Omega\mapsto M_{k}(\mathbb{C})$,

[TABLE]

and we introduce $\mathcal{O}_{n}$-mesurable random variables $(W_{n})$ defined for all $n\in\mathbb{N}$ as

[TABLE]

With a small abuse of notation we identify cylinder sets and their bases, and extend this identification to several associated objects. In particular we identify $O_{n}\in\mathcal{M}^{\otimes n}$ with $\pi_{n}^{-1}(O_{n})$, a function $f$ on $\mathcal{M}^{\otimes n}$ with $f\circ\pi_{n}$ and a measure $\mu^{\otimes n}$ with the measure $\mu^{\otimes n}\circ\pi_{n}$. Since $\mu$ is not necessarily finite, note that we can not extend $(\mu^{\otimes n})$ into a measure on $\Omega$.

Let $\nu$ be a probability measure over $({\mathrm{P}}({\mathbb{C}}^{k}),\mathcal{B})$. We extend it to a probability measure $\mathbb{P}_{\nu}$ over $({\mathrm{P}}({\mathbb{C}}^{k})\times\Omega,\mathcal{J})$ by letting, for any $S\in\mathcal{B}$ and any cylinder set $O_{n}\in\mathcal{O}_{n}$,

[TABLE]

From relation (1), it is easy to check that the expression (9) defines a consistent family of probability measures and, by Kolmogorov’s Theorem, this defines a unique probability measure ${\mathbb{P}}_{\nu}$ on ${\mathrm{P}}({\mathbb{C}}^{k})\times\Omega$. In addition, the restriction of $\mathbb{P}_{\nu}$ to $\mathcal{B}\otimes\{\emptyset,\Omega\}$ is by construction $\nu$.

We now define the random process $(\hat{x}_{n})$. For $(\hat{x},\omega)\in{\mathrm{P}}({\mathbb{C}}^{k})\times\Omega$ we define $\hat{x}_{0}(\hat{x},\omega)=\hat{x}$. Note that for any $n$, the definition (9) of ${\mathbb{P}}_{\nu}$ imposes

[TABLE]

This allows us to define a sequence $(\hat{x}_{n})_{n\in\mathbb{N}}$ of $(\mathcal{J}_{n})$-adapted random variables on the probability space $({\mathrm{P}}({\mathbb{C}}^{k})\times\Omega,\mathcal{J},{\mathbb{P}}_{\nu})$ by letting

[TABLE]

whenever the expression makes sense, i.e. for any $\omega$ such that $W_{n}(\omega)x\neq 0$, and extending it arbitrarily to the whole of $\Omega$. The process $(\hat{x}_{n})$ on $(\Omega\times{\mathrm{P}}(\mathbb{C}^{k}),\mathcal{J},\mathbb{P}_{\nu})$ has the same distribution as the Markov chain defined by $\Pi$ and initial probability measure $\nu$.

Let us highlight the relation between ${\mathbb{P}}_{\nu}$ and density matrices. To that end, let

[TABLE]

By linearity and positivity of the expectation, $\rho_{\nu}\in\mathcal{D}_{k}$. Note that, conversely, for a given $\rho\in\mathcal{D}_{k}$ there exists $\nu$ (in general non-unique) such that $\rho_{\nu}=\rho$. For example, if a spectral decomposition of $\rho$ is $\rho=\sum_{j}p_{j}\pi_{x_{j}}$ then necessarily $\sum_{j}p_{j}=1$, so that $\nu=\sum_{j}p_{j}\delta_{\hat{x}_{j}}$ is a probability measure on ${\mathrm{P}}({\mathbb{C}}^{k})$, and it satisfies the desired relation (11).

This relation motivates the following definition of probability measures over $(\Omega,\mathcal{O})$. For $\rho\in\mathcal{D}_{k}$ and any cylinder set $O_{n}\in\mathcal{O}_{n}$, let

[TABLE]

In particular, for any $S\in\mathcal{B}$ and $A\in\mathcal{O}$,

[TABLE]

The following proposition elucidates further the connections between ${\mathbb{P}}_{\nu}$ and ${\mathbb{P}}^{\rho_{\nu}}$.

Proposition 2.1**.**

The marginal of $\mathbb{P}_{\nu}$ on $\mathcal{O}$ is the probability measure $\mathbb{P}^{\rho_{\nu}}$.

Moreover, if ($\phi$-Erg) holds, ${\mathbb{P}}^{\rho_{\nu_{a}}}={\mathbb{P}}^{\rho_{\nu_{b}}}$ for any two $\Pi$-invariant probability measures $\nu_{a}$ and $\nu_{b}$.

Proof.

By construction it is sufficient to check the equality of the measures on cylinder sets. Let $O_{n}\in\mathcal{O}_{n}$; from the definition of $\mathbb{P}_{\nu}$, and trace and integral linearity we have

[TABLE]

The equality between the marginal of ${\mathbb{P}}_{\nu}$ on $\mathcal{O}$ and ${\mathbb{P}}^{\rho_{\nu}}$ follows.

If $\nu$ is an invariant measure, on the one hand

[TABLE]

On the other hand,

[TABLE]

so that $\rho_{\nu}$ is a fixed point of $\phi$.

Hence if ($\phi$-Erg) holds, $\rho_{\nu}$ is the unique fixed point of $\phi$ in $\mathcal{D}_{k}$. Hence $\rho_{\nu_{a}}=\rho_{\nu_{b}}$ and ${\mathbb{P}}^{\rho_{\nu_{a}}}={\mathbb{P}}^{\rho_{\nu_{b}}}$ holds. ∎

In the following we use the measure ${\mathbb{P}}^{\mathrm{ch}}={\mathbb{P}}^{\frac{1}{k}{\rm Id}_{{\mathbb{C}}^{k}}}$ associated to the operator ${\rm Id}_{\mathbb{C}^{k}}/k\in\mathcal{D}_{k}$ as a reference measure. Since for any $\rho\in\mathcal{D}_{k}$ there exists a constant $c$ such that $\rho\leq c\frac{{\rm Id}_{\mathbb{C}^{k}}}{k}$ we have

[TABLE]

The Radon–Nykodim derivative will be made explicit in Proposition 2.2. To that end, we use a particular $(\mathcal{O}_{n})$-adapted process. We define a sequence of matrix-valued random variables:

[TABLE]

and extend the definition arbitrarily whenever ${\operatorname{tr}}(W_{n}^{*}W_{n})=0$. The latter alternative appears with probability 0. Indeed {\mathbb{P}}^{\mathrm{ch}}\big{(}{\operatorname{tr}}(W_{n}^{*}W_{n})=0\big{)}=0 and then by the absolute continuity of ${\mathbb{P}}^{\rho}$ with respect to ${\mathbb{P}}^{\mathrm{ch}}$ we have {\mathbb{P}}_{\nu}\big{(}{\operatorname{tr}}(W_{n}^{*}W_{n})=0\big{)}={\mathbb{P}}^{\rho_{\nu}}\big{(}{\operatorname{tr}}(W_{n}^{*}W_{n})=0\big{)}=0 for any measure $\nu$. The key property of $M_{n}$, that we establish in the proof of Proposition 2.2, is that it is a $(\mathcal{O}_{n})$-martingale w.r.t. ${\mathbb{P}}^{\mathrm{ch}}$.

From the existence of a polar decomposition for $W_{n}$, for each $n$, there exists a unitary matrix-valued random variable $U_{n}$ such that

[TABLE]

This process $(U_{n})$ can be chosen to be $(\mathcal{O}_{n})$-adapted.

The key technical results about $M_{n}$ needed for our proofs are summarized in the following proposition.

Proposition 2.2**.**

For any probability measure $\nu$ over $({\mathrm{P}}({\mathbb{C}}^{k}),\mathcal{B})$, $(M_{n})$ converges $\mathbb{P}_{\nu}\operatorname{-}\mathrm{a.s.}$ and in $L^{1}$-norm to an $\mathcal{O}$-measurable random variable $M_{\infty}$. The change of measure formula

[TABLE]

holds true for all $\rho\in\mathcal{D}_{k}$.

Moreover, the measure $\mu$ verifies (Pur) if and only if $M_{\infty}$ is ${\mathbb{P}}_{\nu}\operatorname{-}\mathrm{a.s.}$ a rank one projection for any probability measure $\nu$ over $({\mathrm{P}}({\mathbb{C}}^{k}),\mathcal{B})$.

Proof.

We start the proof by showing that $M_{n}$ is a ${\mathbb{P}}^{\mathrm{ch}}$-martingale. Recall that for all $n\in\mathbb{N}$ and all $O_{n}\in\mathcal{O}_{n}$,

[TABLE]

From the definition of $W_{n}$, Eq. 8,

[TABLE]

This implies that for an arbitrary $\mathcal{O}_{n}$-measurable random variable $Y$

[TABLE]

where the second equality follows from the stochasticity condition (1), $\int v^{*}vd\mu(v)={\rm Id}_{{\mathbb{C}}_{k}}$. This shows that $(M_{n})$ is a $(\mathcal{O}_{n})$-martingale w.r.t. ${\mathbb{P}}^{\mathrm{ch}}$. Since the sequence $(M_{n})$ is composed of positive semidefinite matrices of trace one, its coordinates are a.s. uniformly bounded by $1$. Therefore, the martingale property implies the $L^{1}$ and a.s. convergence of $(M_{n})$ to a $\mathcal{O}$-measurable random variable $M_{\infty}$. Now note that for any $\rho\in\mathcal{D}_{k}$,

[TABLE]

This way, the convergence of $(M_{n})$ implies the change of measure formula.

We now prove the last part of the Proposition. Using the martingale property one can see that for all $n\in\mathbb{N}$, and any fixed $p\in\mathbb{N}$,

[TABLE]

Since $(M_{n})$ is bounded and almost surely convergent, the Lebesgue dominated convergence Theorem implies that the term $V_{n}^{p}$ is convergent when $n$ goes to infinity. Hence we get that

[TABLE]

which yields that

[TABLE]

At this stage we use the polar decomposition of $(W_{n})$, Eq. (14), to write

[TABLE]

Then we get an expression for the conditional expectation, c.f. the first part of the proof,

[TABLE]

Since $(U_{n})$ are unitary matrices, ${\mathbb{P}}^{\mathrm{ch}}\operatorname{-}\mathrm{a.s.}$, there exists a subsequence along which $(U_{n})$ converges to some random unitary operator $U_{\infty}$. Taking the limit as $n$ goes to infinity along this subsequence in the above expression, yields

[TABLE]

This implies that

[TABLE]

${\mathbb{P}}^{\mathrm{ch}}$-almost surely and for $\mu^{\otimes p}$-almost all $(v_{1},\ldots,v_{p})$. Note that the notion “${\mathbb{P}}^{\mathrm{ch}}$-almost surely” stands for $M_{\infty}$ and $U_{\infty}$ whereas “$\mu^{\otimes p}$-almost all” stands for the $v_{i}^{\prime}s$. The product is zero, only if at least one of the terms in the product is zero. It turns out, however, that both cases can be described by a single condition: the product vanishes only if there exists $\lambda$ such that

[TABLE]

Denoting by $\pi_{\infty}$ the orthogonal projector onto the range of $M_{\infty}$, the condition is equivalent to $\pi_{\infty}U_{\infty}^{*}v_{1}^{*}\ldots v_{p}^{*}v_{p}\ldots v_{1}U_{\infty}\pi_{\infty}=\lambda\pi_{\infty}$. Finally, it follows that

[TABLE]

for $\mu^{\otimes p}$-almost all $(v_{1},\ldots,v_{p})$. Since $U_{\infty}\pi_{\infty}U_{\infty}^{*}$ is an orthogonal projector, the condition (Pur) implies that $\operatorname{rank}(M_{\infty})=\operatorname{rank}(U_{\infty}\pi_{\infty}U_{\infty}^{*})=1$. Since $M_{\infty}$ is a trace one, positive semidefinite matrix this means that $M_{\infty}$ is a rank one projector.

For the converse implication, assume that $M_{\infty}$ is ${\mathbb{P}}^{\mathrm{ch}}$-almost surely a rank one projection but that (Pur) does not hold. Then there exists $\pi$, a rank two orthogonal projector, s.t. for all $n\in{\mathbb{N}}$,

[TABLE]

$\mu^{\otimes n}$-almost everywhere. Since $\mu^{\otimes n}$-almost everywhere $M_{n}\propto W_{n}^{*}W_{n}$, we get that

[TABLE]

$\mu^{\otimes n}$-almost everywhere. Thus, $\pi M_{\infty}\pi\propto\pi,$ and under our assumption that $\operatorname{rank}M_{\infty}=1$ ${\mathbb{P}}^{\mathrm{ch}}\operatorname{-}\mathrm{a.s.}$ and $\operatorname{rank}\pi=2$ this implies that $\pi M_{\infty}\pi=0$, ${\mathbb{P}}^{\mathrm{ch}}$-almost surely. On the other hand for all $n\in\mathbb{N}$ we have ${\mathbb{E}}^{\mathrm{ch}}(M_{n})={\rm Id}_{{\mathbb{C}}^{k}}$, and the $L^{1}$ convergence implies that ${\mathbb{E}}^{\mathrm{ch}}(M_{\infty})={\rm Id}_{{\mathbb{C}}^{k}}$. Then, ${\mathbb{E}}^{\mathrm{ch}}(\pi M_{\infty}\pi)=\pi$ which contradicts $\pi M_{\infty}\pi=0$ ${\mathbb{P}}^{\mathrm{ch}}\operatorname{-}\mathrm{a.s.}$ ∎

By the polar decomposition, the rank of $W_{n}$ is equal to the rank of $M_{n}$ and the proposition thus implies that $W_{n}\rho_{0}W_{n}^{*}/{\operatorname{tr}}(W_{n}\rho_{0}W_{n}^{*})$ approaches the set of pure states for any $\rho_{0}\in\mathcal{D}_{k}$ if and only if (Pur) holds. This is the result of Maassen and Kümmerer [15] mentioned in the introduction. Though $M_{n}$ is not used in [15], the proof relies on similar ideas.

We are now in the position to show that the Markov chain $(\hat{x}_{n})$ is asymptotically an $\mathcal{O}$-measurable process. This is expressed in the following lemma. Whenever (Pur) holds, we denote by $\hat{z}\in{\mathrm{P}}({\mathbb{C}}^{k})$ the $\mathcal{O}$-measurable random variable defined by

[TABLE]

Recall that $d(\cdot,\cdot)$, defined by Eq. (7), is our metric on ${\mathrm{P}}({\mathbb{C}}^{k})$.

Lemma 2.3**.**

Assume (Pur) holds. Then for any probability measure $\nu$ on $({\mathrm{P}}({\mathbb{C}}^{k}),\mathcal{B})$,

[TABLE]

Proof.

We start the proof by showing that for any $\nu$

[TABLE]

Let $\hat{x}$ be fixed and recall from Proposition 2.2 that (Pur) implies $M_{\infty}=\pi_{\hat{z}}$. In order to show (18), it is enough to show that $\hat{x}$ is almost surely not orthogonal to $\hat{z}$. From Eq. 13 and the change of measure formula in Proposition 2.2,

[TABLE]

Hence the event $\{d(\hat{x},\hat{z})=1\}$ has $\mathbb{P}_{\nu}$-measure [math], and (18) follows from the almost sure convergence of $M_{n}$ to $\pi_{\hat{z}}$.

Now using the polar decomposition, Eq. (14), and the fact that proportionality of vectors means equality of their equivalence classes in ${\mathrm{P}}(\mathbb{C}^{k})$, we have

[TABLE]

The first part of the proof then yields,

[TABLE]

∎

The uniqueness of the invariant measure which is the first part of Theorem 1.1 follows as a corollary.

Corollary 2.4**.**

Assume (Pur) and ($\phi$-Erg). Then the Markov kernel $\Pi$ admits a unique invariant probability measure.

Proof.

For an invariant measure $\nu$, the random variable $\hat{x}_{n}$ is $\nu$-distributed for all $n\in\mathbb{N}$. In particular, \mathbb{E}_{\nu}\big{(}f(\hat{x}_{n})\big{)} is constant for any continuous function $f$. On the other hand Lemma 2.3 and Lebesgue’s dominated convergence theorem imply that

[TABLE]

Hence

[TABLE]

Assume now that there exist two invariant measures $\nu_{a}$ and $\nu_{b}$. Since $U_{n}\cdot\hat{z}$ is $\mathcal{O}$-measurable, Proposition 2.1 implies

[TABLE]

Then Eq. 19 applied with $\nu=\nu_{a}$, resp. $\nu=\nu_{b}$ gives

[TABLE]

which means that $\nu_{a}=\nu_{b}$ and the uniqueness is proved. ∎

Assuming only (Pur) we can actually completely characterize the set of invariant measures.

Proposition 2.5**.**

Assuming (Pur) there exists a set $\{F_{j}\}_{j=1}^{d}$ of mutually orthogonal subspaces of ${\mathbb{C}}^{k}$ such that, for each $j\in\{1,\ldots,d\}$ there exists a unique $\Pi$-invariant probability measure $\nu_{j}$ supported on ${\mathrm{P}}(F_{j})$, and the set of $\Pi$-invariant probability measures is the convex hull of $\{\nu_{j}\}_{j=1}^{d}$.

The subspaces $F_{j}$ are the ranges of the extremal fixed points of $\phi$ in $\mathcal{D}_{k}$. This is shown in the proof of this Proposition, that is provided in Appendix B.

Remark 2.6**.**

Assuming ($\phi$-Erg) only, the chain might or might not have a unique invariant probability measure. Indeed if $\operatorname{supp}\mu\subset SU(k)$ Assumption (Pur) is trivially not verified and, as proved in Appendix C, the uniqueness of the invariant measure depends on the smallest closed subgroup of $SU(k)$ containing $\operatorname{supp}\mu$. To illustrate this point, in the same appendix, we study two examples with $\mu$ supported on and giving equiprobability to two elements of $SU(2)$ such that ($\phi$-Erg) holds. In the first example $\Pi$ has a unique invariant probability measure whereas in the second example $\Pi$ has uncountably many mutually singular invariant probability measures.

3. Convergence

We now turn to the proof of the second part of Theorem 1.1, namely the geometric convergence in Wasserstein distance of the process $(\hat{x}_{n})$ towards the invariant measure. We first recall a 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 Kantorovich–Rubinstein duality Theorem as

[TABLE]

where ${\rm 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$.

The proof of Eq. 6 consists of three parts. In the first part we show a geometric convergence in total variation of ${\mathbb{P}}^{\rho}$ to ${\mathbb{P}}^{\rho_{\mathrm{inv}}}$ under the shift $\theta(v_{1},v_{2},\ldots)=(v_{2},v_{3},\ldots)$. In the second one we show a geometric convergence of the chain $(\hat{x}_{n})$ towards an $\mathcal{O}$-measurable process $(\hat{y}_{n})$. Finally, we combine these results to prove Eq. 6.

3.1. Convergence for $\mathcal{O}$-measurable random variables

Let us first discuss the origin of the integer $m$ in Eq. 6. Let $(E_{1},\ldots,E_{\ell})$ be an orthogonal partition of a $\operatorname{supp}\mu$-invariant subspace, i.e.  a family of mutually orthogonal subspaces such that $E_{1}\oplus\ldots\oplus E_{\ell}$ is a $\operatorname{supp}\mu$-invariant subspace. We say that $(E_{1},\ldots,E_{\ell})$ is a $\ell$-cycle of $\phi$ if $vE_{j}\subset E_{j+1}$ for $\mu$-a.a. $v$ (with the convention $E_{\ell+1}=E_{1}$).222As suggested by its name, the notion of cycle for $\phi$ depends only on $\phi$ and not on the specific measure $\mu$ leading to $\phi$ [20, 18, 5]. The set of $\ell\in{\mathbb{N}}$ for which there exists an orthogonal partition $(E_{1},\ldots,E_{\ell})$ is non-empty (as it contains $1$) and bounded (as necessarily $\ell\leq k$).

Definition 3.1**.**

The largest $\ell\in{\mathbb{N}}$ such that there exists a $\ell$-cycle of $\phi$ is called the period of $\phi$. We denote this period by $m$.

Remark 3.2**.**

• •

The above definition of period for $\phi$ is similar to that of the period of a $\varphi$-irreducible Markov chain. It is obvious that if $(E_{1},\ldots,E_{\ell})$ is an $\ell$-cycle of $\phi$ then it is also an $\ell$-cycle of $\Pi$. However, the Markov chain defined by $\Pi$ is not $\varphi$-irreducible in general. Hence the results of **[17]** on the period of $\varphi$-irreducible Markov chains do not apply and the characterization of the period of $\Pi$ remains an open problem.

• •

The above definition shows that the union $\bigcup_{j=1}^{m}E_{j}$ is invariant by $\mu$-a.a. $v$. Hence, the strong irreducibility assumption discussed at the end of the introduction implies that $m=1$.

The following result is a reformulation of the Perron–Frobenius theorem of Evans and Høegh-Krohn, [7]. The original formulation in [7] makes the additional assumption that $E={\mathbb{C}}^{k}$ in ($\phi$-Erg). For the present extension see e.g. [20]. In the following statement, and in the rest of the article, for $X$ an operator on ${\mathbb{C}}^{k}$ we denote $\|X\|_{1}={\operatorname{tr}}|X|$ (all statements are identical with a different norm, but this choice will spare us a few irrelevant constants).

Theorem 3.3**.**

Assume ($\phi$-Erg) holds. Then there exists a unique $\phi$-invariant element $\rho_{\mathrm{inv}}$ of $\mathcal{D}_{k}$ with range equal to the minimal invariant subspace $E$. In addition, there exist two positive constants $c$ and $\lambda<1$ such that, with $m$ defined in Definition 3.1, for any $\rho\in\mathcal{D}_{k}$ and for any $n\in{\mathbb{N}}$,

[TABLE]

Proof.

Theorem 4.2 in [7] implies that $\rho_{\mathrm{inv}}$ is the unique $\phi$-invariant element of $\mathcal{D}_{k}$, that the eigenvalues of $\phi$ of modulus one are exactly the $m$-th roots of unity, and that they are all simple. The statement follows, with $\lambda$ the modulus of the largest non-peripheral eigenvalue. ∎

Recall that $\theta$ is the left shift operator on  $\Omega$, i.e.

[TABLE]

The main result of this section is the following proposition. As announced it concerns the speed of convergence in total variation (expressed in terms of expectation values).

Proposition 3.4**.**

Assume ($\phi$-Erg) holds. Then there exist two positive constants $C$ and $\lambda<1$ such that for any $\mathcal{O}$-measurable function $f$ with essential bound $\|f\|_{\infty}$, any $\rho\in\mathcal{D}_{k}$ and all $n\in{\mathbb{N}}$,

[TABLE]

Proof.

We claim that for any bounded $\mathcal{O}$-measurable function $f$,

[TABLE]

It suffices to prove the relation for all $\mathcal{O}_{l}$-measurable functions for some integer $l$. For such a function,

[TABLE]

which is equal to ${\mathbb{E}}^{\phi(\rho)}(f)$ as claimed.

Applying Eq. 22 multiple times and using the change of measure of Proposition 2.2 we obtain

[TABLE]

for any $\mathcal{O}$-measurable function $f$. Using $|{\operatorname{tr}}(M_{\infty}A)|\leq\|A\|_{1}$ for $A=A^{*}$ (remark that $M_{\infty}\in\mathcal{D}_{k}$ by construction) we then obtain

[TABLE]

and Theorem 3.3 yields the proposition with $C=ck$. ∎

3.2. Convergence to a ${\mathcal{O}}$-measurable process

Let us introduce two relevant processes: for all $n\in{\mathbb{N}}$, let

[TABLE]

and

[TABLE]

Both random variables $\hat{y}_{n}$ and $\hat{z}_{n}$ are $\mathcal{O}_{n}$-measurable.

The random variable $\hat{z}_{n}$ corresponds to the maximum likelihood estimator of $\hat{x}_{0}$. Note that the $\operatorname{argmax}$ may not be uniquely defined. We can, however, define it in an $\mathcal{O}_{n}$-measurable way. The following results will not be affected by such a consideration, and we will not discuss such questions in the sequel. It follows from the definition of $\hat{z}_{n}$ that

[TABLE]

We recall that $z_{n}$ is a vector representative of the class $\hat{z}_{n}$.

Concerning $\hat{y}_{n}$, it can be seen as an estimator of $\hat{x}_{n}$ given the maximum likelihood estimation of $\hat{x}_{0}$. The following proposition establishes consistency of this estimator, we show the geometric contraction in the mean of $(\hat{x}_{n})$ and $(\hat{y}_{n})$. In fact we prove a slightly more general statement that the estimator based on the first $n$ outcomes can be replaced by an estimator based on outcomes in between $l$ and $l+n$. We will prove the almost sure contraction in Proposition 4.4.

Proposition 3.5**.**

Assume (Pur) holds. Then there exist two positive constants $C$ and $\lambda<1$ such that for any probability measure $\nu$ over $({\mathrm{P}}({\mathbb{C}}^{k}),\mathcal{B})$,

[TABLE]

holds for all non-negative integers $l$ and $n$.

In order to prove Proposition 3.5 we study the largest two singular values of $W_{n}$. As is customary in the study of products of random matrices, we make use of exterior products. We recall briefly the relevant definitions: for $p\in{\mathbb{N}}$ and $p$ vectors $x_{1},\ldots,x_{p}$ in ${\mathbb{C}}^{k}$ we denote by $x_{1}\wedge\ldots\wedge x_{p}$ the alternating bilinear form (y_{1},\ldots,y_{p})\mapsto\det\big{(}\langle x_{i},y_{j}\rangle\big{)}_{i,j=1}^{p}. Then, the set of all $x_{1}\wedge\ldots\wedge x_{p}$ is a generating family for the set $\wedge^{p}{\mathbb{C}}^{k}$ of alternating bilinear forms on ${\mathbb{C}}^{k}$, and we can define a hermitian inner product by

[TABLE]

and denote by $\|x_{1}\wedge\ldots\wedge x_{p}\|$ the associated norm. It is immediate to verify that our metric $d$, defined by (7) satisfies

[TABLE]

For $A$ an operator on $\mathbb{C}^{k}$, we write $\wedge^{p}X$ for the operator on $\wedge^{p}{\mathbb{C}}^{k}$ defined by

[TABLE]

Obviously $\wedge^{p}(AB)=\wedge^{p}A\wedge^{p}B$, so that $\|\wedge^{p}(AB)\|\leq\|\wedge^{p}A\|\|\wedge^{p}B\|$. From e.g.  Chapter XVI of [16] or Lemma III.5.3 of [4], we have in addition for $1\leq p\leq k$

[TABLE]

where $a_{1}(A)\geq\ldots\geq a_{k}(A)$ are the singular values of $A$, i.e.  the square roots of eigenvalues of $A^{*}A$, labelled in decreasing order.

Our strategy to prove Proposition 3.5 is to bound the right hand side of Eq. 26 by a submultiplicative function $f:\mathbb{N}\to\mathbb{R}_{+}$ and then use the Fekete’s lemma. We will show that the function

[TABLE]

have these desired properties. The following lemma establishes an exponential decay of this function.

Lemma 3.6**.**

Assume (Pur). Then there exist two positive constants $C$ and $\lambda<1$ such that

[TABLE]

Proof.

First, we prove $\lim_{n\to\infty}f(n)=0$. To prove it, we express the function $f(n)$ using the process $W_{n}$ as

[TABLE]

By definition the eigenvalues of $M_{n}^{\frac{1}{2}}$ are the singular values of $W_{n}/\sqrt{{\operatorname{tr}}(W_{n}^{*}W_{n})}$. Since by Proposition 2.2, $M_{n}$ converges ${\mathbb{P}}^{\mathrm{ch}}\operatorname{-}\mathrm{a.s.}$ to a rank one projection,

[TABLE]

Using Eq. 29 we then conclude that

[TABLE]

Since $\|\wedge^{2}W_{n}\|\leq\|W_{n}\|^{2}\leq{\operatorname{tr}}(W_{n}^{*}W_{n})$, the expression (31) and Lebesgue’s dominated convergence theorem imply $\lim_{n\to\infty}f(n)=0$.

Second, remark that the function $f$ is submultiplicative. Indeed, for $p,q\in\mathbb{N}$ we have

[TABLE]

and the submultiplicativity follows.

By Fekete’s subadditive Lemma, $\frac{\log f(n)}{n}$ converges to $\inf_{n\in{\mathbb{N}}}\frac{\log f(n)}{n}$, which is (strictly) negative (and possibly equal to $-\infty$) since $f(n)\to 0$. Then there exists $0<\lambda<1$ such that $f(n)\leq\lambda^{n}$ for large enough $n$, and the conclusion follows. ∎

We are now in position to prove Proposition 3.5.

Proof of Proposition 3.5.

The Markov property of $(\hat{x}_{n})$ implies that

[TABLE]

Provided inequality (26) is established for $l=0$, the right hand side of the previous equality can be bounded by $C\lambda^{n}$. It is hence sufficient to prove the inequality for $l=0$.

The case $l=0$ follows from Lemma 3.6 if for any $n\in{\mathbb{N}}$ and any probability measure $\nu$,

[TABLE]

To obtain this inequality, note that from the definitions of $\hat{x}_{n}$, $\hat{y}_{n}$ and $\hat{z}_{n}$, we have that

[TABLE]

holds ${\mathbb{P}}_{\nu}$-almost surely. To get the first inequality we used $\|W_{n}z_{n}\|=\|W_{n}\|$, and $\|x_{0}\wedge z_{n}\|=d(\hat{x}_{0},\hat{z}_{n})\leq 1$. In addition, by definition of ${\mathbb{P}}_{\nu}$,

[TABLE]

which is $f(n)$. Therefore (33) holds and Lemma 3.6 yields the proof. ∎

3.3. Convergence in Wasserstein metric

The remainder of Section 3 is directly devoted to the proof of the second part of Theorem 1.1.

Proof of Eq. 6.

We are supposed to prove that

[TABLE]

is exponentially decaying in $n$. The expression in the supremum on the right hand side is unchanged by adding an arbitrary constant to $f$. This freedom allows us to restrict the supremum to functions bounded by $1$, i.e. $\|f\|_{\infty}\leq 1$.

Let $f\in{\rm Lip}_{1}({\mathrm{P}}(\mathbb{C}^{k}))$ be such a function. Our strategy is to approximate $\hat{x}_{mn+r}$ by $\hat{y}_{mp}\circ\theta^{mq+r}$ with $p=\lfloor\frac{n}{2}\rfloor$ and $q=\lceil\frac{n}{2}\rceil$ so that in particular $p+q=n$. Using telescopic estimates and the invariance of $\nu_{\mathrm{inv}}$ we then have

[TABLE]

We bound the terms on the right hand side using Proposition 3.5 for the first two terms and Proposition 3.4 for the last term. To this end let $C$ and $\lambda<1$ be such that bounds in both these propositions hold true. Since $f$ is $1$-Lipschitz continuous we have

[TABLE]

Proposition 3.5 then implies that

[TABLE]

and similarly with $\nu$ replaced by $\nu_{\mathrm{inv}}$. Regarding the last term in the above telescopic estimate we have by Proposition 3.4,

[TABLE]

where we used the constraint $\|f\|_{\infty}\leq 1$ discussed at the beginning of the proof.

Putting these estimates together we get

[TABLE]

and this concludes the proof of Eq. 6 and therefore of Theorem 1.1. ∎

4. Lyapunov exponents

In this section, we study the almost sure stability exponents. We will always assume ($\phi$-Erg) with the additional assumption that the unique minimal $\operatorname{supp}\mu$-invariant subspace $E$ is ${\mathbb{C}}^{k}$.

Remark 4.1**.**

Assuming $E={\mathbb{C}}^{k}$ amounts to saying that $\phi$ has no transient part. Without this assumption, we would have to take into account the almost sure Lyapunov exponent corresponding to the escape from the transient part. See [3] for a precise account of these ideas.

The relevance of this assumption will stem from the following straightforward inequalities: if $\rho$ is any element of $\mathcal{D}_{k}$ then one has

[TABLE]

and if $\rho$ is faithful (i.e. definite positive), then

[TABLE]

In particular, under the assumption that ($\phi$-Erg) holds with $E={\mathbb{C}}^{k}$, thus $\rho_{\mathrm{inv}}>0$ and for any $\rho\in\mathcal{D}_{k}$, we have

[TABLE]

Let us start by proving the following lemma which concerns ergodicity of $\theta$ w.r.t. the measure ${\mathbb{P}}^{\rho_{\mathrm{inv}}}$.

Lemma 4.2**.**

Assume that ($\phi$-Erg) holds. Then the shift $\theta$ on $(\Omega,\mathcal{O})$ is ergodic with respect to the probability measure ${\mathbb{P}}^{\rho_{\mathrm{inv}}}$.

Proof.

Let $A$, $A^{\prime}$ in $\mathcal{O}_{\ell}$. From the definition of ${\mathbb{P}}^{\rho_{\mathrm{inv}}}$, for $j$ large enough, {\mathbb{P}}^{\rho_{\mathrm{inv}}}\big{(}A\cap\theta^{-j}(A^{\prime})\big{)} equals

[TABLE]

and the Perron–Frobenius Theorem 3.3 implies

[TABLE]

for $\mu^{\otimes l}$-almost all $(v_{1},\ldots,v_{l})$ so that

[TABLE]

which proves the ergodicity. ∎

Now we can state our result concerning Lyapunov exponents.

Proposition 4.3**.**

Assume that ($\phi$-Erg) holds with $E={\mathbb{C}}^{k}$, and that (Pur) holds. Assume $\int\|v\|^{2}\log\|v\|^{2}{\rm d}\mu(v)<\infty$. Then there exist numbers

[TABLE]

such that for any probability measure $\nu$ over $({\mathrm{P}}(\mathbb{C}^{k}),\mathcal{B})$:

1. (1)

for any $p\in\{1,\ldots,k\}$,

[TABLE] 2. (2)

$\gamma_{2}-\gamma_{1}<0$* with $\gamma_{2}-\gamma_{1}$ understood as the limit of $\frac{1}{n}\log\frac{\|\wedge^{2}W_{n}\|}{\|W_{n}\|^{2}}$ whenever $\gamma_{1}=-\infty$,* 3. (3)

[TABLE]

Proof.

Let us start by proving (1). Note that $n\mapsto\log\|\wedge^{p}W_{n}\|$ is subadditive by definition. The existence of the ${\mathbb{P}}^{\rho_{\mathrm{inv}}}\operatorname{-}\mathrm{a.s.}$ limits $\lim_{n\to\infty}\frac{1}{n}\log\|\wedge^{p}W_{n}\|$ then follows from ${\mathbb{E}}^{\rho_{\mathrm{inv}}}(\log\|V\|^{2})\leq\int\|v\|^{2}\log\|v\|^{2}d\mu(v)<\infty$, ${\mathbb{P}}^{\rho_{\mathrm{inv}}}\circ\theta^{-1}={\mathbb{P}}^{\rho_{\mathrm{inv}}}$ and a direct application of Kingman’s subadditive ergodic Theorem (see e.g. [19]). The fact that these limits are ${\mathbb{P}}^{\rho_{\mathrm{inv}}}\operatorname{-}\mathrm{a.s.}$ constant comes from the $\theta$-ergodicity of ${\mathbb{P}}^{\rho_{\mathrm{inv}}}$ proved in Lemma 4.2. Since by Eq. 34 any ${\mathbb{P}}^{\rho}$ is absolutely continuous with respect to ${\mathbb{P}}^{\rho_{\mathrm{inv}}}$, Proposition 2.1 and the $\mathcal{O}$-measurability of $\|\wedge^{p}W_{n}\|$ imply the convergence holds ${\mathbb{P}}_{\nu}\operatorname{-}\mathrm{a.s.}$ The numbers $\gamma_{j}$ are uniquely defined, by defining $\sum_{j=1}^{p}\gamma_{j}$ as the ${\mathbb{P}}^{\rho_{\mathrm{inv}}}\operatorname{-}\mathrm{a.s.}$ limit $\lim_{n\to\infty}\frac{1}{n}\log\|\wedge^{p}W_{n}\|$ and imposing the rule that $\gamma_{j+1}=-\infty$ if $\gamma_{j}=-\infty$. This convention and (29) impose that $\gamma_{j+1}\leq\gamma_{j}$ for $j=1,\ldots,k-1$.

Concerning (2), recall the quantity $f(n)$ defined in Eq. 30. Then Eq. 31 and the inequality ${\operatorname{tr}}\,W_{n}^{*}W_{n}\leq k\|W_{n}\|^{2}$ give

[TABLE]

Jensen’s inequality implies

[TABLE]

so that by Lemma 3.6 and Fatou’s lemma, $\log\lambda\geq\gamma_{2}-\gamma_{1}$ with $\lambda\in(0,1)$.

Finally for (3), from Proposition 2.2, we have

[TABLE]

Since ${\mathbb{P}}_{\nu}\operatorname{-}\mathrm{a.s.}$ $|\langle x_{0},z\rangle|>0$ the proposition holds.

∎

From this proposition we deduce the following almost sure convergence rate for the distance between the Markov chain $(\hat{x}_{n})$ and the $(\mathcal{O}_{n})$-adapted process $(\hat{y}_{n})$.

Proposition 4.4**.**

Assume (Pur) holds and ($\phi$-Erg) holds with $E={\mathbb{C}}^{k}$. Then for any probability measure $\nu$ on $({\mathrm{P}}({\mathbb{C}}^{k}),\mathcal{B})$,

[TABLE]

Proof.

Identity (27) and the definition of $\hat{z}_{n}$ imply

[TABLE]

Proposition 4.3 then yields the result. ∎

Acknowledgments

T.B. and C.P. would like to thank Y. Guivarc’h for his useful comments at an early stage of this work. Y.P. and C.P. would like to thank P. Bougerol for enlightening discussions about random products of matrices. Y.P. and C.P. would like to thank L. Miclo for relevant discussions regarding Markov chains. The research of T.B., Y.P. and C.P. has been supported by the ANR project StoQ ANR-14-CE25-0003-01 and CNRS InFIniTi project MISTEQ.

Appendix A Equivalence of (Pur) and contractivity

We assume $\operatorname{supp}\mu\subset GL_{k}({\mathbb{C}})$. Recall that $T_{\mu}$ is the smallest closed sub-semigroup of $GL_{k}({\mathbb{C}})$ that contains $\operatorname{supp}\mu$. It is said to be contracting if there exists a sequence $(a_{n})_{n\in\mathbb{N}}\subset T_{\mu}$ such that $\lim_{n\to\infty}a_{n}/\|a_{n}\|$ exists and is a rank one matrix.

Proposition A.1**.**

Assume $\operatorname{supp}\mu\subset GL_{k}(\mathbb{C})$ and $T_{\mu}$ is strongly irreducible. Then $\mu$ verifies (Pur) iff. $T_{\mu}$ is contracting.

Proof.

By Proposition 2.2 the implication (Pur) $\implies$ contractivity follows by taking $(a_{n})$ a convergent subsequence of $(W_{n}(\omega))$ for $\omega\in\operatorname{supp}{\mathbb{P}}^{\mathrm{ch}}$.

We prove the opposite implication by contradiction. Following [12, Lemma 3], under the assumptions of the proposition, $T_{\mu}$ is contracting iff. for any two $\hat{x},\hat{y}\in{\mathrm{P}}({\mathbb{C}}^{k})$ there exists a sequence of matrices $(a_{n})\subset T_{\mu}$ such that

[TABLE]

Now, assume that contractivity holds but (Pur) does not. Namely, that $T_{\mu}$ is contracting but there exists an orthogonal projector $\pi$ of rank $\geq 2$, such that for any $a\in T_{\mu}$,

[TABLE]

Let $x,y\in\operatorname{range}\pi$ be orthonormal vectors. Then $\langle ax,ay\rangle=\langle x,y\rangle=0$, and $\|ax\|,\|ay\|$ are nonzero, so that $d(a\cdot\hat{x},a\cdot\hat{y})=1$. As this holds for any $a\in T_{\mu}$, contractivity cannot hold. This contradiction yields the proposition. ∎

Appendix B Set of invariant measures under assumption (Pur)

A quantum channel is a map $\phi$ on $M_{k}({\mathbb{C}})$ of a form

[TABLE]

where $\mu$ is a measure satisfying the normalization condition (1). The decomposition of quantum channels to irreducible components was derived in [20, 5, 2]. The space $\mathbb{C}^{k}$ is decomposed into orthogonal subspaces, one subspace is transient and in all other subspaces the map has a canonical tensor product structure. We recall these results.

There exists a decomposition

[TABLE]

with the following properties. We denote by $v_{j}$ the restriction of $v$ to $\mathbb{C}^{n_{j}}$.

1. (e1)

All invariant states are supported in the subspace $L=\mathbb{C}^{n_{1}}\oplus\dots\oplus\mathbb{C}^{n_{d}}\oplus 0$, 2. (e2)

The restriction of $v$ to this subspace is block diagonal,

[TABLE] 3. (e3)

For each $j=1,\dots,d$ there is a decomposition $\mathbb{C}^{n_{j}}=\mathbb{C}^{k_{j}}\otimes\mathbb{C}^{m_{j}},\,n_{j}=k_{j}m_{j}$, a unitary matrix $U_{j}$ on $\mathbb{C}^{n_{j}}$ and a matrix $\tilde{v}_{j}$ on $\mathbb{C}^{k_{j}}$ such that

[TABLE] 4. (e4)

There exists a full rank positive matrix $\rho_{j}$ on $\mathbb{C}^{k_{j}}$ such that

[TABLE]

is a fixed point of $\phi$.

It follows from (e3) and (e4) that the set of fixed points for $\phi$ is

[TABLE]

The decomposition simplifies under the purification assumption.

Proposition B.1**.**

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

[TABLE]

Proof.

The statement follows from the discussion preceding the proposition if we show that (Pur) implies $m_{1}=\dots=m_{d}=1$. Assume that one of the $m_{j}$, e.g. $m_{1}$, is greater than $1$. Let $x$ be a norm one vector in $\mathbb{C}^{k_{1}}$. Then $\pi=U_{j}\pi_{\hat{x}}\otimes{\rm Id}_{\mathbb{C}^{m_{1}}}U_{j}^{*}\oplus 0\oplus\dots\oplus 0$ is a projection with rank bigger than $1$, and by Eq. 38 we have

[TABLE]

for $\mu$-almost all $v$. This contradicts (Pur). ∎

It is clear from Eq. 37 that to each extremal fixed point $0\oplus\dots\oplus\rho_{j}\oplus\dots\oplus 0$ corresponds a unique invariant measure $\nu_{j}$ supported on its range $F_{j}$. The converse is the subject of the next proposition.

Proposition B.2**.**

Assume (Pur) holds. Then any $\Pi$-invariant probability measure is a convex combination of the measures $\nu_{j}$, $j=1,\ldots,d$.

Proof.

Let $\nu$ be a $\Pi$-invariant probability measure. Let $f$ be a continuous function. From Lemma 2.3,

[TABLE]

Proposition 2.1 implies

[TABLE]

with $\rho_{\nu}\in\mathcal{D}_{k}$ a fixed point of $\phi$. By Proposition B.1, (Pur) implies that there exist non negative numbers $t_{1},\ldots,t_{d}$ summing up to one such that $\rho_{\nu}=t_{1}\rho_{1}\oplus\cdots\oplus t_{d}\rho_{d}\oplus 0_{M_{D}({\mathbb{C}})}$. From the definition of ${\mathbb{P}}^{\rho_{\nu}}$,

[TABLE]

where we used the abuse of notation $\rho_{j}\equiv 0\oplus\cdots\oplus\rho_{j}\oplus\cdots\oplus 0$. Using Proposition 2.1, it follows that,

[TABLE]

Then Lemma 2.3 and the $\Pi$-invariance of each measure $\nu_{j}$ yield the proposition. ∎

Appendix C Products of special unitary matrices

Proposition C.1**.**

Assume $\operatorname{supp}\mu\subset SU(k)$. Let $G$ be the smallest closed subgroup of $SU(k)$ such that $\operatorname{supp}\mu\subset G$. For any $\hat{x}\in{\mathrm{P}}({\mathbb{C}}^{k})$, let $[\hat{x}]_{G}$ be the orbit of $\hat{x}$ with respect to $G$ and the action $G\times{\mathrm{P}}({\mathbb{C}}^{k})\ni(v,\hat{x})\mapsto v\cdot\hat{x}$. Namely, $[\hat{x}]_{G}:=\{\hat{y}\in{\mathrm{P}}({\mathbb{C}}^{k})\ |\ \exists v\in G\mbox{ s.t. }\hat{y}=v\cdot\hat{x}\}$. Then, for any $\hat{x}$, there exists a unique $\Pi$-invariant probability measure supported on $[\hat{x}]_{G}$, and this unique invariant measure is uniform in the sense that for any $v\in G$ it is invariant by the map $\hat{x}\mapsto v\cdot\hat{x}$.

Corollary C.2**.**

With the assumption and definitions of the last proposition, if $G=SU(k)$, $\Pi$ has a unique invariant probability measure and this probability is the uniform one on ${\mathrm{P}}({\mathbb{C}}^{k})$.

Proof.

The Corollary being a trivial consequence of $G=SU(k)\implies[\hat{x}]_{G}={\mathrm{P}}({\mathbb{C}}^{k})\ \forall\hat{x}\in{\mathrm{P}}({\mathbb{C}}^{k})$, we are left with proving the Proposition.

Let $P_{\mu}$ be the Markov kernel on $G$ defined by the left multiplication: $P_{\mu}f(v)=\int_{G}f(uv)d\mu(u)$. Since $G$ is compact as a closed subset of $SU(k)$, following [1, Proposition 4.8.1, Theorem 4.8.2], the unique $P_{\mu}$-invariant probability measure $\mu_{G}$ on $G$ is the normalized Haar measure on $G$. Since $G$ is compact, Prokhorov’s Theorem implies that for any $u\in G$,

[TABLE]

Let $\hat{x}\in{\mathrm{P}}({\mathbb{C}}^{k})$. Since $\operatorname{supp}\mu\subset G$, for any $\hat{y}\in[\hat{x}]_{G}$, $\Pi(\hat{y},[\hat{x}]_{G})=1$. Then, $[\hat{x}]_{G}$ being compact, there exists a $\Pi$-invariant measure $\nu$ supported on $[\hat{x}]_{G}$.

Let $f$ be a continuous function on $[\hat{x}]_{G}$. Then,

[TABLE]

For each $\hat{y}\in[\hat{x}]_{G}$ let $u_{y}\in G$ be s.t. $\hat{y}=u_{y}\cdot\hat{x}$. The map $v\mapsto vu_{y}\cdot\hat{x}$ being continuous, setting $u=u_{y}$, the weak convergence (40) and Lebesgue’s dominated convergence theorem imply,

[TABLE]

It follows that $\nu$ is the image measure of $\mu_{G}$ by the application $v\mapsto v\cdot\hat{x}$. The left multiplication invariance of the Haar measure $\mu_{G}$ yields the invariance of $\nu$ by the map $\hat{x}\mapsto v\cdot\hat{x}$ for any $v\in G$. ∎

Example C.3**.**

Let $\mu=\frac{1}{2}(\delta_{v_{1}}+\delta_{v_{2}})$ with,

[TABLE]

Then $G=SU(2)$ and the uniform measure on ${\mathrm{P}}({\mathbb{C}}^{2})$ is the unique $\Pi$-invariant probability measure.

Proof.

Following Proposition C.1, it is sufficient to prove that any element of $SU(2)$ is the limit of a sequence of products of $v_{1}$ and $v_{2}$.

Let $\sigma_{1},\sigma_{2},\sigma_{3}$ be the usual Pauli matrices:

[TABLE]

The Pauli matrices being generators of $SU(2)$ in its fundamental representation, for any $u\in SU(2)$, there exist three reals $\theta_{1},\theta_{2},\theta_{3}\in{\mathbb{R}}$ s.t.,

[TABLE]

Especially, $v_{1}=\exp(i\sigma_{3})$ and $v_{2}=\exp(i\sigma_{1})$. Since for any $i=1,2,3$, $\exp(i\theta_{i}\sigma_{i})=\exp(i(\theta_{i}+2\pi)\sigma_{i})$, taking limits of sequences of powers of $v_{1}$ or $v_{2}$, for any $\theta\in{\mathbb{R}}$, both

[TABLE]

are elements of $G$. It remains to show that any $u\in SU(k)$ is a product of elements equal to $\exp(i\theta\sigma_{1})$ or $\exp(i\theta\sigma_{3})$ with $\theta$ real.

Fix $(\theta_{1},\theta_{2},\theta_{3})\in{\mathbb{R}}^{3}$. Then using spherical coordinates in ${\mathbb{R}}^{3}$, there exist $r\in{\mathbb{R}}_{+}$, $\theta\in[0,\pi]$ and $\varphi\in[0,2\pi[$ such that $\theta_{1}=r\cos\theta$, $\theta_{2}=r\sin\theta\cos\varphi$ and $\theta_{3}=r\sin\theta\sin\varphi$. Then by direct computation,

[TABLE]

It follows that as a product of elements of $G$, $e^{i(\theta_{1}\sigma_{1}+\theta_{2}\sigma_{2}+\theta_{3}\sigma_{3})}\in G$, hence $G=SU(2)$ and the example holds. ∎

Example C.4**.**

Let $\mu=\frac{1}{2}(\delta_{v_{1}}+\delta_{v_{2}})$ with,

[TABLE]

Then $G=\{\pm{\rm Id}_{{\mathbb{C}}^{2}},\pm v_{1},\pm v_{2},\pm v_{1}v_{2}\}$. For $z\in{\mathbb{C}}$, let $e_{z}=(1,z)^{\mathsf{T}}$ and $e_{\infty}=(0,1)^{\mathsf{T}}$. With the conventions $\infty^{-1}=0$, $0^{-1}=\infty$ and $-\infty=\infty$, for any $z\in{\mathbb{C}}\cup\{\infty\}$, $[\hat{e}_{z}]_{G}=\{\hat{e}_{z},\hat{e}_{z^{-1}},\hat{e}_{-z},\hat{e}_{-z^{-1}}\}$ and the measure $\frac{1}{4}(\delta_{\hat{e}_{z}}+\delta_{\hat{e}_{-z}}+\delta_{\hat{e}_{z^{-1}}}+\delta_{\hat{e}_{-z^{-1}}})$ is a $\Pi$-invariant probability measure.

The proof of this example is obtained by explicit computation.