Quadratic Open Quantum Harmonic Oscillator

TL;DR

This paper analyzes the dynamics of a quadratic open quantum harmonic oscillator, demonstrating exponential convergence to a unique equilibrium state and providing explicit convergence rates, with implications for two-photon processes.

Contribution

It introduces a detailed analysis of the quantum open system evolution using a specific Lindblad generator related to the $sl_2$ Lie algebra, including explicit convergence rates.

Findings

Initial density matrices evolve to a fully supported state

Convergence to equilibrium is exponentially fast

Explicit convergence rates are computed for various parameters

Abstract

We study the quantum open system evolution described by a Gorini-Kossakowski-Sudarshan-Lindblad generator with creation and annihilation operators arising in Fock representations of the $sl_2$ Lie algebra. We show that any initial density matrix evolves to a fully supported density matrix and converges towards a unique equilibrium state. We show that the convergence is exponentially fast and we exactly compute the rate for a wide range of parameters. We also discuss the connection with the two-photon absorption and emission process.

Equations

$$\mathcal{L}(x)$$

$$\hbox{\rm Dom}(B)=\hbox{\rm Dom}(B^{+})=\hbox{\rm Dom}(M)=\left\{u=\sum_{n\geq 0}u_{n}e_{n}\,\Big{|}\,\sum_{n\geq 0}n^{2}|u_{n}|^{2}<\infty\right\}$$

$$B{e}_n$$

$$B^{+}{e}_n$$

$$M{e}_n$$

$${\omega }_n=n(n+r-1)$$

$$B\otimes A^{+}(g)+B^{+}\otimes A(g).$$

$$G=-\frac{{\lambda }^2}{2}B{B}^{+}-\frac{{\mu }^2}{2}{B}^{+}B-\mathrm{i}\left({\zeta }^{+}B{B}^{+}+{\zeta }^{-}{B}^{+}B\right)$$

$$L_1=\mu B,L_2=\lambda B^{+}.$$

$$G=-\left(\frac{{\lambda }^2}{2}+\mathrm{i}{\zeta }^{+}\right){\omega }_{N+1}-\left(\frac{{\mu }^2}{2}+\mathrm{i}{\zeta }^{-}\right){\omega }_N$$

$$P_t{e}_n={\mathrm{e}}^{-t\left(\left(\frac{{\lambda }^2}{2}+\mathrm{i}{\zeta }^{+}\right){\omega }_{n+1}+\left(\frac{{\mu }^2}{2}+\mathrm{i}{\zeta }^{-}\right){\omega }_n\right)}{e}_n.$$

$$\mathcal{L}\text{raise}1.0pt-(x)[u,v]=⟨Gu,xv⟩+\sum _{\ell =1}^2⟨L_{\ell }u,x{L}_{\ell }v⟩+⟨u,xGv⟩,$$

$$⟨v,{\mathcal{T}}_t(x)u⟩=⟨v,P_t^{\ast }x{P}_{t}u⟩+\sum _{\ell =1}^2{\int }_0^t⟨L_{\ell }{P}_{t-s}v,{\mathcal{T}}_s(x)L_{\ell }{P}_{t-s}u⟩ds,$$

$$\mathcal{L}\text{raise}1.0pt-(f(N))$$

$${\lambda }^2(n+1)(n+r)\left(f(n+1)-f(n)\right)+{\mu }^{2}n(n+r-1)\left(f(n-1)-f(n)\right)$$

$$\left|{\lambda }^2(2r+3)+{\mu }^2(2r-3)\right|,\frac{\left|{\lambda }^2(3r+1)+{\mu }^2(r-1)\right|}{2},{\lambda }^2|r|,$$

$$\mathcal{L}\text{raise}1.0pt-((N+1\mathrm{l}{)}^2)\le -2({\mu }^2-{\lambda }^2)N^3+b(N+1\mathrm{l}{)}^2.$$

$$\mathcal{L}(x)$$

$$\rho =(1-{\nu }^2)\sum _{n\ge 0}{\nu }^{2n}|e_n⟩⟨e_n|$$

$$\sum _{k=0}^n{\lambda }^2{\nu }^{2k}{\omega }_{k+1}\left(|e_{k+1}⟩⟨e_{k+1}|-|e_k⟩⟨e_k|\right)+\sum _{k=1}^n{\mu }^2{\nu }^{2k}{\omega }_k\left(|e_{k-1}⟩⟨e_{k-1}|-|e_k⟩⟨e_k|\right)$$

$$\underset{n\to \infty }{\lim }{\Vert {\mathcal{L}}_{\ast }({\rho }_n)\Vert }_1=\underset{n\to \infty }{\lim }{\lambda }^2{\nu }^{2n}({\omega }_{n+1}+{\omega }_n)=0$$

$$\mathcal{N}(\mathcal{T})=\left\{x\in \mathcal{B}(\mathsf{h})\mid {\mathcal{T}}_t(x^{\ast }x)={\mathcal{T}}_t(x^{\ast }){\mathcal{T}}_t(x),{\mathcal{T}}_t(x{x}^{\ast })={\mathcal{T}}_t(x){\mathcal{T}}_t(x^{\ast }),\forall t\ge 0\right\}.$$

$${\mathcal{T}}_t(x^{\ast }x)=x^{\ast }x={\mathcal{T}}_t(x^{\ast }){\mathcal{T}}_t(x),$$

$$\mathcal{D}\left(\mathcal{T}\right):=\left\{{\mathrm{e}}^{-\mathrm{i}tH}{L}_{\ell }{\mathrm{e}}^{\mathrm{i}tH},{\mathrm{e}}^{-\mathrm{i}tH}{L}_{\ell }^{\ast }{\mathrm{e}}^{\mathrm{i}tH}\mid \ell \ge 1,t\ge 0\right\}.$$

$$XB\subseteq BX,X{B}^{+}\subseteq B^{+}X$$

$$X{B}^{+}B\subseteq B^{+}XB\subseteq B^{+}BX,XB{B}^{+}\subseteq BX{B}^{+}\subseteq B{B}^{+}X,$$

