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: a consequence of the presence of cold dark matter in galactic halos is the formation of caustics . caustics are regions of infinite density in the limit where the dark matter particles have zero velocity dispersion . the leading dark matter candidates , namely axions and weakly interacting massive particles(wimps ) , have very small primordial velocity dispersions ( @xmath0 cm / s for a 100 gev wimp , @xmath1 cm / s for a 10 @xmath2ev axion @xcite ) and may be expected to form caustics with large density enhancements . if the dark matter is made up of wimps , the particles can annihilate producing gamma rays . it may be possible to detect this gamma ray flux by means of future experiments such as glast @xcite . caustics may possibly be revealed by this technique . caustics form because cold dark matter particles occupy a thin , three dimensional hypersurface in phase space @xcite . there are two kinds of caustics in galactic halos : outer and inner @xcite . the outer caustics are topological spheres surrounding galaxies . the inner caustics have a more complicated geometry which depends on the details of the dark matter angular momentum distribution @xcite . when the angular momentum distribution is dominated by a net rotational term , the inner caustics have the appearance of rings @xcite . observational evidence for caustic rings was found by analyzing the rotation curves of external galaxies @xcite and the rotation curve of the milky way @xcite . a triangular feature seen in the infrared astronomical satellite ( iras ) map of the milky way is also suggestive of the presence of a nearby inner caustic ring @xcite . recently @xcite , possible evidence was found for a caustic ring in the galaxy cluster cl0024 + 1654 . in this article , we discuss dark matter annihilation in caustic spheres and caustic rings and present simulated gamma ray sky maps .    the existence of dark matter caustics is sometimes disputed on the basis of large scale cosmological simulations . however , current simulations use particles with mass of order @xmath3 . such massive particles have spurious collisions which destroy fine details in the dark matter phase space hypersurface . thus the absence of caustics in most large scale simulations is probably due to the finite resolution of these simulations and the large particle masses involved . note that discrete flows and caustics are observed in @xmath4-body simulations when special care is taken . the simulations of stiff and widrow @xcite , which increase the resolution in the relevant regions of phase space , do show discrete flows . more recently , vogelsberger et al . @xcite employed a new technique to study small scale structures such as caustics . further progress in numerical techniques may enable the study of caustic properties in galactic halos using @xmath4-body simulations . in cold dark matter cosmology , galactic halos form through a process of  hierarchical clustering " . small halos form first . larger halos , such as that of the milky way , are the outcome of past mergers , the accretion of smaller halos of various sizes , and the accretion of unclustered dark matter . this general scenario must surely be correct without , however , telling us anything precise . the critical issue for the existence of a caustic is the velocity dispersion of the flow that would form it . to be specific , consider the flow of dark matter that the milky way halo is accreting _ today_. if the dark matter in this flow has not clustered on any small scales yet , its velocity dispersion has the primordial value mentioned earlier , of order 0.03 cm / s for a 100 gev wimp . if it has clustered already , its effective velocity dispersion is the velocity dispersion of the structure it has become part of . the largest structures for which we have evidence of recent accretion are small sattelites , such as the sagittarius dwarf , with velocity dispersion of order 10 km / s . so one may conservatively assume that the flow of wimp dark matter accreting onto the galaxy today has variable velocity dispersion between 0.03 cm / s and 10 km / s . dark matter caustics are generically surfaces on one side of which the density diverges , in the limit of zero velocity dispersion , as @xmath5 where @xmath6 is the distance to the surface and @xmath7 is a constant called the  fold coefficient " . for finite velocity dispersion @xmath8 , the dark matter caustics in a galactic halo are smoothed over a distance @xmath9 where @xmath10 is the flow velocity ( of order 300 km / s ) and @xmath11 is the overall size of the halo ( of order 300 kpc ) . the velocity dispersion cuts off the divergence of the density at the caustic , i.e. the @xmath12 behaviour stops for @xmath13 . for the above range of velocity dispersion , @xmath14 ranges from @xmath15 cm to 10 kpc . since outer caustics have size of order hundreds of kpc , they are not erased by the expected velocity dispersion . inner caustics which have size of order 10 kpc are smoothed out only when the flow is entirely composed of sattelites with velocity dispersion of order 10 km / s . as was mentioned , evidence for caustic rings of dark matter in the milky way was found in the existence and distribution of sharp rises in the galactic rotation curve and the presence of a triangular feature in the iras map of the galactic plane @xcite . this evidence implies an upper limit on the velocity dispersion of order 50 m / s , or an upper limit on the smoothing scale @xmath14 of order @xmath16 pc . the wimp annihilation rate per unit surface of caustic is proportional to the integral @xmath17 this equation shows that caustics enhance the annihilation signal only by a logarithmic factor @xcite which varies from @xmath18 to 20 , depending on the velocity dispersion .    in this article , we limit ourselves to the simple cases of a spherically symmetric outer caustic and an axially symmetric inner caustic ring . our symmetry assumptions are made solely for the purpose of computational ease . in previous work @xcite , we verified that caustics are stable under perturbations , in particular perturbations that break any assumed symmetry . moreover , the symmetry assumptions are consistent with the observational evidence mentioned earlier . previous work on the possibility of detecting dark matter caustics by indirect means include @xcite . caustics and cold flows are also relevant to wimp and axion direct detection experiments @xcite . the number of photons received by a detector per unit area , per unit time , per unit solid angle and per unit energy due to wimp annihilation is given by the expression@xcite @xmath19 the term @xmath20 depends only on the particle physics : @xmath21 @xmath22 is the wimp annihilation cross section times the relative velocity averaged over the momentum distribution of the wimps , @xmath23 is the wimp mass , @xmath24 is the number of photons produced per annihilation per unit energy @xmath25 in the annihilation channel @xmath26 , and @xmath27 is the branching fraction of channel @xmath26 . the term @xmath28 is called the emission measure and depends solely on the dark matter distribution . @xmath29    where @xmath30 is the distance measured along the line of sight(@xmath31 ) in the direction of the galactic co - ordinates @xmath32 , and @xmath33 is the dark matter density . since all detectors have a finite resolution , the emission measure is often expressed in terms of @xmath34 , which represents an average over the resolution of the detector . in this article , we limit ourselves to computing @xmath28 .    as an example , let us compute the emission measure due to a smooth isotropic halo with the density profile @xmath35 consider an observer at a distance @xmath36 from the center of the halo , and a line of sight making an angle @xmath37 with the line joining the observer with the center ( we choose @xmath38 to be the direction of the halo center ) . @xmath39 is the core radius and keeps the density finite . we choose @xmath40 gev/@xmath41 when @xmath42 kpc . the emission measure @xmath43 is plotted as a function of angle @xmath37 in fig . [ figure1 ] , for the cases @xmath44 and 4 kpc . outer caustics are topological spheres . they are fold catastrophes with the universal density profile @xmath45 where @xmath46 is the dark matter density , @xmath47 is the distance from the center of the sphere , measured along the radial direction , @xmath11 is the radius of the sphere , @xmath7 is a constant for the particular caustic , and @xmath48 is the unit step function . for the self - similar infall model @xcite of the milky way halo , the radii of the first @xmath49 outer caustics @xmath50 kpc @xcite and the corresponding values of @xmath51 @xmath52 @xcite . we will assume a thickness @xmath53 for the caustic sphere , which sets a density cutoff equal to @xmath54 . consider a spherical outer caustic of radius @xmath11 and thickness @xmath53 . let an observer be located at a distance @xmath36 from the center of the sphere @xmath55 . consider a line of sight passing through the sphere making an angle @xmath37 with the line joining the observer with the center of the sphere(@xmath38 is the direction of the center ) . the emission measure @xmath43 is plotted in fig . [ figure2 ] for @xmath56 , for the two cases @xmath57 and @xmath58 . the maximum value of @xmath59 occurs when the line of sight is tangent to the sphere . [ figure3 ] shows a simulated sky map . let us now consider the case where the observer is within the caustic sphere . we once again assume a sphere of radius @xmath11 and thickness @xmath53 . the observer is located at a distance @xmath36 from the center(@xmath60 ) . as before , @xmath37 is the angle made by the line of sight with the line joining the observer with the center ( @xmath61 gives the direction of the center ) . [ figure4 ] shows the emission measure plotted as a function of angle @xmath37 for the three cases @xmath62 and @xmath63 . [ figure5 ] ( a),(b),(c ) show simulated skymaps .    in both cases , the emission measure for the first few outer caustics of a halo like the milky way is quite small compared to that obtained from the smooth halo of eq . [ smooth ] , for @xmath64 . however , @xmath43 would be larger for caustics of smaller radii and could possibly make a significant contribution to the total flux @xcite . the inner caustics are made up of sections of the higher order catastrophes . as mentioned in the introduction , the inner caustics resemble rings when the dark matter angular momentum distribution is dominated by a net rotational term . here , we assume that this is the case . a caustic ring is a closed tube whose cross section has three cusps . consider a caustic ring of radius @xmath39 , horizontal extent @xmath65 and vertical extent @xmath66 . see @xcite for details . we choose cylindrical co - ordinates @xmath67 . let us define the two dimensionless variables @xmath68 and @xmath69 . for simplicity , we assume axial symmetry about the @xmath70 axis and reflection symmetry about the @xmath71 plane . with these assumptions , we can obtain an analytic expression for the density @xmath33 close to a caustic ring , in the limit @xmath72 . for the special case of @xmath71 , the dark matter density is given by @xcite @xmath73 @xmath74 is a constant for the particular caustic ring . for the self - similar infall model of the milky way halo @xcite , @xmath75 for the first few @xmath76 dark matter flows . for the more general case of non - zero z , the density @xmath77 takes the form @xcite : @xmath78 where @xmath79 are the real roots of the quartic equation @xmath80 note that the assumption of reflection symmetry about the @xmath71 plane implies that the dark matter density is a function of @xmath81 . see appendix [ app ] for a derivation of eqs . [ d1 ] - [ t ]      consider an observer located at a distance @xmath36 from the center of a circular caustic ring with @xmath82 and @xmath83 . [ figure6 ] shows simulated skymaps of the emission measure . the caustic ring is aligned with the galactic plane in ( a ) and inclined to the plane in ( c ) . the emission measure is largest when the line of sight is tangent to the ring . the tricusp shape of the cross section is clearly visible in ( b ) . let us now consider the case where the observer is located within the tricusp shaped tube of a circular caustic ring . the observer is in the plane of the ring at a distance @xmath36 from the center with @xmath84 . we also set @xmath82 . [ figure7 ] shows a simulated skymap of the emission measure . the bright ring and the line at @xmath85 occur when the lines of sight pass close to the cusps . we see that for the inner caustics , the emission measure can be much larger than that obtained from a smooth halo ( except near the center of the halo ) . for a comparison of the annihilation flux with the measured gamma ray background , we refer the reader to @xcite . we have discussed dark matter annihilation in outer and inner caustics and provided simulated skymaps of the emission measure . in section [ section2 ] , we considered spherical outer caustics . when the observer is located outside the sphere , the emission measure @xmath43 is largest when the line of sight is tangent to the sphere , falling off rapidly to zero thereafter ( figs [ figure2],[figure3 ] ) . when the observer is within the sphere , the emission measure only changes gradually with the angle @xmath37 ( figs . [ figure4],[figure5 ] ) .    in section [ section3 ] , we considered circular , symmetric tricusp ring caustics . when the observer is outside the caustic ring , the emission measure is largest when the line of sight is tangent to the ring . looking tangentially , one may hope to identify the tricusp cross section of the ring[fig . [ figure6](b ) ] . when the observer is within the tricusp , the emission measure is largest when the lines of sight pass close to the cusps , resulting in the bright ring seen in fig . [ figure7 ] . we find that for inner caustics , @xmath43 can be significantly larger than that produced by a smooth halo , especially in directions away from the halo center . consider a single cold flow of dark matter particles falling in and out of the halo , forming a caustic . let us consider cylindrical co - ordinates @xmath86 where @xmath87 . for simplicity , we assume axial symmetry about the @xmath88 axis and reflection symmetry about the @xmath89 plane . let @xmath90 be the caustic ring radius , defined as the point of closest approach of the particles in the @xmath89 plane ( the particles with the most angular momentum are in this plane ) . consider a reference sphere such that each particle in the flow passes through the sphere once . we can then assign to each particle in the flow , a two parameter label @xmath91 which serves to identify the particle . we choose @xmath92 where @xmath93 is the polar angle of the particle at the time it crossed the reference sphere . @xmath94 is the time when the particles just above the @xmath89 plane cross this plane . see @xcite for details . we may expand @xmath95 and @xmath88 in a taylor series @xcite , @xmath96 where @xmath97 and the ring radius @xmath39 are constants for a given flow . in terms of these constants , we can express the horizontal ( @xmath65 ) and vertical ( @xmath66 ) extents of the tricusp : @xmath98 let us define the dimensionless variables @xmath99 and @xmath70 @xmath100 using eq . [ pq ] and eq . [ xz1 ] , we may express eq . [ rhoz ] in terms of @xmath99 and @xmath70 @xmath101 where @xmath102 . the density at points close to the caustic is given by @xcite @xmath103 where @xmath104 is the absolute value of the two dimensional jacobian determinant @xmath105 @xmath106 \right |   \label{d2}\ ] ] and the sum is over the individual dark matter flows that exist at that point . we use the self - similar infall halo model with angular momentum @xcite to estimate the mass infall rate @xmath107 @xmath74 is a parameter that characterizes the density of the flow , @xmath10 is the speed of the particles forming the caustic , and @xmath108 is the rotation velocity of the halo . when @xmath109 , we have two possible solutions : @xmath112 corresponding to the two flows that exist at each point . using eq.[d ] , [ d2 ] and [ infall ] , and approximating @xmath113 and @xmath114 we find the dark matter density at points close to the caustic @xmath115              we now consider the general case of computing the dark matter density at points outside the @xmath71 plane . for @xmath120 , neither @xmath121 nor @xmath122 can be zero , so we may express @xmath121 in terms of @xmath122 in eq . [ xz2 ] , to obtain a quartic equation for @xmath122 : @xmath123 the dark matter density is given by @xmath124 where @xmath79 are the real roots of eq . [ quartic ] . the number of real roots is given by the sign of the discriminant @xmath125.\ ] ] there are four real roots if @xmath126 and two if @xmath127 .    to solve eq . [ quartic ] , let us make the substitution @xmath128 , thereby eliminating the cubic term in eq . [ quartic ] : @xmath129 where @xmath130 . we may solve eq . [ nocubic ] by expressing it as the product of two quadratics : @xmath131 @xmath132 solves the cubic @xmath133 ^ 3 - 2 ( x + \frac{1}{2 } ) \left [ c^2 \right ] ^2 + \left [ \left ( x+\frac{1}{2 } \right ) ^2 + 4\gamma \right ] c^2 - x^2 = 0 \label{cubic}\ ] ] and @xmath134 . \label{d}\ ] ] let us define the two variables @xmath135 and @xmath4 @xmath136 @xmath132 has one real root @xmath137 if @xmath138 and three real roots @xmath139 for @xmath140 if @xmath141 . @xmath142 may be obtained using any of the values of @xmath143 . the desired solution of the quartic eq . [ quartic ] is then @xmath144 once the real roots @xmath79 are known , eq . [ dm_density ] and @xmath145 may be used to obtain the dark matter density . arnold , s.f . shandarin , y.b . zeldovich , geophysical and astrophysical fluid dynamics , * 20 * , 111 ( 1982 ) ; s.f . shandarin and y.b . zeldovich , rev . 61 , 185 ( 1989 ) ; p. sikivie and j.r . ipser , phys . b * 291 * , 288 ( 1992 ) ; a. natarajan and p. sikivie , phys . rev . d * 72 * , 083513 ( 2005 )                c.j . hogan , phys . d 64 , 063515 ( 2001 ) ; l. bergstrm , j. edsjo and c. gunnarsson , phys . rev d , * 63 * 083515 ( 2001 ) ; c. charmousis , v. onemli , z. qiu and p. sikivie , phys . rev d**67 * * , 103502 ( 2003 ) ; r. gavazzi , r. mohayaee and b. fort , astron . & astrophys . , * 445 * , 43 ( 2006 ) ; v. k. onemli , phys . rev . d 74 , 123010 ( 2006 ) ; r. mohayaee and s.f . shandarin , mon . not . * 366 * , 1217 ( 2006 ) ;          g. gelmini and p. gondolo , phys . rev d * 64 * , 023504 ( 2001 ) ; a.m. green , phys . rev . d 63 , 103003 ( 2001 ) ; j. d. vergados , phys . rev . d 63 , 063511 ( 2001 ) ; f.s . ling , p. sikivie and s. wick , phys . rev . d 70 , 123503 ( 2004 ) ; l.d . duffy et al , phys . rev . d 74 , 012006 ( 2006 ) ; c. savage , k. freese and p. gondolo , phys . rev . d 74 , 043531 ( 2006 )
Abstract:

Answer: dark matter caustics are small scale , high density structures believed to exist in galaxies like ours . if the dark matter consists of weakly interacting massive particles , these caustics may be detected by means of the gamma rays produced by dark matter particle annihilation . we discuss particle annihilation in outer and inner caustics and provide sky maps of the expected gamma ray distribution .


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: a continuous water jet falling gently downward through a tap and its subsequent breakup into discontinuous drops is an everyday - witnessed common phenomenon . there have been many attempts theoretically as well as experimentally to understand the physical mechanism behind the breaking up of jets into drops instead of continuing to keep reducing its cross - sectional area forever . however , the complete understanding of this macroscopic phenomenon has eluded so far . the scientific study of this phenomenon has a long history . savart ( 1833 ) experimentally observed the instability of liquid jets that clearly indicated that the liquid column ( jet ) undergoes changes before breaking up into disjoint droplets . savart s experimental observations suggested that the liquid column develops necks and bulges in the form of waves which ultimately leads to the formation of disjoint droplets out of the continuous column . this problem of instability of liquid jets was taken up by plateau and then by rayleigh ( 1879 ) and developed a theory which came to be known as plateau - rayleigh theory . rayleigh@xcite showed that perturbations of wavelengths larger than a certain value ( larger than the jet diameter ) grow rapidly with time which ultimately make the column unstable ( rayleigh instability ) against formation of droplets based on surface energy considerations of inviscid liquids . chandrasekhar@xcite later extended the theory to viscous liquids . many experimental investigations have been conducted to examine the validity of plateau - rayleigh theory@xcite . the investigations on the instability of water jet and related phenomena have been extensively reviewed@xcite . the study continues to be of current interest@xcite . the plateau - rayleigh theory , describes how the response of perturbations of certain wavelengths keep growing with time whereas those with lower wavelengths decay with time . it is the fastest growing perturbation that utimately breaks the jet . however , the theory is silent about the origin of the perturbations . the experiment of goedde and yuen , for example , applied external perturbations to study the length of the liquid jet ( measured from the root of the jet ) before it breaks up . however , even if no external perturbations are applied , one obtains a finite deterministic length of the liquid jet . in this case , obviously , the perturbation must have an origin at the root of the jet at the mouth of the nozzle , through which the liquid flows out to form the jet . naturally , the perturbations will have amplitudes of random sizes . as a consequence , the length of the jet should have random values . however , on the average , the length turns out to be a deterministically fixed number depending on the conditions of the experiment . in some recent works , in order to overcome this difficulty of fixing the origin of the perturbations , umemura and co - workers@xcite emphasized the idea that the perturbations get generated and sustained self - consistently . this is based on the observed fact that soon after the jet breaks up ( at one of the points on the lowermost neck ) the tip of the remaining column contracts to make its shape round once again to minimize the surface energy . the tip contraction gives rise to an upstream propagating capillary wave which upon reflection at the mouth of the nozzle moves downstream with doppler modified wavelengths . some of these waves with the right wavelength ( as discussed by rayleigh ) cause the liquid column to break producing another contraction of the tip of the column and the process repeats once again . we set up and conduct an experiment to verify the existence and effect of recoil capillary waves on the length of continuous water jet . this we accomplish by damping the recoil capillary wave using an earlier used method@xcite . the experimental set up ( figs.1,2 ) consists of two water tanks , the first one ( tank 1 ) is placed at a lower height to maintain the constant level of water ( i.e. , to keep same pressure throughout the experiment ) and the second one ( tank 2 ) is placed at a higher height to supply water at a constant rate to the tank 2 . both of these tanks are fitted with one tap each . the tap of the tank 2 is connected to the inlet of tank 1 by a rubber pipe . a glass nozzle is connected to a tap of the tank 1 by a rubber pipe through a flow control valve , which permits us to control the flow rate of water . the excess water in the tank 1 flows out through the level outlet to the lower reservoir tank ( tank 3 ) ( this tank is kept at the bottom level ) . the water from the tank 3 is pumped up to the tank 2 by using a water pump . a transparent rectangular beaker with a level outlet on one of its vertical sides is placed vertically below the nozzle so that the water jet falls directly on the water kept in the beaker . the water level in the beaker is maintained fixed by letting the excess water flow out through the level outlet of the beaker . the beaker is placed on a horizontal platform fitted to a travelling microscope so that the beaker can be smoothly moved vertically up and down and its position measured . a laser - pointer - and - detector arrangement is also fitted to the platform so that the horizontal laser beam incident normally on the vertical surface of the beaker and passes through the path of the water jet and then through the opposite surface of the beaker before it is collected by the detector . the laser - detector arrangement can be moved vertically and fixed as desired and the vertical position of the horizontal laser beam measured . a six digit counter mounted with a 24-hour clock is connected to the detector to count the number of drops and also to pin - point the first initial breaking point of the jet . our experimental set up is similar in essentials to that of goedde and yuen ( compare fig.1 of ref.@xcite ) . the ambient condition of temperature and humidity is controlled and allowed to settle down at desired values before the experiment is conducted and the condition is kept fixed throughout the duration of experiment . the water is issued ( jetted ) vertically downward through a long glass nozzle ( of length larger than about ten times its internal diameter , @xmath0 ) in to the water in the beaker through the stagnant ambient air atmosphere . as long as the jet remains continuous the detector remains quiescent . however , a breakup into a drop ( disjoint of the jet ) is detected as a count . the vertical distance between the nozzle exit and the laser beam gives the breakup length . initially , the laser beam is made to face the continuous jet by moving the platform up and then the platform is gradually lowered so that the beam position distance from the nozzle exit increases . at every position of the laser - detector arrangement the counts are recorded for two minutes for several times and their average calculated to obtain the average number of drop - counts per second . the platform is slowly and gradually lowered by small steps and at each position the experiment is repeated to obtain the average count rate . naturally , the count rate begins from zero , reaches a threshold at which the rate just begins to show nonzero value and then gradually keeps increasing as the platform is slowly lowered in small steps . the process continued till the rate reaches a maximum ( saturation ) value . throughout the above process the water flow rate is kept fixed . the same process is then repeated for several values of flow rates . note that after each change of flow rate the flow and the jet is allowed to become steady before the measurement process is begun .    for our purpose , we perform two distinct kinds ( sets ) of experiments . in one we let the laser beam pass grazing the water level in the beaker ( just about 0.2 cm above the water surface on the beaker ) . in the other the beam is kept at a height of about 1.5 cm above the water surface . in the first set of experiments , if the point of separation between a drop and the tip of the remaining jet is such that before the moving tip recoils it just touches the water surface , the process of jet - tip recoil is prevented and thus the recoil capillary wave propagating up the jet length is damped whereas in the other set of experiments this situation never occurs . the experiment is performed in an enclosure where the air current is minimized and , as stated earlier , the temperature and humidity are kept fixed for all sets of measurements . thus , the enclosure inside the lab is made sure to be stagnant and the external disturbance on the jet is minimized . ) as a function of distance of the laser beam from the nozzle exit at the flow rate of 55.5 cc / min and @xmath1 mm for the two sets of experiments ( laser beam close to the water surface and far from the surface . the inset shows a magnified graph to show the first ( jet ) breakup point.,width=377,height=264 ]    ) as a function of distance of the laser beam from the nozzle exit ( same as fig.3 but only for the set close to the surface ) . it also shows the corresponding unnormalized breakup length distribution.,width=377,height=264 ]    figure 3 shows the average number of counts ( drops ) per second , for a water flow rate of 55.5 cc / min and the nozzle inner diameter of 1.71 mm , as the vertical distance of the laser beam from the lower tip of the nozzle is increased . each point on the graph corresponds to an average of three set of counts for a duration of two minutes each . the counts begin from zero ( indicating continuous jet ) to a saturation value ( indicating the finality of the breakup process ) . from this data we calculate the distribution , @xmath2 , of breakup lengths . the probabilty density of breakup lengths is roughly calculated as the derivative of the count rate with respect to breakup length , fig . 4 . in fig.4 the distribution @xmath2 is not normalized . the normalized distribution @xmath3 will be obtained as @xmath4 . the distribution provides information about the mean value of breakup lengths . the mean values are plotted in fig.5 for the same nozzle of inner diameter 1.71 mm but for different flow rates for the two sets of experiments ( with and without damping ) . the two sets of mean values lie in a band . the approximate equality of the two mean values is understandable , however , because once the water level in the observation beaker is `` well below '' the breakup point , it is as good as being far below and equivalent to a point of the set without damping . the inset of fig.3 gives a magnified diagram of the curves of fig . 3 in the lower range of the distance of laser beam from the nozzle exit . from this graph one can obtain the first ( initial ) breakup lengths in the two cases ( with and without damping ) and are plotted in fig . 5 as a function of water flow rates . the figure clearly shows that the first breakup length with the damped recoil capillary wave is larger than that without the damping for all values of flow rates measured . this indicates that the recoil capillary waves do exist and they aid in making the water jet unstable against breakup . the mean breakup lengths plotted in fig.5 , though not exactly identical to the earlier reported results ( for example , p. 104 of ref.@xcite ) , the qualitative features are very similar , showing various regimes of jet breakup . however , information on the first breakup length is entirely new . it shows that the first breakup lengths peak at a smaller reynolds number ( mean re=1079.3 ) and weber number ( we=7.54 ) than in case of the mean breakup lengths ( re=1392.6 , mean we=12.51 ) . our measurements are done at the temperature of ( 25@xmath5.5)@xmath6c and at relative humidity of 80@xmath53% . however , in order to calculate re and we , we have used the tabulated values of surface tension of water @xmath7nm@xmath8 and coefficient of viscosity @xmath9 kg m@xmath8s@xmath8 .        the inset of fig.3 gives a magnified diagram of the curves of fig . 3 in the lower range of the distance of laser beam from the nozzle exit . from this graph one can obtain the first ( initial ) breakup lengths in the two cases ( with and without damping ) and are plotted in fig . 5 as a function of water flow rates . the figure clearly shows that the first breakup length with the damped recoil capillary wave is larger than that without the damping for all values of flow rates measured . this indicates that the recoil capillary waves do exist and they aid in making the water jet unstable against breakup . we have performed a very simple experiment the results of which may be considered to directly give a proof of the existence of recoil capillary waves and their effect on the jet breakup length . though the experiment is not so precise and sophisticated the qualitative features shown are unmistakable . a high speed photographic measurement may help in arriving at a result with more confidence . we hope to obtain a conclusive result once we have the required facility . lord , j. w. s. rayleigh , proc . london math . soc . , 4 ( 1879 ) . rayleigh , , vol . ii , pp . 351 - 375 , dover publications , new york , 1945 . s. chandrasekhar , , dover publications , new york , 1961 . e. f. goedde , and m. c. yuen , j. fluid mech . , 495 ( 1970 ) . j. eggers , rev . , , 865 ( 1997 ) . j. eggers , and e. villermaux , rep . prog . , 036601 ( 2008 ) . s. p. lin , , cambridge university press , 2010 . s. rajendran , , m. s. thesis ( in mechanical engineering ) , university of cincinnati , 2012 . a. umemura , phys . e , 046307 ( 2011 ) . a. umemura , s. kawanabe , s. suzuki , and j. osaka , phys . e , 036309 ( 2011 ) . mahato , k. singh , and s. taraphdar , solid state physics ( india ) , 167 ( 2001 ) .
Abstract:

