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: double bars  systems where a small , faster - rotating `` inner '' bar forms concentrically within a large , `` outer '' bar  were proposed by shlosman , frank , & begelman ( 1989 ) as a mechanism for efficiently feeding gas into the centers of galaxies and fuelling activity there . though isolated examples of real double - barred galaxies were known earlier ( e.g. , de vaucouleurs 1974 ) , high - resolution imaging has revealed numerous examples in the last decade ; see friedli ( 1996 ) for a review . one observational signature is the twisting of isophotes within a large bar ; elmegreen et al .  ( 1996a ) argued that such twists were rather common in early - type barred galaxies . not all such isophote twists are due to inner bars , however : erwin & & sparke ( 1999 ) showed that twists , and peaks in the ellipticity profiles of isophotes , could also be caused by luminous , stellar nuclear rings or by kiloparsec - scale _ disks _ embedded within bars . here , we discuss some results from a survey designed to find out just how common inner bars  and other structures inside bars  might be , and learn something of their characteristics . a more complete description and analysis will appear in erwin & sparke ( 2001a , b ) . to minimize confusion due to dust , we concentrated on early type ( s0 and sa ) galaxies . we chose all barred ( sb or sab ) s0 and sa galaxies in the ugc , north of @xmath0 , which met the following criteria : nearby ( @xmath1 km / s ) , large and easily resolved ( major axis @xmath2 ) , and low to moderate inclination ( major : minor axis ratio @xmath3 ) . the virgo cluster was excluded , in order to concentrate on field galaxies . this resulted in a total of 38 galaxies ( the `` wiyn sample '' ) ; twenty of these are s0 , with ten s0/a and eight sa galaxies . the diameter limit excludes some low - surface - brightness galaxies , but this is probably the only significant bias . we observed all but two of the galaxies in @xmath4 and @xmath5 with the 3.5 m wiyn telescope in tucson , arizona , between december , 1995 , and march , 1998 . seeing ranged from 0.613 , with a median of 08 . we obtained archival hst images , both wfpc2 and nicmos , for over half the galaxies . for ngc 936 , we obtained @xmath4 , @xmath6 , and @xmath5 images ( 07 ) with the 4.2 m william herschel telescope in december , 2000 . we have also begun compiling a catalog of known double - bar and inner - disk systems , combining galaxies from our survey with other wiyn observations and with galaxies drawn from the literature . restricting ourselves to strong , unambiguous detections , we have taken information from a number of studies , including shaw et al . ( 1993 ) ; wozniak et al . ( 1995 ) and friedli et al . ( 1996 ) ; mulchaey , regan , & kundu ( 1997 ) ; jungwiert , combes , & axon ( 1997 ) , mrquez et al . ( 1999 ) ; and greusard et al . we use the published images and ellipse fits  as well as archival hst images wherever available  to remeasure bar position angles and sizes in a consistent manner , and to eliminate mistaken inner bar / disk identifications due to nuclear rings , strong dust lanes , etc . we re - classify some galaxies as inner - disk rather than double - bar , using criteria explained below . to date , we have measurements for a total of 23 double - bar and 17 inner - disk galaxies ; this set is the `` expanded '' sample . to analyze the wiyn sample , we used a combination of ellipse fits to isophotes , unsharp masking , and color maps to identify and measure a variety of central structures in these barred galaxies : inner bars , inner disks , nuclear rings , nuclear spirals , and near - nuclear off - plane dust . we were particularly careful to discriminate between structures which can produce similar distortions in isophotes and ellipse fits ( figure  1 ) . inner bars are surprisingly common : they occur in at least one - quarter ( possibly as many as 40% ) of the barred s0sa galaxies in the wiyn sample , with no discernible dependence on hubble type or outer - bar strength . one galaxy is actually _ triple_-barred ( ngc  2681 ; see erwin & sparke 1999 ) . a typical inner bar in the wiyn sample is 240750 pc in radius ( median = 400 pc ) , 614% the size of its `` parent '' outer bar ( median = 12% ) , and 38% of its galaxy s @xmath7 ( median = 5% ) . in the expanded sample , we find a wider range in relative and absolute sizes ( 2401360 pc , 623% of outer - bar length , 313% of @xmath7 ) , but largely unchanged median values ( 740 pc , 13% , and 6% , respectively ) . we see no preferred angle between inner and outer bars , in either the wiyn or expanded samples . this agrees with previous results ( buta & crocker 1993 ; wozniak et al . 1995 ) and supports models where inner bars rotate at different speeds than outer bars ( maciejewski & sparke 2000 ) . about half of the inner bars in the wiyn sample are surrounded by nuclear rings  dusty , star - forming , or purely stellar . about one - third of the wiyn sample inner bars are devoid of strong dust lanes , while two are found in the midst of substantial _ off - plane _ gas . structurally , inner bars appear similar to the outer bars ( see figure  2 ) , with the `` flat '' luminosity profiles seen in large bars in early - type galaxies ( elmegreen et al . 1996b ) . finally , the presence or absence of inner bars has _ no _ statistically significant effect on nuclear activity in the wiyn sample : galaxies with only one bar are as likely to host agns as those with two .      here , our classification is based purely on orientation : an elliptical structure ( within a large bar ) which is aligned to within 10 of the galaxy s outer disk is called an inner disk , as long as unsharp masking does not show it to be , e.g. , a nuclear ring instead . clearly , this class will include at least some inner bars with chance alignments , along with unresolved nuclear rings and possibly flattened inner bulges . however , there are several pieces of evidence indicating that our inner disks form a physically distinct class .    in the wiyn sample , inner disks are as common as inner bars in s0 galaxies , but practically absent in later hubble types : 35% of the s0 galaxies had an inner disk , while we found one inner disk in an s0/a galaxy and none in sa galaxies . this contrasts strongly with inner bars , which are equally common for s0 , s0/a , and sa galaxies . inner disks are also too common for them _ all _ to be chance alignments of randomly oriented inner bars . finally , inner disks are systematically larger with respect to their parent bars ( median size = 19% for wiyn sample , 20% for expanded sample ) ; k - s tests indicate that inner bars and disks come from different parent populations at a 99.8% confidence level ( 96% for the wiyn sample ) . curiously , inner disks are not significantly larger than inner bars in absolute size , or as a fraction of @xmath7 ; the difference is due mostly to the bars of inner - disk galaxies being systematically smaller ( median radius = 3.9 kpc or 0.39 @xmath7 ) than the ( outer ) bars of double - barred galaxies ( 5.7 kpc or 0.52 @xmath7 ) . at least some of the structures we label inner disks have radial profiles similar to larger - scale galactic disks , and different from bars . figure  2 shows that the inner disk of ngc 3945 ( the largest and brightest in our sample ) differs dramatically in its luminosity profile from a typical inner bar ; it has an exponential profile with a scale length @xmath8 pc , and is distinct from the inner bulge profile . radial luminosity profiles may be a useful tool for determining the true identity of inner disks . kinematic data will be needed to determine if these inner disks are dynamically cool ; this would also distinguish them from oblate exponential bulges . we find nuclear rings in a quarter of the wiyn sample galaxies . most of these rings are dusty and/or star - forming ; such rings are preferentially found in sa galaxies . two of the s0 galaxies have purely stellar nuclear rings , which may be the remnants of past circumnuclear starburst episodes ( see erwin , vega , & beckman , this volume , for more details ) . we find nuclear spirals in @xmath9% of the galaxies ; two of these have blue , luminous stellar arms , indicating recent star formation , but the other six are seen only in dust absorption . finally , 30% of the s0 galaxies have evidence for off - plane gas in the central kiloparsec , in addition to misaligned h i gas at much larger radii . de vaucouleurs , g. 1974 , in _ formation of galaxies , iau symposium no . 58 _ , ed . j.  r. shakeshaft ( dordrecht : reidel ) , 335 elmegreen , d. , et al . 1996 , , erwin , p. , & sparke , l. s. 1999 , , 512 , l37 erwin , p. , & sparke , l. s. 2001a , b , in prep . friedli , d. 1996 , in _ iau colloquium 157 : barred galaxies _ , eds . r. buta , b. g. elmegreen , & d. a. crocker ( san francisco : asp ) , 378 friedli , d. , et al . 1996 , , 118 , 461 greusard , d. , et al . 2000 , , 145 , 425 jungwiert b. , combes f. , & axon , d. j. 1997 , , 125 , 479 maciejewski , w. , & sparke , l. s. 2000 , , , 313 , 745 mrquez , i. , et al .  1999 , , 140 , 1 mulchaey , j. s. , regan , m. w. , & kundu , a. 1997 , , 110 , 299 shaw , m. , et al . 1993 , , 273 , 31 shlosman , i. , frank , j. , & begelman , m. c. 1989 , nature , 338 , 45 wozniak , h. , et al . 1995 , , 111 , 115
Abstract:

