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: directional detection is one of the research frontiers on wimps and several techniques are being studied and developed . in a recent paper  @xcite , it has been proposed that large arrays of aligned carbon nanotubes ( cnts ) might be of use for the directional detection of dark matter ( dm ) weakly interacting massive particles ( wimps ) .    as a result of a wimp scattering on the surface of a cnt , a carbon ion , with energy up to few tens kev , might be extracted and propagate through multiple collisions within the cnt empty volume , provided that its initial velocity is sufficiently aligned with the cnt axis : a sort of _ ion channeling _ phenomenon . however , in a large aligned cnt array , the ion might be equally channeled by the interstices among different nanotubes , as mentioned in  @xcite . the purpose of this paper is to illustrate the dynamics in these regions and , by means of monte carlo ( mc ) simulations , to provide an evaluation of the efficiency of a large number of cnt arrays as the target of a directional dm detector carbon nanotubes each , are the building blocks of our detector . a large number of them must be packed , to obtain a detector surface of the order of 100 m@xmath1 . ] .    at first , we consider particles satisfying channeling conditions . a continuum repelling potential , from the collective structure of atoms on the cnt surface , can be used to describe the channeling kinematics  @xcite . for nanotube radii much larger than the typical c electron orbital size , the potential can be well approximated by a step barrier . in this approximation , cnts are considered as full solid cylinders , totally reflecting those ions having the right initial conditions , namely a transverse energy lower than the continuum potential barrier ( * ? ? ? the ion motion in the transverse plane is that of a particle bouncing in a 2d lattice of circular obstacles . longitudinal energy and velocity along cnts axes are conserved . ions are assumed to move freely between consecutive reflections and their trajectories are a combination of broken lines  @xcite . given these initial assumptions , we can validate our monte carlo ( mc ) simulations on known results on classical billiards  @xcite .    in a second step , the initial conditions of the extracted ions are computed from the recoiling ion momentum distribution @xcite . when channeling conditions are not met , the scattered ion is allowed to penetrate the cnt and it is then subject to larger energy losses @xcite and angular deviation effects @xcite , due to the interactions with atoms at the cnt surfaces . the results obtained here give more strength to the proposal presented in  @xcite and give important suggestions on how to devise the cnt arrays to maximize the wimp detection . although we are still in the stage of a conceptual design and theoretical considerations , the outcome of this work is another step towards the testing of the prototype envisaged in  @xcite . with really small expected rates , the challenge of dark matter detection experiments is to distinguish an unambiguous signal from a large background . direct detection experiments identify the recoils produced by incident wimps on the detector target nuclei . the signal is subject to an annual modulation , related to the variation of the earth s relative motion with respect to the dm galactic halo . up to now , only the dama experiment @xcite claimed the detection of an annual modulation with an undeniable statistical significance . cresst @xcite and cogent @xcite also reported anomalous results , but their consistency with dama ones is a vexed question .    moreover , these results are in conflict with those from cdms @xcite , xenon10 @xcite , zeplin - ii @xcite and zeplin - iii @xcite , which , on the contrary , found no evidence of annual modulation . directional detection of dark matter was proposed in @xcite , where it was observed that a strong forward - backward asymmetry is expected in the case of an isothermal spherical galactic halo . because of this , directional detection , even if in the case of poor angular resolution , might also provide information on the apparent wimp wind direction . several experiments , including drift @xcite , newage @xcite , mimac @xcite , dmtpc @xcite and news @xcite , are focused on this goal . each of them is based on specific detection techniques . drift , newage , mimac and dmtpc are dm directional detectors using as target materials low pressure gas mixtures . nuclear recoils are identified by measuring the induced ionization with a time projection chamber ( tpc ) technique . dmtpc set a @xmath2 cl upper limit on the spin - dependent wimp - proton cross section of @xmath3 cm@xmath1 for 115 gev wimps @xcite , drift s exclusion curve for spin - dependent wimp - proton interactions reaches 1.1 pb at 100 gev @xcite , while mimac s reaches 5400 pb for 150 gev wimps . news proposes to use nuclear emulsions as both the wimp target and the tracking detector by reconstructing the direction of the recoiled nucleus . finally , there is also a proposal to use dna or rna , so that the path of the recoiling nucleus can be tracked to nanometer accuracy @xcite .    in  @xcite , a new type of directional detector was outlined . it consists of large arrays of cnts which provide the target material for wimp - nuclei collisions and an anisotropic sensitivity to the recoil direction , complemented by a low pressure tpc - based readout technique to detect the emitted nuclei . given the relatively low mass of the c nuclei , such detector would extend the wimp search to low masses ( @xmath0  gev ) . as a first approximation , we consider the cnt system to be a regular array of @xmath5 cylindrical obstacles and assume that the particles can not penetrate them . elastic collisions are also assumed .    [ cols="^,^ " , ]     let @xmath6 be the height of the cylindrical obstacles , @xmath7 the radius of their circular basis ( see fig . [ fig : primitive - cell ] ) and assume the cylinder axes are parallel to the @xmath8 direction ( the _ longitudinal _ direction ) , in the @xmath9 half - space , their bases being on the @xmath10 plane ( the _ transverse plane _ ) . we call , respectively , the @xmath10 plane and its parallel plane at @xmath11 the _ bottom _ and the _ top _ ends of the system . assume that the bottom bases centers of the @xmath12 cylinders form a square regular array on the @xmath10 plane . let @xmath13 be the distance between the centers of the bases of two neighboring cylinders and @xmath14 the minimal width of the interstices ( see fig . [ fig : primitive - cell ] , left panel ) . consider a particle entering the lattice through the bottom plane , thus with longitudinal velocity component @xmath15 , and let @xmath16 be the modulus of the projection of the velocity on the transverse plane . assuming uniformly distributed positions and transverse velocity directions as initial conditions at the bottom end , we estimate the typical time needed by a particle to exit the array from the top end . the 2d motion of the system is equivalent to the motion of a particle inside an `` infinite - horizon '' periodic billiard table : a billiard in which there are corridors running straightforwardly towards the perimeter . the average time spent by a particle inside the billiard can be estimated by using the so called machta - zwanzig approximation @xcite , which holds in the limit @xmath17 ( as for the finite - horizon case ) . we assume therefore that the particle spends in each elementary cell ( see fig . [ fig : primitive - cell ] , right panel ) a random time exponentially distributed with mean @xmath18 and then jumps with equal probability @xmath19 to any of the four neighboring cells ) , see app . [ taur ] . ] . hence , we can estimate the typical total time @xmath20 spent by a particle in the billiard by considering a two - dimensional homogeneous random walk on a square of @xmath5 points in @xmath21 . the walker starts on the square from a uniformly distributed random point . we then evaluate the mean number of hops @xmath22 performed by the walker to touch the external boundary of the square . the typical time needed by the particle to exit the billiard is then given by @xmath23 we evaluate @xmath22 numerically by averaging the total number of jumps performed by each particle in an ensemble of @xmath24 particles . by interpolating the results as a function of the side length @xmath25 , we find @xmath26 . the quadratic dependence on the side length of the square is quite obvious , since the random walk is symmetric , while the precise value of the numerical prefactor is not .    in the following , we briefly discuss the algorithm we used in mc simulations for an infinite horizon billiard . the semi - analytical expression in eq . ( [ e : ele010 ] ) , with @xmath26 , is then compared to mc numerical results . we consider a square primitive cell with side of length @xmath27 and periodic boundary conditions ( fig . [ fig : primitive - cell ] , right panel ) . the cnt transverse section is a centered circle in the cell . each cell is labeled with the indices @xmath28 , which define its position in the @xmath5 lattice . a particle moves in the lattice and propagates through multiple reflections on obstacles . the starting position of the particle , projected onto the transverse plane , is determined by @xmath29 , the initial @xmath30 coordinates within the cell and the initial @xmath31 direction ; these quantities are randomly extracted from uniform distributions . we also assume that the particle trajectory starts outside the nanotube .    as a first step , we check if the ion trajectory intersects any obstacle in the cell : 1 ) if it does , we calculate the reflection angle and determine the new ion trajectory ; 2 ) if it does not , we evaluate the ion position at the boundary of the cell , according to its direction . then , we exploit the periodic boundary conditions to place the particle on the opposite side of the cell , with proper position and transverse velocity ( fig . [ fig : primitive - cell ] , right panel ) . this procedure is repeated until a stop condition is fulfilled , that is either when the ion exits the lattice or it covers a length @xmath32 in the transverse plane . below , we show the mc results for the mean exit time , @xmath20 , as a function of @xmath25 from a lattice with @xmath33 nm and @xmath34 nm . when the exit time @xmath20 from the transverse lattice is less than the time it would take to the ion to exit from the top surface of the cnt array , we are losing a signal count in our ideal dm detector . we evaluate @xmath20 for @xmath35 and we interpolate the results with @xmath36 , where @xmath37 is a free fit parameter . the same is done for a @xmath38 nm and @xmath34 nm lattice . the results are shown in fig . [ fig : n - vari - fit ] .    ) , as functions of @xmath25 , for @xmath39 nm , @xmath40 nm ( circles ) and @xmath41   nm , @xmath42   nm ( squares ) lattices . the point interpolation returns @xmath43 ( continuos line ) and @xmath44 ( dashed line ) , respectively . the transverse velocity of the carbon ion is @xmath45 , where @xmath46 ev is its transverse energy and @xmath47 its mass.,width=207 ]    finally , in fig . [ fig : differentdeltafit ] , we show the ratio @xmath48 of the semi - analytical and numerical results for the average exit time from the transverse lattice , extracted for different values of @xmath49 . as shown in fig . [ fig : differentdeltafit ] , the machta - zwanzig approximation becomes effective once @xmath50 .     between semi - analytical and numerical results for the mean exit time from the transverse lattice , for different values of @xmath49.,width=226 ] in this section , we address the question whether a particle reaches the top end before exiting the lattice through its sides . we will consider arrays with larger lattice spacings , @xmath51 nm , and @xmath52 nm . these are the experimental values obtained with a scattering electron microscope ( sem ) image analysis of a sample of aligned cnts ( fig . [ fig : cnt - photo ] ) .     with respect to the electron beam.,width=226 ] we then calculate the lateral losses through the sides in the case of a regular lattice and in the case of a random system of cnts . in this second case , which is probably a better approximation for a real cnt array , the spacing between tubes is not constant . specifically , we assume a gaussian distribution for the displacements of the cnt center positions around the sites of a regular lattice . the simulated carbon ion has an initial energy @xmath53 kev and @xmath54 ev and it can either exit the cnt lattice from the top , after covering the cnt length @xmath55 m , or from the lateral walls , i.e. exit from the transverse lattice , before the top is reached . in the latter case , a potential signal dm counting is lost . we call @xmath56 the fraction of particles leaving from the external perimeter before reaching the top . we expect a @xmath57 behavior , which can be explained as follows . if @xmath25 is sufficiently large , we expect that most of the particles , which exit from the sides , originate within a thin frame close to the external perimeter . the number of particles generated with random initial positions and momenta on the transverse lattice scales with @xmath12 , whereas the number of particles generated close to the perimeter scales with @xmath25 . hence , we expect the @xmath58 scaling .     on a transverse lattice with parameters @xmath52 nm and @xmath51 nm as a function of the lattice size @xmath25 . @xmath56 stands for the fraction of particles leaving from the external lattice perimeter , before reaching the cnt top . we provide results for a regular ( squares ) and a random gaussian lattice ( circles ) , with zero mean @xmath59 and @xmath60 nm . the curves , @xmath61 and @xmath62 , which fit the regular and random lattice values of @xmath56 , are shown as continuous and dashed lines , respectively.,width=226 ]    the results of our mc simulations for the regular and random lattices are shown in fig . [ fig : exit - distance - plot ] with two curves superimposed , @xmath63 and @xmath64 , which fit the data in the regular and random cases , respectively . it is interesting to observe that random lattices trap more particles than regular ones because infinite horizons are closed due to the random displacement of the obstacles .    in fig . [ fig : etavsh - plot ] , we show another plot for @xmath56 as a function of the nanotube height , @xmath6 , with @xmath65 . it is worthwhile noting that lateral losses , if @xmath25 is sufficiently large , do not change dramatically as @xmath6 grows . thus , relying on these results for @xmath56 as a function of @xmath25 and @xmath6 , we can conclude that , if we consider a cnt array with @xmath25 of the order of @xmath66 or larger and a realistic value of the height , @xmath67 300 @xmath59 m , the amount of lateral losses can be kept negligibly small .     for a lattice with parameters @xmath52 nm and @xmath51 nm as a function of the cnt height @xmath6 . @xmath56 stands for the fraction of particles leaving from the external lattice perimeter , before reaching the cnt top . we provide results for a regular ( squares ) and a random gaussian lattice ( circles ) , with null mean @xmath59 and @xmath68 nm.,width=226 ]    these are explanatory results , taken in the limit of small ion longitudinal energy . an estimation of the lateral losses for ions whose initial conditions are produced by the scattering of a galactic halo wimp on the target nucleus , as computed in @xcite , are discussed in sec . [ interactions ] ( see table [ tab : final - results ] ) . in the previous simulations we considered the special case of ions whose transverse energy is below the potential barrier of the cnts so that reflections from nanotubes occur in the transverse lattice . in the following we will consider the general case of ion initial conditions being produced by the scattering of a galactic halo wimp on the target nucleus , as computed in @xcite . in the case the channeling initial conditions are not met , ions do not get reflected at nanotube s surfaces and are also subject to energy losses @xcite and angular deviations @xcite , due to the interaction with atoms when crossing nanotubes . thus , unlike ref . @xcite , the ion motion is not limited to the internal volume of a carbon nanotube , but ions can wander within the entire cnt array with the possibility of traversing cnt surfaces .    in order to consider a more realistic situation , in addition to _ chiral _ cnts , we also include a @xmath69 fraction of _ armchair _ tubes , which is approximately the fraction of metallic nanotubes in a generic cnt sample @xcite . the potential barrier of chiral cnts is axially symmetric and equal to 390 ev , while the potential of armchair cnts can be approximated as a square pulse shape ( * ? ? ? * fig . 5 ) . for the sake of simplicity , the latter is written as @xmath70 here , @xmath71 and @xmath72 are cylindrical coordinates , with @xmath73 , @xmath74 is a random phase , extracted from a uniform distribution in the @xmath75 $ ] range , and 304/2 is the number of  windows " in the tube wall , namely those portions of the cnt surface for which @xmath76 ev . nuclear stopping prevails on the electronic one @xcite for very slow ions ( with a kinetic energy of a few kev per nucleon ) . the latter can thus be neglected , as a first approximation . following ref . @xcite , with some appropriate modifications , we will discuss the formalism to calculate energy losses . the position dependent nuclear stopping force is given by @xmath77 where @xmath72 is the radial coordinate in the plane transverse to the cnt axis and @xmath78 the nuclear distribution , averaged over a cylindrical surface of radius @xmath72 coaxial with the cnt . it can be approximated as @xmath79 where @xmath80 nm is the transverse thermal vibration amplitude at 300 k and @xmath81 nm@xmath82 the atomic density on the cnt surface @xcite . @xmath83 depends on the interatomic potential chosen to model the ion - nucleus interaction .    . @xmath84 is the impact parameter . here , we neglect possible angular deviation effects.,width=132 ]    our choice is the zbl ( ziegler , biersack and littmark ) universal repulsive potential @xcite , because it provides , in the energy range we consider here , results that are more precise than those of classical atomic models , like thomas - fermi or bohr @xcite . in ref . @xcite , the zbl potential was interpolated with the formula @xmath85 where @xmath86 . here , the indices 1 and 2 refer to the projectile and target nuclei , respectively , and @xmath87 is the reduced energy , defined as @xmath88 where @xmath89 is the interatomic screening radius in lindhard s theory ( @xmath90 is the bohr radius ) , @xmath91 the ion kinetic energy , @xmath92 and @xmath93 ( @xmath94 ) the masses and charges of the projectile and target nuclei . we call @xmath84 the impact parameter with respect to the center of the nanotube ( fig . [ fig : dedkovbroken ] ) . when the ion crosses the cnt wall , it loses an amount of energy @xmath95 where @xmath96 . it is then convenient to average @xmath97 over the impact parameter distribution of an homogeneous beam , @xmath98 . in this way , we obtain a parametrization for the average energy loss in a cnt crossing depending on ion kinetic energy only .      when the ion crosses a nanotube surface , it is also subject to angular deviations due to coulomb scattering with nuclei . we consider two possible regimes , depending on the angle of the incident particle with respect to the cnt surface : 1 . single scattering 2 . multiple scattering . a rough way to distinguish between them is discussed in the following .    . the blue line is the limit ion trajectory , which distinguishes between single and multiple scattering regimes . the orange circles are the screened coulomb fields ( of radius @xmath99 ) of two carbon nuclei in the graphene lattice , belonging to the same unit cell . @xmath100 nm is the carbon - carbon bond length in graphene.,width=188 ]    we introduce here a limit angle @xmath101 . if the incident angle between the ion direction and the cnt surface is larger than @xmath101 , the ion is subject to single ( rutherford ) scattering , otherwise undergoes multiple scattering . the incident angle is equal to @xmath101 when the ion trajectory , as it crosses the nanotube surface , is tangent to the screened coulomb fields of two carbon atoms belonging to the same unit cell . the coulomb fields are described as spheres , with radius @xmath102 . @xmath101 can then be calculated as ( see fig . [ fig : thetams ] ) @xmath103 where @xmath100 nm is the carbon - carbon bond length in graphene . if the particle enters the cnt with an angle larger than @xmath101 , we can calculate the ion - nucleus scattering angle using rutherford scattering . in particular , the scattering angle , @xmath104 , is extracted from the distribution @xmath105 @xmath106 . the probability of carbon - ion  nucleus scattering can be estimated as the ratio between the fraction of surface of the graphene lattice unit cell `` occupied '' by the screened nuclear coulomb field and the total area of the unit cell . with an effective shielding radius of the order of @xmath107 , the probability turns out to be @xmath108 .      when a scattered carbon ion moves grazing a cnt wall ( a graphene layer ) with a direction almost parallel to its axis , multiple scattering might occur . the most appropriate treatment to estimate angular deviations of low energy ions undergoing @xmath109 scatterings , is described in  @xcite . the starting point is the differential scattering cross - section by lindhard , nielsen and scharff ( lns ) @xcite , which describes the dynamics of a projectile incident upon a target atom on the basis of classical scattering theory and the thomas - fermi model of atoms . in the treatment in @xcite further assumptions are considered : a ) atoms are described as spheres with radius @xmath110 ; b ) binary collision events have azimuthal symmetry ; c ) there is negligible energy loss in any single collision ; d ) the scattering angles are small . in this way a simplified expression of the angular distribution of the projectile after multiple scattering can be obtained in terms of reduced variables  @xcite . these are the reduced energy , @xmath87 [ see eq . ( [ eqn : reduced - epsilon ] ) ] , the reduced total scattering angle , @xmath111 , and the reduced thickness , @xmath112 , given in dimensionless units . under these assumptions , the angular distribution of the projectile after multiple scattering can be written as @xmath113      \end{array}\ ] ] where @xmath114 and @xmath115 are functions , calculated by numerical integration and tabulated in ref . @xcite . we consider a random lattice , with parameters @xmath27 = 58 nm , @xmath7 = 5 nm , @xmath116 , @xmath60 nm , size @xmath117 and nanotube height of 300 @xmath59 m . the open end of the cnt array is directed along the wimp wind direction , which is opposite to the sun s direction in the galactic frame ( leading towards the cygnus constellation ) . the particle trajectory is considered to start outside the nanotube internal volume at some given distance from the bottom of the cnt array basis extracted from a uniform random distribution . we define different event categories . 1 ) _ top events _ , when the ions exit from the top end of the array , with energy above a certain detection threshold ( 1 kev ) : these are the ions that can eventually be detected ; 2 ) _ bottom events _ , when the ions reach the cnt substrate ; 3 )  _ dead events _ , when the ions reach energies below the detection threshold ( because of energy loss ) while they are still traveling inside the array ; 4 ) _ side events _ , when the ions exit from the lateral sides of the array . for all of them , the initial kinematic conditions are taken from distributions expected for wimp - nuclei collisions . we define @xmath118 to be the angle between the ion trajectory and the sun s direction . ccccc top & bottom & side & dead & total + 30802 & 1786 & 447 & 66965 & 100000 +      a mc simulation including all the effects described above is performed to study the cnt array efficiency . the latter is defined as the ratio between top events and the total number of scattering events . the results are shown in figs . [ fig : initial - distr - e][fig : final - distr - e - cos ] and table [ tab : final - results ] . the estimated efficiency of the cnt array is 31% .    . we also show the distributions of top , dead and bottom events as a function of the initial conditions.,width=207 ]        when the initial ion energy distribution is peaked at low values of @xmath91 , as in fig . [ fig : initial - distr - e ] , most of the particles do not succeed to exit the array , as their energy sets below the detection threshold before they can leave the lattice . on the other hand , the detection rate heavily depends on the orientation of the cnt , as expected ; in particular , since the recoil distribution is peaked for initial angles @xmath119 , fig . [ fig : initial - distr - cos ] , the number of top events is much larger than that of bottom events . this is an important confirmation of the hypothesis of an anisotropic response of the cnt array to the wimp wind described in @xcite . moreover , almost @xmath19 of top events leave the array because of _ secondary _ channeling effects , which can change the initial particle direction bringing them to satisfy the channeling conditions in the cnt interior .        in fig . [ fig : final - distr - e - cos ] , we show the final ion kinetic energies and directions as the particles reach the top or the bottom of the cnt array . the ions might roughly lose on average 4 kev during their path towards the top of the array . as expected , top events take place when ions have an initial kinetic energy higher than the mean initial one ; moreover , we observe that their mean final angle is closer to @xmath120 than the mean initial one , due to secondary channeling effects .      after providing an estimation of the detector efficiency for a given set of lattice parameters , we want to see how it varies with the lattice spacing at fixed cnt radius , @xmath52 nm . this can give important indications to optimize the cost and size of the detector . for example , if distances between tubes are small , we can arrange more cnts on a given surface .    for small values of @xmath27 , particles undergo many collisions with cnts . therefore , their mean energy loss is larger and the number of dead events increases with respect to the large corridor case . however , a relatively strong decrease of the lattice spacing does not heavily reduce the number of top events . specifically , the efficiencies of a lattice with @xmath121 nm and one with @xmath122 nm are @xmath123 and @xmath124 , respectively .    cc @xmath27 [ nm ] & cnt detector mass [ kg ] + 11 & 11.8 + 30 & 1.6 + 45 & 0.7 + 58 & 0.4 +    given this , we want to find a value for the lattice spacing , which may provide a suitable detector mass and a good efficiency . we assume that the total surface of the detector is 100 m@xmath1 . the detector is made up of 100 panels . each of them has a surface of 1 m@xmath1 and consists of a double  brush " , composed of a large number of adjacent cnt arrays , located in a low pressure gas chamber . including the active detection volume of the chamber , each panel might have a thickness of @xmath0 mm . the detector volume is then roughly @xmath125 m@xmath126 . we assume that cnt array units are closely packed with each other , so that the number of cnts in a panel is @xmath127 . considering a cnt surface density of @xmath128 g@xmath129m@xmath82 and a nanotube height of 300 @xmath59 m , we compute the total cnt mass of the detector for different values of @xmath27 . the results are reported in table [ tab : detector - mass ] . in order to have about 1 kg of cnt mass , we should consider arrays with distances between nanotubes less than 30 nm . these values for @xmath27 correspond to an array efficiency of @xmath123 . we studied the dynamics of low energy ions in a carbon nanotube lattice with monte carlo simulations . our results were used to provide an evaluation of the efficiency of a wimp directional detector , composed of a large number of adjacent cnt arrays . collisions with wimps might struck carbon ions out of their lattice positions and inject them in the array with initial kinematical conditions depending on the relative orientation of the cnt axes with respect to the apparent wimp wind from the cygnus . the smaller the angle between these two directions , the higher the probability of channeling within the single nanotubes or from the interstices among them . the ions whose trajectories exit from the surface with open end nanotubes ( on the opposite side cnts terminations are closed by a substrate as in fig . [ fig : cnt - photo ] ) , are eventually detected and constitute the directional dm signal . the array efficiency depends on several technical specifications , as the length of the cnt array basis in terms of number of nanotubes and the lattice spacing . ions with right initial conditions , which anyway get lost from the lateral sides of the array , are weakening the signal . it is thus worthwhile to estimate the lateral ion losses as a function of the size of the array . we showed that lateral ion losses can be kept sufficiently low if we consider cnt arrays with a number of nanotubes @xmath130 , within a large range of nanotube heights . as can be seen from the results of table  [ tab : final - results ] , calculated with the initial kinematic conditions from the ion recoil distribution ( see @xcite and figs . [ fig : initial - distr - e ] , [ fig : initial - distr - cos ] ) , the expected lateral losses are of the order of @xmath131 , approximately .    therefore , commercially available cnt arrays , with a lateral side of @xmath132 nanotubes , are substantially unaffected by this problem . the detector efficiency should also depend on the lattice spacing , @xmath27 , between nanotubes  in a weak way as we found . because of this , we can consider a small lattice spacing , @xmath133 nm , in order to pack a suitable detector mass , of the order of 1 kg , on a given surface ( 100 m@xmath1 , approximately ) preserving , at the same time , the detection efficiency .    an important outcome of the present study , which meets the expectations of ref . @xcite , is that a cnt array cooperates to effectively provide the required anisotropic response for the directional detection of wimps . in fact , the number of top events estimated is much larger than that of bottom events ; the latter corresponds to the number of top events when the cnt brush is reversed with respect to the direction of the wimp wind . specifically , the ratio between bottom and top events is @xmath134 . therefore , the double cnt brush detector in @xcite could effectively discriminate the wimp signal from the background . in @xcite the analysis was focused on the calculation of ion channeling within the fundamental constituent of the lattice : a single nanotube . it was shown that it is difficult to match the right channeling conditions inside a single cnt immediately after a wimp collision . however the cnt array as a whole  including interstices among nanotubes  can enlarge the acceptance of the detector cooperating at driving the trajectories of not perfectly collimated ions towards the open exit . indeed , we showed here that those ions which can not be channeled in a single nanotube , but have initial kinematic conditions close to the channeling ones , have further chances to be channeled by the array , including the interstitial space between cnts .    in ref . @xcite , the calculated fraction of ions channeled in the nanotubes , immediately after a wimp scattering event , is @xmath135 ; see ( * ? ? ? * fig . 7 ) . the remaining large majority , which is the subject of the present study , consists of those particles scattered outside the single nanotube and/or not satisfying the very strict channeling conditions within a single cnt . the amount of these latter particles , which might be anyway guided by the system of nanotubes towards the top end , can be extracted from table [ tab : final - results ] . therefore , by combining the previous results of ref . @xcite and those of the present paper , we estimate that the average fraction of particles which can in principle be detected grows up to @xmath136 after a wimp scattering event .    to quantify the channeling potentiality of our cnt lattice as a whole with respect to the single nanotube case , we can also introduce an acceptance angle , @xmath137 . specifically , the acceptance angle is calculated in such a way that the integral over the ion recoil distribution up to @xmath137 is equal to the fraction of channeled top events , namely those events which reach the top end satisfying channeling conditions . the computed value for the cnt array is @xmath138 , while that for the single nanotube , calculated by using the findings of ref . @xcite , is @xmath139 . thus , an increase in the detector efficiency of the lattice as a whole , with respect to the single nanotube , can be translated in terms of a larger detector acceptance . such an improvement has obvious consequences on the construction specifications of a dark matter detector , especially in terms of target mass . testing experimentally the single carbon nanotube array unit , with electron or neutron probes for example , will finally tell to which extent the approximation used in our calculations are realistic . we thank m. diemoz for supporting this research activity and m.g . betti and c. mariani for giving access to the cnis facility at sapienza univ . we acknowledge former collaboration on these topics with l.m . capparelli and d. mazzilli . g.c . acknowledges partial support from erc ideas consolidator grant crysbeam g.a . n.615089 . following ref . @xcite , we calculate the residence time , @xmath140 , of a particle in an infinite - horizon billiard with an elementary square cell . see fig . [ fig : mz ] ( left panel ) .      for @xmath145 , the particle motion in the elementary square cell can be regarded as ergodic . thus , the probability that a particle leaves the cell in a time smaller than @xmath143 is given by the ratio between the phase space available to exit , which we call @xmath146 , and the total phase space , @xmath147 . the total volume of the phase space of the cell is @xmath148 , with @xmath149 . the calculation of the phase space available to exit the cell is more complex and will be discussed in the following . we refer to fig . [ fig : mz ] . we assume that , due to the smallness of @xmath150 , the obstacle curvature can be ignored . a particle , which is at distance @xmath141 away from the border , leaves the cell in a time smaller than @xmath143 if @xmath151 and @xmath152 . thus , the portion of velocity space for the particle to exit the cell has an area @xmath153 . we can now integrate over the rectangle @xmath154 of fig . [ fig : mz ] ( right panel ) and obtain @xmath155 where we also have a factor of four because there are four possible exit paths . the particle exits the cell when @xmath156 @xmath157 . the average residence time in the cell is thus @xmath158    l.  m.  capparelli , g.  cavoto , d.  mazzilli and a.  d.  polosa , phys . dark univ . * 9 - 10 * , 24 ( 2015 ) ; l.  m.  capparelli , g.  cavoto , d.  mazzilli and a.  d.  polosa , corrigendum to phys . dark univ . * 9 - 10 * , 24 ( 2015 ) . g. j. alner _ et al . _ , instrum . meth . a * 555 * , 173 ( 2005 ) ; s. burgos _ et al . _ , * 31 * , 261 ( 2009 ) ; j.  b.  r.  battat _ et al . _ [ drift collaboration ] , phys . dark univ . * 9 - 10 * , 1 ( 2014 ) .
Abstract:

Answer: large arrays of aligned carbon nanotubes ( cnts ) , open at one end , could be used as target material for the directional detection of weakly interacting dark matter particles ( wimps ) . as a result of a wimp elastic scattering on a cnt , a carbon ion might be injected in the body of the array and propagate through multiple collisions within the lattice . the ion may eventually emerge from the surface with open end cnts , provided that its longitudinal momentum is large enough to compensate energy losses and its transverse momentum approaches the channeling conditions in a single cnt . therefore , the angle formed between the wimp wind apparent orientation and the direction of parallel carbon nanotube axes must be properly chosen . we focus on very low ion recoil kinetic energies , related to low mass wimps ( @xmath0  gev ) where most of the existing experiments have low sensitivity . relying on some exact results on two - dimensional lattices of circular obstacles , we study the low energy ion motion in the transverse plane with respect to cnt directions . new constraints are obtained on how to devise the cnt arrays to maximize the detection efficiency .


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: galaxy clusters play an important role in precision cosmology that complements other techniques like type ia supernovae luminosity - distance relation measurements @xcite , baryonic acoustic oscillations ( baos ) angular - distance relation measurements @xcite , and observations of the cosmic microwave background radiation ( see @xcite for a review ) . for example , @xcite exploit dark energy s influence on the growth of structure by using 37 moderate - redshift and 49 low - redshift clusters to measure the shape of the galaxy - cluster mass function and its redshift evolution , which constrain the dark - energy density parameter @xmath1 to @xmath2 in a non - flat @xmath3 cosmology and the dark energy equation of state parameter @xmath4 to @xmath5 in a flat cosmology . @xcite have obtained similar results from measurements of the evolving number density of the largest clusters in order to constrain @xmath4 to @xmath6 . strategies like these that compare model predictions to galaxy - cluster sample statistics inevitably confront sample error and sample bias . future surveys expected to gather samples of 1040 thousand galaxy clusters will maximize survey reach while maintaining sufficient observation quality in order to minimize sample error , but they still must grapple with a major source of sample bias , which is scatter in the relationship used to infer cluster mass from an observable mass proxy . an important mass - observable relation for galaxy cluster studies connects dark matter halo mass to the temperature of the intracluster medium ( icm ) inferred from its x - ray spectrum ( its `` x - ray temperature '' @xmath7 ) . in this paper we investigate the possibility of correcting for scatter in this relation using temperature - inhomogeneity in the icm , and discuss challenges that may exist in such a program . a significant amount of uncertainty in the dark - energy constraints obtainable from large cluster surveys derives from uncertainty in scatter about the mean scaling relations obeyed by galaxy clusters bulk properties @xcite . the key galaxy cluster property to measure when trying to constrain dark energy with clusters is the cluster s mass , which can not be directly observed . theoretical considerations predict correlations among halo mass and more readily observed cluster properties , like its galaxy richness , the velocity dispersion of its galaxies , @xmath7 , the sunyaev - zeldovich decrement , the gas mass , and @xmath8 parameter , which is the product of @xmath7 and the gas mass inferred from x - ray observations @xcite . theory also predicts intrinsic scatter in these relations owing to variation in cluster dynamical state ( see , for example * ? ? ? * ; * ? ? ? * ; * ? ? ? one way to deal with intrinsic scatter is to join a cluster model to a cosmological model and simultaneously fit for the parameters of both , leveraging the statistical power of large surveys and `` self - calibrating '' the mass - observable relations . another approach is to combine observables that tend to depart from the expected scaling relations in opposite ways , yielding a new , low - scatter composite observable . an example low - scatter composite observable is @xmath8 , since at a given halo mass , offsets in the measured gas mass at fixed total mass tend to anti - correlate with offsets in the measured temperature @xcite . another family of low - scatter composite observables attempt to measure structural variation directly . mergers , relaxation , and non - adiabatic processes like radiative cooling , star formation , and feedback ought to leave a visible imprint that may allow us to measure and correct for scatter . for example , one might use imaging to quantify resolved morphological substructure . @xcite apply two observationally - motivated structure measures , the power ratios @xcite and the centroid shift @xcite , to a sample of galaxy clusters simulated with  @xcite . they find that cluster structure correlates strongly with bias in mass estimates derived from @xmath7 under the assumption of hydrostatic equilibrium and accounting for cluster structure can be used to correct some of the bias . similarly , @xcite apply the power ratios , centroid shift , and axial ratio @xcite substructure measures to the same sample of simulated clusters used in this paper , and find that cluster substructure correlates with departures from the mean  relationship in the sense that clusters with more substructure tend to have a lower temperature at a given mass , and can be used to refine mass estimates derived from the icm x - ray temperature . @xcite likewise find that greater substructure , as quantified with power ratios , correlates with lower temperature at a given mass . @xcite find a strong correlation between mass - temperature scatter and halo concentration in their sample of simulated clusters , with cooler clusters appearing more concentrated than warmer clusters at similar mass . @xcite observe , however , that line - of - sight projection effects lead to significant uncertainties in morphologically - derived substructure measures . substructure also becomes more difficult to resolve at high redshift . spectral signatures of dynamical state are therefore attractive because they are aspect - independent and redshift - independent . one such spectral signature of dynamical state is the `` temperature ratio '' @xcite , which divides a `` hardband '' spectral - fit temperature by a `` broadband '' spectral - fit temperature . an energy cut applied to a broadband spectrum produces a hardband spectrum and serves to filter out cooler line - emitting components of the icm . @xcite studied the effects of relaxation on the observable properties of galaxy clusters , using a sample of numerically - simulated galaxy clusters generated by @xcite . they found that hardband ( 2.09.0 kev ) x - ray spectral fit temperatures average @xmath920% higher than broadband ( 0.59.0 kev ) temperatures and suggested that this effect may signal the presence of cool , luminous sub - clusters lowering the broadband temperature . @xcite corroborated these findings in simulations that included radiative cooling . figures [ fig : f01 ] and [ fig : f02 ] illustrate the effect . both show a counts spectrum and single - temperature fit for a typical , unrelaxed cluster in our sample , with an aperture set at @xmath10 and a core region out to @xmath11 excised . figure [ fig : f01 ] is for a single - temperature model fit to the broad band , while figure [ fig : f02 ] is for the hard band . note the excess emission relative to the model above 4.0 kev and the deficit below 2.5 kev that arises because the model can not simultaneously fit both a hot component and a cooler , line - emitting component @xcite . figure [ fig : f02 ] shows the same spectrum , but with the residuals for a model fit just to the hard band . in this case the fit is much better over the range from 2.0 to 7.0 kev but under - predicts the emission below 2.0  kev . @xcite suggested that the temperature skewing they observed might indicate a real temperature skewing detectable in real clusters using . @xcite fit single - temperature emission models to the hard band ( 2.07.0 kev ) and broad band ( 0.77.0 kev ) for a large ( n = 192 ) sample of clusters with observations selected from the  data archive . these authors show that for a large , heterogeneous sample of clusters across a broad temperature range , the distribution of  has a mean of 1.16 and an rms deviation @xmath12 , with  tending to be larger in merging systems . they also report that while this signal is significant in the aggregate , the errors for any single measurement and the scatter across all of the measurements together pose challenges for any effort to use  either to select for merging systems or to obtain more accurate mass estimates . in order to meet these challenges it is important also to examine  and similar spectral signatures of dynamical state in a simulation context . like the original study of @xcite , this paper examines the temperature ratio for a sample of simulated clusters and simulated  observations , and as in the subsequent study by @xcite , the simulated clusters analyzed here were generated using a hydrodynamical code with radiative cooling included . this paper is organized as follows . in @xmath13 , we describe the simulated clusters in our sample along with the 2 code for generating their mock x - ray observations . in @xmath14 , we present our analysis methods and define the spectral signatures of dynamical state we examine . in @xmath15 , we report and discuss our results , and @xmath16 summarizes our work . this study is based on an analysis of 118 clusters simulated using the cosmological hydrodynamics tree+sph code gadget-2 @xcite , which were simulated in a standard @xmath17 cold dark matter ( @xmath17cdm ) universe with matter density @xmath18 = 0.3 , @xmath19 = 0.7 , @xmath20 = 0.04 , and @xmath21 = 0.8 . the simulation includes radiative cooling assuming an optically - thin gas of primordial composition , with a time - dependent uv background from a population of quasars , and handles star formation and supernova feedback using a two - phase fluid model with cold star - forming clouds embedded in a hot medium . all but four of the clusters are from the simulation described in @xcite , who simulated a box @xmath22 on a side , with @xmath23 dark matter particles and an equal number of gas particles . the present analysis considers the 114 most massive clusters within this box at @xmath24 , which all have @xmath25 greater than @xmath26 . these are referred to as the b04 sample in the remainder of this paper . by convention , @xmath27 refers to the mass contained in a sphere which has a mean density of @xmath28 times the critical density @xmath29 , and whose radius is denoted by @xmath30    that cluster set covers the @xmath91.5 - 5 kev temperature range , but the @xmath22 box is too small to contain significantly hotter clusters . we therefore supplemented it with four clusters with masses @xmath31 and temperatures @xmath32 kev drawn from a dark - matter - only simulation in a larger @xmath33 box @xcite , referred to in this paper as the d09 sample . the cosmology for this simulation also was @xmath17cdm , but with @xmath21 = 0.9 . these were then re - simulated including hydrodynamics , radiative cooling , and star formation , again with gadget-2 and using the zoomed - initial - conditions technique of @xcite , with a fourfold increase in resolution . this is comparable to the resolution of the clusters in the smaller box . adding these four massive clusters to our sample gives a total of 118 clusters with @xmath25 in the range @xmath34 to @xmath35 . the mean structural properties of massive clusters drawn from a sample with @xmath21 = 0.9 may differ somewhat from those of similar - mass clusters in a @xmath21 = 0.8 universe because they reflect a more advanced state of cosmic evolution . however , these four additional clusters carry minimal statistical weight in the context of the overall sample . they are included primarily to evaluate whether the mean and dispersion of their @xmath36 values are consistent with those of the lower - mass systems . the simulated galaxy clusters in our sample are processed with the _ x - ray map simulator _ version 2 ( 2 ) @xcite to generate x - ray images suitable for standard  reduction techniques . in its first step 2 uses the outputs of the hydrodynamic code to calculate the emissivity of each simulation element within the chosen field of view and to project this onto the image plane . in its second step it convolves the resulting flux with the appropriate response of a given detector . in the case of our simulated observations , the second step applies the response matrix file and ancillary response file for the acis s3 ccd , with a 200 kilosecond exposure time . in order to separate out potential calibration issues from the focus of this particular study , these response matrices implement a constant response over the detector , and in generating the simulated x - ray images , the evolved zero - redshift clusters are shifted to a redshift sufficient to fit @xmath37 within the 16 arcminute field of view . because of this step , though the clusters in the sample are distributed over a range of physical sizes , they all have approximately the same apparent size projected onto the image plane and into the simulated observations .      in this paper , cluster mass refers to @xmath25 , while @xmath7 refers to our estimated `` spectral - fit temperatures '' obtained using @xcite to fit single - temperature mekal plasma models to various energy bands in simulated 2 x - ray observations taken from the numerically - simulated clusters . the steps involved in the process are described in detail in @xmath38 [ sec : analysis ] . figure [ fig : f03 ] shows the mass - temperature relation based on our sample of simulated clusters and estimated spectral fit temperatures , measured in a broad 0.77 kev band within an annular aperture from @xmath39 to @xmath10 as in @xcite . the best fits to the power - law form @xmath40 have the coefficients @xmath41 , @xmath42 for the full sample of clusters . for the subset of clusters with @xmath43kev , the best - fit parameters are @xmath44 , @xmath45 . as is generally the case for simulated clusters , the power - law indices of the mass - temperature relations found here are consistent with cluster self - similarity and the virial theorem @xcite . these relationships have scatter , which we characterize by the standard deviation in log space @xmath46 about the best - fit mass at fixed temperature @xmath47 . we find @xmath48 . 2 kev are plotted in purple open triangles , while clusters with @xmath49 2 kev are plotted in blue open squares , and the four massive d09 clusters in green open circles . mean relations are plotted with solid lines for the whole sample , and dashed lines for clusters with @xmath49 2kev . finally , note that average temperatures are taken within a core - excised annulus whose outer radius is placed at @xmath10 . , scaledwidth=80.0% ] we use 2 to simulate x - ray observations  provided as standard  event files  then subject them to a series of reduction steps using the  interactive analysis of observations package ( ciao ) v4.1 @xcite . all of the event files have at least 100k counts , while those for the most luminous clusters have over half a million counts . simulated observations of this quality provide us with an opportunity to address the question of how much of the scatter in our temperature - structure statistics is intrinsic . consequently , our first reduction step filters each raw event file into a new set of event files by randomly sampling the events it records . the number of counts in each file is its `` count level , '' and in order to cover the space of typical  archival observations , we apply the filtering step to each of the raw event files four times , at count levels 15k , 30k , 60k , and 120k . these count levels roughly map to the range between observations of short duration or of relatively low surface - brightness objects to observations of long duration or of relatively high surface - brightness objects . an example of a filtered event file appears in figure [ fig : f04 ] as a surface - brightness image . it displays a common feature of numerical simulations of galaxy clusters , which is the presence of relatively dense , cool , and metal - rich substructures that continuously undergo mass accretion and have not yet come into thermal equilibrium with the hot icm surrounding them @xcite . these bright point - like spots , or `` cool lumps '' , typically are associated with the dense cool cores of smaller halos that have merged with the primary halo . generally regarded as an unphysical artifact of numerical simulations , at least insofar as their temperatures , densities , and concentrations are concerned , commonly they are excised before further analysis is conducted @xcite .     and [ fig : f02 ] . the aperture is set at @xmath10 and a core region of size @xmath50 is excised . note the surface - brightness peaks in the image . in this study these are detected automatically by the ciao tool wavdetect , which generates source regions identified in this figure by red ellipses . labeled `` cool lumps '' in this analysis , they are subjected to varying degrees of excision , from no masking ( _ left _ ) to full masking ( _ right _ ) , as described in the text . , scaledwidth=100.0% ]    in order to study the effect that excising cool lumps has on measures of temperature substructure , we produce from the originals several new event files whose cool lumps have been excised to an increasing degree . the ciao ` wavdetect ` tool identifies peaks in the photon distribution by correlating an event file s image with a sequence of `` mexican - hat '' wavelet functions of differing scale sizes , measured in pixels , then generates a source list with associated region files . we use its default sequence of wavelet scale sizes in this analysis ( 2 and 4 pixels ) , though we apply the tool four times to each filtered event file , each time incrementing the multiplicative factor by which the source regions are scaled . in ciao ` wavdetect ` , the parameter governing this multiplicative factor is ` ellsigma ` and in our study ranges between 0 to 3 . a value of 0 for ` ellsigma ` is equivalent to `` no masking '' while a value of 3 corresponds to what we call `` full masking '' . the effect is to produce versions of each cluster observation with a range of masking . note that the wavelet scale size is a constant fraction of the cluster size because all the clusters have been redshifted so that @xmath37 fits in the field of view . we finish the extraction phase of our analysis with the following steps . first , for every excised file we apply an aperture of @xmath10 and excise the central @xmath51 region , generating new copies of the event files . this step , including the centering algorithm , matches the procedure in @xcite so that regions measured in the simulated clusters correspond to those in the  archival observations . next , each fully - processed event file file is extracted into a standard pulse - invariant ( pi ) spectral file binned so as to have at least 25 counts in each energy channel . the extraction phase is complete when each original raw event file generates pi spectral files suitable for spectral fitting in . the pi spectral files are then fed into  v12.5.0 for spectral fitting . in order to form the temperature ratio  two fits are performed for each spectrum . the first fit is over the hard band from 2.0@xmath527 kev , for which the @xmath52 factor exists in order to shift the 2 kev cutoff from the observer s frame to the cluster s rest - frame . the simulated clusters occupy a range of redshifts because larger clusters are translated to higher redshifts in order to fit @xmath37 within 16 arc - minutes , when creating artificial observations . the second fit is over the broad band , including all energy channels in the spectrum from 0.77 kev . every fit is made by minimizing the @xmath53 statistic for a single - temperature mekal model multiplied by a warm - absorber ( wabs , to account for galactic absorption ) . the galactic column is fixed to @xmath54 and the metallicity to 0.3 solar , leaving the temperature of the emission component , its h density ( although this makes no difference ) , and its normalization as the only free parameters . we adopt two spectral measurements of temperature structure . the first is the temperature ratio  of @xcite and is found in the following way . for each simulated cluster in our sample , for each count level ( 15k , 30k , 60k , 120k ) , and for each cool lump masking level ( 0 , 1 , 2 , 3 ) , we find a spectral - fit temperature in the hard band and the broad band and form the temperature ratio . emission from cooler , line - emitting metal - rich parts of the icm ought to be excluded from the hard band , so we expect that a  value greater than unity signals the presence of merging sub - clumps . our second measure of temperature , which we call the `` cool residual , '' ( @xmath55 ) compares an observation s actual broadband count rate to a model - predicted count rate for a spectral model fit only over the hard band . again , we determined a spectral - fit temperature for each combination of cluster , count level , and masking level , except that in this case we performed only a hardband fit . since emission from cooler components should be excluded from this band , in general we achieve good fits even when the count level is sufficiently high that broadband fits may formally have large reduced @xmath53 values . from this hardband model fit we estimate the corresponding broadband count rate and calculate its percentage deviation from the actual broadband count rate . while single - temperature systems will have actual count rates that are essentially the same as their model count rates , the introduction of cooler , luminous substructure components should generate an excess broadband count rate relative to the model . this measure of temperature substructure avoids the uncertainties associated with a single - temperature fit to the broad band . having calculated spectral - fit temperatures , , and the cool residual for all of the clusters in our sample , we have the means to study how well temperature structure correlates with departures from the mean  relation , and how well these measures can be used to obtain better mass estimates . we will also examine how our results depend on the removal of the cool lumps . we present in figure [ fig : f05 ] the temperature ratios  for our simulation sample , plotted as a function of the broadband temperature fit . this figure is similar to figure 8 in @xcite , which shows that the mean value of  observed among clusters in the _ chandra _ archive is @xmath56 , with a standard deviation of @xmath57 . we remind the reader of important differences between the two samples being considered . the archive sample of @xcite contains clusters most of which have @xmath58 3 kev , whereas most of the clusters in our sample have @xmath59 3 kev . nevertheless , even with that caveat the temperature ratios of the simulated clusters are distributed in approximately the same way as are those in the real sample of @xcite , with @xmath60 and @xmath61 , provided the cool lumps are not excised ( circles in figure [ fig : f05 ] ) . when the cool lumps are excised ( squares in figure [ fig : f05 ] ) , the mean and variance of  diminish , with @xmath62 and @xmath63 . these values are for the full 120k count data , though the full set of @xmath64@xmath65 for the four ` ellsigma ` values and four count values are tabulated in table [ tbl:01 ] . this table provides mean values for  and its variance for all of the clusters in the sample , and for a subset whose  value is greater than 2 kev . much of the variance derives from the lower - temperature clusters , and as they are removed the variance in  drops significantly , especially when full masking is applied . figure [ fig : f06 ] presents a similar effect for @xmath55 . here , we see that the excision of cool lumps again reduces the mean and variance of the temperature substructure measure , although in this case the effect is more dramatic . evidently , the presumably unphysical cool lumps in simulated clusters may be necessary to reproduce the distribution of quantitative temperature substructure measures found in real clusters . this is a subject which we return to in section [ sec : masking ] . our original motivation for conducting this study was to determine if temperature structure , as quantified by the temperature ratio , correlates with and can be used to correct for departures from the mean mass - temperature relation . in order to test this idea , we define the `` mass offset '' at fixed temperature by the relation @xmath66 \label{eq : deltam}\ ] ] where @xmath67 is the cluster s actual mass , and @xmath68 is the mass predicted from the mean  relation . figure [ fig : f07 ] plots the mass offsets for our simulated clusters as a function of , for the case in which predicted masses are derived from the relation for the broadband spectral fit temperature . values of calculated both with cool lumps excised ( squares ) and without excision ( circles ) appear in this figure . error bars on  are omitted for clarity . while there is some correlation , such that clusters with more temperature structure ( larger ) tend to be more massive than predicted by the mean  relation , the trend is weak and has substantial scatter . excising the cool lumps weakens the trend further . figure [ fig : f08 ] shows the same comparison for the @xmath55 measure instead of , and in this case we again find that excising the cool lumps has an even larger effect for although in both cases accounting for temperature substructure does not greatly reduce scatter in mass offset .     from the mean  relationship for the clusters in combined b04+d09 simulation sample . the temperatures in this relation are our broadband spectral - fit temperatures . squares correspond to simulated 2 observations that are processed by wavdetect with ellsigma set to 3 ( full masking ) , while circles correspond to observations that have ellsigma set to 0 . , scaledwidth=80.0% ]     and mass offset @xmath69 from the mean  relationship for the clusters in combined b04+d09 simulation sample . the temperatures in this relation are our broadband spectral - fit temperatures . squares correspond to simulated 2 observations that are processed by wavdetect with ellsigma set to 3 , while circles correspond to observations that have ellsigma set to 0 . , scaledwidth=80.0% ]      we now examine the decline in  and its variance as cool lumps are excised , beginning with figure [ fig : f09 ] , which focuses on . here , we plot the statistic s standard deviation @xmath70 as a function of the two dimensions along which we adjust our analysis pipeline , with the top panel devoted to the masking strategy , and the bottom panel devoted to the count level . focusing on the top panel , we see that increasing the ` ellsigma ` parameter of the ciao ` wavdetect ` tool from 0 to 3 reduces @xmath71 from 0.15 down to 0.13 , when the observations are relatively `` poor '' ( with @xmath915k counts ) . with higher - quality observations of 60k or 120k counts , the decline is in @xmath71 greater , going from 0.12 down to 0.07 as the ` ellsigma ` parameter rises from 0 to 3 , and @xmath71 for full masking ends up being significantly less than the value of 0.10 observed in the @xcite sample . as wavdetect s ellsigma parameter ranges from 0 to 3 . these are plotted for a family of simulated observations of increasingly higher quality , from approximately 15k counts to approximately 120k counts . ( _ bottom _ ) decline in @xmath71 as simulated 2 observations go from lowest - quality ( approximately 15k counts ) to highest - quality ( approximately 120k counts ) . these are plotted for a family of simulated observations with wavdetect s ellsigma parameter ranges from 0 to 3 . , scaledwidth=90.0% ]    figure [ fig : f10 ] helps show why the temperature ratio  and its variance decline as the cool lumps are removed . it plots the sample average of the relative change in temperature as the ciao ` wavdetect ` ` ellsigma ` parameter is increased . as more of each cool lump is excised , both the broadband spectral fit temperature , and the hardband fit temperature increase . however , the increase is significantly larger for the broadband temperature , as the cool lumps contribution to the flux is already largely excluded from the hardband fits . the top panel in this figure shows this effect for spectra with approximately 120k counts , while the bottom panel is for spectra with approximately 15k counts . .,scaledwidth=90.0% ]    another way of visualizing the effect of more aggressive masking is depicted in figure [ fig : f11 ] . the top panel in this figure is similar to figure [ fig : f06 ] in that it occupies the   plane , with  plotted as a function of for our simulated clusters . arrows illustrate the shift in  and in , as the cool lumps are excised . some of the clusters experience very large shifts in both quantities , while others experience no shifts at all . the former are associated with clusters that have many well - defined and bright cool lumps , while the latter correspond to those clusters that are completely free of cool lumps . the bottom panel in this figure shows just the size of this shift for the  statistic , from which we can see that some clusters indeed have a shift of precisely 0 . again , these are clusters whose simulated x - ray observations are unchanged after applying the ciao ` wavdetect ` tool because it finds no sources to mask out .    . notice that the distribution of most of the points shifts significantly to higher spectral fit temperatures and smaller temperature ratios as cool lumps are excised . notice also that some simulated clusters are unafflicted by cool lumps , so that their net shift is 0 . ( _ bottom _ ) the  component ( y - axis ) of the shift presented in the top panel . note that these are for simulated 2 observations of maximum quality , having approximately 120k counts . , scaledwidth=90.0% ]      various researchers have established a clear correlation between morphological measures of cluster dynamical state and @xmath72 @xcite , a correlation which is less apparent in our study of spectral measures of substructure . similiarly , @xcite found a correlation between the temperature ratio and structure for their  sample , in that merging events are associated with elevated . to probe this issue further we focus on  and compare it to several metrics for cluster substructure . these are the centroid variation @xmath73 , the axial ratio @xmath74 , and the power ratios @xmath75 and @xmath76 ( see * ? ? ? * and references therein ) . the centroid variation @xmath73 measures the skewness of a cluster s two - dimensional photon distribution by calculating the variance in a series of isophotes for the cluster surface - brightness map the axial ratio @xmath74 measures a cluster s elongation , which tends to increase during merger events . the power ratios @xmath75 and @xmath76 decompose a cluster s surface - brightness image into two - dimensional multipole expansions , capturing different aspects of a cluster s geometry . @xmath75 relates to the ellipticity in an image and is similar to the axial ratio @xmath74 , while @xmath76 measures the `` triangularity '' in an image .    in @xcite we calculated these morphological substructure metrics for a superset of the b04 sample of simulated galaxy clusters and compared them to departures from the  relation ( see fig . 8 in that paper ) . here we use the same substructure measures for the 114 b04 clusters for which we are able to calculate  and compare the results . these are presented in figure [ fig : f12 ] . whereas in @xcite there is a clear correlation between morphological substructure and @xmath72 , we find in this study that there is little or no correlation between substructure and for our b04 sample .    , the centroid variation @xmath73 , and the power ratios @xmath75 and @xmath75 . squares correspond to simulated 2 observations that are processed by wavdetect with ellsigma set to 3 ( full masking ) , while circles correspond to observations that have ellsigma set to 0 . , scaledwidth=90.0% ]    in order to understand how temperature structure and morphological structure can be uncorrelated , we looked at four simulated clusters having approximately the same temperature ( @xmath77 ) and occupying the relative extremes in  and in the centroid variation @xmath73 . two were selected for relatively low ( @xmath78 ) but large variation in @xmath73 ( @xmath79@xmath80 ) . two were selected for large  ( @xmath81 ) but again large variation in the morphological structure parameter @xmath82 ( @xmath83@xmath80 ) . these clusters are presented in figure [ fig : f13 ] with surface - brightness contours overplotted and with the associated values for  and @xmath73 . the two clusters in the left column appear more symmetric in their contours , while the two in the right column exhibit noticeable centroid shift . however , the morphologically apparent structure in the lower right cluster is not observed spectrally in the temperature ratio . evidently , neither morphology nor  is a perfect measure of relaxation in our simulation sample , as it contains clusters with small centroid shift @xmath82 and obvious substructure , as in the upper left of this figure . and , it contains clusters with large @xmath73 that are nearly isothermal , as in the lower right of this figure . we used a sample of galaxy clusters simulated with radiative cooling and supernova feedback , along with simulated  x - ray observations of these clusters , to study temperature inhomogeneity as a signature of cluster dynamical state . specifically , we adopted two methods of quantifying temperature inhomogeneity spectroscopically , the temperature ratio @xcite and the cool residual @xmath55 . the former is the ratio of a hardband x - ray spectral - fit temperature to a broadband temperature , and becomes greater than 1 for clusters whose icm contains cool , over - luminous sub - components . the latter is the excess broadband count rate relative to the count rate predicted by a model fit to the hard x - ray band . though our simulated clusters are typically less massive and have lower temperatures than the  archive clusters in @xcite , we find that their temperature ratios  occupy generally the same distribution as the observed clusters . we also looked for an opportunity to combine  and the cool residual with the mean mass - temperature relation to obtain better mass estimates than are achieved just with the scaling relation alone . we find , however , that while both  and the @xmath55 are correlated with offset from the  relation , these correlations are weak , at least for this sample . we conclude that these measures of temperature inhomogeneity are not very effective at reducing scatter in the mass - temperature relation . finally , we note a particular difficulty that arises when trying to use clusters from hydrodynamic simulations to calibrate scatter - correction observables based on temperature inhomogeneity , such as . historically , simulated clusters have tended to exhibit their own kind of `` over - cooling problem '' , in which dense lumps of luminous , cool gas associated with merged sub - halos appear . these cool lumps are often excised from simulated clusters before their global properties are measured , but masking them appears to make the remaining temperature structure of the simulated clusters overly homogeneous . this finding suggests that real clusters may have cooler sub - components , that are more diffuse and less concentrated than in their simulated counterparts . some physical process , perhaps thermal conduction , turbulent heat transport ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) , or a more aggressive form of feedback , prevents cool lumps from forming in real clusters and might not completely eliminate those temperature inhomogeneities . newer simulations incorporate treatments of conduction and agn feedback , as well as more accurate treatments of mixing , and it will be interesting to revisit these temperature inhomogeneity measures in simulated clusters when large samples of such simulated clusters become available . the authors wish to thank stefano borgani for contributing the simulation data on which this project was based . this work was supported by nasa through grants nng04gi89 g and nng05gd82 g , through chandra theory grant tm8 - 9010x , and through chandra archive grant saoar5 - 6016x . some of it was done at the kavli institute for theoretical physics in santa barbara , supported in part by the national science foundation under grant no . nsf phy05 - 51164 . eisenstein , d.  j. , zehavi , i. , hogg , d.  w. , scoccimarro , r. , blanton , m.  r. , nichol , r.  c. , scranton , r. , seo , h. , tegmark , m. , zheng , z. , anderson , s.  f. , annis , j. , bahcall , n. , brinkmann , j. , burles , s. , castander , f.  j. , connolly , a. , csabai , i. , doi , m. , fukugita , m. , frieman , j.  a. , glazebrook , k. , gunn , j.  e. , hendry , j.  s. , hennessy , g. , ivezi , z. , kent , s. , knapp , g.  r. , lin , h. , loh , y. , lupton , r.  h. , margon , b. , mckay , t.  a. , meiksin , a. , munn , j.  a. , pope , a. , richmond , m.  w. , schlegel , d. , schneider , d.  p. , shimasaku , k. , stoughton , c. , strauss , m.  a. , subbarao , m. , szalay , a.  s. , szapudi , i. , tucker , d.  l. , yanny , b. , & york , d.  g. 2005 , astrophysical journal , 633 , 560        fruscione , a. , mcdowell , j.  c. , allen , g.  e. , brickhouse , n.  s. , burke , d.  j. , davis , j.  e. , durham , n. , elvis , m. , galle , e.  c. , harris , d.  e. , huenemoerder , d.  p. , houck , j.  c. , ishibashi , b. , karovska , m. , nicastro , f. , noble , m.  s. , nowak , m.  a. , primini , f.  a. , siemiginowska , a. , smith , r.  k. , & wise , m. 2006 , in ciao : chandra s data analysis system , vol . 6270 , 60                                    perlmutter , s. , aldering , g. , goldhaber , g. , knop , r.  a. , nugent , p. , castro , p.  g. , deustua , s. , fabbro , s. , goobar , a. , groom , d.  e. , hook , i.  m. , kim , a.  g. , kim , m.  y. , lee , j.  c. , nunes , n.  j. , pain , r. , pennypacker , c.  r. , quimby , r. , lidman , c. , ellis , r.  s. , irwin , m. , mcmahon , r.  g. , ruiz - lapuente , p. , walton , n. , schaefer , b. , boyle , b.  j. , filippenko , a.  v. , matheson , t. , fruchter , a.  s. , panagia , n. , newberg , h. j.  m. , couch , w.  j. , & project , t. s.  c. 1999 , astrophysical journal , 517 , 565          rasia , e. , meneghetti , m. , martino , r. , borgani , s. , bonafede , a. , dolag , k. , ettori , s. , fabjan , d. , giocoli c. , mazzotta , p. , merten , j. , radovich , m. , & tornatore , l. 2011 , new journal of physics , in press    riess , a.  g. , filippenko , a.  v. , challis , p. , clocchiatti , a. , diercks , a. , garnavich , p.  m. , gilliland , r.  l. , hogan , c.  j. , jha , s. , kirshner , r.  p. , leibundgut , b. , phillips , m.  m. , reiss , d. , schmidt , b.  p. , schommer , r.  a. , smith , r.  c. , spyromilio , j. , stubbs , c. , suntzeff , n.  b. , & tonry , j. 1998 , astronomical journal , 116 , 1009                    rccccc 15000 & 0 & 1.19 & 0.15 & 1.20 & 0.15 + 15000 & 1 & 1.17 & 0.13 & 1.18 & 0.13 + 15000 & 2 & 1.13 & 0.12 & 1.13 & 0.11 + 15000 & 3 & 1.13 & 0.13 & 1.13 & 0.12 + 30000 & 0 & 1.16 & 0.12 & 1.16 & 0.12 + 30000 & 1 & 1.13 & 0.11 & 1.13 & 0.11 + 30000 & 2 & 1.10 & 0.10 & 1.10 & 0.10 + 30000 & 3 & 1.10 & 0.10 & 1.09 & 0.10 + 60000 & 0 & 1.13 & 0.10 & 1.14 & 0.10 + 60000 & 1 & 1.11 & 0.08 & 1.11 & 0.08 + 60000 & 2 & 1.08 & 0.06 & 1.08 & 0.06 + 60000 & 3 & 1.07 & 0.06 & 1.07 & 0.06 + 120000 & 0 & 1.12 & 0.11 & 1.12 & 0.12 + 120000 & 1 & 1.10 & 0.09 & 1.10 & 0.09 + 120000 & 2 & 1.07 & 0.07 & 1.07 & 0.07 + 120000 & 3 & 1.07 & 0.07 & 1.07 & 0.06 +
Abstract:

Answer: precision cosmology studies based on wide - field surveys of galaxy clusters benefit from constraints on intrinsic scatter in mass - observable relationships . in principle , two - parameter models combining direct measurements of galaxy cluster structural variation with mass proxies such as x - ray luminosity and temperature can be used to constrain scatter in the relationship between the mass proxy and the cluster s halo mass and to measure the redshift evolution of that scatter . one candidate for quantifying cluster substructure is the icm temperature inhomogeneity inferred from x - ray spectral properties , an example of which is , the ratio of hardband to broadband spectral - fit temperatures . in this paper we test the effectiveness of  as an indicator of scatter in the mass - temperature relation using 118 galaxy clusters simulated with radiative cooling and feedback . we find that , while  is correlated with clusters departures @xmath0 from the mean  relation , the effect is modest . 2x - mas