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: frictionless flow is one of the defining features of superfluidity @xcite . for a moving obstacle with velocity @xmath2 in a superfluid , the frictionless nature of the superfluid near the obstacle breaks down when @xmath2 exceeds a certain critical velocity @xmath1 . according to landau s criterion this critical velocity is estimated as @xmath3 $ ] , where @xmath4 is the excitation spectrum , @xmath5 is the planck constant , and @xmath6 is the wave vector , with @xmath7 , see refs . an object moving with a velocity above @xmath1 dissipates energy via the creation of elementary excitations , for example , vortices or phonons . superfluidity was first observed in liquid helium @xmath8 and helium @xmath9 . since then , superfluidity has been studied in quantum gas systems of bosons @xcite , fermions @xcite , as well as of bose - fermi mixtures @xcite . the phenomenon of superfluidity is closely related to the bose - einstein condensation ( bec ) of interacting gases . interestingly , a uniform two - dimensional ( 2d ) system can not undergo the bec transition because the formation of long - range order is precluded by thermal fluctuations @xcite . however , it forms a superfluid with quasi - long range order via the berenzinskii - kosterlitz - thouless ( bkt ) mechanism @xcite . the quasi - long range order of this state refers to the algebraic decay of the single - particle correlation function . the algebraic exponent of this correlation function increases smoothly with temperature . at the critical temperature , the superfluid density of the system undergoes a universal jump of @xmath10 , where @xmath11 is the de broglie wavelength . experiments on 2d bosonic systems , such as a liquid helium film @xcite , and trapped bose gases @xcite have shown indications of the bkt transition . furthermore , a trapped 2d system can form a bec due to the modified density of states @xcite and leads to an interesting interplay of the two phase transitions @xcite . quasi - long range order in 2d bosonic systems can be detected via interference and time - of - flight techniques @xcite . however , as a direct method , superfluidity of ultracold atomic gases was probed using a local perturbation , in particular via laser stirring . for example , superfluidity of 3d becs was tested via laser stirring in refs . @xcite . in the experiment @xcite , thermal relaxation of a perturbed 2d quasicondensate is studied . @xcite reported on stirring a trapped 2d bose gas of @xmath0rb atoms with a blue - detuned laser , moving on a circular path around the trap center . the circular motion ensures that the harmonically trapped 2d gas is probed at a fixed phase - space density . by choosing different radii of the circular motion , the superfluid transition was explored . in this paper , we provide a quantitative understanding of the experiment using a c - field simulation method . we demonstrate that a blue - detuned laser of intensity comparable to the mean - field energy causes dissipation due to the creation of vortex - antivortex pairs . this is in contrast to laser stirring with a red - detuned laser @xcite , where dissipation occurs via phonons @xcite . furthermore , we study the relaxation dynamics of the stirred gas following the stirring process , which shows a slow energy transport between the condensate and the thermal cloud . we identify the origin of this slow relaxation to be vortex recombination and diffusion . we show that this effect can explain quantitatively the shift of the measured critical phase - space density in the experiment . this paper is organized as follows . in sec . [ sec_sim_method ] we describe the simulation method that we use . in sec . [ sec_vc_sim ] we determine the critical velocity @xmath1 of the stirred gas , based on which we identify the superfluid to thermal transition . in sec . [ sec_mechanism ] we discuss the dissipation via vortex pairs . in sec . [ sec_comp_expt ] we compare the simulation results with the experiment . in sec . [ sec_nonequilibrium ] we analyze the relaxation of the stirred gas , and in sec . [ sec_conc ] we conclude . we simulate the stirring dynamics of a weakly interacting 2d bosonic system using the c - field simulation method that we used for a 3d system in ref . we describe this method in the following . we start out with the hamiltonian of the unperturbed system , @xmath12.\end{aligned}\ ] ] @xmath13 and @xmath14 are the bosonic annihilation and creation operator , respectively . the 2d coupling parameter @xmath15 is given by @xmath16 , where @xmath17 is the dimensionless interaction , @xmath18 is the atomic mass , @xmath19 is the 3d @xmath20-wave scattering length , and @xmath21 is the harmonic oscillator length of the confining potential @xmath22 in the @xmath23 direction . @xmath24 is the trap frequency along the @xmath23 direction . @xmath25 describes the external potential , which is a harmonic trap , @xmath26 . @xmath27 is the trap frequency in the radial direction and @xmath28 is the radial coordinate . we introduce a time - dependent term to describe laser stirring , @xmath29 where @xmath30 is the time - dependent stirring potential and @xmath31 is the density operator at the location @xmath32 . the stirring potential is a gaussian with a width @xmath33 and a strength @xmath34 , @xmath35 which is centered at @xmath36 . we move @xmath37 along a circular path as a function of time @xmath38 . we perform numerical simulations by mapping this system on a lattice system , which also introduces a short - range cutoff ; see appendix [ app_sim_heating ] . this short - range cutoff is of the order of the healing length @xmath39 , with @xmath40 being the density . we describe both the equations of motion and the initial state within a c - number representation , which corresponds to formally replacing the operators @xmath13 by complex numbers @xmath41 . furthermore , we approximate the initial ensemble by a classical ensemble , within a grand - canonical ensemble of temperature @xmath42 and chemical potential @xmath43 . we sample the initial states via a classical metropolis algorithm . the simulation setup consists of a disc - shaped 2d circular condensate of @xmath0rb atoms . this choice of the 2d circular condensate is inspired by the experimental setup of ref . @xcite . in the simulations we consider @xmath44 @xmath0rb atoms confined by the harmonic potential in both the radial and transverse directions . the trap frequencies are @xmath45 and @xmath46 . here the scattering length is @xmath47 , which yields @xmath48 . the temperature of the trapped gas is in the range @xmath49 . the simulation parameters that we use , are in the typical range of the experimental parameters of ref . @xcite . for simulations of a quasi- and a pure-2d trap geometry we use a lattice of @xmath50 and @xmath51 sites , with the lattice discretization length @xmath52 , respectively . we choose @xmath53 such that it is smaller than , or comparable to , the healing length @xmath54 and the de broglie wavelength @xmath55 , see ref . the trapped gas is in the pure-2d regime if @xmath56 . when @xmath57 and @xmath43 are comparable to @xmath58 , it is in the quasi-2d regime . after initializing the trapped system at temperature @xmath42 , we switch on the stirring potential described by eq . [ eq_stir_potential ] . in the experiment @xcite the trapped gas is stirred with a blue - detuned laser beam moving on a circular path around the trap center . for the circular motion of stirring we choose @xmath59 , where @xmath60 and @xmath61 are the stirring radius and frequency , respectively . for the stirring potential we use the strength @xmath62 and the width @xmath63 , in accordance with the experiment . the stirring sequence is the following : we linearly switch on the stirring potential over @xmath64 , let it stir the system for @xmath65 , and then switch it off over @xmath64 . this is again inspired by the experimental choices . we repeat this for various stirring velocities @xmath66 by changing both @xmath60 and @xmath61 . by choosing different @xmath60 we stir the different regimes of the trapped gas , the superfluid , the thermal , and the crossover regime .    after stirring we calculate the total energy @xmath67 using the unperturbed hamiltonian in eq . [ eq_hamil ] , where we use @xmath41 instead of @xmath13 . from this energy we determine the equilibrium temperature @xmath68 of the stirred gas . we infer this temperature by numerically inverting the temperature dependence of the equilibrium state , @xmath69 . we elaborate on this in appendix [ app_sim_heating ] . from the temperature difference between the stirred and initial system , the heating @xmath70 is determined . we also calculate the local energy , as well as the vortex and anti - vortex distribution . we define the local energy as @xmath71 , where @xmath72 refers to the nearest neighbor sites . @xmath73 , @xmath74 , and @xmath75 are the complex - valued field , the density , and the trap potential at site @xmath76 , respectively . @xmath77 and @xmath78 are the bose - hubbard parameters , see appendix [ app_sim_heating ] . for the vortex distribution , we calculate the phase winding around the lattice plaquette of size @xmath79 , using @xmath80 , where the phase differences between sites is taken to be @xmath81 $ ] . @xmath82 is the phase field of @xmath41 . we identify a vortex and an antivortex by a phase winding of @xmath83 and @xmath84 , respectively . by counting all vortices and antivortices we determine the total number of vortices . we restrict this counting to the the superfluid region of the gas as we describe below . to study the superfluid behavior we stir a 2d quasicondensate with a repulsive gaussian potential . we prepare a trapped 2d quasicondensate of @xmath85 @xmath0rb atoms at temperature @xmath86 . we show the simulated density profile of the trapped gas in fig . [ fig_heating](a ) . we stir the gas with a circularly moving , repulsive stirring potential at stirring radius @xmath87 . as mentioned in sec . [ sec_sim_method ] , we use the strength @xmath62 and the width @xmath63 for the stirring potential . this strength @xmath34 is well above the local mean - field energy @xmath88 at the stirring location . after stirring we determine the induced heating @xmath89 from the equilibrium temperature @xmath90 of the stirred gas , [ sec_sim_method ] for details . by varying the stirring frequency @xmath61 we determine @xmath91 as a function of stirring velocity @xmath92 . we show @xmath93 determined for various @xmath2 in fig . [ fig_heating](b ) . the induced heating is almost negligible at low @xmath2 , its onset occurs at a velocity @xmath1 , and for @xmath94 it increases rapidly . we quantify the onset of heating using a fitting function , @xmath95 + b,\ ] ] which is discussed in ref . @xcite , with the free parameters @xmath96 , @xmath97 , and @xmath1 . for the simulated heating shown in fig . [ fig_heating](b ) , this function gives a critical velocity of @xmath98 . we compare this critical velocity to the bogoliubov estimate of the phonon velocity @xmath99 at the stirrer location . the bogoliubov velocity is determined by @xmath100 . the observed critical velocity is @xmath101 . this is notably different from the case of an attractive stirring potential , where @xmath102 @xcite . we explain this reduction of @xmath1 for a repulsive stirring potential in sec . [ sec_mechanism ] .    by choosing different radii @xmath60 we explore the various regimes of the trapped gas . we use the same strength @xmath34 and the same width @xmath33 as above . for each @xmath60 , we first determine the induced heating @xmath91 as a function of @xmath2 , and then by using the fitting function given in eq . [ fit ] we determine @xmath1 . we show @xmath1 determined at various @xmath60 in fig . [ fig_heating](c ) . the stirring radii are in the range @xmath103 . for @xmath104 , there is no significant change of @xmath1 . as @xmath60 reaches the crossover regime , @xmath1 is reduced sharply and for @xmath60 above the crossover regime , @xmath1 is zero . according to the bkt prediction in a uniform system @xcite with @xmath105 combined with local density approximation , the crossover regime should occur at @xmath106 . this prediction is in good agreement with the crossover regime identified by the simulated @xmath1 . thus , we clearly identify the superfluid , the crossover , and the thermal regimes by the finite , the sharply decreasing , and the zero critical velocities @xmath1 , respectively . we note that in the crossover region the decrease of @xmath1 is as sharp as the size of the stirrer allows . furthermore , we note that the observed almost constant @xmath1 for @xmath107 can be due to the accelerated circular motion and the large strength of the stirring potential @xcite . the observed critical velocities are in the range @xmath108 . to understand what leads to this reduction of the critical velocity with regard to the phonon velocity , we investigate the time evolution of the phase field @xmath109 of a single realization of the thermal ensemble . we obtain this phase field from the complex field @xmath110 via the phase - density representation @xmath111 . in fig . [ fig_dissipation ] we show the phase evolution of the trapped 2d quasicondensate stirred at @xmath112 . we use the velocity @xmath113 , which is above the steep onset of dissipation related to the breakdown of superfluidity . the phase evolution of the unperturbed gas shows rather weak phase gradients . as stirring is switched on , the phase field around the stirring potential starts to fluctuate . these fluctuations develop into strong phase gradients , which result in the creation of vortex - antivortex pairs . this can be confirmed by calculating the phase winding around each plaquette of our numerical grid , as described in sec . [ sec_sim_method ] . we show the calculated phase winding in fig . [ fig_dissipation ] , where vortices and antivortices are shown as circles and triangles , respectively . this indeed shows the creation of vortex - antivortex pairs during stirring . we recall that the stirring strength @xmath62 is much larger than the mean - field energy @xmath114 at the stirring location , which results in a strong reduction of the density at the stirrer location . this density reduction serves as a nucleation site for the creation of vortex pairs . we note that this mechanism of vortex - pair - induced dissipation is suppressed for an attractive stirring potential , as shown in ref . this scenario of dissipation induced by vortex pair creation is consistent with a recent experiment @xcite . we now compare the results of our simulation with the experiment @xcite . we first show the comparison between the experiment and simulation for the heating as a function of @xmath2 . in the superfluid regime , we stir the quasicondensate at the radius @xmath87 . the simulated density profile is shown in the inset of fig . [ fig_heating_comp_sf](a ) . after stirring we let the stirred gas relax for @xmath115 of relaxation time and then determine the induced heating from the temperature of the wings of the cloud . we fit these wings to the hartree - fock prediction , @xmath116,\ ] ] with the fitting parameters @xmath42 and @xmath117 . this method is adopted according to experiment , in which the temperature of the stirred gas is determined in the same way , following a relaxation of @xmath115 as well . we denote this heating determined from the wing temperature @xmath118 by @xmath119 , with @xmath42 being the initial temperature . we show the simulated @xmath120 and their comparison with the experimental ones for various @xmath2 in fig . [ fig_heating_comp_sf](a ) . the measured and simulated heating are found to be in good agreement if we base the comparison on @xmath120 . we also compare the measured @xmath120 with the simulated @xmath121 determined from the equilibrium temperature @xmath68 of the stirred gas . we show @xmath121 as open circles in fig . [ fig_heating_comp_sf](a ) . they show agreement at low and intermediate velocities @xmath2 , whereas they differ at large @xmath2 . this noticeable difference at large @xmath2 is due to the absence of global thermal equilibrium of the stirred gas . as explained in sec . [ sec_nonequilibrium ] , the stirred gas relaxes by transporting the excess energy between the superfluid in the central part and the thermal cloud in the periphery , which is a slow process . the absence of global thermal equilibrium leads to a smaller wing temperature than the equilibrium temperature . the results shown in fig . [ fig_heating_comp_sf](a ) indicate that the onset of heating occurs at a velocity @xmath1 , and for @xmath122 heating increases rapidly . both in experiment and simulation @xmath1 is determined using the fitting function in eq . [ fit ] . in fig . [ fig_heating_comp_sf](b ) we show the comparison between the experiment and simulation for stirring the thermal region of the trapped gas of @xmath123 atoms . the simulated density profile of the gas is given in the inset of fig . [ fig_heating_comp_sf](b ) . both the measured and simulated heating @xmath120 are in good agreement . the simulated @xmath121 determined from the equilibrium temperature of the system are below the measured @xmath120 at large @xmath2 . as we will explain in sec . [ sec_nonequilibrium ] , this is again due to the absence of global thermal equilibrium . as the stirred thermal cloud has more excess energy than the condensate , the wing temperature is larger than the equilibrium temperature . the results shown in fig . [ fig_heating_comp_sf](b ) indicate that heating occurs at all @xmath2 , which results in a zero @xmath1 . next , we show in fig . [ fig_heating_comp_sf](c ) the comparison between the experiment and simulation for @xmath1 that are determined by stirring the superfluid , the crossover , and the thermal regime . in the experiment @xcite @xmath1 is measured for different configurations of the total number of atoms @xmath124 , the temperatures @xmath42 , and the stirring radii @xmath60 . we compare the measured @xmath1 with the simulated @xmath1 determined by stirring the 2d gas in sec . [ sec_vc_sim ] . we show both the measured and simulated @xmath1 as a function of the dimensionless parameter @xmath125 . the parameter @xmath126 characterizes the degree of degeneracy of the cloud and is the relavant parameter in the sense that the thermodynamic properties of the gas depend only on the ratio @xmath127 @xcite . we refer to @xmath1 determined from the wing temperature @xmath118 and from the equilibrium temperature @xmath68 as @xmath128 and @xmath129 , respectively . both the measured and simulated @xmath128 show good agreement . the measured @xmath128 and the simulated @xmath129 agree in the superfluid and thermal regime , while they differ in the crossover regime . for the measured and simulated @xmath128 , the crossover regime occurs at @xmath130 and @xmath131 , respectively . however , for the simulated @xmath129 , it occurs at @xmath132 . the theoretical prediction for the bkt transition in a uniform gas @xcite with @xmath133 occurs at @xmath134 . this prediction is comparable to the simulated crossover regime identified by @xmath129 , whereas its comparison with the crossover regimes identified by the measured and simulated @xmath128 shows a shift . this shift was observed in ref . @xcite , but could not be explained . we conclude that the experiments of ref . @xcite can be reproduced quantitatively if the wing temperature is used , rather than @xmath68 . this suggests that the system has not relaxed to thermal equilibrium after the waiting time of @xmath115 . we confirm and elaborate on this point and the underlying mechanism in the following sections . we now investigate the relaxation of the system , following the stirring process in the superfluid regime . this includes a discussion of the influence of the confinement of the system in the @xmath23 direction . for strong confinement , the system approaches a purely 2d limit , while it is quasi-2d for intermediate confinement . we first analyze the energy - flow dynamics of a stirred trapped gas in the purely 2d limit , and then compare this dynamics with a quasi-2d gas . for a pure-2d trapped gas , we consider a gas of @xmath135 @xmath0rb atoms , which is strongly confined in the transverse direction by the harmonic potential . the temperature @xmath86 is smaller than the transverse trap energy @xmath136 , so that the gas is in the ground state in this direction . as the width of the condensate in the @xmath23 direction is smaller than the lattice discretization length @xmath53 , we simulate this system using a single @xmath137-layer of lattice only , see sec . [ sec_sim_method ] . we stir the gas at @xmath138 for @xmath139 at a velocity @xmath94 . after that we switch off the stirring potential and let the gas relax . we calculate the local energy @xmath140 of the stirred gas and its final equilibrated local energy @xmath141 , as described in sec . [ sec_sim_method ] . we show the evolution of the excess energy @xmath142 for various relaxation times @xmath143 in fig . [ fig_energyevolution](a ) . @xmath144 is the maximum density of the system . the evolution of @xmath145 after @xmath146 shows that most of the stirring - induced energy resides within the superfluid region . the system then relaxes by transporting this excess energy to the thermal cloud . this process occurs slowly and the system achieves fully equilibration only after about @xmath147 relaxation time , remarkably . we now study this energy - flow dynamics for the case of a quasi-2d gas . we consider the quasicondensate that we use in sec . [ sec_vc_sim ] . the initial temperature of the gas and the harmonic potential in the transverse direction are equal and half of those in the pure-2d case , respectively . the resulting system is a quasi-2d gas . we simulate this system using five @xmath137-layers of lattice in the @xmath23 direction . we stir the gas using the same stirring parameters as for the pure-2d case . we show the evolution of the excess energy @xmath145 of the stirred gas for various @xmath143 in fig . [ fig_energyevolution](b ) . in this case , @xmath140 , @xmath141 and @xmath144 are the column ( i.e. integrated along the @xmath23 axis ) quantities . the evolution of @xmath145 after @xmath146 is similar to the pure-2d gas . again , the system relaxes by transporting the excess energy to the thermal cloud . the equilibration process is slightly slower than for the pure-2d gas but again of the same order of @xmath147 . we note that the measured relaxation times in ref . @xcite are indeed of the order of the relaxation times that we find here .    in the experiment @xcite the waiting time after stirring and before measurement is @xmath148 , which is shorter than the relaxation times that we observe here . this indicates that in the experiment thermal equilibrium between the superfluid and thermal cloud is not fully established , which influences the measured heating . it leads to respectively lower and higher measured heating for stirring the superfluid and thermal parts of the cloud .      to understand what causes this slow relaxation for the system , we now examine the evolution of the density and vortices of the system . we calculate the local density as @xmath74 and vortices as described in sec . [ sec_sim_method ] . we show the density and vortex evolution of a single realization of the stirred pure- and quasi-2d gas after various @xmath143 in fig . [ fig_vortexdynamics ] . for both systems , the density relaxation is hard to recognize , whereas the vortex evolution clearly exhibits decay of vortices . thus , the system relaxes via decay of the induced vortices . vortices can decay via both annihilation of a vortex with an antivortex , and drifting out to the thermal region of the cloud . for the pure-2d gas the number of vortices after @xmath146 is larger than the quasi-2d gas . @xcite reported that vortex annihilation in a pure-2d system is strongly suppressed as compared to a quasi-2d system because the vortex lines are impermeable to tilting @xcite and bending @xcite . so , the suppression of vortex annihilation can be a reason for the long - lived vortices for the pure-2d gas .     as a function of time @xmath38 for the pure- and quasi-2d gas . we fit the vortex decay with the bi - exponential function in eq . [ eq_fit_decay ] and determine the decay times @xmath149 ( see text ) . the fitted curves are shown as the dashed lines . ]    to make a quantitative comparison for the vortex relaxation between the two system , we count the total number of vortices within the superfluid region of the cloud at a detection radius @xmath150 , and average it over @xmath151 realizations . we show the averaged vortex number @xmath152 as a function of time @xmath38 for both systems in fig . [ fig_vortexdecay ] . as stirring is switched on at time @xmath153 , @xmath152 starts to increase approximately linearly . it reaches its maximum @xmath154 at @xmath155 . after the stirring is switched off , it decays approximately exponentially . for both systems the nature of vortex growth and decay are the same , but the rates with which they grow and decay are different . for the pure-2d gas the growth and decay rates are larger and smaller than those for the quasi-2d gas , respectively . the enhanced growth and the suppressed decay rate for the pure-2d gas can be due to the suppression of vortex annihilation , as mentioned above , and a slow vortex drift . we quantify the vortex decay rate using the function , @xmath156 with the free parameters @xmath157 . from the fit , we determine @xmath158 and @xmath159 for pure- and quasi-2d gas , respectively . these decay times are similar to those determined from the mean excess energy in appendix [ app_energy_flow ] . the fast decay @xmath160 and the slow decay @xmath161 are essentially connected to the vortex annihilation and drift lifetime , respectively . for the pure-2d gas @xmath160 and @xmath161 are larger than and equal to those for the quasi-2d gas , respectively . we show in table [ tab : rdet ] @xmath154 and the extracted @xmath149 at varying @xmath162 , for both systems . @xmath160 and @xmath161 increase weakly as @xmath162 is increased . however , the following conclusions are essentially independent of the choice of @xmath162 . overall , @xmath154 and @xmath160 are larger for pure-2d gas than those for quasi-2d gas , respectively , while @xmath161 are similar for both systems . we compare the simulated @xmath149 of quasi-2d gas to the waiting time of @xmath148 in the experiment @xcite . this time is twice as large as the fast decay @xmath160 , whereas it is smaller than the slow decay @xmath161 . this suggests that most vortex recombination processes have occurred at the time of the measurement . however , the vortex drift to the thermal cloud has not occurred , and the system is in a metastable state , not in the equilibrated state . this is the mechanism that is responsible for the difference between the wing temperature and the equilibrium temperature . we note that in ref . @xcite an estimate for the time of a vortex line drifting to the thermal cloud was given . while this estimate was for a three dimensional system , we find that the analytical estimate of ref . @xcite gives a timescale that is consistent with our simulation . we also note that the vortex lifetime is suppressed at high tempearures @xcite .      in this section we elaborate on the relaxation dynamics of the stirred trapped gas that we discuss in sec . [ sec_nonequilibrium ] . we first discuss the energy flow dynamics in sec . [ app_energy_flow ] and then the vortex dynamics in sec . [ app_vortex_dynamics ] .      here we elaborate on the energy flow dynamics for the pure- and quasi-2d trapped system . we stir both systems with the stirring potential at velocity @xmath177 . after stirring we calculate the excess energy @xmath142 as described in sec . [ sec_nonequilibrium ] and by averaging this energy over the superfluid region of the gas , we calculate the mean energy @xmath181 . we show the evolution of @xmath181 for both systems in fig . [ fig_app_energy_decay](a ) . @xmath181 decays approximately exponentially as the excess energy @xmath145 within the superfluid region outflows to the thermal cloud . we quantify the energy decay time using the fitting function in eq . [ eq_fit_decay ] . from the fit , we determine the decay times @xmath182 and @xmath183 for pure- and quasi-2d gas , respectively . these decay times are similar to the ones that we determine from the vortex decay in sec . [ sec_vortex ] . in addition to the stirring at @xmath177 , we show @xmath181 corresponding to stirring at @xmath184 for the pure- and quasi-2d gas in fig . [ fig_app_energy_decay](b ) and ( c ) , respectively .        next , we turn to the vortex relaxation of the stirred gas . we stir the pure- and quasi-2d gas using different stirring radii @xmath60 and stirrer strengths @xmath34 . we calculate the total number of vortices within the superfluid region of the cloud , as described in sec . [ sec_vortex ] , and average it over @xmath186 realizations . we show the averaged vortex number @xmath152 normalized by its maximum vortex number @xmath185 as a function of time @xmath38 for the pure- and quasi-2d gas in figs . [ fig_app_vortex_decay](a ) and [ fig_app_vortex_decay](b ) , respectively . for the pure-2d gas the relaxation of vortices is slower than that for the quasi-2d gas , as shown in sec . [ sec_vortex ] .            c. raman , m. khl , r. onofrio , d. s. durfee , c. e. kuklewicz , z. hadzibabic , and w. ketterle , phys . 83 * , 2502 ( 1999 ) ; r. onofrio , c. raman , j. m. vogels , j. r. abo - shaeer , a. p. chikkatur , and w. ketterle . phys . lett . * 85 * , 2228 ( 2000 ) . i. ferrier - barbut , m. delehaye , s. laurent , a.t . grier , m. pierce , b. s. rem , f. chevy , and c. salomon , science * 345 * , 1035 ( 2014 ) ; m. delehaye , s. laurent , i. ferrier - barbut , s. jin , f. chevy , and c. salomon , phys . 115 * , 265303 ( 2015 ) . v. bagnato and d. kleppner , phys . a , * 44 * 7439 ( 1991 ) ; d. s. petrov , d. m. gangardt , and g. v. shlyapnikov , j. phys . iv * 116 * , 5 ( 2004 ) ; h. holzmann , m. chevallier , and w. krauth , europhys . lett . * 82 * , 30001 ( 2008 ) . we have studied the superfluid to thermal transition of a trapped 2d bose gas of @xmath0rb atoms by stirring it with a repulsive stirring potential on a circular path around the trap center . the superfluid transition was probed by choosing different radii of the circular motion . we have identified the superfluid , the crossover , and the thermal regime by the finite , the sharply decreasing , and the zero critical velocity , respectively . the superfluid region of the gas yields critical velocities that are in the range @xmath163 , where @xmath164 is the phonon velocity . we have demonstrated that the onset of dissipation is due to the creation of vortex - antivortex pairs . the comparison of the simulation with the experiment shows good agreement if the temperature measurement of the experiment is imitated in the simulation , i.e. by extracting the wing temperature . however , we confirm the systematic shift that was observed in experiment , if thermal equilibrium is assumed . we have demonstrated that the absence of thermal equilibrium after the waiting time that was used in experiment is due to a remarkably slow relaxation mechanism : the energy transport across the superfluid to thermal interface occurs only on timescales of seconds . this slow transport mechanism is due to the slow drift of vortices out of the superfluid into the thermal wings of the system . we emphasize that this mechanism is relevant for many on - going experiments in the field of ultracold atoms , and their temperature measurements . furthermore , this effect of suppressed transport across critical interfaces is in itself intriguing , and could be studied in a future cold atom experiment with clarity . we acknowledge support from the deutsche forschungsgemeinschaft through grants no . ma 5900/1 - 1 and no . sfb 925 , the hamburg centre for ultrafast imaging , and from the landesexzellenzinitiative hamburg , supported by the joachim herz stiftung . jd thanks jrme beugnon for many fruitful discussions and acknowledges support from erc ( synergy grant uquam ) . in this section we show how we determine the equilibrium temperature @xmath68 of a stirred gas using the c - field method described in sec . [ sec_sim_method ] . we discretize the continuum hamiltonian in eq . [ eq_hamil ] by the bose - hubbard hamiltonian @xcite on a 2d square lattice , @xmath165 @xmath73 and @xmath74 are the complex - valued field and the density at site @xmath76 , respectively . @xmath166 indicates nearest - neighbor bonds . for a lattice discretization length @xmath53 , the bose - hubbard parameters are related to the continuum parameters via @xmath167 and @xmath168 . the 2d coupling parameter @xmath15 is given by @xmath16 , where @xmath169 is the dimensionless interaction , @xmath18 is the atomic mass , @xmath19 is the 3d s - wave scattering length , and @xmath170 is the harmonic oscillator length of the confining potential @xmath171 in the @xmath23 direction . @xmath24 is the trap frequency in the @xmath23 direction . the harmonic trapping potential is @xmath172 . @xmath27 is the trap frequency in the radial direction and @xmath173 is the radial coordinate . we first initialize the system in a thermal state at temperature @xmath42 via classical monte carlo , and then calculate its energy @xmath174 using the hamiltonian in eq . [ eq_hubbard ] . by varying the temperature of the system @xmath42 while keeping the total number of atoms @xmath124 fixed , we calculate the energy @xmath175 as a function of @xmath42 . in fig . [ fig_app_temp_map ] we show the temperature dependence of the energy @xmath175 for the pure- and quasi-2d gas that are described in sec . [ sec_nonequilibrium ] .    to determine the heating , we first stir the gas with the repulsive stirring potential as described in sec . [ sec_sim_method ] and then after stirring calculate its energy @xmath175 using eq . [ eq_hubbard ] . we numerically invert this energy @xmath175 to the equilibrium temperature @xmath68 using the temperature dependence shown in fig . [ fig_app_temp_map ] . finally , from the temperature difference between the stirred and initial system , the heating @xmath176 is determined .         for the pure- and quasi-2d gas stirred at velocity @xmath177 . we fit the energy decay with the bi - exponential function in eq . [ eq_fit_decay ] to determine the decay times @xmath149 ( see text ) . the fitted curves are shown by the dashed lines . we show @xmath178 normalized by its maximum mean energy @xmath179 corresponding to @xmath177 and @xmath180 for the pure- and quasi-2d gas in panels ( b ) and ( c ) , respectively . ]
Abstract:

Answer: we investigate the superfluid behavior of a two - dimensional ( 2d ) bose gas of @xmath0rb atoms using classical field dynamics . in the experiment by r. desbuquois _ et al . _ , nat . phys . * 8 * , 645 ( 2012 ) , a 2d quasicondensate in a trap is stirred by a blue - detuned laser beam along a circular path around the trap center . here , we study this experiment from a theoretical perspective . the heating induced by stirring increases rapidly above a velocity @xmath1 , which we define as the critical velocity . we identify the superfluid , the crossover , and the thermal regime by a finite , a sharply decreasing , and a vanishing critical velocity , respectively . we demonstrate that the onset of heating occurs due to the creation of vortex - antivortex pairs . a direct comparison of our numerical results to the experimental ones shows good agreement , if a systematic shift of the critical phase - space density is included . we relate this shift to the absence of thermal equilibrium between the condensate and the thermal wings , which were used in the experiment to extract the temperature . we expand on this observation by studying the full relaxation dynamics between the condensate and the thermal cloud .


Instruction: Write an abstract for the following acadamic paper.

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

Paper: in this paper we address the realization and identification problem for a subclass of piecewise - affine hybrid systems . * contribution of the paper * we define the class of _ piecewise - linear systems _ ( abbreviated by * pwl * ) . * pwl*s are a subclass of piecewise - affine hybrid systems . the continuous dynamics of * pwl *  is determined by a finite collection of affine subsystems . however , in contrast to traditional piecewise - affine systems , we allow any change of the continuous state during a discrete - state transition , as long as the new state belongs to the set of designated initial states of the affine subsystem associated with the new discrete state . in addition , we do not impose any specific mechanism for triggering discrete - state transitions . we formulate the realization problem for this system class and partially solve it by providing necessary and sufficient conditions for existence of a realization . we also present conditions for minimality . we show that the outputs of any * pwl *  can also be described by a switched ar model . the main conclusion is that any system can be transformed to a minimal system _ with one discrete state _ while preserving input - output behavior . this means that without further restrictions , the identification problem for such systems is not necessarily interesting . we discuss a number of restrictions on the system structure which avoid this problem . note that the conclusion above is not valid for other classes of hybrid systems . for hybrid systems from @xcite , there examples of input - output maps which provenly can not be realized by systems with one discrete state . in addition , we present an identification algorithm for systems with full observations . this algorithm is illustrated by example of physical and biological relevance . * motivation * the motivation for studying realization theory for * pwl*s is that it provides the theoretical foundations for systems identification . the motivation for investigating identification of * pwl*s is the following . * first , several systems of interest can be modeled by * pwl*s and * pwl*s tend to be convenient for analysis . however the parameters of * pwl *  models are often not directly available and hence they have to be estimated from measurements . * the problem of estimating the network dynamics of complex biological systems can be reduced to the identification problem for * pwl*s .    * identification of * pwl*s and biological networks * below we will elaborate on the relationship between identification of * pwl*s and the estimation of dynamic interactions of complex dynamical systems . numerous vital processes in nature involve complex signaling networks , varying from gene - protein interaction networks @xcite and complex communication in microbes @xcite to the synchronization in the heart of higher animals @xcite . in most cases , these signaling networks do not have a fixed topology , but vary their structure according to their internal states and certain external conditions @xcite . this flexibility allows the organism to exhibit complex behavior and flexible responses to various environmental conditions . the underlying biological problem is to describe these flexible dynamical networks , based on ( partial ) observations of their internal states . examples of potential observations are time series of gene expressions , protein densities , oxygen stress and sugar concentrations .    in almost all cases , there is no explicit mathematical model available that even remotely describes the network dynamics . this is caused by the sheer size of the problem ; the many thousands of genes , rnas , proteins and other types of molecules involved , the complexity and idiosyncrasy of their interactions , and the large amount of noise in biological systems . basically , many agents involved in the network are unknown let alone their mechanism of interaction . mathematically , the processes above can be viewed as nonlinear dynamical systems with unknown parameters . such systems often have several equilibria , and hence their behavior can be seen as a mixture of certain affine systems , where each affine system is obtained by linearizing the original system around one of the equilibria . that is , such systems can be approximated by a * pwl*. in addition to their simplicity , the advantage of * pwl *  approximations is that they neatly capture the interaction among various state components around various equilibria . more precisely , interaction between the @xmath0th and @xmath1th state components can be viewed as the property that the @xmath2th entry of the jacobian of the system around equilibrium is non - zero . since we are often interested exactly in interactions rather than the detailed dynamics , it makes sense to recast the problem of identifying such interactions into the problem of identifying * pwl*s approximations . while the paper does not directly address identification of interactions , we believe that the presented results represent a step towards the solution of that problem . * related work * identification of piecewise - linear ( -affine ) hybrid systems has been subject of intensive research , @xcite . realization theory for hybrid systems was investigated in @xcite . the application of hybrid systems to modelling complex dynamical systems , originating from biology , and their identification has been a subject of intensive research @xcite . the results of this paper are new , to the best of our knowledge . what distinguishes the contribution of this paper from the existing work is * ( a ) * the system class considered , in particular , the emphasis on continuous - time systems , * ( b ) * the emphasis on realization theory and algorithms , * ( c ) * the details of the identification algorithm . * outline of the paper *  [ sect : pwl ] presents the formal definition of the system class of interest .  [ sect : approx ] discusses the relationship between nonlinear systems with complex dynamics and piecewise - affine hybrid systems .  [ sect : real ] presents the results on realization theory , and  [ sect : ident ] presents the identification algorithm and the results of the numerical experiments . * notation * we use the standard notation . we denote by @xmath3 the time - axis . we denote by @xmath4 the @xmath5 identity matrix . we denote the set of natural numbers including zero by @xmath6 . the aim of the section is to define the class of piecewise - linear systems formally . a _ piecewise - linear system _ ( abbreviated as pwl ) is a dynamical system determined by @xmath7 here @xmath8 is the finite set of discrete modes , @xmath9 is the output space , @xmath10 is the state - space of the system in mode @xmath11 , for each @xmath12 @xmath13 , @xmath14 , @xmath15 , and @xmath16 are the parameters of the affine system in mode @xmath11 , @xmath17  is the set of initial states of the affine system in mode @xmath11 . the _ state space _ @xmath18 of @xmath19 is @xmath20 we call @xmath19 _ linear _ , if @xmath21 and @xmath22 for all @xmath11 , otherwise @xmath19 is called _ affine_.    we will use the following short - hand notation @xmath23    informally , the evolution of @xmath19 takes place as follows . as long as the value of the discrete state @xmath24 at time @xmath25 does not change , the continuous state and the continuous output an time @xmath25 change according to the affine system @xmath26 and @xmath27 . the discrete state can change at any time , however , we do not allow consecutive changes of discrete states immediately one after the other . if the discrete state changes to @xmath28 at time @xmath25 , then the new continuous state should satisfy @xmath29 , i.e. it should belong to the set of the designated initial states of the new discrete state @xmath28 . note that we do not specify the mechanism which triggers the change of discrete state . we also do not specify the initial discrete state , it is chosen by an unspecified mechanism . hence , the _ description above allows for several state- and output - trajectories . _    note that the definition of * pwl*s presented above differs from the one used in the literature . in the standard definition it is assumed that after a discrete - state transition , the new continuous state depends on the previous one . in contrast , here we only require that the new continuous state belongs to the designated set of initial states of the affine system associated with the new discrete state . this implies that _ the behavior of the system in discrete state @xmath30 depends only on the affine system associated with that discrete state and it does not depend on the discrete modes which were visited in the past . _ another consequence of the definition is that the model is not predictive . that is , the knowledge of the system parameters and the inputs does not uniquely determine the state- and output - trajectory of the system . in order to define the evolution of * pwl*s formally , we need to introduce the following notation and terminology .    a time event sequence is a strictly monotone sequence @xmath31 such that @xmath32 , @xmath33 and for all @xmath34 , @xmath35 . if @xmath36 then we require that @xmath37 . if @xmath38 , then by convention @xmath39 . the role of time event sequences is to formalize the time instances at which discrete events occur . the restrictions formulated in the definition imply that the set of switching times does not have an accumulation point , i.e. no zeno - behavior can take place . a state - trajectory of @xmath19 is a map @xmath40 such that there exists a time event sequence @xmath41 and a sequence of discrete modes @xmath42 such that for all @xmath43 , @xmath44 , it holds that for all @xmath45 , @xmath46 and @xmath47 the time event sequence @xmath48 is called the sequence of _ switching times _ of the state - trajectory @xmath49 . we _ denote by @xmath50 the set of all state - trajectories of @xmath19_.    an output - trajectory of @xmath19 is a map @xmath51 such that the following holds . there exists a state - trajectory @xmath49 of @xmath19 , such that @xmath52 for all @xmath53 . here , @xmath54 is the _ readout - map _ of @xmath19 , defined as @xmath55 we _ denote by @xmath56 the set of all output - trajectories of @xmath19_.    in the sequel , unless stated otherwise , _ @xmath57 denotes a function @xmath58_. the definition above implies that the external behavior of a * pwl *  is exactly a function of this type . the function @xmath57 is said to be realized by * pwl *  @xmath19 , if @xmath59 is an output - trajectory of @xmath19 , i.e. if @xmath60 . in this case @xmath19 is called a realization of @xmath57 . the dimension of a * pwl *  @xmath19 is a tuple @xmath61 , where @xmath62 if @xmath19 is linear , and @xmath63 , if @xmath19 is affine . that is , the first element of @xmath64 is the number of discrete states , the second element is the number of continuous state components . the additional dimension in the case of affine * pwl*s stems from the need to store the vectors @xmath65 , @xmath66 , @xmath11 . we use the following partial order relation on @xmath67 . we say that @xmath68 is smaller than or equal @xmath69 , denoted by @xmath70 , if @xmath71 and @xmath72 . note that the order relation @xmath73 in @xmath67 is indeed a partial order , it is not possible to compare all elements of @xmath67 . a * pwl *  @xmath19 is said to be a minimal realization of @xmath57 , if for any * pwl *  @xmath74 which is a realization of @xmath57 , @xmath75 . the definition above implies that if @xmath19 is minimal , then it is possible to compare its dimension to the dimension of any other realization of @xmath57 . since we work with a partial order on dimensions , the existence of a minimal system is not at all obvious . find conditions for existence of a * pwl *  realization of @xmath57 . characterize minimal * pwl*s  realizations of @xmath57 . find algorithms for computing a * pwl*s  realization from finite data . as it was indicated before , our motivation for studying the realization problem for * pwl*s is to lay the theoretical foundations for identification of * pwl*s . we present below the formulation of the identification problem . assume that the value of @xmath57 and its derivatives up to order @xmath76 , @xmath77 is measured at time instances @xmath78 . based on the ( possibly noisy ) data @xmath79 , find a * pwl *  realization of @xmath57 . our main motivation for studying identification of * pwl*s is the following . the problem of identifying interaction networks of complex dynamical systems can be reduced to the identification problem of * pwl*s . indeed , consider a partially observed non - linear system of the form @xmath80 where @xmath81 and @xmath82 are sufficiently smooth functions . we are interested in the dynamics of the network structure of . by _ the network structure of at time @xmath25 _ we mean the directed graph with nodes numbered by @xmath83 and where an edge goes from node @xmath0 to node @xmath1 if the @xmath0th component of @xmath84 influences the @xmath1th component of @xmath85 . as the values of @xmath84 changes , so does @xmath86 . hence the network structure described above depends on time . the basic question of interest is how to derive the dynamics of the interaction network based on the observed output @xmath87 . one obvious approach is to estimate @xmath57 and @xmath88 from @xmath87 , but this is feasible only if special restrictions are put on @xmath57 and @xmath88 . notice , however , that if we look at local linearizations of @xmath57 around an equilibrium and the interactions are strong enough , then @xmath89 influences @xmath90 , if the @xmath2th entry of the jacobian of @xmath57 is non - zero . prompted by the considerations above we propose the following approach to modelling and identification of the network structure of systems of the form . we assume that has finitely many equilibrium points . we assume that under influence of noise the system jumps between neighborhoods of these equilibrium points . we assume that the transition from the neighborhood of one equilibrium point to the neighborhood of another equilibrium point takes place very fast . hence , the behavior of the system during the transition can be ignored . we associate a * pwl *  @xmath19 of the form with as follows . let @xmath91 be the number of equilibrium points . for each @xmath11 , let @xmath92 be the corresponding equilibrium point and define @xmath93 here @xmath94 and @xmath95 denote the jacobian of @xmath57 and @xmath88 respectively . let @xmath96 be a suitably chosen neighborhood of @xmath92 . it then follows that if the state @xmath84 of is close enough to @xmath92 , we can approximate as follows @xmath97 as it was noted above , the relevant information regarding the interactions of state variables is already included in @xmath98 . hence , _ for the purposes of network reconstructions , identification of the * pwl*s approximation is sufficient_. below we present some basic results on realization theory of * pwl*s . first we present a characterization of minimality , after that we present conditions for existence of a realization and a realization algorithm . throughout the section , _ @xmath57 denotes a function @xmath58_.      in order to present our first result , we introduce the notion of a _ linear system with state - jumps_.    a linear system with state jumps , abbreviated as ( * lssj * ) is a linear * pwl *  @xmath19 of the form with one discrete state , i.e. @xmath99 . we will identify the * lssj *  @xmath19 with the collection of data @xmath100 , where @xmath101 and @xmath102 . the reason we call a * pwl *  with one discrete state a linear system with state jumps is that the system behaves as a linear system , with the exception that its state occasionally jumps back to one of the initial states . for * lssj*s , we can easily define the concepts of observability and span - reachability . a * lssj *  @xmath103 is called _ observable _ , if @xmath104 is an observable pair , and @xmath19 is called _ span - reachable _ , if @xmath105 . [ lssj : min ] note that @xmath106 is an @xmath107 invariant subspace containing @xmath108 , hence by restricting @xmath107 and @xmath109 to @xmath110 we can transform @xmath19 to a span - reachable * lssj *  @xmath111 while preserving output trajectories , i.e. @xmath112 . by applying linear observability reduction , we can transform the * lssj *  @xmath111 to a * lssj *  @xmath113 which is span - reachable , observable and has the same output trajectories as @xmath19 , i.e. @xmath114 .    [ phas : min ] assume that @xmath57 admits a realization by a * pwl*. then @xmath57 admits a minimal * pwl *  realization @xmath19 such that @xmath19 is a span - reachable and observable * lssj*. a * lssj *  realization of @xmath57 is minimal if and only if it is span - reachable and observable . the theorem above says that in general , the external behavior of any * pwl*s can be represented by a linear system with several initial states . that is , without further restrictions , the realization and identification problems for * pwl*s are equivalent to that of * lssjs*.    the proof of theorem [ phas : min ] relies on the following transformations .    [ pahs2linpahs ] define the _ linear _ * pwls *  @xmath115 associated with an affine @xmath19 as follows . @xmath116 where @xmath117 , @xmath118 , @xmath119 and @xmath120    the output trajectories of @xmath19 and @xmath115 coincide , i.e. @xmath121 , and @xmath122 . a state trajectory @xmath49 of @xmath115 is always a map of the form @xmath123 , @xmath53 such that @xmath124 , @xmath53 is a state - trajectory of @xmath19 . since @xmath125 , the statement of the proposition follows . let @xmath19 be a linear * pwl *  of the form . define the * lssj *  @xmath126 associated with @xmath19 as follows . @xmath127 where @xmath62 and @xmath128    the output trajectories of @xmath19 and @xmath126 are the same , i.e. @xmath129 , moreover , @xmath130 where @xmath131 .    for each @xmath11 , let @xmath132 be the subset of @xmath133 of the form @xmath134 , @xmath135 . for each @xmath136 of the above form , let @xmath137 . it then follows that @xmath138 is a state - trajectory of @xmath126 if and only if @xmath139 is a state - trajectory of @xmath19 , where @xmath140 for some @xmath11 if and only if @xmath141 . since @xmath142 for all @xmath143 , the statement of the proposition follows .    using the two transformations defined above , and remark [ lssj : min ] , we can present the proof of theorem [ phas : min ] . assume that @xmath19 is a * pwl *  realization of @xmath57 . if @xmath19 is affine , then replace @xmath19 with the associated linear * pwl *  @xmath115 , which is also a realization of @xmath57 . hence , we can assume that @xmath19 is already a linear * pwl*. construct then the * lssj *   @xmath126 associated with @xmath19 and apply remark [ lssj : min ] to @xmath126 to obtain a span - reachable and observable * lssj *  @xmath113 such that @xmath144 . it then follows that @xmath145 . hence , since @xmath19 was arbitrary , it is enough to look for minimal realizations among span - reachable and observable * lssj*s . notice that for any two * lssj*s  @xmath146 and @xmath147 , either @xmath148 or @xmath149 , hence , among all the possible * lssj*s  realizations of @xmath57 , there must exist a minimal one . the conditions for existence of a realization will be formulated using the rank of the hankel - matrix of @xmath57 . in order to define the hankel - matrix of @xmath57 , we have to define the notion of markov - parameters . to that end , we need the notion of _ piecewise - analytic functions_.    the map @xmath57 is called _ piecewise - analytic _ , if there exist a finite or infinite number of time instances , @xmath150 , @xmath151 , @xmath152 , @xmath44 for some @xmath153 , such that the following holds . if @xmath154 , let @xmath155 . then we require that @xmath33 and @xmath156 and for each @xmath44 , @xmath152 , @xmath57 is analytic on @xmath157 , but @xmath57 is not analytic on any neighborhood of @xmath158 in @xmath159 . we call the points @xmath160 the points of _ non - analyticity_. we define the set @xmath161 of indices of points of non - analyticity .    the intuition behind the definition is as follows . if @xmath57 has a realization by a * pwl *  @xmath19 , then the only points where @xmath57 is not analytic are the points where the corresponding state - trajectory of @xmath19 switches from one discrete mode to another . in fact , one can show that there always exists a state - trajectory of @xmath19 which yields @xmath57 as output trajectory and which switches only at time instances at which @xmath57 is not analytic . hence , if @xmath57 has a realization by a * pwl * , then @xmath57 is piecewise - analytic and the points of non - analyticity tell us the switching times of a * pwl *  realization of @xmath57 .    in the sequel , _ @xmath57 is assumed to be piecewise - analytic_.    assume that @xmath57 is piecewise - analytic and let @xmath160 be the points of non - analyticity of @xmath57 . for each @xmath162 , define the @xmath0th _ markov - parameter _ @xmath163 of @xmath57 as a sequence @xmath164 @xmath165    it is easy to see that the collection of markov - parameters @xmath166 determines the map @xmath57 uniquely . we use the markov - parameters to define the hankel - matrix of @xmath57 . we define the _ hankel - matrix _ @xmath167 of @xmath57 as the infinite matrix , rows of which are indexed by @xmath168 , and columns of which are indexed by @xmath169 . the entry of @xmath167 indexed by row index @xmath170 and by column index @xmath171 equals @xmath172_{(i , r),(j , l)}=(\mathbf{m}_{l}^f(i+j))_{r}\ ] ] where @xmath173 denotes the @xmath76th entry of markov - parameter @xmath174 . the rank of @xmath167 , denoted by @xmath175 , is the dimension of the linear space spanned by the columns of @xmath167 . [ theo : ex ] the map @xmath57 can be realized by a * pwl *  if and only if it is piecewise - analytic and @xmath176 . moreover , a minimal * lssj *  realization of @xmath57 of dimension @xmath177 can be constructed from @xmath167 .    * only if  * assume that @xmath57 can be realized by a * pwl*. then it can be realized by a * lssj *  @xmath103 . then @xmath178 @xmath179 for some @xmath180 , i.e. @xmath181 are the markov - parameters of the linear system @xmath104 from some initial condition . it then follows from linear systems theory that @xmath182 .    * if  * assume that @xmath183 and fix a basis in the linear span of the columns of @xmath167 . define @xmath109 as the matrix in this basis of the linear map which maps a column to the vector formed by the rows of that column indexed by @xmath184 , in this order . define the matrix @xmath107 as the matrix in this basis of the linear map which maps the column indexed by @xmath171 to the column indexed by @xmath185 . finally , let @xmath108 be the set of coordinate vectors of the columns of @xmath167 indexed by @xmath186 , @xmath187 . then @xmath188 is a * lssj *   realization of @xmath57 . if @xmath19 is another * lssj *  realization of @xmath57 , then from the * only if * part it follows that @xmath189 , i.e. @xmath190 is a minimal realization of @xmath57 .    in algorithm [ alg0 ] we present a kalman - ho - like realization algorithm for * lssj*s  and hence * pwl*s . to this end , we need the following definition . [ def : part ] fix integers + @xmath191 . define the set @xmath192 define the submatrix @xmath193 of @xmath167 as the matrix which is formed by the intersection of the rows of @xmath167 indexed by the elements of @xmath194 and the columns of @xmath167 indexed by the elements of @xmath195 . compute the decomposition @xmath196 such that @xmath197 and @xmath198 and @xmath199 . define @xmath200 as the matrix formed by the columns of @xmath201 indexed by elements of @xmath202 . let @xmath203 the matrix such that the column of @xmath204 indexed by @xmath171 equals the column of @xmath201 indexed by @xmath185 . define @xmath205 as follows . * the @xmath0th row of @xmath109 equals the row of @xmath206 indexed by @xmath207 , @xmath208 . * the matrix @xmath107 is the solution of the equation @xmath209 * the set @xmath108 consists of the columns of @xmath201 indexed by indices of the form @xmath186 , @xmath210 . [ alg0:theo ] if @xmath211 , and @xmath212 , then for some @xmath213 , the * lssj *  returned by algorithm [ alg0 ] is a minimal realization of @xmath57 . the conditions above hold if @xmath57 has a realization by a * lssj *  of dimension at most @xmath214 and with a finite set of initial states . the results on realization theory of * pwl*s indicate that the realization and identification problems for * pwl*s are in general ill - posed in the following sense . realizability by a * pwl*s is equivalent to realizability by a minimal * pwl*s with one single discrete state . hence , there seems to be no intrinsic way to choose discrete states based on the data . since the original intention was use the combination of the dynamics in various discrete modes to explain complex dynamics , the above conclusion is an unpleasant one . note that the main problem is that one can increase the dimension of the continuous state - space to encode discrete states . one way to remedy this is to place restriction on the dimension . to this end , we introduce the notion of @xmath215 realization . a * pwl *  @xmath19 is said to be a @xmath215 * pwl * if @xmath19 is of the form and @xmath216 and for all @xmath11 , @xmath217 . the map @xmath57 is said to admit a @xmath215 realization , if there exists a @xmath215 * pwl *  which is a realization of @xmath57 . the motivation for using @xmath215 realizations is that by choosing an appropriate @xmath218 , we can make sure that the growth in the number of continuous states can not be used to replace discrete states . below we present conditions for existence of a @xmath215 realization of @xmath57 . to this end , we have to introduce the following concept . consider the partitioning @xmath219 of the set of @xmath220 . for each @xmath221 , let @xmath222 be the sub - matrix of hankel - matrix @xmath167 which is formed by the columns of @xmath167 indexed by indices of the form @xmath171 , @xmath223 and @xmath224 . [ theo : pos1 ] the map @xmath57 has a realization by a linear @xmath215 * pwl * if and only if there exists a partitioning @xmath225 of @xmath226 with @xmath227 such that for all @xmath228 , @xmath229 . if the latter condition holds , then @xmath19 can be computed from @xmath167 with @xmath8 and with @xmath230 for all @xmath11 . theorem [ theo : pos1 ] above says that @xmath57 can be realized by a @xmath215 realization if the columns of the hankel - matrix @xmath167 can be divided in @xmath231 clusters such that for each fixed @xmath232 , all columns indexed by @xmath233 end up in the same cluster . moreover , the linear span of the elements of each cluster should have dimension at most @xmath218 .    in order to prove theorem [ theo : pos1 ] , we need the following transformation of a * lssj *  to a @xmath215 realization .    [ lss2pwls ] consider a partition @xmath234 of @xmath220 for some @xmath235 . consider an observable * lssj *  @xmath236 such that @xmath237 is a realization of @xmath57 . then there exists a unique state trajectory @xmath49 of @xmath237 , such that the switching times of @xmath49 coincide with the points of non - analyticity @xmath238 of @xmath57 and @xmath239 for all @xmath53 . define the * pwl *  realization @xmath240 of @xmath57 _ associated with @xmath237 and @xmath241 _ as follows . the * pwl *  @xmath240 is of the form , such that @xmath242 , and    * for all @xmath11 , let @xmath243 be the set of states @xmath244 such that for some @xmath224 , @xmath245 . define @xmath246 set @xmath247 and choose a basis of @xmath132 . * for each @xmath11 , define @xmath98 as the matrix in the basis of @xmath132 of the linear map @xmath248 obtained by restricting @xmath107 to @xmath132 . * let @xmath96 be the set of vector representations in the basis of @xmath132 of the elements of @xmath243 . * let @xmath249 be the matrix representation in the basis of @xmath132 of the linear map @xmath250 obtained by restricting @xmath109 to @xmath132 . if @xmath237 is a * lssj *  realization of @xmath57 , then the * pwl *  @xmath240 is a realization of @xmath57 .    * only if  * assume that @xmath19 is a * pwl *  realization of @xmath57 of the form , such that @xmath216 and @xmath217 for all @xmath11 . then there exists a state trajectory @xmath49 of @xmath19 such that @xmath251 and the switching times of @xmath49 contain the points of non - analyticity @xmath238 of @xmath57 . define @xmath252 . it then follows that @xmath253 is a partitioning of @xmath220 . moreover , for any @xmath254 , @xmath255 is the markov - parameter of the linear system @xmath256 from some initial condition . hence , by using linear systems theory we can show that @xmath257 . * if  * if the conditions of the theorem hold , then @xmath258 . construct the * lssj *  realization @xmath190 of @xmath57 , as in the proof of theorem [ theo : ex ] . apply definition [ lss2pwls ] to obtain a linear * pwl *  @xmath259 . then @xmath19 has @xmath227 discrete states . notice that the space @xmath132 from definition [ theo : pos1 ] is then isomorphic to the column space of @xmath222 , @xmath260 . hence @xmath261 . similarly to the general case , in algorithm [ alg05 ] we state an algorithm for computing a @xmath215 realization of @xmath57 .    apply algorithm [ alg0 ] to @xmath262 and denote the result by @xmath263 . [ alg05:clust ] compute a partitioning @xmath225 of @xmath264 such that @xmath227 and for all @xmath228 , @xmath265 . here @xmath266 is the sub - matrix of @xmath193 formed by the columns @xmath171 where @xmath267 and @xmath268 . return the * pwl *  @xmath269 . [ alg05:theo ] if @xmath57 has a @xmath215 realization of dimension at most @xmath270 , and and with a finite set of initial states , then for some @xmath213 , algorithm [ alg0 ] returns a minimal realization of @xmath57 .      below we show that @xmath57 is realizable by a * pwl *   if and only if @xmath57 satisfies a switched ar model .    a _ switched ar system _ ( abbreviated as * sars * ) is a tuple @xmath271 where @xmath8 , @xmath272 , @xmath273 and @xmath274 for all @xmath275 . the function @xmath57 is said to satisfy the * sars *  model @xmath276 , if the following holds . let @xmath238 be the points of non - analiticity of @xmath57 . for any @xmath162 , and for any @xmath277 , denote by @xmath278 the @xmath279th order right - hand derivative @xmath280 . then we require that for any @xmath281 there exist @xmath282 , such that @xmath283 , @xmath284 the set @xmath285 is called the set of discrete modes of @xmath276 . notice that if @xmath99 , then the * sars *  @xmath276 corresponds to an ar system , and @xmath57 satisfies @xmath276 if it satisfies an ar equation . note , however , that in contrast to trajectories of an ar systems , @xmath57 is not smooth . [ sars : theo ] the function @xmath57 has a realization by a * pwl *  with @xmath91 discrete states if and only if there exists a * sars *  @xmath276 with @xmath91 discrete modes such that @xmath57 satisfies @xmath276 .    * only if  * let @xmath19 be a realization of @xmath57 . without loss of generality we can assume that @xmath19 is linear , it is of the form and @xmath286 for all @xmath11 . then for each @xmath162 , there exists @xmath282 such that on @xmath157 , @xmath57 is an output trajectory of the linear system @xmath287 . from classical theory we then know that there exists matrices @xmath288 such that holds . * if  * assume @xmath57 satisfies a * sars *  @xmath276 of the form . define the linear * pwl *  @xmath289 of the form such that @xmath289 is a realization of @xmath57 as follows . the sets of discrete modes of @xmath289 and @xmath276 coincide , @xmath286 , @xmath290 and the matrices @xmath249 and @xmath291 , @xmath11 are the matrices of the linear system which corresponds to the ar @xmath292 . here we consider the linear system , state vector of which is the regressor @xmath293 . theorem [ phas : min ] and theorem [ sars : theo ] yield the following . the function @xmath57 has a realization by a * pwl *  if and only if it satisfies an ar system , i.e. for some @xmath294 , @xmath295 , @xmath296 below we present an algorithm which computes a realization with full observations based on measurements at finitely many time instances . more precisely , the algorithm solves the following problem . [ ident : full_obs ] fix integers @xmath272 and @xmath295 . consider a piecewise - analytic function @xmath297 and a finite sequence of time instances @xmath298 and assume that @xmath299 are known , i.e. we know that value of @xmath57 and its derivative at time instances @xmath300 . find a * pwl *  @xmath19 of the form with full - observations such that @xmath19 realizes @xmath57 , and @xmath301 , @xmath21 , and @xmath286 for all @xmath11 . * motivation of the identification problem * the motivation for considering the identification problem with full observations is the following .    1 . notice that the identification problem with full observations and the identification problem for * sars*s are equivalent . indeed , a * pwl *  with full observation can be considered as a * sars*s with @xmath302 . conversely , for a * sars *  @xmath276 , the associated * pwl *  @xmath289 from the proof of theorem [ sars : theo ] can be viewed as a * pwl *  with full observations , if the high - order derivatives of @xmath57 can be measured . hence , by theorem [ sars : theo ] , _ the identification problem for * pwl*s is equivalent to the identification problem for * pwl *  with full observations _ , if we assume that high - order derivatives of the output can be measured too . the problem of identification with full observations is still a non - trivial problem . even in this case we have problems with identifiability , see example [ example1 ] . 3 .   the solution of problem [ ident : full_obs ] enables us to prove experimentally the feasibility of approximating complex systems by * pwl*s . this is done by applying the identification algorithm which solves problem [ ident : full_obs ] to simulated trajectories of several well - known complex systems . [ example1 ] consider the * lssj*s @xmath303 and @xmath304 where @xmath305 is the identity matrix and @xmath306 and the remaining parameters are as follows ; @xmath307 and @xmath308 . it then follows that @xmath309 , @xmath53 can be realized both by @xmath146 and @xmath147 , but clearly @xmath310 . in other words , _ realizations of @xmath57 with full observation are not identifiable . _ the reason for this phenomenon is that the state of both @xmath146 and @xmath147 live in the subspace @xmath311 , @xmath312 , hence the difference in @xmath313 and @xmath314 is not visible . there is a wealth of results on identification of * sars * , mostly addressing the discrete - time siso case . many of these results , in particular , the algebraic approach @xcite could probably be adapted to the system class of this paper . investigating such adaptations remains a topic of future research . note that it is not at all clear that the convergence results for the discrete - time case have meaningful counterparts in the continuous - time case . in algorithm [ alg2 ] we present a solution for solving problem [ ident : full_obs ] . algorithm [ alg2 ] is an iteration consisting of the following steps . first , the algorithm starts with a random initialization of @xmath315 . subsequently , in each iteration step , the current estimates @xmath316 are updated by applying first the step [ alg : pwl-1 ] for estimating the updated weights with the fixed system parameters @xmath317 , and then applying step [ alg : pwl-2 ] with the previously updated weights @xmath318 for updating the estimates of the system parameters @xmath317 . the interpretation of the weights is as follows : @xmath319 is one if in @xmath158 the discrete state @xmath11 is active , and it is zero otherwise . the algorithm terminates when an absolute criterion @xmath320 falls below a pre - specified threshold @xmath321 . the iteration fails if after a pre - defined maximum number of iterations @xmath322 the criterion @xmath320 has not yet reached the lower threshold @xmath321 . unfortunately , the conditions under which the algorithm terminates and returns a correct realizations are not known yet . the criterion @xmath320 of the algorithm always decreases or stays the same during the iterations , and the same holds for the cost functions in step [ alg : pwl-1][alg : pwl-2 ] . hence , the criterion @xmath320 and the cost functions in step [ alg : pwl-1][alg : pwl-2 ] will converge to a fixed value . whether this value is a global optimum remains a topic of future research . note that the true parameters of a realization of @xmath57 render @xmath320 zero , but there might be many of them , see example [ example1 ] . initialize the estimates of @xmath323 , @xmath324 , @xmath325 . [ alg : pwl-1 ] solve the optimization for the weights @xmath326               \end{split}\ ] ] where for @xmath327 , @xmath328 , @xmath329 , @xmath11 , @xmath330    note that the optimization problem above is a linear programming problem , hence the optimal values of @xmath331 take values in the set @xmath332 . [ alg : pwl-2 ] using the weights @xmath331 recalculate the optimal values of @xmath333 , @xmath11 by solving the following minimization problem . @xmath334 the optimization problem is equivalent to the following linear least squares problem . collect the unknowns @xmath333 , @xmath11 into the matrix @xmath335 @xmath336 define the matrices @xmath337 and @xmath338 @xmath339 then the solution of can be read off from the entries of @xmath335 , where @xmath340    update @xmath341 using the new values . return @xmath19 of the form with @xmath286 , @xmath301 , @xmath21 and @xmath342 as computed in the last iteration , for all @xmath11 , and @xmath343 for all @xmath11 . we conducted several numerical experiments . they were all performed on a pc with intel pentium 4 cpu 2.8 ghz processor and 1.21 gb ram memory under windows xp , using matlab 7.5 including the optimization toolbox . in line with the definitions above , we use @xmath231 and @xmath218 to quantify the dimension of the * pwl *  of interest and @xmath344 to denote the number of data points . dynamical data @xmath299 was sampled from a given dynamical model @xmath345 at regular time intervals @xmath346 , @xmath347 for some @xmath348 . several initial states were chosen randomly , and a sampled trajectory was obtained for each of them . we then concatenated these sampled trajectories into one time series . notice that concatenation of finite components of two output trajectories of a * pwl *  is itself a finite component of a valid output trajectory . hence , the combined data can be viewed as originating from the measurements of one valid output trajectory of a * pwl*. gaussian noise with zero mean and variance @xmath349 was added to the obtained data . finally , algorithm [ alg2 ] was applied to the resulting data . in the numerical experiments we first tested the performance of the algorithm on artificial * pwl *  systems . furthermore , we studied the algorithm on lorenz systems and on the tyson - novak model for the cell cycle of budding yeast . we generate an artificial * pwl *  system containing @xmath231 affine subsystems , each of the same dimension @xmath350 . the state switching is obtained by partitioning the state space in @xmath231 subsets @xmath351 . the switching is generated by using the following switching law : the discrete mode @xmath0 is active , if the continuous state belongs to @xmath352 . in the absence of noise , and if the subsystems are sufficiently sampled , the estimation algorithm is usually able to find the system parameters , though occasionally it gets stuck in a local minimum . the amount of empirical data required to reconstruct the system parameters , within the bounds set by the noise , increases proportionally with @xmath218 and @xmath231 . figure [ n1 ] below shows the result of the algorithm for @xmath353 , @xmath354 , @xmath355 data points . the voronoi partitioning of state space of the original dataset ( black ) and the reconstructed dataset ( red ) are identical . the reconstructed system parameters and the original system parameters are identical as well .    ] further tests reveal that algorithm [ alg2 ] is able to find the correct parameters even if the regions @xmath351 are disconnected . when applying zero - mean gaussian noise @xmath356 , the accuracy of the reconstruction decreases as @xmath357 increases . figure [ n3 ] below shows the result on a dataset of dimension @xmath353 with @xmath354 subsystems and @xmath355 data points with a snr ( signal - to - noise ratio ) of @xmath358 . two subsystems could not be reconstructed , but the correlation between the three reconstructed systems with the best fitting original subsystems is @xmath359    ]      we applied algorithm [ alg2 ] to lorenz systems @xcite . this example provides insight how the identification algorithm performs on systems that are not piecewise affine by nature . the lorenz attractor can be visualized approximately as consisting of two linear subsystems ` left ' and ` right ' , with two different central attractors . a third region can be defined that involves the transition between the left and right subsystem . in figure [ n4 ] the application of algorithm [ alg2 ] to lorenz systems with @xmath360 data points is presented for @xmath361 . for @xmath362 the two main parts of the lorenz attractor are found . for @xmath363 also the transition region from the left to the right attractor is found . applying more than @xmath364 subsystems does not improve the performance . , title="fig : " ] , title="fig : " ]      we applied algorithm [ alg2 ] to data simulated using the tyson - novak model @xcite for the yeast cell cycle . the cell cycle is the succession of events whereby a cell grows and divides into two daughter cells . the simulated data was used to construct a * pwl *  approximation of the system using algorithm [ alg2 ] . this resulted in a clear partition of the cell cycle dynamics as a function of the selected number of classes @xmath231 . the best qualitative results were obtained for @xmath362 , see figure [ n6-a ] . moreover , for @xmath365 the results did not improve . it was found that the obtained two discrete modes correspond to the phases _ s _ ( dna synthesis ) and _ m _ ( mitosis ) of the cell cycle . the reconstruction remains valid under noise with snr at most @xmath366 .    ] we presented some basic results on realization theory of * pwl*s and a practical identification algorithm . analysis of the correctness and convergence of the presented algorithm remains a topic of future research . bako , l. , mercre , g. , vidal , r. , and lecoeuche . , s. ( 2009 ) . identification of switched linear state space models without minimum dwell time . in _ ifac symposium on system identification , saint malo , france_.                              ma , y. and vidal , r. ( 2005 ) . identification of deterministic switched arx systems via identification of algebraic varieties . in _ hybrid systems : computation and control _ , volume 3414 of _ lncs _ , 449  465 .
Abstract:

Answer: the paper address the realization and identification problem for a subclass of piecewise - affine hybrid systems . the paper provides necessary and sufficient conditions for existence of a realization , a characterization of minimality , and an identification algorithm for this subclass of hybrid systems . the considered system class and the identification problem are motivated by applications in systems biology .    ,    realization theory , minimization , identification , hybrid systems , network topology , gene - protein networks .