Answer: we present results of a high - resolution imaging survey of barred s0sa galaxies which demonstrate that the central regions of these galaxies are surprisingly complex . we see many inner bars  small , secondary bars ( typically less than a kiloparsec in radius ) located inside of , and probably rotating faster than , the large primary bars . these are present in about one quarter to one third of all our sample . in contrast to some theoretical expectations , they do _ not _ seem to enhance agn activity significantly . a third of barred s0 s appear to host kiloparsec - scale _ disks _ within their bars ; but the frequency of such inner disks is much lower in our s0/a and sa galaxies . in addition , we find one example of a _ triple _ barred galaxy , and two cases of purely _ stellar _ nuclear rings  probably the fossil remnants of past circumnuclear starbursts . we comment briefly on results from an ongoing analysis of known double - barred systems , extending to hubble types as late as sbc , and discuss their characteristic sizes and orientations .


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: pentacene as well as other members of the oligoacene family are famous organic semiconductors . single crystals of these materials can be grown in a very high quality , and thus very high charge carrier mobilities could be achieved in e.g. pentacene and tetracene @xcite . these have been exploited in organic field effect transistors in view of fundamental as well as applied aspects @xcite . moreover , organic semiconductors are also of interest for manufacturing large area and flexible devices for light emission or light harvesting @xcite . in this context , the electronic excitations in pentacene play an important role . the optical properties of pentacene and related oligoacenes have been studied previously with the aim to microscopically understand their behavior and gain fundamental insight into the materials itself but also into the entire class of organic semiconductors . one subject of particular interest is the dynamics of excitons in these structures as they are decisive for light absorption and emission properties @xcite . a number of recent studies indicate that the fundamental singlet excitations in oligoacene and other @xmath1-conjugated molecular crystals are of complex nature and can not be rationalized on the basis of molecular frenkel excitons only . instead , a mixture of frenkel and charge transfer ( ct ) excitations has been discussed to be responsible for the observed behavior @xcite . in particular , recent theoretical work suggests that neither the davydov splitting , arising from the interaction of the two inequivalent molecules in the unit cell of oligoacenes , nor the exciton dispersion can be understood without the inclusion of charge transfer excitons @xcite . further , also the spectral shape of the electronic excitation spectrum can depend substantially on the admixture of charge transfer excitations @xcite . moreover , for oligoacenes it has been predicted that the charge transfer contribution to the lowest lying singlet excitons becomes more prominent for larger molecules @xcite . in view of this prediction and due to the fact that the exciton dispersion in organic semiconductors determines the material s behavior in light emitting or absorbing devices , it is highly desirable to know the exciton dispersion in such crystals . for instance , luminescence from exciton states at finite momentum in n , n-dimethylperylene-3,4,9,10-dicarboximide ( meptcdi ) is weak due to the indirect character of the corresponding transition @xcite . previous experimental investigations of the archetype material pentacene have provided insight into the exciton dispersion along some crystal directions , but a complete picture could not be measured due to limitations in the energy resolution @xcite . in this contribution , we present a comprehensive analysis of the energetically lowest lying electronic excitations of pentacene single crystals using electron energy - loss spectroscopy ( eels ) in transmission with unprecedented energy resolution and at low temperature ( 20 k ) . eels is able to determine the energy of electronic excitations at finite momentum transfers throughout the entire brillouin zone of a solid , which can provide valuable insight into the electronic properties of topical materials @xcite . our studies allow the exact analysis of the exciton dispersion in the @xmath0 reciprocal lattice plane of pentacene , where parallel to the @xmath2 direction previous experiments could not probe the singlet exciton at all . moreover , the present investigation provides insight into the complex behavior of the higher energy satellite features , which are discussed in terms of the contribution of vibrations as well as charge transfer excitons . pentacene single crystals with a high quality and typical crystal dimensions of @xmath3 mm length , @xmath4 mm width , and @xmath5 mm thickness were obtained via a directional sublimation of two or three times purified pentacene ( fluka ) . the crystal growth was carried out at temperatures between 280@xmath6 c and 220@xmath6 c in closed , evacuated pyrex ampoules and a horizontal two - zone furnace was used . the growth procedure lasted about four to six weeks .    our eels investigations require thin , single crystalline pentacene films ( thickness @xmath7 100 nm ) which were cut from the largest , flat surface of a single crystal platelet with an ultramicrotome using a diamond knife . these cut films were attached to standard electron microscopy grids@xcite and then transferred into the eels spectrometer.@xcite all eels measurements were carried out using a 170 kev spectrometer thoroughly discussed in a previous publication . @xcite at this high primary beam energy only singlet excitations are possible.@xcite the energy and momentum resolution were 85 mev and 0.03 @xmath8 , respectively . the eels signal , which is proportional to the loss function @xmath9 ) ) , was determined for various momentum transfers , * q * , parallel to the directions of the corresponding reciprocal lattice vectors . moreover , a he flow cryostat allows to choose the sample temperature in the range of 20 k to 400 k.    prior to the investigation of the electronic excitations , our films were thoroughly characterized in - situ using electron diffraction . these investigations clearly document that the films are single crystalline as can be seen from fig . 1 where we depict the electron diffraction profiles parallel to the reciprocal lattice directions @xmath10 and @xmath2 . only the respective bragg diffraction features can be observed . in addition , we have measured these electron diffraction profiles at temperatures of 300 k and 20 k and determined the length of the @xmath10 and @xmath2 reciprocal unit vectors . a comparison to literature data as depicted in fig . 1 ( right panel ) reveals that the crystal structure of our films corresponds to the so - called single crystal polymorph @xcite ( polymorph 1 following the notation of ref . ) . the temperature dependence of the @xmath10 and @xmath2 unit vectors arises from the lattice contraction and is in good correspondence to x - ray diffraction results from the literature as seen in fig . 1 ( right panel ) . note that data for temperatures below 90 k have not been reported so far . finally , the diffraction experiments were used to orient the momentum transfer parallel to selected reciprocal lattice directions in the measurements of the electronic excitation spectrum .     and @xmath11 from our electron diffraction measurements ( square symbols ) in comparison to x - ray diffraction data from the literature . the data for 90 k and 290 k are taken from ref . , while those for 180 k stem from ref . . the very good agreement at 300 k demonstrates that our crystal structure corresponds to the single crystal polymorph ( according to ref . ) ]    in order to properly normalize the data taken along different directions in the reciprocal space , we carried out a kramers - kronig analysis ( kka ) . the raw data first have been corrected by subtracting contributions of multiple scattering processes and elimination of the contribution of the direct beam @xcite . the normalization necessary within the kka analysis has been done using the sum rule for all valence excitations ( for details see ref . ) . the results of the kka analysis show that data for different directions can be properly normalized at 2.8 ev to reveal the relative intensities of the corresponding excitations . furthermore , as molecular crystals often are damaged by fast electrons , we repeatedly checked our samples for any sign of degradation . in particular , degradation was followed by watching an increasing amorphous - like background in the electron diffraction spectra and an increase of spectral weight in the loss function in the energy region below the first excitation feature . it turned out that under our measurement conditions the spectra remained unchanged for about 12 h. samples that showed any signature of degradation were not considered further but replaced by newly prepared thin films . the results from the different films have been shown to be reproducible . ) and for momentum vectors @xmath12 parallel to the reciprocal lattice directions @xmath10 and @xmath2 . ] we start the presentation of our results with a comparison of the loss function for momentum vectors @xmath12 parallel to the reciprocal lattice directions @xmath10 and @xmath2 and a small momentum transfer of 0.08 @xmath8 . at such a small momentum transfer , which is much smaller than the size of the brioullin zone in the @xmath0 plane of @xmath13 0.85 @xmath8 , we probe essentially vertical transitions , the so called optical limit , i.e. a comparison to optical data is possible . following the results of our kka analysis , the data are normalized at 2.8 ev . 2 presents the highly anisotropic electronic response of pentacene in the two directions studied at a sample temperature of 20 k. for oligoacenes , it is well known that the lowest singlet exciton is split into two davydov components , which is a consequence of two non - equivalent molecules in the triclinic unit cell . for pentacene , these two davydov components are spilt by about 0.14 ev , and they are predominantly polarized along the two crystal directions @xmath14 and @xmath15 , respectively @xcite . in general , this is also seen in our data as revealed in fig . however , it is important to realize the our data are taken along the _ reciprocal _ lattice directions @xmath10 and @xmath2 , which due to the triclinic crystal symmetry of pentacene are not parallel to the crystal axes @xmath14 and @xmath15 . in particular , the angle between the crystal axis @xmath15 and its reciprocal counterpart @xmath2 amounts to about 14@xmath6 in pentacene single crystals @xcite . therefore , the spectrum in fig . 2 representing the electronic excitations spectrum along the @xmath2 direction reveals two maxima at 1.84 ev and 1.98 ev , which arise from the two davydov components , while for a momentum transfer parallel to @xmath10 we only observe the lower davydov feature . both curves in fig . 2 also show higher energy satellites starting above 2 ev . the spectral shape and the energy of these satellite features is in very good agreement to the results of optical measurements @xcite . the origin of those will be discussed in more detail below . in comparison the previously published results from eels studies of pentacene single crystals @xcite , the data in fig . 2 allow the identification of the excitonic features for all momentum vectors in the reciprocal @xmath0 plane in detail . ( panel a ) , @xmath16 ( panel b ) , and @xmath2 ( panel c and d ) . the size of the momentum transfer is given in @xmath8.,title="fig : " ]   ( panel a ) , @xmath16 ( panel b ) , and @xmath2 ( panel c and d ) . the size of the momentum transfer is given in @xmath8.,title="fig : " ]   ( panel a ) , @xmath16 ( panel b ) , and @xmath2 ( panel c and d ) . the size of the momentum transfer is given in @xmath8.,title="fig : " ]   ( panel a ) , @xmath16 ( panel b ) , and @xmath2 ( panel c and d ) . the size of the momentum transfer is given in @xmath8.,title="fig : " ]    the momentum dependence of the excitation spectra for three different directions in the reciprocal @xmath0 plane are shown in fig . these data clearly demonstrate the strong momentum dependence of the excitons , i.e. our data provide the exciton band structure . for a momentum transfer parallel to @xmath10 ( fig . 3a ) , we observe a clear upshift of the exciton with increasing momentum transfer . this upshift is accompanied by a decreasing spectral weight . very similar data are observed for a momentum transfer parallel to the diagonal direction ( @xmath16 ) ( fig . 3b ) . the data for @xmath12 @xmath17 @xmath2 are more complex due to the presence of both davydov split exciton components . we therefore present the dispersion data along @xmath2 in two parts . at lower momentum transfers ( up to 0.2 @xmath8 , fig . 3c ) , the lowest energy component disperses to larger energies while the second exciton does not show a significant change . going to larger momenta the upper davydov split component disperses negatively reaching the initial energy of the lower component ( @xmath18@xmath8 ) at about 0.8 @xmath8 ( fig . ( panel a ) , @xmath16 ( panel b ) , and @xmath2 ( panel c ) reciprocal lattice directions . the gray line is intended as a guide to the eye.,title="fig : " ]   ( panel a ) , @xmath16 ( panel b ) , and @xmath2 ( panel c ) reciprocal lattice directions . the gray line is intended as a guide to the eye.,title="fig : " ]   ( panel a ) , @xmath16 ( panel b ) , and @xmath2 ( panel c ) reciprocal lattice directions . the gray line is intended as a guide to the eye.,title="fig : " ]    the exciton dispersion in pentacene is summarized in fig . 4 for the three directions discussed above . in general , our data outline that excitons in single crystalline pentacene can move coherently through the crystal giving rise to a well defined dispersion . the dispersion along @xmath10 and the diagonal ( @xmath16 ) is in perfect agreement to the previous results @xcite . the bandwidth of the exciton bands is about 110 mev . for momentum transfers parallel to @xmath2 , the lower davydov component behaves analogous to the other directions whereas the upper component seems to have a mirror - like dispersion relation . the latter is expected within a picture that relies on local excitons and their predominant nearest neighbor interaction , since the hopping matrix elements between the two different molecules in the unit cell , which thereby are mainly responsible for the exciton movement , are the same as for the davydov splitting @xcite but enter the exciton dispersion with opposite sign . however , we emphasize that previous conclusions , that it is impossible to describe the exciton dispersion using a tight - binding approach based upon nearest neighbor interactions only @xcite , are fully corroborated . within such an approach , which would essentially describe local ( frenkel - like ) excitons , three different transfer integrals ( @xmath19 , @xmath20 , and @xmath21 ) are considered@xcite ( the latter is taken into account as the pentacene crystal unit cell contains two symmetrically inequivalent molecules , see above ) . already the dispersion for momenta parallel to the @xmath10 direction can only be modelled when @xmath19 is negligibly small ( for details see refs . @xmath21 then is fixed to describe this direction . going further to the other two directions presented here , it is not possible to model both with the only left parameter @xmath22 . finally , since the data as presented in this publication are taken at 20 k while previous results have been measured at room temperature ( 300 k ) , the good agreement of the two data sets indicates an exciton dispersion in pentacene that is virtually independent of temperature below room temperature . advanced theoretical considerations of oligoacenes have revealed that the dispersion of the singlet excitons and the observed davydov splitting can only be understood when an admixture of charge transfer excitations to the lowest frenkel - type excitons is taken into account @xcite . moreover , these calculations indicate a 48 % ( ! ) charge transfer character of the energetically lowest singlet exciton in pentacene . as a consequence , it is not surprising that a simple tight - binding approach is insufficient to model the exciton dispersion relations as it starts from a local , molecular ( frenkel - like ) wave function and its nearest neighbor interactions , which is in contrast to the much larger charge transfer exciton wave function . finally , our high resolution data additionally reveal that the spectral width of the excitons is momentum independent throughout the whole measured momentum range . this is in good agreement with estimates of the exciton binding energies in pentacene of the order of 0.5 ev or larger @xcite , i.e. much larger than the exciton band width . in such a case , the processes limiting the exciton lifetime can be expected to be momentum independent . in addition to the well defined exciton features in the range below about 2 ev , there are further excitations visible above this energy for all spectra presented in figs . 2 and 3 . moreover , the momentum dependent behavior in this energy range is rather complex . in general , these excitations above about 2 ev have been discussed in the past on the basis of vibrational satellites and charge transfer excitons . the optical absorption spectrum of individual pentacene molecules in solution @xcite is characterized by a progression of satellite features which most likely are caused by the coupling of the excited electronic state with molecular vibrations . this is similar to other oligoacenes and related molecules . thus , it is reasonable to expect that at least part of the spectral intensity in our spectra above about 2 ev is caused by vibrational satellites . we argue however , that on the basis of vibrational effects the intensity distribution and the observed changes upon momentum variation can not be understood using vibrational degrees of freedom only . the behavior as seen in our data would imply an extremely strong momentum dependence of the electron - vibration coupling , which seems rather unlikely in view of the local ( intra - molecular ) nature of the relevant vibrations . the impact of the strong mixture of frenkel and charge transfer degrees of freedom in pentacene however also changes the higher energy part of our spectra substantially . a 48 % contribution of charge transfer excitons to the lowest singlet exciton feature - as predicted @xcite - in turn also means that the anti - bonding part of this mixed frenkel - charge transfer wave function must have a large frenkel character , as well . moreover , previous experiments using electro - absorption indicate that ( charge transfer ) excitons are located at about 0.3 - 0.5 ev higher than the lowest singlet excitation @xcite . consequently , the higher - lying , anti - bonding component of the mixed wave function also receives from the larger spectral weight of the frenkel - type excitations , and should show up in the corresponding spectra . we therefore attribute the spectral structures above about 2 ev as seen in our data to a superposition of vibrational satellites and excitons of mixed frenkel - ct character , whereas the latter can have a significant momentum dependence , in analogy to the lower component as visualized in fig . 3 . furthermore , also the high energy part of this mixed wave function will couple to the molecular vibrations , which renders the spectral shape even more complicated . since the admixture of frenkel and ct excitons in oligoacenes is expected to depend on the length of the molecules @xcite , a systematic study of the electronic excitation spectra of different oligoacenes might give further insight into the microscopic origin of the various excitations and their interplay . to summarize , we have studied the low energy singlet excitation spectrum of singly crystalline pentacene films at 20 k as a function of momentum transfer using electron energy - loss spectroscopy . our data confirm the substantial anisotropy of the electronic excitation in the @xmath0 reciprocal lattice plane . furthermore , we have determined the full band structure within this plane for the lowest singlet excitons . these excitons clearly disperse with a total band with of about 110 mev . a comparison to previous measurements taken at room temperature indicates that the exciton dispersion is rather temperature independent . following recent theoretical work , the exciton dispersion can only be understood with the inclusion of charge transfer excitons that significantly mix with the molecular - type frenkel excitons . the spectral width of the excitons is momentum independent throughout the entire brioullin zone in the @xmath0 reciprocal lattice plane , i.e. their life time does hardly vary with momentum . inspection of the excitation features somewhat above the exciton bands indicates that they can not be rationalized in terms of vibrational satellites only . instead , there is evidence that also admixed frenkel - ct excitons have to be considered . we are grateful to r. hbel , s. leger and r. schnfelder for technical assistance . financial support by the deutsche forschungsgemeinschaft within the forschergruppe for 1154 ( project kn393/14 ) is gratefully acknowledged .
Abstract:

Answer: we have carried out electron energy - loss investigations of the lowest singlet excitons in pentacene at 20 k. our studies allow to determine the full exciton band structure in the @xmath0 reciprocal lattice plane . the lowest singlet exciton can move coherently within this plane , and the resulting exciton dispersion is highly anisotropic . the analysis of the energetically following ( satellite ) features indicates a strong admixture of charge transfer excitations to the exciton wave function .


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 unprecedented progress of laser technology in the recent decade made available the realization of femtosecond laser sources with relativistic : @xmath0 ( in optical domain ) and ultrarelativistic intensities up to @xmath1  by the chirped pulse amplification technique @xcite . it is expected that the next generation of ultrarelativistic lasers will significantly exceed these intensities in the near future @xcite . the interaction of such powerful electromagnetic radiation with matter opens real possibilities for revelation of ultrafast laser - matter interaction dynamics in supershort time scales with new research fields , such as relativistic optics mourou , rel - rep1,rel - rep2 and attoscience @xcite , as well as many nonlinear electrodynamic effects at the ultrarelativistic interaction with the qed vacuum and accelerator beams @xcite .    at the laser - matter interaction the high - order harmonic generation ( hhg ) by atoms / ions @xcite is a subject of extreme interest , since hhg is one of the most promising mechanisms for creation of coherent xuv sources @xcite and intense attosecond pulses atto3 , in particular , for the control of electrodynamic processes on attosecond time scales . to reach a far x - ray region one needs the atoms or ions with a large nuclear charge and laser fields of ultrahigh intensities at which the nondipole interaction and relativistic effects become essential @xcite . at such intensities the state of ionized electron becomes relativistic already at the distances @xmath2 ( @xmath3 is the wavelength of a laser radiation ) irrespective of its initial state , hence the investigation of laser - atom / ion induced processes such as stimulated bremsstrahlung ( sb ) , above threshold ionization ( ati ) etc . , require relativistic consideration @xcite . specifically , the relativistic drift of a photoelectron due to the magnetic field of a strong electromagnetic wave becomes the major inhibiting factor in the relativistic regime of hhg . due to this drift the significant hhg suppression takes place and diverse schemes have been proposed to overcome this negative effect @xcite . for compensation of the magnetic drift in relativistic domain , in refs . @xcite , scheme5 two counterpropagating laser pulses of linear and circular polarizations were considered , respectively , at which the configuration of a standing wave is formed . the effect of resulting magnetic field of the standing wave is vanished near the stationary maxima at the waves linear polarizations @xcite . while at the circular polarizations of counterpropagating laser pulses one can achieve the fully vanishing of resulting longitudinal magnetic force , responsible for magnetic drift scheme5 . however , the latter takes place at the adiabatic turn on of the waves @xcite , which is not valid for short laser pulses of relativistic intensities . so , a standing wave configuration formed by the two counterpropagating laser beams of linear polarizations is of interest due to the simplicity to realize such field structure providing incomparable large hhg rates in the relativistic regime . at the lengths much smaller than a wavelength of a pump wave the effective field of the standing wave may be approximated by the uniform periodic electric field . note that the consideration of the hhg in the field of a standing wave in the paper @xcite is based on the semiclassical model , and only the first order relativistic effects have been taken into account . the entire relativistic quantum - mechanical consideration of the hhg in a standing wave or in a nonstationary strong electric field on the base of the dirac equation is still lacking .    to overcome the suppression of the hhg in the relativistic regime caused by the magnetic field of a strong laser pulse , in the present work the relativistic hhg by hydrogenlike ions in the uniform periodic electric field , formed by the two linearly polarized counterpropagating laser beams of relativistic intensities , at the distances much less than a laser wavelength , is studied . we propose also a scheme of relativistic hhg in plasma , where a traveling laser pulse in the own frame of reference ( moving with the pulse group velocity in plasma ) is transformed into the uniform periodic electric field , i.e. the wave magnetic field in this case vanishes completely , in contrast to a standing wave configuration formed by the counterpropagating laser beams . so , it is also feasible the effective hhg in plasma with the copropagating ultraintense laser and fast ion beams of the same velocities . note that in refs . @xcite a scheme with counterpropagating relativistic ion and laser beams was considered , where due to the large doppler upshift the high harmonic frequency may appear in the @xmath4-@xmath5 domain .    the organization of the paper is as follows . in sec . ii we present analytic and fully relativistic quantum theory of hhg arising from the dirac equation . in sec . iii with the help of fast fourier transform algorithm we present hhg spectrum for diverse ion and field parameters . finally , conclusions are given in sec . let two linearly polarized plane electromagnetic waves with carrier frequency @xmath6 and amplitude of the electric field @xmath7@xmath8propagating in the opposite directions in vacuum , interplay with the hydrogenlike ions having the charge number of the nucleus @xmath9 . we will assume that @xmath10 , where @xmath11 is the characteristic size of the atomic system and @xmath12 is the wavelength of a pump wave ( for the hhg this condition is always satisfied ) .    at the photoionization of an atom / ion in the strong traveling wave field taking into account the relativistic drift due to the magnetic field , one can expect that the probability of ionized electron recombination with ionic core could be non - negligible only if the electron is initially born with a nonzero velocity oppositely directed to incident laser beam . the probability of tunneling ionization with nonzero initial velocity is given by quantum mechanical tunneling theory @xcite as : @xmath13where @xmath14 is the initial velocity , and @xmath15 is the electric field strength of the wave ( here and below , unless stated otherwise , we employ atomic units ) . however , according to ( [ ion ] ) the ionization probability falls off exponentially if this velocity @xmath14 becomes larger than the characteristic atomic velocities , irrespective of it s direction . since we study the case of superstrong laser fields with @xmath16 ( @xmath17 is the relativistic invariant parameter of the wave intensity ) when the energy of the interaction of an electron with the field over a wavelength becomes comparable to the electron rest energy , the required velocities becomes comparable with the light speed @xmath18 ( @xmath19 @xmath20 ) . hence , the probability ( [ ion ] ) in the such fields is practically zero . therefore , in considering case of a standing wave formed by the laser beams ( [ waves ] ) , a significant input in the hhg process will be conditioned by the ions situated near the stationary maxima of the standing wave . for this points the magnetic fields of the counterpropagating waves cancel each other . since the hhg is essentially produced at the lengths @xmath21 , we will assume the effective field to be : @xmath22where @xmath23 , and @xmath24 is the unit polarization vector . as is clear from the above consideration for hhg , one needs to exclude the relativistic drift of a photoelectron due to the magnetic field of a traveling wave , i.e. the magnetic component of the wave . this can also be achieved in the plasma- like medium with a refractive index @xmath25 . indeed , let a plane , transverse , and linearly polarized electromagnetic wave with a frequency @xmath6 propagates in plasma ( @xmath26 ; @xmath27 is the plasma frequency ) . for the electric and magnetic field strengths we have@xmath28 ,   \label{pet}\]]where suppose that an ion beam copropagates together a laser beam in such plasma with a mean velocity @xmath30 equal to a laser beam group velocity : @xmath31 . to find out the hhg rate it is convenient to solve the problem in the center - of - mass ( @xmath32 ) frame of ions . in this frame the wave vector of the photons is @xmath33 . the traveling electromagnetic wave is transformed in the @xmath32 frame into the field avet - book:@xmath34@xmath35hence , in the frame of ions center - of - mass the problem of hhg in a plasma is reduced to one in vacuum in the field of the same configuration ( [ et ] ) with @xmath36 . when @xmath37 , one needs relativistic ion beams . concerning the ion beams of necessary relativism @xmath38 , note that the relativistic ion beams in arbitrary charge states , with the lorentz factor up to about @xmath39 is supposed to be realized at the new accelerator complex at gesellschaft fr schwerionenforschung ( gsi ) ( darmstadt , germany ) . it is worthy to note that if a laser frequency @xmath6 is close enough to plasma frequency @xmath40 , one can achieve the implementation of proposed scheme for hhg with nonrelativistic ion beams . in this case one can use quasimonoenergetic and low emittance ions bunches of solid densities generated from ultrathin foils -nanotargets by supershort laser pulses of relativistic intensities target .    to find out the relativistic probabilities of hhg in the field ( [ et ] ) we arise from the dirac equation:@xmath41where @xmath42 and @xmath43 are the dirac matrices in the standard representation , @xmath44 are the pauli matrices and @xmath45 is the operator of the kinetic momentum ( @xmath46 ) . without loss of generality one can take the polarization vector @xmath24 aligned with the @xmath47  axis of spherical coordinates . we denote the atomic bound states by @xmath48 , where @xmath49 indicates the set of quantum numbers that characterizes the state @xmath50 . here @xmath51 is the principal quantum number , @xmath52 is the whole moment , @xmath53 is the orbital moment and @xmath54 is the magnetic quantum number .    following the ansatz developed in the ref . @xcite , the time dependent wave function can be expanded as@xmath55here@xmath56are the dirac free solutions @xcite with energy @xmath57 and polarization states @xmath58:@xmath59as an initial bound state wave function @xmath60 we assume the ground - state bispinor wave function @xcite for the hydrogenlike ion with the quantum numbers @xmath61 , @xmath62 , @xmath63 , and @xmath64:@xmath65here @xmath66 is the euler gamma function , @xmath67and @xmath68 are the polar and azimutal angles , @xmath69 and @xmath70 is the energy of the ground state . we assume @xmath71 , and the parameter @xmath72 can take values lying in the interval @xmath73 . in the expansion ( [ exp ] ) we have excluded the negative energy states , since the input of the particle - antiparticle intermediate states will lead only to small corrections to the processes considered . neglecting the depletion of the ground state @xmath74 and the free - free transitions due to coulomb field , the dirac equation for @xmath75 reads as:@xmath76@xmath77where@xmath78is the atomic dipole matrix element for the bound - free transition . the latter can be calculated with the help of integral @xcite:@xmath79where @xmath80 is an associated legendre function of degree @xmath53 and order @xmath81 , @xmath82 is a bessel function of order @xmath83 . thus , for @xmath84 we have :    @xmath85    @xmath86    @xmath87 ,   \label{d1}\ ] ]    @xmath88 ,   \label{d_1}\ ] ]    where@xmath89    the last term in eq . ( [ c ] ) describes spin flip due to free - free transitions and proportional to transverse momentum @xmath90 . as will be seen below , the main contribution to the hhg process is conditioned by the electrons with momentum components @xmath91 . hence , we will neglect the last term in eq . ( [ c ] ) , which is valid for the fields @xmath92 , where @xmath93 is the schwinger field strength . then , from eq . ( [ c ] ) for the probability amplitudes @xmath94 we obtain:@xmath95@xmath96 dt^{\prime \prime } \right\ } ,   \label{cp}\]]where @xmath97 is the vector potential of resulting electric field .    for the harmonic radiation perpendicular to the polarization direction @xmath98 one needs mean value of the @xmath47 component of the electron current density @xmath99 . using eqs . ( [ exp ] ) , ( [ cp ] ) and neglecting the contribution by free - free transitions , we obtain : @xmath100@xmath101@xmath102where @xmath103@xmath104is the relativistic classical action of an electron in the field and @xmath105 .    as in the nonrelativistic case , the hhg rate is mainly determined by the exponential in the integrand of eq . ( [ jz ] ) with exact relativistic classical action . the integral over the intermediate momentum @xmath106 and time @xmath107 can be calculated using saddle - point method . the saddle momentum  is determined by the equation:@xmath108which for momentum components gives @xmath91 , and the @xmath47 component of momentum ( @xmath109 ) is given by the solution of the equation : @xmath110 in contrast to nonrelativistic @xcite and relativistic @xcite cases for a traveling wave , this equation can not be solved analytically and requires numerical solution . integrating the latter over @xmath106 , for the current density we obtain:@xmath111@xmath112where@xmath113note that contribution from spin flip due to bound - free transitions in the saddle - point vanishes : @xmath114 . at @xmath115 , for @xmath116 one can recover nonrelativistic result : @xmath117 . the complex saddle times @xmath118 are the solutions of the following equation : @xmath119which may be expressed by the transcendental equation @xmath120then expressing the saddle time as @xmath121 , with @xmath122 , one can obtain the saddle momentum @xmath123and the imaginary part of the saddle time:@xmath124where @xmath125 . taking into account eq . ( [ ps ] ) , from eq . ( [ sad_r ] ) for the real part of the saddle time we obtain @xmath126as usual , @xmath127 is interpreted as the birth time of the photoelectron which returns at the moment @xmath128 to the core and generates harmonic radiation . the transition dipole moment has singularity at the saddle times and the integration has been made with the help of the formula:@xmath129 ^{\frac{\nu -1}{2 } }   \notag \\ & & \times \frac{\gamma \left ( \nu /2\right ) e^{-\lambda f\left ( x_{0}\right ) } % } { \gamma \left ( \nu \right ) } . \label{asim}\end{aligned}\]]thus , taking into account eqs . ( [ sad_p]-[asim ] ) , we obtain the ultimate formula for current density:@xmath130     and ionization time versus the electron time evolution in the continuum ( return time ) for the various laser intensities . ]    the formula ( [ main ] ) is the analogous to the nonrelativistic formula for the dipole moment in the three step model @xcite . here the summation is carried out over the solutions of eq . ( [ ret ] ) . the tunneling ionization amplitude @xmath131 is : @xmath132@xmath133  the propagation amplitude is given by the expression @xmath134and the recombination amplitude is : @xmath135 , where    @xmath136    @xmath137 . \label{crec}\ ] ] in this section we present some numerical simulations for hhg when relativistic effects become essential . in considering relativistic case the saddle time and energy gain essentially depend on the intensity of the pump wave . the latter leads to the modification of the hhg spectrum compared with nonrelativistic one . in fig . 1 we present the solution of eq . ( [ ret ] ) for the born times @xmath127 which are limited to a quarter of the laser period . the fig . 1 also illustrated the photoelectron energy gain versus the electron s time evolution in the continuum ( return time ) for the various laser intensities . it is well known that the cutoff frequency of hhg is defined by the maximal classical energy gain . for this reason in fig . 2 the maximum energy gain of the photoelectron as a function of the relativistic parameter of the wave intensity @xmath138 is displayed . as we see , the relativistic cutoff essentially differs from the nonrelativistic one for @xmath139 and the shift of the cutoff position to the lower values of the harmonic order for the same laser intensity becomes evident . 1 also reveals the multiplateau character of the harmonic spectrum like to the nonrelativistic one .    to find out harmonic spectrum , the fast fourier transform algorithm has been used . , as a function of the harmonic order , for an ion with @xmath140 , @xmath141 , and frequency @xmath142 ( @xmath143 @xmath144 ) . the gray curve represents the hhg spectrum for a traveling wave . ] ( @xmath145 @xmath144 ) and @xmath146 . ] standing wave intensity @xmath147 ( @xmath148 ) , and frequency @xmath149(@xmath143 @xmath144 ) . the gray curve represents the hhg spectrum for a traveling wave . ] figures 3 and 4 display harmonic emission rate via modulus squared of the fourier transform of the current density @xmath150 as a function of the harmonic order @xmath151 for an ion with @xmath140 . the frequencies and relativistic parameters of intensity are : @xmath152 ( @xmath153 ) and @xmath154 ( @xmath155 ) , respectively . figure 5 displays harmonic emission rate for an ion with @xmath156 and field intensity @xmath147 ( @xmath157 ) with the wave frequency : @xmath158 . for the comparison we have also brought the spectra for a traveling wave which was obtained using the results of ref . @xcite . as we see from these figures , with the increase of the laser intensity the hhg rate with a standing wave field by many orders of magnitude is larger than hhg rate with a traveling wave . besides , the relativistic quantum mechanical harmonic cutoff is shifted to the lower values compared with nonrelativistic one .    , @xmath159 , and @xmath160 ( @xmath161 @xmath162 ) . ]    , @xmath159 , and @xmath163 @xmath164 . ]    . ] figures 6 , 7 , and 8 display harmonic emission rates when the relativistic parameter of the wave intensity is : @xmath165 . in this case the relativistic effects are essential . figure 6 corresponds to @xmath166 and @xmath167 . figures 7 and 8 correspond to @xmath168 with @xmath152 and @xmath154 , respectively . for the last three figures the cutoff position is @xmath169 , which is in good agreement with semiclassical analysis ( see fig . 2 ) . for these parameters the harmonic emission rates with a traveling wave are negligible , meanwhile the relativistic rates with a standing wave field at photons energies approaching to @xmath170 region , are considerable . we have presented a theoretical treatment of the hhg by hydrogenlike ions with a large nuclear charge in the relativistic regime with a standing wave configuration , formed by the two linearly polarized counterpropagating laser beams of relativistic intensities , at the distances much less than a laser wavelength . the investigation of the hhg with the such pump field configuration is conditioned by known fact that in the relativistic regime the significant hhg suppression takes place in the field of a strong traveling wave because of the magnetic drift of a photoelectron along the wave propagation direction . we have proposed a new scheme for relativistic hhg with an ultraintense traveling laser pulse and copropagating ions in a plasma . in this case the relativistic hhg occurs exactly without the magnetic drift of the photoelectron , because the traveling wave in the ions own frame of reference ( moving with the laser pulse group velocity ) in underdense plasma is transformed into the purely uniform periodic electric field ( in contrast to a standing wave  counterpropagating laser beams ) . hence , we expect that the efficient hhg with the current ultraintense laser pulses and forthcoming relativistic ion beams ( at gsi - darmstadt accelerator ) will be feasible in the plasmas of conventional densities .    on the base of the dirac equation we have presented the analytic and fully relativistic quantum theory of hhg . with the help of the dynamic solution of the dirac equation and fast fourier transform algorithm we have calculated the harmonic spectrum . the obtained results have been applied to hydrogenlike ions with a moderate nuclear charge . the harmonic spectrum displays the significant difference compared both to non - relativistic and relativistic spectra with a traveling laser pulse . in particular , the shift of the cutoff position to the lower values of the harmonic order is considerable . in the relativistic regime of interaction the harmonic emission rates for a traveling laser pulse are negligible , meanwhile the relativistic rates with a standing wave field , at photons energies approaching to @xmath170 region , are considerable . esirkepov , s. v. bulanov , k. nishihara , t. tajima , f. pegoraro , v. s. khoroshkov , k. mima , h. daido , y. kato , y. kitagawa , k. nagai , and s. sakabe , phys . lett . * 89 * , 175003 ( 2002 ) ; s. ter - avetisyan , m. schnrer , p. v. nickles , m. kalashnikov , e. risse , t. sokollik , w. sandner , a. andreev , and v. tikhonchuk , phys . lett . * 96 * , 145006 ( 2006 ) ; x. q. yan , h. c. wu , z. m. sheng , j. e. chen , j. meyer - ter - vehn , phys . rev . lett . * 103 * , 135001 ( 2009 ) .
Abstract:

Answer: high - order harmonic generation ( hhg ) by hydrogenlike atoms / ions in the uniform periodic electric field , formed by the two linearly polarized counterpropagating laser beams of relativistic intensities , is studied . the relativistic quantum theory of hhg in such fields , at which the impeding factor of relativistic magnetic drift of a strong wave can be eliminated , is presented arising from the dirac equation . specifically , a scheme of hhg in underdense plasma with the copropagating ultraintense laser and fast ion beams is proposed .


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: for the hybrid monte carlo algorithm ( hmc)@xcite , often used to study quantum chromodynamics ( qcd ) on the lattice , one is interested in efficient numerical time integration schemes which are optimal in terms of computational costs per trajectory for a given acceptance rate . high order numerical methods allow the use of larger step sizes , but demand a larger computational effort per step ; low order schemes do not require such large computational costs per step , but need more steps per trajectory . so there is a need to balance these opposing effects . omelyan integration schemes @xcite of a force - gradient type have proved to be an efficient choice , since it is easy to obtain higher order schemes that demand a small additional computational effort . these schemes use higher - order information from force - gradient terms to both increase the convergence of the method and decrease the size of the leading error coefficient . other ideas to achieve better efficiency for numerical time integrators are given by multirate or nested approaches . these schemes do not increase the order but reduce the computational costs per path by recognizing the different dynamical time - scales generated by different parts of the action . slow forces , which are usually expensive to evaluate , need only to be sampled at low frequency while fast forces which are usually cheap to evaluate need to be sampled at a high frequency . a natural way to inherit the advantages from both force - gradient type schemes and multirate approaches would be to combine these two ideas .    previously , we studied the behavior of the adapted nested force - gradient scheme for the example of the @xmath0-body problem @xcite and would like to learn more about their usefulness for lattice field theory calculations . due to the huge computational effort required for qcd simulations , it is natural to attempt an intermediate step first . we chose the model problem of quantum electrodynamics ( qed ) in two dimensions , the schwinger model @xcite , since it is well - suited as a test case for new concepts and ideas which can be subsequently applied to more computationally demanding problems @xcite . as a lattice quantum field theory , it has many of the properties of more sophisticated models such as qcd , for example the numerical cost is still dominated by the fermion part of the action . the fact that this model , with far fewer degrees of freedom , does not require such large computational effort makes it the perfect choice for testing purposes . we compare the behavior of numerical time integration schemes currently used for hmc @xcite with the nested force - gradient integrator @xcite and the adapted version introduced in @xcite . we investigate the computational costs needed to perform numerical calculations , as well as the effort required to achieve a satisfactory acceptance rate during the hmc evolution . our goal is to find a numerical scheme for the hmc algorithm which would provide a sufficiently high acceptance rate while not drastically increasing the simulation time . the paper is organized as follows . in section 2 we give a short overview of the hmc algorithm and numerical schemes for time integration , which are used in hmc . in section 3 we present the 2-dimensional schwinger model and introduce the idea of the force - gradient approach and the resulting novel class of numerical schemes . section 4 is devoted to the results of a comparison between widely used algorithms and the new approach and section 5 draws our conclusion . in this section we provide a general overview of the hmc algorithm @xcite to introduce our novel integrator . we also present some standard numerical time integrating methods used in hmc , as well state - of - the - art numerical schemes , which we later compare by applying them to the two - dimensional schwinger model .      in the hybrid monte carlo algorithm , the quantum lattice field theory is embedded in a higher - dimensional classical system through the introduction of a fictitious ( simulation ) time @xcite . the gauge field @xmath1 is associated with its ( fictitious ) conjugate momenta @xmath2 , and the classical system is described by the hamiltonian , @xmath3 + { \mathcal{b}}[u],\ ] ] where @xmath4 $ ] and @xmath5 $ ] represent the kinetic and potential energy respectively .    for a given configuration @xmath1 , a new configuration @xmath6 is generated by performing an hmc update @xmath7 , which consists of two steps :    * * molecular dynamics trajectory : * evolve the gauge fields @xmath1 , elements of a lie group , and the momenta @xmath2 , elements of the corresponding lie algebra , in a fictitious computer time @xmath8 according to hamilton s equations of motions @xmath9 since analytical solutions are not available in general , numerical methods must be used to solve the system of eqn .  . the discrete updates of @xmath1 and @xmath2 with an integration step @xmath10 are @xmath11 leading to a first - order approximation at time @xmath12 . since the momenta @xmath2 are elements of lie algebra , we have an additive update of @xmath2 . on the other hand , the links @xmath1 must be elements of the lie group , therefore an exponential update is used for @xmath1 to preserve the underlying group structure . * * metropolis step : * accept or reject the new configuration @xmath13 with probability @xmath14 where @xmath15 .      in this paper we are concerned with numerical time integration schemes , which preserve the fundamental properties of geometric integration , time - reversibility and volume - preservation . all numerical schemes presented below possess these necessary properties . * basic schemes : * well - known , commonly used integration schemes in molecular dynamics are given by    * the leap - frog method , a 3-stage composition scheme of the discrete updates defined above : @xmath16 * and a 5-stage extension widely used in qcd computations : @xmath17    * force gradient schemes : * force - gradient schemes increase accuracy by using additional information from the force gradient term @xmath18 , with @xmath19 defining lie brackets . the 5-stage force - gradient scheme proposed by omelyan et al @xcite is the simplest ; @xmath20 here we also test the modification of the force - gradient method proposed in @xcite , where the force - gradient term @xmath21 is approximated via a taylor expansion . an extension is given by the 11-stage decomposition @xcite , recently implemented as the integrator in the open source code openqcd as one of the standard options @xcite @xmath22where @xmath23 , @xmath24 , @xmath25 and @xmath26 are parameters from equation ( 71 ) in ref . @xcite .    * nested schemes : * qed and qcd problems usually lead to hamiltonians with the following fine structure @xmath27 + { \mathcal{b}}_{1}[u]+ { \mathcal{b}}_{2}[u],\ ] ] where the action of the system can be split into two parts : a fast action @xmath28 such as the gauge action , and a slow part @xmath29 , for example , the fermion action . this allows us to apply the idea of multirate schemes ( an idea known as nested integration in physics literature)@xcite in order to reduce the computational effort . at first we consider the nested version of the leap - frog method @xmath30 where the inner cheaper system @xmath31+{\mathcal{b}}_{1}[u]$ ] is solved by @xmath32 with @xmath33 being a number of iterations for the fast part of the action . our main goal is to compare the above - mentioned methods with more elaborated nested schemes : in @xcite , a similar 5-stage decomposition scheme has been recently introduced : @xmath34    a nested version of , which has been used in  @xcite reads @xmath35 where @xmath36 with @xmath37 and @xmath38 . in the limit @xmath39 we have @xmath40 . note that this approach uses force - gradient information at all levels , i.e. , the high computational cost of high order schemes appears at all levels . one may overcome this problem by using schemes of different order at the different levels without losing the effective high order of the overall multirate scheme . for the latter , we include appropriate force gradient information as we explain in the following for the case of two time levels , where the gauge action plays the role of the fast and cheap part , and the fermionic action plays the role of the slow and expensive part . the reasoning is as follows : if one uses the 5-stage sexton - weingarten integrator of second order for the slow action , and approximates the fast action by @xmath41 leap - frog steps of step size @xmath42 , the error of the overall multirate scheme will be of order @xmath43 . with the use of force gradient information only at the slowest level it is possible to cancel the leading error term of order @xmath44 . as @xmath45 usually holds in the multirate setting , the overall order is then given by the leading error term of order @xmath46 , i.e. , the scheme has an effective order of four . one example for such a scheme for problems of type is given by the 5-stage nested force - gradient scheme introduced in @xcite @xmath47 to summarize , the adapted scheme   differs from the original one   in two perspectives :    * the force gradient scheme for the fast action is replaced by a leap - frog scheme . * only the part @xmath48 of the full force gradient @xmath49 is needed to gain the effective order of four . the numerical schemes - and - are second order convergent schemes . methods - and - have the fourth order of convergence . we do not consider integrators of higher order than four since the computational costs are too high . the schemes of the same convergence order differ from each other by the number of stages ( updates of momenta and links per time step ) . usually methods with more stages have smaller leading error coefficients and therefore have better accuracy , but higher computational costs . we would like to determine which integrator would represent the best compromise between high accuracy and computational efficiency .    we will apply all these numerical integration schemes  to the two - dimensional schwinger model . the most challenging task from the theoretical point of view is to derive the force - gradient term @xmath21 . in the next section we introduce the schwinger model and explain how to obtain the force - gradient term . the 2 dimensional schwinger model is defined by the following hamiltonian function @xmath50 = \frac{1}{2}\sum_{n=1,\mu=1}^{v,2 } p_{n,\mu}^2 + s_g[u ] + s_f[u].}\ ] ] with @xmath51 the volume of the lattice . unlike qcd , where @xmath52 and @xmath53 , for this qed problem , the links @xmath1 are the elements of the lie group @xmath54 and the momenta @xmath55 belong to @xmath56 , which represents the lie algebra of the group @xmath54 . this makes this test example very cheap in terms of the computational time . this together with the fact that the schwinger model also shares many of the features of qcd simulations , makes the schwinger model an excellent test example when considering numerical integrators : a fast dynamics given by the computationally cheap gauge part @xmath57 $ ] of the action demanding small step sizes , and a slow dynamics given by the computationally expensive fermion part @xmath58 $ ] allowing large step sizes . the pure gauge part of the action @xmath59 sums up over all plaquettes @xmath60 in the two - dimensional lattice with @xmath61 and is given by @xmath62 the links @xmath1 can be written in the form @xmath63 and connect the sites @xmath0 and @xmath64 on the lattice ; @xmath65 $ ] , @xmath66 , @xmath67 @xmath68 are respectively space and time directions and @xmath69 is a coupling constant . note that from now on we will set the lattice spacing @xmath70 . the fermion part of the action @xmath71 is given by @xmath72 where @xmath26 is a complex pseudofermion field . here , @xmath73 denotes the wilson  dirac operator given by @xmath74 where @xmath75 are the pauli matrices @xmath76 @xmath77 is the mass parameter and the kronecker delta @xmath78 acts on the pseudofermion field by @xmath79 with @xmath80 the pseudofermion field , a vector in the two - dimensional spinor space taking values at each lattice point @xmath0 . in order to proceed with the numerical integration we need to obtain the force @xmath81 and the force gradient term @xmath21 . the force term @xmath82 with respect to the link @xmath83 is given by the first derivative of the action @xmath84 and can be written as @xmath85 since the numerical schemes  use the multi - rate approach , the shifts in the momenta updates are split on @xmath86 and @xmath87 and we can consider them separately . the force terms @xmath86 and @xmath87 are obtained by differentiation over @xmath54 group elements , which for the schwinger model is the standard differentiation .    the force associated with link @xmath88 from the gauge action is given by @xmath89 the force term of the fermion part is given by @xmath90 \,\ ] ] where vectors @xmath91 and @xmath92 are given @xmath93    for the numerical methods and we need to find the force gradient term @xmath94 with respect to the link @xmath83 . in case of the schwinger model this term reads @xmath95    for simplicity we decompose the force gradient term in four parts @xmath96 this decomposition is also useful since the numerical integrator only uses the term @xmath97 by construction . as shown in @xcite , to obtain the fourth order convergent scheme from the second order convergent method we must eliminate the leading error term , which is exactly represented by @xmath97 . for completeness we discuss all 4 parts below . the @xmath98 part of the force - gradient term is @xmath99 \end{aligned}\ ] ] with the set of plaquettes @xmath100 then by using the vectors @xmath101 defined in we obtain the @xmath102 piece of the force - gradient term given by @xmath103 . } \end{aligned}\ ] ] the second derivative of the fermion action is @xmath104 \xi + } \nonumber \\ & & { 2 \operatorname{re } \chi^\dagger \frac{\partial d}{\partial q_\mu(n ) } ( d^\dagger d)^{-1 } \frac{\partial d^\dagger}{\partial q_\nu(m ) } \chi \ , , } \nonumber \\ & & { = 2 \operatorname{re } \left [ z_{1,m,\nu}^\dagger \frac{\partial d}{\partial q_{\mu}(n ) } \xi +   \chi^\dagger \frac{\partial d}{\partial q_{\mu}(n ) } d^{-1 } w_{2,m,\nu } -   \chi^\dagger   \frac{\partial^2 d}{\partial q_{\nu}(m ) \partial q_{\mu}(n ) } \xi +   \chi^\dagger \frac{\partial d}{\partial q_{\mu}(n ) } d^{-1 } z_{1,m,\nu}\right ] } \nonumber \\ & & { = 2 \textrm{re } \left [ z_{1,m,\nu}^\dagger w_{2,n,\mu } +   w_{1,n,\mu}^\dagger z_{2,m,\nu } -   \chi^\dagger   \frac{\partial^2 d}{\partial q_{\nu}(m ) \partial q_{\mu}(n ) } \xi \right ] } \end{aligned}\ ] ] in terms of the vectors @xmath105 and @xmath92 defined in . now the fields @xmath106 and @xmath107 are given by @xmath108 with @xmath109    in order to calculate @xmath110 and @xmath111 it is possible to perform the summation of @xmath112 before the inversions of @xmath73 and @xmath113 to get @xmath114 and @xmath115 which save @xmath116 additional inversions for the force gradient terms . it follows for the force gradient term @xmath111 @xmath117\ ] ] with @xmath118 + z_1 \right ) \ , . \end{aligned}\ ] ] the expression for @xmath110 can be obtained from the one for @xmath111 by replacing in and the vector @xmath119 with @xmath120 defined in . it is important to mention that the computationally most demanding part of the numerical integration of the schwinger model and quantum field theory in general is the inverse of the dirac operator @xmath121 . every momenta update , which includes fermion action requires 2 inversions of the dirac operator , the addition of the force - gradient term @xmath21 requires 4 more inversions . therefore leap - frog based methods and need 4 computations of @xmath121 per time step ; schemes and 6 times ; force - gradient based methods 8 for and , 10 for and the 11 stage method has 12 inversions of the dirac operator . since we use the multi - rate approach for schemes , and , which leads generally to fewer macro time steps needed than for the standard schemes we expect the integrator will be the most efficient choice among the methods considered . in the next section we present numerical tests of this prediction . in this section we apply the numerical integrators  to compute the molecular dynamics step for the schwinger model when studied with the hmc algorithm . we consider a @xmath122 by @xmath122 lattice with a coupling constant @xmath123 and mass @xmath124 . the parameters were taken from @xcite and correspond to the scaling variable @xmath125 defined in @xcite.we have chosen them to simulate close to the scaling limit with light fermions and also to increase the impact of the fermion part of the action . we use one thermalised gauge configuration . for each integrator and value of the step - size we generate @xmath126 independent sets of momenta and integrate the equations of motion on a trajectory of length @xmath127 . we compute the absolute error @xmath128 and estimate its statistical error from the standard deviation . also the parameter @xmath33 is chosen in such a way to make micro step size to be @xmath129 times smaller than the macro step size @xmath10 . figure [ fig:1 ] presents the comparison between the numerical integrators  . it shows the absolute error @xmath128 versus the step - size of the numerical scheme . here the multi - rate schemes , , and outperform their standard versions as expected . also it is easy to see that the scheme has the best accuracy and the nested force - gradient method just slightly edges the adapted nested force - gradient scheme . figure [ fig:2 ] presents the cpu time , required for the proposed integrators  , versus the achieved accuracy . we can observe that the nested force - gradient method and adapted nested force- gradient method show much better results in terms of a computational efficiency than the integrators and ; and even compared to the 11 stage scheme . here we can see that the modification of proposed in @xcite also performs better than its original version . it shows almost similar computational costs as nested versions of the force - gradient approach - , since it has the same number of @xmath121  ( see table [ tab:1 ] ) . but it is less efficient because the proposed nested approach is more precise . .step - sizes and number of inversions of @xmath73 per step and per trajectory for acceptance rate of 90% [ cols="^,^,^,^,^",options="header " , ]     table [ tab:1 ] shows the number of inversions of the dirac operator @xmath73 , which is needed to reach 90% acceptance rate of the hmc . since @xmath121 is the most computationally demanding part it is important to see how many of these inversions are required per each trajectory . from table [ tab:1 ] it easy to see that the adapted nested force - gradient method and nested force - gradient method need the least number of @xmath121 per trajectory to reach the chosen acceptance rate @xmath130 . we can also claim that methods and have a potential to perform even better with respect to the computational effort in the case of lattice qcd problems , since the impact of the fermion action and the computational time to obtain the inversion of the dirac operator @xmath73 is much more significant . we presented the nested force - gradient approach and its adapted version applied to a model problem in quantum field theory , the two - dimensional schwinger model . the derivation of the force - gradient terms was given and the schwinger model was introduced . nested force - gradient schemes seem to be an optimal choice with relatively high convergence order and low computational effort . also it would be possible to improve the algorithm by measuring the poisson brackets of the shadow hamiltonian of the proposed integrator and then tuning the set of optimal parameters , e.  g. micro and macro step sizes . + in future work we will apply this approach to the hmc algorithm for numerical integration in lattice qcd . here we expect the adapted nested - force gradient scheme to outperform the original one , if we further partition the action into more than two parts , by using techniques to factorize the fermion determinant : less force - gradient information is needed for the most expensive action , and only leap - frog steps are needed for the high frequency parts of the action . this work is part of project b5 within the sfb / transregio 55 _ hadronenphysik mit gitter - qcd _ funded by dfg ( deutsche forschungsgemeinschaft ) . s.  duane , a.d . kennedy , b.j . pendleton , d.  roweth , hybrid monte carlo , phys . b195 ( 1987 ) , pp . e.  hairer , c.  lubich , g.  wanner , geometric numerical integration : structure - preserving algorithms for ordinary differential equations , springer , berlin , 2002 . omelyan , i.m . mryglod , r.  folk , symplectic analytically integrable decomposition algorithms : classification , derivation , and application to molecular dynamics , quantum and celestial mechanics , comput . 151 ( 2003 ) , pp .
Abstract:

Answer: we study a novel class of numerical integrators , the adapted nested force - gradient schemes , used within the molecular dynamics step of the hybrid monte carlo ( hmc ) algorithm . we test these methods in the schwinger model on the lattice , a well known benchmark problem . we derive the analytical basis of nested force - gradient type methods and demonstrate the advantage of the proposed approach , namely reduced computational costs compared with other numerical integration schemes in hmc .


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: graphene has entered the nanoelectronics scenario only recently and it is intensely investigated as one of the most promising candidates to replace silicon as a channel material in nanoscale transistors . even though graphene is a gapless material , a significant energy gap can be opened in different ways , like for example by `` rolling '' it in carbon nanotubes  @xcite or by defining narrow graphene stripes through electron beam lithography  @xcite or chemical reaction  @xcite . however , several unsolved technological problems arise : state of the art technology can not indeed conveniently control the chirality and the nanotube position , as well as define graphene nanoribbons with widths close to 1 - 2 nm with atomically flat edges , in order to obtain a reasonable semiconducting energy gap . recently , theoretical models  @xcite and experiments  @xcite have shown that bilayer graphene has an energy gap controllable by a vertical electric field . one could exploit this property to use bilayer graphene as a channel material for fets , defining an energy gap only when really needed , i.e. when the fet must be in the off state . this capability could open the possibility of patterning a bilayer graphene sheet with lithographic techniques and non prohibitive feature sizes .    in this letter , we want to assess the possibility of fabricating a tunable - gap bilayer graphene ( bg ) fet with independent gate operation , based on atomistic numerical simulations . we assume that all the technological challenges associated to the reliable fabrication of bilayer graphene fets can be solved , and evaluate the potential performance of near ideal device structures . in our opinion , this is one of the most powerful uses of computer simulations in guiding and orienting nanoelectronics research . a semi - analytical model based on the effective mass approximation has been presented in order to compare bilayer against monolayer graphene transistors  @xcite . however , in that work a constant induced energy bandgap is assumed for bilayer graphene , missing the most important material property of the material .    for a physically sound analysis of bg - fets , we have developed a code based on the non - equilibrium green s function formalism ( negf ) , with a tight - binding hamiltonian on a @xmath1 orbital basis set in the real space , which has been included in our in - house device simulator nanotcad vides  @xcite . we will show that bg - fets miss the itrs requirements  @xcite for the @xmath0 ratio ( @xmath2 ) by a large amount , since the induced gap is not sufficient to suppress band - to - band tunneling currents . our approach is based on the self - consistent solution of the poisson and schrdinger equations within the negf formalism and the ballistic transport assumption , as described in  @xcite . the considered bilayer graphene hamiltonian is composed by the two single layer graphene hamiltonians , coupled by the @xmath3=0.35  ev hopping parameters in correspondence of overlaying atoms along the @xmath4-direction  @xcite : the elementary cell is depicted in fig . semi - infinite contacts have been modeled along the @xmath5 direction by means of self - energies , while periodic boundary conditions have been imposed in the @xmath6 direction , with period equal to @xmath7 , where @xmath8  nm is the carbon - carbon bonding distance . the simulated device is a double - gate bg - fet , whose structure is shown in fig . we assume metal gates , and a 1.5  nm layer of sio@xmath9 as gate dielectric . we also assume an air gap of 0.5  nm between the dielectric interface and the position of carbon sites , as evaluated in  @xcite . the channel is 15  nm long , and the inter - layer distance is 0.35  nm . the source and drain extensions are 10  nm long , and are doped with an equivalent molar fraction of fully ionized donors @xmath10 . particular attention has to be posed in the computation of the mobile charge . considering for the sake of simplicity the equilibrium case , the mobile charge on a given site on layer @xmath11 ( @xmath12 ) of the bilayer graphene sheet reads @xmath13 \label{rho}\end{aligned}\ ] ] where @xmath14 is the elementary charge , @xmath15 is the intrinsic ( mid - gap ) fermi level , @xmath16 is the fermi - dirac occupation factor , and @xmath17 is the local density of states on layer @xmath11 .    let us stress the fact that , in order to avoid an unphysical antiscreening behavior , it is necessary to include the effect of dielectric polarization along the direction perpendicular to the graphene bilayer , due to the response of valence electrons to the applied field . we include this effect by introducing a polarization factor @xmath18 defined as @xmath19 de } , \label{alpha}\ ] ] where @xmath20 is the potential difference between the two graphene layers . by using @xmath18 , we can write the total charge ( mobile charge , valence band electrons and ions ) on a given site of layers 1 and 2 as @xmath21 \nonumber \\ \rho^\prime_2&=&\rho_2-q[1 - 2\alpha(v ) ] .\end{aligned}\ ] ]    however , performing the integral from @xmath22 to @xmath15 in ( [ alpha ] ) at each iteration step can be too computationally demanding . to this purpose , we compute the polarization factor as a function of @xmath20 only once , for an infinitely long bg - fet ( solid line in fig .  [ fig1]c ) . as can be seen , @xmath23 can be reasonably approximated by a straight line @xmath24 . in our code , we have used an odd polynomial of 11th order obtained through least mean square fitting ( symbols in fig . [ fig1]c ) .    in fig . [ fig2 ] the transfer characteristics for @xmath25 are shown , for different gate configurations . in particular , in fig . [ fig2]a the transfer characteristics as a function of the top gate voltage @xmath26 for two bottom gate voltages ( @xmath27 and @xmath28 ) are shown . as can be seen , large currents can be obtained above threshold , while the bg - fet shows huge problems to switch off . even after applying really negative gate voltages , the drain - to - source current decreases only by roughly six times . the ratio even worsens when @xmath28 , since in this case a larger amount of positive charge is injected in the valence band , due to band - to - band tunneling , which pins the potential along the channel , and further degrades the gate control over the barrier . such results are close to those obtained in experiments on monolayer graphene devices  @xcite , showing that the tunable gap is too small to effectively suppress band - to - band tunneling in the off state .    in order to increase the vertical electric field between the two layers of graphene , we have driven the device by imposing a fixed voltage between the two gates @xmath29 .    even in this case , device transconductance ( i.e. the slope of the transfer characteristic ) almost remains the same , regardless of the applied @xmath30 . this can be explained by the large accumulation of positive charge in the channel , caused by band - to - band tunneling in correspondence of the drain  @xcite : this charge screens the vertical electric field , limiting the size of the opened bandgap .    in fig [ fig3]b , the transmission coefficient as a function of the transversal wave vector @xmath31 is shown for @xmath32 , @xmath25 and @xmath33 , i.e. , close to minimum achievable current , and in correspondence of the band minimum . different energy intervals are identified with letters and compared to a sketch of the band edge profile in fig . [ fig3]a as an aid to interpretation . as can be seen , energy intervals tagged as @xmath34 , @xmath35 and @xmath36 correspond to the bandgaps in the channel , source and drain extensions , respectively : the gap is not uniform , it is maximum in the channel , where it is almost 100 mev . one can also see the bound states in the valence band in the energy interval @xmath37 and the valence band states in @xmath38 , which are populated by holes for negative voltages , and which eventually pin the potential of the two mono - layers of graphene , i.e. reducing the gate effectiveness in controlling the channel barrier and the bandgap . band - to - band tunneling can in principle be reduced by engineering the electric field between source and channel . here we make two attempts , by introducing a lateral spacer on the gate sides of 2.5  nm , and by introducing a linear doping profile degrading to zero in 2.5  nm under the spacer . both options have a positive but limited effect , as shown in fig . 3c , and do not solve the problem . we have explored the possibility of realizing fets by exploiting the tunable - gap property of bilayer graphene , through a code based on tight - binding negf device simulations . aiming at technology foresight , we have considered the overly optimistic case of an ideal structure and of ballistic transport . our computer simulations show that tuning the gap is not a very promising technique for achieving appropriate switching characteristics of bg - fets , due to the too weak suppression of band - to - band tunneling . a more complete exploration of the design space , in terms of structure and bias , would be required to cast a definitive answer , that is beyond the scope of this letter . the work was supported in part by the ec seventh framework program under project grand ( contract 215752 ) , by the network of excellence nanosil ( contract 216171 ) , and by the european science foundation eurocores program fundamentals of nanoelectronics , through funds from cnr and the ec sixth framework program , under project dewint ( contracteras - ct-2003 - 980409 ) .    1 s.  iijima , `` helical microtubules of graphite carbon '' , _ nature _ , vol . 354 , pp . 56 - 58 , nov . 1991 . z.chen , y - m . lin , m. j. rooks and p. avouris , `` graphene nano - ribbon electronics '' , phys . 228 - 232 , 2007 . x. li , x. wang , l. zhang , s. lee , h. dai , `` chemically derived , ultrasmooth graphene nanoribbon semiconductors '' , science , vol . 1229 - 1232 , 2008 . j. nilsson , a.h . castro neto , f. guinea , n.m.r . peres , `` electronic properties of graphene multilayers '' , phys . 97 , p. 266801 t. ohta , a. bostwick , t. seyller , k. horn , e. rotenberg , `` controlling the electronic structure of bilayer graphene '' , science , vol . 951 - 954 , 2006 . oostinga , h.b . heersche , x. liu , a.f . morpurgo , `` gate - induced insulating state in bilayer graphene devices '' , nature material , vol . 7 , pp.151 - 157 , 2008 . y. ouyang , p. campbell , j. guo , `` analysis of ballistic monolayer and bilayer graphene field - effect transistors '' , appl . 92 , p. 063120 g. fiori , g. iannaccone `` code for the 3d simulation of nanoscale semiconductor devices , including drift - diffusion and ballistic transport in 1d and 2d subbands , and 3d tunneling '' , _ journal of computational electronics _ 63 - 66 , 2005 . international technology roadmap for semiconductor 2007 , available : http://public.itrs.net . g. fiori , g. iannaccone and g. klimeck , `` a three - dimensional simulation study of the performance of carbon nanotube field - effect transistors with doped reservoirs and realistic geometry '' _ ieee trans . . dev . _ , vol . 53 , pp.1782 - 1788 , 2006 . m. c. lemme , t. j. echtermeyer , m. baus , and h. kurz , `` a graphene field - effect device '' , _ electr . 282 - 284 , 2007 . g. fiori , g. iannaccone , `` simulation of graphene nanoribbon field - effect transistors '' , _ ieee electr . 28 , issue 8 , pp . 760 - 762 , 2007 .
Abstract:

Answer: we explore the device potential of tunable - gap bilayer graphene fet exploiting the possibility of opening a bandgap in bilayer graphene by applying a vertical electric field via independent gate operation . we evaluate device behavior using atomistic simulations based on the self - consistent solution of the poisson and schrdinger equations within the negf formalism . we show that the concept works , but bandgap opening is not strong enough to suppress band - to - band tunneling in order to obtain a sufficiently large @xmath0 ratio for cmos device operation . * keyworks : * negf , graphene , bilayer , tight - binding hamiltonian , poisson / schrdinger .