Instruction: Write an abstract for the following acadamic paper.

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

Paper: the isolation of single - layer and few - layer graphene has triggered a huge burst of interest , mainly motivated by the observation of unconventional electronic properties , which stem from the dirac - like low - energy electronic structure of graphene characterized by a gapless conical dispersion [ ] . the huge electronic mobility of free - standing graphene originates from its chiral properties , tightly linked with the lack of electron backscattering phenomena near the fermi level [ ] . a drawback of this characteristic band structure is the absence of an energy gap between the dirac cones , which would be highly desirable for the exploitation of graphene in device applications . so far , however , the effective employment of graphene - based materials in low - energy electronics has been hindered by the difficulty of opening a bandgap without affecting the electronic mobility . understanding the fundamental mechanisms responsible of gap opening in graphene is thus of the highest importance in the perspective of engineering new efficient switch on / off devices .    from the theoretical point of view , the simplest way to open a gap in the conical dirac - like dispersion of a two - dimensional honeycomb material is to induce an inequivalence between the two carbon sublattices a and b. this corresponds to include a @xmath0 term in the dirac - like hamiltonian : @xmath1 + \frac{\delta}{2}\hat{\sigma}_z , \label{hgap}\end{aligned}\ ] ] where @xmath2 is the dirac velocity , @xmath3 is the momentum relative to the k point , and where @xmath4 are @xmath5 pauli matrices defined in the subspace of the two carbon orbitals for unit cell . the term @xmath6 represents an energy off - shift of the dirac point due to finite doping . the total dispersion results thus : @xmath7 , where @xmath8 . note that this model predicts for small @xmath3 a parabolic behavior @xmath9 for both conduction and valence bands with an effective mass @xmath10 proportional to the band gap @xmath11 : @xmath12 . in order to highlight the strict relation between the gap opening and the onset of an effective mass , we define this scenario as `` massive '' ( ms ) gap model , in contrast with a `` massless '' ( ml ) gap model that we will discuss below . it should be worth to note , in addition , that , in the massive gap model , the band dispersion recovers a standard linear behavior ( @xmath13 ) for @xmath14 , so that the extrapolation of the upper band dispersion overlaps the lower band and vice versa . angle - resolved photoemission spectroscopy ( arpes ) is one of the most direct experimental techniques that allow to investigate the energy - momentum dependence of the electronic states in solids . although arpes can be applied also to suspended graphene , as recently demonstrated in ref . [ ] , it is more often applied to study the electronic structure of supported graphene . when graphene is in contact with a solid material finite doping effects , as well as modifications ( screening ) of the many - body interactions , are observed . gap opening is frequently reported , with magnitudes that depend on the growth process and on the substrate [ ] . these observations are commonly interpreted in term of the massive model discussed above . at a closer look , however , the actual evidence of such physical behavior in most of the cases is not assessed . a first issue regards the effective opening of a gap in some epitaxially grown graphene samples , where a spectral anomaly at the dirac point was shown to be associated with the signature of plasmaronic subbands rather than with a simple gap opening [ ] . on the other hand , arpes data for other supported graphene layers do not show the characteristic `` diamond - shaped '' dispersion of the plasmaron at the dirac point [ ] . even in these cases , however , the actual experimental band dispersion presents many inconsistencies with the simple massive dirac gap model , pointing thus towards different mechanism of gap opening .    the arpes data of ref . [ ] , for instance , reported a gap @xmath15 ev which was initially discussed in terms of the massive gap model , and critically revised in ref . two main unconventional features were there pointed out in regards to the arpes data : ( @xmath16 ) in spite of a large gap opening , the conduction and valence bands in ref . [ ] retained a conical shape , in contrast with the expected parabolic behavior ; ( @xmath17 ) the linear extrapolation of the upper and lower bands appeared to be _ misaligned _ , with an energy shift corresponding to the gap .    in order to account for these unconventional features , an alternative scenario was proposed in ref . [ ] in terms of a phenomenological `` massless '' gap model with band dispersion @xmath18 $ ] . the anomalous features ( @xmath16)-(@xmath17 ) could be naturally reproduced assuming a self - energy of the form @xmath19 $ ] , where @xmath20 . this phenomenological model was also shown to provide precise predictions ( e.g. on the density of states ) which could be experimentally checked [ ] . the aim of the present paper is twofold . on one hand we would like to provide experimental evidence of the validity of the massless gap model in cvd grown graphene . on the other hand we discuss in details at the theoretical level the possibility of the appearing of a massless gap as result of a tendency towards a second - order phase transition . we believe that a full understanding and controlling of the unconventional properties of the massless gap in graphene on substrates can open new perspectives in the bandgap engineering in graphene - based materials . before addressing a quantitative analysis of the experimental arpes dispersion , and a deeper discussion about the possible origin of the anomalous bandgap features there observed , we summarize briefly here the comparison between the massive and massless gap models . we first write the non - interacting hamiltonian in the form : @xmath21 where @xmath3 are the momenta relative to the k point . the eigenvalues of eq . ( [ h0 ] ) , @xmath22 , describe the dirac cone dispersion , while the non - trivial dependence of ( [ h0 ] ) on the angle @xmath23 accounts for the chiral properties of the eigenstates . a massive gap is induced if the a and b sublattices are electrostatically inequivalent , for instance as an effect of the substrate . such inequivalence is formally taken into account by the additional term @xmath24 as in eq . ( [ hgap ] ) , resulting in the well - known dispersion @xmath25 two regimes can be identified in this energy dispersion . at high - energies ( @xmath26 ) the effect of the gap is negligible and the band dispersion recovers the normal linear behavior @xmath13 . at low - energies ( @xmath27 ) the band acquires a parabolic shape @xmath28 $ ] , with @xmath12 . it is worth noting that in this regime the chiral structure of the band dispersion is strongly affected by the opening of the gap [ ] . the onset of a massive gap in the low - energy electronic structure of graphene has been discussed not only as due to the interaction with the substrate , but also as a possible effect of a spontaneous excitonic instability induced by the long - range coulomb interaction [ ] . in this latter context , the term ( [ gap ] ) should be regarded as an _ order parameter _ , rather than an intrinsic property of the system . for real graphene systems , however , the dimensionless coulomb coupling constant results to be much smaller than the critical value required for the excitonic instability . on the other hand , in the absence of such excitonic broken symmetry , in the normal state the long - range coulomb interaction has been shown to be accounted , at the leading order , by a self - energy term [ ] : @xmath29    note here the much weaker dependence on @xmath30 , than the bare dirac hamiltonian . indeed , the conventional linear behavior in @xmath30 is here counterbalanced by a logarithmic divergence @xmath31 induced by the long - range coulomb interaction . note finally that in this case the chiral structure is preserved down to @xmath32 .    the massless gap model investigated in this paper , alternative to the massive model described by eq . ( [ gap ] ) , resembles features from both ( [ gap ] ) and ( [ coulomb ] ) , and it can be mathematically described by a self - energy term [ ] : @xmath33 where the off - diagonal terms do _ not _ scale to zero for @xmath32 , but they still preserve the full chiral properties . as mentioned above , the characteristic energy - momentum dispersion results to be @xmath34 . \label{eless}\end{aligned}\ ] ] note that this models accounts in a simple way for the anomalous features described in the introduction , namely : @xmath16 ) the dirac cones appear to be just split by the gap @xmath11 without affecting their conical shape ; @xmath17 ) the perfectly linear behavior of the upper ( lower ) band does not extrapolate onto the linear behavior of the lower ( upper ) band .    as we are going to see , such anomalies can permit to identify in a clear way one gap model rather than the other one in the experimental data . in this section we present a detailed analysis of some recent arpes measurements on single - layer graphene on ir(111 ) decorated by ir clusters and na alkali metals . we show how these measurements display , in correspondence with the opening of a large energy gap [ ] , an anomalous behavior of the dirac cones compatible with the massless gap model . the experiment was performed at the vuv - photoemission beamline on the elettra storage ring in trieste under ultrahigh vacuum conditions . the ir(111 ) crystal was cleaned by cycles of ar@xmath35 ion sputtering and annealing at @xmath36 k. surface order and cleanliness of the sample were checked by low - energy electron diffraction and core level photoemission measurements . graphene was grown by thermal cracking of ethylene ( c@xmath37h@xmath38 ) on the ir substrate held at 1300 k for an exposure of 100 l ( 1l corresponds to an exposure of @xmath39 mbar for 1s ) . under these experimental conditions , graphene had lattice vectors aligned in plane to those of the substrate and displayed a moir pattern originating from the interface lattice mismatch [ ] . tiny ir amounts were evaporated from a current heated thin film plate ( @xmath40 mm width , @xmath41 mm thickness ) at an evaporation rate of about @xmath42 ml / s , as determined by core level photoemission measurements [ ] . ir evaporation was performed at a temperature of 350 k. the deposition of ir results in the nucleation of size selected ir clusters on the hcp regions of the moir supercells ( see fig . [ f - sketch ] for a schematic sketch ) [ ] . the process was finely controlled to saturate the hcp regions of the moir - derived graphene supercells with one ir cluster and to avoid cluster percolation . the optimal ir coverage was determined to be 0.15 ml ( monolayer ) ir , with reference to the density of an ir(111 ) plane [ ] . additional na was evaporated from a commercial getter source at room temperature and adsorbed to fill the residual exposed graphene surface . the resulting system ( na / ir / g ) was characterized by means of angle - resolved photoemission , as reported in ref . [ ] . the valence band photoemission analysis was carried out by a scienta r4000 electron energy analyzer with @xmath43 ev photons , close to the cooper minimum of the ir5d levels . the spectra were collected at 120 k with energy and angular resolution of 30 mev and 0.3@xmath44 , respectively .    the arpes intensity map , along the direction perpendicular to @xmath45k passing through the k point ( p@xmath45k ) is shown in fig . [ f - map]a , where the white and green symbols represent the position of the maxima in the energy - distribution - curves ( edcs ) and in the momentum - distribution - curves ( mdcs ) , respectively . the fermi momentum along the p@xmath45k direction is @xmath46 @xmath47 that , taking into account the trigonal warping of the bands , gives an electron doping of @xmath48 electrons per c atom . from the edc at @xmath49 we estimate a bandgap @xmath50 ev . two remarkable features stand out from the arpes data : ( @xmath16 ) both conduction and valence bands present a linear dispersion throughout the probed @xmath3-space region , except within a very small region @xmath51 @xmath47 . this range is thus much smaller than that predicted by the massive gap model ( @xmath52 @xmath47 ) ; ( @xmath17 ) the linear extrapolations of the upper and lower branched do not match , at variance with the expectations of the massive gap model . we note additionally the lack of any diamond - like features close to the k point , which would be typical trademarks of plasmaronic coupling [ ] . on the other hand , plasmaron signatures in arpes data have been so far reported only for epitaxial graphene on sic , whereas they are absent in graphene grown on metallic substrates , most probably as a consequence of screening effects . such band structure of the na / ir / g system appears to be incompatible with the standard massive gap model , whereas they could be naturally reproduced within the context of the massless gap model .    in order to address at a more quantitative level this issue , we provide here a detailed careful analysis of the arpes data for na / ir / g . before testing the different models , it is necessary to estimate the basic parameters which are _ independent _ of the model itself , e.g. the fermi velocity @xmath2 and the the dirac energy @xmath6 in the absence of the gap . we estimate @xmath53 from the linear behavior of the lower band in the energy window from @xmath54 ev to @xmath55 ev , which yields @xmath56 ev /  . this value is consistent with the speed value derived from the upper band . we found this value was also perfectly compatible with the linear slope of the upper band . a rough estimate @xmath57 ev of the dirac energy in the absence of the gap is obtained as the midgap point at @xmath49 ( @xmath58 ev , @xmath59 ev ) . this estimate can be cross - checked with the average energy between the upper and lower band @xmath60 for generic @xmath3 , displayed in fig . [ f - map]b . the result shows a negligible momentum dependence , resulting in @xmath61 ev . note that the negligible dependence of @xmath62 on @xmath3 signalizes an almost perfect symmetry of the experimental upper and lower bands with respect to the dirac energy . this is compatible with both the massless and massive gap models , but it rules out plasmaronic effects which would result in an asymmetry and in a finite @xmath3 dependence of @xmath62 . the unbiased determination of @xmath6 and @xmath2 allows to test the validity of the different models on the experimental data . [ f - disp ] compares the band dispersion extracted from the momentum and energy distribution curves with the prediction of massless and massive modes . the best agreement is found for the massless model with @xmath63 ev . note the mismatch between the upper and lower dirac cones , that reproduce perfectly the experimental data . the same value of the gap can be indeed also inferred from such misalignment . on the other hand , no satisfactory fit was possible for the massive model . we show here thus the predictions of the massive model for @xmath64 ev , that reproduces the value at k of the upper and lower bands . in the previous section we have shown how a careful analysis of the band anomalies observed by arpes in in na / ir / g points out towards a phenomenology of the graphene bands compatible more with massless gap model rather than with the standard massive one . it is clear that the different models represent a trademark for different underlying physical processes . in this section we address thus this issue , in order to clarify which kind of physical processes can be responsible at the microscopic level for the massless dirac phenomenology , with the final aim to control and engineering such processes .    along this perspective , a brief summary of the microscopic physics related to the massive model can be useful . as mentioned in the introduction , the massive dirac model essentially stems from the presence of a local potential term @xmath0 which makes the two carbon atom sublattices inequivalent . this potential can be naturally induced by the interaction of graphene with the substrate . the term @xmath65 in the hamiltonian appears then as an _ external _ field arising from the environment . alternatively , the onset of a mass term @xmath66 has been discussed within the context of a spontaneous symmetry breaking , where the driving mechanism is the unscreened long - range coulomb interaction [ ] . the term @xmath67 can be viewed thus in this scenario as an anomalous _ self - energy _ term , i.e. as an _ order parameter _ , arising from many - body effects . it should be however mentioned that in this case the mass term @xmath68 appears to be momentum dependent , i.e. divergent for @xmath69 [ ] , rather than @xmath3-independent .    from a general point of view , since the presence of a linear dirac cone in ideal flat graphene is dictated by symmetry reasons , it is clear that the opening of a gap and the splitting of the upper and lower bands , both in the massive and in the massless model , must be associated with some kind of breaking symmetry . we have indeed just seen how the symmetry breaking ( spontaneous or not ) responsible for the massive model is the equivalence between the local potentials of the two sublattices . in ref . [ ] the general constraints that a phenomenological self - energy must satisfied to give rise to a massless dirac model have been discussed . it is on the other hand intriguing to discuss in the present work the conditions that could make possible the appearance of a self - energy ( [ less ] ) as result of a tendency towards a spontaneous symmetry breaking . it is clear that the effective breaking of symmetry , in real materials , can be prevented by many causes . for instance , as we will discuss below , a good candidate for a quasi - breaking of symmetry is the long - range coulomb interaction . in real materials , and in particular in na / ir / graphene , such interaction is expected to be screened resulting in a cut - off of the singular behavior for small exchange momenta @xmath70 . asymptotically , this should correspond to `` regularization '' of the interaction , and thus to a linear behavior of the dispersion ( although with a very steep slope ) for momenta @xmath3 close to the k point . for larger @xmath3 , however , the tendency towards an instability can still be reflected in a quasi - gapped dispersion which can be compatible with the experimental observations . from this point of view , we present the analysis of the onset of a spontaneous symmetry breaking just as an illustrative example for the gross features of the resulting phenomenology . in similar way , we do not want to provide here a quantitative evaluation of the possible critical couplings required to induce a massless gapped case , rather to compare , at a qualitative level , this strength of this instability with the strength of the corresponding instability associated with the spontaneous generation of a mass according the model ( [ hgap ] ) . in order to have thus a direct and clear way to compare the two possible instabilities , we address here this issue within the context of a mean - field solution , keeping in mind that a more compelling high - order analysis is needed for a quantitative investigation .    in order to address this issue in the most convenient way , we introduce the nambu notation where the hamiltonian of ideal free - standing graphene can be written as : @xmath71 where @xmath72 is a global index taking into account both spin and valley ( degenerate ) degrees of freedom @xmath73 , @xmath74 ) , @xmath75 is the spinor in the sublattice basis @xmath76 , and where @xmath77.\end{aligned}\ ] ]    we assume a long - range coulomb interaction which can we written in the momentum space as : @xmath78 where @xmath79 , and where @xmath80/\epsilon_0 \kappa   |{\bf q}|$ ] . the term @xmath81 $ ] in @xmath82 takes into account the subtraction of the positive charged background , while @xmath83 is the dimensionless relative in - plane dielectric constant . the phenomenology of the normal state self - energy arising from such many - body interaction , as well as of several possible _ anomalous _ self - energies associated with different spontaneous symmetry breaking , has been discussed in an extended way using different accurate techniques , as for instance renormalization group ( rg ) or static random - phase approximation ( rpa ) . in general , in order to have an accurate estimate of the _ magnitude _ of the corresponding self - energies ( and of the minimum strength of the interaction in the case of spontaneous symmetry breaking ) , a self - consistent screening of the quasi - particle excitations and of the the effective screened coulomb interaction is required .      neglecting the hartree term , which does not play any role in the present context , and focusing on the exchange ( fock ) term , we can write the generic self - energy as : @xmath84 note that , since the interaction @xmath82 is peaked at small @xmath85 , the intra - valley scattering is dominant , supporting the validity of the low - energy model versus a full tight - binding treatment . note also that at this level of approximation ( unscreened static coulomb interaction ) the self - energy does not depend on the frequency @xmath86 . within this scheme is thus convenient to define an effective mean - field hamiltonian @xmath87 , and to write a recursive equation : @xmath88    since @xmath89 depends on the angular variables @xmath90 , @xmath91 essentially through only @xmath92 , one can see that @xmath93 . apart numerical prefactor , the momentum dependence of @xmath94 on @xmath30 can be caught already by a perturbative analysis . replacing @xmath95 with @xmath96 in the right side term of eq . ( [ heff ] ) , or , equivalently , neglecting @xmath97 in the right side term of eq . ( [ fock ] ) , we obtain : @xmath98 \label{alpha}\end{aligned}\ ] ] where @xmath99 is the coupling constant for suspended graphene , @xmath100 is a high - momentum cut - off limiting the validity of the dirac model , approximately determined by the size of the brillouin zone , and where @xmath101    a careful analysis of ( [ emish ] ) shows that at the leading order @xmath102const . , a well known result . note that , within this context , the self - energy associated with the massless gap model is formally described by @xmath103const . i.e. by @xmath104 . it is clear that such analytical dependence is not spontaneously generated by the self - consistent solution of ( [ heff ] ) in the normal state , and it must be regarded as an order parameter of a possible second order phase transition . in order to address this possibility , we rewrite eq . ( [ heff ] ) as : @xmath105 & = & \hat{h}_0({\bf k } ) , \label{hself}\end{aligned}\ ] ] where @xmath106 $ ] is the functional @xmath106=t \sum_{{\bf p},n } v({\bf k - p})/ [ i\omega_n\hat{i}-\hat{x}]$ ] . broken symmetry phases can be generated when the equation @xmath107 & = & 0 \label{inst}\end{aligned}\ ] ] admits a non - trivial solution @xmath108 . we are here interested only in investigating the instability of the normal state towards a broken symmetry phase . within this context we can expand thus the right side term of eq . ( [ inst ] ) at the _ linear _ order with respect to @xmath109 . at @xmath110 we get thus the susceptibility equations : @xmath111 where @xmath112 $ ] is the energy dispersion in the normal state . within the spirit of a mean - field approach we can neglect in @xmath113 the self - energy many - body effects driven by the coulomb interactions , and approximate @xmath114=\hbar v |{\bf p}|$ ] . ( [ lastgap ] ) can be used to investigate the instability towards both a massive as well as massless gap . in the first case @xmath115 , which clearly breaks the symmetry of the original hamiltonian , while in the second case @xmath116 . since terms @xmath117 are already present in the normal state , it is not straightforwardly apparent at this level the nature of the symmetry breaking . on this regards it should be noticed that in the massless gap model the symmetry breaking is more specifically associated with the anomalous @xmath3-dependence of the order parameter , in similarity with @xmath118 flux phases discussed in different contexts , rather than with the pauli structure .    to make clearer this point , we employ the massless gap model @xmath119 $ ] as an ansatz for the _ right _ side of eq . ( [ lastgap ] ) . it is easy to check that the equations for @xmath120 , @xmath121 result to be decoupled and degenerate . from a careful analysis of the angular variable , it is also straightforward to check that the resulting @xmath109 in the left side of eq . ( [ lastgap ] ) appears to be indeed @xmath122 $ ] , where the prefactor @xmath123 will be determined by the effective self - consistency of eq . ( [ lastgap ] ) . using the explicit expression for the coulomb interaction , we obtain thus : @xmath124    it is remarkably relevant to notice that , contrary to eq . ( [ emish ] ) , the integral over @xmath125 is well behaved for both @xmath32 and @xmath126 . in this limits we can thus properly look for the self - consistent solution @xmath127 . we obtain thus the susceptibility equation : @xmath128 where @xmath129 \simeq 1.13 . \label{chi2}\end{aligned}\ ] ] for free - standing graphene with @xmath130 we can predict thus a finite critical coupling @xmath131 above which the system undergoes spontaneously a transition towards a massless gap phase . we should stress that such value of @xmath132 has been derived within the lowest order approximation , and a more compelling analysis is needed for a reliable quantitative estimate . it is however instructive to compare such value with the critical coupling required to induce , _ at the same level of approximation _ , a spontaneous generation of mass . such instability can be determined by using the ansatz @xmath133\hat{\sigma}_z$ ] in eq . ( [ lastgap ] ) . after few straightforward steps , we obtain the susceptibility equation @xmath134 where @xmath135 corresponding , for @xmath130 , to the well - known result @xmath136 [ ] .    although the critical coupling @xmath137 for the spontaneous generation of massless gap results to be larger than the critical coupling for the massive gap , it remains still of the same order of magnitude of @xmath132 within the same level of approximation . we know however that @xmath132 is significantly increased when higher order renormalization effects ( dynamical screening , self - energy effects ,  ) are taken into account , with possible @xmath132 as large as @xmath138 , of the same order thus of the present estimate of @xmath137 . becomes thus crucial , in order to assess the dominant instability of the normal state , to include high - order effect also for the massless gap case , to permit a quantitative comparison . in addition , it should not be excluded a possible interplay , in the broken symmetry phase , between the two order parameters , where one could favor and make more stable the other one . this scenario is suggested by the simple analysis of the experimental data of fig . [ f - disp ] , where the massless gap model reproduce well the experimental dispersion in a broad energy range , but where a finite parabolic behavior is nevertheless present at very small @xmath30 . we have already discuss how such weak parabolic behavior can not account in any case for the large observed gap . the experimental dispersion could be thus compatible with the possible coexistence of the both qualitatively different gaps , with the massless ( larger ) one determining the high - energy range , and the massive ( smaller ) one accounting for the residual parabolic behavior at small @xmath30 . further analysis in this direction is however needed to investigate this issue .    in conclusion , in this work we have provided , using angle - resolved photoemission spectroscopy , a compelling evidence of the presence of an unconventional gap in cvd grown graphene , properly described by a massless gap model . we have further investigated , on a theoretical ground , the microscopic conditions that can give rise to a spontaneous phase transition towards such massless gapped phase . we have shown that , for suspended graphene , a critical value of the coulomb coupling is needed , with larger but comparable values with the ones required to induce a breaking symmetry towards a massive map . this analysis , stricty valid for suspended graphene , can suggest a possible path to induce such unconventional gap also in cvd graphene , although in this case specific features of these materials ( screening of the coulomb interaction , interference with the moir pattern ) must be taken into account . the fully understanding and controlling of the unconventional properties of the massless gap in graphene can open new perspectives in the bandgap engineering in graphene - based materials . e.c . acknowledges support from the european project fp7-people-2013-cig `` lsie_2d '' and italian national miur prin project 20105zztse . l.b . acknowledges support from the italian miur under the projects firb - hybridnanodev - rbfr1236vv and prin-2012x3yfz2 . financial support was also provided by italian miur project firb - futuro in ricerca 2010-project plasmograph grant no . rbfr10m5bt . et al . _ , b * 84 * , 115401 ( 2011 ) . s.  y.  zhou , g.  -h . gweon , a.  v.  fedorov , p.  n.  first , w.  a.  de  heer , d.  -h . lee , f.  guinea , a.  h.  castro  neto , and a.  lanzara , nat . mat . * 6 * , 770 ( 2007 ) .
Abstract:

Answer: the opening of a gap in single - layer graphene is often ascribed to the breaking of the equivalence between the two carbon sublattices . we show by angle - resolved photoemission spectroscopy that ir- and na - modified graphene grown on the ir(111 ) surface presents a very large unconventional gap that can be described in terms of a phenomenological `` massless '' dirac model . we discuss the consequences and differences of this model in comparison of the standard massive gap model , and we investigate the conditions under which such anomalous gap can arise from a spontaneous symmetry breaking . keywords : graphene , bandgap , dirac cone , angle - resolved - photoemission


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: jordan  brans  dicke ( jbd ) theory @xcite is the simplest extended theory of gravity , depending on one additional parameter , the brans  dicke coupling @xmath5 , as compared to general relativity . as einstein s theory is recovered in the limit @xmath6 , there will always be viable jbd theories as long as general relativity remains so too . as such , it acts as a laboratory for quantifying how accurately the predictions of general relativity stand up against observational tests . the most stringent limits are derived from radar timing experiments within our solar system , with measurements using the cassini probe @xcite now giving a two - sigma lower limit @xmath7 ( improving pre - existing limits @xcite by an order of magnitude ) .    with precision cosmological data now available , particularly on cosmic microwave background ( cmb ) anisotropies from the wilkinson microwave anisotropy probe ( wmap ) @xcite , it has become feasible to obtain complementary constraints from the effect of modified gravity on the structure formation process , as suggested in ref .  @xcite . that paper focussed on the way that @xmath5 alters the hubble scale at matter  radiation equality , which is a scale imprinted on the matter power spectrum , in an attempt to identify how large an effect can be expected . subsequently , the expected total intensity and polarization microwave anisotropy spectra in the jbd theory were computed , and a forecast of the sensitivity to @xmath5 of data from the wmap and planck satellites carried out exploiting a fisher matrix approach @xcite .    in this paper we make a comprehensive comparison of predictions of the jbd theory to current observational data , using wmap and other cmb data plus the galaxy power spectrum as measured by the two - degree field ( 2df ) galaxy redshift survey . we define jbd models in terms of eight parameters , which are allowed to vary simultaneously . our paper is closest in spirit to work by nagata et al . @xcite , who considered a more general model , the harmonic attractor model , which includes jbd as a special case . however their dataset compilation was restricted to the wmap temperature power spectrum . the constraint we will obtain is not competitive with the very stringent solar system bound given above ( though the analysis of ref . @xcite indicates that a limit as high as 3000 might eventually be reached by the measurements of the planck satellite ) , but it is complementary in that it applies on a completely different length and time scale . such constraints can therefore still be of interest in general scalar  tensor theories where @xmath5 is allowed to vary ; for instance nagata et al . @xcite find that in some parameter regimes of the harmonic attractor model the cosmological constraint is stronger than the solar system one . in that regard , our result is most comparable to cosmological constraints imposed on @xmath5 from nucleosynthesis , which give only a weak lower limit of @xmath8 @xcite . the lagrangian for the jbd theory is @xmath9 where the brans  dicke coupling @xmath5 is a constant , and @xmath10 is the brans  dicke ( bd ) field whose present value must give the observed gravitational coupling . we have included factors of @xmath11 to define @xmath12 as dimensionless . the equations for a spatially - flat friedmann universe are @xcite @xmath13 where @xmath14 is the cosmological scale factor , and @xmath15 and @xmath16 are the energy density and pressure summed over all types of material in the universe . the universe is assumed to contain the same ingredients as the wmap concordance model @xcite , namely dark energy , dark matter , baryons , photons and neutrinos . we make the simplifying assumptions of spatial flatness , dark energy in the form of a pure cosmological constant , and effectively massless neutrinos whose density is related to that of photons by the usual thermal argument . the present value of @xmath12 must correctly reproduce the strength of gravity seen in cavendish - like experiments , which requires @xcite @xmath17 where here and throughout a subscript ` 0 ' indicates present value . we will assume that the value of @xmath18 in our solar system is representative of the universe as a whole , though this may not be absolutely accurate @xcite . we also assume that the initial perturbations are given by a power - law adiabatic perturbation spectrum .    when the universe is dominated by a single fluid there are a variety of analytic solutions known @xcite , where @xmath12 is typically constant during a radiation era , slowly increasing during a matter era , and then more swiftly evolving as dark energy domination sets in . however we need solutions spanning all three eras and so will solve the equations numerically , for which we use the integration variable @xmath19 . an example of the evolution is shown in fig . [ bdevolve ] . evolution of the bd field from early in radiation domination to the present . it is just possible to see the evolution of @xmath12 increase as @xmath20 domination sets in . the cosmological parameters are @xmath21 , @xmath22 , and @xmath23 in units of the standard cosmology critical density.,title="fig : " ] +    the basic parameter set we use to build our cosmological models contains the following parameters    = @xmath5   = brans  dicke coupling + @xmath24 present hubble parameter @xmath25 $ ] + @xmath26 baryon density + @xmath27 cold dark matter density + @xmath28 curvature perturbation amplitude + @xmath29 perturbation spectral index + @xmath30 reionization optical depth + @xmath31 galaxy bias parameter , @xmath32    where @xmath33 is the ratio of the ( observed ) galaxy power spectrum to the ( calculated ) matter power spectrum . other parameters are fixed by the assumptions above , and the radiation energy density is taken as fixed by the direct observation of the cmb temperature @xmath34k @xcite . an important subtlety that must be taken into account is that the extra terms in eq . ( [ bdfield ] ) , plus the cavendish - like correction to the present value of @xmath12 , means that the usual relation between the hubble parameter and density , used to define the critical density and hence density parameters , no longer applies . generically , the extra terms require an increase in the present value of @xmath15 to give the same expansion rate , the correction being of order @xmath35 . because of this subtlety , we define the density parameters @xmath36 by dividing by the critical density for the standard cosmology , meaning that the density parameters do nt quite sum to one for a spatially - flat model . operationally , we proceed as follows . we seek a background evolution corresponding to a particular value of @xmath37 and of the present physical matter density . we can assume the initial velocity of the bd field @xmath38 is zero deep in the radiation era , which leaves us two parameters , the early time value of @xmath12 and the value of the cosmological constant , to adjust in order to achieve the required values . this is a uniquely - defined problem , with the necessary values readily found via an iterative shooting method . we carry out the evolution of density perturbations using a modified version of the code defast , based on cmbfast @xcite and originally written to study quintessence scenarios where the dark energy scalar field is minimally @xcite or non - minimally @xcite coupled to the ricci scalar . the architecture of defast is based on the version 4.0 of cmbfast , although there has been a progressive code fork in the subsequent versions . defast takes as input the parameter set described in the previous subsection , and returns the microwave anisotropy spectra ( for temperature and polarization ) and the matter power spectrum . a dynamical and fluctuating scalar field , playing the role of the dark energy and/or the bd field , is included into the analysis together with the other cosmological components , following the existing general scheme @xcite .    in order to bring the model description into the formalism used by defast , we redefine the bd field and coupling according to @xmath39 which brings the lagrangian into the form @xmath40 where @xmath41 is now a canonical scalar field non - minimally coupled to gravity . we implement the cosmological constant in the code by giving @xmath41 a constant potential energy . our calculations include the effect of perturbations , with the initial perturbations in @xmath41 fixed by the requirement of adiabaticity . the correction to the background expansion rate from the dynamics of @xmath41 is the most relevant effect on the cmb power spectrum , appearing as a projection plus a correction to the integrated sachs  wolfe ( isw ) effect , as discussed in detail in ref . @xcite . the data we use are taken from wmap @xcite and the 2df galaxy redshift survey expressed as 32 bandpowers in the range @xmath42mpc @xcite . in order to incorporate the 2df data , the galaxy bias parameter @xmath31 is taken to be a free parameter for which the analytic marginalization scheme of ref . @xcite can be applied . we also consider the effect of including the high-@xmath43 cmb data from vsa @xcite , cbi @xcite , acbar @xcite . our present analysis does not include supernovae data . inclusion of the modification to the luminosity distance from @xmath5 would be straightforward . however the variation of the gravitational coupling @xmath44 means that supernovae can no longer be assumed to be standard candles , and ref . @xcite suggests that the effect from varying @xmath44 dominates . further , inclusion of supernovae data may be particularly susceptible to the possibility that the local value of @xmath12 in the vicinity of the supernova may not match the global cosmological value @xcite . nevertheless , it would be interesting to investigate robust methods for including such data , also in connection with alternative observational strategies @xcite . we carry out the data analysis using the now - standard markov chain monte carlo posterior sampling technique , by modifying the june 2004 version of the cosmomc program @xcite to call defast to obtain the spectra . cosmomc computes the likelihood of the returned model and assembles a set of samples from the posterior distribution . we take full advantage of cosmomc s mpi capabilities by running the code across 19 sun v60x xeon 2.8ghz processors . the metropolis  hastings algorithm is run at a temperature of 1.3 in order to better sample the non - gaussian direction of our posterior distribution which results from the degeneracy between @xmath24 and @xmath45 , both of which have a strong effect on the angular diameter distance . the final chains are then cooled and importance sampled @xcite . it can be noted that for the purposes of posterior sampling , we have parametrized the jbd cosmology using @xmath46 simply because it is more straightforward to obtain the samples we need , while simultaneously suppressing the possibility of jumping to regions with @xmath47 . specifically , we use a flat prior on @xmath48 $ ] where the lower cutoff has been adjusted to the point where the likelihood function is no longer sensitive to the effect of varying the brans  dicke parameter and the @xmath49 model is thereby recovered . as usual , this jeffreys prior , which is defined here as a flat prior on the logarithm of a parameter of unknown scale , has the interesting property of invariance under scale reparametrizations @xcite . for this reason it serves as a reasonable substitute for working with a more desirable physical parameter which could be identified to isolate and give a linear response in the isw effect , mainly responsible for the upper bound on @xmath45 . the optical depth @xmath30 is parametrized using @xmath50}$ ] , where @xmath51 is the fraction of photons that remain unscattered through reionization , since the combination @xmath52 is well constrained by the cmb . the results that we present are based on around 100,000 raw posterior samples , and while the basic constraints can be derived with significantly fewer samples , this large number assures more robust constraints on the derived parameter @xmath5 when we use importance sampling in order to adjust for the change in prior density @xcite . marginalized 1d posterior distributions ( solid lines ) on the base parameters as listed in section ii . also displayed is the mean likelihood of the binned posterior samples ( dotted lines).,title="fig : " ] + turning first to the constraints on the basic parameter set , from figure  [ 1d_constraints ] we note the overall consistency of our results with the current observational picture ( see for example ref . @xcite and a work by two of the current authors ref . @xcite ) , finding the @xmath53 marginalized probability regions to be @xmath54 note that part of our allowed region lies outside the priors assumed by nagata et al . @xcite . as usual for joint analyses of cmb and galaxy power spectrum data , it is unnecessary to impose a further constraint on @xmath24 . marginalized 1d posterior distributions ( solid lines ) on the bd parameter @xmath55 ( upper panel ) . also displayed are the derived importance sampled constraints ( correcting for the change in prior density ) on the more familiar @xmath5 ( lower panel , no smoothing ) . we obtain a @xmath2 marginalized probability bound of @xmath56 , corresponding to a bound on the bd parameter @xmath57 . ]    the primary focus of our study has been to derive constraints on the bd parameter for which , from the outset , we have expected only to find a one - sided bound ; the situation can only become more interesting when both the angular diameter distance and the recombination history become much better probed by the cmb . this expectation is indeed confirmed by the data , as shown in figure  [ 1d_xi_omega ] in which we display the region of highest posterior density . the lower panel detailing the posterior constraint on @xmath5 has been obtained by importance sampling to correct for the change in prior density when changing parameters from @xmath55 to @xmath5 ( we note that the mean likelihood of the binned posterior obtained from sampling @xmath55 performs well for putting a bound on @xmath5 , demonstrating less sensitivity to the details of the prior density ) . we obtain calculate marginalized probability upper bound and the main result of this paper to be @xmath58 the corresponding marginalized probability lower bounds on the bd parameter are found to be @xmath59 this bound is nicely consistent with the expectation for wmap given by the fisher matrix analysis of ref . @xcite . marginalized 2d posterior distribution in the @xmath55@xmath24 plane . the solid lines enclose 95@xmath60 and 99@xmath60 of the probability . under this parametrization there is clearly a geometrical degeneracy.,title="fig : " ] +    we present in figure  [ h0xi ] the 2d posterior constraints in the @xmath55@xmath24 plane , in order to demonstrate the degeneracy and covariance between these two parameters . in a more refined analysis , one could replace @xmath24 with the dimensionless parameter @xmath61 more appropriate to the study of the cmb , where @xmath62 is the sound horizon at recombination and @xmath63 is the angular diameter distance to the last - scattering surface @xcite . finally , in figure  [ cmb_matter ] we display two models , our best - fit @xmath20cdm model with parameters @xmath64@xmath65@xmath66@xmath67 @xmath68 , and a best - fit jbd model with parameters @xmath64@xmath69 , in order to illustrate how the observables change at finite @xmath5 . here the jbd model lies in the vicinity the contour enclosing @xmath53 of the posterior probability distribution and was selected by running a short monte carlo exploration at fixed @xmath70 . note that although in principle the parameter @xmath55 could be extended to @xmath71 , whereby the bulk of the parameter space would be composed of the @xmath49 model , in practice it is reasonable to adjust the lower cutoff to the point where the likelihood function loses sensitivity to the variation of @xmath55 so that the brans  dicke model alone is explored by the mcmc . consequently , the probability contours can reasonably be interpreted to describe the most credible region of the brans  dicke model parameter space .    a comparison between a @xmath20cdm model ( solid line ) and a jbd @xmath20cdm model with @xmath72 ( dashed line ) . the data are the 2df galaxy power spectrum and the models the matter power spectrum convolved with the 2df window functions , and whose overall amplitude is left as a free parameter . detailed parameters are given in section  [ sec : obsconst].,title="fig:",width=283 ] +    our current analysis leaves the bias parameter free , and so constrains only the shape of the matter power spectrum . we note however that the jbd model has a significantly higher amplitude , indeed requiring a modest antibias @xmath73 , which at least in part is due to the more rapid perturbation growth ( @xmath74 during matter domination @xcite ) in the jbd theory . for comparison the @xmath20cdm model has a best - fit bias @xmath75 . this suggests that precision measures of the present - day matter spectrum amplitude , as for instance may become available via gravitational lensing , could significantly tighten constraints . we also note that there is a shift in the location of the baryon oscillations in the matter power spectrum as compared to the @xmath20cdm model ; these are mostly erased by the 2df window function , but future high - precision measurements of those may also assist in constraining @xmath5 . we have carried out the same analysis including also the data from vsa , cbi and acbar in the multipole range @xmath76 . this high-@xmath43 data leads to a slightly tighter bound on the brans  dicke parameter , @xmath77 corresponding to @xmath78 at @xmath2 marginalized probability . however , at the same time inclusion of this new data leads to an unexpectedly large shift in the spectral index , to @xmath79 at @xmath2 marginalized probability , so that the harrison  zeldovich spectrum is only just included ( this statement remains true in the general relativity limit ) . whether this points to some emerging tension in the combined dataset , a harmless statistical fluctuation , or a hint of the breaking of scale - invariance , can be addressed only in the light of the next round of cmb observations . while our constraint on @xmath55 marginalizes over @xmath29 , in the interests of quoting a robust bound we have given as our main result the weaker limit obtained without including the high-@xmath43 data . our ultimate constraint @xmath57 can be compared with that of nagata et al . @xcite , who quote results corresponding to @xmath80 at two - sigma and @xmath81 at four - sigma . the former constraint is much stronger than projected in ref . @xcite , and stronger than one would expect from a nave assessment that the corrections to observables should be of order @xmath35 . if we plotted a model with @xmath82 in our figure  [ cmb_matter ] , it would lie practically on top of the @xmath20cdm model . however their latter constraint is in reasonable agreement with ours , and they do highlight that it is this constraint which corresponds to a sharp ridge of deteriorating chi - squared in their analysis , indicating that their constraint should conservatively be taken as @xmath81 . we have derived a constraint on jordan  brans  dicke gravity from current cosmological observations , including cosmic microwave background ( cmb ) anisotropy data and the galaxy power spectrum data . our main result is to obtain a @xmath2 marginalized probability lower bound on the brans  dicke parameter . this result is complementary to the very strong solar system limit provided by cassini , @xmath83 , as it probes entirely different length and timescales . our analysis is based on a markov chain monte carlo technique varying the basic cosmological parameters and @xmath5 .    at the present precision level , the greatest part of the constraining power comes from the shape of the cmb acoustic peaks , in particular from the first - year observations of wmap . therefore , assuming an extension to four years of the wmap observations , we expect some further improvement on the limit on @xmath5 from cosmology . further help is also expected from other structure formation data , as they improve quality and precision in coming years . in particular we have highlighted that an accurate measure of the present - day matter power spectrum amplitude , for instance from gravitational lensing , may be powerfully constraining when compared to the primordial amplitude from the cmb . a leap forward in this and other contexts is expected from the observations of the planck surveyor probe , to be launched in 2007 . those observations are expected to be cosmic variance limited for the whole spectrum of cmb temperature anisotropy down to the damping tail , and to provide an accurate measurement of the gradient mode of the cmb polarization and its correlation with total intensity up to the sixth acoustic peak @xcite . according to the forecasts of chen and kamionkowski @xcite , the limit on @xmath5 from planck should be around an order of magnitude stronger than that from wmap , and hence vastly stronger than the nucleosynthesis constraint . whether that improvement can be realised from actual planck data , of course , remains to be seen . v.a .  was supported at sussex by a marie curie fellowship of the european community programme human potential under contract hpmt - ct-2000 - 00096 , c.b . and f.p . in part by nasa ltsa grant nng04gc90 g , s.m.l .  in part by the european union cmbnet network at geneva , and a.r.l .  by pparc . s.m.l .  thanks the university of geneva for hosting the computations in this work , and c.b .  and s.m.l . acknowledge visits to sussex supported by pparc . we thank john barrow for useful discussions , and acknowledge the use of the legacy archive for microwave background data analysis ( lambda ) . support for lambda is provided by the nasa office of space science . et al . _ , * 148 * , 161 ( 2003 ) , astro - ph/0302213 ; l. verde _ et al . _ , astrophys . j. supp . * 148 * , 195 ( 2003 ) , astro - ph/0302218 ; g. hinshaw _ _ , astrophys . j. supp . , * 148 * , 135 ( 2003 ) , astro - ph/0302217 .                  a. kosowsky , m. milosavljevic , and r. jimenez , phys . rev . d*66 * , 063007 ( 2002 ) ` astro - ph/0206014 ` ; r. jimenez , l. verde , h. peiris , and a. kosowsky , phys . rev . d*70 * , 023005 ( 2004 ) , ` astro - ph/0404237 ` .
Abstract:

Answer: we use cosmic microwave background data from wmap , acbar , vsa and cbi , and galaxy power spectrum data from 2df , to constrain flat cosmologies based on the jordan  brans  dicke theory , using a markov chain monte carlo approach . using a parametrization based on @xmath0 , and performing an exploration in the range @xmath1 $ ] , we obtain a @xmath2 marginalized probability bound of @xmath3 , corresponding to a @xmath2 marginalized probability lower bound on the brans  dicke parameter @xmath4 .


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: strong gravitational lensing of background objects by massive clusters of galaxies is expected to produce multiple images with wide separations , @xmath0 , and the multiple images of background lensed galaxies have long been known with radii of giant arc systems up to @xmath1 ( e.g. @xcite ) . however , despite a number of searches , primarily in large - area radio surveys ( e.g. @xcite ) , widely - separated multiple images of qsos have not yet been found . this is attributable to the low probability of strong lensing along any given line of sight coupled with the relatively low numbers of known qsos . the widest separation qso lensed systems have separations @xmath2 @xcite . other wide - separation multiple qso systems found to date are more likely to be `` binary qsos '' ( that is , a pair of active black holes cohabiting the same parent host galaxy ) rather than multiple images of the same qso @xcite .    we have searched a catalogue of 22,163 qsos that comprises the recently completed 2df qso survey ( `` the 2qz survey '' ) for candidate wide - separation lensed qso systems and have obtained higher signal - to - noise spectra in an attempt to eliminate obvious non - lensed systems and to identify systems that are likely to be genuine lensed systems . after an incomplete survey of the most likely candidates , there remain six systems which , with varying degrees of confidence , are likely to be lensed . the 2qz survey covers 740 deg@xmath3 in two strips of sky , one centred on the south galactic pole region and the other on the celestial equator in the region of the north galactic pole . the survey photometry is described by @xcite and the initial release of the first 10,000 qsos is described by @xcite . in constructing the 2qz survey , candidate qsos were first selected from @xmath4 , @xmath5 and @xmath6 photographic photometry of morphologically - stellar objects , and objects whose colours lay bluewards of the main stellar locus were selected for low resolution spectroscopy at the two - degree field ( 2df ) facility at the anglo - australian telescope . the initial spectra obtained were of sufficient quality that reliable identifications of qsos and measurements of their redshifts could be made , with an r.m.s . error in redshift measurement of typically 0.0035 @xcite , although it should be noted that the redshift error for an individual object depends on its brightness and emission - line equivalent widths as well as observing conditions . the absolute accuracy of the photometry in each band is estimated as @xmath7 magnitudes at the faint survey limit @xcite . the relative photometric accuracy for objects measured close together on the same photographic plate is better than this , however : @xcite found that , for objects more than two magnitudes brighter than the completeness limit of the photographic plate , relative photometry accurate to @xmath8 magnitudes was typical , with that value rising rapidly to values @xmath9 magnitudes as the plate limit is approached . the resulting catalogue of photographic photometry and spectroscopically confirmed qsos was then searched for pairs or multiple systems whose colours were consistent given the large photometric uncertainties and whose redshifts were consistent given the typical redshift uncertainty . to determine the colour consistency , it was assumed that two qsos , if lensed , should have the same @xmath4 , @xmath5 and @xmath6 magnitudes apart from a magnitude offset arising from differential magnification of the images in a lensed system , and thus assuming negligible effects due to variations in qso colours with time ( see section  [ timedelay ] ) . the statistic used as the measure of colour consistency was @xmath10 where @xmath11 is the difference in apparent magnitude of the two qsos measured in waveband @xmath12 and @xmath13 is the average magnitude difference for that qso pair . redshift consistency was measured by the difference , @xmath14 in 2qz catalogue redshift of two qsos . candidate multiple systems were initially selected if they satisfied the criteria @xmath15 , @xmath16 and angular separation @xmath17 . the colour criterion is close to the expected scatter in this statistic , given the photometric errors at the survey limit , and hence we expect the initial colour selection to select only @xmath18% of the lensed systems . the redshift criterion is a factor two larger than the expected scatter in @xmath14 and this criteria should allow the majority of lensed systems to be selected . the stringent colour criterion was chosen in this initial study in order to maximise the probability of finding lensed systems , and hence further observations should be carried out to establish the true completeness of the selected sample . a further check was made on the consistency of their radio properties by searching for detections in the nvss @xcite . one candidate pair with clearly discrepant radio detections was removed from the sample of candidates . the nvss properties of two remaining pairs are discussed further in section 3 . candidates with @xmath19 were selected from 12307 qsos in the sgp region and from 9856 qsos in the ngp region of the 2qz survey , and 38 pairs of qsos in total met these criteria . of the qsos with confirmed redshifts , there were no candidate systems selected with more than two members , although some of the pairs do have additional @xmath20-selected companions for which 2qz spectroscopic observations either were not carried out or did not lead to a positive identification . our expectation is that the majority of the selected pairs will not prove to be lensed systems . test samples were also constructed in attempt to assess the possible contamination . first , qsos pairs that are more widely separated , but still within one degree of each other and measured on the same photographic plates , were selected with @xmath21 and the distribution of @xmath22 was compared with the @xmath22 distribution of the candidate pairs . no statistically - significant difference was found , although this is probably a reflection of the large errors in @xmath22 arising from the large photographic photometry errors . a second test sample comprising qso pairs with @xmath23 and @xmath16 was selected and the distribution of angular separations , normalised by the total number of pairs selected with @xmath24 , was compared with the candidate lens sample selected to have @xmath15 . this time there was an excess of pairs on arcminute scales in the sample of candidate lensed systems consistent with a minority of the pairs being lensed systems , but because of the small numbers of pairs selected the test was not statistically significant . hence at this stage there is no evidence either for or against the hypothesis that the selected sample contains a subset of lensed systems . further evidence needs to be collected on all members of the sample in order to establish its completeness and contamination .    a key part of that evidence is the follow - up spectroscopy described in the next section . because of limitations of observing time and observing conditions , only a subset of the candidates have had such observations , and only that subset is described in the remainder of this paper . for prioritising the observations , the pairs were ranked in @xmath25 and also in @xmath22 : a composite rank was produced by multiplying these ranks and the top candidates in each region were observed . in the sgp region , 10 pairs were observed , but the spectroscopic follow - up is particularly sparse in the ngp region , where only for one pair have spectra been obtained of sufficient quality to allow some assessment of whether or not it may be lensed . at this stage the sample reported here can not be regarded as either being complete or free from contamination by non - lensed pairs of qsos . we postpone any discussion of the completeness of the sample to later work . we anticipate that in future the candidate qso pairs in the ngp region will also have the benefit of more accurate , five - band , colours from the sloan digital sky survey @xcite , which should lead both to samples of candidates with lower contamination and with better quantified photometric uncertainties .    in the following sections each pair is designated by the iau naming convention applied to the mean position of the pair , in the fk5 ( equinox j2000 ) coordinate system . [ cols="^,^,^,^,^,>,^,^,^,^,^,^",options="header " , ]     [ edr ] if these systems of paired qsos are indeed lensed qsos , then the lenses should be identifiable as clusters of galaxies . as described in section 2 , one of the possible lensed systems , j@xmath26 , appears associated with a galaxy cluster at @xmath27 . ten galaxies in the vicinity have redshifts , measured by the 2df galaxy redshift survey , that are close to this value , with a velocity dispersion of about 600  km  sec@xmath28 ( fig .  [ fchart ] ) . such a velocity dispersion would however produce an image splitting of only about @xmath29 for an isothermal sphere lens , so if this cluster is the lens then either an anisotropic or a more centrally concentrated mass distribution is needed . the total mass required to produce the observed image splitting is @xmath30 m . there is also a weak rosat - detected x - ray source close to the position of the two bright galaxies which lie along the line joining the two qsos , although the origin of the emission is not known . at present it is not clear if sufficient mass is present in this cluster to produce the observed image splitting . one other system , j0150@xmath312921 , may also have a low - redshift cluster associated with it , again based on redshift measurements from the 2df galaxy redshift survey , although in this case only three redshifts are currently known ( fig .  [ fchart ] ) .    in order to estimate the expected distribution of lens redshifts we consider the calculations presented by @xcite in which the lensing clusters are assumed to have spherically - symmetric @xcite mass profiles , and to have an abundance predicted by the @xcite evolving mass function . here we assume a cold dark matter power spectrum @xcite modified to include an approximate correction for baryons @xcite in a @xmath32cdm cosmology with @xmath33 and @xmath34 and a power - spectrum normalisation given by @xmath35 . in this model , it is found that the distribution of lens redshifts is not sensitive to the choice of cosmological parameters , although the abundance of lensing systems is . magnification bias is calculated as discussed by @xcite .    within this model , it is clear that a lensing cluster is expect to have @xmath36 , and for the widest separations the model predicts that any lens should have @xmath37 . however , the total probability predicted for @xmath38 lensed systems is so low that such systems are expected only if there are significant departures from the assumptions of the model , in which case the lens redshift distribution would need to be recalculated . at this stage it can not be regarded as conclusive that these candidate systems are indeed lensed , and we also do not yet have an accurate assessment of the completeness of the lensing search . nonetheless it is interesting to compare the results presented here with the expectations of lensing models . @xcite and @xcite have calculated the expected distribution of angular separations arising from strong lensing by massive clusters of galaxies in the standard cosmological framework that assumes structure grows hierarchically from initial gaussian perturbations in a universe dominated by cold dark matter and dark energy . with cosmological parameters as determined from the wmap observations @xcite the expected number of systems with @xmath39 in the 2qz survey of @xmath40 qsos is about 0.5 @xcite , so it is possible , although not highly likely , that the 2qz survey contains wide - separation lensed systems . however , the distribution is expected to be dominated by angular separations towards the lower end of the range , and in this standard model the expected number of systems with @xmath41 is only @xmath42 @xcite . one possible explanation is that the widest - separation systems presented here are not in fact lensed , despite the similarity of the redshifts and spectra . a second is that the values of cosmological and cluster model parameters are not as assumed in calculating the above numbers : key parameters are @xmath43 , @xmath44 , @xmath45 and also parameters of the model such as the dependence of cluster halo concentration on mass , redshift and cosmology and the cosmic scatter in halo concentration ( for details see @xcite ) . or it may be that the models do not yet comprise the full complexity required to model accurately the statistics : possible deficiencies include the assumption of the @xcite mass profile ( see also @xcite ) and the assumption of a smooth distribution of mass . other departures from the standard cosmological model , such as the existence of non - gaussian fluctuations , could also modify the lensing probabilities .    for these reasons we believe it is imperative to test carry out observational tests of the candidate lensed systems presented here in order to test these models , rather than using the predictions of the models to argue against the interpretation of these systems as being lensed . previous searches have been undertaken for lensed qsos with comparable separations , and it is interesting to investigate whether the detection of any systems in 2qz would be inconsistent with those searches . @xcite searched a sample of size @xmath46 radio - loud qsos without finding any evidence for multiple images in the separation range @xmath47 . @xcite examined a complete sample of 11,670 flat - spectrum radio sources for possible multiple images separated in angle by @xmath48 and found no lensed systems . an incomplete sample of @xmath49 radio sources yielded one possible lensed system . the first step in comparing these results with those from the 2qz survey is to compare sample selection . the range of angular separations probed by the radio surveys is generally smaller than considered here , although had any 2qz pairs with @xmath39 been lens candidates they should have been identified as such . for comparison purposes we shall use the cluster - lensing model to compare surveys covering different ranges in separation , although we have argued above that the overall lensing probabilities predicted by such models may be too low to explain the existence of more than one lensed system in the 2qz survey . the model of @xcite with parameters @xmath50 , @xmath34 , @xmath35 and @xmath51 predicts relative lensing probabilities of @xmath52 , @xmath53 and @xmath54 for the separation ranges @xmath55 , @xmath56 and @xmath57 respectively , assuming the same magnification bias in each case . if the models are indeed deficient in some way , or if the values of parameters are different from that assumed , these relative probabilities will be different , but they suffice for the discussion in this section .    in fact , because of the double radio and optical selection of the first survey @xcite the magnification bias could be a factor 5 smaller for that survey compared with the jvas / class or 2qz surveys ( @xcite ) , so the overall expected relative lensing probabilities in the three surveys are @xmath58 , @xmath53 and @xmath54 for first , jvas / class and 2qz respectively . we also need to make some assumption about sample completeness and to multiply the probabilities by the sample sizes . here we shall adopt a conservative view of the jvas / class result and say that out of 15,000 qsos no lensed systems were found , and we adopt an assumed completeness of 50% for the 2qz survey . at this stage this figure is illustrative only . with these assumptions , the predicted relative numbers of systems in each survey are 0.1 , 1.6 and 1.0 for the first , jvas / class and 2qz surveys respectively . now suppose the entire dataset contains one lensed system : the _ a posteriori _ probability of that system being found in the 2qz survey is 37% . however , if there are more systems in the 2qz survey , the probability of there being none in the other surveys falls to 13 , 5 , 2 % for 2 , 3 or 4 2qz systems . thus , if the surveys are indeed comparable in this way , the existence of two lensed systems within the 2qz survey would not be inconsistent with the other surveys , but a larger number could be . finally , however , we note that another possibly significant factor in calculating the expected lensing probabilities is the redshift distribution of the source population . @xcite , @xcite and @xcite have demonstrated that cluster lensing probability strongly increases with the mean redshift of the source population , and also depends to some extent on the distribution of source redshifts about the mean . the redshift distribution is not well - known for the class and first surveys , but appears likely to have a median value in the range @xmath59 ( @xcite ) . the 2qz median redshift is 1.5 . lensing probability increases by a factor about 3 as the source redshift is increased from @xmath60 to @xmath61 @xcite , so if the typical redshift of sources in the radio surveys is even slightly lower than that of qsos in the 2qz survey this would also lower the relative probability of discovering lensed systems in the radio surveys . our overall conclusion from this comparison is that , observationally , the results from the 2qz and radio surveys may not be incompatible , but clearly further work needs to be done both to establish more accurately the redshift distribution of the jvas / class sources and to establish the reality of the 2qz candidate systems . evidence has been presented for the discovery of the first lensed qsos with multiple images separated on arcminute scales . if these are genuine lensed systems then we should expect the associated lenses to be observable as massive clusters of galaxies . further evidence needs to be collected to establish that these are indeed lensed , however . one such piece of evidence that is independent of assumptions of the nature of the lens would be the existence of common ly-@xmath62 forest absorption in candidate lensed qso spectra . the existence of giant arcs associated with a system would also provide strong evidence for cluster lensing , especially if the host galaxies of the qsos could be detected and shown to be distorted by lensing . detection of the lens itself , either by direct imaging or weak - lensing , would also provide convincing evidence . consideration of the expected distribution of angular separations indicates that , although a number of the identified systems may be lensed , a number of those with the greatest separations may not be . confirmation or discounting of the widest separated systems is particularly important as any such systems would challenge our current understanding of lensing by galaxy clusters .    among the examples found , one candidate lensed system has been found to display broad civ absorption in one image and not the other , which could be due to the formation of such absorption systems on timescales of tens of years . further searches for wide - separation lenses may be partially incomplete because such absorption systems would affect optical colours , spectra and the flux in other wavebands , including possibly the x - ray band . broad - line region variability may also render search methods based on colour and/or spectroscopic comparison partially incomplete unless generous selection criteria are allowed . finally , it may be that optical searches are biased against lenses with separations less than a few tens of arcsec because of obscuration in galaxies associated with the cluster lens . * acknowledgments * + the 2df qso survey was based on observations made with the u.k . schmidt and anglo - australian telescopes . we thank all the present and former staff of the anglo - australian observatory for their work in building and operating the 2df facility . follow - up spectroscopic observations were carried out at the anglo - australian telescope and at the william herschel telescope . the william herschel telescope is operated on the island of la palma by the isaac newton group in the spanish observatorio del roque de los muchachos of the instituto de astrofisica de canarias . we acknowledge the use of the sloan digital sky survey archive for the photometry presented in table  [ edr ] . we are also grateful to the 2df galaxy redshift survey team for providing redshifts of galaxies in the fields shown here , in advance of publication , and would like to thank raylee stathakis for useful discussions and aidan crook for assistance with observing preparation . aml acknowledges the support of the portuguese fundao para a cincia e a tecnologia . 99 bardeen , j.m . , bond , j.r . , kaiser , n. & szalay , a.s . apj 304 , 15 ( 1986 ) . colless , m. et al . , astro - ph/0306581 ( 2003 ) . condon j.j . , cotton w.d . , greisen e.w . , yin q.f . , perley r.a . , taylor g.b . & broderick j.j . 115 , 1693 ( 1998 ) . croom s.m . , smith r.j . , boyle b.j . , shanks t. , loaring n.s . , miller l. & lewis i.j . mnras 322 , l29 ( 2001a ) . croom s.m . , shanks t. , boyle b.j . , smith r.j . , miller l. , loaring n.s . & hoyle f. mnras 325 , 483 ( 2001b ) . dodorico v. , cristiani s. , dodorico s. , fontana a. , giallongo e. & shaver p. , a. & a. 339 , 678 ( 1998 ) . green p.j . , kochanek c. , siemiginowska a. , kim d - w . , markevitch m. , silverman j. , dosaj a. , jannuzi b.t . & smith c. apj 571 , 721 ( 2002 ) . kochanek c.s . , falco e.e . & muoz j.a . , apj 510 , 590 ( 1999 ) . kochanek , c.s . , falco , e.e . , impey , c. , lehar , j. , mcleod , b. & rix h .- w . castles web site : http://cfa-www.harvard.edu/castles ( 2002 ) . kneib , j .- p . , ellis , r.s . , smail , i. , couch , w.j . , sharples , r.m . , apj 471 , 643 ( 1996 ) . li , l .- x . & ostriker , j.p . , apj 566 , 652 ( 2002 ) . lopes , a.m. & miller , l. , astro - ph/0306174 ( 2003 ) . ma , feng , mnras 335 , l99 ( 2002 ) . maoz d. & rix h.w . , apj , 416 , 425 ( 1993 ) . maoz d. , rix h.w . , gal - yam , a. , gould , a. , apj , 486 , 75 ( 1997 ) . mitchell p.s . , ph.d thesis , university of edinburgh ( 1989 ) . navarro , j.f . , frenk , c.s . & white , s.d.m . , apj 490 , 493 ( 1997 ) . ofek e.o . , maoz d. , prada f. , kolatt t. & rix h - w . , mnras 324 , 463 ( 2001 ) . ofek e.o . , rix h - w . , maoz d. & prada f. , mnras 337 , 1163 ( 2002 ) . oguri , m. , lee , j. & suto , y. , astro - ph/0306102 ( 2003 ) . phillips p.m. , browne i.w.a . & wilkinson p.n . , mnras 321 , 187 ( 2001 ) . phillips p.m. , browne i.w.a . , jackson n.j . , wilkinson p.n . , mao s. , rusin d. , marlow d.r . , snellen i. & neeser m. , mnras 328 , 1001 ( 2001 ) . sahu k.c . , shaw , r.a . , kaiser , m.e . , baum , s.a . , ferguson , h.c . , hayes , j.j.e . , gull , t.r . , hill , r.j . , hutchings , j.b . , kimble , r.a . , plait , p. & woodgate , b.e . , apj 492 , l125 ( 1998 ) . sheth , r.k . & tormen , g. mnras 308 , 119 ( 1999 ) . smail , i. , couch , w.j . , ellis , r.s . & sharples , r.m . , ap.j . 440 , 501 ( 1995 ) . small t.a . , sargent w.l.w . & steidel c.c . 114 , 2254 ( 1997 ) . smith r.j . , croom s.m . , boyle b.j . , shanks t. , miller l. & loaring n.s . , submitted to mnras ( 2002 ) . spergel , d.n . , verde , l. , peiris , h.v . , komatsu , e. , nolta , m.r . , bennett , c.l . , halpern , m. , hinshaw , g. , jarosik , n. , kogut , a. , limon , m. , meyer , s.s . , page , l. , tucker , g.s . , weiland , j.l . , wollack , e. & wright , e.l . , apj supplement series 148 , 175 ( 2003 ) . stoughton c. et al . 123 , 485 ( 2002 ) . sugiyama , n. apj supplement 100 , 281 ( 1995 ) . wambsganss , j. , bode , p. & ostriker , j.p . , astro - ph/0306088 ( 2003 ) . zhao , d. , jing , y.p . & brner , g. , apj 581 , 876 ( 2002 ) .
Abstract:

Answer: we report the possible discovery of multiple gravitationally - lensed images of qsos with angular separations on arcminute scales . the qsos were selected from the completed 2df qso survey as having redshifts and optical colours which are statistically consistent . in this paper we present higher - quality optical spectra of the candidates and discuss the likelihood of their genuinely being systems lensed by massive clusters of galaxies . from a comparison of the spectra it appears that up to six pairs of qsos may be lensed multiple images , although the true number may be less than that and further observations should be undertaken to amass more evidence and to detect the lensing clusters . two of the candidates may be associated with low redshift clusters of galaxies . [ firstpage ]    gravitational lenses  qsos .