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: non - newtonian fluids peculiar dynamics has focused a large attention over the last decades . their ability to flow or break @xcite , to exhibit viscous or elastic behavior @xcite , thixotropy effects @xcite or to sustain loads without flowing @xcite opened a wide range of applications , including cosmetics , food industry , environment and biology . among them , micellar fluids differ from polymers by their internal dynamics : they continuously break and restructurate @xcite , and their internal kinetics determines their length distribution - and , hence , their macroscopic properties @xcite . the behavior of bubbles rising through such fluids is complex , due , in particular , to the non - trivial coupling between the bubble and the fluid rheology @xcite . this has prevented , up to now , a complete theoretical description of the system , and favored experimental studies . due to the viscoelastic properties , the bubble shape is generally elongated @xcite , its tail ends by a cusp @xcite , and both its geometry and velocity oscillate during its rising through the fluid @xcite . in shear - thinning fluids , the local perturbation due to a rising bubble - or a falling sphere - creates a negative wake @xcite and a corridor of reduced viscosity @xcite . as a consequence , successive bubbles may interact one with the other @xcite , if the emission period is shorter than the time for the perturbation created by the bubble to vanish . an extensive study of the degassing regimes when a constant air flux is imposed at the bottom of a complex , yield - stress fluid column ( gel or immersed granular material ) has revealed the existence of three different regimes @xcite . on the one hand , at low flow - rate , bubbles rise independently one from the other . on the other hand , at high flow - rate , an open channel connects the bottom nozzle to the fluid surface ; this channel develops instabilities , forming a ` bubble chain ' @xcite . finally , at intermediate flow - rate , the system spontaneously oscillates between the two previous regimes , exhibiting a complex intermittent dynamics @xcite .    ) , via a chamber of volume @xmath0 . @xmath1 is of the order of a typical bubble diameter ( a few millimeters to centimeters ) . ]    in this article , we extend these previous works to the case of a fluid without yield stress , and for a _ thin _ fluid layer , i.e. , when the fluid height is smaller than the typical size of a bubble . we report the existence of a peculiar degassing regime in a thin layer of micellar fluid ( ctab / nasal mixture ) . when a constant air flow is imposed at its bottom , in the bubbling regime , the air is released through periodic series of successive bubbles , hereafter named _ cascades_. when the layer is thin enough , the cascades merely consists in successive opening and closing of the fluid layer above the air injection point . for convenience , these periodic degassing apertures will also be named _ bubbles _ in the following . we measure the overpressure at the base of the fluid column , at the injection point . the main characteristics of the cascades ( periodicity , maximum overpressure ) are analyzed in regards with the fluid rheology . we show that inside a cascade , the bubble properties vary continuously : the mean overpressure drops , whereas the emission time increases . the transition between two regimes inside the cascades is directly linked to the fluid rheology and , in particular , provides a direct measurement of the viscoelastic characteristic time associated with the maxwell fluid . the experimental cell consists of a cylinder made of plexiglass ( diameter 74  mm , height 270  mm ) , filled with the fluid up to a height @xmath1 ( fig .  [ fig : expsetup ] ) . air is injected at constant flow - rate @xmath2 ( from 0.17 to 1.72  ml / s ) through an injection hole ( diameter @xmath3  mm ) at the bottom of the fluid column . the air - injection system consists of a mass - flow controller ( bronkhorst , mass - stream series d-5111 ) connected to a chamber of volume @xmath0 , from which the air flows at the column bottom . the volume of the chamber can be easily tuned by changing the water level into the chamber ( see fig . [ fig : expsetup ] ) . the water itself makes it possible to inject humid air inside the fluid , thus avoiding any drying of the sample over the experimental time . a differential pressure sensor ( 223 bd-00010 ab , mks instruments ) measures the variations of the overpressure @xmath4 inside the chamber , corresponding to the pressure variations at the bottom of the cell . the fluid is a semi - dilute micellar system obtained by a mixing at equimolar concentration , inside pure water , sodium salicylate ( nasal , sigma aldrich ) and hexadecyltrimethylammonium bromide ( ctab , sigma aldrich ) . the mixture of this two chemical components causes the formation of a network of giant entangled micelles , which break down and reform continuously @xcite . on a macroscopic point of view , the fluid exhibits shear - thinning , viscoelastic properties @xcite , which can be tuned by varying the fluid concentration @xmath5 ( from 0.03 to 0.5  mol.l@xmath6 ) . the rheology of these well - controlled mixtures is characterized by rheometer measurements ( see appendix  [ app : rheology ] and [ app : modulii ] ) . the height @xmath1 of the thin layer ranges from 5 to 35  mm . unless specified , the results presented here are for a typical height of 5  mm and a concentration @xmath7  mol.l@xmath6 . when varying the system parameters ( @xmath8 ) , we observe three different degassing regimes through the thin layer : a _ bubbling _ regime , for which bubbles are emitted one after the other ; an _ open channel _ regime , for which the system is able to sustain a stable channel connecting the injection nozzle at the base of the fluid column to the fluid free surface ; and an _ intermittent _ regime , for which the system alternates spontaneously between the bubbling and the open channel regime  this latter pinching off intermittently . similar degassing regimes have been reported in yield - stress fluids : gels @xcite or in immersed granular media @xcite . phase diagram of the different degassing regimes depending on the parameters ( _ v , q _ ) for @xmath7  mol.l@xmath6 ( a ) and @xmath9  mol.l@xmath6 ( b ) [ ( @xmath10 ) bubbles , ( @xmath11 ) intermittence , ( @xmath12 ) open channel ] . dashed lines are eyeleads . the gray zone correspond to the space of parameters where bubble cascades are observed [ @xmath13  mm ] . ]    typically , when increasing the injected air flow - rate @xmath2 , all other parameters being constant , the system goes from the bubbling , intermittent and finally open channel regime ( fig . [ fig : regimes ] ) . when increasing the fluid concentration , it becomes easier for the system to sustain an open channel , and it is necessary to go to higher volumes @xmath0 and smaller flow rate @xmath2 to observe the bubbling regime ( fig . [ fig : regimes]a and b ) @xcite . on the contrary , increasing the fluid column height makes it more difficult to open a channel , and shifts the regime boundaries toward smaller volumes and higher flow rates . note that the precise boundary between two different regimes is difficult to determine , due to the finite acquisition time ( the system may not have the time to switch from bubbles to an open channel , or vice - versa ) .    the formation of an open - channel can be qualitatively explained . indeed , due to the fluid shear - thinning properties ( appendix  [ app : rheology ] ) , when a bubble rises through the fluid , its wake is characterized by a local viscosity smaller than the surrounding fluid . if the flow rate is large enough , the following bubble will rise through a fluid with a smaller effective viscosity  and thus , will rise faster , and so on , until the system is able to sustain an open channel through the fluid column . note that this channel does not resemble a cylinder , but rather a bubble chain , similar to previous observations in yield - stress and non yield - stress fluids @xcite . in a wide range of fluid concentration ( from @xmath14 to @xmath15  mol / l ) , when the fluid layer height @xmath1 is of the order of the size of a bubble ( typically , @xmath16  mm ) , we do not observe , in the bubbling regime , successive bubbles but rather periodic series of bubbles , hereafter named _ bubble cascades_. these cascades are clearly observed in the overpressure signal @xmath4 recorded at the bottom of the fluid column ( fig . [ fig : cascade ] , top ) . the pressure signal corresponding to one cascade exhibits successive rises and drops , each of them corresponding to the pressure increasing in the chamber , followed by a bubble emission ( fig . [ fig : cascade ] middle , a and b ) . after each bubble cascade , the overpressure remains at @xmath17 for a few seconds , before the next bubble cascade ( fig . [ fig : cascade ] middle , c ) . during this time interval , the air escapes continuously through the thin fluid layer . the value of the overpressure is given by the charge loss of the air flow through the hole ( roughly a few pa , fig . [ fig : cascade ] , middle , c ) . this hole then suddenly closes , and the next cascade starts . the sequence consisting of a bubble cascade followed by a hole opened through the fluid layer repeats periodically in time . the main goal of this work is to describe thoroughly the cascades properties .     separating a sequence ` bubble cascade + hole ' , as a function of the air flow - rate @xmath2 . _ inset : _ maximum overpressure @xmath18 reached in the cascades ( see fig .  [ fig : cascade ] ) . ( symbol , @xmath0 [ ml ] ) : ( @xmath10 ,  56 ) , ( @xmath19 ,  71 ) , ( @xmath20 ,  83 ) , ( @xmath12 ,  161 ) . dashed lines indicate the mean value of all the data . ] the general properties of the cascades are reported in figure  [ fig : dtpmax ] . on the one hand , we observe that the time interval between each sequence ( cascade + hole ) is very stable and , over the range of parameters explored , does not depend neither on the injected flow rate @xmath2 , nor on the chamber volume @xmath0 . the same observation is reported for the maximum overpressure reached inside each cascade ( fig . [ fig : dtpmax ] , inset ) . on the other hand , the number of bubbles emitted per cascade , @xmath21 , depends linearly on the injected flow rate ( fig . [ fig : bubnum]a ) . the associated slope , @xmath22 , varies with the chamber volume @xmath0 . for @xmath13  mm and small volumes , we report a linear , decreasing relationship @xmath23 , where @xmath24 is a constant ( fig . [ fig : bubnum]b , black dots ) . for large values of @xmath0 ( @xmath25  ml ) or larger height , however , @xmath22 is independent of the chamber volume . these results can be interpreted as follows .    for small @xmath1 , the layer is of the order of a bubble height , and whenever a bubble is emitted , it pierces the layer . via the open channel thus formed , the overpressurized air escapes the chamber . for small chamber volume @xmath0 , the chamber reservoir quickly empties , and a dependence on @xmath0 is clearly seen . when @xmath0 reaches higher values , it acts as a pressure reservoir and @xmath22 remains roughly constant .    for larger @xmath1 , the air is no longer able to pierce the fluid layer , and bubbles are emitted in the fluid , rise and burst without connecting anymore the injection nozzle to the fluid free surface . in this case , the overpressurized air trapped in the chamber can not escape directly through the fluid , but via the successive bubbles . the bubble size is fixed by the nozzle size and the fluid rheological properties , and @xmath22 is independent of @xmath0 .     as a function of the air flow rate @xmath2 . ( symbol , @xmath0 [ ml ] ) : ( @xmath11 ,  48 ) , ( @xmath10 ,  56 ) , ( @xmath19 ,  71 ) , ( @xmath26 ,  105 ) , ( @xmath12 ,  161 ) . the gray lines correspond to the linear interpolation for each series of experiments ( @xmath0 fixed ) [ @xmath13  mm ] . ( b ) @xmath22 as a function of @xmath0 . the series reported in ( a ) display a linear , decreasing relationship up to @xmath27  ml . for larger fluid layer height , the slope is almost constant . ( symbol , @xmath1 [ mm ] ) : ( @xmath10,5 ) ; ( @xmath19,8 ) ; ( @xmath20,18 ) ; ( @xmath12,27 ) ; ( @xmath11,35 ) . _ inset : _ @xmath22 as a function of the fluid concentration [ @xmath28  ml , @xmath29  ml / s , @xmath13  mm ] . no apparent relationship is found . ] the total volume of gas emitted during a cascade can be written @xmath30 . as the number of bubbles emitted per cascade always depends linearly on @xmath2 , we can write @xmath31 , where @xmath32 is a constant for a given series of experiment ( @xmath5 , @xmath1 and @xmath0 fixed ) . we thus find that the average gas volume emitted per bubble , @xmath33 , is constant . for @xmath34  mm , this constant does not depend anymore on @xmath0 ( fig . [ fig : bubnum]b ) .    inside a cascade , however , the bubble properties ( maximum overpressure and emission duration ) vary continuously ( fig .  [ fig : cascade ] , top and middle ) . in the next two sections , we investigate these variations , and see which informations they bring on the system . during the release of a single bubble , the overpressure @xmath4 exhibits three different stages ( fig .  [ fig : cascade ] , bottom ) . first , we observe a linear pressure increase ( fig . [ fig : cascade ] , bottom , region  1 ) . the overpressure signal then departs from the linear tendency ( fig . [ fig : cascade ] , bottom , region  2 ) until the bubble is emitted ( sudden pressure drop , fig . [ fig : cascade ] , bottom , region  3 ) . in this section , we describe each part of the pressure signal , and show that a simple maxwell model can account for the different observations .      when submitted to a sudden stress ( pressure increase ) , at a short time scale , the ctab / nasal mixture does not flow , and the system is equivalent to a chamber of volume @xmath0 continuously filled by a gas flow @xmath2 . the overpressure is given by @xmath35 where @xmath36  pa denotes the atmospheric pressure . the experimental slope is consistent with this linear pressure increase ( dashed gray line , fig . [ fig : cascade ] bottom , region  1 ) , without any adjustable parameters . in the following section , we estimate the threshold pressure @xmath37 above which the fluid starts flowing when the pressure at its bottom increases .    ) and a dashpot representing the viscous loss ( viscosity @xmath38 ) . @xmath39 and @xmath40 indicate the displacements generated by the rising bubble and the dashpot shift after the @xmath41 bubble , respectively . ]      after a certain time @xmath42 , associated with an overpressure @xmath37 ( point @xmath43 , fig . [ fig : cascade ] , bottom ) , the overpressure departs from its linear increase . at this point , the fluid starts flowing , and a bubble is nucleated and grows at the tip of the injection hole ( fig . [ fig : cascade ] , bottom , region  2 ) . in order to estimate this threshold pressure , we describe the fluid with a maxwell model , consisting of an ideal elastic spring attached to an ideal dashpot ( fig . [ fig : maxwell ] ) . this simple model represents the viscoelastic behavior of energy storage and viscous loss , which can be quantified by the elastic ( @xmath44 ) and viscous ( @xmath45 ) modulus , respectively . these modulii can be considered constant as a function of the applied stress in a given frequency range ( @xmath46  hz ) representative of the frequency of bubbles rising through the fluid . the modulii are estimated from the plateau obtained from oscillation measurements ( see appendix  [ app : modulii ] , fig .  [ fig : moduliisigma ] ) to @xmath47  pa and @xmath48  pa .    the equation describing the opening of a hole of radius @xmath39 due to a bubble rising in the fluid layer can be written : @xmath49 where @xmath50 is a constant . the first term in the right - hand side represents the closing force due to the fluid surface tension @xmath51 , where @xmath52  mn / m . the second term describes the elastic force which tends to shift the spring back to its initial length , where @xmath53 is a constant which can be approximated to the inverse of the nozzle radius @xmath54 . @xmath4 is the pressure inside the chamber ( fig . [ fig : expsetup ] )  and , thus , inside the bubble which starts being generated at the injection nozzle . finally , the last term quantifies the weight associated with the fluid layer height @xmath55 above the newborn bubble . at the limit where the fluid starts flowing , @xmath56 , @xmath57 and the threshold pressure , given by the condition @xmath58 , can be written @xmath59 before the emission of the first bubble in the cascade , @xmath60 and we can estimate @xmath61  pa , which is consistent with the pressure signal measured in the experiments ( fig . [ fig : cascade ] , top ) .     as a function of the bubble number in the cascade . in the main part of the cascade , @xmath62 is a linear decaying function of the bubble number @xmath21 ( symbol , @xmath5  [ mol.l@xmath6 ] ) : ( @xmath64 ) , ( @xmath65 ) [ @xmath28  ml , @xmath29  ml / s ] . ]    as the bubble properties vary continuously inside the cascade , we investigate the evolution of @xmath37 as a function of the bubble number . we find that , in the main part of the cascade , @xmath37 is a linear , decaying function of the bubble number ( fig . [ fig : dpy ] ) . the corresponding slope , constant from one cascade to the other , increases with the fluid concentration ( fig .  [ fig : dpy ] , inset ) . this linear relationship can be explained by a simple heuristic model , based on the key ingredients of the maxwell model ( fig .  [ fig : maxwell ] ) . we consider successive bubbles rising through or piercing the thin fluid layer . when the bubble rises during the characteristic time @xmath66 , @xmath56 ( fig . [ fig : maxwell ] ) and we can write the differential equation describing the temporal evolution of @xmath40 : @xmath67 which gives , right after the bubble rise @xmath68 where @xmath69 is the initial condition from which @xmath40 evolves and @xmath70 the characteristic time associated with the maxwell fluid .    after the bubble rise , we suppose that the fluid layer closes almost immediately , due in particular to the surface tension . the fluid is hence at rest , and @xmath40 relaxes towards @xmath71 during a time @xmath72 , where @xmath73 is the average time between bubble emission . for sake of convenience , we consider here that @xmath73 and @xmath66 are the same for each bubble . during this stage of relaxation , @xmath40 obeys the following equation : @xmath74 which gives , after a time @xmath75 , @xmath76 . we can therefore write the recurrence equation giving the displacement after the @xmath41 bubble : @xmath77 by recurrence , and considering that @xmath78 , we thus get the equation describing the displacement @xmath40 : @xmath79 for @xmath80 and @xmath81 ( both the characteristic times of bubble rising through the fluid layer and bubble formation and emission are smaller than the maxwell time ) , the slope @xmath82 can be approximated to @xmath83 . from eq . [ eq : dpy ] , we thus obtain the following expression for the threshold pressure : @xmath84 the average time for a bubble emission ( formation and growth , regions 2 and 3 in fig .  [ fig : cascade ] , bottom ) inside a cascade , for @xmath7  mol.l@xmath6 , can be estimated to @xmath85  s. the characteristic time @xmath86 is obtained from rheological measurements ( see appendix ) and can be estimated to @xmath87  s , with @xmath88  pa.s and @xmath47  pa . we thus estimate a slope @xmath89 of about @xmath90 , in agreement with the experimental results ( see fig .  [ fig : dpy ] for @xmath7  mol.l@xmath6 ) .      above the threshold overpressure @xmath37 , the fluid starts flowing and a bubble grows at the tip of the injection hole ( fig .  [ fig : cascade ] , bottom , region  2 ) . the overpressure then departs from its linear increase , as a consequence of the volume increase due to the bubble growth .    by deriving the general equation for an ideal gas , and denoting @xmath91 the volume of the bubble hence formed , we get @xmath92 further integration leads to the general expression of the pressure variation in time : @xmath93 where @xmath94  j  k@xmath6  mol@xmath6 is the ideal gas constant and @xmath95  l is the molar volume of air at 25@xmath96c . here we do not develop further the calculation , but note that the departure from the linear trend is linear with the bubble volume , which provides a rough estimation of this latter , from about @xmath97  ml to 10  ml from the beginning to the end of the cascade . finally , the bubble reaches the free surface  we remind here that the fluid layer is of the order of the bubble size  , the gas is suddenly released and the overpressure quickly drops ( fig . [ fig : cascade ] , bottom , region  3 ) . in order to get an estimate of the characteristic time over which the pressure drops , we write , on the one hand , bernoulli s equation to describe the air flowing from the chamber of volume @xmath0 through the opening of diameter @xmath98 , up to the surface of the fluid layer : @xmath99 , which gives the flow velocity through the opening @xcite : @xmath100 on the other hand , by considering the air as an ideal gas , we get @xmath101 , where @xmath0 is the initial gas volume , equal to the chamber volume , and @xmath91 the volume variation corresponding to the volume of the bubble connected to the injection point . the typical time to empty the chamber can be written @xmath102 , where @xmath103 is the volumetric flow - rate through the hole . by using eq . [ eq : v ] , we get the characteristic time for the pressure drop : @xmath104 where @xmath105  kg.m@xmath106 is the air density , and @xmath107  pa the atmospheric pressure . @xmath108 is of the order of a few milliseconds , compatible with the measurements of the pressure drop ( region  3 in fig . [ fig : cascade ] , bottom ) . note that the drastic pressure decrease due to the hole opening is followed by oscillations observed right after the bubble emission , which correspond to the elastic response of the fluid to the sudden stress imposed by the closing of the bubble walls , after the air release @xcite .      here , we develop qualitative arguments to estimate the lifetime of the open channel . we note @xmath109 the hole radius for which the channel remains open at the end of a cascade , and @xmath110 the associated value of the inner displacement from the maxwell model developped above ( fig . [ fig : maxwell ] ) . when the channel from the injection nozzle to the fluid free surface remains open , @xmath111 and we can write , based on eq . [ eq : motion ] : @xmath112 and @xmath113 eq . [ eq : aperture ] provides the value of @xmath110 after the last bubble in the cascade , while eq . [ eq : rstar ] gives the channel radius right after the pressure drop in the system . by substracting eqs . [ eq : aperture ] and [ eq : rstar ] and neglecting the capillary forces when the channel if fully open , we find that the channel remains open when the overpressure @xmath114 . this rough approximation is consistent with the experimental value . indeed , from fig . [ fig : cascade ] , the channel remains open once the overpressure reaches a value close to 50  pa , of the order of @xmath44 for an excitation frequency of the order of the frequency of bubbles rising through the fluid ( @xmath46  hz , see appendix  [ app : modulii ] ) .    after the channel opening , we can write the balance between the main forces at stake , the weight and viscous dissipation , @xmath115 . the characteristic lifetime of the open channel can therefore be roughly estimated to @xmath116 , of the order of a second . experimentally , @xmath117 is of a few seconds . this discrepancy can be explained first , by the rough approximation for the radius dynamic equation ; then , by the effective fluid viscosity which , for micelles , can be larger during elongational flowing than the viscosity measured under shear @xcite . = 147  ml , @xmath1=5  mm , @xmath5=0.1  mol.l@xmath6 , @xmath2=0.65  ml / s ] . _ inset : _ variations of @xmath118 with the fluid characteristic time @xmath119 , determined here as @xmath120 ( see text ) . from left to right : @xmath121 and @xmath122  mol.l@xmath6 . ]    the bubble emission characteristics vary continuously inside the cascade : through time , a bubble inside the cascade is emitted with a smaller pressure drop @xmath4 , and a longer emission time @xmath123 ( fig . [ fig : cascade ] ) . note that for the heuristic model developped in section  [ sec : dpy ] , we considered @xmath124 roughly constant and equal to @xmath125 . in this section , we investigate the variations of both @xmath123 and @xmath4 for the successive bubbles inside one cascade . figure  [ fig : dpvsdt ] displays the mean value of the overpressure @xmath126 , over a bubble release , as a function of the time interval @xmath123 over which the same bubble is emitted , for a series of bubble cascades . all data from different cascades , from the same experimental series , have been superimposed . over different cascades , the data all collapse in the same curve , for a given set of parameters ( @xmath5 , @xmath1 , @xmath2 , @xmath0 ) . the curves remain unchanged when both the volume @xmath0 and the flow rate @xmath2 vary . two different regimes can be distinguished , at short and long time scales , respectively ( fig . [ fig : dpvsdt ] , solid lines ) . the limit between both regimes defines a pressure and time threshold , respectively @xmath127 and @xmath118 . a statistical study over different volumes ( @xmath128  ml ) and flow rates ( @xmath129  ml / s ) gives a constant pressure threshold @xmath130  pa . by comparing the characteristic time @xmath118 to the maxwell characteristic time @xmath119 of the fluid ( see appendix ) , we get the direct relationship @xmath131 this relationship holds true for the three different fluid concentrations where bubble cascades are observed ( fig . [ fig : dpvsdt ] , inset ) . measuring the characteristic time linked with the bubble cascades therefore provides a direct signature on the fluid rheology . injecting air through a thin layer of micellar fluid displays a wide range of dynamic behavior . we report the existence of a peculiar degassing regime , the _ bubble cascades _ regime , for which the air is released via successive bubbles which properties ( maximum overpressure and emission duration ) vary continuously through time . this regime is observed over a wide range of parameters ( air flow - rate @xmath2 , chamber volume @xmath0 and fluid height @xmath1 ) . the cascades repeat periodically in time , separated by a few seconds during which a channel remains open between the injection nozzle and the fluid free surface , through which the air degass continuously . measuring the overpressure at the injection point makes it possible to investigate the different stages : pressure increase in the chamber , bubble formation and fluid flow , bubble emission and pressure drop . we find that the cascade periodicity and maximum overpressure depend neither on @xmath2 nor @xmath0 . the number of bubbles emitted per cascade depends linearly on the injection flow - rate @xmath2 . all the different steps of the overpressure evolution can be explained by a simple heuristic model , following the classical maxwell description of a viscoelastic fluid .    finally , we underline the interest of such simple experiment . in addition to reporting a new phenomenon ( the bubble cascades ) , we point out that measuring the evolution of the overpressure inside the cascades provides a direct insight into the fluid viscoelastic characteristic time , linked with its rheology . further work will concentrate on the microscopic behavior of the micelles . in particular , the time over which the rising bubble shears the fluid ( typically , 0.2  s ) suggests a partial alignment of the micelles ( see fig . [ fig : rheology ] , appendix  [ app : rheology ] ) . visualizing the fluid by birefringence will make it possible to determine if the micelles are aligned through time by the shear flow generated by the successive bubbles , and how this alignment is linked with the local rheology . we also propose to investigate the effective fluid viscosity under elongation , which can be much larger than the shear viscosity , and could explain the discrepancy between the estimated and measured time for the open channel lifetime ( section  [ sec : channel ] ) . the rheology of the ctab / nasal mixtures is characterized by measurements on two different rheometers : c - vor 150 , bohlin instruments and ar1000 , ta instruments . all measurements are performed with a plate - plate geometry . sand paper is glued to the plates in order to prevent any sliding at the walls ( typical rugosity of the order of 1  @xmath132 m ) . the following results do not aim at a full rheological characterization of the samples used in our experiments . they only provide the general mechanical behavior of the different mixtures , and a support for the interpretation proposed in section  [ sec : rheology ] . more detailled informations on the rheology of such systems can be found , for instance , in @xcite .     vs shear rate @xmath133 obtained by increasing @xmath134 [ @xmath7  mol.l@xmath6 ] . _ inset : _ viscosity as a function of shear rate . the arrow indicates the shear - thinning behavior [ c - vor 150 bohlin rheometer , plate - plate geometry , diameter 60  mm , gap 400  @xmath132 m , waiting time @xmath135  s per point ] . ] the semi - dilute solution behaves as a newtonian fluid at low shear rate ( @xmath136  s@xmath6 ) , with a viscosity plateau @xmath88  pa , and exhibits non - newtonian properties for higher shear rate ( fig . [ fig : rheology ] ) . the flow curve is classical for micellar fluids under shear @xcite , with a plateau in @xmath137 vs. @xmath134 between @xmath138  s@xmath6 and @xmath139  s@xmath6 . the first transition in the flow curve ( @xmath140 ) provides an access to the characteristic time associated with the viscoelastic maxwell model ( fig . [ fig : maxwell ] ) , @xmath141  s.    for @xmath142 , the flow is homogeneous and isotropic ; for @xmath143 , shear - bands appear and the flow is strongly inhomogeneous ; for @xmath144 , finally , the flow is homogeneous and nematic , with the micelles aligned in the shear direction @xcite . as the typical time for a bubble to rise up the fluid layer is @xmath85  s ( section  [ sec : dpy ] ) , we expect to be in the stress plateau and , therefore , to observe a partial alignment of the micelles ( section  [ sec : conclusion ] ) . ) and viscous ( @xmath145 ) modulii as a function of the applied stress for the ctab / nasal mixture at @xmath7  mol.l@xmath6 [ oscillation test @xmath146  hz , ar1000 rheometer , plate - plate geometry ] . ] figure  [ fig : moduliisigma ] displays the elastic ( @xmath44 ) and viscous ( @xmath45 ) moduli for the fluid @xmath7  mol.l@xmath6 , for an oscillation test ( @xmath146  hz , typically the period between two bubbles in our experiments ) . we observe that over a large range of applied stress , the elastic and viscous modulii are constant . the average plateau values of @xmath44 and @xmath45 displayed in fig . [ fig : moduliisigma ] provide another estimation of the maxwell time : @xmath147 we find @xmath87  s for @xmath7  mol.l@xmath6 , consistent with the estimation from the previous viscosity measurements ( appendix  [ app : rheology ] ) . the average plateau value strongly increases as a function of the fluid concentration ( fig . [ fig : moduliic ] ) . the values of @xmath119 for different fluid concentrations ( fig . [ fig : dpvsdt ] inset , section  [ sec : dpy ] ) have been estimated from the oscillatory measurements and eq . [ eq : taug ] . when the fluid concentration increases , the characteristic time associated with the maxwell model increases ( see appendix ) . for the same excitation ( bubble rising ) , the fluid behavior gets closer to an elastic solid response at short time ( easier channel opening ) ; at long time , due to the higher viscosity , it is more difficult to close the channel , which is therefore more easily sustained .
Abstract:

