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 physics of normal - metal  superconducting ( ns ) hybrid structures has been of interest for several decades . in the sixties , the scattering of a quasi - particle at a ns interface was first described by andreev . @xcite thus , it is now widely known that an electron approaching to a superconductor from the normal - metal side , with an energy above the fermi level , can be reflected at the interface as a hole with an energy below the fermi level and same momentum , but in opposite direction . a cooper pair is thus formed in the superconductor . this scattering process is known as andreev reflection . more recently , the interest in this subject was renewed in the field of mesoscopic physics , where the electron wave function remains coherent and therefore quantum interference effects play an important role . several effects of phase - coherence andreev scattering are under current investigation . in addition , the effects of disorder and ballistic - chaotic scattering in hybrid structures have been widely studied . for example quantities such as the density of states and the ehrenfest time as well as phenomena such as coherent backscattering , excitation and gap fluctuations , superconducting proximity effect , and resonant tunneling have been investigated in mesoscopic hybrid structures.@xcite also , the effects of andreev scattering on the conductance through quantum dots whose classical dynamics is integrable or mixed ( regular and chaotic ) have been investigated.@xcite    on the other hand , random matrix theory ( rmt ) has been useful to study the statistical properties of many different phenomena in normal metals @xcite as well as in superconductors . in particular , the effects of ns interfaces on the conductance of disordered and chaotic billiards  quantum dots  have received considerable attention.@xcite recently , some statistical properties of the conductance fluctuations in mesoscopic ns heterostructures have been experimentally studied in ref .  [ ] where the explanation of the experimental results was based on rmt of quantum transport through chaotic quantum dots ; it was also pointed out the importance of considering the potential barriers formed at the ns interface .    using a scattering approach to the problem of quantum transport and rmt , here we study the statistical properties of the conductance through a quantum - chaotic dot with a superconducting lead attached to it . in particular , we are interested in the distribution of the conductance . the configurations of one and two normal leads additionally attached to the dot are considered . we take into account the presence of a schottky barrier  which is formed in real ns junctions of transparency @xmath0 in the superconducting lead . thus , for these two setups we calculate the complete distribution of conductances for different values of @xmath0 . as we will see , the conductance distribution is affected by the presence of the superconducting region as well as the scattering at the schottky barrier . this paper is organized as follows . in the next section we briefly review general ideas about electron - hole scattering at ns junctions . we also introduce the scattering matrices associated to the different elements of our experimental setup . in sec . iii , we study the two - terminal configuration . firstly , we calculate the total scattering matrix and the conductance . secondly , we present our results for the distribution of the conductance for the cases of time - reversal symmetry ( @xmath1 ) and broken time - reversal symmetry ( @xmath2 ) . simple analytical expressions are provided for the conductance distribution . iv is devoted to the three - terminal configuration . we calculate the necessary scattering - matrix elements in order to calculate the conductance and present the results for its distribution for both symmetries ( @xmath1 and @xmath2 ) , at different values of the transparency of the schottky barrier . a summary and conclusions of our work are presented in sec . v. we start our study with some general ideas on the scattering mechanism of electrons in a ns junction , as illustrated in fig . [ fig_1 ] . as we mentioned above , in a ns junction an incident electron from the normal side can be reflected by the superconductor as an electron or as a hole . the electron and hole wave functions @xmath3 and @xmath4 , respectively , satisfy the bogoliubov - de gennes equation@xcite given by @xmath5 where @xmath6 is the single electron hamiltonian and @xmath7 is the potential which couples the electron and hole wave functions . we assume that @xmath6 is independent of spin - orbit interactions . on the other hand , far from the ns boundary the superconducting gap @xmath7 goes to zero in the normal - metal region , whereas in the superconducting region the pair potential has an amplitude @xmath8 and phase @xmath9 . near the ns interface , however , @xmath7 depends on the details of the ns junction ; assuming that the width of the ns interface is small compared to the superconducting coherence length and neglecting any pairing interaction in the normal region , the pair potential is usually modeled as a step potential : @xmath10 , @xmath11 being the step function . now we consider that the energy of an incident electron is in the energy gap ( @xmath12 , @xmath13 being the fermi energy ) where there are no propagating wave solutions in the superconducting region , i.e. , the electron is reflected as a hole or electron . let us denote by @xmath14 the reflection matrix describing the andreev reflection . within the andreev approximation , where only andreev scattering processes are taking into account and with the above assumptions for @xmath7 , a solution of eq . ( [ bdeg ] ) in the normal region can be written as@xcite @xmath15\phi_{n}(y , z , e_f+\epsilon ) \nonumber           \\ & + & \left(\begin{array}{c }             0 \\ r^a_{he }           \end{array } \right)\exp[{ik_n^hx}]\phi_{n}(y , z , e_f-\epsilon ) \ .\end{aligned}\ ] ]    in this expansion , @xmath16 is the transverse wave function that carries a current in the @xmath17 direction and @xmath18 labels the transverse mode  open channel , while the wave number @xmath19 is given by @xmath20 and @xmath21 .    in the linear - response limit , i.e. , zero bias voltage ( @xmath22 ) , we have @xmath23 . in this limit , using the continuity of the wave function and its derivative at the interface , the reflection matrix @xmath24 is given by @xmath25 and @xmath26 , where @xmath27 is the @xmath28 unit matrix .     and @xmath29 are related by the scattering matrix @xmath30 , see eq . ( [ matrixsa ] ) . ]      as we have mentioned , we consider attached leads supporting one open channel , therefore from now on we specialize to this case . thus , the @xmath31 scattering matrix for the andreev reflection which relates the coefficients @xmath32 and @xmath29 , see fig . [ fig_1 ] , is given by @xmath33 with @xmath34 where @xmath35 and @xmath36 . we also consider the potential barrier created at the ns interface due to the mismatch of the conduction bands of the two different regions . we model this schottky barrier as a planar barrier characterized by its transparency @xmath0 . it is of relevance to consider the effect of this potential on the transport since , as we will see later , the conductance and therefore its statistical properties are affected depending on the strength of @xmath0 .    on the other hand , since in the normal region electrons and holes are not coupled , the scattering matrix of the schottky barrier @xmath37 can be written as @xmath38 with the scattering matrix for electrons @xmath39 given by@xcite @xmath40 while the scattering matrix for holes @xmath41 satisfies @xmath42 . @xcite      we now introduce the polar representation for the scattering matrix @xmath43 associated to the quantum dot . for electrons , the @xmath44 @xmath45-matrix is written as @xmath46 where @xmath47 and the phases @xmath48 , @xmath49 , @xmath50 , and @xmath51 are independent variables . in presence of time reversal symmetry , @xmath45 is a symmetric matrix and therefore @xmath52 and @xmath53 . as before , since in the normal region the electrons and holes are not coupled , the complete scattering matrix for the quantum dot can be written as @xmath54 with @xmath55 . first we consider a two - terminal setup . as shown in fig . [ fig_2 ] , a normal lead is attached to a normal quantum dot while a second lead with a superconducting region and a schottky barrier of transparency @xmath0 is also attached . both leads support one open channel . in a similar two - terminal setup and for a large number of channels in the leads , the average and variance of the conductance have been studied in ref .  [ ] . we are ready to calculate the full scattering matrix of the system , fig . [ fig_2 ] , since we have introduced already the scattering matrices associated to the components of our hybrid structure in the previous section , i.e. , we only need to combine properly the scattering matrices @xmath30 , @xmath37 , and @xmath43 from eqs . ( [ matrixsa ] ) , ( [ matrixss0 ] ) , and ( [ matrixsd ] ) , respectively . we combine first the scattering matrices associated to the andreev reflection @xmath30 and schottky barrier @xmath37 . the resulting scattering matrix @xmath56 is given by @xmath57 where we have defined the matrices @xmath58 from eq . ( [ matrixss ] ) , @xmath59 , @xmath60 , and @xmath61 , we obtain @xmath62 with @xmath63 . we now combine the above expression for @xmath56 with the scattering matrix of the quantum dot @xmath43 . the resulting scattering matrix @xmath64 is given by @xmath65 after some algebra we write the total matrix @xmath64 as @xmath66 where @xmath67 \ , \nonumber \\ s_{eh }   & = & -ice^{i\left(\phi+\psi-\psi'+\theta-\theta'\right ) } \gamma \tau \ , \nonumber \\ s_{he } & = & -ice^{-i\left(\phi+\psi-\psi'+\theta-\theta'\right ) } \gamma \tau \ , \nonumber \\ s_{hh } & = & 2ce^{-2i(\theta+\theta')}\big[(\gamma-2)\sqrt{1-\tau}+ \nonumber   \\ & & \sqrt{1-\gamma}\left((2-\tau ) \cos{(\psi+\psi')}-i\tau \sin{(\psi+\psi')}\right ) \big ] \ , \nonumber \\\end{aligned}\ ] ] with @xmath68^{-1 } \ .\ ] ]     is introduced in the superconducting lead . ]    finally , the dimensionless conductance @xmath69 due to the presence of the ns interface in our two - terminal setup is given , in terms of elements of the total scattering matrix , by@xcite @xmath70 where @xmath71 is the conductance quantum . we now calculate the distribution of the conductance by assuming that the statistics of the scattering matrix of the quantum dot @xmath43 follows one of the circular ensembles , namely , the circular orthogonal ensemble ( @xmath1 ) or the circular unitary ensemble ( @xmath2 ) .    from eq . ( [ conductancens ] ) , the distribution of the conductance is given by @xmath72 with invariant measure@xcite @xmath73    in the following subsections we provide analytical expressions for the distribution of the conductance for the two - terminal configuration by the evaluation of the integrals in eq . ( [ defpofg ] ) , for the cases @xmath1 and @xmath2 . we recall that for the orthogonal symmetry , @xmath74 and @xmath75 are symmetric scattering matrices and therefore @xmath52 and @xmath53 in the polar representation of eq . ( [ polar ] ) . thus , from eq . ( [ elementsofstotal ] ) we find @xmath76 ^ 2 } \ .\ ] ] before considering the general case of arbitrary @xmath0 let us consider the particular case @xmath77 , i.e. , maximum transparency of the schottky barrier . in this case eq . ( [ gbeta1 ] ) reduces to @xmath78 thus from eqs . ( [ conductancens ] ) and ( [ defpofg ] ) , we obtain the simple expression for the distribution @xmath79^{3/2 } } \ , \ \ \mathrm{for } \ \ \gamma=1,\ ] ] which we have already normalized . we now study the case of arbitrary @xmath0 . therefore we need to perform the integrals in eq . ( [ defpofg ] ) with @xmath80 given by eq . ( [ gbeta1 ] ) . the calculation is direct but somewhat lengthy . we finally get @xmath81^{1/2 } \times \nonumber   \\ & & \ \ \ \ \ \ \ \   \ \   \   \left[e(x , m)-f(x , m ) \right ] \big\vert_{\tau_-}^{\tau_+ } \ , \end{aligned}\ ] ] where @xmath82 @xmath83 and @xmath84 in eq . ( [ pofggeneralbeta1 ] ) are the elliptical integrals of first and second kind , respectively , which have to be evaluated at the limits @xmath85 and @xmath86 given by @xmath87^{-1 } \ .\ ] ] we have verified our predictions for @xmath88 , eqs . ( [ pofgamma1 ] ) and ( [ pofggeneralbeta1 ] ) , by comparing with numerical simulations of a quantum - chaotic dot ( in all our numerical simulations we have considered a bunimovich stadium geometry whose classical dynamics is chaotic ) . in fig . [ fig_3 ] , we show both numerical and analytical results for @xmath89 , 0.5 , and 1 . we observe good agreement between theory and the corresponding numerical simulations . for the two - terminal setup for different values of the transparency @xmath0 . good agreement is observed between theory [ eqs . ( [ pofgamma1 ] ) and ( [ pofggeneralbeta1 ] ) ] and the corresponding numerical simulations ( histograms ) . ]     for the two - terminal setup for different values of the transparency @xmath0 . good agreement is observed between theory [ eqs . ( [ pofggamma1beta2 ] ) and ( [ pofggeneralbeta2 ] ) ] and the corresponding numerical simulations ( histograms ) . ] we now consider the circular unitary ensemble . we recall that this symmetry class is relevant for situations where time - reversal symmetry is absent . from eq . ( [ elementsofstotal ] ) we have @xmath90 ^ 2 } \ . \nonumber \\\end{aligned}\ ] ] as in the previous subsection we start with the special case @xmath77 . thus , from eqs . ( [ conductancens ] ) , ( [ defpofg ] ) , and ( [ gbeta2 ] ) the distribution of conductance is given by @xmath91 the general case of arbitrary @xmath0 involves a direct but lengthy calculation , as in the @xmath1 case . for @xmath2 , however , the distribution gets a simpler expression . after some algebra we find @xmath92\gamma^2}{\sqrt{g}\left [ \sqrt{2}\gamma+\sqrt{g}\left(2-\gamma \right)\right]^3 } \ .\ ] ] in fig . [ fig_4 ] we plot @xmath93 , as given by eqs . ( [ pofggamma1beta2 ] ) and ( [ pofggeneralbeta2 ] ) , and compare with numerical simulations of a quantum - chaotic dot . we choose @xmath89 , 0.5 , and 1 . as for the case @xmath1 , for broken time - reversal symmetry a good agreement between theory and numerics is also observed .    to conclude this section we notice that , for the two - terminal configuration studied here and for both symmetries ( @xmath1 and @xmath2 ) , the distribution of the conductance becomes broader as the value of the transparency of the schottky barrier is increased , see figs . [ fig_3 ] and [ fig_4 ] . this broadening of the distribution implies that large values of @xmath69 ( values of @xmath69 at the tail of the distribution ) become statistically favored as the schottky barrier becomes transparent . we anticipate that the effect of the schottky barrier on the conductance is somewhat different in the three - terminal configuration that we study in the following section . we now consider a more complex situation . as shown in fig . [ fig_5 ] , two normal leads are attached to a normal quantum dot while a third lead with a superconducting material and a schottky barrier is also attached . all the leads support one open channel . as in the previous section , we calculate the total scattering matrix of the hybrid structure . the present case is , however , algebraically more involved than the two - terminal configuration . is introduced in the superconducting lead . ] we define a scattering matrix @xmath64 that relates the incoming and outgoing wave amplitudes in the normal leads ( leads 1 and 2 in fig .  [ fig_5 ] ) : @xmath94 where @xmath64 has the following structure @xmath95 here , the subscripts denote the leads . see fig . [ fig_5 ] . we recall that we are working at energies where there are no propagating waves in the superconducting region . in order to obtain the explicit expressions for the elements of the scattering matrix @xmath64 , we need to combine the scattering matrix of the superconducting region [ eq . ( [ matrixsa ] ) ] , the schottky barrier [ eqs . ( [ matrixss0 ] ) and ( [ matrixss ] ) ] , and the quantum dot . the scattering matrix associated to the quantum dot @xmath43 can be written in two blocks , as in section [ introductionc ] , one for the electrons and the other for the holes . thus , for the three - terminal setup , we write for electrons @xmath96 with a similar expression for the holes . once we have introduced the scattering matrices associated to the different parts of our setup , we can follow the same steps of the previous section to obtain the total scattering matrix . on the other hand , however , the dimensionless conductance @xmath69 for our three - terminal setup is given by@xcite @xmath97 thus , we need only the elements @xmath98 and @xmath99 of the total scattering matrix @xmath64 , eq . ( [ matrixst2leads ] ) , in order to calculate the conductance . after some algebra we find for these two scattering - matrix elements @xmath100 where @xmath101 \ , \\ a^e&= & t'^s_e t_{31}^e\frac{1}{1-r_{33}^h r'^s_h } \ , \\ \label{definitions1 } c^e&= & r^s_e+\frac{t'^s_e r_{33}^e t^s_e}{1-r_{33}^e r'^s_e } \ .\end{aligned}\ ] ] @xmath102 and @xmath103 in eqs . ( [ s21ee ] ) and ( [ s21eh ] ) are defined as @xmath104 and @xmath105 , respectively , by replacing @xmath106 . we can observe that if the third lead were absent , the expression for conductance , eq . ( [ gtwoterminal ] ) , is reduced to the usual landauer - bttiker formula @xmath107 .    for the particular values of @xmath108 and 1 , the expressions for the elements @xmath98 and @xmath99 are simplified . therefore for these two special values of @xmath0 , it is useful to show the expressions of @xmath98 and @xmath99 since they exhibit explicitly the different scattering processes among the three leads attached to the quantum dot that contribute to the conductance .      for the case of zero reflection at the schottky barrier the expressions for @xmath98 and @xmath99 in eqs . ( [ s21ee ] ) and ( [ s21eh ] ) , respectively , are simplified to @xmath109 therefore the conductance is given by @xmath110      for the case of total reflection at the schottky barrier the conductance is given by @xmath112 clearly , there is no electron - hole scattering processes in this case .      in order to calculate the distribution of conductances we need to average over the orthogonal or unitary ensembles ( @xmath1 and @xmath2 ) , as in the previous section . ( [ defpofg ] ) . in the three - terminal configuration , however , the expression for the conductance [ eqs . ( [ gtwoterminal])-([s21eh ] ) ] is more complicated , even at the values of the transparency @xmath113 and 1 , where the expressions for the scattering - matrix elements @xmath98 and @xmath99 are simpler . therefore we are not able to provide analytical expressions for @xmath114 , as for the two - terminal setup . thus we proceed with our analysis by generating numerically an ensemble of conductances from eqs . ( [ gtwoterminal])-([definitions1 ] ) , where the elements of the scattering matrix associated to the quantum dot are generated according to the orthogonal and unitary ensembles .    in presence of time - reversal symmetry ( @xmath1 ) , the distribution of the conductance at different values of the transparency is shown in fig . [ fig_6 ] . at zero transparency , the distribution is described by@xcite @xmath115 ( solid line in fig . [ fig_6 ] ) , which is the known result for a quantum - chaotic dot with two normal single - channel leads , i.e. , there is no effect of the superconducting region on the conductance , as one might expect for a nontransparent schottky barrier . for @xmath116 , however , @xmath88 is affected . from fig . [ fig_6 ] we can see that for intermediate values of the transparency ( @xmath89 and @xmath117 ) the distribution probability is finite for @xmath118 , whereas for null and complete transparency ( @xmath108 and @xmath77 , respectively ) the probability is confined in the interval @xmath119 . therefore multiple scattering in the third lead , where the schottky barrier is present , produces larger values of conductance . when time reversal symmetry is broken ( @xmath2 ) the distribution of the conductance is quite different from the @xmath1 case . in fig . [ fig_7 ] we show @xmath93 at different values of the transparency . first we notice that for @xmath111 , @xmath93 reduces to the known result@xcite @xmath120 ( shown as the solid line in fig .  [ fig_7 ] ) . at intermediate values of the transparency of the schottky barrier ( @xmath89 and @xmath117 ) , we observe again that multiple scattering in the superconducting lead can produce larger values of conductance than for the cases @xmath111 and @xmath77 , i.e. , the distribution @xmath93 develops a tail for @xmath121 at intermediate values of @xmath0 .     for the three - terminal setup for different values of the transparency @xmath0 . notice that @xmath122 for @xmath123 , see right panels . the solid curve for @xmath111 is given by @xmath124 . ]     for the three - terminal setup for different values of the transparency @xmath0 . as for @xmath1 , here @xmath125 for @xmath123 , see right panels . the horizontal line for @xmath111 corresponds to @xmath120 . ] we have studied the effects of a superconducting lead attached to a normal quantum - chaotic dot on the statistical properties of the conductance using random matrix theory and a scattering approach to the transport in normal - metal  superconducting ( ns ) hybrid structures . we have calculated the complete distribution of the conductance in the presence and absence of time - reversal symmetry ( @xmath1 and @xmath2 , respectively ) and for two different configurations : two and three - terminal setups . the conductance distribution has been analyzed as a function of the transparency @xmath0 of the schottky barrier which is modeled as a planar barrier formed at the ns interface . the importance of considering the effect of a potential barrier formed at ns junctions has been already recognized in transport experiments,@xcite however , opposite to the one open channel case that we study here , in these experiments a large number of channels are open . for the two - terminal configuration we provided analytical expressions for the distribution of the conductance . we found that the distribution becomes more and more concentrated at @xmath126 as the value of the transparency is reduced , for both symmetries @xmath1 and @xmath2 . this behavior might be expected considering that the conductance is due to the scattering at the ns interface ( andreev reflection ) , i.e. , at zero transparency of the barrier the situation is reduced to a normal quantum dot with a normal lead attached where only reflection processes at the lead are present . for the three - terminal configuration , the analytical expressions for the total scattering matrix and the conductance are more complicated than for the two - terminal case . we thus obtain the distribution of the conductance from an ensemble of conductances generated numerically according to our analytical expressions for the elements of the total scattering matrix . we found that in presence and absence of time - reversal symmetry , conductance values larger than unity are possible at intermediate values of the transparency of the schottky barrier . as a consequence , the distribution of the conductance is broader than for the cases of zero transparency ( @xmath111 ) and an ideal ns interface ( @xmath77 ) , i.e. , multiple scattering in the superconducting region with the schottky barrier increases the possibility of larger values of the conductance .    to conclude , we have seen that the presence of a schottky barrier in the lead , which is formed in real experiments due to the mismatch of the conduction bands of the normal - metal and the superconductor , has a significant effect at the level of the distribution of the conductance . the progress in the fabrication of nanoscale heterostructures has made possible the experimental study of statistical properties of the conductance through such devices.@xcite these experiments have pointed out the relevance of non - perfect contacts on the study of conductance . we thus think that the effects on the statistics of the conductance through quantum dots due to the presence of andreev scattering and potential barriers at the leads , as those described here , are of experimental interest . this work was supported by the  ramn y cajal " program of the spanish ministry of education and science and fondo social europeo . j. a. m .- b . thanks the depto . de fsica terica for its hospitality , as well as the program ayudas para estancias de investigadores de excelencia of the university of zaragoza for financial support . he also acknowledges support from conacyt mexico ( grant cb-2006 - 01 - 60879 ) and viep - buap . a. h. a. acknowledges hospitality at depto . de fsica terica , universidad de zaragoza , and financial support from the gobierno de aragn ( programa de movilidad de investigadores m1034 - 2007 ) and bifi .
Abstract:

Answer: we investigate the conductance statistics of a quantum - chaotic dot  a normal - metal grain  with a superconducting lead attached to it . the cases of one and two normal leads additionally attached to the dot are studied . for these two configurations the complete distribution of the conductance is calculated , within the framework of random matrix theory , as a function of the transparency parameter of the schottky barrier formed at the interface of the normal - metal and superconducting regions . our predictions are verified by numerical simulations .


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: one can expect that the effective interaction between quasiparticles in strongly correlated electron systems to be very complex . the interaction will depend obviously on the charge , but also more generally on the spin and current carried by the quasiparticles . on the border of long - range magnetic order it is plausible that the dominant interaction channel is of magnetic origin and depends on the relative spin orientations of the interacting quasiparticles . it has been shown that this magnetic interaction treated at the mean field level can produce anomalous normal state properties and superconducting instabilities to anisotropic pairing states . it correctly predicted the symmetry of the cooper state in the copper oxide superconductors@xcite and is consistent with spin - triplet p - wave pairing in superfluid @xmath0 [ for a recent review see , e.g. , ref . one also gets the correct order of magnitude of the superconducting and superfluid transition temperature @xmath1 when the model parameters are inferred from experiments in the normal state of the above systems . there is growing evidence that the magnetic interaction model may be relevant to other materials on the border of magnetism .    thus far the magnetic interaction model has been explored in very simple cases . the most extensively investigated example is that of a nearly half - filled single - band in a square or cubic lattice . these studies have revealed a number of interesting features that are quite in contrast to those expected for conventional phonon mediated pairing . in the latter case , the interaction is local in space , but non - local in time , whereas on the border of magnetism , one expects the interaction to be strongly non - local in both space and time . for nearly antiferromagnetic metals the magnetic interaction is oscillatory in space and superconductivity depends on the ability of the electrons in a cooper pair state to sample mainly the attractive regions of these oscillations . because of the strong retardation in time , the relative wavefunction of the cooper pair must be constructed from bloch states with wavevectors close to the fermi surface . furthermore , the allowed symmetries of the cooper pair wavefunction are restricted by the crystal structure . the possibility of constructing a cooper pair state with maximum probability in the attractive regions of the magnetic interaction can be severely constrained by these requirements . therefore , one expects that the robustness of magnetic pairing to be very sensitive to details of the electronic and lattice structures .    on the border of ferromagnetism , one is not hampered by the oscillatory nature of the magnetic interaction which in the simplest model is attractive everywhere in space and time in the spin - triplet channel . at first sight this would seem to be the most favourable case for magnetically mediated superconductivity . however , the results of the numerical calculations presented in ref . @xcite indicate that the highest mean field @xmath1 for the cases considered is obtained for d - wave pairing in the nearly antiferromagnetic state in a quasi-2d tetragonal lattice . in this particular case , as explained in ref . @xcite , it turns out to be possible to ideally match the cooper pair state to the attractive regions of the magnetic interaction .    on the border of ferromagnetism , magnetic pairing in the spin - triplet state has the disadvantage that only the exchange of magnetic fluctuations polarised along the direction of the interacting spins , i.e. , longitudinal fluctuations , contribute to the quasiparticle interactions . for a spin - rotationally invariant system , both longitudinal and transverse fluctuations contribute to pairing only for a spin - singlet state . another disadvantage of being on the border of ferromagnetism is that for otherwise similar conditions the suppression of @xmath1 due to the self interaction arising from the exchange of magnetic fluctuations is stronger than in the corresponding case on the border of antiferromagnetism . this disadvantage can be mitigated in systems with strong magnetic anisotropy in that the effect of the transverse magnetic fluctuations on the self interaction would be suppressed while the strength of the pairing interaction arising from the longitudinal magnetic fluctuations need not be reduced . this may apply in systems with strong spin - orbit interactions or in the spin - polarised state close to the border of ferromagnetism . these arguments  @xcite have stimulated a new search for evidence of superconductivity on the border of itinerant electron ferromagnetism in cases where spin anisotropy is expected to be pronounced , such as @xmath2 . this search has proved fruitful because it led to the first observation of the coexistence of superconductivity and itinerant electron ferromagnetism in @xmath2@xcite and shortly thereafter in @xmath3@xcite and @xmath4@xcite . the prediction of the simple model presented in ref .  @xcite that magnetic pairing is more robust in the quasi-2d square lattice than in the cubic structure seems to have been borne out by recent experiments . namely , one finds an order of magnitude increase in the maximum @xmath1 and in the range in pressure where superconductivity is observed on the border of metallic antiferromagnetism when the simple cubic lattice of @xmath5@xcite is stretched along one principal axis by the insertion of non - magnetic layers to form the tetragonal compounds @xmath6@xcite ( @xmath7 is @xmath8 , @xmath9 or @xmath10 ) . these systems , albeit quite anisotropic , would not normally be considered to be quasi - two dimensional and it is not clear at first sight that the model calculations carried out in ref . @xcite are directly relevant . the purpose of this paper is to show that for comparable model parameters , the robustness of magnetic pairing increases gradually as one goes from a cubic structure to a more and more anisotropic structure on the border of metallic antiferromagnetism and ferromagnetism . this behaviour of the mean field transition temperature is in stark contrast to that of the `` one - loop '' fluctuations corrections to @xmath1 . the latter corrections typically depend logarithmically on the degree of anisotropy and would be expected to be negligible for materials such as @xmath6 . we do not expect some of the results of ref . @xcite to be generic properties of the magnetic interaction model . we have already stressed that even in simple cases , the robustness of magnetic pairing can be very sensitive to certain details of the lattice and electronic structure . even in the single band problem , many such structures have not yet been extensively studied theoretically . furthermore , we expect that the range of possibilities to be greatly expanded in the presence of more than one partially filled electronic band . most known materials on the border of magnetism crystallise in other than simple cubic or tetragonal structure and have more than one band crossing the fermi level . for these more complex systems , one would not expect the model of ref . @xcite to be directly relevant . for example , the observation of spin - triplet rather than spin - singlet d - wave pairing in some multi - band materials with strongly enhanced antiferromagnetic spin fluctuations , such as @xmath11 and @xmath12 , may not be inconsistent with the idea of magnetic pairing . a detailed study of magnetic pairing in multi - band systems for a range of crystal structures would shed light on the possible forms of superconductivity and the conditions most favorable for their observation . the simple model calculations suggest that anisotropic forms of superconductivity should be a generic property of systems on the border of metallic magnetism . it may seem surprising therefore that there are still so few observations of this phenomenon . in many cases , the multiplicity of bands and , for example , magnetic fluctuations in the non - bipartite lattice may weaken magnetic pairing to such an extent that quenched disorder may completely suppress superconductivity . an illustration of this point is the dramatic collapse of the spin - triplet superconducting transition temperature in @xmath12 in the presence of @xmath13 impurity concentrations as low as 0.1% .    at first sight , the magnetic interaction model is mathematically analogous to the conventional electron - phonon problem with the generalised magnetic susceptibility playing the role of the phonon propagator . one would therefore expect that a simple analytic expression similar to that proposed by mcmillan could be used to represent approximately the @xmath1 calculated numerically via the eliashberg equations . our attempts in this direction have not , however , proved successful@xcite . a recent study suggests that there may be a fundamental reason for inapplicability of the mcmillan - style expression for @xmath1@xcite . on the border of long - range magnetic order , the incoherent part of the electron green function , which is ignored in the simplest treatment , plays a major role in the formation of the cooper pair condensate . the traditional picture in which superconductivity arose from pairing of well defined ( weakly damped ) quasiparticles appears inadequate on the border of metallic magnetism even in the mean field eliashberg treatment . we note that in our model the coupling of the quasiparticles to the magnetic fluctuations is a phenomenological constant to be inferred from normal state properties that formally includes that part of the vertex correction which is local in space and time . calculations have shown that the neglect of vertex corrections that are non - local in space and time is justified at least in some cases of physical interest@xcite . when the magnetic correlation length becomes sufficiently large , however , these neglected non - local vertex corrections ( including superconducting phase fluctuations ) may become important . their effect on @xmath1 and the normal state properties are as yet incompletely understood . we consider quasiparticles in a simple tetragonal lattice described by a dispersion relation    @xmath14    with hopping matrix elements @xmath15 and @xmath16 . @xmath17 represents the electronic structure anisotropy along the z direction . @xmath18 corresponds to the quasi-2d limit while @xmath19 corresponds to the 3d cubic lattice . for simplicity , we measure all lengths in units of the respective lattice spacing . in order to reduce the number of independent parameters , we take @xmath20 and a band filling factor @xmath21 as in our earlier work .    the effective interaction between quasiparticles is assumed to be isotropic in spin space and is defined in terms of the coupling constant @xmath22 and the generalised magnetic susceptibility which is assumed to have a simple analytical form consistent with the symmetry of the lattice . @xmath23    where @xmath24 and @xmath25 are the correlation wavevectors or inverse correlation lengths in units of the lattice spacing in the basal plane , with and without strong magnetic correlations , respectively . let    @xmath26    where @xmath27 parameterises the magnetic anisotropy . @xmath28 corresponds to quasi-2d magnetic correlations and @xmath29 corresponds to 3d magnetic correlations .    in the case of a nearly ferromagnetic metal the parameters @xmath30 and @xmath31 in eq . ( [ chiml ] ) are defined as    @xmath32    where @xmath33 is a characteristic spin fluctuation temperature . note that our definition of @xmath33 may differ from the characteristic spin fluctuation temperature scales used by other authors .    in the case of a nearly antiferromagnetic metal , the parameters @xmath30 and @xmath31 in eq . ( [ chiml ] ) are defined as    @xmath34    as in our previous work@xcite , the band structure and generalized magnetic susceptibility are modeled independently . this choice may be inconsistent when all of the contributions to @xmath35 come from the chosen band . however , it allows us , in principle , to deal with the case where there are other important contributions to the generalized magnetic susceptibility . it has been argued that the latter case is of relevance to the ruthenates@xcite , and most likely the heavy - fermion systems .    a complete description of the model , the eliashberg equations for the superconducting transition temperature and their method of solution can be found in the appendix . we note that the model is fully defined by the phenomenological parameters describing the electronic structure @xmath36 , the generalised magnetic susceptibility @xmath35 and the interaction vertex @xmath22 . in principle , these parameters can be estimated from experimental studies of the normal state . in particular , the resistivity can be used to estimate the dimensionless coupling parameter @xmath37 the value of which is between 10 and 20 for the simplest rpa model for the magnetic interaction . the dimensionless parameters at our disposal are @xmath37 , @xmath38 , @xmath25 and @xmath24 . for comparison with the results of our earlier work@xcite , we take @xmath39 and @xmath40 . in 2d , this @xmath33 corresponds to about @xmath41 for a bandwidth of 1 ev while our choice of of @xmath42 is a representative value . the results of our numerical calculations of the mean field critical temperature @xmath1 as a function of the electronic and magnetic anisotropy parameters @xmath17 and @xmath27 , respectively , are shown in figs 1 and 2 for representative values of the parameters @xmath43 and @xmath37 . figures 1a - c illustrate the results for a nearly anti - ferromagnetic metal and figures 2a - c for a nearly ferromagnetic metal . note that our previous calculations correspond to the quasi-2d case @xmath44 and to 3d case @xmath45 . a glance at figure 1 reveals a clear pattern in the variation of @xmath1 with both anisotropy parameters @xmath17 and @xmath27 . we notice that @xmath1 increases gradually and monotonically as the system becomes more and more anisotropic in either the electronic structure or in the magnetic interaction . in going from 3d to quasi-2d , @xmath46 is found to increase by up to an order of magnitude for otherwise fixed parameters of the model . the increase becomes least pronounced as for small @xmath43 and large @xmath37 . the behavior in the nearly ferromagnetic case , figure 2 , though broadly similar to that of the nearly antiferromagnetic metal , shows some interesting differences . in some cases , the minimum @xmath1 occurs for 3d electronic structure , but quasi-2d magnetic interaction . also , in all cases considered the maximum @xmath1 is obtained for a quasi-2d electronic structure and strongly anisotropic , but not 2d magnetic interactions .      in order to make a comparison with the corresponding electron - phonon problem it is instructive to define a mass renormalization parameter @xmath47 and interaction parameter @xmath48 . we define    @xmath49    where    @xmath50    and    @xmath51    for p - wave spin triplet pairing @xmath52 while    @xmath53    in the case of d - wave spin - singlet pairing @xmath54 . the fermi surface averages are given by    @xmath55    in practice , we compute the fermi surface average with a discrete set of momenta on a cubic or tretragonal lattice and we replace the delta function by a finite temperature expression    @xmath56    where @xmath57 is the fermi function . note that @xmath58 as @xmath59 . we have used @xmath60 and @xmath61 in all of our calculations . the finite temperature effectively means that van hove singularities will be smeared out . note that the fermi surface average that appears in @xmath47 , eq . ( [ lambda1 ] ) plays a role similar to that of @xmath62 in the case of phonon mediated superconductivity . from the definitions of the parameters @xmath63 eqs . ( [ lambda1 ] ) ,  ( [ lambda2 ] ) and our model for @xmath35 eq . ( [ chiml ] ) , we see that @xmath63 are directly proportional to the dimensionless factor @xmath64 . thus we will consider the quantities    @xmath65    which are functions only of @xmath66 , @xmath67 and @xmath43 .    in figures 3 and 4 we show @xmath68 , @xmath69 and the ratio @xmath70 for a representative value of @xmath43 in the case of a nearly antiferromagnetic metal and nearly ferromagnetic metal , respectively . the trends in both cases are the same . @xmath68 and @xmath69 are seen to increase gradually and monotonically in going from 3d to quasi-2d . however , @xmath69 grows faster than @xmath68 so the ratio @xmath70 also increases in going from 3d to quasi-2d . this qualitative trend in the ratio is consistent with the behavior of @xmath1 obtained from the numerical solution of the eliashberg equations . in the ferromagnetic case , however , it fails to reproduce the fact that the minimum @xmath1 is not necessary for a fully 3d system and that the maximum @xmath1 is obtained for strongly anisotropic yet not quasi-2d systems . the results of the calculations for both the nearly ferromagnetic and nearly antiferromagnetic metals show that the robustness of magnetic pairing increases gradually as one goes from a cubic to a more and more anisotropic structure with parameters other than @xmath27 and @xmath17 left unchanged . these results are consistent with our previous findings  @xcite and with the calculations for @xmath71 and @xmath72 presented in ref . @xcite . in an earlier study , nakamura et al . @xcite found that @xmath1 could increase by up to a factor of three in going from 3d to 2d for their choice of model parameters . the effect of anisotropy on @xmath1 for nearly ferromagnetic and nearly antiferromagnetic metals is qualitatively similar . this phenomenon arises from the increase with growing anisotropy of the density of states of both the quasiparticles and of the magnetic fluctuations that mediate the quasiparticle interaction . this effect could be further enhanced in the case of a nearly antiferromagnetic metal by the change in the pattern of the oscillations of the magnetic interaction . it can be seen from figure 5 that the strength of the interaction in the repulsive sites outside of the nodal plane of the @xmath73 state gets reduced while crucially the attraction in the basal plane gets enhanced as one goes from the cubic to a more and more anisotropic tetragonal lattice . this enhancement is the consequence of the increase of the phase space of soft magnetic fluctuations as one goes from a cubic to a quasi - two dimensional structure . since our model potential varies smoothly with the tetragonal distortion , parameterised by @xmath27 in figures 1 and 2 , it is clear that these effects occur gradually with increasing separation between the basal planes . the calculations assume that the maximum magnetic response for a nearly antiferromagnetic metal occurs at the commensurate wavevector defined by @xmath74 and @xmath75 , where @xmath76 and @xmath77 are the lattice constants in the basal plane and along the tetragonal axis respectively , reintroduced here for clarity . the oscillations in the magnetic interaction potential along the tetragonal axis obviously depend on the value of @xmath78 . however , the enhancement of the attraction in the basal plane and the reduction of the interaction elsewhere as one goes from a cubic to a more and more anisotropic lattice do not depend on the particular value of @xmath78 . therefore , we expect the qualitative conclusions of this paper to be independent of @xmath78 . the robustness of the pairing is further enhanced by the gradual change in the electronic band from a 3d to a quasi-2d form ( see eq . ( [ eps ] ) ) . the reduced hopping along the distortion axis , parameterised by @xmath17 in figures 1 and 2 , implies a reduced electronic bandwidth and hence increased density of electronic states . our calculations show that this too leads to a gradual increase in @xmath1 with increasing distortion of the lattice .    in a nearly ferromagnetic metal , one again benefits from the reduction of the electronic band width and the increase of the interaction in the basal plane as one goes from a cubic to a tetragonal lattice ( see figure 6 ) . however , the suppression of the interaction between the basal planes has a less dramatic effect on the border of ferromagnetism than antiferromagnetism because in the latter case one suppresses key repulsive regions of the interaction ( figure 5 ) . these simple arguments explain how the pairing effects of the interaction are strengthened by a tetragonal distortion in our model . however , the same effects also contribute to an enhanced self - interaction which acts to suppress @xmath1 . the relative importance of the pair forming and pair breaking effects of the magnetic interaction can not be inferred by the above physical picture alone . the numerical calculations show that for most cases considered here the pair forming effects dominate . the balance is particularly delicate on the border of ferromagnetism where the suppression of @xmath1 brought about by the self - interaction is pronounced . a physical interpretation of this suppression of @xmath1 is given in ref . the same interpretation may explain , for example , why the maximum of @xmath46 in the nearly ferromagnetic case is for a strongly anisotropic yet not quasi-2d pairing potential ( figure 2 ) . a most striking manifestation of the interplay between the pair - forming and pair - breaking tendency of the magnetic interaction is the breakdown of the mcmillan - style expression for @xmath1 in terms of the parameters @xmath48 and @xmath47 ( see eqs . ( [ lambda1],[lambda2 ] ) . this was noted in ref . @xcite and has been interpreted in ref . @xcite in terms of the important role played by the incoherent part of the green function which is ignored in the simplest treatments , but is included in the present and earlier work@xcite where the full momentum and frequency dependence of the self energy is taken into account . the calculations show that the lattice anisotropy may increase the robustness of magnetic pairing in the mean - field approximation . superconducting phase fluctuations which are not included in this approximation may be expected to suppress @xmath1 in the 2d limit . therefore , in practice , one would think that the most favorable case for magnetic pairing is that of strong but not extreme anisotropy . as noted in the introduction and in the previous two sections , the robustness of magnetic pairing can be very sensitive to certain details of the magnetic interaction and electronic structure . therefore , one should exercise caution in making quantitative comparisons between the results of our calculations and experiment . for instance , one would expect all of the parameters of the model ( not solely @xmath27 and @xmath17 ) to change simultaneously with increasing lattice anisotropy . the changes brought about in going from a cubic to a tetragonal lattice may even be much more complex than considered here . in particular , the number of partially filled bands may itself change . as also mentioned in the introduction , this could have in some cases even more dramatic consequences on superconductivity than the effects taken into account in our simple one - band model . the theoretical framework developed for systems on the border of magnetism can be translated to describe systems on the border of other types of instabilities , such as charge density wave or ferroelectric instabilities . the above given phase space argument to explain the increased robustness of magnetic pairing with increasing lattice anisotropy should carry over in part to these other pairing mechanisms , at least at the one - loop mean - field level ( see , e.g. , ref . @xcite ) . while some understanding of the properties of the magnetic interaction model has been gained over the last few years ( e.g. , the conditions for robust pairing of electrons ) , there are many cases where the predictions of the model have not been worked out . of particular importance is the role of the multiplicity of partially filled bands which may be expected to be the key to understanding exotic superconductivity observed in nearly magnetic materials such as @xmath11 and @xmath12 . we would like to thank a.v . chubukov , p. coleman , s.r . julian , p.b . littlewood , a.j . millis , a.p . mackenzie , d. pines , d.j . scalapino and m. sigrist for discussions on this and related topics . we acknowledge the support of the epsrc , the newton trust and the royal society . we consider quasiparticles on a cubic or tetragonal lattice . we assume that the dominant scattering mechanism is of magnetic origin and postulate the following low - energy effective action for the quasiparticles :          where @xmath83 denotes the three pauli matrices . the quasiparticle dispersion relation @xmath36 is defined in eq . ( [ eps ] ) , @xmath84 denotes the chemical potential , @xmath85 the inverse temperature , @xmath86 the coupling constant and @xmath87 and @xmath88 are grassmann variables . in the following we shall measure temperatures , frequencies and energies in the same units . to get @xmath89 to decay as @xmath92 as @xmath93 , as it should , we introduce a cutoff @xmath94 and take @xmath95 for @xmath96 . a natural choice for the cutoff is @xmath97 . we have checked that our results for the critical temperature are not sensitive to the particular choice of @xmath94 used .              where @xmath101 is the quasiparticle self - energy , @xmath102 the one - particle green s function and @xmath103 the anomalous self - energy . the chemical potential is adjusted to give an electron density of @xmath21 , and @xmath80 is the total number of allowed wavevectors in the brillouin zone . in eq . ( [ gap ] ) , the prefactor @xmath104 is for triplet pairing while the prefactor @xmath105 is appropriate for singlet pairing . only the longitudinal spin - fluctuation mode contributes to the pairing amplitude in the triplet channel . both transverse and longitudinal spin - fluctuation modes contribute to the pairing amplitude in the singlet channel . all three modes contribute to the quasiparticle self - energy . the momentum convolutions in eqs . ( [ sigma],[gap ] ) are carried out with a fast fourier transform algorithm on a @xmath106 lattice . the frequency sums in both the self - energy and linearized gap equations are treated with the renormalization group technique of pao and bickers@xcite . we have kept between 8 and 16 matsubara frequencies at each stage of the renormalization procedure , starting with an initial temperature @xmath107 and cutoff @xmath108.the renormalization group acceleration technique restricts one to a discrete set of temperatures @xmath109 . the critical temperature at which @xmath110 in eq . ( [ gap ] ) is determined by linear interpolation . p. monthoux , phys . b * 55 * , 15261 ( 1997 ) . chubukov , p. monthoux and d.r . morr , phys . b * 56 * , 7789 ( 1997 ) . p. monthoux , phil . mag . b * 79 * , 15 ( 1999 ) . s. nakamura , j.phys . . jpn . * 69 * , 178 ( 2000 ) . n.e . bickers and s.r . white , phys . rev . b * 43 * , 8044 ( 1991 ) .
Abstract:

Answer: we compare predictions of the mean - field theory of superconductivity for nearly - antiferromagnetic and nearly - ferromagnetic metals for cubic and tetragonal lattices . the calculations are based on the parameterisation of an effective interaction arising from the exchange of the magnetic fluctuations and assume that a single band is relevant for superconductivity . the results show that for comparable model parameters , the robustness of magnetic pairing increases gradually as one goes from a cubic structure to a more and more anisotropic tetragonal structure either on the border of antiferromagnetism or ferromagnetism .


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: reliable estimates of heavy - light meson decay constants are important , since they appear in many processes from which fundamental quantities can be extracted @xcite . theoretical investigations have focused on estimating decay constants for the weakly decaying pseudo - scalar meson and its hqet ( heavy quark effective theory ) related vector meson . whereas the decay constant of the weakly decaying pseudo - scalar meson is of paramount importance for determining fundamental quantities , the decay constant of the @xmath2-wave vector meson plays a role in exclusive @xmath3 transitions @xcite and in radiative leptonic decays of heavy - light mesons @xcite .    while those decay constants have been and continue to be studied intensively , the decay constants of the more highly excited heavy - light states have been normally ignored . this note attempts to rectify this situation , by predicting decay constants for many higher - excited resonances . that could be important phenomenologically on several accounts . first , cleo recently observed a significant wrong charm contribution in @xmath4 decays @xcite , @xmath5 governed essentially by the @xmath6 quark transitions . the @xmath7 transitions were overlooked in all previous experimental analyses . under the factorization assumption @xcite , wherein the virtual @xmath8 hadronizes independently to the rest of the system , a quantitative modelling of the @xmath7 transitions can be undertaken once theory provides the decay constants for @xmath9 . second , reliable estimates for decay constants of @xmath10 allow one to test whether color - allowed and color - suppressed decay amplitudes interfere constructively for the @xmath11 modes , as has been seen for the @xmath12 transitions @xcite . third , such estimates enable us to better predict subtle cp violating phenomena . decay constants are defined through matrix elements of vector and pseudo - vector currents between meson states and the vacuum . therefore , in order to calculate them , one has to find a way to evaluate hadronic matrix elements . the mock - meson method @xcite has been frequently used in the literature for that purpose @xcite . in this paper we follow the same approach , and use the mock - meson method in order to obtain expressions for the decay constants of heavy - light mesons , in terms of integrals over momentum - space bound - state wave functions . for numerical estimates we decided to use the simplest relativistic generalization of the schr@xmath13dinger equation @xcite , sometimes called the spinless salpeter equation . the rest of the paper is organized as follows . we begin with a brief description of the mock - meson method in section [ mock ] . our approach is based on the @xmath14-@xmath14 coupling scheme , since it is more appropriate for heavy - light mesons than the usual @xmath15-@xmath2 scheme . expressions for the decay constants of heavy - light meson states are given in section [ dcon ] . the relativistic quark model and our numerical estimates are described in section [ nres ] . there we also comment on recent calculations of pseudo - scalar and vector decay constants @xcite . our conclusions are summarized in section [ conc ] . as already mentioned , the mock - meson approach @xcite has been widely used for calculations of hadronic matrix elements @xcite . the basic idea of the method is simple . the mock meson is defined as a collection of free quarks weighted with a bound - state wave function . the mock - meson matrix elements @xmath16 can then be calculated using full dirac spinors . on the other hand , the physical matrix elements @xmath17 can always be expressed in terms of lorentz covariants with coefficients @xmath18 , which are lorentz scalars . in many simple cases , @xmath17 and @xmath16 will be of the same form . the mock - meson prescription then says that in those cases one should simply take @xmath19 . indeed , this correspondence is exact in the zero - binding limit and in the meson rest frame . away from this limit the mock amplitudes are in general not invariant by terms of order @xmath20 .    in this paper we are primarily concerned with the decay constants of heavy - light @xmath21 mesons . in the @xmath22 limit , heavy quark symmetry tells us that the angular momentum of the light degrees of freedom ( ldf ) in the heavy - light meson decouples from the spin of the heavy quark , and both are separately conserved by the strong interaction @xcite . therefore , total angular momentum @xmath14 of the ldf is a good quantum number . for each @xmath14 there are two degenerate heavy meson states ( @xmath23 ) , which can be labeled as @xmath24 , where @xmath25 . this implies that in the case of heavy - light mesons the @xmath14-@xmath14 coupling is more appropriate than the @xmath15-@xmath2 coupling scheme . for this reason , we first define the ldf states @xmath26 as clebsch - gordan ( cg ) combinations of the eigenstates of orbital angular momentum @xmath27 , and those of the spin of the light quark @xmath28 , with cg coefficients denoted as @xmath29 . combining the ldf states with those of the heavy antiquark @xmath30 ( with cg coefficients @xmath31 ) , we get the @xmath21 mock meson state in its rest frame , @xmath32 in the above expression @xmath33 , @xmath34 is the mock - meson mass , and the color wave function ( subscript c denotes color ) is written explicitly . also , @xmath35 is a normalized momentum wave function , where @xmath36 denotes all other quantum numbers of a state not connected to angular momentum ( e.g. , radial quantum number ) . the factor @xmath37 appears due to our normalization convention for creation and annihilation operators @xcite , @xmath38 , etc . the mock - meson states as given in ( [ mmes ] ) are normalized to @xmath39 . as already observed in @xcite , the mock - meson approach suffers from a number of ambiguities , such as the choice for quark masses , or the definition of the mock - meson mass @xmath34 . in the spirit of the method , the mock - meson mass should be defined as @xmath40 . however , as pointed out in @xcite , the mock - meson mass has been introduced to give the correct relativistic normalization of the meson s wave function , and hence the use of the physical meson mass @xmath41 instead of @xmath42 may be more appropriate . we adopt the same approach , and write @xmath43 . we also note that the heavier the mesons are , the less important it is how the mock - meson mass is defined , since the relativistic effects and binding energies become less significant . as far as quark masses are concerned , the self - consistency of the model requires the use of constituent quark masses . in our error estimates we have included variations of constituent light quark masses over a range of about @xmath44 , and also of heavy quark masses over a range of about @xmath45 , so that we believe that uncertainties introduced by a particular choice of quark masses are being properly taken into account . decay constants of heavy - light mesons are defined through matrix elements of vector @xmath46 and pseudo - vector @xmath47 currents between a meson state and the vacuum . following standard definitions in the literature @xcite , for pseudo - scalar ( p ) , vector ( v ) , scalar ( s ) , and pseudo - vector ( a ) mesons , we write @xmath48 note that in the heavy quark limit @xmath49 and @xmath50 states are degenerate ( @xmath2-waves ) , and so are @xmath51 and @xmath52 ( @xmath0-wave states ) . the spin 2 members of @xmath0-wave ( @xmath53 ) and @xmath1-wave ( @xmath54 ) doublets do not couple leptonically due to conservation of angular momentum .    in order to obtain expressions for decay constants in terms of integrals over momentum - space meson wave functions , we evaluate the matrix elements ( [ fps])-([fv32 ] ) in the meson rest frame using ( [ mmes ] ) . of course , any choice of polarization for spin 1 mesons should yield the same result . as mentioned earlier , current matrix elements between states defined in ( [ mmes ] ) and the vacuum can be evaluated exactly with full dirac spinors . because of spherical symmetry , the momentum - space wave function can be written in the form @xmath55 in the above @xmath56 are the usual spherical harmonics , and @xmath57 is the radial part of the wave function , where @xmath58 denotes @xmath59 henceforth . using ( [ wf ] ) , and keeping track of the relevant cg coefficients , we find that all heavy - light meson decay constants in the mock - meson approach can be written in the form @xmath60 where @xmath61 r_{n0}(p)\ , \label{fi1}\\ f_{v_{1/2}}(p ) & = & \left [ 1 + \frac{1}{3 } \frac{p^{2}}{(m_{q}+e_{q})(m_{\overline{q}}+e_{\overline{q } } ) } \right ] r_{n0}(p)\ , \label{fi2}\\ f_{s}(p ) & = & \left [ \frac{1}{(m_{q}+e_{q } ) } - \frac{1}{(m_{\overline{q}}+e_{\overline{q } } ) } \right]p r_{n1}(p)\ , \label{fi3}\\ f_{a_{1/2}}(p ) & = & \left [ \frac{1}{(m_{q}+e_{q } ) } + \frac{1}{3 } \frac{1}{(m_{\overline{q}}+e_{\overline{q } } ) } \right]p r_{n1}(p)\ , \label{fi4}\\ f_{a_{3/2}}(p ) & = & \left[\frac{2\sqrt{2}}{3 } \frac{1}{(m_{\overline{q}}+e_{\overline{q } } ) } \right]p r_{n1}(p)\ , \label{fi5}\\ f_{v_{3/2}}(p ) & = & \left[\frac{2\sqrt{2}}{3 } \frac{1}{(m_{q}+e_{q } ) } \frac{1}{(m_{\overline{q}}+e_{\overline{q } } ) } \right]p^{2 } r_{n2}(p)\ . \label{fi6}\end{aligned}\ ] ] expressions ( [ fi1 ] ) and ( [ fi2 ] ) were found in @xcite and @xcite , respectively . it is interesting to observe that in the limit @xmath22 ( [ fi])-([fi6 ] ) become @xmath62 with @xmath63 equality of @xmath64 and @xmath65 , and also that of @xmath66 and @xmath67 , as well as vanishing of @xmath68 and @xmath69 , are in the heavy quark limit expected from the heavy quark symmetry . in order to obtain numerical estimates for the decay constants of heavy - light mesons , we consider the simplest and widely used generalization of the non - relativistic schr@xmath13dinger equation @xcite with hamiltonian given by @xmath70 where for @xmath71 we take the qcd - motivated coulomb - plus - linear potential @xcite @xmath72 for the sake of simplicity , we take @xmath73 to be a fixed effective short range coupling constant . the effective string tension of the model can be determined from the requirement that the linear regge structure of the model in the light - light limit agrees with the observed slope of the @xmath74 trajectory @xcite . fixing @xmath75 , other parameters can be chosen so that the model reproduces the observed spin - averaged spectrum of the known heavy - light states . one such set of parameters includes constituent quark masses @xmath76 , @xmath77 , @xmath78 , and @xmath79 , and also @xmath80 , @xmath81 , and @xmath82 . and @xmath83 mesons ( @xmath84 and @xmath85 states ) of about @xmath86 . ] as can be seen in table [ spectrum ] , these parameters yield an excellent description of the observed spin - averaged heavy - light spectrum . we now turn to the discussion of pseudo - scalar and vector decay constants . recently , ref . @xcite used ( [ srqm ] ) with six different potentials , and with current quark masses from @xcite , minimized the hamiltonian with respect to the variational parameter @xmath87 of a single harmonic oscillator ( ho ) wave function , @xmath88 and then used the wave function obtained in this way to get pseudo - scalar and vector decay constants from ( [ fi ] ) , ( [ fi1 ] ) and ( [ fi2 ] ) . however , a single harmonic oscillator ( ho ) basis state is not a suitable approximation for the meson wave function . namely , lattice simulations @xcite show that heavy - light wave functions fall exponentially with large @xmath89 ( @xmath90 ) , and therefore ho wave functions ( @xmath91 ) can not be expected to reflect the correct dynamics of heavy - light mesons . if single basis states are used , a much better choice would be pseudo - coulombic ( pc ) basis states @xcite which fall exponentially with large @xmath89 and appear to be in a good agreement with the lattice data , as can be seen in figure [ wfplt ] . models such as the one we are using here are usually solved by diagonalizing the hamiltonian matrix in a particular ( truncated ) basis , with basis states depending on some variational parameter @xcite . as one increases the number of basis states , the dependence of eigenvalues and eigenfunctions on the variational parameter should vanish for the lowest states . in the case of qcd - motivated potentials the solutions obtained with the pc wave functions converge much more rapidly with an increase in the number of basis states , than those obtained with the ho wave functions . we illustrate that in figure [ bplt ] , by plotting the dependence of energy of the lowest @xmath92 state on the variational parameter for @xmath93 , 5 and 15 basis states , for both pc and ho wave functions . one can clearly see that the lowest @xmath92 ho wave function is not a very good trial wave function in a variational calculation of ( [ srqm ] ) ( with qcd - motivated potentials ) . furthermore , even if one believes that the @xmath93 ho result for a state energy is acceptable ( it is roughly @xmath94 higher than the exact solution , as can be seen in figure [ bplt ] ) , that still does not justify the use of a single ho basis state as a meson wave function . this issue is clearly important in calculations where a correct description of meson dynamics is needed , such as calculations of meson decay constants . results obtained by varying a single ho basis state are thus to be interpreted as non - relativistic estimates of some effective harmonic oscillator potential , and not as the results of a qcd - motivated relativistic quark model . one can now observe that if one uses enough basis states , the choice of basis wave functions should not matter , and pseudo - scalar and vector decay constants should be obtainable from the relativistic quark model considered here . the problem is , however , that the @xmath92 wave function is divergent at the spatial origin @xcite , i.e. , @xmath95 the singularity for @xmath96 is related to the singularity of the short - range coulomb potential . by increasing the number of ( usually finite at @xmath97 ) basis states , one is gradually beginning to see that singularity @xcite . furthermore , from ( [ div ] ) one can see that the degree of divergence highly depends on the choice of @xmath73 . because of that , one can expect that pseudo - scalar and vector decay constants can not be reliably estimated within the model we are considering . in figures [ fpsplt ] and [ fvplt ] we demonstrate the dependence of the pseudo - scalar ( @xmath1-meson ) and vector ( @xmath98-meson ) decay constants on the number of basis states ( @xmath99 ) , for both pc and ho wave functions . as one can see , for small @xmath99 both @xmath64 and @xmath65 are significantly increasing with an increase in @xmath99 . by including enough basis states , the dependence on @xmath99 would eventually vanish . ) the results for @xmath64 are better behaved than those for @xmath65 . for example , @xmath64 obtained with @xmath100 pc states are usually larger than those obtained with @xmath101 by only a few @xmath102 . on the other hand , the same increase in @xmath99 in general leads to increase in @xmath65 by several hundred @xmath102 . ] however , as implied by ( [ div ] ) , both @xmath64 and @xmath65 are quite sensitive to the particular choice of parameters of the model . in our calculations we have observed that results obtained with fixed @xmath99 can vary up to a few hundred @xmath102 . because of that , we were not able to obtain reliable estimates of @xmath64 and @xmath65 from the model considered in this paper . and [ fvplt ] it should be clear that in the model considered here the ratio @xmath103 also can not be determined with reasonable errors . ]    one possible solution of the problem discussed above would be to replace the @xmath104 potential with the one - loop single gluon exchange potential , i.e. , @xmath105 . the @xmath92 solution of ( [ srqm ] ) in that case is still divergent , but the divergence is only logarithmic @xcite . this should lead to much more stable results than the ones shown in figures [ fpsplt ] and [ fvplt ] . these results should also be much less dependent on the specific choice of the model parameters . in fact , such a calculation for @xmath64 ( for @xmath1 , @xmath106 , @xmath4 , and @xmath107 mesons ) was already performed by capstick and godfrey in @xcite using the model of @xcite . the dependence of their results on the number of basis states was not shown , but the authors of @xcite stated that they believed that the model overestimates pseudo - scalar decay constants ( e.g. , for @xmath1 meson they found @xmath108 with uncertainty of 20% ) . even though it is important to investigate what really happens with both @xmath64 and @xmath65 in such a model , we shall not consider it in the present paper . we next discuss the heavy - light @xmath0 and @xmath1-wave decay constants . while we were not able to obtain reliable results from ( [ srqm ] ) and ( [ pot ] ) for the @xmath2-waves , the situation for @xmath0 and @xmath1-waves is completely different . in figures [ fsplt ] , [ fpvplt ] , [ fpv2plt ] and [ fv2plt ] we show the dependence on the number of basis states ( @xmath99 ) , for scalar ( @xmath2 ) , two pseudo - vector ( @xmath109 and @xmath110 ) , and vector ( @xmath111 ) decay constants , respectively . all the results shown are for the @xmath10 mesons . as one can see in those figures , in general only a few basis states are needed for results to become independent of @xmath99 , even though the derivatives of the actual @xmath112 and @xmath113 wave functions are singular at spatial origin @xcite . mev for @xmath66 and @xmath69 as little as 10 pc basis states usually were needed , while to achieve the same accuracy for @xmath67 and @xmath68 requires in general about 50 to 75 pc basis states . ] furthermore , as @xmath99 increases the ho results approach the pc results ( always from below ) which shows that the difference between the two basis sets is slowly vanishing . however , even with 15 basis states ( when the state energy obtained from the model is essentially equal for both pc and ho wave functions ) , we can still see the difference for @xmath67 ( figure [ fpvplt ] ) and for @xmath68 ( figure [ fpv2plt ] ) . this reflects the difference in the wave functions obtained from the two basis sets . the reason why both pc and ho basis states yield almost the same results for @xmath66 ( figure [ fsplt ] ) , even though @xmath84 state is also a @xmath0-wave , is the minus sign in ( [ fi3 ] ) . of course , because of the much more rapid convergence , the @xmath114 results are to be preferred over the ho results . our calculations of @xmath0 and @xmath1-wave decay constants showed that their dependence on the particular choice of the model parameters is significantly smaller than the corresponding dependence of @xmath66 and @xmath65 . we present the results for @xmath10 , @xmath9 , @xmath115 , and @xmath116 mesons in tables [ td ] , [ tds ] , [ tb ] , and [ tbs ] , respectively . through [ tbs ] were obtained with 25 pc basis states , which was more than enough for the accuracy of less than @xmath117 in all cases considered . ] to obtain these results the effective string tension @xmath118 of the model was determined from the observed slope of the @xmath74 trajectory . for a fixed @xmath75 other parameters were obtained from the spectrum of the known heavy - light states . experimental meson masses were used in ( [ fi ] ) only when their quantum numbers were unambiguously determined . else , we used model predictions for the appropriate spin - averaged masses , which are also shown in tables [ td ] through [ tbs ] .    in order to estimate uncertainties introduced by a particular choice of the constituent mass of @xmath119 and @xmath120 quarks , we have varied @xmath75 in the range from @xmath121 to @xmath122 . for a given @xmath75 , by adjusting @xmath123 we have also varied constituent heavy quark masses in the range of about @xmath45 ( e.g. , @xmath124 was varied in the range from about @xmath125 to about @xmath126 ) . we emphasize that a good description of the spin - averaged heavy - light meson spectrum was always maintained . results for the decay constants obtained in this way depend on the assumption for the unknown spin - averaged mass of @xmath127 and @xmath83 mesons ( @xmath84 and @xmath85 states ) . to take into account ambiguities introduced in our results in that way , we have repeated all calculations for this unknown mass in the range from @xmath128 to @xmath129 . errors quoted in tables [ td ] through [ tbs ] reflect the uncertainty due to the unknown @xmath0-wave mass , as well as the uncertainties related to the choice of constituent quark masses discussed above .    as one can see from those tables , in spite of the fact that our calculations are performed for a broad range of constituent quark masses , and also for a wide range of the unknown @xmath0-wave mass , as long as a good description of the observed heavy - light meson spectrum is maintained , the @xmath0 and @xmath1-wave heavy - light decay constants are all predicted rather precisely . it is also interesting to observe that the decay constants of strange @xmath51 , @xmath52 , and @xmath130 states are slightly smaller than those of the corresponding non - strange states . the main reason is ( besides the @xmath131 dependence of ( [ fi ] ) ) the light quark dependence of ( [ fi3 ] ) , ( [ fi4 ] ) and ( [ fi6 ] ) . on the other hand , ( [ fi5 ] ) does not depend on the light quark mass , so that @xmath132 factor in ( [ fi ] ) plays a much more significant role , and as a result @xmath68 for the strange states are larger than the ones for non - strange states . also note that @xmath66 for @xmath133 and @xmath134 are larger than those of the corresponding @xmath127 and @xmath135 states , while it is the other way around in the case of @xmath67 . the reason for this are the minus and plus signs in ( [ fi3 ] ) and ( [ fi4 ] ) , respectively . finally , the fact that @xmath136 and @xmath130 @xmath115 states have decay constants smaller than those of the corresponding @xmath10 states , can be easily explained with the @xmath137 dependence of ( [ fi5 ] ) and ( [ fi6 ] ) . in this paper we have examined decay constants of heavy - light mesons within the mock - meson approach @xcite . we obtained all the relevant expressions in the @xmath14-@xmath14 coupling scheme . for numerical estimates we employed a simple and widely used relativistic quark model @xcite . it is based on a spinless salpeter equation with qcd - motivated coulomb - plus - linear potential . the effective string tension is chosen so that the regge structure of the model in the light - light limit is consistent with experiment , and other parameters are based on the good description of the known spin - averaged heavy - light meson masses .    due to the singular nature of the @xmath138 wave functions at spatial origin @xcite , we were not able to obtain reliable estimates of pseudo - scalar and vector decay constants . on the other hand , even though we have allowed for large variations of input parameters , our results show that the model predicts a rather narrow range for all lowest @xmath0 and @xmath1-wave heavy - light decay constants . such precisely predicted decay constants allow us to estimate the @xmath139 production fractions in @xmath118 decays governed by the @xmath6 transitions under the factorization assumption . quantitative predictions regarding the interference of color - allowed and color - suppressed amplitudes in @xmath11 modes can now be formulated . these and some other consequences of our findings will be discussed elsewhere @xcite . we thank j. f. amundson , e. eichten and d. zeppenfeld for discussions . s. veseli would also like to thank the theory group for hospitality during his visit to fermilab . this work was supported in part by the u.s . department of energy under contracts no . de - ac02 - 76ch03000 and de - fg02 - 95er40896 , and in part by the university of wisconsin research committee with funds granted by the wisconsin alumni research foundation . see , for instance , n. isgur and m. b. wise , phys . d * 41 * , 151 ( 1990 ) ; m. b. wise , phys . d * 45 * , 2188 ( 1992 ) ; g. burdman and j. f. donoghue , phys . b * 280 * , 287 ( 1992 ) ; g. burdman and j. kambor , fermilab report , fermilab - pub-96/033-t ( hep - ph/9602353 ) . .relativistic quark model predictions compared to experimental spin - averaged heavy - light meson masses . parameters of the model are @xmath76 , @xmath77 , @xmath78 , @xmath79 , @xmath80 , @xmath81 , and @xmath82 . the unknown @xmath127 and @xmath83 mesons ( @xmath84 and @xmath85 states ) were assumed to have a spin - averaged mass of 2400@xmath142 . heavy quark symmetry arguments then lead to the spin - averaged mass of 2502 @xmath102 for the corresponding @xmath135 and @xmath143 mesons . [ cols="<,^,^,^,^ " , ]
Abstract:

Answer: we investigate decay constants of @xmath0 and @xmath1-wave heavy - light mesons within the mock - meson approach . numerical estimates are obtained using the relativistic quark model . we also comment on recent calculations of heavy - light pseudo - scalar and vector decay constants .    6.1 in 8.5 in .25 in