$$X(s+N)\subseteq (s+N)X$$

$$f(N)B^{+}{e}_n-B^{+}f(N)e_n={\omega }_{n+1}^{1/2}(f(n+1)-f(n))e_n$$

$$\underset{t\to \infty }{\lim }{\mathcal{T}}_{\ast t}(\eta )=\rho$$

$$X=\left(\sum _{k\ge n}(k+1{)}^{-1}(k+r{)}^{-1}\right)|e_n⟩⟨e_n|$$

$$\mathcal{L}\text{raise}1.0pt-(X)=-|e_0⟩⟨e_0|\le 0.$$

Full text

Ameur Dhahri111Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, I-20133 Milano (Italy), ameur.dhahri@polimi.it, Franco Fagnola222Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, I-20133 Milano (Italy) franco.fagnola@polimi.it and Hyun Jae Yoo333Department of Applied Mathematics, Hankyong National University, 327 Jungang-ro, Anseong-si, Gyeonggi-do 456-749, Korea yoohj@hknu.ac.kr

Keywords: quantum harmonic oscillator, quantum Markov semigroup, Fock representations of the $\mathfrak{sl}_{2}$ algebra, spectral gap.

Subject Classification: 46L55, 82C10, 60J27.

1 Introduction

Models of quantum harmonic oscillators are usually based on commutation relations. The Heisenberg-Weyl algebra commutation relations $[H,A]=-A$, $[H,A^{+}]=A^{+}$, $[A,A^{+}]=\mathrm{1}$, or in terms of position $Q=(A^{+}+A)/\sqrt{2}$ and momentum $P=\mathrm{i}(A^{+}-A)/\sqrt{2}$, $[H,Q]=-\mathrm{i}P$, $[H,P]=\mathrm{i}Q$, $[Q,P]=\mathrm{i}$ are the foundation at the best known one. This model arises, for instance, replacing time derivatives in the classical equation $q^{\prime\prime}=-q$ by commutators with the Hamiltonian operator $H=(P^{2}+Q^{2})/2$ so that we can write it as $[H,[H,Q]]=Q$. If we fix $H$ and define $P=\mathrm{i}[H,Q]$, then the double commutator equation reads as $[H,P]=\mathrm{i}Q$ and if, moreover, we want $H,P,Q$ to be elements of a Lie algebra, the Jacobi identity $[H,[P,Q]]+[P,[Q,H]]+[Q,[H,P]]=0$ implies that $[P,Q]$ commutes with $H$. The most natural choice as $-\mathrm{i}\mathbb{1}$ corresponds to the commutation relations of the Heisenberg-Weyl algebra.

Other choices lead to different models of quantum oscillators (see, for instance, [6] and the references therein) and for some of them it is possible to develop a complete theory describing explicitly spectra of observables, eigenvectors, time evolution, etc. The choice $[P,Q]=-2\mathrm{i}H$ corresponds to the commutation relations of the $\mathfrak{sl}_{2}$ Lie algebra.

This is a three dimensional simple ∗-Lie algebra with basis $\{B^{+},B,M\}$, commutation relations $[B,B^{+}]=M$, $[M,B^{+}]=2B^{+}$, $[M,B]=-2B$ and involution $B^{*}=B^{+}$, $M^{*}=M$. The construction of Fock representations of $\mathfrak{sl}_{2}$ Lie algebra and of the current algebra associated to its central extension motivated a large number of papers extending it in different directions: see Ref.[26] for the case of free white noise; [4] for the connection with quantum Lévy processes; Refs.[1, 2, 10, 11] for the construction of the quadratic Fock functor.

The weak coupling limit (see [5]) of an harmonic oscillator coupled with a reservoir in equilibrium with inverse temperature $\beta>0$ gives rise to a fundamental model of an open quantum system with a lot of deep properties and quantities that can be computed explicitly called in the literature the open quantum harmonic oscillator (see e.g. Ref.[21] and the references therein). If we consider, instead, the formal Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) generator arising in the weak coupling limit of an oscillator based on the Fock representations of $\mathfrak{sl}_{2}$ commutation relations we find

[TABLE]

where $\zeta^{\pm},\lambda,\mu$ are real parameters and $\lambda,\mu>0$. This is called the quadratic open quantum harmonic oscillator because the operators $B,B^{+}$ are the annihilation and creation operators arising in Fock representations of $\mathfrak{sl}_{2}$ and the action of $BB^{+}$ and $B^{+}B$ (see formulae (2)) is quadratic with respect to the level of the system while, for the open quantum harmonic oscillator, it is linear. Constants $\lambda^{2},\mu^{2}$ are related with the inverse temperature $\beta$ by $\lambda^{2}\mu^{-2}=\hbox{\rm e}^{-s\beta}$ for some $s>0$.

This is a simple and natural model, however, contrary to what happens for the open quantum harmonic oscillator, it does not admit explicit solutions except for the formula of the invariant state. As an example, if one looks at the action on the abelian algebra of functions of the number operator, one finds a birth-and-death process with quadratic jump rates for which explicit representations for transition probabilities, to the best of our knowledge (see Ref.[25]), are not known.

In this paper, we first show that the formal GKSL generator with unbounded operators $B,B^{+}$ generates a unique quantum Markov semigroup and we establish the existence of a unique explicit equilibrium state. Then, we study the behaviour of the evolution of states and observables for all values of parameters involved. We prove that any initial state converges towards the unique equilibrium state for the trace norm (Theorem 1). We also prove (Theorem 3) that any initial state $\rho_{0}$, in particular also a pure state, evolves to a faithful state $\rho_{t}$ for all $t>0$. Moreover, we show that, for some special values of a parameter $r$ determining the Fock representation of the $\mathfrak{sl}_{2}$ commutation relations this model is intimately related with the two-photon absorption and emission process studied in [8, 16]. Finally, we show that convergence towards the unique invariant state is exponentially fast (with respect to the Hilbert-Schmidt norm induced by the invariant state) and we also compute the sharp exponential rate for a wide range of parameters (Theorem 6). Our analysis shows, in particular, that the decay rate of off-diagonal terms of density matrices is smaller than the rate of convergence of the diagonal part towards the unique equilibrium state for more and more values of the parameter $r$ as the inverse temperature $\beta$ becomes big, i.e. the reservoir becomes cooler. In other words, at low temperatures, decoherence is slower than relaxation for $r$ away from [math].

