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: interesting wave patterns often form on growing icicles that are covered with a thin layer of flowing water ( see fig.1)@xcite . for many of these patterns , the wavelength has a gaussian distribution centered at approximately 8 mm ; however , despite their common occurrence , there is no quantitative explanation for this wavelength distribution @xcite,@xcite . these waves are associated with the growth of the icicles and the flow of fluid along the icicle . hence , there are several processes occurring that should be considered . these include crystallization from the melt , latent heating at the ice - melt interface , laminar flow with two interfaces ( ice - melt and melt - air ) , evaporation of liquid , and transport of heat through the surrounding air . that the waves tend to encircle the icicle clearly indicates the importance of gravity - induced flow , although the specific interactions between the flow , the ice growth , and the heat flow through both interfaces must be considered .        in the study of crystal growth , such a surface instability is usually explained by mullins - sekerka theory ( hereafter abbreviated as ms ) . ms theory is based on two observations : laplace instability and gibbs - thomson effect . ( for detailed explanation , refer to text books and original paper by mullins and sekerka cited in @xcite . )    to a good approximation , the ice in the icicle has a uniform temperature of 273 k , thus temperature gradients into the ice are insignificant , and the external temperature field is time independent and satisfies laplace s equation . we further assume that the external temperature is below 273 k ; i.e. , the ice is not melting on average . at a convex point , the temperature gradient is higher . because heat flow is proportional to the gradient of temperature , the larger heat flow at a convex point rapidly removes latent heat from the ice surface thus allowing the convex points to increase in size rapidly . conversely , concave regions grow relatively slowly . this phenomenon suggests that short - wavelength fluctuations increase in amplitudes more rapidly than long wavelength fluctuations . we refer to this as the laplace instability . next we explain the gibbs - thomson effect . the surface of a solid object has its own energy per area called surface free energy . if a molecule attaches itself to the surface near a convex point , the surface area increases , resulting in an increase in energy . on the other hand , if a molecule becomes attached to the surface near a concave point , absorption of the molecule makes the surface area smaller . therefore , absorption of molecule at a concave point is more energy - efficient than is absorption at a convex point . for this reason , the melting point depends on the curvature of an object , i.e. , the shape of the surface area . the melting point is lower at a convex surface ( easy to melt ) and is higher at a concave surface ( hard to melt ) . such an effect suppresses the fluctuation and makes surface flat . this is called gibbs - thomson effect , which is opposite to that of laplace instability . laplace instability enhances shorter wavelength fluctuation , and gibbs - thomson effect suppresses shorter wavelength fluctuation . from these two effects , we have fluctuation of specific wavelength mainly . these two effects are incorporated in mullins - sekerka s theory @xcite,@xcite . however , the simple application of this theory can not take place in the case of icicles . the reasons are as follows . : firstly , the water layers on icicles are too thin to produce a laplace instability . because the instability requires that the water thickness be larger than the wavelength of the fluctuation . secondly , the curvature is too small . ( wavelength 1 cm is much larger than 10 @xmath0 which gt effect requires to be effective . ) therefore the nature of wave patterns along icicles can not be explained by naive ms theory . + we neglect the fluid instabilities that require turbulence because there should be no turbulence in these water films @xcite . this is because the layers are only about 100 @xmath0 thick and the flow speeds are about 2 - 4cm / s ; the resulting reynolds numbers are only about 1 and the flow is laminar . the hydrodynamics in thin water layer is also discussed by wettlaufer et.al to explain the premelting dynamics , but the discussion here is essentially different @xcite .    for our analysis , we assume water flow on a ramp ( see fig.2 ) because it is simpler to treat and relevant experiments for this geometry are available @xcite . much of the same processes and relative length scales occur in both system because the water - layer thickness ( @xmath1 ) is much smaller than the radius of icicle ( @xmath2 ) . furthermore , matsuda @xcite , observed wave patterns on such an ice ramp ; for example , at @xmath3 , the wavelength was about 8 mm . degrees.,width=226 ]    liquid flow of a thin water on a flat ramp(benney s liquid film ) produces waves @xcite(see also @xcite ) ; however , these waves travel down the ramp and thus are unlike the case on icicles . due to the explicit calculation , the wavelength of these travelling surface waves is about 1 cm , which agrees with wavelength along icicles , but they move at about 4 - 8 cm / s which is twice the speed of the fluid . benney s wave is caused by gravity and surface tension , but it is unclear how it applies to the standing waves on ice unless the travelling waves can become pinned to a fixed location ; such a pinning mechanism has not yet been proposed . our approach is to assume static flow with small ripples on the ramp surface , and then calculate the growth rate for the ripples by solving the thermal diffusion equation in the background fluid . in section 2 , we discuss the fluid dynamics of a thin water flowing along a ramp , and then in section 3 , we couple the thermal diffusion process to the flow . the thermal diffusion in air is solved in section 4 , and in section 5 we discuss the growth rate of fluctuation on icicles by combining the solutions for thermal diffusion equations in two regions ; air and water . we consider the fluid mechanics of a thin water layer with depth @xmath4 as sketched in fig.3 . over each wavelength , the average depth is @xmath5 . there are two material boundaries , one is the solid - liquid boundary sl , and the other is the air - liquid boundary al . the x - axis is along the ramp and increases in the downhill direction , whereas y- is the outward normal to the ramp . the sl surface is not flat and given by    @xmath6    al surface is given by @xmath7 where @xmath8 means average over wavelength in x - direction .        in the experimental data @xcite , surface velocity of fluid is about 3cm / s , and by using nusselt equation which is shown later @xcite , the average water - layer thickness @xmath5 was around @xmath9 mm . the wavelength @xmath10 is about 1 cm experimentally . we use the parameter @xmath11 which is @xmath12 for the typical experimental values above . in general , having @xmath13 defines the long - wavelength approximation . our key assumptions are as follows . * steady - state(time - independent ) flow . * @xmath14 ( long wavelength approximation ) * incompressible fluid . we prefer to use the following dimension - less variables . for length , @xmath15 the thickness of water layer is thus @xmath16 and the height of the sl and al surfaces are    @xmath17    the characteristic flow velocity in the x - direction is    @xmath18    where @xmath19 is the gravitational acceleration and the viscosity @xmath20 at 0 degrees centigrade . this is a nusselt equation , the theoretically predicted velocity at al surface when @xmath21 , i.e. , for a flat sl surface with uniform thickness @xcite . this unperturbative solution results by equating the gravitational force to the viscous force . the velocity distribution is parabolic in y ; @xmath22 we consider perturbations of this solution . by using the above formula for the speed , we relate the experimentally - determined flow @xmath23 to @xmath5 and , surface velocity @xmath24 . the flow quantity is defined by @xmath25 , where @xmath26 , the mean speed ; r , the radius of icicle . in the experiment @xcite , matsuda used the flow @xmath27 and width of gutter @xmath28 ( @xmath29 corresponds to @xmath30 ) because this produced the clearest waves . this gives @xmath31 . from @xmath32 and @xmath33 ( by setting @xmath3 ) , we get @xmath34 with @xmath35 . on the other hand , his measurement of the surface mean velocity by observing the motion of fine grain was @xmath36 at @xmath3 . hence , we assume @xmath37 with @xmath38 as the experimental surface speed and water - layer thickness .    in the y - direction , characteristic velocity is @xmath39 we write @xmath40 as speed in the x - direction , @xmath41 as that in the y - direction , and @xmath42 as the pressure . dimensionless speeds and pressure are given by @xmath43 other dimensionless constant parameters are , respectively , the reynolds number and the weber number @xmath44 where @xmath45 is the surface tension of liquid water . approximate values of these quantities , @xmath46 , @xmath47 , and @xmath48 predict @xmath49 and @xmath50 . this value of reynolds number indicates laminar flow . the flow components u and v satisfy the steady - state navier stokes equation for incompressible fluids : @xmath51    the incompressibility condition is @xmath52 mass conservation law at the water - air interface requires @xmath53 the boundary conditions at the sl surface are zero fluid velocity : @xmath54    stress balancing condition on al surface ; that is , free surface condition is @xmath55 the index means @xmath56 and @xmath57 is the normal unit vector to the al surface . @xmath58 is the pressure just under the al surface , and @xmath59 is the atmosphere pressure . for the above equations , we approximated surface tension term as @xmath60 by neglecting the second order term in @xmath61 . now we rewrite the equations in dimensionless form . + the incompressible fluid condition is @xmath62 to hold this condition automatically , we introduce dimensionless stream function by @xmath63    in the following , we use the stream function instead of the velocity . also , we drop the * mark on the dimensionless quantities . all the quantities are dimensionless till the end of this section . + now navier stokes equation becomes @xmath64- 2\mu^2   \psi_{xxyy } \nonumber \\ & & -\mu^3 r[\psi_x \psi_{xxy}-\psi_y \psi_{xxx } ] -\mu^4 \psi_{xxxx},\end{aligned}\ ] ] where the indices indicate derivatives with respect to @xmath65 and @xmath66 . the 4th - order derivatives appearing in l.h.s . come from viscous term with taking derivative to cancel out the pressure term in equation ( [ eq : ns ] ) . the pressure is determined from navier stokes equation by using the stream function as follows . @xmath67 or , @xmath68    for the stream function we include only 0-th and 1st orders in @xmath61 , and for pressure , only 0-th order terms are kept . the navier - stokes equation for the stream function and the pressure equations become @xmath64,\\ p_x & = & \frac{1}{\mu}+ \frac{1}{2\mu}\psi_{yyy}-\frac{r}{2}(\psi_y \psi_{xy } - \psi_x \psi_{yy}),\\ p_y & = & - \cot \theta.\end{aligned}\ ] ]    next we consider the boundary conditions up to @xmath69 . @xmath70 where we define @xmath71 as order 1 , because w is @xmath72 in our case . when @xmath73 , flat laminar flow occurs with the solution @xmath74 which is easily shown to satisfy navier stokes equation and all boundary conditions . therefore , we consider the perturbations from this solution . the precise perturbative calculations are given in appendix , and as a result we obtain the height of al surface ,    @xmath75    and the stream function given by @xmath76    in our approximation we have only 0-th order terms in @xmath61 . this is because the first order terms in @xmath61 are also proportional to the small quantity @xmath77 , the amplitude of a small ripple , which makes them effectively second order quantities . the form of stream function is intuitively understood easily . because it is just the modification of unperturbative solution for the velocity to be vanished at non flat sl boundary . now we consider the thermal diffusion process in the fluid that was obtained in the previous section . we make the following assumptions .    1 .   we use the long wavelength approximation up to first order in @xmath61 ( @xmath78 ) . we ignore thermal expansion of the water and thus retain incompressible fluid . heat transport is through steady - state thermal diffusion with flow . note that steady - state is valid because the time scale for temperature change is much longer than the time scale for the ice crystal growth . the heat flow is given by @xmath79 , where @xmath80 is the thermal conductivity of water , @xmath81 the temperature , and @xmath82 is the specific heat of water . the steady - state continuity condition is given by dropping the time derivative . @xmath83 where @xmath84 is the thermal diffusivity of water and the incompressibility condition was used . furthermore , we dropped the term for the thermal energy coming from energy dissipation of fluid @xcite , @xmath85 because this term is much smaller than the other terms . in a mean while we use dimensionless parameters @xmath86 again . from the previous section , @xmath87 where @xmath88 is the dimensionless stream function . then equation ( [ eq : thermal ] ) becomes @xmath89,\label{eq : diff}\ ] ] where we have dropped a @xmath90 term . in experimental data , @xmath91 , @xmath92 , @xmath93 , @xmath94 . this gives @xmath95 in the following sections we drop * mark on the dimensionless quantities again . we start from equation ( [ eq : diff ] ) with stream function ( [ eq : stream ] ) . the temperature at sl boundary will be very nearly equal to the melting temperature @xmath96 at atmospheric pressure . the surface tension for the curvatures in the experiment ( gibbs - thomson effect ) can alter the melting temperature by at most @xmath97 , which can be neglected . hence we expand the solution in powers of @xmath98 . @xmath99 the left and right sides of ( [ eq : diff ] ) become @xmath100 @xmath101    setting the left and right sides of ( [ eq : diff ] ) equal gives @xmath102 the coefficients are found recursively . @xmath103 all the coefficients are determined when @xmath104 is known . the first nine coefficients are the followings by using the definition @xmath105 . @xmath106    because @xmath107 , we consider only through these second - derivative terms . + on the sl surface , @xmath108 on the al surface , @xmath109 a(x ) , \label{ondo}\\ q(al ) & = & -\kappa[1 + \frac{5}{12 } \hat{d } + \frac{239}{10080 } \hat{d}^2 ] a(x),\label{flow}\end{aligned}\ ] ] where @xmath23 is the heating resulting from the temperature gradients . we have omitted the @xmath110 terms . these additional terms are @xmath111 in ( [ ondo ] ) , and @xmath112 in ( [ flow ] ) . to be comparable with @xmath113 term , @xmath114 needs to be about @xmath115 . therefore , this approximation is valid when @xmath116 . to determine @xmath117 , we must consider the temperature and heat flow at the al surface . we consider the thermal diffusion in air to consider the temperature and heat flow at the al surface . we note two points here . first , we can not approximate the icicle system as the ramp . the picture of ramp is a good approximation when we consider inside the thin water - layer , but not good for outside . therefore we treat the icicle as cylindrical object , and consider the thermal diffusion outside . second , we can not use the same dimensionless variable as before since our physical space is outside . therefore we use different dimensional coordinates in this section . the diffusion equation in air is given by @xmath118 let us write down in cylindrical coordinate . @xmath119t(r,\theta , x ) = 0.\ ] ] we assume axial symmetry , and so we have @xmath120 . therefore , we work with @xmath121 t(r , x ) = 0.\ ] ] because the surface oscillates in the x - direction , we assume that the solution has the form , @xmath122 where we have assumed that the icicle is an infinitely long column with small surface fluctuations . @xmath123 satisfies @xmath124f(r ) = 0,\ ] ] and @xmath125 satisfies @xmath126g(r ) = 0.\ ] ] the solution is , @xmath127 where @xmath128 is the 0-th modified bessel function , and a , b and c are the constants . @xmath129 is the mean radius of the icicle including the thickness of water layer @xmath5 . note that the constant c is the order @xmath130 , because it is introduced from the fluctuation on icicle . we define local coordinate @xmath66 by @xmath131 the solution is @xmath132 note that the ramp system is retained by taking the limit @xmath133 with fixing @xmath134 and @xmath135 finite .    @xmath136    where b and c are other constants . hereafter we work with finite @xmath129 case , because the ramp case is always retained by taking the limit as above . the reader might think the appearance of logarithmic term strange , since it diverges at large @xmath137 . but the appearance of such a term is natural for infinitely long axially symmetric source . as real icicles have finite length , this solution is valid only close to the icicle ; far from the icicle , the icicle acts like a point source of heat and we must mach the near - icicle and far - icicle solutions at @xmath138 where l is the length of the icicle . this matching includes the temperature at infinity and partly determines the coefficients a , b and c ; in addition , the mean growth rate of the icicle radius also determine coefficient b , which includes the information of the temperature at infinity . two boundary conditions apply to the al surface : continuous temperature and continuous heat flow across the water - air boundary . comparing these two sets of equations ( [ ondo ] ) and ( [ ondo2 ] ) , ( [ flow ] ) and ( [ flow2 ] ) , the solution requires that @xmath153 , where @xmath154 and @xmath155 are constants with dimension of temperature .      at the sl surface , the growth rate of ice is given by @xmath165    the form of growth rate is different to usual ms theory . in ms theory , @xmath166 for the surface @xmath167 . but now we have another term @xmath168 . to understand its physical meaning , we write the relative growth rate @xmath169 as the growth rate in a reference frame moving with velocity @xmath144 ( @xmath170 ) . @xmath171 the steady - state condition means , @xmath172 where @xmath173 is the hight of sl surface in reference frame . from the steady - state condition , the time scale for growth of fluctuation is very long compared to our observing time scale . by solving the equation : @xmath174 we obtain @xmath175 with relations @xmath176 the essential point is that the fluctuation is not only growing up , but also traveling downwards ( @xmath177 ) .    therefore the amplification factor is determined from f. by using the relation @xmath178 , we have @xmath179 note that @xmath180 . in the case of @xmath181 , meaning a thick icicle , we obtain @xmath182    because @xmath183 , this form is similar to the amplification factor given by mullins - sekerka theory . the amplification factor increases proportional to wave number by thermal diffusion in air as expected for a laplace instability , and it decays by interaction with fluid @xmath184 . the thermal diffusion in thin water flow works just like gibbs - thomson effect , since the fluid changes temperature distribution to be uniform and inhibit laplace instability . from these two effects , we have maximum value for @xmath185 .     versus the dimensionless wave number @xmath114 . laplace instability in air and a gt like effect due to hydrodynamics amplifies perturbations that have a wavelength given by an @xmath114 close to 2.2.,width=188 ]    the maximum amplification factor occurs when @xmath186 , which determines a preferred wavelength . @xmath187 by using @xmath188 with experimental data : @xmath189 we have @xmath190 which well agrees with experimental value of 8 mm . @xcite also , in agreement with observations , this value does not directly depend on external temperature .    using the constraints between @xmath24 , @xmath5 , and flow quantity @xmath23 as , @xmath191    then @xmath23 is related to @xmath192 as @xmath193 where @xmath3 is assumed . therefore if @xmath23 is proportional to @xmath129 for usual icicles , its wave length is uniquely determined . but for thin icicles with @xmath194 , we should include the @xmath195 term that we have neglected . next we consider the travel of fluctuations along the icicle . the traveling phase velocity @xmath196 is the following function of @xmath114 . @xmath197 where @xmath198 is defined by @xmath199 and its form is given in figure 5 . that the ripples travel down an icicle versus @xmath200 , where v is the growth rate of icicle.,width=188 ]    at the 8-mm wavelength , fluctuation travels downwards very slowly with speed @xmath201 , where @xmath144 is the mean growth rate ( @xmath202 ) . for the ambient air temperature of -8 degrees centigrade , this speed is about 1 mm per hour in experiment . because water flow carries heat flow downwards , this kind of motion seems to be natural . actually the speed @xmath203 when @xmath204 ( @xmath205 ) . then we can imagine that such speed will be increasing function of fluid velocity @xmath24 . but this is not true . for ( solid ) surface fluctuation to change its form , it is necessary to emit latent heat non uniformly . like the amplification ( growing ) factor becomes maximum at the specific value of @xmath114 , the traveling speed may have such dependence of @xmath114 . qualitatively this is due to the fact that fluid not only to carry heat downwards ( small @xmath114 region in fig.5 ) , but also make the temperature field uniform and thus suppresses the heat diffusion ( large @xmath114 region in fig.5 ) . we have shown that an icicle covered with thin water flow makes wave - like undulations on the ice during solidification , and the preferred wavelength is determined by a ms - like theory . thermal diffusion in air makes wavelength shorter , i.e. the amplification factor becomes larger for shorter wavelengths . on the other hand , thermal diffusion in thin water flow makes the wavelength larger , i.e. the amplification factor becomes smaller for shorter wavelengths . from these two effects , a specific wavelength emerges with a maximum amplification factor of fluctuation . the thermal diffusion in the thin water flow works just like the gibbs - thomson effect because the water flow makes the temperature distribution more uniform and thus inhibits the laplace instability . this then is one of our main results . our @xmath206 depends on @xmath207 and flow quantity @xmath208 . ( @xmath207 dependence will be observed for flow on inclined gutter experiment : fig2 . for an icicle , @xmath3 should be used . ) @xmath209    where we have used the relative nusselt equation @xcite , @xmath210^{1/3}.\ ] ]    the @xmath207 dependence of the wavelength for the ramp case is discussed by matsuda @xcite experimentally as @xmath211,\ ] ] with @xmath212 . however , the number of plotting data is few , and thus we can not give a definite result experimentally . @xmath213 is 0.3 in our theory . the difference exists , but not a big disagreement qualitatively .    in nature , all icicles have their own flow rates @xmath23 , but almost all have the ripples with the same wavelength . this is explained by our analysis because the wavelength depends on ratio of @xmath23 to @xmath129 , but not only on @xmath23 itself . it is then natural to assume that in nature , @xmath23 is proportional to @xmath129 . the selected wavelength of our analysis does not depend explicitly on external temperature . although the mean growth rate @xmath144 and also the amplification factor increase with a decrease of temperature , the selected wavelength with the maximum amplification factor is independent of icicle growth rate ( [ af ] ) .    we have also shown that surface ripples are expected to travel downwards during icicle growth with a speed of 0.5 times the average normal growth rate of the ice . this should be checked by experiment . our theory may be used to map the waves around mineral stalagmite by changing the diffusion equation from temperature field to a solute density field . furthermore , our diffusion equation in fluid is mathematically similar to the schrdinger equation for a harmonic oscillator having the complex valued potential . therefore , the algebraic method may be possible to be used for analysis . this point will be discussed elsewhere . the authors are grateful to prof . r. takaki , and dr . olivier for the valuable discussion . one of the authors n.ogawa thanks prof . k. fujii for the continual encouragement . the authors would like to thank prof . e. yokoyama , dr . nishimura , and prof . r. kobayashi for the helpful discussions and encouragement , and to ms k. norisue for her help in english corrections . we define the fluctuation fields @xmath214 and @xmath215 by the following relation : @xmath216    then for the fluctuation fields , we have following equations : @xmath217,\label{ns}\\ \tilde{p}_x & = & \frac{1}{2\mu}\tilde{\psi}_{yyy } -\frac{r}{2}[(-y^2 + 2y+\tilde{\psi}_y)\tilde{\psi}_{xy } \nonumber \\ & & -(-2y+2+\tilde{\psi}_{yy } ) \tilde{\psi}_x],\label{px}\\ \tilde{p}_y & = & 0.\label{py}\end{aligned}\ ] ]        and neglect higher order in @xmath61 . the solution for the 0-th ordered stream function and pressure are @xmath222 y \nonumber \\ & & + [ ( 1+\tilde{h})\eta^2 \sin^2 x + \frac{1}{3 } \eta^3 \sin^3 x],\\ \tilde{p}^{(0)}&= & -\frac{w_0}{\sin \theta}(\tilde{h}_{xx}-\eta \sin x)+(\tilde{h}+\eta \sin x)\cot \theta.\label{p0}\end{aligned}\ ] ]      which is obtained by ( [ px ] ) with the help of ( [ p0 ] ) . this expression is consistent with ( [ ns ] ) . after some integrations with the boundary conditions , we obtain @xmath224 ~y^3 + \frac{1}{2}c_1(x ) y^2 + c_2(x ) y + c_3(x).\end{aligned}\ ] ] @xmath225 are determined by the boundary conditions . the result are    @xmath226,\\ c_2 & = & -\frac{1}{2}r(1+\tilde{h})\tilde{h}_x \eta^4 \sin^4 x \nonumber\\ & + & [ \frac{w_0}{\sin \theta}(\tilde{h}_{xxx}-\eta \cos x)-(\tilde{h}_x +   \eta \cos x ) \nonumber \\ & \times&\cot \theta]\eta^2 \sin^2 x   - c_1(x)\eta \sin x,\\ \frac{\partial c_3}{\partial x}&= & -\frac{r}{5}[(1 +   \tilde{h})\tilde{h}_x]_x \eta^5 \sin^5 x \nonumber\\ & -&\frac{r}{2}(1 + \tilde{h})\tilde{h}_x \eta^5 \sin^4 x~\cos x \nonumber \\ & + & \frac{\eta^3}{3}[\frac{w_0}{\sin \theta}(\tilde{h}_{xxxx}+\eta \sin x )   \nonumber \\ & -&(\tilde{h}_{xx}- \eta \sin x)\cot \theta ] \sin^3 x \nonumber \\ & - & \frac{\eta^2}{2}\partial_x c_1(x ) \cdot \sin^2 x - \eta \partial_x   c_2(x ) \cdot \sin x.\end{aligned}\ ] ]          from the fact that @xmath234 starts from @xmath235 in its expansion , we can determine @xmath236 as , @xmath237 + c_3,\\ \partial_x c_3 & = & \eta^3 \sin x ( \frac{w_0}{\sin \theta}+\cot   \theta)\nonumber \\ & \times & [ \frac{4}{3 } \eta \sin^3 x - \eta \sin x + 3\sin^2x -2].\end{aligned}\ ] ]          99 t. tozuka , kagaku * 8 * 218 ( 1938 ) ( japanese ) ; c. knight , j. crystal growth . * 49*,193 - 198 ( 1980 ) ; n. maeno and t. takahashi , low temperature science * 43 * ( 1984 ) , 125 - 138 and 139 - 147 ( japanese ) ; l. makkonen , j. glaciol . * 34*,116 , 64 - 70 ( 1988 ) ; k. szilder and e. p. lozowski , ann . * 19 * , 141 - 145 ( 1994 ) ; n. maeno , l. makkonen , k. nishimura , k. kosugi , t. takahashi , j. glaciol . * 40 * ( 135 ) , 319 - 326 ( 1994 ) . w. w. mullins and r. f. sekerka , j. appl . * 35 * ( 1964 ) , 444 - 451 ; see also  fractal growth phenomena " , tamas vicsek , ( world scientific , singapore , 1989 ) ; for the theory of crystal growth , general text books are for example , ivan v markov,crystal growth for beginners " , world scientific , singaporle,1995 ; a. a. chernov ,  modern crystallography iii " , springer series in solid - state sciences 36 , springer 1984 .    s. r. coriell and robert l. parker , j. app . phys , vol36 , no.2 ( 1965 ) 632 ; s. r. coriell and s. c. hardy , j. app.phys , vol 40 , no.4 ( 1969);s . c. hardy and s. r. coriell , j. crystal growth 3 , ( 1968 ) 569 ; y. furukawa and k. nagashima , applied mathematical sciences * 7 * , no.3 , sept , 28 - 43 ( 1997 ) ( in japanese ) . related references are therein . ashton g. d. , heat transfer fluid mech.inst . , proc.(1972 ) 200 - 213 . ; ashton g. d. and kennedy j. f. , j. hydraul . asce 98 , ( 1972 ) 1603 - 1624 ; carey k. l. geol.surv.prof . pap . 550-b , ( 1966 ) b192 - 198 ; gilpin r. r. , hirata t. , cheng k. c. , j. fluid mech . * 99 * , ( 1980 ) 619 - 640 ; hsu k. s. , ph.d thesis . univ.iowa ( 1973 ) 147pp ; hsu k. s. , locher f. a. , kennedy j. f. , j. heat transfer 101 ( 1979 ) 598 - 602 . ; larsen p. a. , j. boston soc . eng . , 56 ( 1969 ) 45 - 67 . ; thorsness c. b. and hanratty t. j. , aiche j.,25 ( 1979 ) 686 - 697 , 697 - 701 . ; summary is given in : m. epstein and f. b. cheung , ann . fluid mech . , 15 ( 1983 ) 293 - 319 . the amplification rate of a wave ( fluctuation ) depends on its wavelength . only the wave with the largest amplification rate can survive during the development of the fluctuations , other modes will be observed as lower - amplitude noise . the constant rate of @xmath185 means that the growth is exponential in time ; hence , small differences of this constant results in large differences in fluctuation amplitude at long times . in this sense , a mode distribution like that in fig . 4 will be hidden in the external noise and thus not observed in experiments . this situation is the same as in the original ms theory .
Abstract:

