Instruction: Write an abstract for the following acadamic paper.

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

Paper: the yields of the heavy residues @xcite , i.e. , the large remnants of the heavy member of an asymmetric reacting pair of nuclei , are known to comprise a large fraction of the reaction cross section for intermediate energy nuclear collisions . however , there are experimental difficulties in studying the formation of these nuclei in asymmetric reactions using `` normal '' kinematics . the energies of the residues are low ( @xmath415 kev / nucleon @xcite ) and the masses are large . experimental thresholds cause one to miss substantial portions ( @xmath550% ) of the product distributions @xcite . one can overcome the problem of detecting low fragment energies by studying these asymmetric reactions in `` inverse kinematics '' , i.e. , bombard a low mass nucleus with a large mass nucleus . by using a high resolution spectrometer / detector system to observe the projectile - like fragments ( plfs ) from such reactions , one can gain important information about the reaction mechanism , complementary to that obtained in exclusive studies , where only the light reaction partners are observed with high isotopic resolution . pioneering studies by bazin , et al.@xcite , faure - ramstein , et al . @xcite , pfaff , et al . @xcite and hanold , et al . @xcite have shown the utility of this approach in studying the heavy reaction products from kr+x and xe+x mass asymmetric collisions at intermediate energies . in a previous paper @xcite , we reported the results of a similar study of the reaction of 20 mev / nucleon @xmath0au with the light nuclei @xmath6c and @xmath7al .  in these cases , we focussed on _ fusion - like events _ and tried to understand the details of the particle emission / fission competition in the de - excitation of these hot nuclei .  in this paper , we report the results of an extension of the previous work to study the yields and velocities of the heavy residues and fission fragments from the interaction of @xmath0au ( 20 mev / nucleon ) with an intermediate mass nucleus , @xmath1zr .  in these systems , fusion - like collisions are less frequent . previous studies @xcite , @xcite of the collisions of heavier nuclei at intermediate energies have shown binary dissipative collisions to be the dominant reaction mechanism . evidence has been presented @xcite for a sequential decay of one of the initial binary fragments leading to a three ( or more ) body final state along with dynamically emitted imfs from the collision complex @xcite .  in some of the binary encounters the projectile - like fragments have been found @xcite to have very high temperatures ( t@xmath87 mev ) . a variety of phenomenological models have been used , with modest success , to describe these collisions @xcite . the goal of this extension of previous work is to understand quantitatively the reaction mechanism(s ) responsible for residue formation in more symmetric collisions at 20 mev / nucleon . we shall present evidence indicating that , in addition to the traditional evaporation residues from peripheral collisions , some residues are formed in `` hard '' collisions , possibly resulting from very asymmetric fission and/or imf emission . the distributions of these residues , which are very proton - rich , can not be fully described by available phenomenological models of these collisions . although the au + zr reaction has not , to the best of our knowledge , been studied before , there are a large number of studies @xcite , @xcite of the analogous @xmath0au + kr reaction at intermediate energies ranging from 4.9 to 70 mev / nucleon with some studies @xcite at a beam energy of @xmath420 mev / nucleon . what this work adds to this extensive data set is detailed , high resolution information on the properties of the heavy , au - like fragments in these reactions garnered from physical rather than radiochemical measurements . ( the previous studies of the kr + au reaction generally point to the presence of both deep inelastic and fusion - like mechanisms . the fragments of au arise as both fission fragments and evaporation residues @xcite with the residues being formed in collisions involving a larger dissipation of the projectile energy . ) the paper is organized as follows : in section 2 , a brief description of the experimental apparatus , the measurements and the data analysis is given . in section 3 , yield distributions are presented and compared with previous work . also the velocity distributions are discussed . in section 4 , the results of the measurements are compared to modern models of the dissipative features of intermediate energy collisions . finally , conclusions from the present study are summarized in section 5 . the experiment was performed at the national superconducting cyclotron laboratory at michigan state university using the a1200 fragment separator @xcite . this work is a direct extension of the work of souliotis et al . @xcite , where the collisions of 20 mev / nucleon @xmath0au with c and al were studied . it should be noted that use of the a1200 for very asymmetric collisions with heavy beams ( where the grazing angle is small , i.e. 1 - 2@xmath9 @xcite ) , leads to rather efficient collection of the projectile residues ( that can reach 50% ) at this energy regime . however , using the a1200 for more symmetric systems , as for the au+zr reaction in this work ( where the grazing angle is 5.5@xmath9 ) , results in substantial loss of residue cross sections , since the measurements are performed near 0@xmath9 , while the residue cross sections are expected to peak at larger angles , close to the grazing angle . in this case , only a small fraction of the residue cross section can be collected ( of the order of a few percent ) . despite this limitation , we attempted to exploit the possibility of getting isotopically resolved au fragments with this apparatus and tried to extract some information regarding the production mechanism . the experimental setup , calibration procedures and data analysis were similar as in ref . @xcite and are briefly summarized below .    a 20 mev / nucleon @xmath0au beam , produced by the k1200 cyclotron , interacted with a @xmath1zr target of thickness 1.0 mg/@xmath10 . ( the target thickness was such that the maximum energy loss of the beam in traversing the target was 0.26 mev / nucleon @xcite . ) the reaction products were analyzed with the a1200 fragment separator operated in the medium acceptance mode [ with an angular acceptance of 0.8 msr ( @xmath11 mr , @xmath12 mr ) and a momentum acceptance of 3% ] . the primary beam struck the target at an angle of 1.0@xmath9 relative to the optical axis of the spectrometer . the a1200 provided two intermediate dispersive images and a final achromatic image ( focal plane ) . at the focal plane , the fragments were collected in a three - element ( @xmath13e@xmath14 , @xmath13e@xmath15 , e ) si surface barrier detector telescope . the 300 mm@xmath16 si detectors were 50 , 50 and 300 @xmath17 m thick , respectively . time of flight was measured between a parallel plate avalanche counter ( ppac ) and a microchannel plate detector ( mcp ) positioned at the first dispersive image and at the focal plane , respectively , and separated by a distance of 14 m. the ppac at the first dispersive image was also x  y position sensitive and used to record the position of the reaction products . the horizontal position , along with nmr measurements of the a1200 dipole fields , was used to determine the magnetic rigidity @xmath18 of the particles . a series of measurements at overlapping magnetic ridigity settings of the spectrometer were performed in the region 1.5001.750 tesla - meters . this region is below the rididity b@xmath19=1.824 tm of elastically scattered 20 mev / nucleon @xmath0au particles at their equilibrium charge state ( q@xmath20=69.5 ) @xcite . the reaction products were characterized by an event - by - event measurement of de / dx , e , time of flight , and magnetic rigidity . the response of the spectrometer / detector system to ions of known atomic number z , mass number a , ionic charge q and velocity was measured using analog beams as described previously @xcite . the procedure to extract z , q and a from the measured quantities is outlined below and described in more detail in @xcite . the transmission detectors ( @xmath21 and @xmath22 ) were calibrated by directly correlating pulse height with energy loss as calculated using the data of hubert et al . for the stopping ( e ) detector , pulse - height defect ( phd ) corrections were necessary to obtain an accurate reconstruction of the residual energy of the particles . the procedure developed in @xcite was used . the total energy e@xmath23 of the particles entering the spectrometer ( i.e. the sum of the energies deposited in the si detectors plus small corrections due to the presence of the ppacs and mcp ) had an overall resolution ( fwhm ) of 1.1% . the velocity resolution was about 0.5% . the determination of the atomic number z was based on the energy loss of the particles in the first @xmath24 detector and their velocity and was reconstructed using the expression :    @xmath25    where @xmath26 is the velocity of the ion entering the detector and @xmath27e the energy loss . in order to determine the functions @xmath28 , @xmath29 and @xmath30 in the velocity range of interest , we used the data of hubert et al . @xcite  to obtain the coefficients of eq . ( 1 ) for the z range 2585 and in the energy range 1224 mev / nucleon by applying a least - squares fitting procedure at each energy , in steps of 0.5 mev / nucleon . the values of each coefficient at the various energies were then fitted with polynomial functions of velocity . the atomic number z of the particles was reconstructed from the measured @xmath13e and @xmath31 using eq . ( 1 ) with a resulting resolution ( fwhm ) of 0.9 z units for heavy residues ( a=160200 ) and 0.6 z units for fission - like ( a=100140 ) residues . the ionic charge @xmath32 of the particles entering the a1200 was obtained from the total energy e@xmath23 , the velocity and the magnetic rigidity according to the expression : @xmath33 where e@xmath23 is in mev , b@xmath34 in tm , @xmath35 and @xmath36 . the measurement of the ionic charge q had a resolution of 0.7 and 0.5 q units for heavy and fission - like residues , respectively . since the ionic charge must be an integer , we assigned integer values of q for each event by setting appropriate windows on each peak of the q spectrum ( see below ) . using the magnetic rigidity and velocity measurement , the mass - to - charge a / q ratio of each ion was obtained from the expression : @xmath37 now , combining the q determination with the a / q measurement , the mass a was obtained as : @xmath38 ( q@xmath39 is the integer ionic charge determined as above ) with an overall resolution ( fwhm ) of 1.0 and 0.7 a units , again for heavy and fission - like residues respectively . the requirement that q be within @xmath400.3 units of its coresponding integer value was applied to eliminate contributions of adjacent q values to the mass assignement . similarly , the requirement that z be within @xmath400.3 units of its coresponding integer value was applied to avoid contributions of adjacent z values to the distributions of a given element . combination and appropriate normalization of the data at the various magnetic rigidity settings of the spectrometer provided fragment distributions with respect to z , a , q and velocity . correction of missing yields caused by charge changing at the ppac ( positioned at the dispersive image ) was performed as in @xcite , based on the equilibrium charge state prescriptions of baron et . al . @xcite . the distributions were then summed over all values of q. it should be pointed out , as mentioned previously , that the resulting distributions in z , a and velocity are residue yield distributions at a reaction angle of 1@xmath9 with an angular window of 1.1@xmath9 , in the magnetic rigidity range 1.5001.750 tm . the extracted residue yield distributions , sorted by mass number a , atomic number z , and velocity @xmath41 ( @xmath42@xmath43/c ) , are shown in figs . 1a , 1b and 1c , respectively . the yields have been normalized to beam current and target thickness and are given in mb in fig . the a distributions are summed over all values of z and @xmath44 . the z  a distributions are summed over all values of @xmath41 , and finally , the @xmath44a distributions are summed over all values of z.    it is apparent from examining figs . 1a and 1b that the spectrometer acceptance and resolution have distorted the fragment yield distributions . for example , it was not possible to cleanly separate the beam nuclei from quasielastic fragments . therefore , products that were very close in mass number , atomic number and energy to the projectile were either not collected by the spectrometer or rejected in the analysis . as a result , a small subset of quasielastic events is displayed in fig . 1 . ( the quasielastic events shown in fig . 1 are very proton rich relative to the projectile nucleus , which is incorrect . their yields are greatly reduced compared to other fragments which is also incorrect . ) we have chosen therefore to focus on the non - quasielastic events ( a @xmath45 160 ) in the analysis . au with @xmath1zr . ( a ) isobaric yield distribution . the data are shown as solid points , whereas the solid and dashed lines indicate the yield curves from radiochemical measurements of the reactions 35 mev / nucleon @xmath46kr+@xmath0au @xcite and 21 mev / nucleon @xmath47xe+@xmath0au @xcite , respectively . ( b ) yield distributions as a function of z ( relative to the line of @xmath48 stability , z@xmath49 ) and a. highest yield contours are plotted with thicker lines . successive contours correspond to a decrease of the yield by a factor of 2 . the dashed line indicates the expected values for intermediate energy projectile fragmentation @xcite , the dot - dashed line indicates the position of the modified `` evaporation residue attractor line '' @xcite and the solid line indicates the position of the proton dripline for odd z nuclei @xcite . finally , the thick line segments show the z@xmath50 values obtained in the radiochemical measurements of 21 mev / nucleon @xmath47xe+@xmath0au @xcite . ( c ) velocity vs. mass distributions . the definition of the velocity lines is as follows : horizontal full line : beam velocity ; horizontal dashed line : velocity corresponding to 15% linear momentum transfer ( lmt 15% ) . , width=343,height=623 ]    the mass yield curve ( fig . 1a ) , has two humps , the distorted quasielastic one peaking at a=193 and the other extending from a=170 down to a=80 . the second bump in the mass yield distribution consists of two groups of fragments ( as can be deduced with the aid of fig . 1c ) ; one group resulting from modest ( @xmath415% ) momentum transfer events with a=120160 ( heavy residues ) and another group with a=90120 and velocities of forward - moving fission fragments .    the radiochemical mass yield curves of target - like residues ( and fission fragments ) from the reactions 35 mev / nucleon @xmath46kr+@xmath0au @xcite and 21 mev / nucleon @xmath47xe+@xmath0au @xcite are also shown for comparison ( solid and dashed lines , respectively ) . as noted earlier , due to the spectrometer acceptance , the measured yields , are much lower than the radiochemical cross sections and bear only scant resemblance to them .    in fig . 1b , the yield distributions of the reaction products are shown as contour plots of z versus a. the fragment z is given relative to the line of beta - stability , z@xmath51 , which was taken as z@xmath51 = a/(1.98 + 0.0155a@xmath52)@xcite . the observed fragment distributions lie on the proton - rich side of stability ( z  z@xmath53 . shown in fig . 1b , are the centroids of the isobaric charge distributions from the radiochemical measurements of target - like fragments from the reaction 21 mev / nucleon @xmath47xe+@xmath0au @xcite . the radiochemical centroids for near - au residues are close to stability due to the dominance of quasielastic events in this region . in the symmetric fission fragment region a=90120 , the centroids are in fair agreement with the present data , as they both correspond to fission fragments . in the region a@xmath4150 , the discrepancy may be due to a combination of experimental factors ( e.g. limited number of measured nuclides in the radiochemical data and acceptance cut in the present data ) .    also shown in fig . 1b are the centroids of the z distributions seen in intermediate energy projectile fragmentation @xcite . despite the fact that no fission is involved in these reactions , there is some similarity ( but differences in the slope of z vs a ) between this prescription and the centroids of the yield distribution of the present data . charity @xcite has recently discussed the concept of an `` evaporation residue attractor line''(eal ) . the eal is a line in the n  z plane that represents the locus of fragment yields produced by the evaporation of neutrons and charged particles from highly excited nuclei . ( the statistical model used to predict the position of the eal is gemini ) .  in fig . 1b , we show the position of the modified heavy fragment eal for the decay of highly excited heavy nuclei where the primary fragment is near the valley of beta stability .  as shown in this figure , the centroids of the observed yields are similar to the modified eal for this reaction ( a=80160 ) , but the slope of z vs a is not the same . this similarity in charge distributions for the products from the present reaction and the eal for a = 120160 indicates that these fragments are decay products of highly excited nuclei .    what causes the centroids of the contours in the data to become more n - rich compared to these two prescriptions as the fragment mass decreases ? a possible answer is the occurrence of fission which produces more n - rich fragments than fragmentation ( residue production ) . ( the two empirical prescriptions @xcite do not include the effects of fission . ) the argument is that at a=100 , a larger fraction of the fragments arise from fission ( fig . 1c ) than at a=150 , thus producing the observed slope of z vs a. evidence for this comes from comparison with the z distributions measured radiochemically for the reaction of 20 mev / nucleon xe with au , where the fission of the au - like nucleus is more prominent ( fig . the observed centroids of the z distributions for this reaction are more n - rich than the fragments observed in this work . the definitive data comes from comparing fig . 1b and 1c and noting that as fission ( high velocity events ) becomes more prominent , the charge distribution becomes more n - rich . one does note that the contours of the observed yields in fig . 1b lie close to the position of the proton dripline ( plotted for odd z nuclei @xcite ) . some of the observed products have 57 more protons than the stable isobar . to examine this point further , we show the detailed fragment charge distributions for selected a values for non - quasielastic fragments in fig . au with @xmath1zr . the data are shown as solid points . the full - line curves are calculations using the dit+icf code coupled to gemini ( see text ) . , width=423,height=634 ]    for fig . 2 , an attempt has been made to correct the measured yields for the spectrometer angular acceptance , so that we may obtain estimates of fragment production cross sections , despite the fact that only a small fraction of the cross section was collected in the present measurements . correction factors for the angular distributions of the fission - like fragments ( a=90160 ) were obtained assuming isotropic binary decay of the primary reaction product . for simplicity , this product was assumed to be a au nucleus with velocity as this observed for the heavy residues ( which , as already mentioned , correspond to @xmath415% linear momentum transfer ) and was assumed to scatter off the target at an angle somewhat inside the grazing angle . ( the grazing angle is 5.5@xmath54 degrees and the scatter angle was taken to be 4.0@xmath54 . ) typical correction factors were 40 , 60 , 90 , 150 for masses a=160 , 140 , 120 , 100 , respectively . after applying mass - dependent correction factors to the measured yields , the resulting yields were normalized to values of cross sections of non - isomeric independent - yield nuclides from the radiochemical data of the reaction ( 21 mev / nucleon ) @xmath47xe+@xmath0au @xcite . fragment production cross sections resulted with an estimated uncertainty of a factor of 2 ( due to the correction procedure ) and , for several a values , are presented in fig . 2 .    in examining fig . 2 , we are struck by the large number of nuclide yields measured for each mass number a , when compared to typical radiochemical measurements where , at best , the yields of two or three isobars are measured for each a value . the measured distributions appear to be rather broad . the general proton - rich character of these yield distributions indicates that these intermediate energy collisions may be useful for the synthesis of heavy proton - rich nuclei , complementary to the usual practice of synthesizing these nuclei via symmetric compound nuclear reactions involving the pxn exit channel @xcite . an attempt to produce and identify proton emitting nuclei using the au+zr reaction is described in @xcite .    the velocity vs mass distributions ( fig . 1c ) show the general absence of fusion - like collisions and the prevalence of dissipative collisions . as already mentioned , one also sees evidence for the occurrence of fission events for a=90120 with velocities higher than that of the beam , corresponding to fission fragments moving forward in the rest frame of the fissioning nucleus ( the spectrometer settings were such that only forward - moving fragments were collected from the two kinematical solutions ) . in addition , one sees a group of events ( a=120160 ) with rather constant velocity , lower than that of the beam ( by approx . 5% ) corresponding to @xmath4 15% fractional linear momentum transfer ( flmt ) . the existence of a constant fragment velocity , roughly independent of fragment mass , is not consistent with fission , where the fragment velocity should be inversely correlated with the fragment mass ( as it is for a=90120 for the au - like fragment ) @xcite . we associate this group of fragments ( a=120160 ) mainly with residues that have emitted one or more imfs in accord with a recent radiochemical study of target - like residues from the @xmath55pb ( 29 mev / nucleon ) + @xmath0au reaction @xcite . from the velocity of this group of fragments ( relative to that of the beam ) we can estimate an average total kinetic energy loss of @xmath4350 mev . we also note in fig . 1c , that the fragment group a=90120 , along with the fission fragments , contains residues characterized by constant velocity ( as in the a=120160 group ) . it is rather remarkable that residues with less than half the mass of the projectile nucleus are observed in this reaction . to gain insight into the underlying reaction mechanisms , we have performed simulations using an appropriate phenomenological model . the basic features of the reaction mechanism model are pre - equilibrium emission , dissipative deep inelastic transfer ( dit ) and for the most central collisions , an incomplete fusion component ( icf ) . the hybrid model we used that included these features has been described by veselsky @xcite .    for each event corresponding to a given impact parameter ( partial wave ) , the effect of pre - equilibrium emission was calculated using a variant of the exciton model . then the interaction of the projectile and target nuclei are simulated using the tassan - got / stephan model @xcite for deep inelastic transfer assuming stochastic nucleon exchange ( in the region of orbital angular momentum @xmath56 = 300870 ) . for trajectories where the overlap between projectile and target nuclei exceeds 3 fm , it is assumed that incomplete fusion takes place with this mechanism being modeled using a variant of the abrasion - ablation model @xcite . following the creation of the primary fragments by this hybrid mechanism , the statistical de - excitation of the excited primary fragments was simulated using gemini @xcite . the results of this calculation were then filtered by the angular and momentum acceptance of the spectrometer . this model has been used successfully @xcite to describe the residues formed in the 35 mev / nucleon kr + au reaction @xcite and the 20 mev / nucleon au + ti reaction studied by us @xcite . simulations were run for partial waves from @xmath56=0870 . de - excitation of the primary fragment distribution was done using gemini . this statistical deexcitation code uses monte carlo techniques and the hauser - feshbach formalism to calculate the probabilities for fragment emission with z@xmath452 . heavier fragment and fission fragment probabilities are calculated using the transition state formalism of moretto @xcite . in the gemini calculations , we have used lestone s temperature dependent level density parameter @xcite , a fading of shell corrections with excitation energy and we enabled imf emission . finally , fission delay was enabled with parameters as in @xcite . each partial - wave distribution was appropriately weighted and combined to give the overall fragment a , z and velocity distributions ( figs . 3 and 4 ) . au with @xmath1zr . ( a ) isobaric yield distributions . the data are shown as solid points . the upper full line is the result of dit+icf / gemini calculation and the dashed line is from dit / gemini ( see text ) . the lower full line and dashed lines are the results of the same calculations as the upper ones , but with a cut corresponding to the angular and momentum acceptance of the spectrometer . ( b ) charge distributions . the data are shown as contours ; the calculated values from dit+icf / gemini are shown as i ) thick full line : with acceptance cut and , ii ) thin full line : without acceptance cut , and the ones from dit / gemini are shown as a dashed line ( without acceptance cut ) . , width=343,height=423 ]        au with @xmath1zr . , width=317,height=476 ]    au with @xmath1zr . , width=343,height=264 ]    in fig . 3a , we compare the measured and calculated mass yield curves for a number of conditions . the calculations were performed neglecting small impact parameters ( deep inelastic transfer only , dit / gemini ) or including all impact parameters ( dit+icf / gemini ) . the upper dashed line is from dit / gemini calculation and the upper solid line is from dit+icf / gemini calculation ( without filtering with the spectrometer acceptance ) . the corresponding lower curves are from the same calculations in which filtering with the spectrometer acceptance has been peroformed . as we see , the shape of the mass yield curve is significantly altered by the overall spectrometer cut . note however , the similarity of shapes between the measured yield curve and the calculated yield curves filtered with the spectrometer acceptance .    in fig . 3b , the measured z vs a yield distribution is compared with the predicted positions of z@xmath57 , the most probable fragment atomic number for a given a value ( given with respect to the line of @xmath44 stability ) . displayed are the results of dit / gemini ( dashed line , no acceptance filter ) , and of dit+icf / gemini ( no acceptance filter : thin solid line , with spectrometer acceptance filter : thick line ) . we observe that the latter calculation with the acceptance cut does a fair job in describing the z@xmath50 values of the a=160190 residues , but predicts more n - rich fragments , consistent with fission , in the region a@xmath58160 . thus , the heavy residue group a=120160 seen in the data is not predicted by the calculation , due to the prevalence of fission as a deexcitation channel in the gemini calculation . the calculated fragment charge distributions are compared with the data in fig . the calculated centroids of the distributions are too n - rich as noted above , but the predicted widths of the charge distributions agree rather well with the observations .    in fig . 4 , we show the predicted velocity distributions ( which are to be compared to fig . 1c ) . in fig . 4a , we show the unfiltered velocity distribution , while in fig . 4b we show the velocity distribution as filtered by the spectrometer angular and momentum acceptance . qualitatively , some features of the measured distributions are reproduced , i.e. , the appearance of residue events ( a = 160190 ) with mean velocities similar to the observed velocities and the characteristic forward moving fission fragment velocity pattern for the lighter fragments ( a@xmath58120 ) which are also similar to the observations . however , as already pointed out , the predicted velocities of the a=120160 group are inversely correlated with mass and ( along with their average z values ) are consistent with fission , whereas the measured properties of this fragment group are those of residues . we conclude that the hybrid model proposed by veselsky , coupled to a modern statistical de - excitation code , can describe some gross features of the residue distributions , especially those not too far from the projectile nucleus , whereas for lower mass residues it predicts mostly fission . we might then ask what we might learn about residue formation from the simulation model only . to gain insight regarding the primary interaction stage , we show in fig . 5 the simulated properties ( dit+icf / gemini ) of the collision as a function of impact parameter b. very large excitation energies and angular momenta are involved in these collisions . interestingly , the average spin of products in the impact parameter range 68.5 fm is predicted to be higher than the maximum angular momentum for which the symmetric fission barrier vanishes ( approx . 80@xmath59 ) which causes these nuclei to fission .    in fig . 6 , we show the simulated fragment mass distributions arising from collisions with various impact parameters . from these impact parameter sorted mass distributions , we see that the near projectile fragments originate in very peripheral collisions ( quasielastic events ) with ( fig . 5 ) low associated excitation energies and angular momenta . the rest of the fragments are produced by collisions involving a large range of impact parameters ( and thus the fragment mass number a can not be used as an impact parameter trigger ) . these broad mass distributions represent events where the heavy fragment from the collision either fissioned symmetrically ( a=90120 ) or fissioned very asymmetrically or emitted one or more imfs in its de - excitation ( a=120160 ) . as already mentioned , the measured velocity spectra ( fig . 1c ) indicate that these fragments ( a=120160 ) are not the result of symmetric fission . what have we learned from this high resolution experimental study of the residues produced in the reaction 20 mev / nucleon @xmath0au + @xmath1zr ? two findings seem of particular importance . they are : ( a ) the observation of surviving residues ( a=120160 ) which appear to result from `` hard collisions '' . they are relatively neutron - deficient with broad charge distributions and velocity , on average , corresponding to a linear momentum transfer of @xmath4 15% . ( b ) a semi - quantitative description of heavy residue formation in intermediate energy collisions between massive nuclei by the hybrid model of veselsky followed by the statistical deexcitation code gemini was found . however , the calculation was not able to describe the residue group with a=120160 , and gave primarily fission fragments in the whole mass range a@xmath58160 . we gratefully acknowledge the support of the a1200 group and operations staff of michigan state university during the measurements and the help of l. hart in the analysis of the data . one of us ( g.a.s ) is also thankful to m. veselsky for his help in the course of calculations and stimulating discussions . financial support for this work was given , in part , by the u.s . department of energy under grant no . de - fg06 - 88er40402 , no . de - fg03 - 97er41026 , and contract de - ac03 - 76sf00098 and the national science foundation under grant no . phy-95 - 28844 . 00 k. aleklett , m. johansson , l. sihver , w. loveland , h. groening , p.l . mcgaughey , and g.t . seaborg , nucl . a149 ( 1989 ) 591 . w. loveland , k. aleklett , l. sihver , z. xu , c. casey , d.j . morrissey , j.o . liljenzin , m. de saint - simon and g.t . seaborg , phys . c41 ( 1990 ) 973 . d. bazin , d. guerreau , r. anne , d. guillemand - mueller , a.c . mueller , and m.g . saint - laurent , nucl . a515 ( 1990 ) 349 . b. faure - ramstein , f. auger , j.p . wieleczko , w. mittig , a. cunsolo , a. foti , e. plagnol , j. quebert , and j.m . pascaud , nucl . a586 ( 1995 ) 533 . r. pfaff , d.j . morrissey , w. benenson , m. fauerbach , m. hellstrom , c.f . powell , b.m . sherrill , m. steiner , and j.a . winger , phys . c53 ( 1996 ) 1753 . hanold , et al . , phys c52 ( 1995 ) 1462 . souliotis , k. hanold , w. loveland , i. lhenry , d.j . morrissey , a.c . veeck , and g.j . wozniak , phys . c57 ( 1998 ) 3129 . o. granier , et al . a481 ( 1988 ) 109 . h. barz , j.p . bondorf , c.h . dasso , r. donangelo , g. pollarolo , h. schulz , and k. sneppen , phys . c46 ( 1992 ) r42 . a. yokoyama , w. loveland , j.o . liljenzin , k. aleklett , d.j . morrissey , and g.t . seaborg , phys . c46 ( 1992 ) 647 . marchetti , et al . c48 ( 1993 ) 266 . m. aboufirassi , et al . , lpcc 93 - 14 , september , 1993 ; j.f . lecolley , et al . b325 ( 1994 ) 317 . garcia - solis , a.c . mignerey , h. madani , a.a . marchetti , d.e . russ , and d. shapira , phys . c52 ( 1995 ) 3114 . baldwin , et al . 74 ( 1995 ) 1299 . m. morjean , et al . phys . a591 ( 1995 ) 371 . denterria , f. fernandez , e. laguera , m. debeauvais , j. ralarosy , s. jokic , m. zamani , and j. adloff , phys . c52 ( 1995 ) 3179 . j. toke , et al . a583 ( 1995 ) 519 . otto , g.t . seaborg , and m.m . fowler , phys . c17 ( 1978 ) 1071 . m. rajagopalan , l. kowalski , d. logan , m. japlan , j.m . alexander , m.s . zisman , and j.m . miller , phys . c19 ( 1979 ) 54 . d. dalili , et al . , z. phys . a316 ( 1984 ) 371 . d. dalili , et al . , a320 ( 1985 ) 349 . g. rudolf , et al . 57 ( 1986 ) 2905 . d. dalili , et al . a454 ( 1986 ) 163 . r. lucas , et al . a464 ( 1987 ) 172 . h. solbach , h. freiesleben , w.f.w . schneider , d. schull , p. braun - munzinger , b. kohlmeyer , m. marinescu , and f. puhlhofer , nucl . a467 ( 1987 ) 349 . w. loveland , k. aleklett , l. sihver , z. xu , c. casey , and g.t . seaborg , nucl . a471 ( 1987 ) 175c . k. aleklett , w. loveland , m. de saint - simon , l. sihver , j.o . liljenzin , and g.t . seaborg , phys . lett b236 ( 1990 ) 404 . m. zamani , m. debeauvais , j. ralarosy , j.c . adloff , f. fernandez , s. jokic , and d. sampsonides , phys . c42 ( 1990 ) 331 . a. adorno , a. bonasera , m. avinato , m. colonna , a. cunsolo , g.c . dileo , m. ditoro , and f. gulminelli , nucl . a529 ( 1991 ) 565 . l. stuttge , et al . a539 ( 1992 ) 511 . b. lott , et al . , z. phys . a346 ( 1993 ) 201 . peaslee , et al . c49 ( 1994 ) r2271 . w. skulski . , et al . , in advances in nuclear dynamics , w. bauer , a. mignerey , eds . , ( plenum , new york , 1996 ) pp 2130 . j. colin , g. bizard , a. genoux - lubain , c. lebrun , j.f . lecolley , m. louvel , ch . meslin , g. rudolf , and l. stuttge , nucl . a583 ( 1995 ) 449 . b. djerroud et al . , phys . c64 ( 2001 ) 034603 ; w. skulski , et al . , c53 ( 1996 ) r2594 . w.q . shen , et al . c57 ( 1998 ) 1508 . p. moller , j.r . nix , and k.l . kratz , at . data and nucl . data tables 66 ( 1997 ) 131 . c. n. davids , et al . ( 1996 ) 592 . souliotis , physica scripta t88 ( 2000 ) 153 . d. dalili , et al . , nucl phys . a454 ( 1986 ) 163 .
Abstract:

Answer: the yields and velocity distributions of heavy residues and fission fragments from the reaction of 20 mev / nucleon @xmath0au + @xmath1zr have been measured using the msu a1200 fragment separator . a bimodal distribution of residues is observed , with one group , resulting from peripheral collisions , having fragment mass numbers a=160200 while the other group , resulting from hard"collisions , has a=120160 . this latter group of residues can be distinguished from fission fragments by their lower velocities . a model combining deep - inelastic transfer and incomplete fusion for the primary interaction stage and a statistical evaporation code for the deexcitation stage has been used to describe the properties of the product distributions . @xmath2 ,    ,    @xmath3 ,    , and    .    nuclear reactions , incomplete fusion , deep inelastic scattering 25.70.-z , 25.70.jj , 25.70.lm


Instruction: Write an abstract for the following acadamic paper.

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

Paper: the emerging field of spintronics @xcite , also known as magneto- or spin - electronics , has provided a strong motivation towards the fabrication of novel magnetic materials . a strong interest has been focused on the so - called half - metallic ferromagnets , _ i.e. _ metallic compounds for which either the majority or minority spin - channel presents a gap at the fermi level @xmath0 @xcite . such behavior is detected in various perovskite structures @xcite , some oxides like cro@xmath1 @xcite and magnetite ( fe@xmath2o@xmath3 ) @xcite , several heusler alloys @xcite , or some diluted magnetic semiconductors ( dms ) @xcite .    in 2000 , akinaga and collaborators succeeded to grow a thin film of the binary alloy cras on top of gaas @xcite . it was found that cras was actually continuing the zinc - blende ( zb ) structure of gaas crystallizing in a metastable phase . the most astonishing property was the total spin magnetic moment of 3 @xmath4 , which is consistent with the half - metallic character predicted by relevant electronic structure calculations published in the same paper and reproduced with larger accuracy in latter publications @xcite . the thin films of cras have also the great advantage of a high curie temperature @xmath5 , exceeding the 400  k ( limit of the magnetometer used in the experiment ) @xcite . this discovery intrigued the interest of experimentalists and several attempts were made to grow similar metallic compounds as thin films . et al . _ have grown zb crsb films by molecular beam epitaxy and found a @xmath5 similar to the cras case @xcite . moreover , the growth of small zb mnas islands on gaas substrates @xcite has been achieved and cras / gaas multilayers have been fabricated @xcite showing that coherent heterostructures are possible . these half - metallic materials combine high values of @xmath5 with the coherent growth on top of semiconductors and as the know - how on the growth increases , they can be considered as very promising candidates for spintronics applications . several theoretical studies have been devoted to the zinc - blende compounds of transition elements with group - v and vi elements ( respectively , pnictides and chalcogenides ) . the calculations of the total energy and structural optimization show , in agreement with experiment , that the ground state crystal structure of these compounds is not of the zinc - blende type . in most cases the equilibrium crystal structure is of the nias type @xcite . this structure does not lead to half - metallicity @xcite . the zinc - blende structure can only be obtained by the epitaxial growth of few layers on semiconductors . galanakis and mavropoulos studied in detail the magnetic properties and the position of the fermi level in zb v- , cr- and mn - pnictides and chalcogenides as a function of the lattice constant @xcite . several authors have confirmed these results @xcite . some of the binary semiconductors crystallize in the wurtzite structure and a number of studies is devoted to such systems @xcite . of special interest is the case of half - metallic zb mnc compound whose properties differ from the properties of other mn compounds @xcite . also the cases of unsupported clusters @xcite and nanocrystallites @xcite have been studied . an important aspect of the studies is the stability of the half - metallic gap . it has been shown that the tetragonal deformation of the zb lattice does not destroy the gap @xcite and for some systems it even makes the zb structure more stable @xcite . contrary to other half - metallic systems , in the case of the zb transition metal compounds it is possible to retain the half - metallicity both at the surfaces of these materials @xcite and in the case of multilayers with binary semiconductors @xcite . the v and cr atoms have the ability to increase considerably their spin moment as a result of decreasing number of neighbors keeping in this way the half - metallicity of the system . contrary to cr and v , the mn spin moment is saturated in these compounds that leads to the loss of the half - metallicity at the surface @xcite . the spin - orbit coupling has very weak effect on the spin - polarization of the electron states at the fermi level for the @xmath6 atoms up to 5th period of the mendeleev s table @xcite . in the case of mnbi , where bi is a heavy element from the 6th period , the spin - orbit coupling leads to the decrease of the spin - polarization from 100% to 77% @xcite . recently , chioncel and collaborators studied the influence of the correlation effects on the electron structure of cras @xcite . they found that the spin - magnon interaction leads to the appearance of nonquasiparticle states in the spin - minority channel . the states are shown to lie above the fermi level and to be sensitive to the lattice constant . in contrast to the gap properties , the exchange interactions in these compounds attracted much less attention . shirai has shown that in the case of pnictides ( vas , cras and mnas ) the ferromagnetic state is energetically preferable compared to the antiferromagnetic state @xcite . only in the case of zb feas an opposite trend was found @xcite . sakuma demonstrated that the mn compounds isoelectronic to cras ( mnsi , mnge and mnsn ) are half - metallic ferromagnets @xcite . in these compounds , similar to cras , the long - range ferromagnetism is stabilized by the short - range direct mn - mn interaction and an indirect interaction through the @xmath6 atoms . the curie temperature in both cras and the isoelectronic mnz compounds was found to exceed 1000  k @xcite . kbler calculated @xmath5 for vas , cras and mnas and found the largest value of 1041  k in the case of cras compared with 764  k in vas and 671  k in mnas @xcite . taking mnas as an example , it was shown that the loss of the half - metallicity in a compressed lattice leads to a strong decrease of the curie temperature @xcite . sanyal and collaborators have calculated @xmath5 for v- , cr- and mnas compounds assuming the gaas and inas lattice constants @xcite . their mean field approximation ( mfa ) results agree with the calculations of other authors while the values of @xmath5 obtained with monte carlo simulations are @xmath7100 k lower than the mfa values for vas and mnas and @xmath7300 k for cras . note that the mean - field formula for the curie temperature used in references @xcite is derived for a single - sublattice system . therefore the exchange interaction between 3@xmath8 and @xmath6 atoms is not treated consequently . also the formula of kbler @xcite is derived for an one - sublattice case . in general , the treatment of a multiple - sublattice materials needs solving a system of equations with the number of equations equal to the number of sublattices ( section [ sec2 ] ) . below we show that for materials studied in this paper the one - sublattice approximation is well founded . in other systems , however , it leads to a large error @xcite .    in this work we report the theoretical investigation of the magnetic exchange interactions for selected zb compounds . we study two cr compounds , cras and crse , and four mn compounds , mnas , mnsi , mnge and mnc . mnas is isoelectronic to crse whereas other three mn compounds are isoelectronic to cras . we use the calculated exchange parameters to estimate the curie temperature in both the mean field approximation ( @xmath9 ) and random phase approximation ( @xmath10 ) . a special attention is paid to the influence of the position of the fermi level on the magnetic properties . the paper is organized as follows : in section  [ sec2 ] , we discuss the calculational technique . in sections  [ sec3 ] and [ sec4 ] , we present the results for the cr and mn based half - metallic compounds . section [ sec5 ] is devoted to mnc . finally we summarize in section  [ sec6 ] . the calculations are carried out using the augmented spherical waves method ( asw ) @xcite within the atomic  sphere approximation ( asa ) @xcite . the exchange  correlation potential is chosen in the generalized gradient approximation @xcite . a dense brillouin zone ( bz ) sampling @xmath11 the radii of all atomic spheres are chosen equal . the method of the calculation of the exchange constants has been already presented elsewhere @xcite . here , to make the paper reasonably self - contained a brief overview is given . we describe the interatomic exchange interactions in terms of the classical heisenberg hamiltonian @xmath12 in equation [ eq : hamiltonian2 ] , the indices @xmath13 and @xmath14 number different sublattices and @xmath15 and @xmath16 are the lattice vectors specifying the atoms within sublattices , @xmath17 is the unit vector pointing in the direction of the magnetic moment at site @xmath18 . we employ the frozen  magnon approach @xcite to calculate interatomic heisenberg exchange parameters . the calculations involve few steps . in the first step , the exchange parameters between the atoms of a given sublattice @xmath13 are computed . the calculation is based on the evaluation of the energy of the frozen  magnon configurations defined by the following atomic polar and azimuthal angles @xmath19 the constant phase @xmath20 is always chosen equal to zero . the magnetic moments of all other sublattices are kept parallel to the z axis . within the heisenberg model  ( [ eq : hamiltonian2 ] ) the energy of such configuration takes the form @xmath21 where @xmath22 does not depend on * q * and the fourier transform @xmath23 is defined by @xmath24    in the case of @xmath25 the sum in equation ( [ eq : j_q ] ) does not include @xmath26 . calculating @xmath27 for a regular @xmath28mesh in the brillouin zone of the crystal and performing back fourier transformation one gets exchange parameters @xmath29 for sublattice @xmath13 . the determination of the exchange interactions between the atoms of two different sublattices @xmath13 and @xmath14 is discussed in reference @xcite .    first , the curie temperature is estimated within the mean  field approximation ( mfa ) for the case of a multi  sublattice material by solving the system of coupled equations @xcite @xmath30 where @xmath31 is the average @xmath32 component of @xmath33 and @xmath34 . equation [ eq_system ] can be represented in the form of eigenvalue matrix problem @xmath35 where @xmath36 , @xmath37 is a unit matrix and @xmath38 is the vector of @xmath39 . the largest eigenvalue of matrix @xmath40 gives the value of @xmath9 @xcite .    within the random phase approximation ( rpa ) the curie temperature is given by the relation @xcite @xmath41 where @xmath42 is the spin - wave dispersion . equation ( [ eq_rpa ] ) does not take into account the presence of several sublattices and will be used to estimate the curie temperature with account for the interactions within the sublattice of the 3@xmath8 atoms . as we will show these interactions give a leading contribution to the curie temperature . ll|ccccc compound & a( ) & @xmath43 & @xmath44 & @xmath45 & @xmath46 & @xmath47 + cras  a@xmath48}$ ] & 5.65 & 3.24 & -0.31 & -0.03 & 0.10 & 3.00 + cras  a@xmath49 & 5.53 & 3.17 & -0.25 & -0.03 & 0.11 & 3.00 + cras  a@xmath50 & 5.47 & 3.10 & -0.23 & -0.03 & 0.10 & 2.95 + crse  a@xmath51}$ ] & 5.82 & 3.92 & -0.10 & 0.02 & 0.16 & 4.00 + crse  a@xmath49 & 5.62 & 3.84 & -0.04 & 0.02 & 0.18 & 4.00 + crse  a@xmath50 & 5.56 & 3.77 & -0.04 & 0.02 & 0.16 & 3.92 + we begin with the discussion of the cr compounds : pnictide cras and chalcogenide crse . the half - metallicity in these compounds results from the formation of bonding and antibonding hybrides between the @xmath52 states ( @xmath53 , @xmath54 , and @xmath55 ) of the transition metal atom ( here cr ) and the @xmath56 states of the four neighboring @xmath6 atoms ( as or se ) . the hybridizational gap is partly filled by the @xmath57 states of the @xmath58 atom ( @xmath59 and @xmath60 ) . the hybridization takes place in both spin channels . the position of the hybridization gap is different for majority and minority electrons as a consequence of the exchange splitting @xcite . there are in total four occupied minority states , one coming from the @xmath61 states of the @xmath6 atom and three bonding 3__d _ _ states . the minority - spin @xmath57 orbitals lay above the fermi level . thus the total spin moment @xmath62 in @xmath63 is given by the relation @xmath64 , where @xmath65 is the total number of valence electrons in the unit cell @xcite . in the cases of cras ( 11 valence electrons ) and crse ( 12 valence electrons ) the total spin moment should be 3 and 4 @xmath4 , respectively .    in table [ table1 ] we present the atom - resolved and total spin moments . notations void1 and void2 are used for two types of empty spheres used in our asa calculations to describe correctly the charge distribution in the zinc - blende structure . for each system we performed calculations for three different lattice parameters . the first parameter ( case i ) corresponds to the experimental lattice constant of gaas(cds ) for cras(crse ) . for these lattice constants both compounds are half - metals . this is seen in figure [ fig1 ] where the spin - projected density of states ( dos ) is presented . the total spin moment is exactly integer as expected for a half - metallic system . the cr atoms posses a very large spin moment of 3.24 @xmath63 in cras that reaches the value of 3.92 @xmath63 in crse . a larger moment of crse is a consequence of one extra valence electron per formula unit . both as and se atoms have small induced moments antiparallel to the cr moments . the discussion of the spin moments of the @xmath6 atoms and their relative magnitude can be found in reference @xcite . the half - metallic gap of crse is much larger than for cras . this is explained by the larger cr moment in crse that leads to a larger exchange splitting @xmath66 and , as a result , to a higher position of the minority @xmath57 bands . case in table [ table1 ] ( fermi level at the higher - energy edge of the half - metallic gap ) and the bottom to the a@xmath50 case in the same table ( 1% contraction of the lattice constant with respect to the a@xmath49 case ) . ]    in the next step , we compress the lattice simulating a hydrostatic pressure . the contraction of the lattice constant results in a higher energy position of the majority @xmath56 states with a corresponding increase of the fermi energy , @xmath67 . the calculations are performed for two decreased values of the lattice parameter . the first ( case ii ) corresponds to the position of the fermi level at the bottom edge of the minority - spin conduction band ( fig . [ fig1 ] ) . the second parameter ( case iii ) is obtained by a further 1% decrease of the lattice constant . in the latter case the fermi level lies within the spin - minority conduction band . since half - metallicity is lost the total magnetic moment in case iii becomes non - integer ( table [ table1 ] ) . the contraction of the lattice leads to an increase of the hybridization between the cr states and the states of the @xmath6 atoms . as the result , both the cr moment and the induced moment of the @xmath6 atom decrease . in case ii , the changes of different moments compensate resulting in the same integer total moment . the decrease of the cr spin moment causes a small decrease of the width of the gap . next we used the frozen - magnon approach to calculate the exchange parameters for both systems and three values of the lattice constants ( figure [ fig2 ] ) . the exchange interaction between the induced moments of the @xmath6 atoms is very small and the corresponding parameters are not presented .    for the semiconductor lattice parameter ( left panels in figure [ fig2 ] ) the interactions are characterized by the large ferromagnetic coupling of the nearest cr atoms . in cras , this coupling is about twice the coupling in crse . in the zinc - blende structure the nearest surrounding of a cr atom is formed by 4 @xmath6 atoms . the second - nearest coordination sphere contains cr atoms . although the exchange interaction cr-@xmath6 atom is antiferromagnetic it makes the parallel directions of the cr moments energetically preferable . therefore the indirect exchange interaction of the cr moments through the @xmath6 atoms stabilizes further the ferromagnetic configuration of the cr moments ( a similar situation was found in full - heusler alloys @xcite ) . the interaction is short range and involves only the nearest cr atoms . the influence of the lattice compression is very different for the two systems . in cras , the compression leads to a small decrease of the exchange constants revealing the robustness of the ferromagnetism of the compound . on the other hand , in crse the effect of the compression is very strong . already for the a@xmath49 case where the system is still half - metallic the exchange interaction between the nearest cr atoms is halved with respect to the case of the cds lattice constant . the interaction between second - nearest cr atoms strongly increases keeping antiferromagnetic character . this is reflected in a very small calculated curie temperature ( table [ table2 ] ) . when the lattice is further compressed and the fermi level moves into the minority - spin conduction band , the ferromagnetic interaction between the nearest cr atoms almost vanishes whereas the antiferromagnetic interaction between the second nearest cr atoms increases further . this leads to the instability of the ferromagnetic state . the calculated curie temperatures are collected in table [ table2 ] . they reflect the properties of the exchange interactions discussed above . we report two estimations of the curie temperature within the mfa , one with only the cr - cr interactions taken into account and other with both cr - cr and cr-@xmath6 atom interactions considered , and one rpa estimation .    in case i , for both systems the curie temperature is much higher than the room temperature . for cras , @xmath5 is about twice larger than in crse and exceeds 1000 k. the compression of the lattice leads in the case of cras to a very moderate decrease of the curie temperature . in crse , the @xmath5 drops strongly with decreasing lattice parameter becoming negative in case iii revealing the instability of ferromagnetism . two mean - field estimations give very similar results ( table [ table2 ] ) revealing the leading role of the cr - cr interactions . therefore the neglect of the inter - sublattice cr-@xmath6 atom interactions does not result in a substantial error for these compounds . ll|ccc compound & a( ) & @xmath9(cr - cr ) & @xmath9 & @xmath10 + cras  a@xmath48}$ ] & 5.65 & 1576 & 1600 & 1176 + cras  a@xmath49 & 5.53 & 1498 & 1512 & 1011 + cras  a@xmath50 & 5.47 & 1376 & 1387 & 860 + crse  a@xmath51}$ ] & 5.82 & 706 & 709 & 508 + crse  a@xmath49 & 5.62 & 109 & 117 & 24 + crse  a@xmath50 & 5.56 & -248 & -251&-194 + since the intra - sublattice interactions dominate , the use of the rpa formula ( [ eq_rpa ] ) is justified . overall the @xmath5 values within rpa are by 20 - 30% smaller than the mfa ones . the difference between mfa and rpa estimations is explained by a different weighting of the spin - wave excitations . in the mfa all excitations are taken with equal weight whereas in the rpa the weight decreases with increasing energy of the excitation @xcite .    the rpa is expected to provide better estimation of the curie temperature . in section [ sec5 ] , we will discuss the relation between the mfa and rpa values of @xmath5 in more detail . the only experimental information about the magnetic transition temperature of the systems considered is that the curie temperature of cras is well above 400 k @xcite which is in good correlation with our estimations . in the case of cras we can compare the calculated @xmath5 values of 1600 k within mfa and 1176 k within rpa with the results of previous ab - initio studies . sanyal and collaborators calculated the curie temperature within mfa and found a value of 1320 k for the gaas experimental lattice constant and 1100 k for the inas experimental lattice constant @xcite . their monte carlo simulations gave for the same lattice constants 980 k and 790 k , respectively . kbler used a ginzburg - landau approach with a simplified form of the dynamical susceptibility and a relaxation parameter fitted to reproduce the experimental curie temperature of nimnsb . his estimation of @xmath5 of cras gives 1041 k @xcite . finally , sakuma s mean - field calculation for the theoretical equilibrium lattice constant of 5.82   gave @xmath5 about 1400 k @xcite . summarizing , all calculations agree that the curie temperature of cras should exceed or be about 1000 k. a high @xmath5 value that is stable with respect to the variation of the lattice constant makes cras an interesting candidate for spintronics applications .    ll|ccccc compound & a( ) & @xmath68 & @xmath69 & @xmath45 & @xmath46 & @xmath47 + mnas  a@xmath70}$ ] & 5.87 & 4.16 & -0.26 & -0.01 & 0.10 & 4.00 + mnas  a@xmath49 & 5.73 & 4.09 & -0.20 & -0.01 & 0.11 & 4.00 + mnas  a@xmath50 & 5.68 & 4.01 & -0.19 & -0.01 & 0.10 & 3.91 + mnge  a@xmath48}$ ] & 5.65 & 3.43 & -0.39 & -0.08 & 0.04 & 3.00 + mnge  a@xmath49 & 5.61 & 3.40 & -0.37 & -0.07 & 0.04 & 3.00 + mnge  a@xmath50 & 5.55 & 3.35 & -0.34 & -0.07 & 0.04 & 2.98 + mnsi  a@xmath48}$ ] & 5.65 & 3.51 & -0.49 & -0.07 & 0.05 & 3.00 + mnsi  a@xmath49 & 5.52 & 3.42 & -0.42 & -0.06 & 0.06 & 3.00 + mnsi  a@xmath50 & 5.46 & 3.37 & -0.39 & -0.06 & 0.05 & 2.97 + we proceed with the discussion of mn compounds . here we study mnsi and mnge which are isoelectronic with cras and mnas which is isoelectronic with crse .    in table [ table3 ] we present the atom - resolved and total magnetic moments . mnas with the lattice constant of inp , and mnsi and mnge with the lattice constant of gaas are half - metallic and their total spin moments following the relation presented above are 4 @xmath63 for mnas and 3 @xmath63 for mnsi and mnge ( see table [ table3 ] and figure [ fig3 ] ) . the mn spin moment is about 3.5 @xmath63 for mnsi and mnge , and exceeds the 4 @xmath63 in the case of mnas reflecting the presence of one more valence electron in the latter compound . the moments of the @xmath6 atoms are antiferromagnetically coupled to the mn moments . the value of the half - metallic gap is determined by the value of the exchange splitting and the energy position of the @xmath56 states of the @xmath6 atom while the value of the exchange splitting is governed by the value of the mn moment . as a result , the largest gap is obtained for mnas followed by mnsi and the smallest gap is obtained for mnge .     for the description of details . ]    as for the cr compounds the contraction of the lattice increases the energy of the majority @xmath56 states leading to an increased fermi level . for all three compounds the hybridization between the @xmath56 states of the @xmath6 atom and the @xmath52 states of mn increases and the mn spin moment slightly decreases with the lattice contraction . the smaller mn moments and , as a result , smaller exchange splittings lead to slightly smaller gap - widths .    ) , middle panel to case ii and the right panel to case iii . ]    in figure [ fig4 ] we present the calculated exchange constants . the exchange interaction between the induced moments of @xmath6 atoms is very weak and is not presented . there is strong similarity in the behavior of the exchange parameters for the mn compounds considered here and the cr compounds discussed in the preceding section . the properties of mnge and mnsi are analogous to the properties of the isovalent cras while properties of mnas are similar to the properties of isovalent crse . therefore the corresponding parts of the discussion of the cr compounds remain valid also for the mn systems . ferromagnetism is stabilized by both the direct mn - mn exchange interaction and by the indirect interaction through the @xmath6 atoms . the exchange interactions are short - range and the mn - mn exchange constants for mnas are about a factor two smaller than for mnsi and mnge . the compression has little effect on mnsi and mnge in contrast to mnas . again there is a strong similarity between mnas and crse although the dependence of the exchange parameters on the lattice constant in mnas is somewhat weaker . in particular , the antiferromagnetic exchange interaction between second - nearest mn atoms developing with contraction has lower value than in crse . this results in preserving the ferromagnetic ground state of mnas for all lattice parameters considered while in crse the ferromagnetic state becomes unstable ( table [ table2 ] ) . the curie temperatures of the mn compounds are collected in table [ table4 ] . again , as for cr compounds , the mn - mn interaction is dominant and the account for the exchange interaction mn-@xmath6 atom influences the value of the curie temperature weakly . in agreement with the properties of the exchange constants the curie temperatures of mnge and mnsi depend weakly on the lattice parameter in contrast to mnas where this dependence is strong . the curie temperature of mnas at highest contraction considered assumes the value of about @xmath71th of the curie temperature at the semiconductor lattice constant . the rpa gives values of @xmath5 that are about 30 % smaller than the corresponding values within the mfa . for mnsi and mnge , also the rpa value of @xmath5 remains considerably above the room temperature at all lattice constants studied . the stability of the high @xmath5 value with respect to the variation of the lattice constant makes mnge and mnsi promising candidates for spintronics applications . mnsi has an additional advantage of a larger half - metallic gap . next we compare the results of our calculation of @xmath5 with the results of previous calculations where available . for mnas , we can compare our mfa - estimation of 701 k and rpa estimation of 551 k obtained for the inp lattice constant with the curie temperatures reported by sanyal et al @xcite . the mfa - estimation of sanyal et al gave the value of 640 k for the inas lattice parameter which is very close to the parameter used in our calculations . their monte carlo simulation gave a @xmath5 of 530 k which is very close to our rpa estimation . kbler obtained the value of 671 k for mnas with the inas lattice constant . when he compressed the lattice down to the gaas lattice constant of 5.65  the @xmath5 dropped sharply to 210 k @xcite in agreement with our results . finally , sakuma s mfa - estimation of the curie temperature of mnge for a lattice constant of 5.82   gave a value of about 1200 k which is very close to our mfa value of 1234 k @xcite . thus the results of different calculations are in reasonable agreement with one another .    ll|ccc compound & a( ) & @xmath9(mn - mn ) & @xmath9 & @xmath10 + mnas  a@xmath70}$ ] & 5.87 & 679 & 701 & 551 + mnas  a@xmath49 & 5.73 & 400 & 417 & 275 + mnas  a@xmath50 & 5.68 & 201 & 230 & 136 + mnge  a@xmath48}$ ] & 5.65 & 1160 & 1234 & 716 + mnge  a@xmath49 & 5.61 & 1180 & 1241 & 740 + mnge  a@xmath50 & 5.55 & 1184 & 1231 & 695 + mnsi  a@xmath48}$ ] & 5.65 & 1255 & 1402 & 848 + mnsi  a@xmath49 & 5.52 & 1310 & 1401 & 857 + mnsi  a@xmath50 & 5.46 & 1291 & 1367 & 911 + the last heusler compound studied in this paper is mnc . this system was predicted to be half - metallic by qian and collaborators @xcite . first we discuss the results of the calculation with the theoretical equilibrium lattice constant reported by qian et al . the magnetic properties of mnc differ from the properties of the mn compounds discussed above . mnc has 11 valence electrons : 6 of majority - spin and 5 minority - spin electrons . this electron distribution results in a magnetic moment which is much smaller than in the compounds discussed above ( table [ table5 ] ) . in contrast to the other compounds the half - metallic gap is now in the majority - spin channel ( figure [ fig5 ] ) . the majority - spin dos has the following structure : low in energy in the interval [ -0.55 ry,-0.4 ry ] lie @xmath61-states . above , in the interval [ -0.4 ry,-0.2 ry ] there are three bonding hybrids between the @xmath56 states of c and the @xmath52 states of mn . this is seen in the atom - projected dos of mnc ( right panel in figure [ fig5 ] ) where both mn and c states have similar contribution in this energy interval . the group of states next in energy ( interval [ -0.2 ry,-0.05 ry ] ) does not hybridize with the lower states . these are the mn @xmath57 states separated by a gap from the antibonding c @xmath56-mn @xmath52 hybrids . the fermi level lies within this gap in the majority - spin dos . for the minority - spin channel , the position of the fermi level belongs to the region of the @xmath57 states . these states are shifted by the exchange splitting to higher energies with respect to the corresponding majority - spin states .     and a@xmath50 cases in table [ table5 ] . right panel : atom- and spin - resolved dos for the half - metallic mnc and mnsi compounds . ] the comparison of the dos of half - metallic mnc and mnsi ( figure [ fig5 ] ) allows to reveal the reason for the move of the semiconductor - type behavior from the spin - minority channel in mnsi to the spin - majority channel in mnc . first , the two systems differ by the value of the magnetic moment that leads in the case of mnc to a smaller relative exchange shift of the spin - up and spin - down states and to the occupation of a larger number of the spin - down mn 3@xmath8 states . in particular , the minority - spin @xmath57 states are partly occupied in mnc whereas in mnsi these states lie above the half - metallic gap . second , there is difference in the position of the @xmath57 states in the hybridizational gap between bonding and antibonding @xmath52-@xmath56 hybrids which in the case of mnsi results in the gap between occupied majority @xmath57 and empty antibonding hybrid states .    to extend the study of mnc to the states with the fermi level lying at the top - edge of the half - metallic gap ( case ii ) and above it ( case iii ) we need in this case to expand the lattice ( table [ table5 ] ) opposite to all systems considered above where the lattice must be compressed . this difference comes from the move of the semiconductor channel to the majority - spin subsystem . the expansion of the lattice makes the energy bands narrower and decreases the strength of the @xmath58-@xmath56 hybridization . at the same time the mn moment increases increasing the exchange splitting between majority and minority states . these properties lead to increasing occupation of the majority - spin states and the loss of the half - metallicity ( figure [ fig5 ] ) . the calculated exchange constants are presented in figure [ fig4 ] . the ferromagnetism is stabilized by the direct mn - mn interactions and the mn - c ferromagnetic coupling . since the mn moment is much smaller than in other compounds the value of the mn - mn exchange interaction is substantially lower ( figure [ fig4 ] ) . on the other hand , in mnc the first and second nearest mn atoms are characterized by a similar strength of the ferromagnetic exchange interaction while in other systems only the nearest mn atoms interact strongly . the expansion of the lattice leads to decreasing mn - mn exchange interactions . the calculated curie temperatures are presented in table [ table5 ] . the curie temperature of mnc is smaller than in other mn compounds being about 500 k ( it is comparable with the rpa estimation of @xmath5 of mnas ) . the value of the curie temperature decreases with expansion of the lattice showing again that half - metallicity is favorable for high curie temperature .    ll|ccccc compound & a( ) & @xmath68 & @xmath69 & @xmath45 & @xmath46 & @xmath47 + mnc  a@xmath72 & 4.20 & 1.16 & -0.14 & -0.01 & -0.01 & 1.00 + mnc  a@xmath49 & 4.23 & 1.16 & -0.14 & -0.01 & -0.01 & 1.00 + mnc  a@xmath50 & 4.27 & 1.13 & -0.16 & -0.01 & -0.01 & 1.05 + compound & a( ) & & @xmath9 & @xmath10 + mnc  a@xmath72 & 4.20 & & & 520 & 526 & 507 + mnc  a@xmath49 & 4.23 & & & 467 & 476 & 458 + mnc  a@xmath50 & 4.27 & & & 346 & 363 & 327 +    a remarkable feature of the calculation for mnc is a very small difference between the @xmath5 values estimated within the mfa and the rpa . in section [ sec3 ] we briefly discussed the reason behind the different values of the curie temperatures estimated within mfa and rpa . a striking contrast between mnc and the systems considered in the previous sections allows gaining a deeper insight into the physics of the formation of the curie temperature . we will base our analysis on the comparison of two systems : mnc and mnsi . characterizing the relative difference of the mfa and rpa values of the curie temperature by the relation @xmath73/t_\mathrm{c}^\mathrm{rpa}$ ] we get for mnc a small value of 2.5% compared to a large value of 47% in mnsi .    in mfa , the curie temperature is determined by an arithmetic average of the magnon energies while in rpa @xmath5 is determined by the harmonic average of the same quantities . it is an arithmetic property that the mfa estimation is larger than the rpa one or equal to it if all averaged numbers are equal to each other . in terms of magnon energies , @xmath9 is equal to @xmath10 in the case that the magnon spectrum is dispersion - less : the magnon energy does not depend on the wave vector @xmath28 .     and [ 0 q q ] directions in the brillouin zone . the vertical broken line shows the brillouin zone boundary for the [ 0 q q ] direction . ]    in figure  [ fig6 ] we plot the frozen  magnon energies of mnsi and mnc for two selected directions in the reciprocal space . we see that the mnc curves are flat in the second half of the * q * interval showing here a very weak dispersion . on the other hand , the mnsi curves are much closer to a simple cosinusoidal form with considerable dispersion in this part of the * q * interval . comparing now the low-*q * parts of the curves we notice that mnc curves lie higher than the mnsi curves . therefore , the relative contribution of the low - energy magnons to the rpa value of the curie temperature is less important in mnc compared to mnsi . this combination of features of the wave - vector dependencies of the frozen - magnon energies is responsible for a much smaller difference between the rpa and mfa estimations of the curie temperature in the two cases . since the form of the magnon dispersion is reflected on the properties of the inter - atomic exchange interactions it is instructive to compare the patterns of the exchange interactions of the two systems . taking into account that only mn - mn interactions contribute importantly in the curie temperature we will focus on the comparison of these interactions . the analysis of the mn - mn exchange interactions in mnsi and mnc ( figure  [ fig4 ] ) shows a strong difference between the two systems . in mnsi , only the first - neighbor exchange interaction plays important role . all other interactions are much smaller . on the other hand , in the case of mnc the second - neighbor interaction is comparable with the first - neighbor one . the further few interactions though sizable play less important role . the presence of two first coordination spheres with similar values of the exchange interactions , on the one hand , makes the exchange interaction in mnc more long ranged and , on the other hand , it effectively increases the coordination number of the atoms with leading exchange interaction . either of these properties leads to a better accuracy of the mfa estimation . we have presented a detailed study of the stability of the ferromagnetism in six transition - metal compounds crystallizing in the zinc - blende structure . we calculated from first - principles the electronic structure and employed the frozen - magnon approximation to evaluate the exchange parameters and the curie temperature . in the calculation of the curie temperature we used both the mean - field and random - phase approximations . considering the cras and crse compounds we found in the case of cras a very high curie temperature exceeding 1000 k. the ferromagnetism of cras is stable with respect to the compression of the lattice . on the other hand crse possesses about two times lower curie temperature which decreases strongly with lattice compression . the properties of mnsi and mnge are similar to the properties of the isovalent crse while the properties of mnas are similar to the properties of isovalent crse . mnc presents an interesting case that differs from all other systems studied . this system is also half - metallic with the half - metallic gap in the majority - spin channel in contrast to other systems where the half - metallic gap is in the minority - spin subsystem . the curie temperature of mnc was found to be about 500 k.    in all cases considered we found that half - metallicity is favorable for ferromagnetism and leads to an increased curie temperature . we conclude that two compounds , mnsi and cras , possess properties which make them promising candidates for spintronics applications : large half - metallic gap and high curie temperature . the stability of the curie temperature in these systems allows to expect that being grown on the substrate of a binary semiconductor they will preserve the desired magnetic properties .    the financial support of bundesministerium fr bildung und forschung is acknowledged . 99                                                                  zhang m , hu h , liu g , y. cui , liu z , wang j , wu g , zhang x , yan l , liu h , meng f , qu j and li y 2003 * 15 * 5017 ; zhang m , hu h , liu g , liu z , cui y and wu g 2004 _ j. low temp . phys . _ * 135 * 267 ; zhang m , hu h , liu z , liu g , cui y and wu g 2004 * 270 * 32                                                        anderson p w , 1963 , _ theory of magnetic exchange interactions : exchange in insulators and semiconductors _ , in _ solid state physics _ , edited by f. seitz and turnbull ( academic press , new york ) , vol .
Abstract:

Answer: based on first - principles electron structure calculations and employing the frozen - magnon approximation we study the exchange interactions in a series of transition - metal binary alloys crystallizing in the zinc - blende structure and calculate the curie temperature within both the mean - field approximation ( mfa ) and random - phase approximation ( rpa ) . we study two cr compounds , cras and crse , and four mn compounds : mnsi , mnge , mnas and mnc . mnc , mnsi and mnge are isovalent to cras and mnas is isoelectronic with crse . ferromagnetism is particular stable for cras , mnsi and mnge : all three compounds show curie temperatures around 1000 k. on the other hand , crse and mnas show a tendency to antiferromagnetism when compressing the lattice . in mnc the half - metallic gap is located in the majority - spin channel contrary to the other five compounds . the large half - metallic gaps , very high curie temperatures , the stability of the ferromagnetism with respect to the variation of the lattice parameter and a coherent growth on semiconductors make mnsi and cras most promising candidates for the use in spintronics devises .