arxiv: v1 [math-ph] 7 Aug 2013

Transcription

1 Optimal Hamiltonian of Fermion Flows Luigi Accardi 1, Andreas Boukas 2 arxiv: v1 [math-ph] 7 Aug Centro Vito Volterra Università di Roma Tor Vergata via Columbia, Roma, Italia accardi@volterra.mat.uniroma2.it 2 Department of Mathematics and Natural Sciences American College of Greece Aghia Paraskevi, Athens 15342, Greece andreasboukas@acgmail.gr Abstract After providing a general formulation of Fermion flows within the context of Hudson-Parthasarathy quantum stochastic calculus, we consider the problem of determining the noise coefficients of the Hamiltonian associated with a Fermion flow so as to minimize a naturally associated quadratic performance functional. This extends to Fermion flows results of the authors previously obtained for Boson flows. fermion.tex 1 Introduction In the Heisenberg picture of quantum mechanics, the time-evolution of an observable X is described by the operator process j t (X = U (tx U(t where U(t = e ith and the Hamiltonian H is a self-adjoint operator on the wave function Hilbert space. The unitary process U(t and the self-adjoint process j t (X satisfy, respectively, du(t = ihu(tdt, U( = I and dj t (X = i[h,j t (X]dt, j (X = X where [H,j t (X] := Hj t (X j t (XH In the presence of quantum noise, the equation satisfied by U(t is replaced by a Hudson-Parthasarathy quantum stochastic differential equation driven by 1

2 that noise (see [15] and [17] and the corresponding equation for the quantum flow j t (X is interpreted as the Heisenberg picture of the Schrödinger equation in the presence of noise. The emergence of such equations as stochastic limits of classical Schrödinger equations is described in [8]. From the point of view of quantum control theory, the problem of minimizing a quadratic performance functional associated with a Hudson-Parthasarathy quantum flow driven by Boson quantum noise, has been considered in [1], [3]-[5] and [9]-[13]. In this paper we consider the same problem for quantum flows driven by Fermion flows. Our approach is based on the representation of the Fermion commutation relations on Boson Fock space obtained in [6] and [16]. A unified approach to the quadratic cost control of quantum processes driven by Boson, Fermion, Finite- Difference and a wide class of other quantum noises, can be found in [2] with the use of the representation free quantum stochastic calculus of [7]. This paper is structured as follows: In section 2 we describe the main features of Hudson-Parthasarathy calculus and the representation of the Fermion commutation relations on Boson Fock space. In section 3 we obtain the quantum stochastic differential equations satisfied by Fermion flows. In section 4 we obtain a new algebraic formulation of Fermion evolution equations and flows and we describe the equations satisfied by their structure maps. Finally, in section 5 we define the quadratic performance functionals associated with Fermion flows and we obtain the coefficients of the corresponding Hamiltonian for which these functionals are minimized. 2 Quantum Stochastic Calculus The Boson Fock space Γ := Γ(L 2 (R +,C over L 2 (R +,C is the Hilbert space completion of the linear span of the exponential vectors ψ(f under the inner product < ψ(f,ψ(g >:= e <f,g> where f,g L 2 (R +,C and < f,g >= + f(sg(sds where, here and in what follows, z denotes the complex conjugate of z C. The annihilation, creation and conservation operator processes A t, A t and Λ t respectively, are defined on the exponential vectors ψ(g of Γ by A t ψ(g := t g(sds ψ(g A tψ(g := ǫ ǫ=ψ(g +ǫχ [,t] Λ t ψ(g := ǫ ǫ=ψ(e ǫχ [,t] g. 2

3 The basic quantum stochastic differentials da t, da t, and dλ t are defined by da t := A t+dt A t da t := A t+dt A t dλ t := Λ t+dt Λ t. Hudson and Parthasarathy defined in [15] stochastic integration with respect to the noise differentials da t, da t and dλ t and obtained the Itô multiplication table da t dλ t da t dt da t dλ t da t dλ t da t dt da t dt. We couple Γ with a system Hilbert space H and consider processes defined on H Γ. The fundamental theorems of the Hudson-Parthasarathy quantum stochastic calculus give formulas for expressing the matrix elements of quantum stochastic integrals in terms of ordinary Riemann-Lebesgue integrals. Theorem 1. Let M(t := t E(sdΛ s +F(sdA s +G(sdA s +H(sds where E, F, G, H are (in general time dependent adapted processes. Let also u ψ(f and u ψ(g be in the exponential domain of H Γ. Then < u ψ(f,m(tu ψ(g >= t < u ψ(f,( f(sg(se(s+g(sf(s+ f(sg(s+h(s u ψ(g > ds. Proof. See Theorem 4.1 of [15]. Theorem 2. Let and M(t := t E(sdΛ s +F(sdA s +G(sdA s +H(sds M (t := t E (sdλ s +F (sda s +G (sda s +H (sds where E, F, G, H, E, F, G, H are (in general time dependent adapted processes. Let also u ψ(f and u ψ(g be in the exponential domain of H Γ. Then < M(tu ψ(f,m (tu ψ(g >= t {< M(su ψ(f,( f(sg(se (s+g(sf (s+ f(sg (s+h (s u ψ(g > + < (ḡ(sf(se(s+f(sf(s+ḡ(sg(s+h(su ψ(f,m (su ψ(g > + < (f(se(s+g(su ψ(f,(g(se (s+g (su ψ(g >}ds. 3

4 Proof. See Theorem 4.3 of [15]. The connection between classical and quantum stochastic analysis is given in the following: Theorem 3. The processes B = {B t, t } and P = {P t, t } defined by and B t := A t +A t P t := Λ t + λ(a t +A t +λt are identified with Brownian motion and Poisson process of intensity λ respectively, in the sense that their vacuum characteristic functionals are given by and Proof. See Theorem 5 of [13]. < ψ(,e isbt ψ( >= e s2 2 t < ψ(,e ispt ψ( >= e λ(eis 1t. The processes A t, A t satisfy the Boson Commutation Relations [ ] A t,a t := A t A t A t A t = ti. In [16] Hudson and Parthasarathy showed that the processes F t and F t defined on the Boson Fock space by F t := F t := t t J s da s J s da s satisfy the Fermion anti-commutation relations It follows that {F t,f t } := F tf t +F t F t = ti. df t = J t da t (2.1 df t = J t da t. (2.2 Here J t is the self-adjoint, unitary-valued, adapted, so called reflection process, acting on the noise part of the Fock space and extended as the identity on the system part, defined by J t := γ( P [,t] +P (t,+ 4

5 where P S denotes the multiplication operator by χ S and γ is the second quantization operator defined by γ(uψ(f := ψ(u f. The reflection process J t commutes with system space operators and satisfies the differential equation (cf. Lemma 3.1 of [16] dj t = 2J t dλ t (2.3 J = 1. 3 Fermion Evolutions and Flows As shown in [16], see also [6], Fermion unitary evolution equations have the form (( du t = ih L L dt+l W df t LdF t +(1 W dλ t U t (3.1 with adjoint du t = U t (( ih L Ldt L df t +W LdF t +(1 W dλ t (3.2 where, U = U = 1 and, for each t, U t is a unitary operator defined on the tensor product H Γ(L 2 (R +,C of the system Hilbert space H and the noise Fock space Γ. Here H, L, W are in B(H, the space of bounded linear operators on H, with W unitary and H self-adjoint. We identify timeindependent, bounded, system space operators X with their ampliation X 1 to H Γ(L 2 (R +,C. Proposition 1. Let φ t (T,S := V t (T +SJ tv t where T,S are bounded system space operators and V t,v t satisfy the quantum stochastic differential equations dv t = = ( αdt+βdf t +γdf t +δdλ t V t (αdt+βj t da t +γj t da t +δdλ t V t and dv t = Vt (α dt+β df t +γ df t +δ dλ t = Vt (α dt+β J t da t +γ J t da t +δ dλ t 5

6 where α, β, γ, δ are bounded system space operators. Then with dφ t (T,S = φ t (α T +T α+γ T γ,α S +Sα+γ Sγdt (3.3 + φ t (γ S +Sβ γ S(2+δ,γ T +T β +γ T δda t + φ t (β S +Sγ δ Sγ,β T +T γ +δ T γda t + φ t (δ T +T δ +δ T δ, (2S +Sδ +δ S +δ SδdΛ t Proof. Making use of the algebraic rule we find φ (T,S = T +S (3.4 d(xy = dxy +xdy +dxdy dφ t (T,S = dvt (T +SJ tv t +Vt d((t +SJ tv t +dvt d((t +SJ tv t = dvt (T +SJ t V t +Vt T dv t +Vt Sd(J t V t +dvt T dv t +dvt Sd(J t V t. But, by (2.3, the Itô table for the Boson stochastic differentials and the fact that J 2 t = 1 Thus d(j t V t = dj t V t +J t dv t +dj t dv t = 2J t dλ t V t +J t (αdt+βj t da t +γj t da t +δdλ t ( 2J t dλ t αdt+βj t da t +γj t da t +δdλ t V t = (αj t dt+βda t γda t (2+δJ t dλ t V t. dφ t (T,S = Vt (α dt+β J t da t +γ J t da t +δ dλ t (T +SJ t V t +Vt (αdt+βj T t da t +γj t da t +δdλ t V t +Vt S (αj t dt+βda t γda t (2+δJ t dλ t V t (α dt+β J t da t +γ J t da t +δ dλ t +V t ( T αdt+βj t da t +γj t da t +δdλ t V t +Vt (α dt+β J t da t +γ J t da t +δ dλ t S (αj t dt+βda t γda t (2+δJ t dλ t = φ t (α T +T α+γ T γ,α S +Sα+γ Sγdt +φ t (γ S +Sβ γ S(2+δ,γ T +T β +γ T δda t +φ t (β S +Sγ δ Sγ,β T +T γ +δ T γda t +φ t (δ T +T δ +δ T δ, (2S +Sδ +δ S +δ SδdΛ t V t V t 6

7 With the processes U t and U t of (3.1 and (3.2 we associate the Fermion flow and the reflected flow where X is a bounded system space operator. j t (X := U t XU t = φ(x, (3.5 r t (X := j t (XJ t = φ(,x (3.6 Corollary 1. The Fermion flow j t (X and the reflected flow r t (X defined in (3.5 and (3.6 satisfy the system of quantum stochastic differential equations and with dj t (X = j t ( i [H,X] 1 2 (L LX +XL L 2L XL dt (3.7 +r t ([L,X] W da t +r t (W [X,L] da t +j t(w XW X dλ t dr t (X = r t ( i [H,X] 1 2 (L LX +XL L 2L XL dt (3.8 j t ({L,X}W da t j t (W {X,L} 2XL da t r t(w XW +X dλ t j (X = r (X = X. Here, as usual, [x,y] := xy yx and {x,y} := xy +yx. Proof. We replace V t and V t in Proposition 1 by U t and U t. Then, equation (3.7 is a special case of (3.3 for T = X, S = and α = (ih L L β = L W γ = L δ = W 1 while equation (3.8 is a special case of (3.3 for T =, S = X and α,β,γ,δ as above. 4 Generalized Fermion Flows Fermion flows can be formulated and studied in a manner similar to Boson (also called Evans-Hudson flows (cf. [14] and [17]. Let B(H denote the space of 7

8 bounded system operators. We define B(H-valued operations and on B(H B(H by Notice that if is the reflection map, then (T 1,S 1 (T 2,S 2 : = T 1 T 2 +S 1 S 2 (4.1 (T 1,S 1 (T 2,S 2 : = T 1 S 2 +S 1 T 2. (4.2 ρ(x,y := (y,x (T 1,S 1 (T 2,S 2 = (T 1,S 1 ρ(t 2,S 2. We also define the B(H B(H-valued product map on (B(H B(H 2 by x y := (x y,x y = (T 1 T 2 +S 1 S 2,T 1 S 2 +S 1 T 2 (4.3 where x = (T 1,S 1 and y = (T 2,S 2. From the definitions it follows that the -product is associative with unit id := (1, where 1 and are the identity and zero operators in B(H. Let the flow φ t (T,S := Vt (T +SJ t V t be as in Proposition 1 and let x = (T 1,S 1 and y = (T 2,S 2. Then φ t (xφ t (y = φ t (T 1,S 1 φ t (T 2,S 2 = V t (T 1 +S 1 J t V t V t (T 2 +S 2 J t V t = V t (T 1 +S 1 J t (T 2 +S 2 J t V t = V t (T 1T 2 +S 1 S 2 +(T 1 S 2 +S 1 T 2 J t V t = φ t (T 1 T 2 +S 1 S 2,T 1 S 2 +S 1 T 2 = φ t (x y i.e φ t is a homomorphism with respect to the -product. Since φ t is the solution of (3.3, (3.4 this suggests considering flows satisfying quantum stochastic differential equations of the form dφ t (x = φ t (θ 1 (xdt+φ t (θ 2 (xda t +φ t (θ 3 (xda t +φ t (θ 4 (xdλ t (4.4 with φ (x = x id+x id where x = (T,S B(H B(H and for i = 1,2,3,4, θ i : B(H B(H B(H B(H are the structure maps. In the case of (3.3, (3.4 θ 1 (T,S = (α T +T α+γ T γ,α S +Sα+γ Sγ θ 2 (T,S = (γ S +Sβ γ S(2+δ,γ T +T β +γ T δ θ 3 (T,S = (β S +Sγ δ Sγ,β T +T γ +δ T γ θ 4 (T,S = (δ T +T δ +δ T δ, (2S +Sδ +δ S +δ Sδ 8

9 where α = (ih L L β = L W γ = L δ = W 1. The general conditions on the θ i in order for φ t to be an identity preserving -product homomorphism are given in the following: Proposition 2. Let φ t be the solution of (4.4. Then φ t (xφ t (y = φ t (x y φ t (x = φ t (x φ t (id = 1 if and only if the θ i satisfy the structure equations θ 1 (x y +x θ 1 (y+θ 2 (x θ 3 (y = θ 1 (x y (4.5 θ 2 (x y +x θ 2 (y+θ 2 (x θ 4 (y = θ 2 (x y (4.6 θ 3 (x y +x θ 3 (y+θ 3 (x θ 4 (y = θ 3 (x y (4.7 θ 4 (x y +x θ 4 (y+θ 4 (x θ 4 (y = θ 4 (x y (4.8 and, with denoting adjoint and x = (T 1,S 1 x = (T 1,S 1, with Proof. θ 1 (x = (θ 1 (x (4.9 θ 2 (x = (θ 3 (x (4.1 θ 3 (x = (θ 2 (x (4.11 θ 4 (x = (θ 4 (x (4.12 θ 1 (id = θ 2 (id = θ 3 (id = θ 4 (id = (,. (4.13 φ t (x y = φ t (xφ t (y dφ t (x y = dφ t (xφ t (y+φ t (xdφ t (y+dφ t (xdφ t (y which implies that φ t (θ 1 (x ydt+φ t (θ 2 (x yda t +φ t (θ 3 (x yda t +φ t (θ 4 (x ydλ t = φ t (θ 1 (xφ t (ydt+φ t (θ 2 (xφ t (yda t +φ t (θ 3 (xφ t (yda t +φ t(θ 4 (xφ t (ydλ t +φ t (xφ t (θ 1 (ydt +φ t (xφ t (θ 2 (yda t +φ t (xφ t (θ 3 (yda t +φ t(xφ t (θ 4 (ydλ t + φ t (θ 2 (xφ t (θ 3 (ydt+φ t (θ 2 (xφ t (θ 4 (yda t +φ t (θ 3 (xφ t (θ 4 (yda t +φ t(θ 4 (xφ t (θ 4 (ydλ t 9

10 from which collecting the dt,da t,da t and dλ t terms on each side, using the homomorphism property and then equating the coefficients of dt,da t,da t and dλ t on both sides we obtain (4.5-(4.8. Similarly i.e φ t (x = φ t (x dφ t (x = dφ t (x φ t (θ 1 (x dt+φ t (θ 2 (x da t +φ t (θ 3 (x da t +φ t (θ 4 (x dλ t = φ t (θ 1 (x dt+φ t (θ 2 (x da t +φ t (θ 3 (x da t +φ t (θ 4 (x dλ t and (4.9-(4.12 follow by equating the coefficients of dt,da t,da t and dλ t on both sides. Finally φ t (id = 1 dφ t (id = which by the linear independence of dt and the stochastic differentials implies ( Optimal Noise Coefficients As in [2] and [3], motivated by classical linear system control theory, we think of the self-adjoint operator H appearing in (3.1 as fixed and we consider the problem of determining the coefficients L and W of the noise part of the Hamiltonian of the evolution equation (3.1 that minimize the evolution performance functional Q ξ,t (u = [ XU t ξ L LU t ξ 2 ]dt+ 1 2 LU T ξ 2 (5.1 where T [,+, ξ is an arbitrary vector in the exponential domain of H Γ and X is a bounded self-adjoint system operator. By the unitarity of U t, U T,J t, and J T, (5.1 is the same as the Fermion flow performance functional J ξ,t (L,W = [ j t (Xξ j t(l Lξ 2 ]dt+ 1 2 j T(Lξ 2 (5.2 and the reflected flow performance functional R ξ,t (L,W = [ r t (Xξ r t(l Lξ 2 ]dt+ 1 2 r T(Lξ 2 (5.3 WeconsidertheproblemofminimizingthefunctionalsJ ξ,t (L,WandR ξ,t (L,W over all system operators L,W where L is bounded and W is unitary. The motivation behind the definition of the performance functionals (5.1, (5.2 and (5.3 can be found in the following theorem which is the quantum stochastic analogue of the classical linear regulator theorem. 1

11 Theorem 4. Let U = {U t, t } be a process satisfying the quantum stochastic differential equation du t = (F U t +u t dt+ψu t df t +ΦU t df t +ZU t dλ t, U = 1, t [,T] (5.4 with adjoint du t = (U t F +u t dt+u t Ψ df t +U t Φ df t +U t Z dλ t, U = 1 (5.5 where < T < +, the coefficients F, Ψ, Φ, Z are bounded operators on the system space H and u t := ΠU t for some positive bounded system operator Π. Then the functional Q ξ,t (u = [< U t ξ,x 2 U t ξ > + < u t ξ,u t ξ >]dt < u T ξ,u T ξ > (5.6 where X is a system space observable, identified with its ampliation X I to H Γ, is minimized over the set of feedback control processes of the form u t = ΠU t, by choosing Π to be a bounded, positive, self-adjoint system operator satisfying ΠF +F Π+Φ ΠΦ Π 2 +X 2 = (5.7 ΠΨ+Φ Π+Φ ΠZ = (5.8 ΠZ +Z Π+Z ΠZ =. (5.9 The minimum value is < ξ,πξ >. We recognize (5.7 as the algebraic Riccati equation. Proof. Let Then θ t =< ξ,u t ΠU t ξ >. (5.1 dθ t =< ξ,d(u t ΠU t ξ >=< ξ,(du t ΠU t +U t ΠdU t +du t ΠdU t ξ > (5.11 which, after replacing du t and du t by (5.4 and (5.5 respectively, and using (2.1, (2.2 and the Itô table, becomes dθ t =< ξ,u t ((F Π+ΠF +Φ ΠΦdt+(Φ Π+ΠΨ+Φ ΠZJ t da t (5.12 +(ΨΠ +ΠΦ+Z ΠΦJ t da t +(Z Π+ΠZ +Z ΠZdΛ t U t ξ > and by (5.7-(5.9 By (5.1 + < ξ,(u t ΠU t +U t Πu t dtξ >. dθ t =< ξ,u t (Π2 X 2 U t dtξ > + < ξ,(u t ΠU t +U t Πu tdtξ > (5.13 θ T θ =< ξ,u T ΠU T ξ > < ξ,πξ >. (

12 while by (5.13 θ T θ = (5.15 (< ξ,u t (Π 2 X 2 U t ξ > + < ξ,(u t ΠU t +U t Πu t ξ >dt By (5.14 and (5.15 < ξ,u T ΠU T ξ > < ξ,πξ >= (5.16 (< ξ,u t (Π 2 X 2 U t ξ > + < ξ,(u t ΠU t +U t Πu t ξ >dt. Thus Q ξ,t (u = (< ξ,u T ΠU T ξ > < ξ,πξ >+Q ξ,t (u (5.17 (< ξ,u T ΠU T ξ > < ξ,πξ >. Replacing the first parenthesis on the right hand side of (5.17 by (5.16, and Q ξ,t (uby (5.6 we obtain after cancelations Q ξ,t (u = (< ξ,(u t Π 2 U t +u t ΠU t +U t Πu t +u t u t ξ > dt + < ξ,πξ > (5.18 = (u t +ΠU t ξ 2 dt+ < ξ,πξ > which is clearly minimized by u t = ΠU t and the minimum is < ξ,πξ >. Theorem 5. Let X be a bounded self-adjoint system operator such that the pair (ih, X is stabilizable i.e there exists a bounded system operator K such that i H + KX is the generator of an asymptotically stable semigroup. Then, the quadratic performance functionals (5.2 and (5.3 associated with the Fermion flow {j t (X := U t XU t, t } and the reflected flow {r t (X := j t (XJ t, t }, where U = {U t, t } is the solution of (3.1, are minimized by and L = 2Π 1/2 W 1 (5.19 W = W 2 (5.2 where Π is a positive self-adjoint solution of the algebraic Riccati equation i[h,π]+π 2 +X 2 = (5.21 and W 1, W 2 are bounded unitary system operators commuting with Π. It is known (see [18] that if the pair (ih, X is stabilizable, then (5.21 has a positive self-adjoint solution Π. Moreover independent of T. min L,W J ξ,t (L,W = min L,W R ξ,t (L,W =< ξ,πξ > (

13 Proof. Looking at (3.1 as (5.4 with u t = 1 2 L LU t, F = ih, Ψ = L W, Φ = L, and Z = W 1, (5.6 is identical to (5.1. Moreover, equations (5.7- (5.9 become i[h,π]+l ΠL Π 2 +X 2 = (5.23 L Π ΠL W +L Π(W 1 = (5.24 (W 1Π+Π(W 1+(W 1Π(W 1 =. (5.25 By the self-adjointness of Π, (5.24 implies that while (5.25 implies that i.e (5.2. By (5.24 and the fact that in this case [L,Π] = [L,Π] = (5.26 [W,Π] = [W,Π] =. (5.27 Π = 1 2 L L i.e L L = 2Π (5.28 equation (5.23 implies (5.21. Equations (5.26 and (5.28 also imply that which implies (5.19. [L,L ] = (5.29 References [1] L. Accardi, A. Boukas, Control of Elementary Quantum Flows, Proceedings of the 5th IFAC symposium on nonlinear control systems, July 4-6, 21, St. Petersburg, Russia. [2], Quadratic Control of Quantum Processes, Russian Journal of Mathematical Physics, 9 (22, no. 4, 381-4, MR [3], Control of Quantum Langevin Equations, O pen Systems and Information Dynamics, 1 (23, no. 1, 89-14, MR [4], From classical to Quantum Quadratic Control, Proceedings of the International Winter School on Quantum Information and Complexity, January 6-1, 23, , Meijo university, Japan. [5], Control of Quantum Stochastic Differential Equations, Proceedings of the Fourth International Conference on System Identification and Control Problems (SICPRO 5, January 25-28, 25. [6] D. Applebaum, R. L. Hudson, Fermion Ito s Formula and Stochastic Evolutions, C omm. Math. Phys. 96, no.4, 1984 (

14 [7] L. Accardi, F. Fagnola, J. Quaegebeur, A representation free quantum stochastic calculus, J ournal of Functional Analysis, 14, no. 1, 1992 ( [8] L. Accardi., Y. G. Lu, I. Volovich, Q uantum theory and its stochastic limit, Springer Verlag, 22. [9] A. Boukas, Linear Quantum Stochastic Control, Q uantum Probability and related topics, QP -PQ IX, , World Scientific Publishing, River Edge, New Jersey, [1], Application of Operator Stochastic Calculus to an Optimal Control Problem, M at. Zametki, 53 (1993, no. 5, 48-56(Russian. Translation in Math. Notes, 53 (1993, no. 5-6, , MR 96a 817. [11], Operator Valued Stochastic Control in Fock Space with Applications to Noise Filtering and Orbit Tracking, J. Prob. Math. Stat., 16 (1994, no. 1. [12], Stochastic Control of Operator-Valued Processes in Boson Fock Space, Russian J. Math. Phys., 4 (1996, no. 2, , MR 97j [13], Application of Quantum Stochastic Calculus to Feedback Control, Global J. Pure and Applied Math., 1 (25, no. 1. [14] M. P. Evans, Existence of Quantum Diffusions, P robab. Th. Rel. Fields, 81, (1989, [15] R. L. Hudson, K.R. Parthasarathy, Quantum Ito s Formula and Stochastic Evolutions, C omm. Math. Phys., 93 (1984, [16], Unification of Fermion and Boson Stochastic Calculus, C omm. Math. Phys., 14, (1986, [17] K. R. Parthasarathy, An Introduction to Quantum Stochastic Calculus, Birkhauser, Boston, [18] A. J. Pritchard, R. F. Curtain, The Infinite Dimensional Riccati Equation, J. Math. Anal. Appl., 47 (