Answer: quantitatively - unexplained stationary waves or ridges often encircle icicles . such waves form when roughly 0.1 mm - thick layers of water flow down the icicle . these waves typically have a wavelength of 1 cm approximately independent of external temperature , icicle thickness , and the volumetric rate of water flow . in this paper we show that these waves can not be obtained by naive mullins - sekerka instability , but are caused by a quite new surface instability related to the thermal diffusion and hydrodynamic effect of thin water flow .


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: although the cuprate superconductors have been studied for a quarter century , the exact form of the interplay of antiferromagnetic(af ) and d - wave pairing(dsc ) order at low hole doping is not known reliably . experiments on very low doped cuprates show insulating behavior and a critical hole concentration , @xmath0 , for the onset of superconductivity . the presence of a strongly varying potential due to the random distributions of the acceptors is believed to be important at very low hole densities . there are reports of a spin glass phase separating the af and dsc regions of the phase diagram @xcite . a series of neutron scattering experiments on la@xmath1sr@xmath2cuo@xmath3 found stripe order but with a change of orientation from parallel to the cu - o - cu bonds at @xmath4 to a @xmath5 direction at @xmath6 @xcite . neutron scattering experiments on yba@xmath7cu@xmath8o@xmath9 found spin density wave ( sdw ) order in the hole concentration range 0.05@xmath10@xcite . on the other hand coexistence of uniform af order with dsc has been deduced from nmr studies of multilayer hg- and ba@xmath7ca@xmath8cu@xmath3o@xmath11(f@xmath12o@xmath13)@xmath7-cuprates in a substantially larger hole density range@xcite . in this paper we analyse the role of af and dsc interlayer coupling in multilayer cuprates and their effect on the coexistence on af and dsc order .    the observation by kitaoka , mukuda and collaborators @xcite of substantial zero magnetic field shifts on a unique cu - site in the innermost layers led them to claim uniform af order in underdoped layers with @xmath14 . the multilayer cuprates are doped by acceptors in the insulating blocks between the multilayer blocks , but their random potential will be weakened in the inner layers due to screening by the more highly doped metallic outer layers . thus these multilayer cuprates offer the best possibility to reliably determine af and dsc ordering at low doping in a clean single layer . in this paper we illustrate the potential of multilayer systems by considering the case of 4-layer cuprates for which experimental results exist in a wide density range . variational monte carlo calculations for an ideal single layer strong coupling t - j model find a substantial coexistence region of spatially uniform af and dsc order up to a critical hole density @xmath14 @xcite . theoretical support for dsc order at low doping in a single layer also comes from exact diagonalization studies on clusters of the t - j model containing up to 32-sites doped with 2 holes @xcite . when extrapolated to an infinite layer , these point to a robust d - wave cooperon ( bound hole pair ) resonance at low doping . the virtual exchange of cooperons will act as a pairing mechanism for doped holes in near nodal states at low hole densities , in the presence of strong af local ( and probably long range ) order @xcite . the two order parameters to be determined are the sublattice magnetization of af order and the pairing amplitude of dsc ordering . the former is directly measured by the zero field shift on cu - sites in nmr experiments . @xcite in multilayer samples the josephson couplings between the dsc order in neighboring layers lead to the josephson plasma resonances which appear in infrared spectra with the electric field oriented along the c - axis @xcite . the energies of these resonances are proportional to the product of the dsc order parameters . the combination of these two measurements on the same samples can , in principle , be used to determine the values and the dependence of both af and dsc order parameters on the hole density . to obtain information on the interplay of af and dsc order in a single layer , we need to include the effects of interlayer couplings between the respective order parameters in the analysis . note , the interlayer spin coupling disfavors sdw order in adjacent planes with different hole densities and therefore different periodicities . as a result we restrict our attention to the case of commensurate af and dsc order .    in this paper we report on calculations using the renormalized mean field theory ( rmft ) , a method introduced in the early days by zhang and collaborators @xcite , to treat the strong coupling t - j model . we extend the method to multilayer t - j models . rmft calculations@xcite generally agree well with more accurate variational monte - carlo calculations @xcite in the case of a single layer .    to treat multilayer materials we require the values of the interlayer couplings . the heisenberg spin - spin coupling has been directly measured with a strength typically of order a tenth of the intralayer coupling . the interlayer josephson coupling constants are less well known , as will be discussed further below . a second key input is the value of the hole densities in inequivalent layers . kitaoka and collaborators have measured the cu - knight shift in each layer and used the room temperature value as input in an empirical formula to estimate the hole density in each layer . we shall return to these inputs below . we start by defining the t - j model for a single plane . @xmath16 for the @xmath17 layer is given by @xmath18 where @xmath19 is an annihilation operator of a spin @xmath20 electron at site @xmath21 in layer @xmath22 , @xmath23 is the gutzwiller projection operator which enforces the no - double occupancy condition in cuprates , @xmath24 is the hopping integral between site i and site j with nearest neighbor hopping @xmath25 , next nearest neighbor hopping @xmath26 , and @xmath27 for the remaining site pairs . j is the superexchange coupling between nn sites .    the interlayer coupling , @xmath28 , between two neighboring layers @xmath22 and @xmath29 includes a superexchange spin - spin coupling and an interlayer hopping term . the interlayer superexchange coupling is given by @xmath30 and the interlayer hopping term reads @xmath31 where@xcite @xmath32 , with @xmath33 .    in a system with more than two layers , the hole concentrations on two adjacent inequivalent layers are different and the fermi surfaces of the two layers are mismatched . because planar momentum is conserved in interlayer hopping , the mismatch of the fermi surfaces means that direct interlayer single particle hopping and interlayer pairing can be neglected . however a cooper pair can tunnel from the fermi surface of one layer into a pair state off the fermi surface on a neighboring layer , conserving the planar momentum of the individual electrons , this process can be written down in standard perturbation theory as @xmath34 with @xmath35 a characteristic energy . this excited state cooper pair can subsequently relax to the fermi surface of the second layer through interactions such as the j - term in the single plane hamiltonian . this virtual tunneling process through pair states off the fermi energy will be the leading contribution of the interlayer josephson coupling , arising from cooper pair hopping between inequivalent layers leading to an interlayer josephson coupling @xmath36 where @xmath37 is the pair hopping amplitudes , @xmath38 is a characteristic energy related to the intralayer relaxation process , and the prime on the summation denote s that the summation are only in the vicinity of fermi surfaces with cut - off @xmath39 . schematic illustration for the hole carrier charge distribution @xmath40 , electric field @xmath41 ( direction indicated by arrows ) in 4-layer cu - oxides . horizontal lines represent layers . sub - indices @xmath21 and @xmath42 represent inner and outer planes , respectively . the top and bottom dashed lines represent chemically doped charge.,scaledwidth=26.0% ]    in multilayer materials the distribution of hole densities on the individual layers plays an important role . theoretical estimates of the hole densities can be made using a hartree approximation scheme by including the electrostatic energy in the hamiltonian @xcite . for simplicity , we consider each layer as an infinite plane with homogeneous charge . the electric field generated by such a plane is @xmath43 , where @xmath44 is the relevant dielectric constant , @xmath40 is the hole concentration of the layer , and @xmath45 is the lattice constant . for a 4-layer system , we denote @xmath46 and @xmath47 as the hole concentrations of the outer layers and the inner layers , respectively . the average hole concentration is @xmath48 , while the electron concentration of the acceptor layer above and below the 4-layers is @xmath49 . the spatial distribution of electric field is depicted in fig . [ fig : charge ] . we are only interested in the terms relevant to the charge inhomogeneity , so the electrostatic energy between acceptor layer and outer layer can be ignored . then the electrostatic energy per plaquette in the 4-layer case is @xmath50 , where @xmath51 is the distance between two adjacent cuo@xmath7 layers . the hamiltonian for multilayer system then reads @xmath52 the dielectric screening of the interlayer electric field arises from displacements of the ionic layer between the neighboring layers and is not identical to the static c - axis dielectric screening due to optical phonons etc . we adjust the value of @xmath44 to fit the estimated values for the layer densities ( see below ) . in this section , we will use rmft to analyse the multi - layer hamiltonian . the main element of the rmft theory is the so - called gutzwiller approximation , which replaces the gutzwiller projection in the hamiltonian by simple numerical factors . let @xmath53 and @xmath54 be the average values of operator @xmath55 in the gutzwiller projected state and in the unprojected state , respectively , then @xmath56 , with @xmath57 a numerical factor , depending on the hole density and the process associated with the operator @xmath55 . the coexistence of af and dsc order in the t - t-j model for a single layer material has been investigated within rmft recently by k .- y . yang et al . @xcite . here we generalize their calculations to a multilayer system . note that in the rmft the superconducting order parameter is non - zero at wave - vector off the fermi surface . the virtual tunneling process through pair states off the fermi energy when the two layers have different hole concentrations we discussed in section ii may also be taken into account by including the pair scattering on the same layer along with the second order interlayer hopping process conserving momentum of single electrons . below we shall use eq . ( 4 ) for the interlayer tunneling directly to study the multi - layer system described in eq . we now introduce the following mean fields @xmath58 where @xmath59 indicates average with the unprojected state . since we are only interested in the ground state , we can assume that the mean fields @xmath60 are all real . in the presence of af order , the magnetization on the sublattices a and b are opposite , and the mean fields @xmath61 from nnn sites depend on the sublattice and spin . we introduce the following mean fields , @xmath62 and @xmath63 . since @xmath64 , we can further simplify the notation as @xmath65 and @xmath66 . on the other hand , the mean fields @xmath61 for nn sites are independent of the sublattice @xmath67 and spin @xmath20 , so we can set the mean fields @xmath68 . the hole concentration on each layer satisfies the condition . @xmath69    then we consider the gutzwiller factors . in intralayer terms , we adopt the same form of the factors @xmath70 , and @xmath71 as used previously for the single layer t - j model @xcite . for the interlayer hopping and superexchange coupling terms , we make a simple assumption that the gutzwiller factors are @xmath72 and @xmath73 , respectively . the total energy reads @xmath74 where @xmath75 is given by the rmft for single layer t - t-j model @xcite , the electrostatic energy @xmath76 is specified in the previous section , and the interlayer energy @xmath77 takes the form @xmath78 with @xmath79 . note that the magnetic moment @xmath80 and superconducting order parameter @xmath81 are related to the mean fields by the gutzwiller renormalization factors@xcite , @xmath82    these approximations lead to the mean field hamiltonian @xmath83 where @xmath84 and @xmath85 , and @xmath86 . the order parameters and hole concentrations calculated in 4-layer t - t-j model eq . ( 5 ) . model parameters are @xmath87 , @xmath88 , interlayer pair hopping @xmath89 , interlayer spin coupling @xmath90 . upper panel : the hole concentrations for each layer @xmath91 vs average hole concentrations @xmath40 . thick blue curves are for inner layer ( @xmath47 ) , thin red for outer layer ( @xmath46 ) . lower panel : magnetic moment @xmath29 and superconducting order parameter @xmath92 as function of hole doping for each layers @xmath91 . thick blue curves are for inner layer , thin red for outer layer . left panel : dielectric constant @xmath93 , right panel : @xmath94 . dashed curves in ( b ) and ( d ) are for single layer results , plotted for comparison.,scaledwidth=40.0% ]    the mean field energy defined by eqn ( 8) has been minimized numerically for the case of 4-layers . the results depend on the hole densities which in turn depend on the value of the dielectric constant , @xmath44 . in fig . [ fig : orders ] , the two upper figures show the density in each layer as function of the total doping for two choices @xmath95 and @xmath96 . the lower figures show the resulting values of the order parameters in each layer . the single layer result is also shown in the figure as a dashed line for comparison . the sc order parameters of the 4-layer system and single layer system are almost equal and independent of @xmath44 . this is a consequence of the weak josephson coupling between adjacent layers in the multilayer system due to the mismatch of the fermi surfaces and the suppression of the interlayer hopping by the strong onsite coulomb interaction . .[tab : epsilon ] comparison between nmr experimental results on ba@xmath7ca@xmath8cu@xmath3o@xmath11(f@xmath12o@xmath13)@xmath7@xcite and the renormalized mean field theory for 4-layer @xmath97 model . @xmath98 is the content of @xmath99 in the compound . @xmath100 and @xmath101 are the hole concentrations on the outer and inner planes . @xmath102 and @xmath103 are the af magnetic moment on the outer and inner planes , respectively . @xmath104 is the fitted value of the dielectric constant from the hole distribution . @xmath105 and @xmath106 are the josephson plasma frequencies between two inner layers and between inner layer and outer layer , respectively . [ cols="^,^,^,^,^,^,^,^,^,^ " , ]     the interlayer af coupling slightly enhances the afm moment and critical doping @xmath0 of the inner layers and leads to a finite afm moment on the outer layers when their hole concentration is much larger than the single layer critical doping @xmath0 . such a long tail arises because of the lower hole concentration on the inner layers . a smaller @xmath44 leads to larger charge imbalance between the outer and inner layers and so to larger @xmath0 of outer layers , as shown in fig . [ fig : orders ] . if one neglects the long tail of the afm moment of the outer layers , the phase diagram of multi - layer system is very similar with that of a single layer system , independent of the value of the dielectric constant @xmath44 . this indicates that experiments on multilayer system can capture the essential physics of single layer t - j model    we turn now to the relation of the model with experiment . the two order parameters we investigated above are the afm magnetization and the pairing amplitude of dsc ordering . the former has been directly measured by the zero field shift in nmr . in the nmr experiments , shimizu , kitaoka and collaborators @xcite have also used an empirical scaling of the cu knight shift , measured at room temperature , to deduce hole densities of individual layers . note this scaling form has evolved with time . in tab . [ tab : epsilon ] , we show the value of the relative dielectric needed to fit the charge distribution recently determined from nmr experiments on the 4-layer ba@xmath7ca@xmath8cu@xmath3o@xmath11(f@xmath12o@xmath13)@xmath7 cuprate . note , in these oxy - fluoride cuprates the hole density differs substantially from the nominal f- concentration , y . though the dielectric constant required to fit the hole density changes from 87 to 247 as the total hole density decreases from 0.46 to 0.264 , the physics in each layer should not change very much as analyzed above . in table 1 , we also show the afm moments calculated by the rmft and measured by nmr experiments . the former values are considerably larger than the latter , indicating an over estimation of the afm order in the rmft theory . the comparison between theory and experiment is sensitive to the values of the hole density . for example , if we scale the value of the estimated hole density so that the optimal density with the highest tc for a single layer is @xmath107 rather @xmath108 quoted in mukuda et al.@xcite , this leads to an increase in the input values for the hole density in the rmft calculations by a factor of 1.25 . then the calculated values of the af moment shown in tab 1 , are smaller , although still too big , and in better agreement with experiment . turning to the dsc ordering , as mentioned earlier the pairing amplitudes enter into the josephson plasma frequencies observed in the infrared spectra with the electric field oriented along the c - axis . in multilayer cuprates , there are two kinds of josephson plasma frequencies which correspond to inter - multilayer and intra - multilayer josephson couplings respectively . in this paper , we are only interested the latter , i.e. the josephson plasma frequencies related to the josephson coupling between two adjacent cuo@xmath7 layers in same unit cell . the josephson energy between those two layers can be written as @xmath109 , where @xmath110 and @xmath111 are the phase of the superconducting order parameter of the two layers respectively . for the two layers with the same hole concentration , @xmath112 is given in our rmft by , @xmath113 where @xmath114 is the anomalous green s function of layer @xmath22 , @xmath115 is the matsubara frequency , @xmath116 , and @xmath117 indicates the summation is only over the reduced zone . the josephson plasma frequency which corresponds to the josephson coupling between those two layers is proportional to the josephson coupling constant @xcite , i.e. @xmath118 . so we have @xmath119 where @xmath120 is a constant related to @xmath121 and finite frequency dielectric constant @xmath122 . for simplicity , we can assume that @xmath120 is weakly dependent on the hole concentrations and materials . for a 4-layer material , there are two intra - unit cell josephson plasmon frequencies @xmath106 and @xmath105 which correspond to the josephson coupling between two inner layers and between the inner and outer layers respectively . a recent optical experiment on a 4-layer hg - based compound has shown that the josephson plasmon frequencies are 360 @xmath123 for inner layers and 540 @xmath123 for outer and inner layers respectively leading to a ratio of @xmath124@xcite . according to eqn . , @xmath125 . using the estimated hole concentrations of the hg - compound , the ratio calculated with rmft is around @xmath126 which is fairly close to the experimental value . to make further comparison , we calculate the factor @xmath120 with @xmath127 , where the plasmon frequencies are the experimental values for the hg - compound and the @xmath128 is calculated with rmft . using this estimate for @xmath120 allows us to obtain values for the plasma frequencies for the 4-layer ba@xmath7ca@xmath8cu@xmath3o@xmath11(f@xmath12o@xmath13)@xmath7 material . the results are quoted in table [ tab : epsilon ] . in this paper we presented a series of calculations using the rmft to treat the af and dsc order parameters in 4-layer cuprates . the rmft calculations / citestripe - rvb and also the results of single vmc calculations / citevmc , give a more extended region of coexistence of the two order parameters than that observed in single and bilayer cuprates . the extended coexistence region in the 4- layer cuprates agrees with the conclusion of the nmr experiments in the inner layers . in our analysis the extended region comes not from strong interlayer magnetic and josephson coupling , but from the extended coexistence region already present in single layers in the rmft approximation . this analysis suggests that the difference in the density range of coexisting order in multilayer cuprates compared to experiments on single and bilayer cuprates , is due to the better screening of the external potential modulation in the inner layers of multilayer cuprates , when compared to single and bilayer cuprates . the latter are adjacent to the random acceptors , while the metallic outer layers screen the random potential at the inner layers . in principle the combination of zero field nmr and infrared josephson plasma frequencies allows one to determine both af and dsc order parameters . however , in order to go further and directly compare measured and calculated order parameters a better determination of the interlayer josephson coupling constant and its variation with the hole densities in the individual layers , is required . to this end a series of measurements of the josephson plasma frequencies and the af moments by zero field nmr at different hole densities on the same multilayer material would be helpful . 99 for a review see m .- h.julien , physica b * 329 - 333 * , 693 ( 2003 ) s. wakimoto , r. j. birgeneau , m. a. kastner , y. s. lee , r. erwin , p. m. gehring , s. h. lee , m. fujita , k. yamada , y. endoh , k. hirota , and g. shirane , phys . b * 61 * , 3699 ( 2000 ) d. haug , v. hinkov , y. sidis , p. bourges , n. b. christensen , a. ivanov , t. keller , c. t. lin , and b. keimer , new j. phys * 12 * , 105006 ( 2010 ) m. ogata and h.fukuyama , rep . prog . phys . * 71 * , 036501 ( 2008 ) a. l. cheryshev , p. w. leung , and r. j. gooding , phys . rev . b * 58 * , 13594 ( 1998 ) for a recent review see h. mukuda , s. shimizu , a. iyo , and y. kitaoka , j. phys . soc . jpn . * 81 * , 011008 ( 2012 ) a. j. leggett , prog . 36 * , 901 - 930 ( 1966 ) d. van der marel , j. supercond . * 17 * , 559 ( 2004 ) f. c. zhang , c. gros , t. m. rice , and h. shiba , supercond.sci.technol . * 1 * , 36 ( 1988 ) k .- y.yang , w. q. chen , t. m. rice , m. sigrist , and f. c. zhang , new j. phys . * 11 * , 055053 ( 2009 ) g. j. chen , r. joynt , f. c. zhang , and g. gros , phys . b * 42 * , 2662 ( 1990);t . giamarchi and c. lhuillier , phys . b * 43 * , 12943 ( 1991 ) ; a. himeda and m. ogata , phys . rev . b * 60 * , r9935 ( 1999 ) s. chakravarty , a. sudbo , p. w. anderson , and s. strong , science * 261 * , 337 ( 1993 ) ; o. k. andersen , a. i. liechtenstein , o. jepsen , and f. paulsen , j. phys . solids * 56*,1573 , ( 1995 ) m. di stasio , k. a. mueller , and l. pietronero , phys . lett * 64 * , 2827 ( 1990 ) j. y. gan , y. chen , z. b. su , and f. c. zhang , phys . lett . * 94 * , 067005 ( 2005 ) m. ogata and a. himeda , j. phys . soc . jpn . * 72 * , 374 ( 2003 ) s. shimizu , s. iwai , s. tabata , h. mukuda , y. kitaoka , p. m. shirage , h. kito , and a. iyo , phys . b * 83 * , 144523 ( 2011 ) y. hirata , k. m. kojima , s. uchida , m. ishikado , a. iyo , h. eisaki , and s. tajima , physica c : superconductivity 470 , s44 ( 2010 ) t. m. rice , k .- y . yang , and f. c. zhang , rep . . phys . * 75 * , 016502 ( 2012 )
Abstract:

Answer: a more extended low density region of coexisting uniform antiferromagnetism and d - wave superconductivity has been reported in multilayer cuprates , when compared to single or bilayer cuprates . this coexistence could be due to the enhanced screening of random potential modulations in inner layers or to the interlayer heisenberg and josephson couplings . a theoretical analysis using a renormalized mean field theory , favors the former explanation . the potential for an improved determination of the antiferromagnetic and superconducting order parameters in an ideal single layer from zero field nmr and infrared josephson plasma resonances in multilayer cuprates is discussed .


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: starting with the pioneering work of burlaga ( @xcite ) a considerable number of studies have been dedicated to the analysis of the intermittent nature of the solar wind turbulence ( see @xcite for a review ) . in particular , the phenomenon of intermittency , as the inhomogeneity of the energy transfer rate in the turbulent nonlinear cascade , has been identified in fluctuations of plasma variables and of the magnetic field in the entire heliosphere @xcite , in the distant heliosheath @xcite , and near the heliopause @xcite . the level of intermittency in the solar wind turbulence has been mainly determined by the analysis of the anomalous scaling of structure functions @xcite , by considering probability distribution functions ( pdfs ) of fluctuations and departure from gaussianity for smaller scales @xcite , by using the scale behavior of the flatness ( or fourth - order moment ) and fitting of the pdfs @xcite , or by fitting models of turbulence @xcite . + multifractal analysis is a higher - order method used to reveal the intermittent nature of various characteristics of turbulence ( see , for instance , @xcite ) . in particular , the degree of multifractality is a measure of the level of heterogeneity in the local scaling indices ; in other words , it informs about the departure from a strict self - similarity ( e.g. , * ? ? ? thus , it is also a quantitative measure of the intermittency successfully used in many studies of turbulence @xcite . + the analysis of observations at the ecliptic showed that fast solar wind is generally less intermittent than slow wind for both wind speed and magnetic field components @xcite . moreover , @xcite and @xcite studied the intermittent properties of the fast wind in the inner heliosphere from helios data ( between @xmath0 and @xmath1 au ) and concluded that fast wind intermittency increases with the distance from the sun . + considerable effort has been devoted to investigating the evolution of intermittency in interplanetary magnetic fields with distance and latitude in the solar wind ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? in particular , @xcite focused on the different intermittent levels of magnetic field fluctuations measured by _ ulysses _ during solar minimum ( 1994 - 1996 ) and solar maximum ( 2000 - 2001 ) and found a high level of intermittency throughout the both phases of the solar cycle . moreover , they showed that the slow wind has a lower level of intermittency as compared with the fast flow . however , these authors used only the wind speed as a criterion to discriminate between pure coronal fast and typical equatorial slow wind . further analysis of pure polar coronal fast wind at solar minimum between 1994 and 1996 @xcite shows that inertial range intermittency at time scales between 40 and 200 seconds , described as departure of pdfs from gaussianity , increases with increasing radial distance from the sun . + @xcite performed a detailed selection of the so - called  pure " states of the solar wind encountered by _ ulysses _ between 1992 and 1997 . it was shown that slow wind is more intermittent than fast wind and no radial evolution was found for the slow wind . however , the heliocentric radial variations during the time intervals studied by @xcite were rather limited ( @xmath2 au ) , as was the heliolatitude range ( @xmath3 ) . there are very few comparative studies of solar wind intermittent turbulence at solar minimum and solar maximum out of the ecliptic plane . therefore , a systematic analysis of data for a larger range of radial distances and heliolatitudes may allow a better understanding of the nature of intermittency out of the ecliptic plane . + in this study , we use _ ulysses _ data recorded during the last two minima of the solar cycle ( 1997 - 1998 and 2007 - 2008 ) , as well as during the previous maximum ( 1999 - 2001 ) ; the database includes 130 time intervals of  pure " slow and fast solar wind , whose selection procedure is described in the next section . based on these data , we construct a more complete picture of the multifractal scaling and intermittent nature of the turbulence beyond the ecliptic plane . this is the comprehensive study of such a magnitude devoted to the multifractal spectrum of fluctuations of the interplanetary magnetic field strength measured at different heliocentric distances and heliographic latitudes during two solar minima ( 1997 - 1998 , 2007 - 2008 ) and a solar maximum ( 1999 - 2001 ) . the selected data cover a radial range between 1.5 and 5.4 au and a heliolatitude range between @xmath4 and @xmath5 . _ ulysses _ is a joint esa - nasa mission launched in 1990 october . the probe completed roughly three solar orbits before the mission ended in 2009 june . the orbit of _ ulysses _ was periodic ( 6.2 years ) , with perihelion at 1.3 au and aphelion at 5.4 au and a latitudinal excursion of @xmath6 allowing the study of both latitudinal and radial dependence of the solar wind over a wide range of solar activity conditions @xcite . in this study , we focused on the slow and fast solar wind data measured in situ by the _ ulysses _ during two solar minima ( 1997 - 1998 , 2007 - 2008 ) and a solar maximum ( 1999 - 2001 ) . to identify the solar wind type and origin we define a set of criteria and thresholds applied on plasma and magnetic field data . thus , we use five parameters to discriminate between the fast and slow solar wind state : ( 1 ) the radial velocity , ( 2 ) the proton density , ( 3 ) proton temperature parameters measured by the swoops instrument @xcite , ( 4 ) the oxygen ion ratio ( @xmath7 ) from swics @xcite , and ( 5 ) the magnetic compressibility factor calculated from magnetic field measurements by the vhm - fgm magnetometer @xcite . for each of the five parameters , we defined a threshold value as a transition from slow to fast wind:(1 ) @xmath8 km s@xmath9 , ( 2 ) @xmath10 @xmath11 , ( 3)@xmath12 k , ( 4 ) @xmath13 , and ( 5 ) @xmath14 @xcite . we also take into account the solar cycle variation of these thresholds ( e.g. , * ? ? ? * ) and exclude time intervals corresponding to interplanetary transients like shocks and cmes , identified by gosling and forysth ( see http://www.sp.ph.ic.ac.uk/ulysses/shocklist.txt ) and gosling and ebert ( see http://swoops.lanl.gov + /cme_list.html ) . we identified @xmath15 time intervals of fast solar wind ( @xmath16 during solar maximum and @xmath17 during solar minimum ) and @xmath18 time intervals of slow solar wind ( @xmath19 during solar maximum and @xmath20 during solar minimum ) . the multifractality is checked for each data set that contains a few days long time series of the magnetic field strength measured with a @xmath21 hz sampling rate . the multifractal approach is a generalization of fractal analysis used to characterize objects that demonstrate a variety of self - similarities at various scales . a detailed description of the multifractal theory can be found in , e.g. , the works of @xcite or @xcite . one of the basic characteristics of multifractal scaling is the multifractal spectrum @xmath22 , sketched in figure  1 of the paper by @xcite . as a function of a singularity strength @xmath23 , @xmath22 portrays the variability in the scaling properties of the measures . in particular , the width of the spectrum quantifies the degree of multifractality of a given system ( e.g. , * ? ? ? * ; * ? ? ? * ) . + there are several techniques to obtain multifractal spectra from the experimental data . in this study , we use two methods : the partition function technique @xcite and the direct method of determination of a multifractal spectrum proposed by @xcite , both successfully tested in many situations for space plasma turbulence @xcite and briefly described below . + in the first step of the analysis , we construct a multifractal measure @xcite from the first moment of increments of magnetic fluctuations at a scale @xmath24 , @xmath25 where @xmath26 , with @xmath27 from @xmath28 to @xmath29 denoting the data set length as a power of 2 , while @xmath30 denotes the magnetic field strength . next , we decompose the signal in segments of size @xmath24 and associate with each @xmath31th segment a probability measure defined by @xmath32 this quantity can be interpreted as a probability that the portion of fluctuation is transferred to a segment of size @xmath24 . as is typical , at a given position @xmath33 , where @xmath34 is the average solar wind speed , the temporal scales measured in units of sampling time @xmath35 can be interpreted as the spatial scales @xmath36 ( the taylor hypothesis ) next , the scaling of the partition function , @xmath37 , for various box size @xmath24 and moments of order @xmath38 is considered according to @xmath39 where @xmath40 is the scaling exponent . the scaling exponents @xmath40 are extracted from the slopes of the logarithm of the partition function @xmath41 versus @xmath42 ( eq . ( [ e : w : spf ] ) ) . in the next step , @xmath22 is determined by a legendre transformation of the @xmath40 curve with @xmath43= \tau'(q)$ ] and @xmath44 @xcite . in the case of the direct method of determination of a multifractal spectrum , the so - called pseudoprobability measure , @xmath45 , is determined from @xmath46 next , the multifractal spectrum can be calculated directly from the average of the pseudoprobability measure @xmath47 in eq . ( 7 ) @xcite : @xmath48 where @xmath49 indicates ensemble averaging . the corresponding average value of the singularity strength is given by@xcite @xmath50    the scaling range where a linear regression is applied to determine @xmath40 , @xmath51 , and @xmath52 is not known a priori , and it must be chosen before the analysis . in this study , we apply an automatic selection procedure proposed by @xcite . in most cases , we identified a scaling range from 4 s to a few days that can be somewhat related to the existence of the inertial range and with the dissipation effects observed on smaller scales . cases , for which the optimal scaling range could not be determined , have been rejected from a further analysis . + the two points @xmath53 and @xmath54 , at which @xmath55 , quantify the degree of multifractality @xcite : @xmath56 owing to the limited amount of data available from _ ulysses _ in the solar wind , we can only determine the points near the maximum of @xmath22 . therefore , in the last step of the analysis we fit the observations with a theoretical two - scale binomial multiplicative model @xcite to extrapolate the values of @xmath54 and @xmath53 and to obtain the degree of multifractality . + it is worth to noting that the agreement of the multifractal spectrum with the phenomenological intermittent two - scale model confirms that multifractal structures detected at small scales are generated by an intermittent cascade . moreover , the degree of multifractality calculated by using fitted two - scale model parameters ( @xmath57 ; for more details , see , e.g. @xcite ) measures inhomogeneity of the considered fluctuations and their deviation from a strict self - similarity . that is why @xmath58 is also a measure of intermittency , which is in contrast to self - similarity ( * ? ? ? figure  [ f : w : sp ] presents an example of an @xmath22 multifractal spectrum obtained using the experimental values of the magnetic field strength for the fast solar wind measured by _ ulysses _ in 2007 ( from january 22 to 28 ) at @xmath59 au and @xmath60 . the continuous line shows the @xmath22 spectrum derived from the theoretical two - scale model @xcite fitted with the observations ; the theoretical fit allows us to determine the width of the multifractal spectrum as the measure of multifractality ( intermittency ) , @xmath61 , as defined by eq . [ e : w : dm ] and discussed in the previous paragraph . the parameter @xmath58 was computed for all selected time intervals during two solar minima ( 1997 - 1998 , 2007 - 2008 ) and maximum ( 1999 - 2001 ) . the results are shown in figures  [ f : w : rad ] and  [ f : w : lat ] . + figure  [ f : w : rad ] suggests a slow evolution of intermittency in the fast and slow solar wind with the distance @xmath62 from the sun . the fast wind multifractal spectra seem to exhibit a small decrease of intermittency with the distance as @xmath63 . the rate of this decrease has been determined by the least squares linear fit to the log - log values of @xmath58 against @xmath62 to all points both from solar minimum and maximum ( figure [ f : w : rad]a ) . for the slow solar wind ( figure [ f : w : rad]b ) , the evolution of intermittency with distance @xmath62 is more variable but anyway the decrease with distance with rate @xmath64 is observed anyway . + a latitudinal dependence of the multifractality of turbulence is also found , as illustrated in figure  [ f : w : lat ] . the fast solar wind at solar minimum exhibits a decrease in intermittency as the latitude increases ; the smallest values of multifractality are found near solar pole . another significant effect is the existence of a symmetry with respect to the ecliptic plane , observed mainly during solar minima ( 1997 - 1998 and 2007 - 2008 ) , which confirms previous observations @xcite . this is an indication that the fast polar solar wind has similar turbulent properties in the two hemispheres . we can not confirm this trend at solar maximum since the statistics of slow wind intervals in our database is rather poor and the solar wind is much more variable . + the results denoted in figures  [ f : w : rad ] and  [ f : w : lat ] by blue circles correspond to solar minimum conditions ( 1997 - 1998 ) , heliocentric distances of the order of @xmath65 au , and heliographic latitudes smaller than @xmath66 . for this data subset , we observe that the level of multifractality is slightly lower for fast than for slow solar wind , with mean @xmath58 about @xmath67 ( compared to @xmath68 for the slow wind ) . these results are consistent with previous studies performed in the ecliptic plane @xcite and with previous analyses of _ ulysses _ data for the years 1992 - 1997 @xcite , which postulate that intermittency in the slow solar wind is higher than for the fast wind . + results for the second minimum ( 2007 - 2008 ) , marked in figures  [ f : w : rad ] and  [ f : w : lat ] by black circles , seem to suggest somewhat different behavior . however , the data from this time interval have been collected for a wider range of latitudes , @xmath69 , and for a radial distance from the sun between @xmath70 and @xmath71 au . our hypothesis here is that the effect of change on solar wind multifractality with heliographic latitude may interfere with variations of this parameter observed with respect to the distance from the sun . the statistics of `` pure '' slow wind data is very poor ( only three intervals for 2007 - 2008 ) , correlated with the end of ulysses mission ; therefore , a comparison of the variation of multifractality between fast / slow wind and the solar cycle is not possible . + in general , in the case of _ ulysses _ data , it is difficult to separate latitudinal and radial evolution of intermittency . the calculated multifractality as function of both heliocentric distance and heliographic latitude for the fast and slow solar wind are summarized in figures [ f : w : fast_slow]a and [ f : w : fast_slow]b , respectively . one may see that fluctuations of the magnetic field strength both in the fast and slow solar wind exhibit increased intermittent features at small distances from the sun . moreover , once again we see that at small heliolatitudes and at distance of about @xmath72 au , the slow wind is more intermittent than the fast . further , analysis of results presented in figures  [ f : w : rad ] and  [ f : w : lat ] shows a rather similar range of the degree of intermittency determined for the solar maximum and minimum . in other words , we have less dependence on the phases of solar activity than on the type of the solar wind @xcite . + additionally , the multifractal analysis of _ ulysses _ magnetic field data presents a somewhat smaller degree of multifractality , @xmath73 , than those obtained previously for velocity fluctuations , @xmath74 @xcite . + in this letter , we have performed systematic studies of the multifractal and intermittent behavior of fluctuations of the interplanetary magnetic field strength for a broad range of heliographic latitudes , heliocentric distances , and for the maximum and minimum phase of the solar cycle . the obtained results can be summarized as follows : ( 1 ) the slow solar wind during solar minimum ( 1997 - 1998 ) , at distances of @xmath65 au and close to the equatorial plane , exhibits on average a higher level of intermittency than fast solar wind ; ( 2 ) the fast solar wind measured during the second minimum ( 2007 - 2008 ) , at distances of @xmath75 au and latitudes of @xmath76 reveals a decrease of intermittency ; ( 3 ) the level of intermittency during solar maximum present non - orderly behavior ; however , in many cases , we observe the same degree as those determined for data from solar minimum ; and ( 4 ) during solar minima ( 1997 - 1998 and 2007 - 2008 ) , similar intermittent properties of the fast polar solar wind in the two hemispheres are observed . + generally , at small distances from the sun , both in the slow and fast solar wind , we observe the highest degree of multifractality , suggesting that most of the intermittent events have a solar origin . moreover , the level of intermittency seems to decrease also with latitude . this effect can be related to the fact that at radial distances of the order of 1 au and less , the solar wind turbulence is still developing ( see , e.g. , the structure function analysis of @xcite ) ; thus , the solar wind structures that contribute to intermittency have not reached statistical stationarity . at larger radial distances and higher heliolatitudes , the turbulence is decaying as there are fewer  sources " of turbulence than at low latitudes , like , e.g. , velocity stream shears . this hypothesis could be tested by investigation of _ ulysses _ data not included in this study . another effect that could possibly contribute to the decrease of intermittency with the radial distance is the inherent expansion of the structures involved in the intermittent transfer of energy ; their statistics would therefore decrease within the ranges of scales considered in this study . possible tests of this hypothesis could be provided by numerical simulations with an expanding solar wind . + based on the investigations presented in this letter , we conclude that the evolution of `` pure '' slow and fast solar wind turbulence beyond the ecliptic plane is probably insufficient to maintain a high level of intermittency . while previous analyses of the radial evolution of solar wind intermittency were based on probability distribution functions and their moments , in this study we use the degree of multifractalality as a key descriptor . although the two approaches have similarities , they do not capture exactly the same aspects of solar wind fluctuations and scale behavior . in a future study , we will explore the consistency between the two approaches for the same data set . we would like to thank the swoops , swics , and mag instruments teams of the _ ulysses _ mission for providing data . this work was supported by the european community s seventh framework programme ( [ fp7/2007 - 2013 ] ) under grant agreement no . were additionally supported by the national science center , poland ( ncn ) , through grant 2014/15/b / st9/04782 . m.m.e . acknowledges support from the romanian ministry of national education , cncs , uefiscdi , project number pn - ii - id pce-2012 - 4 - 0418 , and from the interuniversity attraction poles programme initiated by the belgian science policy office ( iap p7/08 charm ) .                                          , h. , wawrzaszek , a. , macek , w.m . & chang , t. 2010 , in aip conf . ser . , 1216 , twelfth int . solar wind conf . , ed . m. maksimovic , k. issautier , n. meyer - vernet , m. moncuquet , & f. pantellini ( melville , ny : aip ) , 124
Abstract:

Answer: we study intermittency as a departure from self - similarity of the solar wind magnetic turbulence and investigate the evolution with the heliocentric distance and latitude . we use data from the _ ulysses _ spacecraft measured during two solar minima ( 1997 - 1998 and 2007 - 2008 ) and one solar maximum ( 1999 - 2001 ) . in particular , by modeling a multifractal spectrum , we revealed the intermittent character of turbulence in the small - scale fluctuations of the magnetic field embedded in the slow and fast solar wind . generally , at small distances from the sun , in both the slow and fast solar wind , we observe the high degree of multifractality ( intermittency ) that decreases somewhat slowly with distance and slowly with latitude . the obtained results seem to suggest that generally intermittency in the solar wind has a solar origin . however , the fast and slow streams , shocks , and other nonlinear interactions can only be considered as the drivers of the intermittent turbulence . it seems that analysis shows that turbulence beyond the ecliptic plane evolves too slowly to maintain the intermittency with the distance and latitude . moreover , we confirm that the multifractality and intermittency are at a lower level than in the ecliptic , as well as the existence of symmetry with respect to the ecliptic plane , suggesting that there are similar turbulent properties observed in the two hemispheres .