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: charge - parity - time ( @xmath1 ) symmetry is the most fundamental type of symmetry in quantum field theory 1742 - 6596 - 335 - 1 - 012011,aguilararevalo20131303 , where it holds for all relativistically invariant systems obeying the causality principle . its reduced form , _ viz_. , the @xmath0 symmetry , is almost exact too , save the small violation by weak nuclear forces @xcite . the @xmath2 operator is composed of three factors : parity transformation , @xmath3 , which reverses the coordinate axes ; charge conjugation , @xmath4 , which swaps particles and antiparticles ; and time reversal , @xmath5 . the proof of the @xmath1 and @xmath0 symmetries ( when the latter is relevant ) applies to hermitian hamiltonians ( @xmath6 ) , subject to the condition @xmath7 , which guarantees that the spectrum of the hamiltonian is real . however , one can not deduce from the @xmath1 or @xmath8 symmetry that the respective hamiltonian is necessarily hermitian @xcite . indeed , the consideration of hamiltonians which commute with a reduced symmetry operator , @xmath9 , demonstrates that they may contain an anti - hermitian ( dissipative ) part , provided that it is spatially antisymmetric ( odd ) , while the hermitian one is even @xcite . the spectrum of such a hamiltonian remains purely real up to a critical value of the strength of the anti - hermitian part , at which the @xmath9 symmetry is broken , making the system an essentially dissipative one ( recently , a model with _ unbreakable _ @xmath9 symmetry was found ; it includes defocusing cubic nonlinearity with the local strength growing from the center to periphery @xcite ) .    while in the quantum theory the possibility of the existence of non - hermitian @xmath9-symmetric hamiltonians is a purely theoretical one , such systems have been realized , theoretically muga , muga1,muga2,muga3,muga4,muga5,review and experimentally kip , feng05082011,citeulike:11031904,ruter:09,10.1117/12.807739,physreva.82.010103,2012natur.488 .. 163r , in optics , making use of the fact that the wave - propagation equation , derived in the standard paraxial approximation , is identical to the schrdinger equation in nonrelativistic quantum mechanics . in this context , the spatially symmetric and antisymmetric hermitian and anti - hermitian terms of the hamiltonian are represented , respectively , by even and odd distributions of the refractive index , and of the local gain - loss coefficient in the photonic medium . a @xmath9-symmetric electronic circuit was built too , following similar principles @xcite .    the essential role played by the kerr nonlinearity in optics has suggested the development of models in which the hamiltonian includes a quartic hermitian part too . the nonlinearity gives rise to families of @xmath9 -symmetric solitons , that were investigated in detail in continuous and discrete systems muga5,ptsolitons , ptsolitons1,ptsolitons2,ptsolitons3,ptsolitons4,ptsolitons5 , including @xmath9 -symmetric dual - core couplers couplers , couplers1,coupler - management . models combining the @xmath9 symmetry with quadratic nonlinearity in the dynamical equations ( i.e. , cubic terms in the respective hamiltonians ) were also elaborated vvk1,vvk2,vvk3 .    as a subject of quantum field theory , the @xmath1 and @xmath0 symmetries mainly relate to elementary particles 2012arxiv1201.1594s , srednicki2007quantum , kursunogammalu2013confluence,2013phlb .. 718.1500c . on the other hand , the above - mentioned works on the implementation of non - hermitian @xmath9-symmetric hamiltonians in photonics suggest looking for a possibility to design optical settings that would realize non - hermitian hamiltonians featuring the full @xmath1 symmetry , as well as its @xmath0 reduction . a possibility to implement the former symmetry was recently explored in ref . @xcite , which addressed not optics , but rather a two - component bose - einstein condensate with the spin - orbit coupling between the components , one of which is subject to the action of loss , and the other one is supported by gain . in terms of optics systems , the symmetry of that models is similar to the @xmath9 symmetry in a dual - core waveguide , with a combination of continuous @xmath10  transformation acting in the longitudinal direction , and another @xmath3  transformation which swaps the two cores . a similar symmetry was proposed in ref . @xcite , which put forward a @xmath9-symmetric coupler subject to the action of management " , in the form of periodic simultaneous switch of the signs of the coupling and gain - loss coefficients . the present work aims to offer emulation of the @xmath0 symmetry in a two - component optical system , which , at the phenomenological level , may be considered as a dual - core waveguide with opposite signs of the group - velocity dispersion ( gvd ) in the cores and a phase - velocity mismatch between them , embedded into an active medium . we demonstrate that the system can be derived , without phenomenological assumptions , as a model of the spatial - domain propagation for two fundamental - frequency ( ff ) components with orthogonal polarizations of light , pumped by an undepleted second - harmonic ( sh ) wave in a birefringent medium with the @xmath11 nonlinearity . we  further investigate conditions for persistence and breaking of the @xmath0 symmetry , both analytically and numerically . in particular , the addition of cubic ( kerr , alias @xmath12 ) terms to the system with the active coupling breaks the symmetry at the nonlinear level , but helps to stabilize confined breather states ( oscillatory solitons ) , and gives rise to a family of stable stationary gap solitons . the paper is organized as follows : section 2 introduces the model , in its both forms ( phenomenological and the one based on the @xmath11 interaction ) and reports analytical results . numerical findings for the linear and nonlinear systems are presented in section 3 , and section 4 concludes the paper . at the phenomenological level , we consider the copropagation of optical modes @xmath13 and @xmath14 in a dual - core coupler with opposite gvd signs in the cores ( cf . refs . @xcite , where a similar feature was introduced in different contexts , and also refs . @xcite , where systems with opposite signs of group velocities were considered in the contexts of coupled right- and left - handed waveguides ) , and a phase - velocity mismatch @xmath15 between them . the linear coupling between the cores is provided by cross - gain terms , with strength @xmath16 , which is possible when the coupler is embedded into an active medium , as recently proposed in ref . the model is represented by the following system of propagation equations , which include the kerr nonlinearity too , with respective coefficient @xmath17 ( all the quantities are dimensionless ) : @xmath18here @xmath19 is the propagation distance , @xmath20 is the reduced time agrawal2001nonlinear , the gvd coefficients are scaled to be @xmath21 , and @xmath22 may be transformed to @xmath23 by changing @xmath24 . positive and negative values of @xmath17 can also be transformed into each other by substitution @xmath25 , therefore , in what follows below we consider only @xmath26 . then , rescaling allows one to fix @xmath27 , but we prefer to keep it as a free parameter , the variation of which helps to monitor a transition from the weakly nonlinear system to a strongly nonlinear one . it is relevant to mention that a dissipative discrete system with opposite signs of the discrete dispersion and a wavenumber mismatch between the components was introduced in ref . however , that model included the dissipative coefficient in a single equation , therefore it did not realize the symmetry considered here .    equations ( [ u ] ) and ( [ v ] ) can be derived from the non - hermitian ( complex ) lagrangian , which is usual for @xmath9-symmetric systems:@xmath28 dt   \notag \\ -i\gamma \int_{-\infty } ^{+\infty } \left ( u^{\ast } v+uv^{\ast } \right ) dt , \label{l}\end{gathered}\]]which generates the respective non - hermitian hamiltonian in an obvious way.the total energy , @xmath29 dt\equiv e_{u}(z)+e_{v}(z ) ,   \label{p}\]]is not conserved by eqs . ( [ u ] ) and ( [ v ] ) . instead , the system gives rise to the following energy - balance equations:@xmath30the equality of @xmath31 and @xmath32 , i.e. , the conservation of @xmath33 , means that the linear coupling of the present type causes mutual amplification or attenuation of both components .    the linear version of eqs . ( [ u ] ) and ( [ v ] ) , with @xmath34 , are invariant with respect to the @xmath0 transformation , defined as@xmath35where the swap of @xmath13 and @xmath14 stands for @xmath3 , and the complex conjugation  for @xmath4 ( conserved @xmath36 may be considered as the respective charge ) . it is relevant to compare the system of eqs . ( u ) , ( [ v ] ) and their invariance transformation ( [ cpt ] ) with the previously studied model of the @xmath9-symmetric coupler , which was based on the following equations couplers , couplers1,coupler - management : @xmath37obviously , eqs . ( [ upt ] ) and ( [ vpt ] ) are invariant with respect to transformation @xmath38 , which , in the present context , may be considered as corresponding to the @xmath39 symmetry , the reversal of @xmath19 playing the role of additional @xmath5 . it is relevant too to compare the present model to the system of equations with opposite gvd terms , coupled by the usual conservative terms , rather than by those representing the gain and loss @xcite : @xmath40where @xmath41 is a real coupling constant . the linear version of this system is invariant with respect to the _ anti_-@xmath0 transformation : @xmath42 anti " corresponding to the relative sign flip , cf . ( [ cpt ] ) . the nonlinearity breaks the symmetry of system ( [ u ] ) , ( [ v ] ) , as the opposite relative signs of the gvd and cubic terms in the two equations make it impossible to swap @xmath13 and @xmath14 , which represents the @xmath3 transformation in eq . ( [ cpt ] ) . nevertheless , nonlinear effects are obviously interesting too . it is demonstrated below that the nonlinearity creates solitons in the present system . in this connection , it is relevant to mention recently introduced nonlinear models with alternating gain and loss , which do not obey the condition of the @xmath9 symmetry , but nevertheless support stable solitons @xcite .    a solution to the linear version of eqs . ( [ u]-[v ] ) in the form of plane waves , @xmath43 produces a dispersion relation for the wavenumber and frequency:@xmath44obviously , in the case of @xmath45 the spectrum given by eq . ( [ dr2 ] ) is pure real , provided that@xmath46while in the case of @xmath47 the spectrum always includes an imaginary component . the change of the spectrum from real to a partly imaginary one , with the increase of the gain - loss coefficient , at @xmath48 ( provided that @xmath45 ) implies the breakup of the @xmath0 symmetry , similar to the phase transition which is the generic feature of @xmath9-symmetric systems @xcite . if condition ( [ thr ] ) holds , the spectrum given by eq . ( [ dr2 ] ) features a _ bandgap_,@xmath49      while the system of eqs . ( [ u ] ) , ( [ v ] ) was introduced above phenomenologically , its linear version can be derived , in the spatial domain ( rather than in the temporal one ) , starting from the fundamental propagation model for two ff and one sh  components of light waves , @xmath13 , @xmath50 and @xmath51 , respectively , in the dissipation - free medium with the type - ii @xmath52 interaction @xcite:@xmath53where @xmath54 is the transverse coordinate , and @xmath55 is an ff - sh wavenumber mismatch . then , adopting the usual approximation for parametric down - conversion , of undepleted sh pump , we replace it by a constant , @xmath56 , neglecting eq . ( [ chi2w ] ) , denote @xmath57 , and apply the complex conjugation to eq . ( [ chi2v]):@xmath58these equations differ from the linear version of eqs . ( [ u ] ) , ( [ v ] ) only by the replacement of @xmath20 by @xmath54 .    as for cubic terms , they can be added to eqs . ( [ chi2u]-[chi2w ] ) as ones accounting for the kerr nonlinearity in the @xmath11 waveguide . however , in terms of eqs . ( [ uphys ] ) and ( [ vphys ] ) , the resulting cubic terms will be different from those adopted in eqs . ( [ u ] ) and ( v ) , as the above - mentioned complex conjugation of eq . ( [ chi2v ] ) will produce the cubic term in eq . ( [ vphys ] ) with the sign opposite to that in eq . ( [ v ] ) [ incidentally , in this case the cubic terms do not break the @xmath0 invariance of eqs . ( [ u ] ) and ( [ v ] ) ] . furthermore , because the original components , @xmath13 and @xmath50 , correspond to two orthogonal polarizations of the ff wave , the nonlinear extension of eqs . ( [ uphys ] ) and ( [ vphys ] ) should also include the respective xpm ( cross - phase - modulation ) terms , _ viz_. , @xmath59 and @xmath60 , respectively , assuming that the four - wave mixing terms may be neglected , as usual , due to sufficiently strong birefringence @xcite . in the present work , we focus on the nonlinear terms adopted in eqs . ( [ u ] ) , ( [ v ] ) , while those corresponding to the derivation for the @xmath11 system will be considered elsewhere . the system based on eqs . ( [ u ] ) , ( [ v ] ) can be investigated in an analytical form for broad small - amplitude pulses , with widths ( @xmath61 ) and amplitudes satisfying conditions @xmath62 in this case , the linearized version of the system yields , in the lowest approximation , two different relations between the field components : one solution has@xmath63and another one features@xmath64where @xmath65 and @xmath66 are slowly varying amplitudes , in comparison with @xmath67 . the substitution of these expressions into lagrangian ( [ l ] ) leads to a real effective lagrangians for the slowly varying functions ( its imaginary part cancels out in the present approximation ) , which , in turn , give rise to either one of the two following nonlinear schrdinger ( nls ) equations for the slow evolution:@xmath68where @xmath69 and @xmath70 pertain to @xmath71 and @xmath72 , respectively , while the effective gvd and nonlinearity coefficients are @xmath73note that eq . ( [ slow ] ) implies that the system conserves the total energy in the present approximation , which complies with the fact that the substitution of relations ( [ v - u ] ) and ( [ u - v ] ) into energy - balance equations ( [ dp / dz2 ] ) yields @xmath74 . fundamental solutions to the linear version of eq . ( [ slow ] ) are well known in the form of expanding gaussians ( coherent states , in terms of quantum mechanics ) @xcite:@xmath75where @xmath76 is an initial width , and @xmath77 are arbitrary amplitudes . this result is drastically different from that obtained for broad pulses in the linear version of eq . ( [ k ] ) with the conservative coupling @xcite . the latter system gives rise to an effective equation for a slowly varying function in the form of a single linear schrdinger equation with the gvd term periodically ( in @xmath19 ) changing its sign , thus generating robust oscillating gaussian pulses , rather than expanding ones ( [ gauss ] ) . further , the full nonlinear equations ( [ slow ] ) for @xmath71 and @xmath78 give rise to commonly known solutions for bright and dark solitons , respectively . in particular , the bright nls solitons with an arbitrary ( small ) propagation constant , @xmath79 , are @xmath80    it is also worthy to note that the equations for @xmath71 and @xmath72 are , severally , galilean invariant ones , i.e. , the linear and nonlinear versions of the equation for @xmath71 give rise , severally , to moving gaussians and bright solitons , while the full underlying system ( [ u ] ) , ( [ v ] ) does not feature the galilean invariance .      inside of bandgap ( [ gap ] ) , it is natural to look for stationary gap - soliton solutions @xcite of eqs . ( [ u ] ) , ( [ v ] ) in the form of @xmath81for which eqs . ( [ u ] ) and ( [ v ] ) reduce to a system of stationary equations,@xmath82(recall @xmath83 is fixed ) . unlike the usual coupled - mode system for bragg gratings @xcite , there is no substitution which could reduce eqs . ( [ u ] ) and ( [ v ] ) to a single equation , therefore the system should be solved numerically , in the general case . the stability of the gap solitons should then be tested numerically too . as concerns broad solitons ( [ sol ] ) , they actually correspond to the gap solitons at values of @xmath84 close to the bottom of bandgap ( [ gap ] ) , @xmath85 , cf . a similar relation between the general gap solitons and broad ones in the standard model of nonlinear bragg gratings @xcite . a gap - soliton solution to eqs . ( [ u ] ) and ( [ v ] ) with strong asymmetry between the two components can be found in an approximate form for the limit case of weak gain - loss coupling , @xmath86 , at @xmath87 , i.e. , exactly at the central point of bandgap ( [ gap ] ) . the zero - order approximation ( for @xmath88 ) is@xmath89for @xmath90 ( in this approximation , @xmath87 is not required ) . then , the first - order correction is determined by the linearized equation for @xmath91:@xmath92an exact closed - form solution to eq . ( [ 1 ] ) can be found , by means of the fourier transform , at the center of the bandgap , i.e. , at @xmath87 ( a similar solution was reported , in another context , in ref . \right\ } . \label{vexact}\]]an exact solution to eq . ( [ 1 ] ) is available too in the special case of @xmath94:@xmath95 -}\arctan { \left [ \exp \left ( { % 2}\sqrt{\frac{q}{5}}t\right ) \right ] } \right\ } . \label{3/5}\end{gathered}\]]these solutions exist under exactly the same condition , @xmath45 , which was adopted above . although it may not be immediately obvious , both solutions ( [ vexact ] ) and ( [ 3/5 ] ) are even functions of @xmath20 , exponentially decaying at @xmath96 . these approximate analytical solutions are compared with their numerically found counterparts below , see figs . fig : c_gamma_0.1 and [ fig : c_gamma_0.1_k_3_5 ] . equations ( [ u]-[v ] ) were solved numerically by means of the split - step fourier - transform method @xcite , for four different sets of initial conditions . two of them were taken with a gaussian pulse in either component : @xmath97and two other initial sets are given by eqs . ( [ v - u]-[u - v ] ) and ( [ gauss ] ) with @xmath98 . the interval for the temporal variable was fixed as @xmath99 , to ensure that reflections from its boundaries did not affect the results , moreover all quantities plotted are dimensionless . we start the analysis with the case of weak coupling , @xmath100 . a typical numerical solution for this case , presented in fig . fig : l_amp_gamma_0.1 , shows the expansion of the gaussian launched in the form of initial conditions ( [ ic1 ] ) . the @xmath14 component remains weak as the coupling constant is small and , accordingly , the total energy remains very close to the initial value . a detailed comparison with eqs . ( [ v - u]-[slow ] ) demonstrates that the asymptotic stage of the evolution , for broad pulses , is accurately predicted by the analytical approximation . next , we consider the situation close to the @xmath1-symmetry - breaking threshold ( [ thr ] ) , namely , with @xmath101 for @xmath102 . the respective numerical solution , generated by initial conditions ( [ ic1 ] ) , is displayed in fig . [ fig : l_amp_gamma_0.9 ] , which , naturally , demonstrates strong coupling between the two components and more dramatic evolution . in this case too , the asymptotic stage of the evolution for broad pulses is correctly predicted by the above - mentioned analytical approximation .    to test the symmetry of the system , we have also performed simulations of the evolution starting from initial conditions ( [ ic2 ] ) , with swapped components @xmath13 and @xmath14 . comparison of the respective results , displayed in figs . [ fig : v_amp_gamma_0.9](a ) and [ fig : v_amp_gamma_0.9](b ) , with their counterparts shown above in fig . [ fig : l_amp_gamma_0.9 ] confirms the symmetry . furthermore , the detailed comparison of the real and imaginary parts of the two components in both cases ( not shown here in detail ) exactly corroborates the full @xmath0 symmetry implied by definition ( cpt ) .    for @xmath103 , when the the @xmath0 symmetry of the system is broken , according to eq . ( [ thr ] ) , direct simulations ( not shown here ) demonstrate blowup of solutions , as should be expected above the symmetry - breaking point @xcite .    the analytical approximation for broad pulses , based on eqs . ( [ v - u ] ) and ( [ gauss ] ) , was directly tested by comparing its predictions with the numerically simulated evolution commencing from the initial conditions produced by eqs . ( [ v - u ] ) and ( [ gauss ] ) with @xmath98 and the upper sign in the latter equation . figure [ fig : u_numeric_gamma_0.1 ] shows that the respective analytical and the numerical results are almost identical . the comparison produces equally good results ( not shown here in detail ) if the initial conditions are taken , instead , as per eq . ( [ u - v ] ) and eq . ( [ gauss ] ) with the lower sign , at @xmath98 . on the other hand , for strong coupling , e.g. , for @xmath101 , when @xmath104 is no longer a rapidly oscillating carrier in comparison with slowly varying @xmath71 and @xmath78 , see eqs . ( [ v - u ] ) and ( [ u - v ] ) , the analytical approximation is no longer relevant . the comparison with the numerical results corroborates this expectation ( not shown here in detail either ) . simulations of the nonlinear system , based on eqs . ( [ u ] ) and ( [ v ] ) , were performed by varying the nonlinearity coefficient , @xmath17 , and ( as above ) the coupling coefficient , @xmath16 . the initial conditions were taken in the form of eq . ( [ ic1 ] ) , unless stated otherwise . we start by addressing the weakly coupled system with weak nonlinearity , _ viz_. , the one with @xmath105 and @xmath106 . for @xmath107 and @xmath108 , figs . [ fig : amp_nl_0.1_gamma_0.1](a ) and ( b ) demonstrate that the focusing nonlinearity readily causes self - trapping of a robust oscillating quasi - soliton . thus , the weak nonlinearity , while breaking the @xmath0 symmetry ( see above ) , creates the self - confined modes . next , we increase the strength of the coupling to @xmath109 , keeping the nonlinear term small , with @xmath108 . in this case , fig : amp_nl_0.1_gamma_0.8(a ) and ( b ) demonstrate strong self - focusing of the modes occurring around @xmath110 , followed by the propagation of the confined mode in a sufficiently robust form , although with more conspicuous emission of radiation waves than in the case of @xmath107 , cf . fig : amp_nl_0.1_gamma_0.1(a ) . thus , in this case too , the system tends to form oscillatory quasi - soliton modes .    the formation of these solitons is readily explained by eq . ( [ slow ] ) for @xmath71 . indeed , it is easy to check that the width and amplitude of the emerging solitons satisfy conditions ( [ broad ] ) . the solitons are observed in figs . [ fig : amp_nl_0.1_gamma_0.1 ] and fig : amp_nl_0.1_gamma_0.8 in an oscillatory form , which is different from the stationary solution ( [ sol ] ) , in accordance with the well - knows fact that perturbed nls solitons may feature long - lived vibrations , similar to those observed in these figures @xcite .    for @xmath111 and the same weak nonlinearity , @xmath108 , the simulations demonstrate that amplitudes of both modes , @xmath13 and @xmath14 , diverge after a short propagation distance , which implies that the symmetry breaking takes place in these cases , that are close to the threshold ( [ thr ] ) of the symmetry breaking . the small difference of @xmath101 from the exact linear threshold , @xmath112 , is compensated in this case by the presence of the nonlinearity which , as said above , is also a @xmath0-symmetry - breaking factor . the swap of initial conditions ( [ ic1 ] ) and ( [ ic2 ] ) in the nonlinear system produces a strong effect . indeed , in the above simulations , performed for input ( [ ic1 ] ) , the pulse was launched in component @xmath13 , where the nonlinearity is self - focusing [ see eq . ( [ u ] ) ] , while initial conditions ( [ ic2 ] ) imply that the pulse is launched into component @xmath14 with the self - defocusing cubic term , see eq . ( [ v ] ) . accordingly , in the latter case , the simulations produce the results displayed in fig . missing_v_amp_nl_0.1_gamma_0.8 : instead of the quick self - trapping ( cf . fig . [ fig : amp_nl_0.1_gamma_0.8 ] ) , the pulse launched in the @xmath14 component features slow expansion . an additional difference is that the frequency of oscillations observed in the latter case is approximately half of that observed in fig . [ fig : amp_nl_0.1_gamma_0.8 ] . the increase of @xmath17 at a fixed coupling constant , @xmath16 , enhances the @xmath0-symmetry - breaking effects , and eventually leads to destruction of the quasi - soliton . in particular , for the weakly coupled system , with @xmath107 , the destabilization of the quasi - soliton sets in at critical value @xmath113 , as shown in fig . fig : amp_nl_0.1_gamma_0.1 . in this case , the integral energy slowly grows with @xmath19 , which is followed by blowup at very large values of @xmath19 ( not shown here in detail ) . ( a ) , but at the critical value of the nonlinearity strength , @xmath114 , at which the gradual destruction of the quasi - soliton commences.,scaledwidth=50.0% ]    by means of systematical simulations , we have collected the critical values of @xmath17 , at which the quasi - soliton suffers the onset of the destabilization , eventually leading to the blowup , at increasing values of the coupling constant , @xmath16 . the corresponding dependence between @xmath115 and @xmath16 , shown in fig . ( [ fig : stable_nl_vs_gamma ] ) , naturally demonstrates that the critical strength of the nonlinearity vanishes when @xmath16 approaches the threshold of the symmetry breaking in the linear system , @xmath112 , see eq . ( [ thr ] ) [ recall the normalization is fixed by setting @xmath116 in eqs . ( [ u ] ) and ( [ v ] ) ] .    ) , in the plane of @xmath117 . , scaledwidth=50.0% ]    the analytical approximation based on eqs . ( [ v - u ] ) and ( [ sol ] ) was tested for the broad solitons too . figure [ fig : u_numeric_gamma_0.1_nl ] shows that the respective analytical and the numerical results are almost identical , thus validating the analytical approximation for the nonlinear system . the quasi - solitons considered above are built as breathers , featuring permanent oscillations in both components . on the other hand , eqs . ( [ 0]-[vexact ] ) predict the existence of stationary gap solitons in the same system . to check this possibility in the numerical form , we solved eqs . ( u-[v ] ) by means of the newton s method @xcite . this was done using the approximate analytical solution , given by eqs . ( [ 0 ] ) and ( [ vexact ] ) , as the initial guess . the results are produced here for @xmath118 and different values of the coupling constant , @xmath16 .    in the case of weak coupling , @xmath107 , when it is natural to expect the solutions to be strongly asymmetric , in terms of the two components , the numerical solution at @xmath87 , i.e. , at the center of the bandgap , is very close to its above - mentioned analytical counterpart , as seen in fig . fig : c_gamma_0.1 .     and @xmath119 ( note different scales in the panels ) , obtained from eqs . ( [ u ] ) and ( [ v ] ) for @xmath116 , @xmath120 , @xmath121 , and @xmath87 . the respective analytical solution is given by eq . ( 0 ) and ( [ vexact ] ) . , scaledwidth=50.0% ]    for larger values of @xmath16 , the numerical solution differs from the analytical one , obtained under condition @xmath122 , although the difference remains relatively small for @xmath123 , as shown in fig . [ fig : c_gamma_0.5 ] . the difference becomes essential for @xmath101 [ in fact , very close to the symmetry - breaking threshold ( [ thr ] ) ] , as can be seen in fig . [ fig : c_gamma_0.9 ] .    , but for @xmath124.,scaledwidth=50.0% ]    , but for @xmath125.,scaledwidth=50.0% ]    at @xmath126 , the numerically found solutions are still close to the the strongly asymmetric analytical expressions given by eqs . ( [ 0 ] ) and ( vexact ) for @xmath87 , provided that @xmath127 is small enough , see fig . fig : c_gamma_0.2_k_0.2 for @xmath128 and @xmath129 . however , at larger @xmath84 , such as @xmath130 with the same @xmath128 [ note that @xmath130 falls into the bandgap ( [ gap ] ) in this case ] , the numerical solution for the @xmath91 component strongly deviates from the analytical expression given by eq . ( v ) for @xmath87 , while the @xmath131 component is still close to the simple analytical form ( [ 0 ] ) , see fig . [ fig : c_gamma_0.2_k_0.8 ] .    , but for @xmath132 , and the numerical solution taken for @xmath129 [ recall that analytical expression ( [ vexact ] ) pertains to @xmath87 ] . , scaledwidth=50.0% ]    , but for the numerical solution taken for @xmath130.,scaledwidth=50.0% ]    an analytical solution for strongly asymmetric gap solitons ( corresponding to @xmath134 ) was also obtained above in the form of eqs . ( [ 0 ] ) and ( [ 3/5 ] ) for @xmath135 . for @xmath107 , this solution is virtually identical to its numerically found counterpart , as shown in fig . fig : c_gamma_0.1_k_3_5 .     and @xmath119 , obtained from eqs . ( [ u ] ) and ( [ v ] ) for @xmath116 , @xmath120 , @xmath136 , and @xmath135 . the respective analytical solution is given by eqs . ( [ 0 ] ) and ( [ 3/5 ] ) . , scaledwidth=50.0% ]    finally , the stability of the numerically generated gap solitons was tested by using them as initial conditions in direct simulations of eqs . ( [ u ] ) and ( [ v ] ) . the results , not shown here in detail , demonstrate that all the tested examples of gap solitons are stable , both at the center of the gap , @xmath87 , and off the center , including values of the coupling constant ( such as @xmath101 ) which are close to the symmetry - breaking threshold ( [ thr ] ) . the objective of this work is to introduce a model of a dual - core waveguide which may implement an optical system featuring the @xmath0 symmetry . essential ingredients of the model are opposite gvd signs in the two cores , a phase - velocity mismatch between them , and the linear coupling of the gain - loss type , which is possible if the waveguide is embedded into an active medium , or may be provided by the propagation in the @xmath11 medium of the type - ii ( three - wave ) type , neglecting the depletion of the sh ( second - harmonic ) pump . nonlinear cubic terms , which destroy the symmetry , were considered as well ( in the case of the @xmath11 medium , they are different from those considered here ) . it is predicted in an approximate analytical form and demonstrated numerically that the linear system gives rise to expanding gaussian pulses . relatively weak nonlinearity produces an essential effect , building broad oscillatory quasi - solitons , which are destroyed in direct simulations if the nonlinearity is too strong . further , the analysis predicts a general family of stationary gap solitons in the nonlinear system , that have been also found and checked for the stability in the numerical form , the broad solitary pulses being a limit case of the gap solitons near the bottom edge of the bandgap . the analysis can be continued by considering higher - order modes [ it is well known that linear schrdinger equations ( [ slow ] ) give rise to higher modes in the form of hermite - gauss wave functions ] and interactions between fundamental solitons in the nonlinear version of the system . on the other hand , it was mentioned above that eq . ( [ slow ] ) for @xmath72 suggests the existence of dark solitons in the present system , which is an interesting issue too . still another possibility is to consider the different form of the cubic nonlinearity , corresponding to the underlying @xmath137 system . a challenging extension is to consider a two - dimensional version of the model , which may be based on a dual - core planar waveguide embedded into the active medium .
Abstract:

Answer: we introduce a model of a dual - core optical waveguide with opposite signs of the group - velocity - dispersion ( gvd ) in the two cores , and a phase - velocity mismatch between them . the coupler is embedded into an active host medium , which provides for the linear coupling of a gain - loss type between the two cores . the same system can be derived , without phenomenological assumptions , by considering the three - wave propagation in a medium with the quadratic nonlinearity , provided that the depletion of the second - harmonic pump is negligible . this linear system offers an optical realization of the charge - parity ( @xmath0 ) symmetry , while the addition of the intra - core cubic nonlinearity breaks the symmetry . by means of direct simulations and analytical approximations , it is demonstrated that the linear system generates expanding gaussian states , while the nonlinear one gives rise to broad oscillating solitons , as well as a general family of stable stationary gap solitons .


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: networks of balls and springs or frames of nodes connected by compressible struts provide realistic models for physical systems from bridges to condensed solids . their elastic properties depend on their coordination number @xmath0  the average number of nodes each node is connected to . if @xmath0 is large enough , the networks are elastic solids whose long - wavelength mechanical properties are described by a continuum elastic energy with non - vanishing elastic moduli . if @xmath0 is small enough , the networks have deformation modes of zero energy  they are floppy . as @xmath0 is increased from the floppy side , a critical value , @xmath4 , is reached at which springs provide just enough constraints that the system has no zero - energy  floppy " modes @xcite ( or mechanisms @xcite in the engineering literature ) , and the system is _ isostatic_. the phenomenon of rigidity percolation @xcite whereby a sample spanning rigid cluster develops upon the addition of springs is one version of this floppy - to - rigid transition . the coordination numbers of whole classes of systems , including engineering structures @xcite ( bridges and buildings ) , randomly packed spheres near jamming @xcite , network glasses @xcite , cristobalites @xcite , zeolites @xcite , and biopolymer networks @xcite are close enough to @xmath4 that their elasticity and mode structure is strongly influenced by those of the isostatic lattice . ( a ) section of a kagome lattice with @xmath5 and @xmath6 three - site unit cells . nearest - neighbor bonds , occupied by harmonic springs , are of length @xmath7 . the rotated row ( second row from the top ) represents a floppy mode . next - nearest neighbor bonds are shown as dotted lines in the lower left hexagon . the vectors @xmath8 , @xmath9 , and @xmath10 indicate symmetry directions of the lattice . the numbers in the triangles indicate those that twist together under pbcs in zero modes along the three symmetry direction . note that there are only @xmath11 of these modes . ( b ) section of a square lattice depicting a floppy mode in which all sites along a line are displaced uniformly . ( c ) twisted kagome lattice , with lattice constant @xmath12 , derived from the undistorted lattice by rigidly rotating triangles through an angle @xmath13 . a unit cell , bounded by dashed lines , is shown in violet . arrows depict site displacements for the zone - center , i.e. , zero wavenumber , @xmath14 mode which has zero ( nonzero ) frequency under free ( periodic ) boundary conditions . sites @xmath15 , @xmath16 , and @xmath17 undergo no collective rotation about their center of mass whereas sites @xmath15 , @xmath18 , and @xmath19 do . ( d ) superposed snapshots of the twisted lattice showing decreasing areas with increasing @xmath13 . ] though the isostatic point always separates rigid from floppy behavior , the properties of isostatic lattices are not universal ; rather they depend on lattice architecture . here we explore the the unusual properties of a particular class of periodic isostatic lattices derived from the two - dimensional kagome lattice by rigidly rotating triangles through an angle @xmath13 without changing bond lengths as shown in fig . [ fig : periodic - isostatic ] . the bulk modulus @xmath20 of these lattices is rigorously zero for all @xmath21 . as a result , their poisson ratio acquires its limit value of @xmath22 ; when stretched in one direction , they _ expand _ by an equal amount in the orthogonal direction : they are maximally auxetic @xcite . these modes represent collapse pathways @xcite of the kagome lattice . modes of isostatic systems are generally very sensitive to boundary conditions @xcite , but the degree of sensitivity depends on the details of lattice structure . for reasons we will discuss more fully below , modes of the square lattice , which is isostatic , are in fact insensitive to changes from free boundary conditions ( fbcs ) to periodic boundary conditions ( pbcs ) , whereas those of the undistorted kagome lattice are only mildly so . the modes of both , however , change significantly when rigid boundary conditions ( rgbs ) are applied . we show here that in all families of the twisted kagome lattice , modes depend sensitively on whether fbcs , pcbs or rgbs are applied : finite lattices with free boundaries have floppy surface modes that are not present in their periodic or rigid spectrum or in that of finite undistorted kagome lattices . in the long wavelength limit , the surface floppy modes , which are present in any @xmath23 material with @xmath24 , reduce to surface rayleigh waves @xcite described by a conformally invariant energy whose analytic eigenfunctions are fully determined by boundary conditions . at shorter wavelengths , the surface waves become sensitive to lattice structure and remain confined to within a distance of the surface that diverges as the undistorted kagome lattice is approached . in the simplest twisted kagome lattice , all floppy modes are surface modes , but in more complicated lattices , including ones with uniaxial symmetry , we construct , there are both surface and bulk floppy modes .    arguments due to j.c . maxwell @xcite provide a criterion for network stability : networks in @xmath1 dimensions consisting of @xmath25 nodes , each connected with central - force springs to an average of @xmath0 neighbors , have @xmath26 zero - energy modes when @xmath3 ( in the absence of redundant bonds - see below ) . of these a number , @xmath27 , which depends on boundary conditions , are trivial rigid translations and rotations , and the and the remainder are floppy modes of internal structural rearrangement . under fbcs an pbcs , @xmath28 equals @xmath29 and @xmath1 , respectively . with increasing @xmath0 , mechanical stability is reached at the isostatic point at which @xmath30 . the maxwell argument is a global one ; it does not provide information about the nature of the floppy modes and does not distinguish between bulk or surface modes . the kagome lattice of central force springs shown in fig . [ fig : periodic - isostatic](a ) is one of many locally isostatic lattices , including the familiar square lattice lattice in two dimensions [ fig . [ fig : periodic - isostatic](b ) ] and the cubic and pyrochlore lattices in three dimensions , with exactly @xmath2 nearest - neighbor ( @xmath31 ) bonds connected to each site not at a boundary . under pbcs , there are no boundaries , and every site has exactly @xmath23 neighbors . finite , @xmath25-site sections of these lattices have surface sites with fewer than @xmath32 neighbors and of order @xmath33 zero modes . the free kagome lattice with @xmath34 and @xmath35 unit cells along its sides [ fig .  [ fig : periodic - isostatic](a ) ] has @xmath36 sites , @xmath37 bonds , and @xmath38 zero modes , all but three of which are floppy modes . these modes , depicted in fig . [ fig : periodic - isostatic](a ) , consist of coordinated counter rotations of pairs of triangles along the symmetry axes @xmath8 , @xmath9 and @xmath10 of the lattice . there are @xmath34 modes associated with lines parallel to @xmath8 , @xmath35 associated with lines parallel to @xmath10 , and @xmath39 modes associated with lines parallel to @xmath9 .    in spite of the large number of floppy modes in the kagome lattice , its longitudinal and shear lam coefficients , @xmath40 and @xmath41 , and its bulk modulus @xmath42 are nonzero and proportional to the nearest neighbor ( @xmath31 ) spring constant @xmath43 : @xmath44 , and @xmath45 . the zero modes of this lattice can be used to generate an infinite number of distorted lattices with unstretched springs and thus zero energy @xcite . we consider only periodic lattices , the simplest of which are the twisted kagome lattices obtained by rotating triangles of the kagome unit cell through an angle @xmath13 as shown in figs . [ fig : periodic - isostatic](c ) and ( d ) @xcite . these lattices have @xmath46 rather than @xmath47 symmetry and , like the undistorted kagome lattice , three sites per unit cell . as fig . [ fig : periodic - isostatic](d ) shows , the lattice constant of these lattices is @xmath12 , and their area @xmath48 decreases as @xmath49 as @xmath13 increases . the maximum value that @xmath13 can achieve without bond crossings is @xmath50 so that the maximum relative area change is @xmath51 . since all springs maintain their rest length , there is no energy cost for changing @xmath13 , and as a result , @xmath20 is zero for every @xmath52 , whereas the shear modulus @xmath53 remains nonzero and unchanged . thus , the poisson ratio @xmath54 attains its smallest possible value of @xmath22 . for any @xmath21 , the addition of next - nearest - neighbor ( @xmath55 ) springs , with spring constant @xmath56 ( or of bending forces between springs ) stabilizes zero - frequency modes and increases @xmath20 and @xmath57 . nevertheless , for sufficiently small @xmath56 , @xmath57 remains negative . figure [ fig : auxetic - phased ] shows the region in the @xmath58 plane with negative @xmath57 . phase diagram in the @xmath59 plane showing region with negative poisson ratio @xmath57 . ] we now turn to the linearized phonon spectrum of the kagome and twisted kagome lattices subjected to pbcs . these conditions require displacements at opposite ends of the sample to be identical and thus prohibit distortions of the shape and size of the unit cell and rotations but not uniform translations , leaving two rather than three trivial zero modes . the spectrum @xcite of the three lowest frequency modes along symmetry directions of the undistorted kagome lattice with and without @xmath55 springs is shown in fig . [ fig : kagome - spectrum](a ) . when @xmath60 , there is a floppy mode for each wavenumber @xmath61 running along the entire length of the three symmetry - equivalent straight lines running from @xmath62 to @xmath63 to @xmath62 in the brillouin zone [ see inset to fig . [ fig : kagome - spectrum ] ] . when @xmath64 , there are exactly @xmath65 wavenumbers with @xmath61 along each of these lines for a total of @xmath66 floppy modes . in addition , there are three zero modes at @xmath67 corresponding to two rigid translations , and one floppy mode that changes unit cell area at second but not first order in displacements , yielding a total of @xmath68 zero modes rather than the @xmath69 modes expected from the maxwell count under fbcs . this is our first indication of the importance of boundary conditions . the addition of @xmath55 springs endows the floppy modes at @xmath60 with a characteristic frequency @xmath70 and causes them to hybridize with the acoustic phonon modes [ fig . [ fig : kagome - spectrum](a ) ] @xcite . the result is an isotropic phonon spectrum up to wavenumber @xmath71 and gaps at @xmath63 and @xmath62 of order @xmath72 . remarkably , at nonzero @xmath13 and @xmath60 , the mode structure is almost identical to that at @xmath73 and @xmath74 with characteristic frequency @xmath75 and length @xmath76 . in other words , twisting the kagome lattice through an angle @xmath13 has essentially the same effect on the spectrum as adding @xmath55 springs with spring constant @xmath77 . thus under pbcs , the twisted kagome lattice has no zero modes other than the trivial ones : it is  collectively " jammed in the language of references @xcite , but because it is not rigid with respect to changing the unit cell size , it is not strictly jammed . ( a)phonon spectrum for the undistorted kagome lattice . dashed lines depict frequencies at @xmath60 and full lines at @xmath74 . the inset shows the brillouin zone with symmetry points @xmath63 , @xmath62 , and @xmath78 . note the line of zero modes along @xmath79 when @xmath60 , all of which develop nonzero frequencies for wavenumber @xmath80 when @xmath74 reaching @xmath81 on a plateau beginning at @xmath82 defining a length scale @xmath83 . ( b ) phonon spectrum for @xmath84 and @xmath60 . the plateau along @xmath79 defines @xmath85 and its onset at @xmath86 defines a length @xmath87 . ] to understand the origin of the differences in the zero - mode count for different boundary conditions , we turn to an elegant formulation @xcite of the maxwell rule that takes into account the existence of redundant bonds ( i.e. , bonds whose removal does not increase the number of floppy modes @xcite ) and states in which springs can be under states of self - stress . consider a ring network in two dimensions shown in fig .  [ fig : self - stress ] with @xmath88 nodes and @xmath89 springs with three springs of length @xmath7 and one spring of length @xmath90 . the maxwell count yields @xmath91 zero modes : two rigid translations , one rigid rotation , and one internal floppy mode  all of which are  finite - amplitude " modes with zero energy even for finite - amplitude displacements . when @xmath92 , the maxwell rule breaks down . in the zero - energy configuration , the long spring and the three short ones are colinear , and a prestressed state in which the @xmath90-spring is under compression and the three @xmath7-springs are under tension ( or vice versa ) but the total force on each node remains zero becomes possible . this is called a state of self - stress . the system still has three finite amplitude zero modes corresponding to arbitrary rigid translations and rotations , but the finite - amplitude floppy mode has disappeared . in the absence of prestress , it is replaced by two  infinitesimal " floppy modes of displacements of the two internal nodes perpendicular of the now linear network . in the presence of prestress , these two modes have a frequency proportional to the square root of the tension in the springs . thus , the system now has one state of self stress and one extra zero mode in the absence of prestress , implying @xmath93 , where @xmath94 is the number of states of self stress . this simple count is more generally valid as can be shown with the aid of the equilibrium and compatibility matrices @xcite , denoted , respectively , as @xmath95 and @xmath96 . @xmath95 relates the vector @xmath97 of @xmath98 spring tensions to the vector @xmath99 of @xmath100 forces at nodes via @xmath101 , and @xmath102 relates the the vector @xmath103 of @xmath100 node displacements to the vector @xmath104 of @xmath98 spring stretches via @xmath105 . the dynamical matrix determining the phonon spectrum is @xmath106 . vectors @xmath107 in the null space of @xmath95 , ( @xmath108 ) , describe states of self - stress whereas vectors @xmath109 in the null space of @xmath102 represent displacements with no stretch @xmath104 , i.e. , modes of zero energy . thus the nullspace dimensions of @xmath95 and @xmath102 are , respectively , @xmath94 and @xmath110 . the rank - nullity theorem of linear algebra @xcite states that the rank @xmath111 of a matrix plus the dimension of its null space equals its column number . since the rank of a matrix and its transpose are equal , the @xmath95 and @xmath102 matrices , respectively , yield the relations @xmath112 and @xmath113 , implying @xmath114 . under pbcs , locally isostatic lattices have @xmath2 exactly , and the maxwell rule yields @xmath115 : there should be no zero modes at all . but we have just seen that both the square and undistorted kagome lattices under pbcs have of order @xmath33 zero modes as calculated from the dynamical matrix , which , because it is derived from a harmonic theory , does not distinguish between infinitesimal and finite - amplitude zero modes . thus , in order for there to be zero modes , there must be states of self - stress , in fact one state of self - stress for each zero mode .    in the square lattice under fbcs , @xmath116 and @xmath117 , there are no states of self stress , and @xmath118 zero modes depicted in fig . [ fig : periodic - isostatic](b ) . under pbcs , the dimension of the nullspace of @xmath95 is @xmath119 , and there are also @xmath120 zero modes that are identical to those under fbcs . we have already seen that there are @xmath121 zero modes in the free undistorted kagome lattice . direct evaluations @xcite ( see text s1 ) of the dimension of the null spaces of @xmath95 and @xmath102 for the undistorted kagome lattice with pbcs yields @xmath122 when @xmath64 . the zero modes under pbcs are identical to those under fbcs except that the @xmath123 modes associated with lines parallel to @xmath9 under fbcs get reduced to @xmath34 modes because of the identification of apposite sides of the lattice required by the pbcs as shown in fig . [ fig : periodic - isostatic](a ) . thus the modes of both the square and kagome lattices do not depend strongly on whether fcbs or pbcs are applied . under rbc s , however , the floppy modes of both disappear . the situation for the twisted kagome lattice is different . there are still @xmath124 zero modes under fbcs , but there are only two states of self stress under pbcs and thus only @xmath125 zero modes , as a direct evaluation of the null spaces of @xmath95 and @xmath102 verifies ( see appendix for details ) , in agreement with the results obtained via direct evaluation of the eigenvalues of the dynamical matrix @xcite . all of the floppy modes under fbcs have disappeared . ( a ) ring - network with @xmath126 showing internal floppy mode . ( b ) ring - network with @xmath92 showing one of the two infinitesimal modes . ] an effective long - wavelength energy @xmath127 for the low - energy acoustic phonons and nearly floppy distortions provides insight into the nature of the modes of the twisted kagome lattice . the variables in this theory are the vector displacement field @xmath128 of nodes at undistorted positions @xmath129 and the scalar field @xmath130 describing nearly floppy distortions within a unit cell . the detailed form of @xmath127 depends on which three lattice sites are assigned to a unit cell . figure [ fig : periodic - isostatic](c ) depicts the lattice distortion @xmath14 for the nearly floppy mode at @xmath63 ( with energy proportional to @xmath131 ) along with a particular representation of a unit cell , consisting of a central asymmetric hexagon and two equilateral triangles , with @xmath132 sites on its boundary . if sites @xmath15 , @xmath16 , and @xmath17 are assigned to the unit cell , then the distortion @xmath14 involves no rotations of these sites relative to their center of mass , and the harmonic limit of @xmath127 depends only on the symmetrized and linearized strain @xmath133 and on @xmath14 : @xmath134 , \label{eq : full - elastic}\end{aligned}\ ] ] where @xmath135 is the symmetric - traceless stain tensor , @xmath136 , @xmath137 , @xmath138 , @xmath139 , and @xmath140 . the last term in which @xmath141 is a third - rank tensor , whose only non - vanishing components are @xmath142 , is invariant under operations of the group @xmath47 but not under arbitrary rotations . the @xmath143 term is the only one that reflects the @xmath46 ( rather than @xmath47 ) symmetry of the lattice . there are several comments to make about this energy . the gauge - like coupling in which the isotropic strain @xmath144 appears only in the combination @xmath145 guarantees that the bulk modulus vanishes : @xmath14 will simply relax to @xmath146 to reduce to zero the energy of any state with nonvanishing @xmath144 . the coefficient @xmath78 can be calculated directly from the observation that under @xmath14 alone , the length of every spring changes by @xmath147 , and this length change is reversed by a homogenous volume change @xmath148 . in the @xmath149 limit , @xmath150 , and the energy reduces to that of an isotropic solid with bulk modulus @xmath151 if the @xmath152 and @xmath153 terms , which are higher order in gradients , are ignored . the @xmath153 term gives rise to a term , singular in gradients of @xmath154 , when @xmath14 is integrated out that is responsible for the deviations of the finite - wavenumber elastic energy from isotropy . at small @xmath13 , the length scale @xmath155 appears in several places in this energy : in the length @xmath156 and in the ratios @xmath157 , @xmath158 , and @xmath159 . at length scales much larger than @xmath155 , the @xmath152 and @xmath153 terms can be ignored , and @xmath14 relaxes to @xmath146 leaving only the shear elastic energy of an elastic solid proportional @xmath160 . at length scales shorter than @xmath155 , @xmath14 deviates from @xmath146 and contributes significantly to the form of the energy spectrum . if @xmath15 , @xmath18 and @xmath19 in fig . [ fig : periodic - isostatic](d ) are assigned to the unit cell , then @xmath14 involves rotations relative to the lattice axes , and the energy develops a cosserat - like form @xcite that is a function of @xmath161 rather than @xmath14 .    the modes of our elastic energy in the long - wavelength limit ( @xmath162 ) are simply those of an elastic medium with @xmath24 . in this limit , there are transverse and longitudinal bulk sound modes with equal sound velocities @xmath163 and @xmath164 , where @xmath165 is the particle mass at each node and @xmath166 is the mass density . in addition there are rayleigh surface waves @xcite in which there is a single decay length ( rather than the two at @xmath167 ) , and displacements are proportional to @xmath168 $ ] with @xmath169 for a semi - infinite sample in the right half plane so that the penetration depth into the interior is @xmath170 . these waves have zero frequency in two dimensions when @xmath24 , and they do not appear in the spectrum with pbcs . thus this simple continuum limit provides us with an explanation for the difference between the spectrum of the free and periodic twisted kagome lattices . under fbcs , there are zero - frequency surface modes not present under pbcs .    further insight into how boundary conditions affect spectrum follows from the observation that the continuum elastic theory with @xmath24 depends only on @xmath171 . the metric tensor @xmath172 of the distorted lattice is related to the strain @xmath173 via the simple relation @xmath174 ; and @xmath175/2 $ ] , which is zero for @xmath176 , is invariant , and thus remains equal to zero , under conformal transformations that take the metric tensor from its reference form @xmath177 to @xmath178 for any continuous function @xmath179 . the zero modes of the theory thus correspond simply to conformal transformations , which in two dimensions are best represented by the complex position and displacement variables @xmath180 and @xmath181 . all conformal transformations are described by an analytic displacement field @xmath182 . since by cauchy s theorem , analytic functions in the interior of a domain are determined entirely by their values on the domain s boundary ( the  holographic " property @xcite ) , the zero modes of a given sample are simply those analytic functions that satisfy its boundary conditions . for example , a disc with fixed edges ( @xmath183 ) has no zero modes because the only analytic function satisfying this fbc is the trivial one @xmath184 ; but a disc with free edges ( stress and thus strain equal to zero ) has one zero mode for each of the analytic functions @xmath185 for integer @xmath186 . the boundary conditions @xmath187 and @xmath188 on a semi - infinite cylinder with axis along @xmath189 are satisfied by the function @xmath190 when @xmath191 , where @xmath192 is an integer . this solution is identical to that for classical rayleigh waves on the same cylinder . like the rayleigh theory , the conformal theory puts no restriction on the value of @xmath192 ( or equivalently @xmath193 ) . both theories break down , however , at @xmath194 beyond which the full lattice theory , which yields a complex value of @xmath195 , is needed . figure [ fig : ky - kx](a ) shows an example of a surface wave . at the bottom of this figure , @xmath196 is an almost perfect sinusoid . as @xmath197 decreases toward the surface , the amplitude grows , and in this picture reaches the nonlinear regime by the time the surface at @xmath198 is reached . figure [ fig : ky - kx](b ) plots @xmath199 as a function of @xmath193 obtained both by direct numerical evaluation and by an analytic transfer matrix procedure @xcite for different values of @xmath13 ( text s1 ) . the rayleigh limit @xmath200 is reached for all @xmath13 as @xmath201 . interestingly the rayleigh limit remains a good approximation up to values of @xmath193 that increase with increasing @xmath13 . the inset to fig . [ fig : ky - kx ] , plots @xmath202 as a function of @xmath203 and shows that in the limit @xmath204 ( @xmath205 ) , @xmath199 obeys an @xmath13-independent scaling law of the form @xmath206 . the full complex @xmath207 obeys a similar equation . this type of behavior is familiar in critical phenomena where scaling occurs when correlation lengths become much larger than microscopic lengths . the function @xmath208 approaches @xmath209 as @xmath210 and asymptotes to @xmath211 for @xmath212 . thus for @xmath213 , @xmath214 and for @xmath215 , @xmath216 . as @xmath13 increases , @xmath217 is no longer much larger than one , and deviations from the scaling law result . the situation for surfaces along different directions ( e.g. , along @xmath218 rather than @xmath198 ) is more complicated and will be treated in a future publication @xcite . ( a)lattice distortions for a surface wave on a cylinder , showing exponential decay of the surface displacements into the bulk . this figure was constructed by specifying a small sinusoidal modulation on the bottom boundary and propagating lattice - site positions upward to the free boundary at the top under the constraint of of constant lengths and periodic boundary conditions around the cylinder . distortions near and at the top boundary , which have become nonlinear , are not described by our linearized treatment . ( b ) @xmath219 as a function of @xmath220 for lattice rayleigh surface waves for @xmath221 , in order from bottom to top . smooth curves are the analytic results from a transfer matrix calculation , and dots are from direct numerical calculations . the dashed line is the continuum rayleigh limit @xmath200 . curves at smaller @xmath13 break away from this curve at smaller values of @xmath207 than do those at large @xmath13 . at @xmath222 , @xmath223 diverges at @xmath224 . the inset plots @xmath202 as a function of @xmath225 for different @xmath13 . the lower curve in the inset(black ) is the @xmath13-independent scaling function of @xmath226 reached in the @xmath227 limit . the other curves from top to bottom are for @xmath228 , @xmath229 , @xmath230 , and @xmath231 ( chosen to best present results rather than to match the curves in the main figure ) . curves for @xmath232 are essentially indistinguishable from the scaling limit . the curve at @xmath233 stops because @xmath234 . ] the @xmath46 twisted kagome lattice is the simplest of many lattices that can be formed from the kagome and other periodic isostatic lattices . figures [ fig : otherauxetic ] ( a ) and ( b ) show two other examples of isostatic lattices constructed from the kagome lattice . most intriguing is the lattice with pgg symmetry . its geometry has uniaxial symmetry , yet its long - wavelength elastic energy is identical to that of the @xmath46 twisted kagome lattice , i.e , it is isotropic with a vanishing bulk modulus , and its mode structure near @xmath67 is isotropic as shown in fig .  [ fig : otherauxetic ] ( c ) . thus , this system loses long - wavelength zero - frequency bulk modes of the undistorted kagome lattice to surface modes . however , at large wavenumber , lattice anisotropy becomes apparent , and ( infinitesimal ) floppy bulk modes appear . thus in this and related systems , a fraction of the zero modes of the under fbcs are bulk modes that are visible under pbcs , and a fraction are surface modes that are not . randomly packed spheres above the jamming transition with average coordination number @xmath235 exhibit a characteristic frequency @xmath236 and length @xmath237 and a transition from a debye - like ( @xmath238 ) to a flat density of states at @xmath239 @xcite . the square and kagome lattices with randomly added @xmath55 springs have the same properties @xcite . a general  cutting " argument @xcite provides a procedure for perturbing away from the isostatic limit and an explanation for these properties . however , it only applies provided a finite fraction of the of order @xmath240 floppy modes of a sample with sides of length @xmath241 cut from an isostatic lattice with pbcs are extended , i.e. , have wave functions that extend across the sample rather than remaining localized either in the interior or at the surface the sample . clearly the twisted kagome lattice , whose floppy modes are all surface modes , violates this criterion ; and indeed , the density of states of the lattice with @xmath242 shows debye - behavior crossing over to a flat plateau at @xmath243 . adding next nearest neighbor bonds gives rise to a length @xmath244 and crossover to the plateau at @xmath245 . the pgg lattice in fig . [ fig : otherauxetic](a ) , however , has both extended and surface floppy modes , so its crossover to the a flat plateau occurs at @xmath246 rather than at @xmath247 or @xmath248 . ( a ) kagome - based lattice with pgg space group symmetry and uniaxial @xmath249 point group symmetry . ( b ) lattice with p6 space symmetry but global @xmath250 point - group symmetry . ( c ) density plot of the spectrum of the lowest frequency branch of the pgg uniaxial kagome lattice . the spectrum is absolutely isotropic near the origin point @xmath63 , but it has a zero modes on two symmetry related continuous curves at large values of wave - number . ] our study highlights the rich and remarkable variety of physical properties that isostatic systems can exhibit . under fbcs , floppy modes can adopt a variety of forms , from all being extended to all being localized near surfaces to a mixture of the two . under pbc s , the presence of floppy modes depends on whether the lattice can or can not support states of self stress . when a lattice exhibits a large number of zero - energy edge modes , its mechanical / dynamical properties become extremely sensitive to boundary conditions , much as do the electronic properties of the topological states of matter studied in    quantum systems @xcite . the zero - energy edge modes observed in our isostatic lattices are collective modes whose amplitudes decay exponentially from the edge with a finite decay length , in direct contrast to the very localized and trivial floppy modes arising from dangling bonds . we focussed primarily on high - symmetry lattices derived from the kagome lattice , but the properties they exhibit , namely a deficit of floppy modes in the bulk and the existence of floppy surface modes , are shared by any two - dimensional system with a vanishing bulk modulus ( or the equivalent in anisotropic systems ) . three - dimensional analogs of the twisted kagome lattice @xcite can be constructed by rotating tetrahedra in pyrochlore and zeolite lattices @xcite and in cristobalites @xcite . these lattices are anisotropic . with @xmath31 forces only , they exhibit a vanishing modulus for compression applied in particular planes rather than isotropically , but we expect them to exhibit many of the properties the two - dimensional lattices exhibit . finally , we note that maxwell s ideas can be applied to spin systems such as the heisenberg anti - ferromagnet on the kagome lattice @xcite , and the possibility of unusual edge states in them is intriguing . we are grateful for informative conversations with randall kamien , andea liu , and s.d . guest . this work was supported in part by nsf under grants dmr- 0804900 ( t.c.l . and x.m . ) , mrsec dmr-0520020 ( t.c.l . and a.s . ) and jqi - nsf - pfc ( k. s. ) . in this section , we will provide some of the details for calculating the equilibrium and compatibility matrices @xmath95 and @xmath102 and the dynamical matrix @xmath251 under periodic boundary conditions . the equilibrium matrix @xmath95 relates the @xmath252 dimensional vector of bond tensions @xmath97 to the @xmath100 dimensional vector @xmath99 of forces at the nodes via @xmath101 . we label unit cells by the their position vectors @xmath253 , where @xmath254 and @xmath255 are integers and @xmath256 are the primitive lattice vectors . ( here we take the lattice spacing to be one , i.e. , @xmath257 . ) each unit cell has three sites , labeled @xmath15 , @xmath16 , and @xmath17 as shown in fig . [ fig : lattice1 ] , so we label the @xmath258 sites by @xmath259 , where @xmath260 and @xmath261 is the number of unit cells . forces at the nodes are denoted by @xmath262 and their @xmath263 components as @xmath264 , where @xmath265 . they have to be balanced by stretching forces in the springs located on the bonds . there are six bonds per cell , which we can take to be the six bonds in the distorted hexagon shown in fig . [ fig : lattice1 ] , oriented parallel to the unit vectors @xmath266 for @xmath267 , where @xmath13 is the twist angle of the twisted kagome lattice . when @xmath73 , these vector reduce to the edge vectors of a symmetric hexagon with @xmath268 . the @xmath263 bonds , labeled with @xmath269 , are occupied by springs under tension @xmath270 that exert pulling ( for @xmath271 ) forces , direct along the bond vectors @xmath272 , on nodes . we define vectors @xmath273 . because we have periodic boundary conditions , we can express both @xmath274 and @xmath275 in terms of their fourier components , @xmath276 the equilibrium matrix is diagonal in @xmath277 , so it breaks up into independent @xmath278 blocks for each of the @xmath6 independent vectors @xmath277 . for simplicity , we consider only the case @xmath64 for which @xmath279 and for which the number of sites is @xmath280 . there are four bonds incident on each site . the equations relating @xmath262 to @xmath275 follow from fig . [ fig : lattice1 ] : @xmath281 we map @xmath282 to a single index @xmath165 with @xmath283 and @xmath284 and then fourier transform . the result is that @xmath95 breaks up into independent @xmath278 blocks @xmath285 for each @xmath277 . this leads to @xmath286 where @xmath287 the null space of this matrix is easily calculated with the aid of mathematica .    * case 1 , @xmath73 * : this is the untwisted kagome lattice . when @xmath61 , the null space of @xmath288 is empty unless @xmath277 is perpendicular to one of the primitive lattice vectors : @xmath289 and @xmath290 , in which case there is a single vector , @xmath291 where @xmath292 is the number of cells , for each @xmath192 and @xmath293 in the null space . when @xmath67 , there are three vectors , @xmath294 , @xmath295 , and @xmath296 in the null space . these are the @xmath297 limits of @xmath298 , @xmath299 , and @xmath300 , respectively . we can , therefore , take the set of null - space vectors of @xmath95 to be @xmath298 , @xmath299 , and @xmath299 for the @xmath34 values of @xmath301 , which include @xmath297 . linear combinations of these vectors are also null - space vectors , and we can use them to construct vectors confined to a single line of bonds . for example the vector @xmath302 corresponds to a state in which springs on all @xmath15-bond in cells with centers at @xmath303 for any @xmath254 and on all @xmath11-bonds in cells with centers at @xmath304 are all under the same tension . this is a state of self - stress in which all bonds along a particular straight line parallel to the @xmath8-axis are under tension . similar states of self - stress for any of the @xmath305 lines of bonds can be constructed .    * case 2 , @xmath306 * : this is the twisted lattice . the null space of @xmath288 is empty for _ all _ @xmath307 , and it contains only two vectors when @xmath308 . these correspond to two states of self - stress in which the stress the stress on the six bonds take on both positive and negative values that are identical in all unit cells . the compatibility matrix @xmath102 relates the @xmath100 displacements @xmath103 to the @xmath309 stretches @xmath310 . the displacements are labeled the same way as the forces , i.e. as @xmath311 and the stretches in the same way as the bond tensions , i.e. as @xmath312 . the equations relating @xmath311 and @xmath313 are @xmath314 fourier transforming and using the same @xmath315 notation as for the equilibrium matrix , we obtain @xmath316 where @xmath317 as required by the constraint @xmath318 . finally , the dynamical matrix is @xmath319 , whose components are @xmath320 \nonumber\\ d_{12}({\mathbf{q } } ) & = & = -\frac{1}{2 } \sqrt{3 } \cos ( 2 \alpha ) \nonumber\\ d_{13}({\mathbf{q } } ) & = & = d_{31}^*({\mathbf{q } } ) = - e^{-i \left(\frac{q_x}{2}+\frac{\sqrt{3 } q_y}{2}\right ) } \sin ^2\left(\frac{\pi } { 6}-\alpha \right ) -e^{-i q_x}\sin ^2\left(\alpha + \frac{\pi}{6}\right ) \nonumber\\ d_{14}({\mathbf{q } } ) & = & d_{41}^*({\mathbf{q } } ) = \frac{1}{2 } e^{-i q_x}\cos \left(\frac{\pi } { 6}-2 \alpha \right ) + \frac{1}{2}e^{-i \left(\frac{q_x}{2}+\frac{\sqrt{3 } q_y}{2}\right ) } \cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{15}({\mathbf{q } } ) & = & d_{51}^*({\mathbf{q } } ) = -2e^{-\frac{1}{2 } i \sqrt{3 } q_y } \cos ^2(\alpha ) \cos \left(\frac{q_x}{2}\right ) \nonumber\\ d_{16}({\mathbf{q } } ) & = & d_{61}^*({\mathbf{q } } ) = i e^{-\frac{1}{2 } i \sqrt{3 } q_y } \sin ( 2 \alpha ) \sin \left(\frac{q_x}{2}\right ) \nonumber\\ d_{22}({\mathbf{q } } ) & = & 2 \sin ^2(\alpha ) + \cos ^2\left(\frac{\pi } { 6}-\alpha \right ) + \cos ^2\left(\alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{23}({\mathbf{q } } ) & = & d_{32}^*({\mathbf{q } } ) = \frac{1}{2}e^{-i q_x } \cos \left(\frac{\pi } { 6}-2 \alpha \right ) + \frac{1}{2 } e^{-i\left(\frac{q_x}{2}+\frac{\sqrt{3 } q_y}{2}\right ) } \cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{24}({\mathbf{q}})&= & d_{42}^*({\mathbf{q } } ) = -e^{-i \left(\frac{q_x}{2}+\frac{\sqrt{3 } q_y}{2}\right ) } \cos ^2\left(\frac{\pi } { 6}-\alpha \right ) -e^{-i q_x}\cos ^2\left(\alpha + \frac{\pi}{6}\right ) \nonumber\\ d_{2,5}({\mathbf{q } } ) & = & d_{5,2}^*({\mathbf{q } } ) = ie^{-\frac{1}{2 } i \sqrt{3 } q_y } \sin ( 2 \alpha )   \sin \left(\frac{q_x}{2}\right ) \nonumber\\ d_{2,6}({\mathbf{q } } ) & = & d_{6,2}^*({\mathbf{q } } ) = -2e^{-\frac{1}{2 } i \sqrt{3 } q_y } \sin ^2(\alpha )   \cos \left(\frac{q_x}{2}\right ) \nonumber\\ d_{33}({\mathbf{q } } ) & = & 2 \sin ^2\left(\frac{\pi } { 6}-\alpha \right)+ 2 \sin ^2\left(\alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{34}({\mathbf{q } } ) & = & d_{43}^ * ( { \mathbf{q } } ) = 0 \nonumber\\ d_{3,5}({\mathbf{q } } ) & = & d_{5,3}({\mathbf{q } } ) = -e^{-i \left(\frac{\sqrt{3 } q_y}{2}-\frac{q_x}{2}\right ) } \sin ^2\left(\frac{\pi } { 6}-\alpha \right ) -e^{i q_x}\sin ^2\left(\alpha + \frac{\pi}{6}\right ) \nonumber\\ d_{3,6}({\mathbf{q } } ) & = & d_{6,3}^*({\mathbf{q } } ) = -\frac{1}{2 } e^{i q_x}\cos \left(\frac{\pi } { 6}-2 \alpha \right ) -\frac{1}{2}e^{i\left(\frac{q_x}{2}-\frac{\sqrt{3 } q_y}{2}\right ) } \cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{4,4}({\mathbf{q } } ) & = & 2 \cos ^2\left(\frac{\pi } { 6}-\alpha \right)+ 2 \cos ^2\left(\alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{4,5}({\mathbf{q } } ) & = & d_{5,4}^*({\mathbf{q } } ) = -\frac{1}{2 } e^{iq_x}\cos \left(\frac{\pi } { 6}-2 \alpha \right ) -\frac{1}{2}e^{i\left(\frac{q_x}{2}-\frac{\sqrt{3 } q_y}{2}\right ) } \cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{4,6}({\mathbf{q } } ) & = & d_{6,4}^*({\mathbf{q } } ) = -e^{i \left(\frac{q_x}{2}-\frac{\sqrt{3 } q_y}{2}\right ) } \cos ^2\left(\frac{\pi } { 6}-\alpha \right ) -e^{i q_x}\cos ^2\left(\alpha + \frac{\pi}{6}\right ) \nonumber\\ d_{5,5}({\mathbf{q } } ) & = & \sin ^2\left(\frac{\pi } { 6}-\alpha \right)+ \sin ^2\left(\alpha + \frac{\pi } { 6}\right)+2 \cos ^2(\alpha ) \nonumber\\ d_{5,6}({\mathbf{q } } ) & = & d_{6,5}^*({\mathbf{q } } ) = \frac{1}{2}\cos \left(\frac{\pi } { 6}-2 \alpha \right ) + \frac{1}{2}\cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{6,6}({\mathbf{q } } ) & = & 2 \sin ^2(\alpha ) + \cos ^2\left(\frac{\pi } { 6}-\alpha \right ) + \cos ^2\left(\alpha + \frac{\pi } { 6}\right)\end{aligned}\ ] ] its eigenvalue spectrum can easily be calculated with the aid of mathematica . the results are shown in fig .  3 of the main text . surface modes decay exponentially into the bulk . they are characterized by a wavevector @xmath321 parallel to the plane of the surface , their frequency @xmath322 , their decay length @xmath323 perpendicular to the surface . in lattices with nearest - neighbor forces only , @xmath322 and @xmath323 can be determined by setting @xmath324 , where @xmath277 is the wavevector that appears in the dynamical matrix and @xmath325 is the component of @xmath277 perpendicular to the surface ( @xmath326 is the unit vector perpendicular to the surface ) , setting @xmath327 , requiring @xmath328 = 0,\ ] ] where @xmath329 is the unit matrix , and that the equation of motion of the surface layer be satisfied . in the case of surface modes of with zero frequency , the relation @xmath330=0 $ ] determines @xmath331 as a function of @xmath321 . the surface modes of the twisted kagome lattice have zero frequency , so we need only solve this equation . we consider only the case in which the surface is perpendicular to the @xmath197-direction and set @xmath332 and @xmath333 . we find @xmath334 = \frac{9}{64 } a(q)\left[r+r^{-1 } - 2g(q)\right ] \left[r+r^{-1 } - 2g^*(q)\right ] , \label{eq : d(q , r)}\ ] ] where @xmath335 \\ g(q ) & = & \frac{\cos(q/2)[3 - 2 \cos p + \cos p \cos ( 4 \alpha ) ] + 2 i \sqrt{3 } \sin^3 ( 3 q/2 ) \sin ( 4 \alpha)}{2 - \cos q + ( 1 - 2 \cos q ) \cos ( 4 \alpha)}\end{aligned}\ ] ] solving for @xmath111 in eq . ( [ eq : d(q , r ) ] ) and to a lattice spacing of @xmath336 , we obtain @xmath337 .\ ] ] this is the function that is plotted in fig . 5(b ) in the text . note that @xmath338 has both a real and an imaginary part , but that in the limit @xmath339 , @xmath340 .
Abstract:

Answer: model lattices consisting of balls connected by central - force springs provide much of our understanding of mechanical response and phonon structure of real materials . their stability depends critically on their coordination number @xmath0 . @xmath1-dimensional lattices with @xmath2 are at the threshold of mechanical stability and are _ isostatic_. lattices with @xmath3 exhibit zero - frequency  floppy " modes that provide avenues for lattice collapse . the physics of systems as diverse as architectural structures , network glasses , randomly packed spheres , and biopolymer networks is strongly influenced by a nearby isostatic lattice . we explore elasticity and phonons of a special class of two - dimensional isostatic lattices constructed by distorting the kagome lattice . we show that the phonon structure of these lattices , characterized by vanishing bulk moduli and thus negative poisson ratios and auxetic elasticity , depends sensitively on boundary conditions and on the nature of the kagome distortions . we construct lattices that under free boundary conditions exhibit surface floppy modes only or a combination of both surface and bulk floppy modes ; and we show that bulk floppy modes present under free boundary conditions are also present under periodic boundary conditions but that surface modes are not . in the the long - wavelength limit , the elastic theory of all these lattices is a conformally invariant field theory with holographic properties , and the surface waves are rayleigh waves . we discuss our results in relation to recent work on jammed systems . our results highlight the importance of network architecture in determining floppy - mode structure .