The paper is organized as follows. In Section 2, we introduce the model of the quadratic open quantum harmonic oscillator. The full characterization of invariant states and the asymptotic behaviour of the associated quantum Markov semigroup are studied in Section 3. The close relationship with two-photon absorption and emission process is studied in Section 4. In Section 5, we show that for all initial state the support of the state evolved at any time $t>0$ is full. The rate of the exponentially fast convergence towards the unique invariant state is studied in Section 6.

2 The model

Let $\mathsf{h}$ be the Hilbert space $\mathsf{h}=\ell^{2}({\mathbb{N}})\simeq\Gamma(\mathbb{C})$ with canonical orthonormal basis $(e_{n})_{n\geq 0}$. We consider the operators $B,B^{+},M$ of the Fock representation of the renormalized square of the white noise Lie algebra $B$, $B^{+}$ and $M$ with domain

[TABLE]

defined, on vectors of the canonical orthonormal basis, by

[TABLE]

where $r>0$ is a real parameter (see Ref.[4] Section 3.2 p. 134 for the explanation why $r$ must be non-negative) and

[TABLE]

which is strictly positive for all $n\geq 1$ and satisfies $\omega_{0}=0$.

Note that the domain of $B,B^{+}$ and $M$ coincides with the domain of the number operator $N$ defined by $Ne_{n}=ne_{n}$ for all $n\geq 0$.

We consider the formal Lindblad generator (1) which is of weak coupling limit type (see Refs.[3, 5]) since it arises in the weak coupling limit of a system with Hilbert space $\ell^{2}(\mathbb{N})$ and Hamiltonian $H_{S}$ given by the number operator $N$ coupled to a Boson reservoir in equilibrium with inverse temperature $\beta>0$ and interaction operator

[TABLE]

Constants $\lambda^{2},\mu^{2}$ satisfy $\lambda^{2}\mu^{-2}=\hbox{\rm e}^{-s\beta}$ for some $s>0$.

Moreover (see Section 4) for $r=1/2$ (resp. $r=3/2$) and a suitable choice of the real constants $\zeta^{-},\zeta^{+}$ we find the even (resp. odd) part of the two-photon absorption and emission generator studied in Refs.[8, 16].

Let $G$ be the operator defined on the domain ${\rm Dom}(N)$ of the number operator by

[TABLE]

and let $L_{1},L_{2}$ be the operators defined on ${\rm Dom}(N)$ by

[TABLE]

Clearly $G$ is a function of the number operator $N$, defined on the same domain of $B,B^{+}$ and $M$, since

[TABLE]

with negative real part hence generates a strongly continuous semigroup of contractions $(P_{t})_{t\geq 0}$ on $\mathsf{h}$ explicitly given by

[TABLE]

For every $x\in{\mathcal{B}}(\mathsf{h})$ the formal generator is the sesquilinear form

[TABLE]

for $u,v\in{\rm Dom}(G)={\rm Dom}(N^{2})$. One can easily check that conditions for constructing the minimal quantum dynamical semigroup (QDS) associated with the above $G,L_{1},L_{2}$ ((H-min) in Ref.[12]) hold and this semigroup ${\mathcal{T}}=({\mathcal{T}}_{t})_{t\geq 0}$ satisfies the so-called Lindblad equation

[TABLE]

for all $u,v\in{\rm Dom}(G)$.

A straightforward computation using the CCR (it could be done considering quadratic forms on the linear manifold generated by vectors $(e_{n})_{n\geq 0}$ if one wants to cope with unboundedness of the involved operators but we prefer to simplify the notation) shows that

[TABLE]

Taking $f(n)=(n+1)^{2}$, for $n\geq 1-r$ we easily find

[TABLE]

and, for $0\leq n<1-r$, i.e. $n=0$ we obviously find [math]. Therefore, defining as $b$ the maximum of the three constants

[TABLE]

we have

[TABLE]

As a consequence, if $\lambda\leq\mu$, {\mathcal{L}}\kern-9.0pt\raise 1.0pt\hbox{-} satisfies a well known criterion for conservativity (Ref.[12] Theorem 3.40). Moreover, for $\lambda>\mu$ the formal generator satisfies a simple criterion for nonconservativity, see Ref.[19], Example 2. Then the minimal QDS is Markov (or conservative) if and only if $\lambda\leq\mu$. It follows from conservativity that the minimal QDS is the unique solution of equation (4). Moreover an operator $x\in{\mathcal{B}}(\mathsf{h})$ belongs to the domain of the generator ${\mathcal{L}}$ if and only if the sesquilinear form {\mathcal{L}}\kern-9.0pt\raise 1.0pt\hbox{-}(x) is bounded (see Ref.[12] Prop. 3.33 p.64).

The action of ${\mathcal{L}}$ on the linear manifold ${\cal M}={\rm span}\{|e_{j}\rangle\langle e_{k}|:j,k\geq 0\}$ of finite range operators is given by

[TABLE]

3 Invariant states and asymptotic behaviour

The behaviour of the quadratic open quantum harmonic oscillator and the structure of its invariant states depends crucially upon the parameters $\lambda$ and $\mu$. We begin by considering the case where $\mu>\lambda>0$.

Proposition 1

If $\nu=\lambda/\mu<1$ then the normal state

[TABLE]

is invariant.

Proof. Let ${\mathcal{L}}_{*}$ be the generator of the predual semigroup ${\mathcal{T}}_{*}=({\mathcal{T}}_{*t})_{t\geq 0}$, acting on the Banach space of trace class operators on $\mathsf{h}$. Consider the approximations ${\rho}_{n}=(1-\nu^{2})\sum_{k=0}^{2n}\nu^{2k}|e_{k}\rangle\langle e_{k}|$, of $\rho$ by finite rank operators.

