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: many body perturbation theory ( mbpt ) diagrammatic approaches offer a path to the solution of the quantum many body problem of solids , which complements alternative methods such as qmc ( quantum monte carlo ) and cc ( couple cluster ) . this approach received attention for over half a century but there is no complete understanding of how different selection of diagrams performs for different physical quantities . these insights are important in the search for predictive first principles methods for correlated solids . in this work we investigate these questions in the framework of the homogeneous electron gas ( heg ) in a neutralizing positively charged background . the heg is very useful for testing in a simplified setting the methods presently being developed to study the electronic structure of solids for two reasons : it requires less computational effort and there is a natural benchmark since some properties have also been calculated using qmc methods @xcite . this model describes well the properties of alkali metals . most common uses of diagrammatic approaches are based on non - interacting green functions such as a lda ( local density approximation ) green s function . in this work , however , we are interested in self consistent ( sc ) diagrammatic approaches . hence we do not consider ambiguities related to the choice of the non interacting green function for reference . we consider two classes of methods . one initiated by hedin@xcite , which carries out a perturbative expansion in the fully self consistent ( renormalized ) green function which obeys the dyson equation and in screened interaction w. an alternative philosophy ( qsgw ) uses the lowest order diagrams for the polarizability and the self energy with a green s function which is determined by means of a quasi - particle self consistency condition @xcite . in this work we perform qsgw calculations using a previously introduced linearized approach@xcite that has an advantage of being implementable on matsubara frequencies . an important cornerstone in the application of the self - consistent mbpt - based approximations to the heg is the work by holm and von barth @xcite where the authors applied scgw approximation ( i.e. lowest order diagram in the perturbation series in terms of g and w ) to calculate the total energy and spectra of the heg . their principal conclusions are that scgw severely overestimates the band width but gives the total energy very close to the qmc results . thus , the work @xcite raised the question of how ( if at all ) one can get accurate spectra of the heg using scmbpt approach . progress on this question was made by shirley @xcite who showed that if an accurate w is known ( he used qmc input to evaluate w ) then the self - consistency in the green s function g and the lowest order vertex corrections in the self energy mostly cancel out each other justifying in a certain degree the so called one - shot @xmath0 approaches ( lowest order approach using a non interacting greens function @xmath1 ) without vertex corrections to the self energy . takada @xcite used qmc data to parametrize electron - hole four - point irreducible interaction and performed sc vertex - corrected gw calculations with the vertex defined from a condition that the ward identity ( wi ) is satisfied . in this respect his approach resembles the idea of qsgw method where wi is imposed by construction . the importance of vertex corrections in electron gas studies was reported recently in non - self - consistent calculations.@xcite from other studies on the subject , one should mention the work [ ] , where the spectral properties of the electron gas were studied with ab - initio coupled - cluster theory .    in this study we go beyond the earlier diagrammatically inspired works by removing the important limitation related to the use of monte - carlo data for the parametrization of the screened interaction or the electron - hole four - point irreducible interaction . we examine the fully self - consistent ( in both g and w ) diagrammatic schemes involving diagrams of higher order . we study the relative importance of different diagrams for the polarizability and the self energy . the paper begins with a brief presentation of the self - consistent schemes we use to solve the hedin s equations ( section [ meth ] ) . section [ res ] provides the results obtained and a discussion . the conclusions are given afterwords . detailed account of the vertex - corrected schemes we use in this work has been given before in ref . [ ] . for completeness , below we briefly repeat the essentials of the approach and point out the simplifications in technical implementation in the case of the electron gas . we solve the hedin s equations @xcite self - consistently using different approximations for the three point vertex function @xmath2 . the three point vertex function enters the formally exact expressions for the polarizability and the self energy ( in space - time variables )    @xmath3    @xmath4    where the integration / summation over repeated arguments is understood , and @xmath5 is the spin index . we consider two different types of approximations for @xmath2 . the first type consists in expanding the vertex function in terms of the screened interaction to a specified order . keeping only zero order term ( @xmath6 ) in both @xmath7 and @xmath8 corresponds to the famous gw approximation . we will also consider the expansion of the vertex up to the first order ( @xmath9 ) in both the polarizability and the self energy expressions . this approximation is conserving ( like gw ) as the corresponding p and @xmath8 can alternatively be obtained by differentiating the same @xmath10-functional @xcite . the second type of approximation for @xmath2 consists in solving the bethe - salpeter equation @xmath11 with a certain approximate expression for the functional derivative @xmath12 in ( [ vert_0 ] ) . we will consider two expressions for the kernel @xmath13 in this work . the first is obtained by using the gw form for @xmath8 in the functional derivative and neglecting the derivative of the screened interaction @xmath14 , i.e. @xmath15 ( we will call the corresponding vertex as @xmath16 ) . diagrammatically it corresponds to keeping only the first term on the right hand side of fig.[theta_appr ] . in the second approximate expression for @xmath13 we also are using the gw form for the self energy in the functional derivative but we keep the terms up to the second order in @xmath17 in the derivative @xmath14 ( we will use abbreviation @xmath18 for the corresponding vertex ) . in this case , the obtained vertex function corresponds to keeping all three terms for @xmath13 ( fig.[theta_appr ] ) . it is important to point out that the diagrams resulted from @xmath14 allow the spin flips ( as it is clear from fig.[theta_appr ] ) , the importance of which was pointed out in ref . [ ] .    .,width=7 ]    in the particular case , when g and w have been found self - consistently with @xmath19 and @xmath20 , the vertex @xmath18 yields the physical polarizability in scgw approximation ( defined as a functional derivative of electronic density with respect to the total electric field ) . in other cases , when @xmath18 is evaluated with @xmath8 and/or @xmath7 including additional diagrams , the kernel shown in fig . [ theta_appr ] is only an approximation to the derivative @xmath21 , and as a result , the vertex @xmath18 does nt provide physical @xmath7 anymore . essentially , in a search for an optimal approximation , we have to trade between the accuracy of @xmath8 in the dyson equation for @xmath22 ( which is generally better when we increase the number of diagrams beyond simple @xmath19 ) and the degree of the `` deviation '' of the polarizability from the physical one . .diagrammatic representations of the polarizability and the self energy in sc schemes of solving the hedin equations . arguments in square brackets specify @xmath22 and @xmath17 which are used to evaluate the vertex function . other details are explained in the main text . [ cols="^,^,^",options="header " , ]     in table [ bw_comp ] we compare our results for the band width with those obtained by shirley @xcite , who based the calculations ( partially ) on qmc input . band width was determined as the difference between the pole in the spectral function at * * k**=0 and chemical potential . if we assume that shirley s results are close to the exact ones , we can draw certain conclusions about our approaches . as one can conclude from the table , three schemes ( d , e , and g ) show the best performance with small difference between themselves . common for these three schemes are two facts : they all include the solving of the bethe - salpeter equation for the polarizability ( but slightly differently as it was explained above ) and they all apply first order vertex correction in the self energy . scheme c also involves the solving of bse for the polarizability , but it does nt use vertex correction to the self energy and , as a result , shows worse performance . similarly , the conserving scheme b , which applies first order vertex corrections to the @xmath7 and @xmath8 , shows worse performance , because it misses the effects of bse in the w. nevertheless , scheme b seems to be better than scheme c , demonstrating the importance of vertex corrections in the self energy . the performance of qsgw approach is slightly better than performance of scgw , but it is not competitive with schemes d , e , or g. thus , for the studied range of densities of the electron gas , qsgw can not be considered as a reasonable approximation in terms of its predictive power . .,width=8 ]     obtained with different approximations in comparison with qmc result@xcite . the results based on local field factors ( lff ) are taken from the ref.[].,width=8 ]    . the results based on local field factors ( lff ) are taken from the ref.[].,width=8 ]    summarizing our results for one - electron properties , we show them ( for the purpose of reference ) for the spectrum ( @xmath23=3 ) , quasiparticle renorm - factor z and effective electron mass @xmath24 in figures [ bw ] , [ z_factor ] , and [ m_eff ] correspondingly . the spectrum was defined as the positions of the poles in spectral function ( for every * k * ) . in order to evaluate the renorm - factor and the effective mass , we used the formulae ( with @xmath25 being the fermi momentum and @xmath26 being the chemical potential ) @xmath27 and @xmath28 correspondingly . for z - factor there are qmc data @xcite but the data involve extrapolation to the thermodynamic limit and do not look as accurate enough . instead , we compare our results with the ones based on local field factors ( lff ) @xcite which include certain amount of qmc input . it is interesting that the dependence on @xmath23 of the effective mass is weaker than similar dependence of the z - factor suggesting that the frequency derivative and momentum derivative in eq.([ms ] ) are canceling out considerably .    . comparison with qmc data from ref . schemes c and d give identical results for this quantity.,width=8 ]    in figure [ e1_4 ] the static ( zero - frequency ) inverse dielectric function for 3d electron gas is shown for @xmath29 . as one can see from the graph we are able to improve the agreement with qmc data considerably when using the vertices of increased complexity . clearly , the best dielectric function is obtained from the `` physical '' polarizability ( schemes c and d ) , even if the last corresponds to the scgw approximation ( in the sense that the diagrams are evaluated using g and w from scgw ) . one can also point out the importance of including the diagrams with spin flips ( scheme e results in better dielectric function than scheme g does ) . at the same time scheme e is worse than schemes c / d which reflects the above mentioned fact about trading - additional diagram in the self energy in scheme e violates the requirement of the polarizability to be physical . one can also relate the shortcomings of qsgw approach in the one - electron spectra to the poor description of the screening . as one can see from figure [ e1_4 ] , in the physically important range ( @xmath30=0.5 - 3.0 ) the inverse dielectric function in qsgw approximation is larger than the one from qmc data by a factor 2 - 3 .    . @xmath31 is the compressibility of the non - interacting electron gas . to make the plot , we used the vertex @xmath32 in schemes c , d , and e , and the vertex @xmath33 in scheme g.,width=8 ]     at small @xmath34 for @xmath29 . symbols @xmath35 show the calculated data points , line is drawn for convenience.,width=8 ]    in fig.[compress ] we compare our calculated compressibilities with the results based on the qmc data . we have obtained the qmc compressibility @xmath36 as the derivative of the chemical potential @xmath26 ( @xmath37 ) with respect to @xmath23 : @xmath38 . in our vertex - corrected schemes , we have evaluated the compressibility from the ratio of two limits of the vertex function ( @xmath39 and @xmath40 ) and the effective mass @xmath24 ( see for instance ref . [ ] )    @xmath41    the two limits of the vertex function , @xmath39 and @xmath40 , were defined as follows . in the momentum - frequency representation , the vertex function ( solution of the eq.([vert_0 ] ) ) can be conveniently considered as dependent on fermionic momentum+frequncy ( * k*,@xmath42 ) and on bosonic momentum+frequency ( * q*,@xmath34 ) , i.e. @xmath2(*k*,@xmath42;*q*,@xmath34 ) . in these variables , @xmath43 and @xmath44 , with two vectors , * k * and * q * , assumed to be parallel . quantities @xmath39 and @xmath40 are related to the quasiparticle renormalization factor @xmath45 and the landau fermi liquid parameter @xmath46 : @xmath47 , @xmath48 . in order to demonstrate that the above two limits are well defined numerically , we have plotted the vertex function @xmath32 from scheme d in the fig.[vrt_plot ] . the limit @xmath39 corresponds to the value of the function at @xmath49 exactly ( approximately 4.58 on the graph ) . the limit @xmath40 corresponds to the extrapolation of the function to zero frequency ( approximately 1.31 on the graph ) . for the purpose of reference , the calculated fermi liquid parameter @xmath46 is shown in fig.[f0s ] . again , the limits @xmath39 and @xmath40 have been used to obtain the data for the plot .     of the electron gas calculated with vertex - corrected schemes.,width=8 ] it is clear from the figure [ compress ] that the behavior consistent with the presence of the pole in the compressibility ( at about @xmath50 , `` dielectric catastrophe''@xcite ) can only be obtained based on the vertex from schemes c and/or d. in other words , only the vertex corresponding to the physical polarizability can be useful in reproducing the qmc compressibility . it is interesting , that the additional self - consistency iterations for the vertex function in schemes e and g ( as compared to scheme d ) which only slightly change the one - electron spectra , worsen significantly the quality of the calculated compressibility . also , it is obvious that the first order vertex ( scheme b ) is totally insensitive to the presence of a pole in the compressibility .    ) as a function of bosonic momentum * q * for @xmath29 . fermionic frequency and momentum correspond to their values at the fermi surface . schemes c and d give identical results for this quantity.,width=8 ]    an additional insight on the origin of differences in the dielectric function obtained with approximate methods can be gained when one looks at the static vertex as a function of bosonic momentum ( fig . [ vrt_stat ] ) . it is clear that the range of momenta where the calculated dielectric function shows the largest differences ( @xmath51 ) correlates very well with the range of momenta where the humps in the vertex function show very different heights . .,width=8 ]    in fig.[e_c ] we present the calculated correlation energy of the electron gas as a function of @xmath23 . it was obtained as the difference between the expectation value of the hamiltonian corresponding to the selected level of approximation and the expectation value of the hamiltonian in the hartree - fock approximation . in all vertex - corrected schemes the exchange - correlation part was evaluated as a convolution of the green s fucntion and the self energy . excluding scheme c which misses the vertex correction to the self energy and , as a result , shows rather different from other schemes behavior , one can conclude that vertex corrections do not improve the result obtained in scgw approach . only at @xmath23 less than 2 there is a small improvement . however , the description of the slope of the curve is improved . it is especially easy to see in fully self - consistent schemes b , e , and g. for these three schemes there is an almost @xmath23-independent deviation from qmc data , the origin of which is not clear at the moment . in conclusion we have applied the self - consistent diagrammatic approaches based on the hedin equations to study the properties of the 3d heg . we have found , that the inclusion of the most important diagrammatic sequences can reproduce with a good accuracy the one - electron spectra and dielectric properties of the heg in the range of metallic densities . for the one - electron spectra , the corrections to the polarizability and to the self energy are equally important . for the dielectric properties the vertex correction to the self energy is of secondary importance . in all cases , the important conclusion is that the calculation of the polarizability should follow , as close as possible , its definition as a functional derivative of the density with respect to the total electric field . our conclusions concerning one - electron spectra of 3d electron gas are similar to the conclusions made earlier for the spectra of alkali metals and semiconductors in refs . [ ] , namely , that the best spectra are obtained when the set of diagrams for the polarizability is obtained from bse , whereas the first order vertex correction is applied to the self energy ( schemes d , e , and g ) . our benchmarks quantified the inaccuracy the qsgw approximation to predict one - electron spectra of the electron gas at metallic densities ( approximately 15% error ) . we track this inaccuracy to the poor description of screening in qsgw approach ( with an error up to a factor 2 - 3 in the physically important range of momenta ) . concerning the vertex - corrected schemes , one can advocate using scheme d , which combines good accuracy and computational efficiency ( time - consuming bse has to be solved only once ) . this work was supported by the u.s . department of energy , office of science , basic energy sciences as a part of the computational materials science program . gk was supported by the simons foundation under the many electron problem collaboration .
Abstract:

