it is well known that the fiber optical parametric amplifiers ( fopas ) , employing @xmath0 nonlinearity based four wave mixing ( fwm ) to transfer energy from one or two strong pump fields to a weak signal field , can be used in the conventional optical communication systems for signal - processing applications , such as optical amplification , phase conjugation , optical limiters , 3r generators , and time - domain de - multiplexer etc .
in fact , similar to the @xmath1 crystals based optical parametric amplifiers ( opas )  @xcite , fopas are also candidates for generating non - classical light for quantum communication  @xcite .    in quantum communication
, information can be encoded discretely on each photon or continuously in quadrature - phase amplitudes of an optical field .
the former requires single - photon sources that can be generated with weak nonlinear interaction  @xcite , but the latter relies on strong interaction to generate quantum correlation between the amplitudes of optical fields for entanglement  @xcite .
moreover , continuous variable ( cv ) entanglement offers several advantages over its discrete variable counterpart .
the main advantage is the experimental feasibility in unconditional production of entangled states which in turn enables quantum information protocols to be accomplished unconditionally  @xcite .
the past decade has seen the growing interest in developing nonclassical light sources via the fwm in fibers , because they are compatible with the fiber network  @xcite . in order to enhance the fwm and to suppress the raman scattering  @xcite ,
a pulse pumped fopa is preferred because the peak power of pulsed laser is very high even at a modest average power . however , for the majority of the experiments performed to date , the fwm in fiber is in the low gain regime , and the fiber based nonclassical light sources with vacuum injection are in the domain of a few photons and thus can be only employed for discrete variable information encoding . in order to develop the fiber based source of the cv nonclassical light , the pulse pumped fopa needs to be operated in the high gain regime .    up to now , only a couple of experiments reported the generation of cv nonclassical light via the pulse pumped fwm in fibers @xcite . in 2001 , sharping and co - works realized the pulsed twin beams by using the fopa for the first time @xcite . the amplified signal and
generated idler beams bear a strong quantum correlation in intensity ( photon - number ) fluctuations in the sense that the quantum noise level of their intensity difference is lower than the shot - noise limit ( snl ) . in that experiment ,
the gain of fwm is less than 3 and the observed intensity difference noise of twin beams is only about 1.1 db ( 2.6 db after correction for losses ) lower than the snl .
recently , using the fopa with photon number gain of about 16 db , our group demonstrated the twin beams with the intensity - difference noise below the snl by 3.1 db ( 10.4 db after correction for losses ) @xcite .
the results indicate that compared with its @xmath1 crystal counterparts  @xcite , the pulse pumped high gain fopa is a simple and promising system for generating the cv nonclassical light , which has the potential applications in realizing the quantum communication protocols and in performing the high precision measurement .
there have been a few theoretical papers analyzing the generation of pulsed cv nonclassical light via the high gain @xmath1 crystals based opas  @xcite . for the pulse pumped high gain fopas
, however , the investigation is mostly within the domain of classical characteristics  @xcite and the quantum theory has not been worked out yet .
this is different from the case of continuous wave ( cw ) or quasi - cw pumped high gain fopas  @xcite . in order to understand the high gain fopa for further reducing the intensity difference noise of pulsed twin beams
, we focus on developing a multi - mode quantum theory in this paper .    the rest of the paper is organized as follows .
after briefly introducing the quantum description of a single - mode fopa in sec .
ii , we deduce a multi - mode quantum theory of a pulse pumped non - degenerate fopa in sec .
the evolution of non - degenerate signal and idler beams is described by the bogoliubov transformation in a multi - frequency mode form . according to the hamiltonian of fwm in fibers ,
we derive the four transformation functions of the bogoliubov transformation in the analytical form of infinite series .
the calculation shows that the noise performance of the high gain fopa highly depends on the joint spectral function ( jsf ) . in sec .
iv , we show that the fopa with factorized jsf can be described by a quantum model of single - mode fopa , which is similar to the model given in sec .
ii , because both signal and idler beams can be in single temporal mode , respectively .
the results illustrate the validity of the multi - mode theory . in sec .
v , we focus on analyzing the fopa with a very broad gain bandwidth . in this situation , the jsf is non - factorable .
the calculated results show that although the pulse pumped high gain fopa generally has a noise figure greater than the 3 db quantum limit , the intensity difference noise of the twin beams , weakly depends on the jsf , can be significantly reduced to lower than the snl . to closely resemble the real experimental situation ,
the dependence of the intensity difference noise of twin beams is numerically studied by considering the influences of experimental imperfections , including the quantum efficiencies of detectors , the collection efficiency of twin beams , and the excess noise of input signal .
finally , we conclude in sec .
[ conclusion ] .
we start with the quantum description of single - mode linear fopas , whose pumps are single frequency lasers with negligible depletions . in this section
, we will introduce the bogoliubov transformation in the single frequency mode form and the analytical expressions of key parameters of the fopa , including the photon number gain , the noise figure , and the intensity difference noise of the twin beams .
figure [ fopa](a ) is a typical diagram of the non - degenerate phase - insensitive fopa , which is referred to as the fopa hereinafter for brevity .
when the wavelengths of the pump are properly set , the phase matching of fwm in fiber is satisfied . in the fwm process ,
two pump photons at frequencies @xmath2 and @xmath3 are simultaneously scattered into a pair of signal and idler photons at frequencies @xmath4 and @xmath5 via the @xmath0 nonlinearity , and the energy conservation relation @xmath6 is satisfied . since each generated signal photon
is always accompanied by the birth of an idler photon , the two form correlated twin beams .
a weak input is injected into the fopa from either the signal or the idler field . for convenience , we refer the field with weak injection as the signal , while the other with vacuum input as idler . when the weak signal is amplified by the strong pump field via the fwm process , an idler beam is generated , and the frequency overlap between the amplified signal and generated idler fields is negligible for the non - degenerate fopa .
( a ) a typical diagram of the non - degenerate phase - insensitive fopa .
@xmath7 and @xmath8 , filters ; d1 and d2 , detectors ; @xmath9 and @xmath10 , radio frequency amplifiers ; dds , differential detection system .
( b ) the sketch map of the frequency modes for the single frequency laser pumped fopa .
two pump photons of a single frequency laser uniquely defines the single frequency modes of a pair of signal and idler photons .
( c ) the sketch map of the frequency modes for the pulse pumped fopa . a single frequency signal mode ( dotted arrow )
will be correlated to many idler frequency modes ( enveloping curve in gray ) , and vice versa . since the pump is in multi - frequency modes , the pulsed signal and idler beams
are also in multi - frequency modes ( enveloping curves in black).,scaledwidth=65.0% ]    when the fopa is pumped with a single frequency cw laser , as illustrated in fig .
1(b ) , the two pump photons have the same frequency mode ( @xmath11 ) , a pair of signal and idler photons in the frequency modes @xmath4 and @xmath5 are uniquely defined . in this case , the input - output relation of the single mode fopa is given by the well - known bogoliubov transformation : @xmath12 where @xmath13 and @xmath14 refer the input and output modes , the subscript @xmath15 respectively denote the signal and idler fields , the coefficients @xmath16 and @xmath17 are determined by the parametric gain of fwm and satisfy the relation @xmath18 . for input signal in a coherent state @xmath19 with photon number much greater than 1 , i.e. , @xmath20 it is straightforward to deduce the expressions of all the parameters of the fopa by calculating the quantum average of the corresponding operators over the state @xmath21 , where @xmath22 denotes the vacuum input of the idler field and the superscript @xmath23 denotes the single frequency mode case .
the amplified signal ( generated idler ) beam is then measured by the detector d1(2 ) after passing through the filter @xmath24 to reject the strong pump ( fig .
1(a ) ) . the average photon number @xmath25 and the photon number fluctuation @xmath26 of the signal ( idler ) fields are converted to the dc and ac currents of d1 and d2 , respectively . for identical
d1 and d2 , when the efficiencies of the detector d1(2 ) and filter @xmath24 are ideal , the average photon numbers of the amplified signal and generated idler beams can be written as    [ sig_i ] @xmath27 @xmath28    where @xmath29 is the photon number operator of signal ( idler ) beam . in our theoretical model
, we assume the condition @xmath30 is always satisfied , so the contribution of spontaneous terms @xmath31 in the right hand side of eqs .
( [ sig_i]a ) and ( [ sig_i]b ) are negligible . according to the definition of the photon number gain
 the ratio between the average photon numbers of amplified signal and input signal @xmath32 we have @xmath33 .
