in the currently popular hierarchical picture of structure formation , galaxies are thought to be embedded in massive dark halos .
these halos grow from density fluctuations in the early universe and initially contain baryons in a hot gaseous phase .
this gas subsequently cools and some fraction eventually condenses into stars .
much progress has been made in observationally delineating the global star formation history and the resulting build - up of stellar mass ( e.g. , madau et al .
1998 , brinchmann & ellis 2000 ; dickinson et al . 2003
; bundy , ellis & conselice 2005 ) . however , many of the physical details , particularly the roles played by feedback and cooling , essential for a full understanding of how galaxies form , remain uncertain . models ( e.g. , van den bosch 2002 ; abadi et al .
2003 ) have great predictive power in this area but only by assuming presently - untested prescriptions for these effects .
obtaining further insight into how such processes operate is thus an important next step not only in understanding galaxy evolution , but also in verifying the utility of popular models , as well as the hierarchical concept itself .
one approach towards understanding this issue is to trace how the stellar mass in galaxies forms in tandem , or otherwise , with its dark mass .
the first step in this direction began with studies of scaling relations between the measurable properties of disk galaxies , specifically the relation between luminosity and maximum rotational velocity ( tully & fisher 1977 ) .
studies utilizing roughly a thousand spiral galaxies have revealed a tight correlation between absolute magnitude and the maximum rotational velocity for nearby galaxies ( haynes et al .
1999 ) . the limited data at high redshift suggests the tf relation evolves only modestly , equivalent to at most 0.4 - 1 magnitudes of luminosity evolution to @xmath4 ( vogt et al .
1997 ; ziegler et al . 2002 ; bohm et al .
how the tully - fisher relation evolves with redshift is still controversial , although it appears that fainter disks evolve the most ( bohm et al . 2004 ) , and that selection effects are likely dominating the differences found between various studies .
furthermore , it has been difficult for modelers to reproduce the tully - fisher relation to within 30% ( e.g. , cole et al .
2000 ) , making it an important constraint on our understanding of the physics behind galaxy formation .
unfortunately , any interpretation of the tf relation is complicated by the fact that both luminosity and virial mass might be evolving together .
a more physically - motivated comparison would be between stellar and virial mass .
not only does this relation break potential degeneracies in the tf technique , but it also samples more fundamental quantities . in this paper
we begin this task by investigating the evolution in the fraction of the total mass that is in a stellar form .
this can be accomplished with some uncertainty by contrasting the _ stellar mass _ of a galaxy with its _
halo mass_. we selected disk galaxies for our initial study since these two quantities can be effectively probed observationally for these galaxies with various assumptions ( e.g. , van den bosch 2002 ; baugh et al .
2005 ) .
this paper presents the first study of the near - ir tf relation , as well as a comparison between stellar and halo masses , for 101 disk galaxies within the redshift range @xmath5 drawn mostly from the deep1 redshift survey ( vogt et al .
our goal is to address several questions relating to the mass assembly history of disks .
as our sample is not formally complete in any sense , we can not derive general conclusions concerning the history of _ all _ present - day disks .
however , we can determine whether the disks selected from the deep1 survey in the sampled redshift range are still accreting matter and converting baryons into stellar disks at a significant rate .
we construct stellar mass tully - fisher and stellar mass / halo mass relation for our sample and find that there is little evolution in either from @xmath6 .
this suggests that the dark and stellar components of disk galaxies grow together during this time .
the paper is organized as follows :  2 contains a description of the sample including the fields used and the different data products and a discussion on quality control .
3 describes how various quantities , such as the halo and stellar masses , are derived from the basic data .
 4 presents our results and  5 summarizes our conclusions .
we assume the following cosmology throughout this paper : h@xmath7 = 70 km s@xmath8 mpc@xmath8 , @xmath9 = 0.7 and @xmath10 = 0.3 .
our sample consists of 101 galaxies , 93 of which are drawn from the deep1 survey ( see vogt et al .
2004 ; simard et al . 2002 ; weiner et al .
2004 ) , and 8 of which were obtained independently .
each of these systems has a resolved rotation curve which was obtained with the keck telescope using lris ( oke et al .
1995 ) in the redshift range @xmath11 to @xmath12 .
the additional data presented in this paper , which enables stellar masses to be compared with virial masses , consists of near - infrared imaging described in @xmath132.2 .
full details of the selection criteria for the galaxies are discussed in deep1 papers ( vogt et al . 1996 , 1997 , 2005 )
to which the interested reader is referred . moreover , the necessary assumptions implicit in the derivation are detailed in those papers , to which we refer the reader .
briefly , galaxies were selected morphologically as elongated disks in hubble space telescope ( hst ) f814w ( i@xmath14 ) images with i@xmath15 .
the inferred inclination was chosen to be greater than 30to facilitate a measurement of the rotational velocity .
the optical images used for both photometric and morphological analyses come from hst wide field planetary camera-2 ( wfpc2 ) observations of the groth strip ( groth et al .
1994 ; vogt et al .
2004 ) , the hubble deep field ( williams et al .
1996 ) , and cfrs fields ( brinchmann et al . 1998 ) .
using i@xmath14 images , structural parameters were determined for each galaxy using the gim2d and galfit packages ( simard et al .
2002 ; peng et al .
we fit a two - component model to the surface brightness distribution , assuming a de vaucouleurs law for the bulge , and an exponential for the disk component . based on these fits , the disk scale length , r@xmath16 , and bulge - to - disk ratio ( @xmath17 )
were determined .
the uncertainties in the r@xmath16 values determined through this method are possibly underestimated , as has been explored through careful 2-d fitting ( e.g. jong et al 1996 ) . to be inclusive of possible effects , we will incorporate an additional 30% uncertainty for our overall error on the measured values of @xmath18 .
details of the observations , reductions and extraction of maximal rotational velocities @xmath19 from the lris spectroscopy are presented in vogt et al .
( 1996 , 1997 , 2005 ) .
briefly , each disk galaxy was observed along its major axis , as determined from the hst images ( simard et al .
the maximum velocity , @xmath19 , is determined by fitting a fixed form for the rotation curve scaled according to the i@xmath14 disk scale length @xmath18 .
the assumed rotation curve has a linear form which rises to a maximum at 1.5@xmath20 r@xmath21 , and remains flat at larger radii .
our assumption that 1.5@xmath20 r@xmath16 is the radius where the rotation curves for disks reach their maximum is reasonable , based on a similar behavior for local disk systems of similar magnitudes ( persic & salucci 1991 ; sofue & rubin 2001 ) .
our rotation curves are also visible out to several scale lengths , or roughly 2 - 5 ( e.g. , vogt et al .
1997 ) , adequate for measuring v@xmath22 .
this model form is then convolved with a seeing profile that simulates the conditions under which the observations were taken , and v@xmath23 is determined by iterative fitting .
effects from the width of the slit , slit misalignment with the major axis , and varying inclination were taken into account when performing these fits , and calculating the resulting errors .
typically , @xmath19 is determined to a precision of 10 - 20% due to these various sources ( e.g. , vogt et al .
1997 ) .
the new data we present in this paper consist of deep near - infrared observations of the deep1 extended sample .
precision near - infrared photometry is the critical ingredient for determining stellar masses for our sample ( brinchmann & ellis 2000 ) .
photometry was acquired with a @xmath24 filter with three different instruments : the keck near infrared camera ( nirc , matthews & soifer 1994 ) , the ukirt fast - track imager ( ufti , roche et al . 2003 ) , and the cooled infrared spectrograph and camera for ohs ( cisco , motohara et al .
2002 ) on the subaru 8.4 meter telescope ( see bundy et al . 2004 for a more detailed description of these data ) .
nirc has a field of view of 38  and a pixel scale of 0.15  pixel@xmath8 .
the equivalent numbers for ufti are : 96  field of view , with a pixel scale of 0.091 pixel@xmath8 .
the cisco camera has a field of view of 108with a pixel scale of 0.105  pixel@xmath8 .
the typical depths for these images is @xmath25 = 20.5 - 21 ( vega ) with a typical seeing of @xmath26 0.8 .    in each case , the infrared data were taken with a dither pattern whose step size exceeded the typical size of the galaxies of interest .
the data were reduced by creating sky and flat - field images from sets of disregistered neighboring science frames .
standard stars were observed for calibration purposes during the observations .
some images taken in good seeing but through thin cloud were subsequently calibrated in photometric conditions via shallower exposures taken with the wide field infrared camera ( wirc , wilson et al . 2003 ) on the hale 5 meter telescope .
we measure our photometry in the @xmath25-band , and hst i@xmath14 and v@xmath27 bands , within a scale - length factor , either 1.5r@xmath16 or 3r@xmath16 , both of which are generally large enough to avoid seeing effects from the @xmath24 imaging .
we then extrapolate the total magnitudes , within each band , out to infinite radius by using the fitted parameters for an exponential disk derived in the hst i@xmath14 band . to compare observables over a range in redshift ,
it is necessary to reduce all measures to a standard rest - frame .
galactic extinction corrections were applied using the formalism of schlegel , finkbeiner & davis ( 1998 ) and internal extinction was accounted for according to measured inclinations and luminosities using the precepts of tully et al .
it is debatable whether internal extinction corrections derived for nearby spirals are applicable to higher redshift disks .
our sample is mostly composed of systems with m@xmath28 where extinction could be a cause for concern .
however , direct extinction measurements in moderate redshift disks , determined through overlapping pairs , find a modest overall extinction ( white , keel & conselice 2000 ) .    to derive absolute magnitudes and rest - frame colors , we estimated non - evolutionary @xmath29-corrections .
the k - band @xmath29-corrections were computed from spectral energy distributions determined alongside the stellar mass fits .
we experimented with both the original algorithm developed by brinchmann & ellis ( 2000 ) and also an independent one developed by bundy , ellis & conselice ( 2005 , see @xmath133.2 ) .
both methods agree very well .
the approach is similar to that used by vogt et al .
( 1996 , 1997 ) where @xmath29-corrections were calculated from model seds from gronwall & koo ( 1995 ) based on various star formation histories in the bruzual & charlot ( 1993 ) code .
we also note that our derived k - correction values are very similar to non - evolutionary @xmath29-corrections derived by poggianti ( 1997 ) using seds from poggianti & barbaro ( 1996 ) .
the virial mass of a galaxy - that is , the combined total of the dark , stellar and gaseous components - is perhaps its most fundamental property .
however , obtaining an accurate estimate from observed quantities is necessarily difficult since there is no _ a priori _ agreed model for the relative distributions of the various components .
our approach to this challenge will be to use both analytical techniques and semi - analytical simulations to investigate the relationship between the total halo mass and our dynamical and structural observables . although necessarily approximate and debatable in terms of the assumptions made , we will attempt , where possible , to investigate the uncertainties involved by contrasting the two approaches in the context of our data .    for a simple virialized system , such as a circularly rotating disk , we can place constraints on the total mass within a given radius @xmath30 independent of model assumptions .
if @xmath31 , where we estimate the maximum rotational velocity v@xmath22 is reached , the mass within @xmath30 is given by ,    @xmath32    where @xmath31 .
this assumes that the rotation curve reaches the maximum velocity by 1.5 @xmath20 r@xmath16 , ( persic & salucci 1991 ) , otherwise v@xmath22 should be replaced by v@xmath33 .
there are a few possible approaches for determining total halo masses , some of which require the use of simulations to convert observed dynamical qualities , usually v@xmath22 , into halo masses .
we take a basic approach using eq .
( 1 ) to obtain the total mass with a given radius , and take a suitably large total radius of 100 kpc to measure the total halo mass .
this is often the extent of disk hi rotation curves , and similar to the sizes of dark matter halos ( e.g. , sofue & rubin 2001 ) .
another approach now being used ( e.g. , bohm et al . 2004 ) has been proposed by van den bosch ( 2002 ) . in this case , it is argued from analytic simulations of disk galaxy formation that the quantity @xmath34 gives , on average , the best empirical representation of the virial mass for simulated disks .
the zero - point of the relationship between r@xmath16v@xmath35/g and virial mass is claimed to be independent of feedback and independent of the mass of the halo ( van den bosch 2002 ) .
semi - analytic models based on @xmath36cdm ( e.g. , cole et al . 2000 ; benson et al .
2002 ; baugh et al .
2005 ) suggest , however , that the ratio between v@xmath37r@xmath16/g and halo mass is not as simple as the above formalism implies .
the latest _ galform _ models from baugh et al .
( 2005 ) and lacey ( priv .
communication ) show that equation [ 1 ] and van den bosch ( 2002 ) under - predicts the dark halo mass .
these models show that the relationship between the virial mass at r@xmath16 and the total halo mass changes as a function of mass in the sense that the ratio is higher for lower mass halos . physically this can be understood if high mass halos have a larger fraction of their baryonic mass in a hot gaseous phase that is not traced by the formed stellar mass .
using the _ galform _ models to convert our observables into a total halo mass is however potentially inaccurate as these models can not reproduce the tully - fisher relation to better than 30% .
however , there are reasons to believe that the semi - analytic models put the correct amount of halo mass into their modeled galaxies ( e.g. , benson et al .
this however does not necessarily imply that the v@xmath22 values in these models are able to accurately match the halo masses .
independent determinations of total halo masses are necessary to perform this test . however , the masses from this approach match within 0.5 dex the masses from using eq .
( 1 ) with a suitably large total radius .    with these caveats in mind , we have used the _ galform _ model results to fit the relationship between the virial mass at r@xmath16 ( eq .
1 ) , and the total mass of the halo ( m@xmath38 ) , a ratio which we call @xmath39 .
we fit @xmath40 as a linear function of @xmath41 , such that @xmath42 . using the _ galform _ results we fit @xmath43 and @xmath44 at redshifts @xmath45 .
we find that the functional form of @xmath40 does not change significantly with redshift , with typical values @xmath46 and @xmath47 .
the value of the halo mass m@xmath38 is then given by :    @xmath48    observational and model uncertainties contribute to errors on these virial and halo mass estimates in two ways .
measured scale lengths given by the gim2d and galfit fitting procedure (  2.1 ) , give an average error of 0.12 kpc , although we add an extra error to this to account for systematics seen when performing one component models to disks / bulge systems ( de jong 1996 ) . at the same time , systematic difficulties in the rotation curve analysis (  2.1 ) can arise .
vogt et al .
( 2005 ) discuss these issues in some detail and conclude the average error is @xmath49 kms@xmath8 ( typically @xmath5010 - 20% ) . following the discussion of rotation curve fitting in @xmath132.1 , it is possible we do not measure the true v@xmath22 required for insertion in equations ( 1 & 2 ) .
even a modest underestimate of v@xmath22 would lead to a significant error in @xmath51 . in combination , these measurement uncertainties imply virial and halo masses precise to no better than 30% .
the semi - analytic method for computing total halo masses is likely limited , and thus we also investigate the total halo masses found through the van den bosch ( 2002 ) formalism , and through the use of equation ( 1 ) out to 100 kpc .
we find that the total halo mass through eq .
( 1 ) , eq .
( 2 ) and van den bosch ( 2002 ) are all fairly similar
. we account for these differences in our use of total halo masses in  4.3 .
our procedure for deriving stellar masses follows the multi - color method introduced by brinchmann & ellis ( 2000 ) . combining hst v@xmath27 , i@xmath14 , and near - infrared @xmath24 photometry for a galaxy of known redshift
, we fit a range of template seds synthesized using software by bruzual & charlot ( 2003 ) .
our photometry is done in a matched aperture large enough to avoid seeing problems associated with the @xmath24-band .
these fitted seds constrain the @xmath24-band mass / light ratio .
our code for computing stellar masses only uses models with ages less than the age of the universe at the redshift observed . a @xmath52 analysis normalized by the near infrared @xmath24-band luminosity yields the stellar mass
. the template seds were constructed sampling a range of exponentially - declining star formation rates , metallicities and ages with a salpeter imf .
using other stellar mass functions , such as chabrier , kennicutt or kroupa would result in computed stellar masses that are smaller by @xmath53 0.3 dex .
we further assumed a simple exponentially declining star formation history with @xmath54 values ranging from 0.001 - 15.0 gyr , and metallicities from z=0.005 - 5 in solar units .
typical uncertainties in this method are a factor of three ( brinchmann & ellis 2000 ; papovich , dickinson & ferguson 2001 ; drory et al .
2004 ; bundy , ellis & conselice 2005 ) .
the sed is constrained by the observed colors as measured in an aperture of radius 1.5 @xmath20 r@xmath16 in each band which is optimal in terms of signal - to - noise ( s / n ) .
we assume that there are no color gradients and use the colors measured within this aperture as the global color .
the total @xmath24-band light is measured by extrapolating the 1.5 @xmath20 r@xmath16 flux to infinity , assuming the same exponential fit as measured in the hst i@xmath14 image .
errors on the stellar masses arising from photometric uncertainties can be determined in a monte carlo fashion .
simulated exponential disks of a known magnitude were inserted into the reduced hst and ground based images and photometrically recovered using the tools that were applied to the sample .
the simulated disks were arranged to randomly sample the selected ranges of disk scale - length @xmath18 and inclination .
the derived photometric errors were then input into the stellar mass calculations from which a @xmath55 range was calculated .
the average error in m@xmath56 is 0.47 dex .
the calculation of stellar masses through this technique is limited to a degree by systematics which are difficult to constrain with the current data set , with @xmath24-band magnitude errors only a small source of uncertainty .
the models of the spectral evolution of galaxies depend both on the observed stellar libraries and underlying theory , and some fundamental uncertainties remain ( c.f .
bruzual & charlot 2003 ) .
we expect uncertainties of about 5 - 10% , given the range of possible models ( charlot , worthey & bressan 1996 ) .
this is not important for the conclusions in this paper but might be an issue with larger and more accurate datasets in the future .
we also can not constrain the amount of recent star formation produced in bursts ( e.g. , kauffmann et al . 2003 ; bell & de jong 2001 ) , resulting in a slight systematic overestimate of the stellar mass .
the latter effect only becomes important at large burst fractions and should be @xmath57% for our big spirals ( c.f .
drory et al 2003 ) .
the final systematic uncertainty is the dust correction adopted .
we use the calzetti ( 1997 ) extinction law in our stellar mass calculations , although other extinction laws produce very small differences ( at most 10 - 20% ) in the resulting stellar mass ( papovich et al .
2001 ) . taken together
we are likely to have systematic uncertainties in our stellar mass estimates amounting to @xmath58 dex which is lower than our random uncertainty , thus they will not influence our results .
the tully - fisher ( tf ) relation has been the traditional method for investigating how dark halos and the stellar components of disk galaxies relate . the optical relation has been studied in disks galaxies out to redshifts @xmath4 by e.g. , vogt et al .
( 1996 , 1997 , 2005 ) , ziegler et al .
( 2002 ) and bohm et al .
these investigations have found between 0.4 and @xmath59 magnitudes of rest - frame b - band luminosity evolution in disks between @xmath60 and @xmath61 .
this luminosity evolution is derived by assuming that the slope of the tf relation at high redshift is the same as it is at @xmath60 .
the question of differential evolution in the relation ( ziegler et al . 2002 ; bohm et al . 2004 ) , for example that the faintest galaxies evolve more rapidly , remains an important unknown .
although the purpose of this paper is to move beyond the tf relation , we begin by plotting the k - band tf relation for our sample ( figure  1 ) .
one might expect evolution of the tf relation in the k - band to display a clearer signal than the b - band since the effects of dust are mitigated and passive evolution should be more uniform in its effect across the sample .
when we assume that the slope of the k - band tully - fisher is the same as the local value from verheijn ( 2001 ) we find no significant evolution , as is also found in the b - band ( vogt et al .
we perform these fits using both a downhill simplex amoeba and levenberg - marquardt @xmath62 minimization , both of which give the same results .
we find a fading of 0.04@xmath630.24 magnitudes for systems at @xmath64 compared with the @xmath60 relationship and a brightening of 0.37@xmath630.23 magnitudes for systems at @xmath65 , consistent with no evolution .
the scatter does not evolve significantly ( 1.13 magnitudes for the @xmath660.7 sample , and 0.72 magnitudes for systems at @xmath64 ) .
in all cases , the observed slope and scatter change only slightly when we ignore internal extinction corrections .
although our results are broadly consistent with earlier , smaller samples , the interpretation of any evolutionary signal is complicated in two ways .
first , only a limited range of luminosity and rotational velocity can be sampled at high redshift , leading to great uncertainties given the intrinsic scatter .
second , as luminosity and rotational velocity are indirect measures of the assembly state of the galaxy , both may be evolving in complex ways , that mask actual evolutionary changes .
the first step beyond the tf relation is to compare the stellar mass to the measured maximum velocity - a relation we will call the _ stellar mass - tully fisher relation_. the classical b - band tf relation scales such that @xmath67 .
this coupling becomes even steeper for the local stellar mass tully - fisher relation in nearby disks , @xmath68 ( bell & de jong 2001 ) .
ideally , we seek to measure the all - inclusive baryonic tf relation , but measuring the gas content of high redshift disks is not yet feasible .
we can estimate how much cold gas we are missing in our stellar mass inventory by investigating the gas mass fractions for nearby disk galaxies . through examinations of the luminosities and hi masses for nearby disks , mcgaugh & de blok ( 1997 )
conclude that galaxies which are massive , bright , red , or have a high - surface brightness , have very little gas in comparison to bluer , fainter , lower surface brightness systems .
systems which are brighter than m@xmath69 have gas mass fractions that are typically 0.1 or lower . since our selection finds the most luminous , high surface brightness systems , which are also red , they are the least likely sub - class of disks to have a high gas content .    the stellar mass tf relation is shown in figure  2 where , as before , we divide the sample into two redshift bins , split at @xmath70 .
each panel contains a solid line giving the @xmath60 best fit and a dashed line illustrating the @xmath63 3 @xmath71 uncertainty in this fit ( bell & de jong 2001 ) . as was the case for the conventional tf relations , no significant evolution in the zero - point
is observed .
the bell & de jong ( 2001 ) @xmath72 stellar mass tully - fisher relation can be written as m@xmath73 . by holding the slope of this relationship constant , we find that the zero point is best fit by 0.45@xmath630.12 at @xmath74 and 0.41@xmath630.13 at @xmath64 .
neither of these are however significantly different from the @xmath60 relationship , and are very similar to each other .
this implies that , if growth continues , the stellar and dark components are growing together .
for example , if disk assembly since @xmath75 proceeded only by the addition of stellar mass at a uniform rate of 4 m@xmath76 year@xmath8 , the local zero point would be discrepant at the 4@xmath71 level .
this lack of evolution is important for understanding how disk galaxy formation is occurring ( see  4.4 ) .
moreover , the scatter in the stellar mass tully - fisher relation is similar to that observed in the k - band tully - fisher after converting the k - band magnitude scatter into a luminosity and assuming an average stellar mass to light ratio .
the typical scatter ( in log m@xmath56 units ) in stellar mass for these is 0.65 for disks at @xmath74 and 0.48 for those at @xmath64 .
the final step in our analysis will be an attempt to convert our measured quantities into a comparison of the stellar and halo masses as discussed in @xmath133 .
recognizing the considerable uncertainties involved , figure 3 shows as a function of redshift the ratio , @xmath77 m@xmath56/m@xmath78 , between the stellar masses and virial masses ( eq .
1 ) within 3r@xmath16 for our sample . as expected , the stellar masses are nearly always less than the independently - derived virial masses , indicating that our methods for computing these values are not dominated by large systematic errors .
it can also be seen from this figure that there is wide range of m@xmath56/m@xmath78 values at every redshift from @xmath79 .
the open pentagons show the median values of @xmath80 as a function of redshift .
the dotted horizontal line shows the global baryonic mass fraction @xmath81 = 0.171 derived from wmap results ( spergel et al .
within 3r@xmath16 it appears that many disks have m@xmath56/m@xmath78 values higher than this limit .
it also appears that within 3r@xmath16 the stellar mass of the disk dominates the virial mass .
this indicates that within the visual parts of some of our sample , the disk component accounts for roughly all the mass , which is consistent with a maximal disk interpretation . however , there are clear examples on figure  3 where either the stellar mass is overestimated and/or the virial mass is underestimated .
both of these are possibilities , since our stellar masses are potentially too high from using a salpeter imf , and we might be underestimating the value of v@xmath22 due to seeing , or a lack of depth in the lris observations . there is also a slight bias such that at higher redshifts nearly all galaxies sampled have high v@xmath22 values , which are likely maximal , while at lower redshifts we are sampling systems with lower v@xmath22 values , that have not yet fully formed their stellar masses .
this is likely part of the reason that the stellar mass to virial mass ratio increases slight at higher redshifts .
however , we are also interested in constraining the relationship between the total halo mass , m@xmath38 , and the total stellar mass m@xmath56 . as we discuss in
 3 , it is very difficult to accurately obtain halo masses .
we therefore show in figure  4 total halo masses derived through eq .
( 1 ) at 100 kpc and through the semi - analytical approach ( eq .
although we can not accurately determine total halo masses for individual systems , our main goal is to compare how the ratio of stellar to halo mass changes with redshift .
figure  4 shows this relationship divided into the same redshift bins as in figures 1 & 2 .
to first order , the halo masses and stellar masses of disk galaxies should correlate if star formation is regulated in the same manner in halos of different masses ( steinmetz & navarro 1999 ) .
although there is no reason to expect any particular functional form between these two quantities , there is a reasonably well - fit linear relationship between them .
we find that the zero point and slope of this relation does not change significantly between low and high redshift .
there is also no obvious change in the scatter from high to low redshift ( 0.26 cf .
the most significant outcome of figures 3 and 4 however is the remarkable similarity in trends found at high and low redshift suggesting some disks had completed the bulk of their stellar assembly by @xmath821 or more likely that the stellar and dark masses of galaxies grow together continuously .
figure  5 shows the distribution of @xmath80 = m@xmath56/m@xmath38 for systems more massive than , and less massive than , the average halo mass , m@xmath83 .
the solid line is the universal baryonic mass ratio .
as mentioned earlier , most disk stellar mass fractions are lower than the cosmic ratio , which seems appropriate given our stellar mass inventory is not intended to account for all associated baryons .
although some dispersion is expected , conceivably some fractions are overestimated or underestimated , for example by making incorrect assumptions about the imf , or underestimating v@xmath22 ( and hence m@xmath38 ) by insufficient sampling of the rotation curve .
figure  5 also shows a tentative population of disk galaxies with remarkably low stellar fractions , the most extreme cases occurring in objects with halo masses m@xmath84 .
these galaxies deviate from the @xmath72 stellar mass tf relation by more than 4 @xmath71 .
investigation of the individual systems that lie in this category shows them to be undergoing vigorous star formation as inferred by bluer than average @xmath85 colors .
a weak correlation was found between @xmath80 and @xmath85 . because these systems are blue , they may have kinematic asymmetries in their rotation curves which may raise their v@xmath22 values . however , this is unlikely the effect producing this slight correlation as most of the low @xmath80 systems also have low v@xmath22 values .
we investigated several other possible correlations involving @xmath80 ( for example with the bulge / disk ratio ) but no significant trends were found .
to investigate the implications of our results for the assembly history of stellar mass in disk galaxies since @xmath4 , we return to the behavior of the stellar fraction @xmath80 = m@xmath56(@xmath533r@xmath16)/m@xmath78(@xmath533r@xmath16 ) vs. redshift ( figure  3 ) , and the relationship between total stellar mass ( m@xmath56 ) and halo mass ( m@xmath38 ) at @xmath64 and @xmath74 ( figure 4 )
. trends in these relations may shed light on the physical processes regulating star formation in disk galaxies .
this is only the case however if the value of r@xmath16 for disks does not grow with time , which appears to be the case for the largest disks , based on our own limited data , and in statistical studies of @xmath86 disks ( ravindranath et al .
2004 ; conselice et al .
2004 ; cf .
ferguson et al .
2004 for higher redshifts ) .
first , figure  3 shows that the median stellar vs. virial mass value for our sample does not change significantly with redshift .
we can compare two simple extreme models to this result to determine the likely method by which disk galaxies are forming at @xmath3 . in a ` monolithic collapse ' model , where the dark halo and baryons for a galaxy are in place at high redshift , and there is no change in these mass components with time ( e.g. , eggen , lynden - bell & sandage 1962 ) ,
the average value of m@xmath56/m@xmath78 should increase with time .
this assumes that effective radii distribution for disks does not change with time , which is observed ( ravindranath et al .
there is evolution in the monolithic model only in the sense that baryons are gradually converted to stars over time .
this is the opposite to the observed trend in which the ratio of m@xmath56/m@xmath78 is roughly constant .
on the other hand , in a hierarchical picture , the value of m@xmath56/m@xmath78 remains relatively constant with time .
shown in figure  3 are model predictions from _ galform _ which show that the value of m@xmath56/m@xmath78 remains relatively constant with redshift .
while we do not have a complete sample , we do find that the average and median values of m@xmath56/m@xmath78 for our sample remain constant with redshift within the errors , which is consistent with the hierarchical idea .
the higher average m@xmath56/m@xmath78 value at @xmath86 is likely produced by our selection of the brightest disks at these redshifts ( vogt et al .
the _ galform _ model however does not predict the full range of mass fractions that we find . at all redshifts
there are systems dominated within r@xmath16 by their stellar masses , while some systems have low m@xmath56/m@xmath78 ratios .
this is either produced by an observational bias (  4.3 ) , or the simulations predict too much dark matter in the centers of disks .
figure  4 shows the relationship between total stellar and total halo masses ( see  3.1 ) for our sample , divided into our two redshift ranges .
we also show on figure  4 the _ galform _ semi - analytic model values for this relationship for disk dominated galaxies with the 80% completeness indicated . as in figure  3 ,
no strong evolution in m@xmath56/m@xmath38 ratios is predicted .
as galaxies in the semi - analytic models grow by accreting smaller systems , or intergalactic gas coupled with dark matter , and then converting the newly - obtained gas into stars quickly , the relationship between m@xmath56 and m@xmath38 remains constant .
this is generally what we find , as well as a good agreement with the model predictions .
in fact , the agreement between the total stellar and halo masses with the _ galform _ models is slightly better than the comparison between the stellar and virial masses with 3@xmath20r@xmath16 , which further suggests that the luminous components of disk galaxies are dominated by the stellar mass (  4.3 ) .    in summary , the fact that the ratio m@xmath56/m@xmath38 remains relatively constant in our sample with redshift , and that the stellar - mass tf relation does not evolution , are indications that disk galaxies are forming through the accretion of both dark and baryonic mass .
disk galaxies are undergoing star formation at @xmath3 at a rate of a few solar masses per year .
if disks at @xmath4 contained all the baryonic or halo masses that they have at @xmath60 , we would see an increase in m@xmath56/m@xmath38 with time . because we do not see this trend
, it appears that new stars form out of gas accreted from the intergalactic medium , which is coupled with dark matter at a constant ratio .
this baryonic to dark matter ratio is such that both m@xmath56 and m@xmath38 grow with time , creating a stellar mass to halo mass ratio that is on average relatively constant .
we present the results of a dynamical and structure study of 101 disk galaxies drawn mostly from the deep1 survey with redshifts in the range @xmath820.2 - 1.2 .
new infrared imaging observations are presented which enable us to derive reasonably reliable stellar masses and thereby to construct the stellar mass tully - fisher relation and its redshift dependence . using various formalisms drawn from analytic and semi - analytical models , we attempt to convert our dynamical data to make the first comparisons of the relative fractions of stellar and total mass in our sample . notwithstanding the considerable uncertainties and sample incompleteness
, the results are encouraging and suggest remarkably little evolution in the mix of baryons and dark matter since @xmath821 .
although our sample is not formally complete in luminosity or mass , we explore the degree to which there may have been evolution in the relative distribution of stellar and virial masses and find the following :    1 .
massive disk galaxies exist out to @xmath4 with halo masses as large as 10@xmath87 , roughly as large as the most massive disks in the nearby universe .
these systems also contain a large amount of stellar mass .
at least some disk galaxies are nearly fully mature in their stellar content at @xmath88 .
we confirm earlier studies based on smaller samples and find no significant evolution in the zero point or scatter of the rest - frame @xmath89-band tully - fisher relation out to @xmath12 .
the stellar mass tully - fisher relation out to @xmath12 is likewise largely consistent with the relation found for nearby disks .
we find no significant evolution in our sample after comparing systems at redshifts greater than and less than @xmath90 .
4 .   although there are clearly great uncertainties in estimating total halo masses from our dynamical data , we find the distribution of the ratio of stellar and halo masses remains relatively similar from @xmath60 to @xmath12 .
the stellar fraction observed can be understood if the bulk of the baryons associated with massive disk galaxies have already formed their stars .
a modest number of massive galaxies have very low stellar fractions , consistent with continued star formation as revealed by their blue @xmath91 colors .
these results are in relatively good agreement with @xmath36cdm analytical and semi - analytical models ( benson et al .
2002 ; baugh et al .
2005 ) , suggesting that disk galaxy formation is hierarchical in nature .
our overarching conclusion from this study is that no significant evolution in the stellar mass fraction can be detected in the population of regular massive disks since @xmath821 .
although biases and uncertain assumptions may affect detailed quantities at the @xmath92 dex level , the absence of gross trends is consistent with the conclusion that the bulk of these systems grow at @xmath3 by the accretion of dark and baryonic material .
this conclusion is however tempered by the fact that our selection is affected by a bias that likely changes with redshift .
for example , because our sample was selected by a @xmath93 limit , we are studying the brightest disks at the highest redshifts .    these results , while intriguing , will be examined in much greater detail with the upcoming deep2 redshift survey ( davis et al .
2002 ) , which will combine near - infrared imaging from the hale 5 meter telescope with rotation curves obtained using the deimos spectrograph on the keck ii telescope .
there are several issues that can be better addressed with deeper spectroscopy and a more careful selection of high - z disks . the total number of disks in this much larger sample is expected to be on the order of thousands and will solidify and expand on the preliminary results revealed by this study .
we thank the deep team for generously enabling us to augment their optical sample and catalogs with near - infrared data in order to make this comparison .
we also thank cedric lacey and andrew benson for access to their semi - analytical model simulation results , and xavier hernandez and ken freeman for comments regarding this work .
cjc acknowledges support from a national science foundation astronomy and astrophysics postdoctoral fellowship .
npv is pleased to acknowledge support from nsf grants nsf-0349155 from the career awards program , nsf-0123690 via the advance it program at nmsu , and ast 95 - 29098 and 00 - 71198 administered at ucsc , and nasa stsci grants go-07883.01 - 96a , ar-05801.01 , ar-06402.01 , and ar-07532.01 .
abadi , m.g . ,
navarro , j.f . , steinmetz , m. , & eke , v.r .
2003 , apj , 591 , 499 baugh , c.m . ,
lacey , c.g . ,
frenk , c.s . ,
granato , g.l . ,
silva , l. , bressan , a. , benson , a.j .
, & cole , s. 2005 , mnras , 356 , 1191 bell , e.f .
, & de jong , r.s . 2001 , apj , 550 , 212 benson , a.j . , cole , s. , frenk , c.s . ,
baugh , c.m . , & lacey , c.g .
2000 , mnras , 311 , 793 benson , a.j . , lacey , c.g . , baugh , c.m .
, cole , s. , & frenk , c.s .
2002 , mnras , 333 , 156 bohm , a. , et al .
2004 , a&a , 420 , 97 brinchmann , j. , & ellis , r.s .
2000 , apj , 536 , l77 brinchmann , j. , et al .
1998 , apj , 499 , 112 brinchmann , j. 1999 , phd , cambridge university bruzual , a.g . , &
charlot , s. 2003 , mnras , 344 , 1000 bundy , k. , fukugita , m. , ellis , r.s . , kodama , t. , & conselice , c.j .
2004 , apj , 601 , l123 bundy , k. , ellis , r.s . , & conselice , c.j .
2005 , apj , astro - ph/0502204 calzetti , d. 1997 , aj , 113 , 162 cole , s. , lacey , c.g . , baugh , c.m . , & frenk , c.s .
2000 , mnras , 319 , 168 charlot , s. , worthey , g. , & bressan , a. 1996 , apj , 457 , 625 conselice , c.j .
2004 , apj , 600 , l139 davis , m. , et al .
2002 , spie , 4834 , 161 de jong , r.s .
1996 , a&a , 118 , 557 dickinson , m. , papovich , c. , ferguson , h. c. , & budvari , t. 2003 , apj , 587 , 25 drory , n. , bender , r. , feulner , g. , hopp , u. , maraston , c. , snigula , j. , & hill , g.j .
2004 , apj , 608 , 742 eggen , o.j . , lynden - bell , d. , & sandage , a.r .
1962 , apj , 136 , 748 gronwall , c. , & koo , d.c .
1995 , apj , 440 , l1 groth , e.j . ,
kristian , j.a . ,
lynds , r. , oneil , e.j . , balsano , r. , rhodes , j. 1994 , bass , 185 , 5309 ferguson , h.c . ,
2004 , apj , 600 , l107 haynes , m.p . , et al .
1999 , aj , 117 , 2039 kauffmann , g. , et al .
2003 , mnras , 341 , 33 kinney , a. , calzetti , d. , bohlin , r. , mcquade , k. , storchi - bergmann , t. , & schmitt , h.r .
1996 , apj , 467 , 38 madau , p. , pozzetti , l. , & dickinson , m. 1998 , apj , 498 , 106 mcgaugh , s.s .
, & de blok , w.j.g .
1997 , apj , 481 , 689 matthews , k. , & soifer , b.t .
1994 , in `` infrared astronomy with arrays , the next generation '' , editied by ian s. mclean , 190 , 239 motohara , k. , et al .
2002 , pasj , 54 , 315 oke , j.b .
1995 , pasp , 107 , 375 papovich , c. , dickinson , m. , & ferguson , h.c .
2001 , apj , 559 , 620 peng , c.y .
, ho , l.c . , impey , c.d . , & rix , h - w .
2002 , aj , 124 , 266 persic , m. , & salucci , p. 1991
, apj , 368 , 60 poggianti , b.m .
1997 , a&as , 122 , 399 poggianti , b.m . , & barbaro , g. , 1996 , a&a , 314 , 379 ravindranath , s. , et al .
2004 , apj , 604 , l9 roche , p.f .
2003 , spie , 4841 , 901 schlegel , d.j . , finkbeiner , d. , & davis , m. 1998 , apj , 500 , 525 simard , l. et al .
2002 , apjs , 142 , 1 sofue , y. , & rubin , v. 2001 , ara&a , 39 , 137 spergel , d.n .
2003 , apjs , 148 , 175 steinmetz , m. , & navarro , j.f .
1999 , apj , 513 , 555 tully , r.b . , & fisher , j.r . 1977 , a&a , 54 , 661 tully , r.b . ,
pierce , m.j . , huang , j - s . , saunders , w. , verheijen , m.a.w . , &
witchalls , p.l .
1998 , aj , 115 , 2264 van den bosch , f.c .
2002 , mnras , 332 , 456 verheijen , m.a.w .
2001 , apj , 563 , 694 vogt , n.p . ,
forbes , d.a . ,
phillips , a.c .
, gronwall , c. , faber , s.m . ,
illingworth , g.d .
, & koo , d.c .
1996 , apj , 465 , l15 vogt , n.p . ,
1997 , apj , 479 , l121 vogt , n.p .
2004 , apj , submitted vogt , n.p .
, et al . 2005 , in preparation weiner , b.j . , et al .
2004 , apj , in press williams , r. , et al .
1996 , aj , 112 , 1335 white , r.e . ,
keel , w.c . ,
conselice , c.j .
2000 , apj , 542 , 761 wilson , j.c . ,
2003 , spie , 4841 , 451 ziegler , b.l .
2002 , apj , 564 , l69


ultracold atomic gases offer various opportunities for the study of interacting many - body quantum systems in a well - controlled environment @xcite . for instance , the bose - hubbard model  a simplified description for bosons in a periodic potential  can be studied with bose - einstein condensates confined in optical lattices @xcite .
quantum gases have also gained much attention lately regarding the emergence of an effective space - time @xcite : their low - energy phase fluctuations obey the same covariant field equations as a scalar quantum field in a certain curved space - time .
hence , the study of phonons in this laboratory system might shed some light on aspects of cosmic quantum effects , e.g. , hawking radiation @xcite or the freezing and amplification of quantum fluctuations in expanding spacetimes @xcite .
although the fluctuations in bose - einstein condensates are usually small , it has recently become possible in experiments to go beyond the classical order parameter and resolve signatures of the fluctuations @xcite .    theoretically , the fluctuations in a bose - einstein condensate are usually treated as small perturbations of the mean field .
the solution of the coupled field equations is rather demanding and often requires further approximations , especially for time - dependent condensates .
the hartree - fock - bogoliubov method , see , e.g. , @xcite , permits in principle the self - consistent propagation of the mean field and the quantum correlations for arbitrary time - dependences of the trap potential or interaction strength .
but the scaling of the numerics with system size often limits the actual calculations to a low number of dimensions , certain symmetries , or a short time interval .
thermal condensates might be studied using the projected gross - pitaevskii equation @xcite , where the low - energy part of the fluctuations is expanded into an arbitrary basis and , in view of the large thermal occupation , treated classically ; higher excitations as well as the vacuum effect are omitted . on the other hand , studies in the context of expanding spacetimes , e.g. , @xcite ,
indeed focus on the quantum fluctuations but often assume a homogeneous background or start with the hydrodynamic action , which is only valid on scales longer than the healing length .    in this article
, i will discuss the evolution of the quantum fluctuations in trapped time - dependent bose - einstein condensates .
the linear field operators will be expanded into their eigenmodes thus permitting the diagonalization of the ( initial ) evolution equations .
although basis expansions of the field operator are frequently used ( e.g. , in @xcite ) , these references usually consider the harmonic oscillator eigenfunctions , a large number of which must be used in order to describe the excitations properly . by adopting the eigenmodes ,
a much smaller part of the basis needs to be considered and many more situations will become numerically feasible .
( note , however , that in order to obtain the fluctuation eigenmodes , a relatively large number of oscillator functions must be employed  but they need not be propagated . )    this article is organized as follows .
section [ secfieldequations ] reviews the field equations and their linearization .
the bogoliubov - de gennes equation for the linear quantum fluctuations can be diagonalized by an eigenmode expansion , which will be performed in sec .
[ seceigenmode ] .
however , as soon as trap frequency or interaction strength are varied , off - diagonal terms appear .
this coupling of different modes is negligible for excitations with energies much smaller than the chemical potential even in time - dependent condensates , as will be shown in sec .
[ sectf ] , where the order parameter is treated in the thomas - fermi approximation .
it is then also possible to establish the analogy between phase fluctuations and a massless scalar field in a certain curved spacetime .
if the excitation energies are of the same order as the chemical potential , the coupling of different modes might lead to a population transfer , which will be illustrated in an example in sec .
[ sec1d ] .
in dimensionless units , the field operator @xmath0 of a trapped ( quasi)-@xmath1-dimensional bose - einstein condensate obeys the non - linear schrdinger field equation @xcite @xmath2          \hat \psi \ , ,      \label{fieldequation }      \end{aligned}\ ] ] with @xmath1-dimensional coupling strength @xmath3 . by changing the @xmath4-wave scattering length @xmath5 through feshbach resonances @xcite or by varying the trap frequency @xmath6 ,
an external time - dependence can be prescribed on the condensate .
the gas cloud will adapt to these changes and it will either expand or contract and with it the quasi - particle excitations residing upon it .
a part of this background motion can be accounted for by transforming to new coordinates @xmath7 with scale factor @xmath8 .
the field operator then reads @xcite @xmath9 the phase @xmath10 is chosen such as to generate an isotropic velocity field @xmath11 which ( at least partially ) describes the expansion / contraction of the condensate .
if the scale factor @xmath8 obeys @xmath12 with @xmath13 being the initial value of the coupling strength , a scaled field equation follows @xmath14 \hat \psi \ , ,      \label{scaledheisenberg }      \end{aligned}\ ] ] where trapping and interaction terms have acquired the same time - dependent pre - factor @xmath15 and all other coefficients are time - independent .
( the scale factor @xmath8 on the left hand side might be included into a redefined time @xmath16 , see sec .
[ secspacetime ] . )      for large particle numbers @xmath17 , one might formally expand the field operator into inverse powers of @xmath18 @xcite @xmath19 here , @xmath20 and @xmath21 are the atomic operators .
they commute with the linear @xmath22 and higher - order quantum excitations @xmath23 and thus yield the exact conservation of particle number .
the order parameter @xmath24 in the center of the trap but diminishes towards the edge of the condensate .
insertion of the expansion into the scaled heisenberg equation yields the gross - pitaevskii equation for the classical background @xmath25 @xcite @xmath26 \psi_0 \ , .
\label{gp }      \end{aligned}\ ] ] the linear quantum fluctuations @xmath27 obey the bogoliubov - de gennes equation @xcite @xmath28 \hat \chi          + f^2 g_0 \psi_0 ^ 2 \hat\chi^\dagger \
, ,      \label{bdg }      \end{aligned}\ ] ] and the residual terms comprise the equation of motion for @xmath23 @xmath29 \hat \zeta          + g_0 \psi_0 ^ 2 \hat \zeta^\dagger      \nonumber\\      & \quad + g_0 \left ( 2 \psi_0 \hat \chi^\dagger \hat \chi              + \psi_0^ * \hat \chi^2              + \hat \chi^\dagger
\hat \chi^2          \right ) + { { \cal o } } ( g_0 \hat \zeta ) \ , .
\label{zeta }      \end{aligned}\ ] ] these higher orders @xmath23 must remain small in order for the mean - field expansion to be valid , i.e. , for the linearized equation to be applicable .
this means that the terms involving products of @xmath27 must remain small because they act as source terms for higher orders @xmath23 .    from the gross - pitaevskii equation ,
i can infer when the evolution of the order parameter @xmath25 is solely described by the scale factor @xmath8 : apart from the trivial case @xmath30 , this occurs only when the spatial derivatives can be neglected with respect to the interaction and trapping terms , @xmath31 , i.e. , in the thomas - fermi approximation .
density and phase of the order parameter @xmath32 then assume the form @xmath33 where the heaviside step function @xmath34 is @xmath35 for @xmath36 and @xmath37 elsewhere . in this approximation , the motion of the classical background becomes stationary and the scaled coordinates @xmath38 are co - moving with the condensate .
the bogoliubov - de gennes equation can be tackled by introducing self - adjoint operators @xmath39 with @xmath40 being the phase of the order parameter .
these operators resemble ( relative ) density and phase fluctuations @xmath41 and @xmath42 up to the prefactor @xmath43 .
since this prefactor eventually becomes large near ( and beyond ) the surface of the condensate , the smallness of @xmath44 and @xmath45 can not be ensured .
therefore , i will stick to @xmath46 in the following , but still refer to them as density and phase fluctuations .
they obey @xmath47 \hat \chi_-          \ , = \ , - \mathcal k_+ \hat \chi_+ \ , ,      \nonumber\\      & \frac{1}{2 } \left [ b^2 \frac{\partial}{\partial t }          + \bm v_0 \bm \nabla_x          + \frac{1}{2 } ( \bm \nabla_x \bm v_0 )      \right ] \hat \chi_+
\ , = \ , \mathcal k_- \hat \chi_- \ , ,      \label{eomchipm }      \end{aligned}\ ] ] where the velocity field @xmath48 results from the residual background phase beyond the thomas - fermi approximation .
however , @xmath49 is small for @xmath50 and if @xmath51 does not change too swiftly , because an almost homogeneous phase will develop with only small deviations near the boundary of the condensate . whereas attractive @xmath52 invalidate the thomas - fermi approximation and generally @xmath53 even in the center of the trap .
the differential operators on the right - hand sides @xmath54 generally do not commute for inhomogeneous condensates @xmath55 @xmath56          \ , & = \
,          f g_0 \left\ { ( { \mbox{\boldmath$\nabla$}}_x^2 \varrho_0 )              + 2 ( { \mbox{\boldmath$\nabla$}}_x\varrho_0 ) { \mbox{\boldmath$\nabla$}}_x              \right\ } \ , .      \label{commkpm }      \end{aligned}\ ] ]
in order to define the initial state unambiguously , i will assume that the condensate is at rest before @xmath57 .
then @xmath58 , @xmath59 , @xmath60 , and @xmath61 such that the left hand sides of eqs .
reduce to partial time derivatives and the initial eigenmode equations for @xmath62 follow @xmath63 because @xmath64 for inhomogeneous condensates , cf .
eq .  , density and phase fluctuations of each mode must have different space dependences .
this leads to the expansions @xcite ( i will adopt the sum convention throughout this article for brevity ; any indices appearing only on one side of the equation are to be summed ) @xmath65 of @xmath62 into different eigenmode bases @xmath66 and @xmath67 , see appendix [ appeigenmodes ] for more details on how to obtain @xmath68 .
usually , these two bases are neither orthogonal nor normalized , , but instead can be chosen to be dual to each other @xmath69 note that this condition does not fix the norm of @xmath68 but still permits the multiplication by an arbitrary factor , @xmath70 and @xmath71 .
observables must be unaffected by this ambiguity , see app.[secobservables ] .
insertion of the eigenmode expansion into the linear field equations yields a set of coupled first - order differential equations @xmath72 with time - dependent coefficients .
the symmetric matrices @xmath73 are initially diagonal @xmath74 and @xmath75 . at later times
, they acquire off - diagonal elements because of the different - time commutators @xmath76 \neq 0 $ ] and @xmath77 \neq 0 $ ] when @xmath78 .
the velocity coupling matrix @xmath79 h_m^-      \label{velo }      \end{aligned}\ ] ] is not symmetric , but vanishes for homogeneous phases @xmath80 of the order parameter , e.g. , initially or in the thomas - fermi approximation . for repulsive
@xmath50 and slow variations of interaction strength @xmath81 and trap frequency @xmath6 , the order parameter phase is homogeneous except for small ripples near the boundary of the condensate such that the matrix @xmath82 is usually negligible .    from the evolution equations with the initially diagonal coupling matrices ,
the introduction of bosonic operators @xmath83 and @xmath84 creating or annihilating an initial quasi - particle is straightforward @xmath85 here , a bar shall denote complex conjugation , e.g. , @xmath86 .
the coefficients obey the initial values @xmath87 where the phase has been appropriately chosen and the frequencies @xmath88 .
the upper index of the coefficients @xmath89 and @xmath90 labels the mode , while the lower index denotes the components of this particular mode when expanded in a certain basis , e.g. , the initial eigenfunctions @xmath91 .
since the coupling matrices become non - diagonal even for slow ( adiabatically ) variation of the trap frequency @xmath6 or coupling strength @xmath81 , the initial bases @xmath91 can not represent the eigenmodes at later times .
although @xmath68 might be employed in order to calculate the spatial correlation functions , see appendix [ secobservables ] , the use of these functions might be misleading regarding the correlations between different modes . furthermore , when probing the excitations using , e.g. , the scheme proposed in @xcite , the proper particles defined at the time of measurement will be detected and not the initial ones .
the particle definition in time - dependent background is a non - trivial task , see , e.g. , @xcite .
nonetheless , it is always possible to expand @xmath62 @xmath92 into bases @xmath93 , which are defined such that the coupling matrices and , cf .  , become diagonal at any particular instant @xmath94 @xmath95 of course , the velocity term @xmath96 is then generally non - diagonal and the evolution equations of the different modes will not exactly decouple at this particular instant @xmath94 .
but one should bear in mind that measurement occurs usually in an adiabatic region , where the external parameters are only slowly - varying functions of time .
then , the background phase is approximately homogeneous and the velocity @xmath97 .
hence , quasi - particle creators and annihilators might be introduced analogous to eq .
@xmath98 where @xmath99 .
in the previous sections , i made no approximations except for the linearization and the assumption of an isotropic trap .
the formalism is , in principle , applicable for arbitrary variations of trap frequency @xmath6 and interactions @xmath81 . to this end
, it would be necessary to solve the gross - pitaevskii equation and the linear evolution equations simultaneously .
the numerical solution is complicated by the fact that the coupling matrices and need to be calculated at each time step .
some of the numerical difficulties can be circumvent by adopting the thomas - fermi profile , where density and phase of the background become time - independent ( in the coordinates @xmath38 ) and thus require the calculation of the coupling matrices only once . despite some shortcomings regarding the dynamics of the order parameter , this approximation is usually applicable for repulsive interactions and in the center of the trap , but becomes inaccurate towards the surface of the condensate , where the quantum pressure @xmath100 is relevant .      within the thomas - fermi approximation , the coordinate transformation @xmath101 associated with the scaling transformation
renders the background density time - independent @xmath102 , while the phase becomes homogeneous @xmath103 and thus @xmath104 .
the integrals of the coupling matrices simplify considerably and the evolution equations can be cast into the form @xmath105 i.e. , @xmath106 and @xmath107 can be split into time - independent diagonal parts , cf .
@xmath108 and constant , symmetric , non - diagonal coupling matrices @xmath109 with time - dependent prefactors @xmath110 .
the external variation of trap frequency @xmath6 and coupling strength @xmath81 is solely encoded in the scale factor @xmath8 and the scalar function @xmath111 . note also that the coefficients and thus also the evolution equations for the fluctuations are independent of the initial coupling strength @xmath13 .
the addend @xmath112 appearing in the parentheses of eqs . and can be expressed through the chemical potential @xmath113 , where @xmath114 , cf .
the thomas - fermi equation .
the mode functions @xmath68 of excitations with low energies , @xmath115 , are localized inside the condensate . hence , changing the bounds of the integrals and from infinity to the thomas - fermi radius @xmath116 will not alter these matrix elements significantly . bearing further in mind that for @xmath117 , it follows @xmath118 and .
also because of the restriction to low energies , , and one gets @xmath119 .
hence , the evolution equations of different eigenmodes approximately decouple and i obtain second - order equations of motion for phase @xcite @xmath120          \delta\hat \phi_n \ , = \ , 0 \ , ,      \label{effkfgmode }      \end{aligned}\ ] ] and density fluctuations @xmath121          \delta\hat\varrho_n \ , = \ , 0 \ , ,      \label{effkfgdensitymode }      \end{aligned}\ ] ] where i also introduced proper time @xmath16 .
equation is the evolution equation of a mode of a minimally - coupled massless scalar field in a friedman - lematre - robertson - walker spacetime @xcite , provided the scale factor @xmath122 of the space - time is identified with @xmath123 , cf .
@xcite @xmath124 [ note that the prefactor of the damping term is @xmath125 in eq .  , while it is @xmath1 in a @xmath126-dimensional friedman - lematre - robertson - walker spacetime .
]    the analogy is not restricted to the evolution equations in mode expansion but applies in the low - energy limit of the field equations as well : within the thomas - fermi approximation , @xmath127 for @xmath117 and @xmath61 , and for low excitations , @xmath128 , the phase fluctuations obey a second - order field equation @xcite @xmath129 \delta\hat\phi \ , = \ , 0 \ , ,      \label{effkfg }      \end{aligned}\ ] ] which is similar to that of a minimally coupled scalar field in a friedman - lematre - robertson - walker spacetime .
having established this kinematical analogy , cf .
, some of the concepts of general relativity can be applied to time - dependent bose - einstein condensates .
sonic analogs of horizons @xcite are of particular interest for the study of non - equilibrium effects , because they give a rough estimate whether and when adiabaticity will be violated and the ( quantum ) fluctuations freeze and get amplified , i.e. , ( quasi-)particle production occurs . an effective particle horizon occurs , if a phonon emitted at a time @xmath130 can only travel a finite ( co - moving ) distance , i.e. , if the integral @xmath131 converges to a finite value @xmath132 for @xmath133 .
wavepackets emitted at time @xmath130 at the origin @xmath134 can reach only points within the horizon , @xmath135 , in a finite time .
all other points are concealed by the horizon .
[ for simplicity , i assumed in eq .   an homogeneous sound velocity @xmath136 . ]      in order to point out the analogy of phase fluctuations to cosmic quantum fields
, i formulated the evolution equations and using proper time @xmath137 . on the other hand ,
experiments are usually performed in the laboratory and thus the variations of trap frequency @xmath138 and coupling strength @xmath139 are prescribed in laboratory time @xmath140 . since @xmath137 is a complicated function of @xmath140 , it is not quite obvious whether or not the quantum fluctuations will experience non - adiabatic evolution for a given modulation of @xmath6 or @xmath81 . in laboratory time
@xmath140 , eq
.   reads @xmath141          \delta\hat\varrho_n \ , = \ , 0 \ , ,      \label{effkfgdensitymodet }      \end{aligned}\ ] ] which is the evolution equation of a damped harmonic oscillator with time - dependent coefficients @xmath142 and @xmath143 initially , when @xmath59 , the field modes perform free oscillations . upon the gradual increase of the damping term @xmath144 with respect to the oscillation frequencies @xmath145 , the non - adiabatic evolution of the quantum fluctuations slowly sets in , until they finally freeze and get amplified when both terms , @xmath142 and @xmath145 are of the same order @xcite .
let me discuss the two extremal ways a time - dependent scale factor @xmath8 can be achieved : firstly , only the trap frequency might be varied , while the interaction strength @xmath146 . with instantaneous frequencies @xmath147 ,
adiabaticity can be violated for any finite change @xmath148 if @xmath149 becomes sufficiently large .
the quantum fluctuations can not adapt to the changing background any more , they freeze and get amplified . and , secondly , for static traps , @xmath150 , where @xmath81 is time - dependent .
then , the situations is not so clear because @xmath151 , cf .
hence , only the rapid acceleration of the scale factor @xmath152 will lead to notable changes of the excitation frequencies . on the other hand ,
adiabaticity could be violated by increasing the magnitude of the damping term , @xmath142 .
then , however , a continuous acceleration of @xmath149 is required , because otherwise , if @xmath153 was constant , the system would equilibriate .    as an example for the absence of particle production inside a static trap , @xmath150
, i will consider an exponential sweep of the coupling coefficient @xmath154 with @xmath155 . in this case , eq .
for the scale factor permits an analytic solution @xmath156^{\frac1{2+d } } \exp\left\{\frac{\gamma t}{2+d}\right\ } \ , .
\label{bexpsweepg }      \end{aligned}\ ] ] the coefficients of eq .
become time - independent @xmath157 and the density eigenmodes are just damped harmonic oscillators with solutions @xmath158 where @xmath159 is some residual oscillating function with frequency @xmath145 .
hence , the density - density fluctuations diminish @xmath160 , and , consequently , the phase - phase fluctuations increase .
but this is just the adiabatic evolution , because @xmath161 and thus @xmath162 , cf .
[ secobservables ] .
this means that no quasi - particle production occurs for the dynamics in the hydrodynamic regime , i.e. , for low - energy excitations with order parameter treated in the thomas - fermi approximation .
these findings can also be inferred from the ( absence of an ) effective particle horizon . for the particular shape of the scale factor @xmath149
follows @xmath163 because @xmath164 .
note that the presented solution assumes @xmath165 at all times , especially also for @xmath166 .
hence @xmath167 at some time @xmath168 when @xmath169 , cf .
eqs .   and . on the other hand ,
@xmath170 and @xmath171 for a condensate at rest , see eq .  .
since both solutions ( static initial state and exponential sweep ) can not be matched at @xmath172 such that @xmath139 and @xmath149 are both continuous , the switching on of the exponential sweep would excite breathing oscillations .
these oscillations , however , generally affect particle production , e.g. , through parametric resonance .
the simplest application of the presented formalism consists in a quasi - one - dimensional condensate . in highly anisotropic traps , where the perpendicular trap frequency @xmath173 is much larger than the chemical potential , the motion in the perpendicular directions
is restricted to the ground state and might be integrated out .
an effectively one - dimensional field equation follows , where the interaction strength @xmath174 can be varied through feshbach resonance or by changing @xmath173 .
however , one should bear in mind that the transversal extent of the condensate @xmath175 has to be much larger than the @xmath4-wave scattering length @xmath5 such that the interaction of different atoms can still be described through three - dimensional scattering theory . for simplicity
, i will adopt in this section the thomas - fermi approximation for the order parameter but will permit arbitrary energies for the excitations .
calculated using 100 harmonic oscillator basis functions .
the frequencies equal the chemical potential at the intersection with the black line .
the lowest odd mode converges for @xmath176 to the classical thomas - fermi breathing frequency @xmath177 ( green ) , which follows from eq .  .
, scaledwidth=40.0% ]    in figure [ spectrum1d ] , the excitation frequencies @xmath178 of the lowest modes are plotted versus the chemical potential @xmath179 . for @xmath180 ,
one has the equidistant spectrum of the harmonic oscillator @xmath181 . whereas for high @xmath182
, the frequencies become almost independent of the chemical potential . in particular
, the frequency @xmath183 of the lowest excitation with odd parity tends for @xmath176 towards the thomas - fermi breathing frequency @xmath184 obtained from eq .  .
the discrepancy between these two frequencies @xmath183 and @xmath185 for finite values of the chemical potential @xmath179 hints at shortcomings of the thomas - fermi approximation .
in particular , eq .
does not describe the breathing motion of the background properly and care must be taken when employing eqs .   for the order parameter .
nonetheless , eq .   still predicts the correct order of magnitude of the characteristic response time of the background , @xmath186 , to variations of trap frequency or interaction strength .
hence , it is still possible to discuss several cases , where the difference of @xmath187 and @xmath188 is either small or does not matter : firstly , if the shape of the condensate varies only slowly , i.e. , if @xmath189 and the condensate can adapt to changes of @xmath138 and/or @xmath139 immediately .
secondly , if no breathing oscillations are excited , e.g. , because the condensate expands or contracts , see also eqs .   and . and , thirdly
, for very large chemical potentials , @xmath176 , the thomas - fermi approximation becomes exact .
the quantum fluctuations are in the hydrodynamic regime and their effective evolution equations decouple . in order to address cases where breathing of the background occurs
, it would be necessary to abandon the thomas - fermi approximation and to solve the gross - pitaevskii equation , which could , e.g. , be done by expanding the order parameter @xmath25 into oscillator functions @xcite .
a stationary condensate , i.e. , @xmath58 , can be accomplished through simultaneous variation of trap frequency @xmath138 and interaction strength @xmath139 , cf .
@xmath190 though , one should be aware that this only holds within the thomas - fermi approximation : the instantaneous chemical potential must at all times be much larger than the trap frequency @xmath191 if both were of the same order , the kinetic term , @xmath192 , in the gross - pitaevskii equation becomes important ; no stationary background could be realized even for simultaneous variation of @xmath138 and @xmath139 .      for an exponential sweep @xmath193 the effective second - order equation for the phase fluctuations in the hydrodynamic limit @xmath194 \delta\hat\phi \ , = \ , 0 \ , .
\label{effkfgexp }      \end{aligned}\ ] ] is that of a massless scalar field in a de  sitter spacetime with exponentially growing scale factor @xmath195 , cf .
eq .    which is believed to describe the universe during the epoch of cosmic inflation @xcite
for the time - dependence , the integral is finite and an effective sonic horizon occurs ; the quantum fluctuations freeze and get amplified .    instead of solving eq .   for @xmath196
, i will consider the evolution equation of the density fluctuations @xmath197 @xmath198 \hat x_n^+ \ , = \ , 0 \ , ,      \label{densitysweepeom }      \end{aligned}\ ] ] where @xmath199 .
obviously , all modes undergo the same evolution just at different times .
eq .   can be solved analytically in terms of bessel functions @xcite @xmath200 where @xmath201 .
the hankel functions @xmath202 have the proper asymptotics for early times @xmath203 such that the operators @xmath84 annihilate the initial vacuum state .
the phase fluctuations read @xmath204 from these expressions and , i can infer the correlations of each mode . at late times
follows @xmath205 comparison with the adiabatic values @xmath206 , cf.eq .  , yields the quasi - particle number at late times @xmath207 the occupation number of all modes grows exponentially though at different times @xmath208 , where the shift @xmath209 is determined by the excitation frequencies @xmath210 .
however , one should bear in mind that eq .
is only valid for a limited time before leaving the hydrodynamic regime .      in order to go beyond the effective space - time description and thus the analytical findings -
, i will now consider the full evolution equations .
the sweep rate @xmath211 shall be chosen such that all modes evolve adiabatically at first , i.e. , @xmath212 for all @xmath213 .
when subsequently reducing trap potential and coupling strength , the excitation frequencies @xmath214 decrease and non - adiabatic evolution sets in at different times for each mode .
the fluctuations freeze and get amplified .
figure [ sweep50 ] shows the instantaneous particle numbers of the lowest three modes for initial chemical potential @xmath215 and @xmath211 .
the lowest excitation , @xmath216 , which becomes non - adiabatic first , acquires the largest particle number .
the next two modes , @xmath217 , experience less squeezing , though , remarkably , @xmath218  an unexpected result , which can be explained by the coupling of different modes : the second even mode , @xmath219 , gets populated from the principle excitation , @xmath216 , whereas the coupling matrix elements between @xmath216 and @xmath220 are zero because of different parity .     and the sweep rate @xmath211 . at @xmath221 ,
the chemical potential equals the trap frequency . for the calculation of the eigenmodes @xmath68 and the coupling matrices , i used the harmonic oscillator functions @xmath222 up to @xmath223 , cf .
appendix [ seceigenmode ] .
the lowest 20 modes @xmath224 were then propagated .
the numerical accuracy was set to @xmath225 .
, scaledwidth=40.0% ]
the main objective of this article was the investigation of quantum fluctuations in time - dependent harmonically - trapped bose - einstein condensates with repulsive interactions . to this end , the linear fluctuations were expanded into their initial eigenmodes and the field equations were diagonalized . this diagonal form
, however , persists only as long as the condensate is at rest ; as soon as trap frequency or interaction strength are varied , the coupling of different modes sets in .
( only part of which can be accounted for by transforming to the instantaneous eigenmodes , though the definition of instantaneous eigenmodes is a non - trivial task . )
two regimes were identified : firstly , for energies much smaller than the chemical potential , the coupling of different modes is negligible and an effective space - time metric might be introduced for the phase fluctuations .
this , however , necessitates a redefinition of the time coordinate such that the required change of trap frequency and/or interaction strength for a certain dynamics of this effective space - time is not obvious .
it turned out that the sole variation of the interaction coefficient @xmath81 in a smooth monotonic way is hardly sufficient to render the evolution of the quantum fluctuations non - adiabatic , since the expansion / contraction of the background might compensate for changes of @xmath139 such that the sound velocity ( in comoving coordinates ) remains constant .
breathing oscillations of the background , excited , e.g. , by the sudden change of the interaction coefficient @xmath139 , might still yield a notable amount of quasi - particles . and secondly , if the excitation energy is of the same order as the chemical potential , different quasi - particle modes couple .
the amplification and freezing of the fluctuations and also the coupling of different modes was illustrated in an example , where the trap frequency and interaction strength were exponentially ramped down such that the shape of the condensate remains constant .
for the considered parameters , a quasi - particle number of 0.5 was obtained in the lowest mode , though higher occupation numbers could be achieved by faster sweep rates or starting with a higher chemical potential .
the inversion of the occupation number in the next two modes could be attributed to the inter - mode coupling : although the third excitation , @xmath219 , experiences a much briefer period of non - adiabatic evolution than the second mode , @xmath220 , only the former couples to the lowest mode , @xmath216 , and gets populated from it .
i would like to thank craig m.  savage and uwe r.  fischer for helpful discussions .
this research was supported by the australian research council .
the aim of this appendix is the derivation of the initial eigenfunctions @xmath68 of density and phase fluctuations . to this end , let me expand the field operator into any orthonormal basis @xmath226 of the underlying hilbert space @xmath227 with some operator - valued coefficients @xmath228 @xmath229 lower - case greek indices ( @xmath230 , @xmath231 , ... ) shall denote the components in this arbitrary basis , while lower - case latin indices ( @xmath232 , @xmath213 , ... ) label the initial eigenmodes .
if the condensate is initially at rest , the scale factor and the phase of the order parameter is homogenous @xmath233 such that @xmath61 and thus @xmath234 , cf .. the evolution equations simplify considerably and can be expanded into the basis @xmath235 .
it follows , cf .
@xmath236 where @xmath237 are defined in eq .  .
the matrices @xmath238 and @xmath239 are real and symmetric but do not commute due to @xmath240 \neq 0 $ ] , cf .  .
from eqs .  , i obtain second - order evolution equations for phase and density fluctuations @xmath241 which can be diagonalized by a transformation to the eigenvectors .    because the matrix @xmath242 is not
symmetric , phase and density fluctuations obey different eigenvalue equations @xmath243 where @xmath213 denotes the vectors and @xmath230 merely counts the components in the particular basis representation . for brevity
, i will only discuss the eigenvectors of @xmath244 in the following , though the same applies for those of @xmath245 as well .
the eigenvectors are generally not orthogonal @xmath246 but @xmath247 usually still forms a basis .
( note , however , that the eigenvectors of a non - symmetric matrix not always span the entire vector space .
but if the @xmath248 were no basis , the evolution equation could not be diagonalized .
this would mean that there existed some fluctuations , which have constant losses to some eigenmodes  rather unphysical in view of the stationary initial state considered here .
therefore , i will not discuss this case any further . ) with the @xmath249 forming a basis of @xmath227 , there must exist a dual basis @xmath250 in the space of linear functionals ( the dual ) on @xmath227 , which obeys @xmath251 [ roughly speaking , the elements of @xmath227 become column vectors in the basis expansion , whereas row vectors correspond to the functionals on @xmath227 , i.e. , the elements of the dual . since the eigenvectors @xmath248 form a basis , the matrix with the @xmath248 as columns must be invertible .
the rows of the ( left ) inverse matrix then comprise the elements of the dual basis @xmath250 in the particular basis expansion . ]    after the multiplication of the first of eqs .
with @xmath252 from the left and summation over @xmath230 follows @xmath253 which , because @xmath254 is a basis , implies that the @xmath255 are the left eigenvectors of @xmath256 with the same eigenvalues @xmath257 @xmath258 transposition yields @xmath259 i.e. , the eigenvalue equation of @xmath245 , cf .  .
hence , @xmath260 and @xmath261 and @xmath262 ( i.e. , density and phase fluctuations ) must have the same spectrum @xmath263 .
for simplicity , @xmath264 , which can be achieved by renormalization of @xmath265 , in the following .
note that the spectrum of a real non - symmetric matrix might contain pairs of complex conjugate eigenvalues @xmath257 , @xmath266 , which can be seen when taking the complex conjugate of .
complex eigenvalues @xmath267 are associated with exponentially growing solutions , i.e. , unstable modes .
since i am interested in the quantization of the stationary initial state , i will not discuss this case but instead assume @xmath268 @xmath269 .
( as can be easily verified , the components of the eigenvectors @xmath270 and @xmath271 must be real - valued as well . ) nonetheless , complex eigenvalues might still occur in dynamical situations , e.g. , during the signature change event proposed in @xcite or during ( quantum ) phase transitions @xcite ( see also ref .
@xcite for an illustrative example ) .
the initial evolution equations can be diagonalized by multiplication with the left eigenvectors @xmath272 which leads to the definition of the density and phase fluctuation eigenmodes @xmath224 through @xmath273 where the spatial mode functions @xmath274 follow from comparison of eqs .   with and the duality of the @xmath270 and @xmath271 implies eq .   for the @xmath68
.    in view of
, the initial evolution equations can be transformed to the new basis , cf.eq .
@xmath275 where the transformed matrices @xmath276 have become diagonal .
note that the transformation is not orthogonal because the matrices comprising of the eigenvectors @xmath270 and @xmath271 are not orthogonal .
hence , the commutators are _ not _ preserved , in particular @xmath277 where the latter can be inferred from eq .
@xmath278 since @xmath106 and @xmath107 also commute with their product @xmath279 , which is diagonal , they must be diagonal , too .
the duality condition does not fix the norm of the eigenvectors but still permits the multiplication with an arbitrary factor @xmath280 @xmath281 this renormalization then leads to a stretching / shrinking of the operator - valued coefficients , cf .
@xmath282 since this is merely a basis transformation , the time - evolution of the quantum fluctuations @xmath224 must be unaffected . to see this , recall the definitions and of the coupling matrices : they are the matrix elements of the operators @xmath237 and of @xmath283 with respect to the basis functions @xmath68 and thus acquire additional factors as well .
as expected , all of these factors cancel such that the time - evolution remains unchanged . in particular , the eigenfrequencies are invariant @xmath284 .
furthermore the observables should not depend on the particular normalization of the basis functions .
the relative density - density correlations at time @xmath140 read @xmath285 where i used the instantaneous eigenmode basis @xmath286 , see eq .. obviously , the factors @xmath280 and @xmath287 contributed by @xmath286 and @xmath288 cancel each other and the spatial correlations are independent of the normalization @xmath280 .
similarly , the expression for the spatial phase - phase correlations @xmath289 yields the same result regardless of the employed basis .
hence , the factors @xmath280 can be chosen at will . there
exist , however , several convenient choices for @xmath280 : firstly , the density modes might be normalized to unity , @xmath290 .
this is advantageous if left and right eigenvectors are the same , i.e. , if @xmath244 is symmetric .
the drawback is that if @xmath244 is not symmetric and therefore @xmath291 , only one of the mode functions , @xmath292 , can be normalized to unity , whereas the norm of the @xmath293 follows from eq .  . and ,
secondly , these factors @xmath280 can be fixed by demanding @xmath294 . in this case , the prefactors of initially obey @xmath295 and @xmath296 and both quadratures @xmath297 .
however , one should note that neither of the eigenfunctions @xmath68 is generally normalized to unity , @xmath298 .    of course , the coefficients @xmath299 and @xmath300 do depend on the particular choice of the basis functions @xmath301 and thus also on the factors @xmath280 .
for instance , the adiabatic density - density and phase - phase correlations of a particular mode read @xmath302 the dependence on the factor @xmath280 becomes important regarding the low excitations in the thomas - fermi approximation , cf .
[ sectf ] : the modes decouple and it is not necessary to introduce a new spatial basis @xmath68 at the time of measurement .
one has instead @xmath303 and @xmath304 such that @xmath305 while @xmath306 , i.e. , the phase and density fluctuations apparently increase or decrease even for adiabatic evolution . in view of this @xmath280 ambiguity ,
the ( absolute ) density - density or phase - phase correlations provide no adequate measure for the squeezing ( i.e. , non - adiabaticity ) of a single mode .
the fourier transforms of eqs .   and on the other hand are independent of @xmath280 but do not represent the excitation eigenmodes .
as will be shown in the following , the ( instantaneous ) quasi - particle number measures the relative deviation of density and phase fluctuations from their adiabatic values . in view of the different expansions and of density and phase fluctuations into their initial and adiabatic eigenfunctions , the corresponding creation and annihilation operators @xmath83
, @xmath84 and @xmath307 , @xmath308 , respectively , can be transformed by virtue of a bogoliubov transformation .
for the annihilators @xmath308 follows in particular @xmath309 with the bogoliubov coefficients @xmath310 and @xmath311 .
the first line follows from inversion of eq .   and
the second line can be inferred after transforming to the new basis @xmath301 . since @xmath312 for non - adiabatic evolution ,
the quasi - particle number operator @xmath313 will have a non - zero expectation value as well @xmath314          \right\rangle      \end{aligned}\ ] ] where the commutator @xmath315 = -i$ ] .
noting that the prefactors @xmath316 and @xmath317 are just the adiabatic density - density and phase - phase correlations , see eq .
, the particle number can be rewritten @xmath318 i.e. , the instantaneous particle number gives just the relative deviation of the density and phase correlations from their adiabatic values .
note that expression does not contain the correlations between different modes . to this end
, it would be necessary to evaluate @xmath319 , which is fourth order in the @xmath308 .
this observable can be reduced to expectation values quadratic in in the @xmath308 by virtue of wick s theorem , see , e.g. , @xcite .
another interesting aspect regards the scaling of the correlation functions and with the interaction strength : within the thomas - fermi approximation , see sec .
[ sectf ] , the evolution equations are independent of @xmath13 .
all properties of the linear excitations are determined by the chemical potential @xmath179 and the variations @xmath320 and @xmath321 , in particular the expectation values @xmath322 .
only the normalization factor @xmath323[\mu_0 - v(\bm x')]}$ ] in eqs .   and depends on @xmath13 .
hence , the relative density - density and phase - phase correlations in the thomas - fermi approximation are both proportional to @xmath13 for fixed @xmath179 .
a.  fisher , p.  b.  weichman , g.  grinstein , and d.  s.  fisher , phys .
b * 40 * , 546 ( 1989 ) ; d.  jaksch _ et al .
_ , phys .
lett .  * 81 * , 3108 ( 1998 ) ; d.  jaksch and p.  zoller , ann .
( n.y . ) * 315 * , 52 ( 2005 ) ; r.  schtzhold , m.  uhlmann , y.  xu , and u.  r.  fischer , phys .  rev .
lett .  * 97 * , 200601 ( 2006 ) ; c.  kollath , a.  m.  luchli , and e.  altman , phys .  rev .
* 98 * , 180601 ( 2007 ) .              c.  barcel , s.  liberati , and m.  visser , phys .
a * 68 * , 053613 ( 2003 ) ; c.  barcel , s.  liberati , and m.  visser , int .
j.  mod .
d * 12 * , 1641 ( 2003 ) ; c.  barcel , s.  liberati , and m.  visser , class .  quant .
grav .  * 18 * , 1137 ( 2001 ) ; _ ibid . _ * 18 * , 3595 ( 2001 ) .
dimensionless units can be introduced through the initial trap frequency @xmath324 and the atom mass @xmath232 .
all energies are then measured in units of @xmath324 , the dimensionless time is defined through @xmath325 , and the characteristic length @xmath326 , where @xmath327 .
e.  tiesinga , b.  j.  verhaar , and h.  t.  c.  stoof phys .  rev .
a * 47 * , 4114 ( 1993 ) ; s.  inouye _ et al .
_ , nature * 392 * , 151 ( 1998 ) ; ph .
courteille , r.  s.  freeland , and d.  j.  heinzen , phys .
lett . * 81 * , 69 ( 1998 ) .            as similar expansion of the linear fluctuations is also considered in @xcite . there
, the authors expand the linear field operator @xmath27 of a static / stationary condensate into eigenmodes and not the phase and density fluctuations @xmath62 .
this leads to complex mode functions @xmath328 and @xmath329 , while @xmath68 are real .
the eigenfunctions of density and phase fluctuations @xmath330 and @xmath45 are defined slightly differently than those of @xmath62 @xmath331 noting that @xmath332 the center of the trap , the eigenfunctions @xmath292 and @xmath333 ( and similarly @xmath293 and @xmath334 ) must be approximately proportional for low excitations .
hence , the @xmath335 for @xmath336 indeed represent the proper density and phase eigenmodes @xmath337 and @xmath338 up to a constant prefactor @xmath339 .
l.  e.  sadler _ et al .
_ , nature * 443 * , 312 ( 2006 ) ; h.  saito and m.  ueda , phys .
a * 72 * , 023610 ( 2005 ) ; h.  saito , y.  kawaguchi , and m.  ueda , phys .
lett .  * 96 * , 065302 ( 2006 ) ; a.  lamacraft phys .
lett .  * 98 * , 160404 ( 2007 ) ; m.  uhlmann , r.  schtzhold , and u.  r.  fischer , phys .
* 99 * , 120407 ( 2007 ) .