The operators ${\rho}_{n}$ belong to the domain of ${\mathcal{L}}_{*}$ and we can write ${\mathcal{L}}_{*}(\rho_{n})$ as $(1-\nu^{2})$ times

[TABLE]

Terms in the above summations vanish because $\mu^{2}\nu^{2(k+1)}=\lambda^{2}\nu^{2k}$ for all $k\geq 0$. Moreover

[TABLE]

because $\lambda<\mu$. Since the operator $\mathcal{L}_{*}$ is closed, it follows that $\rho$ belongs to the domain of $\mathcal{L}_{*}$ and $\mathcal{L}_{*}(\rho)=0$.

In order to show uniqueness of the invariant state (6) we begin by recalling that the support projection $p$ of an invariant state with density matrix $\rho$, i.e. the orthogonal projection onto the range of $\rho$, satisfies $\mathcal{T}_{t}(p)\geq p$ for all $t\geq 0$ (see e.g. [14] Theorem II.1). Such projections, called subharmonic, are easily characterized in terms of invariant subspaces of operators $P_{t}$ and $L_{1},L_{2}$ considered in Section 2. A QMS is called irreducible if the only subharmonic projections are the trivial ones $0,\hbox{\rm 1\kern-2.8ptl}$. In this case, it is well-known (see Ref.[18] Lemma 1) that a faithful invariant state, if it exists, is unique because the set of fixed points for the QMS $\mathcal{T}$ is the trivial algebra $\mathbb{C}\hbox{\rm 1\kern-2.8ptl}$. In our framework we can prove the following.

Proposition 2

The QMS $\mathcal{T}$ is irreducible for all $\lambda\leq\mu$. In particular, if $\lambda<\mu$, the state (6) is the unique $\mathcal{T}$-invariant state.

Proof. The range of any non-trivial subharmonic projection determines an invariant subspace for the operators $P_{t}$ for all $t>0$ (see Ref.[14] Theorem III.1). Since these operators are normal and compact, these invariant subspaces are generated by eigenvectors of $P_{t}$. Moreover, knowing the spectral decomposition of $P_{t}$ (it is a function of the number operator!) we infer that they are generated by collections of vectors $(e_{n})_{n\in I}$ for some subset $I$ of $\mathbb{N}$. Invariance of these subspaces for $B$ and $B^{+}$ implies then that they must coincide with the whole of $\mathsf{h}$. This proves that the QMS is irreducible.

If $\lambda<\mu$ the QMS admits the faithful invariant state (6) and so the set of fixed points for the QMS $\mathcal{T}$ is the trivial algebra $\mathbb{C}\hbox{\rm 1\kern-2.8ptl}$. It follows then from Lemma 1 of Ref.[18] that (6) is the unique invariant state. $\square$

Applying the main result of Ref.[9] we can also show convergence towards the invariant state in trace norm. As a preliminary step we prove the following result which is interesting on its own

Proposition 3

If $\nu=\lambda/\mu<1$ the decoherence free subalgebra

[TABLE]

and the fixed point algebra $\mathcal{F}(\mathcal{T})=\left\{x\in\mathcal{B}(\mathsf{h})\mid\mathcal{T}_{t}(x)=x\ \forall t\geq 0\right\}$ are trivial.

Proof. It is well-known that $\mathcal{N}(\mathcal{T})$ is a von Neumann subalgebra of $\mathcal{B}(\mathsf{h})$ (see e.g. Proposition 2.1 (3) of Ref.[9]). Moreover, since the invariant state $\rho$ defined in (6) is faithful, also $\mathcal{F}(\mathcal{T})$ is a von Neumann subalgebra of $\mathcal{B}(\mathsf{h})$. Indeed, if $x$ belongs to ${\mathcal{F}}(\mathcal{T})$, then, by $2$-positivity, $\mathcal{T}_{t}(x^{*}x)\geq\mathcal{T}_{t}(x^{*})\mathcal{T}_{t}(x)=x^{*}x$ and $\hbox{\rm{tr}}(\rho(\mathcal{T}_{t}(x^{*}x)-x^{*}x))=0$ because $\rho$ is invariant. It follows that $\mathcal{T}_{t}(x^{*}x)=x^{*}x$ i.e. $x^{*}x\in\mathcal{F}(\mathcal{T})$.

As a by product, if $x\in\mathcal{F}(\mathcal{T})$, then

[TABLE]

and the same identity holds exchanging $x$ and $x^{*}$, i.e. $\mathcal{F}(\mathcal{T})$ is contained in $\mathcal{N}(\mathcal{T})$.

Thus, it suffices to prove that $\mathcal{N}(\mathcal{T})$ is trivial. To this end, we apply Theorem 4.1 of Ref.[9] characterizing $\mathcal{N}(\mathcal{T})$ as the generalized commutator of the set of unbounded operators

[TABLE]

where $H=\zeta^{+}BB^{+}+\zeta^{-}B^{+}B$. The additional technical domain assumptions that can be easily checked taking as $D$ the linear manifold spanned by finite linear combinations of vectors $e_{n}$ of the orthonormal basis and as operator $C$ the number operator or $(N+\mathbb{1})^{2}$.

If $X$ is an operator in the generalized commutator of (7), then it is, by definition of generalized commutator, bounded and, in particular, it satisfies

[TABLE]

(meaning that $BX$ is an ampliation of $XB$ and $B^{+}X$ is an ampliation of $XB^{+}$). It follows that

[TABLE]

and so, since the difference $BB^{+}-B^{+}B$ is $2N+r\hbox{\rm 1\kern-2.8ptl}$, $NX$ is an ampliation of $XN$ and

[TABLE]

