the coherence property plays an essential role in the wave nature of bose - einstein condensates ( becs ) , and it has been investigated intensively @xcite after the experimental realizations of becs in dilute gases @xcite . for a single bec which is a perfect quantum fluid at zero temperature , the existence of a macroscopic wave function ( or order parameter ) means that spatial coherence is an intriguing property of the condensate . in
the dilute bose - condensed gases , the spatial coherence can be investigated directly by interfering two bose condensates . for two coherently separated becs , because the relative phase of two sub - condensates is locked ,
it is not surprising that there is a clear interference pattern when the two sub - condensates are allowed to overlap . in the celebrated experiment by andrews _
et al _ @xcite , however , high - contrast fringes were observed even for two completely independent condensates at an initial time . in this experiment , to prepare two initially independent condensates , the dilute bose gases were evaporatively cooled in a double - well potential created by splitting a magnetic trap in half with a far - off blue - detuned laser beam ( see also @xcite ) .
in particular , the height of the external potential due to the laser beam is much larger than the chemical potential of two separated condensates , which means that the tunneling current can be safely omitted . after switching off the double - well potential , the two initially independent condensates overlapped and high - contrast fringes
were observed in ref .
this experimental result shows clearly that there is a spatial coherence property after the overlapping between two independent interacting condensates .
although high - contrast fringes were observed for two initially independent condensates , in many literature ( see for example leggett , pethick , stringari and references therein ) , it is shown that there is no interference term in the density expectation value for two initially independent condensates . to solve the contradiction between this result and the observed high - contrast fringes @xcite for two initially independent condensates ,
the observed high - contrast fringes were interpreted with the aid of the high - order correlation function @xmath0 ( which is an oscillation function of @xmath1 ) and quantum measurement process in the present popular viewpoint .
several theories have been proposed to interpret the observed high - contrast fringes for two initially independent condensates such as the stochastic simulations of the photon detection for atoms jav , the expansion of fock state by the linear superposition of coherent states @xcite , and the continuous measurement theory @xcite .    in the present work ,
we calculate the density expectation value for two initially independent condensates by including carefully the interatomic interaction . quite different from the simple derivation in other theories
, it is found that _ for the case of two initially independent condensates , upon expansion , there is a nonzero interference term in the density expectation value when the interatomic interaction and exchange symmetry of identical bosons are both taken into account .
_ to exclude any possible error in the cumbersome derivations , we calculate the density expectation value and evolution equations by both the first - quantization and second - quantization methods , and the same results are obtained . after removing the double - well potential trapping the two initially independent condensates ,
we give the theoretical result of the density expectation value which agrees with the interference pattern observed in @xcite .
to show the interaction - induced coherence process between two initially independent bose condensates when both the interatomic interaction and exchange symmetry of identical bosons are considered carefully , the paper is organized as follows . in sec .
ii , we give the general expression of the density expectation value for two initially independent condensates based on the many - body wave function of the whole system .
it is shown that there is a nonzero interference term when the wave functions of two initially independent condensates are no more orthogonal after their overlapping . in sec .
iii , we give the general expression of the overall energy of the whole system , and the evolution equations are given by the action principle . in sec .
iv , we give a brief proof why the wave functions of two initially independent condensates become non - orthogonal after their overlapping when the interatomic interaction is considered . in this section ,
an effective order parameter is introduced to show further the interaction - induced coherence formation process for two initially independent condensates . in sec .
v , we give the numerical results for the evolution of the density expectation value according to the experimental parameters in ref .
it is shown clearly that our theoretical results of the density expectation value agree with the experimental results of the interference patterns . in sec .
vi , the evolution of the density expectation value is given for different coupling constants .
it is shown that increasing the coupling constant has the effect of enhancing the interference effect for two independent condensates . to show further the physical mechanism of the interaction - induced coherence process , in sec .
vii , we prove that the theoretical results based on the second - quantization method is the same as the results derived from the many - body wave function in the previous sections . in sec .
viii , the role of quantum fluctations and the orthogonality of the many - body wave function are discussed .
finally , we give a brief summary and discussion in the last section .
first , we give a brief introduction to the scheme of observing the interference effect of two separated condensates . as shown in fig .
1 , the double - well potential can be created by superposing a far - off blue - detuned laser beam ( which generating a repulsive optical dipole force for atoms ) upon the magnetic trap . the combined double - well potential is@xmath2 + u_{0}e^{-\left ( x - x_{0}\right ) ^{2}/w_{x}^{2}-\left ( y^{2}+z^{2}\right ) /w_{\perp } ^{2 } } , \label{external potential}\]]where the first term is the trapping potential due to the magnetic trap , while the second term represents the potential due to the laser beam .        for the bose gases confined in this double - well potential
, there are two quite different cases .
( i ) if the height @xmath3 of the external potential due to the laser beam is smaller than the chemical potential of the system , the two sub - condensates can be regarded to be coherently separated after the evaporative cooling . in this situation , the two coherently - separated sub - condensates can be regarded as a single condensate .
the relative phase of two sub - condensates is thus locked .
( ii ) if the height @xmath3 of the laser beam is much larger than the chemical potential of the system so that the tunneling effect can be omitted safely , the two condensates are completely independent after the evaporative cooling .    for two
coherently separated sub - condensates , every atom is described by the following normalization wave function @xmath4 /\sqrt{n},\]]where @xmath5 and @xmath6 are the normalization wave functions accounting for the two sub - condensates .
@xmath7 is the total number of particles . before removing the double - well potential ,
the average particle numbers in each condensate are @xmath8 and @xmath9 , respectively . after removing the double - well potential , the evolution of the wave function @xmath10
can be obtained based on the gross - pitaevskii ( gp ) equation @xcite .
the density expectation value is then @xmath11 + n_{2}\left\vert \phi _ { c2}\left ( \mathbf{r},t\right ) \right\vert ^{2}.   \label{coherent - density}\end{aligned}\]]the second term in the above equation accounts for the interference effect when there is an overlapping between two sub - condensates upon expansion .
for @xmath12 satisfying the gp equation , the nonlinear effects in the interference pattern of two coherently separated sub - condensates were investigated in refs @xcite .    as shown in fig .
1(b ) , if the intensity of the blue - detuned laser beam is sufficiently high so that the tunneling effect can be omitted , the two condensates can be regarded to be completely independent .
in this situation , the number of particles @xmath8 and @xmath9 in each of the two condensates are fixed .
we assume that there are @xmath8 particles described by the wave function @xmath13 in the left well , while there are @xmath9 particles described by the wave function @xmath14 in the right well . to investigate clearly the role of the exchange symmetry of identical particles and interatomic interaction , we calculate in the following the density expectation value @xmath15 directly from the many - body wave function . in sec .
vii , we will also calculate the density expectation value based on the second - quantization method , and the same result is obtained .
we will prove in sec .
iv that in the presence of the interatomic interaction , @xmath16 and @xmath17 will become non - orthogonal after the overlapping between two initially independent condensates .
thus we consider here the following general case from the beginning@xmath18after removing the double - well potential , the two initially independent condensates will overlap , and thus one should consider the indistinguishability of identical particles .
when the exchange symmetry of identical bosons is considered , the many - body wave function is@xmath19 ,   \label{wave - function}\end{aligned}\]]where @xmath20 denotes the @xmath21 permutations for the particles in different single - particle state @xmath13 or @xmath14 .
@xmath22 is a normalization factor to assure @xmath23 .
@xmath22 is determined by the following equation : @xmath24 ^{1/2}=1.\]]in this paper , to give a concise expression for various coefficients such as @xmath22 , we have introduced the rule @xmath25 .
when the non - orthogonal property between @xmath13 and @xmath14 is considered , one should note that the normalization constant @xmath22 is relevant to the parameter @xmath26 . from the form of the many - body wave function ( [ wave - function ] ) , the quantum depletion originating from interparticle interaction is omitted .
thus , this form of quantum state is valid when @xmath27 with @xmath28 and @xmath29 being respectively the scattering length and mean distance between particles .
the role of quantum depletion will be discussed in sec .
viii .    from the above many - body wave function , after straightforward derivations , the exact expression of the density expectation value takes the following form : @xmath30 + c_{d}\left\vert \phi _ { 2}\left ( \mathbf{r},t\right ) \right\vert ^{2 } , \label{density}\end{aligned}\]]where the coefficients are@xmath31 in the above summations,@xmath32    for two independent ideal condensates , before the overlapping between the two condensates , we have @xmath33 . based on the schrdinger equation , it is easy to verify that after the double - well potential separating the condensates is removed , we have @xmath34 at any further time .
thus @xmath35 , and the density expectation value is given by@xmath36 in this situation , the interference term is zero in the density expectation value .    in the presence of the interatomic interaction ,
@xmath37 can be a nonzero value after the two initially independent condensates begin to overlap . for @xmath37 being nonzero
, we see clearly from the coefficient @xmath38 that the interference term in ( [ density ] ) gives a contribution to the density expectation value .
based on eqs .
( [ adi ] ) , ( [ bdi ] ) and ( [ cdi ] ) , we have@xmath39and@xmath40it is easy to understand that when @xmath41 and @xmath42 , @xmath43 can not be omitted , and thus there would be clear interference patterns .
this shows clearly that for large particle number , a small @xmath44 can play important role in the density expectation value when @xmath45 and @xmath46 .    generally speaking , even in the presence of the interatomic interaction
, @xmath44 is much smaller than @xmath47 because @xmath48 is an oscillation function about the space coordinate . as shown above , however
, a nonzero value of @xmath49 can give significant contribution to the density expectation value for large @xmath8 and @xmath9 . in fig .
2(a ) , we give the relation between @xmath50 and @xmath26 for @xmath51 based on eqs .
( [ bd ] ) and ( [ cd ] ) .
the relation between @xmath50 and @xmath52 for @xmath53 is shown in fig .
generally speaking , for @xmath54 and @xmath55 , one has @xmath56 .    ) , ( [ bd ] ) and ( [ cd ] ) , shown in fig .
2(a ) is the relation between @xmath50 and @xmath57 , while shown in fig .
2(b ) is the relation between @xmath50 and @xmath58 for @xmath59 .
it is shown clearly that the interference term can play an important role in the density expectation value for @xmath60 being larger than @xmath61 . ]
for the many - body wave function @xmath62 , in the presence of the interatomic interaction , the evolution equation can be obtained based on the standard quantum mechanical principle .
after the double - well potential is removed , the evolution equation is given by@xmath63where the hamiltonian of the whole system in the first - quantization method is@xmath64here the coupling constant @xmath65 with @xmath28 being the scattering length . in the above expression of the hamiltonian , we have used the well - known two - body pseudopotentials .
it is well known that the action principle is quite useful to derive the gp equation for a single condensate ( see for example ref .
similarly , we consider here the evolution of @xmath62 based on the action principle . to get the evolution equation based on the action principle
, we first give the general expression of the overall energy of the whole system . after removing the double - well potential ,
the overall energy of the whole system is @xmath66after straightforward derivations , the exact expression of the overall energy is given by@xmath67where the kinetic energy @xmath68 is@xmath69 in addition , the interaction energy @xmath70 of the whole system is given by@xmath71 , \label{exact - inter}\end{aligned}\]]where the coefficients are given by@xmath72    by using the ordinary action principle and the energy of the whole system , one can get the following coupled evolution equations for @xmath13 and @xmath73 : @xmath74where @xmath75 and @xmath76 are functional derivatives .
ii , we have shown that the nonzero interference term in the density expectation value origins from the assumption that @xmath37 can be a nonzero value after the overlapping between two initially independent condensates . in the last section ,
we have given the evolution equations about @xmath13 and @xmath14 .
thus , an important question emerges naturally : is it physical for @xmath37 being nonzero with the development of time ?    based on the evolution equations ( [ gp1 ] ) and ( [ gp2 ] ) given in the last section , one can understand easily that @xmath37 becomes nonzero after the overlapping between two initially independent condensates in the presence of the interatomic interaction . before removing the double - well potential , @xmath79 and @xmath80 because there is no overlapping between @xmath13 and @xmath14 .
after removing the double - well potential and at the beginning of the overlapping between two condensates , @xmath81 and @xmath82 . for @xmath83 and @xmath84 ,
the overall energy can be approximated as@xmath85the above overall energy is obtained by setting @xmath86 in eq .
( s - overallenergy ) .    based on eqs .
( [ gp1 ] ) and ( [ gp2 ] ) , the approximate evolution equations for @xmath16 and @xmath17 are then @xmath87and @xmath88where @xmath89and @xmath90    we see that @xmath91 is not equal to @xmath92 for @xmath93 .
this leads to an important result that @xmath37 can be nonzero when there is an overlapping between two initially independent condensates .
one can understand this result further through the following equation which determines the evolution of @xmath37 : @xmath94where @xmath95 is a nonzero function given by @xmath96    the above analyses show clearly that why @xmath37 becomes nonzero after the overlapping between two initially independent condensates for @xmath93 . obviously , for @xmath97 being zero , @xmath34 at any further time because @xmath98 .
our numerical calculations in the following section also show that @xmath37 can be a nonzero value in the presence of the interatomic interaction .    in the presence of the interatomic interaction ,
it is inconsistent to assume that @xmath37 is always zero with the development of time .
if @xmath34 , the evolution equations eqs .
( [ gp1 ] ) and ( [ gp2 ] ) are exact .
however , based on these two evolution equations , @xmath99 will increase from zero after the overlapping .
as shown above , the evolution of @xmath13 and @xmath14 is determined by the expression of the overall energy and the evolution equations ( gp1 ) and ( [ gp2 ] ) .
it seems that there is no exact solution for these evolution equations considering the fact that even there is no exact solution for three - dimensional nonlinear schrdinger equation .
however , for @xmath100 and @xmath101 , we find that there is a quite simple evolution equation by introducing an effective parameter order .
this would also contribute to our understanding of the interaction - induced coherence process for two initially independent condensates .    for the cases of @xmath54 , @xmath102 and @xmath103 ,
first we introduce the following effective order parameter @xmath104 which is given by @xmath105based on this effective order parameter , the density expectation value can be approximated well as@xmath106 + c_{d}^{\prime } \left\vert \phi _ { 2}\left ( \mathbf{r},t\right ) \right\vert ^{2 } ,   \label{app - density}\end{aligned}\]]where@xmath107    to compare with the exact expression of the density expectation value given by eq .
( [ density ] ) , fig .
3 shows the ratio @xmath108 , @xmath109 and @xmath110 .
we see that for @xmath111 , @xmath112 , the approximate density expectation value @xmath113 agrees very well with the exact expression of the density expectation value .    ,
@xmath114 and @xmath115 .
we see that for @xmath116 , @xmath117 and @xmath118 , the density expectation value given by the effective order parameter agrees well with the exact expression of the density expectation value given by eq .
( [ density ] ) . ]    based on the effective order parameter , for @xmath111 , @xmath112 and @xmath119 , we also find that the overall energy of the whole system can be approximated very well as @xmath120where@xmath121and@xmath122    for the cases of @xmath54 , @xmath102 and @xmath103 , it is easy to verify that @xmath123 based on the analogous analyses about the density expectation value . for @xmath54 , @xmath124 and @xmath103 , one can also prove the result of @xmath125
. based on eq .
( eappr - inter ) , @xmath126 can be expanded as : @xmath127 , \end{aligned}\]]where @xmath128to compare with the exact expression of the interaction energy given by eq .
( [ exact - inter ] ) , fig .
4 shows the relation between @xmath129 ( @xmath130 ) and @xmath26 for @xmath51 .
it is shown clearly that for @xmath54 and @xmath55 , @xmath131 , and thus @xmath132 .
in this situation , the overall energy can be approximated well as @xmath133the evolution of the effective order parameter can be then obtained by using the action principle .
based on this approximate energy , after removing the double - well potential , it is quite interesting to note that the evolution of the effective order parameter can be described very well by the ordinary gross - pitaevskii equation : @xmath134     ( @xmath130 ) and @xmath57 for @xmath51 .
we see that for @xmath135 and @xmath136 , @xmath137 , which means that @xmath125 . ]
based on eqs .
( [ app - density ] ) and ( [ s - app - equation ] ) , we see that the emergence of the effective order parameter @xmath138 gives us strong evidence that the coherence is formed in the interaction process between two initially independent condensates , and thus results in the emergence of high - contrast interference fringes .
the effective order parameter and the approximate density expectation value suggest strongly that a full coherence is formed between two initially independent condensates for the cases of @xmath139 and @xmath55 .
now we turn to give the theoretical results of the density expectation value according to the experimental parameters in ref .
@xcite where clear interference patterns were observed for two initially independent condensates . in the experiment of ref .
@xcite , @xmath140 condensed sodium atoms were confined in a magnetic trap with @xmath141 , @xmath142 .
a blue - detuned laser beam of wavelength @xmath143 was focused into a light sheet with a cross section of @xmath144 @xmath145 by @xmath146 @xmath145 .
the long axis of the laser beam was perpendicular to the long @xmath147axis of the condensate . for a laser power of @xmath148 ,
the barrier height is about @xmath149 @xmath150 , which is much larger than the chemical potential @xmath151 @xmath150 . for this laser beam
, the two condensates can be regarded to be independent because they are well separated and the tunneling effect can be omitted . with these experimental parameters and the @xmath152wave scattering length @xmath153 , the initial profile of the two condensates
is shown in fig .
5 . at the initial time
, the overlapping between two condensates can be omitted safely , and thus @xmath33 .
after the double - well potential is removed , the evolution of the density expectation value @xmath154 ( in unit of @xmath155 ) is given in fig . 5 through the numerical calculations of eqs . ( [ gp1 ] ) , ( [ gp2 ] ) and ( [ density ] ) .
we see that there is a clear interference pattern in the density expectation value which agrees with the experimental result . shown in the inset
is the evolution of @xmath156 for these parameters .
we see that @xmath157 increases from zero with the development of time . for the expansion time of @xmath158 ,
the numerical result of @xmath156 shows that @xmath159 . to check further
that the nonzero value of @xmath160 does not originate from numerical error , with the same initial conditions , we have verified in the numerical calculations that @xmath160 is always zero ( smaller than @xmath161 ) if the scattering length is assumed as zero .     for the experimental parameters in ref .
]    based on the evolution of the density expectation value shown in fig . 5 , we see that the overall width ( about @xmath162 @xmath145 ) of the ultracold gases in the @xmath147direction does not increase obviously .
this is due to the fact that the initial density distribution is cigar - shaped , and thus the expansion in the @xmath147direction is very slow , while the expansion in @xmath163 and @xmath164 directions is much quick . in the experiment of ref .
@xcite , one can see clearly that there is no obvious expansion in the long @xmath147axis for cigar - shaped condensate . for two initially independent condensates , when the double - well potential is switched off , one should note that in the regime close to @xmath165 @xmath145 shown in fig .
5 , the ultracold gases expand rapidly in the @xmath147direction because in this regime the ultracold gases have higher kinetic energy .
thus , although the total width of the system in the @xmath147direction does not increase obviously , the rapid expansion in the central regime leads to the overlapping between two initially independent condensates , and results in the interference effect .
the overall width of the system shown in fig .
5 is smaller than the experimental result of about @xmath166 @xmath145 @xcite .
this difference may come from the expansion of thermal cloud in this experiment @xcite .
after @xmath158 expansion , the numerical result in @xcite for two coherently separated condensates also showed that the overall width of the system is about @xmath162 @xmath145 .
we see that the interatomic interaction plays an essential role in the emergence of the interference effect for two initially independent condensates . generally speaking , increasing
the particle number will enhance the effect of the interference term in the density expectation value .
based on eqs .
( [ gp1 ] ) and ( [ gp2 ] ) , increasing the coupling constant @xmath97has the effect of increasing @xmath37 .
together with the relation between @xmath50  and @xmath26  illustrated in fig .
2(a ) , this shows that increasing the interatomic interaction will enhance the effect of the interference term . to show more clearly the interaction - induced coherence process between two initially independent condensates , in this section we consider the density expectation value for different coupling constants .
we consider here the evolution of the density expectation value for one - dimensional case . at @xmath167 , to give a general comparison , the initial wave functions for two independent condensates are assumed to be identical for different coupling constants .
the initial wave functions are respectively given by @xmath168 ,   \label{inital1 } \\
\phi _ { 2}\left ( x_{l},t=0\right ) & = & \frac{1}{\pi ^{1/4}\sqrt{\delta _ { 2}}}% \exp \left [ -\frac{\left ( x_{l}-x_{l2}\right ) ^{2}}{2\delta _ { 2}^{2}}\right ] .
\label{inital2}\end{aligned}\]]in the above wave functions , we have introduced a dimensionless variable @xmath169 with @xmath170 being a length . in the present work
, we assume that @xmath171 and @xmath172 .
for these parameters , at @xmath167 , the two condensates are well separated . in the numerical calculations of the coupled equations given by eqs .
( [ gp1 ] ) and ( gp2 ) , it is useful to introduce the dimensionless variable @xmath173 with @xmath174 , and dimensionless coupling constant @xmath175 .
in addition , the particle number is assumed as @xmath176 . in real experiments ,
interatomic interaction plays very important role in the initial ground - state wave function of the condensates .
however , in principle , one can prepare the state given by eqs .
( [ inital1 ] ) and ( [ inital2 ] ) by adjusting the trapping potential for different coupling constants . in this section ,
the identical initial wave functions for different coupling constants would be helpful in the comparison of the density expectation value for different coupling constants .    for @xmath177 , we consider the evolution of the density expectation value in free space . the evolution of @xmath13 and @xmath14 is obtained based on the numerical calculations of eqs .
( [ gp1 ] ) and ( [ gp2 ] ) . from @xmath178 and @xmath14
, we can get @xmath26 , and thus the density expectation value based on eq .
( [ density ] ) . shown in fig .
6 is the evolution of the density expectation value for different coupling constants .
it is shown clearly that increasing the coupling constant has the effect of enhancing the coherence effect , and results in higher contrast in the interference patterns .
( in unit of @xmath179 ) for different coupling constants @xmath180 . shown in the inset of each figure is the relation between @xmath181 and dimensionless time @xmath182 .
for two ideal condensates shown in fig .
6(a ) , we see that there is no interference pattern even there is an overlapping between two condensates . for the case of @xmath183 shown in fig .
6(b ) , we see that low - contrast interference patterns begin to emerge due to the interaction - induced coherence process . in fig .
6(c ) and fig .
6(d ) , we see that there are high - contrast interference patterns . in particular , in fig .
6(d ) for @xmath184 , two initially independent condensates can be regarded to be fully coherent near @xmath185 . ]
the essential reason for the emergence of the interference term of two initially independent condensates lies in that because of the exchange symmetry of identical bosons and interatomic interaction , the two initially independent condensates become coherent after the overlapping between the two condensates .
the physical mechanism of this interaction - induced coherence can be understood further based on the second quantization method .
thus , in this section , we give the density expectation value and evolution equations based on the second - quantization method .
the merit of the second - quantization method lies in that the indistinguishability of identical bosons is satisfied when the correct commutation relation of the field operator is used .      for two initially independent condensates comprising particle number @xmath8 and @xmath9 ,
the corresponding quantum state is ( see for example refs . @xcite):@xmath186where @xmath187 is a normalization constant to assure @xmath188 .
@xmath189 ( @xmath190 ) is a creation operator which creates a particle described by the single - particle state @xmath13 ( @xmath14 ) in the left ( right ) condensate .
similarly to eq .
( [ wave - function ] ) , the quantum depletion is omitted in this sort of quantum state .    for two initially independent condensates , for the state given by eq .
( initial - state ) , it seems that there is no interference term in the density expectation value with the following simple calculation : @xmath191 in the above equation , the field operator @xmath192 is expanded as @xmath193 with @xmath194 and @xmath195 being the annihilation operators .
one should note that , to get ( [ nd - field(2 ) ] ) from ( nd - center ) , there is an implicit assumption that @xmath194 and @xmath190 are commutative .
this holds when @xmath196 . when @xmath197 = 0 $ ] , it is easy to understand that the interference term ( the last term in eq .
( nd - center ) ) is zero in @xmath15 .
as shown in the previous sections , @xmath13 and @xmath14 will become non - orthogonal in the presence of the interatomic interaction .
the operators @xmath194 and @xmath195 can be written as@xmath198and@xmath199here @xmath192 is the field operator . by using the commutation relations of the field operators @xmath200=0
$ ] and @xmath201=\delta \left ( \mathbf{r}_{1}-\mathbf{r}% _ { 2}\right ) $ ] , it is easy to get the commutation relation@xmath202=\zeta ^{\ast } .
\label{commre}\]]we see that @xmath194 and @xmath203 are not commutative any more for @xmath204 being a nonzero value .
in this situation , it is obvious that one can not get the result ( nd - field(2 ) ) from ( [ nd - center ] ) any more .
this means that one should be very careful to get the correct density expectation value for two initially independent condensates .
it is well - known that the field operator should be expanded in terms of a complete and orthogonal basis set .
generally speaking , the field operator @xmath205 can be expanded as : @xmath206where @xmath13 and @xmath207 are orthogonal normalization wave functions .
assuming that @xmath208 , based on the conditions @xmath209 and @xmath210 , we have @xmath211 and @xmath212 @xmath213 . based on@xmath214we have@xmath215it is easy to get the following commutation relations : @xmath216 & = & [ \widehat{k}^{\dagger } , \widehat{k}% ^{\dagger } ] = 0,[\widehat{k},\widehat{k}^{\dagger } ] = 1 ,   \nonumber \\ \left [ \widehat{a}_{1},\widehat{a}_{1}\right ] & = & [ \widehat{a}_{1}^{\dagger } , % \widehat{a}_{1}^{\dagger } ] = 0,[\widehat{a}_{1},\widehat{a}_{1}^{\dagger } ] = 1 , \nonumber \\
\lbrack \widehat{k},\widehat{a}_{1 } ] & = & [ \widehat{k},\widehat{a}% _ { 1}^{\dagger } ] = 0 .
\label{kcom}\end{aligned}\ ] ]    because @xmath217 and @xmath218 are commutative , it is convenient to calculate the density expectation value @xmath219 by using the operators @xmath217 and @xmath220 . after straightforward derivations ,
the exact expression of the density expectation value is @xmath221 ^{n_{2}}\widehat{a}_{1}^{n_{1}}\left [ \widehat{a}% _ { 1}^{\dag } \phi _ { 1}^{\ast } + \widehat{k}^{\dag } \left ( \phi _ { 2}^{\prime } \right ) ^{\ast } \right ]   \nonumber \\ & & \left [ \widehat{a}_{1}\phi _ { 1}+\widehat{k}\phi _ { 2}^{\prime } \right ] \left ( \widehat{a}_{1}^{\dagger } \right ) ^{n_{1}}\left [ \frac{\widehat{k}% ^{\dagger } } { \beta } + \zeta ^{\ast } \widehat{a}_{1}^{\dagger } \right ] ^{n_{2}}\left\vert 0\right\rangle   \nonumber \\ & = & \alpha _ { d}\left\vert \phi _ { 1}\left ( \mathbf{r},t\right ) \right\vert ^{2}+2\beta _ { d}\times \mathrm{re}\left ( e^{i\varphi _ { c}}\phi _ { 1}^{\ast } \left ( \mathbf{r},t\right ) \phi _ { 2}\left ( \mathbf{r},t\right ) \right ) + \gamma _ { d}\left\vert \phi _ { 2}\left ( \mathbf{r},t\right ) \right\vert ^{2 } , \label{ndensity}\end{aligned}\]]where the coefficients are @xmath222 in addition , the normalization constant @xmath187 is determined by@xmath223    the above density expectation value is obtained based on the second quantization method .
although it seems that the expressions of the coefficients given by eqs .
( [ ads ] ) , ( [ bds ] ) , ( [ cds ] ) and ( normconsts ) are quite different from the results calculated from the many - body wave function , we have proven that @xmath224 , @xmath225 , @xmath226 and @xmath227 . thus , the density expectation value given by eq .
( [ ndensity ] ) is the same as the result calculated from the many - body wave function @xmath62 which satisfies the exchange symmetry of identical bosons .
based on the second quantization method , after removing the double - well potential , the hamiltonian in the second - quantization method is@xmath228after straightforward calculations , the overall energy of the whole system is @xmath229where the kinetic energy @xmath68 is given by@xmath230 in addition , the interaction energy @xmath70 of the whole system is given by@xmath231 , \end{aligned}\]]where the coefficients are given by@xmath232we have proven that the above overall energy is equal to the results based on the many - body wave function , by checking that @xmath233 for @xmath234 , @xmath235 , @xmath236 , @xmath237 .    when the non - orthogonal property between @xmath13 and @xmath14 are considered , because the derivations of the density expectation value and overall energy are quite cumbersome , the same results based on the first - quantization method and second - quantization method give us strong evidence that our derivations are correct .
in both the first - quantization and second - quantization methods of the previous calculations , the quantum depletion originating from the elementary excitations at zero temperature is omitted in the quantum state of the whole system .
based on the bogoliubov theory of the elementary excitations , for the initial quantum state , the number of particles due to the quantum depletion is of the order of @xmath238 and thus the quantum depletion is negligible for bose condensate in dilute gases considered in the present work . with the development of time ,
the role of quantum depletion can also be omitted .
because of the factor @xmath239 ( @xmath240 is the wave number of the elementary excitations ) in the wave function @xmath241 of the elementary excitations , a simple analysis shows qualitatively that @xmath242 and @xmath243 are of the order of @xmath244 with @xmath245 being the spatial size of the system .
this exponential decay of @xmath246 and @xmath247 originates from the integral where there is spatially oscillating phase factor in the wave functions of the elementary excitations and condensates .
thus , the contribution to the effective order parameter and density expectation value due to elementary excitations can be omitted safely .    in the preceding paragraph
, it is shown that in calculating the evolution of the density expectation value , the role of quantum depletion can be omitted
. however , the quantum depletion plays important role in the consistency of our theory .
here we discuss mainly the orthogonality for the quantum state of the whole system , especially about the subtle problem that whether the nonorthogonality between @xmath13 and @xmath14 violates the orthogonality of the whole system which must be satisfied .
this sort of problem about orthogonality exists also for a single condensate . for a single condensate , if the quantum depletion is omitted , the many - body wave function is@xmath248based on the many - body schrdinger equation and the action principle , it is easy to get the following gp equation@xmath249assume that there are two orthogonally initial quantum states @xmath250 and @xmath251 . with the development of time , because of the nonlinearity of the gp equation , @xmath252 and @xmath253 may become non - orthogonal in the presence of interatomic interaction .
it is obvious that @xmath254 ^{n}$ ] .
although for large @xmath255 , this integral can be approximated as zero , based on the consideration of consistency , the physical mechanism of the exact orthogonality between @xmath256 and @xmath257 is an interesting problem .
we think that the effect of quantum depletion omitted in eq .
( many - body - single ) is the essential reason for this sort of unphysical nonorthogonality . for the hamiltonian ( [ hamiltonian - first ] ) of the whole system ,
the term @xmath258 is non - factorable about the coordinate , and thus the exact solution of the many - body wave function is non - factorable too . for a single condensate , the exact many - body wave function
can be assumed as@xmath259here @xmath260 accounts for the non - factorable component , and @xmath261 is a normalization constant . in the hard - sphere approximation , @xmath262 if @xmath263 for any @xmath264 ( @xmath265 is the hard - sphere radius ) , and @xmath266 for otherwise situation .
@xmath267 represents the quantum depletion , and was successfully used to calculate the quantum depletion of superfluid liquid @xmath268 in ref .
the omission of the quantum depletion means that we approximate @xmath260 as @xmath47 . for dilute bose condensed gases , @xmath269 .
thus , the wave function ( [ many - body - single ] ) can describe well the condensate .
when the orthogonality of the whole system is considered , we stress here that the wave function @xmath270 should be used . the non - factorable factor
@xmath260 assures the orthogonality of the many - body quantum state , because @xmath270 is the exact solution of the many - body schrdinger equation .
for dilute bose condensed gases , these analyses lead to two results : when the dynamic evolution of dilute bose condensed gases is considered , omitting the quantum depletion can give us quite good description ; when the consistency especially the orthogonality is considered , one should consider the role of quantum depletion .    it is natural to generalize the above analyses to two initially independent condensates .
the nonorthogonality between @xmath13 and @xmath14 does not mean in any sense the violation of the orthogonality of the quantum state of the whole system .
we stress here again that when orthogonality is considered , we should check whether the quantum state of the whole system satisfies the orthogonality .
when the quantum depletion is considered , the many - body wave function takes the following form@xmath271here @xmath62 is given by eq .
( [ wave - function ] ) .
@xmath272 satisfies exactly the many - body schrdinger equation , and thus this sort of wave function satisfies the orthogonality condition .
assume that @xmath273 is constructed by @xmath274 and @xmath275 , while @xmath276 is constructed by @xmath277 and @xmath278 .
assume further that @xmath274 , @xmath279 , @xmath280 and @xmath278 are orthogonal to each other at the initial time , so that @xmath273 and @xmath276 are orthogonal at the initial time .
omitting the quantum depletion and roughly speaking , the integral between @xmath281 and @xmath282 is of the order of @xmath283 ^{\alpha } \left [ \int \left ( \phi _ { 1}^{a}\right ) ^{\ast } \phi _ { 2}^{b}d\mathbf{r}\right ] ^{\beta } \left [ \int \left ( \phi _ { 2}^{a}\right ) ^{\ast } \phi _ { 1}^{b}d% \mathbf{r}\right ] ^{\gamma } \left [ \int \left ( \phi _ { 2}^{a}\right ) ^{\ast } \phi _ { 2}^{b}d\mathbf{r}\right ] ^{\delta } $ ] ( here @xmath284 ) .
we see that similarly to a single condensate , the orthogonality between @xmath281 and @xmath282 is quite good for large @xmath255 when the quantum depletion is omitted . but to assure the exact orthogonality of the whole quantum state , we must consider the role of quantum depletion . although we do not provide here a fully rigorous proof , these analyses lead to the result that : although for dilute bose condensed gases , the omission of the quantum depletion will not play important contribution to the evolution of @xmath285 and @xmath14 , it s the quantum depletion that assure the orthogonality of the whole quantum state .
up to date , we mainly use the bases @xmath13 and @xmath14 to study the dynamic evolution and density expectation value . considering the nonorthogonality between @xmath13 and @xmath14 in the presence of interparticle interaction ,
it is attractive to discuss the physical picture using the orthogonal bases @xmath13 and @xmath207 introduced in sec . vii .
if we use the orthogonal bases @xmath13 and @xmath286 , the quantum state of the whole system is @xmath287 .
we see that the number of particles in the orthogonal modes @xmath13 and @xmath207 are no more definite .
this quantum state becomes a superposition of different number of particles in the orthogonal modes @xmath13 and @xmath207 .
this result is very natural because the coupling ( interaction ) between two initially independent condensates leads to the coherent transfer of particles between two condensates . in ref .
@xcite , it was shown that for a single condensate , a coherent state description of the bose condensed system is a robust state in the presence of the interactions between the condensate and its environment .
in the present work , our research shows that the interaction between two condensates makes each condensate become a coherent superposition of different particle number . the results in ref .
barnett about the robustness of the coherent state imply that @xmath288 is a robust quantum state even when the coupling with the environment exists .
in summary , upon expansion , we calculate the density expectation value of two initially independent condensates . it is found that there is a nonzero interference term when the interatomic interaction and the exchange symmetry of identical bosons are both considered carefully .
in fact , it is well - known that the interaction plays an essential role in the formation of the order parameter of bose - condensed gases , i.e. the formation of a stable coherent property . here
, we provide an example in which the interaction induces coherent evolution between two initially independent condensates .    in the present popular viewpoint , the high - order correlation function and quantum measurement process are used to interpret the interference patterns observed for two initially independent condensates .
although we find that there is already high - contrast interference patterns in the density expectation value for the experimental parameters in ref .
@xcite , it is still possible that the high - order correlation function and thus the mechanism of the interference due to the measurement process play an important role in the emergence of the interference patterns .
thus , to investigate more clearly the interference pattern due to the measurement process alone , we believe an experimental investigation of two independent ideal condensates would be very interesting , because there is no interference term for ideal condensates in the density expectation value . in the last few years ,
the rapid experimental advances of feshbach resonance where the scattering length can be tuned from positive to negative make this sort of experiment be feasible .
note added : the present paper is an expansion of our two previous unpublished works @xcite .
in addition , after the submission of the present paper , we noticed a paper @xcite by l. s. cederbaum et al .
, where the role of interaction in the emergence of interference patterns for two initially independent bose condensates is stressed .
ketterle w , durfee d s and stamper - kurn d m 1999 in _ making , probing and understanding bose - einstein condensates , _ proceedings of the international school of physics , enrico fermi  course cxl , edited by inguscio m , stringari s and wieman c e ( ios press : amsterdam ) , pp .
67359                                  xiong h w , liu s j and zhan m s 2005 _ preprint _ cond - mat/0507354 xiong h w and liu s j 2005 _ preprint _ cond - mat/0507509 cederbaum l s , streltsov a i , band y b and alon o e 2006 _ preprint _ cond - mat/0607556