Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

CRC Press eBooks

The Quantization of the Motion of the Liquid

Isaak M. Khalatnikov

Published Nov 23, 2018• 0 citations
Publisher

Abstract

A classical liquid may be described by specifying the density ρ and the current vector j , which are defined in the following manner: ( 3-1 ) ρ = ∑ m α δ ( R − r α ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9780429502897/20d72e0b-9d39-4175-acc6-a943547edf8a/content/eq51.tif"/> ( 3-2 ) j = ∑ p α δ ( R − r α ) = ∑ m α v α δ ( R − r α ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9780429502897/20d72e0b-9d39-4175-acc6-a943547edf8a/content/eq52.tif"/>

DOI
10.1201/9780429502897-3

Related papers

Physical prerequisites

In-Gee Kim

The quantum revolution in the beginning of the 20th century changed our concepts of dynamics. The dynamic variables such as position r $ {\mathbf{r}} $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_1.tif"/> and momentum p $ {\mathbf{p}} $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_2.tif"/> are replaced by the corresponding position and momentum operators r ^ $ {\hat{\mathbf{r}}} $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_3.tif"/> and p ^ $ {\hat{\mathbf{p}}} $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_4.tif"/> , respectively When one would like to observe a dynamical variable, say ω , of a particle, one has to introduce a wavefunction ψ , which contains all the dynamical information of particles, and to apply the corresponding dynamical operator Ω ^ $ \hat{\Omega } $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_5.tif"/> to the wavefunction ψ . Then one may obtain the desired dynamical value of the particle as the eigenvalue of the operator, 1.1 Ω ^ ψ = ω ψ . $$ \hat{\Omega }\psi = \omega \psi . $$ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/math1_1.tif"/>

•Mar 14, 2018•0 citations

The Kinetic Equation for Elementary Excitations 13

Isaak M. Khalatnikov

The distribution function of the elementary excitations satisfies the kinetic equation ( 18-1 ) ∂ n ∂ t + ∂ n ∂ r ⋅ ∂ H ∂ p − ∂ n ∂ p ⋅ ∂ H ∂ r = ℐ ( n ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9780429502897/20d72e0b-9d39-4175-acc6-a943547edf8a/content/eq539.tif"/>

CRC Press eBooks•Nov 23, 2018•0 citations

Direct quantization of classical systems using the equation of motion and removing the degeneracy of Lagrangian formalism

Shahram Khosravi

Direct quantization of classical systems using the equation of motion and removing the degeneracy of Lagrangian formalism

Pizhūhish/hā-yi nuvīn-i fīzīk•Sep 1, 2017•0 citations

Velocity-dependent forces and non-hydrodynamic initial conditions in quantum and classical fluids

T. R. Kirkpatrick, D. Belitz

arXiv:2304.06157•The European Physical Journal Special Topics•Dec 1, 2023•3 citations

Quantum gravity with photon interactions, cavity resonators with inhomogeneities, classical and quantum optimal control of fields

Harish Parthasarathy

We wish j t ( X ) to track the noiseless trajectory X d ( t )∈ L ( h ). Assume non-demoltion measurements Y o ( t )= U ( t )* Y i ( t ) U ( t ) are made and the Belavkin filter for π t ( X ) =  j t ( X ) | η t ) , η t = σ ( Y 0 ( s ) : s ≤ t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1006.tif"/> has been constructed as d π t ( X ) = F t ( X ) d t + G t ( X ) d Y o ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1007.tif"/> trajectory estimation error X d ( t )− π t ( X ) is given as feedback to the state equations in the form d j t ( X ) = j t ( θ 0 ( X ) ) d t ⁢ + ⁢ j t ( θ 1 ( X ) ) d A ( t ) + j t ( θ 2 ( X ) ) d A ( t ) * + ⁢ j t ( θ 3 ( X ) ) d Λ ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1008.tif"/> 192 + ⁢ K ( t ) ( X d ( t ) − π t ( X ) ) d t https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1009.tif"/> Alternately, if we assume that X d follows a noiseless trajectory, ie, it evolves according to the equation X d ( t ) , = j t ( 0 ) ( X d ) = U 0 ( t ) * X d U 0 ( t ) , U 0 ( t ) = e x p ( − i t H 0 ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1010.tif"/> then we can give as feedback the estimation error E ( t ) = X d ( t ) − T r 2 ( π t ( X ) ( I ⊗ | ϕ ( u ) > < ϕ ( u ) | ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1011.tif"/> which is a system observable. This error feedback is given to the above state equations, ie, to the Evans-Hudson flow resulting in the dynamics d j t ( X ) = j t ( θ 0 ( X ) ) d t + j t ( θ 1 ( X ) ) d A ( t ) + j t ( θ 2 ( X ) ) d A ( t ) * + j t ( θ 3 ( X ) ) d Λ ( t ) + K ( t ) E ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1012.tif"/> Alternately, we can incorporate this error feedback in the original HP equation as d U ( t ) ⁢ = ⁢ ( - ( i H + ⁢ P + ⁢ K ( t ) E ( t ) ) d t + L 1 d A ( t ) ⁢ + ⁢ L 2 d A ( t ) * ⁢ + ⁢ S d Λ ( t ) ) U ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1013.tif"/> resulting in the Evans-Hudson flow the same as above but with θ 0 (X) replaced by θ 0 ( X )+ iK ( t )[ E ( t ), X ]. Luc-Bouten’s method of control is slightly different. Here, we choose a system observable Z and give an infinitesimal control unitary U c ( t , t ⁢ + ⁢ d t ) ⁢ = ⁢ U ( t ⁢ + ⁢ d t ) * e x p ( - i Z d Y i ( t ) ) U ( t + d t ) ⁢ = ⁢ e x p ( i Z ( t + d t ) d Y o ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1014.tif"/> Note that Y o ( t ) ⁢ = ⁢ U ( t + d t ) * Y i ( t ) U ( t + d t ) , Y o ( t + d t ) ⁢ = ⁢ U ( t + d t ) * Y i ( t + d t ) U ( t + d t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1015.tif"/> and hence, forming the difference, we get d Y 0 ( t ) = U ( t + d t ) * d Y i ( t ) U ( t + d t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1016.tif"/> Also, Z commutes with dY i ( t ) and hence Z ( t + dt )= U ( t + dt )* ZU ( t + dt ) commutes with U ( t + dt )* dY o ( t ) U ( t + dt ). Let ρ c ( t ) denote the controlled state at time t, ie after applying the Belavkin filter and control upto time t to the HP evolved state. Then, we apply the Belavkin filter from t to t + dt , via ρ B ( t + d t ) - ρ c ( t ) + δ ρ B ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1017.tif"/> where δ ρ B ( t ) = L t * ( ρ c ( t ) ) d t + ( M t ρ c ( t ) + ρ c ( t ) M t * - T r ( M t + M t * ) ρ c ( t ) ) ρ c ( t ) ) ( d Y o ( t ) - T r ( ρ c ( t ) ( M t + M t * ) ) d t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1018.tif"/> Finally, the filtered and controlled state at time t + dt is given by ρ c ( t + d t ) - ρ c ( t ) ⁢ + ⁢ δ ρ c ( t ) = https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1019.tif"/>

CRC Press eBooks•Feb 3, 2021•0 citations

Design of quantum gates using electrons, positrons and photons, quantum information theory and quantum stochastic filtering

Harish Parthasarathy

Assume that system 1 evolves according to the dynamics ρ ′ ( t ) ⁢ = i [ H 1 , ρ ( t ) ] + ⁢ θ 1 ( ρ ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math407.tif"/> 78 and system 2 evolves according to σ ′ ( t ) = − i [ H 2 , σ ( t ) ] + ⁢ θ 2 ( σ ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math408.tif"/> where θ 1 , θ 2 are Lindblad maps: θ 1 ( ρ ) ⁢ = ⁢ ( - 1 / 2 ) Σ k = 1 N ( L k * L k ρ + ρ L k * L k - 2 L k ρ L k * ) θ 2 ( σ ) ⁢ = ⁢ ( - 1 / 2 ) ⁢ Σ k = 1 M ⁢ ( P k * P k σ ⁢ + ⁢ σ P k * P k - 2 P k σ P k * ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math409.tif"/> Then calculate the rate of change of the relative entropy between the two states at time t : H ( ρ ( t ) , σ ( t ) ) ⁢ = T r ( ρ ( t ) ( l o g ( ρ ( t ) ) ⁢ − l o g ( σ ( t ) ) ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math410.tif"/> hint: Use the formula for the differential of the exponential map in the form: d d t e x p ( Z ( t ) ) ⁢ = e x p ( Z ( t ) ) g ( a d ( Z ( t ) ) − 1 ( Z ′ ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math411.tif"/> where g ( z ) = z 1 − exp ⁡ ( − z ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math412.tif"/> Equivalently, Z ′ ( t ) ⁢ = g ( α d ( Z ( t ) ) ( e x p ( − Z ( t ) ) . d d t e x p ( Z ( t ) ) ) f https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math413.tif"/> Take ρ ( t )= exp ( Z ( t )) and deduce that d d t l o g ( ρ ( t ) ) ⁢ = g ( a d ( l o g ( ρ ( t ) ) ) ) ( ρ ( t ) − 1 ρ ′ ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math414.tif"/> Now note that g ( z ) ⁢ = 1 1 − z / 2 ! + z 2 / 3 ! + … = - 1 + Σ k ≥ 1 c [ k ] z k https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math415.tif"/> Note that ρ(t) has eigenvalues in the range [0,1]. So log(ρ(t)) has all eigenvalues in the range (−∞,0]. Hence, ad(log(ρ(t)) has eigenvalues in the range ℝ = (−∞,∞). So in order that we be able to substitute ad(log(ρ(t))) for z in the above equation, we require that the Taylor series for g(z) be convergent for all z ∈ ℝ which may not be the case. Note that Tr(ρ ' (t)) = 0 and that for k ≥ 1, T r ( ρ ( t ) a d ( l o g ( ρ ( t ) ) k ( ρ ( t ) − 1 ρ ′ ( t ) ) ) ⁢ - T r ( a d ( l o g ( ρ ( t ) ) ) k ( ρ ′ ( t ) ) ) ⁢ = 0 https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math416.tif"/> since the commutator of two operators has trace zero (provided that the operators are bounded). Thus, we get T r ( ρ ( t ) d d t l o g ( ρ ( t ) ) ) ⁢ = 0 https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math417.tif"/> 79 and so d d t T r ( ρ ( t ) . ⁢ l o g ( ρ ( t ) ) ) ⁢ = T r ( ρ ′ ( t ) . ⁢ l o g ( ρ ( t ) ) ) - T r ( T ( ρ ( t ) ) . ⁢ l o g ( ρ ( t ) ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math418.tif"/> where T ( X ) ⁢ = - i [ H 1 , X ] + θ 1 ( X ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math419.tif"/> Now, T r ( [ H 1 , ρ ] . ⁢ l o g ( ρ ) ) ⁢ = T r ( [ H 1 ⁢ . ⁢ l o g ( ρ ) , ρ ] ) ⁢ = 0 https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math420.tif"/> so d d t ⁢ T r ( ρ ( t ) ⁢ . ⁢ l o g ( ρ ( t ) ) ) ⁢ = T r ( θ 1 ( ρ ( t ) ) ⁢ . ⁢ l o g ( ρ ( t ) ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math421.tif"/>

CRC Press eBooks•Feb 3, 2021•0 citations
A classical liquid may be described by specifying the density ρ and the current vector j , which are defined in the following manner: ( 3-1 ) ρ = ∑ m α δ ( R − r α ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9780429502897/20d72e0b-9d39-4175-acc6-a943547edf8a/content/eq51.tif"/> ( 3-2 ) j = ∑ p α δ ( R − r α ) = ∑ m α v α δ ( R − r α ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9780429502897/20d72e0b-9d39-4175-acc6-a943547edf8a/content/eq52.tif"/>
The quantum revolution in the beginning of the 20th century changed our concepts of dynamics. The dynamic variables such as position r $ {\mathbf{r}} $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_1.tif"/> and momentum p $ {\mathbf{p}} $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_2.tif"/> are replaced by the corresponding position and momentum operators r ^ $ {\hat{\mathbf{r}}} $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_3.tif"/> and p ^ $ {\hat{\mathbf{p}}} $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_4.tif"/> , respectively When one would like to observe a dynamical variable, say ω , of a particle, one has to introduce a wavefunction ψ , which contains all the dynamical information of particles, and to apply the corresponding dynamical operator Ω ^ $ \hat{\Omega } $ https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/inline-math1_5.tif"/> to the wavefunction ψ . Then one may obtain the desired dynamical value of the particle as the eigenvalue of the operator, 1.1 Ω ^ ψ = ω ψ . Ω^ψ=ωψ. \hat{\Omega }\psi = \omega \psi . Ω^ψ=ωψ. https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781315196596/a397f31e-8107-49bb-8694-b978d38168d6/content/math1_1.tif"/>
We wish j t ( X ) to track the noiseless trajectory X d ( t )∈ L ( h ). Assume non-demoltion measurements Y o ( t )= U ( t )* Y i ( t ) U ( t ) are made and the Belavkin filter for π t ( X ) =  j t ( X ) | η t ) , η t = σ ( Y 0 ( s ) : s ≤ t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1006.tif"/> has been constructed as d π t ( X ) = F t ( X ) d t + G t ( X ) d Y o ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1007.tif"/> trajectory estimation error X d ( t )− π t ( X ) is given as feedback to the state equations in the form d j t ( X ) = j t ( θ 0 ( X ) ) d t ⁢ + ⁢ j t ( θ 1 ( X ) ) d A ( t ) + j t ( θ 2 ( X ) ) d A ( t ) * + ⁢ j t ( θ 3 ( X ) ) d Λ ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1008.tif"/> 192 + ⁢ K ( t ) ( X d ( t ) − π t ( X ) ) d t https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1009.tif"/> Alternately, if we assume that X d follows a noiseless trajectory, ie, it evolves according to the equation X d ( t ) , = j t ( 0 ) ( X d ) = U 0 ( t ) * X d U 0 ( t ) , U 0 ( t ) = e x p ( − i t H 0 ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1010.tif"/> then we can give as feedback the estimation error E ( t ) = X d ( t ) − T r 2 ( π t ( X ) ( I ⊗ | ϕ ( u ) > < ϕ ( u ) | ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1011.tif"/> which is a system observable. This error feedback is given to the above state equations, ie, to the Evans-Hudson flow resulting in the dynamics d j t ( X ) = j t ( θ 0 ( X ) ) d t + j t ( θ 1 ( X ) ) d A ( t ) + j t ( θ 2 ( X ) ) d A ( t ) * + j t ( θ 3 ( X ) ) d Λ ( t ) + K ( t ) E ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1012.tif"/> Alternately, we can incorporate this error feedback in the original HP equation as d U ( t ) ⁢ = ⁢ ( - ( i H + ⁢ P + ⁢ K ( t ) E ( t ) ) d t + L 1 d A ( t ) ⁢ + ⁢ L 2 d A ( t ) * ⁢ + ⁢ S d Λ ( t ) ) U ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1013.tif"/> resulting in the Evans-Hudson flow the same as above but with θ 0 (X) replaced by θ 0 ( X )+ iK ( t )[ E ( t ), X ]. Luc-Bouten’s method of control is slightly different. Here, we choose a system observable Z and give an infinitesimal control unitary U c ( t , t ⁢ + ⁢ d t ) ⁢ = ⁢ U ( t ⁢ + ⁢ d t ) * e x p ( - i Z d Y i ( t ) ) U ( t + d t ) ⁢ = ⁢ e x p ( i Z ( t + d t ) d Y o ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1014.tif"/> Note that Y o ( t ) ⁢ = ⁢ U ( t + d t ) * Y i ( t ) U ( t + d t ) , Y o ( t + d t ) ⁢ = ⁢ U ( t + d t ) * Y i ( t + d t ) U ( t + d t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1015.tif"/> and hence, forming the difference, we get d Y 0 ( t ) = U ( t + d t ) * d Y i ( t ) U ( t + d t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1016.tif"/> Also, Z commutes with dY i ( t ) and hence Z ( t + dt )= U ( t + dt )* ZU ( t + dt ) commutes with U ( t + dt )* dY o ( t ) U ( t + dt ). Let ρ c ( t ) denote the controlled state at time t, ie after applying the Belavkin filter and control upto time t to the HP evolved state. Then, we apply the Belavkin filter from t to t + dt , via ρ B ( t + d t ) - ρ c ( t ) + δ ρ B ( t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1017.tif"/> where δ ρ B ( t ) = L t * ( ρ c ( t ) ) d t + ( M t ρ c ( t ) + ρ c ( t ) M t * - T r ( M t + M t * ) ρ c ( t ) ) ρ c ( t ) ) ( d Y o ( t ) - T r ( ρ c ( t ) ( M t + M t * ) ) d t ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1018.tif"/> Finally, the filtered and controlled state at time t + dt is given by ρ c ( t + d t ) - ρ c ( t ) ⁢ + ⁢ δ ρ c ( t ) = https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math1019.tif"/>
Assume that system 1 evolves according to the dynamics ρ ′ ( t ) ⁢ = i [ H 1 , ρ ( t ) ] + ⁢ θ 1 ( ρ ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math407.tif"/> 78 and system 2 evolves according to σ ′ ( t ) = − i [ H 2 , σ ( t ) ] + ⁢ θ 2 ( σ ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math408.tif"/> where θ 1 , θ 2 are Lindblad maps: θ 1 ( ρ ) ⁢ = ⁢ ( - 1 / 2 ) Σ k = 1 N ( L k * L k ρ + ρ L k * L k - 2 L k ρ L k * ) θ 2 ( σ ) ⁢ = ⁢ ( - 1 / 2 ) ⁢ Σ k = 1 M ⁢ ( P k * P k σ ⁢ + ⁢ σ P k * P k - 2 P k σ P k * ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math409.tif"/> Then calculate the rate of change of the relative entropy between the two states at time t : H ( ρ ( t ) , σ ( t ) ) ⁢ = T r ( ρ ( t ) ( l o g ( ρ ( t ) ) ⁢ − l o g ( σ ( t ) ) ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math410.tif"/> hint: Use the formula for the differential of the exponential map in the form: d d t e x p ( Z ( t ) ) ⁢ = e x p ( Z ( t ) ) g ( a d ( Z ( t ) ) − 1 ( Z ′ ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math411.tif"/> where g ( z ) = z 1 − exp ⁡ ( − z ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math412.tif"/> Equivalently, Z ′ ( t ) ⁢ = g ( α d ( Z ( t ) ) ( e x p ( − Z ( t ) ) . d d t e x p ( Z ( t ) ) ) f https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math413.tif"/> Take ρ ( t )= exp ( Z ( t )) and deduce that d d t l o g ( ρ ( t ) ) ⁢ = g ( a d ( l o g ( ρ ( t ) ) ) ) ( ρ ( t ) − 1 ρ ′ ( t ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math414.tif"/> Now note that g ( z ) ⁢ = 1 1 − z / 2 ! + z 2 / 3 ! + … = - 1 + Σ k ≥ 1 c [ k ] z k https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math415.tif"/> Note that ρ(t) has eigenvalues in the range [0,1]. So log(ρ(t)) has all eigenvalues in the range (−∞,0]. Hence, ad(log(ρ(t)) has eigenvalues in the range ℝ = (−∞,∞). So in order that we be able to substitute ad(log(ρ(t))) for z in the above equation, we require that the Taylor series for g(z) be convergent for all z ∈ ℝ which may not be the case. Note that Tr(ρ ' (t)) = 0 and that for k ≥ 1, T r ( ρ ( t ) a d ( l o g ( ρ ( t ) ) k ( ρ ( t ) − 1 ρ ′ ( t ) ) ) ⁢ - T r ( a d ( l o g ( ρ ( t ) ) ) k ( ρ ′ ( t ) ) ) ⁢ = 0 https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math416.tif"/> since the commutator of two operators has trace zero (provided that the operators are bounded). Thus, we get T r ( ρ ( t ) d d t l o g ( ρ ( t ) ) ) ⁢ = 0 https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math417.tif"/> 79 and so d d t T r ( ρ ( t ) . ⁢ l o g ( ρ ( t ) ) ) ⁢ = T r ( ρ ′ ( t ) . ⁢ l o g ( ρ ( t ) ) ) - T r ( T ( ρ ( t ) ) . ⁢ l o g ( ρ ( t ) ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math418.tif"/> where T ( X ) ⁢ = - i [ H 1 , X ] + θ 1 ( X ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math419.tif"/> Now, T r ( [ H 1 , ρ ] . ⁢ l o g ( ρ ) ) ⁢ = T r ( [ H 1 ⁢ . ⁢ l o g ( ρ ) , ρ ] ) ⁢ = 0 https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math420.tif"/> so d d t ⁢ T r ( ρ ( t ) ⁢ . ⁢ l o g ( ρ ( t ) ) ) ⁢ = T r ( θ 1 ( ρ ( t ) ) ⁢ . ⁢ l o g ( ρ ( t ) ) ) https://s3-euw1-ap-pe-df-pch-content-public-p.s3.eu-west-1.amazonaws.com/9781003163626/7a1b5c8e-5b02-4966-b7fd-03fb766469c4/content/math421.tif"/>