for all $s>0$. Left and right multiplying by the resolvent $(s+N)^{-1}$, since the operators $(s+N)^{-1}X$ and $X(s+N)^{-1}$ are bounded, we find $(s+N)^{-1}X=X(s+N)^{-1}$ for all $s>0$. This shows that $X$ commutes with every spectral projection $|e_{n}\rangle\langle e_{n}|$ of the number operator and so it is a function $f(N)$ of the number operator itself. However, since $\omega_{n+1}>0$ for all $n\geq 0$ if $r>0$

[TABLE]

vanishes if and only if $f$ is constant and so the generalized commutator of (7) is trivial. $\square$

We are now in a position to prove the following.

Theorem 1

If $0<\nu=\lambda/\mu<1$ then (6) is the unique invariant state and

[TABLE]

in trace norm for all initial state $\eta$.

Proof. Since $\mathcal{F}(\mathcal{T})=\mathcal{N}(\mathcal{T})=\mathbb{C}\mathbb{1}$ by Proposition 3, the conclusion is immediate from Theorem 3.3 of Ref.[9]. $\square$

We complete the study of the asymptotic behaviour by the following.

Proposition 4

If $\lambda=\mu>0$ the QMS is transient. In particular, it has no invariant state.

Proof. By Theorem 5 Ref. [15], it suffices to find a strictly positive bounded operator $X$ such that $\mathcal{T}_{t}(X)\leq X$ for all $t\geq 0$. Inspired by a result on classical birth and death processes ([22] Theorem 2 and Lemma 1), we consider the operator

[TABLE]

which is clearly bounded since $\sum_{k\geq 0}(k+1)^{-1}(k+r)^{-1}<\infty$ and is a function of the number operator. A straightforward computation shows that

[TABLE]

It follows that $X$ belongs to the domain of $\mathcal{L}$ ([12] Prop. 3.33 p.64) and

[TABLE]

so that $\mathcal{T}_{t}(X)\leq X$ for all $t\geq 0$. Since the QMS $\mathcal{T}$ is transient, it has no invariant state by Proposition 6 of Ref.[15]. $\square$

In the case where there is a faithful invariant state, it is not difficult to show that the quantum detailed balance condition (in most of its quantum formulations as in [17]) holds.

4 Relationship with the two-photon absorption and emission process

The two-photon absorption and emission process quantum Markov semigroup is generated by

[TABLE]

where $a,a^{+}$ are the usual annihilation and creation operators in $\mathsf{h}$, $\lambda,\mu,\xi^{+},\xi^{-}\in\mathbb{R}$, $\lambda\leq\mu$.

This generator has been studied in Ref.[16] for $\xi^{+}=0$, however, this does not change any result of that paper. In particular, it has been proved that the even and odd projections

[TABLE]

are harmonic (i.e. invariant) for the QMS generated by $\mathcal{L}_{\rm tp}$. As a consequence we can consider the restricted semigroups $\mathcal{T}^{\rm e}$ and $\mathcal{T}^{\rm o}$ on the von Neumann subalgebras $p_{\rm e}\mathcal{B}(\mathsf{h})p_{\rm e}$ and $p_{\rm o}\mathcal{B}(\mathsf{h})p_{\rm o}$, identified respectively with $\mathcal{B}(p_{\rm e}\mathsf{h})$ and $\mathcal{B}(p_{\rm o}\mathsf{h})$.

Let $U_{\rm e},U_{\rm o}$ be the unitary operators

[TABLE]

A straightforward computation shows that, if $r=1/2$,

[TABLE]

and, if $r=3/2$, similarly

[TABLE]

so that, in both cases,

[TABLE]

As a consequence, the quadratic open quantum harmonic oscillator generator $\mathcal{L}$, for $\xi^{\pm}=4\zeta^{\pm}$, satisfies

[TABLE]

This shows that the QMS of the quadratic open quantum harmonic oscillator is unitarily equivalent (up to the multiplicative constant $4$) to the QMS obtained by restriction of the two-photon absorption and emission process to the even (resp. odd) states of the number operator for $r=1/2$ (resp. $r=3/2$), for a suitable choice of the parameters $\xi^{\pm},\zeta^{\pm}$. This analogy will serve as an inspiration to study the exponential speed of convergence towards the equilibrium state.

5 Instantaneous spread of state supports

In this section we will show that for all initial state $\eta$ the support of the state $\mathcal{T}_{*t}(\eta)$ at any time $t>0$ is the whole of $\mathsf{h}$.

This property follows from a recent result by Hachicha, Nasroui and Gliouez [20] Theorem 3.3 for QMSs associated with operators $G,L_{\ell}$ (in our case $L_{1}=\lambda B^{+},L_{2}=\mu B$) constructed form generators.

Theorem 2

Suppose that $G$ generates an analytic semigroup in a sector $\Delta=\{z\in\mathbb{C}\mid|\hbox{\rm arg }z|<\theta\}$ with $\theta\in]0,\frac{\pi}{2}]$ and, moreover, that

[TABLE]

for all $k>0.$ For all state $\omega=\sum_{j\in J}\omega_{j}|\psi_{j}\rangle\langle\psi_{j}|,$ with $\omega_{j}>0$ for all $j\in J$ and all $t\geq 0,$ the support projection $\mathcal{S}_{t}(\omega)$ of the state $\mathcal{T}_{*t}(\omega)$ is the closure of linear manifold generated by vectors

[TABLE]

for all $j\in J$, $n\geq 1,$ $m_{1},\cdots,m_{n}\geq 0$ and $\ell_{1},\cdots,\ell_{n}\geq 1$ where $\delta_{G}^{m}$ denotes the $m$-th iterated commutator with $G$ and $\delta_{G}^{0}(L_{\ell})=L_{\ell}$.

In our framework the operator $G$ can be written as

[TABLE]

Note that the self-adjoint operator $N^{2}$ generates a semigroup $(e^{zN^{2}})_{z\in\mathbb{C}-\mathbb{R}_{+}}$ defined on the complex plane without the positive real half axis which is an analytic semigroup in the half plane $\{z\in\mathbb{C}\,\mid\,\Re(z)<0\}$. Thus, by a change of variable, the operator

