Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the idea that the images of distant objects might be distorted and even multiply imaged by the gravitational potentials of intervening objects is one of the earliest predictions of general relativity . gravitational lenses , perhaps better described as gravitational mirages , are now proving powerful tools in cosmological studies . the equations of general relativity can be cast in a form where the gravitational potential acts like an index of refraction . as in the case of terrestrial mirages , variations in this index of refraction cause distortions of bundles of light rays , and occasionally permit viewers to see several such bundles which were emitted in different directions . the result is a multiply imaged ( or strongly lensed ) gravitational lens system . perhaps 50 distant active and normal galaxies are known , at present , to be strongly lensed by mass concentrations lying roughly along what would otherwise be the shortest path between the galaxy and the earth . the most spectacular case , shown in figure 1 , is that of the cluster of galaxies cl0024 + 1654 which produces five readily recognized images of a more distant ring galaxy . in this image the background lensed object is identified by its blue colour , while the cluster member galaxies appear reddish yellow . the cluster distorts each of the five images somewhat differently , but these distortions are easily understood as the effect of the cluster s gravitational potential .    one important use of such systems is to map the gravitating mass in the lensing object . another use is to magnify features in distant objects which would otherwise be too small to be resolved . several star forming galaxies at very high redshift have been studied in unprecedented detail thanks to the resolution of these megaparsec - scale natural telescopes . in addition to studying individual objects , be they sources or lenses , multiply imaged systems can be used to derive global cosmological parameters : the matter density of the universe , @xmath1 ; the cosmological constant , or the vacuum contribution to the energy content of the universe , @xmath2 ; and the expansion rate of the universe , @xmath0 . the first two of these dictate the overall geometry of the universe , i.e. how distances to objects relate to their measured redshifts . because lenses and sources populate a wide range of redshifts , statistical properties of large samples of lensed objects , such as the distribution of their source and lens redshifts and angular separation of images , depend on the geometry of the universe . for example , a non - zero cosmological constant increases the volume of space at high redshifts , hence the abundance of distant lensed systems can be used to estimate @xmath2 .    while statistics of complete samples of lensed objects are needed to measure @xmath1 and @xmath2 , a single lensed system , such as a multiply - split high redshift qso ( quasi - stellar object , or quasar ) can in principle yield the value of @xmath0 . in such a situation the light travel times for the bundles of light from multiple images are different , permitting the observer to see the emitting source at two different moments in time . the resultant time difference , or time delay , depends upon the differences in the path lengths among the multiple light paths . because all of the angles are known in many lensed systems , a difference in path lengths is sufficient to determine all distances , and thereby fix the linear scale of the universe . this `` time delay '' phenomenon and its utility for the determination of @xmath0 was the subject of a recent workshop at the nuffield radio astronomy laboratories , jodrell bank , part of the university of manchester s department of physics and astronomy . virtually everything one might want to calculate for a gravitational lens can be computed from the expression for the light travel time for a photon traveling from a source to an observer ( narayan and bartelmann 1996 , blandford and kundi 1996 ) . as can be seen from the equation in figure 2 , there are two contributions to the light travel time : a geometric contribution , proportional to the square of the difference in angular position between where the source would have appeared and where the image does appear , and a gravitational contribution , which results from the gravitational potential altering the effective index of refraction . for a `` thin '' lens , the gravitational effects can be lumped together as a 2d projected potential by integrating the 3d potential along the line of sight . a 2d plot of this light travel time is called the arrival time surface ; the example in figure 3 shows isochrones . fermat s principle tells us that images of the source will appear wherever the light travel time is a minimum , a maximum or a saddle point . to find these stationary points one takes the gradient of the light travel time with respect to position on the sky and sets it equal to zero . this gives a vector equation , often called `` the lens equation '' , for the positions of the images on the plane of the sky . there may be one , three , five or more solutions to this equation , giving an odd number of images . applying the lens equation to the arrival time surface of figure 3 one can readily identify the positions of five images , marked as dots . the cross is the position of the source in the absence of the lens . the time delay between the images is just the difference in `` height '' of the arrival time surface at the image positions .    taking the gradient of each cartesian component of the lens equation separately , one finds expressions for the partial derivatives of the source position with respect to the image position . assembling these terms ( which include second derivatives of the projected potential ) as a matrix gives the `` back projection matrix '' , telling what size source corresponds to what size image . more useful is its inverse , `` magnification matrix '' , which tells us how a source is distorted . an important point here is that for a potential which is sharply peaked at its centre , the back projection factor is large ( because of the large second derivatives of the potential ) , and the magnification small . while in principle one should see odd numbers of images , in practice the images which ought to appear near the centres of the lensing galaxies ( maxima in the light travel time ) are so strongly demagnified that they are rarely detected . hence virtually all the lenses one sees are classified as doublets or quadruplets , in seeming contradiction of the requirement that there be an odd number of images . figure 3 shows the shapes and sizes of the five images produced by our lens from a small circular source . because lensing conserves surface brightness of the source ( a consequence of liouville s theorem ) , the apparent sizes of the images are directly proportional to their magnifications . note that  as for actually observed galaxy lenses  our lens demagnifies the central image so that it is barely detected . the light travel times are proportional to the distance from the observer to the lens and the distance from the observer to the source and are inversely proportional to the distance from the lens to the source . the relative sizes of these quantities can be determined from the redshifts of the lens and the source , but the overall scaling factor depends upon the hubble constant , the value for which continues to be a matter of some debate . if the difference in light travel times can be measured , and if the gravitational potential can be modelled , one can solve for the hubble constant . some 18 years have passed since the discovery of the first gravitationally lensed quasar by walsh and collaborators ( walsh _ et al . _ 1979 ; walsh , 1989 ) . if measuring distances to such systems were really so straightforward , one might reasonably wonder why questions regarding the cosmological distance scale havent long since been resolved . as one might have guessed from the name of the workshop , not all lenses are equally useful  hence the quest for the `` golden lens '' . perhaps the most frustrating fact is that , despite their reputation , most quasars are surprisingly non - variable ; there are beautifully multiply - imaged quasars which fail to show appreciable flux variations over the timescale expected for a time delay . other strongly lensed systems are so symmetrical that the time delays are both very short and very sensitive to small errors in the measurement of the relevant positions . still other systems are quite promising except for the fact that either the lensing galaxy or the lensed quasar is too faint for a redshift measurement .    there have been sociological hurdles as well . the scheduling of 50 or 100 appropriately spaced observations for a time - delay monitoring program is not always easy , especially with changes in instrument configurations and the traditional scheduling scheme at optical observatories whereby a single program is given a contiguous block of several nights . even if the foregoing problems are overcome and a near - perfect lens system is found with its properties carefully measured , there still remains a challenge of constructing a model for the gravitational potential of the system . the task is not trivial because the observables provide a rather limited set of constraints for the model , leaving us to rely on various assumptions about what a galaxy should look like in general . in addition , it often happens that the main lensing galaxy is embedded in a group or cluster of galaxies , like in the case of q0957 + 561 . such additional mass concentrations around the lens make model building even more difficult , and the resultant @xmath0 even more uncertain .    because both observational and theoretical aspects of the problem are challenging and time - consuming there has been something of a chicken - and - egg problem , with observers reluctant to embark upon a monitoring campaign without first having a good model , and model builders reluctant to make detailed models before time delays have been measured . one has to start somewhere , and the recent class survey goes a long way towards addressing the observational challenges . the past couple of decades have produced a number of large comprehensive surveys in astronomy . all the `` early '' lenses were serendipitously discovered during some of these surveys ; for example the original double quasar q0957 + 561 a , b was discovered during the optical spectroscopic follow - up of jodrell bank 966 mhz sources . the quadruple system q2237 + 0305 , another famous member of the lensing family , was discovered in the center for astrophysics ( cfa ) galaxy redshift survey . only a handful of surveys have been specifically designed to uncover new lenses , or have finding new lenses as one of their main goals . the mit - green bank 5 ghz , parkes - mit - nrao 5 ghz , and the hubble space telescope ( hst ) snapshot surveys are among these .    by far the largest of such surveys is the ongoing class , an acronym which stands for the cosmic lens all sky survey ( myers _ et al . this effort is the follow - on to the jodrell bank vla astrometric survey ( jvas , patnaik 1993 ) . the class collaborators , including astronomers from jodrell bank , the netherlands and the us , have obtained radio images of some 10 000 flat spectrum radio sources using the very large array in socorro , new mexico . only by heavily automating their data reduction were the investigators able to sift through the enormous volume of data gathered to identify the most likely candidates . a dozen or so of these have proven probable lenses . the class candidates have been followed up at higher angular resolution using the merlin interferometric array , and at optical and ir wavelengths using the william herschel telescope on la palma and the keck telescope on mauna kea . optical images of many of these newly discovered lenses have been taken with the hst , most of them clearly showing the galaxy responsible for the lensing . new lenses from the class and jvas surveys have brought new surprises . contrary to the theoretical wisdom of the past decade , most of the galaxy lenses in this sample are classified as spirals or s0s and not ellipticals . the two best examples of spiral galaxy lensing are the objects b0218 + 357 and b1600 + 434 . the disk in the latter case appears to be almost exactly edge - on . spurred by these observations , theorists have now turned to the problem of modelling spirals . unlike elliptical galaxy models whose potential is well represented by some variation on the basic isothermal sphere , the spiral galaxies also have a flattened disk which has to be taken into account because the lensing cross section increases dramatically when disk is turned edge - on . it appears that q0957 + 561 , the first lens to be discovered , was not unique in being a member of a group or cluster of galaxies . several other lenses have companions in close projection to the main galaxy lens , some with confirmed matching redshifts . in hindsight this is not surprising  galaxies are known to cluster , and the larger optical depth presented by close collections of galaxies increases magnification , thereby enhancing the chance that the background object is included in flux limited surveys . the class has also uncovered a few unusual and puzzling lenses . perhaps the most spectacular of these is the new lens b1933 + 504 , which shows a total of 10 radio images . observationally the images can be grouped into those of similar radio spectra and compactness . with the help of analytical models , one concludes that one is seeing three components of an active radio galaxy ; the compact flat spectrum central core and one of the steep spectrum lobes are both quadruply imaged , while the second radio lobe is doubly imaged . the abundance of multiple images makes this system ideal for modellers ; it fails to qualify as a `` golden lens '' because the lack of an optical counterpart , even with the hst , means that the redshift of the source is unknown . the jvas lens b0218 + 357 is an excellent system for time delay measurements . it consists of two components which have been observed to vary in radio flux ( vla 15 ghz ) , the degree of polarization , and the polarization position angle . as the variations can be tracked by all three methods , and are relatively large , 10% , it should be possible to measure the time delay to high fractional accuracy . the current estimate , @xmath3 days , is likely to improve rapidly . the modelling of this system should be helped by the existence of a radio ring , and the fact that the lensing galaxy , a face - on spiral , is clearly visible . the galaxy contains neutral and molecular gas , seen in absorption against the background radio images . however , before b0218 + 357 is modelled in detail two aspects of the systems must be understood : the separation of the two radio images and two optical components are different by 10% ( significant given the size of the error - bars ) , while the flux ratio of the images differs by a factor of nearly 40 between the optical and the radio . dust , known to exist in spiral galaxies , may be in part responsible for these observations . jvas lens b1422 + 231 , a quadruple system with three very bright images and a fourth faint one , presents a serious problem for modellers . the image positions can be easily fitted by a range of models , but when radio flux ratios are taken into account there is no acceptable analytical fit for the galaxy . it is generally thought that radio flux ratios should provide useful additional constraints for the overall potential of the galaxy . this is because the radio - emitting region of the source quasar is substantially larger than the lensing cross - sections of individual stars in the lensing galaxy , so that microlensing by stars should not introduce flux variations in the images above and beyond those attributed to the smooth potential of the galaxy as a whole . moreover , radio flux is not obscured by dust . thus it appears either the models tried so far are inadequate for b1422 + 231 , or that its substructure significantly affects image magnifications . class b1608 + 656 , shown in figure 4 , is another potential `` golden lens '' . the source is unusual because it is not a quasar but a post - starburst radio galaxy . the extended nature of the galaxy gives rise to a broken ring observed in the lensed image . the system is being monitored with the vla , and the time delays are reported to be of the order of a few weeks . a potential problem is that the lensing galaxy does not seem simple : it is either an interacting system , or its image is split by an obscuring dust lane . the present sample of class lenses , though exciting and fascinating , is not yet a complete sample in a statistical sense . this is because the quadruplets and higher multiplicity systems are much easier to confirm as lenses using geometry , without the help of optical counterparts and redshifts . double systems , as well as larger separation lenses with several arcseconds between images , and lenses without obvious galaxies are more difficult to confirm as lenses . as @xmath0 can be determined from an individual lens system , ( or a collection of lenses to take into account modelling and observational errors ) the current class results can already be used for that purpose , while some other cosmological applications ( those that require statistically homogeneous samples ) must wait until the class survey is complete . the progress reported at the jodrell bank meeting was along three principal fronts . the first of these was the discovery of a large number of new lenses , the most interesting of which are described above . the second area in which considerable strides have been made has been in the measurement of time delays . the relatively small amplitude of flux variations demands scrupulous attention to calibration in carrying out observations . by the end of the meeting the number of secure measurements of time delays had doubled from two to four . the best determined time delay is that for the oldest known system , q0957 + 561 . using two year s worth of data obtained in 1995 - 1996 at the apache point observatory in new mexico and a range of different methods to analyze the light curves ( see figure 5 ) , kundi , wambsganss and collaborators ( kundi _ et al . _ 1997 ) found the time delay between a and b to be @xmath4 days , confirming the earlier , less robust result of schild and cholfin ( 1986 ) . the time delays in a quadruple pg1115 + 080 , the second gravitational lens to be discovered , are currently determined to @xmath5 accuracy . this is likely to improve with the addition of more data . the lenses b0218 + 357 and b1608 + 656 now have measured time delays , from radio observations , for which the fractional uncertainty is also likely to improve rapidly .    a third front on which progress was reported was in modelling the gravitational potentials , which make a substantial contribution to the total time delay . the potentials are constrained by the observed pattern of multiple images . but there are few images , so the potentials are poorly constrained . the lens equation shows that the image positions depend upon derivatives of potential rather than on the potential itself . but it is the potential that determines time delays ; reconstructing it from a small number of derivatives is an ill - posed problem .    in a few cases some extra constraints to aid reconstruction are provided by the additional images observed in the system . for example , vlbi images of q0957 + 561 show a core and a radio jet with five blobs , each of which is also doubly imaged . these extra constraints are much needed in this system where the lens consists of a galaxy plus a parent cluster and hence requires more parameters for even minimal modelling . there is an intrinsic degeneracy in the modelling because from the image data alone we do not know how the total mass causing the image separation is split up between the galaxy and the more smoothly distributed cluster mass . this `` mass - sheet degeneracy '' which translates into an uncertainty in @xmath0 , can be broken in a number of ways . one is by observing additional images , such as jet components in q0957 + 561 . another , is by studying the lens itself : a measured velocity dispersion within the main galaxy can be , albeit with caveats , converted into its mass . alternatively , a weak shear map of the background galaxies , or the velocity dispersion of galaxies within the cluster , can be used to estimate cluster mass . it is encouraging that different ways of breaking the degeneracy in q0957 + 561 lens lead to nearly the same estimates of @xmath0 .    a further strategy to surmount the problem of limited constraints is through the modelling of the lensing mass . some groups have adopted `` reasonably '' parameterized galaxy models , while others constructed non - parametric models imposing additional constraints requiring smoothness and a degree of central concentration . the ultimate test is whether modellers using different techniques agree on the predicted time delays . agreement appears to be best in the case of q0957 + 561 , which is also the lens with the best measured time delay . the current status of lens time - delay measurements is given in table 1 , where we have listed separately the fractional uncertainties associated with the time - delay measurements and with the models for each of the four lenses with secure time delays . in the last column we show the derived value of the hubble constant for three of the four systems . agreement is quite good , but with the simplifying assumptions made in modelling , a mere three systems seems hardly adequate . the onus for improving this state of affairs falls equally on those measuring delays and on those making models . gravitational lens time delays are hardly the only means of measuring extragalactic distances and there has been considerable talk recently ( kennicutt 1996 ) about the convergence of the many different methods . the crucial rung on what is known as the `` cosmological distance ladder '' is the period - luminosity relation of cepheid variable stars ( jacoby _ et al . their use as distance indicators dates back to the historic curtis - shapley debate of 1920 ( trimble 1995 ) where cepheids were used , though somewhat erroneously , to argue _ against _ the extragalactic nature of `` spiral nebulae '' . the modern distance ladder ( really a coupled set of ladders ) relies ultimately on parallax measurements , but includes as a key link the calibration of the period - luminosity relation for cepheids . there is considerable debate at present about the extent to which this relation depends upon the abundances of elements heavier than helium , since these abundances can be quite different in other galaxies from the values measured locally in the milky way . moreover the distances to milky way cepheids have become more uncertain rather than less uncertain with the first published results from the hipparcos astrometric satellite .    unlike the distance ladder , the gravitational lensing method measures distances on truly cosmological scales , and does so directly , bypassing all the rungs of the ladder together with their associated errors . one other method can claim that : the sunyaev - zel@xmath6dovich technique basically compares the line of sight extent of the cluster , estimated from the properties of its hot x - ray emitting gas , to the cluster s @xmath0-dependent projected size . systematic uncertainties associated with projection effects and gas clumpiness tend to underestimate @xmath0 .    while at present the gravitational lenses give distances which are no more accurate than the other distance indicators , they are at least subject to a completely different set of systematic errors . the major recognized source of these is lens modelling , including the mass - sheet degeneracy . these systematic uncertainties have not yet been explored mostly because lens models used until recently were all rather similar . the new and better data from the class survey and other campaigns have prompted a more comprehensive look at modelling , a trend which will soon allow us to assess critically the associated errors . the mass - sheet degeneracy can in principle be mostly resolved with independent measurements of the smooth cluster mass around the images .    given that all methods are subject their own systematic uncertainties , it is essential to pursue all the available methods in order to reach a convergence . several requirements need to be met in order to measure @xmath0 from gravitational lensing . first , a multiply imaged lens system has to be found that possesses a number of desirable characteristics , like well resolved images , measurable source variability , easily visible lensing galaxy and optical image counterparts . once a system is found , it must be extensively observed to get well sampled image light curves and accurate astrometric information . several of the new class lenses and several previously known satisfy these requirements to various degrees ; a few appear to be potential `` golden lenses '' . next , a good analytical model for the gravitational potential of the lensing galaxy/ galaxies is needed . because the observed parameters of the lens are not nearly enough to constrain the model properly , many different approaches must be explored , so as to quantify carefully the uncertainties arising from modelling . since several parametric and non - parametric modelling techniques have been published recently a comparison of these will soon be possible .    finally , after random and systematic errors have been assessed , the comparison of model predicted and observed time - delays between images will yield the value of the hubble constant . taking stock of the recent observational and theoretical developments , the prognosis for the gravitational lens method is excellent .                                  .observed and predicted uncertainties in the time delays between images in four gravitationally lensed systems with most secure measured time delays . the best estimate for the @xmath0 is quoted in the first three cases . [ cols="^,^,^,^ " , ]
Abstract:

Answer: gravitational lensing is now widely and successfully used to study a range of astronomical phenomena , from individual objects , like galaxies and clusters , to the mass distribution on various scales , to the overall geometry of the universe . here we describe and assess the use of gravitational lensing as `` gold standards '' in addressing one of the fundamental problems in astronomy , the determination of the absolute distance scale to extragalactic objects . this is commonly parameterized by the hubble constant , @xmath0 , the current expansion rate of the universe . the elegance of the underlying geometrical principle of the gravitational lensing method combined with the recent advances in observations and modelling makes it a very promising technique for measuring @xmath0 .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: ultracold bosonic systems have attracted much interest since the successful observation of bose - einstein condensation in @xmath0rb atoms  @xcite . one of the interesting topics in this field is an optical lattice system  @xcite , which is formed by loading the ultracold atoms in a periodic potential . this gives us clean correlated bosonic systems on various lattices with controllable parameters . in fact , remarkable phenomena have been observed such as the phase transition between mott insulating and superfluid states on cubic  @xcite and triangular lattices  @xcite . recently , the possibility of a supersolid state , where solid and superfluid states coexist , has been discussed . the existence of the supersolid state has experimentally been suggested in the @xmath1he system  @xcite , and has theoretically been discussed in strongly correlated systems such as bosonic  @xcite and fermionic  @xcite systems and bose - fermi mixtures  @xcite . as for a hard - core bosonic model , it has been clarified that lattice geometry as well as intersite correlations play a key role in stabilizing the supersolid state  @xcite . the hard - core bose - hubbard model on a triangular lattice has been discussed in detail  @xcite . in the weak coupling region , the superfluid state is realized . on the other hand , in the strong coupling region , the solid state appears at a one - third ( two - thirds ) filling , where one of three sites is filled ( empty ) in a @xmath2 ordering with wave vector @xmath3  @xcite . between the superfluid and solid states , the supersolid state is realized around half - filling , where a part of particles crystallize and the others form the superfluid state  @xcite . by contrast , no such supersolid states have been studied in the bose - hubbard model on layered triangular lattice . this should be crucial to observe the supersolid state in the realistic optical lattice systems  @xcite . an important point is that the supersolid state on each layer has three - fold degeneracy for particle distributions . therefore , it is desired to discuss how the interlayer couplings affect the stability of the supersolid state .    to clarify this point , we systematically study ultracold bosons in an optical lattice with the layered triangular structure , which is schematically shown in fig . [ fig : layered - model ] ( a ) . we discuss ground - state properties of this model , combining cluster mean - field ( cmf ) theory  @xcite with the density matrix renormalization group ( dmrg ) method  @xcite . the rich phase diagram is obtained where the solid and supersolid states with various particle distributions compete with the superfluid state . we demonstrate that the supersolid state is realized even at half filling , in addition to the well - known supersolid states discussed in a single layered model  @xcite . we also discuss the nature of the supersolid states . the paper is organized as follows . in  [ 2 ] , we introduce the bose - hubbard model for the layered optical lattice system and explain the possible particle distributions for the solid and supersolid states . we briefly outline our theoretical approach . in  [ 4 ] , we study ground - state properties in the layered system and discuss the stability of the supersolid states . a summary is given in the final section . we consider ground - state properties of interacting bosons on the layered triangular lattice . here , we assume sufficiently large onsite interactions . in the case , the system should be described by the following hard - core bose - hubbard model as , @xmath4 where @xmath5 denotes the summation over nearest neighbor sites and @xmath6 is the creation ( annihilation ) operator of a boson at site @xmath7 . the hopping @xmath8 and the intersite interaction @xmath9 when the sites @xmath7 and @xmath10 are located in the same ( distinct ) layer . this model is symmetric when interchanging particles with holes . the symmetric condition is given by @xmath11 , where @xmath12 is the chemical potential and @xmath13 . in the paper , we restrict our discussions to the system with @xmath14 , where the particle density @xmath15 and @xmath16 is the total number of sites .    when the interlayer couplings are zero @xmath17 , the model hamiltonian eq . ( [ hardcore - bose - hubbard ] ) is reduced to the hard - core bose - hubbard model on a single - layer triangular lattice  @xcite . in the limit , it is known that the supersolid states are realized , which are surrounded by the superfluid state in the weak coupling region and solid states with @xmath18 and @xmath19 in the strong coupling region  @xcite . note that at the particle - hole symmetric case , two supersolid states with distinct particle density are degenerate . therefore , when the system is fixed at half filling , the genuine supersolid state is not stable , but the phase separation between two supersolid states appears in the strong coupling region . when @xmath20 , the model is reduced to the one - dimensional bose - hubbard model , where no supersolid state appears  @xcite . in the paper , we study the hard - core bose - hubbard model on the layered triangular lattice to discuss how the supersolid states is realized . first , we define the order parameters for the possible ground states . in the layered system with intersite interactions , various type of particle distributions are naively expected . in fact , in the simple ( classical ) limit with @xmath21 , a phase diagram should be obtained , as shown in fig . [ fig : layered - solid ] . when @xmath22 ( @xmath23 ) , an empty ( fully - occupied ) state appears . when @xmath24 ( @xmath25 ) , the solid state appears with @xmath26 ( @xmath27 ) . if one focuses on a certain layer , the solid state has three - fold degeneracy . this means that the layered system has a large residual entropy in this limit . here , we specify the solid structure with a unit cell composed of six and nine sites as @xmath28- and @xmath29-type solids , respectively ( see fig . [ fig : layered - solid ] ) . note that we did not find any particle distributions with larger periods when the hoppings @xmath30 are introduced . the @xmath28- and @xmath29-type particle distributions are characterized by the following quantities as @xmath31 where @xmath32 on the other hand , at half filling @xmath33 , a different particle distribution appears . [ fig : layered - solid ] shows that the particle density in each layer alternate in the @xmath34-direction . the state is characterized by the following quantity , in addition to @xmath35 , as @xmath36 here , we specify this particle distribution as the `` @xmath37''-type one , for simplicity . in the weak coupling limit , the superfluid ( sf ) state is naively expected . the order parameter is defined as @xmath38    when the system has both @xmath39 and some of the other parameters defined above , there appears the coexistence between the superfluid and particle density wave , and we can say that the supersolid state is realized . for simplicity , we specify the solid ( supersolid ) states as the s-@xmath28 , s-@xmath29 , and s-@xmath37 ( ss-@xmath28 , ss-@xmath29 , and ss-@xmath37 ) , corresponding to the particle distributions . the relation between the ordered states and quantities are summarized in table [ tab : solid - order - parameter ] . .the order parameters of the possible ground states [ cols="^,^,^,^,^,^,^,^",options="header " , ]     to discuss ground - state properties in the hard - core bose - hubbard model , we make use of the cmf method  @xcite . in the method , the original lattice model is mapped to an effective cluster model , where particle correlations in the cluster can be taken into account properly . the expectation values of the intercluster hamiltonian are obtained via a self - consistency condition imposed on the effective cluster problem . this method has an advantage in discussing quantum phase transitions correctly since not only stable states but also metastable states can be treated . therefore , the cmf method has successfully been applied to the quantum spin systems  @xcite and bosonic systems  @xcite .    in the study , we consider the three - leg tube model with the length @xmath40 as the effective cluster model , as shown in fig . [ fig : layered - model ] ( b ) . the effective hamiltonian is explicitly given as , @xmath41\nonumber\\    \label{cmf - hamiltonian}\end{aligned}\ ] ] where the symbol @xmath42 [ @xmath43 denote the summation over nearest neighbor sites in the cluster ( between clusters ) . note that this three - leg tube is too thin to quantitatively discuss quantum phase transitions in the system  @xcite . however , it is known that an effective three - site cluster in a single layer reproduce the reasonable phase diagram . an important point in this study is that the cluster with a large value of @xmath40 is necessary to discuss the stability of the solid and supersolid states with long - period structures . to this end , we here make use of the dmrg technique as an effective cluster solver  @xcite . it is known that this method is powerful for the one - dimensional systems  @xcite . furthermore , by combining the dmrg method with a mean - field theory , the phase transitions in higher dimensions have been discussed for the heisenberg models  @xcite and fermionic hubbard models  @xcite . some details of the cmf+dmrg method are explained in ref . @xcite . in the following , we take the intralayer interaction ( @xmath44 ) as unit of energy and fix the ratio of interactions as @xmath45 . in the section , we discuss ground state properties in the hard - core bose - hubbard model on the layered triangular lattice . we then clarify the role of the interlayer coupling in stabilizing the supersolid state in the system . first , we consider the hard - core bosonic system with a low particle density , fixing the chemical potential as @xmath46 . combining the cmf method with the dmrg , we obtain the order parameters for the system with a fixed ratio @xmath47 , as shown in fig . [ fig : low1 ] . ( open circles ) , @xmath48 ( open squares ) , the superfluid order parameter @xmath39 ( solid circles ) , and the particle density @xmath49 ( solid triangles ) as a function of @xmath50 in the system with @xmath47 and @xmath51 . , width=302 ]    in this case , the particle density in the layer is always uniform in the @xmath34-direction ( @xmath52 ) . when @xmath53 , the particle density is one - third ( @xmath26 ) , the order parameter @xmath35 is finite and the other order parameters are zero . this means that the s-@xmath28 state is realized in the limit . the introduction of @xmath54 little decreases @xmath35 and the s-@xmath28 state still remains . at @xmath55 , the particle density is not changed , but @xmath35 suddenly vanishes and @xmath48 appears instead , where @xmath56 . this implies that the particle distributions are drastically changed , and the s-@xmath29 state is realized . further increase in the hopping induces another first - order phase transition at @xmath57 , where @xmath58 . in the case , both @xmath35 and @xmath39 are induced instead of @xmath48 and the particle number becomes away from the commensurate filling @xmath59 . therefore , we can say that the coexisting state with the @xmath28-type particle distribution ( ss-@xmath28 ) is realized . increasing the hopping @xmath54 , superfluid correlations are enhanced and the @xmath28-type particle distribution becomes obscure , as shown in fig . [ fig : low1 ] . at @xmath60 , the first - order phase transition occurs with small jumps in the order parameters and the particle density , where @xmath61 . in the large @xmath54 region , no spatial distribution appears in the particle density and the genuine superfluid state is realized .    to discuss how the @xmath28-type state competes with the @xmath29-type one , we also show the ground state energy in fig . [ fig : low2 ] .    , s-@xmath29 , ss-@xmath28 and ss-@xmath29 states in the system with @xmath47 and @xmath62 . open circles represent the first - order transition points and solid ones represent the second - order transition points in metastable states . inset shows the energy difference @xmath63 , where @xmath64s@xmath65 and ss-@xmath28 . , width=264 ]    the energy for each state is smoothly varied when the interlayer hopping is varied . we find that these curves intersect each other at two points @xmath66 and @xmath67 , where the first - order quantum phase transitions occur . if we focus on the restricted space specified by the @xmath28-type ( @xmath29-type ) particle distribution , we find the critical point between the solid and supersolid states at @xmath68 , which is shown as a solid circle in fig . [ fig : low2 ] . however , the corresponding energy is higher than the ground - state one , which means that the quantum phase transition occurs in the metastable state . therefore , we could not find critical behavior in the ground state when @xmath47 .    by performing similar calculations for different values of @xmath69 , we obtain the zero - temperature phase diagram as shown in fig . [ fig : lower - filling ] .     and @xmath46 . the solid ( dashed ) lines represent the phase boundary , where the first - order ( second - order ) transition occurs . , width=264 ]    in the weak coupling region , the superfluid state appears with a finite @xmath39 . increasing the interactions ( @xmath70 ) , translational symmetry is broken in addition to @xmath71 symmetry and the @xmath28-type supersolid state is realized . further increase of the interactions yields two solid phases ( s-@xmath28 and s-@xmath29 ) with @xmath26 . we find that the phase transition between the solid and supersolid states with the @xmath28-type particle distribution is of second - order and the others are of first - order . there are two remarkable points in the phase diagram . one of them is that the s-@xmath29 state is stable against the s-@xmath28 and ss-@xmath28 states . since the s-@xmath28 and s-@xmath29 states are degenerate in the classical limit ( @xmath21 ) , the competition should be understood by the strong coupling expansion . the energies for the s-@xmath28 and s-@xmath29 states are expressed by the second - order perturbation theory as @xmath72,\\    \frac{e _ { s-\beta}}{nv } & \sim & -\frac{1}{3}\left [   \frac { 2}{3}\big(\frac{t'}{v}\big)^2 + \frac{6}{2+r}\big(\frac{t}{v}\big)^2 \right ] . \label{perturbation - solid - energy}\end{aligned}\ ] ] this means that the @xmath28-type solid state is mainly stabilized by the intralayer hopping while the @xmath29-type solid state is by the interlayer hopping . the first - order phase transition between two solid states occurs at a certain point @xmath73 this is consistent with the phase boundary in the phase diagram although a quantitative difference appears due to a thin tube structure of the effective cluster . for the above reason , the s-@xmath29 state is realized in the strong coupling region with @xmath74 . the other point is that the ss-@xmath28 state is always stable against the ss-@xmath29 state , in contrast to the competition between the solid states discussed above . this may be explained by considering the nature of the supersolid states . the supersolid state is the coexistence between the solid state with a certain particle distribution and the superfluid state . therefore , the stability of the supersolid state with @xmath75 is roughly discussed in terms of the effective model , which is composed of localized particles and itinerant `` defect '' particles  @xcite . in this case , the defect particles can hop on layered honeycomb lattices , as shown in fig . [ fig : layered - honeycomb ] .     and @xmath29 structures . solid circles represent localized particles . bold solid and dashed lines represent the hoppings @xmath76 and @xmath54 for the defect particles . , width=264 ]    if the particle gains the kinetic energy against the intersite repulsion , particles are spontaneously doped . this induces the coexistence between the bose - einstein condensation for defect particles and the solid state for localized particles , that is , the supersolid state is realized . by using the first - order perturbation , we obtain the single particle energy for each state as @xmath77 where @xmath78 . the relation @xmath79 is always satisfied , implying that the ss-@xmath28 state is more favorable due to a one - dimensional network along the @xmath34-direction . this result is consistent with the fact that ss-@xmath29 state never appears in the ground - state phase diagram ( see fig . [ fig : lower - filling ] ) . we also find that the phase boundary between the ss-@xmath28 and s-@xmath28 states roughly obtained from @xmath80 is consistent with that in the phase diagram ( not shown ) .    in the section , we have clarified that the ss-@xmath28 state is indeed realized at the low particle density . as discussed above , the nature of the ss-@xmath28 state is essentially the same as that of the supersolid state in the single layered model . therefore , we can say that the supersolid state realized in the single layer is stable against the interlayer coupling .      here , we consider the hard - core bosonic system at the symmetric condition with @xmath81 to clarify how the supersolid state is stabilized . fig . [ fig : half1 ] shows the order parameters for the system with a fixed ratio @xmath47 . ( open circles ) , @xmath82 ( solid squares ) , and the superfluid order parameter @xmath39 ( solid circles ) as a function of @xmath50 in the system with @xmath83 and @xmath81 . , width=264 ]    in this case , the system is always half - filled ( @xmath84 ) and @xmath48 is zero . therefore , the @xmath29-type particle distribution is not realized at half filling . when @xmath85 , @xmath35 and @xmath82 are finite and @xmath86 , where @xmath56 . this implies that the s-@xmath37 state is realized . increasing the hopping @xmath54 , the second - order phase transition occurs to the ss-@xmath37 state at @xmath55 , where @xmath39 is induced . a further increase in the hopping @xmath54 increases @xmath39 and decreases @xmath35 and @xmath82 . at the transition point , @xmath82 suddenly vanishes and the phase transition occurs to the ss-@xmath28 state , where the small jumps appear in the curves of @xmath35 and @xmath39 . the first - order transition point is determined as @xmath87 , by deducing the crossing point of the energy curves for competing states ( not shown ) . the increase in the hopping @xmath54 decreases @xmath35 and increases @xmath39 . finally , @xmath35 vanishes and the second - order phase transition occurs to the superfluid state at a critical point @xmath88 , where @xmath89 .    . the solid ( dashed ) lines represent the phase boundary , where the second- ( first-)order phase transition occurs . , width=264 ]    we perform the cmf calculations for several values of @xmath69 and end up with the phase diagram at half filling , as shown in fig . [ fig : half3 ] . in the strong coupling region , the s-@xmath37 state is realized , where alternating behavior appears in the local particle density . the increase in the hoppings @xmath76 and @xmath54 enhances superfluid correlations and the quantum phase transition occurs . the ss-@xmath37 , ss-@xmath28 , and superfluid states are realized , depending on the parameters . when @xmath69 is small and @xmath50 is large , the obtained order parameters strongly depend on the tube length @xmath40 in the cmf+dmrg method , and thereby we could not determine the phase boundaries . what is the most important is that the supersolid states are widely realized between the superfluid and solid states . this is contrast to the results discussed in the single layered model where the supersolid state is not stable at half filling . to clarify how the ss-@xmath37 and ss-@xmath28 states are stabilized at half filling ,     when @xmath90 ( a ) and 0.5 ( b ) . radius of left and right semicircles at each sublattice represent a local particle density @xmath91 and a local superfluid order parameter @xmath39 , respectively . , width=302 ]    we show local particle density and superfluid order parameter in fig . [ fig : sshf - density ] . in the ss-@xmath37 state @xmath92 and @xmath93 , the superfluid order parameter is finite at each site . we also find that most of bosons are located at sublattice a in a certain layer and at sublattices b and c in a different one . this feature recalls us the supersolid state in the single layered model away from half filling . in fact , the ss-@xmath37 state may be regarded as the supersolid state with alternating layered structure , where the supersolid state with @xmath94 is realized in the even ( odd ) layers . therefore , we can say that this ss-@xmath37 state in the layered model is stabilized by the self - doping effect through the interlayer hopping . in the other word , the translational symmetry for the particle density in each layer is spontaneously broken , which stabilizes the ss-@xmath37 state at half filling . on the other hand , a different mechanism stabilizes the ss-@xmath28 state . when @xmath47 and @xmath95 , the ss-@xmath28 state is realized with three distinct sublattices in each layer , in contrast to the ss-@xmath37 state , as shown in fig . [ fig : sshf - density ] . at the sublattices a and a , the local particle density is close to half filling and the other sublattices are almost empty or fully - occupied . this enhances superfluid correlations along the chains through the sublattices a and a. in fact , the order parameter at each sublattice is maximum and its proximity effect stabilizes the superfluid state in the whole system . on the other hand , the intersite interactions @xmath96 stabilize the alternating structure in the local particle density , where b and c ( b and c ) sublattices are almost fully - occupied ( empty ) , as shown in fig . [ fig : sshf - density ] ( b ) . therefore , we can say that the ss-@xmath28 state at half filling is mainly stabilized by the one - dimensional network in @xmath34-direction .    before closing the section , we would like to comment on experiments on the optical lattice . it is known that the supersolid state in a single layered model has three - fold degeneracy . therefore , the interlayer coupling should play a crucial role in observing the supersolid state in the layered triangular lattice . in the study , we have clarified that the supersolid states are widely stable at not only the incommensurate fillings but also half filling . therefore , we can say that ultracold bosons on the layered triangular lattice is one of the most appropriate systems to observe the supersolid states . we have studied ground - state properties in the hard - core bose - hubbard model on a layered triangular lattice . combining cluster mean - field theory with the density matrix renormalization group method , we have discussed how stable the supersolid states realized in a single layer triangular lattice are against the interlayer coupling . we have obtained the rich phase diagram of the system . it has been found that two types of supersolid states are realized even at the commensurate filling . the authors would like to thank i. danshita and d. yamamoto for valuable discussions . this work was partly supported by japan society for the promotion of science grants - in - aid for scientific research grant number 25800193 and the global coe program  nanoscience and quantum physics " from the ministry of education , culture , sports , science and technology ( mext ) of japan .
Abstract:

Answer: we study ground - state properties in a hard - core bose - hubbard model on a layered triangular lattice . combining cluster mean - field theory with the density matrix renormalization group method , we discuss the effect of the interlayer coupling on the supersolid states realized in a single layered model . by examining the distributions for the particle density and superfluid order parameter , the rich phase diagram of the system is obtained . we find that the supersolid states are widely stabilized at a commensurate filling , in contrast to the case of the single layered model . the nature of the supersolid states is also addressed .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: since the establishment of the theory of general relativity by a. einstein , the study of classical gravity has been greatly extended . it provides some fascinating research subjects , especially cosmology and black holes . these two topics are powerful motivations for studying physics in curved spacetime . the quantum effects of matter fields in curved spacetime is also important for studying cosmology and black hole physics . in standard modern cosmology , our universe began with an accelerated expansion , known as inflation . the quantum effect of matter fields in the inflationary universe is believed to give rise to cosmological structures such as galaxies . however , this is merely theoretical scenario lacking any experimental evidence . another example relates to the very mysterious objects known as black holes . according to s. w. hawking s prediction  @xcite , black holes are not completely black but rather they emit black body radiation through particle creation in the curved spacetime . it is very difficult to confirm this surprising phenomenon experimentally due to the extreme low intensity of the emitted radiation .    with such a background , new attempts to discover analogous phenomena in fluid systems have appeared . in a pioneering study , unruh found that hawking radiation is not restricted to gravity fields but that it can also be observed in a perfect fluid with a sonic horizon  @xcite . the basic concept involves identifying a sound wave propagating in a flowing fluid with a scalar field in curved spacetime . in this identification , fluid flow corresponds to curved spacetime and provides the spacetime metric . the sound wave can then be considered as a field in curved spacetime . this analogy has been applied not only to hawking radiation but also to other field theoretical phenomena in curved spacetime , including cosmological particle creation  @xcite . one of the merits of considering such an analogy lies in possibility of observing quantum effects in a laboratory . for the purpose of performing a precise investigation of quantum effects , we should consider quantum fluids . bose - einstein condensates ( becs ) in trapped cold atoms represent one of the best systems in this context  @xcite , for the following three reasons : ( i ) a bec system is quite dilute and weakly interacting , so that its microscopic theory can be easily studied , ( ii ) many physical parameters of becs are experimentally controllable , and ( iii ) various quantum phenomena of becs can be directly observed , in stark contrast with other quantum fluids , such as superfluid helium systems . therefore , becs are one of the most promising quantum fluids for an analogous system .    in this spacetime analogy , spacetime corresponds to the condensate wave function described by the gross - pitaevskii ( gp ) equation  @xcite . furthermore , quanta of the scalar field in spacetime correspond to bogoliubov quasiparticles , which are low energy excitations on condensates and can be described by the bogoliubov - de gennes ( bdg ) equations . bogoliubov quasiparticles not only exist in becs at finite temperatures but they are also spontaneously created through the motion of becs , which is one of the important problems associated with defining a bec and its superfluidity in nonequilibrium bose systems . one of the main subjects of the analogy in bec is to verify the hawking radiation  @xcite . not only the hawking effect , other aspects of black hole physics as quasinormal frequencies or superradiance were also investigated  @xcite . another stream is on cosmological models representing homogeneous expanding universe  @xcite or inflationary universe  @xcite . the analogy in bec was firstly formulated by garay _ et al . _ @xcite and explored by barcel _ _  @xcite . in their formulation , the analogue geometry is read from a second - order differential equation for the perturbed phase of the condensate wave function which obeys the gp equation . in order to connect quantum physics in curved spacetime to the behavior of a realistic quantum fluid , leonhardt _ et al . _ @xcite investigated the hawking effect within the bogoliubov theory of the elementary excitations in bec . more detailed correspondence between the quanta in analogue spacetime and the bogoliubov quasiparticles was discussed by jain _ et al . _ @xcite , giving an analytical expression for particle creation spectrum in terms of the bogoliubov mode functions in the case of homogeneous bec .    in this paper , we formulate the spacetime analogy in terms of the bdg equations as an extension of @xcite , by exactly relating bogoliubov quasiparticles with quantum particles in curved spacetime . it is shown that the number operator of quanta in analogue spacetime is different from that of bogoliubov quasiparticles , unless the corresponding field is normalized correctly . with the correct normalization , the number of quanta in analogue spacetime can be exactly related with that of bogoliubov quasiparticles . furthermore , it is shown that the orthonormal relation between the normalized mode functions exactly corresponds to that of bdg wave functions . we also explicitly define bogoliubov quasiparticle creation as a phenomenon caused by a change in the complete set of bdg wave functions for the motion of a bec , and we define analogue particle creation in terms of bogoliubov quasiparticle creation . our formulation based on the bdg equations is generally applicable and gives a simple formula for particle creation spectrum in terms of the bogoliubov mode functions . in the derivation of the formula , we consider the correctly normalized field so that the formula is relevant quantitatively . our derivation of the formula is different from that given in  @xcite , and consequently the formula is generalized to be applicable in general cases . furthermore , our formulation is independent of whether condensates are homogeneous or inhomogeneous , and the obtained formula is valid even for inhomogeneous condensates . we then perform a numerical simulation for the simplest quantum effect in analogue curved spacetime : particle creation by the analogue expanding universe as an expanding bec . the spectrum obtained is consistent with the thermal maxwell - boltzmann distribution , the temperature of which has already been achieved in actual experiments . in this simulation , we set the initial state for the excitation field to be vacuum just after the condensate starts to expand , and there is no external manipulation after the initial time and the particle creation occurs spontaneously , which should be contrasted to the previous works  @xcite . in this sense , the system is isolated . the particle creation simulated here is not a result of any perturbation but a field theoretical phenomenon caused by a difference in definition of quasiparticles . this paper is organized as follows . in section [ sec : analogue ] , we derive an analogue spacetime metric from the bdg equations . in section [ sec : qft - bdg ] , we construct quantum field theory on the analogue spacetime based on bdg theory . in section [ sec : particle - creation ] , we define bogoliubov quasiparticle creation and give a formula for analogue particle creation . in section [ sec : numerical ] , we show the results of a numerical simulation giving particle creation in an expanding and shrinking bec in a trapping potential . the conclusions are given in section [ sec : summary ] . we start with the hamiltonian describing bosons trapped by an external confining potential @xmath0 given by @xmath1   \psi({\bf x}).\end{aligned}\ ] ] here @xmath2 and @xmath3 are boson field operators in the second quantized form , @xmath4 is the kinetic energy operator , @xmath5 is the particle interaction strength and @xmath6 is the chemical potential @xcite . in the heisenberg picture , the field equation is given by @xmath7 \psi({\bf x } , t).\end{aligned}\ ] ] the density operator and the current density operator are respectively defined by @xmath8.\end{aligned}\ ] ]    in a bec , the field @xmath9 consists of two components , @xmath10 , where @xmath11 describes the condensate @xcite with @xmath12 denoting a grand canonical ensemble average of @xmath13 . the bogoliubov field @xmath14 describes fluctuations about the condensate . the condensate wave function @xmath15 obeys the following gp equation @xcite , @xmath16 \psi({\bf x } , t ) \nonumber \\ & \quad + u_0 |\psi({\bf x } , t)|^2 \psi({\bf x } , t).\end{aligned}\ ] ]    the condensate wave function @xmath15 is rewritten in terms of the amplitude , @xmath17 , and the phase , @xmath18 , as follows , @xmath19 hereafter , the spatial and time dependence of fields , represented by @xmath20 , is implicit for conciseness of notation . the condensate contributions to the density and current density are @xmath21 where @xmath22 is the velocity of the condensate . with these quantities , the gp equation can be rewritten as two real equations @xmath23 the first equation is bernoulli s equation for the condensate current and the second equation represents conservation of the current . density fluctuations and current density fluctuations are described by @xmath24.\end{aligned}\ ] ] here , we take the linearized approximation , that is , we neglect higher order fluctuations , like @xmath25 . the fluctuation fields @xmath26 obey the bdg equations . these fields can be u(1 ) gauge transformed as @xmath27 the field equations are @xmath28 where @xmath29 is defined as @xmath30 performing a short calculation , one obtains @xmath31 with these fields , @xmath32 and @xmath33 take the simple forms @xmath34 where we have introduced @xmath35 with @xmath36 note that the field @xmath37 is the velocity potential in the first order approximation with respect to fluctuation fields .    taking the difference of eqs . ( [ eq : b - dg - tilde-0 ] ) and ( [ eq : b - dg - tilde ] ) , we see that the operators @xmath38 and @xmath39 obey the conservation law @xmath40 therefore , the continuity equation is satisfied within first order terms for fluctuation fields . under the hydrodynamic approximation , in which the density gradient is smooth over the local healing length @xmath41 , the sum of eqs . ( [ eq : b - dg - tilde-0 ] ) and ( [ eq : b - dg - tilde ] ) yields @xmath42 then , eq . ( [ eq : conserv-1 ] ) can be rewritten as @xmath43 where @xmath44 . now , we introduce a symmetric tensor @xmath45 where @xmath46 is a constant and @xmath47 . we can then write eq . ( [ eq : rela - org ] ) as @xmath48 = 0.\end{aligned}\ ] ] if we identify @xmath49 with a spacetime metric through @xmath50 , the field @xmath51 obeys the equation of motion for a minimally coupled scalar field in curved spacetime with the metric ( [ eq : effectivemetric - unruh ] ) , and can be seen to be a field in curved spacetime .    the equation for the field @xmath37 given in ( [ eq : rela - org ] ) was already derived as a perturbation equation of the gp equation , see for example @xcite . in this section , we have rederived it directly from time dependent bdg equations as an extension of analysis given in @xcite . this derivation is physically equivalent to the original derivation of analogue metric given in @xcite . now , the field @xmath37 is found to be a corresponding field in analogue spacetime . in this section , however , we show that the number of quanta described by @xmath37 in the analogue spacetime does not quantitatively correspond to that of quasiparticles in bec , and there is an exact expression for a quantum field in analogue spacetime which is different from @xmath37 . to address this point , we construct quantum field theory on the analogue spacetime by using the bdg theory .    for that purpose , we introduce linear transformation and separate operators and wave functions so that the operators @xmath52 and @xmath53 respectively become annihilation and creation operators , as @xmath54 ,   \quad \tilde{\phi}^{\dagger } = \sum_j [ u_j^\ast \alpha_j^{\dagger } -v_j \alpha_j ] . \label{expansion - operators}\end{aligned}\ ] ] the operators @xmath52 and @xmath53 are set to satisfy the boson commutation relation @xmath55 = \delta_{jk}$ ] . the wave functions @xmath56 and @xmath57 obey the following matrix equations : @xmath58 where @xmath59 . they then satisfy the orthonormal conditions @xmath60 the operators @xmath52 and @xmath53 describe the annihilation and creation of bogoliubov quasiparticles , respectively . here , we identify these with annihilation and creation operators of quanta in the analogue spacetime . the field @xmath51 is expanded in terms of operators by use of the expansion ( [ expansion - operators ] ) . @xmath61 ,   \label{eq : expand - phi}\end{aligned}\ ] ] where @xmath62 is defined as @xmath63 firstly , we check that the coefficient functions @xmath62 satisfies the relativistic wave equation with the metric ( [ eq : effectivemetric - unruh ] ) . ( [ eq : gp - s ] ) leads to @xmath64 with these definitions @xmath65 and @xmath66 the matrix equation ( [ eq : u_v ] ) then gives equations for the functions @xmath67 and @xmath68 as @xmath69 where the hydrodynamic approximation , @xmath70 has been used . from eq.([eq : u+v ] ) , @xmath71 substituting this into eq.([eq : u - v ] ) , we obtain @xmath72 the conservation of the condensate current ( [ eq : conserv-0 ] ) leads to @xmath73 for any function @xmath74 . therefore , we obtain @xmath75 furthermore , by noting that @xmath76 one finds that the functions @xmath77 satisfy @xmath78 = 0 ,   \label{eq : kg - f}\end{aligned}\ ] ] with metric ( [ eq : effectivemetric - unruh ] ) , as expected .    now , we calculate the klein - gordon ( kg ) inner product between functions @xmath77 in the analogue spacetime , which is a conserved product for solutions to the wave equation in curved spacetime . the kg product for functions @xmath79 and @xmath80 is defined as @xmath81 where @xmath82 is a unit normal vector orthogonal to the spacelike surface @xmath83 , that is , @xmath84 is a future - directed timelike vector . @xmath85 is an invariant volume element of the surface @xmath83 . @xmath86 is the induced metric on @xmath83 . if @xmath79 and @xmath80 are solutions of the klein - gordon equation ( [ eq : kg - f ] ) , the value @xmath87 does not depend on the choice of @xmath83 as long as @xmath83 is a cauchy surface @xcite . if we choose a @xmath88-constant surface as @xmath83 , the induced metric is @xmath89 then , we obtain @xmath90 . the normal vector to the surface is @xmath91 by using eqs . ( [ eq : byproduct - conserv ] ) and ( [ eq : u+v ] ) , the kg product between the coefficient functions can be calculated as @xmath92 . \label{eq : f - uv}\end{aligned}\ ] ] using the orthonormal relation ( [ eq : orthonormal - bdg ] ) , we obtain @xmath93 similar calculations show that @xmath94 and @xmath95 . therefore , it is seen that the coefficient functions are normal each other . at the same time , it is also seen that they are not orthonormal . we find that a set of functions @xmath96 is a set of orthonormal mode functions for the field @xmath51 . then the number operator of quanta described by @xmath51 is @xmath97 the number of quanta in analogue spacetime does not correspond to that of bogoliubov quasiparticles . in order to identify these numbers , the field has to be normalized so as to make the coefficient functions orthonormal . it is easily seen that a set of functions @xmath98 is a complete set in quantum field theory in analogue curved spacetime , that is , these functions satisfy @xmath99 then the normalized field is @xmath100   = \frac{\sqrt{m}}{\sqrt{\lambda \hbar } } \hat{\phi}^\prime . \label{eq : def - varphi}\end{aligned}\ ] ] then , the number operator of the field @xmath101 is @xmath102 we have obtained orthonormal mode functions ( [ eq : def - f ] ) and normalized scalar field ( [ eq : def - varphi ] ) in the analogue spacetime from bdg wave functions . equivalently , we have exactly related one bogoliubov quasiparticle with one quantum in the analogue spacetime . now , with this normalized field @xmath101 , quantum field theory in curved spacetime can be applicable to bogoliubov quasiparticles quantitatively . in this section , we discuss a simple example of particle creation and derive a formula connecting particle creation spectrum and bdg wave functions . let us consider a condensate which is static and having no excitation at initial time @xmath103 . for @xmath104 , the condensate evolves dynamically and , at @xmath105 , it becomes static again . then , the corresponding effective spacetime evolves dynamically and , in general , quantum field theory in curved spacetime will predict particle creation . therefore , the analogy implies that , bogoliubov quasiparticles will be created at @xmath105 . we choose a complete set at time @xmath106 ( @xmath107 ) , @xmath108 , as each @xmath109 consists of energy eigenfunctions of the bdg equations @xmath110 and @xmath111 at each time @xmath112 . in this way , we define one particle state of the quantum by using annihilation and creation operators @xmath113 and @xmath114 which diagonalize the bdg hamiltonian .    at the initial and final times @xmath106 ( @xmath115 ) , we expand the field using a complete set at each time as @xmath116,\end{aligned}\ ] ] these two sets are related to each other by linear transformations @xmath117.\end{aligned}\ ] ] therefore , the creation and annihilation operators are related as @xmath118 since the initial condensate has no excitations , the state for @xmath101 is in the initial vacuum denoted by @xmath119 , i.e. , @xmath120 . we consider the heisenberg picture and the state for @xmath101 is not time - dependent . then , the expectation value of the final number operator in the @xmath121-th mode @xmath122 is calculated as @xmath123 if @xmath124 is not zero for any @xmath125 , then particles will be created . the expected total number of particles is thus tr(@xmath126 ) . we should note that solutions to the bdg equations @xmath56 , @xmath57 solve relativistic wave equations ( [ eq : kg - f ] ) through the relation ( [ eq : def - f ] ) . therefore , if the time - dependent bdg equations are solved , then the functions @xmath127 satisfy the equation of motion in time - dependent effective spacetime . therefore , time evolution of the initial complete set in analogue spacetime is solved automatically . then , the coefficient @xmath128 is calculated as the kg product at the final time @xmath129 , @xmath130 in the derivation of this formula , we have used ( [ eq : byproduct - conserv ] ) and ( [ eq : u+v ] ) as in the calculation of ( [ eq : complete - set ] ) . here , it should be noted that , in the calculation of @xmath128 , not the set of functions @xmath77 but @xmath131 has to be considered so that the bogoliubov coefficient @xmath128 is relevant as a particle creation spectrum quantitatively . this expression for the bogoliubov coefficient ( [ eq : formula - b ] ) was obtained in @xcite for particle creation in homogeneous bec . now , the formula ( [ eq : formula - b ] ) has been shown to be valid for particle creation in dynamical condensates even when it is inhomogeneous , because , in our formulation , any simplification for the bogoliubov mode functions is not needed to be assumed . therefore , our formulation is much more powerful . it is powerful enough to show that the formula ( [ eq : formula - b ] ) is valid for general cases . in this way , the spectrum of particle creation ( [ eq : formula - spectrum ] ) can be expressed in terms of the bdg wave functions via the formula ( [ eq : formula - b ] ) .    finally , we comment that the total number of particles @xmath132 is not conserved . @xmath133 in our model , since the conservation is closed only in the gp equation for the condensate particles ; @xmath134 to overcome this unphysical result , we have to account for the quantum backreaction which is the effect of the dynamical development of the bogoliubov field to the bec . in this paper , we assume that the quantum backreaction can be neglected when the number of the bec particles far exceeds the number of bogoliubov quasiparticles . now , we show a numerical simulation giving particle creation in an expanding and shrinking bec in the trapping potential @xmath135 , where @xmath136 , @xmath137 , and @xmath138 are the trapping frequencies of the potential along the @xmath139 , @xmath140 , and @xmath141-axes respectively . we consider the case in which a bec is tightly trapped only in the @xmath139-axis ; @xmath142 , and apply the one - dimensional simulation by neglecting the spatial derivatives of the @xmath140 and @xmath141-axes . we find the stationary solution of the gp equation at the trapping frequency @xmath143 , and then change the frequency to @xmath144 at @xmath145 to make the expanding and shrinking bec . this situation can be interpreted as a cosmological expansion and contraction , because , in the expansion of condensate for example , the sound velocity becomes small and the time interval @xmath146 during which a particle on the condensate can travel over it , @xmath147 becomes long , where @xmath148 is the thomas - fermi edge of the condensate . this implies that the size of the analogue spacetime measured by the sound velocity also expands . furthermore , the sound velocity is not constant with respect to the spatial coordinate , which implies that the analogue spacetime is inhomogeneous . therefore , the analogue spacetime corresponds to an inhomogeneous expanding universe . similarly , a shrinking bec corresponds to an inhomogeneous contracting universe . although the condensate starts to move due to changing the trapping frequency , we are not interested in particle creation caused by such an external manipulation but in spontaneous one as a quantum field theoretical phenomenon . for that purpose , we set vacuum for bogoliubov quasiparticles at time just after changing the frequency , say @xmath149 . this is the initial condition in this simulation . after @xmath149 , therefore , the system is isolated without external manipulation , and particle creation can be regarded as a consequence of time evolving analogue spacetime . the spectrum of particle creation @xmath150 can be calculated as follows . firstly , we calculate the time evolution of @xmath151 and @xmath152 with that of the condensate wave function @xmath153 by using eq . ( [ eq : u_v ] ) , starting from the quasi - steady solution at @xmath154 : @xmath155 then , at @xmath129 , we obtain @xmath151 and @xmath152 after the time evolution and @xmath156 and @xmath157 by using eq . ( [ eq : stationary - u_v ] ) again . finally , we obtain the spectrum of particle creation from eq . ( [ eq : formula - b ] ) .    as the physical parameters , we use @xmath158 kg , @xmath159 hz , @xmath160 , and @xmath161 hz considering a bec experiment of @xmath162rb atoms , the total number of which is set to be @xmath163 . as the coupling constant of atomic interactions @xmath5 , we use @xmath164-wave scattering : @xmath165^{2/5},\end{aligned}\ ] ] with a scattering length of @xmath166 nm . as a numerical method , we use the fully dealiased chebyshev - galerkin method in space @xmath167 with the dirichlet boundary condition and the runge - kutta - gill method in time with @xmath168 , where @xmath169 is the characteristic length defined by the trapping potential . figure [ fig : tfedge](a ) shows the time development of the thomas - fermi edge @xmath170 of the condensate fitted by a harmonic function . dynamics of the condensate is almost periodic , and its period is @xmath171 which is consistent with that of the breathing mode @xmath172 . dependence of the thomas - fermi edge @xmath170 ( a ) and the sound velocity @xmath173 ( b ) at the center on time . here , @xmath169 , @xmath174 are the characteristic length and velocity defined by the trapping potential . ] dependence of the thomas - fermi edge @xmath170 ( a ) and the sound velocity @xmath173 ( b ) at the center on time . here , @xmath169 , @xmath174 are the characteristic length and velocity defined by the trapping potential . ] the sound velocity @xmath175 decreases with an increase in @xmath170 as shown in fig . [ fig : tfedge](b ) , also becoming periodic in time . the decrease in the sound velocity shown in fig . [ fig : tfedge](b ) corresponds to expansion of the analogue spacetime . we select the time at which the size of condensate becomes maximum as the final time @xmath176 , i.e. , @xmath177 , and calculate @xmath110 and @xmath111 ( @xmath107 ) . before calculating the spectrum of the particle creation , we need to verify whether the effect of quantum backreaction can be neglected , as discussed above . figures [ fig - noncondensate](a ) and ( b ) show the time dependence of the total number of non - condensate particles ; @xmath178 and its spectrum at @xmath129 . both figures reveal that the number of bogoliubov quasiparticles is small enough relative the number of bec @xmath179 that the quantum backreaction will not alter the result qualitatively . dependence of the number of the bogoliubov quasiparticles on time ( a ) and its spectrum ( b ) . ] dependence of the number of the bogoliubov quasiparticles on time ( a ) and its spectrum ( b ) . ]    next , we investigate its particle creation . figure [ fig : bnm ] shows the spectrum of the particle creation @xmath150 . in the low energy region @xmath180 corresponding to fig . [ fig : bnm](a ) , the spectrum has the thermal maxwell - boltzmann distribution : @xmath181 with the temperature @xmath182 . this result means that the spectrum of the particle creation in the analogue expanding universe obeys the thermal distribution with some temperature @xmath183 . it is widely believed that the thermal spectrum of particle creation is caused by formation of a static horizon . however , at the time @xmath176 , there is no sonic horizon . this particle creation is quite nontrivial .    to confirm this experimentally , a bec is needed at a temperature lower than @xmath183 , to ensure that it is not be submerged by thermal noise . using realistic parameters , @xmath183 is estimated to be 5.60 nk which has been already achieved in many bec experiments .     the spectrum of particle creation at @xmath184 in the low energy region @xmath185 ( a ) and @xmath186 ( b ) . ] the spectrum of particle creation at @xmath184 in the low energy region @xmath185 ( a ) and @xmath186 ( b ) . ]    on the other hand , in the high energy region @xmath187 shown in fig . [ fig : bnm ] ( b ) , the spectrum exhibits power - law behavior @xmath188 . this power - law behavior , however , seems to be difficult to detect because particle creation in this region is very weak . in the high energy region , the modes are not hydrodynamic and the analogy breaks down . the power - law behavior might be explained by trans - planckian physics in analogy@xcite . note that the particle creation shown in this simulation occurs spontaneously in the sense that the initial state for the excitation field has no excitation at the time just after the changing of the trapping frequency , and there is no external manipulation after the condensate starts to move . this is contrasted to the previous simulation giving homogeneous cosmological analogue @xcite , in which the external manipulation like changing strength of the atomic interaction is required . in conclusion , we have formulated quantum field theory in analogue spacetime based on the bdg equations . in this analogy , quanta in curved spacetimes are explicitly exactly related with bogoliubov quasiparticles on condensates . it has been demonstrated that , the orthonormal relations for bogoliubov modes correspond to that for quanta in effective spacetime with respect to the kg product . we derive a simple formula for the particle creation spectrum in terms of bdg wave functions , which can be applied to simple dynamical evolution whose initial and final condensations are quasi - static . furthermore , we calculate the particle creation in the analogue expanding universe by numerically solving time - dependent bdg equations for an expanding bec . the spectrum obtained is consistent with the thermal maxwell - boltzmann distribution for the temperature @xmath189 nk , which is experimentally accessible . this supports the experimental testing of particle creation in an expanding universe by using a bec .    because we neglected the effect of the quantum backreaction from the bogoliubov field to the bec , we are unable to calculate the precise amount of particle creation . to overcome this deficiency , we are extending the model to incorporate the quantum backreaction and will report on this in near future . furthermore , numerical simulations of particle creations in dynamical bec corresponding to hawking radiation will be reported in future studies .    _ acknowledgments : _ mk acknowledges jsps research fellowships for young scientists ( grant no . mt is supported in part by a grant - in - aid for scientific research from jsps ( grant no . 18340109 ) and by a grant - in - aid for scientific research on priority areas from mext ( grant no . 17071008 ) . hi is supported by a grant - in - aid for scientific research fund of the ministry of education , science and culture of japan ( grant no . 19540305 ) .
Abstract:

Answer: we construct quantum field theory in an analogue curved spacetime in bose - einstein condensates based on the bogoliubov - de gennes equations , by exactly relating quantum particles in curved spacetime with bogoliubov quasiparticle excitations in bose - einstein condensates . here , we derive a simple formula relating the two , which can be used to calculate the particle creation spectrum by solving the time - dependent bogoliubov - de gennes equations . using our formulation , we numerically investigate particle creation in an analogue expanding universe which can be expressed as bogoliubov quasiparticles in an expanding bose - einstein condensate . we obtain its spectrum , which follows the thermal maxwell - boltzmann distribution , the temperature of which is experimentally attainable . our derivation of the analogy is useful for general bose - einstein condensates and not limited to homogeneous ones , and our simulation is the first example of particle creations by solving the bogoliubov - de gennes equation in an inhomogeneous condensate .