the parametric beating of a weak probe field with a prepared raman coherence in a far - off - resonance medium has been extensively studied @xcite .
it has been demonstrated that multimode laser radiation @xcite and incoherent fluorescent light @xcite can be replicated into raman sidebands . since a substantial molecular coherence
can be produced by the two - color adiabatic raman pumping method @xcite , the quantum conversion efficiency of the parametric beating technique can be maintained high even for weak light with less than one photon per wave packet @xcite . to describe the statistical properties of a weak quantum probe and its first - order stokes and anti - stokes sidebands in the parametric beating process , a simplified quantum treatment
has recently been performed @xcite .
it has been shown that the statistical properties of the quantum probe can be replicated into the two sidebands nearest to the input line , in agreement with the experimental observations @xcite .    however , many experiments have reported the observations of ultrabroad raman spectra with a large number of sidebands @xcite . in the experiments with solid hydrogen @xcite , at least two anti - stokes sidebands and two stokes sidebands have been observed . in the experiment with molecular deuterium @xcite ,
a large raman coherence @xmath0 and about 20 raman sidebands , covering a wide spectral range from near infrared through vacuum ultraviolet , have been generated . in rare - earth
doped dielectrics with low raman frequency and long - lived spin coherence , a substantial raman coherence @xmath1 and an extremely large number of sidebands ( about @xmath2 ) can also be generated @xcite .
broad combs of raman sidebands @xcite have been intensively studied because they may synthesize to subfemtosecond @xcite and subcycle @xcite pulses .
the generation of broad combs of raman sidebands has always been examined as a semiclassical problem .
while classical treatments are sufficient for many purposes , a quantum treatment is required when the statistical properties of the radiation fields are important . on the other hand , broad combs of raman sidebands with similar nonclassical properties and different frequencies
may find useful applications for high - performance optical communication .
therefore , it is intriguing to examine the quantum aspects of high - order coherent raman processes .    in this paper
, we extend the treatment of ref .
@xcite to study various quantum properties of multiorder sidebands generated by the beating of a quantum probe field with a prepared raman coherence in a far - off - resonance medium . under the conditions of negligible dispersion and limited bandwidth
, we derive a bessel - function solution for the sideband field operators .
we analytically and numerically calculate various quantum statistical characteristics of the sideband fields generated from a single - mode quantum input .
we show that , with increasing the effective medium length or the raman sideband order , the autocorrelation functions , cross - correlation functions , photon - number distributions , and squeezing factors undergo oscillations governed by the bessel functions
. meanwhile , the normalized autocorrelation functions and normalized squeezing factors of the single - mode probe field are not altered and can be replicated into a broad comb of generated multiorder raman sidebands .
we study the mixing and modulation of photon statistical properties in the case of two - mode input .
we show that the prepared raman coherence and the medium length can be used as control parameters to switch a sideband field from one type of photon statistics to another type , or from a non - squeezed state to a squeezed state and vice versa .
we also discuss two - photon interference in coherent raman scattering .
although the multiorder coherent raman scattering can produce a broad comb of sideband fields with different frequencies , it behaves in many aspects as a beam splitter @xcite .
therefore , in this paper , we also make comparison of this conventional device with our system as and when it is possible .    before we proceed , we note that , in related problems , the generation of correlated photons using the @xmath3 and @xmath4 parametric processes has been studied @xcite .
the correlations between the stokes and anti - stokes sidebands and the possibility of transferring a quantum state of light from one carrier frequency to another carrier frequency ( multiplexing ) have been discussed for resonant systems @xcite .
the paper is organized as follows . in sec .
[ sec : model ] , we describe the model and present the basic equations . in sec .
[ sec : single - mode ] , we study various quantum characteristics of the sideband fields generated from a single - mode quantum input . in sec .
[ sec : two - mode ] , we discuss the quantum properties of the sideband fields generated from a two - mode quantum input . finally , we present the conclusions in sec .  [
sec : summary ] .
we consider a far - off - resonance raman medium shown schematically in fig .  [ fig1 ] .
level @xmath5 with energy @xmath6 is coupled to level @xmath7 with energy @xmath8 by a raman transition via intermediate levels that are not shown in the figure .
we send a pair of long , strong , classical laser fields , with carrier frequencies @xmath9 and @xmath10 , and a short , weak , quantum probe field @xmath11 , with one or several carrier frequencies , through the raman medium , along the @xmath12 direction .
the timing and alignment of these fields are such that they substantially overlap with each other during the interaction process .
the driving laser fields are tuned close to the raman transition @xmath13 , with a small finite two - photon detuning @xmath14 , but are far detuned from the upper electronic states @xmath15 of the molecules .
we assume that all the frequency components of the input probe field are separated by integer multiples of the raman modulation frequency @xmath16 .
the driving fields adiabatically produce a raman coherence @xmath17 @xcite . when the probe field propagates through the medium , it beats with the prepared raman coherence .
since the probe field is weak and short compared to the driving fields , the medium state and the driving fields do not change substantially during this step .
the beating of the probe field with the prepared raman coherence leads to the generation of new sidebands in the total output field @xmath18 .
the frequencies of the sideband fields @xmath19 are given by @xmath20 , where @xmath21 is integer and @xmath22 is a carrier frequency of the input probe field .
the range of @xmath21 should be appropriate so that @xmath23 is positive .
the probe field is taken to be not too short so that the fourier - transformation limited broadening is negligible .
we assume that the prepared raman coherence @xmath17 is substantial so that the spontaneous raman process is negligible compared to the stimulated and parametric processes .
consequently , the quantum noise can be neglected . unlike ref .
@xcite , our model does not require any restriction on the magnitude of the coherence as all raman sidebands are included .
when we take the propagation equation for the classical raman sidebands @xcite and replace the field amplitudes by the quantum operators , we obtain @xmath24 here , @xmath25 and @xmath26 are the dispersion and coupling constants , respectively .
we have denoted @xmath27 , where @xmath28 is the molecular number density .
we take all the sidebands to be sufficiently far from resonance that the dispersion of the medium is negligible . in this case
, we have @xmath29 and @xmath30 .
we write @xmath31 $ ] , where @xmath32 and @xmath33 , and assume that @xmath34 and @xmath35 are constant in time and space .
we change the variables by @xmath36 $ ] . using photon operators , we can write @xmath37 . here , @xmath38 is the quantization length taken to be equal to the medium length , @xmath39 is the quantization transverse area taken to be equal to the beam area , @xmath40 is a bloch wave vector , and @xmath41 and @xmath42 are the annihilation and creation operators for photons in the spectral mode @xmath21 and the spatial mode @xmath40 . for simplicity , we restrict our discussion to the case where each sideband field contains only a single spatial mode ( with , e.g. , @xmath43 )
. then , eq .  ( [ 1 ] ) yields @xmath44 where @xmath45 . for the medium length @xmath38 , the evolution time is @xmath46 .
it follows from eq .
( [ 4 ] ) that the total photon number is conserved in time .
note that eq .
( [ 4 ] ) represents the heisenberg equation for the fields that are coupled to each other by the effective interaction hamiltonian @xmath47 the interaction between the sideband fields via the prepared raman coherence is analogous to the interaction between the transmitted and reflected fields from a conventional beam splitter @xcite .
however , the two mechanisms are very different in physical nature .
the most important difference between them is that the two fields from the conventional beam splitter have the same frequency while the sideband fields in the raman scheme have different frequencies .
in addition , the model hamiltonian ( [ 8b ] ) involves an infinitely large number of raman sidebands , separated by integer multiples of the raman modulation frequency @xmath48 . despite these differences , the model ( [ 8b ] )
can be called the multiorder raman beam splitter . the parameters @xmath49 determine the transmission and scattering coefficients for the fields at the raman beam splitter .
we assume that the bandwidth of the generated raman spectrum is small compared to the characteristic probe frequency @xmath22 . in this case , the @xmath21-dependence of the coupling parameters @xmath50 can be neglected , that is , we have @xmath51 . with this assumption , we find the following solution to eq .
( [ 4 ] ) : @xmath52 here , @xmath53 is the @xmath54th - order bessel function .
the expression ( [ 8 ] ) for the output field operators is a generalization of the bessel - function solution obtained earlier for the classical fields @xcite .
the number of generated raman sidebands is characterized by the effective interaction time @xmath55 or , equivalently , the effective medium length @xmath56 , where @xmath57 the coefficient @xmath58 characterizes the strength of the parametric coupling and is proportional to the prepared raman coherence @xmath34 , that is , to the intensities of the driving laser fields .
the bessel functions @xmath59 are the transmission ( @xmath60 ) and scattering ( @xmath61 ) coefficients for the raman sidebands , similar to the transmission and reflection coefficients of a conventional beam splitter .
the assumption of limited bandwidth requires @xmath62 , that is , @xmath63 @xcite . in what follows we use the expression ( [ 8 ] ) to calculate various quantum statistical characteristics of the sideband fields , namely , the autocorrelation functions , the two - mode cross - correlation functions , the squeezing factors , and the relation between the @xmath64 representations of the output and input states
in this section , we consider the case where the input probe field has a single carrier frequency @xmath22 . in other words , we assume that the sideband @xmath65 is initially prepared in a quantum state @xmath66 and the other sidebands are initially in the vacuum state .
the density matrix of the initial state of the fields is given by @xmath67      we study the autocorrelations of photons in the generated raman sidebands .
we use eq .
( [ 8 ] ) and apply the initial density matrix ( [ 8a ] ) to calculate the normally ordered moments @xmath68 of the photon - number operators @xmath69 .
the result is @xmath70 in particular , the mean photon numbers of the sidebands are given by @xmath71 here , @xmath72 is the photon - number operator for the input field .
the @xmath73th - order autocorrelation functions of the sidebands are defined by @xmath74 . from eqs .
( [ 9 ] ) and ( [ 9a ] ) , we find @xmath75 where @xmath76 is the @xmath73th - order autocorrelation function of the input field .
( first row ) and second - order autocorrelation function @xmath77 ( second row ) , both scaled to their initial values for the probe field , as functions of the effective medium length @xmath56 ( first column ) and the sideband order @xmath21 ( second column ) . in ( a ) and ( c )
, the sideband order is @xmath78 ( solid line ) , @xmath79 ( dashed line ) , and @xmath80 ( dotted line ) . in ( b ) and ( d ) , the effective medium length is @xmath81 . in ( c ) , we have amplified @xmath82 ( dashed line ) and @xmath83 ( dotted line ) by 5 and 10 times , respectively . in ( b ) and ( d ) , the negative side of the @xmath21 axis is not shown because the functions plotted are symmetric in @xmath21 . ]
equations ( [ 9a ] ) and ( [ 9b ] ) indicate that , when we increase the effective medium length @xmath56 or the sideband order @xmath21 , the mean photon number @xmath84 and the autocorrelation function @xmath85 undergo oscillations as described by even powers of the bessel function @xmath86 . such oscillatory behavior is illustrated in fig .
when the sideband order @xmath21 is higher , the onset of @xmath84 occurs later [ see fig .
[ fig2](a ) ] and , hence , so does the onset of @xmath85 [ see fig .
[ fig2](c ) ] . for a fixed @xmath21 , both @xmath84 and @xmath85 reach their largest values at the same optimal medium length @xmath87 , where @xmath88 is the position of the first peak of @xmath89 .
the higher the sideband order @xmath21 , the larger is the optimal length @xmath90 and the smaller are the maximal output values of @xmath84 and @xmath85 [ see figs .
[ fig2](a ) and [ fig2](c ) ] .
figures [ fig2](b ) and [ fig2](d ) show that @xmath84 and @xmath85 are substantially different from zero only in the region where @xmath91 is not too large compared to @xmath56 . for a given @xmath56 , both @xmath84 and @xmath85 achieve their maximal values at @xmath92 .
the normalized @xmath73th - order autocorrelation functions of the sidebands are defined by @xmath93 .
these functions characterize the overall statistical properties , such as sub - poissonian , poissonian , or super - poissonian photon statistics , regardless of the mean photon number . unlike the mean photon number @xmath84 and the autocorrelation function @xmath85 ,
the normalized autocorrelation function @xmath94 does not oscillate when we change the effective medium length @xmath56 or the sideband order @xmath21 .
indeed , with the help of eq .
( [ 9 ] ) , we find @xmath95 where @xmath96 .
equation ( [ 10 ] ) indicates that the generated sideband fields and the probe field have the same normalized autocorrelation functions , which are independent of the evolution time and are solely determined by the statistical properties of the input field . in other words , the normalized autocorrelation functions of the probe field do not change during the parametric beating process and are precisely replicated into the comb of generated sidebands .
such a replication of the normalized autocorrelation characteristics can be called autocorrelation multiplexing . in particular ,
if the photon statistics of the input field is sub - poissonian , poissonian , or super - poissonian , the photon statistics of each of the sideband fields will also be sub - poissonian , poissonian , or super - poissonian , respectively .
this result is in agreement with the experiments on replication of multimode laser radiation @xcite and broadband incoherent light @xcite .
the ability of the raman medium to multiplex the autocorrelation characteristics is similar to the ability of a conventional beam splitter @xcite .
it is not surprising that the normalized autocorrelation functions of the probe field are replicated into the sidebands in the parametric beating process .
such a replication is possible because the medium is far off resonance and the quantum probe field is weak compared to the driving fields . under these two conditions ,
the photon annihilation operators are linearly transformed as described by the linear differential equation ( [ 4 ] ) .
this equation shows that the photon annihilation operators are not mixed up with the creation operators , and therefore the evolution of the annihilation operators is linear with respect to the initial annihilation operators .
we should , however , emphasize here that the replication of all of the normalized autocorrelation functions @xmath97 ( for all orders @xmath73 ) of the input field does not mean the replication of the quantum state @xmath66 .
in fact , the oscillations of the mean photon number @xmath84 and the autocorrelation function @xmath85 indicate that the photon - number distributions and consequently the quantum states of the sidebands evolve in a rather complicated way and are quite different from those for the input field
. a separable state at the input can produce an entangled state @xcite .
in addition , cross - correlations between the sidebands can be generated from initially uncorrelated fields , and a fock state at the input does not produce sideband fields in isolated fock states [ see below ] .
we study the correlations between the generated raman sidebands .
for two different sidebands @xmath54 and @xmath98 ( @xmath99 ) , we have @xmath100 the cross - correlation function for the two sidebands is defined by @xmath101 . using eqs .
( [ 9a ] ) and ( [ 11b ] ) , we find @xmath102 when we extend @xmath103 for @xmath104 , we have @xmath105 . according to eq .
( [ 11a ] ) , the cross - correlation function @xmath103 oscillates when we change the effective medium length @xmath56 or the sideband orders @xmath54 and @xmath98 . such oscillatory behavior is illustrated in fig .
[ fig3 ] .    , scaled with the second - order autocorrelation function of the input field , as a function of ( a ) the effective medium length @xmath56 and ( b ) the sideband order @xmath98 . here
, the sideband order @xmath54 is fixed to @xmath79 . in ( a )
, the sideband order @xmath98 is @xmath78 ( solid line ) , @xmath106 ( dashed line ) , and @xmath80 ( dotted line ) . in ( b ) , the effective medium length is @xmath81 . in ( b ) , the negative side of the @xmath98 axis is not shown because the function plotted is symmetric with respect to the sideband orders . ]
the normalized cross - correlation function is defined by @xmath107 . unlike the function @xmath103
, the normalized function @xmath108 does not oscillate .
indeed , we find the relation @xmath109 as seen from the above relation , the normalized cross - correlation functions @xmath108 for all possible sideband pairs @xmath110 are equal to each other , to the normalized second - order autocorrelation function @xmath111 for each sideband @xmath21 , and to the normalized second - order autocorrelation function @xmath112 of the input field .
when @xmath113 , that is , when the photon statistics of the input field is non - poissonian , we obtain @xmath114 , a signature of cross - correlations between the sidebands .
such correlations are generated although the sidebands are initially not correlated .
in particular , if the input field has a sub - poissonian photon statistics ( @xmath115 ) , anti - correlations between the sidebands ( @xmath116 ) will be generated .
note that the conventional beam splitters also have a similar property @xcite .
the generation of cross - correlations between the sidebands indicates that the quantum states of the generated sidebands are different from that of the input field .
we emphasize that the anti - correlation generation can not be explained by the classical statistics of the fields with positive @xmath64 functions although the sideband dynamics is linear with respect to the field variables .      to get deeper insight into the quantum properties of the generated raman sidebands , we derive the photon - number distributions for the output fields .
the joint photon - number distribution for the output fields is defined by @xmath117 , where @xmath118 is the density matrix of the output state . here
, @xmath119 is the evolution operator . using eqs .
( [ 8 ] ) and ( [ 8a ] ) , we find @xmath120 where @xmath121 is the photon - number distribution of the input field , and @xmath122 is the total photon number .
from eq .
( [ 10a ] ) , the marginal photon - number distribution @xmath123 for the sideband @xmath21 is obtained as @xmath124^k   p_{\mathrm{in}}(n+k ) .
\label{10b}\ ] ] clearly , @xmath123 is in general different from @xmath125 .
thus , despite the replication of the normalized autocorrelation functions , the photon - number distribution of the probe field is not replicated into the sidebands .
we examine several particular cases .
first , we consider the case where the probe field is initially prepared in a coherent state @xmath126 .
this state is characterized by a poisson distribution @xmath127 for the photon number , where @xmath128 . in this case ,
( [ 10a ] ) and ( [ 10b ] ) yield @xmath129 and @xmath130 respectively , where @xmath131 .
thus , the generated sidebands are not correlated , and the marginal photon - number distributions for the individual sidebands remain poisson distributions during the evolution process .
second , we consider the case where the probe field is initially in a fock state @xmath132 . in this case , we have @xmath133 . therefore , eq .  ( [ 10b ] ) yields
@xmath134^{n - n }   \frac{n!}{n!(n - n)!}\ ] ] for @xmath135 , and @xmath136 for @xmath137 .
meanwhile , eq .  ( [ 10a ] ) yields @xmath138 , that is , the joint photon - number distribution is not a product of the marginal photon - number distributions for the individual sidebands .
thus , the sidebands are correlated .
they are not generated in isolated fock states .
finally , we consider the case where the probe field is initially in a thermal state , which is characterized by a boltzmann photon - number distribution @xmath139 . in this case , eq .  ( [ 10b ] ) yields @xmath140 where @xmath131 .
thus , the marginal photon - number distributions for the individual sidebands remain boltzmann distributions during the evolution process .
however , according to eq .
( [ 10a ] ) , we have @xmath138 , a signature of correlations between the generated sidebands .
we examine the squeezing of the field quadratures .
a field quadrature of the @xmath21th mode is defined by @xmath141 .
we say that the @xmath21th mode is in a squeezed state if there exists such a phase @xmath142 that @xmath143 or , equivalently , @xmath144 , where @xmath145 .
the squeezing degree is measured by the quantity @xmath146 .
note that the relation between the squeezing factor @xmath147 and the conventional squeezing parameter @xmath148 is @xmath149 . in terms of the photon operators , we have @xmath150 + [ ( \langle \hat b_q^2\rangle-\langle \hat b_q\rangle^2 )   e^{-2i\varphi}+\mathrm{c.c . } ] .
\label{12c}\ ] ] using eqs .
( [ 8 ] ) and ( [ 8a ] ) , we find @xmath151 and @xmath152 when we insert eqs .
( [ 9a ] ) , ( [ 12a ] ) , and ( [ 12b ] ) into eq .
( [ 12c ] ) , we obtain @xmath153 here , @xmath154 denotes the squeezing factor for the @xmath142-quadrature of the input field . equation ( [ 12 ] ) shows that , if @xmath155 , then @xmath156 .
thus , if the input field is in a squeezed state , then the generated sidebands are also in squeezed states . in other words ,
the squeezing of the input field is transferred to the comb of generated sidebands during the parametric beating process .
the squeezing factor @xmath157 of the sideband @xmath21 is reduced from the input squeezing factor @xmath154 by the factor @xmath158 . unlike the case of linear directional couplers and beam splitters @xcite , the squeezing degree of the probe field
can not be completely transferred to the raman sidebands .
this difference is due to the fact that the linear directional coupler and the beam splitter involve only two output modes while the multiorder coherent raman process involves many more output modes .
note that the phase of the squeezed quadrature of the sideband @xmath21 changes by @xmath159 .
this means that the squeezed quadrature of a generated even - order ( odd - order ) raman sideband is parallel ( orthogonal ) to that of the input field .
we introduce the normalized squeezing factor @xmath160 .
we find the relation @xmath161 where @xmath162 is the normalized squeezing factor for the input field .
equation ( [ 12d ] ) implies that , besides a shift of the quadrature phase angle , the normalized squeezing factor for the input field is replicated into the comb of generated sidebands .
this result can be used to convert squeezing to a new frequency , i.e. , to perform squeezing multiplexing .
the relation @xmath163 indicates that the @xmath56- and @xmath21-dependences of the squeezing factor @xmath157 are similar to those of the mean photon number @xmath164 [ see figs .
[ fig2](a ) and [ fig2](b ) ] .
note that , when the input field is in a coherent state , the sideband fields have no squeezing .
this property is similar to the case of four - wave mixing but is unlike the case of degenerate parametric down - conversion , where perfect squeezing can in principle be obtained . since squeezed states are nonclassical states , the ability of the raman medium to multiplex squeezing from a probe field to its sidebands is a quantum property that can not be described by the classical statistics of the fields with positive @xmath64 functions .
we calculate the quantum state of the output fields for several cases .
first , we consider the case where the input sideband @xmath78 is initially in a coherent state @xmath126 .
the state of the fields at the input is @xmath165 the state of the fields at the output is given by @xmath166 . since @xmath167 , we have @xmath168 . using eq .
( [ 8 ] ) , we find @xmath169 here , @xmath170 is a coherent state of the @xmath21th mode , with the amplitude @xmath171 thus , a probe field in a coherent state can produce sideband fields that are also in coherent states .
such a process can be called coherent - state multiplexing .
this property of the raman medium is similar to the case of conventional beam splitters @xcite .
second , we consider the case where the input sideband @xmath78 is initially prepared in a fock state @xmath132 .
the state of the fields at the input is written as @xmath172 the output state of the fields is given by @xmath173 . with the help of eq .
( [ 8 ] ) , we find @xmath174 here , @xmath175 for @xmath176 , and @xmath177 for @xmath178 .
when @xmath179 and @xmath180 , the output state ( [ 15a ] ) is , in general , a multipartite inseparable ( entangled ) state .    in a particular case where the input state of the probe field is a single - photon state ,
i.e. , @xmath181 , eqs .
( [ 15a ] ) and ( [ 15b ] ) yield @xmath182 here , @xmath183 is the quantum state of a single photon in the sideband @xmath21 with no photons in the other sidebands . in this case , the entanglement between two different sidebands @xmath54 and @xmath98 at the output can be measured by the bipartite concurrence @xmath184 , see @xcite .
finally , we consider the case where the input sideband @xmath78 is initially in an incoherent mixed state @xmath185 with the help of eq .
( [ 15a ] ) , the density matrix of the output state of the fields is found to be @xmath186 to obtain the reduced density matrix @xmath187 for an arbitrary sideband @xmath21 , we trace the total density matrix ( [ 15d ] ) over all sidebands except for the sideband @xmath21 .
then , we find @xmath188 where the marginal photon - number distribution @xmath123 is given by eq .
( [ 10b ] ) .
as seen , the reduced state of each sideband is also an incoherent superposition of fock states .
however , if the input state is different from the vacuum state , then , for @xmath180 , we have @xmath189 , a signature of correlations between the generated sidebands .
moreover , the total density matrix ( [ 15d ] ) of the output fields contains nonzero off - diagonal matrix elements in the fock - state basis .    in a particular case where the initial state of the probe field is a thermal state ,
i.e. , @xmath190(|n\rangle\langle n|)_0 $ ] , the reduced state of each generated sideband is also a thermal state , namely , @xmath191 here , @xmath131 is the mean photon number for the sideband @xmath21 .
the reduced thermal states of the generated sidebands are , however , not isolated from each other .
a far - off - resonance medium with a substantial raman coherence , prepared by two strong driving fields , can efficiently mix and modulate the quantum statistical properties of the sideband fields . to understand this mechanism , we study the case where the input probe field has two carrier frequencies , @xmath22 and @xmath192 , separated by an integer multiple @xmath193 of the raman modulation frequency @xmath48 .
we assume that the raman sidebands @xmath78 and @xmath193 are initially in independent quantum states @xmath66 and @xmath194 , respectively , while the other sidebands are initially in the vacuum state . the density matrix of the initial state of the fields is given by @xmath195 here , @xmath196 .
we study the mixing and modulation of photon statistics of the sideband fields .
when we use eq .
( [ 8 ] ) to calculate the mean photon numbers of the sidebands generated from the initial state ( [ 18 ] ) , we find @xmath197 . \label{19}\end{aligned}\ ] ] furthermore , we find @xmath198 .
\qquad \label{20b}\end{aligned}\ ] ] hence , the second - order autocorrelation function @xmath199 of the sideband @xmath21 is found to be @xmath200 , \nonumber\\ \label{20a}\end{aligned}\ ] ] where @xmath201 , \nonumber\\ \delta_3&=&\langle\hat b_0^\dagger(0)\rangle [ \langle\hat b_\nu^\dagger(0)\hat b_\nu^2(0)\rangle -\langle\hat b_\nu^\dagger(0)\hat b_\nu(0)\rangle\langle \hat b_\nu(0)\rangle].\end{aligned}\ ] ] the first two terms on the right - hand sides of eqs .
( [ 19 ] ) , ( [ 20b ] ) , and ( [ 20a ] ) are the contributions of the individual input sidebands @xmath78 and @xmath193
. the other terms result from the interference between the two interaction channels .     as a function of ( a ) the effective medium length @xmath56 and ( b ) the sideband order @xmath21 in the case of two - mode input .
the input sideband @xmath78 is initially prepared in a fock state , with @xmath80 photons .
the input sideband @xmath79 is initially prepared in a thermal state , with 1 photon in average . in ( a )
, the sideband order is @xmath202 ( solid line ) and @xmath106 ( dashed line ) . in ( b ) , the effective medium length is @xmath203 . ]    unlike the case of single - mode input , in the case of two - mode input , the normalized second - order autocorrelation function @xmath204 depends , in general , on @xmath56 and @xmath21 .
such behavior is illustrated in fig .
[ fig4 ] . when @xmath56 is such that @xmath205 or @xmath206 , we have @xmath207 or @xmath208 , respectively .
consequently , if the two input sidebands have different normalized autocorrelation functions , i.e. , @xmath209 , then , with increasing @xmath56 or @xmath21 , the normalized autocorrelation function @xmath210 will oscillate between the values @xmath211 and @xmath212 [ see fig .
[ fig4 ] ] .
in particular , if the photon statistics of one of the input fields , e.g. , the sideband @xmath78 , is sub - poissonian [ @xmath213 and that of the other input field is super - poissonian [ @xmath214 , then each generated sideband @xmath21 will have complex statistical properties and will oscillate between sub - poissonian [ @xmath215 and super - poissonian [ @xmath216 photon statistics [ see fig .
[ fig4 ] ] .
using the prepared raman coherence @xmath34 or the medium length @xmath38 as a control parameter , we can switch a sideband field from super - poissonian photon statistics to sub - poissonian or vice versa .
similar modulation of photon statistics has been demonstrated in a linear directional coupler @xcite .
we study the mixing and modulation of the squeezing properties of the sideband fields .
when we use eq .
( [ 8 ] ) to calculate the amplitudes @xmath217 and @xmath218 of the sidebands generated from the initial state ( [ 18 ] ) , we find the expressions @xmath219 and @xmath220 we insert eqs .
( [ 19 ] ) , ( [ 20 ] ) , and ( [ 21 ] ) into eq .
( [ 12c ] ) . then , we obtain the squeezing factor @xmath221 where @xmath222 and @xmath223 are the initial squeezing factors of the sidebands @xmath78 and @xmath193 , respectively . as seen , the squeezing factor @xmath147 of the sideband @xmath21 is a superposition of the input squeezing factors @xmath222 and @xmath223 , taken with the quadrature phase shifts @xmath224 and @xmath225 , respectively , and weighted by the factors @xmath158 and @xmath226 , respectively .     as a function of ( a ) the effective medium length @xmath56 and ( b ) the sideband order @xmath21 in the case of two - mode input .
the input sideband @xmath78 is initially prepared in a squeezed vacuum state , with the squeezing parameter @xmath227 and the phase @xmath228 .
the input sideband @xmath79 is initially prepared in a thermal state , with the mean photon number 1 . in ( a ) ,
the sideband order is @xmath202 ( solid line ) and @xmath106 ( dashed line ) . in ( b ) , the effective medium length is @xmath203 . ]    unlike the case of single - mode input , in the case of two - mode input , the normalized squeezing factor @xmath229 varies , in general , with @xmath56 and @xmath21 .
such behavior is illustrated in fig .
[ fig5 ] . when @xmath56 is such that @xmath205 or @xmath206 , we have @xmath230 or @xmath231 , respectively .
consequently , if the normalized squeezing factors @xmath232 and @xmath233 of the two input fields are different , the normalized squeezing factor @xmath234 will oscillate between the values @xmath232 and @xmath233 . in particular , if one of the two input fields , e.g. , the sideband @xmath78 , is squeezed [ @xmath235 and the other input field is not squeezed [ @xmath236 , then , each generated sideband @xmath21 will have complex squeezing properties and will oscillate between a squeezed state [ @xmath237 and a non - squeezed state [ @xmath238 , see fig .
using the prepared raman coherence @xmath34 or the medium length @xmath38 as a control parameter , we can switch a sideband field from a non - squeezed state to a squeezed state or vice versa .
note that a similar result has been obtained for a linear directional coupler @xcite .     and
( b ) mean photon number @xmath239 as functions of the effective medium length @xmath56 in the case where the sideband @xmath78 is initially prepared in a weak squeezed vacuum state and the sideband @xmath79 is initially prepared in a strong coherent state .
the parameters for the initial states of the input sidebands are @xmath227 , @xmath228 , and @xmath240 .
the curves are calculated for the sidebands @xmath78 ( solid line ) , @xmath79 ( dashed line ) , and @xmath241 ( dotted line ) . ]
we analyze a particular case where the sideband @xmath78 is initially in a squeezed vacuum state @xmath242|0\rangle_0 $ ] and the sideband @xmath193 is initially in a coherent state @xmath243 . here
, @xmath244 is a complex number , the modulus @xmath245 characterizes the amount of squeezing , and the phase angle @xmath246 characterizes the alignment of the squeezed vacuum state in phase space . for the input sideband @xmath78
, we have @xcite @xmath247 , @xmath248 , and @xmath249 . for the input sideband @xmath193
, we have @xmath250 , @xmath251 , and @xmath252 .
then , we find from eq .
( [ 19 ] ) that the mean photon number of an arbitrary sideband @xmath21 is @xmath253 we find from eq .
( [ 22 ] ) that the maximal squeezing of the sideband @xmath21 occurs in the @xmath254-quadrature where @xmath255 .
the corresponding value of the squeezing factor is @xmath256 as seen from eq .
( [ 24 ] ) , squeezing can be transferred from the initial squeezed vacuum state of the sideband @xmath78 to the other sidebands .
the squeezing factors of the sidebands are independent of the amplitude @xmath257 of the initial coherent state of the sideband @xmath193 . meanwhile , the mean photon number of each sideband is governed not only by the squeezing parameter @xmath258 of the initial state of the sideband @xmath78 but also by the amplitude @xmath257 of the initial state of the sideband @xmath193 . using this fact
, we can manipulate to get optimized mean photon numbers and squeezing degrees of the sideband fields at the output as per requirement .
in particular , we can convert squeezing from a weak field to a much stronger field . to illustrate this possibility , we plot in fig .  [ fig6 ] the squeezing factor @xmath259 and the mean photon number @xmath239 as functions of the effective medium length @xmath56 for the parameters @xmath227 , @xmath228 , @xmath260 , and @xmath261 . in this case ,
the most negative value of the input squeezing factor @xmath262 is achieved at @xmath263 and is given by @xmath264 , indicating the squeezing degree 86% .
the mean photon number of the input squeezed vacuum state is @xmath265 , rather small .
the solid lines in fig .
[ fig6 ] show that the sideband @xmath78 , initially prepared in a weak squeezed vacuum state , can be significantly enhanced while keeping its squeezing degree substantial . meanwhile , the dashed lines show that , for @xmath266 , the sideband @xmath79 , initially prepared in a strong coherent state , is squeezed by about 29% and has the mean photon number of about 41 .
similarly , the dotted lines show that , for @xmath267 , a generated new sideband @xmath241 is squeezed by about 16% and has the mean photon number of about 39 .
thus , from a weak squeezed field at the input , we can obtain other output squeezed fields that have smaller but still substantial squeezing degrees , much larger mean photon numbers , and different frequencies .
we show the possibility of quantum interference between the probability amplitudes for a pair of photons with different frequencies in the coherent raman process .
we assume that the sidebands @xmath78 and @xmath79 are initially prepared in independent single - photon states .
this initial condition corresponds to the situation where two photons with different frequencies @xmath22 and @xmath268 are incident into the raman medium .
the input state of the fields can be written as @xmath269 the output state of the fields is given by @xmath270 . with the help of eq .
( [ 8 ] ) , we find @xmath271|1_k 1_l\rangle .
\label{16a}\end{aligned}\ ] ] here , the fock state @xmath272 is the state of two photons in the sideband @xmath21 with no photons in the other sidebands , and the fock state @xmath273 is the state in which there is one photon in each of the sidebands @xmath54 and @xmath98 but no photons in the other sidebands .
it follows from eq .
( [ 16a ] ) that the probability for finding two photons in the sideband @xmath21 is @xmath274 the joint probability for finding one photon in each of the sidebands @xmath54 and @xmath98 ( @xmath99 ) is given by @xmath275 ^ 2 .
\label{16d}\ ] ] the probability @xmath276 for having one and only one photon in the sideband @xmath21 is @xmath277 the mean photon number of the sideband @xmath21 is @xmath278 we find the relations @xmath279 , @xmath280 , and @xmath281 , which reflect the symmetry of the generated stokes and anti - stokes sidebands with respect to the two input sidebands 0 and 1 .     for finding one photon in each of the sidebands @xmath78 and @xmath79 as a function of the effective medium length @xmath56 . ]    when we insert @xmath60 and @xmath282 into eq .
( [ 16d ] ) , we obtain the following expression for the joint probability for finding one photon in each of the sidebands @xmath78 and @xmath79 : @xmath283 ^ 2.\end{aligned}\ ] ] this expression shows that the joint probability @xmath284 may become zero at certain values of @xmath56 [ see fig .
[ fig7 ] ] .
this is a signature of destructive interference between two channels that form the state @xmath285 . in the first channel ,
each of the two photons individually transmits through the medium without any changes .
the two - photon probability amplitude for this channel is @xmath286 . in the second channel , both the photons
are scattered from the prepared raman coherence and exchange their sidebands .
the two - photon probability amplitude for this channel is @xmath287 .
since raman scattering produces a phase shift of @xmath288 for each photon , the probability amplitudes for the two channels ( the transmission and scattering of both the photons ) are @xmath289 out of phase .
the interference between the two channels is therefore destructive , yielding the output state @xmath285 with the joint probability @xmath284 given above . when the medium length @xmath38 is such that @xmath290 , the interference between the two two - photon amplitudes becomes completely destructive , and therefore the state @xmath285 is removed from the output state ( [ 16a ] ) .
we denote such a medium length by @xmath291 .
the positions of the zeros of @xmath284 depicted in fig .
[ fig7 ] indicate that the first three values of @xmath291 are given by @xmath292 , 3.11 , and 4.68 .
it is interesting to note that @xmath293 can be determined in an experiment using a single - mode input .
indeed , in the case where a single sideband @xmath78 is initially excited , the mean photon numbers of the generated sidebands are given by eq .
( [ 9a ] ) .
therefore , the effective medium length @xmath293 corresponds to the situation where the probe sideband @xmath78 and its adjacent sidebands @xmath294 have the same mean photon numbers at the output .    there exist literatures on two - photon interference in various systems @xcite .
two - photon interference in coherent raman scattering , described above , is an analogy of two - photon interference at a conventional beam splitter @xcite .
we emphasize that two - photon interference in coherent raman scattering involves copropagating photons with different frequencies in a collinear scheme .
to be more general , we consider the case where an arbitrary number of sidebands is initially excited .
we find that an arbitrary multimode coherent state @xmath295 of the input fields produces a coherent state @xmath296 of the output fields . here
, the output amplitudes @xmath297 are linearly transformed from the input amplitudes @xmath298 as given by @xmath299 consequently , the diagonal coherent - state representation @xmath300 of an arbitrary input quantum state @xmath301 determines the representation @xmath302 of the output state @xmath303 via the equation @xmath304 here , we have introduced the notation @xmath305    if the input state @xmath301 is a classical state @xcite , @xmath300 must be non - negative and less singular than a @xmath14 function , and consequently so must @xmath302 . in this case , the output state @xmath303 is also a classical state .
moreover , since the multimode coherent state @xmath306 is separable and the weight factor @xmath302 is non - negative , the output state @xmath303 is , by definition , separable @xcite .
therefore , a necessary condition for the output fields to be in an inseparable ( entangled ) state or , more generally , in a nonclassical state is that the input field state is a nonclassical state .
a similar condition has been derived for the beam splitter entangler @xcite .
note that , in the case where we use a single - mode input field @xmath65 , prepared in an arbitrary quantum state with the coherent - state representation @xmath307 ( the stokes and anti - stokes sideband fields are initially in the vacuum state ) , eq .
( [ 17 ] ) becomes @xmath308 it has been shown in ref .
@xcite that , when the input field is prepared in an even or odd coherent state , a multipartite entangled coherent state can be generated .
we have studied the quantum properties of multiorder sidebands generated by the beating of a quantum probe field with a prepared raman coherence in a far - off - resonance medium . under the conditions of negligible dispersion and limited bandwidth
, we have derived a bessel - function solution for the sideband field operators .
we have analytically and numerically calculated various quantum statistical characteristics of the multiorder sideband fields .
we have examined the quantum properties of the sideband fields in the case of single - mode quantum input .
we have shown that , when we change the effective medium length or the raman sideband order , the autocorrelation functions , the cross - correlation functions , the photon - number distributions , and the squeezing factors undergo oscillations governed by the bessel functions . when the sideband order is higher , the onset of the sideband generation occurs later and , therefore , so does the onset of the sideband autocorrelation functions . the mean photon number and
the autocorrelation functions of each sideband reach their largest values at the same optimal medium length determined by the first peak of the corresponding bessel function .
the higher the sideband order , the larger is the optimal length and the smaller is the maximal output values of the mean photon number and sideband autocorrelation functions .
meanwhile , the normalized autocorrelation functions and normalized squeezing factors of the probe field are not altered by the parametric beating process .
they are replicated into the comb of generated multiorder sidebands . as the result ,
the single - mode input field and the generated sidebands have identical normalized autocorrelation functions and identical normalized squeezing factors .
in other words , they have similar quantum statistical properties  the same type of photon statistics and the same type of squeezing .
in addition to this resemblance , it has been shown that an input field in a coherent state can produce sideband fields in coherent states .
it has also been shown that , when the input field is prepared in a thermal state , the reduced state of each generated sideband is also a thermal state .
therefore , the multiorder coherent raman process can be used to multiplex the statistical properties of a quantum probe field into a broad comb of different frequencies .    as far as replicating the statistical properties of the input probe into its sidebands
is concerned , the raman medium appears to behave as a linear system .
however , the replication of the normalized autocorrelation functions and normalized squeezing factors of the probe field does not mean the replication of the quantum state at all .
the photon - number distributions and the quantum states of the sidebands evolve in a rather complicated way .
cross - correlations between the sidebands can be generated from initially uncorrelated fields .
an inseparable state can be generated from a separable nonclassical state . although the dynamics of our model system is linear with respect to the field variables , the possibilities of interesting quantum phenomena such as anti - correlation generation , squeezing multiplexing , and entangled - state generation represent the quantum properties that can not be described by the classical statistics of the fields with positive @xmath64 functions .
we have also studied the mixing and modulation of photon statistical properties in the case of two - mode quantum input .
we have shown that the prepared raman coherence and the medium length can be used as control parameters to switch a sideband field from one type of photon statistics to another type , or from a non - squeezed state to a squeezed state and vice versa .
in addition , we can switch nonclassical properties , such as sub - poissonian photon statistics and squeezing , from one frequency to another frequency .
we have also shown an example of quantum interference between the probability amplitudes for a pair of photons with different frequencies .
we have made interesting observations that the multiorder coherent raman scattering behaves in many aspects as a conventional beam splitter and hence can be called a multiorder raman beam splitter .
the two systems have the same underlying physics : the fields are linearly transformed from the input values
. however , the two systems are different in their natures . unlike the conventional beam splitter
, the multiorder coherent raman process can efficiently produce a broad comb of sideband fields whose frequencies are different and are separated by integer multiples of the raman modulation frequency .
the number of generated raman sidebands increases with the effective medium length , which is proportional to the product of the medium length and the prepared raman coherence .
the bessel functions of the effective medium length play a similar role as the transmission and reflection coefficients of a conventional beam splitter .
the ability of the far - off - resonance raman medium to generate a broad comb of fields with similar quantum statistical properties and to switch the quantum statistical characteristics of the radiation fields from one type to another type may find useful applications for high - performance optical communication networks .
in addition , two - photon interference in coherent raman scattering may find various applications for high - precision measurements and also for quantum computation .    finally , we emphasize that the coupling between the raman sidebands can be controlled by the magnitude of the prepared raman coherence , that is , by the intensities of the driving fields . in a realistic far - off - resonance raman medium , such as molecular hydrogen or deuterium vapor @xcite , solid hydrogen @xcite , and rare - earth
doped dielectrics @xcite , a large raman coherence and , consequently , a large number of raman sidebands can be generated by the two - color adiabatic pumping technique .
in such a system , the generation of a broad comb of high - order raman sidebands with nonclassical statistical properties is , in principle , feasible .
therefore , we expect that the coherent raman scattering technique using quantum fields will become a practical and efficient method for a wide range of applications in nonlinear and quantum optics .
v. p. kalosha and j. herrmann , phys .
lett . * 85 * , 1226 ( 2000 ) ; v. p. kalosha and j. herrmann , opt . lett .
* 26 * , 456 ( 2001 ) ; fam le kien , nguyen hong shon , and k. hakuta , phys .
a * 64 * , 051803(r ) ( 2001 ) ; fam le kien , k. hakuta , and a. v. sokolov , _ ibid . _
* 66 * , 023813 ( 2002 ) ; v. kalosha , m. spanner , j. herrmann , and m. ivanov , phys .
lett . * 88 * , 103901 ( 2002 ) ; r. a. bartels , t. c. weinacht , n. wagner , m. baertschy , chris h. greene , m. m. murnane , and h. c. kapteyn , _ ibid . _ * 88 * , 013903 ( 2002 ) .
s. e. harris and a. v. sokolov , phys .
a * 55 * , r4019 ( 1997 ) ; a. v. sokolov , d. d. yavuz , and s. e. harris , opt . lett .
* 24 * , 557 ( 1999 ) ; a. v. sokolov , d. d. yavuz , d. r. walker , g. y. yin , and s. e. harris , phys .
a * 63 * , 051801(r ) ( 2001 ) .
h. p. yuen and j. h. shapiro , ieee trans .
theory * it-26 * , 78 ( 1980 ) ; m. ley and r. loudon , opt . commun . * 54 * , 317 ( 1985 ) ; b. yurke , s. l. mccall , and j. r. klauder , phys . rev . a * 33 * , 4033 ( 1986 ) ; s. prasad , m. o. scully , and w. martienssen , opt .
commun . * 62 * , 139 ( 1987 ) ; z. y. ou , c. k. hong , and l. mandel , opt . commun . *
63 * , 118 ( 1987 ) ; h. fearn and r. loudon , opt . commun . * 64 * , 485 ( 1987 ) ; r. a. campos , b. e. a. saleh , and m. c. teich , phys . rev .
a * 40 * , 1371 ( 1989 ) .
j. torgerson , d. branning , c. monken , and l. mandel , phys .
a * 204 * , 323 ( 1995 ) ; k. mattle , h. weinfurter , p. g. kwiat , and a. zeilinger , phys .
* 76 * , 4656 ( 1996 ) ; d. bouwmeester , j. pan , k. mattle , m. eibl , h. weinfurter , and a. zeilinger , nature ( london ) * 390 * , 575 ( 1997 ) ; t. c. ralph , n. k. langford , t. b. bell , and a. g. white , phys .
a * 65 * , 062324 ( 2002 ) ; t. b. pittman , b. c. jacobs , and j. d. franson , phys .
rev . lett . * 88 * , 257902 ( 2002 ) ; s. p. walborn , a. n. de oliveira , s. padua , and c. h. monken , phys .
rev . lett . * 90 * , 143601 ( 2003 ) .
m. d. lukin , a. b. matsko , m. fleischhauer , and m. o. scully , phys .
82 * , 1847 ( 1999 ) ; a. s. zibrov , a. b. matsko , o. kocharovskaya , y. v. rostovtsev , g. r. welch , and m. o. scully , _ ibid . _
* 88 * , 103601 ( 2002 ) .      z.
y. ou and l. mandel , phys .
lett . * 61 * , 54 ( 1988 ) ; z. y. ou and l. mandel , _ ibid_. * 62 * , 2941 ( 1989 ) ; h. huang and j. h. eberly , j. mod .
optics * 40 * , 915 ( 1993 ) ; m. o. scully , u. w. rathe , c. su , and g. s. agarwal , opt . commun . * 136 * , 39 ( 1997 ) ; a. k. patnaik and g. s. agarwal , j. mod .
optics * 45 * , 2131 ( 1998 ) .