Answer: breakup of water jets under gravity is a common - place phenomenon . the role of surface tension in the instability of water jets was recognized by rayleigh and the theory propounded goes by the name of plateau - rayleigh theory . the necks and bulges down along the jet - length that are created by perturbation waves of wavelengths larger than a certain value keep growing with time and ultimately cause the jet to breakup into drops . the effect of perturbation waves have been investigated experimentally and found to confirm the essentials of the theory . however , there is no unanimity about the origin of these perturbation waves . recently , the idea of recoil capillary waves as an important source of the perturbation waves has been emphasized . the recoil of the end point of the remaining continuous jet at its breakup point is considered to travel upward as a recoil capillary wave which gets reflected at the mouth of the nozzle from which the jet originates . the reflected capillary wave travels along the jet downward with its doppler shifted wavelength as a perturbation wave . we set up an experiment to directly verify the existence and effect of the recoil capillary waves and present some preliminary results of our experiment .


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: power control is an effective tool for interference management in cellular , ad - hoc , and cognitive underlay networks @xcite . the prevailing formulation of power control aims to use a minimum total transmission power to support all links in an interference network at their desired sinr targets . a longstanding issue associated with power control is that the problem often becomes _ infeasible _ , i.e. , it is not possible to simultaneously support all links in the network at their sinr targets . in this case , we must adopt a joint power and admission control approach to _ selectively _ remove some links from the network so that the remaining ones can be simultaneously supported at their desired sinr levels . our goal is to maximize the number of simultaneously supportable links at their required sinr targets while using a minimum total transmission power . theoretically , the joint power and admission control problem is known to be np - hard to solve to global optimality @xcite and to approximate to constant ratio global optimality @xcite , so various heuristic algorithms @xcite have been proposed for this problem . among them , the reference  @xcite proposed a convex approximation - based algorithm for the joint power and admission control problem . instead of directly solving the original np - hard problem , the basic idea of the proposed linear programming deflation ( lpd ) algorithm in @xcite is to approximate the problem by an appropriate convex problem . the solution to the approximation problem  can be used to check the feasibility of the original problem and guide link s removals . the removal procedure is terminated until all the remaining links in the network are simultaneously supportable . the recent work  @xcite developed another lp approximation - based new linear programming deflation ( nlpd ) algorithm for the joint power and admission control problem . in @xcite , the joint power and admission control problem is first equivalently reformulated as a sparse @xmath0-minimization problem and then its @xmath1-convex approximation is used to derive a lp , which is different from the one in @xcite . again , the solution to the derived lp can guide an iterative link removal procedure , and the removal procedure is terminated if all the remaining links in the network are simultaneously supportable . based on the sparse @xmath0-minimization reformulation in @xcite , this paper proposes a new deflation algorithm based on @xmath2 ( @xmath3 ) norm minimization for the joint power and admission control problem . compared to the @xmath1-minimization problem , the @xmath2 norm minimization problem is closer to the original @xmath0-optimization problem . the @xmath6-approximation problem is solved by applying the efficient interior - point algorithm in @xcite to solve its equivalent reformulation . numerical results show that the proposed algorithm compares favorably with the existing approaches @xcite in terms of the number of supported links , the total transmission power , and the cpu time . _ notations : _ we adopt the following notations in this paper . we denote the index set @xmath7 by @xmath8 . lowercase boldface and uppercase boldface are used for vectors and matrices , respectively . for a given vector @xmath9 the notations @xmath10 @xmath11_k$ ] and @xmath12_k|^p~(0\leq p<1)\footnote{strictly speaking , $ \|{\mathbf{x}}\|_p$ with $ 0\leq p<1 $ is not a norm , since it does not satisfy the triangle inequality . however , we still call it $ p$ norm for convenience in this paper.}$ ] stand for its maximum entry , its @xmath13-th entry , and its @xmath2 norm , respectively . in particular , when @xmath14 @xmath15 stands for the number of nonzero entries in @xmath16 finally , we use @xmath17 to represent the vector of an appropriate size with all components being one and @xmath18 to represent the identity matrix of an appropriate size , respectively . consider a @xmath19-link ( a link corresponds to a transmitter - receiver pair ) interference channel with channel gains @xmath20 ( from transmitter @xmath21 to receiver @xmath13 ) , noise power @xmath22 sinr target @xmath23 and power budget @xmath24 for @xmath25 . denote the power allocation vector by @xmath26 and the power budget vector by @xmath27 . treating interference as noise , we can write the sinr at the @xmath13-th receiver as @xmath28    the joint power and admission control problem can be mathematically formulated as a two - stage optimization problem @xcite . specifically , the first stage maximizes the number of admitted links:@xmath29 \mbox{s.t . } & \displaystyle { { \mbox{sinr}}}_k\geq \gamma_k,~k\in\cal s\subseteq\cal k,\\ & \bm{0 } \leq { \mathbf{p}}\leq \bar{\mathbf{p}}. \end{array}\ ] ] we use @xmath30 to denote the optimal solution for problem and call it _ maximum admissible set_. notice that the solution for might not be unique . the second stage minimizes the total transmission power required to support the admitted links in @xmath31 : @xmath32 due to the special choice of @xmath33 power control problem is feasible and can be efficiently and distributively solved by the foschini - miljanic algorithm @xcite . since the developed algorithm in this paper follows the similar idea as the nlpd algorithm in @xcite , we first briefly review the nlpd algorithm in this section . the basic idea of the nlpd algorithm is to update the power and check whether all links can be supported or not . if the answer is yes , then terminate the algorithm ; else drop one link from the network and update the power again . the above process is repeated until all the remaining links are supported . we begin with the introduction of an equivalent normalized channel on which the nlpd algorithm is based . in particular , we use @xmath34 with @xmath35  to denote the normalized power allocation vector , and use @xmath36 with @xmath37  to denote the normalized noise vector . we denote the normalized channel matrix by @xmath38 with the @xmath39-th entry @xmath40 in fact , @xmath41 is the normalized channel gain . it is simple to check that @xmath42 if and only if @xmath43_k \geq 0.$ ]    in @xcite , we reformulate the two - stage joint power and admission control problem and as a single - stage optimization problem @xmath44 where @xmath45 and @xmath46_k$ ] measures the excess transmission power @xcite that the transmitter of link @xmath13 needs in the normalized channel in order to be served with its desired sinr target ( assuming all other links keep their transmission powers unchanged ) . notice that the formulation is capable of finding the maximum admissible set with minimum total transmission power . since problem is still np - hard , we further consider its @xmath1-convex approximation ( equivalent to a lp ) @xmath47 by solving , we know whether all links in the network can be simultaneously supported or not . if not , we drop the link @xmath48_j+|a_{jk}|\left[{\mathbf{q}}_e\right]_k\right)\right\}.\ ] ]    an easy - to - check necessary condition @xmath49 for all links in the network to be simultaneously supported is also derived in @xcite , where @xmath50 @xmath51 and @xmath52 . the necessary condition allows us to iteratively remove strong interfering links from the network . in particular , we remove the link @xmath53 according to the scheme @xmath54 until becomes true . the nlpd algorithm can be described as follows . in this section , we develop a new deflation algorithm based on @xmath6-minimization for the joint control problem and . as seen in section [ secreview ] , the original @xmath0-minimization problem is _ successively _ approximated by the @xmath1-minimization problem in the nlpd algorithm . intuitively , the @xmath2 ( @xmath3 ) norm minimization problem @xmath55 should approximate better than . this is the motivation for the development of the new deflation algorithm based on the @xmath2 norm minimization for the joint power and admission control problem .    comparing the @xmath5 norm minimization problem and the @xmath2 norm minimization problem , we see that problem is convex while problem is nonconvex ( for its objective function is nonconvex with respect to @xmath56 ) . generally speaking , convex problems are relatively easy to solve , while nonconvex optimization problems are difficult to solve . however , not all nonconvex problems are hard since the lack of convexity may be due to an inappropriate formulation . in fact , many nonconvex optimization problems admit a convex reformulation . therefore , convexity is useful but unreliable to test the computational intractability of an optimization problem . a more robust tool is the computational complexity theory @xcite . we show that problem is strongly np - hard . the proof is based on a polynomial time reduction from the max-2unanimity problem , which is shown to be strongly np - hard in @xcite .    [ thm1 ] the @xmath6-minimization problem is strongly np - hard if @xmath57    the complexity result in theorem [ thm1 ] motivates us to approximately solve problem . next , we first give a reformulation of problem , and then propose to use the interior - point algorithm developed in @xcite to solve it . [ thm2 ] the @xmath6-minimization problem can be equivalently reformulated as @xmath58_k^p + \alpha\bar{\mathbf{p}}^t{\mathbf{q}}\\ \mbox{s.t . } & \displaystyle { \mathbf{q}}_e={\mathbf{c}}-{\mathbf{a}}{\mathbf{q}},\\&\mathbf{0}\leq { \mathbf{q}}\leq { \mathbf{e}},~{\mathbf{q}}_e\geq\mathbf{0}. \end{array}\ ] ]    a rigorous proof of theorem [ thm2 ] shall be given in the journal version . here we just shed some light on why it does not harm optimality to restrict @xmath59 notice that problem is equivalent to @xmath60 thus , to show the equivalence of and , it suffices to show that any optimal solution @xmath61 of always satisfies @xmath62 . in fact , assume the contrary that @xmath63 where @xmath64@xmath65_k>0\,\}$ ] , @xmath66@xmath65_k=0\,\}.$ ] then by the balancing lemma ( see ( * ? ? ? * lemma 1 ) ) , we can appropriately reduce the power of links in @xmath67 so that both the first term and the second term in the objective of are strictly decreased . based on theorem [ thm2 ] , by introducing a slack variable @xmath68 to problem , we see that problem is actually equivalent to @xmath69_k^p + \alpha\bar{\mathbf{p}}^t{\mathbf{q}}\\ \mbox{s.t . } & \displaystyle { \mathbf{q}}_e={\mathbf{c}}-{\mathbf{a}}{\mathbf{q}},\,{\mathbf{s}}+{\mathbf{q}}={\mathbf{e}}\\&{\mathbf{q}}\geq\mathbf{0},\,{\mathbf{q}}_e\geq\mathbf{0},\,{\mathbf{s}}\geq\mathbf{0}. \end{array}\]]now , we can apply the interior - point algorithm in @xcite to solve problem . similar to @xcite , we can prove that the potential reduction interior - point algorithm returns an @xmath70-kkt @xcite or @xmath70-global solution of problem in no more than @xmath71 iterations . one may ask why we wish to use interior - point algorithms to solve problem ? the reasons are the following . first , the objective function of problem is _ differentiable _ in the interior feasible region . moreover , we are actually interested in finding a sparse solution @xmath56 of problem ; if we start from a solution , some of whose entries are already zero , then it is very hard to make it nonzero . in contrast , if we start from an interior point , the interior - point algorithm may generate a sequence of interior points that bypasses solutions with the wrong zero supporting set and converges to the true one . this is exactly the idea of the interior - point algorithm developed in @xcite for the nonconvex quadratic programming . the proposed @xmath2 norm minimization deflation ( pnmd ) algorithm is given as follows . there are two unclear points in the pnmd algorithm . one is how to compute the parameter @xmath72 in problem , and the other is which removal strategy will be used in the admission control step . next , we make clear of these two points , i.e. , we shall use @xmath72 given in and the removal strategy in the new deflation algorithm .    in the nlpd algorithm @xcite , the parameter @xmath72 is given by @xmath73 where @xmath74 are two constants , and @xmath75 is determined by the equivalence between problem and the joint problem and , and @xmath76 is determined by the so - called `` never - over - removal '' property . since the @xmath6-minimization problem is closer to the @xmath0-minimization problem , we relax the parameter @xmath72 in to @xmath77where @xmath78 @xmath79 are three constants .    having obtained the solution @xmath80 of problem , we use the removal strategy called smart rule in @xcite to drop the link @xmath53 according to @xmath81 the above operation can be interpreted as removing the link with the largest interference plus noise footprint in the _ normalized _ network . we generate the same channel parameters as in @xcite in our numerical simulations , i.e. , each transmitter s location obeys the uniform distribution over a @xmath82 km  @xmath83  @xmath82 km square and the location of its corresponding receiver is uniformly generated in a disc with radius @xmath84 m ; channel gains are given by @xmath85 where @xmath86 is the euclidean distance from the link of transmitter @xmath21 to the link of receiver @xmath87 each link s sinr target is set to be @xmath88 and the noise power is set to be @xmath89 . the power budget of the link of transmitter @xmath13 is @xmath90 where @xmath91 is the minimum power needed for link @xmath13 to meet its sinr requirement in the absence of any interference from other links . the parameter @xmath2 in problem is set to be @xmath92 and the ones in are set to be @xmath93  and  @xmath94 . the number of supported links , the total transmission power , and the cpu time are the metrics we employ to compare the performance of the proposed pnmd algorithm with that of the lpd algorithm in @xcite , the algorithm ii - b in @xcite , and the nlpd algorithm in @xcite . all figures are obtained by averaging over @xmath95 monte - carlo runs . [ user1 ] , [ time1 ] , and [ power1 ] indicate that the pnmd algorithm can take less cpu time to support more links while with less total transmission power than the existing algorithms ( except the algorithm ii - b ) . as shown in fig . [ power1 ] , the algorithm ii - b transmits the least power among the tested algorithms . this is because the algorithm ii - b supports the least number of links ; see fig . [ user1 ] . in particular , compared to the nlpd algorithm ] , the proposed pnmd algorithm can support ( slightly ) more links with much less total transmission power , and at the same time takes less cpu time . the performance improvement of the proposed pnmd algorithm over the nlpd algorithm is mainly attributed to the @xmath6-approximation problem . the simulation results in fig . [ user1 ] and fig . [ power1 ] show that the admissible set @xmath96 obtained by the proposed pnmd algorithm based on the @xmath6-approximation problem is `` better '' than the admissible set @xmath97 obtained by the nlpd algorithm based on the @xmath1-approximation problem , i.e. , although the cardinality of the two admissible sets @xmath96 and @xmath97 is nearly equal to each other , it takes much less total transmission power to support the links in @xmath96 than to support the links in @xmath98 this is consistent with our intuition that the @xmath99 norm minimization problem is capable of approximating the @xmath0-minimization problem better than the @xmath1-minimization problem and the fact that the maximum admissible set for the joint power and admission control problem may not be unique . in this paper , we have developed a @xmath2 @xmath100 norm minimization deflation algorithm for the joint power and admission control problem . numerical simulations show the proposed algorithm outperforms state - of - the - arts in @xcite in terms of the number of supported links , the total transmission power , and the cpu time . i. mitliagkas , n. d. sidiropoulos , and a. swami , `` joint power and admission control for ad - hoc and cognitive underlay networks : convex approximation and distributed implementation , '' _ ieee trans . wireless commun . 41104121 , dec . 2011 . h.  mahdavi - doost , m.  ebrahimi , and a.  k.  khandani ,  `` characterization of sinr region for interfering links with constrained power , '' _ ieee trans . inf . theory _ , 56 , no .  6 , pp . 28162828 , jun . m. andersin , z. rosberg , and j. zander , `` gradual removals in cellular pcs with constrained power control and noise , '' _ wireless netw . 1 , pp .  2743 , mar . g.  j.  foschini and z. miljanic , `` a simple distributed autonomous power control algorithm and its convergence , '' _ ieee trans . 42 , no .  4 , pp . 641646 , nov . s. a. grandhi , j. zander , and r. yates , `` constrained power control , '' _ wireless personal commun . _ , vol . 1 , no . 4 , pp . 257270 , 1995 . f . liu , y .- h . dai , and z .- q . luo , `` joint power and admission control via linear programming deflation , '' _ ieee trans . signal process . 61 , no . 6 , pp . 13271338 , mar . liu , `` an efficient distributed joint power and admission control algorithm , '' in _ proc . 31th chinese control conference , signal processing and optimization session , _ july , 2012 , pp . liu and e. song , `` sample approximation - based deflation approaches for chance sinr constrained joint power and admission control . '' available online : http://arxiv.org/abs/1302.5973 . d. ge , x. jiang , and y. ye , `` a note on the complexity of @xmath101 minimization , '' to appear in _ math . _ available online : http://www.stanford.edu/  yyye / lpmin_v14.pdf . m. r. garey and d. s. johnson , _ computers and intractability : a guide to the theory of np - completeness_. sf , u.s.a . :  w. h. freeman and company , 1979 . y .- f . liu , y .- h . dai , and z .- q . luo , `` coordinated beamforming for miso interference channel : complexity analysis and efficient algorithms , '' _ ieee trans . signal process . 11421157 , mar . y. ye , `` on the complexity of approximating a kkt point of quadratic proramming , '' _ math . prog . , _ vol . 195211 , 1998 .
Abstract:

Answer: in an interference network , joint power and admission control aims to support a maximum number of links at their specified signal to interference plus noise ratio ( sinr ) targets while using a minimum total transmission power . in our previous work , we formulated the joint control problem as a sparse @xmath0-minimization problem and relaxed it to a @xmath1-minimization problem . in this work , we propose to approximate the @xmath0-optimization problem by a @xmath2 norm minimization problem where @xmath3 , since intuitively @xmath2 norm will approximate @xmath4 norm better than @xmath5 norm . we first show that the @xmath6-minimization problem is strongly np - hard and then derive a reformulation of it such that the well developed interior - point algorithms can be applied to solve it . the solution to the @xmath6-minimization problem can efficiently guide the link s removals ( deflation ) . numerical simulations show the proposed heuristic outperforms the existing algorithms .


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: torque oscillator ( sto ) has attracted much attention due to its potential uses for a microwave generator and a recording head of a high density hard disk drive . the self - oscillation of the sto was first discovered in an in - plane magnetized giant - magnetoresistive ( gmr ) system [ 1 ] . after that , the self - oscillation of the sto has been observed not only in gmr systems [ 2]-[6 ] but also in magnetic tunnel junctions ( mtjs ) [ 7]-[11 ] . the different types of sto have been proposed recently ; for example , a point - contact geometry with a confined magnetic domain wall [ 12]-[14 ] which enables us to control the frequency from a few ghz to a hundred ghz .    recently , kubota _ et al . _ experimentally developed the mgo - based mtj consisting of a perpendicularly magnetized free layer and an in - plane magnetized pinned layer [ 15],[16 ] . they also studied the self - oscillation of this type of mtj , and observed a large power ( @xmath0 @xmath1w ) with a narrow linewidth ( @xmath2 mhz ) [ 17 ] . these results are great advances in realizing the sto device . however , the relation between the current and the oscillation frequency still remains unclear . since a precise control of the oscillation frequency of the sto by the current is necessary for the application , it is important to clarify the relation between the current and the oscillation frequency .    in this paper , we derived the theoretical formula of the relation between the current and the oscillation frequency of the sto consisting of the perpendicularly magnetized free layer and the in - plane magnetized pinned layer . the derivation is based on the analysis of the energy balance between the work done by the spin torque and the energy dissipation due to the damping . we found that the oscillation frequency monotonically decreases with increasing the current by keeping the magnetization in one hemisphere of the free layer . the validity of the analytical solution was confirmed by numerical simulations . we also found that the field - like torque breaks the energy balance , and change the oscillation frequency . the shift direction of the frequency , high or low , is determined by the sign of the field - like torque . this paper is organized as follows . in sec . ii , the current dependence of the oscillation frequency is derived by solving the landau - lifshitz - gilbert ( llg ) equation . in sec . iii , the effect of the field - like torque on the oscillation behaviour is investigated . section iv is devoted to the conclusions . -axis are indicated . ] the system we consider is schematically shown in fig . [ fig : fig1 ] . we denote the unit vectors pointing in the directions of the magnetization of the free and the pinned layers as @xmath3 and @xmath4 , respectively . the @xmath5-axis is parallel to @xmath4 while the @xmath6-axis is normal to the film plane . the variable @xmath7 of @xmath8 is the tilted angle from the @xmath6-axis while @xmath9 is the rotation angle from the @xmath5-axis . the current @xmath10 flows along the @xmath6-axis , where the positive current corresponds to the electron flow from the free layer to the pinned layer . we assume that the magnetization dynamics is well described by the following llg equation : @xmath11 the gyromagnetic ration and the gilbert damping constant are denoted as @xmath12 and @xmath13 , respectively . the magnetic field is defined by @xmath14 , where the energy density @xmath15 is @xmath16 here , @xmath17 , @xmath18 , and @xmath19 are the saturation magnetization , the applied field along the @xmath6-axis , and the crystalline anisotropy field along the @xmath6-axis , respectively . because we are interested in the perpendicularly magnetized system , the crystalline anisotropy field , @xmath19 , should be larger than the demagnetization field , @xmath20 . since the llg equation conserves the norm of the magnetization , the magnetization dynamics can be described by a trajectory on an unit sphere . the equilibrium states of the free layer correspond to @xmath21 . in following , the initial state is taken to be the north pole , i.e. , @xmath22 . it should be noted that a plane normal to the @xmath6-axis , in which @xmath7 is constant , corresponds to the constant energy surface . the spin torque strength , @xmath23 in eq . ( [ eq : llg ] ) , is [ 18]-[20 ] @xmath24 where @xmath25 and @xmath26 are the cross section area and the thickness of the free layer . two dimensionless parameters , @xmath27 and @xmath28 ( @xmath29 ) , determine the magnitude of the spin polarization and the angle dependence of the spin torque , respectively . although the relation among @xmath27 , @xmath28 , and the material parameters depends on the theoretical models [ 20]-[22 ] , the form of eq . ( [ eq : h_s ] ) is applicable to both gmr system and mtjs . in particular , the angle dependence of the spin torque characterized by @xmath28 is a key to induce the self - oscillation in this system .    ) and ( [ eq : frequency ] ) . ] figure [ fig : fig2 ] ( a ) shows the steady state precession of the magnetization in the free layer obtained by numerically solving eq . ( [ eq : llg ] ) . the values of the parameters are @xmath30 emu / c.c . , @xmath31 koe , @xmath32 koe , @xmath33 nm@xmath34 , @xmath35 nm , @xmath36 mhz / oe , @xmath37 , @xmath38 , and @xmath39 , respectively [ 17 ] . the self - oscillation was observed for the current @xmath40 ma . although the spin torque breaks the axial symmetry of the free layer along the @xmath6-axis , the magnetization precesses around the @xmath6-axis with an almost constant tilted angle . the tilted angle from the @xmath6-axis increases with increasing the current ; however , the magnetization stays in the north semisphere ( @xmath41 ) . the dots in fig . [ fig : fig2 ] ( b ) show the dependence of the oscillation frequency on the current . as shown , the oscillation frequency monotonically decreases with increasing the current magnitude . let us analytically derive the relation between the current and the oscillation frequency . since the self - oscillation occurs due to the energy supply into the free layer by the spin torque , the energy balance between the spin torque and the damping should be investigated . by using the llg equation , the time derivative of the energy density @xmath15 is given by @xmath42 , where the work done by spin torque , @xmath43 , and the energy dissipation due to the damping , @xmath44 , are respectively given by @xmath45 ,    \label{eq : w_s}\ ] ] @xmath46 . \label{eq : w_alpha}\ ] ]    by assuming a steady precession around the @xmath6-axis with a constant tilted angle @xmath7 , the time averages of @xmath43 and @xmath44 over one precession period are , respectively , given by @xmath47    \cos\theta ,    \label{eq : w_s_ave } \end{split}\ ] ] @xmath48^{2 }    \sin^{2}\theta . \label{eq : w_alpha_ave}\ ] ] the magnetization can move from the initial state to a point at which @xmath49 . then , the current at which a steady precession with the angle @xmath7 can be achieved is given by @xmath50    \sin^{2}\theta . \label{eq : i_theta } \end{split}\ ] ] the corresponding oscillation frequency is given by @xmath51 . \label{eq : frequency}\ ] ] equations ( [ eq : i_theta ] ) and ( [ eq : frequency ] ) are the main results in this section . the solid line in fig . [ fig : fig2 ] ( b ) shows the current dependence of the oscillation frequency obtained by eqs . ( [ eq : i_theta ] ) and ( [ eq : frequency ] ) , where the good agreement with the numerical results confirms the validity of the analytical solution . the critical current for the self - oscillation , @xmath52 , is given by @xmath53 the value of @xmath54 estimated by using the above parameters is @xmath55 ma , showing a good agreement with the numerical simulation shown in fig . [ fig : fig2 ] ( a ) . the sign of @xmath54 depends on that of @xmath28 , and the self - oscillation occurs only for the positive ( negative ) current for the positive ( negative ) @xmath28 . this is because a finite energy is supplied to the free layer for @xmath56 , i.e. , @xmath57 . in the case of @xmath58 , the average of the work done by the spin torque is zero , and thus , the self - oscillation does not occur . it should be noted that @xmath59 in the limit of @xmath60 . this means the magnetization can not cross over the @xmath61-plane , and stays in the north hemisphere ( @xmath41 ) . the reason is as follows . the average of the work done by spin torque becomes zero in the @xmath61-plane ( @xmath62 ) because the direction of the spin torque is parallel to the constant energy surface . on the other hand , the energy dissipation due to the damping is finite in the presence of the applied field [ 21 ] . then , @xmath63 , which means the damping moves the magnetization to the north pole . thus , the magnetization can not cross over the @xmath61-plane . the controllable range of the oscillation frequency by the current is @xmath64 , which is independent of the magnitude of the applied field .    since the spin torque breaks the axial symmetry of the free layer along the @xmath6-axis , the assumption that the tilted angle is constant used above is , in a precise sense , not valid , and the @xmath6-component of the magnetization oscillates around a certain value . then , the magnetization can reach the @xmath61-plane and stops its dynamics when a large current is applied . however , the value of such current is more than 15 ma for our parameter values , which is much larger than the maximum of the experimentally available current . thus , the above formulas work well in the experimentally conventional current region .    contrary to the system considered here , the oscillation behaviour of an mtj with an in - plane magnetized free layer and a perpendicularly magnetized pinned layer has been widely investigated [ 23]-[26 ] . the differences of the two systems are as follows . first , the oscillation frequency decreases with increasing the current in our system while it increases in the latter system . second , the oscillation frequency in our system in the large current limit becomes independent of the @xmath6-component of the magnetization while it is dominated by @xmath65 in the latter system . the reasons are as follows . in our system , by increasing the current , the magnetization moves away from the @xmath6-axis due to which the effect of the anisotropy field on the oscillation frequency decreases , and the frequency tends to @xmath66 , which is independent of the anisotropy . on the other hand , in the latter system , the magnetization moves to the out - of - plane direction , due to which the oscillation frequency is strongly affected by the anisotropy ( demagnetization field ) . the macrospin model developed above reproduces the experimental results with the free layer of 2 nm thick [ 17 ] , for example the current - frequency relation , quantitatively . although only the zero - temperature dynamics is considered in this paper , the macrospin llg simulation at a finite temperature also reproduces other properties , such as the power spectrum and its linewidth , well . however , when the free layer thickness further decreases , an inhomogeneous magnetization due to the roughness at the mgo interfaces may affects the magnetization dynamics : for example , a broadening of the linewidth . the field - like torque arises from the spin transfer from the conduction electrons to the local magnetizations , as is the spin torque . when the momentum average of the transverse spin of the conduction electrons relaxes in the free layer very fast , only the spin torque acts on the free layer [ 19 ] . on the other hand , when the cancellation of the transverse spin is insufficient , the field - like torque appears . the field - like torque added to the right hand side of eq . ( [ eq : llg ] ) is @xmath67 where the dimensionless parameter @xmath68 characterizes the ratio between the magnitudes of the spin torque and the field - like torque . the value and the sign of @xmath68 depend on the system parameters such as the band structure , the thickness , the impurity density , and/or the surface roughness [ 22],[27]-[29 ] . the magnitude of the field - like torque in mtj is much larger than that in gmr system [ 30],[31 ] because the band selection during the tunneling leads to an insufficient cancellation of the transverse spin by the momentum average .     with ( a ) @xmath69 , ( b ) @xmath70 , and @xmath71 . the current magnitude is @xmath72 ma . ( d ) the time evolutions of @xmath73 for various @xmath68 . ] it should be noted that the effective energy density , @xmath74 satisfying @xmath75 $ ] , can be introduce to describe the field - like torque . the time derivative of the effective energy , @xmath76 , can be obtained by replacing the magnetic field , @xmath77 , in eqs . ( [ eq : w_s ] ) and ( [ eq : w_alpha ] ) with the effective field @xmath78 . then , the average of @xmath79 over one precession period around the @xmath6-axis consists of eq . ( [ eq : w_s_ave ] ) , ( [ eq : w_alpha_ave ] ) , and the following two terms : @xmath80 ,    \label{eq : w_s_flt_ave } \end{split}\ ] ] @xmath81 \\    & -    \frac{2 \alpha \beta \gamma m}{1+\alpha^{2 } }    \frac{\hbar \eta i}{2e \lambda msd }    \left (      \frac{1}{\sqrt{1-\lambda^{2}\sin^{2}\theta } }      -      1    \right ) \\    & \times    \left [      h_{\rm appl }      +      \left (        h_{\rm k }        -        4\pi m      \right )      \cos\theta    \right ]    \cos\theta . \label{eq : w_alpha_flt_ave } \end{split}\ ] ] the constant energy surface of @xmath76 shifts from the @xmath61-plane due to a finite @xmath82 , leading to an inaccuracy of the calculation of the time average with the constant tilted angle assumption . thus , eqs . ( [ eq : w_s_flt_ave ] ) and ( [ eq : w_alpha_flt_ave ] ) are quantitatively valid for only @xmath83 . however , predictions from eqs . ( [ eq : w_s_flt_ave ] ) and ( [ eq : w_alpha_flt_ave ] ) qualitatively show good agreements with the numerical simulations , as shown below .    for positive @xmath68 , @xmath84 is also positive , and is finite at @xmath62 . thus , @xmath85 can be positive for a sufficiently large current . this means , the magnetization can cross over the @xmath61-plane , and move to the south semisphere ( @xmath86 ) . on the other hand , for negative @xmath68 , @xmath84 is also negative . thus , the energy supply by the spin torque is suppressed compared to the case of @xmath69 . then , a relatively large current is required to induce the self - oscillation with a certain oscillation frequency . also , the magnetization can not cross over the @xmath61-plane . ( red ) , @xmath70 ( orange ) , and @xmath71 ( blue ) , respectively . ] we confirmed these expectations by the numerical simulations . figures [ fig : fig3 ] ( a ) , ( b ) and ( c ) show the trajectories of the magnetization dynamics with @xmath69 , @xmath87 , and @xmath88 respectively , while the time evolutions of @xmath73 are shown in fig . [ fig : fig3 ] ( d ) . the current value is 2.0 ma . the current dependences of the oscillation frequency are summarized in fig . [ fig : fig4 ] .    in the case of @xmath89 , the oscillation frequency is low compared to that for @xmath69 because the energy supply by the spin torque is enhanced by the field - like torque , and thus , the magnetization can largely move from the north pole . above @xmath90 ma , the magnetization moves to the south hemisphere ( @xmath91 ) , and stops near @xmath92 $ ] in the south hemisphere , which corresponds to the zero frequency in fig . [ fig : fig4 ] . on the other hand , in the case of @xmath93 , the magnetization stays near the north pole compared to the case of @xmath69 , because the energy supply by the spin torque is suppressed by the field - like torque . the zero frequency in fig . [ fig : fig4 ] indicates the increase of the critical current of the self - oscillation . compared to the case of @xmath69 , the oscillation frequency shifts to the high frequency region because the magnetization stays near the north pole . in conclusion , we derived the theoretical formula of the relation between the current and the oscillation frequency of sto consisting of the perpendicularly magnetized free layer and the in - plane magnetized pinned layer . the derivation is based on the analysis of the energy balance between the work done by the spin torque and the energy dissipation due to the damping . the validity of the analytical solution was confirmed by numerical simulation . we also found that the field - like torque breaks the energy balance , and changes the oscillation frequency . the shift direction of the frequency , high or low , depends on the sign of the field - like torque ( @xmath68 ) . the authors would like to acknowledge t. yorozu , h. maehara , a. emura , m. konoto , a. fukushima , s. yuasa , k. ando , s. okamoto , n. kikuchi , o. kitakami , t. shimatsu , k. kudo , h. suto , t. nagasawa , r. sato , and k. mizushima .                                              in the absence of the applied field , @xmath95 . in this case , @xmath96 given by eq . ( [ eq : i_theta ] ) satisfies @xmath97 . these results mean that , above the critical current , the magnetization immediately moves from the initial state to the @xmath61-plane , and stop its motion . thus , the self - oscillation does not occur .
Abstract:

Answer: the magnetization dynamics of spin torque oscillator ( sto ) consisting of a perpendicularly magnetized free layer and an in - plane magnetized pinned layer was studied by solving the landau - lifshitz - gilbert equation . we derived the analytical formula of the relation between the current and the oscillation frequency of the sto by analyzing the energy balance between the work done by the spin torque and the energy dissipation due to the damping . we also found that the field - like torque breaks the energy balance , and change the oscillation frequency . spintronics , spin torque oscillator , perpendicularly magnetized free layer , the llg equation


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: silicon photomultipliers ( sipms ) are constructed on a single substrate and consist of thousands of microcells operating in the geiger mode . the gain and quantum efficiency of the sipms are comparable to those of the photomultiplier tubes ( pmts ) . they are compact in size and insensitive to magnetic fields . with low - operating voltages , it simplifies the circuit design for electrical safety considerations . therefore , they are a good alternative to pmts @xcite . however , the conventional pmts are available in various sizes that are as large as a few tens of cm in diameter . in contrast , the largest single - chip sipm is less than 1  @xmath0 . for scintillation detectors , the efficiency of scintillation photon collection increases with the area of the photon sensor . moreover , a larger detector volume has a higher sensitivity for radiation detection . because of the small active area of the sipms , it is necessary to use multiple sipms to increase the efficiency of photon collection for a larger scintillator . in this work , scintillation detectors are constructed and tested with a custom 2x2 sipm array . the layout of the sipm array and the geometry of the scintillator were determined by performing geant4 simulations @xcite in which a point source is located 20  cm from the front of the scintillator . the @xmath2-rays are emitted uniformly random in a cone irradiating a circular area inscribed by the front surface . the surface finish of the scintillator is modeled using the `` groundteflonair '' option in the look - up - table ( lut ) @xcite . a thin layer of optical grease , 0.1  mm , is sandwiched between the scintillator and the sipms for transporting scintillation photons . the array has four 6x6 mm@xmath3 sipms arranged in a 2x2 configuration such that the electronics is less complicated . the goal of the simulations was to find the optimum geometry of the scintillator and the layout of the sipms for good @xmath2-ray energy resolution . it is found that the energy resolution of a detector is better for the scintillator with an area comparable to the active area of the sipm array . since the sipms are arranged as a square or rectangle in the array , scintillators with a square or rectangular cross section have a better energy resolution than those with a circular cross section . shown in fig . [ fg : resolgap](a ) is the photon distribution at the exit surface of an 18  mm cubic scintillator . the profile of photon distribution for a 6  mm horizontal strip across the middle of the exit surface is shown in fig . [ fg : resolgap](b ) . as can be seen , more scintillation photons are distributed near the center of the exit surface of the scintillator .    for scintillators with a cross section larger than the active area of the sipms , leaving a gap between the sipms leads to a better energy resolution . in fig . [ fg : resolgap](c ) , the energy resolution for an 18  mm cubic scintillator as a function of the gap size between the active area of the sipms is shown . the best energy resolution is for the gap size between 1 and 2  mm . for larger gap sizes , the resolution gets worse slowly but a shoulder starts to develop on either side of the photopeak resulting in a poorer full - width - at - tenth - maximum . [ fg : spk18gap ] shows the simulated 662  kev @xmath2-ray spectra detected by the 18  mm cubic scintillator and sipm array with the gap size of 0.2 , 1.4 , and 4.0  mm . it can be seen that for the gap size of 4.0  mm a noticeable shoulder appears on the high - energy side of the photopeak and the trough between the photopeak and the compton edge is higher as well . furthermore , the photopeak position shifts to lower channels as the gap size increases . this is due to less photons distributed away from the center of the detection plane as shown in the profile of photon distribution in fig . [ fg : resolgap](b ) . in fig . [ fg : resolgap](d ) , the gap size between the sipms for cubic scintillators with dimensions between 16 and 24  mm to have an optimum resolution is shown . as can be seen , it is necessary to increase the gap size for larger scintillators in order to achieve a better energy resolution . -ray detected by an 18  mm cubic scintillator coupled to a 2x2 sipm array . the gap size between the sipms is ( a ) 0.2 , ( b ) 1.4 , and ( c ) 4.0  mm.,width=312 ] a sipm array with four 6x6  mm@xmath3 sipms @xcite has been constructed following the analysis of the geant4 simulations . three types of scintillators , nai , csi , and cs@xmath4liycl@xmath5:ce@xmath6(clyc ) , have been tested with this sipm array . the scintillator and sipm array were enclosed in a hermetically sealed aluminum container to keep out moisture and ambient light . all the tests were performed in a temperature controlled chamber . the three detectors were tested using a @xmath7cs source with the sipms biased at 27.5  v so that the detectors had the same nominal gain . the output of the preamp was recorded by a waveform digitizer ( struck sis3302 ) operating at 100  mhz . the response of the detectors to the 662  kev @xmath2-ray from @xmath7cs was compared by making a histogram of the integral of the digitized pulses . because the clyc pulse has a long - decay component , an integration time of 20  @xmath8s was used for all the scintillators in order to make an unbiased comparison . consequently , the counting rate was kept low to minimize pulse pileup . as the decay time for nai and csi is shorter , integrating these pulses for such a long time could introduce noise to the integrals . since the comparison was for the detector response , the influence of noise on energy resolution was ignored for the current test . a separate measurement using the proper integration time for each scintillator was carried out to compare the energy resolution . figure  [ fg : lightout ] shows the normalized histogram of 10,000 pulse integrals for the three detectors . the nai detector has the largest output because the sipm response is optimized for the 420  nm scintillation photon . although csi has a larger yield of scintillation photons , the detector response is actually smaller than nai . this is due to the mismatch between the scintillation spectrum of csi and the response of the sipms . lastly , the light output for clyc is almost one half of that for nai . therefore , the pulse integral is the smallest . cs detected by the nai , csi , and clyc detectors . the histograms are obtained by integrating the pulses of the preamp output for 20  @xmath8s.,width=312 ]      the breakdown voltage for the sipms varies individually due to the manufacturing processes . the breakdown voltage also increases with temperature which results in gain variation following temperature changes . in order to optimize the detector resolution , the gain for the sipms was adjusted to be the same within the array . furthermore , the gain for the entire sipm array was stabilized over the temperature range between 20 and 50  @xmath1c @xcite . to check the quality of temperature stabilization , @xmath2-ray sources , @xmath9eu , @xmath10co , @xmath7cs , and @xmath11th , were used . [ fg : det33msrc ] shows the spectra for these sources for the csi detector at 20 , 0 , 20 , and 50  @xmath1c . the centroid of the 662  kev photopeak for @xmath7cs was set to channel 220 for all the temperatures . as can be seen , the location of the photopeaks do not shift for different temperatures because of the stabilization performed for the detector . -rays from @xmath9eu , @xmath10co , @xmath7cs , and @xmath11th sources detected by the csi detector at 20 , 0 , 20 , and 50  @xmath1c.,width=312 ]    to optimize the energy resolution , the preamp output was integrated for 1 , 4 , and 20  @xmath8s for nai , csi , and clyc , respectively . the best resolution obtained was 6.8% for nai , 6.4% for csi , and 7.8% for clyc . it has been reported that an energy resolution of 4% was observed for a 1  @xmath12 clyc scintillator coupled to a pmt . however , when the same clyc scintillator was coupled to sipms the energy resolution was between 6.2% and 8.3% @xcite . the @xmath7cs @xmath2-ray spectrum measured by the csi detector at 20  @xmath1c is shown in fig . [ fg : cs137 ] . the noise is sufficiently low such that the 32  kev peak is clearly visible without contamination . cs @xmath2-rays detected by the csi detector . the gain of the sipms in the array has been matched to achieve an optimum energy resolution of 6.4%.,width=312 ]    the energy resolution of the 662  kev photopeak for csi and clyc measured between 20 and 50  @xmath1c is displayed in fig . [ fg : restemp ] . the variation in energy resolution for the csi detector is small . at high temperatures , the poorer resolution is likely due to noise in the sipms whereas at low temperatures the poorer resolution is attributed to a lower light yield of the scintillator . for the clyc detector , the energy resolution is much worse than that for csi , particularly well below and above room temperature . at 20  @xmath1c , the 662  kev photopeak is barely visible in the @xmath2-ray spectrum . a previous study using a smaller crystal , 1  @xmath13 , obtained mixed results for resolution as a function of temperature @xcite . in that study , sipms from different manufacturers were compared . for the sensl sipm , the variation of energy resolution with temperature is smaller as compared to the present work . for the hamamatsu sipm , resolution as poor as 14% was measured at 10  @xmath1c . since clyc crystals are known to be fragile , it is conceivable that fractures might have developed in the crystal that was used in this work . this makes the use of clyc for portable radiation detectors challenging because these instruments are required to operate between 20 and 50  @xmath1c . in addition to studying the three detectors , a large number of csi detectors have been tested which allow for examining the characteristics of the sipms . fig .  [ fg : distr](a ) shows the distribution of the breakdown voltages for 300 sipms . the standard deviation of the distribution is 0.124  v which is consistent with the manufacturer s specification . [ fg : distr](b ) shows the distribution of the temperature coefficient for the csi detector . the centroid of the distribution is 15  mv/@xmath1c . in contrast , according to the manufacturer , the temperature coefficient for the sipm alone is 21.5  mv/@xmath1c . a @xmath14cf source was used for studying neutron detection by the clyc detector between 20 and 50  @xmath1c . according to the specification , the neutron is expected at 3.2  mev in the @xmath2-ray spectrum @xcite . since @xmath11th has a 2.6  mev @xmath2-peak , it was used along with @xmath14cf to test n-@xmath2 discrimination and as a reference for the energy spectrum . because the largest difference in the pulse shape occurs in the first 2  @xmath8s , a shorter integration time of 4.5  @xmath8s , instead of 20  @xmath8s , was used to speed up the data acquisition . [ fg : clyc1d ] shows the @xmath2-ray spectrum detected by the clyc detector . the photopeak for the @xmath11th 2.6  mev @xmath2-ray appears near channel 2600 and the neutron peak appears near channel 3200 . this agrees with the manufacturer s specification @xcite . as can be seen , the neutron distribution is fairly broad and overlaps with the high - energy tail of the 2.6  mev @xmath2-peak . for this reason , it is necessary to use pulse - shape discrimination ( psd ) to separate neutrons from @xmath2 s . -rays from @xmath11th and neutrons from @xmath14cf detected by the clyc detector.,width=312 ]    the psd is performed by taking the ratio of the pulse integral in the delayed window to the sum of the pulse integral in the prompt and delayed windows , @xmath15 shown in fig . [ fg : cfthpsd20](a ) is the histogram for the psd ratio versus the integral of the preamp pulses . neutrons appears in the group of counts in the upper - right corner . the band of counts located across the center of the histogram are @xmath2-rays . the threshold for detection was set around 1.5  mev to reduce the counting rate . the group of events near the end of the @xmath2 band on the right - hand side is the 2.6  mev @xmath2 peak from @xmath11th . cf and @xmath11th , were present during the measurement . ( b ) histogram of projected psd ratio.,width=312 ]    by projecting the the two - dimensional histogram on to the vertical axis , the histogram of the psd ratio is obtained , as shown in fig . [ fg : cfthpsd20](b ) . the figure - of - merit ( fom ) for psd is defined as the ratio of the difference between the centroid ( @xmath8 ) of the @xmath2 and neutron peaks to the sum of the full - width - at - half - maximum ( @xmath16 ) of the two peaks @xcite , @xmath17 the width for the prompt window and delayed window was varied to search for the best fom . at 20 @xmath1c , the best fom for n-@xmath2 discrimination is 1.9 for the width of the prompt and delayed window of 270 and 850 ns , respectively . the fom increases with decreasing temperature due to the change of pulse shape . at 50 @xmath1c the fom is 1.2 and at 20  @xmath1c the fom is 3.0 , as shown in fig . [ fg : fomtemp ] . although the resolution for clyc is poor at the extreme temperatures , the neutron identification by pulse - shape discrimination works well . the nai and csi scintillation detectors constructed with the custom sipm array have a good energy resolution for the 662  kev @xmath2-ray from @xmath7cs . they are suitable for deploying in a spectroscopic personal radiation detector ( sprd ) such as the flir identifinder r200 . using pulse - shape discrimination , neutrons can be identified without contamination from @xmath2-rays for the clyc detector . however , the resolution for @xmath2-ray detection is poor at extreme temperatures , 20 and 50  @xmath1c . further work is in progress to find solutions for neutron detection in hand - held radiation detection instruments . 99 p.  eckert , h .- c . schultz - coulon , w.  shen , r.  stamen , and a.  tadday , _ nucl . meth . a _ , vol . 620 , no . 2 - 3 , 217 - 226 , apr . 2010 . s.  agostinelli _ et al . _ , _ nucl . instrum . meth . a _ , vol . 3 , pp . 250 - 303 , july 2003 . m. janecek and w.w . moses , _ ieee trans . nucl . 964 - 970 , 2010 . microfc-60035 , sensl technologies ltd . patent pending . rmdinc.com/wp-content/uploads/2014/10/clyc-properties-pdf.pdf k.e . mesick , l.c . stonehill , j.t . morrell , d.d.s . coupland , _ 2015 ieee nss / mic , san diego , ca , 2015 , pp . 1 - 4 . _
Abstract:

Answer: silicon photomultipliers ( sipms ) are a good alternative to photomultiplier tubes ( pmts ) because their gain and quantum efficiency are comparable to pmts . however , the largest single - chip sipm is still less than 1  @xmath0 . in order to use sipms with scintillators that have reasonable sensitivity , it is necessary to use multiple sipms . in this work , scintillation detectors are constructed and tested with a custom 2x2 sipm array . the layout of the sipms and the geometry of the scintillator were determined by performing geant4 simulations . cubic nai , csi , and clyc with 18  mm sides have been tested . the output of the scintillation detectors are stabilized over the temperature range between 20 and 50  @xmath1c by matching the gain of the sipms in the array . the energy resolution for these detectors has been measured as a function of temperature . furthermore , neutron detection for the clyc detector was studied in the same temperature range . using pulse - shape discrimination , neutrons can be cleanly identified without contribution from @xmath2-photons . as a result , these detectors are suitable for deploying in spectroscopic personal radiation detectors ( sprd ) . geant4 , silicon photomultiplier , nai , csi , clyc , pulse - shape discrimination , gain stabilization , temperature , spectroscopic personal radiation detector .