[TABLE]

generates an analytic semigroup in the sector

[TABLE]

which is equivalent to $\Re(z)>0$ for $\zeta^{+}+\zeta^{-}=0$ and

[TABLE]

In any case the semigroup generated by the operator $G_{0}$ in (10) is analytic in the sector

[TABLE]

with the convention $\arctan(+\infty)=\pi/2$. Clearly $\hbox{\rm Dom}(G)=\hbox{\rm Dom}(G_{0})=\hbox{\rm Dom}(N^{2})$ and $G$ is the sum of $G_{0}$ and an operator with domain $\hbox{\rm Dom}(N)=\hbox{\rm Dom}(G_{0}^{1/2})$ where $G_{0}^{1/2}$ is defined as $N$ times a square root of the complex number $(\lambda^{2}+\mu^{2})/2+\mathrm{i}(\zeta^{+}+\zeta^{-})$. It follows then from Corollary 2.4 p.81 of [24] that $G$ generates an analytic semigroup in the sector (11).

We can now prove the following.

Theorem 3

For all initial state $\eta$ the support of the state $\mathcal{T}_{*t}(\eta)$ at any time $t>0$ is the whole of $\mathsf{h}$.

Proof. The assumption (8) obviously holds because $\hbox{\rm Dom}(G^{k})=\hbox{\rm Dom}(N^{2k})$ and $B,B^{+}$ map $\hbox{\rm Dom}(N^{h})$ in $\hbox{\rm Dom}(N^{h-1})$ for all $h$ so that

[TABLE]

For all state $\omega$ we can write its spectral decomposition $\omega=\sum_{j\in J}\omega_{j}|\psi_{j}\rangle\langle\psi_{j}|$, for a collection $(\psi_{j})_{j\in J}$ of orthonormal vectors and $\omega_{j}>0$ for all $j\in J$ and $J$ non-empty. Since $L_{1}=\lambda B^{+},L_{2}=\mu B^{-}$, $\delta_{G}^{0}(B^{\pm})=B^{\pm}$ and $[B,B^{+}]-r\,\mathbb{1}=2N$, by Theorem 2, all vectors

[TABLE]

belong to the support of $\mathcal{T}_{*t}(\omega)$. Write $P_{t}\psi_{j}=\sum_{k\geq 0}z_{jk}(t)e_{k}$ and let $m_{j}$ be the minimum $k$ for which $z_{jk}(t)\not=0$. Since the function $z\to\hbox{\rm e}^{-zN}$ is analytic for $\Re z>0$, and $P_{t}\psi_{j}$ belongs to the support of $\mathcal{T}_{*t}(\omega)$, for all real number $s>0$

[TABLE]

belongs to the support of $\mathcal{T}_{*t}(\omega)$ by (12) for all $s>0$. Taking the limit as $s$ tends to infinity, we conclude that $e_{m_{j}}$ belongs to the support of $\mathcal{T}_{*t}(\eta)$. In the same way, starting from either

[TABLE]

we can conclude that

[TABLE]

belong to the support of $\mathcal{T}_{*t}(\omega)$ for $t>0$ for all $k$. As a consequence, the support of $\mathcal{T}_{*t}(\omega)$ is the whole of $\mathsf{h}$ for all $t>0$.

6 Spectral gap

In this section we discuss the spectral gap of the generator of the semigroup of the quadratic harmonic oscillator. For the purpose we will follow the methods developed in Refs. [7, 8].

6.1 Dirichlet form and spectral gap

Recall the invariant state $\rho$ in (6). Let $L_{2}(\mathsf{h})$ be the space of Hilbert-Schmidt operators on $\mathsf{h}$ with inner product $\langle x,y\rangle=\hbox{\rm{tr}}(x^{*}y)$. Consider the embedding

[TABLE]

Let $T=(T_{t})_{t\geq 0}$ be the strongly continuous contraction semigroup on $L_{2}(\mathsf{h})$ defined by

[TABLE]

Let $L$ be the generator of the semigoup $(T_{t})_{t\geq 0}$. We can check that

[TABLE]

Lemma 1

For $x=\sum_{j,k}x_{jk}|e_{j}\rangle\langle e_{k}|\in\mathcal{M}$,

[TABLE]

with the convention $x_{j-1\,k-1}=0$ if $j=0$ or $k=0$.

The Dirichlet form, defined for $\xi\in D(L)$, is the quadratic form $\mathcal{E}$

[TABLE]

The spectral gap of the operator $L$ is the nonnegative number

[TABLE]

Lemma 2

For $\xi\in D(L)$ with $\xi=\sum_{j,k}\xi_{jk}|e_{j}\rangle\langle e_{k}|$

[TABLE]

In particular, $\mathcal{E}(\xi)=\mathcal{E}(\xi^{*})$, where $\xi^{*}$ is the adjoint operator of $\xi$.

**Proof. ** By Lemma 1 we get

[TABLE]

Rearranging the terms we get the desired expression. $\square$

Proposition 5

If $W=\{\xi\in D(L):\mathcal{E}(\xi)=0\}$ then

[TABLE]

Proof. Since $\mathbb{1}\in\mathrm{Ker}\,\mathcal{L}$, $\rho^{1/2}=\iota(\mathbb{1})\in\mathrm{Ker}\,L$. It is obvious that $\mathrm{Ker}\,L\subset W$. Suppose $\xi=\sum_{j,k}\xi_{jk}|e_{j}\rangle\langle e_{k}|\in W$. By Lemma 2

[TABLE]

and

[TABLE]

Thus $\xi$ must be diagonal and by (14), $\xi_{k+1\,k+1}=\nu\xi_{kk}$. Hence for $k\geq 0$, $\xi_{kk}=\nu^{k}\xi_{00}$, i.e., $\xi=\text{const.}\rho^{1/2}$. This completes the proof. $\square$

Like the model of two-photon absorption and emission process discussed in [8], there are invariant subspaces for the process of quadratic open quantum harmonic oscillator. For $m\in\mathbb{Z}$, define

