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: quantization of classical chaotic systems has long been of interest because of the challenges in identifying quantum signatures of systems that , in a classical limit , exhibit chaotic behavior . various signatures have been identified , such as the spectral properties of the generating hamiltonian  @xcite , phase space scarring  @xcite , hypersensitivity to perturbation  @xcite , and fidelity decay  @xcite , which indicate an underlying chaotic presence in the quantum dynamics . here , we investigate the issue of quantum signatures from the perspective of entanglement : as entanglement is at the heart of quantum mechanics and a crucial resource for quantum information processing  @xcite , the entanglement inherent in quantum chaotic systems could provide a valuable approach to studying decoherence and quantum chaos  @xcite . furthermore , quantum chaos could be seen as an engine for generating entanglement ; indeed quantum chaos could lurk in quantum information processing  @xcite as a deleterious or perhaps even as an advantageous effect . we study entanglement by considering a finite multipartite system , whose collective dynamics obey chaotic hamiltonian dynamics in the classical limit .    previous studies of entanglement in chaotic systems  @xcite have explored bipartite entanglement in pure states , entanglement of qubits in a multiqubit system and average entanglement or entangling power . in this study we explore these different types of entanglement in a single simple system , namely the quantum kicked top  ( qkt )  @xcite . this enables us to compare the degree to which these different types of entanglement display signatures of chaos . an advantage of the dynamics of the qkt is that it obeys a spin algebra symmetry . this spin system can thus be regarded as a composite of distinct spin - half particles thereby admitting a multiqubit interpretation . our system thus allows us to study and compare pairwise entanglement between two qubits as well as bipartite entanglement between the two qubits and the rest of the qubits . the hilbert space for the qkt is finite and the poincar section of the phase space compact , allowing analyses of quantum and classical dynamics uncomplicated by truncation issues . the qkt is well studied and understood thereby simplifying the analysis of the role of entanglement in the system . finally , the qkt possesses a parity symmetry that allows coherent quantum tunnelling to occur for states localized at classical fixed points  @xcite . the hamiltonian evolution may increase the entanglement of the multipartite system , initiated in a collective spin coherent state  @xcite . for this analysis it suffices to employ two entanglement measures . for bipartite entanglement , where the multipartite system is divided into two subsystems , entropy of a subsystem is used to quantify the degree of entanglement between the two subsystems . pairwise entanglement , on the other hand , considers the degree of entanglement between two qubits traced over all remaining qubits and is quantified by the concurrence  @xcite . we present general results for both bipartite and pairwise entanglement in the multipartite qkt and demonstrate that these entanglement measures reveal strong signatures of the classical chaos features corresponding to the onset of chaos and to the edge of chaos  @xcite , which is the boundary between regular and chaotic regimes in the classical phase space . we have studied the behavior of the linear entropy and the concurrence for specific initial states as well as the dynamics of these quantities when averaged over all initial states . we show that the linear entropy increases more rapidly for an initial state centered in a chaotic region of the classical phase space than one centered on a fixed point . this agrees with the behavior of the linear entropy observed in other chaotic systems , supporting the conjecture that classical chaos can enhance bipartite entanglement  @xcite . furthermore , we show that the pairwise entanglement as measured by the concurrence also reveals a dramatic change for a spin coherent state whose mean traverses the edge of chaos on its transit through chaotic and regular regions of the phase space . contrary to the linear entropy , the concurrence rapidly decreases for an initial state located in the chaotic region . while the linear entropy and concurrence can reveal the local chaotic and regular structures in phase space , the entangling power , which is the averaged bipartite or pairwise entanglement , can be used to identify the edge of chaos and quantify the onset of global chaos , much like the classical lyapunov exponent . we show that the entangling power greatly increases as the chaoticity parameter  @xmath0 is increased , and the corresponding classical kicked top makes the transition from predominantly regular to predominantly chaotic behavior . in particular , the behavior of the average linear entropy is qualitatively similar to that of the classical lyapunov exponent , thus revealing a signature of a global feature of the classical chaos . the paper is organized as follows . in sec . ii , we introduce the qkt and its classical dynamics , and introduce bipartite and pairwise entanglement measures . in sec . iii , we study in detail dynamical evolutions of bipartite and pairwise entanglement , and examine the edge of quantum chaos , the onset of quantum chaos via the entangling power . we conclude in sec . the qkt is described by the hamiltonian  @xcite @xmath1 where @xmath2 are spin operators and states are restricted to irrep @xmath3 for which @xmath4 . @xmath5 is the duration between periodic kicks , @xmath6 is the strength of each kick ( which is manifested as a turn by angle @xmath6 ) , and @xmath0 is the strength of the twist . the hamiltonian is an alternative sequence of twists ( @xmath7 term ) and turns ( @xmath8 term ) . the qkt describes a spin system , which can be comprised of multiple systems of lower spins . for @xmath9 the pauli operators for the @xmath10 qubit , a collective spin operator , satisfying the usual su(2 ) algebra , is given by @xmath11 an example of using multiple qubits to simulate the qkt has been presented for trapped ions  @xcite . a standard dynamical description of the qkt is via the floquet operator @xmath12 where the energy is rescaled so that @xmath13 , and @xmath14 are henceforth assumed . the orthogonal eigenstates of @xmath15 , denoted by @xmath16 , which satisfy @xmath17 provide a convenient basis for stroboscopic evolution . an arbitrary state @xmath18 evolves to @xmath19 with @xmath20 . the qkt  @xcite is chaotic in the classical limit . for integrable systems , it is well known from semi - classical theory that the classical actions can be associated with corresponding regular eigenstates of the quantum system with a well defined quantum number . this correspondence breaks down in chaotic systems  @xcite . in quasi - integrable systems with a mixed classical phase space of regular and chaotic regions , some of the eigenstates can still be associated with local actions in the regular regions with corresponding discrete eigenenergies . the remaining eigenstates result in an irregular energy spectrum corresponding to the chaotic region  @xcite . we show here that this underlying regular and chaotic energy spectrum of the floquet eigenstates of the qkt is reflected in the dynamics of the entanglement , depending on whether the initial state is in the regular or chaotic region of the classical phase space . the classical limit of the qkt is obtained by expressing @xmath21 and similar for @xmath22 and @xmath23 and factorizing all moments such as @xmath24 to products of first - order moments . then the classical equations of motion , obtained from the heisenberg operator equations of motion and applying the factorization rule above , are given by  @xcite @xmath25 the stroboscopic evolution described by eq.([map ] ) can be represented in a phase space given by a sphere @xmath26 of unit radius . the classical , normalized angular momentum variables @xmath27 can be parametrized in polar coordinates as @xmath28 where @xmath29 and @xmath30 are the polar and azimuthal angles , respectively . thus the mapping domain is essentially two - dimensional . the stroboscopic dynamics of the classical map is shown in fig .  1 . in the plot , we choose the chaoticity parameter @xmath31 which yields a mixture of regular and chaotic areas of significant size . elliptic fixed points surrounded by the chaotic sea are evident . two such elliptic fixed points have coordinates @xmath32 and @xmath33 as we will see , this phase space structure of the classical kicked top determines behavior of quantum entanglement in the qkt . three hundred stroboscopic trajectories are plotted , each for a duration of 300 kicks.,scaledwidth=40.0% ]      pure - state bipartite entanglement has been calculated in previous studies that explore connections between quantum entanglement and underlying classical chaos  @xcite . quite recently , bettelli and shepelyansky  @xcite studied the behavior of the concurrence in a system exhibiting quantum chaos , and found that the underlying classical chaos leads to exponential decrease of the concurrence down to some residual values . this result shows that the concurrence is very sensitive to the onset of chaos , and understanding entanglement of these systems may help to control quantum chaos and suppress its negative effects in quantum information processing . we consider bipartite entanglement between a pair of qubits selected from a symmetric multi - qubit state and the rest of the system , as well as pairwise entanglement between the two qubits . once we obtain the two - qubit reduced density matrix , the entanglement can be readily calculated . by expressing the reduced density matrix @xmath34 for qubit 1 and 2 in terms of the expectation values of the collective operators , all elements of @xmath34 are conveniently obtained  @xcite . our system governed by the qkt hamiltonian is a composite system , and remains in a pure state at all times if we initially choose a pure state . for pure states , bipartite entanglement is well - defined and can be quantified by entropies of either subsystem . for convenience , we adopt the linear entropy as the entanglement measure , which is defined as @xmath35 where @xmath36 is the reduced density matrix for the first subsystem . the maximum linear entropy for a pure state of a bipartite @xmath37 system , is given by @xmath38 . while we may choose other entropies such as the von neumann entropy as our entanglement measure , the qualitative results are in general independent of choice of entropies for pure states . moreover , the linear entropy and the von neumann entropy are two limiting cases of the rnyi entropy  @xcite , they are thus interrelated and one can be used to estimate the other  @xcite . given our @xmath39-qubit system , we consider another type of entanglement , the pairwise entanglement , i.e. , the entanglement between a pair of qubits . when @xmath40 , the pair of qubits can be in a mixed state . entanglement for a mixed state @xmath41 is quantified by the entanglement of formation . specifically , for a pair of qubits , entanglement is equivalent to the non - positivity of the partially transposed density matrix  @xcite . alternatively , one can use the concurrence  @xcite to quantify the pairwise entanglement . the concurrence is defined as @xmath42with the quantities @xmath43 being the square roots of the eigenvalues in descending order of the matrix product @xmath44 @xmath45 denotes the complex conjugate of @xmath34 . the value of the concurrence ranges from zero for an unentangled state to unity for a maximally entangled state . we present here our studies of the entanglement dynamics of our @xmath39-qubit system governed by the qkt  ( [ top ] ) , with the relevant angular momentum operators @xmath46 given by the collective operators . if we choose the initial pure state to be symmetric under exchange of any qubits , then the state vector at any later time is also symmetric . thus , we can describe the state of the @xmath39-qubit system in terms of the orthonormal basis @xmath47 with @xmath48 . the states @xmath49 are the usual symmetric dicke states  @xcite . state @xmath50 is not entangled , whereas state @xmath51 , the so - called w state  @xcite , is pairwise entangled with concurrence @xmath52 .    to connect the quantum and classical dynamics of the kicked top , we choose the initial state to be the spin coherent state ( scs ) @xmath53 , with  @xcite @xmath54\}.\]]mean of @xmath55 is @xmath56 the initial scs can be rewritten as a multi - qubit product state , and thus exhibits no entanglement ( zero linear entropy and concurrence ) . we start by exploring the dynamics of entanglement for initial states with a mean value in four different regions of the phase space , specifically a fixed point , an integrable ( or kam ) region , a chaotic region , and the border between the integrable and the chaotic region known as ` the edge of chaos '  @xcite . we are also interested in the behavior as the chaoticity parameter @xmath0 is varied . we start with the choice @xmath57 , which exhibits large integrable and large chaotic regions , and we select four states localized in the four regions mentioned earlier .    for convenience we fix @xmath58 and vary @xmath59 : this ` line of latitude ' on @xmath26 includes all four regions we are exploring . an elliptic fixed point arises at @xmath60 , a point in the regular region occurs at @xmath61 , one edge of chaos can be seen at @xmath62 , and a point well in the chaotic sea is located at @xmath63 . the states with means at each of these @xmath64 points in the phase space are chosen to be scss . these states are minimum uncertainty states and are well localized around the four chosen points in phase space .     for @xmath31 and @xmath65.,scaledwidth=40.0% ] we numerically compute bipartite entanglement between two qubits and the other @xmath66 qubits , and display the numierical results in fig .  2 as linear entropy versus time increasing parameter @xmath67 . we observe that the entanglement is enhanced for the initial state centered in the chaotic region after a short time . initially , the linear entropy is zero , and as the dynamics evolve , the entropy increases slowly when the wave packet is centered in the regular region , whereas it exhibits a rapid rise when centered in the chaotic region . the curve with @xmath68 displays the intermediate behavior . furthermore , the entanglement for the state initially with mean at a fixed point displays a periodic modulation that is absent in the evolution of the entanglement for the state initiated in the chaotic region . this periodic modulation is an indicator of the underlying regular classical dynamics and corresponding regular energy spectrum of the floquet eigenstates  @xcite .        ) . the other parameters are the same as those of fig . 2 . , scaledwidth=40.0% ]    figure 3 shows the dynamical behavior of the concurrence . just as in the linear entropy case , we see a rapid change in the concurrence for a state initially centered in the chaotic sea . however , our numerical results suggest that the initial state centered in the fixed point will lead to large pairwise entanglement production , which is opposite to the case of the bipartite pure - state entanglement production , in which the classical chaos enhances the production of bipartite entanglement . of particular interest are the collapses and revivals in the evolution of the concurrence for a state centered on the fixed point . figure 3 ( left - up one ) shows a revival at @xmath69 . additional revivals occur at @xmath70 , @xmath71 and so on . for the chaotic case , we can not observe the revival phenomenon . this quasiperiodic behavior in the regular region indicates that the scs in the regular region has a finite support over the basis set of regular floquet eigenstates .    in fig . 4 we plot the long - time behaviors of the concurrence for the regular case , and observe multiple collapse and revivals . at long times the revivals become sparse , and finally the concurrence reduces to zero . we see that the concurrence is very sensitive to quantum chaos , which is consistent with the observations of bettelli and shepelyansky  @xcite in their studies of concurrence between qubits during the operation of an efficient multiqubit quantum algorithm . however , unlike their system in which the concurrence reached a finite residual value , in our qkt model , the concurrence disappears at very long times ; the difference arises because we assumes symmetrised multi - qubit states , whereas bettelli and shepelyansky allow this symmetry to be broken . the edge of chaos  @xcite is an important issue in the study of quantum chaos . in classical chaos , the edge of chaos is a fractal boundary separating the regular and chaotic regions . however , this fine - grained fractal structure does not translate well into the quantum domain . recently , it was found that the edge of chaos is characterized by a power law decrease in the overlap between a state evolved under chaotic dynamics and the same state evolved under a slightly perturbed dynamics  @xcite . here , we study the edge of quantum chaos from the perspective of entangling powers , which are defined to be either the maximal or the mean entanglement that the evolution operator can generate over all initial states  @xcite . alternatively , given a fixed initial state , we may ask what is the maximal and the mean entanglement that the operator can generate over all time . in general state - averaging and time - averaging are inequivalent and so the two methods yield different results .    in strongly chaotic systems , the two definitions converge due to nearly ergodic dynamics . in this study , we explore both the entanglement averaged over all time as well as the entanglement averaged over initial states . we begin our analysis with the average over time . in practice , for numerical purposes we consider a finite time domain . we define time - averaged entanglement power as the average linear entropy or average concurrence over a time interval @xmath72 ( which should be much longer than other time scales ) as follows : @xmath73 for local unitary operations , the above quantities are necessarily zero .    . the mean linear entropy is plotted for @xmath74 ( solid line ) , @xmath65 ( dashed line ) , and @xmath75 ( dotted line ) ; the mean concurrence is plotted for @xmath75 ( solid line ) , @xmath76 ( dashed line ) , and @xmath77 ( dotted line ) . parameters @xmath78 and @xmath0 are set to @xmath79 and @xmath57 . the time average is over 200 steps , and the subplots below are partly enlarged versions of above ones.,scaledwidth=40.0% ]    we fix the polar angle @xmath79 of the scs as before , and vary the azimuthal angle @xmath30 from @xmath80 to @xmath81 . the center of the scs wave packet thus commences in the chaotic region and passes through two regular islands . figure 5 displays the time - averaged mean linear entropy @xmath82and mean concurrence @xmath83 as a function of the azimuthal angle @xmath59 . when the azimuthal angle goes from @xmath80 to the first regular region , the linear entropy decreases until it reaches a minimum which approximately corresponds to the fixed point @xmath84 . subsequently the mean entropy increases to a flat larger area corresponding to the chaotic region . in contrast to the behavior of the mean linear entropy , the mean concurrence reaches a maximum approximately at the fixed point . when @xmath30 increases from @xmath85 , the mean concurrence first decreases slowly , and then exhibits an abrupt increase to a maximal value . thos turning point sharply displays the edge of chaos . another turning point is obvious from the figure . we also observe that the larger the number of qubits , the wider the regular region . in the large @xmath39 limit , the regular region will coincide with that of corresponding classical chaos of fig .  1 . we further calculate the mean linear entropy and the mean concurrence as a function of @xmath78 and @xmath59 . the contour plots are shown in figs .  4 and 5 . comparing figs .  1 and 4 , we observe that these two figures closely match each other . specifically , the four islands of fig .  4 are evident , reflecting the four stable islands in the classical phase space . comparing figs .  1 and 5 , the four stable islands of fig .  1 closely match those of fig .  5 . thus , we have a good classical - quantum correspondence .     and @xmath78 . the mean is over 200 kicks and other parameters are the same as those in fig . 2.,scaledwidth=40.0% ]     and @xmath78 . the average is over 200 kicks and other parameters are the same as those in fig . 2.,scaledwidth=40.0% ]       ( square line ) , @xmath86 ( diamond line ) , and the lyapunov exponent @xmath87 ( circle line)agaist @xmath0 . the plot of lyapunov exponent corresponds to fig .  1 of ref .  @xcite . the parameter @xmath88 and time average is over 50 kicks.,scaledwidth=40.0% ]    in the previous section , the time - averaged entanglement revealed clearly whether the initial state was in the regular or chaotic region . here , we are concerned with global properties of a chaotic system . for classical systems , the onset of global chaos can be quantified by calculating the global lyapunov exponent . here we define the following entangling power to quantify the onset of global quantum chaos : @xmath89 where @xmath90 is the haar measure . like the global classical lyapunov exponent , this state - averaged entangling power characterizes global properties of the qkt . the entangling power characterizes the entangling capability of the @xmath0-dependent floquet operator . figure 8 shows the entangling powers and the lyapunov exponent versus parameter @xmath0 . we observe that when @xmath91 , the entangling power @xmath92 exhibits a rapid increase and saturates beyond @xmath93 . the rapid increase of the entangling power signifies the onset of quantum chaos , and the saturation implies that global chaos has occurred . between @xmath94 and @xmath95 , when there are still regions of regular islands in the chaotic phase space , a mixture of regular and chaotic behavior is expected . in contract to @xmath92 , @xmath86 becomes very small for @xmath96 , which is also an indicator of onset of quantum chaos . note that @xmath86 has a peak which results from the competition between the entangling power of the qkt floquet operator and the inherent chaos . on one hand , increasing @xmath0 will enhance the entangling power , and on the other hand , the inherent quantum chaos suppresses the pairwise entangling power , thus leading to the peak . for @xmath92 , increase of @xmath0 and the quantum chaos both enhance the linear entropy , and thus no competition exists and no peak appears we have investigated a multi - qubit system whose collective hamiltonian dynamics are chaotic in the classical limit . we studied the particular example of the quantum kicked top , which is a well - studied example of quantum chaos with the advantages of having a finite - dimensional hilbert space ( thereby obviating the need for truncation that arises in infinite - dimensional hilbert spaces ) and of involving only spin operators of no more than quadratic order . although the collective dynamics are well understood , the underlying entanglement of the qubits that collectively make up the quantum kicked top is only just beginning to be understood . here we have developed methods for studying the quantum kicked top , and these methods are applicable to more general systems . we have identified bipartite and pairwise entanglement as two quite distinct measures to determine the entanglement in the system , and we have related the dynamics of these measures of entanglement to chaotic features of the quantum kicked top in the classical limit ; as examples , we have connected the features of the entanglement evolution to local properties such as whether the state is supported predominantly in the regular or chaotic region and also to global properties such as showing that entangling power averaged over states grows similarly to the global lyapunov exponent growth for the classical chaotic system . we have assumed symmetric multi - qubit states throughout , and the entanglement properties studied here reflect this assumption . if the symmetrization condition is broken , different dynamics can be expected . for example , bettelli and shepelyanski  @xcite show a concurrence that reaches a residual steady - state value . they explain this non - zero residue as being a result of symmetry breaking in their system . in contrast , our system exhibits a decay of concurrence to zero . the assumption of symmetric states implies indistinguishability of the qubits . thus even though entanglement may exist in the system , it may not be accessible as a useful tool for quantum information processing due to inherent inability to distinguish between the qubits . an alternative measure of entanglement could be an operational measure that takes into account physical restrictions on accessibility of the entanglement due to symmetries of the system  @xcite .    in summary , our work highlights the connection between entanglement of a multi - qubit state whose collective dynamics is chaotic in the classical limit and introduces valuable methods and measures for studying this entanglement . it would be worthwhile to investigate other quantum chaotic systems using concepts of time - averaging and global entangling power . furthermore , it would interesting to define and compare systems in which the symmetrization results hold to systems where this symmetrization condition is broken . we acknowledge valuable discussions with l.t . stephenson , x.w . hou , h.b . xie , l. yang , and h. zhang . this project has been supported by an australian research council large grant , the hong kong research grants council ( rgc ) , a hong kong baptist university faculty research grant ( frg ) , and alberta s informatics circle of research excellence ( icore ) . arrechi , e. courtens , r. gilmore , and h. thomas , phys . a * 36 * , 5543 ( 1987 ) ; a. perelomov , _ generalized coherent states and their applications _ ( springer - verlag , berlin , 1986 ) ; j.m . radcliffe , j. phys . * 4 * , 313 ( 1971 ) ;                                b. kraus and j.i . cirac , phys . a * 63 * , 062309 ( 2001 ) ; b. kraus , w. dr , g. vidal , j.i . cirac , m. lewenstein , n. linden and s. popescu , z. naturforsch . * 56a * , 91 ( 2001 ) ; j.i . cirac , w. dr , b. kraus , and m. lewenstein , phys . lett . *  86 * , 544 ( 2001 ) ; k. hammerer , g. vidal , and j.i . cirac , phys . a * 66 * , 062321 ( 2002 ) ; g. vidal , k. hammerer , and j.i . cirac , phys . * 88 * , 237902 ( 2002 ) ; w. dr , g. vidal and j.i . cirac , phys . rev . lett . * 89 * , 057901 ( 2002 ) ; x. wang and b.c . sanders , phys . a * 68 * , 014301 ( 2003 ) .
Abstract:

Answer: we explore the dynamics of entanglement in classically chaotic systems by considering a multiqubit system that behaves collectively as a spin system obeying the dynamics of the quantum kicked top . in the classical limit , the kicked top exhibits both regular and chaotic dynamics depending on the strength of the chaoticity parameter  @xmath0 in the hamiltonian . we show that the entanglement of the multiqubit system , considered for both bipartite and pairwise entanglement , yields a signature of quantum chaos . whereas bipartite entanglement is enhanced in the chaotic region , pairwise entanglement is suppressed . furthermore , we define a time - averaged entangling power and show that this entangling power changes markedly as  @xmath0 moves the system from being predominantly regular to being predominantly chaotic , thus sharply identifying the edge of chaos . when this entangling power is averaged over initial states , it yields a signature of global chaos . the qualitative behavior of this global entangling power is similar to that of the classical lyapunov exponent .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: in this paper we report on extensions to previous work on dynamical cavity qed effects in the photon statistics of transmitted light from a driven optical cavity coupled to an ensemble of two level atoms . much work has been done on structural cavity qed effects such as energy level shifts and the modification of spontaneous emission rates . these structural effects can be seen to arise from semiclassical models . in addition work has been done on dynamical effects where the coupling between the cavity field and atoms has a significant effect on the evolution of the system , in particular in the strong coupling regime where a single quantum of energy and hence single quantum fluctuations give rise to nontrivial dynamics . in this regime the field can neither be viewed as mildly perturbed by the atoms ( good cavity limit ) , nor are the atoms mildly perturbed by the field ( bad cavity limit ) . for a review of the work on structural and dynamical effects in cavity qed , see ref . @xcite . the problem of a single two level atom coupled to a single mode field was originally studied by jaynes and cummings and extended to many atoms by tavis and cummings . these models have been extended in recent theoretical work to include spontaneous emission and cavity field decay @xcite , and atomic transit time broadening and detunings @xcite . nonclassical correlations in photon statistics which violate a schwarz inequality have been predicted for this system , including photon antibunching [ defined here as @xmath0 , and sub - poissonian statistics [ @xmath1 . also other effects have been predicted which we refer to as overshoots and undershoots [ @xmath2 where @xmath3 is the normalized second order correlation function ] . examples of these nonclassical correlations from previous weak field results are shown in fig . [ fig : weak ] . figure [ fig : weak](a ) shows photon antibunching and sub - poissonian statistics , ( b ) shows an overshoot violation , and ( c ) shows an undershoot violation . photon antibunching has been seen experimentally in this system by rempe et . al . @xcite . overshoot violations have recently been seen by mielke , foster , and orozco @xcite . in general the theory matches the experiments in terms of qualitative behavior while the quantitative size of the nonclassical effect does not . this has led us to consider complications in the experiments which may be responsible for the discrepancy including deviations from the weak field limit , and dephasing due to the atoms entering and leaving the cavity . the experiments use an atomic beam to introduce atoms into the cavity so that the time of flight across the mode is on the order of ten spontaneous emission lifetimes @xcite . we expect that dephasing due to atomic traversal of the cavity will have a detrimental effect on nonclassical correlations . in addition deviations from the weak driving field limit and interactions with spectator atoms far from the mode waist may be important . these effects are investigated in this paper by numerically solving the master equation for the system and by quantum trajectory simulations . rather than investigating all possible effects at once , we isolate them and try to understand what is most critical . the general outline of this paper is as follows . in section [ sec : model ] we present the physical model of the system under investigation and describe the methods of solution . in section [ sec : strong ] we discuss the photon statistics of the transmitted light outside the weak field limit . section [ sec : dephasing ] presents two models of atomic transit dephasing and the resulting photon statistics . in section [ sec : two atoms ] we include effects of a spectator atom with a coupling which is a fraction of the maximum coupling , and finally we conclude in section [ sec : conclude ] . the system under investigation is an extension of the jaynes - cummings hamiltonian which includes effects of atomic and cavity field decay as well as a coherent driving field . a schematic diagram of the system is shown in fig . [ fig : diagram ] . the field and atomic hamiltonians are given by @xmath4 @xmath5 and the atom - field interaction in the rotating wave approximation is given by @xmath6 the cavity field creation and annihilation operators are @xmath7 and @xmath8 respectively and @xmath9 and @xmath10 are pauli operators for the @xmath11th two - level atom . the atom - field coupling strength is determined by @xmath12 where @xmath13 is the dipole transition matrix element between the two atomic states , @xmath14 is the cavity mode volume , and @xmath15 takes into account the position of the atom in the mode . in previous work it was assumed that the atoms sit at antinodes of the field where the coupling is a maximum . here we allow the atoms to be placed anywhere in the mode so that a range of couplings are allowed for different atoms . the cavity field is driven by a classical laser field with the driving field - cavity field coupling described by the hamiltonian @xmath16 where @xmath17 is the classical laser intensity scaled such that @xmath18 is the photon flux injected into the cavity . throughout we assume the atom , cavity field , and driving field are on resonance ( @xmath19 ) . dissipation in the system gives rise to nontrivial irreversible dynamics . cavity field damping and atomic population and polarization decay are described by superoperators acting on the density matrix of the system which are derived using standard methods @xcite . cavity field damping is described by @xmath20 where @xmath21 is the rate of cavity field damping . atomic population and polarization decay are described by @xmath22 where @xmath23 is the spontaneous emission rate of an atom . the full master equation in the born - markov approximation is then @xmath24 + { \cal l}_f \rho + { \cal l}_a \rho \to { \cal l } \rho .\ ] ] a numerical solution of the master equation is carried out in the fock state basis , and also a quantum trajectory simulation is developed from the master equation . the master equation in the fock state basis is    @xmath25\rho_{n,+;m,+ } \nonumber\end{aligned}\ ] ]    @xmath26    @xmath27\rho_{n,+;m,- } \nonumber\end{aligned}\ ] ]    @xmath28    where @xmath29 and @xmath30 and @xmath31 denote upper and lower atomic states respectively . we have numerically solved the master equation for the steady state for arbitrary driving field by truncating the fock basis at a point where the population of @xmath32 is less than @xmath33 . the second order correlation function @xmath34 is calculated from steady state matrix elements using the quantum regression theorem due to lax @xcite . we have developed a quantum trajectory simulation of this system from the master equation following the formalism of carmichael @xcite . we unravel the master equation into both a piece describing continuous evolution and a set of collapse operators in a way which is based on a simulated photon counting experiment . @xmath35 where @xmath36 is identified as the terms which can be written as commutators or anticommutators and @xmath37 is identified as all terms which can be written as @xmath38 . this particular unraveling is well suited for studies of photon statistics as the @xmath39 s represent quantum jumps due to emission of a photon . the continuous evolution of the system is described by @xmath40 while @xmath41 describes collapse events which punctuate the evolution . we define a closed system hamiltonian and a dissipative hamiltonian from the unraveled master equation as @xmath42 + \left [ h_d,\rho\right ] _ + \\ & = & -\frac{i}{\hbar}\left [ h_a + h_f + h_{af } +   h_l,\rho\right ] - \left [ \left ( a^\dagger a + \sum_j \sigma_+^j\sigma_-^j \right ) , \rho\right ] _ + \nonumber\end{aligned}\ ] ] where @xmath43_+$ ] denotes the anticommutator of @xmath39 and @xmath44 . a non - hermitian hamiltonian which reproduces the continuous evolution of the density matrix is defined as @xmath45 the rest of the master equation enters as collapse operators which are applied at random times when @xmath46 where @xmath47 is a random number between zero and one and @xmath48 the time step size is @xmath49 where @xmath50 is the fastest rate in the problem . in the event that we get two collapse processes in a single time step , we use a random number to choose one of the collapses . the time step is chosen so as to minimize such occurrences . for this system , the collapse operators are @xmath51 and @xmath52 corresponding to emission of a photon from the cavity field and the atom respectively . in section [ sec : dephasing ] , we describe dephasing due to an atom leaving the cavity using another collapse operator . because this unraveling of the master equation is based on photon counting experiments , the calculation of the second order correlation function is carried out quite naturally . the collapse operator @xmath51 corresponds to emission and detection of a photon from the cavity field mode . we calculate @xmath3 by building up a histogram of delay times between photon detections averaged over a long evolution time in a way analogous to experimental measurement . the photon statistics of the transmitted light have already been calculated in the weak field limit using a truncated five state basis where the system has up to two quanta of energy in it @xcite . the three types of nonclassical behavior previously discussed have been seen in subsequent experiments , however these experiments are not strictly in the weak field limit . it is of interest then to calculate the photon statistics for arbitrary driving field and to see to what extent the nonclassical effects persist . it is expected that for a strong enough driving field , the atoms saturate and the nonclassical photon correlations will be washed out because the cavity will basically contain a coherent state which is only mildly perturbed by the presence of the atom . looking at the photon correlations in the weak field limit from the point of view of quantum trajectories , we can interpret the nonclassical effects as resulting from the collapse of the wavefunction . the detection of the first photon emitted from the steady state collapses the wavefunction of the system ( @xmath53 ) and the subsequent time evolution as the system returns to the steady state determines the photon correlations . the second order correlation function is given by the probability of detecting a second photon normalized to the probability of detecting a photon in the steady state . here we start the system in the steady state , collapse the wavefunction , and let it evolve to get @xcite @xmath54 this assumes that there is usually only one photon emitted as the system returns to steady state which is a good approximation in the weak field limit . an example of this is shown in fig . [ fig : collapse ] for the case of the overshoot violation . outside of the weak field limit the photon correlation s are altered for two reasons . most simply , the time evolution following a collapse from steady state will be altered by the stronger driving field . another effect , however , is the presence of multiple collapses before the system returns to the steady state . consider a multiple collapse process . the first photon comes from the steady state and collapses the wavefunction of the system ( @xmath55 ) . now the time evolution occurs as before . however , the second photon collapses the system to a new state which depends on the delay time since the emission of the first photon [ @xmath56 . if a third photon is emitted before the system returns to steady state then its delay time will depend on the details of the evolution from @xmath57 . when averaged over many instances , this process will wash out the nonclassical effects because of the different evolution following different possible states , @xmath57 . an example of this process is shown in fig . [ fig:2photon ] where the conditioned photon number undergoes two collapse events . this figure shows three delay times of @xmath58 , @xmath59 , and @xmath60 with different evolutions of the conditioned photon number resulting in each case . figure [ fig : evolve ] shows the time evolution of @xmath61 following a photon emission from the steady state for a variety of system parameters . the overshoot persists in the evolution of the field following emission of a photon from the cavity for a driving field as large as @xmath62 . the undershoot and sub - poissonian statistics survive for driving fields as large as @xmath63 and @xmath64 respectively . the saturation field strength @xmath65 is the driving field for which @xmath66    the photon statistics of the transmitted field are shown in fig . [ fig : strong ] for the three types of nonclassical effects seen in this system at a variety of driving field intensities . figure [ fig : strong](a ) shows @xmath3 for system parameters ( @xmath67 ) which produce an overshoot violation of the schwarz inequality [ @xmath68 in the weak field limit . at a driving field of @xmath69 the overshoot violation is gone , thus this nonclassical effect is quite dependent on the weak driving field . figure [ fig : strong](b ) shows photon statistics for system parameters ( @xmath70 ) which produce an undershoot violation of the schwarz inequality [ @xmath71 in the weak field limit . here the nonclassical effect disappears at a driving field of @xmath72 showing that this is a more robust effect . figure [ fig : strong](c ) shows photon statistics for system parameters ( @xmath73 ) which produce photon antibunching [ @xmath74 and sub - poissonian statistics [ @xmath1 in the weak field limit . in this case the nonclassical effect persists until @xmath75 where the system shows slight bunching and super - poissonian statistics . ( notice that the nonclassical effects are not as robust as the time evolution of the cavity field would indicate . therefore the destruction of nonclassical effects are in part a result of multiple photon processes . ) for all system parameters the transmitted light becomes super - poissonian as the driving field is increased . we now turn our attention to the effects of atomic traversal of the cavity on the photon statistics . in previous work it has been assumed that the atoms are all fixed at antinodes of the cavity field and so had the maximum coupling @xmath76 . experiments on this system have used atomic beams to send atoms through a cavity . this atomic traversal of the cavity will introduce two new effects . first , the atom / cavity field coupling will depend on the position of the atom in the cavity . one might think that this would destroy the nonclassical correlations . however , the atoms with the largest coupling interact most strongly with the field and are most likely to contribute to the correlations . so the atoms near an antinode will have the largest contribution and other atoms may have little effect on the correlations . this issue will be further addressed in sec . [ sec : two atoms ] . the second effect of atomic traversal is dephasing which occurs when an atom enters or leaves the cavity . it is this effect which we consider in this section . we have used two approaches to model the dephasing due to atomic traversal . the first is to add a term to the master equation which describes nonradiative decay of atomic polarization . @xmath77 this term in the master equation has its origins in collisional processes @xcite and so may or may not accurately describe the dephasing which occurs when an atom leaves the cavity . the second approach uses a quantum trajectory simulation of the system to model the dephasing . in this approach we assume that there is always exactly one atom in the cavity . an atom leaves the cavity and another atom enters the cavity in the ground state at a rate @xmath78 . we can assume the atom enters the cavity in the ground state , but it is not immediately clear how to deal with the state of the exiting atom . this atom is in some superposition of excited and ground state and these states are entangled with the cavity field state . one approach would be to leave the photon number distribution of the cavity field unchanged using a collapse operator which has the following action on the state of the system : @xmath79 however , this is not a consistent application of the quantum trajectories . consider the evolution of the atom after it leaves the cavity . the atom at some later time may emit a photon into the vacuum meaning it was in the excited state when it left the cavity ; or it will never emit a photon , meaning it was in the ground state when it left the cavity . in general , the atom and environment and by entanglement the atom / cavity system will then be described by a density operator . however , we wish to use a pure state to describe the atom / cavity system conditioned on the detection of transmitted photons . to be consistent we must use a pure state to describe the atom after it has left the cavity . we use a collapse operator which picks either the excited state field distribution or the ground state field distribution of the system and then places the new atom in the ground state . this operator has the following action :    @xmath80    this collapse operator is then applied at a gaussian distributed series of times with average @xmath81 and a full width of @xmath81 as this approximates the traversal times of atoms with a maxwell - boltzmann velocity distribution . this model of dephasing differs from collisional dephasing in two important ways . first , it does not enter the deterministic hamiltonian evolution between collapses at all , whereas the collisional dephasing in the trajectory picture would have both a collapse component as well as the decay of coherence in the continuous evolution . second , this dephasing always places the atom into the ground state whereas the collisional dephasing places the atom in the ground or excited state with probabilities determined by the relative populations . now we turn to our results and a comparison of the two types of dephasing . because the traversal dephasing does not affect the deterministic evolution of the system we expect that for a given dephasing rate it will be less destructive of the nonclassical photon statistics than collisional dephasing . in fig . [ fig : dephasing ] we show the second order correlation function with collisional dephasing . all three types of nonclassical effects are quite sensitive to this dephasing with @xmath82 destroying the overshoot [ fig . [ fig : dephasing](a ) ] and undershoot [ fig . [ fig : dephasing](b ) ] while the sub - poissonian statistics survive until @xmath83 in fig . [ fig : dephasing](c ) . our results for the transit dephasing are shown in fig . [ fig : dephasing2 ] . the overshoot in fig . [ fig : dephasing2](a ) is again extremely sensitive to dephasing with classical statistics for @xmath84 while the sub - poissonian statistics in fig . [ fig : dephasing2](b ) are more robust , surviving up to @xmath85 . we now consider the effect of placing two atoms inside the cavity either both at antinodes of the field , or allowing one of the atoms to be arbitrarily placed so that its coupling is in the range @xmath86 to @xmath87 . dephasing is not considered in this section . in the experiments conducted on this system it is likely that there is some effect from spectator atoms which are located away from an antinode of the field and so do not contribute to the nonclassical photon statistics . if there are enough of these atoms they may simply absorb light and emit it out of the cavity thus effectively decreasing the quality of the cavity . as a first step towards understanding the effect of spectator atoms , we place one extra atom in the cavity with a coupling which is some fraction of the original atom s coupling . we use a quantum trajectory simulation to calculate photon statistics .    in fig . [ fig : spectator ] we plot @xmath3 while allowing the coupling of the spectator atom ( @xmath88 ) to vary from @xmath89 to @xmath87 . we see that for all three sets of parameters the photon statistics vary continuously from the single atom result to the two atom result with no qualitative deviation in the photon statistics . this suggests that the spectator atoms can have an observable effect on the statistics , but they do not destroy the nonclassical correlations at least when close to the ideal condition of having a single atom in the cavity at a time . a group of many spectator atoms may be more detrimental to nonclassical correlations . we now consider the case in which both atoms are placed at antinodes of the field . the photon statistics for this system have been solved for an arbitrary number of maximally coupled atoms using a set of symmetrized dicke states to describe the atomic excitation @xcite . this would correspond to an experimental setup where it is not possible to tell which atom spontaneously emitted a given photon . however , there would be experimental situations where the symmetrized states are not valid states with respect to spontaneous emission . here , using quantum trajectories , we consider both a symmetrized collapse and an unsymmetrized collapse for spontaneous emission . the operator for the dicke collapse is @xmath90 so that when a spontaneous emission event occurs , both atoms are collapsed symmetrically . for the non - dicke collapse the atomic operators are used separately so that one atom or the other collapses .    for weak driving field we expect that there will be no difference between these types of collapse because in this limit we only detect photons emitted from steady state . however , for stronger driving fields we begin to detect photons emitted from the collapsed state and the two collapses give a different collapsed state ( @xmath91 or @xmath92 versus @xmath93 ) , and therefore , different photon statistics .    in fig . [ fig : dicke ] @xmath3 is plotted for the two types of collapse for a driving field of @xmath94 . figure [ fig : dicke](a ) shows a significant difference as the nonclassical statistics are completely gone for the non - symmetrized collapse . figure [ fig : dicke](b ) shows no dependence on the type of collapse . figure [ fig : dicke](c ) shows a mild dependence on the type of collapse with a slightly larger value of @xmath95 for the non - symmetrized collapse . we have investigated extensions to previous theoretical work on a driven atoms / cavity system with dissipation . we have calculated the normalized second order correlation function for the transmitted light including effects of arbitrary driving fields , non - radiative dephasing , and arbitrary coupling strength for multiple atoms . we have found that nonclassical field states are easily destroyed by deviations from the weak field limit and by non - radiative dephasing modeled as both collisional dephasing and as atomic transit dephasing . we have also found that allowing two atoms in the cavity with different atom / field coupling strengths does not have a detrimental effect on the nonclassical field . the experiments which have been done on this system have not really in the weak field limit . as @xmath96 the number of counts also goes to zero so it is difficult to carry out experiments in this regime , but this work suggests that it is important . the authors would like to thank howard j. carmichael , luis orozco , and greg foster for many useful conversations .
Abstract:

Answer: we investigate the photon statistics of light transmitted from a driven optical cavity containing one or two atoms interacting with a single mode of the cavity field . we treat arbitrary driving fields with emphasis on departure from previous weak field results . in addition effects of dephasing due to atomic transit through the cavity mode are included using two different models . we find that both models show the nonclassical correlations are quite sensitive to dephasing . the effect of multiple atoms on the system dynamics is investigated by placing two atoms in the cavity mode at different positions , therefore having different coupling strengths .


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 presence of dust grains in the diffuse intergalactic medium is still an open question . if they exist , they should have great impacts not only on our understanding of galaxy evolution but also on interpretation of high - redshift observational data . from x - ray measurements of heavy metal lines , it is evident that clusters and groups of galaxies contain a significant amount of elements processed in galaxies ( e.g. , * ? ? ? * and references therein ) . such elements are likely to have been expelled out of galaxies through ram pressure stripping , galaxy mergers and interactions , blowout by supernovae , etc . dust may also be expelled by these processes or by some independent mechanism such as radiation pressure ( e.g. , * ? ? ? * ; * ? ? ? the x - ray emitting hot gas can collisionally heat the dust @xcite , while it can also destroy small dust grains via sputtering @xcite . dust grains coexistent with the x - ray emitting gas have indeed been detected in a number of elliptical galaxies ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? the temperature , the amount , and the size distribution of the dust in hot plasma can be quite different from those in the normal interstellar medium ( see dwek & arendt 1992 for review on dust  gas interactions ) . once detected firmly , the intracluster or intragroup dust will definitely provide a powerful probe of galaxy evolution and the dust  gas interaction . there have been several suggestions and debates regarding the presence of dust in clusters and groups of galaxies . over five decades ago , @xcite first suggested an enhanced extinction toward the coma cluster of galaxies . more recently , @xcite pointed out that the redshift asymmetries of member galaxies in nearby groups are consistent with the presence of dust in the intragroup medium . the observed oxygen k edge in an x - ray spectrum of the perseus cluster can also be attributed to intracluster dust grains @xcite . the extended submillimeter emission detected in a rich galaxy cluster may partly be due to dust @xcite . there have been extensive searches for the enhanced visual extinction toward high redshift objects behind clusters @xcite and @xcite reported the detection of reddening toward @xmath5 clusters . on the other hand , the intracluster dust is likely to enhance the efficiency of gas cooling @xcite that is in apparent contradiction with recent x - ray observations @xcite . more direct evidence has been searched for in the far - infrared bands . @xcite reported that the far - infrared emission observed with _ iras _ is enhanced by about a factor of 2 in compact groups of galaxies compared with a sample of isolated galaxies . the extended excess emission is also inferred toward abell 262 and abell 2670 @xcite . @xcite , however , pointed out that the results of @xcite are likely to have been overestimated mainly due to the limited spatial resolution of _ @xcite and @xcite reported the detection of extended excess far - infrared emission from the coma cluster , which can be attributed to thermal emission from the intracluster dust . they find such emission only in one out of six clusters they observed with _ infrared space observatory ( iso)_. @xcite combined the _ iras _ maps toward 11,507 galaxy clusters and statistically detected the excess emission , although the contribution of member galaxies is yet to be fully removed . major difficulties in the previous infrared searches are the limited sensitivity and spatial resolution of the detectors . these limitations may have led to significant contamination of individual galaxies as well as galactic cirrus confusion . it will therefore be meaningful to perform infrared observations using current facilities with much improved sensitivity and spatial resolution . @xcite searched for the intracluster dust emission from abell 2029 with _ spitzer _ and reported that the cirrus noise limits the detection of the intracluster dust component .    in this paper , we present the results of our deep observations toward the central region of the coma cluster ( abell 1656 ) with the multiband imaging photometer for _ spitzer _ ( mips ) on board the _ spitzer _ space telescope . coma lies at high galactic latitude ( @xmath6 ) where the cirrus level is expected to be among the lowest . in fact , @xcite and @xcite found with _ iso _ the excess far - infrared emission toward coma based on the brightness ratio @xmath7 , while the s / n at each wavelength is still insufficient to detect the excess directly and the contamination of galaxies is not removed . we attempt to probe the excess emission more directly by approaching the confusion limits of extragalactic sources and galactic cirrus in each band .    throughout the paper , we assume a standard set of cosmological parameters , @xmath8 , @xmath9 , and @xmath10 consistent with the _ wmap _ five year data @xcite . an angular size of 1@xmath11 corresponds to a physical size of 0.46 kpc at the redshift @xmath12 of the coma cluster . the central region of the coma cluster was observed at 24 , 70 and 160 @xmath0 m with the mips in our go3 program on 2007 january 19 . total duration of the observation was 41 ksec . linear scans were performed in a slow scan mode covering the @xmath13 region across the cluster as shown in figure [ fig : pspc ] ( the scan width was reduced to @xmath14 at 70 @xmath0 m due to dead pixels ) . the scan paths were carefully chosen so as to avoid bright galaxies and to lie close to the cluster center ; the scan center was fixed at the mid - point ( @xmath15 ) between two massive central galaxies , ngc4874 and ngc4889 with the position angles of @xmath16 , @xmath17 and @xmath18 at 24 , 70 and 160 @xmath0 m , respectively . the above scan center is @xmath19 northeast of the x - ray center of the cluster in the _ rosat _ pspc image ( @xmath20 ) . this separation is small compared to the extent of the x - ray emission , which is characterized by the core radius of @xmath21 . in total , 32 scans were performed in each band , 16 of which repeatedly covered a @xmath22 region on either side of the scan center . the integration time per pixel after 16 scans was 1600 , 1600 and 160 sec at 24 , 70 and 160 @xmath0 m , respectively . two sets of the 16 scan legs overlapped and the integration time was doubled in the central @xmath23 region , where the x - ray emission is the most prominent . the data were first processed with the mips pipeline , version s17 @xcite to yield the basic calibrated data ( bcd ) . to further improve the images , the pixel - dependent and scan - dependent backgrounds were corrected as follows . in doing so , special care was taken to maintain the extended nonlinear component of the emission . all the results regarding the extended emission are based on unfiltered images . first , we create a median flat field of detector pixels by stacking the normalized bcd frames , excluding the regions within @xmath21 from the cluster center . the normalization is done by dividing each bcd by its median . this is intended to separate any large - scale variations , e.g. , due to the zodiacal light , from the pixel - dependent backgrounds and applied solely to create the flat field . we have checked that the flat field so obtained is properly normalized ; its mean is equal to unity within 0.2% . we then divide the original bcd ( not normalized ) by the normalized median flat field . at 24 @xmath0 m , we further remove scattered light that produces a background modulation as a function of the scan mirror position ; the average offset at each scan mirror position is computed excluding the data within @xmath21 from the cluster center . as mentioned below , these processes will be re - run after the source identification . second , we create for each band a mosaic image using the mopex software package , version 18.2 @xcite developed at the spitzer science center . the pixel size is taken to be a nominal value of @xmath24 , @xmath25 and @xmath26 at 24 , 70 and 160 @xmath0 m , respectively . we exclude the first five bcd frames in each scan leg to omit boost frame transients . we also exclude the frames affected by stimulator flashes at 70 and 160 @xmath0 m . the background matching is applied _ after _ combining the frames in each scan leg ; we only correct the constant offsets among 32 scan legs to assure sufficiently large overlapping area and to preserve any extended feature of the emission . third , we identify point sources in the mosaic images down to the 1@xmath1 detection limit ( @xmath27 , @xmath28 and @xmath29 mjy at 24 , 70 and 160 @xmath0 m , respectively ) using apex @xcite . this limit corresponds to the @xmath30 confusion level in each band @xcite . the regions around the sources are masked out and unused in our analysis . more specifically , the source identification and the masking procedure are done on the basis of median - subtracted images ; a median - filtered image is created and subtracted from the unfiltered image solely for this purpose . to improve the efficiency of source detection against unknown background fluctuations , we vary the size of a median filter box from @xmath31 to @xmath32 pixels . for every source identified in each median - subtracted image , we estimate the radius @xmath33 at which the surface brightness , including that of the side - lobes , drops below a certain threshold noise level using the point response function . we adopt 0.5@xmath1 as the threshold noise level at 24 and 70@xmath0 m , and 1@xmath1 at 160 @xmath0 m , considering the lower uniformity of the data in the latter band . we then return to the _ unfiltered _ image and mask all the pixels within @xmath33 from the sources identified in any of the median - subtracted images . we have checked that the correlation coefficient between the mean brightness of masked pixels ( @xmath34 ) and that of unmasked pixels at @xmath35 is less than 0.2 in any of the median - subtracted images . we thus consider that our masking procedure successfully removes the contribution of resolved sources .    strictly speaking , the source identification and the background correction are not independent of each other . we therefore repeat the above procedure from the flat - field correction excluding the pixels around the point sources in the bcd . after applying the above corrections , we subtract from each unfiltered image the linear gradient over the entire map , so as to remove the variation on scales larger than the galaxy cluster . such a gradient is greatest at 24 @xmath0 m due to the zodiacal light ; it causes the difference of @xmath36 mjy / sr between the map edges which is more than 2 orders of magnitude larger than the average noise level of the image ( table [ tab : images ] ) . by simply subtracting the linear component , we avoid relying on specific models of the zodiacal light emission that are still uncertain at scales probed here . to preserve the average brightness of the entire map , the brightness at the scan center is kept unchanged by the subtraction . figure [ fig : images ] shows the mosaic images after the subtraction of the linear gradient . the point sources are explicitly shown for display purposes . the pixel sizes of the images , the fwhms of the point response function , and the average noise levels are listed in table [ tab : images ] . we have assigned a noise to an individual pixel as the larger of the uncertainty estimated in the mips pipeline and the standard deviation of actual measurements at that pixel . the latter is on average twice as large as the former at 70 and 160 @xmath0 m , and both are comparable at 24 @xmath0 m . this choice should enable us to take account of any residual systematics in each pixel more properly . note that the sensitivity at 160 @xmath0 m is likely limited by galactic cirrus ; the @xmath37 fluctuation level expected at the position of coma is @xmath38 mjy / sr at 70 @xmath0 m and 0.20 mjy / sr at 160 @xmath0 m @xcite .    masking the regions around the detected sources removes 38% , 66% and 64% of the pixels in the 24 , 70 and 160 @xmath0 m images , respectively . we plot the histograms of the pixels with and without masking in figure [ fig : hist ] . in all bands , the histograms before masking are positively skewed due to bright sources . once the sources are excluded , the data follow the gaussian distribution more closely , while there are small excesses at high @xmath1 . these excesses are due to low - significance data affected by noisy detector arrays or low coverage . since the noise assigned to such data is large , their contribution to the mean brightness is negligible . at 160 @xmath0 m , the 1@xmath1 range of the distribution after the source mask , 0.19 mjy / sr , is likely affected by the fluctuations of galactic cirrus emission @xcite . at 24 and 70 @xmath0 m , the 1@xmath1 ranges , @xmath39 and @xmath40 mjy , are comparable to the average noise level listed in table [ tab : images ] .    figures [ fig : prof1][fig : prof3 ] show the surface brightness distribution along the scan path as a function of the distance from the x - ray center of the cluster ( negative values of the distance are assigned to the pixels at lower declination than the x - ray center ) . we plot the brightness at each pixel and their binned average . unless otherwise stated , the error bar indicates the 1@xmath1 statistical error of the mean in each bin . the 70 @xmath0 m data are lacking within @xmath41 from the x - ray center due to dead pixels . we hence adopt fiducially the bin size of @xmath42 ( 100 kpc ) in the present paper so that it fully covers the deficit region while it is still sufficiently smaller than the x - ray core radius of @xmath21 . in all bands , the fluctuations in the brightness , particularly the spikes due to point sources , are significantly reduced after the source mask , yet there is no indication of the central excess .    in the bottom panels of figures [ fig : prof1][fig : prof3 ] , we plot the mean brightness of the removed component , i.e. , the difference between the unmasked and masked brightness . from our linear scans within @xmath43 from the cluster center , we do not find a central concentration of the integrated source intensities either . the average levels of the removed brightness are @xmath44 mjy / sr , @xmath45 mjy / sr and 0.16 mjy / sr at 24 , 70 and 160 @xmath0 m , respectively . we further combine the brightness in the two directions from the cluster center and compare its radial profile with model predictions in figure [ fig : prof_abs ] . the prediction for the intracluster dust emission is based on the theoretical model of @xcite ( see the next section for details ) . for direct comparison with the _ iso _ results , we have adjusted the total amount of dust to reproduce the excess brightness of 0.2 mjy / sr at 120 @xmath0 m within the radius of @xmath46 inferred in @xcite . since we are only interested in the difference of the signal across the cluster , an offset is added to the model brightness to match the mean brightness of the whole map at the edge . at 24 @xmath0 m , the absence of the excess is as expected since the predicted temperature of the intracluster dust is relatively low , @xmath47k ( see also dwek et al . 1990 ) . at 70 and 160 @xmath0 m , the foreground / background fluctuations are comparable to the expected levels of the emission and we do not find a sign of the central excess in either band .    from the observed brightness within @xmath48 relative to the mean of the entire map , we obtain 2@xmath1 upper limits on the central excess as @xmath49 , @xmath50 and @xmath51 mjy / sr , at 24 , 70 and 160 @xmath0 m , respectively . in deriving these limits , we have taken into account the background / foreground fluctuations that dominate over the statistical noise as follows . we have randomly selected a circular region of radius @xmath42 from the map in each band 100 times . from the variance of the mean brightness of selected regions , we have assigned an @xmath37 error of @xmath52 , @xmath53 and @xmath54 mjy / sr at 24 , 70 and 160 @xmath0 m , respectively . @xcite and @xcite found with _ iso _ the excess far - infrared emission toward the center of this cluster based on the surface brightness ratio @xmath7 . although the scan paths and the observed bands are not exactly the same in our observations , we attempt to mimic their analysis by presenting @xmath55 and @xmath56 in figure [ fig : rat ] . while the fluctuation is large , there appears to be a weak sign of a central enhancement in @xmath56 . the enhancement is largely due to a variation in the 160 @xmath0 m map around the cluster center ( fig.[fig : images ] ) ; there are bright extended regions at @xmath57 on each side , presumably produced by far - infrared sources lying on fluctuating cirrus foreground . we have checked that no counterparts are detected in the _ chandra _ x - ray image of this region . while some of the sources are apparent and removed at 70 @xmath0 m , the lower spatial resolution and the larger cirrus fluctuation prevent us from removing them completely at 160 @xmath0 m . in turn , the cluster center appears to be darker than its surrounding at 160 @xmath0 m , resulting in the apparent excess in @xmath56 . figure [ fig : prof3 ] also indicates that the depression at 160 @xmath0 m is not due to over - removal of the sources at the center . we therefore do not consider this to be a detection of the intracluster dust . while there is no direct evidence of the central excess in any of the three bands observed by the mips , we still find a weak sign of excess in the brightness ratio @xmath56 that may resemble the excess in @xmath7 reported by @xcite and @xcite . it is therefore meaningful to explore the consistency of these results .    for this purpose , we subtract from our data the average zodiacal light and the cosmic infrared background levels estimated by the _ spitzer _ background estimator @xcite at the scan center position for the date of observation ; @xmath58 and @xmath59 mjy / sr , @xmath60 and @xmath61 mjy / sr at 70 and 160 @xmath0 m , respectively . we subtract the above values from the whole map since we have already corrected the linear gradient of the map relative to the scan center . we have checked that subtracting position - dependent zodiacal light levels prior to the correction of the linear gradient changes the resulting brightness ratio by no more than 5% . the uncertainties in the absolute calibration are not significant either ; @xmath62 at 70 @xmath0 m and @xmath63 at 160 @xmath0 m @xcite . however , simply because the mid - infrared brightness is dominated by the zodiacal light , the nominal error of @xmath64 in its model @xcite will change the brightness after subtraction at 70 and 160@xmath0 m , @xmath65 mjy / sr and @xmath66 mjy / sr , by @xmath67% and @xmath68% , respectively . keeping these large uncertainties in mind , we investigate the origin of the excess in the band ratio as follows .    figure [ fig : rat_abs ] shows the radial distribution of the brightness ratio @xmath69 after the correction mentioned above . the error bars include only the statistical uncertainties . assuming that the corrected brightness distribution is produced by a combination of a central excess and more extended component , e.g. , galactic cirrus , we apply equation ( 7 ) of @xcite to our observed bands and obtain @xmath70 mjy / sr within the radius of @xmath42 . this value is comparable to the 2@xmath1 upper limit of @xmath71 mjy / sr derived solely from the 70 @xmath0 m data in the previous section . the central excess , if real , should have been nearly canceled at 70 @xmath0 m by the fluctuation of the other component . the interpretation of the central excess for the origin of the observed @xmath56 profile is therefore rather unlikely , although it is not fully rejected . note that equation ( 7 ) of @xcite is based on a key assumption that the excess component contributes only to the variation of the shorter wavelength data . alternatively , it is also possible to interpret the same data with the variation of the longer wavelength data . in the present observation , a plausible explanation may be that the corrected brightness is still dominated by unresolved sources lying on the fluctuating cirrus foreground , i.e. , the absolute values of the signal in both bands are dominated by such sources while the fluctuations at 160 @xmath0 m are mainly produced by underlying cirrus emission . in fact , the value of @xmath72 shown in figure [ fig : rat_abs ] corresponds to the dust temperature of @xmath73 k for the absorption efficiency of @xmath74 @xcite . this temperature is higher than the nominal temperature of galactic cirrus @xmath75 k for which the brightness ratio is lower by more than a factor of 5 . the mean values of @xmath76 mjy / sr and @xmath77 mjy / sr are also higher than the levels of the cirrus emission expected for this region , suggesting the contribution of yet unresolved components .    from the limits on the infrared brightness in each band , we further attempt to derive constraints on the amount of the intracluster dust . model predictions shown in figure [ fig : prof_abs ] @xcite rely on an assumption that the dust is supplied continuously to the intracluster medium from galaxies with the initial size distribution similar to our galaxy , @xmath78 @xcite , and has reached a steady state against sputtering with the size distribution @xmath79 . the galaxy distribution is assumed to trace that of the dark matter given by the navarro  frenk white density profile @xcite , and the gas density follows the isothermal beta model . given these , the curves plotted in figure [ fig : prof_abs ] correspond to the dust surface density of @xmath80 m@xmath81 kpc@xmath82 within the central @xmath42 of the cluster . the 2@xmath1 upper limits on the excess brightness derived in the previous section are 1.3 and 0.37 times the model brightness averaged within the central @xmath42 at @xmath83 m and @xmath84 m , respectively . considering much lower uniformity of the data at 160 @xmath0 m , we take the limit of the @xmath83 m data to obtain @xmath85 m@xmath81 kpc@xmath82 . this limit corresponds to the dust - to - gas - mass ratio of less than @xmath86 within the central 100 kpc , which is nearly 3 orders of magnitude smaller than the galactic value @xcite . the average visual extinction @xmath87 is related to the dust surface density as @xcite @xmath88 where @xmath89 is the extinction efficiency factor , @xmath90 is the dust mass density , and @xmath91 is the grain radius . adopting the standard values for @xmath89 , @xmath90 and @xmath91 quoted in the above equation , our derived limit of @xmath92 m@xmath81 kpc@xmath82 corresponds to @xmath93 within the central 100 kpc . on the other hand , @xcite report the detection of the mean @xmath94 of a few times @xmath95 mag toward sdss qsos behind @xmath96 clusters , and @xcite obtain the limit @xmath97 ( 2@xmath1 ) from reddening of @xmath98 background galaxies for 458 clusters . while these optical measurements are given for much larger ( @xmath99 mpc ) scales , our result is still consistent with them . finally , gravitational lensing can impact on the emission of background sources behind galaxy clusters . as the total surface brightness is conserved by lensing , the bias in the resolved or unresolved brightness is caused solely by the sources that exceed the detection limit as a result of magnification @xcite . the number of such sources is @xmath100 times the number of unlensed sources already lying above the flux limit @xmath101 , where @xmath0 is the magnification and @xmath102 is the index of the unlensed source count , @xmath103 . the fractional change in the resolved brightness has the same sign as and is less than @xmath104 , because the sources contributing to the change have the lowest flux among those resolved . note that the factor @xmath104 is different from the bias in the overall surface density given by @xmath105 ( e.g. , * ? ? ? recent surveys with _ spitzer _ yield the mean values of @xmath106 down to 35 @xmath0jy at 24 @xmath0 m @xcite and @xmath107 down to 1.2 mjy at 70 @xmath0 m @xcite . the depths of these measurements are comparable to those of the present observation . at 160 @xmath0 m , the published counts are currently available for @xmath108 mjy , and an extrapolation using the model of @xcite yields @xmath109 at @xmath110 mjy @xcite . in the case of coma , the einstein radius is @xmath111 and the magnification is expected to vary from @xmath112 at @xmath48 to @xmath113 at @xmath114 for a singular isothermal mass profile . taken together , the above factor in this radius range varies from @xmath115 to @xmath116 at 24 @xmath0 m and 70 @xmath0 m , and @xmath117 to @xmath118 at 160 @xmath0 m . the fractional change in the brightness of the _ removed component _ shown in figures 46 should therefore be smaller than the variations along the scan path . there is in fact no sign of systematic increase in the resolved brightness toward the cluster center in any band . owing to the conservation of the total brightness , the bias in the residual unresolved background should also be small . we hence consider that the impact of gravitational lensing is negligible in the present results . in this paper , we have undertaken a search for the infrared emission from the intracluster dust in the coma cluster of galaxies by the mips on board the _ spitzer _ space telescope . our observations yield the deepest mid and far - infrared images of a galaxy cluster ever achieved with average 1@xmath1 sensitivities per pixel of @xmath119 mjy / sr , @xmath120 mjy / sr and @xmath121 mjy / sr at 24 , 70 and 160 @xmath0 m , respectively . in each of the three bands , we have not detected a signature of the central excess component . our observations yield the 2@xmath1 upper limits on the excess emission within @xmath42 ( 100 kpc ) of the cluster center as @xmath49 mjy / sr , @xmath122 mjy / sr , and @xmath54 mjy / sr , at 24 , 70 , and 160 @xmath0 m , respectively . these results are consistent with the values reported by @xcite for abell 2029 , as well as those derived from a stacking analysis of 11,507 clusters @xcite .    with the improved spatial resolution and sensitivity of the mips , we have obtained much more robust limits on the excess emission from this cluster than previous observations . we still find that the brightness ratio between 70@xmath0 m and 160@xmath0 m shows a marginal sign of the central excess , in qualitative agreement with the finding of @xcite and @xcite by _ iso_. although we can not fully reject the contribution of the intracluster dust , it is more likely due to unresolved infrared sources lying on fluctuating cirrus foreground . combined with the theoretical model of @xcite , we have derived a limit on the surface mass density of the intracluster dust as @xmath92 m@xmath81 kpc@xmath82 . this limit corresponds to the dust - to - gas - mass ratio of less than @xmath123 within the central 100 kpc , which is nearly three orders of magnitude smaller than the galactic value , implying that the dust destruction is very efficient at the center of the cluster . our result can further be converted to the limit on the visual extinction of @xmath124 , which is consistent with the results derived from statistical reddening of background objects behind a large sample of galaxy clusters @xcite . we thank takashi hamana , nobuhiro okabe and norio sekiya for discussions and an anonymous referee for helpful comments . this work is supported in part by grant - in - aid for young scientists by mext ( 18740112 , 19740112 ) . n.o . acknowledges support from the alexander von humboldt foundation . cccc +   + wavelength & 24 @xmath0 m & 70 @xmath0 m & 160 @xmath0 m +   + pixel size [ @xmath11 ] & 2.45 & 4.0 & 8.0 + fwhm [ @xmath11 ] & 6.2 & 18 & 40 + average noise per pixel [ mjy / sr ] & @xmath126 & @xmath120 & @xmath121 +
Abstract:

Answer: we have undertaken a search for the infrared emission from the intracluster dust in the coma cluster of galaxies by the multiband imaging photometer for _ spitzer_. our observations yield the deepest mid and far - infrared images of a galaxy cluster ever achieved . in each of the three bands , we have not detected a signature of the central excess component in contrast to the previous report on the detection by _ infrared space observatory ( iso)_. we still find that the brightness ratio between 70@xmath0 m and 160@xmath0 m shows a marginal sign of the central excess , in qualitative agreement with the _ iso _ result . our analysis suggests that the excess ratio is more likely due to faint infrared sources lying on fluctuating cirrus foreground . our observations yield the 2@xmath1 upper limits on the excess emission within 100 kpc of the cluster center as @xmath2 mjy / sr , @xmath3 mjy / sr , and @xmath4 mjy / sr , at 24 , 70 , and 160 @xmath0 m , respectively . these values are in agreement with those found in other galaxy clusters and suggest that dust is deficient near the cluster center by more than 3 orders of magnitude compared to the interstellar medium .