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 hubbard model@xcite of interacting fermions in a solid is a centerpiece of modern condensed matter physics . it is conventionally defined by restricting the motion of electrons in a crystalline solid to a single band , and simplifying the screened long - range coulomb interactions between electrons to an on - site repulsion :    @xmath0    where @xmath1 , @xmath2 creates a fermion on a site @xmath3 , @xmath4 , the summation in the first term runs over the nearest - neighbor sites of the simple cubic lattice , @xmath5 is the hopping amplitude and @xmath6 is the onsite repulsion . though remarkably simple in appearance , the model has been used to study a wealth of intriguing quantum many - body phenomena that are due to electron correlations in solids , such as interaction - driven insulators @xcite , quantum magnetism and high - temperature superconductivity @xcite . however , despite more than a half of century of intensive investigation , the physics of the model is still not completely understood . the most challenging yet the most interesting regime is the intermediate regime with interaction comparable to half the bandwidth @xmath7 , where @xmath8 is the number of nearest neighbors of a site . this regime offers no small parameter to start a controllable analytic theory . furthermore , exact analytic solutions are only accessible either in one spatial dimension @xcite or an infinite@xcite number of spatial dimensions . substantial progress has been possible with the development of efficient quantum monte carlo ( qmc ) methods ( for a recent review , see ref .  ) accompanied by advances in computer technology . although for generic _ bosonic _ systems virtually any equilibrium property can nowadays be calculated by qmc with a controlled high accuracy @xcite , systematic - error - free simulations of correlated _ fermions _ have been limited to a handful of special cases due to the negative sign problem @xcite . the sign problem manifests itself as an exponential scaling of the simulation time with the system volume and inverse temperature making it practically impossible to obtain any reliable information about the system in the thermodynamic ( td ) limit . although in some cases the sign problem can be completely eliminated by choosing a system - specific representation , its general solution is almost certainly not possible @xcite . a major step toward understanding strongly correlated systems has been the experimental realization of the hubbard model with ultracold atomic gases loaded into optical lattices ( for recent reviews see refs . ) . these systems offer substantial control over the hamiltonian . as a result , these experiments can serve as emulators of quantum many - body systems , which allow the accurate study of a given model in a range of parameters inaccessible by analytic and numeric techniques@xcite . the recent experimental observation of mott physics @xcite in the hubbard model is a major milestone along these lines . the next crucial step would be a realization of the antiferromagnetic ( afm ) transition and the nel state in the hubbard model , which requires a substantial effort in reaching lower temperatures , controlling the equilibration rates , developing new probing techniques , etc . in addition to these inherent challenges , there is also a fundamental problem related to thermometry in ultracold - atom systems since they are insulated from the environment . isolated systems such as these are thus , by their nature , better characterized by entropy , rather than temperature . moreover , probes have to be calibrated in the relevant regime and the results obtained with the setup validated against available benchmarks . for these purposes , reliable and accurate numeric results for fermionic systems are indispensable because they ultimately allow a full quantitative understanding , as was recently demonstrated by the example of a bosonic optical - lattice emulator @xcite . this work provides reliable benchmarks for the realization of the nel state in optical lattices as well as for new theoretical methods . we focus on the special case of the half - filled ( @xmath9 , or equivalently @xmath10 ) three dimensional ( 3d ) hubbard model on a simple cubic lattice . the case of half filling is special due to the @xmath11 symmetry of the hamiltonian , which is ultimately broken by the nel state , making magnetism the leading instability . in the limit @xmath12 , the effects of the interaction can be studied perturbatively . in the opposite limit , @xmath13 , eq . reduces to the afm heisenberg model with @xmath14 . while there is no doubt that in the strongly correlated regime the ground state of the half - filled model is afm , mapping out the finite - temperature phase diagram and studying thermodynamics of the system is extremely challenging @xcite . to this end , we employ the unbiased continuous - time determinantal diagrammatic monte carlo ( ddmc ) technique @xcite , which produces numerically exact ( up to a known statistical error bar ) results for a finite - size system , and which is free from the fermionic sign problem at half filling on bipartite lattices , allowing a reliable extrapolation of results to the thermodynamic limit . we use ddmc to determine several critical properties of the hubbard model with high control and accuracy and compare with high - temperature series expansions@xcite ( htse ) where possible . we study the range of on - site interaction @xmath15 , where the critical temperature @xmath16 of the afm transition is expected to reach its maximum @xcite . results in this regime are vital to optical lattice emulator efforts offering experiments their best chance of observing the afm phase . we obtain the critical temperature @xmath16 and compute important thermodynamic properties of the model  the energy and the entropy , as well as two optical lattice observables : double occupancy and the nearest - neighbor spin - spin correlation function  in the paramagnetic phase as a function of temperature down to @xmath16 . a number of previous works using different unbiased approaches studied @xmath16 ( by the determinant quantum monte carlo ( dqmc)@xcite and the dynamical cluster approximation ( dca)@xcite ) and the thermodynamic properties in question ( by ddmc@xcite , a combination of dca and ddmc@xcite , and by dqmc@xcite ) in this regime . our work improves the accuracy of the previous results at half filling and provides the most accurate to date values of the critical temperature @xmath16 and the entropy at the critical point @xmath17 ( summarized in the table . [ table : tns ] ) . accurate knowledge of @xmath16 is particularly important for determining the critical entropy @xmath17 since @xmath18 is a steep function near the transition , so that the error bar of @xmath17 is mainly due to the uncertainty in @xmath16 . the values of @xmath17 are required for an experimental realization for the afm state . as was noted in ref . , close to the transition , the temperature dependence of the nearest - neighbor spin - spin correlation function @xmath19 ( @xmath20 is the unit vector in the direction @xmath21 ) is significantly more pronounced than that of the double occupancy @xmath22 of a lattice site making measurements of the nearest - neighbor spin - spin correlations @xcite more suitable for accurate thermometry in this regime . our results for the spin - spin correlation function can be used for calibration of such a thermometer .    for over a decade , the dqmc study of the phase diagram of the half - filled 3d hubbard model by staudt _ @xcite has been the main reference for @xmath16 in the correlated regime . representing the state of the art at that time , ref . provides a comprehensive comparison of the dqmc data for the nel temperature with those of preceding qmc simulations and approximate theories , e.g. , dmft @xcite . we shall not reproduce this comparison here and refer the reader to a review of results for @xmath16 predating ref .  . the simulation method of staudt _ et al . _ is based on a discrete hubbard - stratonovich decoupling in the hubbard interaction term which requires discretizing the imaginary - time interval @xmath23 into a finite number of steps of size @xmath24 thereby introducing a systematic error . hence , in addition to the standard extrapolation of the finite - size dqmc data to the td limit , one has to perform an extrapolation with respect to @xmath25 . such a double extrapolation is rather laborious . in practice @xmath24 is usually fixed at a value which is large enough to allow efficient simulations , yet , according to ref . , such that `` the results are not significantly affected by the extrapolation @xmath25 '' . in the absence of an explicit extrapolation , the degree of control over systematic errors can be questioned . this is where our approach is a major improvement over that of staudt _ et al . _ the ddmc technique is formulated directly in continuous imaginary time . therefore finite - size corrections are the only systematic error we have to eliminate in our approach . the cost of the absence of the additional systematic error is the computation complexity of the ddmc , which scales as @xmath26 ^ 3 \omega^3 $ ] versus the linear in @xmath27 \omega^3 $ ] scaling of the dqmc of ref .  , for lattices with @xmath28 sites . as a result , we are limited by the values of interaction @xmath29 , whereas staudt _ _ were able to study the model up to @xmath30 . we have also been able to identify the transition itself with a considerable improvement in accuracy over ref .  . long - range afm order in the system causes divergence of the magnetic structure factor @xmath31 where @xmath32 , and @xmath33 is the afm wave vector , so that @xmath34 is related to the magnetization @xmath35 in the td limit , @xmath36 . in ref . , the transition temperature is found as the point at which @xmath35 starts to noticeably deviate from zero , while the magnetization is obtained from a finite - size extrapolation of @xmath37 with respect to @xmath38 . here the scaling power of the finite - size corrections was used as a fitting parameter . an additional ( indirect ) probe of the transition used by staudt _ _ is the peak in the dependence of the specific heat on temperature . in contrast , we use a much more accurate method to determine the critical point . we find the critical point using a finite - size scaling analysis of @xmath39 in combination with the technique of binder crossings @xcite . this approach allows us to get a reliable and accurate value of @xmath40 by making use of the known scaling law of the magnetic structure factor at the critical point , @xmath41 , where @xmath42 is the anomalous dimension of the particular universality class , so that the quantity @xmath43 becomes scale - invariant at the critical temperature up to non - universal corrections , which we also take into account . since the nel transition is breaking the @xmath11 symmetry of the hubbard hamiltonian at half - filling , the universality class is that of the 3d heisenberg model , which provides the critical exponents for the finite - size scaling analysis .    in the context of previous calculations of the nel temperature , it is worth noting the work by kent _ @xcite , where @xmath16 was found from calculations using dca . although this work did not lead to any improvement of the precision claimed by staudt _ , it was shown that the dca can be used to determine the critical point with a controlled accuracy based on significantly smaller clusters of only @xmath44 lattice sites versus up to @xmath45 sites in ref .   and in this work . in addition , good agreement with the values of @xmath16 of ref .   suggested that potential systematic error of the @xmath25 extrapolation in ref .   is likely to be small . in the range of interaction @xmath46 , our results for @xmath16 agree within the error bars with those of refs .   implicitly confirming this . at smaller @xmath47 , however , we find somewhat lower values of @xmath16 . moreover , at @xmath48 , being unable to reach significantly low temperatures to accurately infer the critical point , we can only claim an upper bound from our finite - size - scaling analysis , @xmath49 , which is already lower than the values claimed in refs .  . the reason for the discrepancy is likely to be the long - range character of correlations at smaller coupling , which can be missed in simple finite - size extrapolation schemes based on data for insufficiently large systems . the thermodynamics of the hubbard model near the nel transition in connection with its experimental realization has been the focus of a number of recent studies @xcite . the dmft results @xcite emphasize the role of double occupancy in detecting the build up of afm correlations . however , its dependence on temperature near the nel transition is relatively flat in the regime where @xmath16 is maximal , as observed in unbiased ( extrapolated to the td limit ) dca calculations @xcite . @xcite obtained the energy and the entropy down to @xmath16 in this regime as well as the equation of state away from half filling , which allowed to get an estimate of the entropy at the transition in a realistic harmonically trapped system . a study@xcite of the system using the same dqmc method as ref .   arrived at an agreement with the results and conclusions of fuchs _ et al._.    the ddmc simulation method used in this work is not capable of capturing thermodynamics of the hubbard model away from half filling due to the pronounced negative sign problem . however , exactly at half filling , it has certain advantages over dca and dqmc . as discussed above , in order to claim unbiased results in the td limit within dqmc one has to resort to a double extrapolation , @xmath50 . ddmc is formulated directly in continuous time allowing a more reliable extrapolation to the td limit . modern efficient solvers for dca also work in continuous time , but the cluster sizes amenable to simulation in the regime of interest are typically less then @xmath51 . however , in dca the clusters are embedded in a self - consistently defined medium , which largely improves the convergence to the td limit . in practice the finite - size dependence for the accessible clusters is notably larger than that in ddmc @xcite . hence , at half filling and @xmath52 , where @xmath16 and @xmath17 are expected to reach their maxima , ddmc currently allows to obtain the most reliable benchmarks for the thermodynamic quantities of interest . the paper is organized as follows . in sec . [ sec : method ] , we discuss the simulation method outlining the general formulation of the ddmc technique and its application to calculating the specific observables in question . section  [ sec : t_n ] is concerned with determining the temperature of the nel transition . [ sec : thermodynamics ] describes the thermodynamics near @xmath16 . here we discuss the extrapolation of the observables to the td limit ( [ subsec : td - limit_extrapolation ] ) , the determination of entropy ( [ subsec : entropy ] ) , and thermometry near @xmath16 ( [ subsec : thermometry ] ) . we summarize the results in sec . [ sec : conclusions ] . the appendix contains tables of the obtained numerical data for the entropy , energy , double occupancy , and spin - spin correlation functions as a function of temperature . we first rewrite the hubbard hamiltonian to a form suitable for numerical simulations by mapping the repulsive model to an attractive model by a particle - hole transformation @xcite . we use the fact that the simple cubic lattice is bipartite and can be split into two interpenetrating sublattices @xmath53 and @xmath54 , so that the hopping term in only connects sites belonging to different sublattices . then we introduce the hole operators for the @xmath55-component : @xmath56 this way , eq .   becomes at half filling @xmath57 where @xmath58 , @xmath59 is the number operator for the attractive model , @xmath60 as appropriate for half filling , and @xmath28 is the total number of sites .    since we only consider half filling , @xmath61 , we follow ref .   and shift the chemical potential according to @xmath62 : @xmath63 at half filling the choice @xmath64 is optimal because it leads to the minimal computational complexity of the simulations ( see below ) , and we use this value of @xmath65 throughout .      the ddmc algorithm works with the weak - coupling expansion series for the finite - temperature partition function for the hubbard model . the latter reads @xmath66 here @xmath67 is the inverse temperature , @xmath68 is the unperturbed partition function , @xmath69 denotes time ordering , and @xmath70/z_0 $ ] denotes the thermodynamic average with respect to the unperturbed hamiltonian @xmath71 . equations  generate the standard feynman diagrams : there are @xmath72 diagrams of order @xmath73 , which can be represented graphically as a collection of @xmath73 vertices connected by single - particle propagators . summing over all the interconnections for a fixed vertex configuration @xmath74 equation takes the form @xmath75 where @xmath76 are @xmath77 matrices with matrix elements given by ( @xmath78 ) @xmath79 ( since we only consider an unpolarized system , we omit the spin index for @xmath80-s and @xmath81-s ) , and @xmath81 being the free - particle green s functions , @xmath82 since @xmath83 , eq .  , is diagonal in momentum space , @xmath84   a^\dagger_{{\mathbf{k}}\sigma } a_{{\mathbf{k}}\sigma } , \nonumber \\ \varepsilon_{\mathbf{k } } = - 2t\sum_{i=1,2,3 } \cos(k_i ) , \label{h_0_k}\end{aligned}\ ] ] the free propagators are calculated by the fourier transform @xmath85 } , \ ; \tau>0 , \label{g_0_ft } \\ g^{(0 ) } ( { \mathbf{k } } , -\tau)=-g^{(0 ) } ( { \mathbf{k } } , \beta-\tau ) , \nonumber\end{aligned}\ ] ] and tabulated before the start of the simulation .    the series ,  serves as a basis for a ddmc simulation of the hubbard model : we set up a random walk in the space of the vertex configurations @xmath86 , eq .  , using the metropolis algorithm @xcite with the weights proportional to @xmath87 , eq .  . since the technique itself is detailed elsewhere @xcite , here we only briefly discuss the specific details of the present implementation . we only stress at this point that since all the terms in the series ,  are positive definite we completely avoid a sign problem . the simplest updating strategy for ddmc simulations consists of adding ( @xmath88 ) and removing ( @xmath89 ) interaction vertices at random positions in @xmath3 and @xmath90 to / from a vertex configuration @xmath86 . however , at half filling and with @xmath64 the series only contains even order terms , hence we employ rank-2 updates , where @xmath91 . using the woodbury - type formulas , both rank-1 and rank-2 updates can be performed in @xmath92 operations . we note that for @xmath93 both even and odd terms are present even at half filling . in this sense the choice of @xmath64 is optimal . the general method for calculating observables in the ddmc simulations uses the standard technique of monte carlo estimators : for an observable @xmath94 we define an estimator @xmath95 such that the average of the latter over the vertex configurations generated by the mc process , @xmath96 , converges to the thermal average @xmath97 , where @xmath98 \;.\ ] ] below we explicitly list estimators for useful observables :    [ [ filling - fraction ] ] filling fraction + + + + + + + + + + + + + + + +    the thermal average of the filling fraction @xmath99 for the spin projection @xmath100 is given by @xmath101 hence the corresponding estimator is@xcite @xmath102 where @xmath103 is a @xmath77 matrix and @xmath104 is a @xmath105 matrix with an extra row and a column corresponding to the extra creation and annihilation operators in . here @xmath3 and @xmath90 are random positions in space and time , respectively . notice that for a half filled model and for @xmath64 , the mc average of the second term in must equal zero . we use this fact to check if a simulation has equilibrated . [ [ kinetic - energy ] ] kinetic energy + + + + + + + + + + + + + +    calculating the kinetic energy for the model requires evaluating the average of @xmath106 , which differs from eq .   only in that the creation annihilation is shifted in space with respect to the creation operator while in eq . both operators reside on the same lattice site . the estimator for the kinetic energy @xmath107 is then@xcite @xmath108 where the matrix @xmath109 differs from @xmath104 by the fact that the creation operator in is shifted in space with respect to the annihilation operator . the extra factor of 2 accounts for the summation over @xmath110 and @xmath111 is the number of bonds of a lattice with @xmath112 sites and periodic boundary conditions ( pbc ) . [ [ interaction - energy ] ] interaction energy + + + + + + + + + + + + + + + + + +    since the series is nothing but an expansion in powers of @xmath113 , the corresponding estimator is readily obtained by a standard trick of considering the hamiltonian @xmath114 and differentiating with respect to @xmath115 . the result is@xcite @xmath116    [ [ afm - structure - factor ] ] afm structure factor + + + + + + + + + + + + + + + + + + + +    calculating the afm structure factor in the particle - hole transformed model requires calculating two independent equal - time density - density correlation functions @xmath117    the estimator for the equal spin density correlation function @xmath118 for @xmath119 is given by @xmath120 where @xmath121 is a @xmath122 matrix with two extra rows and columns corresponding to the extra creation and annihilation operators in eq .  . for @xmath123 , @xmath124 which equals @xmath125 for a half - filled model .    to build the estimator for eq . we proceed similarly to . the resulting estimator for is    @xmath126    where @xmath127 is a @xmath77 matrix constructed by selecting a random vertex , @xmath128 , from a configuration @xmath86 and moving the corresponding row of the matrix @xmath103 to @xmath129 while leaving the corresponding column at @xmath128 . in the paramagnetic phase , @xmath130 , @xmath131 scales to zero exponentially as @xmath132 . in the afm phase , on the other hand , @xmath133 as @xmath132 , where @xmath35 is the sublattice magnetization . right at the critical temperature @xmath134 , where @xmath135 is the anomalous dimension .    in order to locate the transition temperature we thus use the standard finite - size scaling ( fss ) ansatz @xcite @xmath136 where @xmath137 is the correlation length which diverges at the transition as @xmath138 , @xmath139 is a real - valued function which tends to a finite constant as @xmath140 , and the corrections in brackets arise from the leading irrelevant operators ( dots represent the higher - order corrections ) . here the exponent @xmath141 is universal , but the amplitude @xmath142 is not . in accordance with the 3d heisenberg universality class , we take @xmath143 , @xmath144 and @xmath145 @xcite .    the basic idea of using eq . for the fss is as follows : if the corrections - to - scaling ( the 2nd term in brackets in eq .  ) were not present , @xmath146 would be scale independent at the transition point , so that performing the simulations at a series of system sizes @xmath147 and plotting @xmath146 versus temperature , one would observe that all the curves intersect at the same point , @xmath148 . this is what we indeed observe for ( relatively ) large values of @xmath149 at @xmath150 : for @xmath151 our mc results are consistent with @xmath152 in eq .   within statistical errors , see fig . [ fig : tc_u8 ] . the data for the smallest systems of size @xmath153 systematically deviate from the scaling described by eq .   for all values of @xmath47 considered here ( see also below ) , and hence they are omitted from the scaling analysis . we find that the corrections - to - scaling become more pronounced with decreasing @xmath149 : for @xmath154 we find a clear evidence of the shift of the pairwise crossings towards lower temperatures , see fig . [ fig : tc_u5 ] . since simulating larger system sizes is not an option , we employ eq .   including corrections to scaling . the most straightforward way is to follow the evolution of the pairwise crossings with the system size . indeed , expanding eq .   around the crossing of @xmath155 at system sizes @xmath156 and @xmath157 up to the terms linear in @xmath158 we find ( cf ref .  ): @xmath159 we perform the linear fit of the series of crossings @xmath160 versus @xmath161 . then the intercept of the best - fit line yields the nel temperature . this procedure is illustrated in fig . [ fig : omegau5 ] . it is clear from figs . [ fig : tc_u8][fig : omegau5 ] that @xmath153 does not follow the scaling . we attribute it to the effect of the higher - order terms neglected in , and only use @xmath162 in the fitting procedure . we note at this point that the methodology based on binder crossings has a build - in self - consistency check : if , in fact , the criticality were not in the heisenberg universality class , the curves for the magnetic structure factor , rescaled via eq .   would have no reason to cross at a unique point , and thus the whole procedure of  would break down . an equivalent procedure has been suggested in ref .  . again , one expands eq .   up to the linear order in @xmath158 , which leads to @xmath163 which is then used as a four - parameter ansatz for a single nonlinear fit . _ a priori _ , fitting procedures based on and are equivalent and indeed produce consistent results , as illustrated in fig . [ fig : tc_u5 ] . we stress at this point that using requires judicious choice of the temperature range for fitting : including monte carlo points at too high temperatures and/or too small system sizes tends to significantly skew the fit results . in the following we therefore quote the @xmath16-s obtained using eq .  . for @xmath164 , we find the corrections - to - scaling to be larger than those for @xmath154 , see fig . [ fig : tc_u4 ] . in fact , with the accessible systems sizes we are only able to put an upper limit on the nel temperature , @xmath165 . from fig .  [ fig : tc_u4 ] it is clear that @xmath166 and possibly even @xmath167 are simply too small and need to be discarded from the finite - size scaling analysis . our results for the dependence of the nel temperature on @xmath47 are summarized in table  [ table : tns ] and fig . [ fig : tns ] . it is instructive to compare our estimates to the previous unbiased calculations from the literature . while for @xmath168 and @xmath169 our estimates agree with and are more accurate than previous estimates from qmc @xcite and dca @xcite . for smaller values of @xmath149 our estimates are systematically lower . the discrepancy can be traced back to the fss procedure which includes corrections - to - scaling , eq .  : if we were to discard the corrections and identified the binder crossings of @xmath166 and @xmath167 with the nel temperature , such estimates would have agreed with refs . . we therefore conclude that the estimates of @xmath16 presented here are more accurate than results reported to date .
Abstract:

Answer: we study thermodynamics of the 3d hubbard model at half filling on approach to the nel transition by means of large - scale unbiased diagrammatic determinant monte carlo simulations . we obtain the transition temperature in the strongly correlated regime , as well as temperature dependence of energy , entropy , double occupancy , and the nearest - neighbor spin correlation function . our results improve the accuracy of previous unbiased studies and present accurate benchmarks in the ongoing effort to realize the antiferromagnetic state of matter with ultracold atoms in optical lattices .


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 recent years , there have been much interest in dirac materials such as graphene@xcite , topological insulators@xcite , high - temperature @xmath3-wave superconductors@xcite , which are described by the dirac - like equation at low energies . dirac materials show novel physical features because of the chiral band dispersion with spin - momentum locking and the existence of the weyl points at which conduction and valence bands touch , distinguishing themselves from other conventional materials with parabolic energy dispersion at low energies . recently three dimensional ( 3d ) weyl semimetals , which are 3d analogs of graphene , have attracted considerable attention because they are expected to give a broad spectrum of unusual physical properties as two dimensional graphene . in recent experiments of angle - resolved photoemission spectroscopy and scanning tunneling microscopy , several materials such as cd@xmath4as@xmath5@xcite , na@xmath4bi@xcite , nbas@xcite , tap@xcite , zrte@xmath6@xcite and taas@xcite have been identified as weyl semimetals . in addition , several efforts have been made on the realization of weyl semimetals in artificial systems such as photonic crystals@xcite . recently , a new type of weyl semimetals referred to as multi - weyl semimetals has been proposed@xcite . these materials are characterized by double ( triple ) weyl - nodes in which the band dispersion is linear along one direction and quadratic ( cubic ) along the remaining two directions . because of the change in topological nature and the enhancement of the density of states , the anisotropic non - linear energy dispersion and a modified spin - momentum locking structure have important consequences in topological , transport and interaction - induced properties@xcite . despite their broad interest , however , there have been lack of studies on the collective modes of charge oscillations in chiral gapless systems including multi - weyl systems . the goal of this paper is to investigate electronic collective modes of chiral gapless electron - hole systems in 3d and to find wave - vector - dependent plasmon dispersions , considering interband transitions along with chirality of the systems . given the great current interest in gapless linearly dispersing dirac systems and the emergence of multi - weyl semimetals , it is necessary to obtain results for linear and higher - order dispersions in 3d in order to develop intuition about the dispersion and chirality dependence of electronic properties . we consider both gapless semimetals with an isotropic band dispersion and multi - weyl semimetals with an anisotropic band dispersion . we calculate the plasmon modes and the energy loss functions ( spectral strength ) in gapless semimetals , and present analytical and numerical results for long - wavelength plasma frequencies focusing on the effect of chirality . throughout the paper we consider the systems with the finite carrier density , i.e. , either electron or hole doped ( or gated ) systems at zero temperatures . we find many intriguing and unexpected features of the plasmon modes and their landau damping in gapless semimetals . for chiral systems with the isotropic linear and quadratic band dispersion ( @xmath7 in @xmath8 ) the long - wavelength plasmons lie outside the landau damping regions due to the depolarization effects arising from the interband transition and therefore they do not decay by producing electron - hole pairs . for the systems with higher - order dispersions ( @xmath9 ) , the long - wavelength plasma frequencies divided by fermi energy ( i.e. , @xmath10 ) increase as the density decreases and enter the interband single particle excitation ( spe ) region at the critical carrier density , where the dispersion shows a discrete energy jump . for the multi - weyl systems with an anisotropic dispersion ( i.e. , with a linear band dispersion along @xmath11 direction and non - linear in @xmath12-@xmath13 directions ) , one interesting salient feature in the calculated plasma frequency is the density dependence of the plasmon . we find that the density dependence of the long - wavelength plasma frequency along the direction of non - linear dispersion behaves like that of the isotropic linear band model ( @xmath14 ) , while along the direction of linear dispersion it behaves like that of the isotropic non - linear model ( @xmath15 ) . we find that both the plasmons remain undamped over a broad range of density and interaction strength due to the chirality induced red - shift of plasmon modes . we believe that our predictions may be easily observable via inelastic light - scattering spectroscopy@xcite or inelastic electron - scattering spectroscopy @xcite .    the rest of the paper is organized as follows . in sec . ii , we consider semimetals with arbitrary isotropic band dispersion and present analytical and numerical results for long - wavelength plasma frequencies at zero temperature . in sec . iii , we consider multi - weyl semimetals with an anisotropic dispersion and present results for zero - temperature plasma frequencies . in sec . iv , we summarize and conclude our results . before we calculate the collective modes in gapless multi - weyl semimetals with anisotropic energy band dispersions , we first consider the gapless semimetals with the isotropic band dispersions . for the isotropic model we consider the following hamiltonian that describes a chiral gapless system with the symmetric energy dispersion , @xmath16 where @xmath17 is the wave vector , @xmath18 , @xmath19 are pauli matrices acting in the space of two bands involved at the weyl point , and @xmath20 and @xmath21 are material dependent parameters , which have units of energy and momentum , respectively . the corresponding energy dispersions are given by @xmath22 and the eigenfunctions corresponding to the @xmath23 energies are @xmath24 where @xmath25 and @xmath26 . the @xmath27 state with positive energy represents the conduction band , and the @xmath28 state with negative energy represents the valence band . note that the hamiltonian in eq . ( [ isotropic_model ] ) with @xmath14 corresponds to the hamiltonian of weyl semimetals . to understand the effects of chirality on the plasmon properties beyond the energy dispersion of the system , we also consider the following form of a non - chiral gapless hamiltonian : @xmath29 the non - chiral hamiltonian we introduce in eq . ( [ non_chiral_isotropic_model ] ) has exactly the same energy dispersion as the chiral system in eq . ( [ isotropic_model ] ) . unlike the chiral model , however , its eigenstates are completely independent of each other so that the chiral nature of wave functions is absent . it will be shown later that the interband transitions associated with chirality are largely responsible for differences between plasmons in the presence and absence of chirality . plasmons are defined as longitudinal in - phase oscillation of all the carriers driven by the self - consistent electric field generated by the local variation in charge density . to find the full plasmon dispersion at finite wave vectors we need the quantum mechanical many - body theory for the collective motion of all carriers@xcite . within the random phase approximation ( rpa ) , the plasmon dispersion is obtained by finding zeros of the dynamical dielectric function , which is expressed as @xcite @xmath30 where @xmath31 is the coulomb interaction and @xmath32 is a background dielectric constant . the non - interacting polarizability @xmath33 is given by @xmath34 where @xmath35^{-1}$ ] is the fermi distribution function for the band @xmath36 , @xmath37 is the chemical potential , @xmath38 , @xmath39 is the overlap factor , @xmath40 is the angle between @xmath41 and @xmath42 , and @xmath43 denotes the number of degenerate weyl nodes of the system . since we focus on the long - wavelength plasmons in this work , we neglect transitions between different nodes located at different momenta . within this assumption , contributions from other weyl nodes can be taken into account by multiplying the degeneracy factor @xmath43 . it is important to notice that in the non - chiral model described by eq . ( [ non_chiral_isotropic_model ] ) and later by eq . ( [ non_chiral_anisotropic_model ] ) , the overlap factor becomes the kronecker @xmath44 , i.e. , @xmath45 , which excludes the possibility of any interband transitions between the conduction and valence bands . in this sense , the non - chiral model considered here is a two - band model without interband transitions , which is effectively a single - band model with the same energy dispersion . $ ] in ( @xmath46,@xmath47 ) space . ( a ) and ( b ) show the loss functions in the absence of chirality for ( a ) @xmath48 ( parabolic dispersion ) and ( b ) @xmath49 ( cubic dispersion ) , whereas ( c ) and ( d ) show the energy - loss functions in the presence of chirality for ( c ) @xmath48 and ( d ) @xmath49 . the dashed lines represent the boundaries of the intraband and interband electron - hole continua . note that in the presence of chirality , the plasmon energies are shifted downward for all wave vectors compared with those in the absence of chirality . here we use the coupling constant @xmath50 , and the fermi energy @xmath51 for @xmath48 and @xmath52 for @xmath49 , respectively . ] we first consider the loss function of 3d gapless semimetals , which is calculated from the dielectric function ( i.e. , @xmath53 $ ] ) and can be directly measured in experiments such as inelastic electron spectroscopy . the loss function is related to the dynamical structure factor @xmath54 by @xmath55 $ ] , which gives a direct measure of the spectral strength of the various elementary excitations . figure [ fig:1 ] shows the density plots of calculated energy - loss functions in ( @xmath56 ) space for @xmath48 ( parabolic dispersion ) and @xmath49 ( cubic dispersion ) in the absence ( top panels ) and presence ( bottom panels ) of chirality . plasmon dispersions are given by sharp peaks of the energy - loss function , which correspond to the poles of the dielectric function . when both re[@xmath57 and im[@xmath57 become zero ( i.e. , @xmath58 , which defines the plasmon mode ) , the imaginary part of the inverse dielectric function becomes the dirac @xmath44-function , i.e. , @xmath53=w(q ) \delta(\omega-\omega_{\rm p}(q))$ ] with the strength @xmath59/\partial \omega|_{\omega=\omega_{\rm p}(q ) } \right \}^{-1},\ ] ] where @xmath60 is the plasma frequency at a given wave vector @xmath46 . thus , an undamped plasmon shows up as a well - defined @xmath44-function peak in the loss function as indicated by sharp yellow solid lines in fig . [ fig:1 ] . the undamped plasmon mode in general carries most of spectral weights and should be observable in experiments ( i.e. , it is expected that the mode does not decay by electron - hole pairs ) . the dotted lines in fig . [ fig:1 ] represent the boundaries of the intraband and interband electron - hole single particle excitation ( spe ) continua . the electron - hole spe continua show up as weak broad incoherent structure and carries small spectral weight . when the plasmon mode enters the spe continuum at the critical wave vector @xmath61 , the dielectric function @xmath62 has a finite imaginary part and the plasmon mode becomes damped via landau damping . the plasmon mode inside the landau damping region decays by emitting intraband or interband electron - hole pairs , which is now allowed by energy - momentum conservation . the broadened peaks inside the spe regions in fig . [ fig:1 ] indicate the damped plasmons . the plasmon energy scaled by the fermi energy ( @xmath10 ) is strongly dependent on the band structure . one interesting result is the chirality dependence of the plasmon energy , and it is important to note that in the presence of chirality , plasma frequencies are red - shifted . this is due to the depolarization effect on the plasmon modes , arising from interband transitions .     and ( b ) @xmath49 . red ( blue ) lines correspond to the results in the presence ( absence ) of chirality . note that for @xmath48 , @xmath63 diverges in the @xmath64 limit in the absence of chirality , whereas it converges to a finite value in the presence of chirality . for @xmath49 , the chiral plasmon dispersion shows a discrete jump at @xmath65 , where plasmons start decaying into electron - hole pair excitations . here we use the coupling constant @xmath50 . ] figure [ fig:2 ] shows the long - wavelength plasma frequencies for @xmath48 and @xmath49 as a function of fermi energy in the presence ( red ) and absence ( blue ) of chirality . in this figure the boundary of interband spe at @xmath66 is @xmath67 . if the plasma frequency is larger ( smaller ) than @xmath67 , the plasmon decays ( does not decay ) by landau damping . for both @xmath48 and @xmath49 , @xmath10 in the absence of chirality increases monotonically with decreasing density and diverges in the low density limit ( @xmath64 ) . in the presence of chirality , the plasma frequencies are always smaller than those without chirality , and the energy difference between the two plasmons grows as the carrier density decreases . for @xmath48 , as the density goes to zero , the plasma frequency approaches a finite value less than @xmath67 , i.e. , @xmath68 , indicating that the long - wavelength plasmon does not enter the interband electron - hole continuum for the whole range of densities and thus the plasmon is not damped through the landau damping . for @xmath49 , the plasma frequency as a function of fermi energy exhibits a discrete jump at a critical value of @xmath69 , and below the critical value the plasmon energy becomes bigger than @xmath67 and enters into the interband electron - hole continuum decaying via landau damping . note that interband transitions are more suppressed at high densities due to pauli - blocking , which narrows the energy range over which interband transitions can occur . thus , the effect of interband transition on plasmons becomes weaker at high densities , resulting in small depolarization shift in plasmon energy . for this reason , the signature of chirality in the plasmon dispersions is significant at low densities .    to get further insight into the results shown in fig . [ fig:1 ] , here we present the calculated analytic expressions for the leading - order long - wavelength plasma frequencies ( for details , see appendix a ) , @xmath70 in the absence of chirality and @xmath71^{-\frac{1}{2}}\ ] ] in the presence of chirality , where @xmath72 is the coupling constant characterizing the interaction strength . the above results are calculated in @xmath73 limit for @xmath74 . note that for @xmath14 weyl semimetals a linear band dispersion leads to cut - off dependent long - wavelength plasma frequencies@xcite . it is easy to see from eq . ( [ chiral plasma frequency ] ) that the plasma frequency with chirality is red - shifted with respect to that without chirality , showing different density dependence at low densities . this difference originates from the interband transition contribution , which appears in the additional term of the polarizability with the opposite sign of the intraband contribution ( see eq . ( a8 ) in appendix a ) . this indicates that interband transitions and associated chiral nature of wavefunctions contribute to the depolarization of the screening and are responsible for the red - shift of the plasma frequencies . it is interesting that the density dependence of the long wavelength plasma frequencies shown in eq . ( [ nonchiral plasma frequency ] ) can be obtained from the classical plasma frequency @xcite @xmath75 where @xmath76 is the charge carrier density and @xmath77 is the effective mass of the charge carrier . by using the momentum relation @xmath78 ( @xmath79 and @xmath80 are the fermi velocity and the fermi wave - vector , respectively ) and the energy dispersion relation @xmath22 for arbitrary band dispersion , the density dependence of the classical plasma frequency can be calculated as @xmath81 which agrees with the full rpa result in eq . ( [ nonchiral plasma frequency ] ) for the non - chiral case . thus , for @xmath14 we have @xmath82 , i.e. , the plasma frequency scaled by the fermi energy is independent of @xmath83@xcite . for @xmath84 , @xmath85 with @xmath86 , i.e. , the scaled plasma frequency increases as the density ( or fermi energy ) decreases and diverge as @xmath87 . note that the interband transition red - shifts the classical plasma frequencies in the presence of chirality . in the previous section we have discussed the plasmon properties of isotropic gapless semimetals . in this section we explore , within the rpa , the plasmon properties of the anisotropic multi - weyl system whose dispersion is non - linear in the in - plane directions , but linear in the out - of - plane direction . we consider the following hamiltonian that describes a multi - weyl node of order @xmath1 , @xmath88 , h_j^{\rm ch}(\bm{k})= e_0\left[\left({k_{-}\over k_0}\right)^j\sigma_{+ } + \left({k_{+}\over k_0}\right)^j\sigma_{- } + { k_z\over k_0}\sigma_z\right],\ ] ] where @xmath89 , @xmath90 , and @xmath21 and @xmath20 are material dependent parameters . throughout this section the momentum and the energy are normalized by @xmath21 and @xmath20 , respectively . with these normalized quantities , we can write the energy dispersion as @xmath91 , where @xmath92 , @xmath93 is the in - plane momentum corresponding to non - linear dispersion ( @xmath94 ) in the @xmath95 plane and @xmath96 is the out - of - plane momentum along the @xmath11 direction where the dispersion is linear ( @xmath97 ) . the overlap factor @xmath98 is given by @xmath99\ ] ] where @xmath100 and @xmath101 are related to @xmath102 and @xmath103 through the coordinate transformation , respectively ( see eq . ( b1 ) in appendix b ) . for comparison , we also introduce an anisotropic non - chiral model with the same energy dispersion : @xmath104     in the absence [ top panels , ( a ) and ( b ) ] and presence [ bottom panels , ( c ) and ( d ) ] of chirality along the @xmath11-direction [ left panels , ( a ) and ( c ) ] and the in - plane direction [ right panels , ( b ) and ( d ) ] . note that plasmon dispersions along the @xmath11-direction behave similarly as those for the @xmath14 isotropic model and plasmon dispersions along the in - plane direction behave as those for the @xmath48 isotropic model . plasmon frequencies along both directions are red - shifted in the presence of chirality , as in the isotropic model . here we use the fermi energy @xmath51 and the coupling constant @xmath50 . ]    to investigate the plasmon dispersions and their damping we calculate the energy loss function of the system . figure [ fig:3 ] shows the density plots of the calculated energy loss functions for @xmath105 in the absence ( top panels ) and presence ( bottom panels ) of chirality . figure [ fig:3 ] ( a ) and ( c ) show the density plots of the energy loss function in ( @xmath106 ) space for a fixed in - plane wave vector ( @xmath107 ) , and fig . [ fig:3 ] ( b ) and ( d ) show the density plots of the energy loss function in ( @xmath108 ) space for a fixed out of plane wave vector ( @xmath109 ) . the yellow lines with the largest spectral weight in the energy loss correspond to the plasmon dispersions of the system . the two most important salient features of the results shown in fig . [ fig:3 ] are following . first , the plasmon dispersions for both in - plane and out - of - plane directions are red - shifted in the presence of chirality , which is consistent with the isotropic results . second , the plasmon dispersion along the out - of - plane direction ( where the band dispersion is linear ) shows a similar behavior as that in the @xmath14 isotropic system . this is because the longitudinal plasmon oscillations propagating along the out - of - plane direction arise from the collective carriers with linear band dispersion and thus their effective motion is essentially identical to those in gapless semimetals with linear dispersion ( @xmath14 ) . the same argument is also applied to the in - plane plasmon modes which propagate along the in - plane direction , where the carriers have the parabolic dispersion and therefore the plasmon dispersion behaves as that in the @xmath48 isotropic system . [ left panels , ( a ) and ( c ) ] and @xmath110 [ right panels , ( b ) and ( d ) ] . along the in - plane momentum direction [ top panels , ( a ) and ( b ) ] , @xmath10 for both @xmath105 and @xmath110 behaves as that in weyl semimetals . along the @xmath11 direction [ bottom panels , ( c ) and ( d ) ] , for both @xmath105 and @xmath110 @xmath10 converges to a finite value in the @xmath111 limit . here we use the coupling constant @xmath50 . ]    in contrast to the plasmon dispersions , however , the density dependence of the long - wavelength plasmons shows a non - trivial relation between the isotropic and anisotropic models . figure [ fig:4 ] presents calculated long - wavelength plasma frequencies as a function of fermi energy for @xmath105 and @xmath110 in the presence ( red ) and absence ( blue ) of chirality . the in - plane plasma frequencies @xmath10 [ ( a ) for @xmath105 and ( b ) for @xmath110 ] , show the same behavior as the @xmath14 isotropic system , i.e. , the normalized frequencies are weakly dependent on the carrier density over a wide range both with and without chirality . note that this behavior is similar to that in gapless semimetals with linear dispersion . for the @xmath11-direction [ ( c ) for @xmath105 and ( d ) for @xmath110 ] , @xmath10 increases with decreasing fermi energy in the absence of chirality , and diverges in the @xmath64 limit for both @xmath105 and @xmath110 . in the presence of chirality , for @xmath105 and @xmath110 it converges to a finite value in the @xmath64 limit . note that these results are qualitatively similar to the chiral signature of @xmath48 and @xmath49 isotropic systems , respectively , as discussed in sec . ii , though here the energy dispersion along the @xmath11-direction is linear . it is important to notice that this result is against our expectations because the density dependence of the anisotropic plasmons along a direction with a specific energy dispersion does not match with the corresponding plasmons in the isotropic model with the same energy dispersion . it should also be noted that for the chiral case anisotropic plasmons along both the in - plane and out - of - plane directions are significantly red - shifted so that they are outside the electron - hole continua for a broad range of density and coupling strength . for @xmath110 plasmons , this is in sharp contrast to the counterpart of the isotropic case ( @xmath49 ) , where plasmons exhibit a discrete jump merging into the electron - hole continua at a critical fermi energy . here we analyze the previous numerical results in fig . [ fig:4 ] with analytical expressions . similarly as the isotropic case , we can also derive long - wavelength plasma frequencies from the linear response theory in the many - body approach . in the absence of chirality , we obtain plasma frequencies for the anisotropic model to be @xmath113 along the @xmath114 momentum direction and @xmath115 along the @xmath96 momentum direction ( for details , see appendix b ) . the difference in density dependence between the two plasma frequencies can also be understood by the classical derivation of plasma frequency , as in the isotropic case . since the density of states of the anisotropic model in eq . ( [ anisotropic_model ] ) is calculated to be @xmath116 , the total carrier density @xmath76 becomes @xmath117 . since the energy dispersion is anisotropic , it is expected that the effective mass is also anisotropic . thus , electrons would effectively behave as if they have the energy dispersion @xmath94 with @xmath118 along the in - plane @xmath114 momentum direction , and @xmath97 with @xmath119 for the out - of - plane @xmath96 momentum direction . putting @xmath76 and @xmath77 into eq . ( [ classical plasmon ] ) , we can obtain the fermi energy dependence of the plasma frequency in classical limit : @xmath120 along the @xmath114 momentum direction and @xmath121 along the @xmath96 momentum direction . the density dependence of plasmons obtained from both classical and quantum mechanical approaches are consistent with the numerical results in fig . [ fig:4 ] . note that as in the case of the isotropic system , the density dependence of classical plasma frequencies is in good agreement with that of the non - chiral plasmon . in this paper we investigate theoretically electronic collective modes of 3d chiral gapless electron - hole systems and find the wave vector dependent plasmon dispersion . we have calculated long - wavelength plasma frequencies and their density dependence both in semimetals with an isotropic band dispersion and in multi - weyl semimetals with an anisotropic band dispersion . we find that the interband transition associated with chirality leads to the depolarization shift of plasma frequencies irrespective of band dispersion .    for the isotropic parabolic dispersion ( @xmath48 ) , the depolarization shift of the long - wavelength plasmons arises from the interband electron - hole transition and the plasmons lie outside the interband electron - hole continuum ( i.e. , @xmath122 ) . thus , the plasmons do not decay via landau damping . for the cubic dispersion ( @xmath49 ) , we find that the plasma frequency ( @xmath63 ) increases as the density decreases and enters the interband single particle excitation region at a critical carrier density showing a discrete energy jump .    for the anisotropic dispersion , we have calculated the plasma frequencies in a system with a linear band dispersion along one specific direction and non - linear dispersion in the other remaining directions . we find that the density dependence of the long - wavelength plasma frequency along the direction of non - linear dispersion shows a similar behavior as that of the linear band isotropic model ( @xmath14 ) , in which @xmath10 depends weakly on the density over a wide range of fermi energies . on the other hand , the density dependence of the long - wavelength plasma frequency along the direction of linear dispersion shows a similar behavior as that of the isotropic model with the non - linear dispersion . the long - wavelength plasmons along all the directions are undamped due to the chirality induced red - shift of plasma frequencies . our predicted plasmon properties clearly distinguish gapless semimetals from the extensively studied usual parabolic 3d electron systems . we believe that our predictions can be tested in doped gapless semimetals using inelastic light scattering @xcite and electron scattering @xcite spectroscopies . this research was supported by basic science research program through the national research foundation of korea ( nrf ) funded by the ministry of education under grant no . 2015r1d1a1a01058071 ( hm ) and 2014r1a2a2a01006776 ( ehh ) . in this section , we show the derivation of the analytical expressions for the zero - temperature polarization function and plasma frequencies in the @xmath123 limit . the polarization function for the isotropic model can be written as @xmath124 where @xmath125    note that at zero temperature the fermi distribution function is given by @xmath126 , where @xmath127 is the step function [ @xmath128 for @xmath129 and @xmath130 for @xmath131 . we can break the polarization function into two parts due to intraband transitions and interband transitions , respectively : @xmath132 , where @xmath133 and @xmath134 . for convenience , we use the change of variable @xmath135 and introduce dimensionless quantities @xmath136 and @xmath137 , where @xmath138 is the density of states at the fermi energy given by @xmath139 . then we can rewrite @xmath140 and @xmath141 as    @xmath142    and @xmath143    where @xmath144 . in the @xmath123 limit @xmath145 $ ] and @xmath146 . in the presence of chirality the overlap factor is given by @xmath147 , where @xmath148 whereas in the absence of chirality @xmath149 and @xmath150 . by putting all these equations into eq . ( [ eq : intra ] ) and expanding up to the second order in @xmath13 , we can obtain the long - wavelength polarization function . for the intraband part @xmath151 , the polarization functions for the chiral and non - chiral cases are the same up to the second order , given by @xmath152 for the interband part @xmath153 , note that there is no interband contribution for the non - chiral case due to the absence of interband transitions ( @xmath154 ) . for the chiral case , @xmath155 and when we expand the integrand only up to the second order , we can approximate @xmath156 , yielding    @xmath157    where @xmath158 is the hypergeometric function , which can be obtained after substituting @xmath159 . thus the full polarizability in the long wavelength limit ( @xmath123 ) is given by @xmath160 for the non - chiral case whereas @xmath161 for the chiral case . note that the polarization function diverges when @xmath14 because a linear band dispersion leads to the divergence of the polarization function at infinite cut - off@xcite . by putting eq . ( [ nonchiral polarization function ] ) and eq . ( [ chiral polarization function ] ) into the dielectric function in eq . ( 5 ) in the main text , we arrive at the analytic expressions for the leading - order long - wavelength plasma frequencies given in eq . ( 8) and eq . ( 9 ) in the main text . in this section , we show the derivation of the polarization function for the anisotropic case . here we consider only the non - chiral case , where the overlap factor is given by @xmath162 and @xmath154 . note that due to the anisotropic band structure it is no longer convenient to normalize the polarization function with the density of states as in the isotropic case . thus for the anisotropic case we normalize the polarization function by @xmath163 , where @xmath164 , @xmath165 and @xmath166 . in the following , for brevity , the momentum and energy without tilde are considered to be normalized by @xmath21 and @xmath20 , as in the main text . we consider the following coordinate transformation : @xmath167 note that with this coordinate transformation the energy dispersion becomes linear : @xmath168 . the jacobian corresponding to this transformation is given by @xmath169 . for the non - chiral case , there is no interband contribution , thus we consider only the intraband contribution . with a similar procedure as in the previous section , we can write the polarization function for the in - plane momentum @xmath170 along the @xmath13-axis and out - of - plane momentum @xmath171 along the @xmath11-axis as              here we use the relation @xmath181 where @xmath182 is the beta function . note that linear terms in @xmath171 vanish both in eq . ( [ eq : q_parallel ] ) and eq . ( [ eq : q_z ] ) due to the symmetry of the angular integrals . combining eq . ( [ eq : inter ] ) and the dielectric function in eq . ( 5 ) in the main text , we arrive at the analytic expressions for the leading - order long - wavelength plasma frequencies given in eq . ( 15 ) and eq . ( 16 ) in the main text .
Abstract:

Answer: we investigate collective modes in three dimensional ( 3d ) gapless multi - weyl semimetals with anisotropic energy band dispersions ( i.e. , @xmath0 with a positive integer @xmath1 ) . for comparison , we also consider the gapless semimetals with the isotropic band dispersions ( i.e. @xmath2 ) . we calculate analytically long - wavelength plasma frequencies incorporating interband transitions and chiral properties of carriers . for both the isotropic and anisotropic cases , we find that interband transitions and chirality lead to the depolarization shift of plasma frequencies . for the isotropic parabolic band dispersion the long - wavelength plasmons do not decay via landau damping , while for the higher - order band dispersions the long - wavelength plasmons experience damping below a critical density . for systems with the anisotropic dispersion the density dependence of the long - wavelength plasma frequency along the direction of non - linear dispersion behaves like that of the isotropic linear band model , while along the direction of linear dispersion it behaves like that of the isotropic non - linear model . plasmons along both directions remain undamped over a broad range of densities due to the chirality induced depolarization shift . our results provide a comprehensive picture of how band dispersion and chirality affect plasmon behaviors in 3d gapless chiral systems with the arbitrary band dispersion .


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 investigation of the properties of mesons composed of a heavy quark and antiquark ( @xmath5 , @xmath6 , @xmath7 ) gives very important insight into heavy quark dynamics . heavy quarkonia have a rich spectroscopy with many narrow states lying under the threshold of open flavour production . excited states experience different decays among which there are radiative transitions to lower levels . the theoretical analysis shows that many properties of heavy quarkonia , including mass spectra and radiative decay rates , are significantly influenced by relativistic effects . thus their inclusion is necessary for the correct description of the spectroscopy and the determination of quarkonium wave functions . radiative decays are the most sensitive to relativistic effects . some of these decays , which are forbidden in the exact nonrelativistic limit ( so - called hindered transitions ) due to the orthogonality of initial and final meson wave functions , have decay rates of the same order as the allowed ones . in the relativistic description of mesons an important role is played by properties of the confining quark - antiquark interaction , in particular its lorentz structure . thus comparison of theoretical predictions with experimental data can provide valuable information on the form of the confining potential . such information is of great practical interest , since at present it is not possible to obtain the @xmath8 potential in the whole range of distances from the basic principles of quantum chromodynamics ( qcd ) . as it is well known , the growing of the strong coupling constant with distance makes perturbation theory inapplicable at large distances ( in the infrared region ) . in this region it is necessary to account for nonperturbative effects connected with the complicated structure of the qcd vacuum . all this leads to a theoretical uncertainty in the @xmath8 potential at large and intermediate distances . it is just in this region of large and intermediate distances that most of the basic meson characteristics are formed .    at present , a vast set of experimental data is available on the masses and different decays of heavy quarkonia . however , not all states predicted by theory have been observed yet , while the others need confirmation . such missing or unconfirmed states are present both in charmonium ( @xmath9 , @xmath10 , @xmath11 , @xmath12 ) and bottomonium ( spin - singlet @xmath13 and @xmath14 states , @xmath4 states ) . the different possibilities for their experimental observation are proposed and widely discussed in the literature @xcite . radiative transitions from the spin - triplet levels with @xmath15 to these states as well as their subsequent radiative decays play an important role in these proposals . the missing charmonium states can also be searched in @xmath16 meson decays and identified by their radiative transitions @xcite . for this purpose , reliable relativistic predictions for the masses of these states and for the rates of radiative transitions involving them are necessary . the properties of the @xmath0 meson are of special interest , since it is the only heavy meson consisting of two heavy quarks with different flavour . this difference of quark flavours forbids annihilation into gluons . as a result , the excited @xmath0 meson states lying below the @xmath17 production threshold undergo pionic or radiative transitions to the ground pseudoscalar state which then decays weakly . there should be a rather rich set of such narrow states which are considerably more stable than corresponding charmonium or bottomonium states . the cdf collaboration @xcite reported the discovery of the @xmath0 ground state in @xmath18 collisions at fermilab . more experimental data are expected to come in the near future from new hadronic colliders . the purpose of this paper is to give a detailed analysis of mass spectra and radiative transitions in charmonium , bottomonium and @xmath0 mesons with the comprehensive account of the relativistic effects . this will allow one to get valuable information about the lorentz structure of confining quark interactions from the comparison of obtained predictions with available experimental data . on the other hand , it will indicate the processes in which the missing states can be searched for . the paper is organized as follows . in sec . [ rqm ] we describe our relativistic quark model . the expression for the heavy quark - antiquark quasipotential with the account of relativistic ( including retardation effects ) and one loop radiative corrections is given in sec . [ sec : bc ] . there it is applied to the calculation of the charmonium , bottomonium and @xmath0 meson mass spectra . in sec . [ sec : pvdc ] pseudoscalar and vector decay constants of the @xmath0 meson are calculated with the account of relativistic corrections and compared with other theoretical predictions . in sec . [ rad ] the relativistic expressions for the radiative transition matrix elements in the quasipotential approach are given . they are used for the calculation of the decay rates of radiative m1 and e1 transitions in secs . [ m1 ] and [ e1 ] , respectively . the role of relativistic effects in these transitions is investigated . special attention is payed to the influence of the lorentz structure of the quark potential on the relativistic corrections to decay rates . pure vector and scalar potentials as well as their mixture are considered . the obtained results are compared with available experimental data , and the possibilities for searching the missing states in bottomonium are discussed . finally , our conclusions are given in sec . [ sec : conc ] . in the quasipotential approach a meson is described by the wave function of the bound quark - antiquark state , which satisfies the quasipotential equation @xcite of the schrdinger type @xcite @xmath19 where the relativistic reduced mass is @xmath20 and @xmath21 , @xmath22 are given by @xmath23 here @xmath24 is the meson mass , @xmath25 are the quark masses , and @xmath26 is their relative momentum . in the center of mass system the relative momentum squared on mass shell reads @xmath27[m^2-(m_1-m_2)^2]}{4m^2}.\ ] ]    the kernel @xmath28 in eq . ( [ quas ] ) is the quasipotential operator of the quark - antiquark interaction . it is constructed with the help of the off - mass - shell scattering amplitude , projected onto the positive energy states . constructing the quasipotential of the quark - antiquark interaction , we have assumed that the effective interaction is the sum of the usual one - gluon exchange term with the mixture of long - range vector and scalar linear confining potentials , where the vector confining potential contains the pauli interaction . the quasipotential is then defined by @xcite @xmath29 with @xmath30 where @xmath31 is the qcd coupling constant , @xmath32 is the gluon propagator in the coulomb gauge @xmath33 and @xmath34 ; @xmath35 and @xmath36 are the dirac matrices and spinors @xmath37 with @xmath38 . the effective long - range vector vertex is given by @xmath39 where @xmath40 is the pauli interaction constant characterizing the anomalous chromomagnetic moment of quarks . vector and scalar confining potentials in the nonrelativistic limit reduce to @xmath41 reproducing @xmath42 where @xmath43 is the mixing coefficient . the expression for the quasipotential of the heavy quarkonia , expanded in @xmath1 without and with retardation corrections to the confining potential , can be found in refs . @xcite and @xcite , respectively . the structure of the spin - dependent interaction is in agreement with the parameterization of eichten and feinberg @xcite . the quasipotential for the heavy quark interaction with light antiquark without employing the expansion in inverse powers of the light quark mass is given in ref . all the parameters of our model such as quark masses , parameters of the linear confining potential @xmath44 and @xmath16 , mixing coefficient @xmath43 and anomalous chromomagnetic quark moment @xmath40 are fixed from the analysis of heavy quarkonium masses @xcite ( see sec .  [ sec : bc ] ) and radiative decays @xcite ( see secs . [ rad][e1 ] ) . the quark masses @xmath45 gev , @xmath46 gev and the parameters of the linear potential @xmath47 gev@xmath48 and @xmath49 gev have usual values of quark models . the value of the mixing coefficient of vector and scalar confining potentials @xmath50 has been determined from the consideration of the heavy quark expansion for the semileptonic @xmath51 decays @xcite and charmonium radiative decays @xcite . finally , the universal pauli interaction constant @xmath52 has been fixed from the analysis of the fine splitting of heavy quarkonia @xmath53- states @xcite . note that the long - range magnetic contribution to the potential in our model is proportional to @xmath54 and thus vanishes for the chosen value of @xmath52 . in the present paper we will take into consideration the retardation corrections as well as one - loop radiative corrections . the heavy quark - antiquark potential with the account of retardation effects and one loop radiative corrections can be presented in the form of a sum of spin - independent and spin - dependent parts . for the spin - independent part we have @xcite @xmath56\cr\nonumber\\ & & + \frac{1}{2m_1m_2}\biggl(\left\{-\frac43\frac{\bar\alpha_v}{r } \left[{\bf p}^2+\frac{({\bf p\cdot r})^2}{r^2}\right]\right\}_w\cr\nonumber\\ & & -\frac43\frac{\beta_0\alpha_s^2(\mu^2)}{2\pi}\left\{{\bf p}^2    \frac{\ln(\mu r)}{r } + \frac{({\bf p\cdot   r})^2}{r^2}\left(\frac{\ln(\mu r)}{r}-\frac1r\right)\right\}_w\biggr ) \cr\nonumber\\ & & + \left[\frac{1-\varepsilon}{2m_1m_2}-\frac{\varepsilon}{4 } \left(\frac{1}{m_1 ^ 2}+\frac{1}{m_2 ^ 2}\right)\right]\left\{ar\left[{\bf p}^2 -\frac{({\bf p\cdot r})^2}{r^2}\right]\right\}_w\cr\nonumber\\ & & + \left[\frac{1}{4}\left(\frac{1}{m_1 ^ 2}+\frac{1}{m_2 ^ 2}\right)+ \frac{1}{m_1m_2}\right]b{\bf p}^2,\end{aligned}\ ] ] where @xmath57,\\ a_1&=&\frac{31}{3}-\frac{10}{9}n_f,\cr \beta_0&=&11-\frac23n_f.\nonumber\end{aligned}\ ] ] here @xmath58 is the number of flavours and @xmath59 is a renormalization scale . note that for a quantity quadratic in the momenta we use the weyl prescription @xcite : @xmath60    for the dependence of the qcd coupling constant @xmath61 on the renormalization point @xmath62 we use the leading order result @xmath63 in our numerical calculations we set the renormalization scale @xmath64 and @xmath65  gev , which gives @xmath66 for @xmath67 ( charmonium ) ; @xmath68 for @xmath69 ( bottomonium ) ; and @xmath70 for @xmath71 , @xmath72 ( @xmath0 meson ) . the spin - dependent part of the quark - antiquark potential for unequal quark masses with the inclusion of radiative corrections @xcite can be presented in our model as follows : @xmath73 + c\ { \bf s}_1\cdot { \bf s}_2 + d\ { \bf l}\cdot({\bf s}_1-{\bf s}_2),\ ] ] @xmath74\biggr ) -\frac{a}{r}\biggr\}+\frac1{m_1m_2}\frac43\frac{\alpha_s(\mu^2)}{r^3 } \biggl(1 +   \frac{\alpha_s(\mu^2)}{\pi}\biggl[\frac16-\frac{\beta_0}{12}\cr\nonumber\\ & & + \gamma_e\left(\frac{\beta_0}{2}-\frac32\right)+\frac{\beta_0}{2 } \ln(\mu r ) -\frac32\ln(\sqrt{m_1m_2}\,r)\biggr]\biggr)\cr\nonumber\\ & & + \left(\frac{1}{m_1 ^ 2}-\frac{1}{m_2 ^ 2}\right )    \frac{\alpha_s^2(\mu^2)}{2\pi r^3}\ln\frac{m_2}{m_1 }    + \frac12\left(\frac{1}{m_1}+\frac{1}{m_2}\right)^2     ( 1+\kappa)(1-\varepsilon)\frac{a}{r},\\ \cr \label{b } b&= & \frac{1}{3m_1 m_2}\biggl\{\frac{4\alpha_s(\mu^2)}{r^3}\biggl(1 + \frac{\alpha_s(\mu^2)}{\pi}\biggl[\frac{29}{6}-\frac{1}{4}\beta_0 + \gamma_e\left(\frac{\beta_0}{2}-3\right)+\frac{\beta_0}{2}\ln(\mu r ) , \cr\nonumber\\ & & -3\ln(\sqrt{m_1m_2}\,r)\biggr]\biggr ) + ( 1+\kappa)^2(1-\varepsilon)\frac{a}{r}\biggr\},\\ \cr \label{c } c&= & \frac{4}{3m_1 m_2}\biggl\{\frac{8\pi\alpha_s(\mu^2)}{3}\biggl(\left[1 + \frac{\alpha_s(\mu^2)}{\pi}\left(\frac{5}{12}\beta_0-\frac{11}{3 }   -\left[\frac{m_1-m_2}{m_1+m_2}+\frac18\,\frac{m_1+m_2}{m_1-m_2}\right ] \ln\frac{m_2}{m_1}\right ) \right]\delta^3(r)\cr\nonumber\\ & & + \frac{\alpha_s(\mu^2)}{\pi}\left[- \frac{\beta_0}{8\pi}\nabla^2\left(\frac{\ln({\mu } r)+\gamma_e}{r}\right ) + \frac{21}{16\pi}\nabla^2\left ( \frac{\ln(\sqrt{m_1m_2}\,r)+\gamma_e}{r}\right)\right]\biggr)\cr\nonumber\\ & & + ( 1+\kappa)^2(1-\varepsilon)\frac{a}{r}\biggr\},\\ \cr \label{d } d&= & \frac14\left(\frac{1}{m_1 ^ 2}-\frac{1}{m_2 ^ 2}\right)\biggl\ {    \frac43\frac{\alpha_s(\mu^2)}{r^3}\biggl(1 +   \frac{\alpha_s(\mu^2)}{\pi}\biggl[\frac73-\frac{\beta_0}{12}+    \gamma_e\left(\frac{\beta_0}{2}-3\right)+\frac{\beta_0}{2}\ln(\mu r ) \cr\nonumber\\ & & -3\ln(\sqrt{m_1m_2}\,r)\biggr]\biggr ) -\frac{a}{r}-2(1+\kappa)(1-\varepsilon)\frac{a}{r}\biggr\ } + \left(\frac{1}{m_1}+\frac{1}{m_2}\right)^2    \frac{\alpha_s^2(\mu^2)}{2\pi r^3}\ln\frac{m_2}{m_1 } ,     \end{aligned}\ ] ] where @xmath75 is the orbital momentum and @xmath76 , @xmath77 are the spin momenta . for the equal mass case ( @xmath78 ) the second order in @xmath79 contribution of the annihilation diagrams @xmath80 must be added to the spin - spin interaction coefficient @xmath81 in eq . ( [ c ] ) . the correct description of the fine structure of the heavy quarkonium mass spectrum requires the vanishing of the vector confinement contribution . this can be achieved by setting @xmath82 , i.e. , the total long - range quark chromomagnetic moment equals zero , which is in accord with the flux tube @xcite and minimal area @xcite models . one can see from eq . ( [ vsd ] ) that for the spin - dependent part of the potential this conjecture is equivalent to the assumption about the scalar structure of confining interaction @xcite .    to calculate the heavy meson mass spectra with the account of all relativistic corrections ( including retardation effects ) of order @xmath1 and one - loop radiative corrections we substitute the quasipotential which is a sum of the spin - independent ( [ sipot ] ) and spin - dependent ( [ vsd ] ) parts into the quasipotential equation ( [ quas ] ) . then we multiply the resulting expression from the left by the quasipotential wave function of a bound state and integrate with respect to the relative momentum . taking into account the accuracy of the calculations , we can use for the resulting matrix elements the wave functions of eq . ( [ quas ] ) with the static potential ) , also not vanishing in the static limit , is treated perturbatively . ] @xmath83 as a result we obtain the mass formula @xmath84 \rangle + \langle c\rangle \langle{\bf s}_1\cdot { \bf s}_2\rangle + d \langle { \bf l}\cdot({\bf s}_1-{\bf s}_2)\rangle , \ ] ] where @xmath85 the first term on the right - hand side of the mass formula ( [ mform ] ) contains all spin - independent contributions , the second and the last terms describe the spin - orbit interaction , the third term is responsible for the tensor interaction , while the forth term gives the spin - spin interaction . the last term is not zero only for the unequal mass case @xmath86 and leads to the mixing of triplet and singlet meson states with the total angular momentum @xmath87 equal to the orbital momentum @xmath88 .    in table [ charm ] the calculated charmonium mass spectrum is compared with experimental data . for meson states we use the notation @xmath89 , where @xmath90 and @xmath91 is the radial quantum number . our predictions agree with pdg @xcite data within a few mev . our model correctly reproduces both the position of the levels and their fine and hyperfine splittings . in this table we also give the recent belle collaboration data @xcite on pseudoscalar @xmath92 and @xmath93 states observed in @xmath16 decays . the measured mass of the ground spin - singlet state @xmath92 is in good agreement with world averages and predictions of our model , while the radially excited @xmath93 state lies considerably higher than previous experimental indications and most of the theoretical predictions . if these data are confirmed , it will be difficult to accommodate such a small hyperfine splitting @xmath94  mev ( almost four times smaller than 117  mev splitting for the ground state ) in the framework of the quark model . can in principle be influenced by the nearby threshold of the open charm production . ] .charmonium mass spectrum ( in gev ) . [ cols="^,^,^,^,^ " , ]     in tables  [ tab : e1bc ] and [ tab : e1bc2 ] we compare our results for the e1 radiative decay rates of the @xmath0 meson with other quark model predictions @xcite . comparison of the calculations using relativistic @xmath95 and nonrelativistic @xmath96 formulae for decay rates shows that relativistic corrections do not exceed 20% in @xmath0 meson e1 decays . most of the theoretical predictions for e1 transitions between @xmath3 and @xmath2 states of @xmath0 mesons given in table  [ tab : e1bc ] are compatible with each other . the largest differences occur for decays involving @xmath97 and @xmath98 states which are the mixtures of spin singlet and spin triplet states ( [ eq : pmix ] ) due to different mixing angles used by the authors . note that for such transitions there are additional relativistic spin - flip contributions to decay rates ( [ eq : edr ] ) proportional to @xmath99 ( [ eq : chi ] ) which are specific only for @xmath0 mesons . in general , our predictions are closer to the ones of ref .  @xcite .    in table  [ tab : e1bc2 ] we present the e1 radiative decay rates of @xmath0 mesons where either the initial or final state is a @xmath4 wave state . here we find rather large variations in theoretical predictions . the main reason of these distinctions is the difference in values of @xmath4 state masses , which for some states reaches 70  mev ( see table  [ tab : bcm ] ) . since @xmath100 and @xmath101 states of the @xmath0 are rather close , such a difference significantly influences the energy of the emitted photon and thus the decay rates . for decays involving the mixed spin singlet and spin triplet states @xmath97 , @xmath98 ( [ eq : pmix ] ) and @xmath102 , @xmath103 ( [ eq : dmix ] ) the additional relativistic spin - flip contributions ( [ eq : edr ] ) are important , especially for transitions where both initial and final states are mixed states . in this paper the mass spectra and radiative m1 and e1 decay rates of charmonium , bottomonium and @xmath0 mesons were calculated in the relativistic quark model based on the quasipotential approach in quantum field theory . special attention was devoted to the role of the relativistic effects in these processes . since both quarks in the considered mesons are heavy , the @xmath104 expansion was applied . in the mass spectra calculations retardation as well as one - loop radiative corrections were taken into account . we also included the one - loop correction due to the finite charm quark mass to the bottomonium mass spectrum . it was found that this correction is rather small @xcite and its inclusion allows one to obtain an even better fit of the bottomonium excited states with the slightly shifted value of qcd parameter @xmath105 . the calculated mass spectra of charmonium and bottomonium agree with the experimental data within a few mev . comparison of our results for the @xmath0 meson mass spectrum with previous calculations showed that different predictions for ground state and low excitations agree within 30  mev . the pseudoscalar and vector decay constants of @xmath0 meson were calculated using the relativistic wave functions obtained during the mass spectrum calculations . it was found that relativistic effects reduce these constants by approximately 20% and produce the splitting between them of about 70  mev . it was shown that relativistic effects play a significant role in radiative decays of mesons . their form strongly depends on the lorentz structure of the quark - antiquark interaction . the most sensitive are radiative m1 decays , where even for allowed transitions they significantly influence predictions for the rates . an important example is the @xmath106 transition , which is overestimated by a factor of more than two if the nonrelativistic approximation is used . it is argued that the inclusion of relativistic corrections for a pure scalar or vector confining potential is not enough to bring theoretical predictions in accord with experiment . only for the specific mixture of these potentials ( [ vlin ] ) with the mixing coefficient @xmath50 , the agreement can be obtained . for other decay rates this mixing of scalar and vector potentials also gives the best results . the hindered m1 transition rates are dominated by relativistic contributions and are significantly enhanced by them . the comparison of the allowed and hindered m1 rates in bottomonium shows that the latter provide better opportunity for the search of the missing pseudoscalar @xmath107 state of the bottomonium . the analysis of radiative e1 transitions showed that the form of relativistic corrections is less dependent on the lorentz structure of the quark potential than in the case of m1 transitions . however , for some decays , e.g. , @xmath108 , the consideration of the mixed ( [ vlin ] ) vector and scalar potentials ( with the same value of @xmath50 ) is important for bringing decay rates in accord with experimental data . in general , all our predictions for radiative decay rates and branching fractions of charmonium and bottomonium agree with measured values . in the case of the @xmath0 meson radiative e1 decays an important additional relativistic correction to decay rates which causes the flip of the quark spin was found . this contribution ( [ eq : chi ] ) to the radiative e1 decay rate ( [ eq : edr ] ) is specific only for transitions involving mixed states @xmath109 , @xmath110 ( [ eq : pmix ] ) or @xmath111 , @xmath112 ( [ eq : dmix ] ) of @xmath0 and is caused by the difference of the @xmath81 and @xmath113 quark masses . finally , a comparison of various quark model predictions for the radiative m1 and e1 decay rates of @xmath0 has been performed . these radiative transitions along with pionic ones are the main decay channels of the low lying excitations in the @xmath0 meson . the authors express their gratitude to a. martin , m. mller - preussker and v. savrin for support and discussions . we are grateful to j. l. rosner for stimulating correspondence and quoting some of our preliminary results on bottomonium radiative transitions in ref . two of us ( r.n.f and v.o.g . ) were supported in part by the _ deutsche forschungsgemeinschaft _ under contract eb 139/2 - 1 , _ russian foundation for basic research _ under grant no . 00 - 02 - 17768 and _ russian ministry of education _ under grant no . e00 - 3.3 - 45 . 99 e. j. eichten , k. lane and c. quigg , phys . lett . * 89 * , 162002 ( 2002 ) . m. suzuki , phys . d * 66 * , 037503 ( 2002 ) . y .- kuang , phys . d * 65 * , 094024 ( 2002 ) . s. godfrey and j. l. rosner , phys . d * 64 * , 074011 ( 2001 ) ; * 65 * , 039901(e ) ( 2002 ) . s. godfrey and j. l. rosner , phys . d * 66 * , 014012 ( 2002 ) . s. godfrey and j. l. rosner , phys . d * 64 * , 097501 ( 2001 ) . cdf collaboration , f. abe _ d * 58 * , 112004 ( 1998 ) . a. a. logunov and a. n. tavkhelidze , nuovo cimento * 29 * , 380 ( 1963 ) . a. p. martynenko and r. n. faustov , theor . phys . * 64 * , 765 ( 1985 ) [ teor . mat . fiz . * 64 * , 179 ( 1985 ) ] . v. o. galkin , a. yu . mishurov and r. n. faustov , yad . * 55 * , 2175 ( 1992 ) [ sov . phys . * 55 * , 1207 ( 1992 ) ] . d. ebert , r. n. faustov and v. o. galkin , phys . d * 62 * , 034014 ( 2000 ) . e. eichten and f. feinberg , phys . d * 23 * , 2724 ( 1981 ) . d. ebert , v. o. galkin and r. n. faustov , phys . d * 57 * , 5663 ( 1998 ) ; * 59 * , 019902(e ) ( 1999 ) . v. o. galkin and r. n. faustov , yad . fiz . * 44 * , 1575 ( 1986 ) [ sov . 44 * , 1023 ( 1986 ) ] ; d. ebert , r. n. faustov and v. o. galkin , phys . b * 537 * , 241 ( 2002 ) . r. n. faustov and v. o. galkin , z. phys . c * 66 * , 119 ( 1995 ) . n. brambilla , p. consoli and g. m. prosperi , phys . d * 50 * , 5878 ( 1994 ) . s. n. gupta and s. f. radford , phys . d * 24 * , 2309 ( 1981 ) ; ibid . * 25 * , 3430 ( 1982 ) ; s. n. gupta , s. f. radford and w. w. repko , phys . d * 26 * , 3305 ( 1982 ) . j. pantaleone , s .- h . h. tye and y. j. ng , phys . d * 33 * , 777 ( 1986 ) . w. buchmller , phys . b * 112 * , 479 ( 1982 ) . n. brambilla and a. vairo , phys . d * 55 * , 3974 ( 1997 ) . h. j. schnitzer , phys . lett . * 35 * , 1540 ( 1975 ) ; w. lucha , f. f. schberl and d. gromes , phys . rep . * 200 * , 127 ( 1991 ) ; v. d. mur , v. s. popov , yu . a. simonov and v. p. yurov , zh . 78 * , 1 ( 1994 ) [ j. exp . . phys . * 78 * , 1 ( 1994 ) ] ; a. yu . dubin , a. b. kaidalov and yu . a. simonov , phys . b * 323 * , 41 ( 1994 ) ; yu . a. simonov , phys . usp . * 39 * , 313 ( 1996 ) . particle data group , k. hagiwara _ et al . d * 66 * , 010001 ( 2002 ) . belle collaboration , s .- k . et al . _ , phys . 89 * , 102001 ( 2002 ) . c. edwards _ et al . lett . * 48 * , 70 ( 1982 ) . t. a. armstrong _ et al . _ , 69 * , 2337 ( 1992 ) . e705 collaboration , l. antoniazzi _ d * 50 * , 4258 ( 1994 ) . a. h. hoang and a. v. manohar , phys . b * 491 * , 101 ( 2000 ) . a. h. hoang , hep - ph/0008102 , hep - ph/0204299 . m. melles , phys . d * 62 * , 074019 ( 2000 ) . n. brambilla , y. sumino and a. vairo , phys . d * 65 * , 034001 ( 2002 ) . d. ebert , r. n. faustov and v. o. galkin , phys . d * 66 * , 037501 ( 2002 ) . cleo collaboration , s. e. csorna _ et al . _ , hep - ex/0207060 ; d. cinabro _ et al . _ , hep - ex/0207062 . e. j. eichten and c. quigg , phys . d * 49 * , 5845 ( 1994 ) . s. s. gershtein , v. v. kiselev , a. k. likhoded , and a. v. tkabladze , phys . d * 51 * , 3613 ( 1995 ) . l. p. fulcher , phys . d * 60 * , 074006 ( 1999 ) . s. nussinov and m. a. lampert , phys . * 362 * , 193 ( 2002 ) . w. kwong and j. l. rosner , phys . d * 44 * , 212 ( 1991 ) . n. brambilla and a. vairo , phys . d * 62 * , 094019 ( 2000 ) . d. ebert , r. n. faustov and v. o. galkin , mod . a * 17 * , 803 ( 2002 ) ; v. o. galkin , a. yu . mishurov and r. n. faustov , yad 53 * , 1676 ( 1991 ) [ sov . . phys . * 53 * , 1026 ( 1991 ) ] . b. d. jones and r. m. woloshyn , phys . d * 60 * , 014502 ( 1999 ) . r. n. faustov , ann . ( n. y. ) * 78 * , 176 ( 1973 ) ; nuovo cimento a * 69 * , 37 ( 1970 ) . cleo collaboration , a. h. mahmood _ et al . _ , hep - ex/0207057 . v. o. galkin , a. yu . mishurov and r. n. faustov , yad . fiz . * 51 * , 1101 ( 1990 ) [ sov . phys . * 51 * , 705 ( 1990 ) ] . j. l. rosner , phys . d * 64 * , 094002 ( 2001 ) . v. a. novikov _ et al . _ , phys . rep . * 41 * , 1 ( 1978 ) ; r. barbieri , r. gatto and r. kogerler , phys . b * 60 * , 183 ( 1976 ) ; r. barbieri , m. caffo , r. gatto , and e. remiddi , nucl . b * 192 * , 61 ( 1981 ) . w. kwong , p. b. mackenzie , r. rosenfeld , and j. l. rosner , phys . d * 37 * , 3210 ( 1988 ) . r. n. faustov and v. o. galkin , proceedings of the international seminar `` quarks88 '' , tbilisi , 17 - 21 may 1988 ( world scientific , singapore , 1989 ) , p. 624 . y .- p . kuang and t .- m . yan , phys . d * 24 * , 2874 ( 1981 ) ; m. voloshin , yad . fiz . * 43 * , 1571 ( 1986 ) [ sov . * 43 * , 1011 ( 1986 ) ] .
Abstract:

Answer: the mass spectra and electromagnetic decay rates of charmonium , bottomonium and @xmath0 mesons are comprehensively investigated in the relativistic quark model . the presence of only heavy quarks allows the expansion in powers of their velocities . all relativistic corrections of order @xmath1 , including retardation effects and one - loop radiative corrections , are systematically taken into account in the computations of the mass spectra . the obtained wave functions are used for the calculation of radiative magnetic dipole ( m1 ) and electric dipole ( e1 ) transitions . it is found that relativistic effects play a substantial role . their account and the proper choice of the lorentz structure of the quark - antiquark interaction in a meson is crucial for bringing theoretical predictions in accord with experimental data . a detailed comparison of the calculated decay rates and branching fractions with available experimental data for radiative decays of charmonium and bottomonium is presented . the possibilities to observe the currently missing spin - singlet @xmath2 and @xmath3 states as well as @xmath4 states in bottomonium are discussed . the results for @xmath0 masses and decays are compared with other quark model predictions .