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: self - propelled particles convert chemical @xcite , magnetic @xcite , or radiation energy @xcite into directed motion . in this way , kinetic energy is intrinsically introduced into the corresponding system on the single particle level and length scale , in contrast to a process of macroscopic stirring or shaking from outside . far - reaching consequences on the single particle behavior as well as on their collective properties emerge . single artificially generated janus particles can start to self - propel mostly due to phoretic effects when one of the two sides is selectively heated by laser light @xcite or catalyzes a chemical reaction @xcite . in a liquid environment , such single colloidal particles obey rotational as well as translational diffusion . on time scales that are short when compared to the characteristic rotational diffusion time , self - propulsion shows up as directed motion . on significantly longer time scales , the overall diffusion coefficient is considerably enhanced @xcite . also the consequences of a possible additional particle deformability were studied extensively @xcite as it may occur , for example , for self - propelled droplets on surfaces @xcite or in a bulk fluid @xcite . another example of artificial self - propelled particles are granular hoppers that transform vertical vibration energy into directed horizontal motion @xcite .    on the biological side , the conversion of chemical energy into directed motion allows bacteria to individually search for food @xcite or to swim towards or away from illuminating light @xcite . likewise , chemical energy is used by bacteria , amoebae , or tissue cells to crawl on substrates @xcite . propulsion mechanisms are manifold and reach from the beating of a single flagellum @xcite via twisting deformations of whole filamental bacterial bodies @xcite to the synchronization of the motion of thousands of cilia @xcite . when many self - propelled particles act together , collective modes of migration emerge , like swarms that move as single entities @xcite , traveling density bands @xcite , or lanes of counterpropagating particles @xcite . even turbulent states arise that are based on the energy conversion on the single particle length scale @xcite . already without adhesive interactions and in the absence of alignment mechanisms , clusters of self - propelled particles can form in high - density systems @xcite . these clusters appear when there is no interaction that aligns the migration directions of different particles so that they can mutually block their motion . for monodisperse two - dimensional systems , hexagonal packing can arise @xcite . confining walls can support a local clustering or trapping of self - propelled particles at the boundaries @xcite . if the particle interaction is strong enough compared to their self - propulsion , then a formation of active crystals composed of self - propelled particles is to be expected . for this scenario , particle simulations demonstrated the emergence of crystal - like objects that collectively migrate as a single structure @xcite . in contrast , the approach that we recently introduced @xcite and analyze here in breadth is a complementary continuum one . it combines the description of crystalline materials by the phase field crystal model @xcite with ideas from the toner - tu model @xcite to characterize active crystalline systems of self - propelling objects @xcite . to our knowledge , there is no field - theoretical approach preceding ref .  @xcite that includes the translational ordering in active crystals of self - propelled particles . the phase field crystal model is a microscopic field approach introduced by elder et al .  to efficiently model processes in crystalline materials with particle resolution , but on diffusive time scales @xcite . in this way it can reproduce results from molecular dynamics simulations @xcite , yet in a computationally much more efficient way . motivated by dynamical density functional theory @xcite , the phase field crystal model can in principle be viewed as a microscopic theory describing a manifold of solidification and crystallization phenomena . examples are defect motion in strained crystals and multidomain structures @xcite , epitaxial growth @xcite , crystal pinning @xcite , crystal growth of binary alloys @xcite , branching in crystal growth @xcite , or even crack propagation @xcite . quite contrarily , the toner - tu model is a macroscopic hydrodynamic - like continuum description of systems of self - propelling objects @xcite . local spatial  i.e.  translational  order of the single particles is not resolved or explicitly taken into account . the focus in this model is on the orientational order of the migration directions and thus on the macroscopic emergence of collective motion . it was argued that , if locally the migration directions spontaneously order , long - ranged orientational order of the migration directions and therefore global collective motion into one single collective direction emerges even in two spatial dimensions @xcite . our present work starts from an `` active phase field crystal model '' that we recently introduced in a previous study @xcite . it combines ideas from the two outlined theories , the equilibrium phase field crystal model and the toner - tu approach namely . we extend our previous investigation by providing for example more details about inverse active lattice structures and the migration speeds of the collectively migrating textures . our main focus here is on the stability of traveling active crystals . we therefore perform a linear stability analysis of these structures and demonstrate linear stability even on large length scales that are not accessible by direct numerical simulations . furthermore , for particles self - propelling on a substrate and surrounded by a thin fluid film , we include the influence of hydrodynamic interactions : they can disturb or destroy the crystalline and orientational order . as a first step we derive in this section our active phase field crystal model on the basis of dynamical density functional theory for brownian systems @xcite . we consider a system of self - propelled particles in two spatial dimensions . examples would be amoebae or tissue cells crawling on a substrate @xcite or granular hoppers that move horizontally on a vibrating plate and can not jump over each other @xcite . such systems can directly exchange momentum with the ground . the same should also be true for low reynolds number swimmers in a thin film of fluid on a substrate @xcite . we assume that each particle migrates with an effective active drive of magnitude @xmath0 into an individual direction that is characterized by a unit vector @xmath1 . the one - particle density distribution function is therefore given by @xmath2 , with @xmath3 the spatial coordinate and @xmath4 the time . in two spatial dimensions , the orientation of @xmath1 can be parameterized by a single angle @xmath5 . thus the dynamical equation for the one - particle density follows as @xcite @xmath6 in this equation @xmath7 sets the inverse thermal energy , @xmath8 is the translational diffusion tensor , @xmath9 denotes the rotational diffusion constant , and @xmath10 is an equilibrium energy functional that will be specified below . the effect of self - propulsion enters through the last term on the right - hand side via the magnitude of the active drive @xmath0 . our next step is to derive dynamical equations for the particle density and for the polar orientational order of the self - propulsion directions . we denote the average value of the one - particle density distribution function @xmath2 as @xmath11 , i.e.   @xmath12 which is a conserved quantity in a closed system . here , @xmath13 is the spatial volume of the system and @xmath14 is the surface area of the unit sphere . the modulation of the particle density @xmath15 around @xmath11 is obtained by integrating out the orientational degrees of freedom from @xmath2 and subtracting @xmath11 , @xmath16 furthermore , we take into account the first orientational moment of @xmath2 , which we denote as @xmath17 with @xmath18 the dimensionality of the system . from this definition it becomes evident that the polarization field @xmath19 describes the local orientational order of the active driving directions @xmath1 of the self - propelled particles . similarly to the procedure in refs . @xcite , we expand the one - particle density with respect to its orientational dependence as @xmath20    to proceed we insert this expansion into eq . ( [ eqrho ] ) . from the resulting equation , the different fourier modes in terms of the angle @xmath21 are calculated . in this way we obtain dynamical equations for @xmath22 and @xmath23 similar to the ones in refs . @xcite , now , however , including the active drive . applying the approximation of constant mobility and assuming an isotropic translational diffusion @xmath24 , the resulting expressions are considerably simplified . we finally find @xmath25 \partial_{\tau}(\bar{\rho}\tilde{\mathbf{p } } ) & = & \frac{\beta \tilde{d}\bar{\rho}}{\pi}\nabla^2\frac{\delta\tilde{\mathcal{f}}}{\delta(\bar{\rho}\tilde{\mathbf{p } } ) } -\frac{\beta \tilde{d}_r \bar{\rho}}{\pi}\frac{\delta\tilde{\mathcal{f}}}{\delta(\bar{\rho}\tilde{\mathbf{p } } ) } \nonumber\\[.05 cm ] & & { } -v\nabla(\bar{\rho}\tilde{\phi}_1 ) . \label{eqpdft}\end{aligned}\ ] ]    in the following we specify the free energy functional @xmath10 . more precisely , we introduce @xmath26 as a sum of a translational part @xmath27 and an orientational part @xmath28 .    on the one hand , @xmath29\phi + \frac{u}{4}\phi^4 \right\},\ ] ] is the free energy introduced as the phase field crystal model @xcite . here , the order parameter @xmath30 is related to the field @xmath31 in eq . ( [ rhoexpand ] ) via @xmath32 and @xmath33 . the parameters @xmath34 , @xmath35 , and @xmath36 determine the magnitude of the energetic contributions , whereas @xmath37 is a measure for the temperature . minimizing this functional for an equilibrium system leads to fluid - like states at higher temperatures or high values of @xmath38 , hexagonally crystalline states at lower temperatures and intermediate values of @xmath38 , and lamellar states at lower temperatures and low values of @xmath38 @xcite . the characteristic length scale of the hexagonally crystalline and lamellar structures is determined by the magnitude of @xmath39 in eq . ( [ fpfc_phi ] ) .    on the other hand , the orientational part is given by @xmath40 ^ 2 \right\}.\ ] ] the concept behind this free energy is similar to the original approach by toner and tu @xcite . we have @xmath41 due to thermodynamic stability . if @xmath42 and @xmath43 , the system can reduce its free energy by an orientational ordering of the self - propulsion directions , given by a nonzero solution of @xmath23 . this corresponds to an intrinsic tendency towards ordered collective motion . if @xmath44 , orientational ordering of the self - propulsion directions always costs orientational free energy . consequently there is a tendency to avoid collectively ordered motion . this case shows a richer behavior than the other one and will mainly be considered below .    predominantly using the ramakrishnan - yussouff expansion @xcite , the phase field crystal model could be connected to density functional theory @xcite ; for a review see ref . expressions for the phenomenological parameters in terms of microscopic correlation functions could be derived . afterwards the procedure was extended to the case of an additional orientational order @xcite . likewise microscopic expressions for the corresponding phenomenological parameters were listed in these references , and we do not repeat them here .    to finally obtain our dynamical equations , we introduce eqs . ( [ fsum])([ftildep ] ) into eqs . ( [ eqpsi1dft ] ) and ( [ eqpdft ] ) . moreover we apply the rescaling rules @xcite @xmath45 this rescales @xmath32 to @xmath46 . all lengths from now on are measured in units of the characteristic length scale @xmath39 of the phase field crystal structures .    an additional simplification of the notation follows from the further rescaling @xmath47 \tilde{c}_4 & = & \frac{u}{4}\,c_4 , \\[.1 cm ] \tilde{d}_r & = & \tilde{d}q_0 ^ 2\,d_r.\end{aligned}\ ] ] at the end of this procedure we obtain our dynamical order parameter equations in the form @xcite @xmath48\psi_1 % \nonumber\\ % & & \qquad { } + 3\bar{\psi}\psi_1 ^ 2+\psi_1 ^ 3\,\big\ } \nonumber\\ & & { } - v_0\,\nabla\cdot\mathbf{p } , \label{eqpsi1}\\[.1 cm ] \partial_{t}\mathbf{p } & = & \nabla^2\left ( c_1\mathbf{p}+c_4\mathbf{p}^2\mathbf{p } \right ) - d_r\left ( c_1\mathbf{p}+c_4\mathbf{p}^2\mathbf{p } \right)\nonumber\\[.05 cm ] & & { } - v_0\nabla\psi_1 . \label{eqp}\end{aligned}\ ] ] the field @xmath49 is a measure for the local particle density , whereas the field @xmath50 measures the local orientational order of the active driving directions of the self - propelled particles . @xmath51 controls the temperature , @xmath52 the mean density , and @xmath53 the strength of the active drive of the individual particles . the ordering behavior of the active driving directions of the self - propelled particles is determined by @xmath54 and @xmath55 , whereas @xmath56 describes the orientational diffusion of the active drive directions . as noted above , the equilibrium phase field crystal model based on eq . ( [ fpfc_phi ] ) shows hexagonally crystalline structures in a certain range of temperature @xmath51 and mean density @xmath52 . changing the temperature or mean density , hexagonal textures can be transformed into lamellar ones . in the non - equilibrium case considered here , such a transition can also be achieved by increasing the active drive @xmath53 of the individual self - propelled particles @xcite . this scenario is depicted in the upper row of fig . [ snapshots ] . the snapshots shown in the figure were obtained by numerically solving eqs . ( [ eqpsi1 ] ) and ( [ eqp ] ) , starting from random initial conditions . we start with a low value of the active drive @xmath53 and a hexagonal structure in panel ( a ) . for the parameter values listed in the figure caption , density peaks form that are symmetrically located within `` + 1''-defects of the polarization field . the overall structure remains at rest . this changes with increasing @xmath53 , when the symmetry of the density peak positions with respect to the `` + 1''-defects of the polarization field is spontaneously broken . this spontaneous shift of the density peaks out of the centers of the `` + 1''-defects of the polarization field represents the symmetry breaking necessary for a net active propulsion to emerge . as a whole , the structure starts to migrate ( travel ) collectively  ( b ) . further increasing @xmath53 deforms the traveling hexagonal crystal to a traveling rhombic and then to a traveling quadratic crystal  ( c ) . finally it leads to a transition to traveling lamellae  ( d ) .    typically , when thinking of crystalline textures in this way , we would identify each density peak with one particle . however , also the situation of cluster crystals is included by our approach . in that case several particles occupy one lattice site @xcite . such active cluster crystalline structures were observed before in a system of deformable self - propelled particles @xcite and in a modified vicsek model confined between two parallel walls @xcite . furthermore , traveling crystalline textures have been observed previously in three - component reaction - diffusion systems @xcite . when we thus identify each density peak in the top row of fig . [ snapshots ] with a cloud of self - propelled particles , we can similarly interpret the corresponding active inverted textures in the bottom row . in each column of fig . [ snapshots ] the two snapshots can be mapped onto each other by noting a special symmetry relation of eqs . ( [ eqpsi1 ] ) and ( [ eqp ] ) . the equations remain invariant under the simultaneous transformations @xmath57 , @xmath58 , and @xmath59 . consequently the series of patterns in the bottom row of fig . [ snapshots ] follows directly as an inversion of the top row patterns . we will concentrate on the characterization of hexagonal , rhombic , and quadratic crystalline textures in the following . nevertheless the results apply equally for the inverted structures due to the symmetry property of the dynamical equations . active honeycomb textures are formed for example by flagellated marine bacteria @xcite .    to quantify the collective motion of the active crystals , we tracked the migration of each individual density peak in the sample . for this purpose , the center of each density peak was determined at fixed time intervals . on the one hand , this time interval was short enough so that identical density peaks could be identified before and after each migration step . on the other hand , it should be long enough so that the migration step is larger than the mesh size of the calculation grid . it was shown previously that  at least for the finite system sizes investigated  a collectively traveling single crystal develops from the random initial conditions via an intermediate multidomain structure @xcite . at the end of this coarsening process towards the active single crystal , all density peaks migrate collectively into the same direction . the migration speed @xmath60 of the resulting active single crystal follows from a sample average of the peak velocity magnitudes of all density peaks , @xmath61 . here @xmath62 denotes the total number of all density peaks in the sample . the velocities @xmath63 of the individual density peaks , with @xmath64 , follow from the tracking procedure described above . corresponding results are depicted in fig . [ vmvsv0 ] for characteristic values of the parameters @xmath54 and @xmath55 . as mentioned above , these parameters control the spontaneous orientational ordering of the migration directions of the individual self - propelled particles .     of the single crystals as a function of the active drive @xmath53 of the individual particles at @xmath65 . for @xmath66 and @xmath67 ( squares ) the directions of active drive of the individual particles order spontaneously and the structures travel collectively with @xmath68 . a nonzero threshold appears for @xmath69 and @xmath67 ( rhombi ) . at @xmath70 and @xmath71 ( circles ) the directions of active drive of the individual particles do not order spontaneously . then the structures are at rest ( @xmath72 ) below a finite threshold value @xmath73 ; they start to travel collectively with nonzero @xmath74 above the threshold.,width=321 ]    we plot in fig . [ vmvsv0 ] the collective migration speed @xmath60 of the single crystal as a function of the active drive @xmath53 of the individual self - propelled particles .    for @xmath70 the active driving directions of the individual particles do not order spontaneously . at low values of the active drive @xmath53 , the overall hexagonally crystalline structure remains at rest , i.e.  @xmath72 . this corresponds to the situation in fig . [ snapshots](a ) . only at a nonzero threshold value @xmath73 do the crystals start to travel collectively , see also fig .  [ snapshots](b ) . here , the nonzero value of @xmath67 only plays a minor quantitative role , and a similar curve is obtained for @xmath75 . in contrast , there is no threshold when the active driving directions of the individual particles order spontaneously , i.e.  @xmath66 and @xmath67 . the whole crystalline structure migrates at each nonzero value of the active drive @xmath53 of the individual particles . even more , the whole crystal propels collectively at the speed of the individual particles , i.e.  @xmath68 . in addition to this curve of negative @xmath54 , we show in fig . [ vmvsv0 ] an intermediate case of @xmath69 and @xmath67 . there , the nonlinear @xmath55-term leads to a small residual threshold behavior at low values of @xmath53 .    during the remaining part of this paper we will be concerned with the case of @xmath70 and @xmath75 . it shows the richer behavior due to the emergence of resting as well as traveling structures that are separated by the finite threshold value @xmath73 . introducing a hydrodynamic - like continuum theory , toner and tu studied the collective motion of flocks of self - propelling objects @xcite . their focus was mainly on the orientational ordering of the self - propulsion directions , and not on a translational , possibly crystalline order of the individual objects . the theory contains nonlinear convective terms due to self - propulsion . analyzing their model , toner and tu found that long - ranged orientational order of the self - propulsion directions can arise in two spatial dimensions @xcite , in contrast to what is known for orientational ordering in equilibrium systems @xcite . such nonlinear convective terms that could lead to long - ranged orientational order of the self - propulsion directions are absent in our approach . therefore the question arises whether a very large crystal will still migrate as a single structure into one direction , or whether it will break up into domains of different migration directions .    instead of the nonlinear convective terms , however , our approach contains a kind of elastic interaction . the phase field crystal free energy functional @xmath76 of eq . ( [ fpfc_phi ] ) acts in the form of a chemical potential in the dynamical equation ( [ eqpsi1 ] ) for the density . it enforces a locally spatially periodic arrangement of the density peaks with a characteristic wave number . deviations from this order are energetically penalized . it has been demonstrated in particle simulations @xcite that solely elastic interactions between self - propelled particles , without any explicit alignment interaction @xcite , can induce global collective motion into one common migration direction . there the elastic interactions guided the particles to find a global mode of migration , even under the influence of moderate stochastic noise . in this way , elasticity enforced global collective motion . due to the finiteness of the system sizes it is difficult , however , to make a rigid statement about what happens on very large length scales on the basis of the particle simulations . our continuum field approach offers an advantage in this respect . the explicit dynamical equations ( [ eqpsi1 ] ) and ( [ eqp ] ) contain all length scales . in principle , their investigation is not limited to a certain system size . an arbitrary large system can be analyzed . in reality , however , to get an explicit solution , we have to solve the equations on a finite numerical lattice . nevertheless , we can extract from our equations the information of whether an arbitrarily large single crystal with a single migration direction is linearly stable . our procedure will be described in the following . as mentioned before , we confine ourselves for this purpose to the case in which the self - propulsion directions do not spontaneously order , i.e.  @xmath70 ( @xmath75 ) . as a first step we linearize the dynamical equations ( [ eqpsi1 ] ) and ( [ eqp ] ) with respect to a small perturbation @xmath77 and @xmath78 of the order parameter fields , @xmath79 here , @xmath80 and @xmath81 correspond to the solutions obtained from directly solving numerically the dynamical equations ( [ eqpsi1 ] ) and ( [ eqp ] ) . these solutions describe , on the finite size of the numerical calculation grid , a single crystalline structure as illustrated for example in fig . [ snapshots ] . in our case , the converged steady - state solutions @xmath80 and @xmath81 , around which we perturb the system , are spatially inhomogeneous . this significantly complicates the analysis as will be shown below .    from the linearization we obtain @xmath82\delta\psi_1   + 6\bar{\psi}\psi_{10}\delta\psi_1 \nonumber\\ & & { } \qquad+3\psi_{10}^2\delta\psi_1 \big\ }   - v_0\,\nabla\cdot\delta\mathbf{p } , \label{eqdeltapsi1}\\[.1 cm ] \partial_{t}\delta\mathbf{p } & = & \nabla^2 c_1\delta\mathbf{p } - d_r c_1\delta\mathbf{p }   - v_0\nabla\delta\psi_1 . \label{eqdeltap}\end{aligned}\ ] ] since the coefficients of the perturbations in the second equation are constant , the system of equations can be significantly simplified in fourier space . we apply the transformations in the form @xmath83 eqs . ( [ eqdeltapsi1 ] ) and ( [ eqdeltap ] ) are transformed accordingly . then the second equation can be solved for @xmath84 , which is inserted into the first equation . in this way , the system of eqs . ( [ eqdeltapsi1 ] ) and ( [ eqdeltap ] ) is reduced to a single equation on the scalar order parameter field @xmath85 , which reads @xmath86\,\delta\psi_1(\mathbf{q},\omega ) \nonumber\\ & & { } + \frac{6\bar{\psi}}{(2\pi)^3}\int\!d^2q'\int\!d\omega'\ ; \psi_{10}(\mathbf{q}-\mathbf{q}',\omega-\omega')\delta\psi_1(\mathbf{q}',\omega ' ) \nonumber\\ & & { } + \frac{3}{(2\pi)^3}\int\!d^2q'\int\!d\omega'\ ; \psi_{10}^2(\mathbf{q}-\mathbf{q}',\omega-\omega')\delta\psi_1(\mathbf{q}',\omega ' ) \,\big\ } \nonumber\\ & & { } -\frac{q^2v_0 ^ 2}{c_1(d_r+q^2)+i\omega}\,\delta\psi_1(\mathbf{q},\omega ) . \label{stabaneq}\end{aligned}\ ] ]    as we can see , the analysis becomes nonlocal in fourier space . the reason is the spatial modulation of the traveling or resting unperturbed structures , as they are for example depicted in fig . [ snapshots ] . in principle , eq .  ( [ stabaneq ] ) represents an infinite system of equations in which all modes @xmath85 of different @xmath87 and @xmath88 are coupled . a similar situation occurs in the study of microphase - separated diblock copolymer systems ref . there , however , an equilibrium situation is investigated , which allows a different treatment of the problem . fortunately the situation can be simplified by exploiting the regular nature of the unperturbed states @xmath80 . first , the structures form a steady state in the sense that they migrate as a single object with a constant migration velocity @xmath89 . in the co - moving frame , the single crystals are at rest . second , the spatial modulation of the unperturbed structures is periodic . using these observations , eq . ( [ stabaneq ] ) can be split into discrete sets of equations in fourier space , each of them still being infinite , however . a numerical evaluation is required to handle the problem . since the unperturbed structures @xmath80 form a steady state , we can write them in the form @xmath90 this relation translates itself into fourier space in the following way : @xmath91 introducing these expressions into eq . ( [ stabaneq ] ) decouples the modes in the frequencies @xmath88 , which significantly simplifies the problem .    for the second simplification , we use the fact that the unperturbed structures are spatially periodic . they are resting or traveling crystals . in fourier space they are therefore characterized by a reciprocal lattice of a discrete set of wave vectors . translating in reciprocal space the whole structure by a reciprocal lattice vector maps the positions of the reciprocal lattice points onto themselves . in this sense , the reciprocal lattice is closed by itself under such translations . the subtractions @xmath92 in eq .  ( [ stabaneq ] ) correspond to these translations . as a consequence , the set of equations ( [ stabaneq ] ) splits into different decoupled discrete subsets of equations . within one such set , all modes are connected via addition or subtraction of the reciprocal lattice vectors of the unperturbed structure . still each set generally contains infinitely many modes , but the discretization already implies a major reduction .    in the following we investigate three different structures : a traveling hexagonal crystal , a traveling rhombic crystal , and a traveling quadratic crystal . as mentioned above , the textures are obtained in this sequence when increasing the active drive @xmath53 . on the left - hand side , [ ftstructure ] shows snapshots of the steady - state structures in real space , i.e.  @xmath93 . on the right - hand side , the corresponding power spectra @xmath94 in two - dimensional fourier space are included .     and ( a ) @xmath95 for a traveling hexagonal crystal , ( b ) @xmath96 for a traveling rhombic crystal , and ( c ) @xmath97 for a traveling quadratic crystal . only a fraction of the numerical calculation box is shown . right column : corresponding power spectra in two - dimensional fourier space . a distorted hexagonal lattice structure is obtained in fourier space in panel ( d ) from the traveling hexagonal case . panel ( e ) gives the power spectrum for the traveling rhombic crystal . finally , the square pattern in real space transforms itself into a square pattern of the two - dimensional power spectrum ( f ) . the intensity of the peaks in the power spectrum is represented by their brightness and was partially rescaled for better visualization . circles mark the directions of the reciprocal lattice vectors of smallest magnitude @xmath98 . the examples of linear stability analysis against linear dynamical instabilities that are shown in fig . [ logdeterminantdynamic ] were performed in these directions.,width=321 ]    the different textures  hexagonal , rhombic , and quadratic  can be identified by eye from the real - space plots . also the spatial arrangements of the inner intensity peaks on the right - hand side of fig . [ ftstructure ] reflect these symmetries . however , when we take into account the magnitude of the intensities , we see that the collective migration breaks the symmetry . most obviously , there is no four - fold symmetry in the intensity magnitudes of fig . [ ftstructure ] ( f ) for the quadratic texture . likewise , the single real - space density peaks in fig . [ ftstructure ] ( c ) appear a little distorted . the direction of collective migration marks the orientation of the remaining symmetry axis . we can see that the magnitude of the peaks in the power spectrum , coded by their brightness in fig . [ ftstructure ] , significantly decays with increasing wave number . as a consequence we introduce a positive threshold value @xmath99 . @xmath93 is approximated by only those modes that satisfy @xmath100 . the same procedure is applied to @xmath101 , i.e.  the fourier transform of @xmath102 .    for each investigated structure , the above procedure must be performed only once in practice . we numerically iterate eqs . ( [ eqpsi1 ] ) and ( [ eqp ] ) forward in time starting from random initial conditions . when the system has reached the steady state , @xmath103 is used as an input for the unperturbed states @xmath93 and @xmath102 . after performing the fourier transform , the values @xmath104 and @xmath101 that satisfy the above threshold criterion together with their locations @xmath87 in the reciprocal space are stored in an ordered list structure for the further analysis . we also numerically measure @xmath89 during this initial numerical iteration and use this value in our further evaluation . ( it is a challenging task to calculate the collective migration velocity @xmath89 analytically . a corresponding attempt in one spatial dimension that is restricted to the threshold vicinity is included in the appendix . ) our intention is to perform a linear stability analysis of the traveling single crystal at long wave lengths , i.e.  at small wave numbers @xmath105 . direct numerical solution of eqs . ( [ eqpsi1 ] ) and ( [ eqp ] ) showed that the structures are stable on the length scale of the size of the numerical box . but the question of what happens on larger length scales has so far remained unanswered . in other words , we are interested in the question whether the single crystal will break up into a multidomain texture for large system sizes . below , we denote the reciprocal lattice vectors of the unperturbed structure as @xmath106 , with @xmath107 a labeling index . they localize the positions of the spectral points in the right column of fig . [ ftstructure ] . the reciprocal lattice vector of smallest magnitude is called @xmath98 . we probe a square region in reciprocal space bordered at two opposing edges by @xmath108 and @xmath109 . this region covers the first brillouin zone . all wave vectors @xmath87 lying in the first brillouin zone belong to a different discrete coupled subset of equations spanned by the nonlocal terms in eq . ( [ stabaneq ] ) . we can not connect two different wave vectors in this region by adding or subtracting reciprocal lattice vectors . still , these subsets of equations are infinite . we have seen , however , that the coupling strengths @xmath110 and @xmath111 that are responsible for the nonlocal coupling in eq . ( [ stabaneq ] ) decay with increasing distance @xmath112 . therefore we introduce a cutoff distance @xmath113 and only include modes that satisfy @xmath114 in the subset of equations for each wave number @xmath87 . in this way , the discrete subsets of equations become finite and can be handled numerically .    in practice , we sample the square region between @xmath108 and @xmath109 in small wave vector steps . at each @xmath87 in this region , we build up the discrete set of equations described above . for this purpose , we numerically search for all modes @xmath115 with wave vectors @xmath116 that can be obtained by iteratively adding or subtracting the reciprocal lattice vectors that predominantly contribute to @xmath117 and @xmath118 ( those were previously stored in a list when calculating the spectra , see above ) . only modes @xmath115 that satisfy the condition @xmath114 are included . ( [ stabaneq ] ) leads for each wave vector @xmath87 and frequency @xmath88 to a linear system of equations for all modes @xmath115 that are coupled to @xmath85 . in practice , our cut - off values are chosen such that the coupling to several hundred modes is included for each wave vector @xmath87 . the coefficients of this system of equations are calculated and stored in matrix form . for practical reasons , to obtain a matrix of real coefficients , we consider the real and imaginary parts of the perturbations @xmath119 as independent variables . each equation in the system of equations ( [ stabaneq ] ) is split into its real and imaginary part , which doubles the number of equations . finally , the determinant @xmath120 of the resulting matrix is calculated by standard numerical procedures @xcite . at the end of this procedure , a vanishing determinant @xmath120 signals a linear instability of the single crystal with respect to the perturbation @xmath85 .    before continuing with the results , we add a technical remark to explain the appearance of the plots included below . due to the many modes that are coupled in each case , the number of multiplications during the process of obtaining the determinant @xmath120 is very large . in fact , @xmath120 often exceeds the size that can be processed on the computer . to circumvent this problem , we rescale the factors in the final multiplication by an adjusted constant positive number . since we are only interested in the question whether @xmath120 vanishes ( or changes sign ) this is a legal procedure . nevertheless , the absolute values of @xmath120 should not be interpreted any more . furthermore , we calculate and plot the logarithm of the determinant @xmath121 $ ] .    in a first step , we consider static instabilities in the co - moving frame . that is , we set @xmath122 ( and likewise @xmath123 ) in eq . ( [ stabaneq ] ) . thus @xmath88 does not describe an additional degree of freedom but is determined by the wave vector @xmath87 . the results corresponding to the three textures in fig . [ ftstructure ] are depicted in fig . [ logdeterminantco ] .     against instabilities of small wave - numbers in the co - moving frame . the analyzed regions in @xmath87-space contain the first brillouin zones of the textures in fig . [ ftstructure ] , respectively . surface plots show the logarithm of the determinant @xmath124 of the stability matrix in the analyzed @xmath87-regions . it has been rescaled for technical reasons in each panel , so that the absolute values should not be interpreted . below the surfaces , the shadow plots represent colored projections of the surfaces into the plane . we define @xmath125 , @xmath126 the component of @xmath87 parallel to @xmath98 , and @xmath127 its component in the perpendicular direction , where @xmath98 is the reciprocal lattice vector of smallest magnitude . a linear instability would be indicated by a vanishing @xmath124 , however , we observe nonzero values in all probed instances ( and no sign changes ) . we thus conclude that the structures are linearly stable against static instabilities in the co - moving frame.,width=321 ]    the determinant naturally vanishes at @xmath128 and at @xmath129 , which can not be shown in the logarithmic plot and leads to the white wholes in the projected shadow plots below each surface plot . we do not find any other position in the first brillouin zone where the determinant vanishes . the small dips visible in the surfaces of @xmath130 $ ] were found to decrease with increasing size of the numerical calculation box that was used to obtain @xmath117 . they were not observed to indicate a linear instability . we thus conclude that the investigated structures are linearly stable against static instabilities in the co - moving frame .    after excluding static linear instabilities in the co - moving frame , we ask for dynamical linear instabilities . for example , on larger length scales , there might appear domains migrating in different directions . this question adds a further degree of freedom to our analysis , namely the frequency @xmath88 . previously , it was fixed by the relation @xmath122 in the co - moving frame . now we can in principle choose any value independently of the selected wave vector @xmath87 . fortunately , for each perturbation @xmath85 , the coupled system of equations has the same size as above in the analysis within the co - moving frame . the additional degree of freedom @xmath88 does not lead to further couplings . we can understand this result from eqs . ( [ stabaneq ] ) , ( [ psi10co ] ) , and ( [ psi10sqco ] ) . on the one hand , eq . ( [ stabaneq ] ) implies that the modes @xmath131 are coupled via the unperturbed steady - state fields @xmath132 and @xmath133 . on the other hand , due to the steady - state form of @xmath80 and @xmath134 , eqs . ( [ psi10co ] ) and ( [ psi10sqco ] ) contain a rigid relation between @xmath87 and @xmath88 in the @xmath135-functions . in this way , each coupling wave vector @xmath136 in eq . ( [ stabaneq ] ) uniquely fixes one coupling frequency @xmath137 . the complete @xmath87-@xmath88 space is too large to be systematically tested . we therefore particularly focused on the representative directions given by the reciprocal wave vectors of each structure , and the directions perpendicular to these . in neither case did we find a vanishing determinant that would indicate a dynamical linear instability . this provides evidence that the investigated traveling active single crystals are also dynamically linearly stable . a few examples of our results are depicted in fig . [ logdeterminantdynamic ] .     against linear dynamical instabilities . the panels show stability probes performed in @xmath87-space along the direction of the reciprocal lattice vector of smallest magnitude @xmath98 . in addition , the frequency @xmath88 was varied in an interval centered around @xmath138 at each wave vector @xmath87 , with @xmath89 the collective migration velocity . we define @xmath139 , @xmath125 , and @xmath140 . surface plots show the logarithm of the determinant @xmath141 of the stability matrix in the analyzed @xmath87-@xmath142 regions . it has been rescaled for technical reasons in each panel , so that the absolute values should not be interpreted . below the surfaces , the shadow plots represent colored projections of the surfaces into the plane . our examples correspond to the directions indicated by the circles in the power spectra of fig . [ ftstructure ] for ( a ) @xmath95 , ( b ) @xmath96 , and ( c ) @xmath97 . a linear instability would be indicated by a vanishing @xmath141 ( or a sign change ) . apparently , the structures are linearly stable against linear dynamical instabilities.,width=321 ]    one might argue that there are sources of uncertainty in our analysis . for example , we obtain our initial reciprocal lattice vectors from a numerical calculation on a grid of finite size . however , if this approximation were problematic , it would rather show up as an instability of the unperturbed state . in contrast to that , we here find that the system is linearly stable despite such approximations . we are therefore confident that our analysis reflects the real behavior of the structures , and that the active single crystals are linearly stable against perturbations . in the previous section , we have seen that the active single crystals are linearly stable even at very large system sizes . the underlying stabilization mechanism is provided by the effective elastic interactions between the density peaks . it is introduced by the phase field crystal free energy functional that acts like a chemical potential for the density peaks in our dynamical equations ( [ eqpsi1 ] ) and ( [ eqp ] ) . this picture is further supported by an analysis of related particle simulations @xcite .    on the contrary the question arises , whether in real experimental systems there might be a counteracting mechanism that could destroy this long - ranged order . a natural candidate is given by hydrodynamic interactions . it turns out that hydrodynamic interactions can destabilize the system already at relatively small system sizes . this can be studied by direct numerical simulation of our dynamical equations , so that a complicated analysis as in the previous section is not necessary . in the situation that we investigate , the self - propelled particles can still directly exchange momentum with the ground . this requirement is satisfied by particles that directly migrate on a substrate , or approximately by swimmers that propel in a low reynolds number environment close to a surface with no - slip boundary conditions . however , we assume that the density peaks and lamellae can additionally interact with each other through a surrounding background fluid film . our goal in the following is a simple qualitative estimate of the effect of hydrodynamic interactions . we assume that the particles and the fluid have the same mass density , and that the system is incompressible . in this way the dynamical equation of mass continuity is trivially satisfied and does not need to be considered explicitly in the following . the location of the particles is described by the concentration or density field @xmath143 . we introduce a velocity field @xmath144 to parameterize the fluid flow that is induced by the particle motion . on the one hand , additional convective terms appear in the dynamical equations ( [ eqpsi1 ] ) and ( [ eqp ] ) that now read @xmath145 \hspace{-.7cm}\partial_{t}\mathbf{p } & = & \nabla^2\frac{\delta\mathcal{f}}{\delta\mathbf{p } } - d_r\frac{\delta\mathcal{f}}{\delta\mathbf{p } } - v_0\nabla\psi_1 -\mathbf{u}\cdot\nabla\mathbf{p } + \mathbf{\omega}\cdot\mathbf{p}. \label{eqphydro}\end{aligned}\ ] ] here the tensor @xmath146 with components @xmath147 describes rotations due to convection . we do not consider flow alignment of the polarization field in this qualitative picture @xcite . on the other hand , the additional dynamical equation for the velocity field @xmath144 in the rescaled form reads @xmath148 it is supplemented by the incompressibility condition @xmath149 , so that an explicit pressure term is not necessary . we can neglect the convective term on the right - hand side in the low reynolds number regime . when the self - propelled particles migrate on the substrate , they push the surrounding fluid and set it into motion . this effect is included by the force term @xmath150 that we introduce as @xmath151 with a sufficiently large coefficient @xmath152 . again @xmath62 is the number of density peaks in the sample . @xmath153 and @xmath63 denote , respectively , the current positions and velocities of the density peaks for @xmath154 . in practice we track the positions and motion of the centers of the density peaks as described above . at each peak center @xmath153 , the fluid flow velocity @xmath155 is then adjusted to the peak velocity @xmath63 . to mimic the presence of the particles in the fluid , the viscosity @xmath156 varies with the density @xmath49 . we set @xmath157 at positions of lowest density @xmath49 and let it linearly increase by a factor of @xmath158@xmath159 when moving to the regions of high density @xmath49 . due to the smooth density profiles in our numerical samples the effect of the contribution @xmath160 is mostly negligible . overall , the volume elements of high density @xmath49 , which represent the presence of the self - propelled particles , are considered as part of the fluid flow . this approach is similar to the `` fluid particle dynamics '' introduced by tanaka and araki @xcite . the last term @xmath161 in eq . ( [ equ ] ) represents a simple way to model the friction of the thin fluid film with its environment , with @xmath162 the friction constant . this term was used previously to investigate the complex flow behavior in quasi - two - dimensional systems of fluid films @xcite . two limits for the magnitude of @xmath163 are obvious . on the one hand , if @xmath163 is small , the friction between the fluid film and its surfaces is low and the migrating particles can easily set the surrounding fluid into motion . on the other hand , if @xmath163 is very large , we see from eq . ( [ equ ] ) that the fluid remains practically at rest with @xmath164 . then the convective terms in eqs . ( [ eqpsi1hydro ] ) and ( [ eqphydro ] ) are negligible and we reobtain the previous eqs . ( [ eqpsi1 ] ) and ( [ eqp ] ) . thus hydrodynamic interactions are not important in the regime of large @xmath163 . we illustrate the flow behavior at different values of the friction parameter @xmath163 in fig . [ fluidflow ] .        for high @xmath163 , corresponding to strong fluid friction , most parts of the fluid film can not be set into motion as illustrated in panels ( a ) and ( d ) . the fluid flow remains very localized around the density peaks and does not extend between different peaks . consequently the hydrodynamic interactions between the density peaks are low . around the density peaks some `` backflow '' of the fluid oppositely to the peak migration directions occurs . this is stressed in a rescaled version of panel ( d ) , included as panel ( e ) ; a `` zoom '' into the region around one single density peak is provided for illustration . with decreasing fluid friction @xmath163 , the `` backflow '' vanishes and the fluid flow more and more extends between the density peaks , as depicted in panels ( b ) and ( f ) . at very low friction parameters @xmath163 , the whole fluid can be set into motion nearly homogeneously , as shown in panels ( c ) and ( g ) . since the fluid flow here extends between the density peaks , hydrodynamic interactions are strong in this case .    as for fig . [ vmvsv0 ] , we measured again the collective migration speed @xmath60 in the form of the sample - averaged peak velocity magnitude , @xmath61 , with @xmath62 the number of the density peaks and @xmath63 the velocity of the single density peaks , @xmath154 . a result of @xmath60 as a function of the active drive @xmath53 of the individual self - propelled particles is plotted in fig . [ alphadep ] for two different     of the single crystals as a function of the active drive @xmath53 of the individual particles in the presence of hydrodynamic interactions at @xmath165 . for low hydrodynamic interactions , corresponding to a large fluid friction parameter @xmath163 , the collective migration speed is lower . for strong hydrodynamic interactions , corresponding to a low fluid friction parameter @xmath163 , the collective migration speed is enhanced . in the latter case the particles can push each other by setting the fluid between them into motion . the critical value of the self - propulsion velocity @xmath73 remains unchanged by the presence of the hydrodynamic interactions . , width=321 ]    values of the fluid friction parameter @xmath163 . in this case @xmath70 , i.e.  there is no alignment mechanism that would lead to a spontaneous orientational ordering of the self - propulsion directions .    as a first result , we found that , within our numerical resolution , the critical self - propulsion velocity @xmath73 , at which collective motion of the active crystal sets in , remains unchanged by the hydrodynamic interactions . second , the hydrodynamic interactions speed up the collective motion of the active single crystal . we can see this from the two data curves plotted in fig . [ alphadep ] . at high values of @xmath163 the fluid friction is large , and the self - propelling particles can not set the surrounding fluid into motion . however , at low fluid friction for low values of @xmath163 , the particles can set the whole surrounding fluid film into motion , as was also demonstrated in fig . [ fluidflow ] . in this way , one density peak can push the preceding peak via the fluid between them . the density peaks support each other in migration via hydrodynamic interactions . at fixed active drive @xmath53 of the individual self - propelled particles , the collective migration speed @xmath60 increases with decreasing fluid friction @xmath163 and increasing hydrodynamic interactions . however , the hydrodynamic interactions disturb the order of the active single crystals . this is illustrated in fig . [ destabilizedcrystal ] .    . ( a ) negligible hydrodynamic interactions at strong fluid friction ( @xmath166 ) allow the emergence of a perfect traveling active single crystal . ( b ) strong hydrodynamic interactions at low fluid friction ( @xmath167 ) destabilize this texture.,width=321 ]    the figure includes two identical systems that only differ by the fluid friction parameter @xmath163 . at high @xmath163 and therefore little hydrodynamic interaction ( a ) , the sample shows a perfect active single crystal that migrates collectively into one direction . on the contrary , the crystal is broken up into different domains at low @xmath163 corresponding to strong hydrodynamic interactions ( b ) . also a partial transition to lamellar textures is induced by the hydrodynamic interactions . in the case without hydrodynamic interactions this transition would typically occur only at higher magnitudes of the active drive @xmath53 @xcite .    in summary , we can say that  at least within the qualitative picture presented  hydrodynamic interactions have a destabilizing effect on the active crystalline structures considered here . in this paper we focused on the stability of traveling active single crystals . an active phase field crystal model was introduced for this purpose . we outlined its derivation via dynamical density functional theory from a microscopic consideration of self - propelled particles that feature an active drive . besides resting hexagonal crystals , traveling hexagonal , rhombic , and quadratic crystals , as well as traveling lamellar textures were obtained from this description . also the inverted structures like resting and traveling honeycomb textures , as well as inverted traveling rhombic and quadratic lattices were observed . a linear stability analysis did not indicate a linear instability of the collectively migrating single crystals . in particular , long - wave - length instabilities beyond the size of the numerical calculation grids did not show up , indicating linear stability also for large system sizes . we found , however , that hydrodynamic interactions can destabilize the single - crystalline textures . the linear stability analysis was significantly complicated by the spatial modulation of the unperturbed steady state . as a consequence of this spatial modulation , our equations became nonlocal in fourier space . we could reduce the problem to manageable size by exploiting the periodicity of the unperturbed state and its steady - state migration . this made the linear stability analysis possible . the problem of a spatial modulation of the unperturbed state occurred already in a previous equilibrium stability study of micro - phase separated diblock copolymer melts @xcite . there , however , the equilibrium conditions allowed a different strategy to perform the analysis . our results on the stability of the traveling crystals are important when we think of the design of new active materials . it is expected that , for example , features of traveling crystals like phononic properties are different from their passive counterparts . we have seen that the structure of the migrating crystals , i.e.  their hexagonal , rhombic , or quadratic form , can be tuned by the activity of the constituent self - propelled particles . it will be interesting to test whether other features of active crystals , such as their phononic behavior , can be similarly tuned by the properties of the individual active components . we shall investigate this interesting question in the future . the authors thank the deutsche forschungsgemeinschaft for support of this work through the german - japanese project lo 418/15 . in this appendix we briefly analytically address the onset of collective migration . our goal is to derive an expression for the collective migration speed @xmath60 . this is a challenging task . we confine ourselves to one spatial dimension that we denote as the @xmath168-direction . furthermore , we consider the case of @xmath70 and @xmath75 . still our results remain restricted to the very close vicinity of the threshold for the onset of collective motion .    under these circumstances , and using the fact that the structures collectively migrate in a steady - state solution with speed @xmath60 , see eq . ( [ steadystate ] ) , we can rewrite eqs . ( [ eqpsi1 ] ) and ( [ eqp ] ) as @xmath169\psi_1   \label{apppsi1 } \nonumber\\ & & { } \qquad+ 3\bar{\psi}\psi_1 ^ 2+\psi_1 ^ 3\,\big\ }   - v_0\,\partial_z{p } , \\[.1 cm ] -v_m\partial_z{p } & = & c_1\!\left ( \partial_z^2-d_r \right)\!p - v_0\partial_z\psi_1 . \label{appp}\end{aligned}\ ] ]    introducing the abbreviations @xmath170 and @xmath171 , we can write the formal solution of eq . ( [ appp ] ) as @xmath172    \partial_z\psi_1   \nonumber\\ & = & { } -\bar{v}_0\partial_z\left [ \gamma+\bar{v}_m\gamma^2 \partial_z    + \bar{v}_m^2\gamma^3\partial_z^2+\bar{v}_m^3\gamma^4\partial_z^3+ ... \right]\psi_1 ,   \nonumber\\&&\label{solp}\end{aligned}\ ] ] where we have defined @xmath173 via @xmath174 with the abbreviation @xmath175 when the operator @xmath173 is applied to an arbitrary function @xmath176 . we have further used the relation @xmath177\psi(z ' ) \nonumber\\ & = & \partial_z \int dz'\,\gamma(z - z')\psi(z')\nonumber\\ & & \end{aligned}\ ] ] in the last step of eq . ( [ solp ] ) . ( [ solp ] ) is now inserted into eq . ( [ apppsi1 ] ) to obtain a closed equation for @xmath49 : @xmath169\psi_1   \nonumber\\ & & { } \qquad+ 3\bar{\psi}\psi_1 ^ 2+\psi_1 ^ 3\,\big\ }     \nonumber\\ & & { } + v_0\bar{v}_0\partial_z^2\big [ \gamma+\bar{v}_m\gamma^2 \partial_z \nonumber\\ & & { } \qquad+\bar{v}_m^2\gamma^3\partial_z^2    + \bar{v}_m^3\gamma^4\partial_z^3+ ... \big]\psi_1 . \nonumber\\\end{aligned}\ ] ] we rewrite this expression in the form @xmath178\psi_1 \nonumber\\ & & { } \qquad+ 3\bar{\psi}\psi_1 ^ 2+\psi_1 ^ 3\,\big\ }     \nonumber\\ & & { } + v_0\bar{v}_0\big [ \gamma+\bar{v}_m\gamma^2 \partial_z \nonumber\\ & & { } \qquad+\bar{v}_m^2\gamma^3\partial_z^2    + \bar{v}_m^3\gamma^4\partial_z^3+ ... \big]\psi_1 ,   \nonumber\\ \label{gintroduced}\end{aligned}\ ] ] where we have defined @xmath179 with an abbreviation analogous to eq .  ( [ abbrgamma ] ) .    next , we multiply eq . ( [ gintroduced ] ) by @xmath180 and carry out the integral over @xmath168 : @xmath181\psi_1.\qquad \nonumber\\ & & \label{inteq}\end{aligned}\ ] ] here the integrals are taken over the whole domain of size @xmath182 . the periodic boundary conditions at the domain boundaries imply @xmath183 for any integrable functional @xmath184 . our last step consists of expanding the operator @xmath173 on the right - hand side of eq . ( [ inteq ] ) into powers of @xmath185 . on the one hand , we have @xmath186 and on the other hand , @xmath187 with integers @xmath188 and @xmath189 . using these relations , we finally obtain from eq . ( [ inteq ] ) an expression for the collective migration speed at the onset of collective motion : @xmath190 where @xmath191g(z - z')\partial_{z'}\psi_1(z ' ) , \nonumber\\ & & \label{b0 } \\ b_1 & = & \frac{v_0 ^ 2}{c_1 ^ 2}\int dz\,\int dz'\,\int dz''\,\left[\partial_z\psi_1(z)\right ] \nonumber\\    & & { } \qquad\times\gamma(z - z')\gamma(z'-z'')\partial_{z''}\psi_1(z '' ) , \nonumber\\ & & \label{b1 } \\ b_3 & = & \frac{v_0 ^ 2}{c_1 ^ 4}\int dz\,\int dz'\,\int dz''\,\int dz'''\,\int dz''''\ , \nonumber\\    & & { } \qquad\times\left[\partial_z^2\psi_1(z)\right]\gamma(z - z')\gamma(z'-z '' ) \nonumber\\    & & { } \qquad\times\gamma(z''-z''')\gamma(z'''-z'''')\partial_{z''''}^2\psi_1(z '' '' ) . \nonumber\\ & & \label{b3}\end{aligned}\ ] ] from the symmetry of the integrands in these expressions , and since @xmath192 and @xmath173 are positive operators , we expect all coefficients @xmath193 , @xmath194 , and @xmath195 to be positive . this implies a supercritical bifurcation of the collective migration speed at the onset of collective motion , which is in line with our simulations in two spatial dimensions . however , since @xmath196 generally depends on @xmath60 , the coefficient @xmath195 might be renormalized from the @xmath60-dependence of @xmath193 and @xmath194 .    directly at the threshold for the onset of collective motion , we may assume a harmonic functional dependence of the form @xmath197 . when we calculate the coefficients @xmath193 , @xmath194 , and @xmath195 from eqs . ( [ b0])([b3 ] ) , we obtain from eq . ( [ eqvm ] ) for the collective migration speed directly above onset : @xmath198^{\frac{1}{2}}.\ ] ] however , higher harmonics in the expression for @xmath196 and its @xmath60-dependence will affect this result . 97ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop _ _  ( , ) @noop _ _  ( , ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( )
Abstract:

Answer: a recently introduced `` active phase field crystal model '' describes the formation of ordered resting and traveling crystals in systems of self - propelled particles . increasing the active drive , a resting crystal can be forced to perform collectively ordered migration as a single traveling object . we demonstrate here that these ordered migrating structures are linearly stable . in other words , during migration , the single crystalline texture together with the globally ordered collective motion is preserved even on large length scales . furthermore , we consider self - propelled particles on a substrate that are surrounded by a thin fluid film . we find that in this case the resulting hydrodynamic interactions can destabilize the order