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 bounce characteristics of baseballs have a large influence in baseball games ; thus , baseball organizations often establish rules concerning official balls . for example , in major league baseball ( mlb ) , the baseballs are made by tightly winding yarn around a small core and covering it with two strips of white horsehide or cowhide@xcite .    for the estimation of the bounce characteristics , the coefficient of restitution @xmath0 is widely used , which is defined as @xmath1 where @xmath2 and @xmath3 are the speeds of incidence and rebound , respectively , in a head - on collision of a ball with a plane . note that the coefficient of restitution determines the loss of translational energy during the collision . the coefficient of restitution depends on the kind of material and the internal structure of the ball , as well as other factors such as the impact speed@xcite , impact angle@xcite , temperature of the ball@xcite , and humidity at which the balls are stored@xcite . various baseball organizations officially determine the range of the coefficient of restitution of baseballs . for example , the coefficient of restitution of an mlb baseball is required to be @xmath4@xcite . regarding japanese professional baseballs , the nippon professional baseball organization ( npb ) first introduced their official baseball in 2011 , which was used in both pacific and central league . table [ tb1 ] shows a chronological table indicating the average coefficient of restitution for baseballs used in japanese professional baseball games from 2010 to 2013@xcite , along with the annual number of home runs@xcite . clearly , the number of home runs decreased drastically in 2011 compared with 2010 , although the difference in the average coefficient of restitution is only on the order of @xmath5 . the average coefficient of restitution increased in 2013 because the npb made a baseball equipment manufacturer change the specification of the baseballs in order to increase the level of offense in baseball games . .coefficient of restitution of baseballs and number of home runs in japanese professional baseball games . [ cols="<,^,^,^,^,<",options="header " , ]     [ tb2 ]    next , we show the distributions of the horizontal position @xmath6 and the vertical position @xmath7 of the pitched balls in the pitching zone . figure [ fig3 ] shows the distribution of the horizontal position of the pitched balls , where the horizontal axis indicates the positions of the grids in fig . [ fig1 ] , which are numbered from left to right . the open circles indicate the calculated probabilities as a function of the horizontal position . to obtain the distribution function approximating these data , we fitted the mixture distribution of the normal distributions , which is expressed as @xmath8 where @xmath9 and @xmath10 are the normal distributions with the same standard deviation : 0.77 . the means of @xmath9 and @xmath10 are @xmath11 and @xmath12 , respectively . in fig . [ fig3 ] , the black solid curve shows eq . ( [ md2 ] ) with the mixing parameter @xmath13 , whereas the blue and red curves show the first and second terms , respectively , in the right - hand side of eq . ( [ md2 ] ) . the coefficient of determination of the fit is equal to 0.986 , which indicate a good fit . this indicates that the pitchers tended to choose the right and left sides of the strike zone with a similar probability to prevent home runs .    ) with @xmath13.,width=264 ]    on the other hand , fig . [ fig4 ] shows the distribution of the vertical position for the pitched balls . the probability is almost constant , except at @xmath14 , which is the highest position . this indicates that the pitchers tended to choose heights all over the strike zone with a similar probability . the frequency at @xmath14 is relatively small to avoid the risk of long ball hitting . on the other hand , the frequency at @xmath15 is almost same as those at @xmath16 inside the strike zone , which may be attributed to that the balls at @xmath15 are difficult to strike . we performed a batting simulation according to the results presented in the previous section . figure [ fig5](a ) shows our simulation setup , in which a home plate was placed at the origin of a cartesian coordinate system . the initial position of the center of mass of a baseball was set at @xmath1718.44 m , 0 m , 1.8 m@xmath18 according to the official baseball rules in japan and the average height of japanese professional baseball players . the baseball was pitched in the negative direction of the @xmath19 axis with an initial velocity of @xmath20 , where @xmath21 is the unit direction vector of the pitch , which will be defined later . here , @xmath22 is randomly chosen according to the distribution function given by eq . ( [ md ] ) with the parameters shown in table ii . .,width=453 ]    the pitch direction was defined as follows . a ball needs to be thrown from about shoulder height and be launched in an almost horizontal direction to cross the plate at the correct height@xcite . thus , we assume that a pitcher throws a ball toward a point @xmath230 , @xmath24 [ m ] , 1.8 m@xmath18 that is on the @xmath6-@xmath7 plane . a thrown ball travels a curved path to cross the plate at a height less than the initial height due to the gravitational force ( fig . [ fig5](b ) ) . @xmath24 is the random variable selected by the distribution function @xmath25 defined in eq . ( [ md2 ] ) . the unit vector of the pitch direction was defined as @xmath26 , where @xmath27 . we consider the bat to be an uniform cylinder 1 m long and the diameter of the base to be @xmath28 m , in accordance with the official baseball rules . the bat is placed along the @xmath6 axis , and its center of mass is positioned at @xmath290 , 0 , @xmath30 [ m]@xmath18 , where @xmath30 is defined later . for later discussions , we respectively label bases of the bat as @xmath31 and @xmath32 so that the @xmath6 component of the center of @xmath32 is larger than that of @xmath31 . let us assume that a thrown ball crosses the @xmath6-@xmath7 plane between the time @xmath33 and @xmath34 , where @xmath35[s ] is the time step of our simulation ( fig . [ fig5](b ) ) . first , we calculate @xmath36 using the @xmath7 component of the intersection between the @xmath6-@xmath7 plane and the line segment connecting @xmath37 and @xmath38 , where @xmath37 is the position of the center of mass of a thrown ball at the time @xmath33 . next , we determine @xmath30 as follows : @xmath39 where @xmath40 represents the random numbers chosen from the normal distribution with a mean of @xmath41 and a standard deviation of @xmath42 m , which is the diameter of an official baseball in japan . figure [ sys ] shows a schematic of a collision between a ball and a bat from the viewpoint of the home base . we assume the bat swings around an axis passing through the center of b. the bat swings with the angular velocity @xmath43 with @xmath44 rad / s , which is the typical value of bat swings . here @xmath45 is the angle between the direction of @xmath46 and @xmath47 , the value of which is randomly chosen from the range @xmath48 . when a ball collides with the bat , we calculate the vector @xmath49 , which is defined by the vector from the center of @xmath32 to the foot of the perpendicular line passing through the center of the ball to the central axis of the bat . the velocity of the bat on the collision is defined by @xmath50 . basically , a pitched ball obeys the equation of motion @xmath51 where @xmath52 is the mass of the ball . the three terms on the right - hand side of eq . ( [ eqofmo ] ) represent the drag force , the magnus force , and the gravitational force , respectively . the drag force @xmath53 generated by the resistance from the air decreases the ball speed . @xmath53 is expressed as @xmath54 where @xmath55 , @xmath56 , @xmath31 , and @xmath57 are the drag coefficient , the density of air , the cross section of a baseball , and the unit vector in the direction of the velocity , respectively . we use @xmath58 , corresponding to a baseball pitched at a high speed of approximately 40.2 to 44.7 m / s@xcite . our choice of the value will be discussed in later section . in addition , we use @xmath59 kg / m@xmath60 , which is the value at 0 c@xmath61 and 1 atm .    on the other hand , the magnus force @xmath62 generated by the backspin of the baseball enhances the flight of the ball after it is thrown and batted@xcite . @xmath62 is expressed as @xmath63 where @xmath64 and @xmath65 are the lift coefficient and the unit vector perpendicular to @xmath57 , respectively . we used @xmath66 , which is a typical value for a spinning ball@xcite . the gravitational acceleration is indicated by the vector @xmath67 m / s@xmath68 ) . we numerically solved eq . ( [ eqofmo ] ) using the velocity verlet scheme@xcite . a thrown ball changes its flight direction when it collides with a bat . we assume that the ball collides with the bat when all of the following conditions are fulfilled :    1 .   the distance between the center of mass of the ball and the central axis of the bat is less than the sum of the radius of the ball and that of the bat . 2 . the @xmath6 component of the ball s position , @xmath69 , satisfies @xmath70 m @xmath71 m.    when these conditions are fulfilled , the ball is reflected by the averaged repulsive force during @xmath72 , @xmath73 where @xmath74 , @xmath75 and @xmath76 are the reduced mass of the ball and the bat , the normal unit vector that is perpendicular to the tangential plane between a colliding ball and a bat , and the relative velocity of the ball to the bat , @xmath77 , respectively .    in eq . ( [ collision ] ) , we use the coefficient of restitution between a ball and a bat @xmath78 , which is defined by @xmath79 where @xmath80 with the moment of inertia @xmath81 of the bat around the center of b@xcite . the derivation of eq . ( [ collision ] ) and the algorithm for the collision are summarized in appendices a and b , respectively . when a batted ball falls on the ground , we calculate the distance @xmath82 between the point of landing and the home plate . and @xmath83 , respectively.,width=340 ]    figure [ fig6 ] shows the probability densities of the flying distance @xmath82 of batted balls ( @xmath84 and @xmath83 ) , each of which was calculated from 1,000 samples . the highest peak position shifts from 75 m to 125 m with increasing coefficient of restitution of a ball . the total frequencies in @xmath85 , which can be regarded as the number of home runs , increases with increasing coefficient of restitution indicating an increase in the home - run probability . thus , we investigate the relationship between the home - run probability and @xmath0 by changing the value of @xmath0 from @xmath86 to @xmath87 .     and home - run probability @xmath88 . solid curve shows eq . ( [ pf_hr ] ) with @xmath89 m and @xmath90 m. , width=340 ]    figure [ fig7 ] shows the relationship between the home - run probability and @xmath0 . here , each data point was calculated from 1,000 samples . we consider a sample with @xmath85 m as a home run and define the home - run probability as the ratio of the number of home runs to 1,000 . as shown in fig . [ fig7 ] , the home - run probability increases with @xmath0 , which is intuitively understood .    here , we estimate the functional form of @xmath88 using a simple theoretical argument . let us suppose that a projectile is launched from the ground with a launch speed @xmath91 and a launch angle @xmath92 under the air friction . the range of the trajectory under the air friction can be estimated by a linear function of @xmath91 ( ref . 2 ) , so that we roughly estimate the range of the trajectory as a linear function of the coefficient of restitution @xmath0 as @xmath93 with the constant @xmath94 . here we assume that the speed of the thrown ball upon the collision with the bat is almost constant . we assume that the probability density of the flying distance of a batted ball which will land at approximately @xmath95 m can be approximated by the normal distribution as @xmath96.\ ] ] by integrating eq . ( [ pdf_hr ] ) from @xmath95 m to infinity , we obtain the probability function @xmath97,\label{pf_hr}\end{aligned}\ ] ] where @xmath98 is the error function . we assume both @xmath94 and @xmath99 as fitting parameters . the solid curve in fig . [ fig7 ] corresponds to eq . ( [ pf_hr ] ) drawn with @xmath100 m and @xmath101 m , which indicates a good estimation . figure [ fig3 ] shows the distribution of the horizontal index @xmath6 of pitched balls at the pitching zone , which can be approximated by the mixture distribution of the normal distributions . because of the method of division of the pitching zone in the database that we used , the number of points was too scarce to identify a more detailed distribution . to obtain a more detailed distribution of the position of pitched balls , we can use more accurate data from pitchf / x@xcite . however , from the viewpoint of pitcher strategy , it may be desirable to aim at the edges of the strike zone ; thus , we may find two peaks at the right and the left edges of the strike zone in the real distribution . figure [ fig6 ] shows the probability density distributions @xmath102 of the distance @xmath82 , each of which has some peaks . for comparison with real data , a japanese group attempted to draw a histogram of the distance of batted balls using data from ultimate zone rating ( uzr)@xcite in 2009 and 2010 . their results show that the histogram has two peaks around @xmath103 and @xmath104 m , which is qualitatively close to our simulation results at @xmath83 . however , their results show that the probability of finding a ball at the distance @xmath105 is very low , although our simulation results show that the probability density at @xmath83 is almost same in the same range . notice that the group constructed the histogram using the positions at which the fielders caught or picked up the batted balls , which are sometimes different from the points of landing . thus , their histogram can be considered as the mixture of the histograms of the position of pickup by the infielders and the outfielders , which may cause the discrepancy between their results and our results . in the uzr analysis , the probability of finding a batted ball between 91.44 and 121.92 m was 0.179 . in our simulation , on the other hand , the home - run probability @xmath88 was @xmath106 at @xmath107 , which is the closest value to the averaged coefficient of restitution in 2010 . the discrepancy between the two values is attributed to that the data in 2009 and 2010 were included in the uzr analysis .    finally , our model is a simplified model based on simple rigid - body collisions , which ignores the rotation of a batted ball . our model includes the lift coefficient @xmath66 as a fixed value , although @xmath64 can vary depending on the rotation of the ball@xcite . for example , recent experiments have demonstrated that @xmath64 strongly depends on the spin parameter , @xmath108 , where @xmath109 and @xmath110 are the angular velocity and the radius of a ball , respectively@xcite . on the other hand , @xmath55 also varies according to reynolds number , although our model includes the drag coefficient @xmath111 as a fixed value . especially , the drag force on a baseball is known to drop sharply at speeds typical for thrown and batted balls@xcite . to improve our model , we will need to incorporate @xmath64 and @xmath55 dependant on the motion of thrown and batted balls . in addition , constructing a model based on elastic collisions in which the friction between the ball and bat is considered will yield more accurate results . in this paper , we analyzed real data for the speed and the course of pitched balls in professional baseball games in japan . our results show that the distribution of the ball speed can be approximated by the mixture distribution of two normal distributions . in addition , we found that the horizontal position of pitched balls in the pitching zone obeys the mixture distribution of two normal distributions , each of which has a peak at the edge of the strike zone . we simulated collisions between baseballs and a bat , where the statistics of the pitching obey our analyzed results . we finally obtained the probability density distribution @xmath102 of the distance @xmath82 between the home plate and the landing point of the batted balls to calculate the home - run probability as a function of the coefficient of restitution . by using a simple theoretical argument with the assumption that @xmath102 around @xmath95 m can be approximated by the normal distribution , we quantified the home - run probability as a function of the coefficient of restitution @xmath0 of baseballs . as stated in the previous section , we will obtain more accurate results in future works by improving our model . developing a methodology to calculate the home - run probability may yield useful information for designing baseballs . in ref.1 , the author derived the rebound velocity of a ball in a head - on collision with a bat@xcite . in this appendix , we derive the rebound velocity of a batted ball in a three dimensional oblique collision . note that we ignore the rotation of a ball and a bat in our argument . figure [ figa1 ] shows a ball of mass @xmath52 colliding with a bat of mass @xmath112 . here we denote the colliding velocities of the ball and the bat as @xmath113 and @xmath114 , respectively . assuming that the ball experiences the averaged force @xmath115 from the bat during the duration @xmath72 , we can describe the rebound velocities of the ball and the bat , @xmath116 and @xmath117 , as @xmath118 from the definition of the averaged force and newton s third law of motion . note that the total momentum of the system is conserved after collision , @xmath119 .     and a ball of mass @xmath52 . the velocities of the centers of mass of the ball and bat are denoted by @xmath113 and @xmath114 , respectively . @xmath75 is the normal unit vector perpendicular to the tangential plane . , width=226 ]    by subtracting eq . ( [ a2 ] ) from eq . ( [ a1 ] ) and introducing the reduced mass @xmath120 defined by @xmath121 , we obtain @xmath122 where @xmath76 and @xmath123 are the relative velocities of the ball to the velocity of bat before and after collision , respectively . by introducing the normal unit vector @xmath75 perpendicular to the tangential plane of the bat and ball , the scalar projection of @xmath123 onto @xmath75 is calculated as @xmath124 where we used @xmath125 . ( [ a4 ] ) in the definition of the coefficient of restitution @xmath78 between ball and bat , @xmath126 we obtain the equation @xmath127   \left [ \beta + { \tilde { \bf v } } \cdot { \bf n } ( 1+\tilde{e})\right ] = 0,\label{qd}\ ] ] where @xmath128 . the two solutions of the quadratic equation eq . ( [ qd ] ) are respectively written as @xmath129 in the two solutions , only @xmath130 corresponds to the averaged force in the collision because @xmath131 becomes @xmath41 ( @xmath132 ) when @xmath133 , which will cause the penetration of the ball into the bat . thus , the averaged force is calculated as @xmath134 will penetrate into the bat at the time @xmath135 , which is unrealistic . to obtain the normal force acting on the ball , we need to calculate the time @xmath136 when the ball touches the bat . , width=226 ]    in appendix a , we explained the way to calculate the averaged force acting on the ball during collision . however , it is difficult to obtain the unit normal vector @xmath75 in naive calculation because the ball can penetrate into the bat during the simulation time step , which is due to the nature of the finite difference approximation of derivatives . in this appendix , we explain the algorithm of the collision between a bat and a ball to avoid penetration used in our simulation . figure [ figa2 ] shows a schematic of a collision between a ball and a bat , where a ball at the time @xmath33 is going to collide with the bat . when the first condition of collision ( see section iii.c ) is fulfilled , the ball has penetrated into the bat at the time @xmath135 . thus , we put the ball back to the previous position at the time @xmath33 to determine the remaining time @xmath137 before collision . from the condition that the distance between the center of the bat and the center of the ball at the time @xmath136 equals to the sum of the radius of the bat and that of the ball , we obtain a quadratic equation for @xmath137 . here we assume that the ball travels linearly with a constant velocity during @xmath137 to simplify the calculation . among the two solutions for the quadratic equation , we choose the positive one for @xmath137 which has a physical meaning .    from the position of the ball at the time @xmath136 , we can calculate all the variables such as @xmath75 , @xmath138 , @xmath114 and @xmath78 in eq . ( [ collision ] ) . the ball reflects from the bat according to the following algorithm . first , we put the ball back to the position at the time @xmath33 . next , we apply the half of the averaged force to the ball during @xmath139 so that the impulse from the bat becomes constant . with this two - step time evolution , the ball bounces from the bat without penetration . d. frenkel and b. smit , _ understanding molecular simulation _ , 2nd edition ( academic press , san diego , 2002 ) . a. m. nathan ,  characterizing the performance of baseball bats " , am . j. phys . , * 71 * , 134143 ( 2003 ) .
Abstract:

Answer: in baseball games , the coefficient of restitution of baseballs strongly affects the flying distance of batted balls , which determines the home - run probability . in japan , the range of the coefficient of restitution of official baseballs has changed frequently over the past five years , causing the number of home runs to vary drastically . we analyzed data from japanese baseball games played in 2014 to investigate the statistical properties of pitched balls . in addition , we used the analysis results to develop a baseball - batting simulator for determining the home - run probability as a function of the coefficient of restitution . our simulation results are explained by a simple theoretical argument .


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: qed in three space - time dimensions has become a testing laboratory for certain nonperturbative aspects of qft like chiral symmetry breaking  @xcite , bound states  @xcite or confinement  @xcite and for various approximation schemes . due to the dimensionality of the coupling constant ( as @xmath0 ) it is a superrenormalizable theory so it is free of infinite renormalization ambiguities typical for four - dimensional theory . it can be formulated in two inequivalent versions : with two- and four - component fermions  @xcite . the properties of the theory are different in these two cases . to the investigation of both versions much attention has been payed over the last twenty years . the work has been especially concentrated on nonperturbative solutions of dyson - schwinger ( ds ) equations with different approximations incorporated in the theory as quenched approximation , rainbow approximation , @xmath1 expansion and various models of vertex function  @xcite . particularly often the multi - flavor theory with the limit @xmath2 has been used , since it avoids infrared problems , which usually become troublesome in lower number of dimensions . the other important point in this analysis has been the unpleasant gauge dependence of the nonperturbative results which constitute the common problem of the approximated studies based on ds equations  @xcite . as it is well known , ds equations constitute an infinite set of relations involving green s functions which form the ` inverted ladder ' type structure : the @xmath3-point functions depend on @xmath4-point ones and so on up to infinity . this set can not be solved without the truncation of such a hierarchy . but even such truncation , which turns the infinite set of equations into only a couple of them , leads to the system which is far from being trivial and requires further simplifications . usually this truncation is accomplished by performing certain assumption for the vertex function , which is chosen most often in the form satisfying ward - takahashi identity . naturally this identity fixes only longitudinal part of the vertex leaving the transverse part a subject of further discussion and improvements  @xcite .    in our previous paper  @xcite we applied a method elaborated earlier in qed@xmath5  @xcite to qed@xmath6 with two - component fermions . it consists on the following five steps :    1 . a certain infrared form of two basic propagators , @xmath7 and @xmath8 ( suggested by perturbative calculations or other methods ) as dependent on a couple of unknown parameters is assumed , 2 . the fermion propagator is represented in the spectral form with one known spectral density @xmath9 , which is in general possible in the infrared domain , 3 . salam s form of the vertex function  @xcite , together with @xmath7 and @xmath8 , is substituted into the first two of the set of ds equations , 4 .   from these two equations the set of self - consistent equations for parameters is derived , 5 . the obtained equations are solved numerically or analytically . this method proved to be relatively effective both in qed@xmath5 , where all parameters were correctly found without the necessity of infinite renormalization , and in spinor version of qed@xmath6 , where obtained results stay in general agreement with other works . in the present paper we would like to extend its application to four - component qed@xmath6 , without chern - simons term , i.e. when gauge bosons ( ` photons ' ) remain massless in spite of interaction . the common effect of various simplifications of ds equations is the gauge dependence of the results ( even of the physical observables ) . the sources of this undesirable behavior are the approximations made to the propagators and to the vertex . it seems therefore valuable to test various possible approaches with regard to that particular feature and much work has already been done in this direction ( see the references above ) . this point lies in the scope of interest of the present work too . this paper is organized as follows . in the next section we define the model itself and give the resulting set of ds equations . in section  [ sec : assump ] the infrared green s function in question are formulated up to several unknown parameters . in section  [ sec : solu ] we substitute these green s functions into ds equations and obtain the set of relations for the introduced parameters . in the last section we present analytical and numerical results and some conclusions . the model is defined through following lagrangian density : @xmath10 where @xmath11 is the gauge parameter . the quantities @xmath12 and @xmath13 denote here the _ bare _ fermion mass and the _ bare _ coupling constant respectively . as mentioned in the introduction , the latter for @xmath14 is a quantity with the dimensionality of @xmath0 . that means that in the quantum theory higher terms of perturbation expansions have better ultraviolet momentum dependence in loop integrations and the model is superrenormalizable .    as already told , in the present paper we deal with four - component fermion field , choosing the following representation for gamma matrices used also in four - dimensional qed : @xmath15    there are two other gamma matrices , which anticommute with all above , and which can serve for defining chiral transformations : @xmath16    the fermion mass term , that was chosen in  ( [ eq : lagr ] ) in the form @xmath17 , breaks the chiral symmetry defined by any of the matrices  ( [ eq : gamma35 ] ) . there is , however , the possibility of the other choice for the mass term  @xcite : @xmath18 , where @xmath19=\left(\begin{array}{cc}\openone   &   0\\ 0 & -\openone     \end{array}\right)\ ; . \label{eq : tauma}\ ] ] this term does not break chiral symmetry , since matrix @xmath20 commutes with both @xmath21 and @xmath22 , but it does violate parity symmetry ( contrary to the choice made in  ( [ eq : lagr ] ) ) .    for matrices  ( [ eq : gammas ] ) we have ordinary relations : @xmath23=4g^{\mu\nu}\ ; , \nonumber\\ \mathrm{tr}\,\left[\gamma^\mu \gamma^\nu \gamma^\rho \gamma^\sigma\right]=4\left(g^{\mu\nu}g^{\rho\sigma}-g^{\mu\rho}g^{\nu\sigma}+g^{\mu\sigma}g^{\nu\rho}\right)\ ; , \label{eq : toga}\end{aligned}\ ] ] where we choose for the metric tensor : @xmath24    the trace of the product of an odd number of gamma matrices equals zero as in four dimensions ( this was not the case in the spinor representation , where it was proportional to the antisymmetric tensor @xmath25 ) . additionally we have the identity : @xmath26    using the lagrangian density  ( [ eq : lagr ] ) , one can derive in the standard way  for instance through feynman path integral  the dyson - schwinger equations for propagators  @xcite . for boson propagator we obtain the relation @xmath27,\nonumber\\ \label{eq : dysph}\end{aligned}\ ] ] which may be given the graphical form shown in figure  [ fig : dspf ] . the propagator @xmath8 on the right hand side is not integrated over three - momenta and therefore it can fully be represented through fermion functions ( and the vertex ) .    the ds equation for the fermion propagator does not allow for such a separation , and has the form @xmath28\ ; , \nonumber\end{aligned}\ ] ]    its graphical representation is shown in figure  [ fig : dsef ] . while introducing the fermion self - energy @xmath29 , it can be rewritten in a simpler manner @xmath30 where @xmath31    higher equations are not considered in our present approach . the vertex function @xmath32 is not taken from ds equations , but postulated in the form proposed by salam  @xcite and used afterwards in the so called ` gauge technique '  @xcite . the eventual extension of our method on higher green s functions will be considered elsewhere . the free propagator of massless vector boson has the standard form : @xmath33    as is well known , the interaction with fermions does not change the longitudinal part of @xmath8 . this is guaranteed by the following ward - takahashi ( wt ) identity @xmath34 consequently the nonperturbative propagator may be then written as @xmath35 with certain unknown function @xmath36 . in the infrared domain we assume this function to have a taylor expansion starting from @xmath37 , since for bispinor fermions the boson field remains massless  @xcite and we expect @xmath38 . this is not the case for massive two - component fermions , for which the topological boson mass is generated . taking the first two terms of this expansion , we have @xmath39 where @xmath40 , together with the renormalization constant @xmath41 , will be determined from consistency conditions . from the unitarity we expect the value of @xmath41 to satisfy the restriction @xmath42  @xcite . this expectation was actually confirmed in our previous works on spinor qed@xmath6  @xcite and qed@xmath5  @xcite .    with the above assumptions the inverse of @xmath43 may be written as @xmath44    the free fermion propagator has the usual form : @xmath45 due to the masslessness of the gauge boson and the absence of the mass gap for emission of soft photons , the pole at mass @xmath12 should , in the full propagator , turn into a branch point at @xmath46 , where @xmath47 is a physical mass . this is expected also on the basis of general considerations on the analytical structure of fermion propagator  @xcite and of confinement , which prohibits @xmath7 from having a simple pole . the numerical results , performed in @xmath1 expansion with bare vertex and in euclidean space , suggesting the existence of _ complex _ singularities rather than _ real _ may as well constitute the effect of the coarse approximations made while solving ds equations : the approximations that are inevitable in any nonperturbative approach . such singularities may be a signal of confinement , but are not prerequisite . in the schwinger model , massless electrodynamics in two space - time dimensions , which exhibits confinement of fermions , the _ real _ singularity in the infrared domain in the fermion propagator has been found ( it has the form of  @xmath48 )  @xcite . our infrared assumption for @xmath7 is then : @xmath49 we will see later that it will turn out to be self - consistent .    the values of the exponent @xmath50 and of mass renormalization constant @xmath51 will be established by the requirement of consistency . the presence of @xmath50 in denominator improves the ultraviolet behavior of loop integrals ( we assume that @xmath52 , to be verified _ a posteriori _ ) .    for our further purposes @xmath7 has to be written in the spectral form @xmath53 with the one spectral density : @xmath54\ ; , \label{eq : ro}\end{aligned}\ ] ] sufficient to define the infrared behavior of the propagator . @xmath55 is here the heaviside step function . the particular form of @xmath9 was given in our previous works  @xcite . using @xmath9 , one can write the vertex function , which we need to put into ds equations  ( [ eq : dysph ] )  and  ( [ eq : dyse ] ) . as mentioned in the introduction we use the salam s vertex  @xcite ( with slight and obvious modification resulting from the form  ( [ eq : sp ] ) of @xmath7 , which in our case contains two terms ) : @xmath56 \ ; . \label{eq : verp}\end{aligned}\ ] ] it automatically guarantees the compliance with the wt identity : @xmath57 inserting  ( [ eq : sp ] ) and  ( [ eq : verp ] ) into the right hand side of  ( [ eq : dysph ] ) we obtain    @xmath58    where we have rewritten this equation for inverse propagator , which is easier to handle , because in this case @xmath43 decouples from other functions . the vacuum polarization tensor may be defined as the spectral integral @xmath59 where @xmath60 denotes the usual perturbative tensor : @xmath61    it may be evaluated in the standard way with the use of identities  ( [ eq : toga ] ) , and for instance by performing wicks rotation and introducing feynman parameters . passing back to minkowski space , we have @xmath62    for @xmath63 the @xmath64 integral is well defined . the transversality of @xmath60 is a favorable consequence of using the vertex function in the form  ( [ eq : verp ] ) , satisfying the wt identity  ( [ eq : kve ] ) . performing the spectral integral over @xmath65 in  ( [ eq : pote ] ) with the use of  ( [ eq : ro ] ) , similarly as it was done in  @xcite , we obtain : @xmath66 the integral over feynman parameter @xmath64 leads to the hypergeometric ( gauss ) function @xmath67 and the ds equation  ( [ eq : invpp ] ) may be given the form @xmath68\ ; . \nonumber\end{aligned}\ ] ]    after the substitution of the expression  ( [ eq : dinv ] ) for the left hand side and cancellation of the tensor structures , we are left with the scalar equation , for which we require the adjustment the first two terms of the taylor expansion in @xmath37 . in that way we get two equations for @xmath41 and @xmath40 : @xmath69 where we expanded the gauss function for small momentum according to the formula : @xmath70 and @xmath71 is the euler function . these are two of the set of four equations for parameters of the model , to be solved in section  [ sec : rec ] .      to get other two equations for parameters , we put  ( [ eq : sp ] ) , ( [ eq : verp ] ) and  ( [ eq : dprfull ] ) together with  ( [ eq : dexp ] ) into the ds equation  ( [ eq : dyse ] ) . we obtain , after some simplifications @xmath72_a+[\sigma(p)s(p)]_b+[\sigma(p)s(p)]_c\ ; , \nonumber\end{aligned}\ ] ] where we divided fermion self - energy into pieces coming from different tensor structures in @xmath73 : @xmath74 and @xmath75 from the transverse part and again from the gauge - dependent longitudinal part . after small rearrangement they are :    @xmath76_a=\label{eq : siga}\\   & & iz_3\kappa^2e^2_0\int dm\rho(m)\int\frac{d^3k}{(2\pi)^3}\left [ \gamma^{\mu}\frac{1}{\not\ ! p\ ; + \not\ ! k - m+i\varepsilon}\gamma_{\mu}\frac{1}{(k^2+i\varepsilon)(k^2-\kappa^2+i\varepsilon)(\not\ ! p - m+i\varepsilon ) } - ( m\rightarrow m)\right]\nonumber\\ & & \!\!\!\!\![\sigma ( p)s(p)]_b=\label{eq : sigb}\\ & & iz_3\kappa^2e^2_0\int dm\rho(m)\int\frac{d^3k}{(2\pi)^3 } \left [ \not\ ! k\frac{1}{\not\ ! p\ ; + \not\ ! k - m + i\varepsilon)(k^2+i\varepsilon)^2(k^2-\kappa^2+i\varepsilon)}- ( m\rightarrow m)\right]\nonumber\\ & & \!\!\!\!\![\sigma ( p)s(p)]_c=\label{eq : sigc}\\   & & \frac{ie^2_0}{\lambda}\int dm\rho(m)\int\frac{d^3k}{(2\pi)^3 } \left [ \not\ ! k\frac{1}{(\not\ ! p\ ; + \not\ ! k - m + i\varepsilon)(k^2+i\varepsilon)^2 } - ( m\rightarrow m)\right]\nonumber\ ; .\end{aligned}\ ] ]    a comment should be made here . to avoid technical difficulty while performing wick s rotation in the above momentum integrals , we analytically continued the value of @xmath77 to imaginary values on the upper half plane of complex @xmath77 . now the deformation of the integration contour as is required by passing into euclidean space , is not disturbed by the inappropriate location of poles since all singularities have the ` feynman ' position . after performing the integrals , we will come back to real values of @xmath77 . the procedure was discussed in  @xcite . with this trick each of the above momentum integrals can be performed in an ordinary way known from perturbation theory . it can easily be seen that all integrals are finite without the need of any regularization . omitting details of this standard calculation , we find ( in minkowski space )    @xmath78\ ; , \nonumber\\ i_b=&&\!\!\!\!\!ie^2_0\int\frac{d^3k}{(2\pi)^3 }   \not\ ! k\frac{1}{(\not\ ! p\ ; + \not\ ! k - m + i\varepsilon)(k^2+i\varepsilon)^2(k^2-\kappa^2+i\varepsilon)}=\frac{e_0 ^ 2}{8\pi\kappa^2}\int\limits_0 ^ 1dx\bigg[\left(1-x\,\frac{\not\ ! p(\not\ ! p+m)}{\kappa^2}\right)\times\label{eq : ib}\\ & & \times\left(\frac{1}{(m^2x - p^2x(1-x))^{1/2 } } -\frac{1}{(m^2x - p^2x(1-x)+\kappa^2(1-x))^{1/2}}\right)+\frac{\not\ ! p(\not\ ! p+m)}{2}\frac{x(1-x)}{(m^2x - p^2x(1-x))^{3/2}}\bigg]\ ; , \nonumber\\ i_c=&&\!\!\!\!\ ! ie^2_0\int\frac{d^3k}{(2\pi)^3 }   \not\ ! k\frac{1}{(\not\ ! p\ ; + \not\ ! k - m + i\varepsilon)(k^2+i\varepsilon)^2}=-\frac{e_0 ^ 2}{8\pi}\int\limits_0 ^ 1 dx\bigg[\frac{1}{(m^2x - p^2x(1-x))^{1/2}}\nonumber\\ & & + \frac{\not \ ! p+m)}{2}\frac{x(1-x)}{(m^2x - p^2x(1-x))^{3/2}}\bigg]\ ; . ic}\end{aligned}\ ] ]    now the contributions to @xmath79 may be written as @xmath80_a = \label{eq : ssa}\\ & & z_3\kappa^2\int dm\rho(m)\left[i_a\frac{1}{\not\ ! p - m+i\varepsilon}-(m\rightarrow m)\right]\ ; , \nonumber\\ & & \!\!\!\!\ ! [ \sigma ( p)s(p)]_b = \label{eq : ssb}\\ & & z_3\kappa^2\int dm\rho(m)\left[i_b-(m\rightarrow m)\right]\ ; , \nonumber\\ & & \!\!\!\!\ ! [ \sigma ( p)s(p)]_c = \frac{1}{\lambda}\int dm\rho(m)\left[i_c-(m\rightarrow m)\right]\label{eq : ssc}\ ; .\end{aligned}\ ] ] all the above spectral integrals can been found similarly as in  @xcite , therefore we omit the technicalities .    in the infrared domain , when @xmath81 , the left hand side of the ds equation  ( [ eq : hsf ] ) , after substituting  ( [ eq : spr ] ) , contains two types of singular terms ( the only terms that are important ) : @xmath82 therefore , to avoid lengthy expressions we will not give the results for  ( [ eq : ssa ] ) , ( [ eq : ssb ] ) , and  ( [ eq : ssc ] ) in their full complexity , because it is sufficient for our goal to only pick out from them the identical singular terms . we find : @xmath80_a \approx \frac{e_0 ^ 2z_3}{8\pi}\bigg\{\bigg[\frac{1}{m^2}\,{\cal i}_1(m^2,\kappa^2)-\frac{1}{m^3}\bigg]\times\nonumber\\ & & \times\frac{\not\ ! p(\not\ ! p+m)}{2}\left({\cal i}_2(m^2,\kappa^2)-\frac{1}{m^3}\right)\nonumber\\ & & -\frac{3}{4\beta m}\bigg]\frac{1}{(1-p^2/m^2)^{\beta}}\bigg\}\ ; , \label{eq : ssa1}\end{aligned}\ ] ] @xmath80_b \approx \frac{e_0 ^ 2z_3}{8\pi}\bigg[\frac{1}{m^3}\,\frac{\not\ ! p+m)}{(1-p^2/m^2)^{\beta+1}}\label{eq : ssb1}\\ & & -\bigg(\frac{\not\ ! p+m)}{2m^3}\ , \frac{1-\beta}{\beta}-\frac{1}{4\beta m}\bigg)\frac{1}{(1-p^2/m^2)^{\beta}}\bigg]\ ; , \nonumber\end{aligned}\ ] ] @xmath80_c \approx -\frac{e_0 ^ 2}{8\pi\lambda}\bigg[\frac{\not\ ! p+m)}{m^3}\,\frac{1}{(1-p^2/m^2)^{\beta+1}}\nonumber\\ & & + \bigg(\frac{\not\ ! p(\not\ ! p+m)}{2m^3}\,\frac{\beta-1}{\beta}+\frac{1}{4\beta m}\bigg)\frac{1}{(1-p^2/m^2)^{\beta}}\bigg]\label{eq : ssc1}\ ; , \end{aligned}\ ] ] where @xmath83 refers to diverging terms , when @xmath81 . the functions @xmath84 and @xmath85 have the following form : @xmath86 if the solution is to be self - consistent close to the fermion mass shell , the divergent terms on both sides of ds equation  ( [ eq : hs ] ) must be identical . equating them , and making use of the fact that up to finite terms we have @xmath87 we derive the following two relations for unknown parameters @xmath88 and @xmath50 : @xmath89\ ; , \label{eq : eq3}\\ 1 & \!\!\ ! = & \!\!\ ! \frac{e_0 ^ 2 z_3}{8\pi}\left[-\kappa^2{\cal i}_3(m^2,\kappa^2)+\frac{3}{4\beta m}\left(\frac{1}{\lambda z_3}-2\right)\right]\ ; .\nonumber\\ \label{eq : eq4}\end{aligned}\ ] ] the function @xmath90 is defined as follows @xmath91\nonumber=\\ & & \!\!\!\!\int\limits_0 ^ 1 dx\frac{(1-x)(x+2)}{(m^2x^2+\kappa^2(1-x))^{3/2}}\ ; . \label{eq : ii3}\end{aligned}\ ] ] the four equations we have obtained , i.e  ( [ eq : eq1 ] ) , ( [ eq : eq2 ] ) , ( [ eq : eq3 ] ) and  ( [ eq : eq4 ] ) are sufficient to determine all parameters . we rewrite them with the use of _ renormalized _ quantities : fermion mass @xmath92 , gauge coupling constant @xmath93 and gauge parameter @xmath94 . besides , it is useful to introduce a dimensionless parameter @xmath95 . after executing the parametric integrals in @xmath84 and @xmath96 and performing the reverse analytical continuation in @xmath77 , we get    @xmath97,\label{eq : deleq}\\ 1 = & & \!\!\!\!\!\frac{\zeta}{2}\bigg[\frac{3}{4\beta } \left(\frac{1}{\lambda_r}-2\right)-2\,\frac{4m^2-\kappa^2}{4m^2+\kappa^2}\nonumber\\ & & \!\!\!\!\!-\frac{\kappa^2}{2m^2}\,\ln ( 4m^2/\kappa^2 + 1)\bigg]\ ; , \label{eq : jedeq}\\ z_3 = & & \!\!\!\!\ ! 1-\frac{2\zeta}{3}\frac{\gamma(\beta+1/2)}{\sqrt{\pi}\gamma(\beta+1)}\ ; , \label{eq : z3eq}\\ \frac{\kappa^2}{m^2 } = & & \!\!\!\!\ ! \frac{15   } { 2\zeta}\frac{\sqrt{\pi}\gamma(\beta+1)}{\gamma(\beta+3/2)}\ ; . \label{eq : kapeq}\end{aligned}\ ] ]    this set of equations for @xmath98 , @xmath50 , @xmath41 and @xmath99 may be solved numerically for certain chosen values of renormalized gauge parameter @xmath100 and the results plot as functions of parameter @xmath101 . let us consider the case of weak coupling , when @xmath101 is small . since the applicability of our method requires @xmath52 , then from equation  ( [ eq : kapeq ] ) we deduce that @xmath99 should be large . it may be justified by the elementary estimation given below . first we rewrite following expression in terms of beta function @xmath102 and its integral representation : @xmath103 and next use the inequalities valid for @xmath52 : @xmath104 now , from  ( [ eq : kapeq ] ) it becomes obvious that for small @xmath101 the left hand side must be large . this is in agreement with our expectations concerning the subsequent terms in the taylor expansion  ( [ eq : dexp ] ) . but for @xmath105 the first term in square brackets on the right hand side of  ( [ eq : jedeq ] ) dominates over all other , which may be easily verified . the expected positivity of @xmath50 requires then considering only gauges for which @xmath106 . this is not surprising , since @xmath50 is a strongly gauge dependent quantity . for our numerical calculations we have then chosen the values of the gauge parameters from that range .    in figure  [ fig : graph1 ] we show the dependence of parameters @xmath77 and @xmath50 on @xmath101 for exemplary value of @xmath107 ( solid lines ) . in figure  [ fig : graph2 ] we similarly plot the dependence of two other parameters : fermion mass renormalization constant @xmath88 and gauge field renormalization constant @xmath41 . it is nice to observe that we have @xmath108 and @xmath109 , as we expected . the dashed lines in these figures represent approximate solutions of the set ( [ eq : deleq])-([eq : kapeq ] ) , as described below .    in the weak coupling regime , assuming that we do not choose the gauge parameter @xmath100 approaching zero ( but still to be less than @xmath110 ) , the equations may be given the following approximated form : @xmath111\ ; , \label{eq : deleqa}\\ 1&\!\!\ ! = & \!\!\ ! \frac{\zeta}{2}\left[\frac{3}{4\beta } \left(\frac{1}{\lambda_r}-2\right)-\frac{12m^2}{\kappa^2}\right]\ ; , \label{eq : jedeqa}\\ z_3&\!\!\ ! = & \!\!\ ! 1-\frac{2\zeta}{3}\ ; , \label{eq : z3eqa}\\ \frac{\kappa^2}{m^2}&\!\!\ ! = & \!\!\ ! \frac{15 } { \zeta}\ ; . \label{eq : kapeqa}\end{aligned}\ ] ] the terms @xmath112 in the first two equations arise from the expansion of @xmath84 and @xmath96 for large @xmath40 , but in fact they may be omitted , since they are of order @xmath101 , as results from the last equation . as already told , the solutions of these simplified equations are drawn in figures  [ fig : graph1 ] and  [ fig : graph2 ] as dashed lines . on the plots for @xmath77 and @xmath88 these curves are not visible since they are almost identical with the full solutions , but also on the other two they do not deviate from the ` exact ' results too much . it is interesting to observe , if and how the values of the parameters depend on the gauge . this is shown in figures  [ fig : graph3 ] and  [ fig : graph4 ] which are performed for the following values from the range @xmath115 $ ] : @xmath116 . the landau gauge can not be used because it would lead to negative value of @xmath50 . the particular stress deserves the observation that the gauge dependence of the parameter @xmath77 is extremely weak . this means that the renormalized gauge boson propagator ( strictly speaking its transverse part ) is practically gauge independent , as it should be . please note that in other approximation schemes applied to ds equations in qed@xmath6 one obtains the gauge dependent value of the polarization scalar  @xcite . the dependence of @xmath41 on @xmath100 is relatively weak for small @xmath101 too . these results are worth noticing , since the full gauge independence should appear in the _ exact _ theory , and one ought not to expect too much from the approximated model , where infrared forms of green s functions are postulated in a simple form . gauge dependence of the ` physical ' fermion mass obtained in our work is relatively strong but it is a common feature of nonperturbative calculations in this model  @xcite . 99 r. d. pisarski , phys . * d 29 * , 2423(1984 ) . t. appelquist _ et al _ , phys . lett . * 55 * , 1715(1985 ) . t. appelquist _ et al _ , phys . rev . * d 33 * , 3704(1986 ) . w .- h . kye and j.k . kim , phys . rev . * d 50 * , 5398(1994 ) . kondo and p. maris , phys . lett . * 74 * , 18(1995 ) ; phys . d 52 * , 1212(1995 ) . a. bashir , a. huet and a. raya , phys d 66 * , 025029(2002 ) . a. bashir and a. raya , nucl * b 709 * , 307(2005 ) . a. bashir _ et al . _ , phys . rev . * c 78 * , 055201(2008 ) . a. bashir _ et al . _ , few - body syst . * 46 * , 229(2009 ) . a. bashir , a. raya and s. sanchez - madrigal , phys . * d 84 * , 036013(2011 ) . t. w. allen and c. j. burden , phys . * d 53 * , 5842(1996 ) , _ ibid . _ * d 54 * 6567(1996 ) . t. matsuyama and h. nagahiro , univ . educ . bf 50 , 1(2001 ) . y. hoshino , arxiv:0706.1603 . m. gopfert and g. mack , commun . phys . * 82 * , 545(1982 ) . c. j. burden and c. d. roberts phys d 44 * , 540(1991 ) . c. j. burden , j. praschifka and c. d. roberts , phys d 46 * 2695(1992 ) . p. maris , phys . d 52 * , 6087(1995 ) . v. sauli , acta phys . * 2 * , 443(2009 ) . a. raya , arxiv : hep - th/0404138 . a. bashir and m. de j. anguiano galicia , few - body syst . * 0 * , 1(2008 ) . a. bashir and a. raya , few - body syst . * 41 * , 185(2007 ) . a. bashir , a. raya and s. sanchez - madrigal , j. phys . * a 41 * , 505401(2008 ) . j. s. ball and t .- w . chiu , phys . * d 22 * , 2542(1980 ) . r. delbourgo , r. zhang , j. phys * a 17 * , 3593(1984 ) . c. n. parker , j. phys . * a 17 * , 2873(1984 ) . d. c. curtis , m. r. pennington , phys . rev . * d 42 * , 4165(1990 ) . a. bashir and a. raya , phys . rev . * d 64 * , 105001(2001 ) . t. radoycki , arxiv:1212.6197 . t. radoycki and i. biaynicki - birula , phys . rev . * d 52 * , 2439(1995 ) . a. salam , phys . rev . * 130 * , 1287(1963 ) . for instance c. itzykson and j .- b . zuber , _ quantum field theory _ , mcgraw - hill , new york 1980 . r. delbourgo , n. cim . * 49 a * , 484(1979 ) . g. thompson , r. zhang , phys . rev . * d 35 * , 631(1987 ) . s. weinberg , _ the quantum theory of fields _ , cambridge university press , new york 1995 . p. maris , phys . rev . * d 50 * , 4189(1994 ) . k. stam , j. phys . * 9 * , l229(1983 ) . t. radoycki , eur . j. * c 6 * , 549(1999 ) .
Abstract:

Answer: quantum electrodynamics in three dimensions in the bispinor formulation is considered . it is shown that the dyson - schwinger equations for fermion and boson propagators may be self - consistently solved in the infrared domain if on uses the salam s vertex function . the parameters defining the behavior of the propagators are found numerically for different values of coupling constant and gauge parameter . for weak coupling the approximated analytical solutions are obtained . the renormalized gauge boson propagator ( transverse part ) is shown in the infrared domain to be practically gauge independent .


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: xmm - newton galactic plane survey ( xgps ) has revealed a large population of x - ray sources ( hands et al . 2004 ) . in order to identify the nature of the brightest sources detected in xgps , an optical campaign has recently been carried out ( motch et al . 2010 ) . in this process , three sources were identified as promising cataclysmic variables ( cvs ) .    the x - ray source designated as xgps-9 ( = xgps - i j183251 - 100106 ) in motch et al . ( 2010 ) is one of the newly found cv candidates in xgps . two optical objects are found within its x - ray error circle , while the fainter one ( @xmath2 ) is suggested to be the promising counterpart . its cv nature was identified through optical spectroscopy which unambiguously shows the strong h , he  i and he  ii emission along with a blue continuum ( motch et al . these features are known to be typical for a magnetic cv .    apart from the optical identification , motch et al . ( 2010 ) have also reported a brief x - ray analysis of archival xmm - newton data obtained by the observations on 15 march 2002 ( obs . ids 0135741601 , 0135744401 ) . the authors found that an absorbed single component thermal plasma model is able to describe the data . however , as xgps-9 located at an off - axis angle @xmath3 in these short exposures ( @xmath4  ks and @xmath5  ks respectively ) , the poor statistic and the degraded angular resolution leave the spectral parameters unconstrained . while the x - ray light curve of xgps-9 ( fig .  8 in motch 2010 ) clearly shows the variability of this source , these short observations preclude the search for any periodic behavior of this system .    in order to tightly constrain its x - ray properties , we have performed a follow - up investigation of xgps-9 with a deep xmm - newton observation . the results of this analysis are presented in this paper . xgps-9 was observed by the * e*uropean * p*hoton * i*maging * c*amera ( epic ) on board xmm - newton , which consists of two metal oxide semiconductor ccd detectors ( mos1 and mos2 ) and a pn ccd detector , on 8 - 9 october 2009 ( obs . i d : 0605480101 ) . while the mos1/2 cameras were operated in full frame mode with a a temporal resolution of 2.6  s , the pn camera was operated in extended full frame mode with a temporal resolution of 199.1  ms . medium filters were used to block optical stray light in all epic cameras . the aimpoint of this observation is ra=@xmath6 dec=@xmath7 ( j2000 ) . using tasks _ emproc _ and _ epproc _ of the xmm science analysis software ( xmmsas version 11.0.0 ) , we have reprocessed all the epic data with the updated instrumental calibration . we subsequently selected only those events for which the pattern was between @xmath8 for both mos cameras and @xmath9 for the pn camera in @xmath10 kev . we further cleaned the data by accepting only the good times when sky background was low for the whole camera ( @xmath11  counts / s , @xmath12  counts / s and @xmath13  counts / s for mos1 , mos2 and pn respectively . ) . after removing all events potentially contaminated by bad pixels , the effective exposures are found to be 55.9  ks , 56.6  ks and 41.3  ks for mos1 , mos2 and pn respectively . we have also utilized the optical monitor ( om ) data for our investigation . xgps-9 has been observed by om in standard imaging mode with three different filters : u ( effective wavelength 3440 ) , uvw1 ( 2910 ) and uvm2 ( 2310 ) . the exposures of om with u , uvw1 and uvm2 are 4.4  ks , 4.4  ks and 2.5  ks respectively . for data reprocessing , flat - fielding and source detection , we have utilized the metatask _ omichain_.    this observation was originally intended for spectro - imaging investigation of the supernova remnant kes  69 ( see figure  [ xmm_img ] ) , which has its x - ray properties poorly constrained ( cf . yusef - zadeh et al . 2003 , bocchino et al . 2012 ) . a detailed analysis of the remnant emission from kes  69 will be published elsewhere ( seo et al . in preparation ) as this is out of scope in this paper . in this observation , xgps-9 is an serendipitous source located at an off - axis angle of @xmath14 ( figure  [ xmm_img ] ) . with the aid of the xmmsas task _ edetect_chain _ , we have run source detection individually on each epic data set . the mean x - ray position of xgps-9 determined from this observation is ra=@xmath15 dec=@xmath16 ( j2000 ) with a resultant uncertainty of @xmath17 by combining the statistical errors inferred from each camera in quadrature . this x - ray position is marked as a black cross in figure  [ xmm_img ] .    for timing analysis , the point source was extracted within a circular region of a @xmath18 radius centered at the x - ray position from each camera , which corresponds to an encircled energy fraction of @xmath19 . 758 counts , 682 counts and 1277 counts were obtained from mos1 , mos2 and pn respectively . all the photon arrival times are subsequently corrected to the solar system barycenter with the xmmsas tool _ barycen _ by adopting the updated ephemeris jpl de405  m error in earth s position . this corresponds to a timing error of @xmath20  @xmath21s which is negligible in our investigation of orbital period . ] and the aforementioned mean x - ray position . in order to search for any periodicity , we have applied the techniques of both epoch - folding and the lomb - scargle periodogram ( scargle 1982 ) on the event list and the binned light curve respectively . the epoch - folding method was performed with the _ efsearch _ subpackage of heasoft . the peak value of the @xmath22 spectrum with 50 bins for the arrival time of merged mos and pn events is at @xmath23s . the uncertainty of epoch - folding is estimated according to the empirical formula of leahy ( 1987 ) . on the other hand , we have binned all of the events into a light curve with 150  s resolution to calculate the lomb - scargle periodogram . the periodogram demonstrates a periodicity at @xmath24s which is consistent with the result of epoch - folding . the uncertainty was estimated by @xmath25 times of monte carlo simulation . the results of the periodicity search with epoch - folding and lomb - scargle method are shown in fig . [ efsearch_ls ] .    in order to identify the nature of this periodicity , we have also folded the photon events based on the period we detected . the 5337.7s period was adopted to fold the light curve as it is directly resulted from epoch - folding . in addition , we have divided the photons into soft band ( 0.3 - 2.5  kev ) and hard band ( 2.5 - 12  kev ) for computing the hardness ratio : ( hard - soft)/(hard+soft ) . the folded light curves of both bands and the hardness ratio are shown in fig . [ efold_hardness ] , where the epoch zero is set at mjd  55113.0383 . the minimum at phase zero is particularly interesting and leads us to further investigate it . this feauture resembles those resulted from the eclipses of x - ray emitting regions in polars ( e.g. mukai et al . 2003 ; vogel et al . we notice there is residual x - ray emission at the minimum . to investigate if the residual is caused by the background , we estimated the count rates within a number of source - free circular regions of 15 " radius around xgps-9 . the average background count rate estimated from the merged dataset is @xmath26  counts / s which is consistent with the minimum of the folded light curve in 0.3 - 12  kev . this suggests the residual x - rays are likely the background contamination . therefore , the feature can possibly be a total eclipse of the x - ray emitting region by the companion star .    to further characterize this feature , we attempted to estimate the eclipse width and to obtain an ephemeris . in view of the relatively low statistics , these properties can not be properly constrained by fitting the eclipse profile . instead , visual inspection was adopted for the measurements reported in this paper . the folded light curves shown in figure  [ efold_hardness ] suggests the ingress and egress are rather sharp which allow us to estimate the eclipse duration with the flat bottom minimum . this leads to an approximate width of @xmath27  s , where the uncertainty is estimated from half of the bin size in figure  [ efold_hardness ] .    for defining the mid - eclipse phase , we adopted the mid - point of this feature in figure  [ efold_hardness ] . the timing of the x - ray eclipse can be described by the following linear ephemeris :    @xmath28    we have also checked these results by adopting different bin size ( 32 , 64 , 128 bin / period ) , we found that both the eclipse width and the ephemeris are consistent with the aforementioned results within the tolerence of uncertainties . to investigate if the periodic signal comes from the synchronous rotation , we further examined the existence of any marginal periodic signal in a larger frequency range using the @xmath22-test with 32 bins and fourier resolution . because the spin period of a typical white dwarf is about hundreds of seconds ( cf . ritter & kolb 2003 ) and the signal of orbital modulation seriously contaminates the periodogram when the searching range is larger than 2.5@xmath29  s@xmath30 , we only considered a periodic signal of spin within 100400  s. within this selected range , the largest peak value can be found at 164.7(1 )  s with the chance probability of @xmath31 . however , this result is not significant if we took into account the number of trials in the searching range . we also tried the decomposition methods ( e.g. the empirical mode decomposition proposed by huang et al . 1998 ) to resolve another strong signal . however , we did not yield any other possible periodicity except for the binary modulation . therefore , we conclude that no spin period can be detected from the current _ xmm - newton _ observation and it may give a hint of the synchronized spin and orbital periods . motch et al . ( 2010 ) have briefly reported the x - ray spectral properties of xgps-9 . they found that a one - component thermal plasma spectrum with a column absorption of @xmath32  @xmath33 can provide a reasonable fit to the data . however , due to the limited photon statistic , the plasma temperature and the x - ray flux in their study remains to be unconstrained . with a much longer integration time of the data adopted in our investigation , we are able to provide a tighter constraint on its spectral properties . for the spectral analysis , we have extracted the spectra within a circular region with a @xmath34 radius in each mos camera and a @xmath35 radius in pn camera around the x - ray position respectively . this corresponds to an encircled energy fraction of @xmath36 in all epic cameras . spectra extracted from each camera were grouped so as to have at least 50 counts per spectral bin . the background spectra were sampled from the low count circular region of @xmath34 radius centered at ra=@xmath37 , dec=@xmath38 ( j2000 ) in each camera . response files were computed by using the xmmsas tasks _ rmfgen _ and _ arfgen_. for the spectral analysis , we used xspec ( version 12.6.0 ) with @xmath39 statistics adopted for all the fittings . all the quoted uncertainties of the spectral parameters are @xmath40 for 1 parameter of interest .    based on the optical properties , motch et al . ( 2010 ) suggested xgps-9 a magnetic cv . to explore this scenario , we examine its x - ray spectrum by fitting an absorbed thermal bremsstrahlung model , which accounts for the x - ray emission from the shock - heated gas in the accretion column . we found that this single component model can describe the data reasonably well with a goodness - of - fit of @xmath39=49.40 for 54 d.o.f .. the spectral fit yields a column density of @xmath41  @xmath33 , a plasma temperature of @xmath42  kev and an unabsorbed flux of @xmath43  erg  @xmath33  s@xmath30 in the energy range of @xmath44  kev . the comparison between the best - fit model and the data is shown in figure  [ xmm_spec ] . we have also computed the confidence contours to investigate the relative parameter dependence between the plasma temperature and the column absorption , which are plotted in figure  [ contour ] .    as x - ray spectrum of cvs can possibly show the presence of line emission , we have also examined the observed spectrum of xgps-9 with xspec model mekal which is a code that models the plasma in collisional ionization equilibrium with the line emission built in . we found that the fit yields a comparable goodness - of - fit with the thermal bremsstrahlung model ( @xmath39=50.11 for 54 d.o.f . ) . the best - fit spectral parameters are also similar ( @xmath45  @xmath33 , @xmath46  kev ) .    for magnetic cvs , part of the thermal bremsstrahlung emission can possibly be compton reflected by the white dwarf surface ( see matt et al . this so - called compton reflection component is expected to be accompanied with a fluorescent k@xmath47 iron line ( de martino et al.2008 ) . for searching evidence for the line feature , we have added a gaussian component to the thermal bremsstrahlung model and fixed the line energy at 6.4  kev . to begin with , we have allowed both normalization and the line width as free parameters . however , we found that the line width can not be constrained ( i.e. the fitted width is essentially zero ) . this led us to fix it at a reasonable value of @xmath48  ev ( see balman 2011 ; de martino et al . this results in a best - fit parameter set of @xmath49  @xmath33 , @xmath50  kev , and a line flux of @xmath51  photons  @xmath33  s@xmath30 with @xmath39=49.11 for 53 d.o.f .. in view of the large uncertainty of fitted line flux and the insignificant improvement of the goodness - of - fit , we conclude that there is no compelling evidence of fe  k@xmath47 line can be found in this observation . we report a @xmath40 upper limit on the line flux of @xmath52  photons  @xmath33  s@xmath30 .    while a thermal bremsstrahlung / mekal can model the spectrum well , as mentioned by motch et al . ( 2010 ) , the spectral hardening can possibly be resulted from the reflection component and/or warm absorption ( see schwarz et al . 2009 ; staude et al . 2008 ) . therefore , the plasma temperature inferred from this simple model can be unphysically high . we have also investigated if any other single component model can provide any acceptable phenomenological description . we notice that a simple absorbed power - law can also result in a comparable goodness - of - fit ( @xmath39=49.83 for 54 d.o.f . ) , and hence we can not distinguish it from the thermal bremsstrahlung / mekal model statistically . the power - law spectral fit yields @xmath53  @xmath33 , a photon index of @xmath54 and an unabsorbed flux of @xmath55  erg  @xmath33  s@xmath30 ( @xmath44  kev ) . the spectral steepness and the flux are fully consistent with those inferred from a brief analysis limited to the hard band ( i.e. @xmath56  kev ) as reported recently ( cf . tab .  1 in bocchino et al . 2012 ) . we have checked the robustness of all the spectral results quoted in this paper by incorporating background spectra sampled from different source - free regions around xgps-9 . we found that the spectral parameters inferred from the fittings with different adopted background are consistent within @xmath40 uncertainties .    we have also attempted to constrain the possible contribution from the soft component by using om data . soft component from a cv can be originated from the optically - thick blackbody like emission of the white dwarf surface . the temperature can range from few tens to hundred of ev ( evans & hellier 2007 ) . searching for the optical / uv counterpart in om data results in non - detection in all three bands , which place a limiting flux density of @xmath57jy , @xmath58jy and @xmath59jy in u , uvw1 and uvm2 respectively . we adopted the column density inferred from the x - ray spectral fitting ( i.e. @xmath60  @xmath33 ) to perform the extinction - correction ( cf . predehl & schmitt 1995 ; cardelli et al . 1989 ) . with the de - reddened limiting flux densities , we place an upper bound for any soft component contribution by assuming a low temperature white dwarf with @xmath61  ev ( see figure  [ om_sed ] ) . the blackbody gives a limiting flux of @xmath62  erg  @xmath33  s@xmath30 in @xmath44  kev , which is @xmath63 of the total x - ray flux in this energy range . utilizing the archival _ xmm - newton _ data , we have detected the periodicity from a new cv candidate xgps-9 . through a detailed timing analysis , only a single periodic signal of @xmath0  hrs can be found . this suggests a possible scenario of synchronous rotation , which favors the interpretation that the system is a magnetic cv of polar type .    for a polar , the white dwarf is accreting matter from its late - type main sequence companion via roche - lobe overflow which will be driven by the strong magnetic field of the white dwarf ( with a surface magnetic field strength of few tens of million gauss ) to form a quasi - radial flow towards the magnetic pole . the accretion flow will eventually form a strong shock above the stellar surface . the shock - heated plasma will be cooled via bremsstrahlung radiation in x - ray regime with typical temperature of few tens of kev , which is consistent with the best - fit plasma temperature ( i.e. @xmath42  kev ) inferred in our spectral analysis . it is interesting to compare the temporal behavior of xgps-9 and 2xmmp  j131223.4 + 173659 , which is a eclipsing polar serendipitously identified in 2xmm catalog ( vogel et al . the properties of these two systems are remarkably similar . first , both objects show a narrow and deep minimum in their phase - folded light curves ( i.e. depression in phase @xmath64 for the case of xgps-9 ) . this feature is interpreted as the eclipse of x - ray emission in the polar region of the white dwarf when the companion star passes through our line - of - sight towards this region . also , the broad depression ( i.e. phase @xmath65 in fig .  [ efold_hardness ] ) , which is likely due to the self - occultation by the white dwarf itself , is observed in both system . the coexistence of these two features indicates that the spin of the white dwarf and the binary orbital periods are synchronized . furthermore , a shallow minimum in the soft band light curve ( in phase @xmath66 in fig . [ efold_hardness ] ) is found before the eclipse . when the x - ray emitting region is obscured by the accretion stream , a pre - eclipse dip can be resulted from the increased photoelectric absorption which makes the dip more prominent in the soft band ( see fig . [ efold_hardness ] ) . we further attempted to place a constraint on the orbital inclination with our estimated eclipse width ( i.e. @xmath67  s ) . assuming the companion star fills its roche lobe and x - ray emitting region is close to the surface of the white dwarf , the duration of the total eclipse of the white dwarf can be expressed as a function of mass ratio and the inclination ( cf . fig .  2 in horne 1985 ) . the eclipse width corresponds to a phase interval of @xmath68 which we considered as an upper limit of the eclipse width as a longer observation might reveal the falling / rising edges in the eclipse profile . we further assume that the companion star is a m - dwarf with mass @xmath69 . for a typical white dwarf ( i.e. @xmath70 ) , this implies a mass ratio of @xmath71 . with these assumptions , the upper bound of the orbital inclination can be placed at @xmath72 .    for further investigation of xgps-9 , dedicated optical / infrared observations will be useful for confirming its source nature . a differential photometric study can enable a comparison with the x - ray light curve reported in this paper . this can provide information for further constraining the accretion geometry . also , the free electrons in the shock region can spiral in the magnetized ionized gas . this can lead to the emission of cyclotron radiation in infrared regime . therefore , infrared spectroscopy and polarimetry are encouraged for determining the magnetic field of the white dwarf directly . we would also like to point out that the current data can not provide a conclusive constraint on the plasma temperature of a thermal bremsstrahlung spectrum at a few tens of kev , which is far above the energy coverage of xmm - newton . this is reflected by its large relative uncertainty reported in this work ( i.e. @xmath73 ) . for a more constraining spectral results , hard x - ray observation is certainly needed . as xgps-9 resides in rather complex region , the recently commenced mission nustar and the upcoming one ( e.g. astro - h ) , which equipped with hard x - ray focusing optics , will provide the desirable instruments for the further investigations of this interesting eclipsing polar . balman , s. 2011 , arxiv:1112.3893 bocchino , f. , et al . 2012 , a&a , 541 , 152 cardelli , j. a. , clayton , g. c. , & mathis , j. s. 1989 , apj , 345 , 245 de martino , d. , et al . 2008 , mmsai , 79 , 246 evans , p. a. & hellier , c. 2007 , apj , 663 , 1277 hands , a. d. p. , warwick , r. s. , watson , m. g. , & helfand , d. j. 2004 , mnras , 351 , 31 horne , k. 1985 , mnras , 213 , 129 huang , n. e. et al . lond . a , 454 , 903 leahy , d. a. 1987 , a&a , 180 , 275 matt , g. , et al . 1998 , a&a , 334 , 857 motch , c. et al . 2010 , a&a , 523 , a92 mukai , k. , hellier , c. , madejski , g. , patterson , j. , & skillman , d. r. 2003 , apj , 597 , 479 predehl , p. , & schmitt , j. h. m. m. 1995 , a&a , 293 , 889 ritter , h. , & kolb , u. 2003 , a&a , 404 , 301 scargle , j. d. 1982 , apj , 263 , 835 schwarz , r. et al . 2009 , a&a , 496 , 833 staude , a. et al . 2008 , a&a , 486 , 899 vogel , j. et al . 2008 , a&a , 485 , 787 yusef - zadeh , f. , wardle , m. , rho , j. , & sakano , m. 2003 , apj , 585 , 319
Abstract:

Answer: we report the results from a detailed analysis of an archival xmm - newton observation of the x - ray source xgps - i j183251 - 100106 , which has been suggested as a promising magnetic cataclysmic variable candidate based on its optical properties . a single periodic signal of @xmath0  hrs is detected from all epic cameras on board xmm - newton . the phase - averaged x - ray spectrum can be well - modeled with a thermal bremsstrahlung of a temperature @xmath1  kev . both x - ray spectral and temporal behavior of this system suggest it as a eclipsing cataclysmic variable of am herculis ( or polar ) type .