detecting a very weak signal which is almost inaccessible within the classical ( i.e. , non - quantum ) regime is one of the most important subjects in quantum information science .
a strong motivation to devise such an ultra - precise sensor stems from the field of gravitational wave detection @xcite .
in fact , a variety of linear sensors composed of opto - mechanical oscillators have been proposed @xcite , and several experimental implementations of those systems in various scales have been reported @xcite .
it is well known that in general a linear sensor is subjected to two types of fundamental noises , i.e. , the _ back - action noise _ and the _ shot noise_. as a consequence , the measurement noise is lower bounded by the _
standard quantum limit ( sql ) _
@xcite , which is mainly due to the presence of back - action noise .
hence , high - precision detection of a weak signal requires us to devise a sensor that evades the back - action noise and eventually beats the sql ; i.e. , we need to have a sensor achieving _ back - action evasion ( bae)_. in fact , many bae methods have been developed especially in the field of gravitational wave detection , e.g. , the variational measurement technique @xcite or the quantum locking scheme @xcite .
moreover , towards more accurate detection , recently we find some high - level approaches to design a bae sensor , based on those specific bae methods .
for instance , ref .
@xcite provides a systematic comparison of several bae methods and gives an optimal solution .
also systems and control theoretical methods have been developed to synthesize a bae sensor for a specific opto - mechanical system @xcite ; in particular , the synthesis is conducted by connecting an auxiliary system to a given plant system by _ direct - interaction _
@xcite or _ coherent feedback _
@xcite .    along this research direction
, therefore , in this paper we set the goal to develop a general systems and control theory for engineering a sensor achieving bae , for both the coherent feedback and the direct - interaction configurations .
the key tool used here is the _ geometric control theory _
@xcite , which had been developed a long time ago .
this is indeed a beautiful theory providing a variety of controller design methods for various purposes such as the non - interacting control and the disturbance decoupling problem , but , to our best knowledge , it has not been applied to problems in quantum physics .
actually in this paper we first demonstrate that the general synthesis problem of a bae sensor can be formulated and solved within the framework of geometric control theory , particularly the above - mentioned disturbance decoupling problem . +
this paper is organized as follows .
section  [ sec:2 ] is devoted to some preliminaries including a review of the geometric control theory , the general model of linear quantum systems , and the idea of bae .
then , in section  [ sec:3 ] , we provide the general theory for designing a coherent feedback controller achieving bae , and demonstrate an example for an opto - mechanical system . in section  [ sec:4 ] , we discuss the case of direct interaction scheme , also based on the geometric control theory .
finally , in section  [ sec:5 ] , for a realistic opto - mechanical system subjected to a thermal environment ( the perfect bae is impossible in this case ) , we provide a convenient method to find an approximated bae controller and show how much the designed controller can suppress the noise . + for a matrix @xmath0 , @xmath1 , @xmath2 , and @xmath3 represent the transpose , hermitian conjugate , and element - wise complex conjugate of @xmath4 , respectively ; i.e. , @xmath5 , @xmath6 , and @xmath7 .
@xmath8 and @xmath9 denote the real and imaginary parts of a complex number @xmath10 . @xmath11 and @xmath12 denote the zero matrix and the @xmath13 identity matrix . @xmath14 and
@xmath15 denote the kernel and the image of a matrix @xmath4 , i.e. , @xmath16 and @xmath17 .
let us consider the following _ classical _ linear time - invariant system : @xmath18 where @xmath19 is a vector of system variables , @xmath20 and @xmath21 are vectors of input and output , respectively .
@xmath22 , and @xmath23 are real matrices . in the laplace domain ,
the input - output relation is represented by @xmath24 where @xmath25 and @xmath26 are the laplace transforms of @xmath27 and @xmath28 , respectively .
@xmath29 is called the _ transfer function_. in this subsection , we assume @xmath30 .
now we describe the geometric control theory , for the disturbance decoupling problem @xcite .
the following _ invariant subspaces _ play a key role in the theory .    * definition 1 : *   let @xmath31 be a linear map .
then , a subspace @xmath32 is said to be @xmath4-_invariant _ , if @xmath33 .    * definition 2 : *   given a linear map @xmath31 and a subspace im@xmath34 @xmath35 , a subspace @xmath32 is said to be @xmath36-_invariant _ , if @xmath37 .    *
definition 3 : *   given a linear map @xmath31 and a subspace ker@xmath38 @xmath35 , a subspace @xmath39 is said to be @xmath40-_invariant _ , if @xmath41 .
* definition 4 : *   assume that @xmath42 is @xmath40-invariant , @xmath43 is @xmath36-invariant , and @xmath44 .
then , @xmath45 is said to be a @xmath46-pair .    from definitions 2 and 3 , we have the following two lemmas .
* lemma 1 : *   @xmath32 is @xmath36-invariant if and only if there exists a matrix @xmath47 such that @xmath48    * lemma 2 : *   @xmath32 is @xmath40-invariant if and only if there exists a matrix @xmath49 such that @xmath50    the disturbance decoupling problem is described as follows .
the system of interest is represented , in an extended form of eq .
, as @xmath51 where @xmath52 is the disturbance and @xmath53 is the output to be regulated .
@xmath54 and @xmath55 are real matrices .
the other output @xmath28 may be used for constructing a feedback controller ; see fig .
[ fig_feedbackframework ] .
the disturbance @xmath52 can degrade the control performance evaluated on @xmath53 . thus it is desirable if we can modify the system structure by some means so that eventually @xmath52 dose not affect at all on @xmath53 . to @xmath56 is zero for all @xmath57 , for the modified system . or equivalently , the controllable subspace with respect to @xmath52 is contained in the unobservable subspace with respect to @xmath53 . ] this control goal is called the disturbance decoupling . here
we describe a specific feedback control method to achieve this goal ; note that , as shown later , the direct - interaction method for linear quantum systems can also be described within this framework .
the controller configuration is illustrated in fig .
[ fig_feedbackframework ] ; that is , the system modification is carried out by combining an auxiliary system ( controller ) with the original system ( plant ) , so that the whole closed - loop system satisfies the disturbance decoupling condition . the controller with variable @xmath58
is assumed to take the following form : @xmath59 where @xmath60 , @xmath61 , @xmath62 , and @xmath63 are real matrices .        then , the closed - loop system defined in the augmented space @xmath64 is given by @xmath65            & = \left[\begin{array}{cc }                 a+bd_{\scriptscriptstyle k}c & bc_{\scriptscriptstyle k}\\                 b_{\scriptscriptstyle k}c & a_{\scriptscriptstyle k}\\            \end{array}\right ]            \left[\begin{array}{c }                 x\\                 x_{{\scriptscriptstyle k}}\\            \end{array}\right ]           + \left[\begin{array}{c }                e\\                o\\            \end{array}\right ] d , ~ \nonumber \\          z&=\left[\begin{array}{cc }                h & o\\           \end{array}\right ]          \left[\begin{array}{c }                x\\                x_{{\scriptscriptstyle k}}\\           \end{array}\right].\end{aligned}\ ] ] the control goal is to design @xmath66 so that , in eq .  , the disturbance signal @xmath52 dose not appear in the output @xmath53 : see the below footnote . here
, let us define @xmath67 ,   \label{eq_extended_matrixae}\end{aligned}\ ] ] @xmath68 , @xmath69 , @xmath70 , and @xmath71 .
then , the following theorem gives the solvability condition for the disturbance decoupling problem .
* theorem 1 : *   for the closed - loop system , the disturbance decoupling problem via the dynamical feedback controller has a solution if and only if there exists a @xmath46-pair @xmath72 satisfying @xmath73 note that this condition does not depend on the controller matrices to be designed . the following corollary can be used to check if the solvability condition is satisfied .
* corollary 1 : *   for the closed - loop system , the disturbance decoupling problem via the dynamical feedback controller has a solution if and only if @xmath74 where @xmath75 is the maximum element of @xmath36-invariant subspaces contained in @xmath76 , and @xmath77 is the minimum element of @xmath40-invariant subspaces containing @xmath78 . these subspaces can be computed by the algorithms given in appendix  a.    once the solvability condition described above is satisfied , then we can explicitly construct the controller matrices @xmath79 . the following intersection and projection subspaces play a key role for this purpose ; that is , for a subspace @xmath80 , let us define @xmath81                   \in \mathcal{v}_{{\scriptscriptstyle e } } \right\ } , ~ \\               \mathcal{v}_{{\scriptscriptstyle p}}:&=                           \left\{x \in \mathcal{x}~\bigg| \left[\begin{array}{c }                                      x\\                                        x_{{\scriptscriptstyle k}}\\                                 \end{array}\right ]                   \in \mathcal{v}_{{\scriptscriptstyle e}},~                   \exists { x}_{{\scriptscriptstyle k } } \in                        \mathcal{x}_{{\scriptscriptstyle k } } \right\}.\end{aligned}\ ] ] then
, the following theorem is obtained :    * theorem 2 : *   suppose that @xmath82 is a @xmath83-pair .
then , there exist @xmath84 , @xmath85 , and @xmath63 such that @xmath86 and @xmath87 hold , where @xmath88
.    moreover , there exists @xmath89 with @xmath90 , and @xmath91 has an invariant subspace @xmath92 such that @xmath93 and @xmath94 . also , @xmath95 satisfies @xmath96 where @xmath97 is a linear map satisfying @xmath98 .
in fact , under the condition given in theorem 2 , let us define the following augmented subspace @xmath99 : @xmath100           \bigg|~ x \in \mathcal{v}_{2 } \right\}.\end{aligned}\ ] ] then , @xmath93 and @xmath94 hold , and we have @xmath101         & = \left[\begin{array}{cc }                     a+bd_{\scriptscriptstyle k}c & bc_{\scriptscriptstyle k}\\                     b_{\scriptscriptstyle k}c & a_{\scriptscriptstyle k}\\            \end{array}\right ]            \left[\begin{array}{c }                     x\\                     nx\\            \end{array}\right ]   \nonumber\\       & = \left[\begin{array}{c }                     ( a+bf)x\\                     n(a+bf)x\\            \end{array}\right ]     \in \mathcal{v}_{{\scriptscriptstyle e}},\end{aligned}\ ] ] implying that @xmath102 is actually @xmath91-invariant .
now suppose that theorem  1 holds , and let us take the @xmath46-pair @xmath82 satisfying eq .  .
then , together with the above result ( @xmath103 ) , we have @xmath104^{\top } \subseteq \mathcal{v}_{{\scriptscriptstyle e } }   \subseteq { { \rm ker}\hspace{0.07cm}}[h~o]$ ] .
this implies that @xmath52 must be contained in the unobservable subspace with respect to @xmath53 , and thus the disturbance decoupling is realized .
here we describe a general linear quantum system composed of @xmath105 bosonic subsystems .
the @xmath106-th mode can be modeled as a harmonic oscillator with the canonical conjugate pairs ( or quadratures ) @xmath107 and @xmath108 satisfying the canonical commutation relation ( ccr ) @xmath109 .
let us define the vector of quadratures as @xmath110^{\top}$ ] .
then , the ccrs are summarized as @xmath111 .
\end{aligned}\ ] ] note that @xmath112 is a @xmath113 block diagonal matrix .
the linear quantum system is an open system coupled to @xmath114 environment fields via the interaction hamiltonian @xmath115 , where @xmath116 is the field annihilation operator satisfying @xmath117 .
also @xmath118 is given by @xmath119 with @xmath120 .
in addition , the system is driven by the hamiltonian @xmath121 with @xmath122 .
then , the heisenberg equation of @xmath123 is given by @xmath124 where @xmath125 is defined by @xmath126         = \left[\begin{array}{cc }               ( \hat{a}_j+\hat{a}_j^*)/\sqrt{2 } \\               ( \hat{a}_j-\hat{a}_j^*)/\sqrt{2 } i \\
\end{array}\right].\ ] ] the matrices are given by @xmath127 and @xmath128 with @xmath129^\top \in \mathbb{r}^{2 \times 2n}$ ] .
also , the instantaneous change of the field operator @xmath125 via the system - field coupling is given by @xmath130 summarizing , the linear quantum system is characterized by the dynamics and the output , which are exactly of the same form as those in eq .
( @xmath131 in this case ) .
however note that the system matrices have to satisfy the above - described special structure , which is equivalently converted to the following _ physical realizability condition _
@xcite : @xmath132      the opto - mechanical oscillator illustrated in fig .
[ fig_opticalsensor ] is a linear quantum system , which serves as a sensor for a very weak signal .
let @xmath133 and @xmath134 be the oscillator s position and momentum operators , and @xmath135 represents the annihilation operator of the cavity mode .
the system hamiltonian is given by @xmath136 ; that is , the oscillator s free evolution with resonant frequency @xmath137 plus the linearized radiation pressure interaction between the oscillator and the cavity field with coupling strength @xmath138 .
the system couples to an external probe field ( thus @xmath139 ) via the coupling operator @xmath140 , with @xmath141 the coupling constant between the cavity and probe fields . the corresponding matrix @xmath142 and vector @xmath143
are then given by @xmath144,~~             c_{1}=\sqrt{\frac{\kappa}{2 } }                  \left[\begin{array}{c }                     0 \\ 0 \\ 1 \\ i \\                  \end{array}\right].\end{aligned}\ ] ] the oscillator is driven by an unknown force @xmath145 with coupling constant @xmath146 ; then the vector of system variables @xmath147^\top$ ] satisfies @xmath148 where @xmath149 , \ ,             b=\sqrt{\gamma}\left[\begin{array}{c }                        0 \\                        1 \\                        0 \\                        0 \\                 \end{array}\right ] , \nonumber \\
c_1&=-b_1^{\top }                = \left[\begin{array}{cccc }                       0 & 0 & \sqrt{\kappa } & 0\\                       0 & 0 & 0 & \sqrt{\kappa}\\                  \end{array}\right ] , \nonumber \\
\hat{w}_{1 }                   & = [ \hat{q}_{1 } , \hat{p}_{1}]^{\top } , ~~           \hat{w}_{1}^{\rm out }                    = [ \hat{q}_{1}^{\rm out } , \hat{p}_{1}^{\rm out}]^{\top}. \label{eq_matrix_originalplant}\end{aligned}\ ] ] note that we are in the rotating frame at the frequency of the probe field .
these equations indicate that the information about @xmath150 can be extracted by measuring @xmath151 by a homodyne detector .
actually the measurement output in the laplace domain is given by @xmath152 where @xmath153 , @xmath154 , and @xmath155 are transfer functions given by @xmath156 thus , @xmath151 certainly contains @xmath150 .
note however that it is subjected to two noises .
the first one , @xmath157 , is the back - action noise , which is due to the interaction between the oscillator and the cavity .
the second one , @xmath158 , is the shot noise , which inevitably appears .
now , the normalized output is given by @xmath159 and the normalized noise power spectral density of @xmath160 in the fourier domain @xmath161 is calculated as follows : @xmath162 the lower bound is called the sql .
note that the last inequality is due to the heisenberg uncertainty relation of the normalized noise power , i.e. , @xmath163 .
hence , the essential reason why sql appears is that @xmath164 contains both the back - action noise @xmath165 and the shot noise @xmath166 .
therefore , toward the high - precision detection of @xmath150 , we need bae ; that is , the system structure should be modified by some means so that the back - action noise is completely evaded in the output signal ( note that the shot noise can never be evaded ) .
the condition for bae can be expressed in terms of the transfer function as follows @xcite ; i.e. , for the modified ( controlled ) sensor , the transfer function from the back - action noise to the measurement output must satisfy @xmath167 equivalently , @xmath164 contains only the shot noise @xmath158 ; hence , in this case the signal to noise ratio can be further improved by injecting a @xmath158-squeezed ( meaning @xmath168 ) probe field into the system .
there are two schemes for controlling a quantum system via feedback .
the first one is the _ measurement - based feedback _
@xcite illustrated in fig .
[ fig_mfvscf ]  ( a ) . in this scheme
, we measure the output fields and feed the measurement results back to control the plant system . on the other hand , in the _ coherent feedback _
scheme @xcite shown in fig .
[ fig_mfvscf ]  ( b ) , the feedback loop dose not contain any measurement component and the plant system is controlled by another quantum system .
recently we find several works comparing the performance of these two schemes @xcite .
in particular , it was shown in @xcite that there are some control tasks that can not be achieved by any measurement - based feedback but can be done by a coherent one .
more specifically , those tasks are realizing bae measurement , generating a quantum non - demolished variable , and generating a decoherence - free subsystem ; in our case , of course , the first one is crucial .
hence , here we aim to develop a theory for designing a coherent feedback controller such that the whole controlled system accomplishes bae .
as discussed in section [ sec:2 - 1 ] , the geometric control theory for disturbance decoupling problem is formulated for the controlled system with special structure ; in particular , the coefficient matrix of the disturbance @xmath52 is of the form @xmath169^\top$ ] and that of the state vector in the output @xmath53 is @xmath170 $ ] .
here we consider a class of coherent feedback configuration such that the whole closed - loop system dynamics has this structure , in order for the geometric control theory to be directly applicable .
first , for the plant system given by eqs .   and , we assume that the system couples to all the probe fields in the same way ; i.e. , @xmath171 this immediately leads to @xmath172 . next , as the controller , we take the following special linear quantum system with @xmath173 input - output fields : @xmath174 where the matrices @xmath175 satisfy the physical realizability condition . note that , corresponding to the plant structure , we assumed that the controller couples to all the fields in the same way , specified by @xmath176 . here
we emphasize that the number of channels , @xmath114 , should be as small as possible from a viewpoint of implementation ; hence in this paper let us consider the case @xmath177 .
now , we consider the coherent feedback connection illustrated in fig .
[ fig_3i3o_general ] , i.e. , @xmath178        where @xmath179 and @xmath180 are @xmath181 unitary matrices representing the scattering process of the fields ; recall that the scattering process @xmath182 with @xmath183 the phase shift can be represented in the quadrature form as @xmath184             = s(\theta )                \left[\begin{array}{cc }                   \hat{q } \\
\hat{p } \\
\end{array}\right ]              = \left[\begin{array}{cc }                     \cos \theta & -\sin \theta \\
\sin \theta & \cos \theta \\                 \end{array}\right ]                 \left[\begin{array}{cc }                   \hat{q } \\                   \hat{p } \\                \end{array}\right ]
. \label{eq_phaseshifter}\ ] ]        here we apply the coherent feedback scheme elaborated in section  [ sec:3 - 2 ] to the opto - mechanical system studied in section  [ sec:2 - 3 ] .
the goal is , as mentioned before , to determine the controller matrices @xmath95 such that the closed - loop system achieves bae . here
, we provide a step - by - step procedure to solve this problem ; the relationships of the class of controllers determined in each step is depicted in fig .
[ fig_controllerconstraint ] .
-phase shifter corresponding to eq .  .
, width=264 ]      * ( ii ) * in the next step we apply theorem  1 to check if there exists a feedback controller such that the above closed - loop system achieves bae ; recall that the necessary and sufficient condition is eq .  , i.e. , @xmath206 , where now @xmath207 \right\ } , \\         \mathcal{h}&={{\rm ker}\hspace{0.07cm}}h={\rm span } \left\ {                      \left[\begin{array}{c }                           1\\                           0\\                           0\\                           0\\                      \end{array}\right ] ,                      \left[\begin{array}{c }                           0\\                           1\\                           0\\                           0\\                      \end{array}\right ] ,                      \left[\begin{array}{c }                           0\\                           0\\                           1\\                           0\\                      \end{array}\right ] \right\}. \end{aligned}\ ] ] to check if this solvability condition is satisfied , we use corollary 1 ; from @xmath208 and @xmath209 , the algorithms given in appendix  a yield @xmath210 implying that the condition in corollary 1 , i.e. , @xmath211 , is satisfied .
thus , we now see that the bae problem is solvable , as long as there is no constraint on the controller parameters .
next we aim to determine the controller matrices @xmath95 , using theorem  2 .
first we set @xmath212 and @xmath213 ; note that @xmath82 is a @xmath83-pair .
then , from theorem 2 , there exists a feedback controller with dimension @xmath214 .
moreover , noting again that @xmath196 , there exist matrices @xmath84 , @xmath85 , and @xmath215 such that @xmath216 these conditions lead to @xmath217 , ~~                g=\left[\begin{array}{cc }                          0&g_{\scriptscriptstyle 12}\\                          -\frac{g}{\sqrt{\kappa}}&g_{\scriptscriptstyle 22}\\                          g_{\scriptscriptstyle 31}&g_{\scriptscriptstyle 32}\\                          0&\sqrt{\kappa}\\                       \end{array}\right ] ,   \\
n&=\left[\begin{array}{cccc }                          n_{\scriptscriptstyle 11 } & n_{\scriptscriptstyle 12 }                                & 0 & n_{\scriptscriptstyle 14 } \\
n_{\scriptscriptstyle 21 } & n_{\scriptscriptstyle 22 }                               & 0 & n_{\scriptscriptstyle 24}\\                       \end{array}\right ] , \end{aligned}\ ] ] where @xmath218 , and @xmath219 are free parameters
. then the controller matrices @xmath95 can be identified by eq .   with the above matrices @xmath220 ; specifically , by substituting @xmath221 , @xmath222 , and @xmath223 in eq .
, we have @xmath224 which yield @xmath225 where @xmath226 is the right inverse to @xmath215 , i.e. , @xmath227 .    *
( iii ) * note again that the controller has to satisfy the physical realizability condition , which is now @xmath228 and @xmath229 .
these constraints are represented in terms of the parameters as follows : @xmath230 where @xmath231 and @xmath232 .
this is one of our main results ; the linear controller achieving bae for the opto - mechanical oscillator can be fully parametrized by eq .   satisfying the condition .
we emphasize that this full parametrization of the controller can be obtained thanks to the general problem formulation based on the geometric control theory .        *
( iv ) * in practice , of course , we need to determine a concrete set of parameters to construct the controller . especially here
let us consider a _ passive system _ ; this is a static quantum system such as an empty optical cavity .
the main reason for choosing a passive system rather than a non - passive ( or _ active _ ) one such as an optical parametric oscillator is that , due to the external pumping energy , the latter could become fragile and also its physical implementation must be more involved compared to a passive system @xcite . now the condition for the system @xmath233 to be passive is given by @xmath234 and @xmath235 ; the general result of this fact is given in theorem  3 in appendix  b. from these conditions , the system parameters are imposed to satisfy , in addition to eq .  , the following equalities : @xmath236 there is still some freedom in determining @xmath237 , which however corresponds to simply the phase shift at the input - output ports of the controller , as indicated from eq .  .
thus , the passive controller achieving bae in this example is unique up to the phase shift . here
particularly we chose @xmath238 and @xmath239 .
then the controller matrices satisfying eqs .   and are determined as @xmath240,~~          c_{\scriptscriptstyle k}=-b_{\scriptscriptstyle k}^{\top }                      = \left[\begin{array}{cc }                           \frac{g}{\sqrt{\kappa } } & 0 \\                           0 & \frac{g}{\sqrt{\kappa } } \\
\end{array}\right ] .   \ ] ] as illustrated in fig .
[ fig_3i3obae ] , the controller specified by these matrices can be realized as a single - mode , 2-inputs and 2-outputs optical cavity with decay rate @xmath241 and detuning @xmath242 .
in other words , if we take the cavity with the following hamiltonian and the coupling operator ( @xmath243 is the cavity mode ) @xmath244 then to satisfy the bae condition the controller parameters @xmath245 must satisfy @xmath246 summarizing , the above - designed sensing system composed of the opto - mechanical oscillator ( plant ) and the optical cavity ( controller ) , which are combined via coherent feedback , satisfies the bae condition .
hence , it can work as a high - precision detector of the force @xmath150 below the sql , particularly when the @xmath247-squeezed probe input field is used ; this fact will be demonstrated in section  [ sec:5 ] .
in this section , we study another control scheme for achieving bae . as illustrated in fig . [ fig_mankei_scheme ]
( a ) , the controller in this case is directly connected to the plant , not through a coherent feedback ; hence this scheme is called the _ direct interaction_. the controller is characterized by the following two hamiltonians : @xmath248 where @xmath249^{\top}$ ] is the vector of controller variables with @xmath250 the number of modes of the controller .
@xmath251 is the controller s self hamiltonian with @xmath252 . also @xmath253 with @xmath254
, @xmath255 represents the coupling between the plant and the controller . note that , for the hamiltonians @xmath251 and @xmath253 to be hermitian , the matrices must satisfy @xmath256 and @xmath257 ; these are the physical realizability conditions in the scenario of direct interaction .
in particular , here we consider a plant system interacting with a single probe field @xmath197 , with coupling matrices @xmath258 and @xmath259 .
then , the whole dynamics of the augmented system with variable @xmath260^{\top}$ ] is given by @xmath261 where @xmath262 , ~~        b_{{\scriptscriptstyle e}}=\left[\begin{array}{c }                         b\\                         o \\
\end{array}\right ] ,   \nonumber   \\        c_{{\scriptscriptstyle e}}&=\left[\begin{array}{cc }                         c & o   \\                  \end{array}\right ] , ~~         b_{{\scriptscriptstyle e}}=\left[\begin{array}{cc }                          b^{\top } & o\\                       \end{array}\right ] ^{\top}. \end{aligned}\ ] ] note that @xmath263 , @xmath264 , and @xmath265 are the same matrices as those in eq .  .
also , comparing the matrices and , we have that @xmath266 , which thus leads to @xmath267 and @xmath268 in theorem  2 .
now , again for the opto - mechanical system illustrated in fig .
[ fig_opticalsensor ] , let us aim to design the direct interaction controller , so that the whole system achieves bae ; that is , the problem is to determine the matrices @xmath269 so that the back - action noise @xmath270 does not appear in the measurement output @xmath271 .
for this purpose , we go through the same procedure as that taken in section  [ sec:3 - 3 ] .    *
( i ) * because of the structure of the matrices @xmath263 and @xmath264 , the system is already of the form , where the geometric control theory is directly applicable .    *
( ii ) * because we now focus on the same plant system as that in section [ sec:3 - 3 ] , the same conclusion is obtained ; that is , the bae problem is solvable as long as there is no constraint on the controller matrices @xmath269 .
the controller matrices can be determined in a similar way to section [ sec:3 - 3 ] as follows .
first , because the @xmath46-pair @xmath82 is the same as before , it follows that @xmath272 , i.e. , @xmath273 .
then , from theorem  2 with the fact that @xmath267 and @xmath268 , we find that the direct interaction controller can be parameterized as follows : @xmath274 the matrices @xmath47 , @xmath49 , and @xmath215 satisfy @xmath275 , and @xmath98 , which lead to @xmath276 , ~~                g=\left[\begin{array}{cc }                          0&g_{\scriptscriptstyle 12}\\                          -\frac{g}{\sqrt{\kappa}}&g_{\scriptscriptstyle 22}\\                          g_{\scriptscriptstyle 31}&g_{\scriptscriptstyle 32}\\
0&0\\                       \end{array}\right ] , \nonumber   \\                n&=\left[\begin{array}{cccc }                          n_{\scriptscriptstyle 11}&n_{\scriptscriptstyle 12 }                                & 0&n_{\scriptscriptstyle 14}\\                          n_{\scriptscriptstyle 21}&n_{\scriptscriptstyle 22 }                                & 0&n_{\scriptscriptstyle 24}\\                       \end{array}\right],\end{aligned}\ ] ] where @xmath218 , and @xmath219 are free parameters .    * ( iii ) * the controller matrices have to satisfy the physical realizability conditions @xmath256 and @xmath257 ; these constraints impose the parameters to satisfy @xmath277 where @xmath231 and @xmath232 .
+ equations , , and provide the full parametrization of the direct interaction controller .    * ( iv ) * to specify a set of parameters , as in the case of section [ sec:3 - 3 ] , let us aim to design a passive controller . from theorem  4 in appendix  b , @xmath278 and
@xmath279 satisfy the condition @xmath280 and @xmath281 , which lead to the same equalities given in eq .  .
then , setting the parameters to be @xmath238 and @xmath239 , we can determine the matrices @xmath278 and @xmath282 as follows : @xmath283 , ~~                  r_2=r_1^\top                     = \left[\begin{array}{cccc }                               0&0&g&0 \\                               0&0&0&g \\
\end{array}\right ] .
\end{aligned}\ ] ] the controller specified by these matrices can be physically implemented as illustrated in fig .
[ fig_mankei_scheme ]  ( b ) ; that is , it is a single - mode detuned cavity with hamiltonian @xmath284 , which couples to the plant through a beam - splitter ( bs ) represented by @xmath285 .
* remark : * we can employ an active controller , as proposed in @xcite . in this case
the interaction hamiltonian is given by @xmath286 , while the system s self - hamiltonian is the same as above ; @xmath284 . that is , the controller couples to the plant through a non - degenerate optical parametric amplification process in addition to the bs interaction . to satisfy the bae condition
, the parameters must satisfy @xmath287 . note that this direct interaction controller can be specified , in the full - parameterization , , and , by @xmath288
we have demonstrated in sections [ sec:3 - 3 ] and [ sec:4 ] that the bae condition can be achieved by engineering an appropriate auxiliary system and connecting it to the plant .
however , in a practical situation , it can not be expected to realize such perfect bae due to several experimental imperfections .
hence , in a realistic setup , we should modify our strategy for engineering a sensor so that it would accomplish _ approximate bae_. then , looking back into section [ sec:2 - 3 ] where the bae condition , @xmath289 , was obtained , we are naturally led to consider the following optimization problem to design an auxiliary system achieving the approximate bae : @xmath290 where @xmath291 denotes a valid norm of a complex function .
in particular , in the field of robust control theory , the following @xmath292 norm and the @xmath293 norm are often used @xcite : @xmath294 that is , the @xmath292 or @xmath293 control theory provides a general procedure for synthesizing a feedback controller that minimizes the above norm . in this paper
, we take the @xmath292 norm , mainly owing to the broadband noise - reduction nature of the @xmath292 controller .
then , rather than pursuing an optimal quantum @xmath292 controller based on the quantum @xmath292 control theory @xcite , here we take the following geometric - control - theoretical approach to solve the problem .
that is , first we apply the method developed in section [ sec:3 ] or [ sec:4 ] to the idealized system and obtain the controller achieving bae ; then , in the practical setup containing some unwanted noise , we make a local modification of the controller parameters obtained in the first step , to minimize the cost @xmath295 .    as a demonstration , here we consider the coherent feedback control for the opto - mechanical system studied in section [ sec:2 - 3 ] , which is now subjected to the thermal noise @xmath296 . following the above - described policy , we employ the coherent feedback controller constructed for the idealized system that ignores @xmath296 , leading to the controller given by eqs .   and , illustrated in fig .
[ fig_3i3obae ] . the closed - loop system with variable @xmath297^{\top}$ ] , which now takes into account the realistic imperfections
, then obeys the following dynamics : @xmath298 where @xmath299 .   \ ] ] @xmath263 , @xmath264 , and @xmath265 are the same matrices given in eq .  .
@xmath296 is the thermal noise satisfying @xmath300 , where @xmath301 is the mean phonon number at thermal equilibrium @xcite .
note that the damping effect appears in the @xmath302 component of @xmath303 due to the stochastic nature of @xmath296 . also , again , @xmath304 and @xmath305 are the decay rate and the detuning of the controller cavity , respectively . in the idealized
setting where @xmath296 is negligible , the perfect bae is achieved by choosing the parameters satisfying eq .  .
the measurement output of this closed - loop system is , in the laplace domain , represented by @xmath306 the normalized noise power spectral density of @xmath307 is calculated as @xmath308 the coefficient of the back - action noise is given by @xmath309 our goal is to find the optimal parameters @xmath310 that minimize the @xmath292 norm of the transfer function , @xmath311 .
norm @xmath312 versus the coupling constant @xmath304 and the detuning @xmath305.,width=264 ]    the system parameters are taken as follows @xcite : @xmath313 mhz , @xmath314 mhz , @xmath315 khz , @xmath316 mhz , @xmath317 , and the effective mass is @xmath318 kg . we then have fig .
[ fig_h2norm ] , showing @xmath319 as a function of @xmath304 and @xmath305 .
this figure shows that there exists a unique pair of @xmath320 that minimizes the norm , and they are given by @xmath321 mhz and @xmath322 mhz , which are actually close to the ideal values .
[ fig_noisepower ] shows the value of eq .   with these optimal parameters @xmath320 , where the noise floor @xmath323 is subtracted .
the solid black line represents the sql , which is now given by @xmath324 then the dot - dashed blue and dotted green lines indicate that , in the low frequency range , the coherent feedback controller can suppress the noise below the sql , while , by definition , the noise power of the autonomous ( i.e. , uncontrolled ) plant system is above the sql .
moreover , this effect can be enhanced by injecting a @xmath158-squeezed probe field ( meaning @xmath325 and @xmath326 ) into the system .
in fact the dashed red line in the figure illustrates the case @xmath327 ( about 9 db squeezing ) , showing the significant reduction of the noise power .
the main contribution of this paper is in that it first provides the general theory for constructing a back - action evading sensor for linear quantum systems , based on the well - developed classical geometric control theory .
the power of the theory has been demonstrated by showing that , for the typical opto - mechanical oscillator , a full parametrization of the auxiliary coherent - feedback and direct interaction controller achieving bae was derived , which contains the result of @xcite .
note that , although we have studied a simple example for the purpose of demonstration , the real advantage of the theory developed in this paper will appear when dealing with more complicated multi - mode systems such as an opto - mechanical system containing a membrane @xcite .
another contribution of this paper is to provide a general procedure for designing an approximate bae sensor under realistic imperfections ; that is , an optimal approximate bae system can be obtained by solving the minimization problem of the transfer function from the back - action noise to the measurement output .
while in section  [ sec:5 ] we have provided a simple approach based on the geometric control theory for solving this problem , the @xmath292 or @xmath293 control theory could be employed for systematic design of an approximate bae controller even for the above - mentioned complicated system .
this is also an important future research direction of this work .
the set of @xmath36-invariant subspaces has a unique maximum element contained in a given subspace @xmath330 .
this space , denoted by @xmath331 , can be computed by the following algorithm : @xmath332    similarly , the set of @xmath40-invariant subspaces has a unique minimum element containing a given subspace @xmath333 , and this space , denoted by @xmath334 , can be computed by the following algorithm : @xmath335
note that a passive system is constituted only with annihilation operator variables ; a typical optical realization of the passive system is an empty optical cavity .
moreover , @xmath345 is already satisfied and @xmath344 leads to @xmath346 .
this is the reason why it is sufficient to consider the constraints only on @xmath347 and @xmath348 .
then the goal here is to represent the conditions @xmath349 in terms of the coefficient matrices of eq .  .
for this purpose , let us introduce the permutation matrix @xmath350 as follows ; for a column vector @xmath351^{\top}$ ] , @xmath350 is defined through @xmath352^{\top}$ ] .
note that @xmath350 satisfies @xmath353 .
then , the coefficient matrices of the above two system representations are connected by @xmath354 where @xmath355.\end{aligned}\ ] ] @xmath356 , @xmath357 , and @xmath358 have the same forms as above .
then , we have the following theorem , providing the passivity condition in the quadrature form :    * theorem 3 : *   the system is passive if and only if , in addition to the physical realizability condition , the following equalities hold : @xmath359 _ proof : _
let us first define @xmath360 $ ] , which leads to @xmath361 .
then , we can prove @xmath362 the condition in the right hand side is equivalent to @xmath363 and @xmath364 , which thus leads to @xmath343 .
also , from a similar calculation we obtain @xmath365    let us next consider the passivity condition of the direct interaction controller discussed in section  [ sec:4 ] .
the setup is that , for a given linear quantum system , we add an auxiliary component with variable @xmath366 , which is characterized by the hamiltonians .
the point is that these hamiltonians have the following equivalent representations in terms of the vector of annihilation operators @xmath367 and @xmath368 : @xmath369 \mathscr{r}_{\scriptscriptstyle k }                      \left[\begin{array}{c }                             \hat{a}_{\scriptscriptstyle k }   \\                             \hat{a}_{\scriptscriptstyle k}^{\sharp } \\
\end{array}\right ] , \nonumber \\            \hat{h}_{\rm int}&=\frac{1}{2}\left (                      \left[\begin{array}{cc }                            \hat{a}^{\dag } & \hat{a}^{\top }
\\                      \end{array}\right ] \mathscr{r}_1                      \left[\begin{array}{c }                             \hat{a}_{\scriptscriptstyle k }   \\                             \hat{a}_{\scriptscriptstyle k}^{\sharp } \\
\end{array}\right ]                   + \left[\begin{array}{cc }                            \hat{a}_{\scriptscriptstyle k}^{\dag }                                  & \hat{a}_{\scriptscriptstyle k}^{\top } \\
\end{array}\right ] \mathscr{r}_2                      \left[\begin{array}{c }                             \hat{a }   \\
\hat{a}^{\sharp } \\
\end{array}\right ] \right).\end{aligned}\ ] ] the matrices @xmath370 , and @xmath371 are of the same forms as those in eq .  .
note that they have to satisfy the physical realizability conditions @xmath372 and @xmath373 . now we can define the passivity property of the direct interaction controller ; that is , if the hamiltonians does not contain any quadratic term such as @xmath374 and @xmath375 , then the direct interaction controller is passive .
the formal definition is given as follows :    * definition 6 : *   the direct interaction controller constructed by hamiltonians is said to be passive if , in addition to the physical realizability conditions @xmath372 and @xmath373 , the matrices satisfy @xmath376 and @xmath377 .
* theorem 4 : *   the direct interaction controller constructed by hamiltonians is passive if and only if , in addition to the physical realizability conditions @xmath256 and @xmath257 , the following equalities hold : @xmath378      c. m. caves , k. s. thorne , r. w. p. drever , v. d. sandberg , and m. zimmermann , on the measurement of a weak classical force coupled to a quantum - mechanical oscillator
. i. issues of principle , rev .
* 52 * , 341 ( 1980 )              c. l. latune , b. m. escher , r. l. de matos filho , and l. davidovich , quantum limit for the measurement of a classical force coupled to a noisy quantum - mechanical oscillator , phys .
a * 88 * , 042112 ( 2013 ) .
t. corbitt , y. chen , e. innerhofer , h. m@xmath379ller - ebhardt , d. ottaway , h. rehbein , d. sigg , s. whitcomb , c. wipf , and n. mavalvala , an all - optical trap for a gram - scale mirror , phys .
lett . * 98 * , 150802 ( 2007 ) .
m. underwood , d. mason , d. lee , h. xu , l. jiang , a. b. shkarin , k. brkje , s. m. girvin , and j. g. e. harris , measurement of the motional sidebands of a nanogram - scale oscillator in the quantum regime , phys .
a * 92 * , 061801 ( 2015 ) .
h. j. kimble , y. levin , a. b. matsko , k. s. thorne , and s. p. vyatchanin , conversion of conventional gravitational - wave interferometers into quantum nondemolition interferometers by modifying their input and/or output optics , phys .
d * 65 * , 022002 ( 2001 ) .
y. liu , s. shankar , n. ofek , m. hatridge , a. narla , k. m. sliwa , l. frunzio , r. j. schoelkopf , and m. h. devoret , comparing and combining measurement - based and driven - dissipative entanglement stabilization , phys .
x * 6 * , 011022 ( 2016 ) .