[TABLE]

One can easily check the following properties.

1. (1)

$\xi\in\mathcal{G}_{m}\Leftrightarrow\xi^{*}\in\mathcal{G}_{-m}$. 2. (2)

The linear spaces $\mathcal{G}_{m}$ are orthogonal in $L_{2}(\mathsf{h})$, and

[TABLE] 3. (3)

Each $\mathcal{G}_{m}$ is invariant for the generator $L$ and so also for the semigroup $T$. 4. (4)

Each $\mathcal{G}_{m}$ is isometrically isomorphic to the space $l^{2}(\mathbb{N})$ of square summable sequences.

By mimicking the proof of [8, Proposition 4], we can show the following.

Proposition 6

[TABLE]

where $A_{m}=\begin{cases}\inf\{\mathcal{E}(\xi):\,\|\xi\|=1,\,\,\xi\in\mathcal{G}_{m}\},&\mathrm{for}\,\,m\neq 0\\ \inf\{\mathcal{E}(\xi):\,\|\xi\|=1,\,\,\xi\in\mathcal{G}_{0},\,\,\xi\perp\rho^{1/2}\},&\mathrm{for}\,\,m=0\end{cases}$.

We will now study separately off-diagonal and diagonal minima.

6.2 Off-diagonal minimum problems

Fix $m>0$. For $\xi\in\mathcal{G}_{m}$ we can write for some sequence $y=(y_{j})_{j\geq 0}$ in $l^{2}(\mathbb{N})$,

[TABLE]

where $y_{j}=\xi_{j\,j+m}$. Then

[TABLE]

Proposition 7

For any $m\geq 1$,

[TABLE]

The lower bound is attained by $\mathcal{E}(\xi)$ for

[TABLE]

and $y_{0}\geq 0$ such that $\|\xi\|^{2}=1$.

Proof. We fix $m\geq 1$ and $\xi\in\mathcal{G}_{m}$. From the formula (15), without loss of generality we may assume $y_{j}\geq 0$ for all $j\geq 0$ because $|a-b|\geq||a|-|b||$ for any complex numbers $a$ and $b$. Then we can rewrite

[TABLE]

Since

[TABLE]

where $\theta_{j}>0$ is such that

[TABLE]

i.e.

[TABLE]

In this way, we get the inequality

[TABLE]

It follows that

[TABLE]

The above lower bound for the Dirichlet form is attained if and only if the Schwarz inequalities (18) are equalities namely $\mu\theta_{j}^{1/2}y_{j+1}=\lambda\theta_{j}^{-1/2}y_{j}$ so that

[TABLE]

for all $j\geq 1$. Iterating we find

[TABLE]

Since

[TABLE]

we find $\sum_{j\geq 0}y_{j}^{2}<\infty$ and so the lower bound is a minimum.

Minimizers can be written in an explicit form. First note that

[TABLE]

Iterating

[TABLE]

$\square$

The function $m\to(m+r-1)+\nu^{2}(m-r+1)$ is a positive and increasing function of $m$ for $m\geq 1$ and so we have

Theorem 4

The off-diagonal gap is

[TABLE]

6.3 Diagonal minima

For any $f$ in $l^{2}(\mathbb{N})$, let us denote by $M_{f}\in\mathcal{B}(\mathsf{h})$ the multiplication operator by $f$, $M_{f}=\sum_{j}f_{j}|e_{j}\rangle\langle e_{j}|$. Then we get

[TABLE]

The (formal) explicit expression for $\mathcal{A}$ is given by

[TABLE]

The invariant measure for this classical birth and death process is

[TABLE]

From now on, whenever there is no confusion we write simply $f$ for $M_{f}\in\mathcal{B}(\mathsf{h})$. By Lemma 2 we can compute

[TABLE]

Proposition 8

For any positive sequence $(a_{n})_{n}$, define the (strictly positive) constant

[TABLE]

Then $\mathrm{gap}(\mathcal{A})\geq B(\nu)^{-1}$.

Proof. We follow the proof of [8, Proposition 7] with a change $(2u+1)(2u+2)$ by $\omega_{u+1}$, which amounts to consider the birth rate $b_{j}=\lambda^{2}\omega_{j+1}$ instead of $b_{j}=\lambda^{2}(2j+1)(2j+2)$. $\square$

The following was proven in [8, Lemma 8].

Lemma 3

Take a positive summable sequence $(a_{n})_{n\geq 0}$ and define the positive decreasing tail sequence $(A_{k})_{k\geq 0}$ by $A_{k}=\sum_{n\geq k}a_{k}$. Then

[TABLE]

Thus the computation of the spectral gap relies on how we choose the sequence $(a_{n})_{n}$. Here we propose the following choice. (cf. [8, Lemma 9])

Lemma 4

Let $a_{n}=\frac{1}{(n+r+1)(n+r)}$. Then the following properties hold.

1. (i)

$A_{k}=\sum_{n\geq k}a_{n}=\frac{1}{k+r}$.

** 2. (ii)

$\sum_{k\geq 0}\nu^{2k}a_{k}=2\nu^{-2(r+1)}\int_{0}^{\nu}\frac{s^{2r-1}}{1-s^{2}}(\nu^{2}-s^{2})ds$.

** 3. (iii)

$\sum_{k\geq 0}\nu^{2k}A_{k}=2\nu^{-2r}\int_{0}^{\nu}\frac{s^{2r-1}}{1-s^{2}}ds$.

Proof. The item (i) is trivial. To prove (ii), we write

[TABLE]

Let

[TABLE]

Then, $b(\nu):=a(\nu)\nu^{2r}$ satisfies

[TABLE]

Since $b(0)=0$, we get

[TABLE]

From this we easily get the desired expression. For (iii), we see by (i) that the value we are looking for is $a(\nu)$, which we computed in the above. $\square$

We can now find a lower bound for the diagonal minimum of the Dirichlet form on vectors $\xi$ orthogonal to $\rho^{1/2}$.

Theorem 5

[TABLE]

