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: positronium ( ps ) , the bound state of electron and positron , was extensively investigated since its discovery @xcite contributing to the development of bound state qed ( see e.g. @xcite for a review on the current status of this field ) . furthermore , this system provided stringent limits on possible deviation from the standard model that could indicate new physics @xcite ( see also @xcite for a comprehensive review of former experiments ) . moreover , in the field of materials science ps found various applications due to its unique properties ( see e.g. @xcite for modern reviews on this subject ) . one recent example is the characterization via positron annihilation lifetime spectroscopy ( pals ) of low - k dielectrics that are potential candidates for the next generation of integrated circuits @xcite . the motivation of the work presented in this paper is to understand if mesoporous silica with an interconnected pore network could be used for producing a high fraction of positronium at low temperatures . this would open the door for a new generation of experiments in fundamental research . cold positronium could be used to improve the precision of spectroscopic studies of ps @xcite or to perform the first spectroscopy of the ps@xmath2 molecule @xcite . furthermore , it could provide an alternative to the methods that are used for anti - hydrogen formation @xcite-@xcite . as it was suggested sometime ago @xcite-@xcite , anti - hydrogen could be formed using charge - exchange of ps with anti - protons . this was demonstrated for the charge conjugate reaction @xcite and in the atrap collaboration resonant charge exchange collisions of positrons with rydberg cs atoms were used to form rydberg ps that via charge exchange with the anti - protons produced anti - hydrogen in rydberg states @xcite . recently , two experiments @xcite were proposed to perform an anti - gravity test using this process . in both experiments , one of the main issues is that a high fraction of ps at low temperatures should be available . another interesting application would be the possibility to perform an experiment in order to confirm the interpretation of the recent dama / libra annual modulation signal @xcite as generated by mirror type dark matter @xcite . a step further would be to achieve bose - einstein condensation of ps @xcite . this would allow the exploration for the first time of the effects of the collective properties of a matter - antimatter system .    in mesoporous silica , ps is produced by injecting positrons into the film and the distribution of the implantation depth follows a makhovian profile @xcite . in the following we solely consider the long lived triplet spin state ( called ortho - positronium with 142 ns lifetime , o - ps ) because the singlet spin state ( called para - positronium , p - ps ) has a very short lifetime of 125 ps and can be considered as annihilating in the target . the o - ps ( for simplicity we will refer to it as ps ) that diffuses into the pores loses its kinetic energy via scattering . if the pores are interconnected , the ps has a probability to tunnel from one pore to another . a fraction of the ps reaches the film surface and exits into vacuum . a classical model of thermalization process was developed by nagashima et al . their calculations reproduce very well the behavior for sio@xmath2 aerogel with pore sizes of about 100 nm . however , a classical approach is not expected to give reliable predictions for ps confined in few nm pores because quantum mechanical effects become relevant . as a matter of fact , in this regime the de broglie thermal wavelength of ps is comparable with the size of the pores . recently , mariazzi et al . @xcite considered phonon scattering to reproduce the thermalization process of ps in a box ( closed porosity ) showing that the minimal energy is not that of the lowest accessible level because the momentum phonons can exchange is fixed . in the same paper , they pointed out that this is not the case in rectangular channels because in one direction the side of the potential well tends to infinity ( z - axis ) . therefore , the magnitudes of the @xmath3 and @xmath4 momentum are quantized but the @xmath5 tends to be a continuum and the minimal energy that ps can reach is given by the ground state in the x - y components . thus measuring the ps emission energy provides a method to distinguish between the two different pores architectures ( see section [ sec : discussion ] ) .    in previous studies of ps emission in vacuum using time of flight ( tof ) , many interesting effects , such as the emission from the surface of different materials @xcite-@xcite were investigated . recently , this technique was applied to study mesoporous and hybrid silica films @xcite-@xcite , to evaluate the continuity barrier @xcite and the effect of thermalization for pore surfaces decorated with different groups @xcite . however , the influence of the temperature on ps emission in vacuum from mesoporous thin films was never studied in detail . to our knowledge , only a work of mills et al . @xcite with fumed silica revealed a dependence with the temperature of the ps emission in vacuum . the fraction of ps in the low energy tail was estimated but no quantitative estimate on the value of the minimal emission energy was obtained . very recently cassidy et al . @xcite measured the emission energy of ps from mesoporous films using doppler spectroscopy . this is a different technique from that used in the present study . results consistent with the one presented in this paper ( though not sample temperature dependencies ) have been obtained . in this paper , we study the ps yield and emission energy in vacuum as a function of the target temperature for two different kinds of mesoporous thin films with the same tetraethoxysilane ( teos ) mineral source for the silica network skeleton precursor : ctacl - teos and f127-teos . the density of the c sample is approximately 1.2 g/@xmath6 and of the f sample is 1.5 g/@xmath6 . both samples were spin - coated on glass similar to the ones we measured in @xcite . the c samples are prepared via a sol - gel process using cetyl trimethyl ammonium chloride ( ctacl ) cationic surfactants as the organic pore generator ( porogen ) agent @xcite . a pure aqueous method is used . the ctacl / teos molar ratio for the films prepared is 0.22 . after deposition , the ctacl - teos / glass samples are treated at 130 @xmath7c and stored in air . the f samples use non - ionic pluronic f-127 triblock copolymer ( eo@xmath8po@xmath9eo@xmath8 ) as surfactant and were prepared in the same way as described in @xcite . both samples were calcinated for 15 minutes at 450 @xmath7c in air immediately before the e@xmath10 measurements . the recorded x - ray diffraction patterns indicate no symmetry in the pore organization . the ethz slow positron beam used for these measurements is described in greater detail in @xcite . the positrons flux is 25000 e@xmath10/s . the slow positron beam is stopped in the sio@xmath2 target . the positrons can either form positronium , i.e. o - ps or p - ps , or annihilate into @xmath11 s . the detection with an mcp ( micro - channel plate ) of the secondary electrons ( se ) emitted when the positrons hit the target serves for tagging the positronium formation . the se leave the target accelerated to 1 - 11 kev by the same voltage applied to the target relative to the grounded transport tube that is used to implant the positrons in the positronium converter . the se are then transported by a magnetic field in the backward direction , as shown in figure [ traj ] . the electrons move along the magnetic field line in spirals and are deflected to the mcp region by the @xmath12 filter . the tagging efficiency varies from 70% to 30% in the energy range 1 - 10 kev . the samples are mounted on a cryocooler head to allow the possibility of varying their temperature in the range of 50 - 400 k.      the start time t@xmath13 for the detectors is triggered by the mcp detecting the se emitted when the positrons hit the target . in both pals and tof detectors the stop is given by one ( or more ) annihilation photons depositing some energy in the calorimeter ( ecal ) . both ecal are composed of bgo crystals with hexagonal shape , 61 mm external diameter and 200 mm length . the time resolution of the system mcp - ecal was measured to be around 5 ns ( fwhm ) . the typical energy resolution of the crystals is about 25 - 30@xmath14 ( fwhm ) . the ecal for the tof is placed behind a lead slit at a distance z from the target that can be varied ( as sketched in figure [ traj ] ) . the detector is screened by four half cylinders of lead surrounding the beam pipe . the thickness of the shielding is 70 mm , the width of the slit is set to 5 mm and its position with respect to the target can be adjusted . for this measurement in order to maximize the signal to background ratio the center of the slit was placed at a distance of 18 mm from the target . the main contribution to the background is given by photons coming from direct positrons and p - ps annihilations in the target ( so called prompt peak ) . some of those photons can be detected in the crystals after compton scattering in the lead shield . the determination of the slit position with respect to the target was done by scanning the slit position in 0.1 mm steps and recording the maximum of the 511 kev annihilation peak . the scans were performed 30 min after the temperature was set on the sample . this technique was very important in order to correct the slit position as a function of the temperature . in fact , at 50 k the contraction of the cryocooler head was measured to be @xmath15 mm in agreement with the prediction of finite elements calculations performed with the comsol package @xcite . the pals detector was designed to have a large acceptance to provide a uniform efficiency for detecting the ps emitted in vacuum @xcite . in order to design the detectors and interpret the data , simulations served as a powerful tool . in the monte carlo ( mc ) simulations of our setup , the 3d eb - fields were calculated with the comsol multi - physics program and the positron / electron trajectories in the beam were simulated with geant4 . the simulation of the photon detection in the apparatus was based on the same package @xcite . new classes were written in order to simulate the o - ps production , propagation in the beam pipe , reflection on pipe walls that were assumed to be knudsen - like ( ps is reflected isotropically ) , pick - off effect or decay ( ps is not a standard particle in geant4 ) . the events for the o - ps@xmath16 process were generated taking into account the decay matrix element @xcite . the geometries of the beam transport pipe , photon detector , positron tagging system and its material were coded into simulations . the results were cross - checked with our experimental measurements for both photon detection @xcite and particle transport in the eb - fields @xcite . the measurements were taken in a clean vacuum of @xmath17 mbar . to avoid water contamination of the film during the cooling down , the target was kept at room temperature for an hour using a heater before lowering its temperature . the cooling cycles were repeated and the data confirmed the reproducibility of the results . the signals from the bgo s photomultipliers are split to record both energy and timing with a charge - to - digital converter ( qdc caen v792 ) and a time - to - digital converter ( tdc caen v775 ) . a cut on the energy deposited in the bgo between @xmath18 kev was applied to optimize the signal to background ratio suppressing compton scattering events in the collimator from direct and p - ps annihilations in the target .      in this section we present the results we obtained from the pals measurements . the spectra are analyzed using the lt9 @xcite program . the lifetimes of the decay components and its fractions are resolved by fitting the pals spectra . the program finds 3 exponentials convoluted with a 5 ns fwhm resolution function of the spectrometer . the shortest exponential ( less than 4 ns ) is originated by direct positron and p - ps annihilations . it is disregarded because we are interested only in contribution of o - ps . we define as ( @xmath19 , @xmath20 ) , and ( @xmath21 , @xmath22 ) the intensities and the lifetimes of the two longer exponentials . to determine the yield @xmath23 of ps emitted in vacuum and the lifetime @xmath24 in the pores of the film we used a model of ps escaping in vacuum @xcite . according to this model the lifetime @xmath24 in the pores of the film is defined as : @xmath25 the yield @xmath23 of ps emitted in vacuum was calculated according to : @xmath26 where @xmath27 is the escape rate of free ps into vacuum . we extracted the films thicknesses for the two samples by fitting the total ps yield as a function of the positron implantation energy to the makhovian profile . these thicknesses are reported in table [ tbl1 ] in which we also present the values of @xmath24 and @xmath28 calculated with the expressions above . we use the results of the fits of the pals spectra at 6 kev for the c sample and 10 kev for the f sample at 50 k and 300 k. these implantation energies will serve as a reference for the rest of the paper . we choose those values in order to maximize the amount of thermalized ps . at these energies the majority of the positrons are still implanted within the films ( no significant drop of the total ps yield @xmath29 is observed , see fig . [ palsdata ] ) and the emission energy is in the constant region ( see next section ) . with @xmath24 one can calculate the effective pore size @xmath30 in the films applying the gidley et al . @xcite extension of the tau - eldrup model @xcite-@xcite ( we will call it hereafter rte model ) . as one can see , at lower temperature the lifetime in the film increases as predicted by the rte model . this can be understood in the following way : the overlap of the ps wave function with the volume contained within a distance @xmath31 nm @xcite from the walls for which the annihilation rate is assumed to increase is less for ps confined in the pores occupying the ground state than for ps in excited states . since at lower temperatures the population of the ground state is higher the pick - off rate decreases . however , the measured lifetimes are lower than what is expected by the calculations using the rte model . for the c sample the discrepancy is less than 10% but for the f sample the measured value is 30% lower suggesting that the pick - off rate increases . deviations from the rte model were already observed in previous measurements using the 2 to 3@xmath32 ratio technique @xcite and in a measurement using pals @xcite . different factors could be responsible for that as pointed out in those papers and deserves further studies . since the parameter @xmath33 of the rte model was calibrated at room temperature we use the lifetimes at 300 k to determine the effective pore size that ps experiences . we report the results at the end of the paper in table [ tbl3 ] for both rectangular channels and cubic boxes . interestingly , we found that the yield of ps emitted in vacuum can be considered as independent of the temperature of the target ( see fig . [ palsdata ] ) .    ) and total yield of ps ( @xmath34 ) for f ( upper plot ) and c samples ( lower plot ) at room temperature and 50 k.__,scaledwidth=50.0% ]    this may seem to contradict the fact that the lifetime in the films increases with the temperature . however , our measurements indicate that the escape rate in vacuum is smaller at lower temperature ( see table [ tbl1 ] ) , explaining the observation that the yield is constant . a possible explanation of this effect can be found considering that ps tunnels from one pore to the other . in this case , the tunneling probability will decrease with temperature explaining our measurements . .[tbl1 ] film thickness ( z ) , ps lifetime in the pores @xmath35 and escape rate @xmath28 at 50 k and 300 k. [ cols="^,^,^,^,^,^",options="header " , ]      figure [ tof_fd ] shows the acquired tof spectra at different implantation energies for the f sample . for every run the time spectra of each crystal are calibrated finding the position of the peak arising from the annihilations in the target that defines t@xmath13 . first we extract the mean energy of ps emission in vacuum applying the analysis method proposed in @xcite as follows . the background due to the annihilations in the target is subtracted by fitting the measured time spectra with the resolution function of our detector determined using targets ( aluminum and kapton ) in which the ps formation is negligible . the tof spectra are corrected for the ps decays and the time spent in front of the detector with the factor @xmath36 . the maximum of the peak distribution defines the mean ps emission energy in the direction perpendicular to the film surface . let us note that , neglecting the reflection of ps in the beam pipe , with the time - of - flight method one measures only the mean of the energy component perpendicular to the collimator that we define as the z - axis ( we define it as @xmath37 ) . the triangles in fig . [ opsevsv ] , represent it as a function of the implantation energy . as a function of the implantation voltage at a target temperature of 300 k. the triangles represents the energy extracted with the maximum analysis method ( ma ) , the squares are after the subtraction of the non - thermalized part ( maximum analysis method after correction of the spectra , ma - cs).__,title="fig:",scaledwidth=50.0% ]   as a function of the implantation voltage at a target temperature of 300 k. the triangles represents the energy extracted with the maximum analysis method ( ma ) , the squares are after the subtraction of the non - thermalized part ( maximum analysis method after correction of the spectra , ma - cs).__,title="fig:",scaledwidth=50.0% ]    one can see that starting from 3 kev for the c sample ( this was confirmed by other measurements on similar samples @xcite ) and 4 kev for the f sample the value of the energy emitted in vacuum tends to be constant . clearly , the emission energy of ps calculated in this way is not the minimal energy due to the confinement in the pores because one has different contributions given by the convolution of the emission energy with the implantation profile . in order to isolate the thermalized part from the tof spectra , one has to subtract the contributions of non - thermalized ps components ( see fig . [ components ] ) . to estimate these non - thermalized contributions , we suppose that the shape of their distribution , nt(t ) , is represented by the tof spectrum obtained at a low implantation energy . to select the implantation energy ( from the ones we had measured ) and the scaling factor of nt(t ) for the subtraction , we relied on the mc . we used the values for which the best fit between the mc and the spectra obtained after the subtraction of nt(t ) was achieved ( for more details see @xcite ) . we found that the best fits were obtained for 2 kev with l@xmath38=200 nm in the case of the c sample , and 3 kev and l@xmath39=350 nm for the f sample . for a given implantation energy e@xmath40 , nt(t ) is scaled down by the fraction of positrons implanted at depths smaller than l@xmath38 and l@xmath39 . this was determined by using a makhovian profile . the results of this analysis are shown as the squares in fig . [ opsevsv ] . a parabolic fit is used to determine the position of the maximum . the statistical error of the fit is typically @xmath419 ns for the f and @xmath416 ns for the c sample . the uncertainty on the determination of the slit position of @xmath42 mm results in a systematic error of the order of @xmath43 mev in the determination of the ps mean emission energy @xmath37 ( for implantation energies above 3 kev ) . the subtraction procedure described above introduces a systematic error that we estimated analyzing the data using different values of l@xmath44 nm and l@xmath45 nm . the estimated error is @xmath46 mev for the c and @xmath47 mev for the f sample . thus the combined statistical and systematic error is at a level of @xmath412.9 mev for the f and @xmath413.0 mev for the c sample .    the results for the c and the f samples at room temperature and at 50 k are shown in figs . [ opsevsv_rt]-[opsevsv_lt ] . for implantation energies higher than 4 kev for the c and 5 kev for the f sample the values @xmath37 of the mean emission energy are constant . in the c sample , @xmath37 reaches its constant value at lower implantation voltages because the pore size is smaller than in the f sample ( see fig . [ opsevsv_rt ] ) . those values are higher than the thermal energy that ps will have if it would thermalize at the temperature of the film . as expected in the presence of confinement in the pores , the mean emission energy is higher for the c sample with pore sizes smaller than the f sample .     as a function of the positron implantation energy for c and f samples at 300 k.__,scaledwidth=50.0% ]     for implantation energies higher than 2 kev at 50 k and 300 k for the c sample . lower plot : ps mean energy as a function of the implantation energies higher than 4 kev at 50 k and 300 k for the f sample.__,scaledwidth=50.0% ]    the tof technique measures @xmath37 , the mean energy of the ps atoms in the z - direction . to find the mean emission energy of ps in vacuum , @xmath37 should be multiplied by a factor @xmath48 that takes into account the angular distribution . assuming that the ps is emitted mono - energetically and isotropically from the surface , one can calculate that @xmath48 is equal 2 . in this estimation , the reflection of ps in the beam pipe and the detector acceptance , i.e. the fact that a fraction of events decaying before or after the collimator aperture are detected , are not taken into account . therefore , to determine @xmath48 considering these effects we used the mc simulation we described in section [ sec : mc ] . as shown in the upper plot of fig . [ mc ] , a satisfactory agreement between the data and the mc ( adding the spectra at 2 kev that takes into account the non - thermalized ps ) is achieved . we attribute the difference between 40 and 100 ns to the approximation used in the subtraction method where the contribution of the non - thermalized ps is underestimated since only a spectrum of a defined energy is used for this correction . the fact that the ps is emitted with an angular spread is clearly supported by the data . as one can see in the lower plot of fig . [ mc ] for ps emitted with no angular spread the data are not reproduced . the physical interpretation is that in the films we studied the pores have no organization thus they are expected to be randomly aligned . the value of @xmath48 estimated with the mc is 1.7 . this is consistent with the expectation of the analytical result and the values reported in previous experiments @xcite . a detailed scan shows that the mean ps energy ( @xmath49 ) decreases with the sample temperature down to a minimum level ( see fig . [ opsenergyvst ] ) . for the c sample this value is basically constant ( 73@xmath415 mev ) in the range of temperature in which we performed our measurement . this can be understood by the fact that in this sample the confinement energy is much higher than the thermal energy at room temperature ( @xmath50 mev ) thus almost all ps is in the ground state . for the f sample there is a weak dependence on the temperature . due to the bigger pore size compared to the c sample , the energy of the ground state is only twice the thermal energy at room temperature . therefore , the probability to find the ps occupying an excited state is higher . as expected , this probability decreases with the temperature thus the minimal energy reaches its constant value of 48@xmath415 mev . the time that ps spends in the films before being emitted in vacuum was not considered in our determination of the emission energy . the measurements presented in fig . [ opsenergyvst ] are in a regime in which a classical approach is not expected to give reliable results . some theoretical work to develop a full quantum mechanical picture of the emission process is required to address this problem ( as pointed out in @xcite as well ) .      to understand the behavior of the value of the minimal energy as a function of the film temperature we present a simple model of ps in thermodynamic equilibrium at a temperature t in rectangular boxes . the expectation value @xmath51 of the hamiltonian operator for ps confined in 1 dimensional infinite well in contact with a reservoir at a temperature @xmath52 is given by : @xmath53 where z is the partition function defined as @xmath54 where @xmath30 is the dimension of the well , @xmath55 is the ps mass , @xmath56 is the principal quantum number and @xmath57 and @xmath58 are the planck and the boltzmann constants . to calculate the mean value @xmath59 of the energy for the 3d case we can use @xmath60 where to calculate @xmath61 and @xmath62 one can substitute the pore side length @xmath30 in eq . [ eq : z ] with @xmath63 and @xmath64 . for ps confined in rectangular pores , we thus obtain : @xmath65 for the case of a cubic box one can set @xmath66 .    to compare the prediction of the pore size that one can extract from the tof measurements with the pals results we fit the data using eq . [ eq : h ] . to construct the function used for the fit we kept only the first 50 terms of the sum . this is very conservative , the probability for ps to occupy a state higher than n@xmath6710 for the kind of target and the temperatures we used in this study is already negligible . the solid lines in fig . [ opsenergyvst ] represent the fit to the data where the pore side length ( @xmath30,@xmath63,@xmath64 ) are left as free parameters . the fits were repeated assuming cubic box pores and the results are shown as the dashed lines in fig . [ opsenergyvst ] . we used migrad from the minuit package @xcite as a minimization procedure . as one can see , the fit to the data suggests that the pores of both samples are better modeled as rectangular pores than cubic boxes . as proposed in @xcite , to compare the values obtained from the fit with the ones extracted from the pals measurements , one has to add twice the parameter @xmath33 ( see previous section ) . the pore side lengths obtained in this way are reported in table [ tbl3 ] . in particular , the fit supports the idea that the pores are better modeled by rectangular channels . the pore side length in one direction ( @xmath64 ) obtained from the fit is much longer than the substrate thickness . since , the same values for the side length @xmath30 and @xmath63 are obtained fitting the data assuming rectangular channels instead of rectangular boxes we do not report this value in table [ tbl3 ] as it has no physical meaning .     with the pore side lengths @xmath30,@xmath63,@xmath64 left as free parameters . the dashed lines were obtained fitting with eq . [ eq : h ] with a single side length free @xmath68 ( cubic box pores).__,scaledwidth=50.0% ]    ccc + & tof & pals + @xmath69 [ nm ] & 3.3@xmath410.1 & 4.2@xmath410.5 + @xmath70 [ nm ] & ( 2.7@xmath410.8 , 2.7@xmath410.8 ) & ( 3.3@xmath410.5 , 3.3@xmath410.5 ) + @xmath71 [ mev ] & 73@xmath415 +   + & tof & pals + @xmath69 [ nm ] & 4.1@xmath410.1 & 6.4@xmath410.5 + @xmath70 [ nm ] & ( 2.9@xmath410.8 , 3.6@xmath411.8 ) & ( 4.3@xmath410.5 , 5.3@xmath410.5 ) + @xmath71 [ mev ] & 48@xmath415 +    according to the quantum mechanical model for ps thermalization of mariazzi et al . @xcite the index @xmath56 in the sum of eq . [ eq : z ] could differ from 1 for cubic pores if the level separation of two close ps energy levels is higher than the maximum momentum that a single phonon can exchange . we performed fits with different @xmath72 but the results did not improve , supporting the idea that the pores are not very well modeled by cubic boxes . to summarize , for the c sample with the tof data the best fit was obtained for pores modeled as square channels of 2.7@xmath410.8 nm side length while the lifetime method gives 3.3@xmath410.5 nm . the best fit to the data for the f sample was obtained for rectangular channels with a cross section of 2.9@xmath410.8 nm by 3.6@xmath411.8 nm that has to be compared with 4.3@xmath410.5 by 5.3@xmath410.5 of the pals measurement . this pore size was extracted by applying the rte model to reproduce the measured lifetime of 74 ns assuming the ratio @xmath73 to be the same as in the tof results . both values obtained with the tof measurement are systematically lower than the pals results . nevertheless , considering the approximation used in our model , the assumptions that the pores can be treated as rectangular channels and the uncertainty in the determination of the pore size we conclude that our results ( summarized in table [ tbl3 ] ) are in reasonable agreement with expectations . in this paper , we show that the yield of ps emitted in vacuum measured with the pals technique is independent of the temperature of the mesoporous thin films . the lifetime in the films increases with a decrease in temperature but the measured values are lower than expected by the rte model . this suggests that another source of pick - off which depends on the temperature should be invoked . further studies are needed to investigate the origin of this effect . the escape rate in vacuum decreases explaining the observation that the yield of ps emitted in vacuum is the same at 50 and 300 k.    furthermore , we show that due to quantum mechanical confinement in the pores the ps emission energy into vacuum has a minimal value . the minimal energy is higher for the sample with smaller pores and this constant value is reached at a lower positron implantation energy . our results are in fair agreement with a model of ps confined in rectangular channels in thermal equilibrium with the sample . the measured minimal energy for the c sample was found to be 73@xmath415 mev while for the f sample it is 48@xmath415 mev . these experimental results provide a solid ground to understand how to produce ps at lower temperature and could serve to develop a more realistic model to interpret the data . this work was supported by the swiss national foundation , the eth zurich and the cnpq ( brazil ) . p.c . would like to thank c. lenz cesar , p. lotti and b. barbiellini for the very useful discussions . m. deutsch , phys . 82 , 455 ( 1951 ) . s. karshenboim , int . j. mod . phys . a 19 , 3879 ( 2004 ) . a. badertscher , p. crivelli , w. fetscher , u. gendotti , s. gninenko , v. postoev , a. rubbia , v. samoylenko , d. sillou , phys . rev . d 75 , 032004 ( 2007 ) [ arxiv : hep - ex/0609059 ] . a. rubbia , int . phys . a 19 , 3961 ( 2004 ) [ arxiv : hep - ph/0402151 ] . s.w.h . eijt , b. barbiellini , a.j . houtepen , d. vanmaekelbergh , p.e . mijnarends , a. bansil , phys . ( c ) 4 , 3883 ( 2007 ) , arxiv:0705.1113v1 . d. w. gidley , h. g. peng , r. s. vallery annual review of materials research 36 , 49 ( 2006 ) . d. w. gidley , w. e. frieze , t. l. dull , a. f. yee , e. t. ryan , and h. m. ho , phys . rev . b 60 , r5157 ( 1999 ) . d. w. gidley , w. e. frieze , t. l. dull , j. sun , a. f. yee , c. v. nguyen , and d. y. yoon , appl . 76 , 1282 ( 2000 ) . s.  chu , a.  p.  mills , j. l. hall , phys . 52 , 1689 ( 1984 ) . m.  s.  fee , a.  p.  mills , s.  chu , e.  d.  shaw , k.  danzmann , r.  j.  chichester and d.  m.  zuckerman , phys . 70 , 1397 ( 1993 ) . d. b. cassidy , h. w. k. tom and a. p. mills , aip conf . 1037 , 66 ( 2008 ) . m. amoretti , et al . [ athena collaboration ] , nature ( london ) 419 , 456 ( 2002 ) . g. gabrielse , n. bowden , p. oxley , a. speck , c. h. storry , j. n. tan , m. wessels , d. grzonka , w. oelert , g. schepers , t. sefzick , j. walz , h. pittner , t. w. hnsch , and e. a. hessels , phys . 89 , 213401 ( 2002 ) . g. b. andresen , [ alpha collaboration ] , phys . b 685 , 141 ( 2010 ) . j. w. humberston , m.  charlton , f.  m.  jacobsen and b.  i.  deutch , j. phys . b 20 , l25 ( 1987 ) . m. charlton , phys . a 143 , 143 ( 1990 ) . b. i. deutch , m. charlton , m. h. holzscheiter , p. hvelplund , l. v. jorgensen , h. knudsen , g. laricchia , j. p. merrison and m. r. poulsen , hyperfine interaction 76 , 153 ( 1993 ) . j. p. merrison , h. bluhme , j. chevallier , b. i. deutch , p. hvelplund , l. v. jorgensen , h. knudsen , and m. r. poulsen , m. charlton , phys . 78 , 2728 ( 1996 ) . c. h. storry , a. speck , d. le sage , n. guise , g. gabrielse , d. grzonka , w. oelert , g. schepers , t. sefzick , h. pittner , m. herrmann , j. walz , t. w. hnsch , d. comeau , e. a. hessels , phys . 93 , 263401 ( 2004 ) . a. kellerbauer et al . [ aegis collaboration ] , nucl . instr . and meth . b 266 , 351 ( 2008 ) . p. perez , l. liszkay , m. f. barthe , j .- boilot , m. carty , c. corbel , p. crivelli , a. curtoni , o. delferrierre , p. desgardin , m. etienne , u. gendotti , m. hassan , a. mohri , t. ohdaira , j .- m . rey , a. rubbia , h. saitoh , y. sauce , r. suzuki , a. walcarius , and y. yamazaki , aip conf . 1037 , 35 ( 2008 ) . r. bernabei et al . [ dama collaboration ] , eur . j. c 56 , 333 ( 2008 ) . [ arxiv:0804.2741 [ astro - ph ] ] . r. foot , phys . d 78 , 043529 ( 2008 ) [ arxiv:0804.4518 [ hep - ph ] ] . a.  badertscher , a.  belov , p.  crivelli , s.  n.  gninenko , j.  p.  peigneux , a.  rubbia and d.  sillou , int . j.  mod . phys .  a 19 , 3833 ( 2004 ) . p. m. platzmann and a. p. mills jr . , phys b 49 , 454 ( 1994 ) . s. valkealahti and r. nieminen , appl . phys . a 32 , 95 ( 1983 ) . m. tuomisaari , r. h. howell , and t. mcmullen , phys . b 40 , 2060 ( 1989 ) . y. nagashima et al . , phys . rev . b 58 , 12676 ( 1998 ) . a. p. mills , jr . , e. d. shaw , r. j. chichester , and d. m. zuckerman , phys . rev . b 40 , 2045 ( 1989 ) . r. s. yu , t. ohdaira , r. suzuki , k. ito , k. hirata , k. sato , y. kobayashi , and j. xu , appl . 83 , 4966 ( 2003 ) . k. ito , r. s. yu , k. sato , k. hirata , y. kobayashi , t. kurihara , m. egami , h.arao , a. nakashima , and m. komatsu , j. appl . 98 , 094307 ( 2005 ) . c. he , m. muramatsu , n. oshima , t. ohdaira , a. kinomura , and r. suzuki , phys . lett . a 355 , 73 ( 2006 ) . m. p. petkov , m. h. weber , k. g. lynn , and k. p. rodbell , appl . 79 , 3884 ( 2001 ) . r. suzuki , t. ohdaira , y. kobayashi , k. ito , r. s. yu , y. shioya , h. ichikawa , h. hosomi , k. ishikiriyama , h. shirataki , s. matsuno , and j. xu , mater . forum 445 - 446 , 224 ( 2004 ) . j. xu , a. p. mills , jr . , and c. case , appl . 87 , 54105 ( 2005 ) . h. k. m. tanaka , t. kurihara , and a. p. mills , jr . b 72 , 193408 ( 2005 ) . c. he , t. ohdaira , m. muramatsu , n. oshima , a. kinomura , r. suzuki , t. oka and y. kobayashi phys . b 75 , 195404 ( 2007 ) . d. b. cassidy , p. crivelli , t. h. hisakado , l. liszkay , v. e. meligne , p. perez , h. w. k. tom , and a. p. mills jr . a 81 , 012715 ( 2010 ) . l. liszkay , c. corbel , p. perez , p. desgardin , m .- f . barthe , t. ohdaira , r. suzuki , p. crivelli , u. gendotti , a. rubbia , m. etienne , and a. walcarius , appl . lett . , 92 , 063114 ( 2008 ) . y. cohen , k. landskron , n. tetretault , s. fournier - bidoz , b. hatton , g. a. ozin , adv . 15 , 593 ( 2005 ) . w. yantasee , y. lin , x. li , g. e. fryxell , t. s. zemanian and v. v. viswanathan , analyst , 128 , 899 ( 2003 ) . n. alberola , t. anthonioz , a. badertscher , c. bas , a.s . belov , p. crivelli , s.n . gninenko , n.a . golubev , m.m . kirsanov , a. rubbia and d. sillou , nucl . instr . method a 560 , 224 ( 2006 ) . p. crivelli , u. gendotti , l. liszkay , p. perez , a. rubbia , mat . forum 607 , 251 ( 2009 ) . agostinelli et al [ geant4 collaboration ] , nucl . instr . method a 506 , 250 - 303 ( 2003 ) . a. ore and j. l. powell , phys . 75 , 1696 ( 1949 ) . a. badertscher , p. crivelli , m. felcini , s.n . gninenko , n.a . goloubev , p. nedelec , j.p . peigneux , v.postoev , a. rubbia and d. sillou phys . p. crivelli ,  a new search for invisible decays of ortho - positronium " , phd thesis , no.16617 , eth zurich , switzerland ( 2006 ) . l. liszkay , c. corbel , l. raboin , j .- boilot , p. perez , a. brunet - bruneau , p. crivelli , u. gendotti , a. rubbia , t. ohdaira and r. suzuki , appl . 95 , 124103 ( 2009 ) . j. kansy , nucl . instrum . meth . a 374 , 235 ( 1996 ) . l. dull , w. e. frieze , d. w. gidley , j. sun , and a. f. yee , j. phys . b 105 , 4657 ( 2001 ) . s. j. tau , chem . 56 , 5499 , ( 1972 ) . m. eldrup , d. lightbody amd j. n. sherwood , chem . phys . lett . 63 , 51 , ( 1981 ) . t. goworek , k. ciesielski , b. jasinska , and j. wawryszczuk , chem . 272 , 91 ( 1997 ) . k. ito , h. nakanishi and y. ujihira j. phys . b 103 , 4555 - 4558 ( 1999 ) . s. mariazzi , l. toniutti , n. patel , r.s . brusa , applied surface science 255 , 191 - 193 ( 2008 ) . c. g. fischer , m. h. weber , c. l. wang , s. p. mcneil , and k. g.lynn , phys . rev . b 71 , 180102r ( 2005 ) . s. thraenert , e. m. hassan , d. enke , d. fuerst and r. krause - rehberg , phys . status solidi c 4 , 3819 ( 2007 ) . u. gendotti , ph.d thesis unpublished , eth zurich , switzerland ( 2010 ) . f. james , m. roos , minuit , function minimization and error analysis , cern d506 ( long writeup ) .
Abstract:

Answer: in this paper , we present measurements of the ortho - positronium emission energy in vacuum from mesoporous films using the time of flight technique . we show evidence of quantum mechanical confinement in the mesopores that defines the minimal energy of the emitted ps . two samples with different effective pore sizes , measured with positron annihilation lifetime spectroscopy , are compared for the data collected in the temperature range 50 - 400 k. the sample with smaller pore size exhibits a higher minimal energy ( @xmath05 mev ) , compared to the sample with bigger pores ( @xmath15 mev ) , due to the stronger confinement . the dependence of the emission energy with the temperature of the target is modeled as ortho - positronium being confined in rectangular boxes in thermodynamic equilibrium with the sample . we also measured that the yield of positronium emitted in vacuum is not affected by the temperature of the target .


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: it is widely believed that quarks , the basic building blocks of strongly interacting matter , lose their confinement at some large density or temperature , at or near the chiral phase transition . many models , for example the instanton liquid model , suggest that the unconfined quarks exert an attractive force on each other . the bcs mechanism predicts that an arbitrarily weak attractive force leads to an instability of the fermi - surface , giving rise to a colour superconductor @xcite . this new phase of quark matter has been the subject of many papers in recent years ; see ref . @xcite for a detailed list of references . the interplay of the order parameter with the flavour symmetry ( either su(2 ) or su(3 ) , since only the lightest quarks are relevant ) can generate many interesting phases . these include the two - flavour colour superconductor ( 2sc ) which pairs quarks of different flavour and colour and leaves one colour unpaired @xcite ; the colour flavour locked phase ( cfl ) which is possible in the three - flavour case @xcite ; and the crystalline colour superconducting phase @xcite . in general , a robust colour superconducting gap of the order of @xmath0 has been predicted by calculations which use different interactions and different superconducting states . the density at which the maximum pairing occurs is several times nuclear density . an important question is what happens in finite systems . if we ever succeed in creating such a phase in the lab , it must occur in finite , colour - singlet , states . equally , if we have a nugget of quark matter at the centre of a neutron star , we expect it to be colour neutral , since it originates from the compression of a colour neutral object and is surrounded by colour neutral matter . the techniques with which to discuss such phenomena have long been used in nuclear physics , where the superfluid nature of finite nuclei should not be described by the number - non - conserving bcs state , but by a projected one @xcite . the additional complication here is the requirement of colour neutrality , since the cooper pairs in 2sc superconductivity are colour anti - triplets . hence an additional projection onto a colour - singlet state is required . in contrast the cooper pairs in nuclear physics already carry zero angular momentum .    in this work we will use the model of ref . @xcite and work at zero temperature . we limit our study to the simplest state , the 2sc . we model the effects of finite size by confining the quarks to a cubic box with antiperiodic boundary conditions . as a result boundary effects are not included . after projection the exact equations for the superconducting gap are very involved , so in this work we confine ourselves to a one - parameter variation . with these limitations , we find that the superconducting gap persists after , though not always before , projection in systems of finite size . finite size effects in superconducting matter have also been considered in ref . @xcite , where it was found that the stability of small lumps of cfl matter was enhanced . in that paper however the emphasis was on thermodynamic effects ; shell structure and projection onto states of good quantum numbers were not considered . the paper is organised as follows : in section [ sec : model ] we describe the model and fix the notation ; in section [ sec : numb_proj ] we apply the projection over baryon number to the 2sc and obtain an expression for the expectation value of the energy in this projected state ; in section [ sec : colour_proj ] we apply the projection over the colour singlet to the 2sc and generalise the expression given in the previous section . in section [ sec : res ] we give the numerical results , and finally , in section [ sec : concl ] we state our conclusions . the reader will find technical details and some useful formulae in the appendices [ app:1 ] and [ app:2 ] . we wish to investigate the effects of finite size , definite baryon number and colour neutrality on the colour superconducting state . as mentioned in the introduction , a variety of models have been so far used to study the emergence of colour superconductivity . we follow ref . @xcite and use a standard njl model with massless quarks carrying two flavours and three colours . we leave for future work the study of more complex scenarios .    following ref . @xcite we model the attractive interaction responsible for the pairing of quarks at large densities with a contact four - quark interaction whose form is suggested by the instanton - induced interaction . the hamiltonian then has the form @xmath1 with    @xmath2    where flavour and colour indices are latin and greek respectively , and the tensor @xmath3 is given by @xmath4    at zero density the quarks acquire a dynamical mass , through the generation of a chiral condensate . this mass drops as the density increases , reaching zero at the chiral phase transition . at higher densities the colour condensate which signals the onset of colour superconductivity forms .    in order to regularise the contributions stemming from large momenta , the authors of ref . @xcite have supplemented the interaction ( [ eq : hint1 ] ) with a form factor , to be associated with each quark of a given momentum . we choose to use a form factor of the form @xmath5    the parameter @xmath6 determines the sharpness of the cut - off , and has been taken to be 10 . this leaves two free parameters , @xmath7 and the coupling constant @xmath8 in the interaction . these are fixed at zero density by requiring the pion decay constant to have its physical value , and by choosing the dynamically generated quark mass to be 400  mev ; this gives @xmath9  mev and @xmath10  mev@xmath11 . the modes of the fields are quantised in the standard way , and we define the creation operators @xmath12 , @xmath13 , @xmath14 for particles , antiparticles and holes relative to the fermi sea , @xmath15 , which can thus be written as @xmath16 for all three colours are equal . in fact the properties of the bcs state are independent of the initial fermi momentum for the condensing quarks . however colour projection becomes involved if @xmath17 is different for the spectator quarks . the fermi sea is the natural choice of ground state at low to intermediate density . at large densities the interaction ( [ eq : hint1 ] ) is responsible for the condensation of cooper pairs of quarks , and we get a new form for the wave function , @xmath18 where @xmath19 \nonumber \\ \prod_{\absum\kvec}\ ! & \left [ \cos\theta_b^r(k ) +   \epsilon_{\alpha \beta 3 } \sin\theta_b^r(k )   e^{i \xi_b^r(k ) } \hat{b}_{r 1\alpha}^\dagger(\kvec ) \hat{b}_{r2\beta}^\dagger(-\kvec ) \right ] \nonumber \\ \prod_{\absum\kvec}\ ! & \left [ \cos\theta_c^r(k ) +   \epsilon_{\alpha \beta 3 } \sin\theta_c^r(k )   e^{i \xi_c^r(k ) } \hat{c}_{r 1\alpha}^\dagger(\kvec ) \hat{c}_{r2\beta}^\dagger(-\kvec ) \right ] \nonumber \\ & \equiv\hat{g}_{l a } \,\hat{g}_{l b } \,\hat{g}_{l c } \label{eq : gl}\end{aligned}\ ] ] and similarly for @xmath20 , with the substitutions @xmath21 . the sum denoted @xmath22 is taken over @xmath23 and @xmath24 only ; the bcs angles @xmath25 and @xmath26 are to be determined variationally .    the bcs wave function @xmath27 of eq . ( [ eq : wf1 ] ) describes the condensation of pairs of quarks with the same helicity , but antisymmetric in flavour and in colour , as can be seen by inspection of eq . ( [ eq : gl ] ) . the condensate carries a colour-@xmath28 order parameter . only two colours , by convention @xmath29 and @xmath30 ( or red and green ) , condense , leaving the third colour inactive . this corresponds to a spontaneous symmetry breaking . the bcs angles are obtained by minimizing the thermodynamic potential @xmath31 the phases @xmath32 are fixed by requiring that the expectation value of @xmath33 is maximally attractive ; this gives @xmath34 minimising with respect to the functions @xmath35 one finds the gap @xmath36 and the angles @xmath37 for particles , antiparticles and holes respectively to be given by the solutions of the coupled equations @xmath38 and @xmath39 \ . \label{eq : gap}\end{aligned}\ ] ]    . the angles are in units of @xmath40 . the solid , dotted and dashed lines correspond to particle , antiparticle and hole contributions respectively.,width=264 ]    in fig . [ fig : paperfig2 ] we display the solutions for the bcs angles @xmath41 , @xmath42 and @xmath43 in the infinite volume for a density corresponding to a chemical potential @xmath44 . the solid , dotted and dashed lines correspond to the particle , antiparticle and hole contributions respectively . the solutions quickly fall off at momenta larger than the cutoff @xmath7 as a result of the action of the form factor . pairing is most effective at the fermi surface , @xmath45 . the bcs wave function ( [ eq : wf1 ] ) does not have a definite baryon number ( the condensate is a superposition of states with different number of pairs ) and it is not in a colour singlet state . however , the number of particles is fixed on average by the lagrange multiplier @xmath46 in eq . ( [ eq : gibbs ] ) . because the fluctuations in the baryon number behave as @xmath47 , they become negligible as a fraction of the baryon number when the condensate contains a large number of particles .        in fig . [ fig : paperfig3 ] we show the quark density as a function of the chemical potential for a fermi gas of quarks ( solid line ) and for the two - flavour bcs condensate ( dashed line ) . at small density the condensate contains more particle pairs than antiparticle or hole pairs ; as the density grows , the reverse is true , since the particle states start to be suppressed from the form factor and more ( fermi ) hole states become available .    in this paper we are interested in the effect of finite size on the colour superconductor . to do this we consider the quarks to be in a cubic box with antiperiodic boundary conditions . we start with a fermi sea in which all levels with @xmath48 are occupied , and consider the projected and unprojected bcs states built on this . for simplicity we only consider filled shells .     for a box of size @xmath49 . dashed line : expectation value of the number of red and green particles in the bcs condensate as a function of @xmath46 . at the intersections between the two curves the number of pairs of particles in the bcs state equals on average the number of pairs of antiparticles and holes . , width=264 ]    in fig . [ fig : paperfig6 ] we have plotted the quark density as a function of @xmath46 for the two active colours in the fermi gas ( continuous line ) and in the bcs condensate ( dashed line ) , in a finite box of size @xmath50 . notice that the dashed curve does not depend on @xmath17 ; in fact only quantities involving the third ( inactive ) colour will depend on @xmath17 . because we work with equal values of @xmath17 ( that is the same filled shells ) for all three colours , at the intersection between the two curves the box will contain on average the same number of red , green and blue particles . these are the points at which we will project out colour singlets .    in the following sections we study the occurrence of colour superconductivity in systems of finite size and with a fixed number of particles . we first project the bcs wave function onto states of definite baryon number ( section [ sec : numb_proj ] ) and then onto colour - singlet states ( section [ sec : colour_proj ] ) . there are various ways to solve the problem of picking a component of fixed particle number from a bcs state . the most efficient ones are based on the calculation of a contour integral ; we follow here the method of residues of ref . this will be used to project the colour superconducting state considered in eq . ( [ eq : wf1 ] ) onto a state with a fixed number of particles . we start by defining the rotated bcs state @xmath51 where @xmath52 is the number operator , @xmath53 is the number of particles in the fermi sea and @xmath54 is a complex number on the unit circle with phase @xmath55 . the rotated state can be created directly from the fermi sea by the rotated operators @xmath56 : @xmath57 where    @xmath58 \nonumber \\ & \times & \prod_{\absum\kvec } \left [ \cos\theta_b^r(k ) +   \zeta^ *   \epsilon_{\alpha \beta 3 }   \sin\theta_b^r(k )   e^{i \xi_b^r(k ) }   \hat{b}_{r 1\alpha}^\dagger(\kvec )   \hat{b}_{r2\beta}^\dagger(-\kvec ) \right ] \nonumber \\ & \times & \prod_{\absum\kvec }   \left [ \cos\theta_c^r(k ) +   \zeta^ *   \epsilon_{\alpha \beta 3 }   \sin\theta_c^r(k )   e^{i \xi_c^r(k ) }   \hat{c}_{r 1\alpha}^\dagger(\kvec )   \hat{c}_{r2\beta}^\dagger(-\kvec ) \right ] \nonumber \\ & \equiv &   \prod_{\absum\kvec } \hat{g}_{l a \alpha \beta}(\zeta,\kvec )    \prod_{\absum\kvec } \hat{g}_{l b \alpha \beta}(\zeta,\kvec )    \prod_{\absum\kvec } \hat{g}_{l c \alpha \beta}(\zeta,\kvec ) ,    \label{eq : glz}\end{aligned}\ ] ]    and similarly for @xmath59 . because the creation of antiparticles or holes lowers the baryon number , the corresponding terms in @xmath60 contain @xmath61 or , equivalently , @xmath62 . the number - projected wave function is produced from an appropriate superposition of rotated states : @xmath63 where @xmath64 is the number of pairs . the only operators whose expectation values in the number - projected state are non - vanishing are those which are particle - number conserving and so commute with the number operator @xmath52 . this enables us to write @xmath65    a convenient way to calculate these matrix elements is to find the thouless operator @xmath66 which maps @xmath67 into @xmath68 @xcite : @xmath69 the explicit expression for @xmath70 is given in appendix [ app:1 ] . there we also give the definitions of the quasi - particle operators @xmath71 which annihilate the state @xmath72 .    the thouless operator can be expanded in a series of terms containing zero , two , four ,  quasiparticle creation operators : @xmath73 depending on the form of the interaction , only a limited number of these terms will contribute to the matrix elements of the hamiltonian . luckily , the interaction ( [ eq : hint1 ] ) selects only the first two terms in the expansion . for left - handed particles these are @xmath74    the function @xmath75 can be expanded as a taylor - laurent series about @xmath76 , with real coefficients @xmath77 . these are related to the normalisation constant @xmath78 in eq . ( [ eq : wf2 ] ) , by @xmath79 , and to the expansion of the bcs state @xmath80 in terms of states of definite particle number : @xmath81 since @xmath80 has unit norm , the sum of the @xmath77 is unity .    from eq . ( [ eq : expect ] ) we obtain the expectation value of the hamiltonian in the number projected state : @xmath82    explicitly , we obtain for the non - interacting part of the hamiltonian : @xmath83 \right.\nonumber\\ & + 2 \,\theta(k - k_f ) \sin^2 \!\theta_a(k)\ , { \cal i}_{a , n}(\theta_a(k ) )   \nonumber \\ & + \left . 2\sin^2 \!\theta_b(k ) \,{\cal i}_{b , n}(\theta_b(k ) )   \right\ } \ . \label{eq : h02}\end{aligned}\ ] ] in this expression the first term describes the hole contribution , the second the particle one , and the last the antiparticle . the first term ( @xmath77 ) in the square brackets describes the spectator colour ; the second term reduces to twice this in the absence of pairing ( when @xmath84 ) . the interaction hamiltonian leads to the more complicated form    @xmath85    the volume of the box is @xmath86 . the contour integrals are included in the definitions of @xmath87 and @xmath88 ( see appendix a ) ; @xmath89 is the angle between @xmath90 and @xmath91 . the form factors in the expression ( [ eq : hint2 ] ) suppress the contributions coming from large momenta . the expectation value of the number operator is obtained from eq . ( [ eq : h02 ] ) by dropping the single - particle energy @xmath92 which appears immediately after the summation sign and swapping the sign of the antiparticle ( @xmath93 ) term . these expressions ( [ eq : h02 ] ) and ( [ eq : hint2 ] ) reduce to the expectation values in the unprojected bcs state if the integrals @xmath94 and @xmath95 and the constant @xmath77 are set to one . the infinite - volume gap equations ( [ eq : solu ] ) and ( [ eq : gap ] ) are easily obtained by minimising @xmath96 with respect to the angles @xmath35 , if the sums are replaced by integrals . however in the finite volume projected case , the equations are much more complicated . for a start , in the cubic box there is no reason to expect the angles to depend only on the length and not on the direction of @xmath97 . this introduces a dependence on the specific geometry that one might be quite happy to ignore . however even then the dependence of the @xmath94 s and @xmath95 s on the @xmath35 s renders the equations very complicated . we do not attempt to solve these . instead we assume that the form of the @xmath35s is given as before by eq . ( [ eq : solu ] ) , and minimise the expression for @xmath98 obtained from eqs  ( [ eq : h02 ] ) and ( [ eq : hint2 ] ) with respect to the gap @xmath36 only .    this process of restricting the bcs angles to have the same form as in the unprojected state will presumably introduce fewest errors where the projected state has a high overlap with the unprojected state . thus although in principle we can project a state of a given baryon number from any bcs state , in practice we choose values of @xmath99 where the expectation value of the number of pairs in the unprojected bcs state has the value desired in the projected state ( see fig .  [ fig : paperfig6 ] ) .    in the bcs state there are non - vanishing expectation values for diquark condensates such as @xmath100 . however as this operator changes the quark number by two units , it does not commute with the number operator @xmath52 and the expectation value in the number - projected state will vanish . none - the - less , the square of the diquark condensate operators do not vanish , and they act as order parameters in this case . we now turn to the issue of colour projection . as we have pointed out before , the colour - superconducting state which forms at large densities must be a colour singlet . however , the state in eq . ( [ eq : wf1 ] ) does not fulfill this requirement , given the form of the operators @xmath101 . for this reason we want to project the colour - singlet state from eq . ( [ eq : wf1 ] ) . this requires us to integrate over the group manifold : @xmath102 is the volume element on the group manifold and @xmath103 is the number projected bcs state . @xmath104 is a unitary operator which performs a rotation in the @xmath105 space . following @xcite we parametrize an element @xmath106 of su(3 ) in the form @xmath107 where @xmath108 are the gell - mann matrices , and write the volume element as @xmath109 to obtain the realization of these rotation operators when acting on our quarks states , we replace the @xmath110 s with the transition operators @xmath111 , which are associated with the su(3 ) group . we can simplify the algebra considerably by noting that @xmath113 \ , \end{aligned}\ ] ] where @xmath29 , @xmath30 and @xmath114 are the values taken by the the colour index . since we are working with states built on fermi seas which have equal numbers of red , green and blue quarks , we see that @xmath115 annihilates only the number - projected bcs state @xmath116 which contains zero net pairs ( _ i.e. _  the number of pairs of holes and antiparticles is the same as the number of pairs of particles ) . we therefore restrict our projection to this state . thus the only operator which acts non - trivially is @xmath117 , which is non - diagonal in the colour space and connects one of the two paired colours with the remaining one . as a result we are able to reduce the integration over the eight different angles @xmath118 in eq . ( [ eq : rotation ] ) down to an integration over the single angle @xmath55 for the rotation generated by @xmath117 . the residual volume element is therefore @xmath119 where @xmath120 is the angle associated with @xmath117 .    as before , we want to project states with zero pairs from bcs states which have zero pairs on average , so we are restricted to working only at certain values of @xmath46 . these are the values at which the two curves of fig . [ fig : paperfig6 ] cross . for small boxes , this is a significant restriction on the densities we can consider .    just as in the case of number projection particle - number - conserving operators have non - zero expectation values , so when we project onto a colour singlet only colour - singlet operators will have non - vanishing expectation values . such operators commute with the colour rotation , and so we can write ( cf . ( [ eq : expect ] ) ) @xmath121        ( recall that the first numerical index on the creation operators is the flavour index , and the second is colour ) . this clearly reduce to the original form eq . ( [ eq : glz ] ) for @xmath126 . by looking at the expressions ( [ eq : rot2 ] ) we see that , by acting with @xmath127 on the fermi sea , one can create pairs of particles in the colours @xmath128 or in the colours @xmath129 . the probability of each of these events will depend on the magnitude of the angle @xmath55 ; clearly , the two extrema of the interval ( i.e. @xmath126 and @xmath130 ) correspond to pure pairing in the @xmath128 and @xmath129 colours respectively . once again , calculations are most conveniently done with the aid of a thouless operator , @xmath131 , which fulfills the relation @xmath132 the expectation value ( [ eq : expct_vl ] ) therefore reads @xmath133 the explicit form of the operator @xmath131 is given in appendix [ app:2 ] . we can expand @xmath134 as a laurent series in @xmath54 , with coefficients @xmath135 . the normalization constants @xmath136 are related to integrated coefficients @xmath137 as in the pure number - projection case the @xmath78 are to the laurent coefficients @xmath77 of @xmath138 . the forms of the expectation values of @xmath139 , @xmath140 and @xmath52 are as given in the previous section , in eqs  ( [ eq : h02 ] ) and ( [ eq : hint2 ] ) , but setting @xmath141 , and with the replacement of @xmath77 with @xmath142 and of the integrals @xmath94 and @xmath95 with the new integrals @xmath143 and @xmath144 which are given in appendix [ app:2 ] . once again a one - parameter variation with respect to the gap @xmath36 is performed .        we can now give the explicit expression for the thouless operator @xmath131 , which can be derived in a similar fashion to @xmath70 . following the same approach as before we first look for an operator , @xmath161 , that is at most bilinear in the quark operators and that it can be applied to @xmath162 to give the colour rotated operator @xmath163 . we will simply state the results for the component corresponding to left handed particles : @xmath169 where @xmath170 and @xmath171 the last term pairs the third colour with one of the other two . this term does not to contribute to the expectation value of the hamiltonian .      in analogy with sec . [ app1:sub3 ] we now define the following objects : @xmath172 where @xmath173 is given by eq . ( [ eq : red ] ) . these reduce to the integrals @xmath95 of eq . ( [ eq : app1j ] ) if @xmath55 is set to zero and the colour - group integral is suppressed . this approximation is particularly useful in the numerical calculation , since it allows us to calculate the integrals analytically . with this simplification , the integrals in eq . ( [ app2:j ] ) can be written in a closed form as : @xmath179 , \nonumber \\ j_{b , n}(\theta_1,\theta_2 ) & = &   \sum_m   d_m(0 )   \left [ i_1    \lambda_{n , m}(\theta^{\rm c}_1,\theta^{\rm c}_2 ) \right . \nonumber \\ & + & \left . \frac{i_2}{2 }    \left ( \sin^2 \!\theta_1   \sin^2 \!\theta_2   \delta_{n , m}(\theta^{\rm c}_1,\theta^{\rm c}_2 ) - \cos^2\!\theta_1   \cos^2\!\theta_2   \delta_{n , m+2}(\theta^{\rm c}_1,\theta^{\rm c}_2)\right ) \right ] , \nonumber \\ j_{c , n}(\theta_1,\theta_2 ) & = & \sum_m   d_m(0 )   \left [ i_1    \lambda_{n , m}(\theta^{\rm c}_1,\theta_2 ) \right . \nonumber \\ & + & \left . \frac{i_2}{2 }    \left ( \cos^2\!\theta_1   \sin^2 \!\theta_2   \delta_{n , m+2}(\theta^{\rm c}_1,\theta_2 ) +   \sin^2 \!\theta_1   \cos^2\!\theta_2   \delta_{n , m}(\theta^{\rm c}_1,\theta_2 ) \right . \nonumber \\ & + & \left . 2   \sin^2 \!\theta_1   \sin^2 \!\theta_2 \ \delta_{n , m+1}(\theta^{\rm c}_1,\theta_2 ) \right ) \right ]   , \end{aligned}\ ] ] where @xmath180 and we have defined @xmath181 and @xmath182 the parameter @xmath177 is defined in eq . ( [ eq : alpha ] ) . the integrals of eq . ( [ eq : lam ] ) can be also evaluated analytically . for a box of size @xmath49 and for different densities corresponding to the intersections between the two curves in the previous figure . the number of particles corresponding to each curve is indicated . right plot : distribution of the rescaled laurent coefficients @xmath145 , as functions of the rescaled variable @xmath146.,width=302 ]        first however we look at the effect of finite size on the unprojected bcs state . as the size of the box grows the fractional fluctuations in the baryon number will drop as @xmath147 , and we expect that the value of the gap and of the energy should converge toward their values in the infinite volume . this behaviour is evident in fig . [ fig : fsize1 ] , which shows the superconducting gap and total energy per particle . convergence is obtained quickly and continuously . now we turn to the effects of number and colour projection . as we have explained in section [ sec : colour_proj ] , the projection over colour is carried out on bcs states having the same total baryon number as the underlying fermi gas they are built on . this restricts us to certain values of @xmath46 , as we saw from fig . [ fig : paperfig6 ] . in eq . ( [ eq : thouless1 ] ) we introduced the function @xmath138 whose laurent coefficients @xmath77 are related to the amount of the projected state with @xmath64 pairs contained in the unprojected state ( see eq . ( [ eq : comp ] ) ) .    fig . [ fig : paperfig7 ] shows the distribution of the laurent coefficients for the function @xmath138 at these points.because we are working at our chosen values of @xmath46 , the coefficients peak at @xmath141 , and we see that the unprojected bcs state has a substantial overlap with the projected state with zero pairs . the distributions are easily approximated with gaussian functions .    in fig . [ fig : gap_l3 ] we show the superconducting gap as a function of chemical potential for a range of box sizes . the dashed line is the finite - volume unprojected gap , and for comparison the continuous line shows the result for the infinite - volume case . we remind the reader that because we perform only a limited variation after projection , as explained above , the projection can only be done for certain discrete values of @xmath46 for which the average net number of pairs in the unprojected bcs state is zero . at these points we show results for the unprojected bcs ( circles , which must coincide with the dashed line ) , and for the number - projected ( triangle up ) and combined number - colour projected bcs ( triangle down ) . the most striking finite - volume effect is the disappearance of the gap for values of @xmath46 which are far from shells . these regions are confined to very low @xmath46 for boxes above 3  fm in length . in these regions the effects of projection on the gap can be striking , as can be seen even at 7  fm . however in all other regions , where the unprojected gap does not disappear , projection has very limited effect . in general the results for colour and number projection lie above those for number projection alone ; the latter show almost no difference from the unprojected case unless the unprojected gap vanishes . the effect of finite size on the energy of the condensate has already been shown in fig . [ fig : fsize1 ] . the effect of projection is to lower the energy by less than 1  mev even for the smallest boxes considered . in the approach to the infinite - volume limit , shell effects for the unpaired quarks vanish more slowly than finite - size or projection effects on the condensate . in this work we have studied the effects of finite size over the colour superconducting state , by considering quark matter at finite density and zero temperature enclosed in a cubic box . we have considered the 2sc superconducting state of ref . @xcite and we have calculated the effects of the projection over definite baryon number and over a colour singlet state . the most significant finite size effect is the vanishing of the gap for values of the chemical potential between widely spaced shells . for these cases , projection has a significant effect , restoring a sizeable gap . otherwise the effects of projection are slight , though in all cases each projection increases the gap to some degree . this behaviour points to the interesting possibility that the interplay between the chiral and the superconducting phase can also change in a finite system . this hypothesis can be verified by extending the present calculation , allowing a competition between the two phases . this topic is currently under investigation @xcite . there are other aspects which require future investigation as well : it would be interesting to extend our results to more general form of pairing interactions , rather than the simple `` instanton - motivated '' interaction of eq.([eq : hint1 ] ) . for technical reasons , the number projection which was carried out exactly in the present work , would require more effort in a general case , where the interaction contributions will contain not only terms which mix the helicity , such as @xmath148 , but also terms which connect fermions with the same helicity , such as @xmath149 . it would be more appropriate in such cases to apply approximate number projection techniques , which have been developed in the contest of nuclear theory @xcite . another aspect of interest , which we plan to investigate in the future , is the effect of the colour projection on different superconducting states , and in particular on the colour - flavour locked phase ( cfl ) , which was originally studied in @xcite . because in this state the colour and flavour degrees of freedom are `` locked '' , the colour projection will affect both colour and flavour , and could possibly yield larger effects than the one observed in the case of the 2sc . we also remark that in the present analysis we have not performed a full variation after projection . instead we have followed a simpler approach , assuming for the solution the same functional form as in the unprojected case and obtaining the numerical solution by performing a minimization of the total energy with respect to the gap @xmath36 . this work was supported by the uk epsrc . the authors would like to thank dr . m. alford for stimulating discussions and suggestions . p.a . would like to thank prof . walecka for useful discussions on the topic . notice that these operators fulfill canonical anticommutation relations . the inverse equations , which allow us to express the original operators in terms of the quasi - particle ones , are also useful and read : @xmath151      we now consider the operator defined in eq .  ( [ eq : thouless1 ] ) . for the part corresponding to @xmath152 we look for solutions satisfying @xmath153 \hat{g}_{l a \alpha \beta}^\dagger(1,\kvec ) & = &   \hat{g}_{l a \alpha \beta}^\dagger(\zeta,\kvec ) , \nonumber \\ \left [ w_1^b(\zeta ) + \epsilon_{\alpha \beta 3 }   w_2^b(\zeta )    \hat{b}^\dagger_{r 1 \alpha}(\kvec )   \hat{b}^\dagger_{r 2 \beta}(-\kvec ) \right ] \hat{g}_{l b \alpha \beta}^\dagger(1,\kvec ) & = &   \hat{g}_{l b \alpha \beta}^\dagger(\zeta,\kvec ) , \nonumber \\ \left [ w_1^c(\zeta ) + \epsilon_{\alpha \beta 3 } w_2^c(\zeta )    \hat{c}^\dagger_{r 1 \alpha}(\kvec )   \hat{c}^\dagger_{r 2 \beta}(-\kvec ) \right ] \hat{g}_{l c \alpha \beta}^\dagger(1,\kvec ) & = &   \hat{g}_{l c \alpha \beta}^\dagger(\zeta,\kvec ) .\end{aligned}\ ] ] these have the forms @xmath154      finally the thouless operator is : @xmath156 where @xmath157 , \nonumber \\ \hat{s}_{r b}(\zeta ) & = &   \prod_{\absum\kvec } \left [ 1 + \left(\zeta^*-1\right )   \epsilon_{\alpha \beta 3 }    \tan\theta_b^r(k )   e^{i \xi_b^r(k ) }   \hat{b}^\dagger_{r 1 \alpha}(\kvec )    \hat{b}^\dagger_{r 2 \beta}(-\kvec ) \right ] , \nonumber \\ \hat{s}_{r c}(\zeta ) & = & \prod_{\absum\kvec } \left [ 1 + \left(\zeta^*-1\right )    \epsilon_{\alpha \beta 3 } \ \tan\theta_c^r(k )   e^{i \xi_c^r(k ) } \hat{c}^\dagger_{r 1 \alpha}(\kvec )   \hat{c}^\dagger_{r 2 \beta}(-\kvec ) \right ] , \nonumber \\ \hat{s}_{r a}(\zeta ) & = & \prod_{\absum\kvec } \left [ 1 + ( \zeta-1 )   \epsilon_{\alpha \beta 3 }    \tan \theta_a^r(k )   e^{i \xi_a^r(k ) } \ \hat{a}^\dagger_{r 1 \alpha}(\kvec )   \hat{a}^\dagger_{r 2 \beta}(-\kvec ) \right ] , \nonumber \\ \hat{s}_{l b}(\zeta )   & = & \prod_{\absum\kvec } \left [ 1 + \left(\zeta^*-1\right )   \epsilon_{\alpha \beta 3 }    \tan\theta_b^l(k )   e^{i \xi_b^l(k ) }   \hat{b}^\dagger_{l 1 \alpha}(\kvec )    \hat{b}^\dagger_{l 2 \beta}(-\kvec ) \right ] , \nonumber \\ \hat{s}_{l c}(\zeta )   & = & \prod_{\absum\kvec } \left [ 1 + \left(\zeta^*-1\right )   \epsilon_{\alpha \beta 3 } \ \tan\theta_c^l(k )   e^{i \xi_c^l(k ) } \hat{c}^\dagger_{l 1 \alpha}(\kvec )   \hat{c}^\dagger_{l 2 \beta}(-\kvec ) \right ]   .\end{aligned}\ ] ]    to apply wick s theorem it is preferable to express the equations above directly in terms of the quasi - particle operators : @xmath158 | \psi \rangle , \nonumber \\ \hat{s}_{r a}(\zeta ) & ~ & \hskip-.5 cm   | \psi \rangle   = \prod_{\absum\kvec } \left [ \left ( \cos^2\!\theta_a^r(k ) + \zeta   \sin^2 \!\theta_a^r(k ) \right )   \right . \nonumber \\ & + & \left . ( \zeta-1 )   \epsilon_{\alpha \beta 3 } \sin\theta_a^r(k ) \cos\theta_a^r(k )   e^{i \xi_a^r(k ) }   \hat{\alpha}_{r 1 \alpha}^\dagger(1,\kvec ) \hat{\alpha}_{r 2 \beta}^\dagger(1,\kvec )   \right ] | \psi \rangle , \nonumber \\ \hat{s}_{r b}(\zeta ) & ~ & \hskip-.5 cm   | \psi \rangle   = \prod_{\absum\kvec } \left [ \left ( \cos^2\!\theta_b^r(k ) + \zeta^ *    \sin^2 \!\theta_b^r(k ) \right )   \right . \nonumber \\ & + & \left . \left(\zeta^*-1\right )   \epsilon_{\alpha \beta 3 }     \sin\theta_b^r(k ) \cos\theta_b^r(k ) e^{i \xi_b^r(k ) }   \hat{\beta}_{r 1 \alpha}^\dagger(1,\kvec )   \hat{\beta}_{r 2 \beta}^\dagger(1,\kvec )   \right ] | \psi \rangle , \nonumber \\ \hat{s}_{l b}(\zeta ) & ~ & \hskip-.5 cm   | \psi \rangle   = \prod_{\absum\kvec } \left [ \left ( \cos^2\!\theta_b^l(k ) + \zeta^ *   \sin^2 \!\theta_b^l(k ) \right )   \right . \nonumber \\ & + & \left . \left(\zeta^*-1\right )   \epsilon_{\alpha \beta 3 }     \sin\theta_b^l(k ) \cos\theta_b^l(k ) e^{i \xi_b^l(k ) }   \hat{\beta}_{l 1 \alpha}^\dagger(1,\kvec )   \hat{\beta}_{l 2 \beta}^\dagger(1,\kvec )   \right ] | \psi \rangle , \nonumber \\ \hat{s}_{r c}(\zeta ) & ~ & \hskip-.5 cm   | \psi \rangle   = \prod_{\absum\kvec } \left [ \left ( \cos^2\!\theta_c^r(k ) + \zeta^ *   \sin^2 \!\theta_c^r(k ) \right )   \right . \nonumber \\ & + & \left . \left(\zeta^*-1\right )   \epsilon_{\alpha \beta 3 } \sin\theta_c^r(k ) \cos\theta_c^r(k ) e^{i \xi_c^r(k ) }   \hat{\gamma}_{r 1 \alpha}^\dagger(1,\kvec )   \hat{\gamma}_{r 2 \beta}^\dagger(1,\kvec )   \right ] | \psi \rangle , \nonumber \\ \hat{s}_{l c}(\zeta ) & ~ & \hskip-.5 cm   | \psi \rangle   = \prod_{\absum\kvec } \left [ \left ( \cos^2\!\theta_c^l(k ) + \zeta^ *   \sin^2 \!\theta_c^l(k ) \right )   \right . \nonumber \\ & + & \left . \left(\zeta^*-1\right )   \epsilon_{\alpha \beta 3 }     \sin\theta_c^l(k ) \cos\theta_c^l(k ) e^{i \xi_c^l(k ) }   \hat{\gamma}_{l 1 \alpha}^\dagger(1,\kvec )   \hat{\gamma}_{l 2 \beta}^\dagger(1,\kvec )   \right ] | \psi \rangle .\end{aligned}\ ] ]
Abstract:

Answer: in this paper we study the effect of finite size on the two - flavour colour superconducting state . as well as restricting the quarks to a box , we project onto states of good baryon number and onto colour singlets , these being necessary restrictions on any observable `` quark nuggets '' . we find that whereas finite size alone has a significant effect for very small boxes , with the superconducting state often being destroyed , the effect of projection is to restore it again . the infinite - volume limit is a good approximation even for quite small systems .