Answer: we report the experimental study of the degassing dynamics through a thin layer of shear - thinning viscoelastic fluid ( ctab / nasal solution ) , when a constant air flow is imposed at its bottom . over a large range of parameters , the air is periodically released through series of successive bubbles , hereafter named _ cascades_. each cascade is followed by a continuous degassing , lasting for several seconds , corresponding to an open channel crossing the fluid layer . the periodicity between two cascades does not depend on the injected flow - rate . inside one cascade , the properties of the overpressure signal associated with the successive bubbles vary continuously . the pressure threshold above which the fluid starts flowing , fluid deformation and pressure drop due to degassing through the thin fluid layer can be simply described by a maxwell model . we point out that monitoring the evolution inside the cascades provides a direct access to the characteristic relaxation time associated with the fluid rheology .


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: several arguments based on the quantum chromodynamics ( qcd ) show that the hadronic matter makes transition from the hadron phase to the quark - gluon one at sufficiently high temperature and/or large baryon - number density . in particular , recent intensive studies on the quark matter suggest that , in low temperature , it becomes a condensed state of quark pairs called the color superconductor ( csc ) , where the color @xmath4 symmetry is spontaneously broken  @xcite . in model calculations , the most attractive quark - quark interaction induced by gluon exchange appears in the total spin singlet ( @xmath5 ) and the color anti - triplet ( @xmath6 ) channel , so that the csc phase is considered to be a state where the quark pairs with these quantum numbers are condensed coherently in a similar mechanism with the bcs theory  @xcite .    generally , there exist two limiting cases in the fermion - pair condensation : weak- and strong - coupling regions in the strength of the two - fermion correlation . in the weak - coupling region , which appears in the high - density quark matter because of the asymptotic - freedom of qcd  @xcite , the coherence length @xmath7 becomes larger than the interquark distance @xmath8 , and the quark matter becomes the weakly - interacting fermion system where the mean - field picture are applicable . the strong - coupling region appears in the moderate- or low - density quark matter ( especially in the neighbor of the deconfinement transition point ) . large pair - fluctuations develop  @xcite and the coherence length is suggested to be very small ; @xmath9 , even around the moderate - density region [ @xmath10 from the analysis of the improved ladder approximation  @xcite . this implies that , for low - density region [ @xmath11 , @xmath12 and the strong quark - pair correlations are described by bosonic degrees of freedom , which we call `` diquark '' in this paper . thus , for the description of this diquark region of the quark matter , the incorporation of two - particle correlations beyond mean - field levels is essentially important . several methods have been proposed to incorporate the correlation effects . a typical one is that of nozires and schmitt - rink , where thermodynamic potentials are calculated in the normal phase using the in - medium ladder approximation , and the fermion number density is shown to include a contribution from the bosonic part through the two - body in - medium scatterings  @xcite . the transition temperature is obtained from the zero - frequency pole of the two - particle correlation function ( thouless criterion for second - order phase transitions ) . the quasi - chemical equilibrium theory is an another approach developed by schafroth , butler , and blatt  @xcite for the fermion - pair condensates . using a cluster expansion technique for the partition function , the fermion - pair condensates are shown to be in chemical equilibrium states between the isolated fermions and the bosonic difermions identified with the in - medium resonant / bound states of two fermions . this method was firstly invented to describe the experimental phenomena of electric superconductivity of metals in 1940 , before the birth of bcs theory . in modern sense , this quasi - chemical equilibrium picture corresponds to the strong coupling superconductor ( @xmath13 ) , and the transition temperature , which is obtained as the bec-@xmath14 of difermions under the equilibrium , is shown to be equivalent with that in nozires - schmitt - rink theory in the strong - coupling limit ( @xmath15 ) .    in the strong - coupling region of the quark matter , we can expect that the similar quark - pair fluctuations develop and the csc is described by the diquark condensate in equilibrium with isolated quarks . the diquark correlation effects have been proposed by jaffe and low  @xcite in relation with the phenomenology of the low - energy hadron - hadron collisions and the multi - hadron spectroscopy . also , roles of composite diquarks in the quark - gluon plasma have been studied by ekelin , et al . @xcite , especially for the hadronization process in the high - energy @xmath16 collision ; they considered the space - time developments of the quark - gluon plasma as an adiabatic process , and discussed the possibility of the diquark bec rather at high-@xmath17 stages in that process . also the diquark correlations have been discussed to study the stability of diquark matter but without regard to diquark bec  @xcite .    in this paper , we make a general formulation of the quasi - chemical equilibrium theory for the relativistic many - fermion system , and apply it to the strong - coupling quark matter with composite diquarks . our special interest is on diquark bec that should be the csc in the strong - coupling region , and we clarify the csc states as diquark bec at finite temperature and density within the quasi - equilibrium theory in less model - dependent way . in sec . ii , we formulate the quasi - chemical equilibrium theory for the difermion formation and the particle - antiparticle pair - creation processes in the relativistic many - body system . we also show some interesting properties for the bec of difermions which are characteristic in the relativistic system . in sec . iii , applying the method to the quark - diquark matter , we show @xmath14 of the diquark bec , a phase diagram on quark / diquark masses , and asymptotic behaviors of them . we also discuss @xmath14 of the color superconductor using phenomenological values of parameters . in sec . iv , effects of the interdiquark interactions are discussed in a mean - field approximation . v is devoted to summary and outlook . throughout this paper , we use the natural units : @xmath18 . to develop the quasi - chemical equilibrium theory , we consider molecular formation / annihilation process in a fermionic matter : @xmath19 where @xmath20 are fermions with masses @xmath21 , and @xmath22 is a composite boson of @xmath23 and @xmath24 with mass @xmath25 ; the composite boson @xmath22 should be a bound state ( @xmath26 ) or a resonance state ( @xmath27 ) . we call it `` molecule '' in both cases .    here we take the quasi - particle picture for the fermionic matter where the system consists of quasi - fermions ( @xmath23 and @xmath24 ) and quasi - bosons ( @xmath22 ) . the two - body interactions between fermions bring about the two - body correlations in the matter ; their major effects are for the binding / resonance energy of the composite boson @xmath22 and the residual interactions among quasi - particles . this quasi - particle interactions are generally regarded to be weak , and we neglect it in the first part of this paper . in sec . iv , we introduce the quasi - particle interactions and discuss their effects within a mean - field approximation . the equilibrium condition for the process ( [ equilib1 ] ) is given by @xmath28 where @xmath29 are the chemical potentials of the particles @xmath23 , @xmath24 and @xmath22 , which should be functions of the temperature @xmath17 and the particle - number densities @xmath30 . for free uniform gases , they are given by the bose / fermi statistics : @xmath31 where @xmath32 . using the relativistic energy - momentum dispersion relations : @xmath33 @xmath34 , eqs . ( [ bose1],[fermi1 ] ) become @xmath35 in the above equations , the relativistic bose / fermi functions @xmath36 and @xmath37 are defined by @xmath38 where we use scaled variables @xmath39 and @xmath40 , and the jttner variable @xmath41 that is defined by @xmath42  @xcite . the boson chemical potential @xmath43 satisfies @xmath44 from the positive norm condition of @xmath45 , and the boson density ( [ bose2 ] ) has a singularity at this upper bound , where a phase transition occurs to the bose - einstein condensate . in the bec region ( @xmath46 ) , the boson density consists of the condensed and thermal parts , @xmath47 , which are given by @xmath48 and @xmath49  @xcite .    in the process ( [ equilib1 ] ) , the particle - number conservations for @xmath20 give constraints for the densities : @xmath50 where @xmath51 are the total number densities of @xmath20 , which consist of isolated fermions and those included in the composite boson @xmath22 . solving eqs . ( [ equilib2 ] ) and ( [ densitycons ] ) with ( [ bose2],[fermi2 ] ) , we can obtain the densities @xmath52 in equilibrium at temperature @xmath17 . the energy scale in the quark - diquark system is so large that the antiparticle degrees of freedom are generally important . here we develop the quasi - chemical equilibrium theory for the particle - antiparticle system . the fermion - antifermion ( quark - antiquark ) or boson - antiboson ( diquark - antidiquark ) pairs can be created through the radiation process mediated by the gauge boson ( gluon ) ; @xmath53 thus the equilibrium condition for the above process is given by @xmath54 , where @xmath55 are the chemical potentials of the particle @xmath56 , the antiparticle @xmath57 and the gluon @xmath58 . we put @xmath59 because the gluon is a massless particle , so that the chemical potentials of @xmath56 and @xmath57 have a relation : @xmath60 .    in the process ( [ equilib3 ] ) , the difference between the particle and antiparticle densities should be conserved : @xmath61 where @xmath62 corresponds to the baryon - number density of quark or diquark degrees of freedom . as given in ( [ bose2],[fermi2 ] ) , the boson / fermion number - densities of the particle and antiparticle at the temperature @xmath17 are given by @xmath63 where @xmath36 and @xmath37 have been defined in ( [ bose3],[fermi3 ] ) . ( [ bose2 ] ) shows that the boson chemical potentials are found to have the upper bound ; @xmath64 . from the relation @xmath65 , they are also found to have the lower bound : @xmath66 where @xmath67 ( @xmath68 ) corresponds to @xmath69 ( @xmath70 ) . furthermore , the upper and lower bounds of @xmath43 in ( [ conditionmud1 ] ) correspond to the boson bec ( @xmath71 ) and antiboson bec ( @xmath72 ) respectively , which means that bosons and antibosons can not be condensed to the bec state at the same time in equilibrium states . the transition temperature @xmath14 of the boson bec ( the antiboson bec ) is determined by ( [ bdensitycons ] ) with @xmath71 ( @xmath73 ) : @xmath74_{\mu_b = \pm m_b , t = t_c}.   \label{ttba}\ ] ]      here we consider a general system of many kinds of fermions @xmath75 and molecular bosons @xmath76 which are bound / resonant states of two fermions . the fermions are assumed to be transmuted through some processes ; @xmath77 . as is same with the previous simple case , the total number density of fermions @xmath78 is to be conserved in molecular formation / annihilation process of @xmath76 , and is given by summing all species of fermions @xmath79 , which also amounts to summation of isolated fermions @xmath80 and those included in composite bosons @xmath81 : @xmath82 to obtain the equilibrium condition for the above condition , we take the helmholtz free - energy density of the system : @xmath83 where @xmath84 is the density - independent part , and the last term has been added to expose the constraint ( [ subtract1 ] ) with the lagrange multiplier @xmath85 . differentiating eq . ( [ mixed_energy1 ] ) with respect to @xmath80 and @xmath81 , one obtains for @xmath86 @xmath87 this means that the chemical potentials of fermions and bosons take the same values @xmath88 and @xmath43 in an equilibrium state , which obey the condition ( [ equilib2 ] ) ; @xmath89 . the range of the chemical potential is @xmath90 for any @xmath76 and , without loss of generality , the boson @xmath91 is assumed to be the one with the smallest mass ( most - stable composite ) ; @xmath92 . then , the range of @xmath43 becomes @xmath93 it shows that , at the upper bound of the chemical potential ( @xmath94 ) , the bec occurs for the boson @xmath91 , but there exists no bec singularity for the other bosons because @xmath95 . in summary , the bec occurs only for the boson with the smallest mass below the transition temperature @xmath14 , and the other bosons with larger masses ( less - stable composites ) can not be condensed to the bec states at any temperature and density ( one - bec theorem ) . the critical temperature @xmath14 for the @xmath91 s bec can be determined from the condition at @xmath96 : @xmath97 } _ { 2\mu_f=\mu_b \rightarrow m_{b_1 } }        = n_{ft}.   \label{upper3}\ ] ]    below @xmath14 , the density of @xmath91 consists of the thermal and condensed parts : @xmath98 -n_{b_1}^{(th ) } , \label{b1c}\end{aligned}\ ] ] where @xmath99 is the total fermion number , and @xmath100 are given by ( [ bose4],[fermi4 ] ) again . in the case of the system with antiparticle degrees of freedom , antibosons can not be condensed to bec state if the baryon - number density is positive . this can be also understood by one - bec theorem . the ground state contributions of bosons ( @xmath91 ) and antibosons ( @xmath101 ) are explicitly given by @xmath102 which means that the antiboson is effectively heavier than boson by @xmath103 and bec occurs only for bosons @xmath91 by using one - bec theorem ( in case that the baryon - number density is negative , @xmath104 , only @xmath105 can then be condensed to bec state at the lower bound @xmath106 ) . in this section , we apply the relativistic quasi - chemical equilibrium theory developed in the previous section to the color superconductivity of quark matter in bcs - bec crossover regime @xcite , where quark - pair fluctuations develop dynamically to create diquarks ( composite bosons ) . a possibility of quark - quark correlations has also been suggested in low - energy hadron physics @xcite , and its extrapolation to the deconfinement phase can be considered to become diquark degrees of freedom in dense region @xcite .      before going to actual calculations we should evaluate the applicability of quasi - chemical equilibrium theory for interacting fermionic systems , which was originally discussed by schafroth , et al . with respect to electron gas in superconducting metals @xcite :    particles of compositions ( fermions @xmath107 and composite molecules @xmath22 ) should be quasi - particles ,    correlations between @xmath108 and @xmath109 are unimportant ,    higher - order correlations are not stronger than the pair - correlations considered ,    any localizations to liquid droplets with large volume exclusion ( gas - liquid transition ) do not occur before the formation of quasi - molecules ,    composites sufficiently obey the bose statistics as expected with small size of composites relative to the interparticle distance . now we discuss the condition 1 ) for quark matter . the lifetime of a particle relates to various decay processes in medium , and to recognize quark as quasi - particle its lifetime should be at least longer than inverse of the qcd energy scale , that is , @xmath110 fm . the lifetime of quarks tends to be shorter near the fermi surface , which can be roughly estimated from the imaginary part of quark self - energy , e.g. , @xcite ; @xmath111 mev@xmath112 . in our model these quarks with short lifetime near fermi surface due to strong pair fluctuations are described as bosonic molecules , and residual quarks in the deeper interior of the fermi sphere have explicitly longer lifetime than the qcd scale . thus , quarks ( antiquarks ) being outside of pairs due to medium effects can be identified as quasi - particles . as for the diquark composites , its lifetime as quasi - particle is expected to be longer than that of quarks in crossover regime , as shuryak and zahed suggested the existence of various binary bound states in the whole regions just above the deconfinement temperature because qcd coupling strength is found not to be so small as one has believed so far @xcite . thus , diquarks should also be identified in medium . concerning other conditions 2 ) and 3 ) , higher - order correlations rather than pair ones are not treated in this sec . this treatment corresponds to the gaussian - type approximations like nozires and schmitt - rink approach @xcite . the effects of higher - order correlations , especially interdiquark interactions , will be discussed in sec . iv with sincere generalization of quasi - chemical equilibrium theory . condition 4 ) is also discussed in the end of sec . iv .    with regard to condition 5 ) , the density dependence of coherence length @xmath7 , characterising typical size of a cooper pair , relative to interquark distance @xmath8 is estimated by using the improved ladder approximation  @xcite , which suggests that @xmath13 for @xmath113 . this means that @xmath15 for lower density region @xmath114 concerned , and our quasi - chemical equilibrium approach with bosonic modes is expected to work well for these strong couplig regimes of quark matter . the critical coupling strength for the formation of diquark bound states is also calculated by using effective four - fermi interaction model with gaussian approximation  @xcite , and critical value @xmath115 ( @xmath116 is ultraviolet cutoff scale ) and its relavance for real quark matter is suggested therein . because of smooth crossover , this critical value of coupling strength corresponding to the unitary limit in cold atomic systems brings no sudden changes for physical quantities . the essential quantity for the behavior of bec - bcs crossover is always the ratio of the coherence length to the interparticle distance . although above quantitative argumets come from some specific model calculations , recent analyses on qcd experiments , e.g. , relativistic heavy ion collision  @xcite and lattice qcd  @xcite , also suggest that even in deconfinement phase qcd holds strong coupling at almost whole region of the phase diagram relevant to possible observations . we expect that our present approach works well in such strong coupling regimes in more general way , especially around the region close to the phase transition line where correlations between quarks should develop largely and the strong coupling picture becomes essential for the pairing problem . in this context three - quark correlations in addition to pairing will be important issue around this region , which can also be treated in our approach .    here we comment on the bec - bcs crossover in beta - equilibriated charge - neutral quark matter , where mismatches of fermi surfaces exist in different flavors . even in the weak coupling regime , various phases have been proposed in that matter in relation with the chromomagnetic instability problem : the gapless superconductor  @xcite , the crystalline color superconductor  @xcite , the gluon condensed state  @xcite , and so on . the strong coupling limit has also been discussed of the two - flavor quark matter with fermi - surface mismatch  @xcite . these considerations suggest that the cooper - pair formation should still be favored in the case if the energy gain by the pairing exceeds the energy differences of the mismatched fermi surfaces . thus we can expect the crossover into the diquark bec in the quark matter with finite mismatches for sufficiently large coupling strength and , at the temperature around @xmath14 , the mismatch effects are expected to be small with large pairing energy . in low temperature region , however , the fermi - surface mismatches makes the unpaired fermions remain with some flavor , and they should form new fermi surfaces ; the calculation of the profile changes of them around the crossover region should be a very interesting problem in future . in the present paper , we take the two - flavor quark matter with same fermi surfaces to discuss typical features of diquark bec ; however , our approach can be extended to the case of the finite fermi- surface mismatches and their effects in the diquark bec and crossover will be discussed in another paper .      in what follows , we take u , d - quark matter around the moderate density region in qcd , and consider the chemical equilibrium between quarks and composite diquarks in the 2sc - type pairing state ( color antisymmetric @xmath6 , flavor antisymmetric , and total - spin @xmath5 )  @xcite . the other types of diquarks can also be created in the quark matter using less attractive channels , e.g. , diquarks in color-@xmath6 , flavor symmetric and the @xmath117 pairing state  @xcite . but , according to the one - bec theorem derived in the previous section , the diquark bec can occur only for the smallest mass ( most stable ) diquark , i.e. in the 2sc - diquark state of u and d quarks , so that we consider only the 2sc - diquarks . generally the diquarks which do not participate in the bec would have the effect of lowering the critical temperature , through the total - number conservation , in comparison with the case of one kind of diquark . we also assume that the quark and diquark masses do not depend on the temperature and density , but the present formalism is also applicable to the case with effective masses dependent on these parameters .    in the quark - diquark matter , the baryon - number density @xmath118 corresponds to the conserved total fermion - number density @xmath78 in ( [ bdensitycons ] ) or ( [ subtract1 ] ) , where @xmath119 is the regularized baryon - number density in quark degrees . in the chemical equilibrium between quarks and 2sc - diquarks , it is given by @xmath120 where @xmath121 are the number densities of quark and antiquark , and @xmath122 are those of diquark and antidiquark . the numerical factors in ( [ qcd_cons1 ] ) come from the degeneracy of the internal degrees of freedom : the color - triplet ( 3 ) , the flavor - doublet ( 2 ) , the spin - doublet ( 2 ) for the quark - density term , and the compositeness ( 2 ) , the color - anti - triplet ( 3 ) , the flavor - singlet ( 1 ) , the spin - singlet ( 1 ) for the diquark - density term . ( [ bose4],[fermi4 ] ) gives the densities @xmath121 and @xmath122 in equilibrium at temperature @xmath17 : @xmath123 where @xmath124 are a quark / antiquark and a diquark / antidiquark masses , and @xmath60 ( @xmath125 ) has been used for the chemical potentials of antiquark and antidiquark . when the diquark ( antidiquark ) bec occurs for @xmath126 , the diquark ( antidiquark ) density consists of the condensed and thermal parts : @xmath127 ( @xmath128 ) , which are obtained in the same way as eqs . ( [ b1nc],[b1c ] ) . especially these condensates @xmath129 and @xmath130 are the superfluid density of csc . if the baryon - number density is positive , diquark bec can only be expected as mentioned in the previous section .    the chemical potentials @xmath131 and @xmath132 should obey the quark - diquark equilibrium condition : @xmath133 which corresponds to ( [ equilib2 ] ) . the set of eqs . ( [ qcd_cons1]-[qcd_chemical_eq4 ] ) determines a composition of the baryon - number density in the quark - diquark quasi - chemical equilibrium system .      as shown in the previous section , diquark bec occurs when the diquark chemical potential reaches the upper limit @xmath134 . in the quark - diquark system considered here , the conditions for @xmath14 are @xmath135 } _ { t = t_c } ,    \label{qcd_critical_mixed1}\end{aligned}\ ] ] where we have used eqs . ( [ qcd_cons1],[qcd_chemical_eq4 ] ) . for any fixed values of the quark and diquark masses @xmath124 , and for the baryon - number density @xmath119 in ( [ qcd_critical_mixed1 ] ) , @xmath14 is obtained by solving these equations . [ fig1 ] shows dependence of @xmath14 on @xmath136 and @xmath137 for the baryon - number densities @xmath138 where @xmath139 is the normal nuclear density . these figures show that @xmath14 increases with the baryon - number density ; especially it can be read from the shrinkage of no - bec region ( area in the front side in the @xmath140,@xmath141-plane ) . this effect can be understood that @xmath14 goes up with increase of the diquark density together with the baryon density , as seen in the behavior of @xmath14 in the ideal bose gases ( see appendix  [ appb ] ) .    [ cols="^,^ " , ]     here we should comment that the other interactions except the diquark - diquark one might give important effects , e.g. , a possible quark - diquark interaction should generate baryonic collective modes in relation with three - quark correlations . as for these higher order correlations , we can discuss these effects by extending our method , referring to higher - order cluster expansion technique  @xcite . finally we should discuss a possibility of the transition of quark matter into classical liquid with large volume exclusion . schafroth , et al . @xcite pointed out that the quasi - chemical equilibrium theory can not be applied to the interacting fermionic system which undergoes the condensation into the classical liquid before the formation of difermion molecules . to evaluate it for quark matter we now estimate the de boer parameter @xmath142  @xcite which gives the ratio of the quantum energy fluctuation to the potential depth @xmath143 : @xmath144 where @xmath145 is a range of interaction and @xmath146 is a particle mass . large @xmath142 corresponds to large quantum fluctuation relative to the potential depth @xmath143 . taking typical values of several parameters for quark matter ; @xmath147 , @xmath148 , @xmath149 , @xmath142 is roughly estimated to be @xmath150 which is very large relative to that of the quantum liquid @xmath151 ( @xmath152 ) . this large quantum fluctuation implies that quark matter can not be condensed into classical liquid and it must be a quantum gas or , rather , quantum liquid . as qualitative feature of liquid @xmath151 can be well described by free bose gas model calculations , we expect , in this sense , the quasi - chemical equilibrium picture should be applicable for quark matter . we have generally developed the quasi - chemical equilibrium theory for the system of isolated fermions and composite bosons of fermion pairs in the relativistic framework , and derived the conditions for occurrence of the composite - boson bec . especially a general property of bec in a multi - component system is suggested ; only the lightest composite bosons can be condensed to the bec state ( one - bec theorem ) . an application has been done for the diquark condensates in quark matter at moderate densities , and only 2sc - type diquarks are treated by virtue of the one - bec theorem . in this framework we have discussed the color superconductivity as the diquark bec in a possible bcs - bec crossover regime . we found the mass phase diagrams which determine the occurrence of bec at a certain temperature for various mass values of quarks and diquarks . with phenomenological values of model parameters , we obtain @xmath14 as a few tens mev , which is lower than those by model calculations within mean - field approximations , @xmath153 . in any case , @xmath14 has a large dependence on quark and diquark masses , and it varies from @xmath154 to @xmath155 along with @xmath156 ( strong coupling limit ) .    at the high-@xmath17 limit , antiparticles produced through the pair - creation processes suppress the dissociation of diquarks into quark pairs , and the ratio of the number densities are found to approach a universal value @xmath157 , which reflects that quantum statistics still remains even in the high-@xmath17 limit . we have also discussed effects of the diquark - diquark interaction , and found within the mean - field approximation that a possible repulsive interaction lowers @xmath14 by @xmath3 in comparison with that of non - interacting case through a chemical - potential renormalization . the other correlation effects beyond the mean - field approximation , e.g. , quark - diquark and three - quark correlations , should also be important in the bcs - bec crossover problem . combining the method of the beth - uhlenbeck approach  @xcite with the present quasi - chemical equilibrium theory , we should discuss the correlation effects and the crossover problem from less model - dependent point of view . the study along this line is now ongoing and will be presented near future . the authors thank t. kunihiro , t. hatsuda , h. suganuma , t. tatsumi , and j.  p. blaizot for many useful discussions . k.n . is thankful to k. tsumura and t. kojo for their critical suggestions and encouragements . in this appendix , we derive some formulas on the critical temperature of the bec transition for single bose gas in the three limiting cases : 1 ) nonrelativistic limit , 2 ) ultrarelativistic limit with a mass correction , and 3 ) ultrarelativistic limit with negative - charge particles ( antibosons ) .    here we consider the function @xmath36 defined in ( [ bose3 ] ) , the denominator of which can be expanded as @xmath158 where @xmath159 . the integral in the above expansion is evaluated to be @xmath160 d\theta   \nonumber\\            & = & \frac{1}{4 } \left\ { k_3(l\betat ) -k_1(l\betat ) \right\ }              = \frac{1}{l\betat } k_2(l\betat ) ,   \label{apa2}\end{aligned}\ ] ] where @xmath161 is the modified bessel function : @xmath162 and we have used a formula : @xmath163  @xcite . substituting eq . ( [ apa2 ] ) into ( [ apa1 ] ) , we obtain @xmath164 the transition temperature @xmath165 is obtained from eq . ( [ bose2 ] ) with @xmath166 , so that @xmath167 using the asymptotic expansion of @xmath161 for @xmath168  @xcite : @xmath169 we obtain the asymptotic expansion of ( [ apa5 ] ) for large @xmath170 ( small @xmath14 ) : @xmath171 where @xmath172 is the thermal de  broglie wave length and the @xmath173-function is defined by @xmath174 . taking the first term in the right - hand side of ( [ apa7 ] ) , we obtain the nonrelativistic limit for @xmath14  @xcite : @xmath175 to obtain the ultrarelativistic limit , we use the asymptotic expansion of @xmath176 with positive integer @xmath177 for @xmath178 : @xmath179 taking the first term of the above expansion and using it for ( [ apa5 ] ) , we obtain the density formula for small @xmath170 : @xmath180 where @xmath181 is the appell function : @xmath182 using the expansion @xmath183 for @xmath184 ( high-@xmath14 limit ) in eq . ( [ apa10 ] ) , we obtain the ultrarelativistic limit for @xmath14 : @xmath185   \nonumber\\         & = & \frac{\zeta(3)}{\pi^2 } t_c^3           + m_b \frac{\zeta(2)}{\pi^2 } t_c^2 , \label{apa12}\end{aligned}\ ] ] where @xmath186 and @xmath187 . it should be noted that the first term in the above equation does not include @xmath25 ; in the ultrarelativistic limit , the boson mass is neglected and the energy - momentum dispersion relation is approximated to be @xmath188 . the second term in eq . ( [ apa12 ] ) corresponds to the small mass correction .    in summary , we obtain the formulas of @xmath14 of the relativistic bec : @xmath189    finally we consider the boson - antiboson system , where the similar calculation with ( [ ttba ] ) gives the expansion for baryon - number density : @xmath190 the formula for @xmath14 is obtained by putting @xmath191 and @xmath192 . using eq . ( [ apa6 ] ) , we obtain the asymptotic expansion for large @xmath170 : @xmath193 where @xmath194 , defined by ( [ apa11 ] ) , corresponds to the antiboson contributions . in the nonrelativistic limit ( low @xmath14 ) , @xmath195 can be neglected because @xmath196 , so that we obtain the same formula as eqs . ( [ apa7],[apa8 ] ) where @xmath45 is replaced with @xmath197 . this means that the antiboson degree of freedom is suppressed and we can put @xmath198 in this limit .    in the ultrarelativistic limit ( small @xmath170 ) , using ( [ apa9 ] ) , eq . ( [ apa13 ] ) becomes @xmath199 where @xmath200 has been used . substituting @xmath201 and @xmath202 , we obtain the ultrarelativistic formula for @xmath14 : @xmath203 comparison with eq . ( [ apa12 ] ) shows that the leading terms from @xmath204 cancels out so that the result depend on @xmath25 in spite of the ultrarelativistic limit .    in summary , we obtain the formulas of @xmath14 for the boson - antiboson system : @xmath205 first we calculate the asymptotic expansion of fermion - number density ; calculations go in the same way as for the boson density formula ( [ apa13 ] ) : @xmath206 where the extra factor @xmath207 comes from the factor @xmath208 in the denominator of the fermion - density ( [ fermi3 ] ) reflecting fermi statistics ( bose statistics gives the extra factor @xmath209 which is embedded in eq . ( [ apa13 ] ) ) . expanding the function @xmath210 and taking the leading - order term in ( [ apa9 ] ) , we obtain the fermion - density formula for small @xmath211 : @xmath212 where @xmath213 has been used to give the factor @xmath214 in ( [ apb2 ] ) . in case of bose statistics , factor @xmath215 from the analysis of eq .  ( [ apa13 ] ) gives @xmath216 in ( [ apa15 ] ) . the difference between the factors @xmath216 and @xmath214 in eqs . ( [ apa15],[apb2 ] ) tells that the effect of quantum statistics still remains even in high-@xmath17 limit , which is discussed in eq . ( [ ratio2 ] ) . now we apply the boson and fermion density - formulas for high-@xmath17 region , ( [ apa15 ] ) and ( [ apb2 ] ) , to quark - diquark system , and evaluate the high-@xmath17 limit of the baryon - number density ratio . using eqs .  ( [ apa15],[apb2 ] ) , the baryon - density ratio of isolated quarks to diquarks is given by @xmath217 @xmath218 where @xmath219 are the baryon numbers of quark and diquark , and @xmath220 are the multiplicity factors for quark and diquark in the present system , which represent the degenerate degrees freedom ( @xmath221 ) . @xmath222 is the composite factor indicating that a diquark is constructed by @xmath223quarks ( in the present case , @xmath224 ) , and @xmath225 is the statistical factor which reflects the statistical difference between quarks ( fermions ) and diquarks ( bosons ) . m. alford , j.  a. bowers , and k. rajagopal , phys . rev . * d63 * , 074016 ( 2001 ) ; j.  a. bowers and k. rajagopal , phys . rev . * d66 * , 065002 ( 2002 ) ; e.  v. gorbar , m. hashimoto , and v.  a. miransky , phys . * 96 * , 022005 ( 2006 ) .
Abstract:

Answer: bose - einstein condensation ( bec ) of composite diquarks in quark matter ( the color superconductor phase ) is discussed using the quasi - chemical equilibrium theory at a relatively low density region near the deconfinement phase transition , where dynamical quark - pair fluctuations are assumed to be described as bosonic degrees of freedom ( diquarks ) . a general formulation is given for the diquark formation and particle - antiparticle pair - creation processes in the relativistic framework , and some interesting properties are shown , which are characteristic for the relativistic many - body system . behaviors of transition temperature and phase diagram of the quark - diquark matter are generally presented in model parameter space , and their asymptotic behaviors are also discussed . as an application to the color superconductivity , the transition temperatures and the quark and diquark density profiles are calculated in case with constituent / current quarks , where the diquark is in bound / resonant state . we obtained @xmath0 for constituent quarks and @xmath1 for current quarks at a moderate density ( @xmath2 ) . the method is also developed to include interdiquark interactions into the quasi - chemical equilibrium theory within a mean - field approximation , and it is found that a possible repulsive diquark - diquark interaction lowers the transition temperature by @xmath3 . + pacs numbers : 12.38.-t , 12.38.mh , 12.38.lg , 03.75.nt