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: it has been realised over the past decade that the supermassive black hole ( smbh ) at the centre of a galaxy bulge plays a major role in galaxy evolution . the process commonly considered is the negative feedback effect which may determine the final stellar mass of a galaxy so that it ends up with the locally observed universal bulge stellar mass over smbh mass ratio ( e.g. cattaneo et al . 2009 , fabian 2009 ) , m@xmath2/m@xmath3 700 ( 500 in marconi & hunt 2003 , mclure & dunlop 2001 or ferrarese et al . 2006 , 700 in kormendy & gebhardt 2001 and 830 in mclure & dunlop 2002 ) . indeed the correlation that exists between the mass of the central smbh of local galaxies and their bulge luminosity ( kormendy & richstone 1995 , magorrian et al . 1998 ) , stellar mass or velocity dispersion ( gebhardt et al . 2000 , ferrarese & merritt 2000 ) , suggests the existence of a physical mechanism connecting the activity of a galaxy nucleus and the bulding of its host galaxy stellar bulge . the most massive galaxies exhibiting the largest over - abundance of @xmath4-elements with respect to the solar value , this means that they formed their stars the fastest ( thomas et al . 2005 ) and indeed a downsizing of the typical star forming galaxy has been observed in deep surveys , i.e. the comoving star formation rate activity of the universe is dominated by more and more massive galaxies going back in lookback time . if negative feedback may provide a way to `` kill '' galaxies to reproduce the downsizing of galaxies , it does not necessarily explain why most stars in massive ellipticals formed rapidly . hence positive feedback must have been at play early on shortly after the formation of the proto - galaxy . major mergers have long been thought to be the best candidates to explain those rapid processes by first triggering a starburst then feeding the central smbh and activating the qso phase ( sanders et al . 1988 ) , but this scenario has been subject to debate , in particular since it was found that agn activity in the nearby universe was not shown to correlate with galaxy pairs ( li et al . 2008 ) . since the masses of smbh and of stars in galaxies are correlated , it is natural to consider the alternative possibility that both processes are related , with one possibly affecting the other . if radio jets are shown to be able to trigger star formation in galaxies , as suggested by early studies of distant galaxies but more or less abandoned since then , they may represent an interesting candidate process to understand the universal mass ratio since the most massive smbh will have the largest accretion rates , hence also most powerful radio jets , possibly triggering more efficiently star formation . because radio jets are difficult to identify , especially in the distant universe , they have not yet been considered as a major actor in the process of galaxy formation . however , they must have been ubiquitous in the past to explain the ubiquity of smbh in the center of local galaxies . together with the following pieces of evidence , this suggests that their role may have been underestimated until now : ( i ) the discovery of a jet - induced ultra - luminous infrared galaxy ( ulirg ) , one of the strongest starbursts ever observed , in the system he0450@xmath02958 , ( ii ) the presence of massive extended emission line regions ( eelr ) surrounding quasars and mostly radio galaxies , ( iii ) evidence that radio jets were more abundant in the past , ( iv ) the presence of large concentrations of molecular gas with an offset with respect to their neighboring qso , possibly resulting from the impact of radio jets , ( v ) theoretical arguments suggesting that radio jets may trigger star formation efficiently . our aim here is not to firmly state that we are convinced that galaxy formation can not be understood without accounting for the role of radio jets , since existing data and models are too sparse to either infer or reject such statement . however , in the near future , new observatories such as alma , evla , e - elt or the jwst will provide crucial observations that should allow us to better address this issue , hence we wish here to bring back to the front this mechanism as a possible driver not only of star formation , but maybe even galaxy formation . evidence for jet - induced star formation has been found in various objects and environments , either far away from the host galaxy , such as in the lobe of radio jets ( e.g. van breugel et al . 1985 ) , or inside the host galaxy , resulting in the so - called radio - optical alignment ( e.g. mccarthy et al . 1987 , rees 1989 ) . mccarthy et al . ( 1987 ) noticed that both the stellar continuum and size of the emission line region of 3cr radio galaxies at @xmath50.6 were highly elongated in parallel with their radio jets . this so - called radio  optical alignment was interpreted as evidence that the radio jets interact with the interstellar medium and stimulate large - scale star formation in the host galaxy ( see also rees 1989 , rejkuba et al . 2002 , oosterloo & morganti 2005 ) . other mechanisms than jet - induced star formation have been suggested in the literature to explain the radio - optical alignment effect such as the scattering of light from the central agn ( e.g. dey & spinrad 1996 ) or the nebular continuum emission from warm line - emitting regions ( dickson et al . 1995 ) . however , many of these objects show clear evidence of star formation . this is , in particular , the case of 4c 41.17 ( @xmath1=3.8 ) , which rest - frame uv continuum emission is aligned with the radio axis of the galaxy , unpolarized and showing p cygni - like features similar to those seen in star - forming galaxies ( dey et al . the most dramatic events , with star formation rates as high as 1000 m@xmath6 yr@xmath7 , are associated with very luminous radio galaxies at redshifts up to z @xmath8 4 often located at the center of proto - clusters of galaxies ( dey et al . 1997 , bicknell et al . 2000 , zirm et al . 2005 ) . even in the closest radio galaxy centaurus a , the reality of jet - induced star formation has been the center of a debate . early on , blanco et al . ( 1975 ) noticed the presence of extended gaseous filaments with bright knots exhibiting strong h@xmath4 and uv emission or loose chains of blue compact objects aligned with the radio jets of the central active nucleus . the origin of these optical filaments has been subject to a long debate regarding the main source of their ionization until the young and massive stars inhabiting these optical knots were individually resolved showing that jet - induced star formation must indeed have taken place in this system ( rejkuba et al . if there is indeed a mechanism through which radio jets may induce star formation in quasars or radio galaxies , then this mechanism may also be at play inside their scaled - down version , i.e. microquasars . the relativistic jets of the microquasar grs 1915@xmath9105 , which show apparent superluminal motion , are aligned with two iras sources themselves associated with radio knots . one of the radio source was resolved with the vla showing a non thermal extension pointing also in the direction of the microquasar on one side and the iras source on the other side ( see fig .  [ fig : grs1915 ] ) . both iras sources appear to be h ii regions ionized by late o or early b stars ( rodriguez & mirabel 1998 ) .        inversely , radio jets are known to be efficient in producing a negative feedback on the environment of the massive radio galaxies located at the center of the so - called `` cooling flow '' clusters ( see fabian 2009 , cattaneo 2009 ) . there , even when its cooling time gets lower than 1 gyr , the intracluster medium ( icm ) temperature does not fall below a given threshold which implies that non gravitational heating is taking place . the presence of buoyant radio bubbles brings direct evidence that on these large scales , radio jets can prevent cooling from happening if they are long - lived . at the smaller scale of individual galaxies , there is no such strong direct evidence for negative feedback and the process must be different since radio jets are collimated , hence should escape galaxies without affecting the bulk of their ism . for this reason , a second mode has been considered , called the radiative or wind mode ( see discussion in fabian 2009 ) . yet observationally , active nuclei are often found associated with star formation and both effects may not be exclusive in the sense that positive feedback may occur first at early stages in galaxy evolution followed by the quenching of star formation at a later stage .    the discovery of a case of a ulirg in which the intense star formation activity appears to be triggered by the impact of a radio jet ( elbaz et al . 2009 and sect . [ sec : he0450 ] ) suggests that , at least in some cases , positive feedback from radio jets can be a powerful mechanism . this is a scaled up version of the prototypical example for jet - induced star formation , minkowski s object , and in both cases , no evidence has been found for the presence of stars older than the starburst which suggests that the process may be responsible for the triggering of the formation of a whole galaxy . the fact that at @xmath1@xmath82 , the comoving space density of radio galaxies with powerful jets was 1000 times higher than at the present epoch ( willott et al . 2001 ) , suggests that the process taking place in these locally rare cases may have been an important actor in the process of galaxy formation and evolution . relic traces of this activity may be seen in the presence of the extended emission line regions ( eelrs ) surrounding qsos , and radio galaxies , with masses which may be as large as that of a massive galaxy , i.e. 10@xmath10 m@xmath6 ( fu & stockton 2008 ) . those gas reservoirs may represent the seeds for new galaxies if they are shocked by a radio jet as it may have happened in he0450@xmath02958 and such events may have been more common in the past . indeed , hydrodynamical simulations of radiative shock - cloud interactions indicate that for moderate gas cloud densities ( @xmath11 1 @xmath12 ) such as those observed in eelrs , cooling processes can be highly efficient and result in more than 50% of the initial cloud mass cooling to below 100 k , hence leading to star formation ( fragile et al . the mere existence of these eelrs itself is suggestive of the impact of radio jets on galaxies since their abundance ratios are similar to those observed in the host galaxies of their neighboring agn , which fu & stockton ( 2008 ) interpret as evidence that they were expelled out of the host galaxy by the radio jets . because they are collimated , radio jets may not affect very strongly the host galaxy of the active nucleus that produces them but they may produce an impact on the local environment . we present here two elements favoring this interpretation :    - the pair of ulirgs in the he0450@xmath02958 system , discussed in the next section , provide evidence for a major event in the formation process of a galaxy related to radio jets ( see sect . [ sec : he0450 ] ) . - several agns have been found with neighboring molecular gas clouds , which may be the result of the impact of a radio jet on its environment ( see sect .  [ sec : co ] ) . he0450@xmath02958 is a nearby luminous ( @xmath13@xmath14@xmath15 ) radio quiet quasar located at a redshift of @xmath1=0.2863 ( canalizo & stockton 2001 ) and presently the only known quasar for which no sign of a host galaxy has been found ( magain et al . this object has been the center of a debate and various scenarios have been proposed to explain the absence of detection of a host galaxy . since dust extinction could potentially explain it , we imaged he0450@xmath02658 with visir , the very large telescope imager and spectrometer in the infrared ( visir ) at the eso - vlt . indeed , a bright source was detected by iras at the approximate location of he0450@xmath02958 within the position error bar of 5 arcmin which , at the luminosity distance of the qso ( @xmath1=0.2863 ) , would translate into the total infrared luminosity of an ultra - luminous infrared galaxy ( ulirg , l@xmath16=l(81000@xmath17m)@xmath1810@xmath19 l@xmath6 , de grijp et al . 1987 , low et al . 1988 ) . in a first analysis , a unique mid - ir source was detected with visir associated with the qso , but combined with hst  nicmos imaging at 1.6@xmath17 m , it was not possible to disentangle a local dust torus from a whole dust - obscured galaxy ( jahnke et al . 2009 ) . however , a re - analysis of the visir data led to the detection of a second mid - ir source in the field , associated with the 7 kpc distant companion galaxy to the qso ( elbaz et al . 2009 ) . 2958 from 0.3 to 400@xmath17 m in the rest - frame ( red line ) . the plain green and blue lines represent the median and 68 percentile reddest seds from the atlas of 47 quasar seds of elvis et al . ( 1994 ) , normalized to the visir measurement at 11.3@xmath17 m ( observed ; 8.9@xmath17 m rest - frame ) associated to the qso ( filled blue circle ) . filled blue star for the companion galaxy . * _ ( b ) _ * visir image with contours from nicmos - f160w at 1.6@xmath17 m ( in yellow , 1.2@xmath17 m in the rest - frame , jahnke et al . 2009 ) and atca co 10 ( dark blue lines , from papadopoulos et al . 2008 ) . * _ ( c ) _ * atca 6208 mhz radio continuum ( dark blue lines , from feain et al . 2007 ) emission associated with the radio jets overlayed on the visir image.,title="fig:",width=211 ] 2958 from 0.3 to 400@xmath17 m in the rest - frame ( red line ) . the plain green and blue lines represent the median and 68 percentile reddest seds from the atlas of 47 quasar seds of elvis et al . ( 1994 ) , normalized to the visir measurement at 11.3@xmath17 m ( observed ; 8.9@xmath17 m rest - frame ) associated to the qso ( filled blue circle ) . filled blue star for the companion galaxy . * _ ( b ) _ * visir image with contours from nicmos - f160w at 1.6@xmath17 m ( in yellow , 1.2@xmath17 m in the rest - frame , jahnke et al . 2009 ) and atca co 10 ( dark blue lines , from papadopoulos et al . * _ ( c ) _ * atca 6208 mhz radio continuum ( dark blue lines , from feain et al . 2007 ) emission associated with the radio jets overlayed on the visir image.,title="fig:",width=134 ] 2958 from 0.3 to 400@xmath17 m in the rest - frame ( red line ) . the plain green and blue lines represent the median and 68 percentile reddest seds from the atlas of 47 quasar seds of elvis et al . ( 1994 ) , normalized to the visir measurement at 11.3@xmath17 m ( observed ; 8.9@xmath17 m rest - frame ) associated to the qso ( filled blue circle ) . filled blue star for the companion galaxy . * _ ( b ) _ * visir image with contours from nicmos - f160w at 1.6@xmath17 m ( in yellow , 1.2@xmath17 m in the rest - frame , jahnke et al . 2009 ) and atca co 10 ( dark blue lines , from papadopoulos et al . * _ ( c ) _ * atca 6208 mhz radio continuum ( dark blue lines , from feain et al . 2007 ) emission associated with the radio jets overlayed on the visir image.,title="fig:",width=134 ]    the visir image provides the first direct evidence that the iras `` source '' is in fact made of two well separated objects , the qso and its companion galaxy . the infrared sed of the `` system '' ( fig . [ fig : im]a ) peaks around @xmath850@xmath17 m ( rest - frame ) which indicates that the bulk of the far infrared light arises from star formation ( the qso and companion galaxy are both contained in the single iras measurement ) and not from and active nucleus that would drop beyond @xmath820@xmath17 m ( netzer et al . 2007 , minimum temperature of @xmath8200 k typically ) . the ir sed of the companion galaxy , which produces the bulk of the far - ir light , is very well fitted with a standard star forming template sed ( red plain line , from the library of chary & elbaz 2001 ) . the total infrared luminosity of the companion galaxy is l@xmath16@xmath82@xmath2010@xmath19 l@xmath6 , which corresponds to a sfr@xmath8340 m@xmath6yr@xmath7 ( conversion factor of kennicutt 1998 ) . hence he0450@xmath02958 appears to be a composite system made of a pair of ulirgs ( both galaxies have similar l@xmath16 ) , where both mechanisms , star formation and qso activity , are spatially separated by 7 kpc in two distinct sites . the fact that the starburst is associated with the companion galaxy is reinforced by the finding that the molecular gas traced by the co molecule appears to avoid the qso but to peak at the location of the companion galaxy ( dark blue contours in fig . [ fig : im]b ) close to the peak mid - infrared emission measured with visir ( papadopoulos et al . 2008 ) . this both suggests that there is a large amount of dust in this galaxy and that there is a large gas reservoir ( m(h@xmath21)@xmath222.3@xmath2010@xmath10 m@xmath6 ) to fuel an intense star formation event . no evidence for the presence of an old stellar population is found in the vlt - fors optical spectrum of the companion , although the radiation of such stars could be diluted in that of the youngest population . the age of the dominant stellar population is about 40200 myr suggesting that the galaxy was recently born . the stellar mass of the galaxy ( m@xmath23@xmath22[56]@xmath2010@xmath10 m@xmath6 ) is consistent with this timescale for the sfr derived from its ir luminosity . the companion galaxy coincides with one of the two radio lobes on both sides of the qso he0450@xmath02958 ( fig . [ fig : im]c ) , suggesting that it is hit by a radio jet emitted by the qso . evidence for shock excitation of the ism in the companion galaxy , as traced by an high [ nii]/h@xmath4 emission line ratio ( letawe et al . 2008b ) , confirms that the association of the jet and the galaxy is physical and not due to a projection effect . the visir image itself shows the presence of a bridge of mid - ir emission between the qso and the companion fig . [ fig : jets]-left . . altogether , this suggests that it is the radio jet that is triggering the starburst in the companion . as a result , he0450@xmath02958 is comparable to minkowski s object , a proto - typical case for jet - induced star formation ( van breugel et al . minkowski s object is a newly formed galaxy , with a stellar population age of only 7.5 myr and a stellar mass of 1.9@xmath2010@xmath24 m@xmath6 ( croft et al . 2006 ) , aligned with the radio jet of ngc 541 , an fr i galaxy ( fanaroff & riley 1974 ) located in the cluster abell 194 . the companion galaxy of he0450@xmath02958 might be a scaled - up version of minkowski s object , i.e. a massive galaxy whose formation was jet - induced . 2958 system with visir 11.3@xmath17 m image ( blue ) and hst in the r - band ( yellow ) . * _ ( right ) _ * composite image showing the details of minkowski s object ( croft et al . 2006 ) . optical image of the field of abell 194 , overlaid with radio continuum ( purple ) , hi data ( dark blue ) , and h@xmath4 data ( light blue).,title="fig:",width=190 ] 2958 system with visir 11.3@xmath17 m image ( blue ) and hst in the r - band ( yellow ) . * _ ( right ) _ * composite image showing the details of minkowski s object ( croft et al . optical image of the field of abell 194 , overlaid with radio continuum ( purple ) , hi data ( dark blue ) , and h@xmath4 data ( light blue).,title="fig:",width=289 ]    evidence for the presence of large amounts of matter in the surrounding of radio quasars , in the form of massive gas clouds exhibiting strong emission lines with high excitation levels , has been gathered since the early 60s . these extended emission lines regions ( eelrs , see stockton & mackenty 1983 ) can reach masses of several 10@xmath10 m@xmath6 ( fu & stockton 2008 , stockton et al . 2007 ) . spectro - imaging and narrow - band filter imaging in the fields of radio quasars demonstrated that eelrs were made of clouds of gas , whose strong emission lines were excited by the hard radiation of the quasar . fu & stockton ( 2008 ) found that low metallicity quasar host galaxies happened to be surrounded by eelrs with similar abundance ratios , which they interpreted as evidence that the eelrs had been expelled from the host galaxy itself . the double - lobed morphologies of extragalactic radio sources shows that the relativistic jets not only couple strongly with the ism of the host galaxies , but are capable of projecting their power on the scales of galaxy haloes and clusters of galaxies ( up to @xmath81 mpc , tadhunter 2007 ) . several eelrs were found in the surroundings of he0450@xmath02958 following the axis of the radio jets . they have been detected up to 30 kpc away from the qso ( letawe et al . the n - w blob of gas , photoionized by the radiation of he0450@xmath02958 , is itself located in the direction of the radio jet opposite to the one pointing at the companion galaxy . the thermodynamic state of eelrs is uncertain since the densities one can infer from their optical spectra are of the order of n@xmath25=100300 @xmath12 and t@xmath25=10,00015,000 k ( fu & stockton 2008 ) , which imply that they should be gravitationally unstable and have already started forming stars . hence fu & stockton ( 2008 ) modeled these regions using two components , a low density component , where most of the mass is locked at densities of order n@xmath261 @xmath12 , and a higher density component responsible for the observed emission lines . hydrodynamical simulations of radiative shock - cloud interactions indicate that for moderate gas cloud densities ( @xmath11 1 @xmath12 ) , cooling processes can be highly efficient and result in more than 50% of the initial cloud mass cooling to below 100 k ( fragile et al . hence , the companion galaxy of he0450@xmath02958 could have formed from a seed eelr that was hit by the radio jet . this process could explain why radio jets do not always produce the same effect . here we search for other cases such as he0450@xmath02958 , where a large mass of molecular gas is found offset with respect to the qso and associated with newly formed stars which formation could be induced by a radio jet . such an offset between molecular gas and qsos may be a common feature in distant radio sources as suggested by klamer et al . ( 2004 ) who produced a systematic search in @xmath1@xmath113 co emitters for an agn offset with respect to the co concentrations and for a connection with radio jets . out of the 12 @xmath27 3 co emitters that they studied , six showed evidence of jets aligned with either the co or dust , five have radio luminosities above 10@xmath28 w hz@xmath7 and are clearly agns , and a further four have radio luminosities above 10@xmath29 w hz@xmath7 indicating either extreme starbursts or possible agns . in the following , we discuss three proto - typical cases classified with increasing redshift from @xmath1@xmath142.6 to 6.4 that may be compared with he0450@xmath02958 .    at @xmath1=1.574 , 3c18 is radio loud quasar with recently formed radio jets , as inferred from their small physical size . a large mass of molecular gas ( @xmath30@xmath8(3.0@xmath310.6)@xmath2010@xmath10 m@xmath6 ) , inferred from co 21 emission , is found associated to the quasar with a positional offset of @xmath820 kpc ( willott et al . 2007 ) .    at @xmath1= 2.6 , txs0828@xmath9193 is a radio galaxy with a neighboring co gas concentration of @xmath81.4@xmath2010@xmath10 m@xmath6 located 80 kpc away and shows no evidence for an underlying presence of stars or galaxy ( nesvadba et al . an upper limit of 0.1 mjy at 24@xmath17 m , from the mips camera onboard spitzer , was obtained by nesvadba et al . ( 2009 ) , who concluded from this limit that no major starburst with more than several hundred solar masses per year could be taking place associated with this co concentration . we used the library of sed templates of local galaxies from chary & elbaz ( 2001 ) to convert this mid infrared measurement into an upper limit for the total ir luminosity at this position of l@xmath32@xmath87.8@xmath2010@xmath19 l@xmath6 , which would translate to a maximum sfr of @xmath81340 m@xmath6 yr@xmath7 using the kennicutt ( 1998 ) conversion factor for a salpeter imf . hence there is still room for a large amount of star formation in this object that may be in a similar stage than he0450@xmath02958 but at a much larger distance .    at @xmath1= 4.695 , br 1202@xmath00725 is the first high - redshift quasar for which large amounts of molecular gas was detected ( omont et al . the co map presents two well separated emission peaks which coincide with radio continuum emission interpreted as evidence for the presence of radio jets by carilli et al . the radio jets themselves would be too faint to be detected at the sensitivity level of the radio image but marginal evidence for variability suggest that the radio emission is not due to star formation . omont et al . ( 1996 ) suggested that the double co emission might be due to gravitational lensing , but no optical counterpart of the qso is found associated with the second source and the co ( 21 ) line profiles are different for the two components ( carilli et al . klamer et al . ( 2004 ) suggested that stars may have formed along the radio jets , providing both the metals and dust for cooling and `` conventional '' star formation to take place afterwards .    in the case of he0450@xmath02958 , the mid infrared emission is spatially associated with the offset co concentration suggesting that similar systems might be found not only through co imaging but also infrared imaging . sdss160705 + 533558 , located at a redshift of @xmath1=3.65 , might be a distant analog of he0450@xmath02958 in that respect since it also presents a positional offset between the maximum sub - millimeter ( from the submillimeter arra , sma ) and optical emission ( clements et al . 2009 ) .    a key test for the role of radio jets in the formation of galaxies will be the detection of either molecular gas or mid infrared emission with a positional offset with respect to their neighboring qso . this test will be fulfilled with the advent of alma , the jwst or project instruments such as metis ( mid infrared e - elt imager and spectrograph ) for the elt ( extremely large telescope ) . bicknell , g.v . , sutherland , r.s . , van breugel , w.j.m . 2000 , apj 540 , 678 blanco , v.m . , graham , j.a . , lasker , b.m . , osmer , p.s . 1975 , apj 198 , l63 canalizo , g. , stockton , a. 2001 , apj 555 , 719 carilli , c. l. , kohno , k. , kawabe , r. , et al . 2002 , aj 123 , 1838 cattaneo , a. , faber , s. m. , binney , j. et al . 2009 , nature 460 , 213 chary , r.r . , elbaz , d. 2001 , apj 556 , 562 clements , d. l. , petitpas , g. , farrah , d. et al . 2009 , apj 698 , l188 croft , s. , van breugel , w. , de vries , w. et al . 2006 , apj 647 , 1040 de grijp , m. h. k. , lub , j. , miley , g. k. 1987 , a&as 70 , 95 dey , a. , spinrad , h. 1996 , apj 459 , 133 dey , a. , van breugel , w. , vacca , w.d . , antonucci 1997 , apj 490 , 698 dickson , r. , tadhunter , c. , shaw , m. , clark , n. , morganti , r. 1995 , mnras 237 , l29 elbaz , d. , jahnke , k. , pantin , e. , le borgne , d. , letawe , g. , a&a 507 , 1359 elvis , m. , wilkes , b.j . , mcdowell , j.c . , et al . 1994 , apjss 95 , 1 fabian 2009 , these proceedings of iau symposium 267 ( arxiv:0912.0880 ) fanarof , b.l . , riley , j.m . 1974 , mnras 167 , 31 feain , i.j . , papadopoulos , p.p . , ekers , r.d . , middelberg , e. 2007 , apj 662 , 872 ferrarese , l. , merritt , d. 2000 , apj 539 , l9 ferrarese , l. , ct , p. , dalla bont , e. et al . 2006 , apj 644 , 21 fragile , p.c . , murray , s.d . , anninos , p. , van breugel , w. 2004 , apj 604 , 74 fu , h. , stockton , a. 2008 , apj 677 , 79 gebhardt , k. , bender , r. , bower , g. , et al . 2000 , apj 539 , l13 jahnke , k. , elbaz , d. , pantin , e. , et al . 2009 , apj 700 , 1820 kennicutt , r.c.jr 1998 , araa 36 , 189 klamer , i.j . , ekers , r.d . , sadler , e.m . , hunstead , r.w . 2004 , apj 612 , l100 kormendy , j. , richstone , d. 1995 , araa 33 , 581 kormendy , j. , gebhardt , k. 2001 , in aip conf.proc.586 , relativistic astrophysics , ed . h. martel & j. c. wheeler , p.363 letawe , g. , magain , p. , courbin , f. 2008a , a&a 480 , 69 letawe , y. , magain , p. , letawe , g. , courbin , f. , hutsemkers , d. 2008b , apj 679 , 983 li , c. , kauffmann , g. , heckman , t.m . , white , s.d.m . , ping , y.p . 2008 , mnras 385 , 1915 low , f.j . , cutri , r.m . , huchra , j.p . , kleinmann , s.g . 1988 , apj 340 , l1 magain , p. , letawe , g. , courbin , f. , et al . 2005 , nature 437 , 381 magorrian , j. , tremaine , s. , richstone , d. , et al . 1998 , aj 115 , 2285 marconi , a. , hunt , l.k . 2003 , apj 589 , l21 mccarthy , p.j . , van breugel , w. , spinrad , h. , djorgovski , s. 1987 , apj 321 , l29 mclure , r.j . , dunlop , j.s . 2001 , mnras 327 , 199 mclure , r.j . , dunlop , j.s . 2002 , mnras 331 , 795 nesvadba , n.p.h . , neri , r. , de breuck , c. , et al . 2009 , mnras 395 , l16 omont , a. , petitjean , p. , guilloteau , s. , et al . 1996 , nature 382 , 428 ooosterloo , t.a . , morganti , r. 2005 , a&a 429 , 469 papadopoulos , p.p . , feain , i.j . , wagg , j. , wilner , d.j . 2008 , apj 684 , 845 rees , m.j . 1989 , mnras 239 , 1 rejkuba , m. , minniti , d. , courbin , f. , silva , d.r . 2002 , apj 564 , 688 rodriguez , l.f . , mirabel , i.f . 1998 , a&a 340 , l47 sanders , d. , soifer , b.t . , elias , j.h . 1988 , apj 325 , 74 stockton , a. , mackenty , j.w . 1983 , nature 305 , 678 tadhunter , c. 2007 , newar 51 , 153 thomas , d. , maraston , c. , bender , r. , mendes de oliveira , c. 2005 , apj 621 , 673 van breugel , w. , filippenko , a.v . , heckman , t. , miley , g. 1985 , apj 293 , 83 willott , c. j. , rawlings , s. , blundell , k. m. , lacy , m. , eales , s. a. 2001 , mnras 322 , 536 willott , c. j. , martnez - sangre , a. , rawlings , s. 2007 , apj 133 , 564 zirm , a.w . , overzier , r.a . , miley , g.k . , et al . 2005 , apj 630 , 68
Abstract:

Answer: we discuss evidence that quasars , and more generally radio jets , may have played an active role in the formation stage of galaxies by inducing star formation , i.e. through positive feedback . this mechanism first proposed in the 70 s has been considered as anecdotic until now , contrary to the opposite effect that is generally put forward , the quenching of star formation in massive galaxies to explain the galaxy bimodality , downsizing and the universal black hole mass over bulge stellar mass ratio . this suggestion is based on the recent discovery of an ultra - luminous infrared galaxies , i.e. an extreme starburst , which appears to be triggered by a radio jet from the qso he0450@xmath02958 at @xmath1=0.2863 , together with the finding in several systems of an offset between molecular gas and quasars , which may be explained by the positive feedback effect of radio jets on their local environment .


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 plasmonics has recently emerged as a promising platform for studying strong light - matter interactions.@xcite graphene is a novel two - dimensional ( 2d ) material @xcite with unique electronic and optical properties @xcite and a wide range of applications such as biosensors,@xcite ultrafast lasers , @xcite optical modulators , @xcite and photodetectors . @xcite clean graphene samples are characterized by long electron scattering times and much lower , compared to metals , ohmic losses due to relatively weak electron - phonon interaction . @xcite charge - doped graphene with fermi energy , @xmath0 , in the range 0.2 - 0.6 ev exhibits a stable in - plane plasmon in the infrared frequency range @xcite with gate - tunable wavelength , @xmath1 , well below radiation wavelength at the same frequency . @xcite    a large local density of states ( ldos ) of graphene plasmons as compared to that of surface plasmon polaritons ( spp ) on metal surfaces ensures very efficient plasmon excitation by a local probe such as an excited dye molecule or an exciton in semiconductor quantum dot ( qd ) situated at a close distance to the graphene sheet . @xcite recent measurements @xcite of nonradiative energy transfer from a dye molecule to graphene pointed to a strong coulomb coupling between electronic excitations in these two systems with transfer rates exceeding the radiative decay rate by a factor of @xmath210@xmath3 . plasmons can also provide efficient coupling between several emitters situated near a graphene sheet , e.g. , by enhancing superradiance from two excited qds @xcite or facilitating long - distance energy transfer between a donor and an acceptor . @xcite these studies explored weak coupling regime between graphene plasmons and excitons , i.e. , the energy spectra of these excitations are not significantly altered by the optical interactions between them .    here we demonstrate that a _ strong _ coupling between excitons and graphene plasmons can be realized in a hybrid graphene - semiconductor system giving rise to a propagating exciton - plasmon state with a mixed 2d dispersion ( hereafter referred to as exciton - plasmariton@xcite ) . strong exciton - plasmon coupling effects have been recently studied between surface plasmon - polaritons ( spp ) or localized surface plasmons in metal structures and excitons in dye j - aggregates , @xcite individual dye molecules , @xcite semiconductor quantum wells ( qw ) , @xcite or qds and nanoscrystals . @xcite in these systems , very large rabi splittings ( relative to vacuum rabi splitting @xcite ) were reported in the range of 100 - 200 mev for dispersion of spp coupled to molecular excitons @xcite and in the range of 200 - 450 mev for absorption / transmission spectra of excitons coupled to localized plasmons ( plexcitons ) , @xcite while relatively small splittings ( @xmath4 mev ) were observed for spp coupled to qw excitons.@xcite the optical excitation energies in the above experiments were in the range 1 - 3 ev to match those of plasmon resonances in metal nanostructures . in contrast , stable graphene plasmons have significantly lower energies ( below @xmath5 ev in highly doped samples ) and , to the best of our knowledge , strong coupling regime in graphene has not yet been explored . here we show that a strong exciton - graphene plasmon coupling can be achieved in a hybrid structure comprised of doped graphene monolayer separated by a thin spacer from a narrow gap semiconductor qw . the proposed structure is schematically shown in fig . [ fig : fig1](a ) . the tunneling barrier between qw and graphene , spacer i , prevents a photoinduced charge transfer from qw to graphene . in this case , quasiparticles in graphene ( plasmons ) and qw ( excitons ) interact with each other only via direct coulomb coupling . the qw rests on a dielectric substrate separated by spacer ii . as examples , we chose four narrow gap semiconductor qws : inas , insb , and hgcdte ( hct ) with two different cd concentrations corresponding to hct bulk bandgaps of @xmath6 ev and 0.3 ev . @xcite in fig . [ fig : fig1](b ) , the exciton dispersion relations for all qws are plotted together with the plasmon dispersion . in principle , both exciton and plasmon energies as well as the coulomb coupling between them are subject to complex dielectric screening determined by spacer layer widths @xmath7 and @xmath8 , qw width @xmath9 , and dielectric constants of all the materials involved . for example , by varying the parameters of the layered structure , the coulomb coupling strength can be tuned in a wide range . in this paper , we assume a simplified model of dielectric screening by introducing a single effective dielectric constant , @xmath10 , and using it as a free parameter . for example , fig . [ fig : fig1](b ) shows the graphene plasmon and qw exciton dispersions for effective dielectric constant of @xmath11 and the graphene charge - doping level of @xmath12 ev . at this doping level , graphene plasmon is landau damped at @xmath13 nm@xmath14 , the onset of damping marked by red cross in fig . [ fig : fig1](b ) . it can be seen that , for the chosen semiconductor qws , the plasmon and excitons dispersion curves intersect at @xmath15 well below the damping onset . at smaller @xmath10 , exciton energies decrease ( due to stronger electron - hole interaction ) and plasmon energies increase , resulting in the exciton - plasmon resonance at even smaller @xmath15 . for small separations between graphene and qw , we thus expect a mixed exciton - plasmon state to form due to strong coulomb coupling between its constituent excitations . we show that the dynamics of exciton - plasmaritons can be accurately described by a simple hamiltonian similar to exciton - polariton hamiltonian with exciton - plasmon coupling constant strongly affected by the system parameters . we explore the strong coupling regime by performing numerical calculations involving full graphene response functions and compare the results to our model . the calculated rabi splitting between the exciton - plasmariton upper and lower dispersion branches strongly depends on system parameters and can reach values @xmath250 - 100 mev for small graphene - qw separation @xmath7 . the paper is organized as follows . in sec . [ sec : hamiltonian ] we present our model for coupled qw exciton and graphene plasmon and derive the exciton - plasmariton hamiltonian and its energy spectrum within plasmon - pole approximation . in sec . [ sec : numerics ] we perform numerical calculations of exciton - plasmariton energy spectrum using graphene full response functions and compare the results to those of our model . technical details of calculations are described in appendices . ( a ) schematics of the device . ( b ) plasmon and exciton dispersion relations are shown by thick turquoise and thin lines , respectively . the bulk bandgap of hct is given in parentheses in the legend . the qw width and the fermi energy in graphene are @xmath16 nm and @xmath12 ev , respectively . plasmon - exciton resonances and the onset of landau damping are marked by black and red crosses , respectively . graphene electron - hole transition continua are shown in grey . , width=288 ] we consider a qw of width @xmath9 separated from a monolayer graphene by a sufficiently high potential barrier of thickness @xmath7 . the system hamiltonian is @xmath17 , where @xmath18 and @xmath19 are many - body hamiltonians for electronic excitations in qw and graphene , respectively , and @xmath20 describes the coulomb interaction between them . in the absence of free electron population in qw , the main contribution to @xmath21 comes from the interaction of interband polarization in qw with the electric field due to charge density fluctuations in graphene ( see appendix  [ app : ham ] ) @xmath22 where @xmath23 is the dipole transition operator between conduction ( @xmath24 ) and valence ( @xmath25 ) bands and @xmath26 are bloch functions ( @xmath27).@xcite the electric field @xmath28 acting on exciton in qw is generated by density fluctuations in graphene @xmath29,\ ] ] where @xmath30 is 3d coordinate , @xmath31 and @xmath32 being in - plane and out - of - plane coordinates , respectively , and integration is carried over graphene plane @xmath33 . the in - plane fourier transform of the coulomb potential , @xmath34 , is given by @xmath35 , where @xmath36 is the in - plane momentum and @xmath37 is the effective dielectric constant . electron density operator and the average electron density in graphene are denoted by @xmath38 and @xmath39 , respectively , . our goal is to derive effective exciton - plasmon hamiltonian . to this end , let us first recast @xmath21 in second - quantized electron - hole form . in a standard manner,@xcite the bloch function splits into envelope and periodic parts , @xmath40 , where the latter defines the interband dipole matrix element averaged over unit cell volume @xmath41 , @xmath42 , while the envelope functions @xmath43 describe `` free '' electrons in qw , where @xmath44 is the size - quantization wave function in qw and @xmath45 is the normalization area ( hereafter set to unity ) . performing 2d fourier transform , we obtain @xmath46 where @xmath47 , @xmath48 are electron ( annihilation ) and hole ( creation ) operators defined in a usual way through their conduction and valence band counterparts,@xcite and @xmath49 is the fourier transform of graphene density operator , @xmath50 being electron operator in graphene ( as usual , the point @xmath51 is excluded ) . the matrix element @xmath52 is evaluated as @xmath53u_{cv}(q),\ ] ] where ( for @xmath54 ) @xmath55 is overlap between gpaphene coulomb potential and qw excitations across the structure . we now introduce the exciton creation operator in a standard manner @xcite as @xmath56 where @xmath57 , @xmath58 ( @xmath59 ) are the relative electron and hole effective masses , and @xmath60 is the 2d fourier transform of eigenfunction @xmath61 of the wannier equation ( @xmath62 is the exciton center of mass momentum and @xmath63 is its quantum state ) . using the orthogonality relation @xmath64 , eq . ( [ eq : hint_eh ] ) takes the form @xmath65 where @xmath66 characterizes the strength of exciton coupling to charge density fluctuations in graphene , and @xmath67 and @xmath68 are , respectively , projections of @xmath69 onto @xmath32-axis and vector @xmath62 ( note that @xmath70 ) .    other terms of the hamiltonian that contain qw electron and hole operators can be expressed in excitonic variables using standard projection technique . the free exciton hamiltonian has the usual form @xcite @xmath71 , where @xmath72 is exciton dispersion ( see below ) while higher - order terms including exciton - exciton scattering and exciton scattering by graphene density fluctuations are negligible for low exciton densities we consider ( see appendix  [ app : ham ] ) . the hamiltonian of qw exciton interacting with graphene density fluctuations then takes the form @xmath73 the energy spectrum of exciton - plasmatiton can be obtained from exciton green function `` dressed '' due to the interactions with graphene density fluctuations . in the rotating wave approximation , the non - interacting exciton green function can be written as @xmath74 , where @xmath75 is free exciton damping . the interaction with graphene excitations results in the appearance of the exciton self - energy in the green function , @xmath76 where , in the lowest order in qw - graphene coupling , the spectrum is determined by the density - density correlation function in graphene , @xmath77 : @xmath78 . it is important to note here that even though the spectra of electronic excitations in graphene and qw are assumed isotropic by themselves , this is not necessarily the case for exciton - plasmaritons due to the possible anisotropy in the coupling constant , @xmath79 . this anisotropy originates from the directionality of @xmath69 and is , therefore , maximized when @xmath69 is parallel to the graphene plane [ see eq . ( [ eq : t_coupl ] ) ] .    within the random phase approximation ( rpa ) , the density - density correlation function in graphene is given by @xmath80 $ ] , where @xmath81 is the polarization bubble in graphene and @xmath82 .    in the following we restrict ourselves to the lowest @xmath83 excitonic state in a narrow qw ( @xmath84 ) characterized by wave function @xmath85 where @xmath86 is a 2d bohr radius . the exciton energy is @xmath87 , where @xmath88 is the excitation energy of the lowest qw exciton . here , @xmath89 is the bulk bandgap energy , @xmath90 is exciton reduced mass , @xmath91 is exciton binding energy , and the second term in the r.h.s . is correction to the bulk bandgap energy due to the quantum confinement in qw . for a symmetric qw , the transverse part of envelope functions of the lowest qw subband is @xmath92 . then the coulomb overlap @xmath93 can be explicitly evaluated as @xmath94 in the long - wave limit , the graphene polarization bubble has the form @xmath95 , so that @xmath77 is dominated by the plasmon pole @xmath96 where @xmath97 and @xmath98 are , respectively , plasmon energy dispersion and damping rate ( in energy units ) . combining eqs . ( [ eq : t_coupl ] ) and ( [ eq : pi_plpole ] ) , we obtain the exciton self - energy as @xmath99 where the exciton - plasmon coupling , @xmath100 , is given by @xmath101 with the above self - energy , the poles of eq . ( [ exciton - green ] ) determine two exciton - plasmariton energy branches , @xmath102,\ ] ] where @xmath103 and @xmath104 . to gain more insight , we note that eq . ( [ spectrum ] ) is the energy spectrum of a simple two - level hamiltonain where both the energies of the two levels , @xmath105 and @xmath106 , and the coupling between them , @xmath107 , are dependent on @xmath62 . this observation can be made precise by casting the graphene hamiltonian in terms of plasmon normal modes , @xmath108 , where @xmath109 is a plasmon creation operator . after expanding the graphene charge density operator in eq . ( [ ham - int ] ) over these normal modes , @xcite the hamiltonian ( [ ham2 ] ) of qw excitons interacting with graphene plasmons takes the form @xmath110\biggr ] , \end{aligned}\ ] ] where exciton - plasmon coupling @xmath111 is defined by eq . ( [ coupling ] ) . the imaginary part of @xmath105 and @xmath106 has to be interpreted in a usual manner as a decay constant of a metastable state . the corresponding lifetimes of an exciton and plasmon are given by @xmath112 and @xmath113 , respectively . the hamiltonian ( [ ham - final ] ) is similar to that describing exciton - polaritons and can be brought to canonical form by a standard bogolubov - hopfield transformation @xcite yielding , near the resonance , two dispersion branches eq . ( [ spectrum ] ) which , in general , are anisotropic in @xmath36-plane due to the dependence of exciton - plasmon coupling @xmath111 on the exciton polarization . note that if @xmath62 is parallel to the projection of @xmath69 onto the @xmath114-plane , then the coupling is maximized and @xmath115 , so that @xmath116 . this condition is assumed fulfilled hereafter . at resonance , @xmath117 , corresponding to momentum @xmath118 , the rabi frequency @xmath119 has the form @xmath120 and strongly depends on system parameters . all the analytical considerations above were based on the plasmon dispersion @xmath121 , eq . ( [ eq : pi_plpole ] ) . however , this simple dispersion relation becomes inaccurate in the high-@xmath10 dielectric environment due to non - local effects . it turns out that the dispersion relation could be analytically corrected to account for such effects ( see appendix  [ app : high - k ] for detail ) . below we compare eq . ( [ spectrum ] ) , augmented by this corrected dispersion relation , to full numerical calculations of the energy spectrum of the system . numerical calculations were performed for four narrow - gap semiconductors , insb ( @xmath122 ev ) , @xcite inas ( @xmath123 ev ) , @xcite and hgcdte ( hct ) with two different cd concentrations corresponding to @xmath6 ev and @xmath124 ev . @xcite the electron and hole effective masses for these materials are taken from refs .  . for each material , the dispersion of qw exciton intersects plasmon dispersion in graphene doped to @xmath125 ev before the landau damping onset ( see fig .  [ fig : fig1 ] ) . unless otherwise noted , full rpa response functions for graphene were employed in numerical calculations ( see appendix  [ app : polariz ] for detail ) and standard materials parameters were used . the interband matrix element @xmath126 was estimated as @xmath127 where @xmath128 is the kane energy , @xmath129 is the kane momentum matrix element and @xmath130 is the inter band transition momentum . for semiconductor materials considered in this work we have @xmath131 ev , @xmath132 ev and @xmath133 ev . @xcite     surface plot of the logarithm of spectral function for the exciton - plasmariton , @xmath134 , for ( a ) hct with @xmath135 ev and ( b ) inas . the effective dielectric constant is set to @xmath136 ( panel a ) and @xmath137 ( panel b ) . other system parameters are set to @xmath12 ev , @xmath138 nm and @xmath139 nm . to make the excitonic `` ridge '' broader and more visible , we substituted @xmath140 ps in this figure . , width=288 ]    figure  [ fig : fig2 ] shows 3d log plot of exciton spectral function @xmath141 ( in arbitrary units ) as defined by eq .  ( [ exciton - green ] ) for htc and inas qws of width @xmath16 nm , separated from the graphene by @xmath138 nm thickness spacer , and effective dielectric constants of the structure @xmath136 and @xmath137 , respectively . plasmon decays time is taken @xmath142 ps . @xcite exciton lifetime can vary in a wide range depending on specific parameters of the system , e.g. , qw thickness.@xcite here we assume a very conservative value of @xmath143 ps . since @xmath144 , the specific value for the exciton lifetime is largely irrelevant . with the chosen parameters , both exciton - plasmariton branches are seen to be well - separated in fig . [ fig : fig2 ] .    in figs . [ fig : fig3 ] and [ fig : fig4 ] we show density plot of @xmath145 vs. momentum @xmath15 for inas and htc qws for two different effective dielectric constants . thin black lines , superimposed on top of the density plots , show the dispersion relations of exciton - plasmariton dispersion branches calculated from eq .  ( [ spectrum ] ) . we see that with increasing @xmath10 , the splitting between exciton - plasmariton branches is strongly reduced while the resonance position is shifted towards higher @xmath15 . both effects can be traced to exciton - plasmon coupling which is directly affected by increase of dielectric constant , @xmath146 . at the same time , increasing @xmath10 reduces the exciton binding energy and hence shifts the resonance position to higher @xmath15 which , in turn , leads to a decrease of @xmath93 resulting in an even stronger reduction of @xmath147 . density plot of the spectral function for the exciton - plasmariton in inas . the effective dielectric constant is set to @xmath148 ( panel a ) and to @xmath137 ( panel b ) . other system parameters are set to @xmath12 ev , @xmath138 nm and @xmath139 nm . , width=288 ]     density plot of the spectral function for the exciton - plasmariton in hct ( @xmath135 ev ) . the effective dielectric constant is set to @xmath149 ( panel a ) and to @xmath136 ( panel b ) . other system parameters are set to @xmath12 ev , @xmath138 nm and @xmath139 nm . , width=288 ]    since the _ exciton _ spectral function is plotted , only the exciton - like branch is expected to be seen at vanishing exciton - plasmon coupling . indeed , the intensity of plasmon - like branches is seen to decrease for larger @xmath10 . as can be seen in all the panels , the analytical curves for the dispersion ( thin black lines ) agree reasonably well with the numerical results , seen as ridges in density plots . therefore , we will rely on analytical model in our further analysis . strength of exciton - plasmon coupling versus ( a ) effective dielectric constant , ( b ) fermi energy , ( c ) thickness of spacer i and ( d ) thickness of the quantum well . the system parameters are set to @xmath136 , @xmath12 ev , @xmath138 nm and @xmath16 nm except for that varied in each panel ( e.g. , @xmath10 in panel a ) . red crosses mark the onset of landau damping and dashed lines correspond to short - living plasmons that live within the electron - hole sector and , therefore , are subject to landau damping . , width=336 ]    in fig . [ fig : fig5 ] we plot the rabi splitting @xmath150 given by eq . ( [ rabi ] ) as a function of various system parameters . panel ( a ) shows the dependence of rabi splitting on the effective dielectric constant of the structure , @xmath10 . in addition to aforementioned reduction of @xmath150 for larger @xmath10 , it has a maximum at @xmath151 followed by its sharp decrease for smaller @xmath10 . the latter is caused by the reduction of the _ total _ exciton energy @xmath152 as the exciton _ binding _ energy @xmath153 increases , which leads to a sharp drop of @xmath150 [ see eq . ( [ rabi ] ) ] . such a strong effect is caused by the cubic energy dependence of plasmon oscillator strength [ see eq . ( [ eq : pi_plpole ] ) ] . the rise of @xmath150 as a function of fermi energy , shown in panel ( b ) , is due to the resonance shift to smaller momenta ( @xmath154 ) which leads to increase of the coulomb overlap @xmath93 . it should be noted that behavior of @xmath155 is rather complicated and depends on qw width @xmath9 as well as separation between graphene and qw planes @xmath7 [ see eq . ( [ eq : ucv ] ) ] . it can be seen that for insb and htc qws with comparable @xmath156 the dependence is linear , while for inas ( larger @xmath156 ) and htc ( 0.2 ev ) ( smaller @xmath156 ) qws there are slight deviations from linearity . the overall magnitude is determined by the prefactor of @xmath157 originating from the graphene coulomb potential which suppresses plasmon - exciton coupling for small @xmath0 corresponding to large @xmath156 . the same prefactor determines the dependence of @xmath150 on spacer thickness @xmath7 in panel ( c ) . indeed , the exponential decay of plasmon - exciton coupling strength with spacer thickness causes @xmath150 to vanish for @xmath7 in the range from 7 nm to 24 nm for qws used , signaling the end of strong coupling regime . finally , panel ( d ) shows the dependence of @xmath150 on qw width @xmath9 . rabi splitting decreases as @xmath155 in eq . ( [ eq : ucv ] ) decreases for large @xmath9 . for smaller @xmath9 , however , the @xmath9-dependence of @xmath150 can become non - monotonic . the reason for this is that at small @xmath9 the exciton energy becomes higher due to the qw size - quantization , so the plasmon and exciton dispersion intersect at higher energies ( and momenta ) leading to a less efficient coupling . finally , for @xmath158 nm the plasmon dispersion enters the electron - hole excitation continuum before it meets the exciton dispersion . in this case , the graphene plasmon is short - lived and , therefore , the excition - plasmariton is no longer a stable excitation . in summary , we studied strong coupling between a graphene plasmon and a two - dimensional exciton in a narrow gap semiconductor quantum well separated from graphene by a potential barrier . we developed an analytical model based on hamiltonian similar to that for exciton - polaritons that accurately describes mixed exciton - plasmon states ( exciton - plasmaritons ) and calculated energy dispersion curves for several semiconductor materials . authors would like to thank g. khodaparast for a helpful discussion . work at lanl was performed under the nnsa of the u.s . doe at lanl under contract no . de ac52 - 06na25396 . work at jsu was supported by the nsf under grant no . the interaction hamiltonian between the qw and graphene is @xmath159 where @xmath160 is the electron field annihilation operator in qw . integral is evaluated over the volume of qw . the electrostatic potential due to charge density fluctuations in graphene is given by @xmath161 where the integral is evaluated over graphene . in the second quantization representation , electron field operators are expanded as @xmath162 where @xmath163 is the electron annihilation operator in band @xmath164 with ( quasi)momentum @xmath165 in the first brillouin zone and @xmath166 is the corresponding bloch function . substituting this definition of the field operator into eq . ( [ eq : hqwg ] ) , the hamiltonian takes the form @xmath167 where @xmath168 in turn , bloch functions can be expand into the basis of wannier functions , @xmath169 , as @xcite @xmath170 where @xmath171 is the total number of unit cells and the sum runs over lattice sites @xmath172 . wannier functions are normalized as @xmath173 . rewriting @xmath174 via wannier functions one obtains @xmath175 since the potential @xmath176 is smooth on atomic scale , we expand it within each unit cell as @xmath177 . the first term , after substituting into ( [ c2 ] ) yields @xmath178 where @xmath179 is the fourier transform of @xmath176 . to evaluate the second expansion term , we neglect the overlap of localized wannier wavefunctions at different lattice sites so that diagonal terms vanish . for non - diagonal terms we obtain @xmath180 where @xmath181 is the valence to conduction band transition dipole moment . using the definition of wannier functions , eq . ( [ eq : wannier ] ) , it can be rewritten in a more standard form @xmath182 where @xmath183 is the conduction ( valence ) bandedge bloch function , and the integration is assumed over the unit cell .    substituting @xmath184 and @xmath185 into eq . ( [ h2 ] ) we obtain @xmath186 with @xmath187 the latter expression can be rewritten in a more intuitive form @xmath188 where @xmath189 is the transition dipole moment operator of the qw unit cell . the operator of electric field due to charge density fluctuations in graphene is @xmath190 . the hamiltonian @xmath191 describes interaction of _ interband _ qw excitations ( excitons ) with the electric field of graphene . it represents the starting point of this work . in contrast , @xmath192 describes interactions of _ intraband _ density fluctuations in qw with those in graphene . after recasting in exciton basis , this term describes qw exciton scattering on graphene plasmons . at low exciton densities and at low temperatures these processes are weak and not considered here . the bare density correlation function ( polarization bubble ) , is calculated within the dirac electrons approximation as @xcite @xmath193}{\sqrt{\omega^{2}-v_{f}^{2}q^{2}}}\nonumber \\ + \left.\frac{\left[g(\delta_{-})+i\pi\right]\theta\left[{\rm re}\left\ { \delta_{-}\right\ } + 1\right]-g(\delta_{+})}{\sqrt{\omega^{2}-v_{f}^{2}q^{2}}}\right],\label{eq : full_pop}\end{gathered}\ ] ] where @xmath194 and @xmath195 . the square roots are chosen to yield positive real parts and the imaginary part of the logarithm is taken in @xmath196 $ ] range . fermi velocity and fermi level ( the latter determines the extent of graphene charge - doping ) are denoted by @xmath197 and @xmath198 , respectively . within the dirac electrons approximation , the density correlation function is insensitive to the sign of the fermi level , so in all the expressions here and in the main text @xmath198 has to be understood as @xmath199 . the two important limiting forms of the density correlation function are ( i ) the long wavelength limit ( @xmath200 , @xmath201 ) , and ( ii ) the static limit ( @xmath202 , @xmath203 ) . the long wavelength limit is given by @xmath204 the static limit of the bare density correlation function is obtained as @xmath205 the naive substitution @xmath206 to account for in - graphene scattering losses in eq . ( [ eq : full_pop ] ) ( @xmath98 is the electron scattering rate in energy units ) is inaccurate in a general case ( especially if @xmath98 is not small ) , since it does not preserve the particle conservation requirement . to correct for this , the more accurate mermin procedure is adopted , yielding @xcite @xmath207 the full ( or  dressed " ) density correlation function , which accounts for screening in graphene , is obtained within the random phase approximation as @xmath208 where @xmath209 is the two - dimensional fourier transform of the coulomb potential within the graphene s plane , @xmath210 . the effective dielectric constant of the environment @xmath10 is in general determined by the specific geometry of the graphene - based device and bulk dielectric constants of its material constituents . for example , in the simplest case of the interface between two homogeneous materials with dielectric constants @xmath211 and @xmath212 , the effective dielectric constant at the interface is given by @xmath213 . @xcite    the plasmon dispersion relation , @xmath214 , is found by requiring the real part of the denominator of eq . ( [ eq : pop_rpa ] ) to vanish . the taylor expansion of the denominator around this point ( up to leading terms in both real and imaginary parts ) leads to the possibility of approximating the full density correlation function within the so called plasmon pole approximation as ( for @xmath215 ) @xmath216 where @xmath217 is the amplitude of plasmon pole , and @xmath218 is the plasmon energy - dissipation rate . the coefficients of the taylor expansion of the denominator of eq . ( [ eq : pop_rpa ] ) are @xmath219 \right|_{\omega=\omega_q},\nonumber \\ b_q&=-v_q{\rm im}\left[\pi_{\gamma}(q,\omega_q)\right].\end{aligned}\ ] ] in the low-@xmath15 limit ( i.e. , @xmath220 ) , the density correlation function in the plasmon pole approximation can be obtained purely analytically by ( i ) substituting eqs . ( [ eq : pop_lowq ] ) and ( [ eq : pop_stat ] ) into eq . ( [ eq : mermin ] ) , and ( ii ) using the so obtained @xmath221 to evaluate the taylor expansion coefficients @xmath222 and @xmath223 at @xmath224 , i.e. , where the real part of the denominator of eq . ( [ eq : pop_rpa ] ) vanishes . the first step produces @xmath225 the second step produces eq . ( [ eq : plpole ] ) with @xmath226 thus , resulting in the expected @xmath227 plasmon dispersion . @xcite    it turns out that for the specific case considered here , i.e. , the plasmon pole approximation in the long wavelength limit , the same analytical expression for @xmath228 could have been obtained in the limit of small @xmath98 by using the substitution @xmath206 instead of the more general mermin s procedure . it has to be emphasized , however , that such an agreement is not general and hard to foresee . therefore , the more accurate mermin s procedure has to be favored over more approximate methods of introducing the finite scattering rate into the density correlation function.@xcite plasmon dispersion in graphene remains @xmath227 only at @xmath229 , and deviates from this simple relation at higher @xmath15 due to non - local effects originating from interaction of plasmons with intra - band electron - hole pair excitations . therefore , in high-@xmath37 environment , the dispersion of plasmons in graphene can deviate from @xmath230 at not very large @xmath15 . to correct for this , we have to find the position of a plasmon pole more accurately than that given by eq . ( [ eq : plpole_apprx ] ) . to this end , we expand the bare density correlation function of graphene at @xmath231 , which yields @xmath232.\ ] ] it is straightforward to check that at @xmath233 this expression reduces to eq . ( [ eq : pop_lowq ] ) . on the other hand , the limit @xmath202 yields eq . ( [ eq : pop_stat ] ) . the equation for the plasmon pole , @xmath234 is easily solved with respect to @xmath235 , yielding the corrected plasmon pole position as @xmath236 where @xmath237 . the amplitude of the plasmon pole becomes @xmath238 it turns out however that the correction to amplitude is rather small , so the original uncorrected amplitude , given in eq . ( [ eq : plpole_apprx ] ) can be used . 56ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) link:\doibase    10.1103/physrevb.87.075447 [ * * ,   ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) in  @noop _ _ ,   ( ,  ,  )  chap . @noop _ _ ,  ed . ( ,  ,  ) @noop * * ,   ( ) @noop * * ,   ( ) in  @noop _ _ ,  vol .  ,  ( ,  ,  ) p.   @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop _ _ ,  ed . ( ,  ,  ) @noop * * , ( ) @noop * * ,   ( )
Abstract:

Answer: we study strong coupling between plasmons in monolayer charge - doped graphene and excitons in a narrow gap semiconductor quantum well separated from graphene by a potential barrier . we show that the coulomb interaction between excitons and plasmons result in mixed states described by a hamiltonian similar to that for exciton - polaritons and derive the exciton - plasmon coupling constant that depends on system parameters . we calculate numerically the rabi splitting of exciton - plasmariton dispersion branches for several semiconductor materials and find that it can reach values of up to 50 - 100 mev .


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: this paper presents a formal statement and an assisted proof of a jordan curve theorem ( jct ) discrete version . in its common form , the theorem says that the complement of a continuous simple closed curve ( a jordan curve ) @xmath0 in an affine real plane is made of two connected components whose border is @xmath0 , one being bounded and the other not . the discrete form of jct we deal with states that in a finite subdivision of the plane , breaking a ring @xmath1 of faces increases by @xmath2 the connectivity of the subdivision . it is a weakened version of the original theorem where the question of bound is missing . however , it is widely used in computational geometry and discrete geometry for imaging , where connection is the essential information @xcite . in fact , we only are in a combinatoric framework , where any embedding is excluded , and where bounding does not make sense .    in computational topology , subdivisions are best described by map models , the most general being _ hypermaps _ we propose a purely combinatorial proof of jct based on this structure . the hypermap framework is entirely formalized and the proofs are developed interactively and verified by the coq proof assistant @xcite . using an original way to model , build and destruct hypermaps , the present work brings new simple constructive planarity and connectivity criteria . it proposes a new direct expression of jct and a simple constructive proof with algorithmic extensions . it is also a large benchmark for the software specification framework we have been developing in the last fifteen years for map models used in geometric modeling and computer imagery @xcite . the useful coq features are reminded and the whole process is described , but the full details of the proofs are omitted . section [ rw ] summarizes related work . section [ ma ] recalls some mathematical materials . section [ fh ] proposes basic hypermap specifications . section [ pc ] proves constructive criteria of hypermap planarity and connectivity . section [ fr ] inductively specifies the rings and their properties . section [ jct ] proves the discrete jct . section [ cl ] concludes . the jct is a result of classical plane topology , first stated by c. jordan in 1887 , but of which o. veblen gives the first correct proof in 1905 . in 1979 , tutte proposes operations and properties of combinatorial maps , _ e.g. _ planarity and euler s formula , defines rings and proves a discrete jct @xcite . our theorem statement is comparable , but our framework is modeled differently and all our proofs are formalized and computer - assisted .    in 2003 , g. bauer and t. nipkow specify planar graphs and triangulations in isabelle / isar to carry out interactive proofs of euler s formula and of the five colour theorem @xcite . however , they do not approach the jct . in 2005 , a. kornilowicz designs for the mizar project a semi - automated classical proof of a continuous form of jct in an euclidean space @xcite . in 2005 also , on his way towards the proof of the kepler conjecture in the flyspeck projet , t. hales proves the jct for planar rectangular grids with the hol light system , following the kuratowski characterization of planarity @xcite .    in 2005 always , g. gonthier _ _ prove the four colour theorem using coq . plane subdivisions are described by hypermaps , and euler s formula is used as a global planarity criterion @xcite . a local criterion , called _ hypermap jordan property _ , is proven equivalent . the main part of this work is the gigantic proof of the four colour theorem with hypermaps and sophisticated proof techniques . the hypermap formalization is very different from ours and it seems that jct is not explicitly proven there . finally , since 1999 , we carry out experiments with coq for combinatorial map models of space subdivisions @xcite . [ hd ] a _ hypermap _ is an algebraic structure @xmath3 , where @xmath4 is a finite set whose elements are called _ darts _ , and @xmath5 , @xmath6 are permutations on @xmath4 . if @xmath7 , @xmath8 is the @xmath9-_successor _ of @xmath10 , @xmath10 is the @xmath9-_predecessor _ of @xmath8 , and @xmath10 and @xmath8 are said to be @xmath9-_linked_.    in fig . [ fig : uhmap1 ] , as functions @xmath5 and @xmath6 on @xmath11 are permutations , @xmath12 is a hypermap . it is drawn on the plane by associating to each dart a curved arc oriented from a bullet to a small stroke : @xmath13-linked ( resp . @xmath2-linked ) darts share the same small stroke ( resp . bullet ) . by convention , in the drawings of hypermaps on surfaces , @xmath9-successors turn counterclockwise around strokes and bullets . let @xmath12 be a hypermap . ( orbits and hypermap cells)[oc ] + ( 1 ) let @xmath14 be @xmath15 functions in @xmath4 . the _ orbit _ of @xmath16 for @xmath14 is the subset of @xmath4 denoted by @xmath17 , the elements of which are accessible from @xmath10 by any composition of @xmath14 . + ( 2 ) in @xmath18 , @xmath19 is the @xmath13-orbit or _ edge _ of dart @xmath10 , @xmath20 its @xmath2-orbit or _ vertex _ , @xmath21 its _ face _ for @xmath22 , and @xmath23 its _ ( connected ) component_.    in fig . [ fig : uhmap1 ] the hypermap contains @xmath24 edges ( strokes ) , @xmath25 vertices ( bullets ) , @xmath25 faces and @xmath26 components . for instance , @xmath27 is the edge of dart @xmath26 , @xmath28 its vertex . faces are defined , through @xmath29 , for a dart traversal in counterclockwise order , when the hypermap is drawn on a surface . then , every face which encloses a bounded ( resp . unbounded ) region on its left is called _ internal _ ( resp . _ external _ ) . in fig . [ fig : uhmap1 ] , the ( internal ) face of @xmath30 is @xmath31 and the ( external ) face of @xmath32 is @xmath33 . let @xmath34 and @xmath35 be the numbers of darts , edges , vertices , faces and components of @xmath18 . ( euler characteristic , genus , planarity ) + ( 1 ) the _ euler characteristic _ of @xmath18 is @xmath36 . + ( 2 ) the _ genus _ of @xmath18 is @xmath37 . + ( 3 ) when @xmath38 , @xmath18 is said to be _ planar_.    for instance , in fig . [ fig : uhmap1 ] , @xmath39 and @xmath40 . consequently , the hypermap is non planar . these values satisfy the following results :    @xmath41 is an even integer and @xmath42 is a natural number .    a non empty connected @xmath43 i.e. with @xmath44 @xmath43 planar hypermap satisfies @xmath45 . when @xmath46 , the _ representation _ of @xmath18 on an _ orientable closed _ surface is a mapping of edges and vertices onto points , darts onto open oriented jordan arcs , and faces onto open connected regions . it is an _ embedding _ when every component of @xmath18 realizes a partition of the surface . then , the genus of @xmath18 is the minimum number of _ holes _ in an orientable closed surface where such an embedding is possible , thus drawing a subdivision , or a polyhedron , by hypermap component @xcite . for instance , all the components of the hypermap in fig . [ fig : uhmap1 ] can be embedded on a torus ( @xmath2 hole ) but not on a sphere or on a plane ( @xmath13 hole ) . when a ( planar ) hypermap component is embedded on a plane , the corresponding subdivision has exactly one unbounded ( external ) face . but a non planar hypermap can never be embedded on a plane : in a drawing on a plane , some of its faces are neither internal nor external , _ e.g. _ @xmath47 in fig . [ fig : uhmap1 ] . conversely , any subdivision of an _ orientable closed _ surface can be modeled by a hypermap . in fact , the formal presentation which follows is _ purely combinatorial _ , _ i.e _ without any topological or geometrical consideration .    -.3 cm      to state the version of jct we will prove , we need the concepts of _ double - link _ , _ adjacent faces _ and _ ring of faces _ in a hypermap @xmath48 . ( double - link and adjacent faces ) + ( 1 ) a _ double - link _ is a pair of darts @xmath49 where @xmath8 and @xmath50 belong to the same edge . + ( 2 ) the faces @xmath51 and @xmath52 of @xmath18 are said to be _ adjacent by the double - link _ @xmath49 if @xmath8 is a dart of @xmath51 and @xmath50 a dart of @xmath52 . we choose a face adjacency _ by an edge _ rather than _ by a vertex _ as does w.t . tutte @xcite . in fact , due to the homogeneity of dimensions @xmath13 and @xmath2 in a hypermap , both are equivalent . ( ring of faces ) [ rof ] + a _ ring of faces _ @xmath1 of length @xmath15 in @xmath18 is a non empty sequence of double - links @xmath53 , for @xmath54 , with the following properties , where @xmath55 and @xmath56 are the edge and face of @xmath57 : + ( 0 ) _ unicity : _ @xmath55 and @xmath58 are distinct , for @xmath59 and @xmath60 ; + ( 1 ) _ continuity : _ @xmath56 and @xmath61 are adjacent by the double - link @xmath53 , for @xmath62 ; + ( 2 ) _ circularity _ , or _ closure : _ @xmath63 and @xmath64 are adjacent by the double - link @xmath65 ; + ( 3 ) _ simplicity : _ @xmath56 and @xmath66 are distinct , for @xmath59 and @xmath60 . this notion simulates a jordan curve represented in dotted lines in fig . [ fig : uring4bis ] on the left for @xmath67 . then , we define the _ break along a ring _ , illustrated in fig . [ fig : uring4bis ] on the right . ( break along a ring ) [ bar ] + let @xmath1 be a ring of faces of length @xmath15 in @xmath18 . let @xmath68 , for @xmath69 , be a hypermap sequence , where the @xmath70 are recursively defined by : + ( 1 ) @xmath71 : @xmath72 = @xmath5 ; + ( 2 ) @xmath73 : for each dart @xmath74 of @xmath4 : @xmath75 = if @xmath76 then @xmath77 else if @xmath78 then @xmath57 else @xmath79 + then , @xmath80 is said to be obtained from @xmath18 by a _ break along _ @xmath1 . finally , the theorem we will prove in coq mimics the behaviour of a cut along a simple jordan curve of the plane ( or of the sphere ) into two components :    let @xmath18 be a planar hypermap with @xmath35 components , @xmath1 be a ring of faces in @xmath18 , and @xmath81 be the break of @xmath18 along @xmath1 . the number @xmath82 of components of @xmath81 is such that @xmath83 . -.3 cm      in coq , we first define an inductive type dim for the two dimensions at stake :    inductive dim : set:= zero : dim | one : dim .    all objects being typed in coq , dim has the type set of all concrete types . its _ constructors _ are the constants zero and one . in each inductive type , the generic equality predicate @xmath84 is built - in but its decidability is not , because coq s logic is intuitionistic . for dim , the latter can be established as the lemma :    lemma eq_dim_dec : forall i j : dim , i = j+  i = j . once it is made , its proof is an object of the sum type \{i = j}+\{`~`i = j } , _ i.e. _ a function , named eq@xmath85dim@xmath85dec , that tests whenever its two arguments are equal . the lemma is interactively proven with some tactics , the reasoning being merely a structural induction on both i and j , here a simple case analysis . indeed , from each inductive type definition , coq generates an _ induction principle _ , usable either to prove propositions or to build total functions on the type . we identify the type dart and its equality decidability eq@xmath85dart@xmath85dec with the built - in nat and eq@xmath85nat@xmath85dec . finally , to manage exceptions , a nil dart is a renaming of 0 :    definition dart:= nat . definition eq_dart_dec:= eq_nat_dec . definition nil:= 0 . -.3 cm      the hypermaps are now approached by a general notion of _ free map _ , thanks to a free algebra of terms of inductive type fmap with @xmath26 constructors , v , i and l , respectively for the _ empty _ ( or _ void _ ) map , the _ insertion _ of a dart , and the _ linking _ of two darts :    inductive fmap : set:= v : fmap | i : fmap->dart->fmap | l : fmap->dim->dart->dart->fmap . for instance , the hypermap in fig . [ fig : uhmap1 ] can be modeled by the free map represented in fig.[fig : hmap2 ] where the @xmath13- and @xmath2-links by l are represented by arcs of circle , and where the orbits remain open . again , coq generates an induction principle on free maps . observers _ of free maps can be defined . the predicate exd express that a dart exists in a hypermap . its definition is recursive , which is indicated by fixpoint , thanks to a pattern matching on m written match m with .... the attribute @xmath86struct m@xmath87 allows coq to verify that the recursive calls are performed on smaller fmap terms , thus ensuring termination . the result is false or true , basic constants of prop , the built - in type of propositions . note that terms are in prefix notation and that @xmath85 is a place holder :    fixpoint exd(m : fmap)(z : dart)struct m : prop:= match m with v = > false | i m0 x = > z = x exd m0 z | l m0 _ _ _ = > exd m0 z end . the decidability exd@xmath85dec of exd directly derives , thanks to a proof by induction on m. then , a version , denoted a , of operation @xmath88 of definition [ hd ] completed with nil for convenience is written as follows , the inverse a@xmath851 being similar :    fixpoint a(m : fmap)(k : dim)(z : dart)struct m : dart:= match m with v = > nil | i m0 x = > a m0 k z | l m0 k0 x y = > if eq_dim_dec k k0 then if eq_dart_dec z x then y else a m0 k z else a m0 k z end . predicates succ and pred express that a dart has a k - successor and a k - predecessor ( not nil ) , with the decidabilities succ@xmath85dec and pred@xmath85dec . in hypermap m of fig . [ fig : hmap2 ] , a m zero 4 = 3 , a m zero 5 = nil , succ m zero 4 = true , succ m zero 5 = false , a@xmath851 m one 2 = 1 . in fact , when a k - orbit remains open , which will be required in the following , we can obtain its top and bottom from one of its dart z. then , we can do as if the k - orbit were closed , thanks to the operations ca and ca@xmath851 which _ close _ a and a@xmath851 , in a way similar to operation @xmath89 of w.t . tutte @xcite . for instance , in fig . [ fig : hmap2 ] , top m one 1 = 3 , bottom m one 1 = 4 , ca m one 3 = 4 , ca@xmath851 m one 4 = 3 .    finally , _ destructors _ are also recursively defined . first , d : fmap->dart->fmap deletes the latest insertion of a dart by i. second , b , b@xmath85:fmap->dim->dart->fmap break the latest k - link inserted for a dart by l , forward and backward respectively . cm      preconditions written as predicates are introduced for i and l :    definition prec_i(m : fmap)(x : dart):= x < > nil /   exd m x. definition prec_l(m : fmap)(k : dim)(x y : dart):= exd m x /  exd m y /   succ m k x /   pred m k y /  ca m k x < > y.    if i and l are used under these conditions , the free map built necessarily has open orbits . in fact , thanks to the closures ca and ca@xmath851 , it can always be considered as a true hypermap exactly equipped with operations @xmath88 of definition [ hd ] . it satisfies the _ invariant _ :    fixpoint inv_hmap(m : fmap):prop:= match m with v = > true | i m0 x = > inv_hmap m0 /  prec_i m0 x | l m0 k0 x y = > inv_hmap m0 /  prec_l m0 k0 x y end . such a hypermap was already drawn in fig . [ fig : hmap2 ] . fundamental proven properties are that , for any m and k , ( a m k ) and ( a@xmath851 m k ) are _ injections _ inverse of each other , and ( ca m k ) and ( ca@xmath851 m k ) are _ permutations _ inverse of each other , and are closures . finally , traversals of faces are based on function f and its closure cf , which correspond to @xmath29 ( definition [ oc ] ) . so , in fig . [ fig : hmap2 ] , f m 1 = nil , cf m 1 = 5 . properties similar to the ones of a , ca are proven for f , cf and their inverses f@xmath851 , cf@xmath851 .    -.3 cm      testing if there exists a path from a dart to another in an orbit for a hypermap permutation is of prime importance , for instance to determine the number of orbits . the problem is exactly the same for @xmath5 , @xmath6 or @xmath29 ( definitions [ hd ] and [ oc ] ) . that is why a _ signature _ sigf with formal parameters f , f@xmath851 and their properties is first defined . next , a _ generic module _ ( or _ functor _ ) mf(m : sigf ) , the formal parameter m being a _ module _ of type sigf , is written in coq to _ package _ generic definitions and proven properties about f and f@xmath851 . among them , we have that each f - orbit of m is _ periodic _ with a positive smallest uniform period for any dart z of the orbit . the predicate expo m z t asserts the _ existence of a path _ in an f - orbit of m from a dart z to another t , which is proven to be a _ decidable equivalence_. note that most of the properties are obtained by _ noetherian induction _ on the length of iterated sequences of f - successors , bounded by the period . appropriate modules , called ma0 , ma1 and mf , are written to instantiate for ( ca m zero ) , ( ca m one ) and ( cf m ) definitions and properties of f. so , a generic definition or property in mf(m ) has to be prefixed by the module name to be concretely applied . for instance , mf.expo m z t is the existence of a path from z to t in a face . in the following , mf.expo is abbreviated into expf . for instance , in fig . [ fig : hmap2 ] , expf m 1 5 = true , expf m 5 3 = false . finally , a binary relation eqc stating that two darts belong to the same component is easily defined by induction . for instance , in fig . [ fig : hmap2 ] , we have eqc m 1 5 = true , eqc m 1 13 = false . we quickly prove that ( eqc m ) is a _ decidable equivalence_.    -.3 cm      we now count cells and components of a hypermap using the coq library module zarith containing all the features of z , the integer ring , including tools to solve linear systems in presburger s arithmetics . the numbers nd , ne , nv , nf and nc of darts , edges , vertices , faces and components are easily defined by induction . euler s characteristic ec and genus derive . the genus theorem and the euler formula ( for any number ( nc m ) of components ) are obtained as corollaries of the fact that ec is even and satisfies 2 * ( nc m ) > = ( ec m ) @xcite . remark that - > denotes a functional type in set as well as an implication in prop :    definition ec(m : fmap ) : z:= nv m + ne m + nf m - nd m. definition genus(m : fmap ) : z:= ( nc m ) - ( ec m)/2 . definition planar(m : fmap ) : prop:= genus m = 0 . theorem genus_theorem : forall m : fmap , inv_hmap m - > genus m > = 0 . theorem euler_formula : forall m : fmap , inv_hmap m - > planar m - > ec m / 2 = nc m. a consequence of the previous theorems is a completely constructive _ criterion of planarity _ , when one correctly _ links _ with l at dimensions @xmath13 or @xmath2 , _ e.g. _ for @xmath13 :    theorem planarity_crit_0 : forall ( m : fmap)(x y : dart ) , inv_hmap m - > prec_l m zero x y - > ( planar ( l m zero x y ) < - > ( planar m /  (   eqc m x y expf m ( ca_1 m one x ) y ) ) ) . so , the planarity of m is preserved for ( l m zero x y ) _ iff _ one of the following two conditions holds : ( 1 ) x and y are not in the same component of m ; ( 2 ) x@xmath851 = ( ca@xmath851 m one x ) and y are in the same face of m , _ i.e. _ the linking operates _ inside the face _ containing y. fig . [ fig : planarity ] illustrates @xmath13-linking inside a face , giving two new faces , and between two ( connected ) faces , giving a new face , thus destroying planarity . finally , after a long development , we prove the expected _ planarity criterion _ , when breaking a link with b , at any dimension , _ e.g. _ for @xmath13 :    lemma planarity_crit_b0 : forall ( m : fmap)(x : dart ) , inv_hmap m - > succ m zero x - > let m0 : = b m zero x in let y : = a m zero x in ( planar m < - > ( planar m0 /  (   eqc m0 x y expf m0 ( ca_1 m0 one x ) y ) ) ) . such a lemma is easy to write / understand as a _ mirror form _ of the @xmath13-linking criterion , but it is much more difficult to obtain . it would be fruitful to relate these constructive / destructive criteria with the static one of g. gonthier @xcite . finally , some useful results quickly characterize the effect of a link break on the _ connectivity _ of a planar hypermap . for instance , when @xmath13-breaking x , a disconnection occurs _ m y x0 :    lemma disconnect_planar_criterion_b0:forall ( m : fmap)(x : dart ) , inv_hmap m - > planar m - > succ m zero x - > let y : = a m zero x in let x0 : = bottom m zero x in ( expf m y x0 < - >  eqc ( b m zero x ) x y ) . -.3 cm      since an edge is always open in our specification , when doing the backward break of a unique @xmath13-link from y or y , we in fact realize a double - link break , as in definition [ bar ] . so , we choose to identify a double - link by the unique dart , we called x , where the @xmath13-link to be broken begins . in fact , with respect to the face f _ on the left of _ the double - link in the ring , there are two cases , depending on the position of x and its forward @xmath13-link , as shown in fig . [ fig : identface ] ( a ) and ( b ) . we decided to distinguish them by a boolean b. then , a double - link is coded by a pair ( x , b ) . so , we implicitely identify each ring face f by the double - link coding on its right in the ring . in fig . [ fig : identface ] ( a ) , face f is identified by ( x , true ) and contains y:= a m zero x , whereas in fig . [ fig : identface ] ( b ) , face f is identified by ( x , false ) and contains x0:= bottom m zero x. these modeling choices considerably simplify the problems . indeed , in closed orbits , a true double - link break would entail @xmath90 applications of b followed by @xmath90 applications of l , and would be much more complicated to deal with in proofs .    -.3 cm      first , we inductively define linear lists of pairs of booleans and darts , with the two classical constructors lam and cons , and usual observers and destructors , which we do not give , because their effect is directly comprehensible :    inductive list : set : = lam : list | cons : dart*bool - > list - > list . such a list is composed of couples ( x , b ) , each identifying a face f : if b is true , f is represented by y:= a m zero x , otherwise by x0:= bottom m zero x ( fig . [ fig : identface ] ) . in the following , bl m l breaks all the @xmath13-links starting from the darts of list l in a hypermap m. now , we have to model the conditions required for list l to be a ring of hypermap m. translating definition [ rof ] , we have four conditions , called pre@xmath85ringk m l , for k @xmath91 , which we explain in the following sections . finally , a predicate ring is defined by :    definition ring(m : fmap)(l : list):prop:=  emptyl l / pre_ring0 m l /  pre_ring1 m l /  pre_ring2 m l /  pre_ring3 m l.    -.3 cm      the predicate distinct@xmath85edge@xmath85list m x l0 saying that the edges of l0 are distinct in m from a given edge of x , pre@xmath85ring0 m l is defined recursively on l to impose that all edges in l are distinct : condition ( 0 ) of definition [ rof ] . it also imposes that each dart in l has a @xmath13-successor , in order to have well defined links , which is implicit in the mathematical definition , but not in our specification whith open orbits . fixpoint pre_ring0(m : fmap)(l : list)struct l : prop:= match l with lam = > true | cons ( x , _ ) l0 = > pre_ring0 m l0 /  distinct_edge_list m x l0 /  succ m zero x end . -.3 cm      then , we define adjacency between two faces identified by xb = ( x , b ) and xb = ( x , b ) , along the link corresponding to xb :    definition adjacent_faces(m : fmap)(xb xb:dart*bool):= match xb with ( x , b ) = > match xb with ( x,b ) = > let y : = a m zero x in let y:= a m zero x in let x0 : = bottom m zero x in let x0:= bottom m zero x in if eq_bool_dec b true then if eq_bool_dec b true then expf m x0 y else expf m x0 x0 else if eq_bool_dec b true then expf m y y else expf m y x0 end end . this definition is illustrated in fig . [ fig : adjacency ] for the four possible cases of double - link codings . so , the predicate pre@xmath85ring1 m l recursively specifies that two successive faces in l are adjacent : condition ( 1 ) in definition [ rof ] :    fixpoint pre_ring1(m : fmap)(l : list)struct l : prop:= match l with lam = > true | cons xb l0 = > pre_ring1 m l0 /   match l0 with lam = > true | cons xb l = > adjacent_faces m xb xb end end . -.3 cm      the predicate pre@xmath85ring2 m l specifies that the last and first faces in l are adjacent : condition ( 2 ) of circularity in definition [ rof ] :    definition pre_ring2(m : fmap)(l : list):prop:= match l with lam = > true | cons xb l0 = > match xb with ( x , b ) = > let y : = a m zero x in match l0 with lam = > let x0 : = bottom m zero x in expf m y x0 | cons _ l = > let xb:= last l0 in adjacent_faces m xb xb end end end .    -.3 cm      the predicate specifiying that the faces of m identified by xb and xb are distinct is easy to write by cases on the booleans in xb and xb. the predicate distinct@xmath85face@xmath85list m xb l0 expressing that the face identified by xb is distinct from all faces of list l0 entails . then , the predicate pre@xmath85ring3 m l says that all faces of l are distinct : condition ( 3 ) in definition [ rof ] :    fixpoint pre_ring3(m : fmap)(l : list)struct l : prop:= match l with lam = > true | cons xb l0 = > pre_ring3 m l0 /  distinct_face_list m xb l0 end . the general principle of the jct proof for a hypermap m and a ring l is a structural induction on l. the case where l is empty is immediatly excluded because l is not a ring by definition . thus the true first case is when l is _ reduced to one element _ , _ i.e. _ is of the form cons ( x , b ) lam . then , we prove the following lemma as a direct consequence of the planarity criterion planarity@xmath85crit@xmath85b0 and the criterion face@xmath85cut@xmath85join@xmath85criterion@xmath85b0 :    lemma jordan1:forall(m : fmap)(x : dart)(b : bool ) , inv_hmap m - > planar m - > let l:= cons ( x , b ) lam in ring m l - > nc ( bl m l ) = nc m + 1 .    when a ring l1 contains _ at least two elements _ , we prove that the condition ` ~`expf m y x0 must hold with the first element ( x , b ) of l1 ( in fact , conditions @xmath92 and @xmath93 are enough ) :    lemma ring1_ring3_connect : forall(m : fmap)(x x:dart)(b b:bool)(l : list ) , inv_hmap m - > let l1:= cons ( x , b ) ( cons ( x,b ) l ) in let y:=a m zero x in let x0:= bottom m zero x in planar m - > pre_ring1 m l1 - > pre_ring3 m l1 - >  expf m y x0 .    in this case , thanks to disconnect@xmath85planar@xmath85criterion@xmath85b0 ( section [ pc ] ) , the lemma entails that the break of the first ring link does never disconnect the hypermap . then , after examining the behavior of pre@xmath85ringk , for k @xmath91 , we are able to prove the following lemma which states that the four ring properties are preserved after the first break in l :    lemma pre_ring_b : forall(m : fmap)(l : list ) , inv_hmap m - > planar m - > let x : = fst ( first l ) in let y : = a m zero x in let x0 : = bottom m zero x in let m1 : = b m zero x in  expf m y x0 - > ring m l - > ( pre_ring0 m1 ( tail l ) /  pre_ring1 m1 ( tail l ) /  pre_ring2 m1 ( tail l ) /  pre_ring3 m1 ( tail l ) ) . the most difficult is to prove the part of the result concerning pre@xmath85ringk , for k @xmath91 . the four proofs are led by induction on l in separate lemmas . for pre@xmath85ring0 , the proof is rather simple . but , for the other three , the core is a long reasoning where @xmath90 , @xmath26 or @xmath94 links are involved in input . since each link contains a boolean , sometimes appearing also in output , until @xmath95 cases are to be considered to combine the boolean values . finally , from jordan1 and pre@xmath85ring@xmath85b above , we have the expected result by a quick reasoning by induction on l , where links are broken one by one from the first :    theorem jordan : forall(l : list)(m : fmap ) , inv_hmap m - > planar m - > ring m l - > nc ( bl m l ) = nc m + 1 .    it is clear that , provided a mathematical hypermap @xmath18 and a mathematical ring @xmath1 conform to definitions [ hd ] and [ rof ] , we can always describe them as terms of our specification framework in order to apply our jct . conversely , given a hypermap term , some mathematical rings can not directly be written as terms . to do it , our ring description and our jct proof have to be slightly extended . however , that is not necessary for the _ combinatorial maps _ ( where @xmath5 is an involution ) terms , for which our ring specification and our jct formalization are _ complete_. this is more than enough to affirm the value of our results . we have presented a new discrete statement of the jct based on hypermaps and rings , and a formalized proof assisted by the coq system . our hypermap modeling with _ open orbits _ simplifies and precises most of known facts . it also allows to obtain some new results , particularly about hypermap construction / destruction , connection / disconnection and planarity . this work involves a substantial framework of hypermap specification , which is built _ from scratch _ , _ i.e. _ exempt from any proper axiom . it is basically the same as the one we have designed to develop geometric modelers via algebraic specifications @xcite . so , we know how to efficiently implement all the notions we formally deal with . the coq system turned out to be a precious auxiliary to guide and check all the process of specification and proof . the preexistent framework of hypermap specification represents about 15,000 lines of coq , and the jct development about 5,000 lines , including about 25 new definitions , and 400 lemmas and theorems . note that all results about the dimension 0 were actually proven , but some planarity properties about dimension 1 , which are perfectly symmetrical , have just been admitted . however , the jct formal proof is complete . so , we have a solid foundation to tackle any topological problem involving orientable surface subdivisions . extensions are in 2d or 3d computational geometry and geometric modeling by introducing embeddings @xcite , and computer imagery by dealing with pixels @xcite or voxels .
Abstract:

Answer: this paper presents a formalized proof of a discrete form of the jordan curve theorem . it is based on a hypermap model of planar subdivisions , formal specifications and proofs assisted by the coq system . fundamental properties are proven by structural or noetherian induction : genus theorem , euler s formula , constructive planarity criteria . a notion of ring of faces is inductively defined and a jordan curve theorem is stated and proven for any planar hypermap .    jean - franois dufourd