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: asking for the time evolution of a system means solving the schrdinger equation @xmath0 . restricted to the case of closed systems - which is done from now on - this equation can formally be solved by @xmath1 with @xmath2 for the sake of simplicity and @xmath3 . the whole dynamics is thus governed by the unitary propagator @xmath4 which maps any given `` initial state '' @xmath5 at time @xmath6 to its propagated counterpart @xmath7 at time @xmath8 . quantum mechanical determinism is expressed by the fact that @xmath5 can also be regarded as a `` final state '' into which a state @xmath9 evolves by satisfying @xmath10 ( see fig . [ fig : evol1 ] ) . @xmath11{images / evol1.eps}\ ] ]    since equation  ( [ prop1 ] ) is valid for all times @xmath12 , giving @xmath13 , one is led to the operator identity @xmath14 in other words : the set of evolution operators @xmath15 forms a unitary group @xmath16 with the group property  ( [ group1 ] ) , giving any propagator @xmath17 its inverse element @xmath18 . the physical meaning of this mathematical statement is profound : in quantum mechanics there is no microscopic arrow of time . any state @xmath19 can be propagated to the future as well as to the past , so that there is a complete time symmetry .    in order to explain an arrow of time despite these facts , there are several ideas a few of them are :    1 . any measurement on the system interrupts the deterministic evolution of the state by applying random projections ( plus normalization ) instead of unitary operations ( the `` collapse of the wave function '' ) . this process can not be reversed . real physical systems can never be regarded as closed . there is always an interacting environment with many unobservable degrees of freedom , inducing non - unitary time evolution of the observed system . 3 .   time symmetry in classical mechanics is only broken with respect to macroscopic properties of a system ( e.g. temperature or entropy ) . considering macroscopic observables in quantum mechanics , this effect also shows up ( coarse graining argument ) . all these ideas have been worked out to elaborate physical models with strong experimental support and all of them break time symmetry ( see e.g.  @xcite ) . so one could be content with the situation and regard those models as fairly good explanations of an arrow of time , which is indeed an undeniable part of human experience . a group of scientists around _ ilja prigogine _ ( nobel prize in chemistry 1977 ) , referred to as the `` brussels school '' , is not satisfied with the situation and has worked out a quantum mechanical formalism especially applying to unstable systems ( i.e. systems involving decay ) , where time symmetry is claimed to be broken at the microscopic level ( see e.g.  @xcite ) . @xmath11{images / evol2.eps}\ ] ]    this fact is expressed by the splitting of the unitary group @xmath16 into two semigroups @xmath20 and @xmath21 of propagators for negative and positive times , such that past and future states can no longer be regarded as equivalent counterparts . this leads to a so - called `` intrinsic irreversibility '' , implying a microscopic arrow of time ( see fig . [ fig : evol2 ] ) . this paper is separated into two different parts . the first part is dedicated to a review of the rigged hilbert space formalism worked out by the brussels school . it is shown how it works , how it applies to physical systems and what its advantages over standard quantum mechanics are . the second part is devoted to the role of time symmetry within the rigged hilbert space formalism . apart from its practical advantages , the brussels school points out that a microscopic arrow of time can be implemented here . it will be shown that such an implementation is possible but not mathematically necessary . writing down the formal solution  ( [ prop1 ] ) of a closed system s schrdinger equation is not yet really satisfying . how can expectation values of observables be numerically calculated ? to this aim one has to find a spectral decomposition of the system s hamiltonian @xmath22 , thus one has to solve the time - independent schrdinger equation @xmath23 . since @xmath22 is a self - adjoint operator , its spectral decomposition is complete and orthogonal and may therefore be used as a spectral decomposition of the unity operator @xmath24 . furthermore the spectrum of @xmath22 is real and consits in general of a discrete part @xmath25 and a continuous part @xmath26 . in special cases there are other parts , e.g. singular continuous and dense spectra , but these cases are excluded here . with this restriction , using the dirac bra - ket notation , the spectral decomposition of the unity operator reads in general @xmath27 where @xmath28 , @xmath29 and @xmath30 is the energy density . as soon as the energy decomposition of the initial state @xmath31 , generally given by @xmath32 is known , the propagated state is also known and reads @xmath33 the expectation value of an observable @xmath34 can be calculated for any time @xmath8 using @xmath35 so the fundamental key to the dynamics of a system is the spectral decomposition of its hamiltonian . all dynamical properties are monitored by this decomposition and here is where the new formalism enters . there is a class of systems where the brussels schools s spectral decomposition of the hamiltonian differs from the usual decomposition : systems envolving decay , called _ unstable systems_.      the basic concept of the brussels school s formalism is to find a different spectral decomposition of the hamiltonian . in contrast to the method of _ complex scaling _ ( see  @xcite ) , the hamiltonian is not modified and remains self - adjoint ; it only acquires a different representation involving generalized eigenvectors contained in an extended distribution space . the modified decomposition can most easily be achieved by using the hamiltonian s _ resolvent _ @xmath36 the spectrum of @xmath22 is defined by the irregularities of its resolvent . the discrete eigenvalues are the singularities of @xmath37 , i.e. those points in the complex plane where the operator @xmath38 is not everywhere defined on the hilbert space @xmath39 . the continuous eigenvalues are those points where @xmath38 is an unbound ( and thus non - continuous ) operator on @xmath39 . anywhere else on the complex plane the resolvent is an analytic function mapping complex numbers onto the set of linear bounded operators . as @xmath22 is self - adjoint , any irregularities of @xmath38 take place on the real axis , where the discrete eigenvalues form the point spectrum and the continuous eigenvalues define a cut in the complex plane ( see fig . [ fig : spectrum ] ) . @xmath40{images / spectrum.eps}\ ] ]    it is possible to find an analytic continuation of @xmath38 from both sides of the cut , each leading to a different operator @xmath41 with @xmath42 understood here and in the following . note that this continuation is not possible at singularities embedded in the continuous spectrum . note further that on the real axis outside the continuous spectrum both operators @xmath43 and @xmath44 coincide with @xmath38 . the two continued resolvents are known as the _ retarded _ and the _ advanced _ _ green operator _ , respectively . the resolvent can be used to construct eigenvectors . the residuum of the resolvent at a discrete eigenvalue yields the dyadic product of the corresponding eigenvector : @xmath45 note that both dyadic products project onto the same subspace corresponding to @xmath46 . if the discrete eigenvalue is isolated , i.e. not embedded in the continuous spectrum , both limits from above and below the real axis coincide . in the following we will assume that all discrete eigenvalues are isolated and we therefore have @xmath47 a discrete eigenvector @xmath48 of @xmath22 can be obtained by applying the residual resolvent to any hilbert vector @xmath31 being not orthogonal to the eigenspace of @xmath46 . : @xmath49 where @xmath50 is an adequate normalization constant , such that @xmath51 . in analogy to the discrete case the residuum of the resolvent at one of the continuous eigenvalues yields the dyadic product of the corresponding continuous eigenvector : @xmath52 note that both dyadic products are equivalent representations of the @xmath53-operator @xmath54 , defined by @xmath55 where @xmath56 is the continuous part of the total hamiltonian @xmath22 . since the continuous eigenvalues are not isolated , there are always pairs of eigenvectors . they correspond to retarded and advanced solutions of the continuous eigenvalue equation and are known as _ incoming _ and _ outgoing waves_. both solutions can equivalently be used in the spectral decomposition and they can be constructed by @xmath57 where @xmath58 is a suitable normalization constant , such that @xmath59 .    using dyadic products , the retarded and advanced green operators can be written as @xmath60 where @xmath61 is any of the two dyadic products @xmath62 and the principal value distribution @xmath63 is defined by @xmath64 the difference of the two green operators gives the dyadic product of the continuous eigenvector : @xmath65 altogether , the spectral decomposition of the unity operator may be expressed by the exclusive use of the resolvent : @xmath66 this unity decomposition will be used as the starting point of the brussels schools s unity spectral decomposition . the resolvent can be analytically continued to the complex plane beyond the real axis entering the _ second riemann sheet _ : @xmath67 on the one half of the complex plane where the continuation started ( the first riemann sheet ) , both operators coincide with the original resolvent . on the other half ( the second riemann sheet ) it might happen that new singularities appear , and that in a pairwise manner : if @xmath68 is a singularity of @xmath69 in the lower half of the complex plane , then @xmath70 is a singularity of @xmath71 in the upper half ( see fig . [ fig : riemann ] ) . are exactly those with a resolvent having complex poles on the second riemann sheet .    @xmath72{images / riemann.eps}\ ] ]    the second riemann sheet poles act as complex resonances and cause the decay of certain states . consider the continuous part of the unitary decomposition  ( [ unity ] ) : @xmath73 the following manipulations of @xmath74 are reducing its domain . as a consequence one is forced to use a smaller _ test space _ , enlarging at the same time the distribution space and inducing this way a certain _ gelfand tripel _ or _ rigged hilbert space _ where the modified completeness relation remains valid . we delay this discussion to section  [ rhs ] . if one analytically deforms the integration path along the continuous spectrum @xmath26 to a path @xmath75 leading through the lower half of the complex plane , one enters the second riemann sheet and crosses the singularities of @xmath69 . let the set of these singularities be defined by @xmath76 and be enclosed by @xmath77 . in the course of the deformation circular integrations around the singularities are split off , each leading to a residuum of @xmath69 at these points ( see fig . [ fig : spectrum2 ] ) . @xmath40{images / spectrum2.eps}\ ] ]    the modified operator now reads @xmath78 introducing a _ complex distribution _ @xmath79 , defined by @xmath80 where @xmath81 and @xmath82 being the analytic continuation of @xmath83 to respectively the upper and lower half of the complex plane , one can rewrite @xmath84 as @xmath85 in analogy to the standard case one defines dyadic products of _ generalized eigenvectors _ of @xmath22 by @xmath86 for the discrete case and , using  ( [ dyad2 ] ) , @xmath87 for the continuous case , where @xmath61 is any of the two dyadic products @xmath62 . hence we have @xmath88 the generalized eigenvectors appear as pairs of right and left eigenvectors . the discrete eigenvectors are called _ gamov vectors _ and can be constructed by : @xmath89 where @xmath31 is a suitable vector and @xmath90 is a normalization constant , such that @xmath91 . by adjunction one sees that @xmath92 . furthermore we have @xmath93 . these objects are left and right eigenvectors of the hamiltonian @xmath22 corresponding to complex eigenvalues : @xmath94 since the gamov vectors have zero norm , they all coincide in hilbert space with the null vector , which is not allowed to be an eigenvector . also , since @xmath22 is self - adjoint , it can not have any eigenvectors corresponding to complex eigenvalues . concluding , the gamov vectors can not be elements of the hilbert space . in fact they are distributions in a suitably chosen distribution space ( see section  [ intrinsic ] ) . the continuous eigenvectors can be constructed in two ways : @xmath95 by application to test vectors with analytic extensions to the complex plane one sees that @xmath96 hence there are two equivalent bases given by @xmath97 both bases can equally be used to span the hilbert space @xmath39 . arbitrarily we select @xmath98 and obtain a _ retarded _ unity decomposition @xmath99 in complete analogy one can modify the spectral decomposition by deformation of the integration path to @xmath81 in the _ upper _ half of the complex plane , entering the domain of @xmath71 with singularities at @xmath100 . one then obtains the two bases @xmath101 such that @xmath102 arbitrarily choosing @xmath103 one obtains an _ advanced _ unity decomposition : @xmath104 which is just the adjoint of the retarded unity decomposition , so @xmath105 . note , however , that both decompositions are equivalent representations of the unity operator , thus in the sense of an operator identity we have @xmath106 . the unity decompositions obtained in the preceding section can be used for different representations of vectors and operators . using the retarded decomposition the hamiltonian reads @xmath107 the adjoint of the above expression yields the hamiltonian in the _ advanced _ decomposition . since both decompositions are equivalent representations , one sees that @xmath108 , i.e. the hamiltonian remains self - adjoint . the propagator in the retarded decomposition reads @xmath109 since the complex numbers @xmath68 are located in the lower half of the complex plane , the contribution from the gamov vectors decays in the future . if an initial hilbert state @xmath31 has gamov components , they will disappear for @xmath110 . so if the state has no rotating components , i.e. @xmath111 , the whole state will disappear in the future with the decay rates of the gamov contributions given by the imaginary part of the contributing @xmath68 : @xmath112 with @xmath113 , @xmath114 and @xmath115 . each contribution from a gamov state also induces a certain energy shift , given by the real part of the @xmath68 . the advantage of the extended spectral decomposition is getting clear now : not only the stable frequencies can immediately be read from the spectrum , but the decay rates and the energy shifts , too . by additionally checking the amplitudes @xmath116 between the initial state and the gamov vectors one has immediate access to the dissipative dynamics of the system . apart from the contribution of the gamov vectors one has to consider the contribution from the generalized continuous eigenvectors . let us call it the _ background contribution_. since the decomposition involves the complex distribution @xmath117 , a test vector @xmath118 has to be introduced . the background contribution then reads : @xmath119 where @xmath120 and @xmath121 . due to @xmath122 the integration path is transformed into the curve @xmath75 below the complex eigenvalues @xmath68 . @xmath123 where @xmath124 and @xmath125 are the analytic continuations of respectively @xmath126 and @xmath127 to the lower half of the complex plane ( provided their existence ) . the curve @xmath75 is arbitrary as long as @xmath128 surrounds @xmath129 , i.e. , @xmath75 is placed _ below _ the poles @xmath68 , which all have negative imaginary parts . for @xmath130 , due to the exponential factor @xmath131 , the background contribution thus becomes more and more neglectable in comparison with the pole contributions . for initial states with high decay rates , i.e. pole contributions with small imaginary parts , the background contribution may be neglected already for small times @xmath130 . this approximation is called the _ pole approximation _ and it is essentially equivalent to the _ wigner - weisskopf approximation _ ( see  @xcite ) . the time evolution of the initial state may be approximated by an exponential decay and its energy spectrum by a breit - wigner distribution . this approximation is very popular , e.g. in quantum optics , and indeed in most cases it is a very good one . obviously the above argument does not hold for negative times . here , the contributions from the gamov vectors diverge . in several publications of the brussels school ( see e.g.  @xcite or _ a. bohm _ in  @xcite ) this is an argument in favour of a broken time symmetry , since for negative times the _ advanced _ decomposition  ( [ unity_ad ] ) is the better choice . however , as will be shown below , this is a rather intuitive reasoning and does not enforce a break of time symmetry . let us see what happens with the retarded decomposition of a hilbert state for negative times . of course the gamov contributions now diverge , but they are compensated by the background contribution . for @xmath132 the background integral  ( [ bg2 ] ) gains influence , since the integration path @xmath75 is placed _ below _ the complex eigenvalues @xmath68 . the same argument , that for increasing @xmath130 the background contribution can be neglected in comparison with the pole contributions , is turned upside down for @xmath133 . the smaller imaginary part of @xmath75 induces a higher divergence than the contribution of the poles . by construction their divergences cancel each other , since the sum of both yields the contribution of the standard continuous eigenvectors ( see eq . ( [ unity_2rep ] ) ) , so a divergence of an initial hilbert state backwards in time is prevented . the divergence of the gamov contributions is even more than compensated : if a hilbert state decays in the future , it will also decay in the past . to see this one can use the advanced decomposition  ( [ unity_ad ] ) to investigate the system s evolution to the past and one would find a completely symmetric situation with majorizing pole contributions decaying backwards in time . the situation is different if a gamov vector @xmath134 is taken as the initial state . now there is only one pole contribution @xmath68 decaying in the future and diverging in the past . but gamov vectors are only `` existing '' in the spectral representation of hilbert states and they themselves do not refer to elements of physical reality , similiar to all other `` generalized eigenstates '' appearing in spectral decompositions . as constituted by the axiomatic foundations of quantum mechanics a physically realizable state is represented by a hilbert vector ( actually a family of hilbert vectors normalized to unity and differing only in phase ) . gamov states are no hilbert vectors , so their properties are not expected to be physically meaningful .    concluding , both decompositions are equivalent , but the use of the retarded decomposition is more _ intuitive _ for @xmath130 . however , this should not be taken as the reason for a broken time symmetry . a more convincing argument in favour of a broken time symmetry is the splitting of _ test spaces_. as we have seen in  ( [ bg2 ] ) , analytic continuations of the test functions are needed to have the decompositions well - defined . in section  [ rhs ] we will get back to this point .      besides the advantages pointed out in section  [ advantages ] we now face another important advantage of the brussels school s formalism . perturbation theory can be applied even in cases where standard quantum mechanics fails . those cases are featured by the phenomenon of instability , which is strongly connected to the sensitivity of a system towards small perturbations . let the system s hamiltonian be of the form @xmath135 where @xmath136 is the perturbation operator , switched on by the perturbation parameter @xmath137 $ ] , and @xmath138 is the unperturbed part of the hamiltonian , whose spectral decomposition is known and reads @xmath139 to simplify the discussion let the spectrum of @xmath138 be non - degenerate , so @xmath140 for @xmath141 . as long as the perturbation strength @xmath142 is not explicetely needed , it is set to 1 and @xmath143 is written as @xmath22 . the spectral decomposition of @xmath22 is yet unknown and reads in general @xmath144 where the sum is taken over the discrete spectrum @xmath25 and the integration is performed over the continuous spectrum @xmath26 of @xmath22 . the continuous spectrum @xmath26 is not affected by the perturbation ( provided the perturbation operator @xmath136 is @xmath138-compact , which is assumed here ) and there is an identity mapping of perturbed and unperturbed continuous eigenvalues , so in both cases the integration is performed over @xmath26 . the discrete case is different . switching to _ brioullin - wigner _ perturbation theory let us define pairs of orthogonal projection operators @xmath145 furthermore let the perturbed eigenvector @xmath146 of @xmath22 be not normalized to 1 but rather let @xmath147 , such that @xmath148 . applying these projections to the perturbed eigenvalue equation @xmath149 one obtains the _ brioullin - wigner formulas _ for the shifted eigenstate @xmath146 and the shifted energy @xmath46 : @xmath150 iteration leads to the corresponding perturbation series @xmath151 take a closer look at the above formulas  ( [ briou1 ] ) and  ( [ briou2 ] ) . they are only valid if the expression @xmath152 is well - defined . here _ resonance _ appears as a major obstacle . there are two different types of resonance making a perturbation expansion of eigenstates and eigenvalues impossible :    1 .   * discrete resonance * : the eigenvalue @xmath153 may be shifted in the course of the perturbation to a neighboured eigenvalue of @xmath138 ( see fig .  [ fig : shift]a ) . expression  ( [ proj1 ] ) is no longer defined and the corresponding perturbation series diverges on account of the summation term . * continuous resonance * : if the eigenvalue @xmath153 lies next to the continuous spectrum of @xmath138 , it may be shifted there ( see fig .  [ fig : shift]b ) . now the integral term in  ( [ proj1 ] ) causes a divergence of the perturbation series .    in both cases the perturbation expansion is limited to a finite convergence radius of the perturbation parameter @xmath142 . outside its range the shifted eigenvalue disappears from the spectrum of @xmath154 and the mapping from the unperturbed to the perturbed spectrum is no longer one - to - one ( see fig .  [ fig : vanish ] ) . @xmath11{images / shift.eps}\ ] ]    now regard the case of discrete unperturbed eigenvalues _ embedded _ in the continuous spectrum . here the convergence radius of @xmath142 reduces to zero , i.e. the perturbation series are no longer analytic in @xmath142 . an arbitralily small perturbation induces continuous resonance . as a consequence embedded states get unstable and start to decay . decaying states can not be part of the spectral decomposition of @xmath154 , they disappear from there and can not be reconstructed for @xmath155 , i.e. @xmath156 . hence perturbation theory is not aplicable . @xmath157{images / vanish.eps}\ ] ]    there is a way out : one has to include unstable states , i.e. gamov vectors , into the spectral decomposition of @xmath154 . in terms of perturbation theory the gamov vectors can be constructed by modified brioullin - wigner equations : @xmath158 where @xmath159 are the singularities of @xmath43 on the lower half of the complex plane and @xmath160 are those of @xmath44 on the upper half . there is a one - to - one mapping between the disappearing discrete eigenvalues @xmath161 and the complex eigenvalues @xmath68 . iteration leads to the corresponding perurbation series @xmath162 the generalized continuous eigenstates are @xmath163 with the complex distribution @xmath164 defined in  ( [ compl_dist ] ) . the standard continous eigenvectors can be constructed using the _ lippmann - schwinger equation _ @xmath165 or its iteration , the _ born series _ : @xmath166 take a closer look at the construction of gamov states in eq . ( [ briou3 ] ) . the divergences due to continuous resonance are prevented by introducing a small complex shift in the denominator and analytic continuation to the second riemann sheet . as a consequence , the energy eigenvalues are shifted to respectively the lower and upper complex plane , inducing exponential decay and breit - wigner energy distribution of the corresponding gamov states . the expressions are still analytic in @xmath142 , so perturbation theory can be applied . this is a major advantage of the brussels school s formalism . the friedrichs model ( k. friedrichs 1948 , see  @xcite ) is the simplest model involving instability . here one single discrete state is coupled to a continuum of states by an interaction hamiltonian @xmath136 . a physical analogon to this model is e.g. the _ wigner - weisskopf model _ ( wigner and weisskopf 1930 , see  @xcite ) of a twolevel atom coupled to the vacuum field . another example is the _ auger effect _ of autoionisation ( see e.g.  @xcite ) . the total hamiltonian is of the form @xmath167 using the _ second resolvent identity _ @xmath168 with the resolvent @xmath169 of the free hamiltonian @xmath138 and the resolvent @xmath170 of the total hamiltonian . the resolvent @xmath37 can be obtained by applying it to the free states @xmath171 and @xmath172 and doing some rearrangement : @xmath173\nonumber \\      & & \left[\langle1|+\int_0^\infty d\omega\,\frac{w^*(\omega ) }      { z-\omega}\langle\omega|\right]\nonumber \\      & & + \int_0^\infty d\omega\frac1{z-\omega}|\omega\rangle\langle\omega|\quad.\end{aligned}\ ] ] the complex function @xmath174 has analytic extensions @xmath175 from above and below to the real axis . the zeros of @xmath176 on the second riemann sheet define the complex poles of the resolvent : @xmath177 the pole @xmath178 can be obtained by numerically calculating the complex zero of @xmath179 . in first order one gets @xmath180 the gamov vectors can be obtained using eqns . ( [ g_eigenvector1a ] ) and  ( [ g_eigenvector1b ] ) : @xmath181_+}|\omega\rangle\right]\\      \langle f_1 ^ -|&=&\frac1{c_1}\lim_{e\rightarrow z_1}(e - z_1)\langle1|r_+(e)\\      & = & \frac1{\sqrt{\eta_+{}'(z_1)}}\left[\langle1|+\int_0^\infty d\omega\ ,      \frac{w^*(\omega)}{[z_1-\omega]_+}\langle\omega|\right]\quad,\end{aligned}\ ] ] where the complex distribution is defined by @xmath182_+ }      : = \lim_{e\rightarrow z_1}\int_0^\infty d\omega\ ,      \frac{\varphi(\omega)}{e-\omega+i\epsilon}\quad.\ ] ] the standard outgoing states can be constructed using eq . ( [ eigenvector2 ] ) : @xmath183\quad.\end{aligned}\ ] ] using eq . ( [ g_eigenvector2 ] ) the generalized continuous eigenvectors read @xmath184\quad,\\      \langle f^-(\omega)|&=&\langle u^+(\omega)|\\      & = & \langle\omega|+\frac{w(\omega)}{\eta_-(\omega)}\times\nonumber \\      & & \times\left[\langle1|+\int_0^\infty d\omega'\,\frac{w^*(\omega')}{\omega-\omega'-i\epsilon }      \langle\omega|\right]\quad,\end{aligned}\ ] ] where the complex distribution @xmath164 is defined by  ( [ compl_dist ] ) with the curve @xmath75 leading from @xmath185 to @xmath186 below @xmath178 .    in the pole approximation , i.e. neglecting the background contribution from the continuous eigenvectors , the initial state @xmath171 decays with a decay rate @xmath187 given by two times the imaginary part of @xmath178 , @xmath188 which is in agreement with fermi s golden rule . when dealing with infinite dimensions , the hilbert space appears to be too small in the following sense : position and momentum operator have no hilbert eigenvectors and dirac s bra - ket formalism is mathematically not justified . to correct this problem , the hilbert space must be enlarged to include _ improper states _ like @xmath53-functions and plane waves . the framework to enable these features is given by the theory of the _ gelfand triplet _ , also called the _ rigged hilbert space _ ( rhs ) . @xmath11{images / leaving.eps}\ ] ]    the concept is the following ( comp . e.g.  @xcite ) . take the hilbert space @xmath189 of squared integrable functions and find a suitable test space @xmath190 with the properties    1 . @xmath191 is a locally convex topological vector space , 2 . @xmath191 is complete with respect to its own topology , 3 . @xmath191 is dense in @xmath192 ,    such that the topological dual @xmath193 of @xmath191 , consisting of all linear continuous functionals on @xmath191 , called _ distributions _ , includes the desired improper states . the hilbert space itself is his own topological dual , i.e. @xmath194 . for any test space @xmath191 being smaller than @xmath39 , the gelfand triplet @xmath195 makes it possible to extend the bracket @xmath196 to all pairs @xmath197 , @xmath198 by @xmath199 on pairs of hilbert vectors the bracket conicides with the scalar product in @xmath192 . the `` kets '' and `` bras '' are now identified by linear functionals on the opposite dual space , i.e. @xmath200 where @xmath201 or @xmath193 . an operator @xmath34 can be extended to the distribution space @xmath193 using @xmath202 hence @xmath34 can only be extended to @xmath193 if its adjoint @xmath203 is defined on and does not lead out of the test space @xmath191 : @xmath204    the test space of standard quantum mechanics is the schwartz space @xmath205 of rapidly decreasing functions , defined by @xmath206 the schwartz space is a physically convenient test space , since here position and momentum operator are everywhere defined , they fulfill  ( [ extend_op ] ) , they are self - adjoint and have their eigenvectors among the _ tempered distributions _ in @xmath207 . furthermore , they keep their real spectrum , which makes the rhs @xmath208 be called a _ tight rigged _ hilbert space . another popular test space is the space @xmath209 of @xmath210-functions with compact support . but @xmath209 is not tight rigged , since the momentum operator obtains a complex spectrum with complex plane waves @xmath211 as eigenvectors . now we turn to the rhs of the brussels school . as already stated in section  [ deform ] the definition of generalized eigenvectors is only possible if the test functions can be analytically continued to whether the lower or the upper half of the complex plane . the retarded right eigenvectors @xmath134 and @xmath212 involve complex distributions continuating the test functions @xmath213 to the lower half of the complex plane , denoted by @xmath214 . there is a class of functions fulfilling this condition : the hardy class @xmath215 from below ( for detailed information on hardy classes see  @xcite ) . if @xmath216 is in @xmath215 then its complex conjugate @xmath217 is in the hardy class @xmath218 from above , whose members can be analytically continuated to the upper half of the complex plane , denoted by @xmath219 . in detail : the function @xmath220 is in @xmath221 if and only if it has an analytic continuation to @xmath222 , such that there is a @xmath223 and @xmath224 thus , restricted to the real line , hardy functions are @xmath225-functions . an equivalent definition can be given using the _ paley - wiener theorem _ : the hardy spaces @xmath221 are formed by the inverse fourier transforms of @xmath225-functions with support on respectively the positive and negative semiaxis @xmath226 : @xmath227 for physical convenience it is necessary , that the test functions are also in the schwartz class @xmath228 . furthermore , their energy representation can be restricted to the positive semiaxis , since energy is bounded from below by zero . altogether , the test spaces of the brussels school are defined by @xmath229 the topological duals @xmath230 then include the complex distributions used in the generalized eigenvectors . the formalism of the brussels school can be applied in the rhs @xmath231 since there are _ two _ rhs s , one has to investigate how to deal with them . the retarded unity decomposition  ( [ unity_ret ] ) contains the dyadic products @xmath232 . applying test vectors @xmath233 and @xmath234 to both sides one obtains @xmath235 the above expression only makes sense , if @xmath127 and @xmath236 are both in @xmath237 , thus @xmath238 must be in @xmath239 . hence the instructions for computing brackets while using the retarded unity decomposition are : on the left hand you have to use @xmath239-vectors and on the right hand @xmath237-vectors . however , since both test spaces are dense in @xmath39 , you would not have to care too much , if there was not a significant subtlety : time symmetry is broken . to see this we investigate the time evolution of a test vector @xmath240 : @xmath241 since @xmath220 is in @xmath215 , there is an analytic continuation @xmath242 to @xmath214 and so there is one for its propagated counterpart , given by @xmath243 obviously , condition  ( [ hardy_def ] ) is violated for @xmath244 and @xmath133 , thus time evolution leads out of the test space for @xmath133 . as a consequence test functions in @xmath237 can only be propagated to the future . analogically , test functions in @xmath239 can only be propagated to the past . using condition  ( [ extend_op ] ) one finds in addition that the extension of the unitary propagator @xmath17 to the distribution spaces @xmath245 and @xmath246 is only possible for respectively @xmath130 and @xmath133 . so retarded gamov vectors @xmath134 in @xmath245 `` propagate '' ( as part of the spectral decomposition of hilbert vectors ) only to the future and advanced gamov vectors @xmath247 in @xmath246 propagate only to the past . the unitary propagator group @xmath248 hence splits into a semigroup @xmath249 propagating to the future ( and applying to vectors in @xmath237 and @xmath245 ) and a second one @xmath250 propagating to the past ( and applying to vectors in @xmath239 and @xmath246 ) . this feature is called by the brussels school the _ intrinsic irreversibility _ of quantum mechanics and is identified with a microscopic arrow of time . is the special choice of the two rhs s _ necessary _ for the correct mathematical implementation of gamov states ? if so , one could regard the splitting of the unitary time evolution group into two semigroups as enforced by the mathematical concept . this point of view is taken by _ a. bohm _ and _ n.l . harshman _ in  @xcite :    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ `` philosophizing alone would not be enough to take a semigroup instead . to arrive at the semigroup , we start from the empirically desirable properties of gamov resonance states and let mathematics determine the path . '' _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    the above mentioned properties are given by @xmath251 with @xmath252 . starting from these conditions , following the reasoning of section  [ intrinsic ] , one is led to the two rhs s with broken time symmetry . however , the argumentation is not valid if a rigged hilbert space can be found , where gamov vectors are included and time symmetry is preserved . let us start with the space @xmath209 of @xmath210-functions with compact support . as already mentioned , this space is a perfect test space fulfilling the requirements given in section  [ rhs2 ] . furthermore we have @xmath253 . now regard the space @xmath254 of inversely fourier transformed @xmath209-functions : @xmath255 since the inverse fourier transform is a one - to - one mapping on the schwartz space @xmath205 , we have @xmath256 let @xmath257 be a @xmath209-function . its inverse fourier transform @xmath258 reads @xmath259 where @xmath260 is the compact support of @xmath257 . the analytic continuation of @xmath220 to the complex plane exists and is an entire function , given by @xmath261 for all @xmath262 ( see  @xcite , theorem 7.22 ) . since the inverse fourier transform is a linear continuous one - to - one mapping , the transformed topology @xmath263 of @xmath209 , @xmath264 transfers all topological properties from @xmath209 to @xmath254 . thus , @xmath265 is a topological vector space fulfilling all requirements of a test space in @xmath266 . since @xmath254 contains functions being analytic in the whole complex plane , the topological dual @xmath267 contains the retarded gamov vectors defined by  ( [ req1 ] ) and  ( [ req2 ] ) . furthermore , it contains also the _ advanced _ gamov vectors @xmath247 with eigenvalues @xmath70 in the upper half of the complex plane , as well as any other generalized eigenvector used in the extended spectral decompositions  ( [ unity_ret ] ) and  ( [ unity_ad ] ) . because of relation  ( [ z_schwartz ] ) all @xmath254-functions are @xmath205-functions , so @xmath254 is already a physical convenient test space and standard quantum mechanics is applicable , too . since @xmath268 the space @xmath254 does neither include any of the test spaces @xmath269 nor is it part of it , i.e. @xmath270 but of course @xmath271 . so @xmath254 is a different test space , whose topological duals include any of the complex distributions needed for a comfortable treatment of decaying systems . however , there is no intrinsic irreversibility . let the @xmath209-function @xmath257 have a compact support bounded by the finite interval @xmath272 $ ] . its inverse fourier transform @xmath220 can be propagated in time yielding @xmath273 , whose fourier transform , @xmath274 has a compact support bounded by the finite interval @xmath275 $ ] , thus it remains compact for all @xmath276 and @xmath277 stays in @xmath209 , so @xmath278 stays in @xmath254 . hence the propagator @xmath17 does not lead out of the test space for any @xmath276 . consequently , the unitary propagator group @xmath279 can be extended to the distribution space @xmath267 without any splitting . the support of @xmath280-functions can be restricted to the positive semiaxis , due to the positivity of energy , and the rhs @xmath281 is thus suitable for the presented concepts . any desired unity decomposition including gamov vectors , retarded or advanced , can equivalently be used . as far as hilbert states are concerned , no unphysical behaviour is predicted and any calculated physical quantity is numerically identical to the one calculated by standard methods .      apart from pure mathematics there is a physical reasoning presented by members of the brussels school who try to derive the special choice of the two hardy classes from the assumption of _ causality _ ( see e.g. schulte and twarock in @xcite or antoine , bohm and harshman in @xcite , from where the citations below are taken ) . this is done the following way . starting point is a special assumption of _ causality _ which is called the `` preparation@xmath282registration arrow of time '' :    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ `` time translation of the registration apparatus relative to the preparation apparatus makes sense only by an amount of @xmath283 . '' _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    this is translated into the language of quantum mechanics to a so - called `` quantum mechanical arrow of time '' :    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ `` a state @xmath284 must be prepared before an observable @xmath285 can be measured in that state , i.e. @xmath286 must be prepared during @xmath287 . '' _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    interpreting measurement formally as a scattering experiment , the prepared state @xmath288 is identified with an incoming state and the measured state @xmath289 with an outgoing state ( see fig .  [ fig : measurement ] ) . survival probabilities of the initial ( i.e. prepared ) state are thus calculated by transition probabilities of the incoming and outgoing state : @xmath290 using the energy decomposition this implies @xmath291 using the paley - wiener theorem  ( [ paley - wiener ] ) one is led to the desired conclusion @xmath292 hence , founding on the given formulation of causality the only correct choice of a rigged hilbert space is the choice of the two rhs s of the brussels school . in other words : standard quantum mechanics has not yet found a microscopic arrow of time because it has been incomplete in its axiomatic foundation . @xmath293{images / measurement.eps}}\ ] ]      the above argumentation is correct , but it can not be taken as a _ derivation _ of the arrow of time . as is obvious even by the choice of words ( `` before '' , `` after '' , etc . ) , the given `` causality '' assumption is already based on the arrow of time , hence it is logically impossible to derive the latter from the former . causality itself does not require an arrow of time , since neither newton s law nor schrdinger s law contradict causality , although both theories are perfectly time - symmetric . yet , the `` preparation@xmath282registration arrow of time '' refers explicitely to _ open subsystems _ , namely the preparation apparatus and the measuring apparatus , which are coupled together to form a composite system . though one can not generally infer from statements about the time evolution of open subsystems to statements about the evoution of closed composite systems , especially not when the subsystems are involved in a measurement process . hence postulating a `` preparation@xmath282registration arrow of time '' for _ all _ systems equals introducing a microscopic arrow of time _ by hand _ into quantum mechanics . once the arrow of time has been put into the mathematics , it is not surprising that it will reappear eventually . the situation can be visualized by fig . [ fig : circular ] . @xmath294{images / circular.eps}\ ] ]    thus , the special choice of the two rhs s must be regarded as being _ equivalent to a postulation of a global microscopic arrow of time_. the split rhs s are a mathematical mirror image of the empirically motivated splitting between past and future . as already stated in the introduction , it is possible to explain an arrow of time due to time - asymmetric evolution of _ macroscopic _ observables and due to _ local _ observations on subsystems . but respecting a _ global microscopic arrow of time _ , there is a major difference between standard and rigged hilbert space quantum mechanics . in standard quantum mechanics there is no way to implement a global microscopic arrow of time except perhaps by applying brute force and splitting the unitary evolution group axiomatically into two semigroups . in rigged hilbert space quantum mechanics one has the choice of either time - symmetric or time - asymmetric microscopic evolution , due to the choice of the test spaces . hence , in the rhs formulation it is possible to implement a global microscopic arrow of time in contrast to standard quantum mechanics . so one could argue that nature _ obviously _ has taken the time - asymmetric choice of two test spaces instead of one . but what does `` obviously '' mean ? one might mention our everyday experiences or global time - asymmetric boundary conditions like the big bang or the second law of thermodynamics . though these experiences and boundary conditions are all related to subsystems or far from being microscopic . the `` universe '' is a construction from observations on subsystems on a _ very _ macroscopic scale . altogether , these time - asymmetric boundary conditions all base on _ local macroscopic observations _ and hence , as already stated , they do not give rise to a real problem with quantum mechanics , which is a set of microscopic laws . so the choice of global microscopic time - symmetry or time - asymmetry remains merely axiomatic and can not be derived from other physical postulates . the quantum mechanical formalism , developped by the brussels school , offers a comfortable treatment of decaying systems , where standard quantum mechanics requires a lot more effort . the generalized spectral decompositions involve gamov states and can be interpreted in a physically intuitive manner . they directly yield important physical quantities like decay rates and energy level shifts . perturbation expansions become possible in cases , where standard quantum mechanics fails . in contrast to standard quantum mechanics , where the implementation of a microscopic arrow of time is not possible , the rigged hilbert space formalism offers the possibility to incorporate such an arrow by the special choice of two different gelfand triplets . due to this choice , the unitary time evolution group splits into two semigroups , each one for respectively the future and the past . however , it is shown that the use of gamov states does not require this splitting . as a counter example , one single gelfand triplet is constructed , where the formalism of the brussels school can be applied and where time symmetry is preserved . furthermore , it is pointed out that the physical arguments in favour of the special choice of two gelfand triplets are not independent from the arrow of time and thus can not be used for a derivation of the latter . hence even within the framework of the rigged hilbert space formalism , the global microscopic arrow of time remains a postulate . i would like to thank f. petruccione , h .- breuer , m. bordemann , and m. wilkens for stimulating discussions and editorial advice and the referee for valuable comments . i also thank all other people talking with me about the topic of this paper and contributing in many ways to its completion .
Abstract:

Answer: despite the fact that the fundamental physical laws are symmetric in time , most observed processes do not show this symmetry . especially the phenomenon of decay seems to involve a kind of irreversibility that makes the definition of a microscopic arrow of time possible . such an intrinsic irreversibility is incorporated within the rigged hilbert space quantum mechanics of the brussels school , contrasting to the statements of standard quantum mechanics . as shown in this paper , the formalism bears significant advantages in the description of decaying systems , however the breaking of time symmetry can be avoided .    2


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: in the past decades , raman spectroscopy@xcite techniques were successfully applied to carbon compounds , such as graphite ( see ref .   and references therein ) and carbon nanotubes.@xcite upon the discovery of graphene,@xcite raman spectroscopy has proven to be a powerful tool to identify the number of layers , structure , doping , disorder , and to characterize the phonons and electron - phonon coupling.@xcite so far , most of the attention was focused on the position and width of the raman peaks .    here we present a detailed calculation of the _ intensities _ of the multiphonon raman peaks in graphene . raman scattering involves an electron - hole pair as an intermediate state ; we show that the multiphonon raman peaks are strongly sensitive to the dynamics of this electron - hole pair . thus , raman scattering can be used as a tool to probe this dynamics . writing the low - energy hamiltonian of the interaction of electrons with the crystal vibrations and the electromagnetic field from pure symmetry considerations , we describe the system in terms of just a few independent coupling constants , considered to be parameters of the theory . the electron scattering rate is introduced phenomenologically as another parameter . the results of the present calculation are used to extract information about these parameters from the raman peak intensities , measured experimentally . as shown below , the raman intensities strongly depend on the electron scattering rate ; moreover , the electron - phonon and electron - electron contributions to this rate can be separated . this is especially important as there are very few techniques giving experimental access to electron scattering rates , which , in turn , determine the transport properties of graphene samples . besides , the quasiclassical character of the process imposes a severe restriction on the electron and hole trajectories which can contribute to the two - phonon raman scattering : upon the phonon emission the electron and the hole must be scattered backwards . this restriction results in a significant polarization memory : it is almost three times more probable for the scattered photon to have the same polarization as the incidend photon than to have the orthogonal polarization . also , the raman intensities depend on electron - phonon coupling constants ; to reproduce the experimental results , one has to take into account renormalization of these coupling constants by electron - electron interaction . this renormalization is missed by local or semi - local approximations to the density - functional theory , typically used for the _ ab initio _ calculation of the coupling constants . since graphene is a non - polar crystal , raman scattering involves electronic excitations as intermediate states : the electromagnetic field of the incident laser beam interacts primarily with the electronic subsystem , and emission of phonons occurs due to electron - phonon interaction . the matrix element of the process can be schematically represented as    @xmath0    here @xmath1 is the initial state of the process ( the incident photon with a given frequency and polarization , and no excitations in the crystal ) , @xmath2 is the final state ( the emitted photon and @xmath3  phonons left in the crystal ) , while  @xmath4 , @xmath5 , label the intermediate states where no photons are present , but an electron - hole pair is created in the crystal and @xmath6  phonons have been emitted . @xmath7  and  @xmath8 , @xmath5 are the energies of these states , and @xmath9 is the inverse lifetime of the electron ( hole ) due to collisions . @xmath10  and  @xmath11 stand for the terms in the system hamiltonian describing interaction of electrons with the electromagnetic field and with phonons , respectively . in the calculations we do not include the phonon broadening , assuming the phonon states to have zero width . first , this approximation is consistent with the available experimental information : the phonon width is about @xmath12 at most,@xcite while the electronic broadening is at least an order of magnitude higher ( see the discussion below , sec .  [ sec:2intensities ] ) . second , this approximation is irrelevant provided that we calculate the integrated intensities of the raman peaks , since they are determined by the total spectral weight of the phonon state which does not depend on the phonon broadening . ( color online ) . schematic representation of the role of electron dispersion ( dirac cones , shown by solid lines ) in the one - phonon ( a , b ) and two - phonon raman scattering ( c , d ) . vertical solid arrows represent interband electronic transitions accompanied by photon absorption or emission ( photon wave vector is neglected ) , dashed arrows represent phonon emission , the horisontal dotted arrow represents the impurity scattering.,width=264 ]    the photon wave vector is negligible , so momentum conservation requires that the sum of the wave vectors of emitted phonons must vanish ( provided that the impurity scattering is neglected ) . for the same reason raman scattering on one intervalley phonon must be impurity - assisted [ process  ( b ) in fig . [ fig : res ] , giving rise to the so - called @xmath13  raman peak ] . @xmath13  peak is absent in the experimental raman spectrum of graphene,@xcite showing that impurity scattering is indeed negligible in these samples . looking at the intermediate electronic states involved in the raman scattering ( fig . [ fig : res ] ) , we notice that for one - phonon scattering [ processes ( a ) ,  ( b ) ] at least one intermediate state must be virtual , since energy and momentum conservation can not be satisfied simultaneously in all processes . thus , at least one of the factors in the denominator of eq . ( [ ramanmatrixelement= ] ) must be of the order of the phonon frequency  @xmath14 [ for the impurity assisted scattering one of the electron - phonon matrix elements in the numerator of eq . ( [ ramanmatrixelement= ] ) should be replaced by the electron - impurity matrix element ] . for the two - phonon scattering [ process  ( c ) ] all intermediate states can be real , so that all energy mismatches in the denominator of eq . ( [ ramanmatrixelement= ] ) can be nullified simultaneously and the result is determined by the electron scattering rate  @xmath9 . we emphasize the qualitative difference between the _ fully resonant _ process  ( c ) and the double - resonant@xcite process  ( b ) , where one intermediate state is still virtual . we also note the analogous difference between the two - phonon processes ( c )  and  ( d ) in fig . [ fig : res ] : only the process  ( c ) is fully resonant , the other one involves an energy mismatch of @xmath15 . as a result , its amplitude will be smaller by a factor @xmath16 . obviously , these arguments can be extended to all multi - phonon processes with odd and even number of phonons involved : in order to annihilate radiatively , the electron and the hole must have opposite momenta ; if the total number of emitted phonons is odd , the electron and the hole must emit a different number of phonons , which is incompatible with energy conservation in all processes . our main focus will thus be on even - phonon processes , as their intensities are determined by the electronic scattering . the full resonance picture presented above assumes the mirror symmetry between the electron and the hole spectra . the electron - hole asymmetry can be included as a correction to the dirac spectrum : the electron and the hole energies can be written as @xmath17 , where @xmath18  is the momentum counted from the dirac point , @xmath19  is the dirac velocity , and @xmath20  is the asymmetry parameter . the energy scale  @xmath21 , quantifying the role of the asymmetry in the raman scattering is defined as @xmath22 , where @xmath23  and  @xmath24 are the frequencies of the incident and the scattered photon ( the details are given in secs . [ sec : warping ] ,  [ sec : intdev ] ) . namely , the arguments of the previous paragraph hold if @xmath25 . in the opposite case , it is @xmath21 that determines the smallest value of the denominators in eq . ( [ ramanmatrixelement= ] ) . we will always assume that both @xmath26 .    in the real space , the typical size of the region of space probed by the electron - hole pair in the fully resonant two - phonon raman scattering , is @xmath27 . for the doubly resonant defect - induced one - phonon scattering , the inverse energy mismatch  @xmath28 determines the time duration of the process by virtue of the uncertainty principle , so the length scale of the process in the real space is @xmath29 . although this length scale is much shorter than @xmath30 , it is still much greater than the lattice constant or the electron wavelength @xmath31 . the fully resonant raman scattering , where the energy is conserved in each of the elementary scattering processes , admits a simple quasiclassical description , described qualitatively in this subsection , and justified rigorously in secs .  [ sec:2raman],[sec:4raman ] . let us denote by  @xmath32 the energy of the electron and the hole in the photoexcited pair . initially , it is given by the half of the excitation frequency  @xmath23 , @xmath33 . after the emission of @xmath3  phonons it is decreased by  @xmath34 ; assuming @xmath35 , we neglect this decrease in the qualitative considerations . thus , during all the time taken by the raman scattering , electron and hole can be viewed as wave packets of the size @xmath36 , propagating across the crystal along classical trajectories . the electron and the hole are created in the same region of space of the size @xmath36 around some point  @xmath37 at the moment of the arrival of the excitation photon . at this initial moment they have opposite momenta @xmath38 , and opposite velocities @xmath39 ( if the electron - hole asymmetry is taken into account , the two velocities will have slightly different magnitude ) , so they move along the straight lines , their positions being @xmath40 , @xmath41 . after a typical time @xmath42 they undergo some scattering processes ( e.  g. , phonon emission ) , where their momenta and ( generally speaking ) energies are changed . each such elementary scattering process occurs during a short time @xmath43 . thus , the trajectories of the electron and the hole after their creation are represented by broken lines , with the typical segment length @xmath44 ( the electron mean free path ) . the crucial point is that in order to recombine radiatively and contribute to raman signal , the electron and the hole should meet again within a spatial region of the size @xmath36 , and have opposite momenta . the latter condition automatically implies that the number of the phonons emitted by the electron and the hole is the same . these considerations are illustrated by fig . [ fig : trajectories ] .    in the presence of a significant electron - hole asymmetry , @xmath45 , the described picture is modified . namely , one of the segments of either the electron or the hole trajectory has the length @xmath46 instead of @xmath30 , the corresponding time travelling being restricted by the phase mismatch rather than by collisions . ( color online ) . ( a ) an example of a quasiclassical electron - hole trajectory contributing to the four - phonon raman scattering . ( b , c )  trajectories with emission of two phonons , not contributing  ( b ) and contributing  ( c ) to the two - phonon raman scattering . in all pictures the lightning represents the incident photon which creates the pair . the solid lines denote the free propagation of the electron and the hole . the flash represents the radiative recombination of the electron - hole pair . the dashed lines denote the emitted phonons.,title="fig:",width=302 ]   ( color online ) . ( a ) an example of a quasiclassical electron - hole trajectory contributing to the four - phonon raman scattering . ( b , c )  trajectories with emission of two phonons , not contributing  ( b ) and contributing  ( c ) to the two - phonon raman scattering . in all pictures the lightning represents the incident photon which creates the pair . the solid lines denote the free propagation of the electron and the hole . the flash represents the radiative recombination of the electron - hole pair . the dashed lines denote the emitted phonons.,title="fig:",width=302 ]   ( color online ) . ( a ) an example of a quasiclassical electron - hole trajectory contributing to the four - phonon raman scattering . ( b , c )  trajectories with emission of two phonons , not contributing  ( b ) and contributing  ( c ) to the two - phonon raman scattering . in all pictures the lightning represents the incident photon which creates the pair . the solid lines denote the free propagation of the electron and the hole . the flash represents the radiative recombination of the electron - hole pair . the dashed lines denote the emitted phonons.,title="fig:",width=302 ] for the single - phonon raman peaks the commonly accepted notations are `` @xmath47 '' for the peak at @xmath48  corresponding to emission of an optical phonon with zero wave vector , and `` @xmath13 '' for the defect - induced peak at @xmath49 corresponding to emission of an optical phonon with the wave vector near @xmath50  or  @xmath51 points of the brillouin zone.@xcite sometimes one also distinguishes the so - called @xmath52-peak at @xmath53 . this peak is also defect - induced , and corresponds to emission of an optical phonon with a small wave vector @xmath54 . as mentioned above , in the present work we study only the clean graphene , hence @xmath13  and  @xmath52 peaks are of no interest to us . unfortunately , there is no single commonly accepted system for labelling of the multiphonon raman peaks . the strong peak at @xmath55 , corresponding to emission of two phonons with the opposite wave vectors near the @xmath50  and  @xmath51 points , was historically called  @xmath56 ( as it is not defect - induced ) ; sometimes it is denoted by @xmath57 or by  @xmath58 to stress that it is the second overtone of the @xmath13  peak . the peak at @xmath59 corresponding to emission of two phonons with two opposite wave vectors near the @xmath60  point is sometimes called @xmath61 , @xmath62 or @xmath63 . the latter notation reflects the fact that the frequency of this peak is not exactly the double of that of the @xmath47  peak , but rather the double of the defect - induced @xmath52  peak .    in the following we use the notation @xmath64 to denote the peak corresponding to emission of @xmath3  phonons with wave vectors within @xmath65 from the @xmath60  point and of @xmath66  phonons with wave vectors within @xmath65 from the @xmath50  or  @xmath51 points . for multiphonon peaks this nomenclature is unambiguous . thus , the peaks at @xmath55 and at @xmath59 will be called @xmath67 and @xmath68 , respectively .      in the clean graphene the only one - phonon raman process , allowed by the momentum conservation corresponds to the emission of the @xmath69  optical phonon with zero wave vector and frequency @xmath48 . for this process the situation turns out to be drastically different from that described by the qualitative considerations of sec . [ sec : fullresonance ] . as shown in sec . [ sec:1raman ] , if one approximates the electron spectrum by the dirac cones , the numerator of eq . ( [ ramanmatrixelement= ] ) vanishes due to high symmetry of the low - energy electronic dirac hamiltonian ( as compared to the microscopic symmetry of the crystal ) . thus , the main contribution to the raman amplitude comes from the regions of the electronic brillouin zone far from the dirac points . as a consequence , the typical energy mismatch in the denominator of eq . ( [ ramanmatrixelement= ] ) is of the order of the whole electronic bandwidth . thus , the raman process , responsible for the @xmath48  peak , is completely off - resonant and the picture shown in fig . [ fig : res ] for the process  ( a ) is wrong . as a result , the intensity of the peak is expected to be insensitive to most external parameters : polarization , electron concentration , degree of disorder , etc . to characterize this intensity , one has to introduce an additional parameter into the theory which has no simple relation to the parameters of the low - energy effective hailtonian . the resulting intensity of the peak is given by eq . ( [ ig= ] ) ; it is proportional to the fourth power of the excitation frequency which is the standard result for raman scattering when the difference between the frequencies of the incident and scattered photons is small . this dependence also agrees with the experimental results of ref .  . note that the results described above do not hold for the defect - induced peak at @xmath49 . for this peak the double resonance picture,@xcite shown in fig . [ fig : res ] , process  ( b ) , is fully adequate .      as the phonons are emitted by electrons with momentum @xmath70 , the largest possible phonon momentum is @xmath71 , corresponding to the electron and hole backscattering ( for the @xmath67  peak at 2700  @xmath72 we count the phonon momenta from the @xmath50  and  @xmath51 points ) . it would be natural to expect that any pair of phonons with opposite momenta @xmath73 and @xmath74 can be emitted , the only exception being @xmath75 which is prohibited by symmetry@xcite and the nearby ones which are suppressed due to the smallness of the matrix elements . these arguments would predict the width of the peak to be of the order of @xmath76 , where @xmath77  is the @xmath50  phonon group velocity ; besides , the shape of the peak would be strongly asymmetric : a sharp cutoff on the high - energy side at the frequency @xmath78 due to the resonance restriction , and a smooth drop - off towards zero at @xmath79 due to the matrix element suppression . the phonon dispersion can be deduced from the dependence of the frequency of the impurity - assisted one - phonon @xmath13  peak in graphite on the excitation energy  @xmath23 : @xmath80.@xcite thus , for @xmath81  ev these arguments give the width of the @xmath67  peak to be about @xmath82  @xmath72 . however , the experimentally observed width is only about 30  @xmath72 at @xmath83  ev , and its shape is quite symmetric.@xcite    the observed small width of the peak is explained by the quasiclassical picture , presented in sec .  [ sec : realspace ] . if upon the emission of phonons the electron and the hole are scattered by an arbitrary angle , as shown in fig . [ fig : trajectories](b ) , they will not be able to meet at the same spatial point in order to recombine radiatively and contribute to the two - phonon raman peak . only if the scattering is backwards , this event is possible , as illustrated by fig . [ fig : trajectories](c ) . this condition fixes the wave vectors of the emitted phonons to be @xmath84 . the small deviations of the scattering angle from  @xmath85 are restricted by the quantum diffraction , and the width of the two - phonon raman peaks , instead of being @xmath86 , is determined by a much smaller energy scale ( see the discussion below ) . the dominance of the electron and hole backscattering manifests itself in the polariazion memory of the raman signal . if the incident light is linearly polarized , the probability of excitation of the electron - hole pair with a given direction of momenta is proportional to @xmath87 , where @xmath88  is the angle between the electric field vector of the light and the momenta . thus , upon backscattering and radiative recombination , the probability to detect a photon of the same polarization as the original one is @xmath89 , and that of the orthogonal polarization is @xmath90 . averaging over  @xmath88 , we obtain the ratio of intensities for the detection of polarization parallel and perpendicular to that of the incident light to be @xmath91 . this ratio may be slightly decreased due to a finite aperture ( see the discussion in sec .  [ sec : ramanprobab ] ) . the calculation of the intensities of the two - phonon raman peaks is performed in sec . [ sec:2raman ] . the explicit expressions for the intensities of the @xmath67  and  @xmath68 peaks , obtained under the assumption of dirac spectrum for the electrons , are represented by eqs . ( [ i2kend= ] ) and  ( [ i2gend= ] ) . both are proportional to @xmath92 , where @xmath9  is the electron ( hole ) inelastic scattering rate . if the latter is smaller that the electron - hole asymmetry  @xmath21 , then , according to the arguments of sec .  [ sec : fullresonance ] , it is @xmath21 that restricts the energy denominators from below . formally , this results in the replacement  ( [ i2ktrig= ] ) in both eqs . ( [ i2kend= ] ) ,  ( [ i2gend= ] ) . numerically , @xmath93 ( see , e.  g. , ref .  ) , so the relative correction to eq . ( [ i2kend= ] ) for small @xmath21 can be estimated as @xmath94 ^ 2/2\sim-{10}^{-4}(\omega_{in}/2\gamma)^2 $ ] . the total electronic broadening @xmath9 was measured by time - resolved photoemission spectroscopy to be 20  mev in ref .   and @xmath95  mev in ref . ( all values taken for @xmath96  ev ) . a recent angle - resolved photoemission spectroscopy ( arpes ) measurement gives a significantly larger value for @xmath97 ( ref .  ) . thus , the case @xmath98 seems to be more relevant for the description of experiments , than the opposite one . the raman matrix element corresponding to emission of two phonons with given wave vectors @xmath99 and @xmath100 is given by eq . ( [ mq2ph= ] ) for @xmath98 . from this dependence one can deduce the lineshape of the two - phonon peaks @xmath101 : @xmath102^{3/2}},\ ] ] where @xmath103 , @xmath104 is the central frequency for each peak , and @xmath105  is the group velocity of the corresponding phonon . thus , the full width at half maximum ( fwhm ) of each peak is given by @xmath106 for @xmath107 the dependence of the raman matrix element on  @xmath108 is described by eq . ( [ warpedbackscattering= ] ) . the lineshape corresponds to _ two _ peaks separated by @xmath109 . experimentally , one sees just one @xmath67 peak with the fwhm about 30  @xmath72 at the excitation frequency @xmath110.@xcite this corresponds to an unrealistically large value of @xmath111 . most likely , this indicates that two - phonon peaks are broadened by other mechanisms , not taken into account in the present work . in particular , eq . ( [ di2dw= ] ) neglects ( i )  the broadening of the phonon states , and ( ii )  the anisotropy of the phonon dispersion ( trigonal warping of the phonon spectrum ) . a detailed study of these effects would require introduction of additional parameters into the theory , so we prefer to postpone such study for the future work . it is worth emphasizing again that the _ integrated _ intensity of the peaks , which is the main focus of the present study , does not depend on these details .    in view of the results of the present paper it is worth mentioning the experimental measurements of the intensity @xmath112 as a function of doping . while in ref . no significant dependence was observed , ref . , where higher doping levels were reached , shows quite a strong dependence of @xmath113 on doping . the intensity @xmath114 of the off - resonant single - phonon @xmath48 peak should not depend on doping ( although the phonon width does exhibit such a dependence , the total spectral weight of the phonon state , determining the integrated intensity of the peak , must be preserved ) . at the same time , the intensity  @xmath112 , if determined by the electron inelastic lifetime , should be sensitive to the concentration of carriers . indeed , in the intrinsic graphene at low temperatures the photoexcited carriers do not participate in electron - electron collisions , as the phase space volume is restricted.@xcite as the carriers are added to the system , the electron - electron collisions become possible , thus the total  @xmath115 increases , and the intensity  @xmath112 is decreased , in qualitative agreement with the observation of ref .  . the motivation to study the four - phonon raman process comes from the following picture for the fully resonant processes . the incident photon creates an electron and a hole  real quasiparticles which can participate in various scattering processes . if the electron emits a phonon with a momentum  @xmath99 , the hole emits a phonon with the momentum @xmath100 , and after that the electron and the hole recombine radiatively , the resulting photon will contribute to the two - phonon raman peak . if they do not recombine at this stage , but each of them emits one more phonon , and they recombine afterwards , the resulting photon will contribute to the four - phonon peak , _ etc . _ three - phonon processes , not being fully resonant , are not interesting in this context . besides phonon emission and radiative recombination , electron and hole are subject to other inelastic scattering processes , which can also be viewed as emission of some excitations of the system . in principle , raman spectrum should also contain the contrubution from these excitations , which are left in the system after the radiative recombination of the electron and the hole . the key point is that for real quasiparticles , the probability to undergo a scattering process  @xmath116 is determined by the ratio of corresponding scattering rate @xmath117 to the total scattering rate @xmath118 , not by the history . this probability determines the relative _ frequency - integrated _ intensity of the corresponding feature in the raman spectrum . thus , the ratio of integrated intensity @xmath119 of the raman peak corresponding to @xmath120 @xmath50  phonons to that for @xmath121 @xmath50  phonons ( @xmath122 ) must be proportional to @xmath123 , where @xmath124 is the rate of emission of each of the two @xmath50  phonons , and the square comes from the phonon emission by the electron and the hole . this conclusion depends weakly on the relation between @xmath115  and  @xmath21 , only through a logarithmic factor .    in the doped graphene , the most obvious competitor of the phonon emission is the electron - electron scattering : the optically excited electron can kick out another one from the fermi sea , i.  e. , to emit another electron - hole pair . thus , raman spectrum should contain contribution from electron - hole pairs ; however , their spectrum extends all the way to the energy of the photo - excited electron ( optical energy ) in a completely featureless way . thus , it can not be distinguished from the parasitic background which is always subtracted in the analysis of raman spectra , and can not be seen in the raman spectrum directly . however , assuming @xmath125 , where @xmath126 is the rate of emission of phonons from the vicinity of the @xmath60  point of the first brillouin zone , and @xmath127 is the electron - electron collision rate , one can extract the value of @xmath128 , relative to phonon emission rates from the experimental data . more precisely , in this way one obtains the rate of all inelastic scattering processes where the electron loses energy far exceeding the phonon energy . note that arguments leading to @xmath129 are not specific for graphene ; in fact , this is nothing but breit - wigner formula , applied once for the electron and once for the hole . multi - phonon raman scattering has been studied in wide - gap semiconductors both experimentally  @xcite ( up to ten phonons were seen in the raman spectra of cds ) , and theoretically  @xcite . in a wide - gap semiconductor an optically excited electron does not have a sufficient energy to excite another electron across the gap , so the electron - electron channel is absent . in addition , interaction with only one phonon mode is dominant , so the ratios of subsequent peaks are represented by a sequence of fixed numbers . the simple band structure ( one valley for cds in contrast to two valleys for graphene ) allowed a calculation of the whole sequence . a more complicated electronic band structure in graphene makes it problematic to calculate the whole sequence , so we restrict ourselves to the calculation of  @xmath130 for the most intense four - phonon peak . this calculation is performed in sec . [ sec:4raman ] . its result depends , besides the relation between @xmath115  and  @xmath21 , also on their relation to the energy scale @xmath131 , characterizing the phonon dispersion . in sec . [ sec:2intensities ] we have already discussed this energy scale ; for @xmath132 we have @xmath133 . the meaning of this energy scale is the difference between the energies of the electron and the hole after each of them has emitted two phonons with almost arbitrary momenta ( the only restriction is that the sum of all four phonon momenta must vanish ) . if @xmath134 , which seems to be the case ( see the discussion in the previous subsection ) then this difference can be neglected , and the intensity of the @xmath135  peak is given by eq . ( [ i4gammadispless= ] ) for @xmath98 [ which is likely to be the case relevant for most experiments , and which was reported in the short paper by the author ( ref .  ) ] , and by eq . ( [ i4dehdispless= ] ) for @xmath107 ; the polarization memory is lost in both these cases , @xmath136 . in the case @xmath137 the intensity  @xmath130 is given by eq . ( [ i4kdisp= ] ) , and a significant polarization memory is expected , up to @xmath138 .    a thorough experimental study of the intensity @xmath130 ( in particular , its dependence on doping ) is still lacking . our prediction for the case @xmath139 , which we believe to be the experimentally relevant one , is@xcite @xmath140      in the calculations of the raman intensities , described above , electron - phonon coupling constants entered as parameters of the theory , without any assumptions about their values , except for the relations fixed by the symmetry of the crystal . in particular , the two - phonon peak intensities @xmath112  and  @xmath141 are determined by two independent dimensionless coupling constants which we denote @xmath142  and  @xmath143 [ see eq . ( [ lambdaphonon= ] ) for the definition ] . a simple estimate of the coupling constants can be obtained from the tight - binding nearest - neighbor model of the graphene crystal . in this model the only parameter characterizing the electron spectrum is the nearest - neighbor electronic matrix element  @xmath144 , and the electron - phonon interaction is characterized by its change with the bond length , @xmath145 . in this model we obtain @xmath146 to be given by the inverse ratio of the corresponding phonon frequencies , about  1.2 ; the same result up to a few percent is obtained from the density - functional theory ( dft ) calculations of ref .  . at the same time , by comparing the experimentally measured intensities of the different two - phonon peaks , and using the result of our calculation performed in sec . [ sec:2raman ] , we can independently extract the ratio of the coupling constants . according to the data of ref .  , @xmath147 , which gives @xmath148 . to explain this discrepancy we first analyzed the effect of the electronic trigonal band warping , which affects @xmath112 and @xmath141 differently . the corresponding calculation is done in sec . [ sec : warping ] . for @xmath81 ev , we estimate the relative contribution of the warping term as @xmath149 for  @xmath112 and @xmath150 for  @xmath141 , which is far too little to account for the observed ratio @xmath151 . we are thus led to the conclusion that the observed ratio @xmath151 must be due to the difference of the coupling constants , not accounted for by the dft calculation . a similar conclusion about the insufficiency of the dft calculation of the electron - phonon coupling constants has been drawn in ref .  , where an attempt was made to explain the experimental data obtained by arpes@xcite using the results of the dft calculation . at the same time , we should note that the dimensionless coupling constant @xmath143 for the phonons near the @xmath60  point , as calculated by dft@xcite ( @xmath152 ) , agrees reasonably well with the measured one : the measurements of the linear in the wave vector  @xmath108 term in the phonon dispersion ( kohn anomaly due to electron - phonon interaction ) , @xmath153 , give @xmath154 ( see ref .  ) ; the measurements of the dependence of the phonon frequency @xmath14 on the electron fermi energy  @xmath155 : @xmath156 give @xmath157 ( ref .  ) @xmath158 ( ref .  ) . we show that the difference between the ratio @xmath148 extracted from the raman peak intensities , and @xmath159 obtained by the dft calculation,@xcite is due to the part of coulomb interaction between electrons , not picked up by the dft when local approximations are used for the exchange - correlation functional , such as the local density approximation ( lda ) or the generalized gradient approximation ( gga ) , namely , logarithmic renormalizations.@xcite coulomb interaction has been known to be a source of logarithmic renormalizations for dirac fermions.@xcite coulomb renormalizations in graphene subject to a magnetic field have been considered in ref .  , coulomb effect on static disorder has been studied in refs .  . essentially , the idea of the renormalization of the coupling constants is that the matrix element of the electron - phonon interaction should be taken not between the non - interacting electronic states , but between the states dressed by the coulomb interaction . if the typical electronic energy in the problem is  @xmath32 ( @xmath160 in the case of raman scattering ) , the renormalization is determined by the coulomb interaction at all length scales from the shortest ones ( lattice constant ) to the electron wavelength @xmath31 . it is this long - range part of the exchange and correlation that is missed by the local approximations in the dft calculation , which take into account correctly only the short - range correlations ( at the distances of the order of the lattice constant ) .    in sec . [ sec : coulombrg ] we calculate the renormalization of the dimensionless electron - phonon coupling constants ( a preliminary account of this work was given in the short paper@xcite ) , and show that the coupling constant  @xmath143 for the phonons near the @xmath60  point is not renormalized ( hence the agreement between the value of  @xmath143 calculated by the dft and measured in the experiments , as mentioned above ) , while the coupling constant  @xmath142 for the phonons near the @xmath161  point , which is responsible for the @xmath67  raman peak , is enhanced by the coulomb interaction . this enhancement depends on the electronic energy , as shown in fig . [ fig : rgflow ] . for the electronic energy @xmath162 this enhancement is in quantitative agreement with the measured ratio @xmath151 , provided that the screening of the coulomb interaction by the substrate is weak . the dependence of the enhancement on the electronic energy translates into the dependence of @xmath151 on the excitation frequency , which can be checked experimentally . similarly , as the coulomb interaction is screened by the substrate with a dielectric constant  @xmath163 ( its high - frequency value ) , the dependence of @xmath151 on  @xmath163 can also serve as an experimental check of the theory . we also show in sec . [ sec : phononrg ] that the electron - phonon coupling itself is a source of logarithmic renormalizations . however , due to the smallness of the coupling constants this effect is much weaker than the effect of the coulomb interaction .      in sec . [ sec : hamiltonian ] the low - energy hamiltonian of the interaction of electrons with the crystal vibrations and the electromagnetic field is written from pure symmetry considerations . in sec . [ sec : symmetry ] the symmetry of the graphene crystal is reviewed . in sec . [ sec : electrons ] the symmetry considerations are used to write the electronic part of the hamiltonian . [ sec : adddirac ] is dedicated to the symmetry analysis of the dirac part of the electron hamiltonian , whose symmetry is significantly higher than the symmetry of the crystal . [ sec : parphonons ] is dedicated to the symmetry analysis of the in - plane crystal vibrations . in secs . [ sec : eopt ] and  [ sec : eac ] we write the hamiltonian of interaction of electrons with the optical and acoustical vibrations , respectively . [ sec : perpphonons ] is dedicated to the symmetry analysis of the out - of - plane vibrations of the graphene crystal . in sec . [ sec : emfield ] the hamiltonian of the interaction of electrons with the electromagnetic field is written .    sec .  [ sec : ramangeneral ] describes the general scheme of the calculation of raman peak intensities using the standard perturbation theory . in sec . [ sec : propagators ] the green s functions are introduced , and in sec .  [ sec : ramanprobab ] the general expression for the raman scattering probability is derived . in sec . [ sec : inelastic ] we discuss the electron inelastic scattering , and calculate the electronic self - energy due to the electron - phonon coupling .    in sec . [ sec:1raman ] one - phonon raman scattering is discussed , and it is shown that the calculation of the one - phonon peak intensity can not be performed within the low - energy theory . in sec . [ sec:2raman ] the two - phonon raman peak intensities are calculated , first , under the assumption of the dirac electron spectrum ( sec . [ sec:2ramandirac ] ) , and then taking into account the trigonal band warping and electron - hole asymmetry ( sec .  [ sec : warping ] ) . in sec . [ sec:4raman ] the intensity @xmath130 of the most intense four - phonon peak is calculated . [ sec : coulombrg ] ,  [ sec : phononrg ] are dedicated to the renormalization of the electron - phonon coupling constants due to coulomb interaction and due to the electron - phonon interaction , respectively . since the typical energy of the incident photon ( about @xmath164  ev ) is much smaller than the @xmath85-electron bandwidth ( @xmath165  ev ) , one can expect the low - energy excitations to play the dominant role . in this section we employ standard symmetry analysis@xcite to fix the form of the low - energy hamiltonian . we prefer not to choose any specific basis and use algebraic properties . the carbon atoms of graphene form a honeycomb lattice with two atoms per unit cell , labeled @xmath166  and  @xmath167 ( fig . [ fig : lattice ] ) , the distance between nearest neighbors being @xmath168   . three out of four electrons of the outer shell of each carbon atom form strong @xmath169  bonds with its three nearest neighbors , and represent no interest to us . the remaining @xmath85  orbitals ( one per each carbon atom ) give rise to the half - filled @xmath85  band . the honeycomb lattice with two atoms ( @xmath166  and  @xmath167 ) per unit cell . tripled unit cells are shown by dashed hexagons.,width=302 ]    .irreducible representations of the groups @xmath170 and @xmath171 and their characters.[tab : c6vc3v ] [ cols="^,^,^,^,^,^,^",options="header " , ] polarization operator with arbitrary matrix vertices @xmath172 , @xmath173 , @xmath174 ( we denoted @xmath175 ) , can be conveniently calculated in the coordinate representation:@xcite @xmath176 if @xmath177  matrices are also present in the vertices , tracing of them is trivial : the matrices in the two vertices must coincide for the trace not to vanish , then their product is equal to the unit matrix . it is convenient to switch to the imaginary time @xmath178 , @xmath179 . then the green s function in the coordinate representation can be found by using the analogy with the 3d coulomb problem . namely , we introduce the third dimension @xmath180 , and notice that @xmath181 since the inverse of the 3d laplacian is the coulomb potential , we obtain @xmath182 using auxiliary relations ( here @xmath183 ) ,                    99 l.  i.  mandelshtam and g.  s.  landsberg , z.  phys . * 50 * , 169 ( 1928 ) ; c.  v.  raman and k.  s.  krishnan , nature * 121 * , 501 ( 1928 ) ; _ ibid . _ * 121 * , 619 ( 1928 ) . s. reich and c. thomsen , phil . a * 362 * , 2271 ( 2004 ) . j. kastner , t. pichler , h. kuzmany , s.  curran , w.  blau , d.  n.  weldon , m.  dlamesiere , s.  draper , and h.  zandbergen , chem . . lett . * 221 * , 53 ( 1994 ) . a. m. rao , e. richter , s.  bandow , b.  chase , p.  c.  eklund , k.  w.  williams , m.  menon , k.  r.  subbaswamy , a.  thess , r.  e.  smalley , g.  dresselhaus , and m.  s.  dresselhaus , science * 275 * , 187 ( 1997 ) . k. s. novoselov , a.  k.  geim , s.  v.  morozov , d.  jiang , y.  zhang , s.  v.  dubonos , i.  v.  grigorieva , and a.  a.  firsov , science * 306 * , 666 ( 2004 ) . a. c. ferrari , j.  c.  meyer , v.  scardaci , c.  casiraghi , m.  lazzeri , f.  mauri , s.  piscanec , d.  jiang , k.  s.  novoselov , s.  roth , and a.  k.  geim , phys . lett . * 97 * , 187401 ( 2006 ) . a. gupta , g. chen , p. joshi , s.  tadigadapa , and p.  c.  eklund , nano lett . * 6 * , 2667 ( 2006 ) . d.  graf , f.  molitor , k.  ensslin , c.  stampfer , a.  jungen , c.  hierold , and l.  wirtz , nano lett . * 7 * , 238 ( 2007 ) . j.  yan , y.  zhang , p.  kim , and a.  pinczuk , phys . rev . lett . * 98 * 166802 ( 2007 ) . s.  pisana , m.  lazzeri , c.  casiraghi , k.  s.  novoselov , a.  k.  geim , a.  c.  ferrari , and f.  mauri , nature mat . * 6 * , 198 ( 2007 ) . c. casiraghi , s. pisana , k.  s.  novoselov , a.  k. geim , and a.  c.  ferrari , appl . * 91 * , 233108 ( 2007 ) . a.  das , s.  pisana , s.  piscanec , b.  chakraborty , s.  k.  saha , u.  v.  waghmare , r.  yiang , h.  r.  krishnamurhthy , a.  k.  geim , a.  c.  ferrari and a.  k.  sood , arxiv:0709.1174 . c.  thomsen and s.  reich , phys . lett . * 85 * , 5214 ( 2000 ) . l. g. canado , a.  jorio , and m.  a.  pimenta , phys . rev . b * 76 * , 064304 ( 2007 ) . r. p. vidano , d.  b.  fischbach , l.  j.  willis and t.  m.  loehr , solid state commun . * 39 * , 341 ( 1981 ) . y.  wang , d.  c.  alsmeyer , and r.  l.  mccreery , chem . * 2 * , 557 ( 1990 ) . i.  pcsik , m.  hundhausen , m.  koos , and l.  ley , j.  non - cryst . solids * 227230 * , 1083 ( 1998 ) . m.  j.  matthews , m.  a.  pimenta , g.  dresselhaus , m.  s.  dresselhaus , and m.  endo , phys . b * 59 * , r6585 ( 1999 ) . r. saito , g.  dresselhaus , and m.  s.  dresselhaus , _ physical properties of carbon nanotubes _ ( imperial college press , london , 1998 ) . s. xu , j.  cao , c.  c.  miller , d.  a.  mantell , r.  j.  d.  miller , and y.  gao , phys . rev . lett . * 76 * , 483 ( 1996 ) . g.  moos , c.  gahl , r.  fasel , m.  wolf , and t.  hertel , phys . rev . * 87 * , 267402 ( 2001 ) . a.  bostwick , t.  ohta , t.  seyller , k.  horn , and e.  rotenberg , nature physics * 3 * , 36 ( 2007 ) . j. gonzlez , f.  guinea , and m.  a.  h.  vozmediano , phys . lett . * 77 * , 3589 ( 1996 ) . r.  c.  c.  leite , j.  f.  scott , and t.  c.  damen , phys . * 22 * , 780 ( 1969 ) . m.  v.  klein and s.  p.  s.  porto , phys . * 22 * , 782 ( 1969 ) . r.  m.  martin and c.  m.  varma , phys . * 26 * , 1241 ( 1971 ) . r.  zeyher , solid state commun . * 16 * , 49 ( 1975 ) . d.  m.  basko , phys . b * 76 * , 081405(r ) ( 2007 ) . s.  piscanec , m.  lazzeri , f.  mauri , a.  c.  ferrari , and j.  robertson , phys . * 93 * , 185503 ( 2004 ) . j.  maultzsch , s.  reich , c.  thomsen , h.  requardt , and p.  ordejn , phys . lett . * 92 * , 075501 ( 2004 ) . m.  calandra and f.  mauri , phys . b * 76 * , 205411 ( 2007 ) . a.  a.  abrikosov and s.  d.  beneslavskii , zh . fiz * 59 * , 1280 ( 1970 ) [ sov . jetp * 32 * , 699 ( 1971 ) ] . j.  gonzlez , f.  guinea , and m.  a.  h.  vozmediano , mod . b * 7 * , 1593 ( 1994 ) ; nucl . b * 424 * , 595 ( 1994 ) ; j.  low temp . phys . * 99 * , 287 ( 1994 ) . j.  gonzlez , f.  guinea , and m.  a.  h.  vozmediano , phys . b * 59 * , r2474 ( 1999 ) . i. l. aleiner , d. e. kharzeev , and a. m. tsvelik , phys . b * 76 * , 195415 ( 2007 ) j.  ye , phys . b * 60 * , 8290 ( 1999 ) . t.  stauber , f.  guinea , and m.  a.  h.  vozmediano , phys . b * 71 * , 041406(r ) ( 2005 ) . m.  s.  foster and i.  l.  aleiner , arxiv:0802.0283 . d.  m.  basko and i.  l.  aleiner , phys . b * 77 * , 041409(r ) ( 2008 ) . secs . [ sec : hamiltonian ] and  [ sec:1raman ] of the present paper were written in coauthorship with i.  l.  aleiner . m.  lax , _ symmetry principles in solid state and molecular physics _ ( dover , new york , 2001 or wiley , new york , 1974 ) . p.  r.  wallace , phys . rev . * 71 * , 622 ( 1947 ) . h.  suzuura and t.  ando , phys . * 89 * 266603 ( 2002 ) . l.  d.  landau and e.  m.  lifshits _ theory of elasticity _ ( pergamon press , 1970 ) ,  25 . i.  a.  lukyanchuk and a.  m.  bratkovsky , arxiv:0707.0466 . m.  lazzeri , s.  piscanec , f.  mauri , a.  c.  ferrari , and j.  robertson , phys . rev . b * 73 * , 155426 ( 2006 ) . m.  lazzeri and f.  mauri , phys . * 97 * , 266407 ( 2006 ) . a.  h.  castro neto and f.  guinea , phys . b * 75 * , 045404 ( 2007 ) . a.  grneis , r.  saito , t.  kimura , l.  g.  canado , m.  a.  pimenta , a.  jorio , a.  g.  souza filho , g.  dresselhaus , and m.  s.  dresselhaus , phys rev . b * 65 * , 155405 ( 2002 ) . v.  b.  berestetskii , e.  m.  lifshits , and l.  p.  pitaevskii , _ quantum electrodynamics _ ( pergamon press , 1982 ) ,  79 . j. maultzsch , s.  reich , and c.  thomsen , phys . rev . b * 70 * , 155403 ( 2004 ) . @xmath194 in graphite was considered long ago [ f.  guinea , j.  phys . c : solid state phys . * 14 * , 3345 ( 1981 ) ] , but the logarithmic divergence was not noted in that work . i.  l.  aleiner and k.  b.  efetov , phys . * 97 * , 236801 ( 2006 ) . e.  mccann , k.  kechedzhi , v.  i.  falko , h.  suzuura , t.  ando , and b.  l.  altshuler , phys . * 97 * , 146805 ( 2006 ) . e. i. blount , physical review * 126 * , 1636 ( 1962 ) . the idea of this calculation was suggested to the author by m.  s.  foster .
Abstract:

Answer: we present a detailed calculation of intensities of two - phonon and four - phonon raman peaks in graphene . writing the low - energy hamiltonian of the interaction of electrons with the crystal vibrations and the electromagnetic field from pure symmetry considerations , we describe the system in terms of just a few independent coupling constants , considered to be parameters of the theory . the electron scattering rate is introduced phenomenologically as another parameter . the results of the calculation are used to extract information about these parameters from the experimentally measured raman peak intensities . in particular , the raman intensities are sensitive to the electron scattering rate , which is not easy to measure by other techniques . also , the raman intensities depend on electron - phonon coupling constants ; to reproduce the experimental results , one has to take into account renormalization of these coupling constants by electron - electron interaction .