Answer: a few approximate schemes to solve the hedin equations @xcite self - consistently introduced in ( phys . rev . b 94 , 155101 ( 2016 ) ) are explored and tested for 3d electron gas at metallic densities . we calculate one electron spectra , dielectric properties , compressibility , and correlation energy . considerable reduction in the calculated band width ( as compared to the self consistent gw result ) has been found when vertex correction was used for both the polarizability and the self energy . generally , it is advantageous to obtain the diagrammatic representation of the polarizability from the definition of this quantity as a functional derivative of electronic density with respect to the total field ( external plus induced ) . for the self energy , the first order vertex correction seems to be sufficient for the range of densities considered . whenever it is possible , we compare the accuracy of our vertex - corrected schemes with the accuracy of the self - consistent quasi - particle gw approximation ( qsgw ) , which is less expensive computationally . we show that qsgw approach performs poorly and we relate this poor performance with inaccurate description of the screening in the qsgw method ( with an error comprising a factor 2 - 3 in the physically important range of momenta ) .


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: ooperative communication is one of the key parts of the future communication protocols , as the deployment of relays can improve link equality , extend coverage range and mitigate inference . in general , there are various relay strategies which are casted into three main categories , i.e. , amplify - and - forward ( af ) , decode - and - forward ( df ) and compressed - and - forward ( cf ) . among these relaying strategies , af strategy which has the lowest complexity is most suitable for practical implementation . it is also well - established that adopting multiple antennas has a potential to improve overall wireless system performance . in order to reap both benefits promised by cooperative communication and multi - input multi - output ( mimo ) systems , linear transceiver design for af mimo relaying systems has been widely researched in @xcite . generally , speaking there are two main kinds of criteria for transceiver design : capacity maximization and mean - square - error ( mse ) minimization . joint design of relay forwarding matrix and destination equalizer for minimizing mse is discussed in @xcite and @xcite . furthermore , joint design of source precoder relay forwarding matrix and destination equalizer for minimizing mse is investigated in @xcite . the capacity maximization transceiver design has been discussed in @xcite .     in most of previous works , channel state information ( csi ) is assumed to be perfectly known . however , this assumption can not be met in practice . channel estimation errors are always inevitable and drastically degrades system performance . it is well known that robust designs can reduce or mitigate the negative effects introduced by imperfect csi . this is also the motivation of our work . in this paper , we jointly optimize source precoder matrix and relay forwarding matrix for mutual information maximization under channel estimation errors . based on the properties of matrix - monotone functions , the optimal structure of robust transceivers is derived . then , an iterative water - filling solution is proposed . finally , the numerical result shows the performance advantage of the proposed robust design . the following notations are used throughout this paper . boldface lowercase letters denote vectors , while boldface uppercase letters denote matrices . the notation @xmath0 denotes the hermitian of the matrix @xmath1 , and @xmath2 is the trace of the matrix @xmath1 . the notation @xmath3 is the hermitian square root of the positive semi - definite matrix @xmath1 , such that @xmath4 and @xmath3 is also a hermitian matrix . for a rectangular diagonal matrix @xmath5 , @xmath6 denotes the main diagonal elements are in decreasing order and @xmath7 denotes the main diagonal elements are in increasing order . for two hermitian matrices , @xmath8 means that @xmath9 is a positive semi - definite matrix . the symbol @xmath10 represents the @xmath11 largest eigenvalue of @xmath1 . in our work , a dual - hop af mimo relay system is investigated , in which there is one source with @xmath12 antennas , one relay with @xmath13 receive antennas and @xmath14 transmit antennas , and one destination with @xmath15 antennas . because of deep fading , the direct link between the source and destination is not taken into account . at the first hop , the source transmits data to the relay . the received signal at the relay is denoted as @xmath16where @xmath17 is the mimo channel matrix between the source and the relay , and @xmath18 is the precoder matrix at the source . the vector @xmath19 is the @xmath20 data vector transmitted by the source with the covariance matrix @xmath21 . furthermore , @xmath22 is the additive gaussian noise vector with correlation matrix @xmath23 .    at the relay , the received signal @xmath24 is multiplied by a forwarding matrix @xmath25 . then the resultant signal is transmitted to the destination . the received signal @xmath26 at the destination can be written as @xmath27where @xmath28 is the mimo channel matrix between the relay and the destination , and @xmath29 is the additive gaussian noise vector at the second hop with covariance matrix @xmath30 . in order to guarantee the transmitted data @xmath19 can be recovered at the destination , it is assumed that @xmath12 , @xmath31 , @xmath14 , and @xmath15 are greater than or equal to @xmath32 @xcite . when channel estimation errors are taken into account , the dual - hop channels read as @xmath33where @xmath34 and @xmath35 are the channel estimates and @xmath36 and @xmath37 are the corresponding estimation errors with zero - mean gaussian distributed entries . additionally , the estimation errors are independent with each other as the channels are separately estimated . referring to estimation errors , the following widely used kronecker structure is adopted @xcite @xcite @xmath38 where the entries of @xmath39 and @xmath40 are identical and independent distributed ( i.i.d . ) with zero mean and unit variance . the column correlation matrices ( @xmath41 and @xmath42 ) and the row correlation matrices ( @xmath43 and @xmath44 ) are determined by training sequences and channel estimators @xcite . to the best of the authors knowledge , for a general case even for a point - to - point mimo system , there is no closed - form solution . in this paper , we focus on the case with column correlations only i.e. , @xmath45as this case corresponds to a practical linear minimum mean square error ( lmmse ) channel estimator @xcite .      at the destination , a linear equalizer @xmath46 is adopted to detect the data vector @xmath19 . the mean - square - error ( mse ) matrix is @xmath47 , where the expectation is taken with respect to random data , channel estimation errors , and noise . in @xcite , it is shown that @xmath48where matrices @xmath49 and @xmath50 are defined as@xmath51it is obvious that @xmath49 is the covariance matrix of the received signal at the relay . using linear bayesian theory , the lmmse equalizer at the destination equals to @xmath52based on which the mse matrix in ( [ mse_final ] ) is rewritten as @xmath53    capacity is one of the most important and widely used performance metrics for transceiver designs . given the received pilots in channel estimation as @xmath54 and @xmath55 , the channel capacity is denoted as @xmath56 , which is the conditional mutual information based on known @xmath54 and @xmath55 @xcite . to the best of our knowledge , the exact capacity for mimo channels with estimation errors at both ends is largely open even for point - to - point mimo systems @xcite . to proceed , a common logic is to derive and use bounds i.e. , lower bound or upper bound . since we aim to maximize channel capacity , lower bound is more meaningful than upper bound . in appendix  [ appedix:5 ] it has been proved that @xmath57 this is a widely - established lower bound @xcite and becomes tighter as estimation errors are smaller .    as a result , the robust transceiver design for maximizing mutual information is formulated as @xmath58 based on the definition of @xmath49 in ( [ r_x ] ) , @xmath59 is a function of @xmath18 . in order to simplify the analysis , we define a new variable@xmath60based on which @xmath61 and the two constraints involved in ( [ prob : opt ] ) become independent . meanwhile , the mse matrix in ( [ mse_matrix_orig ] ) is correspondingly rewritten as @xmath62finally , the optimization problem for the robust design becomes @xmath63 in the following , the optimal solutions of ( [ opt : capacity ] ) will be discussed in details . in our work , we investigate the optimization problem ( [ opt : capacity ] ) from matrix - monotone function viewpoint . the idea of utilizing the properties of matrix - monotone functions to design mimo transceivers has been address in @xcite . in this paper , we extend this idea to robust transceiver designs for a af mimo relaying system .    * definition 1 : * a matrix - monotone function is defined as @xmath64 which maps a matrix variable from a subsect of positive semi - definite matrices to a real number . if @xmath64 is a monotonically decreasing matrix - monotone function on positive semi - definite matrices , it satisfies @xmath65on the other hand , when @xmath64 a monotonically increasing matrix - monotone , it means @xmath66 is a monotonically decreasing matrix - monotone function @xcite .    in the following we focus our attention on a kind of optimization problems with a decreasing matrix - monotone function as objective , which is formulated as @xmath67 solving the optimization problem ( [ opt_x ] ) , two important solutions are derived and are the basis for the following derivations .    * * the optimal solution of ( [ opt_x ] ) satisfies @xmath68/\eta_x = p.\end{aligned}\]]defining the unitary matrix @xmath69 and rectangular diagonal matrix @xmath70 based on the following singular value decomposition @xmath71the optimal solution of the optimization problem ( [ opt_x ] ) has following structure @xmath72\},\end{aligned}\]]where @xmath73 is an unitary matrix and @xmath74 is a rectangular diagonal matrix with real diagonal elements .    _ proof : _ see appendix  [ appedix:3 ] . @xmath75    * * for a complex matrix @xmath76 and a positive define matrix @xmath77 , based on following eigen - decomposition @xmath78when the objective function is @xmath79 the unitary matrix @xmath73 in ( [ conslusion_1 ] ) equals to @xmath80 defining @xmath81@xmath82 , @xmath83 in ( [ conslusion_1 ] ) has the structure of @xmath84,\end{aligned}\]]where @xmath85 is a @xmath86 diagonal matrix .    _ proof : _ see appendix  [ appedix:4 ] . @xmath75      in this section the structure of optimal @xmath87 is derived . based on the matrix inversion lemma , the mse matrix in ( [ mse_matrix ] ) can be rewritten as @xmath88 based on which for any given @xmath18 the optimization problem with respect to @xmath87 becomes as @xmath89    defining unitary matrices @xmath90 and @xmath91 based on the following singular value decomposition @xmath92 we have the following eigen - decomposition @xmath93together with the following singular value decomposition @xmath94and based on * conclusions 1 * and * * 2 * * ] , the optimal @xmath87 has the following structure @xmath95\},\end{aligned}\]]where @xmath96 is a @xmath97 diagonal matrix . meanwhile , for the optimal @xmath87 the following constraint is fulfilled @xmath98/\eta_f = p_r.\end{aligned}\ ] ]      in the following , it will be proved that given the structure of @xmath87 the optimization problem for @xmath18 is the same as that for @xmath87 . using the optimal structure of @xmath87 , we have @xmath99 using the following substitution @xmath100and the matrix inversion lemma again , the mse matrix can be reformulated as    @xmath101    therefore , the optimization problem with respective to @xmath18 is equivalent to @xmath102based on the definitions of @xmath103 and @xmath104 , it can be derived that @xmath105 has the following eigen - decomposition @xmath106 together with following singular value decomposition , @xmath107and using * conclusions 1 * and * * 2**. ] , the optimal @xmath18 has the following structure @xmath108\},\end{aligned}\]]where @xmath109 a @xmath97 diagonal matrix . considering that as there are no constraints on @xmath110 , @xmath110 can be an arbitrary @xmath111 unitary matrix . based on * conclusion 1 * , it can be concluded that the optimal @xmath18 satisfies @xmath112/\eta_p = p_s.\end{aligned}\ ] ]    substituting ( [ structure_p_capacity ] ) into ( [ u_1 ] ) , it can be derived that @xmath113 and then the optimal structure of @xmath87 is @xmath114 based on the optimal structure given by ( [ structure_p_capacity ] ) and ( [ structure_f_capacity ] ) and with regard to the fact that @xmath115 and @xmath116 are determined by @xmath109 and @xmath96 , respectively , the left unknown variables are only @xmath109 and @xmath96 . based on ( [ cons_f ] ) and ( [ cons_p ] ) , the optimization problem ( [ opt : capacity ] ) also equals to @xmath119/\eta_f= p_r \nonumber \\   & \ \ \ \ \ \ \ \ { \rm{tr}}[{\bf{p}}{\bf{p}}^{\rm{h}}(\alpha_1p_s{\boldsymbol \psi}_{sr}+\sigma_{n_1}^2{\bf{i}})]/\eta_p = p_s.\end{aligned}\]]furthermore , with the following diagonal matrices @xmath120and substituting ( [ structure_p_capacity ] ) and ( [ structure_f_capacity ] ) into ( [ opt : capacity_a ] ) , the optimization problem ( [ opt : capacity_a ] ) can be rewritten as @xmath121with respective to the fact the problem ( [ opt : capacity_b ] ) is inherently non - convex and difficult to solve , an iterative water - filling solution is proposed in this paper . when @xmath122 s are fixed , @xmath123 s can be computed as    @xmath124    where @xmath125 is the lagrange multiplier which makes @xmath126 . on the other hand , when @xmath123 s are fixed @xmath122 s can be computed as    @xmath127    where @xmath128 is the lagrange multiplier which makes @xmath129 .    * special cases * : several existing algorithms can be considered as special cases of our proposed solution . @xmath130 when csi is perfectly known and @xmath131 , the proposed solution for @xmath25 reduces to that in @xcite . @xmath130 when csi is perfectly known , the proposed solution for @xmath18 and @xmath25 reduces to that given in @xcite . @xmath130 when the second hop channel is an identity matrix and noiseless , the proposed solution for source precoder design reduces to that given in @xcite . in this section , simulation results are presented to assess the performance of the proposed algorithm and for the purpose of comparison , the algorithm based on the estimated channel only ( without taking the channel errors into account ) @xcite . in the following , we consider an af mimo relay system where the source , relay and destination are equipped with same number of antennas , i.e. , @xmath132 . the channels @xmath17 and @xmath28 are randomly generated according to i.i.d . gaussian distribution .    to estimate the channels , a practical lmmse estimation algorithm is adopted @xcite . for the training sequence @xmath133 , a famous exponential correlation matrix is used to describe the correlation matrix of @xmath133 , i.e. , @xmath134 where @xmath135_{ij}=\rho^{|i - j|}$ ] . as a result , @xmath136 and @xmath137 where @xmath138 is the signal - to - noise ratio ( snr ) in channel estimation process @xcite .    in the simulation ,    . ]    for data transmission process the snr at relay is defined as @xmath139 , and the snr at destination is defined as @xmath140 . for simplicity , it is also assumed that @xmath141 . each point in the following figure is an average of 10000 independent channel realizations . [ fig:1 ] shows the sum rates of different algorithms including the proposed robust design and its counterpart based on estimated csi only when @xmath142 . it can be seen that the performance of the proposed robust design is always better than that of the design based estimated on csi only . furthermore , as the channel estimation snr decreases the performance gain of the robust design becomes larger . robust mutual information maximization transceiver design for dual - hop af mimo relay systems was investigated . with gaussian distributed channel errors the precoder at the source and forwarding matrix at the relay were jointly designed . the structures of the optimal solutions were derived first , which differentiates our work from the existing works . for the unknown diagonal matrices , a well - known iterative water - filling solution was proposed . the simulation result demonstrated the performance advantage of our robust design . denoting @xmath54 and @xmath55 as the received pilots in the separate dual hop channel estimations , the capacity between the source and destination equals to @xmath143where @xmath144 denotes the conditional entropy of @xmath19 when @xmath54 and @xmath55 are known and @xmath145 is the conditional entropy of @xmath19 when @xmath26 , @xmath54 and @xmath55 are known . as @xmath54 and @xmath55 are independent with @xmath19 , @xmath146 . the second term on the right hand side of ( [ capacity_equ ] ) will be discussed in the following .    denoting @xmath147 as the conditional mean and based on the definition of covariance matrix , the conditional covariance matrix satisfies @xmath148where @xmath149 is an arbitrary estimate of @xmath19 including lmmse estimate . it should be pointed out that as csi is not perfectly known , @xmath150 . notice that @xmath151 is the conditional mean and thus we have @xmath152based on which the right hand side of ( [ covariance_matrix ] ) also equals to @xmath153 substituting ( [ covariance_matrix_a ] ) into ( [ covariance_matrix ] ) , we have @xmath154it is also well - known that with fixed covariance matrix , gaussian distribution has the maximum entropy . therefore , it is concluded that @xmath155based on which , an lower bound of the capacity ( [ capacity_equ ] ) is @xmath156the final equality comes from the fact @xmath157 . as @xmath64 is a matrix monotonically decreasing function , it can be proven that for the optimal solution the power constraint is always active , i.e. , @xmath158 @xcite . as a result , we have the following relationship . @xmath159from ( [ eta_x ] ) , the constraint of the optimization problem ( [ opt_x ] ) equals to @xmath160/\eta_{x}=p,\end{aligned}\]]based on which the optimization problem ( [ opt_x ] ) is equivalent to @xmath161/\eta_{x}=p.\end{aligned}\]]then , defining a new variable @xmath162the optimization problem ( [ opt : app_3_1 ] ) is further reformulated as @xmath163    for any given @xmath164 , based on the following singular decompositions @xmath165 there exists a matrix @xmath166 satisfying @xmath167where @xmath74 is a diagonal matrix with the same rank as @xmath168 and @xmath169 is a scalar which makes @xmath170 hold . based on lemma 12 in @xcite , the following inequality holds @xmath171together with the fact that @xmath64 is a matrix monotonically decreasing function , the following inequality holds @xmath172therefore , it is concluded that the optimal @xmath164 has the structure given by ( [ x_opt_str ] ) . furthermore , based on the definition of @xmath164 ( [ app_equ_x ] ) , the optimal @xmath173 has the following structure @xmath174where @xmath74 is a diagonal matrix . in ( [ x_structure ] ) @xmath175 is unknown either . in order to solve @xmath175 , substitute the structure of @xmath173 in ( [ x_structure ] ) into the definition of @xmath175 in ( [ opt_x ] ) , and then we get the following equation @xmath176+\sigma_{n}^2.\end{aligned}\]]this is a simple linear function of @xmath177 , and @xmath175 can be easily solved to be @xmath178\}.\end{aligned}\ ] ] the objective function in ( [ objective ] ) can be reformulated as @xmath179where the second equality is based on the fact that @xmath180 . using the matrix inequality that for two positive semi - definite matrices @xmath181 and @xmath182 i.e. , @xmath183 @xcite , we directly have @xmath184\nonumber \\ & -\sum_i{\rm{log}}[\lambda_i({\bf{x}}^{\rm{h}}{\bf{h}}^{\rm{h } } { \bf{h}}{\bf{x}}/\eta_x+{\bf{i}})].\end{aligned}\]]together with the optimal structure given by * conclusion 1 * , in order to make the equality in ( [ case_b ] ) hold the following equation holds @xmath185 in light of the fact that power is never loaded to the eigenchannels with zero magnitudes @xcite , the diagonal matrix @xmath83 has the following structure @xmath84,\end{aligned}\]]where @xmath85 is a @xmath86 diagonal matrix and @xmath186 .      o. munoz - medina , j. vidal , and a. agustin , `` linear transceiver design in nonregenerative relays with channel state information , '' _ ieee trans . signal process . 6 , pp . 25932604 , june 2007 . y. rong , x. tang , and y. hua , `` a unified framework for optimizing linear nonregenerative multicarrier mimo relay communication systems , '' _ ieee trans . signal process . 48374851 , dec . 2009 .        c. xing , s. ma , y .- c . wu , and t .- s . ng , `` transceiver design for dual - hop non - regenerative mimo - ofdm relay systems under channel uncertainties , '' _ ieee trans . signal process . _ , vol.58 , no.12 , pp.6325 - 6339 , dec . 2010 .    c. xing , s. ma , and y .- c wu , `` robust joint design of linear relay precoder and destination equalizer for dual - hop amplify - and - forward mimo relay systems , '' _ ieee trans . signal process . 4 , pp . 22732283 , apr . 2010 . d. p. palomar , j. m. cioffi , and m. a. lagunas , `` joint robust tx - rx beamforming design for multicarrier mimo channels : a unified framework for convex optimization , '' _ ieee trans . on signal processing _ , 9 , pp . 23812399 , sep . 2003 .
Abstract:

Answer: in this paper , we investigate the robust transceiver design for dual - hop amplify - and - forward ( af ) mimo relay systems with gaussian distributed channel estimation errors . aiming at maximizing the mutual information under imperfect channel state information ( csi ) , source precoder at source and forwarding matrix at the relay are jointly optimized . using some elegant attributes of matrix - monotone functions , the structures of the optimal solutions are derived first . then based on the derived structure an iterative waterfilling solution is proposed . several existing algorithms are shown to be special cases of the proposed solution . finally , the effectiveness of the proposed robust design is demonstrated by simulation results .


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 of the most important questions in particle physics is : where is the higgs boson ? the lep exclusion limits , @xmath0 gev , constraints from precision electroweak data , @xmath1 gev  @xcite , and recently also the tevatron limits , leave about 40 gev window for the standard model ( sm ) higgs boson . this 40 gev window is very interesting and there are several suggestive hints or coincidences related to it . first of all , this window overlaps with the range of higgs masses in which the standard model can be a consistent theory all the way to the grand unification scale or the planck scale , @xmath2 gev . another interesting coincidence is that this window also overlaps with the range of higgs masses predicted in the minimal supersymmetric model ( mssm ) , @xmath3 gev . hence , it is expected that the higgs boson is somewhere in this window and most of the effort is focused on discovery strategies related to this possibility . however there are also several compelling hints that the sm - like higgs boson is below the lep exclusion limits . first of all , electroweak symmetry breaking ( ewsb ) in simple supersymmetric ( susy ) extensions of the standard model , with superpartners near the electroweak ( ew ) scale , generically predicts the higgs boson not heavier than about 100 gev . non - observation of the higgs boson at lep resulted in the  fine - tuning " problem in these models  @xcite . second of all , the best fit to precision electroweak data is achieved for the higgs mass of 87 gev  @xcite . the third hint comes from the lep data : the largest excess of higgs like events at lep corresponds to the higgs mass of 98 gev . it is an interesting coincidence that natural ewsb in susy models , precision ew data , and the largest excess of higgs like events point to the same region . this supports the idea that the higgs boson is light , somewhere near 100 gev , and we missed it at lep . how can this be ? the basic idea is very simple : if the sm - like higgs boson decays in a different way than the higgs boson in the standard model then the usual experimental limits do not apply . such a higgs can be light , as predicted from susy , gives better agreement with precision electroweak data , and can even explain the largest excess of higgs like events at lep  @xcite . in theories beyond the sm the higgs sector is usually more complicated and there are typically many other higgses in addition to the sm - like higgs boson . for example , there are five higgs bosons in the mssm : light and heavy cp even higgses , @xmath4 and @xmath5 , the cp odd higgs , @xmath6 , and a pair of charged higgs bosons , @xmath7 ; seven in the next - to - minimal supersymmetric model ( nmssm ) , which contains an additional singlet field : three cp even higgs bosons , @xmath8 , two cp odd higgs bosons , @xmath9 , and a pair of charged higgs bosons ; and there are many simple models with even more complicated higgs sectors . usually we explore parameter space in which the extra higgses are somewhat heavy  the so called decoupling limit . the decoupling limit is not the only possibility . one of the extra higgses can be light , for example , the singlet cp odd higgs in the nmssm . if it is sufficiently light , the sm - like higgs boson can ( and typically would ) dominantly decay into a pair of cp odd higgses and eventually , depending on the mass of the cp odd higgs boson , into * four * b quarks , * four * @xmath10 leptons , * four * c quarks , * four * @xmath11 leptons , * four * electrons , * four * light quarks , or gluons  @xcite . limits on these 4-body final states are weaker than limits on the sm higgs boson ( decaying into @xmath12 ) and currently @xmath13 leptons , c quarks , and light quarks or gluons final states allow the sm - like higgs boson at @xmath14 gev or even lighter  @xcite , as is predicted from the best theoretically motivated region of the parameter space in supersymmetric theories , and it also gives much better agreement with precision electroweak data . in addition , the subleading decay mode of the higgs boson , @xmath15 , with branching ratio of @xmath16 can completely explain the largest excess ( @xmath17 ) of higgs - like events at lep in the @xmath18 final state ( for a reconstructed mass @xmath19 gev )  @xcite . another , perhaps even more interesting variation of the above nmssm scenario is the scenario with a * doublet - like * cp odd higgs bellow the @xmath12 threshold . for small @xmath20 , @xmath21 , this scenario is the least constrained ( and only marginally ruled out ) in the mssm , and thus easily viable in simple extensions of the mssm  @xcite . surprisingly , the prediction from this region is that all the higgses resulting from two higgs doublets : @xmath4 , @xmath5 , @xmath6 and @xmath7 could have been produced already at lep or the tevatron , but would have escaped detection because they decay in modes that have not been searched for or the experiments are not sensitive to . the heavy cp even and the cp odd higgses could have been produced at lep in @xmath22 but they would avoid detection because the dominant decay mode of @xmath5 , @xmath23 , has not been searched for . the charged higgs is also very little constrained although it could have been pair produced at lep or appeared in decays of top quarks produced at the tevatron . the dominant decay mode of the charged higgs in this scenario is @xmath24 with @xmath25 or @xmath26 . in addition , the charged higgs with properties emerging in this scenario and the mass close to the mass of the @xmath27 boson could explain the @xmath28 deviation from lepton universality in @xmath27 decays measured at lep  @xcite as pointed out in ref . @xcite . the main idea is simple , and it can be realized in a variety of other models . for example , in specific little higgs models the sm - like higgs boson can dominantly decay into four c quarks  @xcite or four gluons  @xcite . four body final states of the higgs boson can also occur in composite higgs models  @xcite , and models for dark matter  @xcite among others . for a review of other scenarios and references , see ref . @xcite . the situation can be even more complicated if there are several lighter higgses . in this case the sm - like higgs boson can cascade decay into the lightest one , and eventually , depending on the mass of the lightest higgs boson , the sm - like higgs would decay into a large number of light jets , or a large number of muons , or a large number of electrons , or large number of soft photons  @xcite . such events would be quite spectacular . some of these signatures were recently studied in ref . @xcite . additional level of complexity arises when more higgs boson share the coupling to the z boson and there is not a single sm - like higgs boson , see _ e.g. _  @xcite . the simplest models allowing the sm - like higgs boson at @xmath14 gev are still those with 4-body decay modes . the nmssm scenario with @xmath29 has been recently studied and searched for at a variety of experiments . in the rest of this talk we will focus on this scenario . the first constraints on these scenarios came from cleo  @xcite . the light cp odd higgs can be produced in upsilon decays , @xmath30 , and the predicted branching ratio from the nmssm model typically varies between @xmath31 for the cp odd higgs mass close to @xmath32 threshold , and @xmath33 for the cp odd higgs mass close to the upsilon mass  @xcite . at present , the strongest constraints come from babar  @xcite that sets limits @xmath34 ( slightly varying around this value for different masses of the cp odd higgs boson ) . in order to satisfy these limits it is typically required that @xmath35 gev for @xmath36 . for smaller @xmath20 these limits become weaker as @xmath26 channel becomes more important . more details about the impact of these searches on the nmssm parameter space can be found in ref . @xcite . searches for @xmath37 were performed and are in progress at the tevatron  @xcite . these searches are not sensitive yet to the sm - like higgs boson at 100 gev . nevertheless it is interesting to note that it is advantageous to search for a subleading decay mode of one of the cp odd higgs bosons , @xmath38 , as suggested in ref . the rate for @xmath39 is suppressed by a factor of @xmath40 compared to @xmath41 , but these events are much cleaner and one can reconstruct the mass of the cp odd higgs boson . this search mode is very promising at the lhc where one expects about 500 events in 1 fb@xmath42 of data ( for 14 tev center of mass energy )  @xcite . the most important constraints on this scenario come from the recent search at aleph  @xcite . naively this search rules out the sm - like higgs that decays into four @xmath10 leptons up to 107 gev . note however , that this search places limits on @xmath43 and thus it is very sensitive to @xmath44 , which in the nmssm is never equal ( or even close ) to 100% . for large @xmath20 , depending on the mass of the cp odd higgs boson , @xmath44 is between 0.9 and 0.7 because of the sizable branching ratio of @xmath45 . in addition , for smaller @xmath20 a new decay mode , @xmath46 , becomes important . folding in realistic branching ratios of @xmath47 one finds that aleph limits allow the sm - like higgs boson at 100 gev ( or slightly lighter ) for any @xmath36 only when @xmath48 gev , and generically for @xmath49  @xcite . recently we have been working on new ways to search for these scenarios . one question was whether we can directly observe the light cp odd higgs ( without relying on producing heavier higgses that decay into the cp odd higgs ) . the direct production cross section of the light cp odd higgs boson in gluon fusion channel is , depending on the mass of the cp odd higgs boson and @xmath20 , between 1 and 100 nb at the tevatron , and between 10 and 1000 nb at the lhc  @xcite . in spite of the large production cross section , picking up the signal on huge background is a serious problem . if one looks for the dominant decay mode , @xmath50 , then it is basically hopeless , however searching for the subleading decay mode , @xmath38 , is actually very promising . we still have to deal with huge background from drell - yan production of @xmath11-pairs , and semileptonic b and c decays . nevertheless , based on existing analyses we showed that current data sets would allow cdf and d0 to improve on limits from babar , especially for masses of the cp odd higgs boson close to the mass of the upsilon  @xcite . at the lhc this is a very promising search mode and already the first fb@xmath42 of data might provide an evidence for , or the discovery of the light cp odd higgs . more details , and predictions for the integrated luminosity needed for the discovery as a function of the mass and couplings of the light cp odd higgs boson can be found in ref . @xcite .    looking for subleading decay modes of the light cp odd higgs might be also advantageous in searches for the charged higgs . if the cp odd higgs boson has a significant dublet component than the charged higgs is generically light , typically lighter than the top quark . depending on the mass of the charged higgs and @xmath20 the @xmath51 can go up to 40% for @xmath52 and @xmath53 gev dropping very fast with increasing @xmath20 and increasing the mass of the charged higgs . the dominant decay mode of the charged higgs is @xmath54 and it was not previously searched for ( till recently ) . the cdf recently performed the search and set the upper limit on @xmath55 to about 10%  @xcite constraining only a small region of the parameter space . the search for the charged higgs decaying into @xmath56 will be relatively easy at the lhc which is the top quark factory . we might again look for a subleading decay mode @xmath38 that we can not afford at the tevatron because the rate would be too small . in addition , we can look for events in which one of the @xmath27 bosons decays into @xmath57 , resulting in 3-muon events with properties that easily stand out from the background . we expect about 30 events of this type with 1 fb@xmath42 of data at the lhc  @xcite .    in conclusion , four body decay modes : @xmath13 , @xmath58 , @xmath59 , and @xmath60 are the simplest possibilities allowing the sm - like higgs boson at @xmath14 gev . for the @xmath61 scenario in the nmssm , searching for dominant decay modes typically requires many tens of fb@xmath42 of data . however , searching for the subleading decays modes is very promising with early data at the lhc . especially searches for * @xmath62  @xcite , * @xmath63  @xcite , * @xmath64  @xcite    could lead to an evidence or discovery already with 1 fb@xmath42 of data . i would like to thank j. gunion , h.d . kim , e. lunghi and a. raval for collaboration on projects this talk is based on . b.  bellazzini , c.  csaki , a.  falkowski and a.  weiler , arxiv:0906.3026 [ hep - ph ] . b.  gripaios , a.  pomarol , f.  riva and j.  serra , jhep * 0904 * , 070 ( 2009 ) [ arxiv:0902.1483 [ hep - ph ] ] . j.  mardon , y.  nomura and j.  thaler , arxiv:0905.3749 [ hep - ph ] . w.  love _ et al . _ [ cleo collaboration ] , phys . * 101 * , 151802 ( 2008 ) [ arxiv:0807.1427 [ hep - ex ] ] . r.  dermisek , j.  f.  gunion and b.  mcelrath , phys . d * 76 * , 051105 ( 2007 ) [ arxiv : hep - ph/0612031 ] .
Abstract:

Answer: we review the motivation for hidden higgs scenarios and discuss the light cp odd higgs scenario in the nmssm as an example . we summarize experimental constraints including recent limits from babar and aleph . the main part of the talk is the discussion of dominant decay modes of the standard model like higgs boson , and related decay modes of the charged higgs and heavy cp even higgs bosons , in these scenarios with the focus on signatures and prospects for the lhc . examples include the direct production of a light cp odd higgs boson , and a light charged higgs boson in top quark decays .


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 crab pulsar wind nebula is powered by its central pulsar born from a supernova explosion in 1054 . it is a very luminous source in almost all wavelengths , from radio to the very high energy ( vhe ) @xmath0-ray bands . the broadband non - thermal emission spectrum is well modeled by a synchrotron - inverse compton ( ic ) scenario with the transition from the synchrotron to ic component occurring at a few hundred mev @xcite . the overall emission seems to be steady , so that the crab nebula has been adopted as a standard candle in high energy astrophysics to calibrate observations from different instruments .    however , detailed images of the crab nebula in optical and x - ray bands indicated dynamical structures at small scales . observations by the hubble space telescope revealed wisps and knots in the nebula , and resolved some substructures of the highly variable wisps @xcite . x - ray observations by rosat and chandra uncovered a jet - torus structure of the inner nebula @xcite , and the equatorial ring is found moving outwards with a speed of @xmath3 @xcite . at higher energies in the @xmath0-ray band , detailed structures of the nebula generally can not be resolved . however , comptel and egret observations indicated that the synchrotron radiation in @xmath4 mev is variable on a timescale of @xmath5 yr .    in september 2010 , agile collaboration reported a @xmath0-ray flare above @xmath6 mev from the direction of the crab nebula , which lasted for about 3 days ( atel # 2855 ; tavani et al . the flux during the flare period is about @xmath7 times higher than the average one . this flare was soon confirmed by the fermi / lat collaboration ( atel # 2861 ) , and archive search of the fermi / lat data revealed another flare in february 2009 @xcite . there were many other simultaneous or follow - up measurements for the flare in september 2010 in x - ray , optical , infrared and radio bands . however , no significant flux enhancement was discovered . in the vhe @xmath0-ray energies , argo - ybj collaboration claimed the detection of a flux enhancement around tev , with a possibly longer duration ( atel # 2921 ) . however , magic and veritas observed the source with a shorter duration during the flare phase and did not reveal any enhancement in flux ( atel # 2967 , # 2968 ) . it was proposed that the flares was due to synchrotron emission of ultra - relativistic electrons with energies up to @xmath2pev @xcite . considering the fact that the rest frame synchrotron radiation in a magnetic field dominated acceleration regime can not exceed @xmath8 mev due to the fast synchrotron cooling of high energy electrons @xcite , the observed gev emission implies doppler shift of the radiation region or other acceleration mechanisms instead of shock acceleration @xcite . @xcite proposed that the @xmath0-ray variability ( flare ) originated from the `` inner knot '' of the crab nebula ( `` knot 1 '' as defined in @xcite ) , with mildly doppler - boosted emission . the instability of the termination shock may cause the variability as revealed by magnetohydrodynamic ( mhd ) simulations . @xcite suggested that electrons are accelerated in a region behind the shock , and the variability was attributed to changes in the maximum energy of accelerated electrons , electron spectral index or the magnetic field .    in this letter we employ a statistical approach to model the @xmath0-ray variability of the nebula . fermi / lat observations have shown that the low - energy synchrotron component is variable on monthly time scale , while the high - energy ic component seems to be stable @xcite . since the synchrotron @xmath0-rays are produced by the highest energy electrons , these observations indicate that fluctuations at the high - energy end of the electron distribution might be responsible for the variability . it is natural to expect that events that can generate the highest energy electrons are rarer , and therefore would lead to the largest fluctuation . therefore the variability and flares in the sub - gev @xmath0-rays can be simply due to the statistical fluctuation of the highest energy electrons achievable in the electron accelerators . lower energy electrons do not suffer from significant fluctuations since many more accelerators can contribute to them simultaneously , which gives rise to a `` steady - state '' emission in both the lower energy synchrotron component and the higher energy ic component . we build a physically possible model to investigate the above picture in more details . we assume that electrons are accelerated in a series of knots with a power - law energy spectrum @xmath9 with the maximum electron energy @xmath10 being proportional to the size of the @xmath11th knotpc , which just corresponds to the size of the emission region as indicated by the day - scale flare . therefore we assume the largest size of the knots is @xmath12 , and for all other knots smaller than @xmath12 the scaling relation between @xmath10 and @xmath13 can well hold . ] . we further assume that the size of knots has a power - law distribution @xmath14 . adding the contribution from all these knots together , noticing that the normalization of the synchrotron spectrum of the knots scales as @xmath15 , we get the total electron spectrum as @xmath16 where @xmath17 is the heaviside step function , @xmath18 is the maximum energy of the largest knot . if @xmath19 , the above integral of the last step is approximately proportinal to @xmath20 given the maximum energy of the largest knot @xmath21 ; if @xmath22 the integral is nearly a constant determined by @xmath18 . therefore the electron spectrum is a power law @xmath23 with power - law index @xmath24 and a cutoff determined by the largest knot . the maximal energy of synchrotron emission in the magnetic field dominated acceleration regime is @xmath8 mev @xcite ( defined by equating the synchrotron cooling time scale to the gyro - period , the minimum of the acceleration time scale ) . for smaller accelerators , the maximum energy is also limited by the size of the accelerator . in order to account for the observations in the @xmath1 mev range , we therefore further employ a mild doppler boost factor . we assume a gaussian distributed lorentz factor of all the knots with a mean value @xmath25 and a standard deviation @xmath26 . such a lorentz factor is consistent with the upper limit of the typical velocity of the jet @xcite . we assume that the angle @xmath27 between the knot motion and the line - of - sight is randomly distributed . the doppler factor @xmath28 ( where @xmath29 is the velocity of the knot in unit of @xmath30 ) therefore is distributed within the range @xmath31 for a @xmath32 distribution of @xmath33 . only knots with large enough size and large enough doppler factor can contribute to high energy synchrotron radiation to account for the flares observed by the fermi / lat . we realize such a picture by a monte - carlo simulation . we directly simulate the synchrotron spectrum in each knot instead of the electron spectrum . the synchrotron spectrum of each knot in its co - moving frame is adopted as @xmath34 , where @xmath35 , @xmath36 is the synchrotron spectral index . where @xmath37 is the power - law index of the overall electron distribution . the magnetic field in each knot is not explicitly used in the calculation . the magnetic field will affect the cooling break of the low energy synchrotron spectrum as well as the high energy cooling time scale . to keep the basic framework of the present study , we may expect the magnetic field not to vary significantly with respect to the knot size .    the energy loss of electrons in the knots can be important and needs to be considered . for instantaneous injection of accelerated electrons , the probability of observing electrons at a given energy is proportional to the corresponding synchrotron cooling time @xmath38 . for low energy electrons with cooling time longer than the age of the crab nebula @xmath39 yr , these electrons will survive and accumulate . therefore we can construct a probability function @xmath40 , describing the effect of cooling on the electron spectrum . the cooling time of electrons producing synchrotron photon with energy @xmath41 is @xmath42 yr . for an average magnetic field @xmath43 mg , a mildly doppler factor @xmath44 , electrons producing photons with an energy @xmath45 ev have a cooling time comparable to @xmath46 . the observed spectrum of crab nebula has a break near @xmath47 ev ( see below fig . [ fig : f1 ] ) .    in this work we mainly focus on the high energy part of the spectra , e.g. , from x - ray to @xmath0-ray band , then @xmath48 . assuming that the accelerated electrons have a spectrum index @xmath49 , the electron spectrum index in the synchrotron cooling dominated regime is then @xmath50 . here we adopt @xmath51 which corresponds to @xmath52 , as required by fitting the radio - optical data . the maximum value of @xmath53 is taken as @xmath54 mev in the co - moving frame ( corresponding to the largest knot ) . this value is derived according to a fit to the average fermi / lat spectrum of the synchrotron component @xcite , given the lorentz factor distribution . the power - law index of the size distribution is adopted as @xmath55 . adding all the contribution from the knots together we can get a total synchrotron spectrum @xmath56 , which can reproduce the observed optical - mev band data of the crab nebula ( fig . [ fig : f1 ] ) . note that these parameters , i.e. , @xmath33 , @xmath53 , @xmath57 and @xmath58 are not uniquely determined . what we adopt in this work is an illustration of the model . once the data from optical to @xmath0-ray band are reproduced , the discussion below is basically unchanged if a different set of parameters is taken . for each knot we randomly generate an angle @xmath27 and a lorentz factor @xmath33 based on the assumed gaussian distribution , and calculate the doppler factor @xmath59 . then the frequency of synchrotron emission from each knot is shifted from @xmath60 to @xmath61 , and the flux is enhanced to @xmath62 . the contribution from all the knots are summed up to obtain the total emission spectrum .     realizations . the thick blue line denote the fit to the multi - wavelength steady state observational data , which can be understood as the average result of many realizations in our model . the black - dashed curve denotes the contribution from a maximum size knot with slightly specific parameters ( see the text ) , which is shown to be able to explain the large flare in september 2010 . the references of the observational data are : radio - optical @xcite , integral @xcite , comptel , fermi / lat @xcite , magic @xcite , hegra @xcite , hess , tibet - as@xmath0 @xcite , and fermi / lat flares @xcite . ]    a series of simulated synchrotron spectra are shown in fig . [ fig : f1 ] , which are compared with the observational data of the crab nebula . note that for the integral data we adopt a flux normalization factor @xmath63 in order to be consistent with the extrapolation of comptel data . in a @xmath64 factor was adopted to match the xmm - newton data . the difference of flux normalization may be due to discrepancies in calibration of different detectors . the fact that the integral data may contain emission from the crab pulsar may also contribute to the difference . for comparison we show the fitting results to the broad band data with the blue thick line , with a broken power - law electron spectrum . this result is similar to the model invoking two population electrons as introduced in @xcite and . it is shown that the simulated synchrotron spectra have very small fluctuations in the x - ray energy band , but fluctuate significantly at higher energies ( @xmath65 mev ) . the average of the simulated @xmath0-ray flux would be consistent with the long term fermi / lat and comptel data . the fluctuations can be responsible for the observed variabilities and flares . the very large flare at september 2010 that is not well reproduced by these simulations can be regarded as an event with very small probability . such an event may be due to a large knot with specific parameters different from what used in the model . the black - dashed line in fig . [ fig : f1 ] shows a possible reproduction of the large flare , produced by a knot with maximum size , doppler factor @xmath66 , cutoff energy of synchrotron radiation in its comoving system @xmath67 mev and an additional normalization factor @xmath68 of the flux . detailed modeling of the knot formation and its evolution is expected to better address this issue . we pick out two realizations and show the synchrotron sky - maps at 10 kev ( left ) and 100 mev ( right ) in fig . [ fig : f2 ] , respectively . to generate a sky - map , we designate an @xmath69 plane with @xmath70 pixels . we randomly assign a position for each knot in the @xmath69 plane , with the largest knot having a size of @xmath71 pixels . we then sum up contributions of all the knots located in each pixel and obtain the sky - map . it is shown that in the x - ray band ( left panels ) the skymaps are relatively smooth with some hot spots due to several large knots . these hot spots may correspond to the chandra and hubble knots of the crab nebula ( atel # 2882 , # 2903 ) . the total flux of x - ray emission , however , is statistically steady from time to time . the width of the total x - ray flux distribution at @xmath72 kev is estimated to be of the order @xmath73 , which is consistent with the upper limit of @xmath74 obtained by swift / bat ( atel # 2858 ) . the @xmath1 mev @xmath0-ray sky - map could differ significantly from each other . in some cases , one or several large knots would dominate the total flux . the predicted energy - dependent spatial scales due to fluctuations may not be ready to test with the available high energy telescopes , but higher angular resolution observations in the future will help . -ray flare , while the lower one shows a realization with relatively weak @xmath0-ray emission . ]    according to our model , the ic component should not show significant variability except in the high energy regime ( @xmath75 tev ) , depending on the maximal electron energies . for example , for a magnetic field @xmath76 mg , the maximal lorentz factor of electron can be as high as @xmath77 . therefore only the ic photons with energies higher than several hundred tev would have relatively large fluctuations . this conclusion is different from that of @xcite , who suggested that the ic @xmath0-ray spectra also vary significantly at multi - tev energies . the upper panel shows the results with weekly bin and the bottom panel for monthly bin . the vertical axes in each panel indicate the average value of the simulations , which is normalized to the observational value @xmath78 @xmath79 s@xmath80 @xcite . the fermi / lat observational data with monthly bin are also plotted in the lower panel . ]        to investigate fluctuations of emission in detail , we show in fig . [ fig : f3 ] the light curves of the simulated synchrotron fluxes in the @xmath1 mev energy band . each bin of the light curves represents an independent realization , and the length of the time bin is proportional to the number of knots . we normalize the time scale so that over a year there is an enhancement in flux by at least a factor of 5 - 6 , which may be responsible for a flare , and plot the weekly and monthly bin lightcurves in the upper and lower panels of fig . [ fig : f3 ] , respectively . for comparison we also plot the fermi / lat observed monthly light curve @xcite in the lower panel of fig . [ fig : f3 ] . it can be seen that the predicted fluctuation level of the simulated results is very similar to that observed by fermi / lat . in order to show this more clearly , we plot the histogram of the monthly bin flux distribution for the simulated results along with that of fermi / lat observation in the left panel of fig . [ fig : f4 ] . these two distributions are indeed similar . in the right panel of fig . [ fig : f4 ] the histogram of flux distribution for weekly bin is also plotted . a concentration toward low - flux is evident , which is distinct from a symmetric distribution around the average value . this is because in the short time scale limit , the distribution of the gamma - ray fluxes is determined by the distributions of the knots ( power - law ) and the doppler factors . only the knots with both large enough sizes and high enough doppler factors can produce large gamma - ray fluxes . therefore the probability to give large fluxes is very small , and the flux distribution concentrates towards low - flux values . in the limit of long time scale , the distribution should be gaussian as a result of central limit theorem . such a prediction can be tested by a detailed analysis of the fermi / lat data . in summary , we propose a statistical picture to explain the observed @xmath0-ray variability and flares of the crab nebula , using the fluctuations of the highest end of the electron spectra . the electrons are thought to be accelerated in a series of knots , with a size distribution @xmath14 and a distribution of the doppler factor . the maximal energy of the electrons in the co - moving frame is assumed to be proportional to the size of the knots . thus the rare knots with large sizes and high doppler boosts may generate electrons up to @xmath2pev , and hence , be responsible for the observed @xmath0-rays flares at fermi / lat . on the other hand , the low energy electrons are generated by many smaller knots . the average effect smooths the fluctuation , so that the low energy synchrotron component and the ic component do not change significantly . this scenario can naturally explain both the variability in the mev - gev band , and the relatively steady emission in lower ( optical , x - ray ) and higher ( tev ) energies . the expected variability of the monthly bin fluxes above 100 mev are well consistent with that observed by fermi / lat . the two large @xmath0-ray flares can also be naturally accounted for without additional assumptions .    in the simulation , we do not consider the detailed structure ( e.g. jet , torus ) of the nebula . the knots are assumed to have a uniform and isotropic distribution . considering a more realistic morphology of the nebula would not affect our conclusion noticeably , as long as the distributions of the knot sizes and doppler factors are similar to those introduced here . in more general terms , these knots may be interpreted as individual electron acceleration events distributed throughout the nebula . we point out that at low energies contributions from large amount of knots are equivalent with previous studies that electrons are accelerated continuously inside the whole nebula . although the model does not discuss the absolute time scales of flares and @xmath0-ray fluctuations from the first principle , they may be related to timescales of developing mhd instabilities at different spatial scales . we would like to thank zhuo li , shuang - nan zhang , yi - zhong fan , fang - jun lu and you - jun lu for helpful discussions . this work is partially supported by the natural sciences foundation of china grants 10773011 , 11075169 , 10633040 and 10921063 , the 973 project grants 2010cb833000 and 2009cb824800 , the nsf grant ast-0908362 and nasa grants nnx10ad48 g and nnx10ap53 g at unlv .
Abstract:

Answer: a statistical scenario is proposed to explain the @xmath0-ray variability and flares of the crab nebula , which were observed recently by the fermi / lat . in this scenario electrons are accelerated in a series of knots , whose sizes follow a power - law distribution . these knots presumably move outwards from the pulsar and have a distribution in the doppler boost factor . the maximal electron energy is assumed to be proportional to the size of the knot . fluctuations at the highest energy end of the overall electron distribution will result in variable @xmath0-ray emission via the synchrotron process in the @xmath1 mev range . since highly boosted larger knots are rarer than smaller knots , the model predicts that the variability of the synchrotron emission increases with the photon energy . we realize such a scenario with a monte - carlo simulation and find that the model can reproduce both the two @xmath0-ray flares over a period of @xmath2 year and the monthly scale @xmath0-ray flux fluctuations as observed by the fermi / lat . the observed @xmath0-ray spectra in both the steady and flaring states are also well reproduced .


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: two types of galaxies account for the vast majority of star formation at @xmath2 : the lyman break galaxies ( lbgs ; steidel et al .  1996 , 2003 ) , and the ultra / hyperluminous infrared galaxies detected at submillimeter and millimeter wavelengths ( blain et al . 2002 and references therein ; bertoldi et al . as both selections identify strong star formation , it is tempting to assume they are sampling a common population . however , current observations do not show a large overlap between the two samples . only a small percentage of lbgs are detected with scuba , even though their high star formation rates predict they should be far - infrared ( fir ) luminous ( chapman et al . 2000 ; peacock et al . the reverse is true as well  the majority of ultraluminous infrared sources ( ulirgs ; @xmath3 ) are so highly extincted that they would be missed in rest - frame uv - selected surveys for high redshift galaxies like the lbgs ( meurer et al . 1999 ; goldader et al . this difference may imply that the two methods select galaxies at different stages of evolution , or with intrinsically different physical characteristics . the spitzer space telescope gives us the first opportunity to study both populations across a wide range of mid - infrared ( mir ) wavelengths . recently , six extremely bright lbgs were discovered in sloan digital sky survey ( sdss ) spectra ( bentz , osmer , & weinberg 2004 ; hereafter bow ) . the spectra of the six sources are typical of lbgs , with bright rest - frame uv continua and clearly detected stellar and interstellar absorption lines . they appear 2 - 3 magnitudes brighter than the most luminous lbgs in the steidel et al .  sample . none of the six show spectral features that would be expected for agn , such as broad or high ionization emission lines . bow find no morphological or environmental indicators which suggest lensing . by contrast , the only comparably bright known lbg is ms1512-cb58 ( yee et al . 1996 ; hereafter cb58 ) ; it is strongly lensed by a dense foreground cluster ( seitz et al . 1998 ) , and has an obvious arclike shape which is discernible from the ground .    in this paper , we present the first spitzer detection of an lbg at 16  with the spitzer infrared spectrograph ( irs ; houck et al . 2004 ) `` blue '' peak up filter which is centered at 16 . `` blue '' peak up array . one of the sdss lbgs , j134026.44 + 634433.2 ( hereafter s - lbg-1 ) , was selected as a science verification target for spitzer . the object has a spectroscopic redshift of z=2.79 . throughout the paper , we will adopt a flat , @xmath4-dominated universe ( @xmath5 km s@xmath6  mpc@xmath6 , @xmath7 ) . we obtained 16  imaging of s - lbg-1 on 14 november 2003 with the irs onboard the spitzer space telescope . between exposures , we moved the telescope in a four point dither pattern , with separations of 4 . at each point in the pattern , two nod positions were observed , and two cycles were taken at each nod position . a total of 16 exposures were taken , with integration times of one minute each . the data were reduced using the irs pipeline at the spitzer science center ( see chapter 7 of the spitzer observer s manual ) . individual frames were registered based on the reconstructed pointing ( accurate to better than 1 ) and refined with centroiding of objects . bad pixels were masked before mosaicing . we have not removed the small ( @xmath8% ) geometric distortion in the images . we have several times redundancy for each pixel on the sky , and we remove faint latent images caused by bright sources in pu observations prior to our program .    photometric calibration utilized routine observations of a bright standard . we measured the flux in 16  images of hd 46190 to calculate the zeropoint by comparison with template spectra based on normalized kurucz ( 1993 ) models ( cohen et al . a curve of growth was used to estimate total flux for a point source from aperture photometry . we also obtained supporting near - infrared ( nir ) data in the j- and k - bands with the wirc instrument ( wilson et al . 2003 ) on the hale 200-inch telescope . the observations were taken on 27 january 2004 . the field was observed for five dithered exposures in each filter . j - band exposures were two coadds of one minute each , and k - band were two coadds of 30 seconds each . the night was not photometric , but conditions were stable over the few minutes of the observations . the wirc field of view contains @xmath9  2mass stars , which we used to calibrate the photometry and estimate the extinction . the final , registered 16  image ( figure 1 ) has an rms noise of 15  @xmath10jy in an aperture with 5  radius , and 25  @xmath10jy in an aperture of 9  radius . s - lbg-1 is detected , with a flux of @xmath1  mjy . seven other objects are detected ( @xmath11  or more ) in the field , one of which is very bright in the mir and quite faint at optical wavelengths . of these , two have no counterpart in the sdss or wirc imaging . five additional sources are marginally detected ( @xmath12 ) , without optical counterparts . however , given the faintness of these sources it is not surprising that deeper optical / nir data will be needed to identify them . the reddest object , j134025.02 + 634358.3 ( hereafter red-1 ) , is discussed further in section [ sec : ero ] . photometric data from irs and wirc are given in table 1 . the wirc imaging shows no objects in the field that were not detected in sdss . there is no indication that s - lbg-1 is resolved , but it is not expected to be . the 16  imaging is diffraction limited , and the point spread function ( psf ) has a full width at half maximum ( fwhm ) of 4 . s - lbg-1 is compact in both the sdss and wirc imaging , which have fwhm of @xmath13  and @xmath14 , respectively . at the redshift of s - lbg-1 , z=2.79 , one arcsecond corresponds to @xmath15   kpc . the 16  filter extends from 13.3 to 18.6  at fwhm , or 3.5 to 4.9  in the rest frame . this wavelength range encompasses br@xmath16  but does not include the polycyclic aromatic hydrocarbon ( pah ) feature centered at 3.3  ( tokunaga et al . the filter samples the rayleigh - jeans tail of the stellar photospheric continuum emission emission , which can be considerable if the galaxy contains a substantial population of old stars . alternatively , 4  marks the beginning of the rise in emission from warm dust , heated either by very young stars or by an active nucleus . the spectral energy distribution ( sed ) of s - lbg-1 is consistent with a strongly starbursting galaxy ( figure 2 ) . in the figure we show a direct comparison with the spectrum of ngc 5253 ( wu et al . 2002 ) , a `` benchmark starburst '' ( calzetti et al . 1999 ) with regions that have undergone intense star formation in the last 100 myr ( see calzetti et al . 1997 and the reference therein ) , and a substantial older stellar population . it has mir emission from dust heated by massive , young stars ( crowther et al . the sed of s - lbg-1 is redder than ngc 5253 in the uv , but the two match quite well if one magnitude of additional extinction is applied . the rest - frame optical color of s - lbg-1 ( the observed j- and k - band ) is bluer than that of ngc 5253 , perhaps implying a lower fraction of old stars . ellingson et al .  ( 1996 ) showed that the optical / nir sed of cb58 could be fit with a recent ( 10 myr ) secondary burst of star - formation added to a large ( 85% by mass ) stellar population with an age of @xmath17   gyr . the young age of the source was later confirmed by uv absorption line studies ( pettini et al . 2002 ) . a marginal detection of cb58 with isocam ( bechtold et al . 1998 ) , at rest - frame 1.8 and 3 , favored a higher percentage of younger stars . taken at face value , the cb58 sed is not consistent with s - lbg-1 unless cb58 rises steeply beyond the rest - frame k - band . the inconsistency could result from a higher dust content in s - lbg-1 and likely a dissimilar spatial distribution of dust , which may not be surprising given the different star formation states of the sources . the intrinsic ( corrected for lensing magnification ) sfr of cb58 is @xmath9  ( teplitz et al .  2000 ) , compared to the @xmath18  in s - lbg-1 ( without correcting for extinction ) . the s - lbg-1 sed is also consistent with a relatively dust - free seyfert galaxy ( figure 2 ) . however , bow see no high excitation lines in the uv . it is also possible that the observed 16  flux of s - lbg-1 indicates a `` buried '' agn that is not seen in the uv . bow discuss the possibility that the s - lbg s are bal qso s , but conclude that they would be unique amongst that class of objects if they were . nonetheless , the luminosity of s - lbg-1 may favor the possibility of agn activity dominating the mir flux . in addition , shorter wavelength spitzer observations of s - lbg-1 are needed to conclusively rule out a very steep mir slope , given the low rest - frame 3  flux of cb58 . laurent et al . ( 2000 ) show that the presence of an agn is revealed in the mir by an excess of emission in the 3 - 6  range . this has been attributed to hot dust emission heated to near sublimation temperatures ( @xmath18  k for silicates and @xmath19  k for graphite ) by the accretion disk . the presence of this continuum has been detected in a number of nearby galaxies hosting an active agn including ngc1068 and ngc4151 ( le floch et al . 2000 , alonso - hererro et al . 2003 ) , as well as centaurus a ( mirabel et al . 1999 ) . however , such excess usually has a positive slope in the 3 - 6  range , making it challenging to identify in distant systems because the old stellar population of the bulge of the galaxy may be missinterpreted as an excess of thermal emission . furthermore , there has been evidence that as the luminosity of dust enshrouded ir galaxies increases beyond @xmath20  , so does the probability that an agn contributes substantially to the heating of the dust ( sanders et al . 1988 ; veilleux , kim , & sanders 1999 ; tran et al . 2001 ) . if the brightest lbgs harbor buried agn , it might indicate that the typical percentage of agn in lbg searches , @xmath21% ( steidel et al . 2002 ) , is underestimated . in that case , the inferred   of lbg systems , and thus the large correction to their inferred contribution to the global density of star formation , would be overestimated .      the prodigious star formation rates of s - lbg-1 ( @xmath22  ; bow ) implies that it must be generating enough uv radiation to put it in the same luminosity class as ulirgs ( @xmath23 ) and perhaps enough to be comparable to the hyper - lirgs found among the scuba galaxies ( @xmath24 ; blain et al . it is not certain , however , how much of this radiation is absorbed and reradiated by dust . the optical to mid - infrared sed of s - lbg-1 is not consistent with arp 220 , given the latter s high extinction in the uv . however , the sed over our sampled wavelength provides only a tenuous indication of the bolometric luminosity . if the s - lbgs in general are not fundamentally different than ulirgs , then one must ask why is nt their large uv luminosity absorbed and reradiated by dust . lbgs are found to be 5 - 20 times underluminous ( less dusty ) in the ir for their uv slopes ( baker et al .  2001 ; van der werf et al .  2001 ) , while ulirgs are equally overluminous ( goldader et al . one possibility is that they are at a different stage of evolution . star formation in @xmath25  ulirgs is triggered by merger activity ( flores et al . 1999 ) . in such systems , the dust is dynamically mixed throughout and can easily absorb most of the uv radiation emitted from actively star - forming hii regions . if the lbgs have a less homogeneous distribution of dust , more uv radiation may escape . another possibility is that the extreme uv radiation field of the most luminous lbgs heats the dust to higher temperatures than in typical starbursts , suppressing the fir emission relative to the uv ( baker et al . this would have the effect of shifting the peak of the reradiated luminosity closer to the spitzer wavelengths . finally , a significantly lower dust content  either the result of extremely low metallicity , or the destruction of grains by the uv radiation field ( as suggested by bow )  would allow greater uv brightness .      in the sdss image , bow find no evidence to suggest strong lensing . one might expect that if the lens were highly extincted , it could be detected in the near- or mid - ir . however , we find no evidence of a massive source close in projection to s - lbg-1 that could be identified as a gravitational lens . the 16  psf does not allow us to identify close ( @xmath26 ) companions , but such objects would likely have been seen in the k - band if they were present . the critical radius is a few arcseconds even for a massive elliptical ( blandford & narayan 1992 ; eisenhardt et al . the brightest mir source in the field , red-1 , lies 30  away . red-1 is somewhat unusual in both brightness and color . it has a flux density of 3 mjy at 16 , but is faint in the optical ( @xmath27 ) . the surface density of such sources is quite low . iso observations of the elais field detected @xmath28   sources down to 3 mjy per square degree , or one per @xmath29  pu fields ( la franca et al . down to 1 mjy , the source counts are only ten times higher ( elbaz et al .  1999 and the reference therein ) . few 3 mjy sources are as red as red-1 . la franca et al . find only 18%  of 15  sources brighter than 1 mjy have @xmath30   magnitudes fainter than 23 .    the morphology of red-1 is hard to quantify . it is faint and has low surface brightness in both the optical and nir . it is detected only at the limit of the wirc and sdss images . it shows no central central concentration and is clearly extended , covering several square arcseconds . the 16  filter samples the @xmath31  pah feature and very small grain continuum emission at @xmath32  and the @xmath33   pah features at @xmath34 . at the higher redshift , red-1 would fall into the hyper - lirg luminosity class , but at @xmath32  it would be only @xmath35  ( chary & elbaz 2001 ) . in addition , the object is not extremely red in r - k as would be expected for a @xmath34  ulirg . we have detected an extremely luminous lyman break galaxy in the first spitzer 16  image of a uv - selected source at high redshift . the mir data point is consistent with the large star formation rate inferred from the uv continuum . the sed confirms that the star formation is likely not extincted enough for the object to be considered a ulirg . we find no evidence for a strong lens in the field . future spitzer observations will be crucial to understanding the connection between vigorously star forming lbgs and ulirgs . deep imaging with irac ( fazio et al . 2004 ) and mips ( rieke et al .  2004 ) will potentially detect moderate luminosities lbgs . the 16   imaging capability of the irs is a powerful additional mode for the study of these objects . the detection of mir emission in s - lbg-1 demonstrates that the brightest lbgs will be observable by irs spectroscopy , allowing more detailed studies of their dust properties . such observations will open a new window onto a class of objects that may be different from the infrared luminous sources that will be more commonly studied with spitzer . we thank l. yan , d. frayer , and j. colbert for helpful suggestions . this work is based in part on observations made with the spitzer space telescope , which is operated by the jet propulsion laboratory , california institute of technology under nasa contract 1407 . support for this work was provided by nasa through an award issued by jpl / caltech . bechtold , j. , elston , r. , yee , h. k. c. , ellingson , e. , & cutri , r. m.  1998 , in `` the young universe : galaxy formation and evolution at intermediate and high redshift '' . eds s. dodorico , a. fontana , and e. giallongo . asp conference series ; vol . 146 , p.241
Abstract:

Answer: we present the first rest - frame @xmath0  detection of a lyman break galaxy . the data were obtained using the 16  imaging capability of the spitzer infrared spectrograph . the target object , j134026.44 + 634433.2 , is an extremely luminous lyman break galaxy at z=2.79 first identified in sloan digital sky survey spectra ( bentz , osmer , & weinberg 2004 ) . the source is strongly detected with a flux of @xmath1  mjy . combining spitzer and sdss photometry with supporting ground - based j- and k - band data , we show that the spectral energy distribution is consistent with an actively star - forming galaxy . we also detect other objects in the spitzer field of view , including a very red mir source . we find no evidence of a strong lens amongst the mir sources .