Instruction: Write an abstract for the following acadamic paper.

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

Paper: 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 .
Abstract:

Answer: the dynamical response of the abrikosov vortex lattice in the presence of an oscillating driving field is calculated by constructing an analytical solution of the time - dependent ginzburg - landau equation . the solution is steady - state , and work done by the input signal is dissipated through vortex cores , mainly by scattering with phonons . the response is nonlinear in the input signal , and is verified for consistency within the theory . the existence of well - defined parameters to control nonlinear effects is important for any practical application in electronics , and a normalised distance from the normal - superconducting phase - transition boundary is found to be such a parameter to which the response is sensitive . favourable comparison with nbn experimental data in the optical region is made , where the effect is in the linear regime . predictions are put forward regarding the suppression of heating and also the lattice configuration at high frequency .