Instruction: Write an abstract for the following acadamic paper.

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

Paper: the rayleigh - taylor instability @xcite , @xcite , @xcite is a classical hydrodynamic instability @xcite with relevance in such diverse fields as plasma physics , astrophysics , meteorology , geophysics , inertial confinement fusion and granular media , for a review see , e.g. , @xcite . generically this instability develops if a layer of liquid is superimposed to an immiscible and less dense liquid such that the potential energy of the system can be reduced by interchanging the liquids . consequently , the initially plane interface between the liquids becomes unstable and the characteristic dimples and spikes develop resulting finally in a stable layering with the lighter fluid on top of the heavier one . a quantitative experimental investigation of the rayleigh - taylor instability requires reliable control of the initial condition . standard procedures like suddenly removed partitions between the fluids @xcite , @xcite or quickly turning the experimental cell upside down @xcite , @xcite clearly produce unpredictable initial perturbations . it is much more convenient to use some additional mechanism which first _ stabilizes _ the unstable layering of the liquids and may later be switched off instantaneously . it is well known that the rayleigh - taylor instability may be suppressed , e.g. , by vertical oscillation of the system @xcite , @xcite and by appropriate temperature gradients @xcite . the first mechanism is likely to induce uncontrolled initial surface deflections when stopped , in the second one it is difficult to abruptly switch off the stabilization .    in the present paper we investigate the possibility to stabilize a potentially rayleigh - taylor unstable system involving a magnetic fluid by external magnetic fields . we will show that for experimentally relevant parameter values moderate fields strengths which can easily be switched on and off are sufficient to achieve the desired stabilization . magnetic fluids are suspensions of ferromagnetic nano - particles in carrier liquids with the hydrodynamic properties of newtonian fluids and the magnetic properties of super - paramagnets @xcite , @xcite . it is well known that a magnetic field parallel to the plane interface between a ferrofluid and a non - magnetic fluid suppresses interface deflections with wave vector in the direction of the field @xcite . this may be used to stabilize the rayleigh - taylor instability in two - dimensional situations where the interface is line , as e.  g. in a hele - shaw cell . in the more natural three - dimensional setting we are interested in here a static magnetic field parallel to the undisturbed interface is not sufficient to stabilize the flat interface since perturbation with wave vectors _ perpendicular _ to the magnetic field will still grow as in the absence of the field . we therefore propose to use a spatially homogeneous magnetic field _ rotating _ in the plane of the undisturbed interface and determine appropriate values of the field amplitude and rotation frequency . an alternative possibility is to use a static _ inhomogeneous _ magnetic field with the magnetic force counterbalancing gravity @xcite . this method was used in @xcite to investigate the 2-d rayleigh - taylor instability in a hele - shaw cell .    before embarking upon the detailed analysis we would like to mention three characteristic features of our method . firstly , we will show that a rotating magnetic field is unable to suppress _ all _ possible unstable modes of the system . in fact it can only stabilize surface deflections with wavenumber modulus larger than some threshold value . perturbations with very long wavelength are , however , not a serious problem in real experiments because these are suppressed automatically by the finite geometry of the sample . secondly , it is well - known that in analogy with the faraday instability an oscillating magnetic field will induce _ new _ instabilities at wave numbers which without field were stable @xcite , @xcite . in our analysis we keep track of these unstable modes and determine the magnetic field strength such that no new instabilities may occur . for suppressing these new modes viscous losses in the liquids will be decisive which is the reason why the viscosities of the two liquids will be consistently taken into account in the analysis . finally , due to the dispersed magnetic grains ferrofluids have usually comparatively high densities . we therefore specialize to the case in which the upper , heavier layer is formed by the magnetic fluid . this should be the typical situation in experiments . nevertheless a similar analysis with analogous results is possible for the reverse situation with the ferrofluid at the bottom of the system superimposed by an even denser non - magnetic liquid . the paper is organized as follows . in section ii we collect the basic equations and boundary conditions . in section iii we linearize these equations around the reference state of a plane interface between the liquids . section iv contains the floquet theory to determine the boundaries separating stable from unstable regions in the parameter plane . after shortly discussing two approximate treatments of the fluid viscosities in section v we present the results of our analysis in section vi . finally section vii contains some discussion . we consider a ferrofluid with density @xmath0 superimposed on a non - magnetic fluid of lower density @xmath1 , see fig . [ rosensweig ] . both layers are assumed to be infinite in horizontal as well as in vertical direction . the densities and the respective viscosities @xmath2 and @xmath3 are taken to be constant . the liquids are immiscible and the interface between them is parametrized by @xmath4 . we will study the stability of a flat interface which we take as the @xmath5-@xmath6-plane of our coordinate system , the undisturbed interface is hence given by @xmath7 . in the absence of a magnetic field this situation is unstable due to the rayleigh - taylor instability @xcite .    in the presence of an external magnetic field @xmath8 the magnetic fluid builds up a magnetization @xmath9 which is assumed to be a linear function of the field , @xmath10 , where @xmath11 denotes the susceptibility related to the relative permeability by @xmath12 . both liquids are subject to the homogeneous gravitational field acting in negative @xmath13-direction and to the interface tension @xmath14 acting at their common interface . the magnetic fluid is additionally influenced by the magnetic force density @xmath15 resulting from the externally imposed spatially homogeneous magnetic field @xmath16 rotating with constant angular frequency @xmath17 in the @xmath5-@xmath6-plane . , density @xmath0 and viscosity @xmath3 is superimposed on a non - magnetic fluid with density @xmath1 and viscosity @xmath2 . the normal vector @xmath18 on the interface @xmath19 points into the non - magnetic liquid . the vector @xmath20 denotes the gravitational acceleration , @xmath14 is the interface tension . the system is subjected to a horizontal rotating magnetic field @xmath21 with angular frequency @xmath17 . ] the time evolution of the system is governed by the following set of equations . the incompressibility of both liquids gives rise to the continuity equations for the velocity fields @xmath22 @xmath23 with @xmath24 where here and in the following the lower ( i.e. non - magnetic ) fluid parameters are denoted with superscript @xmath25 and the higher ( magnetic ) one with superscript @xmath26 . the hydrodynamic equations of motion are the navier - stokes equations @xmath27 with @xmath28 denoting the acceleration due to gravity and the stress tensors @xmath29 given by @xmath30      & + b_{\ell}^{(j)}h_{m}^{(j ) } + \eta^{(j)}\left({\partial}_\ell v_{m}^{(j ) }        + { \partial}_mv_{\ell}^{(j)}\right ) \ ; . \end{split}\ ] ] here @xmath31 denotes the pressure in each liquid , and @xmath32 is the respective magnetic induction . note that also the stress tensor for the non - magnetic liquid contains contributions from the magnetic field which , however , are divergence free and therefore do not give rise to a force density in the lower fluid .    for values of @xmath17 relevant to the present investigation radiative effects are negligible and the magnetic field has to obey the magneto - static maxwell equations @xmath33 \nabla\times{\bf h}^{(j ) } & = 0    \end{split}\quad.\ ] ] in view of the second equation it is convenient to introduce scalar magnetic potentials @xmath34 according to @xmath35 . these potentials then fulfill the laplace equations @xmath36 the above equation have to be supplemented by appropriate boundary conditions . far from the interface the velocities must remain bounded , @xmath37 and the magnetic field must be equal to the externally imposed field , @xmath38 to formulate the boundary conditions at the interface @xmath4 we define the normal vector @xmath18 by @xmath39 the normal component of the stress tensor has to fulfill @xmath40 where @xmath41 is the local curvature of the interface and the symbol @xmath42 denotes here and in the following the difference in the value of the respective quantity slightly above and slightly below the interface . the tangential components of the stress tensor have to be continuous , @xmath43 for all vectors @xmath44 perpendicular to @xmath18 . the motion of the interface is related to the velocity fields in the liquids by the kinematic condition @xmath45 finally , at the interface the normal component of @xmath46 and the tangential component of @xmath8 have to be continuous which gives rise to the following boundary conditions for the magnetic potentials at the interface @xmath47      { \partial}_n\phi^{(1 ) } & = \mu_r{\partial}_n\phi^{(2 ) }    \end{split } \quad .\ ] ] the main purpose of the present work is to investigate whether the rayleigh - taylor instability can be suppressed with the help of a rotating magnetic field . we will hence study the linear stability of the reference state with a flat interface , @xmath48 , in dependence on the magnetic field strength @xmath49 and the angular frequency @xmath17 . the reference solution of the basic equations is given by @xmath50     & \phi_0 = -h_0\left(\cos(\omega t)x + \sin(\omega t)y\right )    \end{split}\qquad .\ ] ] to investigate its stability we introduce as usual small perturbations @xmath51      & \phi^{(j)}=\phi_0+\phi^{(j ) }    \end{split } \ ; , \ ] ] and linearize the basic equations in these perturbations as well as in the interface deflection @xmath19 . we will denote the components of the perturbed velocity vectors by @xmath52 . it is convenient to introduce dimensionless quantities according to @xmath53      t & \rightarrow      \left(\frac{\sigma}{g^3\triangle\rho}\right)^{1/4}t\\[1ex ]      p^{(j ) } & \rightarrow\sqrt{\triangle\rho g \sigma}\;p^{(j)}\\[1ex ]      h^2 & \rightarrow \left(\frac{2}{\mu_0}\frac{\chi +          2}{\chi^2}\sqrt{\triangle\rho         g\sigma}\right ) h^2\\[1ex ]      \rho^{(j ) } & \rightarrow\triangle\rho\;\rho^{(j)}\\[1ex ]      \eta^{(j ) } & \rightarrow\left (      \frac{\triangle\rho\sigma^3}{g}\right)^{1/4}\;\eta^{(j)}\ ; ,    \end{split}\ ] ] where @xmath54 is the critical wave number of the rayleigh - taylor instability for @xmath55 . the linearized set of basic equations ( [ konti],[nse],[laplace ] ) then reads @xmath56     \label{nselin }      \rho^{(j)}{\partial}_t\delta{\bf v}^{(j ) } & = -\nabla\delta p^{(j ) }      + \eta^{(j)}\nabla^2\delta{\bf v}^{(j)}\\              & - 2\frac{\chi+2}{\chi}\nabla\left({\bf h}_0^{(j)}\cdot               \nabla\phi^{(j)}\right ) \\[1ex ]      \label{laplacelin }    \nabla^2\phi^{(i ) } & = 0 \ ; .\end{aligned}\ ] ] in order to eliminated the pressure it is convenient to consider the @xmath13-component of the curl curl of eq.([nselin ] ) which is of the form @xmath57 where we have introduced the kinematic viscosities @xmath58 .    from the boundary conditions ( [ vinfinity ] ) and ( [ maginfinity ] ) we find @xmath59 and @xmath60 the boundary conditions at the interface simplify under linearization . generally we may replace the interface position @xmath4 by @xmath61 to linear order in @xmath19 . therefore the symbol @xmath62 has now the more specific meaning @xmath63 . from ( [ kinematic ] ) we then get @xmath64 implying @xmath65 moreover , the continuity of the flow field @xmath66 together with ( [ kontilin ] ) gives rise to @xmath67 from ( [ tangentialtensor ] ) we find @xmath68 = 0\ ; .\ ] ] finally , linearization of ( [ normaltensor ] ) together with ( [ nselin ] ) yields @xmath69 = \\[1ex ]      & \quad\nabla_\perp^2\left(\left(1 + \nabla_\perp^2\right)\zeta +      2\frac{\chi+2}{\chi}{\bf h}_0\cdot\nabla\phi^{(2)}\big|_{z=0}\right ) \ ; . \end{split}\ ] ] where the horizontal laplace operator is defined by @xmath70 .    the magnetic boundary conditions ( [ magcon ] ) acquire the form @xmath71      & m_0\big(\cos(\omega t){\partial}_x\zeta + \sin(\omega t){\partial}_y\zeta\big ) +      { \partial}_z\big(\mu_r\phi^{(2)}-\phi^{(1)}\big)=0\ ; . \end{split}\ ] ]    to find a solution of the set of linearized equations ( [ laplacelin ] ) , ( [ ccnselin ] ) together with their boundary conditions we may exploit their translational invariance and have to keep in mind their explicit time dependence induced by the second boundary condition ( [ magconlin ] ) for the magnetic field problem . an appropriate ansatz is therefore given by @xmath72       w^{(j)}(x , y , z , t )    \end{array}\right)=    \left (    \begin{array}{c }      \hat\zeta(t)\\[1ex ] \hat w^{(j)}(z , t )    \end{array}\right)e^{i(k_xx + k_yy)}\ ] ] and @xmath73 with the abbreviation @xmath74 eq . ( [ ccnselin ] ) acquires the form @xmath75 moreover the ansatzes ( [ ansatzes2 ] ) and ( [ ansatzes3 ] ) already fulfill ( [ laplacelin ] ) , ( [ maginfinitylin ] ) and the first of the boundary conditions ( [ magconlin ] ) . the second one yields @xmath76 which gives rise to @xmath77 the boundary conditions ( [ winfinity],[kinematiclin]-[nse+normaltensorlin ] ) assume the form @xmath78 \label{kinematiclinf }    { \partial}_t\hat \zeta & = \hat w^{(j)}\big|_{z=0}\\    \label{velocity1f }    \triangle\ , \hat w & = 0\\[1ex ]    \label{velocity2f }    \triangle\ , { \partial}_z \hat w & = 0\\[1ex ] \label{tangentialtensorlinf }   \triangle\big[\eta \left(k^2 + { \partial}_z^2\right)\hat w\big ] & = 0\end{aligned}\ ] ] and , using also ( [ h1 ] ) , @xmath79=\\[1ex ]            & \big(-1 + h_0 ^ 2(1+\cos(2\omega t))k +            k^2\big)k^2\hat\zeta\ ; .    \end{split}\ ] ] we now invoke floquet theory @xcite to solve this system of linear differential equations with time periodic boundary conditions for the amplitudes @xmath80 and @xmath81 . in order to analyze the stability of the flat interface we employ the following floquet ansatz for the time dependence of the interface perturbation amplitude @xmath82 and the @xmath13-components of the velocity @xmath83 : @xmath84 where @xmath85 is the floquet exponent . here @xmath86 is a real number and negative @xmath86 describes stable situation whereas positive @xmath86 signals an instability of the reference state . the imaginary part @xmath87 of the floquet exponent is either zero or one and distinguishes between harmonic ( @xmath88 ) and subharmonic ( @xmath89 ) response of the system @xcite . plugging ( [ floquetansatz ] ) into ( [ nselinf ] ) we find @xmath90 where @xmath91 eq.([dglw ] ) has the solution @xmath92 where the constants @xmath93 can be determined with the help of the boundary conditions ( [ winfinity2]-[tangentialtensorlinf ] ) . as a result the amplitude of the @xmath13-component of the velocity may be expressed in terms of the interface amplitude @xmath94 according to @xmath95       \tilde w_n^{(1 ) } & = \left(\tilde a_ne^{kz } + \tilde        c_ne^{q_n^{(1)}z}\right)\tilde\zeta_n \ ; ,    \end{split}\ ] ] where    @xmath96        \tilde b_n & = \frac{\triangle\eta k^3 + \left(\eta^{(1)}q_n^{(1 ) }          + \eta^{(2)}q_n^{(2)}\right)k^2 +        \left(\eta^{(1)}\left(k+q_n^{(1)}\right ) +           \eta^{(2)}\left(k+q_n^{(2)}\right)\right)q_n^{(2)2 } +          2\eta^{(1)}q_n^{(1)}q_n^{(2)}k}{\eta^{(2)}\left(k + q_n^{(2)}\right ) +          \eta^{(1)}\left(k + q_n^{(1)}\right)}\nu^{(2)}\\[1ex ]        \tilde c_n & = \frac{\left(\eta^{(1)}k            + \eta^{(2)}q_n^{(2)}\right)\left(q_n^{(2)2 } - k^2\right ) }         { \left(k - q_n^{(1)}\right)\left(\eta^{(2)}\left(k+q_n^{(2)}\right ) +          \eta^{(1)}\left(k+q_n^{(1)}\right)\right ) } 2k\nu^{(2)}\\[1ex ]        \tilde d_n & = -\frac{\left(\eta^{(1)}q_n^{(1 ) }            + \eta^{(2)}k\right)\left(k + q_n^{(2)}\right ) }         { \eta^{(2)}\left(k + q_n^{(2)}\right ) + \eta^{(1)}\left(k +            q_n^{(1)}\right)}2k\nu^{(2 ) } \quad .      \end{split}\ ] ]    finally , using these results in ( [ nse+normaltensorlinf ] ) we find a relation of the form @xmath97\right\ } k^3\ , e^{(\alpha+i(\beta      + 2n))\,\omega t}=0 \ ; .\ ] ] since this equation has to hold for all values of @xmath98 all coefficients in the curly brackets must vanish separately . we therefore end up with an infinite homogeneous system of linear equations for the amplitudes @xmath99 in which the off - diagonal terms arise due to the time dependence in ( [ nse+normaltensorlinf ] ) . nontrivial solutions for the @xmath99 require that the determinant of the coefficient matrix vanishes which yields the desired relation between the parameters of the problem , @xmath49 , @xmath100 and @xmath86 . in the present investigation we are mainly interested in the stability boundaries in the parameter plane . we therefore specialize to the case @xmath101 and find for the coefficients @xmath102 in ( [ finaleq ] ) @xmath103+\\[1ex ] & \frac{1}{k}\big[-            3\left(\eta^{(2)}(k\tilde b +            q_n^{(2)}\tilde d ) + \eta^{(1)}(k\tilde a + q_n^{(1)}\tilde            c)\right ) + 1 - k^2\big]\ ; . \end{split}\ ] ]    to exploit the solvability condition for a numerical determination of the stability boundaries we have to truncate the infinite system of linear equations at some finite value @xmath104 of @xmath105 . comparing the results for different values of @xmath104 the accuracy of the procedure may be estimated . for the results presented in section vi we have used @xmath106 , i.e. we have included 39 terms , @xmath107 . before presenting explicit results of our analysis for experimentally relevant parameter values it is instructive to consider two limiting cases for which alternative approaches are available . let us first discuss the situation of ideal liquids , @xmath108 . using ( [ nselinf ] ) , ( [ kinematiclinf ] ) and ( [ velocity1f ] ) we may then express @xmath83 in terms of @xmath82 . plugging the result into the boundary condition ( [ nse+normaltensorlinf ] ) we obtain the following mathieu equation for the amplitude of the surface deflection @xmath109 : @xmath110       & + h_0 ^ 2 k^2\cos(2\omega t)\big)=0\ ; . \end{split}\ ] ] from the standard stability chart of the mathieu equation @xcite we are now able to determine the threshold for the amplitude @xmath49 of the external field necessary to stabilize interface deflections with wavenumber modulus @xmath100 . however , since most ferrofluids are rather viscous this theory will not adequately describe the experimental situation . it is possible to approximately incorporate the influence of viscosity by assuming that the dominant contribution to viscous damping originates far from the interface in the bulk of the fluids where the flow field is identical to the one of ideal liquids @xcite . one may then derive a damped mathieu equation for the amplitudes of the interface deflection of the form @xmath111       & + h_0 ^ 2 k^2\cos(2\omega t)\big)=0 \ ; ,    \end{split}\ ] ] where the damping constant @xmath112 is given by @xmath113 since the damped mathieu equation may be mapped on the undamped one @xcite we may again employ the stability chart of the mathieu equation to discuss the stabilization of a surface deflection mode with wave vector modulus @xmath100 in the presence of small damping . in the following section we compare the results of these approximate estimates with those of our complete treatment . stability chart in the plane spanned by the dimensionless wavenumber @xmath100 and the magnetic field amplitude @xmath49 for the flat interface between a ferrofluid and a non - magnetic fluid in a rotating magnetic field . the dimensionless angular frequency of the magnetic field is @xmath114 corresponding to @xmath115 , the remaining parameters are as given in the main text . the full lines separating the white regions of stability from the gray regions of unstable combinations derive from a numerical solution of ( [ finaleq ] ) . for comparison also the results of the inviscid theory building on ( [ mathieu ] ) and of the approximate treatment of viscosity related to ( [ mathieudamped ] ) are included as long and short dashed lines respectively . the threshold values of the magnetic field and wave number are @xmath116 and @xmath117 corresponding in physical units to @xmath118 and @xmath119 . ]    in this section we display detailed results of our analysis for a typical experimental combination of a ferrofluid and an immiscible non - magnetic fluid which has been used in a related experimental investigation @xcite . the fluid parameters are as follows : @xmath120 , @xmath121 , @xmath122 , @xmath123 , @xmath124 , and @xmath125 . for the capillary length @xmath126 we then obtain @xmath127 . the dimensionless magnetic field amplitude @xmath128 corresponds to a field of @xmath129 , @xmath130 corresponds to a field rotating with frequency @xmath131 . same as fig . [ hk-1 ] for an angular frequency of @xmath132 ( @xmath133 ) . due to the higher frequency the influence of the viscosity is stronger giving rise to the threshold values @xmath134 and @xmath135 corresponding to @xmath136 and @xmath137 . ]    in figs . [ hk-1 ] and [ hk-2 ] we show the regions of instability of the flat interface in the @xmath100-@xmath49 plane . for @xmath55 all perturbations are unstable for which the modulus of the wave vector is smaller than 1 ( in our dimensionless units , cf . ( [ dimlos ] ) ) , which is the well - known trademark of the rayleigh - taylor instability . increasing @xmath49 from zero the interval of unstable wave numbers shrinks and hence more and more long - wave perturbations may be stabilized . however , if @xmath49 gets larger than a threshold value @xmath138 the parametric excitation due to the time - dependent magnetic field gives rise to new instabilities at higher wave - numbers . since these additional unstable modes are clearly unwanted @xmath49 must remain below this threshold value @xmath138 . correspondingly there is a threshold @xmath139 for the wavenumber modulus such that perturbations with @xmath140 can not be stabilized with the help of the magnetic field . as we will detail in section vii these modes have to be stabilized by lateral boundary conditions . we note that with decreasing @xmath17 the tongues of instability move closer together and come nearer to the @xmath100-axis implying @xmath141 and @xmath142 for @xmath143 . it is clearly seen from the figures that the stability regions are strongly influenced by the viscosity of the liquids . in the inviscid theory the tongues of instability all reach the @xmath100-axis implying that any rotating magnetic field would induce new unstable modes at values of @xmath100 that were stable in the absence of the field . therefore a complete suppression of the rayleigh - taylor instability would be impossible . it is also apparent that for realistic parameter combinations the phenomenological inclusion of viscosity in the theoretical description as discussed in the previous section may give results which significantly differ from the complete theory . this is similar to the analysis of the faraday instability performed in @xcite .    figs . [ ohc ] and [ omegakc ] display the dependence of the threshold values @xmath138 and @xmath139 on the angular frequency @xmath17 of the field . clearly @xmath138 increases and @xmath139 decreases with increasing @xmath17 as exemplified also by a comparison between fig . [ hk-1 ] and [ hk-2 ] . for the parameters considered an increase in @xmath17 beyond @xmath144 does not significantly reduce @xmath139 any more . threshold value @xmath138 of the magnetic field amplitude @xmath49 as function of the angular frequency @xmath17 of the field for the parameters given in the main text . the displayed interval @xmath145 corresponds in physical units to @xmath146 . the inset shows a blow - up of the steep increase of @xmath138 for small values of @xmath17 . ]     similarly to fig . [ ohc ] the threshold value @xmath139 of the wavenumber modulus @xmath100 is shown as function of the angular frequency @xmath17 of the field . ]    finally , fig . [ hckc ] combines figs . [ ohc ] and [ omegakc ] and shows the relation between the two threshold values @xmath138 and @xmath139 . threshold value @xmath139 of the wave number @xmath100 versus threshold value @xmath138 of the magnetic field amplitude @xmath49 , again for the special set of parameters given in the main text . ] in the present paper we have investigated the possibility to stabilize a layering of a ferrofluid and a non - magnetic fluid which were potentially unstable due to the rayleigh - taylor instability by a spatially homogeneous magnetic field rotating in the plane of the undisturbed interface . special emphasis was put on an exact treatment of the influence of the viscosities by starting from the complete set of navier - stokes equations for both liquids . our results show that this approach is for experimentally relevant parameter values superior to both the inviscid theory and to a standard phenomenological procedure to include viscous effects using the inviscid flow field . the trademark of the rayleigh - taylor instability is a band of unstable wave numbers extending from @xmath147 up to a threshold value @xmath139 which in the absence of magnetic effects is given by the capillary wavelength @xmath148 . the main result of the present investigation is that @xmath139 may be reduced roughly by a factor of ten with the help of a rotating magnetic field of experimentally easily accessible amplitude and frequency . as expected the stabilization works best for ferrofluids with high susceptibility @xmath11 which , however , have also high densities increasing in turn @xmath149 .    in order to provide a clean initial condition for an experimental study of the onset of the rayleigh - taylor instability one has also to stabilize the modes with @xmath150 . one way to accomplish this suppression may be to use the boundary condition of a finite geometry , i.e. by pinning the contact line between the liquids at the boundary of the sample . in this way all long wave - number modes with @xmath151 are stabilized . here @xmath152 is determined by the linear extension @xmath153 of the sample and roughly given by @xmath154 . modes with @xmath155 are suppressed by surface tension . if one is able to temporarily stabilize the remaining modes by the rotating magnetic field , i.e. if one is able to realize @xmath156 the flat interface is stable . switching off the magnetic field at a given time all modes with @xmath157 will become unstable . since it is easily possible to realize values of @xmath139 significantly smaller than @xmath158 , the wavenumber with largest growth rate in the absence of the field , the ensuing rayleigh - taylor instability should closely resemble the case without lateral boundary conditions . to be precise it should be emphasized that our theoretical analysis is for infinite layers only and does not take into account the influence of lateral boundary conditions . however , the relevant values of @xmath138 and @xmath139 will only marginally be modified .    for an order of magnitude estimate let us consider a cylindrical vessel of diameter @xmath159 . pinning the contact line at the boundary the instability of modes with dimensionless wavenumber @xmath160 will be suppressed . on the other hand a rotating magnetic field with amplitude @xmath161 and frequency @xmath162 realizes @xmath163 ( cf . [ hk-2 ] ) . switching off the magnetic field all modes with @xmath157 will become unstable . for the above example this includes the first 8 cylindrical modes which should allow a rather accurate study of the rayleigh - taylor instability . we hope that our theoretical study may stimulate experimental work along these lines . 99 lord rayleigh , scientific papers ( cambridge university press , cambridge , england , 1900 ) * 2 * , pp . 200 - 207 g. i. taylor , proc . london , ser . a * 201 * , 192 ( 1950 ) d. j. lewis , proc . london , ser . a * 202 * , 81 ( 1950 ) s. chandrasekhar , _ hydrodynamic and hydromagnetic stability _ , ( oxford university press 1961 ) d. h. sharp , physics * d12 * , 3 ( 1984 ) p. f. linden and j. m. redondo , phys . fluids * a3 * , 1269 ( 1991 ) s. b. dalziel , p. f. linden , and d. l. youngs , j. fluid . mech . * 399 * , 1 ( 1999 ) m. s. plesset and c. g. whipple , phys . fluids * 17 * , 1 ( 1974 ) a. lange , m. schrter , m. a. scherer , and i. rehberg , eur . phys . j. b * 4 * , 475 ( 1998 ) g. h. wolf , prl * 24 * , 9 ( 1970 ) g. h. wolf , z. physik a * 227 * , 291 ( 1969 ) j. m. burgess , a. juel , w. d. mccormick , j. b. swift , and h. l. swinney , prl * 86 * , 1203 ( 2001 ) r. e. rosensweig , _ ferrohydrodynamics _ , cambridge university press , cambridge , ( 1985 ) m. d. cowley and r. e. rosensweig , j. fluid mech . * 30 * , 671 ( 1967 ) . g. pacitto , c. flament , j .- c . bacri , and m. widom , pre * 62 * , 7941 ( 2000 ) p. carls , z. huang , g. carbone , and c. rosenblatt , prl * 96 * , 104501 ( 2006 ) r. e. zelazo , and j. r. melcher , j. fluid mech . * 39 * , 1 ( 1969 ) d. w. jordan , and p. smith , _ nonlinear ordinary differential equations _ , 3th edition , oxford university press , oxford ( 1999 ) k. kumar , and l. s. tuckerman , j. fluid . mech . * 279 * , 49 ( 1994 ) ; k. kumar , proc . lond . a * 452 * , 1113 1996 y. fautrelle and a. d. sneyd , j. fluid . mech . * 375 * , 65 ( 1998 ) t. mahr , and i. rehberg , europhys . lett . * 43 * , 23 ( 1998 ) m. abramowitz , i.a . stegun , _ handbook of mathematical functions _ , ( dover publications , inc . , new york , 1972 ) h. lamb , _ hydrodynamics _ , 6th edition , ( dover , new york , 1932 ) ,  348 l. d. landau and e. m. lifschitz , _ hydrodynamik _ , ( akademie - verlag , berlin , 1991 ) ,  25
Abstract:

Answer: the rayleigh - taylor instability of a magnetic fluid superimposed on a non - magnetic liquid of lower density may be suppressed with the help of a spatially homogeneous magnetic field rotating in the plane of the undisturbed interface . starting from the complete set of navier - stokes equations for both liquids a floquet analysis is performed which consistently takes into account the viscosities of the fluids . using experimentally relevant values of the parameters we suggest to use this stabilization mechanism to provide controlled initial conditions for an experimental investigation of the rayleigh - taylor instability .


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: in the standard calculations of @xmath4 oscillations in the medium [ 1,2 ] the interaction of particles @xmath5 and @xmath6 with the matter is described by the potentials @xmath7 ( potential model ) . @xmath8 is responsible for loss of @xmath6-particle intensity . we consider the @xmath0 transitions [ 3,4 ] in a medium followed by annihilation : @xmath9 where @xmath10 are the annihilation mesons . for the process ( 1 ) potential model is used as well [ 5 - 12 ] .    in [ 10,11 ] it was shown that one - particle ( potential ) model mentioned above does not describe the process ( 1 ) and thus total neutron - antineutron transition probability : the process ( 1 ) probability is @xmath11 ( see ( 15 ) ) , whereas the potential model gives @xmath12 , where @xmath13 is the annihilation width of @xmath14 in the medium ( see eq . ( 6 ) of ref . [ 5 ] or eq . ( 1.6 ) of ref . [ 9 ] or eq . ( 16 ) of present paper ] ) . in the potential model the effect of final state absorption ( annihilation ) acts in the _ opposite _ ( wrong ) direction , which tends to the additional _ suppression _ of the @xmath0 transition . since the annihilation is the main effect which defines the speed of process ( 1 ) , the potential model should be rejected . this is because the unitarity condition is used for the essentially non - unitary @xmath1-matrix [ 10,11 ] . the interaction hamiltonian contains the antineutron optical potential @xmath15 and @xmath16 plays a crucial role . the @xmath1-matrix should be _ unitary_.    more formally , the basic equation @xmath17 , @xmath18 follows from the unitarity condition @xmath19 . however in the potential model the @xmath1-matrix is _ essentially _ non - unitary @xmath20 , resulting in @xmath21 because @xmath22 is extremely small : @xmath23 [ 10,11 ] . the above - given basic equation is inapplicable in this case . the potential model was developed in [ 5 - 8 ] . in more recent papers the verious details of the model have been refined , in particular the parameters of optical potential . we do nt dwell on these papers since the heart of the problem is in the non - hermiticity of the optical potential . in [ 13 ] we have proposed the model based on the diagram technique for direct reactions which does not contain the non - hermitian operators . subsequently , this calculation was repeated in [ 14,15 ] . however , in [ 16 ] it was shown that this model is unsuitable because the @xmath0 transition takes place in the propagator which is wrong . the neutron of the nucleus is in the bound state and so the neutron line entering into the @xmath0 transition vertex should be the wave function , but not the propagator , as in the model based on the diagram technique . for the problem under study this fact is crucial . it leads to the cardinal error for the process in nuclei . the @xmath0 transitions in the medium and vacuum are not reproduced at all . if the neutron binding energy goes to zero , the result diverges ( see eqs . ( 18 ) and ( 19 ) of ref . [ 13 ] or eqs . ( 15 ) and ( 17 ) of ref . so we abandoned this model [ 16 ] ( for more details , see [ 17 ] ) .    in [ 10 ] the model which is free of drawbacks given above has been proposed ( model * b * in the notations of present paper ) . however , the consideration was schematic since our concern was only with the role of the final state absorption in principle . in sect . 2 this model as well as the model with bare propagator are studied in detail . the corrections to the models ( sect . 3 ) and time - dependence ( sect . 4 ) are considered as well . in sect . 5 we sum up the present state of the investigations of the @xmath0 transition problem .    also we want to attract the attention to the problem under consideration because of the following reasons : 1 ) for the lower limit on the free - space @xmath0 oscillation time the range of uncertainty is too wide ( see eq . 2 ) the problem involves a number of questions which are of independent interest . the basic material is given in sects . figure 1 and analysis of amplitudes @xmath24 and @xmath25 are of major importance . first of all we consider the antineutron annihilation in the medium . the annihilation amplitude @xmath24 is defined as @xmath26 here @xmath27 is the hamiltonian of the @xmath14-medium interaction , @xmath28 is the state of the medium containing the @xmath14 with the 4-momentum @xmath29 ; @xmath30 denotes the annihilation products , @xmath31 includes the normalization factors of the wave functions . the antineutron annihilation width @xmath13 is expressed through @xmath24 : @xmath32    for the hamiltonian @xmath27 we consider the model @xmath33 where @xmath34 is the effective annihilation hamiltonian in the second quantization representation , @xmath35 is the residual scalar field . the diagrams for the model ( 4 ) are shown in fig . the first diagram corresponds to the first order in @xmath34 and so on . consider now the process ( 1 ) . the qualitative process picture is as follows . the free - space @xmath0 transition comes from the exchange of higgs bosons with the mass @xmath36 gev [ 4 ] and so the subprocess of @xmath0 conversion is scarcely affected by a medium effects . > from the dynamical point of view this is a momentary process : @xmath37 s. the antineutron annihilates in a time @xmath38 . we deal with two - step process with the characteristic time @xmath39 . the neutron wave function is @xmath40 here @xmath41 is the neutron 4-momentum ; @xmath42 , where @xmath43 is the neutron potential . the interaction hamiltonian has the form @xmath44 @xmath45 here @xmath46 is the hamiltonian of @xmath0 conversion [ 6 ] , @xmath47 is a small parameter with @xmath48 , where @xmath49 is the free - space @xmath0 oscillation time ; @xmath50 . in the lowest order in @xmath46 the amplitude of process ( 1 ) is _ uniquely _ determined by the hamiltonian ( 6 ) : @xmath51 @xmath52 @xmath53 , @xmath54 . here @xmath55 is the antineutron propagator . the corresponding diagram is shown in fig . the annihilation amplitude @xmath24 is given by ( 2 ) , where @xmath56 . since @xmath24 contains all the @xmath14-medium interactions followed by annihilation including antineutron rescattering in the initial state , the antineutron propagator @xmath55 is bare . once the antineutron annihilation amplitude is defined by ( 2 ) , the expression for the process amplitude ( 8) _ rigorously follows _ from ( 6 ) . for the time being we do not go into the singularity @xmath57 . transition in the medium followed by annihilation . the antineutron annihilation is shown by a circle . * b * same as * a * but the antineutron propagator is dressed ( see text),title="fig : " ]    one can construct the model with the dressed propagator . we include the scalar field @xmath35 in the antineutron green function @xmath58 @xmath59 , where @xmath60 is the antineutron self - energy . the process amplitude is @xmath61 @xmath62 ( see fig . the block in the square braces shown in fig.1 corresponds to the vertex function @xmath25 . the models shown in figs . 2a and 2b we denote as the models * a * and * b * , respectively . in both models the antineutron propagators do nt contain the annihilation loops since the annihilation is taken into account in the amplitudes @xmath24 and @xmath25 ; the interaction hamiltonians @xmath63 and unperturbed hamiltonians are the same . if @xmath64 , the model * b * goes into model * a*. in this sense the model * a * is the limiting case of the model * b*.    we consider the model * b*. for the process width @xmath65 one obtains @xmath66 @xmath67 where @xmath68 is the annihilation width of @xmath14 calculated through the @xmath25 ( and not @xmath24 ) . the normalization multiplier @xmath69 is the same for @xmath65 and @xmath68 . the vertex function @xmath25 is unknown . ( we recall the antineutron annihilation width @xmath13 is expressed through the amplitude @xmath24 . ) for the estimation we put @xmath70 this is an uncontrollable approximation . the time - dependence is determined by the exponential decay law : @xmath71 equation ( 15 ) illustrates the result sensitivity to the value of parameter @xmath60 .    on the other hand , for @xmath0 transitions in nuclear matter the potential model gives the inverse @xmath13-dependence [ 6 - 11 ] @xmath72 where @xmath15 is the antineutron optical potential . the wrong @xmath13-dependence is a direct consequence of the inapplicability of the model based on optical potential for the calculation of the total process probability [ 11 ] . ( the above - mentioned model describes the probability of finding an antineutron only . )    comparing with ( 15 ) , one obtains @xmath73 where the values @xmath74 mev and @xmath75 mev have been used . the parameter @xmath60 is uncertain . we have put @xmath76 mev only for estimation . the model * b * leads to an increase of the @xmath0 transition probability . the lower limit on the free - space @xmath0 oscillation time @xmath77 increases as well . let @xmath77 and @xmath78 be the lower limits on the free - space @xmath0 oscillation time obtained by means of eqs . ( 15 ) and ( 16 ) , respectively ; @xmath79 is the oscillation time of neutron bound in a nucleus . the relationships between @xmath77 , @xmath78 and @xmath79 are @xmath80 @xmath81 where the well - known equation @xmath82 has been used . for estimation we take @xmath74 mev , @xmath75 mev and @xmath83 s [ 18 ] . ( the limit @xmath83 s was derived [ 18 ] through eq . ( 20 ) from experimental bound on the neutron lifetime in oxygen @xmath84 yr obtained by super - kamiokande collaboration [ 18 ] . ) equation ( 18 ) gives @xmath85    if @xmath64 , eq . ( 18 ) diverges : @xmath86 . this circumstance should be clarified ; otherwise the model under consideration can be rejected . as we will see later , the correct formulation of the problem ( on the finite time interval ) leads to finite result for the model * a * , which justifies our approach . we return to the model shown in fig . we use the basis @xmath87 . the results do not depend on the basis . a main part of existing calculations have been done in @xmath88 representation . the physics of the problem is in the hamiltonian . the transition to the basis of stationary states is a formal step . it has a sense only in the case of the potential model @xmath89const . , when the hamiltonian of @xmath14-medium interaction is replaced by the effective mass @xmath90 because the hermitian hamiltonian of interaction of the stationary states with the medium is unknown . since we work beyond the potential model , the procedure of diagonalization of mass matrix is unrelated to our problem .    the amplitude ( 8) diverges @xmath91 ( see also eq . ( 21 ) of ref . [ 16 ] . ) these are infrared singularities conditioned by zero momentum transfer in the @xmath0 transition vertex . ( in the model * b * the effective momentum transfer @xmath92 takes place . )    for solving the problem the field - theoretical approach with finite time interval [ 19 ] is used . it is infrared free . if @xmath93 , the approach with finite time interval reproduces all the results on the particle oscillations ( see sect . 5.2 of ref . this is the test of above - mentioned approach . however , our purpose is to describe the process ( 1 ) by means of hermitian hamiltonian because the potential model describes the absorption wrongly .    for the model * a * the process ( 1 ) probability was found to be [ 12,16 ] @xmath94 where @xmath95 is the free - space @xmath0 transition probability . owing to annihilation channel , @xmath96 is practically equal to the free - space @xmath0 transition probability . if @xmath97 , eq . ( 23 ) diverges just as the modulus ( 22 ) squared does . if @xmath64 , eq . ( 15 ) diverges quadratically as well . the explanation of the @xmath98-dependence is simple . the process shown in fig . 2a represents two consecutive subprocesses . the speed and probability of the whole process are defined by those of the slower subprocess . if @xmath99 , the annihilation can be considered instantaneous . so , the probability of process ( 1 ) is defined by the speed of the @xmath0 transition : @xmath100 .    distribution ( 23 ) leads to very strong restriction on the free - space @xmath0 oscillation time [ 12,16 ] : @xmath101    if @xmath64 , @xmath102 rises quadratically . so @xmath77 and @xmath103 can be considered as the estimations from below and above , respectively . we show that for the @xmath0 transition in medium the corrections to the models and additional baryon - number - violating processes ( see fig . 3 ) can not essentially change the results . first of all we consider the incoherent contribution of the diagrams 3 . in fig . 3a a meson is radiated before the @xmath0 transition . the interaction hamiltonian has the form @xmath104 in the following the background neutron potential is omitted . the neutron wave function is given by ( 5 ) , were @xmath105 and @xmath106 . for the process amplitude @xmath107 one obtains @xmath108 @xmath109 where @xmath110 is the 4-momentum of meson radiated , @xmath111 is the amplitude of antineutron annihilation in the medium in the @xmath112 mesons . as with model * a * , the antineutron propagator @xmath113 is _ bare _ ; the @xmath14 self - energy @xmath114 . ( the same is true for figs . 3b-3d . ) if @xmath115 , the amplitude @xmath107 increases since @xmath116 , @xmath117 ( the limiting transition @xmath115 for the diagram 3a is an imaginary procedure because in the vertex @xmath118 the real meson is escaped and so @xmath119 . ) the fact that the amplitude increases is essential for us because for fig . 2a @xmath120 . due to this @xmath57 and @xmath121 . let @xmath122 and @xmath123 be the widths corresponding to the fig . 3a and annihilation width of @xmath14 in the @xmath124 mesons , respectively ; @xmath125 . taking into account that @xmath123 is a smooth function of @xmath126 and summing over @xmath124 , it is easy to get the estimation : @xmath127    the time - dependence is determined by the exponential decay law : @xmath128 comparing with ( 15 ) we have : @xmath129 . so for the model * b * the contribution of diagram 3a is negligible . for the model * a * the contribution of diagram 3a is inessential as well . indeed , using eqs . ( 30 ) and ( 23 ) we get @xmath130 where we have put @xmath131 . consequently , if @xmath132 and @xmath133 ( see ( 23 ) ) then the contribution of diagram 3a is negligible . for the @xmath0 transition in nuclei these conditions are fulfilled since in this case @xmath134 mev and @xmath135 yr , where @xmath136 is the observation time in proton - decay type experiment [ 18 ] . in fact , it is suffice to hold condition ( 33 ) only because it is more strong .    in the calculations made above the free - space @xmath0 transition operator has been used . this is impulse approximation which is employed for nuclear @xmath137 decay , for instance . the simplest medium correction to the vertex ( or off - diagonal mass , or transition mass ) is shown in fig . 3b . in this event the replacement should be made : @xmath138 @xmath139 , where @xmath140 is the correction to @xmath47 produced by the diagram 3b . for the model * a * the limit becomes @xmath141 obviously , the @xmath142 can not change the order of magnitude of @xmath103 since the @xmath143 operator is essentially zero - range one . the free - space @xmath0 transition comes from the exchange of higgs bosons with the mass @xmath36 gev [ 4 ] . since @xmath144 ( @xmath145 is the mass of @xmath146-boson ) , the renormalization effects should not exceed those characteristic of nuclear @xmath137 decay which is less than 0.25 [ 20 ] . so the medium corrections to the vertex are inessential for us . the same is true for the model * b*.    consider now the baryon - number - violating decay @xmath147 [ 21 ] shown in fig . it leads to the same final state , as the processes depicted in figs . 2 , 3a and 3b . denoting @xmath148 , for the decay width @xmath149 one obtains @xmath150 @xmath151 . the parameter @xmath152 corresponding to the vertex @xmath147 is unknown and so no detailed calculation is possible . the baryon - number - violating conversion @xmath153 in the medium [ 21 ] shown in fig . 3d can not produce interference , since it contains @xmath154-meson in the final state . for the rest of the diagrams the significant interferences are unlikely because the final states in @xmath155 annihilation are very complicated configurations and persistent phase relations between different amplitudes can not be expected . this qualitative picture is confirmed by our calculations [ 22 ] for @xmath156-nuclear annihilation . it is easy to verify the following statement : if the incoherent contribution of the diagrams 3a-3c to the total nuclear annihilation width is taken into account , the lower limit on the free - space @xmath0 oscillation time @xmath2 becomes even better . to summarise , the contribution of diagrams 3 is inessential for us . the non - trivial circumstance is the quadratic time - dependence in the model * a * : @xmath157 . the heart of the problem is as follows . the processes depicted by the diagrams 2b and 3 are described by the exponential decay law . in the first vertex of these diagrams the momentum transfer ( figs . 3a-3c ) , or effective momentum transfer ( figs . 2b , 3d ) takes place . the diagram 2a contains the infrared divergence conditioned by zero momentum transfer in the @xmath0 transition vertex . this is unremovable perculiarity . this means that the standard @xmath1-matrix approach is inapplicable [ 12,16,19 ] . in such an event , the other surprises can be expected as well . from this standpoint a non - exponential behaviour comes as no surprise to us . it seems natural that for non - singular and singular diagrams the functional structure of the results is different , including the time - dependence . the opposite situation would be strange . the fact that for the processes with @xmath120 the @xmath1-matrix problem formulation @xmath158 is physically incorrect can be seen even from the limiting case @xmath159 : if @xmath160 ( see ( 6 ) ) , the solution is periodic . it is obtained by means of non - stationary equations of motion and not @xmath1-matrix theory . to reproduce the limiting case @xmath161 , i.e. the periodic solution , we have to use the approach with finite time interval .    if the problem is formulated on the interval @xmath162 , the decay width @xmath13 can not be introduced since @xmath163 , @xmath164 . this means that the standard calculation scheme should be completely revised . ( we would like to emphasise this fact . ) the direct calculation by means of evolution operator gives the distribution ( 23 ) . formally , the different time - dependence is due to @xmath110-dependence of amplitudes . we consider eq . ( 26 ) , for instance . if @xmath110 decreases , the amplitude @xmath107 increase ; in the limit @xmath115 it is singular ( see ( 28 ) ) . the point @xmath120 corresponds to realistic process shown in fig . the @xmath98-dependence of this process is the consequence of the zero momentum transfer . the more physical explanation of the @xmath98-dependence is as follows . in the hamiltonian ( 25 ) corresponding to fig . 3a we put @xmath165 . then the virtual decay @xmath166 takes place . the first vertex of the diagram 3a dictates the exponential decay law of the overall process shown in fig . similarly , in the hamiltonian ( 6 ) corresponding to fig . 2a , we put @xmath159 . then the free - space @xmath0 transition takes place which is quadratic in time : @xmath167 . the first vertex determines the time - dependence of the whole process at least for small @xmath13 . we also recall that even for proton decay the possibility of non - exponential behaviour is realistic [ 23 - 25 ] . in both models the antineutron propagators do nt contain the annihilation loops since the annihilation is taken into account in the amplitudes @xmath24 and @xmath25 . the alternative model containing the full in - medium propagator has been considered in [ 26 ] . the sole physical distinction between models * a * and * b * is the definition of antineutron annihilation amplitude ; or , similarly , the non - zero antineutron self - energy in the model * b * which is conditioned by residual scalar field . however , it leads to the fundamentally different results .    if @xmath64 , @xmath168 diverges quadratically . this circumstance should be clarified ; otherwise the model * b * can be rejected . the calculation in the framework of the model * a * gives the finite result , which justifies our approach from a conceptual point of view and consideration of the model * a * at least as the limiting case . in fact , the model * a * seems quite realistic in itself . indeed , we list the main drawbacks of the model * b*.      \2 ) the diagram 2b means that the annihilation is turned on upon forming of the self - energy part @xmath59 ( after multiple rescattering of @xmath14 ) . this is counter - intuitive since at low energies [ 27,28 ] @xmath170 where @xmath171 and @xmath172 are the cross sections of free - space @xmath155 annihilation and @xmath155 scattering , respectively . the inverse picture is in order : in the first stage of the @xmath14-medium interaction the annihilation occurs . this is obvious for the @xmath0 transitions in the gas . the model * a * reproduces the _ competition _ between scattering and annihilation in the intermediate state [ 26 ] . \3 ) the time - dependence is a more important characteristic of any process . it is common knowledge that the @xmath173-dependence of the process probability in the vacuum and medium is identical ( for example , exponential decay law ( 15 ) ) . in the model * a * the @xmath173-dependencies in the vacuum and medium coincide : @xmath174 and @xmath175 . the model * b * gives @xmath176 , whereas @xmath175 . there is no reason known why we have such a fundamental change . \4 ) if @xmath177 , the model * a * reproduces all the well - known results on particle oscillations in contrast to the model * b * ( see sect . 5.2 of ref . [ 12 ] ) . in other words , the model * a * reproduces all the results of potential model . ( recall that our purpose is to describe the process ( 1 ) by means of hermitian hamiltonian because the potential model describes the absorption wrongly . ) the model * a * is free of drawbacks given above . the physics of the model is absolutely standard . for instance , for the processes shown in fig . 3 the antineutron propagators are bare as well . the same is also true for the diagram technique for direct reactions [ 13 - 15 ] . besides , the process amplitude ( 8) is derived _ directly _ from the basic eqs . ( 6 ) and ( 2 ) which are obvious .    however , there is fundamental problem in the model * a * : the singularity of the amplitude ( 22 ) . the approach with finite time interval gives the finite result , which justifies the models * a * and * b * at least in principle . nevertheless , the time - dependence @xmath174 and limit ( 24 ) seem very unusual . the corresponding calculation contains too many new elements . due to this we view the results of the model * a * with certain caution . besides , due to the zero momentum transfer in the @xmath0-transition vertex , the model is extremely sensitive to the @xmath60 . the process under study is _ unstable_. the small change of antineutron self - energy @xmath178 , or , similarly , effctive momentum transfer in the @xmath0 transition vertex converts the model * a * to the model * b * : @xmath179 with @xmath180 . ( this is because the process amplitude is in the peculiar point ( see ( 22 ) ) owing to zero momentum transfer . for the processes with non - zero momentum transfer the result is little affected by small change of @xmath110 . ) although we do nt see the specific reasons for above - mentioned scenario , it must not be ruled out . this is a point of great nicety .    finally , the values @xmath181 and @xmath182 yr are interpreted as the estimations from below ( conservative limit ) and from above , respectively . so the realistic limit @xmath2 can be in the range @xmath183 recall that for the free - space @xmath4 oscillations the @xmath4 transition probability is extremely sensitive to the difference of masses @xmath184 as well . for the neutron in the bound state ( @xmath0 transition in finite nucleus ) the result is the same [ 29 ] , what is consistent with the qualitative process picture : two - step process ( dynamic @xmath0 conversion [ 4 ] , annihilation ) with the characteristic time @xmath185 . the fact that process amplitude is in the peculiar point ( unremovable perculiarity ) is the basic reason why the range ( 38 ) is very wide . the estimation from below @xmath181 exceeds the restriction given by the grenoble reactor experiment [ 30 ] by a factor of 14 and the lower limit given by potential model by a factor of 5 . at the same time the range of uncertainty of @xmath2 is too wide . further theoretical and experimental investigations are desirable .    99 m.l . good , phys . rev . * 106 * , 591 ( 1957 ) . e. d. commins and p. h. bucksbaum , _ weak interactions of leptons and quarks _ ( cambridge university press , 1983 ) . v. a. kuzmin , jetf lett . * 12 * , 228 ( 1970 ) . r. n. mohapatra and r. e. marshak , phys . * 44 * , 1316 ( 1980 ) . p. g. h. sandars , j. phys . g * 6 * , l161 ( 1980 ) . m. v. kazarnovsky _ _ , pisma zh . . fiz . * 32 * , 88 ( 1980 ) . k. g. chetyrkin , m. v. kazarnovsky , v. a. kuzmin and m. e. shaposhnikov , phys . b * 99 * , 358 ( 1981 ) . j. arafune and o. miyamura , prog . theor . phys . * 66 * , 661 ( 1981 ) . w. m. alberico _ et al . _ , nucl . phys . a * 523 * , 488 ( 1991 ) . v. i. nazaruk , mod . a * 21 * , 2189 ( 2006 ) . v. i. nazaruk , eur . j. a * 31 * , 177 ( 2007 ) . v. i. nazaruk , eur . j. c * 53 * , 573 ( 2008 ) . v. i. nazaruk , yad . fiz . * 56 * , 153 ( 1993 ) . l. a. kondratyuk , pisma zh . exsp . theor . fiz . * 64 * , 456 ( 1996 ) . v. kopeliovich and i. potashnikova , eur . j. c * 69 * , 591 ( 2010 ) . v. i. nazaruk , phys . c * 58 * , r1884 ( 1998 ) . v. i. nazaruk , arxiv : 1101.2310 [ hep - ph ] . k. s. ganezer , presented at the international workshop on search for baryon and lepton number violations , berkley , september 2007 ( http://inpa.lbl.gov/blnv/blnv.htm ) . v. i. nazaruk , phys . b * 337 * , 328 ( 1994 ) . b. buch and s. m. perez , phys . 50 * , 1975 ( 1983 ) . j. basecq and l. wolfenstein , nucl . b * 224 * , 21 ( 1983 ) . v. i. nazaruk , phys . b * 229 * , 348 ( 1989 ) . g. n. fleming , phys . b * 125 * , 187 ( 1983 ) . k. grotz and h. v. klapdor , phys . c * 30 * , 2098 ( 1984 ) . p. m. gopich and i. i. zaljubovsky , part . nuclei * 19 * , 785 ( 1988 ) . v. i. nazaruk , eur . j. a * 39 * , 249 ( 2009 ) . t. kalogeropouls and g. s. tzanakos , phys . d * 22 * , 2585 ( 1980 ) . g. s. mutchlev _ et al . _ , d * 38 * , 742 ( 1988 ) . v. i. nazaruk , arxiv : 1004.3192 [ hep - ph ] . m. baldo - ceolin _ et al . _ , z. phys . c * 63 * , 409 ( 1994 ) .
Abstract:

Answer: the new model of @xmath0 transitions in the medium based on unitary @xmath1-matrix is considered . the time - dependence and corrections to the model are studied . the lower limit on the free - space @xmath0 oscillation time @xmath2 is in the range @xmath3 . -0.6 in    11.30.fs ; 13.75.cs    keywords : diagram technique , infrared divergence , time - dependence    * e - mail : nazaruk@inr.ru


Instruction: Write an abstract for the following acadamic paper.

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

Paper: the secondary stars in cvs have been recognised to follow approximately the mass  radius relation of main sequence field stars [ echeverra , 1983 , patterson , 1984 , warner 1995 , smith & dhillon , 1998 ] . a detailed comparison is complicated by the lack of information both on the masses of the secondary stars in cvs and of the masses and radii of main sequence field stars . it is advantageous , therefore , to restrict the comparison to quantities readily observable in cvs , as the orbital period @xmath0 and the spectral type @xmath1 of the secondary star [ beuermann et al . , 1998 , henceforth bbkw98 , smith & dhillon , 1998 ] . this approach requires to construct an equivalent diagram for field stars predicting the orbital period of a cv at which a field star of given spectral type would fill the roche - lobe of the secondary . note that this prediction is almost independent of the mass of the white dwarf primary . figure 1 shows the spectral types of the secondaries in cvs with known orbital period . the data points are based on published optical / near - ir spectra quoted in the compilation of ritter & kolb [ 1998 ] with some recent spectral type determinations added [ e.g. , schwarz et al . , 1999 , smith et al . , 1999 , thomas et al . , 1999 , thorstensen et al . , 1999 ] . purely photometric spectral type assignments have not been included . we have also not included spectral types derived from ir spectra only .    constructing a spectral type period diagram for field stars requires knowledge of their masses and radii . unfortunately , both quantities are not generally well known . in this situation it is important to note that recent progress in the construction of stellar models [ baraffe et al . , 1998 , henceforth bcah98 ] combined with the _ nextgen _ model atmospheres of hauschildt et al . [ 1999 ] have led to a substantially improved definition of the lower main sequence [ leggett et al . , 1996 , henceforth l96 ] . in fact , as noted by bbkw98 , the theoretical and observational radii of m - stars [ l96 ] agree at the few % level for stars of given absolute @xmath2-band magnitude @xmath3 . using the observed radii , we have established a radius scale which allows to estimate the radii of field stars of given @xmath3 [ beuermann et al . , 1999 , henceforth bbh99 ] . for the present purpose , we adopt the same stellar models which accurately reproduce the radii to estimate the masses , too . specifically , we use the theoretical relation between mass and absolute @xmath2-band magnitude , @xmath4 [ bcah98 ] to calculate the above defined ` orbital period ' for a field star of known @xmath3 and @xmath1 . we do not follow the approach of clemens et al . [ 1998 ] who adopted the observational relation between mass and absolute magnitude of henry & mccarthy [ 1993 ] because this relation represents a mean for stars of different age and metallicity and we want to study how differences in metallicity affect the location of stars in the @xmath5 diagram . figure 2 shows the resulting @xmath5 diagram for field stars of spectral types k / m supplemented by the sun . the main sequence of stars with near - solar metallicity is delineated by eight young disk ( yd ) field stars from l96 ( @xmath6 ) , the yd binary yy gem ( @xmath7 ) , the sun ( @xmath8 ) , and further 56 stars ( @xmath9 ) with radii from [ bbh99 ] . the effects of decreasing metallicity are indicated by the binary cm dra ( @xmath7 ) and by four old disk stars ( @xmath10 ) and four halo stars ( @xmath11 ) from l96 . also shown in fig . 2 are the models curves from the bcah98 stellar models , namely the zams model for solar metallicity ( solid line ) , and models for stars aged @xmath12 yrs with 1/3 solar metallicity ( long dashes ) and 1/30 solar metallicity ( short dashes ) . for the model stars , the spectral type is assigned by converting the theoretical colour @xmath13 to @xmath1 [ bbkw98 ] . the solid curve is expected to agree with the locus of the late - type yd stars . the slight displacement is probably due to remaining errors in the transformations used . comparison with fig . 1 indicates that the ( zero age ) main sequence of low - mass field stars with near - solar metallicity coincides with the locus of cv secondaries of the earliest spectral types at any given value of @xmath0 . the upper left in fig . 1 is devoid of cvs and shows that none of the systems included in the figure contains a secondary with metallicity substantially lower than the sun ( larger than solar metallicities are not excluded but remain unproved at present ) . the low space density of popii cvs is in agreement with population studies by stehle et al . [ 1999 ] . the secondaries in many cvs with @xmath14h have a later spectral type than expected for zams secondary stars of near - solar metallicity . these secondaries are expanded over and have lower masses than roche - lobe filling zams stars . the two principal causes for this expansion are the loss of thermal equilibrium due to the on - going mass transfer and nuclear evolution prior to the onset of roche - lobe overflow . kolb & baraffe [ 1999 ] have computed corresponding evolutionary sequences which nicely explain the observed behaviour and of which first results were presented in bbkw98 . the paths included in fig . 2 refer to secondaries starting as zams stars of 1@xmath15 and evolving under mass loss rates of @xmath16@xmath15yr@xmath17 ( solid curve ) and @xmath18@xmath15yr@xmath17 ( dashed curve ) . the point on these paths where the secondary becomes fully convective is indicated by the dotted curve [ kolb & baraffe , 1999 ] . late spectral types in cvs with orbital periods between 3 and 5 h can be explained by this scenario . the evolutionary models suggest that whenever the secondary becomes convective it has reached @xmath19@xmath15 and @xmath20m4.5 . if angular momentum loss drops rapidly at this point the secondary enters the period gap and re - appears below the gap with the same mass and spectral type @xmath21m4.5 . the comparatively late spectral types in cvs with orbital periods larger than 5 h can not be explained by the loss of thermal equilibrium due to mass loss but are consistent with nuclear evolution of the secondary star prior to the onset of mass transfer . two evolutionary paths are included in fig . 2 for secondaries starting mass transfer at @xmath22@xmath15 with a rate of @xmath16@xmath15yr@xmath17 and a central hydrogen fraction which is reduced to 0.16 ( dot - dashed curve ) or is practically exhausted ( dot - dot - dashed curve ) [ bbkw98 , kolb & baraffe , 1999 ] . early evolutionary calculations of a similar type were performed by whyte & eggleton [ 1980 ] . the model calculations of kolb & baraffe [ 1999 ] indicate that stars driven out of thermal equilibrium stay at approximately the same effective temperature and spectral type as an undisturbed star of the same mass . their results allow to estimate the masses of the secondary stars and thereby their roche radii from the observed spectral types . this is a prerequisite for distance determinations of cvs using the surface brightness method [ bailey 1981 ] . the roche - lobe filling secondary stars experience gravity darkening and will not be of uniform surface brightness . more seriously and more difficult to model are the variations in surface brightness caused by irradiation . marsh [ 1990 ] showed that the strengths of the tio bands and the nai@xmath23 absorption - line doublet decrease on the hemisphere of the m5 secondary star in ht cas which faces the primary . naively , one might expect that heating of the secondary star produces an earlier spectral type and a corresponding increase in the tio band strength . observation indicates an opposite result , a decreased flux of the tio features on the illuminated hemisphere , which suggests that a major change in the atmospheric structure takes place as a result of heating . schwope et al . [ 1999 ] present a nice tomographic picture of the non - uniform appearance of the secondary star of qq vul in the nai lines . similar results have been obtained e.g. for am her [ southwell et al . , 1995 ] , illustrating the complications caused by the non - uniform appearance of the secondary star for dynamic mass determinations . this non - uniformity must be taken into account when applying the surface brightness method for distance measurements . barnes & evans [ 1976 ] found that the visual surface brightness of giants and supergiants ( the only stars for which directly measured radii are available , apart from the few eclipsing binaries ) is a function of colour only and independent of luminosity ( gravity ) . bbh99 derived the surface brightness for the m - dwarfs studied by l96 and demonstrated that it deviates from that of giants . the implied small gravity dependence is in perfect agreement with recent model calculations for dwarfs and giants [ bcah98 , hauschildt et al . , 1999 ] . the dwarf / giant difference reaches a maximum for @xmath24m0 ( @xmath25=1.7 , @xmath26=3.5 ) . bbh99 also showed that the surface brightness of dwarfs in the @xmath2-band depends on metallicity in agreement with the predictions of bcah98 . since the well - observed cvs have near - solar metallicities ( figs . 1 and 2 ) , consideration of the gravity and metallicity dependencies allows improved surface brightness - colour relations to be established which are valid , e.g. , for cvs of the galactic disk population . 3 shows the resulting surface brightness @xmath27 in the @xmath2-band as a function of colour @xmath26 which can be fit by the linear relationship given in the figure . this relation differs from the original relation of bailey [ 1981 ] which was widely used in cv research . bailey s calibration of @xmath27 depended on the barnes - evans relation for giants and is , therefore , expected to be low by some 0.3 - 0.4 mag near spectral type m0 or @xmath26@xmath28 . the systematic difference between the two relations is primarily due to the gravity dependence of @xmath27 . the often - cited near constancy of @xmath27 for m - dwarfs does not exist . the scatter in our data is substantially reduced over that in bailey s diagram , and in the similar one by ramseyer [ 1994 ] , because we considered only single dwarfs with near - solar metallicity and avoided colour transformations between different photometric systems . low - metallicity m - dwarfs have @xmath27 values lower by @xmath29 mag . there is a systematic uncertainty in the derived surface brightness of about 0.1 mag carried over from the remaining uncertainty in the radius scale [ bbh99 ] . the secure identification of cv secondary stars requires spectroscopic observations in the optical / near - ir spectral regions . m - stars are best detected by their pronounced tio bands in the red part of the optical spectrum which display a variation in band strength ratios with spectral type [ wade & horne , 1988 , marsh , 1990 ] . we have calibrated the _ absolute strength _ of the flux depression at 7165 relative to the quasi - continuum at 7500 vs. spectral type for a selection of field stars with near - solar metallicities and show the result in fig . the quantity @xmath30 is the tio flux depression expressed as a surface brightness , with @xmath31 the distance and @xmath32 the stellar radius . the spectral fluxes @xmath33 and @xmath34 represent averages over @xmath35 and @xmath36 , respectively . different from the surface brightness @xmath27 above , @xmath37 is given in physical units . a second order polynomial fit to the logarithm of @xmath37 is @xmath38 where @xmath39 is the spectral subtype for m - dwarfs , i.e. @xmath40 for an dm5.5 secondary , and @xmath41 for k7 , preceding m0 . ( 1 ) is applicable to all cvs which have secondaries with near - solar abundances . for metal - poor m - dwarfs of the old - disk and halo populations , the surface brightness @xmath37 is reduced over that of eq . ( 1 ) . in extreme subdwarfs the @xmath42 flux depression disappears and our definition of @xmath37 becomes meaningless . hence , eq . ( 1 ) needs revision should spectral types become available for the secondaries in pop ii cvs , e.g. , in globular clusters . @lclcccccc name & @xmath43 & @xmath1 & @xmath44 & @xmath2 & @xmath45 & @xmath46 & @xmath47 & @xmath48 +   + v834cen & 1.69 & m5.5 & 1.78 & @xmath49 & 0.13&@xmath50 & @xmath51 & + z cha & 1.79 & m5.5 & 1.27 & @xmath52 & 0.13&@xmath53 & @xmath54 & + am her & 3.09 & m4 + & 21.0 & @xmath55 & 0.27 & @xmath56&@xmath57&@xmath58 + u gem & 4.25 & m4 + & 37.0 & @xmath59 & 0.41&@xmath60 & @xmath61 & @xmath62 + [ 0.5ex ] + bl hyi & 1.89 & m5.5 & 1.13 &  & 0.16 & @xmath63 & & + rx0203 & 1.89 & m2.5 & 1.10 &  & 0.45 & @xmath64 & & + uz for & 2.11 & m4.5 & 0.89 &  & 0.20 & @xmath65 & & + tt ari & 3.30 & m3.5 & 2.20 &  & 0.30 & @xmath66 & & + ix vel & 4.66 & m2 ? &  & 10.7&0.50 & & @xmath67 & @xmath68 + ae aqr & 9.88 & k4 &  & @xmath69 & 0.50 & & @xmath70 & @xmath71 +   + ht cas & 1.77 & m5.4 & 0.98 & @xmath72 & 0.13 & @xmath73 & @xmath74 & + ss cyg & 6.60 & k4 &  & @xmath75 & 0.70 &  & @xmath76 & @xmath77 + v1309 ori & 7.98 & m0.5 & 1.09&@xmath78 & 0.45 & @xmath79 & @xmath80 & + the surface brightness method [ bailey , 1981 ] allows the measurement of distances of cvs if the @xmath2-band flux of the secondary star or its spectral flux in the tio@xmath81 band can be measured . to be sure , the secondaries in cvs display a non - uniform surface brightness and care should be taken if only a single spectrum or flux measurement are available . ideally , the full orbital modulation is needed in order to select the most appropriate view on the secondary star for comparison with the appropriate surface brightness of field stars . in general , the @xmath2-band fluxes are not seriously affected by irradiation , while spectral features as the tio bands and the nai absorption lines are greatly diminished in strength on the heated face of the secondary star [ e.g. marsh , 1990 , schwope et al . , 1999 ] . both methods have , therefore , advantages and disadvantages . the advantage of using the tio band strength lies in its easy identification in observed spectra . we have used both the @xmath2-band brightness of the secondary and its tio band strength to derive the distances of selected cvs and show some results in table 1 ( for a more detailed discussion see beuermann & weichhold , 1999 . note that the error in @xmath82 has not been considered and the distances scale as @xmath83 ) . ideally , both methods should yield the same distances which should agree within errors with the trigonometric parallax . this is , in fact , the case for of am her and u gem ( parallaxes by c. dahn , private communication , harrison et al . , 1999 ) . we have confidence also in the distances derived for v834 cen and z cha which have no trigonometric parallaxes yet . there are some hopefully clear cases which represent cvs with well - defined tio features . obtaining independent distance information for these systems , e.g. , from the @xmath2-band fluxes of the secondary stars is clearly desirable . there are also controversial cases , however , which indicate some of the pitfalls of the methods . for ht cas , the distance derived from infrared photometry ( berriman et al . , 1987 ) is probably wrong , as noted already by ramseyer [ 1994 ] . another important case is that of ss cyg which has @xmath84 in quiescence . if that ir flux is entirely due to the k4 secondary , as assumed by bailey [ 1981 ] , our calibration yields @xmath85pc . on the other hand , the secondary contributes only about @xmath86% of the visual flux which implies that it is has @xmath87 and raises the distance to @xmath88pc , consistent with the recent hst fgs parallax of 166 pc [ harrison et al . , 1999 ] . note that assuming a k5 secondary would reduce @xmath31 again to 133pc . finally , the 8-h am herculis binary v1309 ori has an ( evolved ) m0@xmath89m1 secondary and a tio distance of 745 pc . the optical flux and spectral type imply that the secondary accounts for most of the @xmath2-band flux which is inconsistent with the statement of harrop - allin et al . [ 1997 ] of a contribution as low as @xmath90% based on @xmath2-band spectroscopy of the illuminated face of the secondary . the discrepancy is likely due to problems associated with the interpretation of illuminated stellar atmospheres . we conclude that reliable distances of cvs can be obtained from both the @xmath2-magnitude and the tio band strength of the secondary if its flux contribution can be unequivocally determined , illumination effects are taken properly into account , and the viewing direction on to the secondary is known . measuring accurate distances to cvs allows to derive basic quantities as their absolute magnitudes and the mass transfer rates [ warner , 1987 ] . much of our understanding of the physics of cvs [ warner , 1995 ] is based on the tedious derivation of such quantities . i thank isabelle baraffe and ulrich kolb for many discussions and for allowing me to show their evolutionary paths in fig . 1 prior to publication , hans ritter for pointing out some errors in my earlier compilation of spectral types of secondary stars , and marc weichhold for the initial work on this subject . baraffe i. , chabrier g. , allard f. , hauschildt p.h . , 1998 , a&a 337 , 403 ( bcah98 ) . barnes t.g . , evans d.s . , 1976 , mnras 174 , 489 . bailey j. , 1981 , mnras 197 , 31 . berriman , g. , kenyon , s. , boyle , c. , 1987 , aj 94 , 1291 beuermann , k. , weichhold , m. , 1999 , a&a , in preparation beuermann k. , baraffe i. , hauschildt , p. , 1999 , a&a , in press ( bbh99 ) . beuermann k. , baraffe i. , kolb u. , weichhold m. , 1998 , a&a 339 , 518 ( bbkw98 ) . clemens j.c . , , gizis j.e . , obrien m.s . , 1998 , apj 496 , 392 . echeverra j. , 1983 , rev.mex.astron.astrof . harrison t.e . , mcnamara b.j . , szkody p. , mcarthur , b.e . , benedict , g.f . , klemola , a.r . , gilliland , r.l . , 1999 , apj 515 , l93 . harrop - allin , m.k . , cropper , m. , potter , s.b . , dhillon , v.s . , howell , s.b . , 1997 , mnras 288 , 1033 . hauschildt p.h . , allard f. , baron e. , 1999 , apj 512 , 377 . ( _ nextgen _ models ) . henry t.j . , mccarthy jr . w. , 1993 , aj 106 , 773 . kolb , u. , baraffe i. , 1999 , submitted to mnras . leggett s.k . , allard f. , berriman g. , dahn c.c . , hauschildt p.h . , 1996 , apjs 104 , 117 ( l96 ) . marsh t. , 1990 , apj 357 , 621 . patterson j. , 1984 , apjs 54 , 443 . ramseyer t. , 1994 , apj 425 , 243 . ritter h. , kolb u. , 1998 , a&a 129 , 83 . schwarz r. , schwope a.d . , beuermann k. et al . , 1998 , a&a 338 , 465 ( polar rxj0203 + 29 ) . schwope a.d . , beuermann k. , cataln , m.s . , metzner , a. , smith , r.c . , steeghs , d. , 1999 , a&a in press ( qq vul ) . smith d.a . , dhillon v.s . , 1998 , mnras 301 , 767 . smith r.c . , hawkins , n.a . , mehes , o. , 1999 , this conference . southwell k.a . , still m.d . , smith r.c . , martin j.s . , 1995 , a&a 302 , 90 . stehle r. , kolb u. , ritter h. , 1999 , mnras in press . thomas , h .- c . , beuermann , k. , burwitz , v. , reinsch , k. , schwope , a.d . , 1999 , a&a , submitted ( polar rxj1313 - 32 ) . thorstensen , j.r . , taylor , c.j . , kemp , j. , 1999 , this conference . wade r.a . , horne k. , 1988 , apj 324 , 411 . warner , b. , 1987 , mnras 227 , 23 . warner b. , 1995 , cataclysmic variable stars , cambridge astrophysics series 28 , ( cambridge : cup ) . whyte c. , eggleton p. , 1980 , mnras 190 , 801 .
Abstract:

Answer: recent theoretical and observational progress has substantially improved the definition of the lower main sequence and established a new basis for a comparison of main sequence stars and the secondaries in cvs . the evolutionary sequences of kolb & baraffe [ 1999 ] imply that the secondaries in many cvs are expanded compared with main sequence stars of the same mass as a consequence of unusually high mass transfer rates and/or pre - cv nuclear evolution . we show that the location of the secondaries of all well - studied cvs in the spectral type period diagram implies that they are consistent with having near - solar metallicities . we show , furthermore , that the surface brightness of k / m stars depends on gravity and metallicity and present new barnes - evans relations valid for dwarfs of near - solar metallicity and the secondaries in cvs of the galactic disk population . distances derived by the surface brightness method agree with recent measurements of the trigonometric parallaxes of a few selected systems .