Instruction: Write an abstract for the following acadamic paper.

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

Paper: the experimental realization of quantum condensation of quasiparticles with finite lifetimes has opened the possibility to study situations where the steady states of the condensate are not just controlled by energetics , but involve also consideration of pumping and decay , leading to steady states with quasiparticle currents . the exploration of the possible behavior exhibited in these systems is driven particularly by recent experimental progress in microcavity polaritons@xcite and magnons@xcite ( as well as more exotic realizations in helium@xcite ) . microcavity polaritons are the quasiparticles that result from strong coupling between photons confined in planar semiconductor microcavities and excitons confined in quantum wells . the photons are confined by bragg reflectors to a two - dimensional cavity , and for small in - plane momentum have an effectively quadratic dispersion , with a mass of the order of @xmath0  @xmath1 of the free electron mass . from strong coupling between these photon states and exciton states in quantum wells , new normal modes arise : polaritons . these are hybrid light  matter particles , which have both a light effective mass , as well as sufficient interactions to allow a quasi - equilibrium particle distribution to be established . however , the distribution is only quasi - equilibrium , since the photon component is not perfectly confined , so polaritons have a finite lifetime , and to maintain a steady state population , continuous pumping is required . even when the energy distribution of polaritons looks reasonably thermal , this does not mean effects of finite lifetime may be neglected , as one may have thermalized energy distributions on top of steady state flows driven by the pumping and decay . experimentally , it is observed that above a threshold pumping strength , an accumulation of low energy polaritons is accompanied by : a significant increase in temporal coherence , spatial coherence that extends over the entire cloud of polaritons@xcite , the appearance of quantized vortices@xcite , suggestions of changes to the excitation spectrum@xcite , and a robustness of polariton propagation to disorder@xcite . of these , the observation of spontaneous vortices@xcite in polariton condensates provides a strong hint that the steady states of the system are affected by the flow of particles . for a more general introduction to microcavity polaritons , their condensation , and the physical systems studied there exist several reviews or books , e.g. , refs .  . the aim of this article is to investigate the rich varieties of behavior that come from combining spatial profiles driven by particle flow with the dynamics due to the spin degree of freedom that microcavity polaritons also possess . the effects of this spin degree of freedom have been previously considered in equilibrium : due to the weak attractive interaction between opposite spin polarizations , a condensate is expected to be linearly polarized@xcite , i.e. , to have an equal density of left- and right - circularly polarized components . applying a magnetic field can then lead to a phase transition , as the density of one circular polarization increases , and the other decreases , leading through an elliptically polarized phase to a single circularly polarized condensate@xcite . in the absence of further symmetry breaking , the linear or elliptically polarized phase has two gapless modes , corresponding to two broken symmetries : the overall polariton phase , and the direction of linear polarization ( i.e. , orientation of elliptical polarization ) . the transition to a circular polarized state is thus a phase transition , as such a phase has only a single broken symmetry , and a single gapless mode .    assuming no other symmetry breaking terms , the linearly polarized condensate can have independent vortices of the left- and right - circularly polarized components@xcite ; such polarization - dependent phase winding has recently been observed@xcite . these two vortices are not completely independent . even in the absence of any further symmetry breaking , there is a short - range density - density interaction between vortices of opposite polarizations@xcite , however such a term does not depend on the circulations of the two vortices . if one also takes account of te - tm splitting of the electromagnetic modes@xcite , this leads to a term that splits linear polarization states according to whether the polarization is parallel or perpendicular to the polariton momentum ; as such , this provides a circulation dependent interaction between vortices of the two polarizations@xcite . in addition , real materials are expected to possess a small splitting between linear polarization states due to asymmetry of the quantum - well interfaces in non - centro - symmetric crystals@xcite ; such a splitting can also be controlled by applying electrostatic fields@xcite or stress@xcite . such terms will again induce interactions between vortices of the left- and right - circular polarization , and if strong enough , will lead to a pinning of the polarization of a polariton condensate , as has been observed in experiment @xcite . considering the splitting of linear polarizations as a phase - locking term between the two circular polarization components , and combining this with a magnetic field that favors one or other circular polarization , one has a josephson problem@xcite , where the energy favoring linear polarization is a josephson coupling , and the magnetic field leads to a bias between the two fields . @xcite have considered the interplay between these spinor josephson oscillations , and josephson coupling between two different spatial modes of a trapped polariton condensate , showing that complex behavior can arise in this four - mode system without pumping and decay . a similar problem has also been studied in the context of spinor condensates of cold atoms in double well potentials@xcite . if one does include pumping and decay then even the two - mode system ( i.e. , just the dynamics of the spin component ) can show a rich variety of behavior @xcite . due to dissipation , the dynamics settles on an attractor , which may either correspond to a phase - locked ( i.e. , linearly polarized ) state , or may correspond to a limit cycle , in which the phase difference between the two polarization components continually increases , with the cyclic nature arising from @xmath2 periodicity in the phase difference . there also exist parameter ranges where there is bistability , so that which attractor the system settles on depends on the initial conditions . in this article , we will both review this two - mode dynamics , and then extend to the case of including many spatial modes in addition to the spin dynamics .    to describe multiple spatial modes in a trapped , pumped , decaying condensate , a mean - field approach to this problem leads to a complex gross  pitaevskii equation@xcite . it is also possible to include fluctuations beyond mean - field theory within such a formalism by quantum stochastic approaches@xcite . as well as describing the spatial structure of nonequilibrium condensates , the complex gross  pitaevskii equation has also been applied in a wide range of areas , see for example refs .  . one application which is very closely connected to microcavity polaritons is the dynamics of lasers propagating through strongly nonlinear materials , where combinations of cubic and quintic nonlinearities can produce a state with a preferred density , sometimes referred to as `` liquid light''@xcite . the results in this article suggest that consideration of polarization dynamics in such materials may also provide interesting results .    by considering the spinor complex gross  pitaevskii equation to describe the spinor condensate in an harmonic trap , several new classes of attractors are seen to occur in addition to those present for the two - mode system . there are limit cycles describing small oscillations of the phase difference between the two components , and limit cycles with @xmath3 periodicity of the phase difference . in addition , the different dynamical attractors correspond to different spatial profiles , in which the presence or absence of vortex lattices can be influenced by the applied magnetic fields .    as well as describing the steady state profiles , the complex gross  pitaevskii equation can describe the normal modes for small fluctuations about such steady states . in the absence of a spatial trap , one notable consequence of including pumping and decay is that the long - wavelength modes of a dissipative condensate are modified , to become diffusive@xcite ; similar results are seen both for incoherent pumped condensates , and for optical parametric oscillation , where scattering of pumped polaritons into signal and idler states is modified by bosonic enhancement due to the population of the signal@xcite . from this standpoint , another purpose of this article is to address how the combination of pumping , decay , and symmetry - breaking terms modify the spectrum of the spinor condensate . while a diffusive mode for the overall polariton density and phase always exists , the long - wavelength modes for spin modes can show a range of behaviors : diffusive , linearly dispersing , or gapped , dependent on the balance of decay and symmetry - breaking terms . the remainder of this article is arranged as follows . in section ii we introduce the model of spinor condensates as a system of coupled gross  pitaevskii equations with pumping and decay terms . by neglecting the spatial variations of polarized components we discuss the bifurcation diagrams for a two - mode system in section iii . the normal modes of the homogeneous system are obtained in section iv . in section v we study the stability of cross - polarized vortices . the detailed study of the dynamics of the full trapped system is performed in section vi . in section vii we discuss how various regimes can be detected in experiments . the conclusions in section viii summarize our findings . the model we consider in this article consists of the complex gross - pitaevksii equation ( cgpe ) as described in ref .  , taking into account the two possible polarizations of polaritons , written in the basis of left- and right - circular polarized states , denoted by @xmath4 . in addition to the terms discussed in ref .  , three new parameters arise from considering the spin degree of freedom in an applied magnetic field . firstly , in addition to an interaction with the total polariton density @xmath5 . secondly , there is a magnetic field term @xmath6 , which may naturally arise@xcite due to asymmetry at the quantum - well interfaces@xcite , or may be induced by electric fields@xcite or stress@xcite . put together , these yield the coupled cgpe : @xmath7    \psi_l    +     j_1 \psi_r,\end{gathered}\ ] ] with the analogous equation for @xmath8 following by replacing @xmath9 and @xmath10 . in this article , we will either neglect @xmath11 ( in sections  [ sec : review - two - mode ] , [ sec : norm - modes - extend ] ) , or consider a harmonic trap , @xmath12 . in either case , we will rescale energies and lengths in terms of the harmonic - oscillator energy @xmath13 , and length @xmath14 such that @xmath15 . by additionally rescaling the density such that @xmath16 the cgpe can be written as : @xmath17    \psi_l    +     j \psi_r\end{gathered}\ ] ] with the definitions @xmath18 . if considering a harmonic trap , then @xmath19 , otherwise we will take @xmath20 . the dimensionless parameter @xmath21 describes the extent of spin anisotropy of the interaction ; using the estimate in ref . we take @xmath22 . from estimates of the polarization splitting in refs . , we take @xmath23 . using the estimates of trap frequency given in ref . , we will take @xmath24 , and we may write the dimensionless splitting @xmath25 ; in the numerical results , we will consider a variety of values of @xmath26 around this value . the spatially extended , harmonically trapped system is unstable unless pumping is restricted to a finite spot@xcite . for our studies of this system , we therefore consider a circular pumping spot with a cutoff radius @xmath27 : @xmath28 . for constant pumping strength , @xmath27 determines whether the system forms vortices spontaneously@xcite . in order to provide a basis from which to understand the behavior observed in the trapped system , it is first useful to review the results that occur in the `` two - mode model '' , where spatial variation of each component is neglected , so the coupled complex gross  pitaevskii equation reduces to two equations for the variables @xmath29 . this model was studied by @xcite , where the conditions for the existence of a steady state ( i.e. , synchronized ) solution were found . as @xmath30 increased , @xcite found that a synchronized solution could persist for some range of @xmath30 , with an increasing phase and density difference between the two modes , until at some critical @xmath30 the steady state vanished . near the critical @xmath30 , the existence of bistability was noted . this section will both summarize these previous results , as well as analyze the behavior of the desynchronized solution , showing that the dynamics is analogous to that of the damped driven pendulum , or current - biased josephson junction@xcite .    to discuss the dynamics , it is convenient to reparametrize using : @xmath31 and write coupled equations for @xmath32 . ( the global phase @xmath33 does not affect the dynamics . ) in visualizing the dynamics , one may consider a bloch vector , defined by @xmath34 , in terms of the above variables . since there is pumping and decay , the length of the bloch vector is not conserved , hence there is a dynamical equation for @xmath35 . the coupled equations have the form : @xmath36 writing the equations in this form firstly allows one to understand the basic effects of pumping and decay on the dynamics , and , secondly , will be shown below to reduce to the damped driven pendulum in a relevant limit .    for a steady state , it is clear that the condition @xmath37 defines a bloch surface , @xmath38 . one may then study the dynamics of points that start on this surface to characterize the nature of the transition as @xmath30 increases . as shown in fig . [ fig : attraction - basin ] , each starting point @xmath39 is attracted either to a fixed point , or to a limit cycle . the fixed point is the synchronized solution ; the limit cycle is the desynchronized solution , in which the phase difference @xmath40 between the two spin components is changing  this behavior is a cycle since the dynamics are periodic under @xmath41 . as @xmath30 increases , there is a lower critical @xmath42 below which the all initial conditions on the bloch surface flow to the synchronized solution , and an upper critical @xmath43 which corresponds to the previously found@xcite point beyond which synchronized solutions can no longer exist . , as indicated , and shows a colormap according to the final state found by starting from the given initial conditions @xmath44 . points that flow to the fixed point are colored gray , others are white . the fixed point is marked by a red cross , and the limit cycle ( when it exists ) by the magenta line.,width=307 ]    the behavior seen in fig . [ fig : attraction - basin ] can be understood by noticing that the choice of parameters used puts one in the `` josephson regime '' of the josephson equations@xcite , i.e. , @xmath45 , and that in addition the typical dynamics obey @xmath46 . in this case , eq . ( [ eq:6 ] ) simply reduces to @xmath47 . then , by eliminating @xmath48 from eq . ( [ eq:4],[eq:5 ] ) , one finds : @xmath49 this equation describes a damped driven pendulum , or alternatively a current - biased josephson junction . , due to our choice of sign of @xmath26 , `` gravity '' for the pendulum acts to drive @xmath50 ) . the behavior of eq . ( [ eq:7 ] ) is well known ( see , e.g. , ref .  ) ; for large damping @xmath51 there is a simple transition at @xmath52 between fixed points for @xmath53 and limit cycles above . for smaller damping ( explicitly , for for @xmath54 ) , there is a range of bistability , where both limit cycles and fixed points may be found@xcite . the inset of fig . [ fig : ft - two - mode ] shows the bifurcation diagram determined from eq . ( [ eq:4][eq:6 ] ) , and for comparison , the value @xmath43 and the critical damping that result from the approximations leading to eq . ( [ eq:7 ] ) .    as well as explaining the classes of behavior seen , eq . ( [ eq:7 ] ) may be used to explain the critical behavior near @xmath42 . for @xmath51 sufficiently small that bistability exists , then near @xmath42 , the period of the limit cycle grows like @xmath55|$ ] . in the context of the current problem , this period relates to the average chemical potential difference between the two components , @xmath56 . for large @xmath30 , one eventually reaches @xmath57 . the evolution of @xmath58 is illustrated in fig . [ fig : ft - two - mode ] , which plots the spectral weight @xmath59 , choosing initial conditions so the limit cycle is obtained for all @xmath60 .    , showing bifurcation of frequencies as @xmath30 is increased , and illustrating critical behavior @xmath61 . initial conditions for each @xmath30 are chosen so that a limit cycle ( i.e. , desynchronized solution ) is found for all @xmath62 . inset : regions of stability of the synchronized and desynchronized solutions found from eq . ( [ eq:4][eq:6 ] ) . the dotted line marks the approximate @xmath43 appropriate in the josephson regime , eq . ( [ eq:7 ] ) , and the @xmath63 marks the known bifurcation point@xcite in that regime.,width=307 ] before considering the interplay of spin dynamics with the spatial profile of a trapped condensate , one may first consider a simpler problem involving the interplay of spin and spatial dynamics  the finite momentum normal modes of the pumped , decaying spinor condensate . the normal modes of the spinor condensate without pumping and decay were discussed in @xcite . in the equilibrium case , for small @xmath64 , there is an elliptical condensate ( i.e. , a finite density of both spin components ) and there are two gapless linear modes , describing excitations of the global phase , and the relative phase of the two components . when @xmath64 becomes large enough to cause a transition to a circularly polarized state , only a single gapless mode survives , describing phase modes of the condensed component . on the other hand , the presence of pumping and decay is known to replace the linear dispersion of the phase modes with a diffusive behavior at small momentum@xcite . the aim of this section is to see how these two effects are combined . the result is that while the global phase mode remains diffusive , the real part of the relative phase mode can be either gapped , linearly dispersing , or diffusive according to the value of @xmath30 chosen . in the following , we will first present numerical results for the normal - mode frequencies , and then discuss how these can be straightforwardly interpreted in the same josephson regime as discussed above .    for the purpose of numerically calculating the normal modes , it is simplest to consider the bogoliubov parametrization of fluctuations at some wavevector @xmath65 ( allowing for decay ) as : @xmath66 substituting this into eq . ( [ eq:2 ] ) , and linearizing in @xmath67 yields the secular equation @xmath68 = 0 $ ] , where in the basis @xmath69 , the matrix @xmath70 is @xmath71 noting that for plane waves , the kinetic energy term @xmath72 , the matrix elements are : @xmath73 and similarly with @xmath74 . the normal modes calculated this way are shown in figs . [ fig : uniform - stability ] ,  [ fig : uniform - spectrum ] . figure  [ fig : uniform - stability ] shows the modes at @xmath75 ( lower panels ) as well as the value of @xmath30 , and the densities @xmath76 ( upper panels ) corresponding to a solution with a given value of @xmath40 . [ this is to make use of the fact that the steady state of eqs . ( [ eq:4][eq:6 ] ) can be found explicitly for @xmath77 as a function of @xmath40 ( see ref .  ) . ]    , between the two circular polarization components . top left : densities of components in steady state solution . top right : detuning corresponding to given solution . bottom : real ( left ) and imaginary ( right ) parts of frequency for small fluctuations about the steady state . the gray shaded region is unstable to small fluctuations . the crosses mark the conditions used for the finite @xmath65 spectra in fig . [ fig : uniform - stability ] . plotted for @xmath78.,width=307 ]    in fig . [ fig : uniform - stability ] , the range of @xmath40 is restricted to @xmath79 $ ] , since the range @xmath80 $ ] is equivalent to this , after swapping @xmath74 and @xmath81 . within this range , only values @xmath82 correspond to stable solutions . as only the modes at @xmath75 are shown , there is always a zero mode corresponding to global phase rotations . the other three modes divide into an overdamped mode ( largest imaginary part ) , and a pair of modes which can be either overdamped or underdamped , as seen by the bifurcation of the real part at @xmath83 . three values of @xmath40 ( corresponding to three different applied detunings @xmath30 ) are chosen to illustrate the overdamped , critically damped , and underdamped cases , and the normal modes with finite @xmath65 are shown for these points in fig . [ fig : uniform - spectrum ] . in the overdamped case ( for @xmath40 near @xmath84 ) , both the relative - phase and global - phase modes are diffusive at small @xmath65 . when underdamped , the relative - phase mode always has a nonzero real frequency . when critically damped , the real part of the relative - phase mode has a linear dispersion . we will next discuss how this behavior can be understood in the josephson regime .    . in all cases the total density modes are diffusive , while the spin mode changes from diffusive ( top row , @xmath40 nearest @xmath84 ) , to linearly dispersive ( middle row ) to underdamped spin oscillations ( bottom row , @xmath40 nearest @xmath85).,width=307 ]    to understand the modes qualitatively , it is clearest to work in the basis of @xmath86 . extending eqs . ( [ eq:4][eq:6 ] ) to include spatial gradients via a madelung transformation , one finds the equations : @xmath87 where we have introduced the shorthand @xmath88 , with @xmath89 . when considering fluctuations about the uniform state , the gradient terms simplify , as expressions such as @xmath90 are necessarily second order in the fluctuations , and so may be neglected . we will make the same approximations as led to eq . ( [ eq:7 ] ) ( i.e. , assume @xmath91 ) , and use the steady - state solution in this limit ( @xmath92 ) to simplify the equations for fluctuations . if one writes the fluctuations as @xmath93 then since these are real variables , the normal modes have the time dependence @xmath94 , giving a secular equation : @xmath95 = 0,\ ] ] where @xmath96 indicates an expansion in powers of momentum , to understand the small @xmath65 dispersion . the first two terms are : @xmath97    the nature of the @xmath75 normal modes is immediately clear : @xmath33 has no restoring force , and so has a zero frequency oscillation ; @xmath98 has a damping rate @xmath99 , and so describes a decaying mode ; @xmath48 and @xmath40 are mixed , and have modes with frequencies given by : @xmath100 noting that @xmath101 is negative , and that the prefactor is the same as that of @xmath102 in eq . ( [ eq:7 ] ) , one may write : @xmath103 where @xmath104 is a `` plasma frequency '' describing the restoring force as a function of angle . the normal modes are @xmath105 . the transition between under- and overdamping occurs because for @xmath106 the restoring force is large so @xmath107 , while as @xmath108 , the restoring force vanishes , and so there is an intermediate value of @xmath40 ( and hence @xmath30 ) at which @xmath109 , describing the critical damping . note that whatever the choice of parameters , an overdamped regime will always exist since @xmath110 must always tend to zero at @xmath111 . however if @xmath51 is sufficiently large the underdamped regime may vanish . the approximate expressions for the eigenvalues @xmath112 match the general form seen in the bottom panel of fig . [ fig : uniform - stability ] , it is apparent that the full problem has some additional variation of damping rate with @xmath40 , not captured by the approximations used in the above expressions . one may now explain the linear dispersion at this critical damping as arising from degenerate perturbation theory , allowing a @xmath113 perturbation to give a @xmath114 splitting . setting @xmath109 , and writing @xmath115 , then expanding eq . ( [ eq:11 ] ) to leading order in @xmath116 gives : @xmath117,\ ] ] hence describing modes @xmath118 . note that although these modes have a linear dispersion of the real part as a function of @xmath65 , they have a lifetime that remains finite as @xmath119 . as such , this may provide an example where linear dispersion of a given mode in a condensed system need not imply superfluidity of the associated density . further work is required to determine the current - current response function in this system , which will determine whether there is a difference of transverse and longitudinal response for spin currents , however superfluidity of the spin current is not expected here , since the spin orientation is locked by the josephson term @xmath26 . the linear dispersion arises only when the damping matches frequency of this phase - locking term , producing critical damping . when an harmonic trap is introduced , the existence of steady state currents can destabilize the thomas  fermi like density profile , and lead to spontaneous rotating vortex lattices@xcite . before addressing how this instability interacts with the polarization degree of freedom , we first discuss the simpler question of the dynamics and stability of individual vortices of opposite polarization . one should note first that a single vortex in a pumped decaying condensate already has a more complicated structure than the vortex without pumping and decay  the vortex becomes a spiral vortex , with both radial and azimuthal currents ; see appendix  [ sec : spir - vort - single ] for further discussion .    considering a spinor condensate system with vortices in both polarization components , there is an interplay between the weak attractive coupling that invites the vortices to have their cores aligned irrespective of their circulation , the josephson coupling that discourages the alignment of vortices of opposite circulation and the detuning that stabilizes such alignment . to illustrate these possibilities we consider several examples . the notation @xmath120 refers to condensates with vortices of topological charge @xmath121 ( @xmath122 ) in left ( right ) condensate . from eq . ( [ eq:2 ] ) , with @xmath123 , the following set of stability scenarios are found :    @xmath124 : : :    all @xmath125 and @xmath126 vortex    complexes are dynamically stable . @xmath127 : : :    solutions @xmath128 are stable , @xmath125    and @xmath129 are unstable . depending on the strength of    @xmath26 , vortices may start precessing around the center of    the trap , move beyond the boundary forming either a pair of    rarefaction waves with opposite velocities or a complex    @xmath128 , or disappear at the condensate s edge and    bringing in two aligned pairs of vortices @xmath128 and    @xmath130 . [ jjj ] illustrates these possibilities for    @xmath131 and    @xmath132 . @xmath133 : : :    for a given @xmath26 , any sufficiently large    @xmath30 allows the vortex complexes @xmath129    and @xmath125 to stabilize . fig . [ oppvort ] shows such    stabilized @xmath129 complex for    @xmath134     for @xmath135 and @xmath136 ( top row ) , @xmath137 ( second row ) , @xmath138 ( left bottom ) and @xmath139 ( right bottom ) . the initial state is @xmath140 where @xmath141 , @xmath142 and @xmath143 is the vortex core parameter ; see appendix  [ sec : spir - vort - single ] . depicted are the density plots and streamlines of @xmath144 ( left panels in first and second row , bottom row ) and @xmath8 ( right panels in first and second row ) . luminosity is proportional to the magnitude of the density for density plots and to the velocity for streamlines . the size of the pumping is shown as a circle of radius @xmath145 . for @xmath136 the final state consists of two precessing vortices of opposite circulation . for @xmath137 the system evolves into a pair of two gray solitons . for @xmath138 the vortex of the negative circulation in @xmath8 leaves and a vortex of positive circulation enters the cloud . for @xmath139 both vortices leave the condensate and two pairs of vortices of opposite circulation enter . for @xmath138 and @xmath139 @xmath146 , so only @xmath144 is shown . , width=307 ]     for @xmath147 and @xmath137 . depicted are the density plots and streamlines of @xmath144 ( left ) and @xmath8 ( right ) . luminosity is proportional to the magnitude of the density for density plots , and velocity is shown by streamlines . the size of the pumping is shown as a circle of radius @xmath145.,width=307 ]    the stationary state shown in fig . [ jjj ] for @xmath137 is the vorticity - free , but not radially symmetric , solution of eq . ( [ eq:2 ] ) . these are analogous to the rarefaction solitary waves of the nonlinear schrdinger equation discovered by jones and roberts @xcite . in trapped 2d condensates these waves were also found @xcite : they form as two vortices of opposite circulation disappear at the boundary of the condensate . in the conservative gpe these waves propagate with velocities exceeding the velocity of any vortex pair . in spinor damped / driven condensates two such solutions induce flow in opposite directions forming a stationary complex . the plots of the real and imaginary parts of @xmath144 and @xmath8 shown in the left panel of fig . [ grey ] , and the density plots of fig . [ jjj ] can be compared to the bottom panel of fig . 9 and the left panel of fig . 10 of ref .  . it also follows from simulations that @xmath148 , so the found stationary state satisfies a one component ginzburg  landau equation @xmath149\psi + j \psi(x ,- y ) . \label{gl_new}\end{aligned}\ ] ] this also suggests a way to generate solitary waves in spinor condensates . the dark soliton ( obtained by phase imprinting , for instance ) will undergo a transverse snake instability and form two stationary pairs of vortices of opposite circulation , whereas starting with a @xmath129 complex one will obtain stationary rarefaction pulses . ( solid lines ) and @xmath8 ( dashed lines ) of the rarefaction wave complex as a stationary solution of eq . ( [ eq:2 ] ) with @xmath150 for @xmath151 and @xmath152 ( left ) and @xmath153 ( right ) plotted along the axis of density symmetry . the functions are antisymmetric with respect to the origin @xmath154 for @xmath152 , but not for @xmath155.,width=307 ]    as @xmath30 increases , the rarefaction waves in two components lose their antisymmetry and the density minima move away from the center as the right panel of fig . [ grey ] illustrates . for intermediate @xmath30 , complexes combining a rarefaction wave in one component with a precessing vortex in the other component emerge ; for even larger @xmath30 two vortices of opposite circulation move around in the condensate . finally , for even larger @xmath30 , their cores overlap and vortices move to the center of the condensate . these possibilities are illustrated on figs . [ oppvort ] and [ possibledelta ] . the cross - polarized vortices deform the condensate to a slightly oblate form as seen in the density plots of fig . [ oppvort ] . the streamlines and the density plots indicate that the vortices coexist with rarefaction waves that rotate in the counterclockwise direction . ( top row ) and @xmath156 ( bottom row ) for the outcome of instability of cross - polarized vortices obtained by numerical integration of eq . [ eq:2 ] with @xmath150 , @xmath151 and various @xmath30 . the initial state is the same as in fig . luminosity is proportional to density . the size of the pumping is shown as a circle of radius @xmath145 . for @xmath153 two rarefaction pulses form a stationary complex . for @xmath157 a vortex of negative circulation in the @xmath98 component is coupled to a rarefaction pulse in the @xmath158 component . for @xmath159 the two vortices of opposite circulation do not align their cores . for @xmath157 and @xmath159 the complexes precess around the center with the individual vortices moving along an epitrochoid ; see also fig . [ motionvort ] . , width=307 ]    the vortex trajectories with pumping , decay and a spin degree of freedom are nontrivial . for a one - component conservative gpe , a single vortex moves perpendicularly to the background density gradient due to the magnus force@xcite ( the speed of the motion is however nonuniversal , and depends on the global condensate shape @xcite ) . for a two - component conservative gpe with vortices in both components , there is an additional advection of each vortex by the flow pattern of the other component . including also pumping and decay , the trajectories of the vortices are yet more complicated , as illustrated in fig . [ motionvort ] for a @xmath129 complex with @xmath160 . both vortices move along trajectories closely resembling epitrochoids , and the distance between the vortices varies quasi - periodically with time . similarly complicated cycloid trajectories of vortices are known for two - layer fluids with one vortex in each layer  such behavior has been seen for example in models of tropical vortices@xcite . obtained by numerical integration of eq . [ eq:2 ] with @xmath150 , @xmath151 and @xmath159 , and with pumping in a radius @xmath161 . these parameters are the same as for the density profiles in the right hand column of fig . [ possibledelta ] . the trajectory of the vortex of positive ( negative ) circulation in the left ( right ) component is shown as thick black ( thin blue ) line . the start ( end ) points of the time interval plotted are shown as red ( green ) filled circles . note that the vortex trajectories are not closed.,width=153 ] we study the time - dependent problem defined by eq . ( [ eq:2 ] ) with a harmonic trapping potential @xmath150 by numerical integration , using a fourth - order , finite - difference approximation in space and a fourth - order runge  kutta method in time . pumping is restricted to a circular spot of radius @xmath27 centered at the bottom of the trap , as described at the end of section  [ sec : model - spin - polar ] . in all of the following , we take the experimentally realistic parameters@xcite @xmath162 and @xmath163 . with this choice of parameters , the critical radius can be estimated from the thomas  fermi approximation for the problem without polarization@xcite as @xmath164 . this section will present and discuss the transition between synchronized and desynchronized states , and the interplay with the instability to vortex lattice formation . section  [ sec : critical - delta ] first addresses the value of @xmath30 for which a transition occurs , as one changes the size of the pumping spot , and coupling @xmath26 . section  [ sec : attractors ] will then discuss in greater detail the nature of the new attractors near the critical @xmath30 that occur when one includes spatial degrees of freedom . let us start with the simplest trapped problem , taking @xmath165 and @xmath166 . we study the behavior of the system as the detuning @xmath30 is increased in the following way : starting at @xmath152 from a gaussian initial state , we let the system evolve until a steady solution is reached . we then increase @xmath30 in steps , taking for each new @xmath30 the final state at the previous @xmath30 as the new initial state . for each value of @xmath30 , the chemical potentials of the two polarization components are calculated ( see appendix  [ app : numerical - precession ] for details ) . the result is shown in the left panel of fig . [ fig : chempotj1r3-j1r4 ] . the components stay synchronized , sharing a common chemical potential , up to @xmath167 . as the detuning is increased further , the components desynchronize and ac josephson oscillations occur between them ( fig . [ fig : acjosephson ] ) . note that as the components desynchronize , the chemical potentials show an oscillatory time dependence . [ fig : chempotj1r3-j1r4 ] shows a time average in this case .    . for solutions with time - dependent chemical potentials , a time average is shown.,width=307 ]     in the limit of large time . left : @xmath168 as a function of time . right : the chemical potentials of the two polarization components . ( time given in units of @xmath169 . blue and red color denote @xmath158 and @xmath98 components , respectively.),width=307 ]    in terms of a suitably defined , spatially averaged phase difference @xmath40 between the components ( see appendix  [ app : numerical - precession ] ) , the desynchronization transition corresponds to a transition from a fixed point to a limit cycle in the @xmath170 plane , in direct analogy to the two - mode problem ( see fig . [ fig : r3progression ] and compare with fig . [ fig : attraction - basin ] ) .     and @xmath166 . left : synchronized solution for small detuning . right : desynchronized solution in the limit of large detuning ( corrected for numerical precession).,width=307 ]    the right panel of fig . [ fig : chempotj1r3-j1r4 ] shows the chemical potentials of the two polarization components as the detuning is increased stepwise in a system pumped in a spot with radius @xmath171 . the progression from the synchronized solution at small detunings to the desynchronized solution in the limit of large detunings is now more complicated than for a small pumping spot . @xmath171 is too small a pumping radius for vortices to form spontaneously in the corresponding one - component system@xcite . accordingly , the spinor system with zero detuning develops circularly symmetric densities , identical in both components ( but with different phases for the two components ) . for small enough detunings , the system stays synchronized , the densities adjusting to accommodate the common , constant chemical potential ( left - most panels of fig . [ fig : r4progression ] ) .     and @xmath165 . top : phase portraits for the phase difference @xmath40 ( a numerical precession has been removed from the @xmath172 plot ) . middle and bottom : total densities of the @xmath158 and @xmath98 components , respectively.,width=307 ]    as the detuning increases beyond @xmath173 the rotational instability of the high - angular - momentum modes  @xcite reappears , and a four - vortex lattice forms in both components . at this point , the smaller component is sufficiently depressed by the detuning that the effective critical radius is smaller than the pumping radius . the suppressed component becomes rotationally unstable and drags the larger component along . this transition to a vortex - lattice solution causes the ( common ) chemical potential to drop , as shown in fig . [ fig : chempotj1r3-j1r4 ] . while the two components remain on average phase locked , the drop in chemical potential at the formation of the vortex lattice is accompanied by a time - dependence of the chemical potentials , which oscillate together around a common mean . due to the amplitude of the oscillations being slightly larger in the larger component , the fixed point in the @xmath170 plane turns into a small limit cycle , but with only small variations of @xmath40 unlike the @xmath2 periodic cycles in the desynchronized phase . ( middle panels of fig . [ fig : r4progression ] . ) finally , as @xmath30 becomes large , the components desynchronize completely . this situation is shown in the right - hand panels of fig . [ fig : r4progression ] . the difference in space - averaged phase now traces out a limit cycle winding through the full @xmath2 . the overall progression from synchronized solutions at small @xmath30 to desynchronized solutions with josephson oscillations holds very generally for different values of the josephson coupling @xmath26 and also for both small ( no vortices ) and large ( vortex lattice ) pumping - spot sizes . the mechanism in both cases is the same : the synchronized solution is upheld by adjusting the densities and a steady interconversion current forms . after desynchronization , the densities revert to profiles as for the single component condensate , with the average density set by the balance of pumping and decay , but with an additional time - dependent interconversion current . [ fig : densityprof ] shows density profiles of typical synchronized solutions with and without vortices , and also a snapshot of a desynchronized solution . -axis for representative examples of synchronized and desynchronized solutions . left : synchronized solution without vortices ( @xmath174 ) . middle : synchronized solution with vortex lattice ( @xmath175 ) . right : desynchronized solution without vortices ( @xmath176 ) . dashed and solid lines indicate left and right - polarized components , respectively . @xmath165 for all three examples.,width=307 ]    we also performed calculations where we reset the initial conditions to a gaussian density profile , with equal phase of the two polarization components for each value of @xmath30 . based on the results of the two - mode model shown in fig . [ fig : attraction - basin ] , such initial conditions would be expected to find the limit cycle whenever it exists , whereas the stepwise increase of @xmath30 discussed above is intended to follow the fixed point . in these calculations with resetting of initial conditions , desynchronized solutions generally develop at smaller @xmath30 , as one should expect if there is a region of coexistence of synchronized and desynchronized solutions . [ fig : deltac ] shows the highest @xmath30 yielding synchronized solutions for different @xmath27 and different @xmath26 , both for stepwise increased @xmath30 ( top left ) and calculation directly from gaussian initial conditions ( top right ) . comparing the two plots gives an estimate of the region of coexistence . an example ( @xmath177 ) showing the region of bistability is given in the bottom left panel . note that the difference between stepwise increased @xmath30 and reset to gaussian initial condition for each @xmath30 vanishes for large detunings . this need not be the case for larger @xmath27 or larger @xmath26 , where more than one desynchronized , meta - stable state ( distinguished by numbers of vortices ) may be possible at high @xmath30 . yielding synchronized solutions plotted as a function of @xmath27 for @xmath178 and @xmath179 . left : @xmath30 increased stepwise . right : resetting initial conditions for each @xmath30 . bottom left : bifurcation of chemical potentials ( plotted as in fig . [ fig : chempotj1r3-j1r4 ] ) for the two approaches , indicating the region of bistability for the example @xmath177 . bottom right : sketch of regions of stability of synchronized and desynchronized solutions for @xmath166 . compare with inset of fig . [ fig : ft - two - mode ] . , width=307 ]    our calculations suggest that the spatially extended system allows for a rich variation of behaviors . notably , interconversion due to the josephson coupling between the components need not be uniform in space . the interconversion rate is easily calculated from the time - dependent local phase difference @xmath180 as @xmath181 familiar from the josephson effect as found in any textbook@xcite . note that eq . ( [ eq:5 ] ) describes the same physics for the two - mode model . [ fig : interconvmap ] shows a map of the interconversion rate @xmath182 at two close points in time in the large - time limit of the solution directly from gaussian initial conditions with @xmath183 . this solution is desynchronized and exhibits josephson oscillations . however , at any given time there is interconversion in both directions at different points in space .     at two different , but close , times . @xmath183.,width=307 ]    calculations directly from gaussian initial conditions at large detuning with a large pumping spot ( @xmath184 ) indicate that solutions are possible in which the polarization components develop counter - rotating vortex lattices . this leads to rapid density modulations , particularly around the edge of the cloud , and a corresponding pattern of interconversion in opposite directions . in this case , the josephson oscillations are suppressed and the total ( scaled ) mass @xmath168 exhibits rapid , small - amplitude oscillations . these effects are shown in fig . [ fig : counter - rot ] .    , chemical potentials and interconversion rate @xmath182 for a desynchronized solution with counter - rotating components . bottom : densities of the @xmath158 and @xmath98 polarization components , respectively . ( note that the gray scale has been adjusted for each component separately to accentuate the density modulations . ) @xmath185.,width=307 ]      the spatial degree of freedom allows a number of behaviors of the @xmath170 phase portrait that are not possible in the two - mode problem , which only allows fixed points ( synchronized solutions ) and limit cycles with winding number @xmath186 ( desynchronized solutions ) . in the spatially extended system , these behaviors are exemplified in fig . [ fig : r3progression ] for a system without vortices . both these two classes of attractor can also be seen when a vortex lattice exists . in addition , the spatial degree of freedom gives rise to several new classes . we have already noted in fig . [ fig : r4progression ] an example of a synchronized limit cycle . this can be distinguished from the desynchronized limit cycle by the winding number , @xmath187 , over one period . the synchronized limit cycle has winding number @xmath188 , and the desynchronized cycles have winding number @xmath186 . we also find an example of limit cycle with winding number @xmath189 , shown in fig . [ fig : phaseportraitzoo ] . this solution also exhibits another behavior not possible in the two - mode model : a retrograde loop . this solution appears in calculation directly from a gaussian initial condition at @xmath190 , which is barely above @xmath191 for @xmath177 . for larger @xmath30 the loop quickly becomes a cusp and then disappears .     and retrograde loop . left panel shows a blow - up of the loop . a numerical precession has been removed from these plots . ( @xmath192 . ) bottom left : possibly quasi - periodic behavior similar to a limit cycle with winding number @xmath188 ( @xmath193 from stepwise increase ) . bottom right : possible chaotic attractor ( @xmath194 from stepwise increase).,width=307 ]    the bottom panels of fig . [ fig : phaseportraitzoo ] show two other behaviors found when stepwise increasing the detuning in a system with a strong josephson coupling ( @xmath195 ) . both these solutions are basically synchronized ( the time average of the chemical potential is the same for both components ) . however , the large detuning causes the chemical potentials to differ at most instances in time , resulting in behaviors similar to the limit cycles with winding number @xmath188 , but which appear to have a chaotic attractor and/or quasi - periodic behavior . to directly observe rotating vortex lattices in experiments would require time - resolved measurements on timescales of the order of the trap frequency , which would be challenging with current experimental configurations . it is however possible to see signatures of a vortex lattice in the momentum- and energy - resolved photoluminescence spectrum , which can be directly measured in the far field . the spectral weight is given by the modulus squared of the fourier transform of the wavefunction : @xmath196    as an illustration of this , fig . [ fig : spectra ] shows the spectral weight as a function of @xmath197 . the vortex lattice as well as the desynchronization transition can be seen in the dispersion spectrum . [ fig : spectra ] shows the spectra of the two polarization components for two different solutions . the top panels show a synchronized solution that exhibits a vortex lattice without a central vortex . each ring of vortices in the lattices shows up as a side band in the spectrum . if there were no vortices , only the bottom band would be present . as a contrast , the bottom panels show a desynchronized solution that has a vortex lattice with a central vortex . the presence of the central vortex means that the spectral weight vanishes at @xmath198 . note also that the desynchronization causes the spectra of the two components to be shifted relative to each other , whereas in the synchronized case , the two spectra are identical .     for two vortex - lattice solutions . left and right panels show left and right polarized components , respectively . top : synchronized solution without a central vortex . ( @xmath199 . ) bottom : desynchronized solution with a central vortex . ( @xmath200.),width=307 ] the interplay between the effects of pumping and decay in setting the density profile , and the dynamics of the polarization degree of freedom lead to a rich variety of possible nonequilibrium steady states , as well as dynamical attractors for the polarized polariton condensate in a harmonic trap . even neglecting spatial currents , the two - mode model shows nontrivial behavior as a function of applied magnetic field : with strong damping there is a simple transition between synchronized and desynchronized states , whereas for weak damping , a region of bistability exists , in which the state established depends on the initial conditions .    allowing for spatial fluctuations in a homogeneous system , the plane - wave modes on top of the synchronized solution show a second class of transition : if the phase locking term for the spin degree of freedom is large enough , then at weak magnetic fields , this term will provide sufficient restoring force for spin fluctuations to give underdamped global spin oscillations . as the magnetic field increases , the restoring force for such spin waves decreases , eventually vanishing when the synchronized solution becomes unstable . before this instability occurs , there is a transition between underdamped and overdamped spin oscillations , and at this transition , the spin wave energies have a linear dispersion vs momentum .    introducing vortices , the phase locking term favors co - alignment of vortices in the two spin polarizations , but at sufficiently large @xmath30 anti - aligned vortex pairs may become stable . the phase - locking term also makes it possible to obtain a stationary complex of rarefaction waves . the experimental realization of such solitary vorticity - free waves would demonstrate another aspect of superfluid behavior@xcite in the incoherently pumped polariton system . considering the transition between synchronized and desynchronized states in the inhomogeneous system , for sufficiently small pump spots , the transition is similar to the two - mode model . for larger pump spots however , an instability to vortex - lattice formation may preempt desynchronization , driven by the density imbalance required to sustain a synchronized solution . for both large and small pumping - spot radius , the phase portraits near the critical detuning can be more complicated , showing synchronized limit cycles with winding number 0 , desynchronized limit cycles with winding number 2 , and chaotic behavior .    in conclusion , the results presented here illustrate the wide variety of dynamical behavior that can arise in spinor polariton condensates , and suggest that experimental efforts to investigate such behavior should be feasible . this appendix will discuss spiral vortices in the one - component version of eq . ( [ eq:2 ] ) , i.e. , neglecting @xmath26 and @xmath30 . it was shown in ref .   that below the critical radius of pumping for the formation of a vortex lattice there exist states with one or a few spiral vortices . by combining vorticity with pumping and decay , one has both radial and azimuthal supercurrents , both of which modify the density profile , and so these currents interact . this is shown in fig .  [ vortices ] , which shows both the radial velocity ( main figure ) for vortex solutions , and the density profile ( inset ) . the vortex solution requires vanishing density at the origin , this makes the region at small radii a region of net gain , and leads to an outward - flowing current . thus , in the vortex solutions , there exists both a local maximum of radial current , and a point at finite radius ( of the order of the healing length ) at which this current vanishes . , of the vortex solutions of eq . ( [ eq : gpe_rescale ] ) for @xmath201 . the numbers next to the lines indicate the winding numbers of the vortices ( 1 and 2 ) and the ground state ( 0 ) . the tf approximation of the gradient flow of the ground state for small @xmath51 and @xmath202 @xmath203 , is given by the dashed line . the inset compares the density , @xmath204 , of a vortex solution for @xmath205 ( @xmath206 , gray ( red ) line ) with the density of the gpe vortex ( @xmath207 ) with the same number of particles ( @xmath208 , @xmath209 , black line ) . both vortex solutions have topological charge @xmath186.,width=307 ]    in polar coordinates @xmath210 the wavefunction of a spiral vortex of topological charge @xmath211 in the one - component polariton condensate takes the form @xmath212 $ ] , with the equations governing @xmath213 and @xmath214 given by @xmath215 around the center of the trap @xmath216 and @xmath217 can be found recursively in the form of the power series @xmath218 to the leading order @xmath219 showing that the stronger the pumping the larger the outward velocity is in the vortex core . only vortices of topological charge @xmath186 are dynamically stable . for these , @xmath220 , where the numerical value of the vortex core parameter @xmath221 was first calculated by pitaevskii @xcite for the straight line vortex of uniform gross  pitaevskii condensate . we note that around and above critical detuning , solutions to eq . ( [ eq:2 ] ) may exhibit irregular density profiles and irregular vortex lattices . this happens most frequently when solutions are found starting from gaussian initial conditions , but it may also happen for large @xmath27 when @xmath30 is increased stepwise . as an example , we show a solution obtained from gaussian initial conditions at @xmath222 in fig .  [ fig : irregular - density ] . this solution shows an irregular vortex lattice , reminiscent of turbulent behavior . in order to numerically find the chemical potential of a condensate with a rotating vortex lattice , it is useful to use : @xmath223 along with the equation : @xmath224 in many cases , eq . ( [ eq:15 ] ) with @xmath225 is sufficient to determine @xmath226 . however , when there is a vortex at the center of the cloud , the numerator and denominator both vanish , whereas one of @xmath227 will give a well defined value . in such cases , it is then necessary to use eq . ( [ eq:16 ] ) and the integral for @xmath228 to eliminate @xmath229 .    for plotting the phase portraits , it is necessary also to calculate the associated phase , @xmath33 such that @xmath230 . this can similarly be defined as : @xmath231\right\}.\ ] ] however , when there is a central vortex , the same difficulty will arise . to avoid this problem , two possible approaches are used . in some most cases it is sufficient to calculate @xmath232 , and numerically integrate this equation to find @xmath233 . numerical errors in evaluating @xmath234 can however introduce unphysical precession into the plot of @xmath235 . to avoid this , in those cases where the nature of the portrait is simple , we use @xmath236 , with @xmath237 adjusted to remove the spurious precession . this method has been used in figs . [ fig : r3progression ] , [ fig : r4progression ] ( right hand panel ) and [ fig : phaseportraitzoo ] .    where the phase portrait involves finer structure , such as the middle panel of fig . [ fig : r4progression ] , the phase portrait has instead been obtained by : @xmath238\right.\\      & -      \left.\ln\left[\int e^{i s \vartheta } \psi_l(\mathbf{r})\,d^2r\right ]    \right\ } ,    \end{split}\ ] ] with @xmath239 according to whether there is a central vortex . since the two components are corotating in this case , such a definition satisfies : @xmath240 as required .
Abstract:

Answer: quasiparticles in semiconductors  such as microcavity polaritons  can form condensates in which the steady - state density profile is set by the balance of pumping and decay . by taking account of the polarization degree of freedom for a polariton condensate , and considering the effects of an applied magnetic field , we theoretically discuss the interplay between polarization dynamics , and the spatial structure of the pumped decaying condensate . if spatial structure is neglected , this dynamics has attractors that are linearly polarized condensates ( fixed points ) , and desynchronized solutions ( limit cycles ) , with a range of bistability . considering spatial fluctuations about the fixed point , the collective spin modes can either be diffusive , linearly dispersing , or gapped . including spatial structure , interactions between the spin components can influence the dynamics of vortices ; produce stable complexes of vortices and rarefaction pulses with both co- and counter - rotating polarizations ; and increase the range of possible limit cycles for the polarization dynamics , with different attractors displaying different spatial structures .