moreover , according to the definition of the intensity noise of the amplified signal and generated idler fields , we have    [ sig_delta_i ] @xmath34 @xmath35    after using the snl of signal ( idler ) beam @xmath36 to normalize the intensity noise @xmath37 , we obtain @xmath38 indicating the normalized intensity noise of the amplified signal beam @xmath39 is equal to that of the generated idler beam @xmath40 . notice that they are both above the snl because of amplification .    as for the noise figure ( nf ) of the fopa , according to the definition  @xcite @xmath41 with @xmath42 and @xmath43 respectively denoting the signal to noise ratio of the input and amplified signal beams , the nf of the single mode fopa is given by @xmath44 which shows that the parametric amplification process adds the excess noise to the amplified signal beam and the 3 db quantum limit of nf can be achieved in the high gain limit .
in the ideal case , the photocurrents of d1 and d2 are highly correlated and can be subtracted from each other to yield quantum - noise reduction , i.e. , @xmath45 . in practice , the measured noise reduction is very sensitive to the collection efficiency of twin beams and the response function of differential detection system ( dds ) , consisting of d1 and d2 followed by the radio frequency ( rf ) amplifiers @xmath9 and @xmath10 , and the subtractor ( see fig .
1(a ) ) . for brevity , we assume throughout this paper that the transmission efficiencies of the filters @xmath7 and @xmath8 at their central wavelengths are perfect .
in this situation , the collection efficiencies of the signal and idler beams at the single frequency @xmath4 and @xmath5 are perfect , so we only need to analyze the influence of dds with non - ideal response .
the response function of the dds is described by @xmath46 and @xmath47 , which are determined by the quantum efficiencies and electrical circuit of d1 and d2 , respectively .
here @xmath48 denotes the radio frequency .
the dc response of the dds , @xmath49 , is associated with the quantum efficiency through the relation @xmath50 , where @xmath51 is the electron charge , @xmath52 is the quantum efficiency of d1 ( d2 ) , and @xmath53 is the reduced planck constant ; while the ac response , @xmath54 ( @xmath55 ) , is not only determined by the response of d1 ( d2 ) @xmath56 , but also determined by that of the rf amplifier @xmath57 . for convenience , the ratio @xmath58 between the ac response @xmath46 and @xmath47 , defined as @xmath59 which can be adjusted by changing the electronic gain of the amplifier @xmath60 or @xmath61 .    in practice ,
the quantum efficiencies of d1 and d2 are not perfect , i.e. , @xmath62 and @xmath63 . modeling the non - ideal detectors as the beam splitters , which couple the vacuum mode @xmath64 and @xmath65 to the signal and idler fields , the field operator incident on d1 and d2
are given by    [ eq : sm_operator ] @xmath66 @xmath67    in this situation , the measured photon numbers of the amplified signal and idler beams in eq .
( 3 ) are rewritten as    [ sm_isi ] @xmath68 @xmath69    and the intensity noise of individual signal and idler beams in eq .
( 5 ) is rewritten as    [ deltai_sm ] @xmath70 @xmath71    for the operator of the measured intensity difference of the twin beams @xmath72 the normalized intensity difference noise of the twin beams is expressed as @xmath73 where @xmath74 is the corresponding snl of the detected signal ( idler ) beam and @xmath75 is the quantum correlation term . substituting eqs .
( [ sm_isi])-([deltai_sm ] ) into eq .
( [ eq : electrical_gain ] ) , we obtain @xmath76    to illustrate the factors influencing the value of @xmath77 , we first plot @xmath78 as a function of the gain @xmath79 by assuming the ratio @xmath58 in eq . ( 11 ) is equal to 1 and by varying the efficiencies @xmath80 and @xmath81 , as shown in fig .
when the efficiencies of both d1 and d2 are @xmath82 or @xmath83 , @xmath84 decreases with the increase of @xmath79 ; for a given @xmath79 , @xmath84 obtained for @xmath85 is lower than that for @xmath86 .
when efficiencies of d1 and d2 are @xmath87 and @xmath88 , respectively , @xmath84 is higher than that for @xmath89 in the high gain regime , however , @xmath84 can be slightly lower than that for @xmath90 when the gain is within the regime of @xmath91 .
therefore , for the case of @xmath92 , @xmath84 decreases with @xmath79 and increases with the decrease of @xmath52 ; while for the case of @xmath93 , @xmath84 does not always decrease with @xmath79 : after obtaining the minimum @xmath84 at a certain gain , @xmath84 will start to increase with @xmath79 .
( a ) the noise reduction @xmath84 as a function of the gain @xmath79 under different quantum efficiencies for the cases of ( a ) @xmath94 and ( b ) @xmath95 . , scaledwidth=55.0% ]    eq .
( [ eq : r_sm_eta ] ) shows that for the given values of @xmath80 , @xmath81 and @xmath96 , @xmath84 can be minimized when the ratio @xmath58 takes the optimized value @xmath97 we then plot @xmath84 as a function of @xmath79 by varying the efficiencies @xmath80 and @xmath81 for the case of @xmath95 . as shown in fig .
2(b ) , no matter @xmath80 and @xmath81 are equal or not , @xmath84 always decreases with the increase of @xmath79 .
moreover , for a fixed value of @xmath79 , our calculation indicates that @xmath84 decreases with the increase of @xmath98 and @xmath99 .
we note that it is impossible to obtain the ideal noise reduction @xmath100 for the input signal with photon number @xmath101 , even if the detectors d1 and d2 are perfect .
to clearly illustrate this point , let s look at the expressions of @xmath78 when the condition @xmath102 is fulfilled .
for the case of @xmath103 , eq . ( [ eq : r_sm_eta ] )
is rewritten as @xmath104 while for the case of @xmath94 , eq . ( [ eq : r_sm_eta ] )
is rewritten as @xmath105 both eqs .
( [ r_sm_r ] ) and ( [ r_t_sm ] ) show the intensity difference noise of twin beams is always less than the snl ( 1 or @xmath106 db ) , however , the condition @xmath107 is always satisfied .
moreover , a comparison between eqs .
( [ r_sm_r ] ) and ( [ r_t_sm ] ) indicates that in the high gain limit , the noise reduction can be improved up to 3 db by optimizing the ratio @xmath58 .
when the fopa is pumped by a strong pulsed field having multiple frequency components , such as a mode - locked laser with a pulse width of a few pico - seconds , a single frequency signal mode is correlated to many idler frequency modes , and vice versa ( see the dotted arrow and the corresponding enveloping curve in fig.1 ( c ) ) . in this situation , both amplified signal and generated idler beams are pulsed fields with broad bandwidth , as illustrated by the spectrally broadened curves in fig .
therefore , it is necessary to deal with the multi - frequency modes . in this section ,
we first extend the single mode bogoliubov transformation to the multi - frequency domain .
then , according to the hamiltonian of fwm in fiber , we deduce the four transformation functions of the multi - mode bogoliubov transformation , from which the general expressions of the key parameters of a pulse pumped high gain fopa can be derived
.      an optical field propagating along the fiber ( denoted as the z direction ) can be quantized by using one dimensional approximation if the polarization mode is well defined .
therefore , the pulsed field operator in optical fiber is written as @xcite @xmath108 where the positive frequency field operator @xmath109 and negative frequency operator @xmath110 satisfy the relation @xmath111^\dagger$ ] , and @xmath109 is given by @xmath112 here the annihilation operator at frequency @xmath113 , @xmath114 , satisfies the commutation relation @xmath115=\delta(\omega-\omega')$ ] .
note that @xmath116 , so the distance in propagation is equivalent to a delay in time .    for the pulse
pumped fopa , the strong pump pulses remain classical , but the signal and idler fields are quantized . extending bogoliubov transformation in eq .
( [ single_trans ] ) to the multi - frequency domain and using eq .
( 22 ) , we write the input - output relation as  @xcite :    [ eq : multimode_trans ] @xmath117 @xmath118    where the unitary evolution operator @xmath119 is determined by the hamiltonian of the fopa(see later for the explicit form ) , the footnotes @xmath120 denote the frequency ranges of the signal and idler fields , and the four green functions @xmath121 are referred to as the transformation functions . to ensure the satisfaction of the commutation relations of the operators @xmath122 and @xmath123 , the transformation functions should be constrained by the relations :    [ eq : fwm_opa_unitarity_conditions ] @xmath124    the weak input signal pulses , synchronized with the pump pulses , is ideally in a multi - mode coherent state @xmath125d\omega_s \big\ }   | 0 \rangle,\ ] ] where @xmath126 is the complex amplitude , and the frequency distribution function @xmath127 satisfies the normalization condition @xmath128 so that the photon number of the input signal is @xmath129 , which is the same as eq .
( [ input - sig ] ) .
in general , the bandwidth of the input signal is much smaller than that of the central frequency .
so the quasi - monochromatic approximation applies , and the frequency integral range from @xmath106 to @xmath130 can be treated as from @xmath131 to @xmath130 . for the sake of brevity , all the integral ranges from @xmath131 to @xmath130 will be omitted hereinafter .
in contrast to the single mode fopa , wherein the filters @xmath132 and @xmath133 with ideal transmission efficiency at the central frequencies result in perfect collection efficiency of twin beams , the collection efficiency of the broad band multi - mode twin beams is associated with the spectra of @xmath7 and @xmath8 . in this case , the field operators incident on the detectors d1 and d2 are given by    [ eq : vacuum_mixing ] @xmath134 @xmath135    where the complex function @xmath136 ( @xmath137 ) describes the spectrum of the filter @xmath138 . eq . ( [ eq : vacuum_mixing ] ) implies that the measured quantum noise property of the pulse pumped fopa is not only influenced by the quantum efficiencies of d1 and d2 , but also influenced by the spectral functions @xmath139 and @xmath140 .
for the ideal detectors , we should have @xmath102 , while for the perfect collection of the twin beams , we should have @xmath141 .
similar to the single mode fopa , the key parameters of the pulse pumped fopa can be obtained by calculating the quantum average of the corresponding operators over the state @xmath142 . in the next subsection
, we will derive the analytical forms of @xmath121 , which are the key for deriving the formulas of the operators of twin beams ( see eq .
( [ eq : multimode_trans ] ) ) .      the unitary operator in eq .
( [ eq : multimode_trans ] ) is determined by the hamiltonian of parametric process .
so let s begin with the hamiltonian of the co - polarized pulse pumped fwm in single - mode optical fibers  @xcite @xmath143,\ ] ] where @xmath144 is a constant determined by experimental details and the units of quantized optical fields , @xmath0 is the 3rd - order real nonlinearity , @xmath145 and @xmath146 are the classical pump fields ,    [ eq : ea ] @xmath147 and @xmath148    are the quantized negative - frequency field operator of the signal and idler beams , respectively .
the gaussian shaped strong pump pulses propagating along the fiber can be written as @xcite @xmath149 where @xmath150 , @xmath151 and @xmath152 are the bandwidth , central frequency , and wave vector of the pump , respectively , and @xmath153 is related to the peak power through the relation @xmath154 .
the nonlinear coefficient @xmath155 is expressed as @xmath156 , where @xmath157 denotes the effective mode area and @xmath158 is the speed of light in vacuum . after substituting eqs .
( [ eq : ea ] ) and ( [ eq : pump ] ) into eq .
( [ eq : hamiltonian_0 ] ) and changing @xmath159 in eq .
( [ eq : hamiltonian_0 ] ) to @xmath160 , we carry out the integral over the whole fiber length @xmath161 ( from -l/2 to l/2 ) and arrive at the hamiltonian in the time dependent form @xmath162 where @xmath163 is the phase mismatching term , and the term @xmath164 is originated from the self - phase modulation of pump .    to obtain the formula of the unitary evolution operator @xmath165
, we need to carry out the integration over the time and over all the possible combinations of @xmath2 and @xmath3 within the pump bandwidth .
since the time integral gives rise to the @xmath166 function in @xmath167 to guarantee the energy conservation at single - photon level , we arrive at @xmath168\big\},\ ] ] where the so called two - photon joint spectral function ( jsf ) of fwm has the form of  @xcite @xmath169 which is determined by the pump envelop function @xmath170 and phase matching function @xmath171 , and is referred to as the probability amplitude of simultaneously finding a pair of signal and idler photons within the frequency range of @xmath172 and @xmath173 , respectively .
the coefficient @xmath174 in eq .
( [ eq : nondegererate_evolution_operator ] ) determines the gain of fwm , and the coefficient @xmath175 in eq .
( [ jsf_def ] ) is a constant used to ensure the satisfaction of the normalization condition @xmath176 . according to eq .
( [ jsf_def ] ) , the jsf is generally asymmetry since @xmath177 is frequency dependent . therefore
, @xmath178 is generally not satisfied .    to obtain the formula of the transformation functions , we rewrite the input - output relation by substituting eqs .
( [ eq : nondegererate_evolution_operator ] ) into eqs.([eq : multimode_trans ] ) and using baker - hausdorff lemma : @xmath179\nonumber\\   & & \ \
+ \frac{1}{2!}[\hat{a}_j(\omega_j'),[\hat{a}_j(\omega_j'),\hat{b}]]+\frac{1}{3!}[\hat{a}_j(\omega_j'),[\hat{a}_j(\omega_j'),[\hat{a}_j(\omega_j'),\hat{b}]]]+ ... \ ,
~~ \label{lieseries }     \end{aligned}\ ] ] where the operator @xmath180 is expressed as @xmath181.\ ] ] then the transformation functions can be written in the form of an infinite series    [ eq : trans_intergral ] @xmath182\nonumber\\ & & [ \psi^*(\omega_2,\omega_1)\psi^*(\omega_4,\omega_3)\psi^*(\omega_6,\omega_5 ) \cdot\cdot\cdot \psi^*(\omega_s',\omega_{2n-1 } ) ] \}\end{aligned}\ ] ] @xmath183\nonumber\\&&[\psi(\omega_s,\omega_1)\psi(\omega_2,\omega_3)\psi(\omega_4,\omega_5 ) \cdot\cdot\cdot \psi(\omega_{2n},\omega_i ' ) ] \}\end{aligned}\ ] ] @xmath184\nonumber\\ & & [ \psi^*(\omega_1,\omega_2)\psi^*(\omega_3,\omega_4)\psi^*(\omega_5,\omega_6 ) \cdot\cdot\cdot \psi^*(\omega_{2n-1},\omega_i ' ) ] \}\end{aligned}\ ] ] @xmath185\nonumber \\ & & [ \psi(\omega_1,\omega_i)\psi(\omega_3,\omega_2)\psi(\omega_5,\omega_4 ) \cdot\cdot\cdot \psi(\omega_s',\omega_{2n } ) ] \}.\end{aligned}\ ] ]    eq .
( [ eq : trans_intergral ] ) clearly shows that the transformation functions satisfy the relations    [ eq : signal_idler_correlation ] @xmath186 @xmath187 @xmath188    indicating that there is a correlation between the signal and idler beams .
it is worth noting that bogoliubov transformation in multi - frequency modes ( eq ( [ eq : multimode_trans ] ) ) with the transformation function shown in eq .
( [ eq : trans_intergral ] ) is also suitable for describing the single mode fopa .
for example , if the jsf in eq .
( [ jsf_def ] ) takes the limit of a single frequency laser pumped case , i.e. , @xmath189 , with @xmath190 and @xmath191 denoting the single frequencies of the signal and idler beams , the input - output relation of eq ( [ eq : multimode_trans ] ) would transform into @xmath192 with @xmath193 and @xmath194 , which is exactly the same as eq .
( [ single_trans ] ) . hence , eqs .
( [ eq : multimode_trans ] ) and ( [ eq : trans_intergral ] ) can be viewed as a general description of the input - output relation of an fopa .      for the measurements of the signal and idler fields ,
in addition to taking the quantum efficiencies , @xmath80 and @xmath81 , and the spectrum of filter @xmath24 into account , we need to consider the response times of the detectors d1 and d2 as well .
when the response time of each detector is much longer than the pump pulse duration , but is much shorter than the period of two adjacent pump pulses , the average photon number of the amplified signal ( generated idler ) beam per pulse can be expressed as @xmath195 where @xmath196 with @xmath197 is the photon number operator of the pulsed field . the time integral ranges in eqs .
( [ i_j ] ) and ( [ field - operator ] ) are omitted , because they can be treated as from @xmath131 to @xmath130 .    with the filter operator incident on d1 and d2 ( eq .
( [ eq : vacuum_mixing ] ) ) , we substitute eqs .
( [ eq : multimode_trans ] ) , ( [ eq : coh_state ] ) , ( [ field - operator ] ) into eq .
( [ i_j ] ) , and obtain the expressions of the detected photon number    [ eq : signal_idler_average photon number ] @xmath198 @xmath199    where    [ spectra ] @xmath200 and @xmath201    are the power spectra of the signal and idler beams . the first terms in the right hand side of eqs . (
[ eq : signal_idler_average photon number]a ) and ( [ eq : signal_idler_average photon number]b ) are originated from the stimulated emission of input signal , while the second terms are from the spontaneous emission . since we have assumed that the photon number of the weak input signal per pulse is much greater than 1 , the spontaneous emission terms are negligible .    using eq .
( [ eq : is ] ) and the definition of photon number gain in eq .
( [ def_g ] ) , we find the expression of the photon number gain of the pulse pumped fopa is @xmath202 for the case of @xmath203 and @xmath204 .    to compute the intensity noise of the amplified signal ( generated idler ) beam @xmath205 we substitute eqs .
( [ eq : multimode_trans ] ) , ( [ eq : coh_state ] ) , ( [ eq : vacuum_mixing ] ) and ( [ i_j])-([field - operator ] ) into eq .
( [ eq : fluctuation_general ] ) , and arrive at the detailed expression    [ eq : photon_number_fluctuation_fwm_opa ] @xmath206 @xmath207    with    [ eq45 ] @xmath208
@xmath209 @xmath210 @xmath211    and    [ eq46 ] @xmath212\end{aligned}\ ] ] @xmath213.\end{aligned}\ ] ]    accordingly , the normalized intensity noise of signal ( idler ) beam is @xmath214 moreover , we have @xmath215 by substituting eqs .
( [ eq : signal_idler_average photon number ] ) and ( [ rsandri ] ) into eq .
( [ snr - out ] ) . according to the definition in eq .
( [ eq : nf1 ] ) , we obtain the formula of noise figure @xmath216    eqs .
( [ eq : photon_number_fluctuation_fwm_opa])-([rsandri ] ) indicate that the intensity noises of the individual beams @xmath217 are not only related to detection efficiencies , @xmath80 and @xmath81 , but also to the spectral functions of filters , @xmath139 and @xmath140 .
when @xmath80 and @xmath81 are very low , i.e. , @xmath218 ( @xmath137 ) , the intensity fluctuation @xmath26 in eq .
( [ eq : photon_number_fluctuation_fwm_opa ] ) is dominated by the term @xmath219 ( @xmath137 ) , while the contribution of the term @xmath220 ( @xmath221 ) is negligible . on the other hand ,
when the bandwidth of the filters @xmath7 and @xmath8 are much narrower than that of input signal , @xmath26 is also dominated by the term @xmath219 ( @xmath137 ) .
since the non - ideal detector induces the vacuum noise to the individual signal or idler field , the former is easy to understand .
however , the understanding of the latter is not so straightforward .
we think this is because , when the signal and idler beams are in multi - frequency mode , the filter @xmath7 ( @xmath8 ) with bandwidth narrower than that of the amplified signal ( generated idler ) beam can also be viewed as a loss , which introduce vacuum noise to the signal ( idler ) beam as well . hence , the term @xmath219 ( @xmath137 ) in eqs .
( [ eq : h_vac_s ] ) and ( [ eq : h_vac_i ] ) is originated from the loss induced vacuum .
if the detection and collection efficiencies of the signal and idler beam are ideal , we would have @xmath222 ( @xmath137 ) .    finally , we derive the expression of the intensity difference noise of the twin beams by using the definition @xmath223 where the operator @xmath224 is the photon number difference with a weight factor of ac response ratio @xmath58 . using eqs .
( [ eq : multimode_trans ] ) , ( [ eq : vacuum_mixing ] ) , ( [ i_j])-([field - operator ] ) and ( [ eq : fluctuation_general ] ) , eq .
( [ difference - fluctuation_general ] ) is transformed into @xmath225 - ( i_s - ri_i)^2\nonumber\\ = & & \delta i_s^2+r^2\delta i_i^2 - 2r(h_{si1}+h_{si2})|\alpha|^2 ,   \end{aligned}\ ] ] where @xmath226 and @xmath227 are positive terms originated from the quantum correlation between signal and idler twin beams .
it is obvious that @xmath228 is less than the sum of intensity noise of individual beams @xmath229 due to the correlation of twin beams .
consequently , similar to eq .
( [ eq : electrical_gain ] ) , the general expression of the normalized intensity difference noise for pulsed twin beams @xmath230 |\alpha|^2 } { i_s + r^2 i_i } \ ] ] is obtained by substituting eqs .
( [ eq : photon_number_fluctuation_fwm_opa ] ) and ( [ eq : h_si1])-([eq : h_si2 ] ) into eq .
( [ eq : it_h ] ) .
in addition , @xmath231 can be further minimized if the ac response ratio @xmath58 takes the optimized value @xmath232 ^ 2 } } { 2i_i(h_{si1}+h_{si2})}\nonumber\\ & & + \frac{i_i(h_{s1}+h_{s2}+h_{vs})-i_s(h_{i1}+h_{i2}+h{vi})}{2i_i(h_{si1}+h_{si2})}.\end{aligned}\ ] ]    from the formulas of the photon number gain , noise figure and intensity difference noise in eqs .
( [ gdef ] ) , ( [ nf2 ] ) and ( [ eq : r_t2 ] ) , one sees that the key parameters of the fopa are determined by the transformation functions in eq .
( [ eq : trans_intergral ] ) , which highly depend on the jsf in eq .
( [ jsf_def ] ) . in the following sections , we will focus on studying the noise characteristics of the pulse pumped fopas with two typical kinds of jsf .
we first study the fopa with a spectrally factorable jsf . under this condition ,
both the signal and idler modes are in single temporal mode , so we expect to recover the results for single mode fopa in sec .
[ single_mode ] . for the sake of brevity , in this section
, we assume the quantum efficiencies of d1 and d2 are perfect , and the collection efficiency of the twin beams is ideal .
a factorable jsf , obtained by properly regulating the pump pulses and tailoring the dispersion of optical fibers , is written as  @xcite @xmath233 where @xmath234 and @xmath235 with the normalization condition of @xmath236 and @xmath237 are the probability amplitude of finding a pairs of signal and idler photons within the frequency range of @xmath238 and @xmath239 , respectively .
in fact , looking at eq .
( [ eq : factotized jsf ] ) from another point of view , we find @xmath234 and @xmath235 can also be viewed as the gain spectra in signal and idler fields , respectively .    with the jsf in eq .
( [ eq : factotized jsf ] ) , the transformation functions in eq .
( [ eq : trans_intergral ] ) are simplified to    [ eq : fac_trans_func ] @xmath240    note that in this section , we will use the superscript @xmath241 to represent the result of an fopa with the factorable jsf .    with the jsf in eq .
( [ eq : factotized jsf ] ) , the average photon numbers of the amplified signal and generated idler beams per pulse in eq.([eq : signal_idler_average photon number ] ) are simplified to    [ eq : factorized_signal_avg_ph_number ] @xmath242    where the coefficient @xmath243 is determined by the matching between the spectra of signal injection @xmath244 ( see eq .
( [ eq : coh_state ] ) ) and that of the gain in signal band @xmath245
. therefore , we obtain the photon number gain of the fopa @xmath246 by substituting eqs .
( [ eq : factorized_signal_avg_ph_number ] ) and ( [ input - sig ] ) into eq .
( [ def_g ] ) .    to show the factors influencing the noise performance of the fopa with factorable jsf , we first substitute the simplified transformation function in eq .
( [ eq : fac_trans_func ] ) into eqs.([eq45])-([eq46 ] ) and ( [ eq : h_si1])-([eq : h_si2 ] ) . for the case of @xmath247 , the expressions of the terms in eqs .
( [ eq45])-([eq46 ] ) and ( [ eq : h_si1])-([eq : h_si2 ] ) , which determine @xmath248 ( @xmath137 ) in eq .
( [ eq : photon_number_fluctuation_fwm_opa ] ) and @xmath249 in eq .
( [ eq : it_h ] ) , are simplified to @xmath250 @xmath251 @xmath252 @xmath253 consequently , the normalized intensity noise of individual signal and idler beams are @xmath254 and @xmath255 respectively . moreover , according to eq .
( [ eq : factorized_signal_avg_ph_number ] ) , the noise figure in eq .
( [ nf2 ] ) and the normalized intensity difference noise of twin beams in eq .
( [ eq : r_t2 ] ) can be respectively simplified to @xmath256 and @xmath257    equations ( [ rs_fac])-([eq : rt_fac ] ) show that except for the normalized intensity noise of idler beam @xmath258 , the remaining three parameters , @xmath259 , @xmath260 and @xmath261 , are associated with the matching coefficient @xmath262 .
for the case of @xmath263 , if we respectively replace the terms @xmath264 and @xmath265 with @xmath266 and @xmath31 , the analytical expressions of the key parameters ( eqs .
( [ rs_fac])-([eq : rt_fac ] ) ) will be the same as those of the single frequency mode fopa in sec .
[ single_mode ] .
the results indicate that our multi - mode theory of the pulse pumped fopa are valid .
an marked difference between the fopas described by single temporal mode and single frequency mode is the spectra of signal and idler beams .
the spectra of the amplified signal and generated idler beams obtained by substituting eq .
( [ eq : fac_trans_func ] ) into eq.([spectra ] ) are written as    [ spectrum_fa ] @xmath267 and @xmath268    respectively . from eq .
( [ spectrum_fa ] ) , one sees that the spectrum of generated idler is always the same as the spectrum of the gain in idler field , however , different from the fopa pumped with a single frequency laser , the spectra of the injected signal and amplified signal beams are different unless the condition @xmath263 is fulfilled .    to further understand the difference between the fopas described by single temporal mode and single frequency mode
, we then compare the noise performance of the two cases .
we find that for a given value of @xmath269 ( or @xmath270 ) , the noise figure and the intensity difference noise of the twin beams in eqs .
( [ eq : nf_fac ] ) and ( [ eq : rt_fac ] ) are worse than that predicted by eqs ( [ nf_s_sm ] ) and ( [ r_t_sm ] ) unless @xmath263 . in practice
, it is very difficult to maintain the spectral matching condition @xmath271 in the high gain regime due to the self - phase modulation and cross - phase modulation induced spectral broadening  @xcite .
hence , it is very challenging to experimentally realize a pulse pumped fopa capable of being described by a single mode theory .
for the conventional optical communication systems , an fopa with an extremely broad gain bandwidth is desirable .
when the gain bandwidth of the fopa is much broader than the bandwidths of the pump and input signal , the jsf can be simplified by assuming the perfect phase matching condition @xmath272 . in this case , the jsf in eq .
( [ jsf_def ] ) is rewritten as  @xcite @xmath273 which is obviously non - factorable . in this section ,
we focus on numerical investigation of the quantum noise performance of this kind of fopa .    with the jsf in eq .
( [ eq : jsf_pm ] ) , the transformation functions in eq.([eq : trans_intergral ] ) are reformulated as    [ trans_unfac ] @xmath274    where @xmath275 , and the superscript @xmath276 represents the result for the fopa with broad gain bandwidth . for the sake of convenience , the frequency variables @xmath277 and @xmath278 used in previous sections will be replaced with @xmath113 and @xmath279 hereinafter .
we first study the factors influencing the key parameters of the fopa by assuming the detectors are perfect , and collection efficiency of the twin beams is ideal , i.e. , @xmath102 and @xmath280 . without loss of the generality , we assume the spectral function of the pulsed input signal can be described by the gaussian function @xmath281 where @xmath282 denotes the bandwidth . from eqs .
( [ trans_unfac ] ) and ( [ eq : guassian_input_signal_spectrum ] ) , we can rewrite the power spectra of the amplified signal and generated idler beams in eq .
( [ spectra ] ) as    [ spectra_unfac ] @xmath283    for all the terms of the infinite series in eqs .
( [ spectra_unfac]a ) and ( [ spectra_unfac]b ) , whose order corresponds to @xmath284 or @xmath285 , their bandwidths are larger than that of the input power spectrum @xmath286 .
moreover , since the bandwidth of the term in infinite series increases with the order determined by the integers @xmath287 and @xmath288 , the bandwidths of both signal and idler twin beams increase with @xmath289 .    using eq .
( [ spectra_unfac ] ) , the average photon numbers of amplified signal field and generated idler beams in eq .
( [ eq : signal_idler_average photon number ] ) can be rewritten as    [ i_unfac ] @xmath290    where @xmath291 describes the ratio of the pump bandwidth to the input signal bandwidth .
for all the numerically simulated results presented hereinafter , the coefficient @xmath289 is within the range of @xmath292 , and the series in eqs .
( [ trans_unfac ] ) , ( [ spectra_unfac ] ) and ( [ i_unfac ] ) are calculated up to the order of @xmath293 and @xmath294 . to ensure the correctness of this truncation approximation
, we numerically verify the general relation @xmath295 ( @xmath129 ) , which is the inherent nature of signal and idler twin beams originated from the energy conservation of fwm . after calculating the series in eq .
( [ i_unfac ] ) to the order of @xmath293 and @xmath294 , we find that the difference between the calculated result of @xmath296 and  1 " is less than @xmath297 , showing the validity of the approximation .
figures [ spectrum](a ) and [ spectrum](b ) show the normalized power spectra of the amplified signal and generated idler beams .
in the plots , eqs .
( [ spectra_unfac]a ) and ( [ spectra_unfac]b ) are calculated by varying the value of @xmath289 under the condition @xmath298 . as a comparison
, we also show the power spectrum of input signal @xmath299 in fig .
[ spectrum](a ) .
it is obvious that the bandwidths of both signal and idler twin beams are broader than that of the input signal and they increase with @xmath289 , showing the spectrum broadening effect .
moreover , in the low gain regime , the spectrum of amplified signal is mainly determined by the input signal , so the bandwidth of generated idler beam , which is the convolution of the spectra of the input signal and pump , is greater than that of the signal beam . however , in the high gain regime , the bandwidths of signal and idler twin beams are about equal because both the spectra of signal and idler beams are the convolution of the spectra of pump and its counterparts .
the normalized power spectra of ( a ) amplifed signal and ( b ) generated idler beams at different @xmath289 under the condition @xmath298 . as a comparison ,
the power spetrum of input signal @xmath286 is plotted in ( a ) as well.,scaledwidth=48.0% ]    figure [ fig : photon_number_gain_no_filtered ] plots the dependence of the photon number gain @xmath300 upon the pump power for the different ratio @xmath291 . in the calculation , @xmath301 is changed by varying the peak pump power @xmath302 .
one sees that at a fixed peak power @xmath302 , @xmath270 increases with the decrease of @xmath303 because the overlap of temporal mode , which is required for maximizing the gain of fwm , is improved by decreasing the bandwidth of pump .
when @xmath303 is less than @xmath304 , further decreasing the ratio to @xmath305 does not result in an obvious increase in @xmath270 , because no more space is left for improving the temporal mode overlapping .
the result in fig .
[ fig : photon_number_gain_no_filtered ] implies the pump with a longer pulse duration ( narrower bandwidth ) gives a better mode overlapping .
however , it is worth noting that the pump with a short pulse duration helps to obtain a high gain under the condition of low average pump power , which is of practical importance .
photon number gain @xmath270 versus the peak pump power for the different ratio @xmath291 . in the calculation , @xmath272 , @xmath306 , the parameter @xmath307 , and @xmath308.,scaledwidth=33.0% ]    to investigate the quantum noise performance of the fopa with broad gain bandwidth , we first deduce the formulas of terms @xmath220 ( @xmath137 , @xmath309 ) and @xmath219 ( @xmath137 ) , which determine @xmath248 ( @xmath137 ) in eq .
( [ eq : photon_number_fluctuation_fwm_opa ] ) . by substituting eqs .
( [ trans_unfac ] ) and ( [ eq : guassian_input_signal_spectrum ] ) into eqs .
( [ eq45])-([eq46 ] ) , we obtain @xmath310 @xmath311 @xmath312 @xmath313    with the above , we then numerically study the dependence of the normalized intensity noise of individual beams @xmath314 ( @xmath137 ) upon the ratio @xmath303 .
figure [ fig : relative_intensity_noise_no_filtered](a ) plots @xmath315 and @xmath316 versus @xmath270 for the cases of @xmath317 , @xmath318 , and @xmath319 , respectively . from fig .
[ fig : relative_intensity_noise_no_filtered](a ) , one sees that for a fixed value of @xmath270 ( @xmath303 ) , both @xmath315 and @xmath316 increase with @xmath303 ( @xmath270 ) .
moreover , for the given values of @xmath303 and @xmath270 , @xmath316 is greater than @xmath315 .
this is different from the results of @xmath320 ( @xmath137 ) for a single mode fopa [ see eq .
( [ r_cw ] ) ] , but has been experimentally confirmed in references @xcite and @xcite .
( a ) the normalized intensity noise of the individual signal/ idler beams @xmath321 ( @xmath137 ) and ( b ) noise figure nf versus the photon number gain @xmath270 for the different ratio @xmath291 . in this calculation , @xmath272 , and @xmath247.,scaledwidth=33.0% ]
we also compute the noise figure @xmath322 in eq .
( [ nf2 ] ) . figure
[ fig : relative_intensity_noise_no_filtered](b ) plots @xmath322 versus @xmath270 for @xmath317 , @xmath318 , and @xmath319 , respectively .
one sees that for a fixed ratio @xmath303 , @xmath322 increases with @xmath270 in the low gain regime , and approaches a constant in the high gain regime . while for a fixed @xmath270 , the value of nf increases with the ratio @xmath303 .
it is obvious that for the pulsed pump case with @xmath323 , @xmath322 in the high gain regime is always greater than the well known 3db - limit of a single mode fopa .
the normalized intensity difference noise of the twin beams @xmath231 versus photon number gain @xmath270 for @xmath324 , @xmath325 and @xmath326 when @xmath58 is set to ( a ) @xmath94 and ( b ) @xmath327 , respectively.,scaledwidth=55.0% ]     the optimized value @xmath328 and the photon number ratio @xmath329 as a function of @xmath270 . in this calculation
, the phase mismatching term of the fopa is @xmath272 , and the detection efficiencies and the filters satisfy @xmath247.,scaledwidth=40.0% ]    to study the factors influencing the normalized intensity difference noise of the twin beams @xmath231 in eq .
( [ eq : r_t2 ] ) , we need to calculate @xmath330 and @xmath331 as well . by substituting eqs .
( [ trans_unfac ] ) and ( [ eq : guassian_input_signal_spectrum ] ) into eqs .
( [ eq : h_si1 ] ) and ( [ eq : h_si2 ] ) , we find the following relation @xmath332 in figs . [
fig : twin_beams_no_filtered](a ) and [ fig : twin_beams_no_filtered](b ) , we plot @xmath231 as a function of @xmath270 for different values of @xmath303 when the ac response ratio @xmath58 takes the values of @xmath94 and @xmath327 , respectively .
[ fig : twin_beams_no_filtered](a ) shows @xmath231 is independent on @xmath303 and always decreases with the increase of @xmath270 ; while fig .
[ fig : twin_beams_no_filtered](b ) demonstrates that @xmath231 still decreases with the increase of @xmath270 , but @xmath231 can be further decreased and depends upon the ratio @xmath303 . for a fixed gain @xmath270
, @xmath231 decreases with @xmath303 , which is different from the case of @xmath94 .
moreover , when @xmath303 is very small , for example , @xmath324 , the difference between the @xmath231-values with @xmath94 and @xmath327 is about 3 db in the high gain limit ; when @xmath303 is very big , for example , @xmath326 , @xmath231-values with @xmath94 and @xmath327 are almost the same .    the results in figs .
[ fig : twin_beams_no_filtered](a ) and [ fig : twin_beams_no_filtered](b ) not only indicate the noise reduction of @xmath231 can be improved by properly reducing the ratio @xmath303 , but also imply that @xmath328 might depend upon @xmath303 and @xmath328 might be very close to 1 when the value of @xmath303 is greater than 10 .
this is shown in fig .
[ fig6 ] , in which we plot @xmath328 as a function of @xmath270 for @xmath324 , @xmath325 and @xmath326 .
in addition , we also plot the photon number ratio @xmath329 versus @xmath270 in fig .
one sees that @xmath328 is always within the range of @xmath329 and @xmath333 . for a fixed value of @xmath270
, @xmath328 decreases with the increase of @xmath303 .
as we expected , @xmath328 obtained for the case of @xmath326 is about 1 if @xmath270 is greater than @xmath334 .
moreover , in the high gain limit , both @xmath328 and @xmath329 approach @xmath333 , which is irrelevant to @xmath303 .      in practice , the ideal condition @xmath247 can never be fulfilled , the measured intensity difference noise @xmath231 is closely related to the loss of the fopa system .
so we need to study the influence of the quantum efficiencies ( @xmath80 and @xmath81 ) and collection efficiency of twin beams . in this subsection
, we assume the filters placed at the output port of the fopa ( see fig .
1 ) have a gaussian shaped spectrum @xmath335 where @xmath336 and @xmath337 denote the central frequency and bandwidth of the filter @xmath338 @xmath339 , respectively .    using the transformation functions , input signal spectrum and the spectrum of @xmath338 @xmath339 in eqs .
( [ trans_unfac ] ) , ( [ eq : guassian_input_signal_spectrum ] ) and ( [ eq : filter_shape_guassian ] ) , respectively , the measured photon numbers of the amplified signal and generated idler beams in eq .
( [ eq : signal_idler_average photon number ] ) are rewritten as    [ eq : perfect_phase_matching_filtered_signal_ave_pho_num ] @xmath340\   \}^{-1/2 } \\ & & i'^{(b)}_i=\eta_i \sum_{n_1=0}^{\infty } \sum_{n_2=0}^{\infty } \frac{|g'|^{2n_1 + 2n_2 + 2}}{(2n_1 + 1)!(2n_2 + 1 ) ! } \nonumber\\&&\ \ \ \ \ \{\ s^2[(4n_1 + 2)p^2 + 1][(4n_2 + 2)p^2 + 1 ] \nonumber\\&&\ \ \ \ \ \ \ + [ 2(n_1+n_2 + 1)p^2 + 1]\   \}^{-1/2},\end{aligned}\ ] ]      for a fixed value of @xmath351 , if the bandwidth of filters @xmath24 is much narrower than that of the input signal , namely , @xmath352 , we have @xmath353 . in this case , the collection efficiency of the twin beams is very low due to the multi - mode nature of the signal and idler fields  @xcite .
since the narrow band filters @xmath24 introduce the vacuum noise into the individual signal and idler beams , for a given @xmath289 , the value of @xmath354 @xmath339 in eqs ( [ hvacs_f])-([hvaci_f ] ) approaches to @xmath355 and is dominant over the terms @xmath356 ( @xmath357 ) and @xmath358 ( @xmath357 ) ( eqs . ( [ hs1_f])-([hi2_f ] ) and ( [ hsi1_f])-([hsi2_f ] ) ) .
hence , the normalized intensity noise and normalized intensity difference noise approach to the snl , i.e. @xmath359 and @xmath360 .    in general , the collection efficiency of the twin beams is not only associated with the ratio @xmath361 , but also depends on the ratio @xmath303 . to reduce the detrimental effect of the collection efficiency on @xmath231 , both @xmath361 and @xmath303 should be small enough unless the condition @xmath362 or @xmath305 is satisfied  @xcite .
to illustrate this point , we plot the contour of the normalized intensity difference noise @xmath231 as a function of @xmath303 and @xmath361 in fig .
[ contour](a ) . in the calculation , we assume @xmath363 , @xmath94 , and the quantum efficiencies of d1 and d2 are perfect . for each combination of @xmath361 and @xmath303 , we compute the terms in eqs .
( [ eq : perfect_phase_matching_filtered_signal_ave_pho_num])-([hsi2_f ] ) and substitute the results into eq .
( [ eq : r_t2 ] ) .
[ contour](a ) shows that @xmath231 depends on both @xmath361 and @xmath303 .
we notice that with the increase of @xmath303 , the dependence of @xmath231 upon @xmath361 becomes stronger .
when @xmath303 is less than @xmath364 , at a fixed @xmath303 , @xmath231 increases with @xmath361 , but the noise reduction condition of @xmath365 is always achievable .
however , when @xmath303 is greater than @xmath364 , to obtain the noise reduction of @xmath365 , the value of @xmath361 should be smaller than a certain value , which decreases with the increase of @xmath303 .
in real experiment , it is necessary to use the filter @xmath24 with a certain bandwidth to prevent the strong pump background from reaching the detectors d1 and d2 .
so , to improve the noise reduction of twin beams , the pump with a pulse duration longer than that of the input signal , i.e. , @xmath366 , is desirable .
contour plot of ( a ) the normalized intensity difference noise @xmath231 , ( b ) @xmath367 , and ( c ) @xmath368 as a function of the ratios @xmath303 and @xmath361 . in this calculation
, @xmath369 , @xmath363 , and @xmath94 , and @xmath370 , scaledwidth=85.0% ]    figure [ contour](a ) also shows that the contour @xmath371 divides the plot into two zones according to @xmath372 or @xmath373 .
for the former case , there is a one - to - one correspondence between the combination of @xmath361 , @xmath303 , and @xmath231 ; while for the latter case , each @xmath231 maps two values of @xmath361 at a fixed value of @xmath303 : one is very close to @xmath106 , the other is larger . to better understand fig .
[ contour](a ) , we also plot the contours of the sums of the positive terms and negative terms in eq .
( [ eq : r_t2 ] ) in figs .
[ contour](b ) and [ contour](c ) , which respectively illustrate the influence of collection efficiency upon the intensity noise of individual beams and upon the correlation of twin beams . for both figs .
[ contour](b ) and [ contour](c ) , the contour decreases with the increase of @xmath361 ( @xmath303 ) at a fixed value of @xmath303 ( @xmath361 ) . however , since the correlation of twin beams highly relies on the collection efficiency , the dependence of the contour on @xmath361 and @xmath303 in fig .
[ contour](c ) is stronger than that in fig .
[ contour](b ) .
in particular , when @xmath303 increases , the value of @xmath374 in fig .
[ contour](c ) rapidly decreases with the increase of @xmath361 , which is responsible for the multivalued mappings in the zone of @xmath373 ( see fig .
[ contour](a ) ) .    the normalized intensity difference noise of the twin beams @xmath231 versus @xmath270 for ( a ) @xmath375 , ( b ) @xmath325 , respectively .
for each setting of @xmath303 , @xmath361 in eq.([eq : perfect_phase_matching_filtered_signal_ave_pho_num ] ) and ( [ hs1_f])-([hvaci_f ] ) is @xmath376 , @xmath377 , and @xmath378 , respectively . in this calculation , @xmath369 .
, scaledwidth=66.0% ]    to further demonstrate the dependence of the noise reduction of @xmath231 upon the photon number gain @xmath270 , we plot @xmath231 as a function of @xmath270 for @xmath376 , @xmath377 and @xmath378 . in figs .
[ fig : twin_beams_filter_bandwidth](a ) and [ fig : twin_beams_filter_bandwidth](b ) , we assume @xmath102 and @xmath94 , and @xmath303 is fixed at @xmath324 and @xmath325 , respectively . for each combination of @xmath361 and @xmath303 , we compute the terms in eqs .
( [ eq : perfect_phase_matching_filtered_signal_ave_pho_num])-([hvaci_f ] ) by changing @xmath289 ( @xmath270 ) and substitute the results into eq.([eq : r_t2 ] ) .
it is obvious that for @xmath303 and @xmath361 with fixed values , @xmath231 decreases with the increase of @xmath270 .
moreover , at a fixed value of @xmath270 and @xmath361 ( @xmath303 ) , @xmath231 increases with the increase of @xmath303 ( @xmath361 ) due to the decreased collection efficiency of twin beams .    the normalized intensity difference noise of the twin beams @xmath231 versus g for different setting of @xmath80 and @xmath81 .
for the plots in ( a ) , @xmath379 ; while for the plots in ( b ) , r is optimized for minimizing the measured @xmath231 . in this calculation , @xmath318 and @xmath380.,scaledwidth=66.0% ]    having understood the influence of the collection efficiency of twin beams on @xmath231 , we then study the influence of the detection efficiency on @xmath231 when the values of @xmath303 and @xmath361 are fixed . without loss of generality , we assume @xmath325 and @xmath377 .
when the ac response ratio is set to @xmath94 , we plot @xmath231 as a function of @xmath270 by varying @xmath80 and @xmath81 , as shown in fig .
[ fig : twin_beams_filtered_effeciency](a ) .
it is obvious that for the case of @xmath92 , @xmath231 decreases with the increase of @xmath52 at a fixed value of @xmath270 .
however , for the case of @xmath93 , @xmath231 is not a monotonically decreasing function of @xmath270 and @xmath52 . depending on @xmath270 ,
the correlation of photon currents originated from the quantum correlation of the signal and idler beams might become weak or strong due to the unbalanced detection efficiencies . as shown in fig .
[ fig : twin_beams_filtered_effeciency](a ) , @xmath231 for the case of @xmath87 and @xmath88 is lower than that for @xmath90 when @xmath270 is less than 18 ; however , @xmath231 becomes even higher than that for @xmath89 when @xmath270 is greater than 30 .
making comparison between fig .
[ fig : twin_beams_filtered_effeciency](a ) and fig .
2(a ) , which are obtained for fopa with non - factorable jsf and with single frequency pump , respectively , we find there are two differences : ( i ) for @xmath270 , @xmath80 and @xmath81 with certain values , @xmath231 in fig . [
fig : twin_beams_filtered_effeciency](a ) is larger than that in fig .
2(a ) because the collection efficiency of the pulsed twin beams is smaller ; ( ii ) the difference between the curve for @xmath87 and @xmath88 and that for @xmath90 ( @xmath89 ) in fig . [
fig : twin_beams_filtered_effeciency](a ) is more apparent than that in fig .
therefore , for the pulse pumped fopa , the measured @xmath231 can be minimized by properly adjusting the efficiency @xmath80 or @xmath81  @xcite . however , this method will introduce extra loss , and the minimized @xmath231 can only be obtained for @xmath270 with a specified value .
if the ratio @xmath58 is adjustable , we can minimize @xmath231 by optimizing the value of @xmath58 without introducing the vacuum noise . in this case , as shown in fig .
[ fig : twin_beams_filtered_effeciency](b ) , for fixed values of @xmath80 and @xmath81 , @xmath231 always decreases with the increase of @xmath270 ; while for a fixed value of @xmath270 , @xmath231 always decreases with the increase of @xmath98 and @xmath99 .
hence , for both the signal and idler channels , the higher the efficiencies are , the better the noise reduction is .      finally , we analyze the influence of the excess noise of input signal on @xmath231 , since the noise of input signal pulses originated from the mode - locked fiber lasers is often higher than the snl  @xcite .
when the excess noise is included , the input signal can be viewed as a mixture of coherent states , whose density operator and average photon number are written as @xmath381 and @xmath382 respectively , where @xmath383 is a complex random variable and @xmath384 is the classical probability density .
consequently , the formulas of the average photon number of amplified signal and generated idler beams , @xmath385 and @xmath386 , can be obtained by replacing the term @xmath387 in eq .
( [ eq : signal_idler_average photon number ] ) with @xmath388 .
we then calculate the quantum noise @xmath26 in eq .
( [ eq : fluctuation_general ] ) and @xmath389 in eq .
( [ difference - fluctuation_general ] ) by taking the quantum average over the density operator @xmath390 .
after some algebra , we get the normalized intensity noise of the amplified signal ( generated idler ) beam @xmath391 and the normalized intensity difference noise of twin beams @xmath392 where @xmath393 , and the inequality @xmath394 holds for the input signal with excess noise . in eqs .
( [ rj - excess ] ) and ( [ rt - excess ] ) , the first terms @xmath395 and @xmath396 in the right hand sides are the corresponding noise for input signal in a pure coherent state with @xmath397 ( see eqs .
( [ rsandri ] ) and ( [ eq : r_t2 ] ) ) ; while the second terms in the right hand sides are originated from the excess noise of input signal .
equations ( [ rj - excess ] ) and ( [ rt - excess ] ) indicate that the side effect of the excess noise can be eliminated under the condition @xmath398 , and the experiment in ref .
@xcite has verified this for the case of @xmath94 . in the low gain regime ,
the value of @xmath398 is different from @xmath328 for minimizing @xmath231 ( see fig . [
fig : twin_beams_no_filtered](b ) ) , so the excess noise the input signal is deleterious for the noise reduction of @xmath231 .
however , in the high gain limit , we have @xmath399 . in this case , eliminating the side effect of excess signal noise and minimizing @xmath231 can be achieved simultaneously .
therefore , a high gain fopa is desirable for improving the noise reduction of twin beams .
in conclusion , we have developed the multi - mode quantum theory for analyzing the noise characteristics of a non - degenerate phase insensitive fopa pumped by a pulsed laser .
the calculation shows that both the noise figure and the intensity difference noise of twin beams depend on the jsf . in the high gain regime ,
the noise figure is generally greater than the 3db quantum limit unless the jsf is factorable and the spectrum of the injected signal well matches the gain spectrum in signal band ; while the intensity difference noise can be significantly less than the snl . to closely resemble the experiments , the quantum noise of twin beams generated by a broadband fopa
is numerically studied by taking the real experimental conditions into account .
in addition to the influences of the quantum efficiency of detectors and the excesses noise of input signal , the influence of the collection efficiency of the twin beams , determined by the ratio of the pump bandwidth to signal bandwidth @xmath303 and the ratio of the signal bandwidth to filter bandwidth @xmath361 , are carefully analyzed as well .
our theoretical investigation is of practical importance , it not only serves as a guideline for optimizing the measured noise reduction of the twin beams , but also helps for understanding the quantum noise preformance of pulse pumped fopa .    instead of using the bloch - messiah reduction or schmidt decomposition of the jsf , which has been exploited to theoretically analyze the pulsed cv nonclassical light via the high gain spontaneous parametric down conversions in refs .
@xcite , our transformation functions of the bogoliubov transformation are expressed in the infinite series of the gain coefficient @xmath269 .
the accuracy of our numerical simulation can be controlled by the two parameters : one is the precision of the jsf determined by the step size of the frequency range , the other is the order of the series determined by the gain of fwm .
so our method has a greater flexibility .
we believe our theory can be extended to a more generalized opa , including the degenerate opa and the phase - sensitive opa .
moreover , our theory can also be used to study the quantum noise of other quantities , such as the noise correlation between the quadrature amplitudes of the signal and idler beams .
however , it is worth pointing out that the fopa in our model is simplified as a linear amplifier , and we only considered the fwm process . therefore , our calculation only quantitatively explains the results obtained under the condition of negligible pump depletion  @xcite . to model the fopa system
more accurately , other effects like raman effect , higher - order fwm and gain saturation effect should be included as well  @xcite .
this work was supported in part by the state key development program for basic research of china ( no .
2010cb923101 ) and by the specialized research fund for the doctoral program of higher education of china ( no .
20120032110055 ) .