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 .


superconducting films are candidate substances for the improvement of electronics technology in a myriad of applications .
while the low resistance is very attractive in this regard , it has proved difficult to control the nonlinear behaviour of such materials in response to electromagnetic field@xcite . when a magnetic field is strong enough to penetrate into a superconductor in the form of quantised magnetic flux tubes , the vortex state obtains as a mixed state of superconducting phase punctuated by the vortices themselves .
vortices are surrounded by a supercurrent and can be forced into motion by the current resulting from an applied field .    as a topological defect
, a vortex is not only stable under perturbations@xcite but can not decay . the collection of vortices in a type - ii superconductor forms what is called vortex matter , and it is this which determines the physical properties of the system rather then the underlying material properties , in particular driving phase transitions@xcite . in the mixed state ,
a superconductor is not perfect ; it exhibits neither perfect diamagnetism nor zero electrical resistance .
the transport current @xmath0 generates a lorentz force @xmath1 on the vortex and forces it into motion , dissipating energy .    in reaching thermal equilibrium ,
energy is transferred via interactions between phonons and quasiparticle excitations .
small - scale imperfections such as defects scatter the quasiparticles , affecting their dynamics . in _
superconductors , impurities are plentiful and vortices experience a large friction .
this implies a fast momentum - relaxation process .
in contrast is the _ clean _
limit , where impurities are rare and no such relaxation process is available .
it is in this situation of slow relaxation that the hall effect appears .    generally , the @xmath2-@xmath3 phase diagram@xcite of the vortex matter has two phases . in the _ pinned phase _
vortices are trapped by an attractive potential due to the presence of large - scale defects , thus resistivity vanishes .
this phase contains what are known as glass states .
there is then the _ unpinned phase _ in which vortices can move when forced and so a finite resistivity appears .
this phase is also known as the flux - flow region and can be of two types .
one type is a liquid state where vortices can move independently ; the other type is a solid state in which vortices form a periodic abrikosov lattice@xcite resulting from their long - range interacton .
one model for the transition between the pinned and unpinned phases appears in @xcite .    in the unpinned phase ,
the system is driven from equilibrium and experiences a relaxation process .
there are several ways to describe such a system .
a microscopic description@xcite invoking interactions between a vortex and quasi - particle excitations at the vortex core provides a good understanding of friction and sports good agreement with experiments in the sparse - vortex region @xmath4 .
there is also a macroscopic description , the london approach , where vortices are treated either as interacting point - like particles or an elastic manifold subject to a pinning potential , driving force and friction@xcite . in the small - field region , vortices behave as an array of elastic strings .    in the dense - vortex region @xmath5 , where the magnetic field is nearly homogeneous due to overlap between vortices , ginzburg - landau ( gl ) theory , which describes the system as a field , provides a more reasonable model . in dynamical cases , time - dependent gl ( tdgl )
theory is appropriate@xcite ; in gl - type models , additional simplification can come from the lowest landau level ( lll ) approximation which has proven to be successful in the vicinity of the superconducting - normal ( s - n ) phase transition line @xmath6 .
this has been pursued in the static case@xcite ( without driving force ) and in the dynamic case with a time - independent transport current@xcite .
it may be noted that in the glass state , zero resistance within the lll approximation can not be attained@xcite .    based on tdgl theory , we will study the dynamical response of a dense vortex lattice forced into motion by an alternating current induced by an external electromagnetic field .
vortices are considered which are free from being pinned and thermally excited , which in addition to thermal noise would produce entanglement and bending .
we assume the vortices can transfer work done by an external field to a heat bath .
experimentally , a low - temperature superconductor far away from the clean limit is the best candidate for attaining these conditions .
we do not consider thermal fluctuation effects specific to high - temperature superconductors . in a dissipative system driven by a single - harmonic electric field @xmath7 , long after its saturation time we can expect the system to have settled into steady - state behaviour , where the vortices are vibrating periodically with some phase .
the tdgl model in the presence of external electromagnetic field is analysed and solved in . the dynamical s - n phase transition surface @xmath8 is located in @xmath9-space .
this surface coincides with the mean - field upper - critical field @xmath10 in the absence of the applied field , and with the phase - transition surface in the presence of the constant driving field considered by hu and thompson@xcite .
we will provide an analytical formalism for perturbative expansion in the distance to @xmath11 , valid in the flux - flow region .
the response of vortex matter forced into motion by the transport current is studied in .
the current - density distribution and the motion of vortices are treated in . in analysing the vortex lattice configuration in
, a method is utilised whereby the heat - generation rate is maximised .
next are discussed power dissipation , generation of higher harmonics , and the hall effect .
an experimental comparison is made in with far - infrared ( fir ) measurement on nbn .
finally , some conclusions are made in .
let us consider a dense vortex system prepared by exposing a type - ii superconducting material to a constant external magnetic field @xmath12 with magnitude @xmath13 .
we also select the @xmath14 axis of the superconductor to be in the @xmath15 direction .
let the superconductor carry an alternating electric current along the @xmath16 direction , generated by an electric field @xmath17 as shown in .
such a system when disturbed from its equilibrium state will undergo a relaxation process . for our system ,
the tdgl equation@xcite is a useful extension of the equilibrium gl theory .    in the dense - vortex region of the @xmath2-@xmath3 phase diagram ,
vortices overlap and a homogeneous magnetic field obtains . describing the response of such a system by a field ,
the order parameter @xmath18 in the gl approach , is more suitable than describing vortices as particle - like flux tubes , as is done in the london approach@xcite .
a strongly type - ii superconductor is characterised by its large penetration depth @xmath19 and small coherence length @xmath20 , @xmath21 .
the difference between induced magnetic field and external magnetic field is @xmath22 . in the vicinity of the phase - transition line
@xmath10 vortices overlap significantly , and @xmath23 making @xmath24 small . in this case
, the magnetic field may be treated as homogeneous within the sample .
we will have in mind an experimental arrangement using a planar sample very thin compared with its lateral dimensions .
since the characteristic length for inhomogeneity of electric field@xcite @xmath25 is then typically large compared with sample thickness , this implies that the electric field may also be treated as homogeneous throughout@xcite , eliminating the need to consider maxwell s equations explicitly .    in equilibrium
, the gibbs free energy of the system is given by@xcite @xmath26&=&\int\td { \bf r } \bigg\ { \frac{\hbar^2}{2m_{ab}}|{\bf d}\psi|^2 + \frac{\hbar^2}{2 m_c}|\partial_z \psi|^2 \nonumber\\ & & { \ \ \ \ } -\alpha ( \tnst - t)|\psi|^2+\frac{\beta}2 |\psi|^4 \bigg\}\end{aligned}\ ] ] where @xmath27 is the critical temperature at zero field .
covariant derivatives employed here preserve local gauge symmetry and are two - dimensional ; @xmath28 and @xmath29 .
governing the dynamics of the field @xmath18 is the tdgl equation @xmath30 this determines the characteristic relaxation time of the order parameter .
microscopic derivation of tdgl can be found in @xcite in which the values of @xmath31 , @xmath32 and @xmath33 are studied . in the macroscopic case ,
these are viewed simply as parameters of the model . at microscopic scale
, disorder is accounted for by @xmath33 , the inverse of the diffusion constant ; the relation of @xmath33 to normal - state conductivity is discussed in .    in standard fashion , @xmath34
while @xmath35 .
our set of equations is completed@xcite by including ampre s law , writing for the total current density @xmath36 as we shortly make a rescaling of quantities , we have written @xmath37 subscripts here for clarity .
the first term is the normal - state conductivity .
the second term can be written using a maxwell - type equation relating the vector - potential with the supercurrent , @xmath38 .\end{aligned}\ ] ]    this is a gauge - invariant model ; we fix the gauge by considering the explicit vector potential @xmath39 and @xmath40 , corresponding to an alternating transport current .
each vortex lattice cell has exactly one fluxon .
we do not assume the electric field and the motion of vortices are in any particular direction relative to the vortex lattice , by way of rendering visible any anisotropy .    for convenience ,
we define some rescaled quantities .
the rescaled temperature and magnetic field are @xmath41 and @xmath42 .
@xmath43 denotes the mean field upper - critical field , extrapolated from the @xmath27 region down to zero temperature .    in the @xmath44-@xmath45 plane of the crystal we make use of _ magnetic length
_ @xmath46 .
we define @xmath47 where @xmath48 .
the scale on the @xmath14-axis is @xmath49 with @xmath50 .
the coordinate anisotropy in @xmath15 is absorbed into this choice of normalisation , as can be seen in .
the order parameter @xmath18 is scaled by @xmath51 .
the time scale is normalised as @xmath52 .
therefore , frequency is @xmath53 .
note that @xmath54 is then inversely proportional to @xmath45 .
the amplitude of the external electric field is normalised with @xmath55 so that @xmath56 .
after our rescaling the tdgl equation takes the simple form @xmath57 where the operator @xmath58 is defined as @xmath59 with our specified vector potential , covariant derivatives are @xmath60 , @xmath61 and @xmath62 .
we define @xmath63 for convenience .
the tdgl equation is invariant under translation in @xmath15 , thus the dependence of the solution in the @xmath15 direction can be decoupled .
@xmath58 is not hermitian ; @xmath64 where the conjugation is with respect to the usual inner product , defined below .
we will make extensive use of the eigenfunctions of @xmath58 and @xmath65 in what follows .
the eigenvalue equation @xmath66 defines the set of eigenfunctions of @xmath58 appropriate for our analysis ; this can be seen in .
the convention is that @xmath67 when and only when @xmath68 . taking corresponding in @xmath58 and switches it in the resulting @xmath69 to get the ` corresponding ' eigenfunction @xmath70 for @xmath65 . ]
eigenfunctions of @xmath65 to be @xmath71 , the orthonormality @xmath72 may be chosen , so long as @xmath73 . shown in ,
crystal structure determines linear combinations of these basis elements with respect to @xmath74 ; the resulting @xmath75 functions are then useful for expansion purposes below .
the inner product is @xmath76 where the brackets @xmath77 denote an integral over space and time . to define averages over only time or space alone ,
we write @xmath78 or @xmath79 respectively .    ;
( b ) typical small angle , plotted for @xmath80 . in the static case ,
both @xmath81 and @xmath80 are two particular angles will correspond to the same energy .
the applied constant magnetic field @xmath82 is along the @xmath83 direction and time - dependent oscillating electric field @xmath84 is in the @xmath16 direction .
@xmath85 is the the apex angle of the two defining lattice vectors .
the two vectors for the rhombic vortex lattice are @xmath86 and @xmath87 where @xmath88 .
the motion of vibrating vortices , indicated by the red arrow , is disccused in .
, title="fig:",width=151 ] ; ( b ) typical small angle , plotted for @xmath80 . in the static case ,
both @xmath81 and @xmath80 are two particular angles will correspond to the same energy .
the applied constant magnetic field @xmath82 is along the @xmath83 direction and time - dependent oscillating electric field @xmath84 is in the @xmath16 direction .
@xmath85 is the the apex angle of the two defining lattice vectors .
the two vectors for the rhombic vortex lattice are @xmath86 and @xmath87 where @xmath88 .
the motion of vibrating vortices , indicated by the red arrow , is disccused in .
, title="fig:",width=151 ]      states of the system can be parametrised by @xmath89 . by changing temperature @xmath90 , a system with some fixed @xmath91
may experience a normal - superconducting phase transition as temperature passes below a critical value @xmath92 .
such a point of transition is also known as a _
bifurcation point_.    the material is said to be in the normal phase when @xmath18 vanishes everywhere ; otherwise the superconducting phase obtains , with @xmath18 describing the vortex matter . because of the vortices , the resistivity in the superconducting phase need not be zero .
the s - n phase - transition boundary @xmath92 separates the two phases . to study the condensate
, we will use a bifurcation expansion to solve .
we expand @xmath18 in powers of distance from the phase transition boundary @xmath93 .      as in the static case
, we can locate the dynamical phase - transition boundary by means of the linearised tdgl equation@xcite .
this is because the order parameter vanishes at the phase transition , and we do not need to consider the nonlinear term . the linearised tdgl equation is written @xmath94    of the eigenvalues of @xmath58 , only the smallest one @xmath95 , corresponding to the highest superconducting temperature @xmath93 , has physical meaning .
the s - n phase transition occurs when the trajectory in parameter space intersects with the surface @xmath96 where the lowest eigenvalue is calculated in ; @xmath97 utilising a @xmath45-independent frequency @xmath98 and amplitude @xmath99 of input signal , we write @xmath100 in the absence of external driving field , @xmath101 , the phase - transition surface coincides with the well - known static - phase transition line @xmath102 in the mean - field approach . with time - independent electric field at @xmath103 , where the vortex lattice is driven by a fixed direction of current flow ,
the dynamical phase - transition surface coincides with that proposed in @xcite , but with a factor of @xmath104 .
this amplitude difference is familiar from elementary comparisons of dc and ac circuits .    in the above equation
, we can see that in the static case @xmath101 , the superconducting region is @xmath105 .
in addition when @xmath106 , the superconducting region in the @xmath45-@xmath90 plane is smaller than the corresponding region in the static case , as can be seen in ( a ) .    finally , increasing frequency will increase the size of the superconducting region , as in ( b ) ; in the high - frequency limit , the area will reach its maximum , which is the superconducing area from the static case . as with any damped system ,
response is diminished at higher frequencies .
the superconducting state does not survive at small magnetic field ; for example at @xmath107 in ( a ) , the material is in the normal state over most of the @xmath2-@xmath3 phase diagram . later in this paper
we will consider interpretation of this phenomenon .
in particular , when discussing energy dissipation in , we will see that the main contribution to the dissipation is via the centre of the vortex core . at small magnetic field , since there are fewer cores to dissipate the work done by the electric field , the superconducting state is destroyed and the order parameter vanishes .     as a function of @xmath45 for various @xmath99 at @xmath108 and
( b ) @xmath93 as a function of @xmath99 for various @xmath109 at @xmath110 .
the straight line in ( a ) is the @xmath101 curve and corresponds to the mean - field phase - transition line @xmath102 .
states above each line are normal phase while the region below each line is superconducting .
@xmath99 suppresses the superconducting phase as shown in ( a ) , while @xmath109 removes this suppression effect , as shown in ( b ) .
, title="fig:",width=158 ]   as a function of @xmath45 for various @xmath99 at @xmath108 and ( b ) @xmath93 as a function of @xmath99 for various @xmath109 at @xmath110 .
the straight line in ( a ) is the @xmath101 curve and corresponds to the mean - field phase - transition line @xmath102 .
states above each line are normal phase while the region below each line is superconducting . @xmath99 suppresses the superconducting phase as shown in ( a ) , while @xmath109 removes this suppression effect , as shown in ( b ) . ,
title="fig:",width=158 ]      that the vortex matter dominates the physical properties of the system is especially pronounced in the pinning - free flux - flow region . here
we solve by a bifurcation expansion@xcite .
since the amplitude of the solution grows when the system departs from the phase transition surface where @xmath111 , we can define a distance from this surface as @xmath112 and expand @xmath18 in @xmath113 .
the tdgl in terms of @xmath113 is @xmath114 where @xmath115 is the operator @xmath58 shifted by its smallest eigenvalue .
@xmath18 is then written @xmath116 and it is convenient to expand @xmath117 in terms of our eigenfunctions of @xmath118 @xmath119    in principle , all coefficients @xmath120 in can be obtained by using the orthogonal properties of the basis , which are explained in . inserting @xmath18 from equation into tdgl equation , and collecting terms with the same order of @xmath113
, we find that for @xmath121 @xmath122 and for @xmath123 @xmath124 for @xmath125 @xmath126 and so on .
observing , the solution for the equation is @xmath127 where @xmath128 is a particular linear combination of all eigenfunctions with the smallest eigenvalue .
the coefficient of @xmath129 can be obtained by calculating the inner product of @xmath130 with , @xmath131 in the same way , the coefficient of the next order @xmath132 , can be obtained by finding the inner product of @xmath133 with the @xmath125 equation , @xmath134 the inner product of @xmath135 on gives the coefficient for @xmath136 @xmath137 and @xmath138 the solution of tdgl is then @xmath139    in this paper we will restrict our discussion to the region near @xmath93 where the next - order correction can be disregarded ; @xmath140    we would like to emphasise that our discussion at this order is valid in the vicinity of the phase - transition boundary and in particular for a superconducting system without vortex pinning . in such a system , vortices move in a viscous way , resulting in flux - flow resistivity ; no divergence of conductivity is expected .
our results based on were calculated at @xmath129 order , where only the lowest eigenvalue @xmath141 of the tdgl operator @xmath58 makes an appearance .    the next - order correction is at order @xmath132 , and there is now a contribution from higher landau levels . from the symmetry argument in @xcite ,
as long as the hexagonal lattice remains the stable configuration for the system , the next - order contribution comes from the sixth landau level with a factor @xmath142 . even in the putative case of a lattice deformed slightly away from a hexagonal configuration ,
the next contributing term is @xmath143 , since in our system the lattice will remain rhombic .
the vortex lattice has been experimentally observed since the 1960s and its long - range correlations have been clearly observed@xcite with dislocation fraction of the order @xmath144 .
remarkably , the same techniques can be used to study the structure and orientation of moving vortex lattice with steady current@xcite , and with alternating current in the small - frequency regime@xcite . in this subsection
, we will discuss the configuration of the vortex lattice in the presence of alternating transport current in the long - time limit .
in the dynamical case , the presence of an electric field breaks the rotational symmetry of an effectively isotropic system to the discrete symmetry @xmath145 .
in contrast , a rhombic lattice preserves at least a symmetry of this kind along two axes , and the special case of a hexagonal lattice preserves sixfold symmetry .
the area of a vortex cell is determined by the quantised flux in the vortex , which is @xmath146 in terms of our rescaled variables .
as shown in , we choose a unit cell @xmath147 defined by two elementary vectors @xmath148 and @xmath149 .
we will first construct a solution for an arbitrary rhombic lattice parameterised by an apex angle @xmath85 .
consideration of translational symmetry in the @xmath150 direction leads to the discrete parameter @xmath151 . in
we show that in the long - time limit the lowest - eigenvalue steady - state eigenfunctions of @xmath58 must therefore combine to form @xmath152{2 \sigma}\sum_{l=-\infty}^{\infty }          e^{i \frac{\pi}2 l(l-1 ) }          e^{i k_l ( x- v \sin\omega \tau/\omega ) } u_{k_l}(y,\tau )      .\ ] ] here @xmath128 is normalised as @xmath153 the function @xmath154 is given by @xmath155 ^ 2                  }      , \end{aligned}\ ] ] with @xmath156              }      .\end{aligned}\ ] ] in analogy with a forced vibrating system in mechanics , a phase @xmath157 and a reduced velocity @xmath158 have been introduced for convenience in .
the zero electric - field limit , large - frequecy limit and zero frequency limit are consistent with previous studies concluded in .    in our approximation ,
the @xmath159 in is a time - independent quantity from and @xmath160 where @xmath161 .    in the small signal limit @xmath162
, @xmath159 reduces to the abrikosov constant .
the abrikosov constant with either @xmath80 or @xmath163 minimises the gl free energy in the static state@xcite . to be more explicit
, @xmath159 can be expanded in terms of the amplitude of input signal . in powers of @xmath164 , @xmath165 and we find it convenient to write in terms of @xmath166 .
the first term in @xmath159 is the abrikosov constant @xmath167 for hexagonal lattices @xmath168 , whereas for a square lattice @xmath169 .
the next term in @xmath159 is @xmath170 with a coefficient @xmath171 we see that at high frequency , the correction in higher order terms of @xmath166 can be disregarded .
in this section we discuss the current distribution and motion of vortices , energy transformation of the work done on the system into heat , nonlinear response and finally the hall effect .
in addition to the conventional conductivity attributable to the normal state , there is an overwhelming contribution due to the superconducting condensate in the flux - flow regime , tempered only by the dissipative properties of the vortex matter . in this section
we will examine the supercurrent density to investigate the motion of the vortex lattice .
we consider a hexagonal lattice in a fully dissipative system ; the non - dissipative part known as the hall effect will be discussed in .
the supercurrent density @xmath172 is obtained by substitution of the solution into .
@xmath173 and @xmath174 where @xmath175 and @xmath154 is given in . observing , we conceptually split the current into two components .
one part is the circulating current surrounding the moving vortex core as in the static case ; we refer to this component as the _ diamagnetic current_. the other part which we term the _ transport current _ is the component which forces vortices into motion .    :
( a ) @xmath176 ; vortex cores move to the right .
( b ) @xmath177 ; vortex cores move to the left .
vortices are drawn back and forth as the direction of the transport current density alternates .
the magnitude of the current density has maximal regions which tend to circumscribe the cores ; the maxima in these regions move in the plane and their manner of motion can be described as leading the motion of the vortices by a small phase .
the average current in a unit cell leads the motion of the vortex in time by a phase of @xmath178 .
, title="fig:",width=158 ] : ( a ) @xmath176 ; vortex cores move to the right .
( b ) @xmath177 ; vortex cores move to the left .
vortices are drawn back and forth as the direction of the transport current density alternates .
the magnitude of the current density has maximal regions which tend to circumscribe the cores ; the maxima in these regions move in the plane and their manner of motion can be described as leading the motion of the vortices by a small phase .
the average current in a unit cell leads the motion of the vortex in time by a phase of @xmath178 .
, title="fig:",width=158 ]    the diamagnetic current may be excised from our consideration by integrating the current density over the unit cell @xmath147 ; that is , we consider @xmath179 .
we have @xmath180 and @xmath181 with our conventions , the transport current is along the @xmath16-direction . considering the lorentz force between the magnetic flux in the vortices and the transport current , we expect the force on the vortex lattice to be perpendicular .
we identify the locations of vortex cores to be where @xmath182 .
the velocity of the vortex cores turns out to be @xmath183 along the @xmath150-direction .
note that vortex lattice moves coherently .
the vortex motion the electric field with a phase @xmath184 which increases with frequency and reaches @xmath178 asymptotically .
the maximum velocity of vortex motion @xmath185 decreases with increasing frequency .    in
we show the current distribution and the resultant oscillation of vortices .
as anticipated , the transport current and the motion of vortices are perpendicular as the vortices follow the input signal .
the current density diminishes near the core ; it is small there compared to its average value .    in steady - state motion , since the vortices move coherently in our approximation , the interaction force between vortices is balanced as in the static case .
since the system is entirely dissipative , the motion that the vortices collectively undergo is viscous flow .
the vortex lattice responds to the lorentz driving force as a damped oscillator , and this is the origin of the frequency - dependent response .      in static case
the system is described by the gl equation . solving this equation , which is but with zero on the left - hand side ,
will select some lattice configuration .
the global minima of the free energy correspond to a hexagonal lattice , while there may be other configurations producing local minima . in the static case
the lattice configuration can be determined in practice by building an ansatz from the linearized gl solution@xcite and then using a variational procedure to minimise the full free energy .    in the dynamic case
, there is no free energy to minimise ; we must embrace another method of making a physical prediction regarding the vortex lattice configuration .
let us follow @xcite and take as the preferred structure the one with highest heat - generation rate .
though we have at present no precise derivation , our physical justification of this prescription is that the system driven out of equilibrium can reach steady - state and stay in condensate only if the system can efficiently dissipate the work done by the driving force .
therefore , whatever the cause , the lattice structure most conducive to the maintenance of the superconducting state will correspond to the maximal heat generation rate .
the heat - generating rate@xcite is @xmath186          \big\ }      .\end{aligned}\ ] ] @xmath159 is given explictly in and is the only parameter involving the apex angle @xmath85 of the moving vortex lattice . here
@xmath159 plays the same rle as the abrikosov constant @xmath187 in the static case .
corresponding to maximising the heat - generating rate , the preferred structure can be obtained by simply minimising @xmath159 with respect to @xmath85 .
this shows from the current viewpoint of maximal heat - generation rate that vortices are again expected to move coherently .    in or it
is seen that the moving lattice is distorted by the external electric field but this influence subsides at high frequency .
numerical solution for minimising @xmath85 shows that while near the high - frequency limit there remain two local minima for @xmath159 with respect to @xmath85 , the solution near @xmath81 is favoured slightly over that at @xmath188 as the global minimum .
this is as presented in ( a ) .
the two minima tend to approach each other slightly as the frequency begins to decrease further .
in an experimental setting , this provides an avenue for testing the empirical validity of the maximal heat generation prescription , in particular in terms of the direction of lattice movement@xcite .
we put forth the physical interpretation that at high frequency the friction force becomes less important , and the distortion is lessened .
since interactions dominate the lattice structure the system at high frequency will have many similarities with the static case .
energy supplied by the applied alternating current is absorbed and dissipated by the vortex matter , and the heat generation does not necessarily occur when and where the energy is first supplied . in ,
we show an example of this transportation of energy by the condensate . on the left
is shown a contour plot of the work @xmath189 done by the input signal ; points along a given contour are of equal power absorption . on the right of
is shown the heat - generating rate@xcite , @xmath190 .
the periodic maximal regions are near the vortex cores in both patterns .     and @xmath191 : ( a ) work @xmath192 and ( b ) heat generating rate @xmath193 .
the vortex cores are denoted by ` @xmath194 ' in both figures , shown in the @xmath150-@xmath16 plane .
the maximum displacements of vortex cores are shown by the arrow .
the maximal region around the core in ( a ) is elongated by the current .
the similar horizontal broadening around the core in ( b ) is caused by the vortex motion .
energy is transported ; maxima in ( a ) and ( b ) do not coincide . , title="fig:",width=158 ]   and @xmath191 : ( a ) work @xmath192 and ( b ) heat generating rate @xmath193 .
the vortex cores are denoted by ` @xmath194 ' in both figures , shown in the @xmath150-@xmath16 plane .
the maximum displacements of vortex cores are shown by the arrow .
the maximal region around the core in ( a ) is elongated by the current .
the similar horizontal broadening around the core in ( b ) is caused by the vortex motion .
energy is transported ; maxima in ( a ) and ( b ) do not coincide . ,
title="fig:",width=158 ]    in ( b ) , one can see that the system dissipates energy via vortex cores . from a microscopic point of view , cooper pairs break into quasiparticles inside the core ; these couple to the crystal lattice through phonons and impurities to transfer heat .
the interaction between vortices and excitation of vortex cores manifests as friction@xcite .
the power loss of the system averaged over time and space is @xmath195 .
@xmath196      \end{aligned}\ ] ] where @xmath197 is a bessel function of the first kind .
( upper panel ) and expansion parameter @xmath113 ( lower panel ) as a function of frequency @xmath109.,width=340 ]    in is shown the power loss and also @xmath113 as a function of frequency .
@xmath113 is proportional the density of cooper pairs , and can be thought of as an indication of how robust is the superconductivity .    as frequency increases , while @xmath113 tends to an asymptotic value , @xmath198 achieves a maximum and then decreases ; this maximum is due to fluctuations of order parameter caused by the input signal . in a fully dissipative system as considered here , the maximum of each curve is not a resonance phenomenon but is instead caused by fluctuation of the order parameter resulting from the influence of the applied field .    a parallel may be drawn between what we have observed in this section and the suppression of the superconductivity by macroscopic thermal fluctuations commonly observed in high - temperature superconductors . in our case ,
the vortices in a high-@xmath199 superconductor undergo oscillation due to the driving force of the external field .
we may think of this as being analogous to the fluctuations of vortices due to thermal effects alone in a low-@xmath199 superconductor .
although the method of excitation is different , the external electromagnetic perturbation in the present case essentially plays the same rle as the thermal fluctuations in low-@xmath199 situation .
finally , we point out that @xmath113 seems to be an appropriate parameter for determining the amount of power loss .
generically , it seems that for points deeped inside the superconducting region , that is at large @xmath113 compared with its saturation value at high @xmath109 , the power loss due to the dissipative effects of the vortex matter becomes suppressed .
we suggest the possibility that this effect , which is navely intuitive , is in fact physical and more widely applicable than merely the present model .
the practical application of superconducting materials is dependent on how well one can control the inherent nonlinear behaviour . in this section we will focus on the generation of higher harmonics in the mixed state , in response to a single - frequency input signal .
the periodic transport current @xmath200 is an odd function of input signal , and it turns out that the response motion also contains only odd harmonics . from we can calculate the fourier expansion for transport current .
@xmath201      , \ ] ] where the fourier coefficient @xmath202 is @xmath203      e^{-i(2n+1)\theta      } .\ ] ]
we see the response goes beyond simple ohmic behaviour and the coefficients are proportional to @xmath113 .
experimentally , one way of measuring these coefficients is a lock - in technique@xcite which is adept at extracting a signal with a known wave from even an extremely noisy environment .    to make contact with more standard parameters and satisfy our intuition ,
we expand the first two harmonics in terms of @xmath164 . the fundamental harmonic , @xmath204 expanded in powers of @xmath205 is @xmath206 where @xmath207 .
the first term is the ohmic conductivity denoted as @xmath208 , and is reminiscent of drude conductivity for free charged particles .
this is not an unexpected parallel , since the cooper pairs in a superconducting system can be imagined to behave like a free - particle gas .
taking this viewpoint , in the small - signal limit , the ratio @xmath209 gives the relaxation time of the charged particles .
subsequent higher - order corrections all contain @xmath54 in such a way that their contributions are suppressed at large @xmath54 .
the coefficient of the @xmath210 harmonic expanded in powers of @xmath205 is @xmath211 which decreases quickly with increasing @xmath54 .    in
, we show the generation of higher harmonics for three different states in the dynamical phase diagram . for each harmonic labeled by @xmath212 , @xmath213 as a function of @xmath109 has the same onset as @xmath113 .
we can see that @xmath213 reaches a maximum and then starts to decay while @xmath113 saturates .
the coefficients of harmonics with @xmath214 decay to zero in the high @xmath109 limit , where the state is well inside the superconducting region .
we pointed out in and reaffirm here that @xmath113 plays a significant rle in determining the extent of nonlinearity in the system . in turn ,
the parameter which controls this is @xmath54 .
when @xmath54 is large , @xmath113 is brought closer to its saturation value @xmath215 , causing the higher harmonics to be suppressed , and also lessening distortion of the vortex lattice . finally , for a given harmonic
, @xmath213 is generally smaller when @xmath215 is smaller ; this can be seen by comparing ( a ) and ( c ) of .
one might point out that the nonlinear behaviour is decreased at , for example , large @xmath54 .
nevertheless , we view the parameter @xmath216 as more intrinsic to the system , rather than simply characterising the input signal .    a limited parallel
can be drawn between the effect of thermal noise in high-@xmath199 superconducting systems , and the effect of the electromagnetic perturbation in our present case .
it seems that in either case the fluctuation influence can be reduced by moving the state deeper inside the superconducting region .
[ r][r]@xmath217 [ r][r]@xmath218 [ r][r]@xmath219 [ r][r]@xmath220      in contrast to the fully dissipative system we have considered , in this section we will discuss an effect caused by the non - dissipative component , namely the hall effect . in a clean system ,
vortices move without dissipation ; a transverse electric field with respect to current appears .
the non - dissipative part is subject to a gross - pitaevskii description , using a type of nonlinear schrdinger equation@xcite , with a non - dissipative part to the relaxation @xmath33 from . the fully dissipative operator @xmath58 in our previous discussion can be generalised by using a complex relaxation coefficient @xmath221 .
we thus define @xmath222 the ratio @xmath223 is typically on the order of @xmath224 for a conventional superconductor , and @xmath225 for a high-@xmath199 superconductor@xcite .
the hall effect is small here . in normal metals ,
the non - dissipative part gives the cyclotron frequency .
if @xmath226 is the relaxation time of a free electron in a dirty metal , then for typical values of @xmath227 the hall effect becomes negligible . because the supply of conducting electrons is limited , the transverse component increases at the expense of the longitudinal component as the mean free path of excitations grows
it is equivalent to an increase in the imaginary part of the relaxation constant at the expense of the real part .
the eigenvalues and eigenfunctions of @xmath228 can be obtained easily by replacing the @xmath164 in previous results with @xmath229 , @xmath54 with @xmath230 and @xmath231 with @xmath232 .
the transport current along the @xmath150-direction is no longer zero in the presence of the non - dissipative component ; it is propotional to @xmath233 .
the frequency - dependent hall conductivity can be obtained from the first - order expension in @xmath164 , @xmath234 while the hall contribution in the @xmath16 direction is expected to be negligible , as it is of the order of @xmath235 .    in principle
, the crossover between non - dissipative systems and dissipative systems can be tuned using the ratio @xmath233 . in a non - dissipative system , which is the clean limit , the hall effect is important and taking account of the imaginary part of tdgl is necessary . on the contrary , in a strong dissipative system where excitations are in thermal equilibrium via scattering ,
the tdgl equation gives satisfactory agreement .
far - infrared spectroscopy can be performed using monochromatic radiation which is pulsed at a high rate , known as fast far - infrared spectroscopy .
this technique sports the advantage of avoiding overheating in the system , making it a very effective tool in observing the dynamical response of vortices .
in particular , one can study the imaginary part of conductivity contributed mainly from superconducting component .    in
is shown a comparison with an nbn experiment measuring the imaginary part of conductivity .
the sample has the gap energy @xmath236  mev .
the resulting value of @xmath237 is larger than the value expected from bcs theory@xcite .
we consider frequency - dependent conductivity in the case of linearly polarised incident light with a uniform magnetic field along the @xmath15 axis .
the theoretical conductivity contains both a superconducting and a normal contribution .
the total conductivity is obtained from the total current as in where the normal - part conductivity in the condenstate is the conductivity appearing in the drude model .    according to our previous discussion ,
the nonlinear effect of the input signal on nbn is unimportant in the thz region , which corresponds to @xmath238 .
an approximation where the flux - flow conductivity includes only the term @xmath239 from , and hall coefficient @xmath240 from is shown in and the agreement with experiment is good .
the nave way in which we have treated the normal - part contribution is essentially inapplicable to the real - part conductivity .
this is because the real - part conductivity contains information about interactions with the quasi - particles inside the core , making further consideration necessary@xcite .
[ r][r]@xmath241  k and @xmath242  t calculated using of @xcite .
normal - state conductivity is @xmath243  @xmath244 and relaxation time of an electron @xmath245  fs are taken from @xcite .
the theoretical curve has one fitting parameter @xmath246.,title="fig:",width=226 ]
the time - dependent ginzburg - landau equation has been solved analytically to study the dynamical response of the free vortex lattice .
based on the bifurcation method , which involves an expansion in the distance to the phase transition boundary , we obtained a perturbative solution to all orders .
we studied the response of the vortex lattice in the flux - flow region just below the phase transition , at first order in this expansion .
we have seen that there are certain parameters which can be tuned using the applied field and temperature , providing a feasible superconducting system where one can study precise control of nonlinear phenomena in vortex matter .    under a perturbation by electromagnetic waves
, the steady - state solution shows that there is a diamagnetic current circulating the vortex core , and a transport current parallel to the external electric field with a frequency - dependent phase shift and amplitude .
vortices move perpendicularly to the transport current and coherently .    using a technique of maximising the heat - generation rate
, we showed that the preferred structure based on energy dissipation is a hexagonal lattice , with a certain level of distortion appearing as the signal is increased or the frequency is lowered .
energy flowing into the system via the applied field is dissipated through the vortex cores .
we showed that the superconducting part may be thought of as having inductance in space and time .
we have written transport current beyond a simple linear expression .
a comparison between different harmonics of three different states in our four - dimensional parameter space indicated that the nonlinearity becomes unimportant at high frequency and small amplitude , and the influence of the input signal is decreased when the system moves deeper inside the superconducting region , away from the phase - transition boundary .    to observe the configuration of moving vortices , techniques such as muon - spin rotation@xcite , sans@xcite , stm@xcite and others@xcite seem to be promising options . to provide the kind of input signal considered here
, methods such as short - pulse fir spectroscopy as used in @xcite might be applied .
the coefficient we defined in corresponds to conductivity .
we have also seen that a simple parametrisation by complex quantities like conductivity and surface impedence is insufficient to capture the detailed behaviour of the system ; in performing experiments , it should be kept in mind that the nonlinearity can be measured in terms of more appropriate variables as we have shown .
we have viewed the forcing of the system by the applied field to be somewhat analogous to thermal fluctuations , in the sense that they both result in vibration of the vortex lattice .
hence , the influence of the electromagnetic fluctuation is stronger at the nucleation region of superconductivity than deep inside the superconducting phase . besides , since at high frequency
the motion of vortices is limited , the influence from electric field is suppressed , as is the hall effect .    fruitful discussions with b.  rosenstein , v.  zhuravlev , and j.  kolek are greatly appreciated .
the authors kindly thank g.  bel and p.  lipavsk for critical reading of the manuscript and many useful comments .
the authors also have benefitted from comments of a.  gurevich .
we thank j.  r.  clem for pointing out to us reference @xcite .    `
nsc99 - 2911-i-216 - 001 `
we follow @xcite to estimate the coefficient @xmath33 which characterises the relaxation process of the order parameter .
@xmath33 is the inverse of the diffusion coefficient for electrons in the normal state . for a strongly - scattering system , as in the dirty limit@xcite
, the ratio between the relaxation times of order parameter @xmath247 and the vector potential ( or current ) @xmath248 is @xmath249 by definition of the thermal critical field @xmath250 and @xmath251 , as in @xcite , we know the ratio of the two parameters is @xmath252 the coherence length at zero temperature can be written in terms of @xmath43 as @xmath253 and in terms of effective mass as @xmath254 . as a result ,
@xmath255 with @xmath33 , we can retrieve the experimental quantities from the rescaled ones used in calculation . using the @xmath37-subscripted original variables , we write the electric field @xmath256 , and the frequency @xmath257 . the current density @xmath258 . in the case of linear response
, we have @xmath259 where @xmath260 .
we consider the linearised time - dependent ginzburg - landau equation , which has been written in our chosen gauge .
we wish to find the set of eigenfunctions of @xmath58 corresponding to the lowest eigenvalue .
based on knowledge of the solution in the static case , we solve the now time - dependent problem by making the following ansatz@xcite . the electric field along the @xmath16-direction breaks rotational symmetry in the @xmath150-@xmath16 plane , so we write @xmath261    after substitution of @xmath262 for @xmath263 in , comparison of coefficients of powers of @xmath16 gives the following differential equations in @xmath231 .
@xmath264 the solutions are @xmath265 .\ ] ] for a steady - state solution , @xmath266 , we have @xmath267 . @xmath268 as with @xmath269 , we have here @xmath270 .
@xmath271 where @xmath272 is a normalisation constant .
the resulting eigenvalue is @xmath273 now , although in a more realistic treatment of the system , one may introduce some boundary condition restricting @xmath274 to a certain set of values , here we simply select the smallest eigenvalue @xmath275 available to us by setting @xmath274 to zero .
thus equipped with the set of eigenfunctions corresponding to our lowest eigenvalue , we deem them to be the first elements of our basis , labeled by @xmath141 and @xmath74 .
these eigenfunctions of @xmath58 are @xmath276 with @xmath277 ^ 2                  }      , \end{aligned}\ ] ] and @xmath278              } .\end{aligned}\ ] ]    in the same way , the corresponding eigenfunctions of @xmath279 can be obtained .
@xmath280 with @xmath281 where @xmath282              }      \end{aligned}\ ] ] according to our normalisation condition .
the lowest eigenvalue of @xmath279 is @xmath283 .    in the @xmath284 limit
, the system reduces to the case of constant electric field .
the eigenfunctions and eigenvalues are then consistent with those obtained by hu and thomson@xcite . in the limit of zero electric field ,
the eigenfunctions and eigenvalues reduce to those of the lowest landau level static - state solution@xcite .
this is also the same in the @xmath285 limit .
99 d.  e.  oates , _ overview of nonlinearity in hts : what we have learned and prospects for improvement _ , j.  supercond .
novel magn .
* 20 * , 3 ( 2007 ) ; j.  c.  booth , s.  a.  schima , d.  c.  degroot , _ description of the nonlinear behavior of superconductors using a complex conductivty _ , ieee trans .
appl .  supercond .
* 13 * , 315 ( 2003 ) ; a.  v.  velichoko , m.  courier , j.  lancaster , a.  porch , _ nonlinear microwave properties of high @xmath199 thin films _ , supercond .
technol .  * 18 * r24 ( 2005 ) m.  tinkham , _ introduction to superconductivity _ , mcgraw - hill , new york ( 1996 ) d.  charalambous , e.  m.  forgan , s.  ramos , s.  p.  brown , r.  j.  lycett , d.  h.  ucko , a.  j.  drew , s.  l.  lee , d.  fort , a.  amato , u.  zimmerman , _ driven vortices in type - ii superconductors : a muon spin rotation study _ , phys .  rev .
b * 73 * , 104514 ( 2006 ) b.  rosenstein , d.  li , _ ginzburg - landau theory of type ii superconductors in magnetic field _ , rev .
phys .  * 82 * , 109 ( 2010 ) and references therein a.  a.  abrikosov , zh .  eksperim , i teor .
fiz .  * 32 * , 1442 ( 1957 ) a.  gurevich , e.  h.  brandt , _ ac response of thin superconductors in the flux - creep regime _ , phys .
b * 55 * , 12706 ( 1997 ) n.  b.  kopnin , _ theory of nonequilibrium superconductivity _ , university press , oxford ( 2001 ) m.  w.  coffey , j.  r.  clem , _ unified theory of effects of vortex pinning and flux creep upon the rf surface impedance of type - ii superconductors _ , phys .
lett .  * 67 * , 386 ( 1991 ) j.  i.  gittleman , b.  rosenblum , _ radio - frequency resistance in the mixed state for subcritical currents _ , phys .
* 16 * , 734 ( 1966 ) c .- r .
hu , r.  s.  thompson , _ dynamic structure of vortices in superconductors .
ii . @xmath4 _ , phys .
b * 6 * , 110 ( 1972 ) ; r.  s.  thompson , c .- r .
hu , _ dynamic structure of vortices in superconductors _ , phys .
* 27 * , 1352 ( 1971 ) g.  w.  crabtree , d.  o.  gunter , h.  g.  kaper , a.  e  koshelve , g.  k. leaf , v.  m.  vinokur , _ numerical simulation of driven vortex systems _ , phys .
b * 61 * , 1146 ( 2000 ) ; d.  y.  vodolazov , f.  m.  peeters , _ rearrangement of the vortex lattice due to instabilities of vortex flow _ , phys .
b * 76 * , 014521 ( 2007 ) d.  li , a.  malkin , b.  rosenstein , _ structure and orientation of the moving vortex lattice in clean type - ii superconductors _ , phys .
b * 70 * , 214529 ( 2004 ) b.  rosenstein and v.  zhuravlev , _ quantitative theory of transport in vortex matter of type - ii superconductors in the presence of random pinning _ , phys .
b * 76 * , 014507 ( 2007 ) j.  b.  ketterson , s.  n.  song , _ superconductivity _ , university press , cambridge ( 1998 ) l.  l.gorkov , g.  m.  eliashberg , zh .  eksp .
.  fiz . *
54 * , 612 ( 1968 ) ; a.  schmid , phys .
materie * 5 * , 302 ( 1966 ) g.  lasher , _ series solution of the ginzburg - landau equations for the abrikosov mixed state _
* 140 * , a523 ( 1965 ) d.  li , b.  rosenstein , _ lowest landau level approximation in strongly type - ii superconductors _ , phys .
b * 60 * , 9704 ( 1999 ) p.  kim , z.  yao , c.  a.  bolle , c.  m.  lieber , _ structure of flux line lattices with weak disorder at large length scales _ , phys .
b * 60 * , r12589 ( 1999 ) d.  charalambous , p.  g.  kealey , e.  m.  forgan , t.  m.  riseman , m.  w.  long , c.  goupil , r.  khasanov , d.  fort , p.  j.  king , s.  l.  lee , f.  ogrin , _ vortex motion in type - ii superconductors probed by muon spin rotation and small - angle neutron scattering _ , phys .
b * 66 * , 054506 ( 2002 ) a.  t.  fiory , _ quantum interference effects of a moving vortex lattice in al films _ , phys .
* 27 * , 501 ( 1971 ) m.  o.  sonnaillon , f.  j.  bonetto , _ a low - cost , high - performance , digital signal processor - based lock - in amplifier capable of measuring multiple frequency sweeps simultaneously _ , rev .
inst .  * 76 * , 024703 ( 2005 ) y.  ikebe , r.  shimano , m.  ikeda , t.  fukumura , m.  kawasaki , _ vortex dynamics in a nbn film studied by terahertz spectroscopy _ , phys .
b * 79 * , 174525 ( 2009 ) p .- j .
lin , p.  lipavsk , _ time - dependent ginzburg - landau theory with floating nucleation kernel : far - infrared conductivity in the abrikosov vortex lattice state of a type - ii superconductor _ , phys .
b * 80 * , 212506 ( 2009 ) a.  m.  troyanovski , j.  aarts , p.  h.  kes , _ collective and plastic vortex motion in superconductors at high flux density _ , nature , * 339*,665(1999 ) t.  h.  johansen , _ gallery of abrikosov lattices in superconductors _
, ` http://www.fys.uio.no/super/vortex ` c.  f.  lo , _ propagator of the general driven time - dependent oscillator _ , phys .  rev .
a * 47 * , 115 ( 1993 ) b.  rosenstein , private communication , unpublished .