Proof. We choose the sequence $(a_{n})_{n}$ as in Lemma 4. Let us define a function $V:\mathbb{N}\to\mathbb{R}$ by

[TABLE]

We can see that the function $V(u)$ is increasing. In fact, regarding $V(u)$ as a function on the interval $[0,\infty)$, we differentiate it. With a little computation we see that

[TABLE]

Therefore we get by using Lemma 4 (iii)

[TABLE]

The result now follows by Proposition 8. $\square$

Remark. The lower bound of theorem 5 can be written in a closed form by introducing the Lerch function $\Phi:\left\{\,z\in\mathbb{C}\,\mid\,|z|<1\,\right\}\times\mathbb{N}\times]0,+\infty[\to\mathbb{C}$

[TABLE]

By Theorem 5, we have $\mathrm{gap}(\mathcal{A})\geq\mu^{2}\left(\Phi(\nu^{2},1,r)\right)^{-1}$. Moreover, by comparing the off-diagonal explicit minimum and the diagonal lower bound we find the following.

Theorem 6

For all $\nu,r$ such that $\left(2\nu^{2}+r(1-\nu^{2})\right)\Phi(\nu^{2},1,r)\leq 2$ we have

[TABLE]

In particular, if $r\geq 2\nu^{2}/(1-\nu^{2})$ the above identity holds.

Proof. The first claim follows immediately by comparing the diagonal lower bound of Theorem 5 and the off-diagonal minimum of Theorem 4. By the elementary inequality,

[TABLE]

we have $\Phi(\nu^{2},1,r)\leq\left(r(1-\nu^{2})\right)^{-1}$ and so the identity (19) holds, in particular, if $\left(r(1-\nu^{2})\right)^{-1}\left(2\nu^{2}+r(1-\nu^{2})\right)\leq 2$, i.e. $r\geq 2\nu^{2}/(1-\nu^{2})$. $\square$

The graph in Figure 1 shows the values of $\nu,r$ for which the identity $\left(2\nu^{2}+r(1-\nu^{2})\right)\Phi(\nu^{2},1,r)=2$ holds. Clearly, for pairs $(\nu,r)$ lying above (or on) the graph the spectral gap is given by (19).

The exact value of the spectral gap, for pairs $(\nu,r)$ lying below the graph, could be also the diagonal minimum whose exact value is not known and we are unable to compute. This will be clear studying upper bounds for the diagonal spectral gap.

6.4 Upper bound

In this section we discuss the upper bound of the diagonal spectral gap. By definition, any value $\mathcal{E}(\xi)$ with $\|\xi\|=1$, with $\xi$ orthogonal to $\rho^{1/2}$ is an upper bound for the gap. Thus to get a good upper bound we need to choose a vector $\xi$ cleverly.

Looking at the explicit form of the off-diagonal minimizer that we get for $m=1$, we consider the first order polynomial $f_{j}=j-c$ where $c$ is a constant chosen in such a way that $j\to\pi_{j}^{1/2}f_{j}$ is orthogonal to $\pi_{j}^{1/2}$ i.e. $\sum_{j}\pi_{j}(j-c)=0$. It is worth noticing here that this choice yields the desired minimizer of the Dirichlet form on vectors orthogonal to $\rho^{1/2}$ for the usual harmonic oscillator (see [7]). Using the identities

[TABLE]

one computes $c=\nu^{2}/(1-\nu^{2})$. Then, considering $\xi_{j}=\pi^{1/2}_{j}(j-c)$, one finds

[TABLE]

and so

[TABLE]

As a consequence, one has the explicit upper bound

[TABLE]

which is twice the spectral gap in good cases by Theorem 6.

In order to show that the spectral gap converges to [math] as $r\to 0$, we find another upper bound computing the value of the Dirichlet form for another vector $\xi$ suggested by our choice of the sequence $(A_{k})_{k}$ in Lemma 4.

Theorem 7

[TABLE]

where $c(\nu,r)=(1-\nu^{2})\sum_{j=0}^{\infty}\frac{\nu^{2j}}{j+r}=(1-\nu^{2})\Phi(\nu^{2},1,r)$. In particular, for all $\nu>0$ fixed, $\mathrm{gap}\,L$ tends to [math] as $r\to 0$.

Proof. We consider $A_{k}=(k+r)^{-1}$ as in Lemma 4, and define $f_{k}=A_{k}-c=(k+r)^{-1}-c(\nu,r)$, where the constant $c$ is such that $\sum_{j}\pi_{j}f_{j}=0$. We compute

[TABLE]

On the other hand, $\|f\|_{\pi}^{2}=(1-\nu^{2})\sum_{j\geq 0}\left((j+r)^{-1}-c(\nu,r)\right)^{2}\nu^{2j}.$ Therefore,

[TABLE]

We now check that the above upper bound tends to [math] as $r\to 0$. To this end note that we can write the denominator as

[TABLE]

It follows that, if we multiply the denominator by $r^{2}$ in the limit as $r\to 0$ we get $1-(1-\nu^{2})=\nu^{2}$. On the other hand, if we multiply also the numerator by $r^{2}$ we get

[TABLE]

This completes the proof. $\square$

Conclusion. Theorem 7 shows that, for $r$ near [math], $\mathrm{gap}\,L$ is the diagonal gap whose exact value is not known. Moreover, it tends to [math] as $r\to 0$. The exact result of Theorem 6 holds in the white region above the graph in Figure 1. It is worth noticing here that the range of values of $r$ for which our exact result holds becomes closer and closer to the half-line $(0,+\infty)$ as $\nu^{2}=\hbox{\rm e}^{-s\beta}$ goes to [math], namely the inverse temperature $\beta$ goes to $+\infty$. This confirms the intuition that quantum (off-diagonal) effects prevail over the classical (diagonal) ones when the temperature is small.

Acknowledgement

A. Dhahri acknowledges support by Fondo professori stranieri DHG9VARI01 Politecnico di Milano. The research by H. J. Yoo was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2016R1D1A1B03936006).