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 development of increasingly accurate atomic clocks has led to many advances in technology and tests of fundamental physics . in the search for the next generation of clocks and frequency standards , there has been considerable interest in using alkaline earth species because of their narrow intercombination lines in the optical spectrum  @xcite . in order to achieve a high level of short - term stability and long - term reproducibility and accuracy on the clock transition , it is desirable to have a large number of cold atoms located in a well - characterized trap for an improved signal - to - noise ratio ( @xmath0 ) and reduced systematic errors associated with atomic motion . single ion - based systems do effectively eliminate doppler and other motion - related systematic errors when the single ions are confined in the lamb - dicke regime  @xcite , although the achievable @xmath0 is limited by single quantum absorbers . for neutral atoms it is important that changes in the level structure due to the trapping potential do not alter the relevant clock transition frequency . such a scheme has been proposed by trapping alkaline earth atoms in three - dimensional optical lattices tuned to a  magic " wavelength where the relevant states for the clock transition experience exactly the same level shift  @xcite . the @xmath1 forbidden transition ( @xmath2 nm ) in @xmath3sr  @xcite is in particular a promising candidate for a lattice - based optical clock transition because of the long lifetime of the excited state ( @xmath4 s ) and the insensitivity of the @xmath5 states to the polarization state of the trapping light . already there have been efforts towards the cooling and trapping of @xmath3sr  @xcite , and recently this transition was directly observed and measured for the first time  @xcite . calcium , another alkaline earth atom that has been studied extensively as a frequency standard  @xcite , may be a candidate for optical lattice clocks as well .    in the case of @xmath6 independent atoms , one benefits from a @xmath7 improvement in @xmath0 in spectroscopy . however , atoms trapped in an optical lattice can interact with each other and can not truly be considered independent . each optically excited atom represents essentially a point dipole whose radiated electromagnetic field can affect other atoms . these atom - atom interactions can manifest themselves as shifts in the observed transition frequencies . because of the spatial ordering of atoms in a lattice and the potentially high atomic density , it is possible that such interactions may produce very large frequency shifts . one might expect then that dipole - dipole interactions can be much more severe here than in , for example , atomic fountains , and thus could place serious limits on the accuracy of an optical lattice clock if not accounted for . on the other hand , it might be possible to design lattice geometries where this shift is reduced or cancelled . although the trapping lasers are constrained to operate at the  magic " wavelength , the lattice geometry can be altered by changing the relative orientations of the trapping beams , whose degrees of freedom can be characterized by a set of variables @xmath8 .    in this paper , we investigate theoretically the dipole - dipole interaction - induced shifts in the clock transition frequency recovered by ramsey spectroscopy . we show that by varying the lattice geometry we can quantitatively control the clock frequency shift and even reduce the shift to zero . in particular , we give an analytical equation that can be solved giving configurations @xmath9 where constructive interference causes the line shift to be very large . in these  bad " lattice configurations , the magnitude of the line shift scales approximately like @xmath10 . quite generally we propose that by tuning the parameter space @xmath8 to lie in between two of these bad configurations , one can find  good " configurations where the shift is cancelled . the mechanism of cancellation is associated with the destructive interference of contributions to the shift from different atoms in the lattice . it is important to emphasize that the present mechanism and theoretical treatment differ considerably from the conventional approaches used to treat dipole - dipole line broadening and shifts . in the case of atoms in a spatially - ordered lattice geometry , long - range effects are important , and the usual methods involving binary collisions of nearest neighbors  @xcite are not applicable . these long - range effects include , in particular , interference of the far - field dipole radiation produced by the excited atoms , a phenomenon similar to bragg scattering in a crystal . this paper is organized as follows . in sec . [ sec : eqofmotion ] we derive equations describing the evolution of an atomic system with dipole - dipole interactions . these equations are derived assuming that the atoms are in the lamb - dicke regime , with one atom or less per lattice site . in sec .  [ sec : ramsey ] we give a brief review of ramsey spectroscopy and solve for the dipole - dipole induced line shift using perturbation theory . we find that the shift can be qualitatively understood in terms of the classical interaction energies between oscillating dipoles . there is a contribution to the shift that is zeroth order in the interrogation time @xmath11 , which is due to imperfections in the ramsey pulses . even with perfect pulses , however , one finds a shift that is first order in @xmath11 that results from spontaneous decay of the atoms . [ sec : generalization ] discusses how our result for the line shift can be generalized for systems with imperfect filling of the lattice sites and for multilevel atoms . in the case of imperfect filling , one can calculate the mean value of the frequency shift as well as some nonzero variance , due to the uncertainty of how the lattice is filled . in sec . [ sec : latticedesign ] , we derive an equation that can be solved giving lattice configurations where the shift is large due to constructive interference . we derive an approximate scaling law for the shift in these  bad " configurations and discuss how the line shift can be reduced by choosing an appropriate lattice design . in sec . [ sec : example ] we demonstrate these results numerically for one specific lattice configuration . to treat the problem of interacting atoms in a lattice , we consider @xmath6 two - level atoms in the lamb - dicke limit with polarizability along the @xmath12-axis . a simple model of the system consists of treating the atoms as point dipoles , and we further assume that there is one or less atom per lattice site . this corresponds to a mott - insulator state for bosons or a normal state for fermions . in principle , to solve exactly the problem of interacting atoms one would start from the full atom - field hamiltonian and take into account not only all the atomic degrees of freedom but the continuum of electromagnetic field modes . to simplify the theoretical treatment , we effectively eliminate the field in the standard way using the born - markov approximation  ( see appendix ) . this is valid provided that the atomic system evolves slowly on timescales of the correlation time @xmath13 , which is of the order @xmath14 where @xmath15 is the linear size of the system and @xmath16 is the speed of light . as a result of eliminating the field , one finds an effective equation of motion for the density matrix @xmath17 of the atomic system . atom - atom interactions then appear through an effective hamiltonian @xmath18 as well as through a non - hermitian operator @xmath19 : = + [ ] . here , @xmath20 is the atomic hamiltonian for a non - interacting system . writing out all the terms in detail , & = & ^_a -^_abg(k_ab ) + & & [ eq : dpdt]-^_a , bf(k_ab)(\{,}-2 ) -^_a(- ) , where f ( ) & = & ( ^2+(3 ^ 2 - 1)(- ) ) , + [ eq : fg]g ( ) & = & -(^2+(3 ^ 2 - 1)(+ ) ) , and @xmath21 is the angle that @xmath22 makes with the @xmath12-axis .    the first term on the right - hand side of eq .  ( [ eq : dpdt ] ) corresponds to @xmath20 . @xmath23 is the pauli matrix of atom @xmath24 corresponding to the population difference between the excited and ground states , and @xmath25 is the resonance frequency of the dipole transition . the second term corresponds to @xmath18 . here , @xmath26 is the spontaneous decay rate of the excited state of a single , isolated atom , where @xmath27 and @xmath28 is the dipole matrix element between the ground and excited states . @xmath29 is the atomic raising operator on atom @xmath24 , and @xmath30 is the lowering operator on atom @xmath31 . one then sees that the effect of dipole - dipole interactions is an exchange of excitation between pairs of atoms . the strength of interaction is modified by a function @xmath32 that depends on the distance and orientation between two dipoles . it is to be understood that @xmath33 in the functions @xmath34 and @xmath35 . we see that both short - range , near - field  @xmath36 and long - range , far - field  @xmath37 dipole interactions are included in our formalism and are treated on equal footing . the third term on the right side of eq . ( [ eq : dpdt ] ) corresponds to @xmath19 and also is due to atom - atom interactions . it also depends on @xmath38 and has a position dependence described by @xmath39 . the non - hermitian nature of this term is evident through the anticommutator . physically @xmath19 describes the processes of both independent and cooperative decay . finally , we have also added phenomenological dephasing through the @xmath40 term , which in particular includes the effects of a finite , short - term linewidth of the laser interrogating the clock transition .    from eq . ( [ eq : dpdt ] ) , one can derive equations of motion for any atomic operators . for our particular application of ramsey spectroscopy , we find it necessary to solve for the coherence @xmath41 and the two - atom correlation @xmath42 . furthermore , to remove the rapid oscillations due to @xmath25 , it is convenient to work in the frame rotating with the interrogating laser frequency @xmath43 , where & = & -i- + ^_ba(f(k_ab)-ig(k_ab))[eq : ddtpa ] , + & = & -(i+)- -(f(k_ab)+ig(k_ab ) ) -f(k_ab ) + & & + ^_ja , b(f(k_bj)-ig(k_bj ) ) -^_ja , b(f(k_aj)+ig(k_aj ) ) + & & -^_ja , b(f(k_aj)-ig(k_aj))[eq : ddtzapb ] , and @xmath44 . in principle , to solve for the atomic system exactly , equations of motion for higher - order correlations are needed , although typically some approximation is used to truncate the resulting hierarchy of equations . in general these equations can describe light - matter interactions in an optically dense medium , including radiation trapping , level shifts , and superradiance  @xcite . we now analyze the effects of atom - atom interactions on ramsey spectroscopy . starting from the ground state @xmath45 of the system , suppose that one applies a strong probe pulse with the interrogating laser , given in the rotating frame by the hamiltonian h=^_ai(e^i_a - e^-i_a ) , where @xmath46 is the rabi frequency . for simplicity , we have made a plane - wave assumption about the probing laser , taking @xmath47 to be in the positive @xmath48 direction . we also assume that @xmath49 , and suppress the subscript in future calculations . we can do this if the phase error @xmath50 over the length of the sample incurred by making this assumption is small . applying this pulse for a time @xmath51 evolves the system through the unitary operator u=^_a^ (    cc & e^ikx_a + -e^-ikx_a &    ) . the state vector immediately following this pulse is given by = ^_a^(+e^ikx_a ) , where @xmath52 denotes the excited state of the clock transition .    in ramsey spectroscopy , one lets @xmath53 evolve for an interrogation time @xmath11 to @xmath54 . during this time the coherence between the ground and excited states acquires some time - dependent phase that depends on @xmath55 . after a time @xmath11 , one then applies a second pulse corresponding to the inverse unitary operation @xmath56 , and then measures the signal corresponding to @xmath57 , averaged over the final system @xmath58 . @xmath59 corresponds to the total population inversion . because of the second pulse , @xmath59 will now depend on @xmath55 and @xmath11 , allowing one to extract information about the resonance line . formally , we can rewrite @xmath59 as & = & ( ^_au^(t)u ) + & = & ^_a -2^_ae^ikx_a[eq : signaltilde ] , where the averages denoted above apply to @xmath54 , the system immediately before the second pulse .    in the case of @xmath6 non - interacting , independently decaying atoms , & = & -1+e^-t(1- ) , + e^ikx_a & = & e^-(i+/2+/2)t , which gives a corresponding signal [ eq : independentsignal]=-n((1-e^-t)+e^-t^22+e^-(+)t/2 ^ 22 ) . one can see that there is a peak in the signal around @xmath60 . determination of this peak allows one to find the frequency of the transition . one can note two important points about @xmath59 . the contrast in @xmath59 with respect to @xmath55 is maximized when a  perfect " @xmath61 pulse is applied , _ i.e. _ , when @xmath62 . furthermore , the contribution to @xmath59 due to @xmath63 in eq . ( [ eq : signaltilde ] ) is independent of @xmath55 and thus plays no role in determination of the resonance line . thus , one is motivated to define an effective signal @xmath64 that consists of the part of @xmath59 that is actually used to determine the line : [ eq : signal ] s=-2^_ae^ikx_a . the equation above states that from a theoretical standpoint , determination of the resonance line by measuring the population inversion after the second ramsey pulse is equivalent to measuring the real part of @xmath65 directly before the second pulse .      solving for @xmath64 exactly in the presence of dipole - dipole interactions appears to be quite a difficult task . since all interactions are proportional to @xmath38 , our approach is to solve for @xmath64 as a perturbative expansion in @xmath66 . in particular , we solve for the coherence @xmath41 in eq . ( [ eq : ddtpa ] ) to second order in @xmath38 . this requires solving @xmath42 in eq . ( [ eq : ddtzapb ] ) to first order in @xmath38 . the solution for the coherence with appropriate initial conditions is & = & e^-ikx_a(2)e^-(i+/2)t[eq : pafull ] , where a_ab & = & -e^ik(x_a - x_b)(2)+e^ik(x_a - x_b)+(f(k_ab)+ig(k_ab ) ) -f(k_ab)(2 ) + & & -^_ja , b(f(k_bj)-ig(k_bj))e^ik(x_a - x_j)^2(2 ) , + b_ab & = & e^ik(x_a - x_b)^2(2)^_ja , b(f(k_aj)+g(k_aj ) ) , + c_a & = & ^_ba(f(k_ab)-ig(k_ab))e^ik(x_a - x_b ) . ( [ eq : pafull ] ) is correct to every order of @xmath40 . ( [ eq : signal ] ) and ( [ eq : pafull ] ) can be evaluated numerically for a given lattice configuration and number of atoms . to illustrate the general features of the shift , however , we now make the following simplifications . we expand eq . ( [ eq : pafull ] ) to lowest order in @xmath40 . we also assume that the ramsey pulses are nearly perfect @xmath61 pulses , _ i.e. _ , @xmath67 . we then keep terms like @xmath68 but ignore terms like @xmath69 with these simplifications , [ eq : pa ] e^-ikx_a(2)e^-it(1- + -_a ) , where _ a & = & ^_ba(f(k_ab)-ig(k_ab))e^ik(x_a - x_b ) ( ( 2)+ ( 1+(f(k_ab)+ig(k_ab))e^-ik(x_a - x_b ) .. + & & + .. ^_ja , b(f(k_aj)+g(k_aj ) ) ) ) the first three terms in the parentheses of eq . ( [ eq : pa ] ) are a result of expanding the @xmath70 term that appears in the result for independent atoms , given in eq . ( [ eq : independentsignal ] ) . this is just the decay of the signal one would get from independent spontaneous emission . the last term in the parentheses is a correction due to atom - atom interactions . plugging this result into eq . ( [ eq : signal ] ) , we find that s-^22((t)(n- + -^_a_a ) -(t)^_a_a ) . because of the antisymmetric @xmath71 term now appearing in @xmath64 , one immediately sees that dipole - dipole interactions introduce a shift @xmath72 in the ramsey fringes , which can be found by solving @xmath73 . suppose that the inequalities @xmath74 are satisfied . under these conditions , a simple expression for @xmath72 results : & & ^_a^_ba(g(k_ab)-f(k_ab ) ) + & & ( + ( 1+^_ja , b(f(k_aj)+g(k_aj))))[eq : shift ] . the shift given by eq . ( [ eq : shift ] ) yields a simple interpretation . in anticipation of future analysis , we write @xmath72 as [ eq : shiftsymbolic ] = ^_a^_ba_ab(+_a ) , where _ ab & = & g(k_ab)-f(k_ab ) , + & = & , + _ a & = & 1+^_ja , b(f(k_aj)+g(k_aj ) ) . we will see that @xmath75 is a dimensionless quantity corresponding to the classical interaction energy between two oscillating dipoles , @xmath76 is a parameter characterizing the error in the ramsey pulses , and @xmath77 is a dimensionless quantity characterizing cooperative decay of the system .    to show the meaning of the @xmath75 term in eq . ( [ eq : shiftsymbolic ] ) , consider the interaction between a classical , oscillating dipole at @xmath78 excited with phase @xmath79 and the field incident on it due to a classical , oscillating dipole at @xmath80 excited with phase @xmath81 . we assume that both dipoles are oriented along the @xmath12-axis and that their magnitudes @xmath28 are determined from the relation @xmath82 . the classical interaction energy between dipole @xmath24 and the incident field is given by @xmath83 . the field at @xmath78 due to dipole @xmath31 is  @xcite : [ eq : classicalfield ] e_z(_a)=e^i(kx_b - t)e^ikr ( + ( 3 ^ 2 - 1)(- ) ) , where @xmath84 . now using the definitions in eq . ( [ eq : fg ] ) , the interaction energy can readily be rewritten as : u_ab & = & -(_a_b^*(_a))=(g(k_ab)-f(k_ab ) ) + & = & _ ab . one then sees that this indeed corresponds to the first term of eq . ( [ eq : shiftsymbolic ] ) .    although the @xmath75 term in eq . ( [ eq : shiftsymbolic ] ) resembles a classical interaction energy , the terms in parentheses reflect the quantum mechanical nature of the system . there is a contribution to the shift that is zeroth order in the interrogation time @xmath11 and proportional to @xmath76 . one notes that for perfect @xmath61 ramsey pulses , @xmath85 . thus , the zeroth order shift is due to error in the ramsey pulse this can be understood by considering eq . ( [ eq : ddtpa ] ) . one sees that the interaction terms only influence evolution of the coherence through the term @xmath42 . for a perfect @xmath61 pulse , this term is initially zero , and in this case , interactions can not affect the measurement at short times . this effect is due to the nature of dipole - dipole interactions : these interactions can not influence the coherence @xmath41 of an atom when it is in an equal superposition of the ground and excited states . even if a perfect @xmath61 pulse is applied , there is an additional contribution to the shift that is first order in @xmath11 and whose strength is given by @xmath77 . the intuition behind this is also straightforward . even if @xmath42 is initially zero , decay of the excited state will eventually evolve @xmath63 away from zero and back towards its equilibrium value of @xmath86 . once @xmath42 is nonzero , interactions can influence evolution of @xmath41 . the rate of decay is characterized by @xmath77 . the first term in @xmath77 is the contribution from independent decay of the atom back to the ground state . the second contribution involves a sum over other atoms and represents a correction due to the fact that the decay process may in fact be cooperative ( e.g. , superradiance ) . one can easily verify that the contribution from atom @xmath87 is proportional to @xmath88 : ( _ a_j^*(_a))f(k_aj)(kx_aj ) + g(k_aj)(kx_aj ) . this reflects the well - known result that the atomic inversion @xmath63 is driven by the dipole component in quadrature with the incident field . experimentally , knowing the exact number of atoms in the lattice and achieving a filling factor of one atom per lattice site are difficult tasks . most likely , one can experimentally determine the density @xmath89 of atoms in the lattice , such that the probability of occupation at any particular site @xmath24 is @xmath90 , where @xmath91 is the volume of a unit cell . it is straightforward to modify eq . ( [ eq : shift ] ) to the case of imperfect filling . for simplicity , we only consider the shift that is zeroth order in @xmath11 . this shift can be written & = & ^_a^_ba ( g(k_ab)-f(k_ab ) ) + & = & ^_g(k_ab ) + & = & ^_0u ( ) , where @xmath92 , @xmath93 denotes the set of direct lattice vectors , and @xmath94 is the number of pairs of atoms separated by @xmath95 . in the derivation above we have utilized the fact that @xmath96 to cancel the sum of @xmath97 . in a realistic scenario , one neither knows @xmath94 nor @xmath6 exactly . in this case , one must solve instead for the ensemble average @xmath98 and the variance @xmath99 . for large @xmath6 , one can safely pull the factor of @xmath6 out of the ensemble average : . with this simplification , & = & ^_0u()n ( ) + & = & ^_0u()d()(+)v[eq : avgshift ] , + & = & ( ) ^2^_,0u()u(^)n()n(^ ) + & & -(^_0u()n())^2 + & = & ( ) ^2[eq : shiftvariance ] . each term in eq . ( [ eq : shiftvariance ] ) has a clear meaning . to calculate the average shift in eq . ( [ eq : avgshift ] ) , one must evaluate @xmath100 . to do this , one must perform a sum over @xmath101 of the probability that the sites @xmath101 and @xmath102 are both occupied . to find the variance , one must calculate quantities like @xmath103 , and thus the probability that the sites @xmath101 , @xmath102 , @xmath104 , and @xmath105 are all occupied . when these four points are distinct , the probability is simply a product of the probabilities of each point being occupied . this is untrue when one or more of the points overlap . the terms in eq . ( [ eq : shiftvariance ] ) represent corrections due to these overlaps . the product @xmath106 , for example , is due to the overlap of @xmath101 and @xmath104 . our results derived thus far are for the case of two - level atoms . this is the relevant case of study for the @xmath107 forbidden transition proposed for optical lattice clocks , where the simple level structure makes it easier to cancel the relative stark shift in the clock transition . nonetheless , our results can be generalized to more complicated level structure , such as the case of an atom with a single ground and multiple excited states . a simple argument shows that eq . ( [ eq : shift ] ) remains correct to the lowest nontrivial order in @xmath108 . if multiple excited states are present in addition to the one that is initially excited , the equation of motion for @xmath109 in eq . ( [ eq : pa ] ) will contain additional terms like @xmath110 , where the subscript  clock " refers to the clock transition and  other " refers to other excited state levels . initially , @xmath111 and thus t. consequently , at short times , evolution of @xmath109 will be dominated by the clock transition . thus , if imperfections in the ramsey pulse constitute the major source of shift , the multiple excited states will contribute an additional source of shift that is first order in @xmath108 . if decay of the clock excited state constitutes the major source , the multiple excited states will contribute a shift that is of order @xmath112 . eq .  ( [ eq : shift ] ) or eqs .  ( [ eq : avgshift ] ) and  ( [ eq : shiftvariance ] ) can be evaluated numerically for a given lattice configuration and number of atoms . to extract the key features of the shift , we note that the zeroth - order shift in @xmath11 in eq . ( [ eq : shift ] ) essentially consists of adding together the classical dipole interaction energies @xmath113 . for a generic configuration of atoms , the amplitudes of the dipole fields incident on a given dipole @xmath24 tend to interfere . for certain configurations , it is possible that the field amplitudes will add constructively along some direction of propagation @xmath47 . near these configurations one will expect large shifts to result . the condition for constructive interference between radiated dipole fields is similar to that of bragg scattering in a crystal , and is readily found to occur when ||=k_0 , where @xmath114 and @xmath115 is a reciprocal lattice vector . this condition can be rewritten as [ eq : bigshift ] ||^2=2k_0g_x.numerical results indicate that peaks in the line shift do indeed occur when condition ( [ eq : bigshift ] ) is nearly satisfied .    one can easily derive an approximate scaling law for the line shift in these resonant configurations . we define a dimensionless parameter @xmath116 related to the density of atoms by @xmath117 . @xmath116 characterizes the spacing between neighbors in the lattice . in a resonant configuration , the electric fields add constructively , and the total electric field experienced by an atom is approximately e~d^3~_0^ldr~ , where @xmath15 is the linear size of the system . for @xmath6 total atoms , @xmath118 . then ~. experimentally , one has freedom to choose the orientations of the trapping laser beams that form the lattice . the control parameters can be parameterized by a set of variables @xmath8 , which will also determine the reciprocal lattice vectors @xmath119 . one can then find solutions @xmath9 of eq . ( [ eq : bigshift ] ) corresponding to configurations with large line shifts . in the parameter space between two sets of solutions @xmath9 , one can numerically find configurations where the shift is significantly reduced .    in the case of imperfect filling of lattice sites , it will be important to account for not only the mean shift but the variance as well . for large numbers of atoms , eq . ( [ eq : shiftvariance ] ) can not be evaluated exactly without extensive computational resources . with a small filling factor @xmath120 , however , we can estimate that the major contribution to the variance results from the @xmath121 terms , while the @xmath122 terms remain negligible . in this diffuse limit , [ eq : diffusevariance ] ( ) ^22^_0u()^2d()(+)v . in this case , one readily finds that the variance @xmath99 scales like @xmath123 . for @xmath124 , the variance increases with @xmath125 due to the increasing uncertainty of whether a pair of sites will both be occupied , but decreases with @xmath6 due to the decreasing fractional uncertainty in the total number of pairs of atoms separated by a vector @xmath95 . as an illustration of our results , we consider @xmath3sr atoms trapped in a lattice formed by six interfering beams , as shown in fig .  [ fig : sixbeam ] . for @xmath3sr , the  magic " wavelength of the trapping lasers is roughly @xmath126  @xcite , and one can vary the angle @xmath21 between the propagation vectors of the trapping beams . the resulting lattice is tetragonal , with lattice constants of @xmath127 , @xmath128 , and @xmath129 along @xmath48 , @xmath130 , and @xmath12 , respectively . the lattice constants are plotted in fig .  [ fig : latticeconstants ] . the corresponding basis of the reciprocal lattice has lengths @xmath131 , @xmath132 , and @xmath133 . we have ignored the effect of atomic back - action  @xcite on the trapping fields , whereby scattering of light by the atoms introduces phases that might modify the lattice constants . such effects are expected to be stronger in red - detuned lattices , where atoms lie in the antinodes of the potential , and with increasing atomic density . taking into account this back - action does not modify our results , except that now the lattice constants must be solved self - consistently  @xcite . using eq . ( [ eq : bigshift ] ) we can find values of @xmath21 where constructive interference causes the shifts to be large . we focus on two specific solutions , @xmath134 and @xmath135 . for our system we consider @xmath136 atoms in a spherical distribution with uniform density @xmath137 for @xmath138 , and zero density for @xmath139 . the relationship between the density and filling fraction @xmath125 is given by @xmath140 , where @xmath91 is the volume of a unit cell . the critical value @xmath141 is determined by the equation r_0=()^1/3 . we first consider a perfectly filled lattice consisting of @xmath142 atoms . for simplicity we calculate the line shift to zeroth order in the interrogation time @xmath11 . in fig . ( [ fig : perfect ] ) , we plot the quantity @xmath143 as a function of @xmath21 . peaks in the shift are clearly visible at the points @xmath144 that were calculated analytically . it should be noted that the line shift can be very large in one of these  bad " configurations . even in the limit of short interrogation times , one can see that shifts of order @xmath145 are possible . this is perhaps a surprising result , and occurs because the spatial ordering of the atoms allows the interactions to behave constructively at these points . for longer interrogation times , one expects this line shift to become even larger , since the constructive interference in these configurations also leads to superradiant decay and thus a large contribution to the shift that is first order in @xmath11 . one also sees that away from these  bad " points , the shift is strongly suppressed and even becomes zero for one particular value of @xmath21 . we next consider a partially filled lattice consisting of @xmath146 atoms and a filling factor of @xmath147 . [ fig : imperfect ] gives the quantity @xmath148 as a function of @xmath21 . the mean shift @xmath98 still exhibits peaks at the points @xmath144 , and vanishes near @xmath149 . at this  good " point , one can use eq . ( [ eq : diffusevariance ] ) to estimate the variance in the expected shift . within this diffuse approximation , we find that the variance ( 2)3.110 ^ -3 . experimentally , there will be additional sources of error that result from not knowing @xmath137 perfectly , errors in the configuration of the trapping lasers , and from the effects of atomic back - action on the lattice constants . nonetheless , it appears that the error due to dipole - dipole interactions can be made quite small by appropriately designing the lattice . we have derived an expression for the line shift measured in ramsey spectroscopy due to dipole - dipole interactions . we find that the lattice geometry strongly affects the magnitude of the shift , and is peaked in lattice configurations where the interactions between atoms add constructively . because of the spatial ordering in the lattice , the shift can be quite large in these resonant configurations . by tuning the lattice between two of these configurations , one can reduce the dipole - induced line shift to nearly zero .    while the resonant configurations might be bad for clock applications , it might be worthwhile to study these configurations further . the dipole - dipole couplings in an optical lattice offer the possibility of strong , constructive interactions that can be dynamically tuned by changing the lattice geometry . this might be useful for applications such as quantum information processing and might have interesting consequences for studying phenomena such as superradiance and for probing the superfluid - mott insulator transition . we gratefully acknowledge a. srensen , a. andr , and r. walsworth for many helpful discussions . this work was supported by the nsf ( career program and graduate fellowship ) , a. sloan foundation , the david and lucille packard foundation , nist , and onr . in this appendix we derive eq . ( [ eq : dpdt ] ) starting from the full atom - field hamiltonian . for a more detailed derivation and discussion , one can also see  @xcite . the hamiltonian for the atom - field system is h = h_0+v , where h_0 & = & h_+h _ + & = & ^_j(_0 ) + ^_,_(a^_a_+ ) and v & = & -^_j_j + & = & -^_j^_d(+ ) ( _ ) ( _ a_e^i_j+h.c . ) . @xmath150 refers to the internal state of atom @xmath87 , @xmath25 is the atomic transition frequency , and @xmath151 is the frequency of the mode of the electromagnetic field with wave vector @xmath47 and polarization @xmath152 . @xmath104 is the position of atom @xmath87 , and @xmath153 is the polarization direction of the dipoles . we now consider the evolution of the atoms@xmath154field density matrix @xmath155 in the interaction picture . the equation of motion is given by [ eq : atomfield ] & = & , + ( t ) & = & e^ih_0t / ve^-ih_0t/. we can integrate eq . ( [ eq : atomfield ] ) once and substitute the result back into itself . we then trace out the field degrees of freedom to obtain an equation of motion for the atomic density matrix @xmath17 alone : = -_f^t_0d . to make the above equation more useful , we employ the born - markov approximation , replacing @xmath156 with @xmath157 . physically , this amounts to assuming that correlations in the field are negligible , that the field can always be approximated by a vacuum state , and that the correlation time of the atom - field system is much shorter than any atomic timescales . these assumptions safely allow us to extend the time integral to infinity , so that [ eq : bornmarkov ] = -_f^_0d . writing out @xmath158 and @xmath159 , where @xmath160 , give ( t ) & = & ^_a , , -d _ ( _ ) ( e^i_0t+e^-i_0 t ) ( a_e^i(_a-_t ) + a^_e^-i(_a-_t ) ) , + ( t^ ) & = & ^_b,,^ -d_(_^ ) ( e^i_0t^+e^-i_0t^ ) + & & ( a_e^i(^_b-_t^ ) + a^_e^-i(^_b-_t^ ) ) . when we substitute the expansions above into eq . ( [ eq : bornmarkov ] ) and perform the trace , the only nonzero terms will be those associated with @xmath161 , @xmath162 , @xmath163 , and @xmath164 , where @xmath165 and @xmath166 . making these simplifications and replacing the sum on @xmath47 by an integral gives & = & -^_0d^_0k^2dkd^_a , b , d^2 _ ^2(_)^2 . the angular integral is tedious but straightforward and gives d^_e^i_ab(_)^2= 4(^2+(1 - 3 ^ 2)(- ) ) . the time integral can be evaluated using the formula ^_0cde^-ic(kk_0)=(k_0k)i , where @xmath167 denotes the principal value . the delta function above eventually yields the non - hermitian component of the evolution , while the principal value yields the coherent atom - atom interactions . performing all integrals , making the replacement @xmath168 , and changing back to the schrodinger picture yields eq . ( [ eq : dpdt ] ) . the results above were derived treating the atoms in the lattice in first quantization . starting instead from second quantization , one can see that the results derived in first quantization are appropriate in the limit of tight confinement of the atoms , when the overlap between atoms at different sites can be ignored .    in second quantization , the atomic wave function @xmath169 can be expanded in terms of wannier functions : ( ) = ^_nib_ni_n(-_i ) . here @xmath170 denotes the internal state of the atom , @xmath171 is the band index , and @xmath172 labels the lattice sites @xmath101 . to lowest order , we assume that only the harmonic oscillator ground - state wave function is relevant . the atomic hamiltonian is then given in second quantization by _ a=^_ijb_i^b_j_(-_i ) ( + u()+e_)_(-_j ) . in the limit of tight confinement , the overlap integrals for @xmath173 can be ignored , leading back to the atomic hamiltonian used in first quantization . again , in the limit of tight confinement , the overlap integrals for @xmath173 can be ignored , and furthermore the term @xmath174 can be replaced with @xmath175 . in this limit this hamiltonian is equivalent to the electric dipole hamiltonian used in first quantization . hall , m. zhu , and p. buch , j. opt . b * 6 * , 2194 ( 1989 ) . r.j . 85 * , 2462 ( 2000 ) . h. katori , in _ proceedings of the 6th symposium on frequency standards and metrology _ , p. gill , ed . ( world scientific , singapore , 2002 ) , pp . 323 - 330 . t. mukaiyama _ et al . _ , lett . * 90 * , 113002 ( 2003 ) . et al . _ , lett . * 90 * , 193002 ( 2003 ) . et al . _ , j. opt b * 20 * , 968 ( 2003 ) . i. courtillot _ _ , physics/0303023 . t. udem _ et al_. , phys . * 86 , * 4996 ( 2001 ) . g. wilpers _ et al_. , phys . * 89 , * 230801 ( 2002 ) . see , for example , n. allard and j. kielkopf , rev . phys . * 54 * , 1103 ( 1982 ) . m. fleischhauer and s.f . yelin , phys . a * 59 * , 2427 ( 1999 ) . jackson , _ classical electrodynamics _ ( john wiley , new york , 1999 ) , p. 411 . i.h . a * 52 * , 1394 ( 1995 ) . m. weidem@xmath176ller _ et al . a * 58 * , 4647 ( 1998 ) . agarwal , _ springer tracts in modern physics _ , vol . 70 ( springer , berlin , 1974 ) . m. gross and s. haroche , phys . rep . * 93 * , 301 ( 1982 ) .
Abstract:

Answer: motivated by the ideas of using cold alkaline earth atoms trapped in an optical lattice for realization of optical atomic clocks , we investigate theoretically the perturbative effects of atom - atom interactions on a clock transition frequency . these interactions are mediated by the dipole fields associated with the optically excited atoms . we predict resonance - like features in the frequency shifts when constructive interference among atomic dipoles occur . we theoretically demonstrate that by fine - tuning the coherent dipole - dipole couplings in appropriately designed lattice geometries , the undesirable frequency shifts can be greatly suppressed .


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: blazars are most extreme and powerful sources among active galactic nuclei ( agns ) , whose broadband emission is mainly dominated by non - thermal components produced in a relativistic jet pointing toward us @xcite . the blazars are traditionally sub - divided as flat spectrum radio quasars ( fsrqs ) and bl lacertae objects ( bl lacs ) according to their emission - line features , where the bl lacs have weak or no emission lines ( e.g. , equivalent width , ew , of the emission line in rest frame is less than @xmath6 ) while fsrqs have stronger emission lines ( @xmath7 , e.g. , * ? ? ? the different emission - line properties of fsrqs and bl lacs may be triggered by accretion mode transition , where standard cold accretion disk ( e.g. , * ? ? ? * ) exists in fsrqs while advection dominated accretion flow ( adaf , e.g. , * ? ? ? * ; * ? ? ? * ) exists in bl lacs ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? the spectral energy distribution ( sed ) of blazar generally exhibits a two - bump structure . the lower energy bump is commonly ascribed to synchrotron radiation and peaks at infrared to x - ray band , while the second one , attributed to inverse compton ( ic ) scattering , peaks at mev - gev band . the location of the peak for the lower energy bump in the sed , @xmath8 , is also used to classify the sources as low ( lsp , e.g. , @xmath9hz ) , intermediate ( isp , e.g. , @xmath10hz ) and high - synchrotron - peaked ( hsp , e.g. , @xmath11hz ) blazars ( e.g. , * ? ? ? * ; * ? ? ? in contrast to bl lac objects , fsrqs are essentially all lsp blazars . origin of the seed photons for the ic scattering in blazars is an open issue , which include the synchrotron photons ( ssc process , * ? ? ? * ; * ? ? ? * ) and/or external photons ( ec process ) , where the external photons possibly originate from accretion disk ( e.g. , * ? ? ? * ) , broad line region ( blr , e.g. , * ? ? ? * ; * ? ? ? * ) , and/or molecular torus ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) . normally , the seds of hsp blazars ( e.g. , bl lacs ) seem to be consistent with the pure ssc models ( e.g. , * ? ? ? * ; * ? ? ? * ) , while the lsp blazars often require an ec component to explain their @xmath1-ray spectra ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? no ec process in bl lacs may be caused by the disappearance of the blr ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) and/or torus(e.g . * ; * ? ? ? * ) in the low - power sources . identifying the origin of external seed photons in fsrqs is also a diagnostic for the location of the @xmath1-ray emitting region ( e.g. , * ? ? ? for example , the external seed photons will dominantly originate from the cold disk if the @xmath1-ray emitting region is less than several hundreds schwarzschild radius ( @xmath12 ) from the central black hole ( bh ) , while the soft photons will mainly come from the blr if the @xmath1-ray emitting region is larger than several hundreds @xmath12 and less than several thousands @xmath12 . if the @xmath1-ray emitting region stay outside the blr and less than @xmath13 , the most abundant seed photons are ir photons coming from the molecular torus . the cosmic background photons may be important if the gamma - ray emitting region is much larger than @xmath13 ( see figures 2 - 3 in * ? ? ? the jet matter content is also a mystery , which possibly include a ` normal ' plasma of protons and relativistic electrons ( an @xmath5 jet ) , a ` pair ' plasma consisting only of relativistic electrons and positrons ( an @xmath4 jet ) , or combination of them . discriminating between these possibilities is crucial in understanding the physical processes of jet occurring close to the central engine . there are several methods were proposed to explore this issue . the first one is based on synchrotron self - absorption ( ssa ) arguments combined with total kinetic power of jets ( e.g. , * ? ? ? * ; * ? ? ? * ) , where the measurements on the flux and size of radio core lead to constraint on number density of relativistic electrons and magnetic field strength in the jet . using this method , @xcite concluded that the core of m 87 is probably dominated by an @xmath4 plasma . the second one is the constraint from circular polarization in the radio core , which is due to the faraday conversion of linear to circular polarization caused by the low - energy electrons . @xcite suggested that the jet of 3c 279 is composed mainly by @xmath4 pairs based on its small minimum lorentz factor of electrons ( @xmath14 ) that derived from the circular polarizations of radio core . the third constraint comes from the cocoon dynamics of the jets . @xcite suggested that the jet of cygnus a should be dominated by mixture compositions of @xmath15-@xmath16-@xmath17 by comparing the partial pressures of electrons and protons with the observational pressure of cocoon that estimated from the interaction between jet and interstellar medium . the fourth approach is the constraint from the absence of bulk - compton emission in fsrqs ( e.g. , * ? ? ? * ; * ? ? ? * ) , where the electrons in the jet traveling with a bulk lorentz factor @xmath18 will interact with the photons produced by accretion disk and/or broad emission lines and will radiate at soft x - rays . @xcite and @xcite argue against the pure @xmath4 pair jets based on the absence of observational feature of the soft x - ray excess . the jet power is important for understanding jet formation , disk - jet relation , jet kinetic luminosity function , and jet feedback etc . ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? the x - ray cavities in galaxy clusters and giant elliptical galaxies provide a direct measurement of the mechanical energy released by the agn jets through the work done on the hot , gaseous halos surrounding them ( i.e. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? the jet power of fr ii radio galaxies and radio quasars can be estimated directly from the measurements of the hotspot size and equipartition magnetic field strength along with some assumptions ( i.e. , * ? ? ? @xcite estimated the jet power from the number of emitting particles , magnetic field strength and jet velocity that constrained from the radio and x - ray data using the standard ssc theory ( see also * ? ? ? * ; * ? ? ? * ; * ? ? ? modeling the multi - waveband seds of blazars also returns the particle number , magnetic field strength and doppler factor of the jet , where the kinetic jet power can also be derived with some reasonable assumption ( i.e. , * ? ? ? * ; * ? ? ? * ; * ? ? ? it should be noted that the last two methods suffer the uncertainties of the jet matter content and the low - energy cutoff , @xmath2 , of electrons , where the parameter @xmath2 is coupled with the parameter of electron number density ( i.e. , * ? ? ? . the value of @xmath2 can be strongly constrained in some lsp blazars ( e.g. , fsrqs ) if their soft x - ray emission is dominated by the ssc process ( i.e. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) .    in this work , we explore the issue of possible external seed photon fields ( e.g. , blr or ir torus ) for a sample of lsp blazars by a @xmath0 test in the sed modeling . the parameter of @xmath2 can be constrained for these lsp blazars with good quality of soft x - ray data in the sed modelings . we further investigate the issue of jet power and jet matter content if the @xmath2 is known . the sample is described in section 2 . section 3 contains the detailed information of our model . we present the results in section 4 . the last section was devoted to discussion and conclusion . throughout this work , we assume the following cosmology : @xmath19 , @xmath20 and @xmath21 . 0.60 mm    [ cols= " < , < , < , < , < , < , < , < , < , < , < , < , < , < , < , < , < " , ]     in figure 5 , we present the relation between the power in bulk motion of electrons , protons and magnetic field(@xmath22 ) and the 151 mhz radio luminosity @xmath23 for the lsp blazars , where the 151 mhz radio fluxes are selected from _ _ ned _ _ ( see table 2 ) . for comparison , we also present the @xmath24 relation for a sample of fr i ( empty stars ) and fr ii ( empty squares ) radio galaxies ( see * ? ? ? * for more details ) , where the jet kinetic power of fr is and fr iis are estimated from the x - ray cavities and cocoon dynamics respectively . we can clearly find that the jet kinetic power of blazars is systematically higher than that of fr i / iis at given 151 mhz radio luminosity , even though fr iis are normally regarded as same as lsp blazars in unification scheme except their jet viewing angles . we speculate that it may be caused by the assumption of one proton per emitting electron in the jet , and the jet power will be reduced if the jet also include a fraction of @xmath4 pairs .    to evaluate the possible positrons in the jet , we further assume @xmath25 and @xmath26 in the model , where @xmath27 corresponds to the normal jet with pure @xmath5 plasma while @xmath28 corresponds to the jet with pure @xmath4 pairs . in this case , the number density of protons and the jet power carried by protons will be reduced by a factor of @xmath29 , since that positrons emit with the same kind of energy spectrum as electrons . it is possible to derive the value of @xmath30 if the jet power is known from an independent method . to do this , we calculate the jet kinetic power of these lsp blazars from their 151 mhz radio luminosities using the relation of @xmath24 for fr iis @xcite , @xmath31 where the jet kinetic power of fr iis are estimated from their cocoon dynamics . @xmath32 derived from equation ( 7 ) for each lsp blazar is listed in table 2 . by setting @xmath33 , we derive the @xmath30 values for 21 sources , which are listed in column ( 10 ) of table 2 . we noted that other nine sources have @xmath34 , which may be caused by the uncertainties in jet power estimation from @xmath24 or by the uncertainties in the sed modeling . the jet power carried by protons is most possibly negligible in these nine blazars and we will simply assume @xmath35 in them . the distribution of @xmath30 is shown in figure 6 , where the median @xmath30 value is 0.08 ( average @xmath30 value is 0.22 ) . the jet power for cold protons , @xmath36 , electrons / positrons , @xmath37 , poynting flux , @xmath38 , and radiation , @xmath39 , were presented in figure 7 when considering the possible @xmath4 pairs , where the total jet power estimated from @xmath23 is also shown in the bottom panel ( dashed lines represent median values ) .    in figure 8 , we compare the doppler factor , @xmath40 , that derived from our sed fittings with the variability doppler factor , @xmath41 , that estimated from the variability brightness temperature by assuming the intrinsic brightness temperature is limited to an equipartition value ( @xcite , see also @xcite for more details ) . they are positively correlated ( pearson correlation coefficient @xmath42 and chance probability @xmath43 ) and roughly consistent with each other . the simultaneous or quasi - simultaneous _ plank _ , _ swift _ , _ fermi _ and ground - based observations for a sample of blazars provide a great opportunity to explore their jet properties through modeling their seds . our model roughly become transparent at submm waveband for the typical jet parameter of lsp blazars . the fairly good correlation of the variability between mm and optical signals ( e.g. , * ? ? ? * ) or between @xmath44ray and mm signals ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) also support this scenario . therefore , the emission from mm to @xmath44ray may originate from a more or less similar region , and it should be reasonable to fit these data with the one - zone homogeneous model . from the fitting results , we find that the first hump from mm to optical band can be well reproduced by the synchrotron emission , while the second hump can be reproduced by ssc and ec emission for these lsp blazars ( see figures a1 - 6 ) . it should be cautious that the 1@xmath45 uncertainty of model parameter in table ( 1 ) is derived by setting other parameters to be their best - fit values , which may be underestimated if we allow all the parameters to be free ( particularly in the case of the possible degeneracy in the model ) . the observational data are not very good in a few sources ( e.g. , only 8 points in j1911 - 2006 , or data is absent in a certain waveband ) , which will not affect our main conclusion in a statistical sense . the consistency of the doppler factor from the sed fitting with that derived from independent method suggest that our sed modeling should be reasonable , even though there is no reason that the radio doppler factor should be the same as that derived from the @xmath1-ray emitting region ( e.g. , * ? ? ? * ) .      for lsp blazars , the @xmath1-ray emission is mainly contributed by the ec process , and two main candidates of external seed fields have been proposed . the first one is the blr if the @xmath1-ray emitting region stay inside the blr , which can easily explain short timescales of variability reported in some @xmath1-ray observations of blazars ( e.g. , * ? ? ? * ; * ? ? ? * ) , even though the short timescales of variability do not necessarily imply the short distance to the bh . the second candidate is the molecular torus if the @xmath1-ray emitting region stay outside the blr and up to @xmath46 ( e.g. , * ? ? ? * ; * ? ? ? * ) . in this work , we explore this issue through the sed modelings , where the external seed photons from blr or molecular torus are considered respectively . we find that the soft photons may dominantly come from ir torus in these lsp blazars based on the @xmath0 tests ( see figure 1 ) , which may be caused by the kn effect . assuming a typical seed photon frequency @xmath47 , the compton scattering by the electron with a lorentz factor @xmath1 yields @xmath48 within the thomson regime in observational frame ( assuming @xmath49 ) . combining the thomson scattering condition @xmath50 , we obtain @xmath51 . if seed photons come from blr , we get @xmath52 hz for @xmath53 hz ( see also * ? ? ? this implies that the ic emission would be significantly suppressed by kn effect if @xmath54hz , which will lead to a very steep spectrum . however , we have @xmath55 hz @xmath56 gev if main seed photons come from dust torus with @xmath57 hz . therefore , the @xmath1-ray spectrum of ec model with ir seed photons will not become very steep spectrum at observational band . this can roughly explain why the fittings are systematically better for ec model with ir seed photons ( see also * ? ? ? * ; * ? ? ? @xcite proposed that the ir external field may play an important role through analyzing the ratio of ec to synchrotron luminosity for a sample _ fermi _ bright blazars . @xcite found that the bright blazars favor the seed photons from molecular torus over those from the blr based on the lack of the bulk - compton and the klein - nishina features in the broadband spectra . our modeling results give an indirect evidence that the @xmath1-ray emitting region may locate outside the blr and within the molecular torus , which roughly correspond to several thousand to @xmath58 from the bh . we note that the location of @xmath1-ray emitting region should be determined by the energy dissipation or particle acceleration processes within the relativistic jet that responsible for the generation of nonthermal electrons . @xcite explored the velocity field of m87 jet and found that the jet matter is mainly accelerated at region of @xmath46 , and the jet is decelerated at larger radius . @xcite s result suggest that most of jet energy was dissipated at @xmath59 , which may correspond to the @xmath1-ray emitting region of lsp blazars ( e.g. , @xmath60 ) . the detection of @xmath1-rays with rest - frame energy above 20 gev ( e.g. , * ? ? ? * ; * ? ? ? * ) suggest that the @xmath1-ray emitting region should be located outside the highly opaque ( @xmath615 - 10 ) blr , which would not permit photons of such high energies to escape ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? based on the simultaneous flares at mm and @xmath1-ray waveband , it was also suggested that the @xmath1-ray emitting region should be far away from the blr region where the high - energy photons should be produced in a region which is already transparent to the radiation at mm wavelength ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) . through the correlation of mm with @xmath1-ray light curves and direct ultrahigh - resolution 7 mm imaging with the very long baseline array ( vlba ) , @xcite argued that the location of @xmath1-ray emitting region should be @xmath6214 pc from the bh in the jet of oj 287 , which roughly correspond to @xmath63 if considering the bh mass @xmath64 ( e.g. , * ? ? ? * ) . these independent evidences support our conclusions that the @xmath1-ray emitting region may stay outside the blr ( several thousands @xmath12 ) and the dominant seed photons for ec should be dominated by the dusty torus rather than blr . it is well known that the minimum electron lorentz factor ( or the low - energy cutoff ) , @xmath2 , play a crucial role in estimating the powers in particles ( e.g. , @xmath65 ) . normally , this parameter is poorly constrained due to its synchrotron radiation will be self - absorbed . however , these low - energy electrons would instead contribute to the low - energy part of ssc emission . the x - ray emission of lsp blazars is known to be dominated by the ssc ( i.e. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) , which provide a possibility to constrain @xmath2 through modeling their multiwavelength seds ( see figure 3 ) . in this work , we select a sample of lsp blazars with good soft - x - ray data and find that @xmath2 ranges from 5 to 160 ( with a median of 55 ) through their sed modelings . it is interesting to note that @xmath2 value is not sensitive to the possible external seed photon field because it was mainly constrained from the ssc process ( see figure 4 ) . the @xmath2 values around several tens for most of lsp blazars in our fittings are quite consistent with that reported in some of recent works even most of their fittings are evaluated by  eyeball " ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? from a theoretical perspective , @xcite derived the minimum electron lorentz factor @xmath66(m_p / m_e)\varepsilon_e\delta\gamma$ ] ( @xmath67 ) based on the electron energy distribution and jump conditions for a relativistic shock , where @xmath68 is the fraction of shock energy goes into the electrons , @xmath17 is the index of power - law electron distribution ( e.g. , @xmath69 ) and @xmath70 is lorentz factor difference for two colliding shells . we can derive @xmath71 for typical value of @xmath72 for relativistic shock acceleration ( e.g. , * ? ? ? * ) , @xmath73 in gamma - ray bursts and blazars ( e.g. , * ? ? ? * ; * ? ? ? * ) and @xmath74 in blazars , which can roughly explain @xmath75 several tens as constrained from the observations . it should be noted that the lower limits of @xmath2 are not well constrained for several sources ( e.g. , j0457 - 2324 , j1517 - 2422 , j1642 + 3948 , j2202 + 4216 and j2203 + 3145 ) , which is caused by the poor observational data or its x - ray emission is contributed by both ssc and ec .    with the constrained of @xmath2 , we calculate the jet power for each source . we find that the jet kinetic power of lsp blazars is systematically higher than that of fr iis at given 151 mhz luminosity ( see figure 5 ) if assuming one proton per emitting electron in the jet , where these two types of agns are assumed to be intrinsically same except the jet viewing angle ( e.g. , * ? ? ? * ; * ? ? ? one possible reason is that the jet include some positrons , which will reduce the jet power . the x - ray cavities / bubbles in galaxy clusters ( or giant galaxies , i.e. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) and cocoons associated with fr i / ii galaxies ( i.e. , * ? ? ? * ; * ? ? ? * ) provide the possibility to estimate the jet power independently . assuming the jet power of lsp blazars is intrinsically the same as that of fr iis at given 151 mhz luminosity , we find that the jet of most of lsp blazars should be dominated by a pair plasma , where the median ratio of proton to electron is around 0.08 or the number density of positron is around 10 times higher than that of proton . it should be noted that @xmath76 may also caused by some fraction of kinetic jet power will be converted into radiation before it is dissipated in hotspots on large scales . we find that the total radiative jet power normally occupy @xmath77 several percent of jet kinetic power in these lsp blazars , where most of jet radiation come from the @xmath44ray emitting region . after subtracting the radiative power derived from our modeling ( as an approximation for total radiative power ) , we find that @xmath30 value increases a little bit with a median value of 0.22 ( average value is 0.29 ) . it means that the number density of positrons is still around 4 times higher than that of protons . the assumption of full nonthermal electrons in jet may be too simple since that only a small fraction of thermal electrons are accelerated into power - law distribution and most of electrons remains in thermal pool . if this is the case , our conclusion will be strengthened since that the radiation from these thermal electrons is much less than that of nonthermal electrons but the jet kinetic power will increase and the @xmath30 value will decrease ( i.e. , more positrons are needed ) . it should be noted that there are still large uncertainties in estimation of jet power through @xmath78 relation , and , therefore , our results should be only statistically meaningful . based on the possible effect of anisotropic external seed photon field in jet comoving frame , @xcite set the upper limit of @xmath79 on the pair to proton number ratio ( see also * ? ? ? * ) . @xcite claimed that the jets contain more @xmath4 pairs than protons based on the absence of bulk - compton emission in fsrqs , but that jets are still dynamically dominated by protons . based on the vlbi observations and theory of ssa , it was also found that the jet should be dominated by @xmath4 plasma ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? our conclusion roughly consistent with these results even they are derived from different methods .      in this work , we employ the one - zone homogeneous leptonic jet model and @xmath0-minimization procedure to fit the simultaneous or quasi - simultaneous multi - waveband seds for a sample of lsp blazars , where the external seed photons originated from ir torus and blr are considered respectively . our main results are summarized below . \1 ) the sed fitting with external seed photon from ir torus is systematically better than that from blr . this result suggests that the @xmath1-ray emitting region of these lsp blazars most possibly stay outside the blr . \2 ) with the good quality soft x - ray data combined with other multi - wavelength observations , we find that minimum electron lorentz factors , @xmath2 , range from 5 to 160 for these lsp blazars with a median value of 55 , which are not affected by the possible uncertainties of external seed photons . \3 ) assuming one - to - one ratio of proton and electron in jet , we find that the jet power estimated from the fitting parameters is much higher than that of fr ii galaxies at given 151 mhz radio luminosity even though they are assumed to be intrinsically the same in the unification scheme . therefore , we propose a mixture composition of @xmath80 in the jet of these lsp blazars . the number density of @xmath4 pairs should be several times higher than that of @xmath5 pairs if assuming that the jet power of lsp blazars is the same as that of fr iis at given 151 mhz radio luminosity . we appreciate the anonymous referee for his / her helpful suggestions and comments , which help to clarify and improve our work . this work is supported by the nsfc ( grants 11103003 , 11133005 , 11103060 and 11233006 ) , national basic research program of china ( 2009cb824800 ) , and new century excellent talents in university ( ncet-13 - 0238 ) . abdo , a. , et al . 2009 , apj , 697 , 934 abdo , a. et al . 2010 , , 716 , 30 ackermann , m. , et al . 2010 , , 721 , 1383 agudo , i. , et al . 2011 , , 735 , 10 asada , k. , nakamura , m. , doi , a. , nagai , h. , inoue , m. 2014 , , 781 , 2 allen , s. , dunn , r. , fabian , a. c. , taylor , g. b. , & reynolds , c. s. 2006 , , 372 , 21 aleksi c. j. , et al . , 2011 , , 730 , l8 arbeiter , c. , pohl , m. , schlickeiser , r. 2002 , , 386 , 415 bai , j.  m. , liu , h.  t. , & ma , l.  2009 , , 699 , 2002 bednarz , j. , & ostrowski , m. 1998 , phys . lett . , 80 , 3911 bentz , m.  c. , peterson , b.  m. , netzer , h. et al . 2009 , apj , 697 , 160b , r.  d. , znajek , r.  l. 1987 , , 322 , 650 blazejowski , m. , et al . 2000 , , 545 , 107 blumenthal , g. r. , & gould , r. j. 1970 , rev . phys . , 42 , 237 bttcher , m. , basu , s. , joshi , m. , et al . 2007 , , 670 , 968b bttcher , m. , & chiang , j. 2002 , , 581 , 127 brzan , l. ; rafferty , d. , mcnamara , b. , wise , m. , & nulsen , p. 2004 , , 607 , 800 brzan , l. , mcnamara , b. , nulsen , p. , carilli , c. , & wise , m. 2008 , , 686 , 859 brown , a.  m.  2013 , , 431 , 824 cao , g. & wang , j. c. 2013 , mnras , 436 , 2170c cao , x. 2003 , apj , 599 , 147 cao , x. 2003 , apj , 724 , 855 cao , x. 2014 , apj , 783,51 cao , x. & spruit , h. c. 2013 , apj , 765 , 149 cavagnolo , k. w. , et al . 2010 , , 720 , 1066 celotti , a. , & fabian , a.c . 1993 , , 264 , 228 celotti a. , & ghisellini g. 2008 , , 385 , 283 chen , l. , & bai , j.  m.  2011 , , 735 , 108 chen , l. , bai , j.  m. , zhang , j. , liu , h.  t. 2010 , research in astronomy and astrophysics , 10 , 707 cleary , k. , lawrence , c. r. , marshall , j. a. , hao , l. , meier , d. 2007 , , 660 , 117 dammando , f. , antolini , e. , tosti , g. , et al . 2013 , , 431 , 2481d dermer , c.d . , & schlickeiser , r. 1993 , , 416 , 458 donea a .- c . , protheroe r. j. , 2003 , aph , 18 , 377 dutka , m. s. , ojha , r. , pottschmidt , k. , et al . 2013 , , 779 , 174d dunn , r. j. h. , fabian , a. c. , & celloti , a. 2006 , , 372 , 1741 elitzur , m. & ho , l. c. 2009 , , 701 , 91 fabian , a. c. , et al . 2006 , , 318 , 65 fan , j .- h . , yang , j .- h . , liu , y. , zhang , j .- y . 2013 , research in astronomy and astrophysics , 13 , 259 fan , j .- h . , huang , y. , he , t .- m . , yang , j. h. , hua , t. x. , liu , y. , wang , y. x. 2009 , pasj , 61 , 639 fan , z. , cao , x. , & gu , m. 2006 , , 646 , 8 foschini , l. , ghisellini , g. , tavecchio , f. , bonnoli , g. , stamerra , a. 2011 , , 530 , 77 fossati , g. , buckley , j.  h. , bond , i.  h. , et al . 2008 , , 677 , 906 ghisellini , g. , celotti , a. , fossati , g. , maraschi , l. , & comastri , a. 1998 , mnras , 301 , 451 ghisellini , g. , & tavecchio , f. 2008 , , 387 , 1669 g ghisellini g. 2012 , , 424 , 26 ghisellini , g. , & celotti , a. 2001 , , 379 , 1 ghisellini , g. , tavecchio , f. 2009 , , 397 , 985 ghisellini g. & tavecchio f. 2010 , , 409 , 79 ghisellini , g. , tavecchio , f. , & ghirlanda , g. 2009 , , 399 , 2041 ghisellini , g. , tavecchio , f. , foschini , l. , et al . 2010 , mnras , 402 , 497 g giommi , p. , et al . 2012 , , 541 , 160 godfrey , l. e. h. , & shabala , s. s. 2013 , , 767 , 12 gu , m. , cao , x. , & jiang , d. 2009 , , 396 , 984 . gu , q. , & huang , j. 2002 , , 579 , 205 hirotani , k. 2005 , , 619 , 73 ho , l. c. 2008 , ara&a , 46 , 475 h'onig , s. f. & beckert , t. 2007 , , 380 , 1172 hovatta , t. , valtaoja , e. , tornikoski , m. , lahteenmaki , a. 2009 , , 498 , 723 jorstad , s. g. , marscher , a. p. , larionov , v. m. et al . 2010 , apj , 715 , 362j kaspi , s. , brandt , w.  n. , maoz , d. et al . 2007 , apj , 659 , 997k kataoka , j. et al . 2008 , , 672 , 787 kino , m. , kawakatu , n. , & takahara , f. 2012 , , 751 , 101 krawczynski , h. , et al . 2004 , , 601 , 151 knigl , a.,1981 , , 243 , 700 . kushwaha , p. , sahayanathan , s. , singh , k. p. 2013 , mnras , 433 , 2380k laor , a. 2003 , , 590 , 86 lei , w .- h . , wang , d .- x . , & ma , r .- y . 2005 , , 619 , 420 lei , w .- h . , wang , d .- x . , zou , y .- c . , zhang , l. 2008 , chja&a , 8 , 404 len - tavares , j. , valtaoja , e. , giommi , p. , et al . 2012 , apj , 754 , 23l li , s .- l . , & cao , x. 2012 , , 753 , 24 liang , e. w. , liu , h. t. 2003 , 340 , 632l liu , y. & jiang , d .- 2006 , , 637 , 669 liu , h. t. ; bai , j. m. 2006 , , 653 , 1089 lubow , s. h. , papaloizou , j. c. b. , pringle , j. e. 1994 , mnras , 267 , 235 malmrose , m. p. , marscher , a. p. , jorstad , s. g. , et al . 2011 , apj , 732 , 116 m marscher , a. p. , jorstad , s. g. , larionov , v. m. , et al . 2010 , apj , 710l , 126 m maraschi , l. , tavecchio , f. , 2003 , , 593 , 667 maraschi , l. , celotti , a. , & ghisellini , g. 1992 , , 397 , 5 marscher , a. p. , & gear , w.k . 1985 , , 298 , 114 mastichiadis , a. , & kirk , j. g. 1997 , , 320 , 19 mehta , k. t. , georganopoulos , m. , perlman , e. s. , et al . 2009 , apj , 690 , 1706 m narayan , r. , & mcclintock , j. e. 2012 , , 419 , 69 narayan r. , yi i. , 1994 , apj , 428 , l13 narayan r. , yi i. , 1995 , apj , 444 , 231 orienti , m. ; koyama , s. ; dammando , f. et al . 2013 , mnras , 428 , 2418o pacciani l. , et al . , 2012 , mnras , 425 , 2015 padovani , p. & giommi , p. 1995 , , 444 , 567 paliya , v. s. , stalin , c. s. , shukla , a. , sahayanathan , s. 2013 , , 768 , 52p panaitescu , a. , & kumar , p. 2001 , , 560 , l49 potter , w. j. , & cotter , g. 2013 , mnras , 431 , 1840p poutanen j. , stern b. , 2010 , , 717 , l118 reynolds , c. s. , fabian , a. c. , celotti , a. , rees , m. j. 1996 , mnras , 283 , 873 rafferty , d. , mcnamara , b. , nulsen , p. , & wise , m. w. 2006 , , 652 , 216 rybicki , g. b. , & lightman , a. p. 1979 , radiative processes in astrophysics ( new york : wiley - interscience ) , 393 sahayanathan , s. & godambe , s. 2012 , mnras , 419 , 1660s sambruna , r. , et al . 1999 , , 515 , 140 sari , r. , piran , t. , & narayan , r. 1998 , , 497 , 17 savolainen , t. ; homan , d. c. ; hovatta , t. ; et al . 2010 , a&a , 512a , 24s shakura n. i. , sunyaev r. a. , 1973 , a&a , 24 , 337 sikora , m. , moderski , r. , madejski , g. 2008 , , 675 , 71 sikora , m. & madejski , g. 2000 , , 534 , 109 sikora , m. , begelman , m. c. , & rees , m. j. 1994 , , 421 , 153 sikora , m. & begelman , m. c. 2013 , , 764 , 24 sikora , m. , stawarz ,  . , moderski , r. et al . 2009 , apj , 704 , 38s sokolov , a. , & marscher , a. p. 2005 , , 629 , 52 spruit , h. c. 2010 , lecture notes in physics , 794 , 233 tavecchio , f. , pacciani , l. , donnarumma , i. , stamerra , a. , isler , j. , macpherson , e. , & urry , c. m. 2013 , , 435 , 24 tavecchio , f. , maraschi , l. , ghisellini , g. , kataoka , j. , foschini , l. , sambruna , r. m. , & tagliaferri , g. 2007 , , 665 , 980 tavecchio , f. , maraschi , l. , sambruna , r  m. , urry , c.  m. 2000 , , 544 , 23 t tchekhovskoy , a. , narayan , r. , & mckinney , j. c. 2011 , , 418 , 79 tran , h. d. 2001 , , 554 , l19 urry , c. m , padovani , p. 1995 , pasp , 107 , 803 valtonen , m. j. , ciprini , s. , lehto , h. j. 2012 , , 427 , 77v vercellone , s. , striani , e. , vittorini , v. et al . 2011 , , 736l , 38v vovk , ie . , neronov , a. 2013 , , 767 , 103v xu , y .- d . , cao , x. , & wu , q. 2009 , , 694 , 107 wang , j .- m . , staubert , r. , ho , l. c. 2002 , , 579 , 554 wardle , j. f. c. , homan , d. c. , ojha , r. , & roberts , d. h. 1998 , nature , 395 , 457 wehrle , a. e .. et al . 2012 , , 758 , 72 wu , q. , cao , x. , ho , luis c. ; wang , d .- x . 2013 , , 770 , 31 wu , q. , cao , x. , & wang , d. x. 2011 , , 735 , 50 wu , q. , yuan , f. , & cao , x. 2007 , , 669 , 96 wu , q. & cao , x. 2006 , pasp , 118 , 1098 yan , d. , zeng , h. , & zhang , l.  2014 , , 439 , 2933 zhang , j. , sun , x .- n . , liang , e .- w . , lu , r .- j . , lu , y. , zhang , s .- n . 2014 , submitted to apj zhang , j. , liang , e .- w . , zhang , s .- n . , & bai , j.  m.  2012 , , 752 , 157
Abstract:

Answer: we fit the ( quasi-)simultaneous multi - waveband spectral energy distributions ( seds ) for a sample of low - synchrotron - peaked ( lsp ) blazars with a one - zone leptonic model . the seed photons that predominantly come from broad line region ( blr ) and infrared ( ir ) molecular torus are considered respectively in external compton process . we find that the modeling with ir seed photons is systematically better than that with blr photons based on a @xmath0 test , which suggest that @xmath1-ray emitting region most possibly stay outside the blr . the minimum electron lorentz factor , @xmath2 , is constrained from the modeling for these lsp blazars with good soft x - ray data ( ranges from 5 to 160 with a median value of 55 ) , which plays a key role in jet power estimation . assuming one - to - one ratio of proton and electron , we find that the jet power for lsp blazars is systematically higher than that of fr ii radio galaxies at given 151 mhz radio luminosity , @xmath3 , even though fr iis are regarded as same as lsp blazars in unification scheme except the jet viewing angle . the possible reason is that there are some @xmath4 pairs in the jet of these blazars . if this is the case , we find the number density of @xmath4 pairs should be several times higher than that of @xmath5 pairs by assuming the jet power is the same for lsp blazars and fr iis at given @xmath3 .