the ability to coherently transfer atoms out of a bose - einstein condensate ( bec ) in a trapped atomic gas is of fundamental importance in the creation of coherent matter wave generators , often called `` atom lasers '' .
theoretical interest has developed toward characterizing the out - coupling stage of such an atom laser , with the hope of developing efficiency criteria for this process .
all this work has been done in the zero temperature limit , using either a non - interacting gas coupled to a laser mode @xcite or the modified gross - pitaevskii dynamics of the condensate order parameter @xcite .
a complete understanding of the out - coupling process necessarily must take into account the presence of the thermal cloud composed of noncondensate atoms .
recently , the role that the noncondensate atoms play in the out - coupling process has been discussed by japha , choi , burnett and band @xcite . in their work ,
a `` pair - breaking '' contribution was obtained , in which a correlated pair of atoms take part in the output process , with one output atom appearing simultaneously with the creation of an excitation in the trapped bose gas .
out - coupling experiments are only a small class of experiments that can be performed in which a bec is coupled to another gas of atoms ( that may or may not be bose - condensed ) .
the mit group @xcite has performed work in which a bec is scattered from an optical bragg grating .
they show that the excitation spectrum of the scattered atoms gives a direct experimental measure of the dynamic structure factor of the initial bose - condensate .
since the detailed form of the dynamic structure factor @xmath2 is dependent upon the many - body correlations of the bose - condensate , this kind of experiment provides a means of probing the unique correlations induced by a bose - condensate @xcite .    despite the large amount of work that has been performed with regard to the coupling of bose - condensed gases with another ( possibly bose - condensed ) gas ,
there still does not exist in the literature a clean formulation of such a coupling .
even though refs.@xmath3@xcite have sparked interest in this issue , the formalism used in these papers does not show in a transparent way which correlations present in a bose - condensed gas can be probed by such experiments .
moreover , the case in which both trapped gases are bose - condensed was not considered .    in this paper , we set up a general linear response theory approach for treating the coupling of two bose gases at finite temperatures , which gives a formally exact expression for the tunneling current .
the general expressions we derive  within linear response theory  are independent of the specific approximation used for the many - body states of the two uncoupled , trapped bose gases . in this sense , our results are the analogue of showing how the inelastic neutron scattering cross - section is related to the dynamic structure factor @xcite and the general theory of tunneling in superconductors @xcite .
the work closest to ours is by meier and zwerger @xcite , who applied a similar formalism in discussing coherent josephson tunneling between two uniform bose - condensed gases . in sections [ section : linear response theory - general ] and [ section : tunneling currents in terms of correlations functions ] ,
we give a detailed account of the many - body theory needed to find formal expressions for the tunneling currents in terms of single - atom correlation functions .
we apply our theoretical formalism to a class of experiments involving raman scattering of atoms between trapped bose gases .
we calculate the tunneling current expression at finite temperatures using the bogoliubov - popov approximation for the single - particle excitations in the trapped bose gases @xcite .
we discuss coherent tunneling of atoms between two bose - condensed gases , as well as the out - coupling of atoms from a bec in a trap to the vacuum ( i.e. ,  to free particle states ) . in this latter case , as expected , we show that our more general formalism reproduces the results of refs.@xmath3@xcite in the appropriate limit . in section [ sec : lda ] , we give an explicit calculation of the out - coupling current based on treating the noncondensed atoms within a simple local density approximation ( lda ) , such as used in refs .
@xcite .
in essence , we show that besides the dominant contribution from tunneling of atoms out of the condensate , one finds a much smaller contribution which arises from the atoms tunneling out of the noncondensate .
this latter contribution , however , is of special interest since it is directly related to the single - particle correlation function of the noncondensate atoms .
one component has weight proportional to the bogoliubov amplitude @xmath4 , where @xmath5 is the bose - distribution function for an excitation of energy @xmath6 , and increases with temperature . the other , more interesting , component has weight proportional to the bogoliubov amplitude @xmath7 .
this latter process is formally related to the `` negative energy pole '' of the bogoliubov spectral density .
this so - called `` pair - breaking '' term @xcite is of interest as a unique signature of the unusual correlations induced by a bose condensate .
it corresponds to the out - coupling of an atom and the simultaneous creation of a bogoliubov excitation in the trapped gas .
this contribution already is present at zero temperature and , in contrast to the tunneling current from the condensate , contributes to an output tunneling current of atoms with a wide spectrum of energies . as the temperature increases
, the weight of the @xmath8 term decreases while the @xmath9 term starts to build up ( we recall that as the condensate density goes to zero , @xmath10 and @xmath11 ) .
the analogue of this pair - breaking process was discussed earlier in the context of a bose - condensed gas of excitons @xcite .
in contrast to our work , a proposal has recently been put forward @xcite to measure these bogoliubov amplitudes using an optical bragg scattering technique . in this approach
, a two - photon bragg scattering beam is used to excite quasiparticles out of the condensate and then a second bragg beam is used to probe the altered momentum distribution of the sample .
recent work at mit has beautifully verified these predictions @xcite .
unlike this situation , in the raman scattering experiment the energy absorbed by the atoms imparts a momentum kick in addition to causing a change of hyperfine state , allowing the scattered atom to no longer be trapped .
our results give another fairly direct method of measuring the bogoliubov quasiparticle amplitudes . furthermore , since the weight of the out - coupling current is proportional to the bogoliubov amplitude @xmath1 , it is significantly amplified for larger values of the s - wave scattering length @xmath12 such as can be achieved by working close to a feshbach resonance .
we hope our work further stimulates experimental interest in this kind of out - coupling experiment , even though the current from the noncondensate atoms is quite small compared to that from the condensate .
our goal is to calculate the current of atoms tunneling between two atomic gases coupled through a weak external field , which can be treated as a time - dependent perturbation to a time - independent many - body hamiltonian using linear response theory ( see , for example , refs.@xmath3@xcite ) .
the total hamiltonian in the schrdinger representation is@xmath13 where @xmath14 is the hamiltonian for the many - body systems , @xmath15 and @xmath16 are the atomic hyperfine level energies of an atom in gas 1 and gas 2 respectively , @xmath17 is the atomic number operator for gas 1(2 ) , given by@xmath18 the external probe field @xmath19 is explicitly time - dependent .
we assume that the two many - body systems do not interact with one another in the absence of @xmath20 , so that @xmath21 .
this amounts to assuming that the two gases do not overlap , _
_ _ the coupling term of the kind @xmath22 can be neglected . in the absence of @xmath19 ,
the number of atoms in the two systems is fixed by their chemical potentials @xmath23 and @xmath24 .
the tunneling perturbation @xmath19 couples the two many - body systems together by introducing a mechanism by which an atom can tunnel between the two systems .
this perturbation has terms which create an atom in one system while destroying an atom in the other system , and vice - versa .
an example of the kind of experiment described by the tunneling perturbation @xmath19 is the out - coupling phenomenon considered in refs .
@xcite . in this experiment ,
an external electromagnetic ( em ) field causes transitions between two atomic hyperfine levels , one of which is not trapped .
we label the quantum states corresponding to these hyperfine levels by @xmath25 and @xmath26 . assuming that hyperfine level @xmath25 is lower in energy than hyperfine level @xmath26
, energy will be absorbed from ( emitted to ) the external field facilitating the transition @xmath27 ( @xmath28 ) . in more general experiments
, however , the transition can take place between two atomic hyperfine levels that are both trapped .
an example of such a process is a stimulated raman transition between the magnetic hyperfine levels @xmath29 and @xmath30 ( for a detailed discussion , see ref .
we will return to the specifics of a bec coupling experiment in section [ section : coupling between two trapped gases ] when we calculate the tunneling current between two trapped bose - condensed gases .    for the moment it suffices to write down the semiclassical form of the tunneling perturbation @xmath20 arising from the electric - dipole and rotating wave approximations @xcite , @xmath31 here
@xmath32 is the effective energy transfer from the laser fields , @xmath33 is the momentum transfer from the laser fields , and @xmath34 ( taken to be real ) is the interaction strength that is associated with the rabi frequency of the transition between atomic states .
we have also introduced a factor @xmath35 that models the adiabatic switching on of the interaction at @xmath36 . in our final results ,
we take the limit @xmath37 .
we use units in which @xmath38 and @xmath39 , unless otherwise specified .
we next outline the formalism for calculating the tunneling current between two bose gases ( resulting from the action of @xmath20 ) .
the state of the combined system evolves in time according to @xmath40 , where the evolution operator @xmath41 is unitary and satisfies the equation of motion@xmath42 this equation of motion can be formally solved to give @xmath43 , where @xmath44 denotes the time - ordering operator . in the heisenberg picture , operators are written in terms of the corresponding schrdinger picture operators as @xmath45 , while the state of the system is given by the initial state of the system at @xmath36 , namely @xmath46 .
we use the subscript @xmath47 to denote an operator or state in the heisenberg picture with respect to the full hamiltonian @xmath48 .
the heisenberg equation of motion for any operator @xmath49 is then given as@xmath50 .\ ] ]    we define the atomic current _ from _ gas 1 _ to _ gas 2 in terms of the thermal equilibrium ensemble average of @xmath51@xmath52 we have introduced a negative sign in this expression for the atomic current so that the current flow from gas 1 to gas 2 is positive . the ensemble average in ( [ current definition ] ) is given by @xmath53 \hat{u}(t)\right\ } , \label{ensemble average for n1dot 1 } \end{aligned}\ ] ] where the initial equilibrium density operator @xmath54 describes the state of the system at @xmath36 before the coupling @xmath20 is turned on . since @xmath55 = \left [ \hat{v}(t),\hat{n}_{1}\right ]   $ ] , the average in ( [ ensemble average for n1dot 1 ] ) takes the form@xmath56 \hat{u}(t)\right\ } .\ ] ] we now turn to the task of evaluating this ensemble average within linear response theory @xcite .
we first write the average ( [ ensemble average for n1dot 2 ] ) in an interaction picture , in which @xmath20 is taken as a perturbation to the many - body hamiltonian @xmath57 of the uncoupled gases .
later we will introduce an interaction picture involving @xmath58 .
with respect to @xmath59 , we define the interaction representation as @xmath60 in this interaction picture , we also define the time - dependent density operator as @xmath61 using these definitions , we now can write ( [ ensemble average for n1dot 2 ] ) as @xmath62 \mathrm{e}^{-i\hat{h}_{0}^{\prime } t}\right\ } .\ ] ] evaluating the commutator @xmath63   $ ] in ( [ ensemble average for n1dot 3 ] ) , we obtain@xmath64.\ ] ]    the density operator in the interaction picture defined in ( [ density matrix in interaction picture ] ) satisfies the liouville equation @xmath65 , \ ] ] which has a formal solution@xmath66 .\ ] ] iterating ( [ liouville equation in integral form ] ) to first order in @xmath20 , we obtain @xmath67 here@xmath68\ ] ] is the perturbation to the initial equilibrium density operator @xmath54 to first order in @xmath20 , the linear response approximation to the density matrix .    within this first - order approximation for @xmath69 the ensemble average of an operator @xmath49 is given by @xmath70 where @xmath71 the atomic current given by ( [ current definition ] ) can thus be written in linear response theory as@xmath72 where @xmath73 is the equilibrium current and
@xmath74 is the linear response current . substituting our approximation for @xmath69 of ( [ linear response theory approximation to density matrix in interaction picture ] ) together with ( [ v ] ) for @xmath20 into ( [ ensemble average for n1dot 4 ] ) ,
we obtain , after some algebra , explicit expressions for the equilibrium current @xmath75 and the linear response current @xmath76}\left\langle \left [ \hat{\psi } ^{\dagger } _ { 1h_{0}^{\prime } } ( \mathbf{r}^{\prime } , t^{\prime } ) \hat{\psi } _ { 2h_{0}^{\prime } } ( \mathbf{r}^{\prime } , t^{\prime } ) , \hat{\psi } _ { 1h_{0}^{\prime } } ( \mathbf{r},t)\hat{\psi } ^{\dagger } _ { 2h_{0}^{\prime } } ( \mathbf{r},t)\right ] \right\rangle _ { eq}\nonumber \\   &   & + \gamma ^{2}\mathrm{e}^{i[\mathbf{q}\cdot ( \mathbf{r}+\mathbf{r}^{\prime } ) -\omega ( t+t^{'})]}\left\langle \left [ \hat{\psi } _ { 1h_{0}^{\prime } } ( \mathbf{r}^{\prime } , t^{\prime } ) \hat{\psi } ^{\dagger } _ { 2h_{0}^{\prime } } ( \mathbf{r}^{\prime } , t^{\prime } ) , \hat{\psi } _ { 1h_{0}^{\prime } } ( \mathbf{r},t)\hat{\psi } ^{\dagger } _ { 2h_{0}^{\prime } } ( \mathbf{r},t)\right ] \right\rangle _ { eq}\biggr \}.\label{linear response current 1 } \end{aligned}\ ] ] the expressions in ( [ equilibrium current 1 ] ) and ( [ linear response current 1 ] ) are the basis of the rest of this paper .
one sees that there are two different contributions to the linear response current in ( [ linear response current 1 ] ) , corresponding to phases depending on the times through @xmath77 or @xmath78 .
the former will be seen shortly to correspond to a _
normal _ current of single atoms between the two gases , while the latter corresponds to an _ anomalous _ current of atoms between the two gases .
we next want to express the ensemble averages in ( [ equilibrium current 1 ] ) and ( [ linear response current 1 ] ) in terms of thermal green s functions given by standard techniques in finite temperature many - body theory @xcite . to this end
, we first express the time - dependent field operators in ( [ equilibrium current 1 ] ) and ( [ linear response current 1 ] ) in the usual `` heisenberg picture '' used in many - body theory , namely the interaction picture in which the free hamiltonian is taken to be the many - body `` grand canonical hamiltonian '' @xmath79 using the identity @xmath80 , it is easily verified that @xmath81 where we define @xmath82 . using ( [ relation between destruction operators in h0prime and k0 interaction pictures ] ) , it is easy to rewrite ( [ linear response current 1 ] ) in terms of the @xmath83 operators .
we also note that the grand canonical ensemble average of a product of operators of the two gases factorizes @xmath84 where @xmath85 .    below the transition temperature of a bose gas ,
the new phase of matter is described by a symmetry - broken state in which the equilibrium density operator does _ not _ commute with the number operator : @xmath86 \neq 0 $ ] ( the state breaks the system s @xmath87 symmetry ) .
thus , the equilibrium ensemble averages are truly over the @xmath88-ensemble introduced by bogoliubov @xcite , which is symbolically defined as @xmath89 where @xmath90 is the symmetry - breaking grand canonical equilibrium density operator .
the symmetry - breaking hamiltonian is defined as @xcite @xmath91 @xmath92 is the broken symmetry grand canonical ensemble partition function .
we will use the notation @xmath93 to denote equilibrium ensemble averages , with the understanding that they are over a broken symmetry state . an equivalent way of stating this is that atomic field operators can be decomposed into two parts : @xmath94 where @xmath95 is the symmetry - breaking order parameter field and @xmath96 is the fluctuation field operator associated with the noncondensate atoms . in
( [ order - parameter ] ) , @xmath97 is the local condensate density while @xmath98 is the condensate phase . in our formalism
based on linear response from thermal equilibrium , the bose order parameter in ( [ order+fluctuation - field ] ) is that appropriate to thermal equilibrium .
a phase @xmath99 would arise if we were dealing , for example , with a condensate with a vortex .
we also note that the usual time - dependence of the bose - order parameter ( @xmath100 ) is absent because we are working in the @xmath101 interaction representation .
we do not need to assume that the inverse temperature is the same for both gases . indeed
, because the two gases are assumed to be uncoupled in the absence of the perturbation @xmath20 , no mechanism exists to force both gases to reach a common thermodynamic equilibrium .
therefore , the equilibrium ensemble averages in ( [ equilibrium current 1 ] ) , ( [ linear response current 1 ] ) and ( [ ensemble average decomposition ] ) may be viewed as being over two independent gases .
the equilibrium tunneling current in ( [ equilibrium current 1 ] ) is first order in the perturbation @xmath20 , whereas the linear response current in ( [ linear response current 1 ] ) is second order in @xmath20 . using ( [ relation between destruction operators in h0prime and k0 interaction pictures ] ) , ( [ equilibrium current 1 ] ) can be written as @xmath102 where @xmath103 , @xmath104 . introducing ( [ order - parameter ] ) into ( [ equilibrium current 2 ] ) , the equilibrium current is reduced to@xmath105.\ ] ] the magnitude of this equilibrium current is seen to be proportional to the overlap of the two condensate density profiles and thus serves as a probe of the condensate surface structure .    in a non - bose condensed gas ,
the ensemble averages @xmath106 and @xmath107 in ( [ equilibrium current 2 ] ) vanish . as noted in ( [ order - parameter ] ) and the discussion preceding it
, these ensemble averages do _ not _ vanish in a bose - condensed system .
thus , the equilibrium current ( [ equilibrium current 3 ] ) is a direct manifestation of bose broken - symmetry .
this expression for the equilibrium current is a more general form of the _ linearized _ `` internal '' josephson current @xcite for bose - condensed gases .
our expression ( [ equilibrium current 3 ] ) differs from previous results in several ways .
firstly , there is the phase associated with the momentum @xmath33 transferred by the two laser beams .
this will introduce momentum kicks ( bragg scattering ) to the atoms in the two condensates @xcite .
secondly , we have allowed for the case of a spatially varying condensate phase difference , @xmath108 , such as would arise if one of the condensates had a vortex .
the quantity @xmath109 in ( [ equilibrium current 3 ] ) is analogous to the voltage in the standard discussion of the josephson current in superconductors @xcite
. we will show later ( see ( [ deltaiaa ] ) and ( [ deltaiab ] ) ) that this quantity also enters into the phase of the anomalous component of the linear response current in ( [ linear response current 1 ] ) as a second harmonic .
we may interpret the atomic energy level splitting of the atomic states @xmath110 in addition to the chemical potential difference @xmath111 of the two bose gases as a negative voltage impeding the atomic current flow _ from _ gas 1 _ to _ gas 2 .
the energy difference between the laser photons acts as a forward biased voltage serving to negate this reverse bias . in the regime
where @xmath112 , we recover what is termed the d.c .  josephson effect in superconducting systems @xcite . in this regime , as with superconductors , a non - zero difference between the phases of the two condensate order parameters drives a d.c .
atomic current .
we next turn to the linear response current in ( [ linear response current 1 ] ) , which is second order in the tunneling perturbation @xmath20 .
we first want to write it in terms of thermal green s functions by introducing correlation functions of the atomic field operators . as in the previous section ,
we rewrite our time - dependent interaction picture field operators in terms of a heisenberg picture defined with respect to @xmath101 , rather than @xmath59 ( see ( [ relation between destruction operators in h0prime and k0 interaction pictures ] ) ) .
we consider the first term in the expression of the linear response current in ( [ linear response current 1 ] ) .
we refer to this term as the normal current , for reasons which shall become clear .
the normal current is then given by ( @xmath113 ) @xmath114.\label{normal current definition } \end{aligned}\ ] ] we have evaluated the commutators in ( [ linear response current 1 ] ) making use of ( [ ensemble average decomposition ] ) to write the ensemble averages as a product of two - point correlation functions such as @xmath115 in subsequent discussion , we will refer to @xmath116 as @xmath117 for notational simplicity .
the normal current given in ( [ normal current definition ] ) is independent of the time @xmath118 .
we make this more explicit by rewriting the integral and noting that factors of @xmath35 factor can be neglected in the limit @xmath37 . introducing the ( heaviside ) step function @xmath119 ,
the normal current may be written in the explicitly time - independent form @xmath120.\label{normal current 1 } \end{aligned}\ ] ] we observe that , as expected for a current , the second term in ( [ normal current 1 ] ) breaks time - reversal symmetry ; the current in the negative time direction is the negative of that in the positive time direction .
the correlation functions appearing in ( [ normal current 1 ] ) may now be written in terms of their fourier transforms @xmath121 and hence ( [ normal current 1 ] ) becomes @xmath122.\label{normal current 2 } \end{aligned}\ ] ] the integral over the time @xmath123 in ( [ normal current 2 ] ) is trivially given by @xmath124 in the physical limit @xmath37 , we may use the identity @xmath125 where @xmath126 denotes the principal value of an integral .
therefore , the normal current ( [ normal current 2 ] ) may be written in the form @xmath127 where we have introduced the two contributions @xmath128\label{deltaina } \end{aligned}\ ] ] and @xmath129}{\omega ^{\prime \prime } -\omega ^{\prime } -\omega + \delta \varepsilon + \delta \mu } .\label{short time normal current } \end{aligned}\ ] ]    we next show how the correlation functions in ( [ deltaina ] ) and ( [ short time normal current ] ) may be converted to spectral densities , which in turn can be directly calculated from time - ordered single - particle green s functions @xcite .
thus , the problem is reduced to the application of well - known techniques from many - body theory to calculate these green s functions @xcite .
we first rewrite the correlation functions in terms of the order parameter fields and the single - particle correlation functions of the noncondensate field operators .
the diagonal correlation functions above become @xmath130 and@xmath131 using the detailed balance identity @xcite@xmath132 in addition to the bose identity @xmath133 one may easily write @xmath134 in terms of a spectral density ( see p.@xmath37 of ref.@xmath3@xcite ) @xmath135 here @xmath136 is the bose distribution function @xmath137 and the single - particle spectral density is @xmath138 \right\rangle _ { eq}.\ ] ] note that for an inhomogeneous system , this spectral density is in general not real .    in turn , @xmath139 is given by the difference of the retarded and advanced green s functions , @xmath140 and @xmath141 respectively , for the noncondensate quantum field operators of a bose gas ( see p.@xmath3150 of ref.@xmath3@xcite ) : @xmath142 where @xmath143 is the fourier transform of @xmath144 \right\rangle _ { eq}\ ] ] and @xmath145 is the fourier transform of @xmath146 \right\rangle _ { eq}.\ ] ] the point of all this formalism is that standard finite - temperature many - body theory @xcite most naturally determines thermal green s functions at imaginary bose matsubara frequencies @xmath147 .
fortunately , one can determine @xmath148 by a unique analytic continuation @xmath149 to the real frequency axis .
this procedure thus gives the spectral density @xmath150 .
we may express ( [ deltaina ] ) and ( [ short time normal current ] ) in terms of the condensate order parameters and single - particle noncondensate spectral densities .
we find@xmath151 and@xmath152 in these expressions we have introduced the normal current densities @xmath153 which corresponds to atomic tunneling of single particles between the condensate and noncondensate components of the two gases .
in contrast , @xmath154 is associated with atomic tunneling of single atoms between the noncondensates of the two gases .
we observe that in both @xmath155 and @xmath156 , atoms make only two types of transitions between the states of the two gases .
the first transition process involves tunneling of an atom between the condensate of one gas and the noncondensate of the other .
it is represented in @xmath155 and @xmath156 by terms of the form @xmath157 .
this contributes at all temperatures including @xmath158 .
the second transition process involves the tunneling of an atom between the _ noncondensates _ of the two gases .
it appears in @xmath155 and @xmath156 as a term of the form @xmath159 .
this kind of transition can occur at zero temperature only when one of the gases is bose - condensed .
finally we consider the `` anomalous current '' given by the second term in ( [ linear response current 1 ] ) in terms of the condensate order parameters and noncondensate spectral density functions .
we follow the identical procedure used to find ( [ deltaina 1 ] ) and ( [ deltainb 1 ] ) for the contributions @xmath155 and @xmath156 to the normal current
. we obtain finally@xmath160.\label{anomalous current definition } \end{aligned}\ ] ] clearly this current may be written entirely in terms of the so - called anomalous correlation functions involving creation or destruction of _ pairs _ of atoms , hence the name anomalous current .
using @xmath161 , we first rewrite ( [ anomalous current definition ] ) as @xmath162\bigg \}.\label{anomalous current 1 } \end{aligned}\ ] ] this shows that the anomalous current is explicitly _ time - dependent _ , in contrast to the normal quasiparticle current given by ( [ normal current 2 ] ) .
it clearly involves a second harmonic of the first - order josephson current in ( [ equilibrium current 3 ] ) .    following the same procedure we used in writing the normal current in terms of the condensate order parameters and noncondensate spectral density functions of the of the two gases ,
the anomalous current may be rewritten as @xmath163 where we have introduced the quantities @xmath164\label{deltaiaa } \end{aligned}\ ] ] and @xmath165.\label{deltaiab } \end{aligned}\ ] ] the quantities @xmath166 and @xmath167 in the above expressions for @xmath168 and @xmath169 are given by @xmath170 and @xmath171 where we have introduced the total anomalous current density @xmath172 , which is written as @xmath173 here @xmath174 is the anomalous noncondensate single - particle spectral density function , analogous to the noncondensate single - particle spectral density function @xmath139 defined in ( [ spectral - density - definition ] ) .
the first and second terms in ( [ total anomalous tunneling current ] ) correspond to a tunneling current between the condensate of one gas and the noncondensate of the other gas . the third term in ( [ total anomalous tunneling current ] )
represents a tunneling current between the noncondensates of the two gases .
one can rewrite @xmath166 and @xmath167 such that the contributions from the condensate - noncondensate and noncondensate - noncondensate tunneling processes are made more explicit .
thus , we have ( see also ref .
@xcite ) @xmath175 and @xmath176 where we have introduced the condensate - noncondensate current density ( the analogue of ( [ condensate - noncondensate current density ] ) ) @xmath177 and the noncondensate - noncondensate current density ( the analogue of ( [ noncondensate - noncondensate current density ] ) ) @xmath178    we note that @xmath169 in ( [ deltaiab ] ) is the current most directly associated with the _ coherent _ tunneling responsible for the josephson effect in superconductors .
this current is in - phase with the josephson current in ( [ equilibrium current 3 ] ) but involves a second harmonic with a frequency @xmath179 .
it clarifies the role of the noncondensate in the bec josephson effect . on the other hand
, the @xmath168 current in ( [ deltaiaa ] ) is directly related to a _ dissipative _ interference current , just as one finds for the josephson effect in superconductors @xcite .
in this section , we use the results of section [ section : tunneling currents in terms of correlations functions ] to calculate the second order tunneling current involved in a raman scattering experiment .
the spectral density functions are calculated in the finite temperature bogoliubov - popov approximation @xcite .
the josephson tunneling current has already been discussed and is given by ( [ equilibrium current 3 ] ) .
a two - photon bragg scattering experiment involves two laser beams at different frequencies and with different propagation vectors @xcite .
however , we are interested in a different type of two - photon process in which an atom changes its hyperfine state .
this is a raman process and is depicted in fig.@xmath3[fig : raman - figure ] .
we now discuss how the effective frequency @xmath32 and wavevector @xmath33 in ( [ v ] ) are defined in such a raman scattering experiment @xcite . assuming that each laser beam has equal intensity @xmath180 , the combined intensity of the two beams at space - time point @xmath181 is given by @xmath182 where @xmath32 is the frequency difference between the two laser beams@xmath183 the wavevector difference between the two laser beams is @xmath184 where @xmath185 is the wavevector of laser beam @xmath186 .
we may take @xmath187 such that the standard bragg scattering relation ( see , for example , p.@xmath3103 of ref.@xmath3@xcite ) holds and @xmath188 is given by @xmath189 where @xmath119 is the angle between the two counter - propagating beams .
finally the coupling parameter @xmath34 in ( [ v ] ) is related to the single - photon rabi frequency @xmath190 ( assumed to be the same for both hyperfine level transitions ) and the laser beam detuning @xmath191 from the intermediate level @xmath192 ( see fig .  [
fig : raman - figure ] ) according to @xcite @xmath193 where @xmath194 is the so - called two - photon  rabi frequency  " .
the many - body hamiltonian ( in the absence of a tunneling perturbation ) is @xmath21 , where @xmath195 \hat{\psi } _ { 1}(\mathbf{r})+\frac{g_{1}}{2}\hat{\psi } _ { 1}^{\dagger } ( \mathbf{r})\hat{\psi } _ { 1}^{\dagger } ( \mathbf{r})\hat{\psi } _
{ 1}(\mathbf{r})\hat{\psi } _ { 1}(\mathbf{r})\right\ } .\ ] ] in this expression for @xmath196 , @xmath197 is the trapping potential and @xmath198 is the s - wave scattering pseudo - potential with @xmath199 being the s - wave scattering length for atoms in gas @xmath200 .
the expression for @xmath201 has an analogous form .
we recall that the spectral densities in section [ section : tunneling currents in terms of correlations functions ] are defined in terms of the grand canonical hamiltonian @xmath202 .
thus the excitation energies of gas @xmath203 are measured relative to the chemical potential @xmath204 .      in the standard static bogoliubov - popov approximation for excitations in a trapped bose gas @xcite
, the noncondensate field operator @xmath96 can be written as@xmath205 where @xmath206 and @xmath207 are the bogoliubov quasiparticle amplitudes , @xmath208 is the energy of bogoliubov state @xmath186 and @xmath209 ( @xmath210 ) annihilate ( create ) a quasiparticle in state @xmath186 , respectively .
these annihilation and creation operators are defined to satisfy the boson commutation relations@xmath211 = \left [ \hat{\alpha } _
{ m}^{\dagger } , \hat{\alpha } _ { n}^{\dagger } \right ] = 0;\qquad \left [ \hat{\alpha } _ { m},\hat{\alpha } _ { n}^{\dagger } \right ] = \delta _ { mn}.\ ] ] the
bogoliubov quasiparticle amplitudes must satisfy the finite temperature coupled bogoliubov - popov equations of motion @xcite@xmath212 in order for @xmath101 to be reduced to the diagonal form @xmath213 in ( [ bogoliubov - equations - of - motion ] ) , we have introduced the differential operator@xmath214 where @xmath215 is the total static density .    using the commutation relations ( [ bogoliubov - commutation - relations ] ) , one finds that the single - particle spectral density function can be expressed in terms of bogoliubov quasiparticles as@xmath216,\ ] ] which has the fourier transform@xmath217.\ ] ] we note the appearance in ( [ bogoliubov - spectral - density - frequency ] ) of negative energy poles ( @xmath218 ) in addition to the usual positive energy poles ( @xmath219 ) .
these are a result of unique processes induced by the bose condensate .
the two poles follow immediately from ( [ bogoliubov - quasi - particle - quantization ] ) , which shows how the destruction of an atom involves a linear superposition of the destruction ( amplitude @xmath220 ) and creation ( amplitude @xmath221 ) of bogoliubov excitations .    similarly , the off - diagonal or anomalous single - particle spectral density function of relevance to the anomalous current densities @xmath222 in ( [ condensate - noncondensate anomalous current density ] ) and @xmath223 in ( [ noncondensate - noncondensate anomalous current density ] ) is given by @xmath224.\ ] ] the spectral density functions can be shown to satisfy @xmath225 and @xmath226 .      using the results in ( [ bogoliubov - spectral - density - frequency ] ) and ( [ bogoliubov spectral density for anomalous correlations ] )
, we can now write the normal and anomalous linear response currents in terms of bogoliubov excitations . in section [ section : tunneling currents in terms of correlations functions ] , we found that two distinct physical processes contribute to the normal current .
the current density of atoms tunneling between the condensates and noncondensates of the two gases is given by ( [ condensate - noncondensate current density ] ) , or @xmath227.\label{condensate - noncondensate normal current density in bogoliubov approximation } \end{aligned}\ ] ] this type of process involves both condensate atoms and quasiparticles .
the current density associated with the tunneling of atoms between the noncondensates of the two gases has a more complicated form given by ( [ noncondensate - noncondensate current density ] ) , or @xmath228.\label{noncondensate - noncondensate normal current density in bogoliubov approximation } \end{aligned}\ ] ]    as in the case of the normal current , two different physical processes contribute to the anomalous current in ( [ anomalous current definition ] ) .
the first process involves pairs of condensate atoms tunneling to noncondensate states of the other bose gas .
this anomalous current density is given by ( [ condensate - noncondensate anomalous current density ] ) , namely @xmath229.\label{condensate - noncondensate anomalous current density in bogoliubov approximation } \end{aligned}\ ] ] the expression for the anomalous current density between condensate 1 and noncondensate 2 has a similar form .
the final process involves quasiparticles of both gases and is given by ( [ noncondensate - noncondensate anomalous current density ] ) , namely@xmath230.\label{noncondensate - noncondensate anomalous current density in bogoliubov approximation } \end{aligned}\ ] ]    as with the josephson current given in ( [ equilibrium current 3 ] ) , one sees from ( [ condensate - noncondensate normal current density in bogoliubov approximation])-([noncondensate - noncondensate anomalous current density in bogoliubov approximation ] ) that the magnitude of the normal and anomalous quasiparticle currents is determined by the spatial overlap of the bogoliubov wavefunctions described by @xmath231 and @xmath232 for the two systems .    when one uses ( [ condensate - noncondensate normal current density in bogoliubov approximation])-([noncondensate - noncondensate anomalous current density in bogoliubov approximation ] ) in the expressions for the normal ( see ( [ deltaina 1 ] ) and ( [ deltainb 1 ] ) ) and anomalous ( see ( [ deltaiaa ] ) and ( [ deltaiab ] ) ) currents , our results are essentially equivalent to expressions discussed in ref .
the main differences are that meier and zwerger approximated the two bose gases as being uniform ( or homogeneous ) and used the single - particle spectral densities at @xmath158 .
in this section , we consider the special case of atoms out - coupling to the vacuum as discussed in refs .  @xcite
. more precisely , we first deal with the case when gas 2 is a uniform non - interacting gas . the corresponding spectral density function then has the well - known form @xcite@xmath233 where@xmath234 is the the excitation energy of an atom in gas @xmath235 . because the gas is assumed to be in the classical limit , we have @xmath236 .    the only terms in ( [ condensate - noncondensate normal current density in bogoliubov approximation ] ) which are finite reduce to @xmath237 and in ( [ noncondensate - noncondensate normal current density in bogoliubov approximation])@xmath238v_{1j}(\mathbf{r}^{\prime } ) v^{*}_{1j}(\mathbf{r})\delta \left ( \omega + e_{1j}\right ) \biggr \}\delta ( \tilde{\varepsilon } _ { \mathbf{k}2}-\omega ^{\prime } ) .\label{outcoupling current density noncondensate - noncondensate } \end{aligned}\ ] ] we have used the bose identity @xmath239 $ ] in writing the last term in ( [ outcoupling current density noncondensate - noncondensate ] ) . the condition that @xmath240 results in negligible tunneling from gas 2 back to gas 1 .
the current density @xmath241 in ( [ outcoupling current density condensate - noncondensate ] ) corresponds to the _ destruction _ of a condensate atom in gas @xmath200 and the _ creation _ of an excitation in gas @xmath235 with energy @xmath242 .
the first term in the density @xmath243 of ( [ outcoupling current density noncondensate - noncondensate ] ) represents what might be called the `` quantum evaporation '' of an atom , _ destroying _ a bogoliubov quasiparticle with energy @xmath244 and _ creating _ an atom in gas @xmath235 with @xmath245 . the second term in ( [ outcoupling current density noncondensate - noncondensate ] ) corresponds to _ creating _ an atom in gas @xmath235 with energy @xmath245 and simultaneously _ creating _ a quasiparticle excitation in gas @xmath200 with energy @xmath246 $ ] .
this second out - coupling channel in ( [ outcoupling current density noncondensate - noncondensate ] ) is a direct consequence of the unusual condensate - induced correlations in a bose gas .
it corresponds to the `` pair - breaking '' process discussed by the oxford group @xcite .
the analogue of this process appears in the predicted lineshape of the recombination photons emitted in the decay of an exciton in a bose - condensed gas of optically - excited excitons @xcite .
all three out - coupling contributions will be calculated in section [ sec : lda ] .    using ( [ outcoupling current density condensate - noncondensate ] ) and ( [ outcoupling current density noncondensate - noncondensate ] ) , the out - coupling current @xmath247 is@xmath248v_{1j}(\mathbf{r}^{\prime } ) v^{*}_{1j}(\mathbf{r})\delta \left ( k^{2}/2m-(\omega -\delta \varepsilon + \mu _ { 1})+e_{1j}\right ) \biggr ] \biggr \}.\label{bogoliubov outcoupling current } \end{aligned}\ ] ] we refer to the second term as the `` @xmath0 contribution '' and the third term as the `` @xmath1 contribution '' .
it is useful to describe the main features of the out - coupling current spectrum arising from the three terms in ( [ bogoliubov outcoupling current ] ) , as a function of the energy @xmath245 of the out - coupled atoms , for fixed @xmath32 and * @xmath33*. the energy conserving delta functions in ( [ bogoliubov outcoupling current ] ) determine the energies of the out - coupling current peaks for the different processes .
there is a delta function peak at @xmath249 from atoms coming directly from the condensate , where @xmath23 is the condensate energy of a trapped atom in the hyperfine state @xmath25 and @xmath250 is the hyperfine state energy .
relative to this condensate pole , the peak associated with the quantum evaporation term is shifted by a positive amount @xmath251 , while the peak of the `` pair - breaking '' channel is shifted by a negative amount @xmath252 .
these shifts follow from the fact that the quantum evaporation channel arises from the positive energy pole ( related to @xmath0 ) and the pair - breaking channel corresponds to the negative energy pole ( related to @xmath1 ) in the single - particle spectral density .
we next consider the behavior of the three out - coupling current contributions in ( [ bogoliubov outcoupling current ] ) as a function of temperature . for @xmath253 , since @xmath254 , ( i.e.  no thermally excited quasiparticles ) the second channel ( @xmath0 ) has zero weight .
thus , at @xmath253 , only the first and third processes occur . when gas 1 is at a finite temperature ( @xmath255 ) , the second channel ( @xmath0 ) begins to contribute to the out - coupling current , taking weight away from the pair - breaking channel ( @xmath1 ) .
all these features are shown explicitly in the numerical results given in section [ sec : lda ] ( see figs .  [
fig : condensate - noncondensate out - coupling current for q=3d3d2/aho and a=3d3d5nm]-[fig : total out - coupling current for q=3d3d10/aho a=3d3d5 nm t=3d3d0.2tbec ] ) .
the out - coupling current derived in ref .
@xcite is given as the rate of atoms excited out of gas 1 with energy @xmath242 and may be written as ( where we use a notation more consistent with our own)@xmath256v_{1j}(\mathbf{r}^{\prime } ) v^{*}_{1j}(\mathbf{r})\delta \left ( k^{2}/2m-(\omega -\delta \varepsilon + \mu _ { 1})+e_{1j}\right ) \biggr ] \biggr \},\label{oxford - result } \end{aligned}\ ] ] where @xmath257 is the wavefunction for the out - coupled excitation with energy @xmath242 .
we can make a direct comparison between the integrands of ( [ bogoliubov outcoupling current ] ) and ( [ oxford - result ] ) . indeed taking @xmath258 ( a free atom with momentum @xmath259 ) reduces ( [ oxford - result ] ) to the integrand of ( [ bogoliubov outcoupling current ] ) .
thus our result for the out - coupling current reduces to that derived in refs.@xmath3@xcite , where the out - coupled atoms were in free - particle states .
in this section , we evaluate the out - coupling current in ( [ bogoliubov outcoupling current ] ) and estimate the magnitude and temperature - dependence of the quantum evaporation and pair - breaking effects discussed in section [ sec : case of out - coupling to a non - bose - condensed gas ] . in previous work @xcite
, these noncondensate processes were found to be a few percent , relative to the dominant contribution of atoms tunneling out of the high density condensate .
it has also been found @xcite that in current bose gas experiments , the contribution of the noncondensate atoms to the dynamic structure factor @xmath2 at @xmath260 is about @xmath261 of the contribution from the condensate fluctuations ( the latter has been studied in refs .
@xcite at @xmath158 ) .
this small amplitude is , of course , a direct reflection of the strongly peaked condensate at the centre of the trap compared to the broad , low density profile of the noncondensate atoms .
however , while small , the tunneling current from the noncondensate atoms gives direct information about the bose - induced correlations and thus is a worthwhile goal for future bragg scattering and raman experiments .
the condensate current in ( [ bogoliubov outcoupling current ] ) can be rewritten in the form@xmath262 where @xmath263 is the fourier transform of the bose order parameter @xmath264 and we have introduced the abbreviation @xmath265 . it is also useful to introduce a `` detuned frequency '' @xmath266 the noncondensate contribution in ( [ bogoliubov outcoupling current ] ) is more conveniently written in terms of the trapped gas bogoliubov spectral density in ( [ bogoliubov - spectral - density - frequency ] ) , @xmath267 in the rest of this section , for notational simplicity , we shall drop the index `` 1 '' on the quantities associated with the trapped gas .
equations ( [ condensate outcoupling current ] ) and ( [ noncondensate outcoupling current ] ) summarize our main results for the tunneling currents into the vacuum , showing very clearly which properties of the trapped bose gas they depend on .      in the standard thomas - fermi approximation ,
the condensate density profile is given by @xcite @xmath268 where the chemical potential is@xmath269 and @xmath270 the thomas - fermi radius @xmath271 defines the size of the condensate , with @xmath272 the total number of particles in the bose - condensate at temperature @xmath273 and @xmath274 is the oscillator length .
the @xmath158 thomas - fermi results have been extended to finite temperature by adjusting @xmath275 according to the noninteracting trapped gas result @xcite @xmath276 where @xmath277 is the noninteracting trapped gas bose - condensation temperature . using ( [ thomas - fermi condensate density ] ) for the order parameter @xmath278 , the fourier transform in ( [ order parameter fourier
transform definition ] ) can be calculated analytically @xmath279 where @xmath280 is the ordinary bessel function of order 2 .
this result is the same as previously obtained @xcite , taking into account that we use harmonic oscillator units ( length is measured in units of @xmath274 and energy is measured in terms of @xmath281 ) .
writing the integral in ( [ condensate outcoupling current ] ) in terms of spherical coordinates , we obtain using ( [ order parameter fourier transform ] ) @xmath282 where @xmath283 .
the result ( [ lda condensate outcoupling contribution in harmonic oscillator units ] ) may be easily generalized to asymmetric harmonic traps , having trapping potential @xmath284 by introducing the coordinate transformation@xmath285 where @xmath286 .
the same procedure as used above follows through with the same final expressions , apart from @xmath287 being replaced by @xmath288 .
in addition , any wavevector @xmath259 in our tunneling expressions such as ( [ condensate outcoupling current ] ) needs to be interpreted as @xmath289 in an anisotropic trap @xcite .
the total out - coupling current from the condensate given by ( [ condensate outcoupling current ] ) , as a function of the detuning parameter @xmath290 in ( [ definition of detuning parameter ] ) , is shown by the solid lines in figs .
[ fig : condensate - noncondensate out - coupling current for q=3d3d2/aho and a=3d3d5 nm ] and [ fig : total out - coupling current for q=3d3d2/aho a=3d3d5 nm t=3d3d0.2tbec ] .
in all our numerical calculations , for the two - photon rabi coupling we have used @xmath291 . note the logarithmic scale for the out - coupled atom current .
the current magnitude only decreases as a function of increasing temperature by about @xmath292 . from fig .
[ fig : condensate - noncondensate out - coupling for q=3d3d10/aho and a=3d3d5 nm ] , it can be seen that the maximum condensate tunneling current amplitude is sensitive to the magnitude of the momentum transfer @xmath33 .
the decrease in the maximum tunneling current with increasing values of @xmath188 is a result of the fact that the momentum distribution of atoms in the condensate is peaked at zero momentum , decreasing for increasing momentum . for an asymmetric trap ,
this means that the directions of strongest confinement ( largest trap frequency ) will out - couple most efficiently .
we observe also that the peak of the condensate tunneling current shifts toward positive @xmath290 as @xmath188 increases .
the out - coupling current ( [ noncondensate outcoupling current ] ) from the noncondensate atoms will be estimated using a simple local density approximation ( lda ) .
the lda has also been recently used to calculate @xmath2 in trapped gases @xcite .
we first introduce a new coordinate system defined in terms of the relative coordinate@xmath293 and the center - of - mass coordinate@xmath294 writing the single - particle spectral density function @xmath295 in ( [ noncondensate outcoupling current ] ) in terms of these new coordinates , we note that for large condensates and quasiparticle energies much larger than @xmath281 , the spectral density will vary slowly as a function of @xmath296 .
it is thus plausible to approximate @xmath139 by its homogeneous form , but adjusted so that at each point @xmath296 the appropriate condensate density @xmath297 is used .
thus , in the lda , the single - particle spectral density function in ( [ noncondensate outcoupling current ] ) is approximated by@xmath298 where @xmath299 denotes the _ uniform _ gas expression for the spectral density function for a condensate density @xmath97 .    in the bogoliubov - popov approximation ,
the lda form ( [ lda approximation for spectral density ] ) for the spectral density in energy - momentum space is @xcite@xmath300.\label{lda bogoliubov spectral density in momentum space } \end{aligned}\ ] ] thus , in the lda the quasiparticle energies become @xmath301 while the bogoliubov quasiparticle mode amplitudes are given by@xmath302 we have introduced the free atom energy @xmath303 and@xmath304 plays the role of an  effective  " local chemical potential , modifying the chemical potential @xmath305 as a result of the inhomogeneity of the trap .    within the lda , the noncondensate contribution ( [ noncondensate outcoupling current ] ) to the out - coupling current
is given by@xmath306 where as before we have introduced the variable @xmath307 for notational simplicity .
using the identity@xmath308 where @xmath309 , we may evaluate the integral over @xmath296 in ( [ lda noncondensate outcoupling current version 1 ] ) assuming a spherically symmetric trap .
after some algebra , one obtains the following expression for the noncondensate contribution to the out - coupling current@xmath310 where we have introduced the lda spatially averaged spectral density function in energy - momentum space@xmath311 where@xmath312 and harmonic oscillator units are used .
note that the noncondensate out - coupling current in ( [ lda noncondensate outcoupling current in harmonic oscillator units ] ) may be generalized to an asymmetric trap in the same manner as discussed earlier for the condensate out - coupling current ( see remarks below ( [ coordinate transform for asymmetric trap to symmetric trap ] ) ) .
the formula in ( [ lda noncondensate outcoupling current in harmonic oscillator units ] ) has a simple physical interpretation .
the integrand gives the number of emitted atoms with momentum @xmath259 per second .
the first term in ( [ lda energy - momentum space spectral density function ] ) is associated with the positive energy pole in ( [ bogoliubov - spectral - density - frequency ] ) , with amplitude @xmath0 , and only occurs for positive frequency @xmath313 .
the second term in ( [ lda energy - momentum space spectral density function ] ) arises from the negative energy pole in ( [ bogoliubov - spectral - density - frequency ] ) , with amplitude @xmath1 , and only occurs for negative frequency @xmath314 .
these first two terms are associated with quasiparticles in the region of the ( thomas - fermi ) condensate , i.e.  for @xmath315 . the third term , in contrast ,
is associated with quasiparticles lying outside of the condensate region .
these quasiparticles comprise a hartree - fock gas in our simple model calculation .
the dominant weight in ( [ lda energy - momentum space spectral density function ] ) comes from small @xmath316 and @xmath32 .
the quasiparticles with @xmath0 and @xmath1 weights in the region where there is a condensate are the dominant contribution . for @xmath317 ,
one also finds a contribution from the normal hartree - fock gas of atoms existing outside of the condensate .
as the temperature increases toward @xmath318 , this contribution from the normal hartree - fock gas becomes more significant relative to the bogoliubov quasiparticles .
of course , for @xmath319 , it is the sole contribution .
we have evaluated the integral appearing in ( [ lda noncondensate outcoupling current in harmonic oscillator units ] ) numerically .
some results , as a function of the detuning parameter @xmath290 for various temperatures , are shown in figs .
[ fig : condensate - noncondensate out - coupling current for q=3d3d2/aho and a=3d3d5 nm ] and [ fig : total out - coupling current for q=3d3d2/aho a=3d3d5 nm t=3d3d0.2tbec ] .
note the logarithmic scale for the out - coupled atom current .
as a function of detuning , the out - coupling current separates into two regions . from ( [ lda noncondensate outcoupling current in harmonic oscillator units ] ) and
( [ lda energy - momentum space spectral density function ] ) , we observe that for @xmath320 ( or @xmath321 ) , we have @xmath322 for all values of the wavevector @xmath259 . in this region , we find a contribution from the @xmath0 bogoliubov modes , in addition to a smaller contribution from the atoms coming from outside the condensate region , at finite temperatures .
this contribution to the out - coupling current @xmath323 represents the `` quantum evaporation '' process discussed at the end of section [ sec : case of out - coupling to a non - bose - condensed gas ] , being associated with positive energy poles of the single - particle noncondensate green s function .
it is only present at finite temperatures . in fig .
[ fig : total out - coupling current for q=3d3d2/aho a=3d3d5 nm t=3d3d0.2tbec ] , we have plotted the _
total _ out - coupling current as a function of the detuning @xmath290 at an intermediate temperature of @xmath324 .    however , for @xmath325 ( or @xmath326 ) , @xmath327 can be less than zero for a restricted set of wavevectors @xmath259 . in this case , there are two contributions to the out - coupling current @xmath323 .
one arises from the @xmath1 bogoliubov modes and hence represents the so - called pair - breaking process associated with the negative energy poles of the single - particle noncondensate spectral density function ( i.e. ,  a tunneling atom creates a bogoliubov excitation in the trapped gas ) .
this is the only process which contributes to the noncondensate out - coupling current at zero temperature . at finite temperature ,
a second process arises which is associated with quantum evaporation from high momentum states .
this involves the positive energy poles of the single - particle noncondensate spectral density function .
this second contribution decays more quickly as the detuning is increased , compared with the @xmath1 contribution .
we note that these two regions will be doppler shifted toward positive @xmath290 by the nonzero momentum kick @xmath33 from the laser fields .
the overall behavior of our lda results in figs .
[ fig : condensate - noncondensate out - coupling current for q=3d3d2/aho and a=3d3d5nm]-[fig : total out - coupling current for q=3d3d10/aho a=3d3d5 nm t=3d3d0.2tbec ] is similar to that reported in refs .
@xcite , which were based on solving a one - dimensional gross - pitaevskii equation for an isotropic trap .
we have also evaluated the noncondensate contribution to the out - coupling current for @xmath328 , with the results shown in figs .
[ fig : condensate - noncondensate out - coupling for q=3d3d10/aho and a=3d3d5 nm ] and [ fig : total out - coupling current for q=3d3d10/aho a=3d3d5 nm t=3d3d0.2tbec ] .
we observe that the point at which there is no longer any contribution from the @xmath1 bogoliubov modes shifts increasingly toward positive @xmath290 as @xmath188 is increased .
however , unlike for the case of @xmath329 , the contribution at finite temperatures from the @xmath0 bogoliubov modes dominates that from the @xmath1 bogoliubov modes for the range of @xmath290 considered .
this results from the fact that the higher momentum quasiparticles have greater spectral weight .
we have plotted the _
total _ out - coupling current at a temperature of @xmath324 in fig .
[ fig : total out - coupling current for q=3d3d10/aho a=3d3d5 nm t=3d3d0.2tbec ] .
because of the strength of the condensate contribution , some comment is required about the observability of the noncondensate contribution in the total out - coupling current .
since the contribution from the condensate is smallest at high @xmath188 , we compare the condensate and noncondensate contributions at a momentum kick of @xmath328 ( see fig .
[ fig : condensate - noncondensate out - coupling for q=3d3d10/aho and a=3d3d5 nm ] ) . since the out - coupling current as a function of @xmath290 is much broader for the noncondensate contribution than the condensate contribution , there are regions in which the strong contribution from the condensate will not dwarf that from the noncondensate . at @xmath260 , the maximum contribution from the quantum evaporation process
is about @xmath330 of the maximum contribution from the condensate ( see fig .
[ fig : condensate - noncondensate out - coupling for q=3d3d10/aho and a=3d3d5 nm ] ) .
the small but interesting @xmath1 bogoliubov contribution is perhaps more visible as a `` distinct contribution '' in the limit of very low temperatures ( see the @xmath158 results in figs .
[ fig : condensate - noncondensate out - coupling current for q=3d3d2/aho and a=3d3d5 nm ] and [ fig : condensate - noncondensate out - coupling for q=3d3d10/aho and a=3d3d5 nm ] ) .    the unique features associated with the out - coupled atoms coming from the noncondensate
can also be enhanced by working with larger values of the s - wave scattering length @xmath12 .
in particular , one can `` tune '' the scattering length @xmath12 significantly by the use of feshbach resonances @xcite .
the amplification of the out - coupling current of atoms associated with the negative energy poles of the single - particle noncondensate green s function could be looked for in trapped gases near a feshbach resonance .
out - coupling experiments may thus provide important information about the bose - condensed gas excitation spectrum if we can increase the s - wave scattering length substantially . in this case , of course , we expect observable changes from the simple bogoliubov - popov approximation we introduced in section [ section : coupling between two trapped gases ] for the single - particle spectral density .
such many - body effects can thus be probed `` directly '' using the out - coupling kind of experiment we have been discussing .
this paper has studied the weak coupling of two trapped dilute bose gases using a tunneling hamiltonian approach .
when considering tunneling between two trapped condensates , we obtained a coherent josephson - type time - dependent contribution given in ( [ equilibrium current 3 ] ) , with a frequency @xmath331 ( see also ref .
the expression in ( [ equilibrium current 3 ] ) , which corresponds to coherent transfer of atoms between two condensates , is first order in the coupling @xmath34 .
in addition , we found that two types of quasiparticle currents emerge from treating a perturbation of the kind shown in ( [ v ] ) to second order : a static , time - independent current given by ( [ deltaina ] ) and ( [ short time normal current ] ) involving only regular correlation functions of the form @xmath332 , and a time - dependent current given by ( [ deltaiaa ] ) and ( [ deltaiab ] ) involving a second harmonic with frequency @xmath333 .
the latter is associated with anomalous correlation functions of the kind @xmath334 , and is reminiscent of the josephson current carried by cooper pairs in superconductors @xcite .
these tunneling currents were written in terms of single - particle spectral density functions that contain all the details of the microscopic correlations of the bose - condensed gases . in this manner ,
the tunneling dynamics is cleanly _ separated _ from the specific approximation used for the many - body description of the bose - condensed gas .
this separation was one of the main goals of our work .
our paper further serves to unify and extend the work of ref .
@xcite on josephson - type and other kinds of tunneling phenomena in trapped bose gases , as well as recent work @xcite dealing with out - coupling from a single bose - condensed gas .
it provides a transparent way of expressing the various contributions to the tunneling processes involving an inhomogeneous trapped bose gas , before any specific approximation is introduced .
the present paper thus provides a foundation upon which one can study the details of different tunneling processes within various many - body approximations for the dynamics of a bose gas .
our work re - emphasizes the potential usefulness of the kind of out - coupling experiments discussed in refs.@xmath3@xcite as a probe of the many - body physics of a bose - condensed gas , including systems of reduced dimensionality .
as we have noted in the text , our present study contrasts with the work reported in refs .
the latter deal with the density fluctuation spectrum while our present paper ( see also ref .
@xcite ) probes the single - particle green s function of a trapped bose gas , which is perhaps a more basic measure of the condensate - induced correlations between atoms .    for illustration ,
we have used our formalism in sections [ section : coupling between two trapped gases ] and [ sec : case of out - coupling to a non - bose - condensed gas ] to find a general expression for the out - coupling atom current using the bogoliubov - popov quasiparticle approximation .
our results exhibit the same three processes obtained by the oxford group @xcite .
the first process is the tunneling of an atom out of the condensate of the trapped gas .
the second process is the `` quantum evaporation '' of an atom from the trapped gas , _ destroying _ a quasiparticle excitation in the trapped gas and _ creating _ an atom in the out - coupled gas .
the third process , referred to as a `` pair - breaking '' process in ref .
@xcite , involves an atom tunneling out of the trapped gas concurrently with the _ creation _ of a quasiparticle excitation in the trapped gas ( see also ref .
@xcite ) .    in section [ sec : lda ]
, we have numerically evaluated the relative contributions of these three processes as a function of the laser detuning parameter @xmath335 , for various values of the wavevector difference between the two laser beams @xmath336 and the temperature @xmath273 .
this kind of experiment would be especially promising as a probe of the correlations in a bose gas in the case of a large scattering length @xmath12 made possible by working near a feshbach resonance .
this enhances the tunneling contribution from the noncondensate atoms and moreover requires the use of an improved single - particle correlation function which includes quantum depletion of the condensate ( such as the second order beliaev approximation ) .
the tunneling currents from the noncondensate atoms are admittedly small and will require high precision experiments .
we argue that the unique information that can be obtained justifies the effort . in future work
, we plan on using the formalism set up in this paper to discuss the case when the trapped bose condensate is in a vortex state .
in addition , information about the single - particle correlation functions in a one dimensional bose gas @xcite is already present in the out - coupled tunneling current
.    we would like to thank j.e .
williams for discussions of josephson tunneling in trapped bose gases , as well as m.  imamovic - tomasovic and t.  nikuni for discussions concerning the lda .
the research was supported by nserc of canada .
out - coupling atom current for @xmath337 and @xmath338 .
the detuned frequency @xmath339 is defined in ( [ definition of detuning parameter ] ) .
the plot on the left is at @xmath340 .
the plot on the right is at @xmath341 . in figs .
[ fig : condensate - noncondensate out - coupling current for q=3d3d2/aho and a=3d3d5nm]-[fig : total out - coupling current for q=3d3d10/aho a=3d3d5 nm t=3d3d0.2tbec ] , we take @xmath342 atoms , @xmath343 and @xmath344 . ]
out - coupling atom current for @xmath337 and @xmath338 .
the plot on the left shows each component contributing to the out - coupling atom current .
the plot on the right shows the total out - coupling atom current . in both plots @xmath345 . ]
out - coupling atom current for @xmath337 and @xmath346 .
the plot on the left shows each component contributing to the out - coupling atom current .
the plot on the right shows the total out - coupling atom current . in both plots @xmath345 . ]