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: while new algorithms  @xcite are pushing the frontier of quantum computation , their practical implementation requires novel ideas for the design of qubits which can cope with the contradictory constraints of scalability and long decoherence time . solid state electronics offers scalability but involves macroscopic elementary blocks which usually interact strongly with the environment and , thus , suffer from a short decoherence time . in contrast , atomic scale designs ( trapped atoms , photons in cavities , nuclear spins ) can be easily decoupled from the environment but are difficult to combine into large and complex devices . the ideal system for a qubit implementation is a macroscopic device with two quantum states decoupled from the environment which can be manufactured and combined with others by conventional techniques . this ideal is closely approximated by the superconducting phase qubit ( scpq ) that consists of a frustrated josephson junction circuit @xcite with two states distinguished _ mostly _ by their superconducting phases . the simplest example of such a circuit is a _ low inductance _ squid loop consisting of three ( or more ) josephson junctions frustrated by a magnetic field that creates a flux @xmath2 through the loop@xcite , see fig .  1 . the two states of the qubit are defined by the phase changing clockwise or counterclockwise around the loop . the advantage of an  all phase " design is that the superconducting phase is not directly coupled to the classical environment ; furthermore , superconductors themselves have a gap such that the environment has a low density of low energy modes potentially leading to decoherence and the charge fluctuations are screened such that the electric coupling to the environment is negligible . the magnetic coupling to the environment can not be completely avoided as the phase changing around the frustrated loop is coupled to an induced current , however , the latter can be made small if the josephson junctions are weak and the loop inductance is small . in particular , we shall show that the interaction of these currents with the environment leads to a small decoherence even if the environment includes potentially dangerous normal wires . finally , we note that the coupling with the frustrating magnetic field present in the simplest design can be eliminated in designs involving @xmath3-wave superconductors@xcite or @xmath4junctions@xcite . the first problem to be addressed in the development of a scpq is the implementation of a convenient probe which tests whether the device produces coherent rabi oscillations . it is obviously difficult to test a system which is decoupled from the environment and it is even more difficult to check that it remains in the phase coherent state without disturbing it . here we suggest an experiment to probe the state of such qubits and to determine their decoherence time . the idea is to measure the andreev conductance between the ground and a normal wire in the geometry of a fork with the two prongs connecting to two separate islands of the qubit ( fig .  1 ) . if the qubit phase is fixed ( classical regime ) , the low-@xmath5 conductance exhibits a contribution from processes in which electrons from the normal metal are reflected as holes from the qubit boundary , diffuse to the bulk superconductor boundary and are reflected from it as electrons . such a phase coherent electron diffusion in the normal metal prongs of the `` fork '' leads to periodic conductance oscillations ( with period @xmath6 ) as a function of the magnetic flux penetrating through the @xmath7 region of the fork.@xcite this is due to the fact that the magnetic field controls the superconducting phase difference between the two ns contacts of the fork [ ( 1 ) and ( 2 ) in fig .  1 ] and thus influences the electron interference picture . a similar mechanism controls the experiment discussed in the quantum regime when the phase of the qubit fluctuates . qualitatively , at small bias voltage the quantum phase acquired by the reflected electrons fluctuates and the resulting contribution to the conductivity averages out . at larger voltages the phase fluctuations become slow compared with the electron tunneling time and the contribution to the conductivity is restored . below we describe the specific device parameters and derive quantitative formulas describing this effect . a similar proposal  to probe the josephson splitting of levels in a superconducting single cooper - pair transistor by measuring the andreev conductance  was put forward in ref .  . consider the squid loop with four junctions shown in fig .  1 . two degenerate states naturally appear in such a loop if the flux @xmath8 of the external field through the qubit loop is exactly @xmath9 . indeed , in a gauge with @xmath10 , @xmath11 the classical minima of the josephson energy correspond to the phase drops @xmath12 or @xmath13 around the loop ; in these states the phase of the island @xmath14 ( see fig . 1 ) is @xmath15 ( relative to the phase at point @xmath12 ) ; below we refer to these states as @xmath16 and @xmath17 . in order to reduce the parasitic coupling to the environment , the inductance @xmath18 of the loop should be sufficiently small so that @xmath19 , @xmath20 is the critical current . in the absence of any charging energy in the junctions the system prepared in one of these states will stay put forever ; quantum effects are due to the small but finite charging energy @xmath21 determined by the capacitances @xmath22 of the junctions , which should be smaller than the josephson energy @xmath23 . the tunnelling frequency between the two classically degenerate ground states is estimated as @xmath24 , where @xmath25 is the josephson plasma frequency of the contacts ; the coefficient @xmath26 was found by a numerical evaluation of the saddle - point trajectory in the space of the three phase variables @xmath27 , @xmath28 , and @xmath29 . we assume values of @xmath20 in the range @xmath30 a , and capacitances of order of a few ff , resulting in @xmath31 ghz and a tunneling frequency @xmath32 in the range of few ghz . our choice of parameters is limited by the following constraints : ( a ) the system should be sufficiently small to fulfill the condition @xmath33 , but not too small , since ( b ) we need @xmath34 ; finally , ( c ) @xmath35 should not be too small , in order to keep the decoherence time @xmath36 much longer than @xmath37 . once tunnelling is taken into account , the true eigenstates become @xmath38/\sqrt{2}$ ] and @xmath39/\sqrt{2}$ ] and are separated by the energy gap @xmath40 ( here and below we shall use the  spin " notation to describe the _ classically _ degenerate states ) . the deviation of the magnetic flux @xmath41 through the loop from the point @xmath42 removes the degeneracy of the states @xmath16 and @xmath17 . thus the intrinsic hamiltonian of our qubit written in the basis @xmath43 is @xmath44 , where @xmath45 $ ] and @xmath46 . all operations on the qubit can be performed if one is able to vary the effective fields @xmath47 and @xmath48 . the variation of @xmath47 can be achieved by changing the flux @xmath8 in the loop , while the variation of @xmath48 can be implemented by two means : the variation of the gate potential applied to the island @xmath14 changes this field smoothly ( c.f . , ref .  ) , while switching a capacitor @xmath49 in parallel blocks the tunneling abruptly . the coupling of a two - level quantum system to the environment is usually described by the caldeira - leggett  spin - boson " model where the  spin " ( i.e. , two - level system , tls ) is linearly coupled to an ensemble of oscillators.@xcite actually , our tls originates from the dynamics of a continuous phase variable @xmath50 with a potential energy @xmath51 strongly favoring values @xmath50 around @xmath52 . then the term that describes the coupling to the environment can be written , upon integration over environmental modes , in the form of a non - local imaginary - time action @xmath53 ^ 2 $ ] . while in a superconducting system we can hope the low energy spectrum of the environment to be gapped providing favorable conditions for a long coherence time , here , we analyze the ( worst case ) situation where the coupling of low lying modes survives to produce a finite ohmic dissipation with @xmath54 , where @xmath55 is the generalized  friction coefficient " in the equation of motion @xmath56 . the strength of the environmental coupling is determined by the dimensionless parameter @xmath57 ; consider the symmetric system with @xmath58 : for @xmath59 the tunnelling between the states @xmath16 , and @xmath17 is suppressed and the system possesses two classically degenerate ground states . for intermediate coupling @xmath60 tunnelling is finite but _ incoherent _ : at time scales @xmath61 ( with @xmath62 the effective tunneling frequency depending on @xmath63 ) the tls is described by the classical master equation for the occupation probabilities @xmath64 of the two states with a transition rate @xmath65 vanishing at @xmath66 . at @xmath67 this equation takes the form @xmath68 and the  magnetization " @xmath69 monotonically decays to zero . here , we are interested in the case of very weak coupling @xmath70 , where the tls dynamics is weakly perturbed by the environment  @xcite , i.e. , the decoherence time @xmath71 shall be much longer than the period of coherent oscillations @xmath72 . in this case , the magnetization exhibits damped oscillations @xmath73 with a damping rate @xmath74 and a frequency @xmath62 determined by the parameter @xmath75 via  @xcite @xmath76 ,   \qquad \gamma = \omega / q = \omega \tan[\pi\alpha/2(1-\alpha ) ] , \label{lesage}\ ] ] with the renormalized tunnelling frequency @xmath77 ; here , @xmath78 is an upper cutoff for the frequencies of the environmental modes . a measurement of @xmath79 , particularly its quality factor @xmath80 , then allows us to test the degree of quantum coherence in the device .    the direct determination of @xmath79 in a time - domain experiment like the one reported in ref . @xcite ( where a superconducting  single cooper - pair " device was studied ) is a demanding task . as an alternative one may attempt to find the correlation function @xmath81 through a measurement@xcite of the _ rf _ absorption ; indeed , the latter provides the imaginary part @xmath82 of the response function @xmath83 , where @xmath84 \rangle/2i$ ] is the antisymmetric part of the spin - spin correlation function . standard relations then allow us to extract the correlator @xmath85 from the measured absorption @xmath82 : the symmetrized correlation function @xmath86 relates to the correlator via @xmath87 $ ] ( note that @xmath88 at @xmath89 ) . below we discuss yet another experiment giving us access to @xmath85 : we will see that the nonlinear andreev conductance @xmath90 through the qubit s island @xmath14 can be expressed via the fourier - transform @xmath91 of its correlator . scanning magnetic field and voltage in a transport measurement then appears to give a much easier and direct access to the operation of the qubit than the more usual real - time and _ rf_-techniques . being not aware of exact results available for @xmath92 or @xmath93 , we make use of an approximate analysis  @xcite to arrive at an idea of the relevant dependencies in these quantities : for @xmath94 the absorption characteristic is close to a combination of lorentzians : @xmath95 -\gamma/[\gamma^2 + ( \omega + \omega)^2]$ ] , implying a symmetrized correlator of the form @xmath96 $ ] . in the time domain , @xmath97 $ ] , where the last term is due to the non - analyticity of @xmath98 at @xmath99 . we first concentrate on the subgap conductance through a nis contact ( # 1 in fig . 1 ) with a low normal - state conductance @xmath100 and discuss the effect of fluctuations in the phase @xmath101 on the island @xmath14 . consider an ideal scpq operating at the degeneracy point ( @xmath58 ) with the ground and excited states @xmath102 and @xmath103 separated in energy by @xmath104 . the operator @xmath105 creating cooper pairs changes the relative sign of the states @xmath106 and @xmath107 , thus latexmath:[$c^\dagger|\,\!0\rangle =    ground state . therefore @xmath102 is a state with a well - defined number of cooper - pairs _ modulo _ 2 , implying that the usual andreev process transmitting a single cooper pair through the interface is forbidden at low voltage @xmath109 , while at voltages above the gap @xmath110 the andreev conductance will be large . higher - order processes with simultaneous tunnelling of two cooper pairs are allowed , but their rate is of the order of @xmath111 and thus is negligible at @xmath112 . overall , we then expect an andreev conductance @xmath90 similar to the one of a usual nis junction but with the superconducting gap @xmath113 replaced by the two - level separation @xmath104 . on the other hand , for a fixed phase on the island we can expect to observe the usual andreev current across a nis junction . the above consideration defines the simplest example of a situation where the quantum fluctuations of the superconducting phase suppress the andreev conductance of the sin contact . general expressions describing these effects have been derived in ref . ( c.f . , appendix b ) . here , we discuss the result for the nonlinear subgap current @xmath114 for the case where the normal metal is relatively clean ( with a dimensionless conductance @xmath115 ) such that all corrections of relative order @xmath116 can be neglected , @xmath117-\tanh[(e+e'/2-ev)/2t]\right ] . \nonumber\end{aligned}\ ] ] here , @xmath118 is the real part of the cooperon amplitude in the normal metal ( see below ) and @xmath119 is the fourier - transformed autocorrelation function @xmath120 of the superconducting order parameter @xmath121 $ ] on the island @xmath14 . on long time scales @xmath122 the oscillatory modes near the minima of the josephson potential are irrelevant ; the values of @xmath50 in these minima differ by @xmath13 , thus @xmath123 is equivalent to the ( sought for ) tls correlation function @xmath85 discussed in the previous section . the real part of the cooperon amplitude is given by @xmath124 where the sum goes over the cooperon eigenmodes @xmath125 , @xmath1 denotes the volume , @xmath126 the density of states ( for a single projection of spin ) , and @xmath127 the diffusion constant of the normal metal . assuming a wire geometry with @xmath128 , @xmath129 , and @xmath130 denoting its length , width and thickness , eq . ( [ cooperb ] ) reduces to the expression @xmath131^{1/4}$ ] , where the resistance @xmath132 per unit length of the wire is measured in units of @xmath133 . the term @xmath134 describes the electron dephasing present at all @xmath135 due to electron - electron interactions and external noise  @xcite . for a short wire , @xmath136 ; in the end we arrive at the simple result @xmath137 . an expression for @xmath138 analogous to eqs . ( [ main ] ) has been derived earlier  @xcite for a similar physical problem . let us analyze the result ( [ main ] ) in the limit @xmath139 which applies to the regime @xmath140 . at the same time we keep the electron dephasing rate @xmath141 nonzero , since experimentally  @xcite an apparent saturation at @xmath142 has been found . differentiating eq . ( [ main ] ) by @xmath1 at @xmath89 we find @xmath143 ( values @xmath144 are excluded from the integral since @xmath145  at  @xmath89 ) .    consider first of all the degenerate case with @xmath67 in the absence of any dephasing within the qubit , i.e. , @xmath146 and @xmath147 . the evaluation of the integral in ( [ zerot ] ) is trivial and we obtain the result @xmath148 in agreement with the above qualitative arguments . in the following we restrict the discussion to long wires @xmath149 ; then the conductivity above threshold shows a square - root singularity which is smeared by a nonzero dephasing rate @xmath150 . the nonzero conductance below threshold is due either to a finite dissipation rate @xmath151 or a finite temperature . at zero temperature , a finite @xmath74 produces the sub - threshold ( @xmath152 ) conductance @xmath153^{1/2 } ,               } & $ \delta \equiv \hbar\omega - 2ev \ll \hbar\omega/2 $ , } \ ] ] where we have assumed @xmath154 and have made use of the form of @xmath155 discussed at the end of sec .  2 . the numerical coefficient @xmath156 depends on the shape of the correlator @xmath155 and is of order unity ; for the lorentzian shape discussed above we find @xmath157 . a finite phase breaking time @xmath158 in the normal wire changes this result at low voltages @xmath159 : @xmath160 $ ] . finally , a finite temperature also produces a finite sub - threshold conductance which , however , can be well separated ( at @xmath161 ) from the above behavior due to its exponential @xmath162 $ ] rather than algebraic dependence @xmath163 and @xmath164 at low and high voltages . note , that these results for the suppression in the conductance @xmath90 due to phase fluctuations refer to terms of order @xmath165 and hence our analysis is valid only for low @xmath166 . the above discussion clarifies the appearance of a threshold behavior in the degenerate case @xmath58 , where the differential conductance @xmath90 sharply drops below the gap @xmath110 and vanishes at @xmath167 . for @xmath168 the correct ground state is close to a state with a well - defined phase @xmath50 and the conductance attains its semiclassical value @xmath169 . as a consequence , a sharp minimum ( at @xmath42 corresponding to @xmath67 ) is expected in @xmath114 when changing the magnetic flux @xmath8 through the qubit . this qualitative conclusion survives the presence of a finite electron decoherence rate @xmath150 in the normal - metal , however , the precise form of @xmath90 is modified . a short electron decoherence time @xmath170 is frequently observed in thin wires at lowest temperatures as relevant for the experiments on qubits ( c.f . , ref .  ) . these observations are poorly understood theoretically and thus it is more practical to measure @xmath158 independently in the same experimental setup . such a measurement can be carried out via an analysis of the magnetic field dependent andreev conductance in the  fork " geometry as we are now going to discuss . consider the interference experiment for the device shown in fig.1 . the first stage of the proposed experiment is to measure the amplitude @xmath171 in the oscillations of the andreev current @xmath172 as a function of magnetic field @xmath173 with the aim to extract the electron decoherence time @xmath158 in the normal metal ( the superscripts @xmath174 refer to the contributions from the contacts @xmath175 , see fig .  1 ) . the total phase determining the interference current @xmath171 sums up to @xmath176 . away from the degeneracy field @xmath177 , the ( semiclassical ) phase @xmath50 in our device is determined by the minimization of the josephson energy and we obtain @xmath178 , where @xmath179 is the fractional part of @xmath180 . choosing a geometry with @xmath181 , we find for @xmath182 $ ] an interference current @xmath183 , whereas for @xmath184 $ ] the sign changes : @xmath185 . the amplitude  @xcite of the interference component @xmath186 depends on the precise geometry of the  fork " and on the electron dephasing length @xmath187 which should be longer than the distance @xmath188 between the contacts # 1,2 ( otherwise @xmath186 will be exponentially small ) . detailed calculations and measurements of @xmath186 have been presented by the saclay group  @xcite ; they found a discrepancy by a factor 1.8 between the best fit ( with @xmath158 as fitting parameter ) to their data and the calculated expression for @xmath186 ( c.f . , ref . , p.  190 ; this discrepancy may be due to the use of a simplified description of the electron decoherence using a factor @xmath189 in the real time expression for the cooperon  @xcite ) . second , we test for the coherent time evolution of the qubit : in the geometry with @xmath181 , the semiclassical current amplitude @xmath190 is periodic with a period @xmath191 , showing upward cusps at @xmath192 , see fig .  2 , dotted line . these cusps are due to the penetration of single flux quanta @xmath6 into the qubit loop at @xmath193 . at these fields , the interference contribution @xmath194 vanishes and the semiclassical andreev current is given by the sum @xmath195 ( in the generic case @xmath196 the component @xmath171 does not vanish at @xmath193 and the flux penetration leads to a jump in the measured @xmath197 ) . in addition , at @xmath198 the quantum fluctuations in the phase @xmath50 become large and suppress the andreev current @xmath199 to the  active " island g of the scpq , see eq . ( [ zerot ] ) . hence , close to the points @xmath193 a narrow dip is superimposed upon the above - mentioned cusp in the @xmath200 dependence , see fig .  2 , solid line . the relative width @xmath201 of this dip is determined by the relation @xmath202 and provides the coherence gap of the qubit . in addition , we can measure the decoherence rate @xmath74 via an analysis of the conductance @xmath203 at the degeneracy points @xmath193 .  in the limit @xmath204 and @xmath205 the entire conductance is determined  by  the current through the contact # 2 , @xmath206 . the variation of the full differential conductance @xmath207 at low voltages @xmath208 then can be used to extract the value of @xmath74 . note that @xmath118 for the fork differs from the form used in the expression for the subgap conductance , c.f . , ref .. the detailed analysis of @xmath209 for the different relations between @xmath62 , @xmath74 , and @xmath158 will appear in a future publication .    finally , any method measuring the scpq coherence time is useful only if the presence of the measurement circuit does not , by itself , lead to a decoherence rate comparable to the  intrinsic " one . in our setup there are two sources of additional phase @xmath50 decoherence due to the presence of the normal - metal wire nearby : ( i ) the direct coupling of the qubit to the environment via the non - zero andreev conductance , and ( ii ) the inductive coupling between _ ac _ supercurrents in the qubit loop and normal currents in the andreev fork . we consider them one by one : ( i ) comparing the dissipative action for the andreev conductance ( c.f . , ref .   for example ) with the general relation for @xmath210 and the definition of the parameter @xmath57 one finds the contribution @xmath211 , where @xmath212 is the semiclassical andreev conductance via the contact # 1 . thus , a classical subgap resistance of the sn contacts in the m@xmath62 range will produce a value @xmath213 . ( ii ) we calculate the friction coefficient @xmath214 via the energy dissipation rate @xmath215 . the main contribution to @xmath216 is due to the _ ac _ magnetic field applied to the conducting wire ( with thickness @xmath217 nm , width @xmath218 nm and a few microns length @xmath219 ) , which gives @xmath220 , where @xmath221 is the imaginary part of the magnetic susceptibility of a thin - film wire in a magnetic field @xmath222 in the low - frequency limit ( for frequencies @xmath223 and a conductivity @xmath224 cm of the wire the effective skin depth is much longer than the wire s width @xmath218 nm ) . the magnetic field outside the scpq loop ( of size @xmath225 m ) is less than @xmath226 , thus @xmath227 . combining all factors , one finds @xmath228 . for a scpq with @xmath229 and @xmath230 this results in a negligible value @xmath231 .    in conclusion , we have derived expressions for the subgap @xmath232 characteristics , eqs . ( [ main ] ) and ( [ zerot ] ) , which can be used for the determination of the intrinsic correlation function @xmath233\,\exp[-i\phi(0)]\rangle$ ] of a superconducting phase qubit ( scpq ) , once the electron dephasing time @xmath158 that enters the cooperon , eq . ( [ cooperb ] ) , is known . as the latter is difficult to control theoretically at the ultralow temperatures of the experiment we propose an interference experiment probing the scpq with an  andreev fork " , allowing to determine @xmath158 within the same setup and then use its value to extract @xmath36 . in the absence of decoherence in the scpq , @xmath234 , the current @xmath114 is zero ( at @xmath89 ) below the threshold @xmath235 , irrespective of the value of @xmath158 . thus the scpq decoherence time @xmath236 can be measured even if it is longer than @xmath158 . in practice , the most important limitation of this technique seems to be due to the non - vanishing temperature of experiment which produces a nonzero sub - threshold conductance even at @xmath234 . however , in the range @xmath237 and at voltages such that @xmath238 the behavior of @xmath90 is exponentially close to the zero-@xmath5 one ; therefore even values of @xmath239 can be extracted from the measured @xmath90 curves . we thank g.  lesovik , j.  mooij , m.  skvortsov , and c.  van der wal for helpful discussions and the ` zentrum fr theoretische studien ' at eth - zrich for hospitality and financial support . the research of m.v.f .  was supported by the rfbr grant 98 - 02 - 19252 and by the russian ministry of science via the programm  physics of quantum computing " .
Abstract:

Answer: we propose a new method to measure the coherence time of superconducting phase qubits based on the analysis of the magnetic - field dependent _ dc _ nonlinear andreev current across a high - resistance tunnel contact between the qubit and a dirty metal wire and derive a quantitative relation between the subgap @xmath0-@xmath1 characteristic and the internal correlation function of the qubit . pacs numbers : 03.67.lx , 85.25.dq


Instruction: Write an abstract for the following acadamic paper.

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

Paper: the coupling between dynamical systems can give rise to a number of collective phenomena such as synchronization , phase locking , phase shifting @xcite , amplitude death @xcite , phase - flip @xcite , hysteresis @xcite , riddling @xcite and so on @xcite . since communication between the individual systems is mediated by signals that can have a finite transmission time , many studies account for this by the introduction of time  delay in the coupling @xcite . a number of recent studies have examined the manner in which delay coupling can affect the collective dynamics , particularly since time  delay makes the systems effectively infinite dimensional @xcite .    when conservative systems are coupled via time  delayed interactions , then there are additional considerations . to start with , the system becomes explicitly non  conservative and thus the nature of the dynamics changes drastically : in the uncoupled system , the phase flow preserves volumes @xcite , but in the coupled system there can be attractors . this issue is of added interest when studying hamiltonian systems where there can be a hierarchy of conserved quantities @xcite . studies of coupled hamiltonian systems have largely examined the case of instantaneous coupling @xcite which does not affect the hamiltonian structure .    in the present work we study time  delay coupled hamiltonian systems and examine the effect of interaction on the nature of the dynamics . we consider the following examples , that of diffusively coupled harmonic oscillators that models delay  coupled pendulums , for instance , and coupled hnon  heiles oscillators . in the absence of coupling , in the former case the motion is periodic , while in the latter case the dynamics can be ( quasi)periodic or chaotic . in both instances we find that the effect of introducing dissipation is to cause the oscillatory dynamics to be damped to a fixed point , namely we find that there is the so  called amplitude death ( ad ) @xcite as has been seen in delay  coupled dissipative dynamical systems @xcite . although the major effect of the coupling is to make the overall system dissipative , there are differences from the case when non - hamiltonian systems are coupled . when the dynamics is decaying to a point attractor , there is an abrupt transition in the relative phases of the oscillatory transient motion . this is the phase  flip transition that has been seen in a number of other systems @xcite . here , however , there are special values of the time  delay when the coupling term effectively vanishes : the underlying hamiltonian structure then becomes apparent . our main results are presented in sections ii and iii where we consider the cases of coupled harmonic oscillators and coupled hnon heiles systems respectively . we show how ad is reached and the nature of the phase - flip transition in both cases . since the uncoupled systems are hamiltonian , it is possible to define an energy , and while this quantity has been studied in coupled feedback oscillators @xcite as a tool to determine onset of ad , its variation in the ad regime itself has not been examined . in this work we do energy analysis in the ad region and find that the energy dissipation is non monotonic as a function of the coupling , decaying faster prior to the phase  flip transition and slower subsequently . the paper concludes in section iv with a summary and discussion of the results . the simplest system we consider is that of diffusively coupled harmonic oscillators . we consider the following equations of motion @xmath0=0 \label{eq : shm}\end{aligned}\ ] ] where @xmath1=1 , 2 and @xmath2 and @xmath3 and @xmath4 represent the position and the velocity of the @xmath5th oscillator , and @xmath6 is the intrinsic frequency . we take the oscillators to be identical @xcite , @xmath7 . the parameters @xmath8 and @xmath9 represent coupling strength and time delay respectively . in the absence of delay , @xmath10 causes the systems to synchronize completely . due to the simple dynamics of the system no other significant behavior is observed . when , delay is finite the coupling quenches oscillations leading to ad . stability analysis of eq . ( [ eq : shm ] ) around the fixed point , namely the origin , gives the characteristic equations @xmath11taking the roots of the jacobian to be @xmath12 , the condition for marginal stability is @xmath13 , and substituting this condition in eq . ( [ eq : ch ] ) we get @xmath14 where @xmath15 is the critical value at which @xmath13 and @xmath16 is the time - period of the uncoupled oscillators . shown in fig . [ fig : ho1](a ) are the first two lyapunov exponents @xcite of the system as a function of @xmath9 for fixed coupling strength @xmath8 = 1 . the largest le is zero only at @xmath17 ( marked in fig . 1(a ) as b and d ) and remains negative for all other values of @xmath9 , implying that the system is driven to ad except when the delay is an integral multiple of half the time period . further , at these critical delay values the system oscillates at the frequency of the uncoupled system namely @xmath18 , and the parameter space is divided into multiple ad regions by the critical delay values which are independent of the coupling strength @xmath8 . in contrast , in non  hamiltonian systems ad islands are separated by finite range of delay values @xcite where the coupling function need not vanish . hence , in those systems the reappearance of oscillations after ad depends both on the coupling strength and on the delay , whereas in coupled hamiltonian systems we find that this happens only due to delay .    in the ad regime(s ) points of discontinuity in the slope of the largest le ( marked in fig . 1(a ) as a and c ) indicate the change in the relative phases of oscillation . this is the so  called phase - flip transition @xcite , and the difference in the phases of the coupled subsystems changes by @xmath19 ; see ( fig.[fig : ho1]b ) . as in other cases where this phenomenon has been observed , there is simultaneously a discontinuous change in the oscillation frequency @xcite , as shown in fig.[fig : ho1](c ) .    ) as a function of the time  delay @xmath9 for fixed @xmath10 = 1 ; ( b ) phase difference between the oscillators around point a indicated in ( a ) . inset figures show trajectories @xmath20 and @xmath21 as a function of time for the in  phase and the out  of  phase dynamics at @xmath9 = 1 and 2 respectively , namely on either side of the phase  flip transition , and ( c ) the frequency of oscillation of the subsystems . ]    the uncoupled conservative hamiltonian systems are made dissipative through the time  delay coupling , and this is also reflected in the fact that the sum of all the les remains negative for all @xmath9 . defining the energy of the individual oscillators as @xmath22 where @xmath23 is the instantaneous frequency of oscillation , the approach to the fixed point can be seen to be at an exponential rate in the ad regime @xcite as can be seen in fig . [ fig : ho2](a ) . the exponential decay is however modulated , the oscillatory behavior being due to coupling ( see the inset in the figure ) . in order to capture the dynamics , we define a decay constant as @xmath24 where @xmath25 represents the @xmath26-th maxima in the energy time series of the @xmath5th oscillator . the averages @xmath27 and @xmath28 are performed on @xmath26 and 100 initial conditions respectively . since the oscillatory behavior is modulated by exponential decay , @xmath29 is also an exponentially decaying function of @xmath26 , at rate @xmath30 . this rate can be measured experimentally and its variation with @xmath9 is shown in fig . [ fig : ho2](b ) ; the variation mirrors that of the frequency change at phase  flip , suggesting that energy decays more rapidly before the transition than after it .    ) , as a function of time for @xmath31 . the marked box is expanded in the inset . ( b ) the decay rate , @xmath32 as a function of @xmath9 ( see eq . ( [ eq : decay ] ) ) . errorbars ( marked in red ) are calculated for 100 initial conditions . ]    at the critical delay values @xmath33 ( see points b and d in fig . [ fig : ho1](a ) ) the largest lyapunov exponent is zero and thus the motion is periodic . [ fig : ho3](a ) shows orbits for five different initial conditions at the critical point ; these resemble invariant curves as in conservative systems . however , in the vicinity of the critical delay values the largest lyapunov exponent is near zero and the motion appears periodic after an initial transient . the time series of one such periodic orbit is shown in fig . [ fig : ho3](b ) . also , as can be seen in the inset , there is very slow decrease in amplitude . this occurs since delay is not strictly @xmath34 . away from delay @xmath34 , the dissipation is more pronounced . the rate of decrease of the amplitude can be quantified through the measure    @xmath35    where @xmath36 is the @xmath26-th maxima of @xmath3 . this is plotted in fig . [ fig : ho3](c ) as a function of @xmath9 in the vicinity of @xmath37 , namely the point b in fig . [ fig : ho1](a ) . at @xmath33 the rate of decrease of amplitudes approaches zero and hence the orbits are almost periodic . similar behavior is observed at point d of fig.[fig : ho1](a ) . the reason for reappearance of oscillatory motion is straightforward . when the delay is a multiple of half the natural period of oscillation , then the coupling term effectively vanishes since @xmath38 and the system effectively becomes conservative . clearly when this occurs , each initial condition gives rise to an invariant curve ( or , in this case , a nearly invariant curve ) . the better the equality above is realized , the more long lived the transients are . there is the phase  flip transition at a , c and also higher values of @xmath9 . at each transition a phase difference of @xmath19 is introduced between the oscillators resulting in anti  phase synchronization at @xmath39 and in  phase synchrony as @xmath40 and so on . hence , the coupling function becomes zero at these critical delays . the consecutive oscillatory states alternate between having the oscillations being in  phase or out  of  phase , as shown in fig . [ fig : ho3](d ) ( out  of  phase at @xmath37 , namely at b ) and fig . [ fig : ho3](e ) ( in  phase at @xmath41 at the point marked @xmath42 ) respectively . note that the phase - relation between oscillators is independent of the initial conditions .    ; ( b ) time series of position variable @xmath20 as a function of time at @xmath43 . the inset figure shows decrease in amplitude of @xmath20 ; ( c ) the averaged distance @xmath44 with errorbars , as function of delay @xmath9 near @xmath37 . the errorbars are calculated for 100 initial conditions ; ( d ) out - of - phase motion and ( e ) in - phase motion in relative phase plane @xmath45 at @xmath37 and @xmath46 respectively ( corresponding to points b and d of fig . [ fig : ho1](a ) ) . ] in order to examine the dynamics when the uncoupled systems are capable of exhibiting chaotic motion , we examine the behavior of two non - integrable hnon - heiles systems @xcite , @xmath47    as is well  known , in the uncoupled case ( @xmath10=0 ) the system has both regular and irregular behavior largely depending on the total energy as well as on the initial condition @xcite . shown in fig . [ fig : hh0 ] are the poincar maps for two different initial conditions , one leading to regular motion ( red points ) , while another leads to chaotic dynamics ( black points ) , at the same energy @xmath48 = 0.13 just below the dissociation limit @xmath48 = 1/6 .    ) , at energy @xmath49 . red ( outer ) and black ( inner ) dots represent regular and chaotic motion respectively . ]      in the absence of delay @xmath51 , eq . ( [ eq : hh ] ) reduces to the case of simple diffusive coupling . the effect of increasing the coupling strength , namely @xmath10 , is to induce simplicity to the resulting collective dynamics . shown in fig.([fig : hhic ] ) are the the fraction of initial conditions @xmath52 leading to quasiperiodic motion . we take 100 pairs of random initial conditions from the bounded region of phase space ( fig . [ fig : hh0 ] ) . in one case , when the initial condition pairs of quasiperiodic motions are taken , the collective dynamics due to interaction remains quasiperiodic ( solid line ( black circles ) : @xmath53 ) . however , if the initial motion is chaotic then the resulting dynamics becomes quasiperiodic only after certain value of coupling strength ( dashed line ( red triangles ) : @xmath54 ) . similar behavior is observed when initial conditions of mixed chaotic and quasiperiodic motions are used ( dotted line ( green stars ) : @xmath55 ) . these results indicate that for small coupling strength the hamiltonian structure still exists , but for larger values of coupling the collective dynamics becomes quasiperiodic ; in this sense coupling induces simplicity in such coupled systems .    ) , chaotic ( @xmath54 ) and mixed chaotic and quasiperiodic @xmath56 motions respectively . averages have been taken over a sample of 100 initial conditions . ] 0.2 cm [ fig : hhic ]      , and ( b ) the oscillation frequency , for coupling strength @xmath8 = 1 . the inset shows trajectories @xmath57 and @xmath58 for @xmath9 = 1.4 and 1.8 respectively . ( c ) the decay constant @xmath32 as a function of the delay @xmath9 . ]    when the two systems are coupled in presence of delay , @xmath9 ( initial conditions taken from either the regular or irregular motion ) the largest lyapunov exponent quickly becomes negative . this is shown as a function of the delay in fig . [ fig : hh1](a ) : almost as soon as the delay is switched on the oscillators are driven to ad state . the phenomenology of this higher dimensional system is very similar to that of the coupled harmonic oscillators : the largest le has the same characteristic shape  there is a change in the slope at point a , where it is clear that the phase - flip transition occurs . shown in fig . [ fig : hh1](b ) is the common frequency of oscillations as a function of delay which changes discontinuously at @xmath591.65 . since the phase of the oscillators is not clearly defined in such systems , we infer the phase relation from the time  series ( inset figures ) of the two systems in the neighborhood of the point of discontinuity . the phase difference changes from 0 to @xmath19 along with the frequency jump as in the simpler 1dimensional harmonic system .    here also energy decreases exponentially in this ad region . we define the energy of individual systems in the usual manner @xcite @xmath60 and quantify the energy dissipation as in eq . ( [ eq : en ] ) by a quantifier @xmath30 . the variation of this quantity with @xmath9 can be seen in fig . [ fig : hh1](c ) . it confirms that energy dissipation is faster before phase flip transition and slower thereafter .    when the largest lyapunov exponents approach zero ( at points marked b and d in fig . [ fig : hh1](a ) ) the motion is oscillatory , decaying very slowly to the fixed point . critical delays are at half the natural time period of oscillation , and since the natural frequencies of the oscillators is equal to ` 1 ' , the time period is 2@xmath19 . the poincar sections of representative trajectories at @xmath43 are shown in fig . [ fig : hh2](a ) which are quasiperiodic ( cf . section iia ) . in the vicinity of this point , the rate of decay can be computed numerically as discussed earlier , and the quantifier @xmath61 , eq . ( [ eq : dis ] ) defined to locate the oscillatory state also has a minimum for at @xmath62 . at successive critical points ( @xmath43 and @xmath63 ) the quasiperiodic motions alternate in the nature of the phase synchrony ; see fig . [ fig : hh2](d , e ) . unlike the case of harmonic oscillators , though , the coupling term does not quite vanish , so the emergence of oscillations is not as pronounced in this case .    0.8 cm   for three different initial conditions ; ( b ) time series of @xmath20 ; ( c ) variation of @xmath44 as a function of delay . the error bars are computed from a large number of initial conditions . ( d ) the relative phase in @xmath64 plane to show anti - phase motion at @xmath65 and ( e ) in - phase motion at @xmath66.,title="fig : " ] in the present paper we have explored the effects of incorporating time  delay in the coupling between hamiltonian systems . the coupled system effectively becomes dissipative , and in the absence of other attractors , the total system displays amplitude death . however the coupling can effectively vanish for specific values of the time  delay , at which delay the the system naturally appears to be conservative . we find that this behaviour differs from that of delay  coupled non  hamiltonian systems , where ad islands are separated by finite ranges of delay values where the coupling function need not vanish . orbits near these critical delay values reflect both dissipative and conservative behavior : different initial conditions give rise to different ( seemingly ) invariant curves which are very long lived transients . energy dissipation in the ad regime is found to be associated with the phase  flip transition , and the damping is faster prior to the transition , and slower after it .    the finite velocity with which signals are transmitted gives rise to intrinsic delays in the coupling , and as such this is germane in both dissipative as well as conservative systems . nevertheless , the effect of delay has been explored to a limited extent in conservative systems . the results presented here may have more general applicability in coupled systems with other conservation laws such as in ecological contexts @xcite . * acknowledgment : * gs gratefully acknowledges the support of the csir , india . ap and rr acknowledge the research support of the department of science and technology ( dst ) . ap also acknowledges dst for du - dst - purse grant for financial support . this work was started when ap was a sabbatical visitor at the mpi - pks in dresden , germany , and he acknowledges their kind hospitality . y. yamaguchi and h. shimizu , ; k. bar - eli , . d. g. aronson , , ; p. c. mathews and s. h. strogatz , ; g. b. ermentrout , ; v. resmi , g. ambika and r. e. amritkar , ; v. resmi , g. ambika , r. e. amritkar and g. rangarajan , . g. saxena , a. prasad , and r. ramaswamy , phys . reports * 521*,205 ( 2012 ) . a. prasad , . a. prasad , j. kurths , s. k. dana and r. ramaswamy , ; a. prasad , s. k. dana , r. karnataka , j. kurths , b. blasius and r. ramaswamy , chaos * 18 * , 23111 , ( 2008 ) ; c. masoller , m. torrent and j. o. garcia , ; y. chen , j. xiao , w. liu , l. li and y. yang , ; j. m. cruz , j. escalona , p. parmananda , r. karnatak , a. prasad , and r. ramaswamy , phys . rev . e. * 81 * , 046213 ( 2010 ) ; b. m. adhikari , a. prasad , m. dhamala , chaos * 21 * , 023116 ( 2011 ) ; r. karnatak , p. nirmal , a. prasad and r. ramaswamy , phys . e. * 82 * , 046219 ( 2010 ) ; p. nirmal , r. karnatak , a. prasad , j. kurths and r. ramaswamy , . s. kim , s. h. park and c. s. ryu , ; a. prasad , l. d. iasemidid , s. sabesan and k. tsakalis , pramana * 64 * , 513 ( 2005 ) . h. g. schuster and p. wagner , . k. konishi , m. ishii and h. kokame , ; k. konishi , h. kokame and k. hirato , ; k. konishi , ; k. konishi , k. senda and h. kokame , ; k. konishi , h. kokame and n. hara , ; l. b. le , k. konishi and n. hara , nonlinear dyn . * 67 * , 1407,(2012 ) . d. v. ramana reddy , , ; s. h. strogatz , nature ( london ) * 394 * , 316 ( 1998 ) ; f. m. atay , ; f. m. atay , j. jost , ; f. m. atay , ; g. saxena , a. prasad , and r.ramaswamy , phys . rev . e. * 82 * , 017201 ( 2010 ) . w. zou , c. yao and m. zhan , ; w. zou , j. lu , y. tang , c. zhang and j. kurths , ; w. zou , y. tang , l. li and j. kurths , ; w. zou , d. v. senthilkumar , y. tang and j. kurths , . m. tabor , _ chaos and integrability in nonlinear dynamics : an introduction _ ( john wiley & sons , new york , 1989 ) . v. i. arnold , _ mathematical methods of classical mechanics _ ( springer , new york , 1978 ) . d. h. zanette and a. s. mikhailov , ; a. hampton and d. h. zanette , . the flow is integrated using the runge - kutta @xmath67 order scheme with integration step @xmath68 where @xmath69 is fixed . as we increase the value @xmath70 ( checked upto @xmath71 ) the decrease in amplitude in fig . 3(b ) becomes slower . lyapunov exponents are calculated using the method as given in ref . phase and frequency are numerically calculated as per ref . @xcite .    understanding the transient dynamics is essential in describing ad since this state is asymptotically featureless . transients can be significant in applications that are restricted to finite times , as for example in ecology @xcite .
Abstract:

Answer: hamiltonian systems , when coupled _ via _ time  delayed interactions , do not remain conservative . in the uncoupled system , the motion can typically be periodic , quasiperiodic or chaotic . this changes drastically when delay coupling is introduced since now attractors can be created in the phase space . in particular for sufficiently strong coupling there can be amplitude death ( ad ) , namely the stabilization of point attractors and the cessation of oscillatory motion . the approach to the state of ad or oscillation death is also accompanied by a phase  flip in the transient dynamics . a discussion and analysis of the phenomenology is made through an application to the specific cases of harmonic as well as anharmonic coupled oscillators , in particular the hnon - heiles system .


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: this paper is the second of two detailing an observational study of spiral structure for galaxies in the _ spitzer _ infrared nearby galaxies survey ( sings - kennicutt et al . this work also builds directly on an earlier study of m81 presented in kendall et al . ( 2008 ) . in the first paper ( kendall et al . 2011 , hereafter paper i ) , the methods and results from individual galaxies were presented together with a comparison between those galaxies in which grand design spiral structure was detected and those that lacked such structure . in this paper we focus on the subset of @xmath1 galaxies ( henceforth the ` detailed sample ' ) in which it has been possible to characterise the properties of the spiral structure in detail and present an analysis of any trends within this sample . observational results have been crucial to the advances in the study of spiral structure , and the launch of the _ spitzer _ space telescope provided a new resource in its high quality infrared data . this paper aims to build on previous studies in the optical and infrared , and provides a new dataset against which numerical simulations of spiral structure can be tested . there have been numerous observational studies of spiral structure , some key results are as follows : schweizer ( 1976 ) analysed six galaxies using optical data to produce profiles of surface brightness against azimuth . from these data schweizer concluded that there must be an old disc component in the spiral arms of each of the galaxies studied and that the spiral arms reach amplitudes as large as 30 per cent relative to the disc . other studies have confirmed this result ; the measured ratio of arm to disc mass varies with galaxy and the exact method and wavelength used , but as a rule the contrast between arm and inter - arm regions in the optical is in the range 20 to 100 per cent ( e.g schweizer ( 1976 ) ; elmegreen & elmegreen ( 1984 ) ; grosbol et al . ( 2004 ) ) , with more recent studies suggesting that an upper limit of 40 to 50 per cent is appropriate in the majority of cases ( note that contamination from star forming regions can cause this quantity to be over - estimated particularly in the early studies which used the _ i _ band rather than near infrared ( nir ) data ) . the amplitude of spiral arms is potentially important since it might be expected to directly affect the global star formation rate ( sfr ) in the galaxy : larger contrasts between the stellar arm - interarm regions are expected to cause the gas to shock more strongly , become more dense , and thus increase the star formation ( bearing in mind that according to the kennicutt - schmidt law ( kennicutt 1998 ) , star formation varies as @xmath2 . ) this association between star formation and spiral structure would appear to be borne out in cases like m81 and m51 , where the vast majority of star formation lies on well defined spiral arms . some studies have attempted to correlate the rates of star formation with the strength of the stellar spiral arms : for example , cepa & beckman ( 1990 ) found that star formation is triggered preferentially by spiral arms ( with a ` non - linear dependency ' ) . seigar & james ( 1998 ) found a link between sfr measured from far - infrared luminosity and _ k _ band arm strength , and found that the sfr increases with arm strength . similarly , seigar & james ( 2002 ) reported that the spiral arm strength ( in _ k _ ) correlates with h@xmath3 ( a sfr tracer ) suggesting that stronger potential variations ( and associated shocks ) lead to a larger sfr ( up to a limiting threshold , above which the sfr is constant ) . nevertheless , the picture does not always appear to be this simple : studies of flocculent galaxies showing underlying m=2 spiral structure in the mass density as traced by near infrared emission ( elmegreen et al . 1999 , thornley 1996 ) clearly show that star formation is not always preferentially located on two well defined spiral arms . more recently elmegreen et al ( 2011 ) analysed nir imaging for a sample of nearby galaxies from the spitzer s@xmath4 g survey ( sheth et al . 2010 ) and found that the distinctions between grand design and flocculent arm structures seen in the visible extend to the nir structure as well . another property of spiral galaxies that has been extensively studied is the pitch angle ( i.e. the angle between spiral features and the tangential direction ) : kennicutt & hodge ( 1982 ) measured pitch angles and compared them against theoretical predictions from both qsss ( quasi - stationary spiral structure - lindblad 1964 ; lin & shu 1964,1966 ) and sspsf ( stochastic self - propagating star formation - mueller & arnett 1976 ; gerola & seiden 1978 ) theories and found qualitative trends that were in agreement with both theories but poor detailed agreement in both cases . similarly , kennicutt ( 1981 ) compared pitch angle against rsa type ( largely dependent on disc resolution ) , and found a good correlation on average , but with a large degree of scatter . a correlation between pitch angle , _ i _ , and bulge - to - disc ratio was also found , where a larger pitch angle ( more open spiral pattern ) is associated with a smaller bulge fraction ; again there is a large scatter in the pitch angle at given bulge to disc ratio . in addition the maximum rotational velocity was found to correlate strongly with pitch angle , with galaxies with large rotational velocities having more tightly wound spirals .    however , there are indications that the appearance of spiral arms in the nir and optical do not always correlate ( an extreme case being the optically flocculent galaxies with nir arms already noted above ) . some studies using nir data ( e.g. seigar & james ( 1998 ) have not found the same trends as kennicutt . block & wainscoat ( 1991 ) found evidence that the spiral arms are less tightly wound in the nir than the optical . likewise , block et al . ( 1994 ) compare the nir and optical morphology for a sample of galaxies and , for a number of spirals , found significant differences between the pitch angles of the arms in the optical and nir . they found that galaxies generally appear to be of an earlier hubble type when viewed in the nir ( _ k _ band ) than in the optical , the traditional wavelength for such classifications . on the other hand , more recent studies have found a generally good agreement between the pitch angles measured in the nir and b band ( seigar et al . 2006 , davis et al . 2012 ) ; in the few cases where the two values are significantly different , there is a mild preference for the nir pitch angles to be larger . turning now to correlations with galactic parameters , seigar & james ( 1998 ) found no evidence of a correlation between nir pitch angle and bulge fraction or hubble type . seigar et al . ( 2005 ) however reported a convincing connection between pitch angle in the nir and the morphology of the galactic rotation curve , with open arms correlating with rising rotation curves , while pitch angles decline for flat and falling rotation curves . the goal of this paper is to use the new measurements of the properties of the underlying _ stellar _ spiral arms from paper i to re - examine the degrees to which the strengths and shapes of these arms are driven by the structural and environmental properties of their host galaxies . this is done mainly by examining correlations between the amplitude and pitch angles of the stellar arms with measures of galaxy morphology ( independent of spiral structure ) and mass , and with a measure of tidal perturbation from nearby companion galaxies . the same data are also used here to compare the shapes of the stellar arms with those traced by young gaseous and stellar components , and to constrain the degree to which the response of the gas discs and the local star formation rates are influenced by the amplitude of the spiral perturbation in the stellar discs . the data and methods used in this work are described in detail in paper i. a brief summary is as follows ; 3.6 and 4.5@xmath5 m nir data from the _ spitzer _ infrared array camera ( irac ) are used in conjunction with complementary optical data in order to trace the stellar mass . in addition , irac 8@xmath5 m data are used to trace the gas response , with particular emphasis on shocks . the data were analysed using azimuthal profiles ( defined from fitting isophotal ellipses to the axisymmetric components of the galaxies ) or radial profiles in the case that difficulties were encountered using azimuthal profiles . azimuthal profile data were further examined by calculating the fourier components in order to examine the amplitude and phase of individual modes , particularly the m=2 component . this fourier analysis was not possible for radial profile data ; instead the total amplitude and phase were analysed . in this paper we consider only the @xmath1 galaxies that constituted the so - called ` detailed sample ' discussed in paper i. we list various properties of these galaxies ( which we go on to correlate with the properties of their spiral structure ) in table 1 . the methods developed to isolate , identify , and parametrise the massive stellar arms are not well suited to strongly barred galaxies ( which can distort arm shapes very significantly from a logarithmic shape ) , so this paper focuses mainly on normal and weakly barred galaxies ; the role of bars in spiral structure is another important question awaiting further study ( see elmegreen et al . 2011 ) . [ cols="<,<,<,<,<,<,<,<,^",options="header " , ]         despite these caveats , figure 17 demonstrates that the amplitude of spiral responses does indeed correlate positively with the tidal @xmath6 parameter ( in fact figure 17 shows a more convincing correlation than any of the other plots we present here ) . the fact that @xmath7 galaxies have @xmath8 and yet display a significant spiral pattern implies that tidal interaction _ is clearly not a pre - requisite _ for inducing spiral structure ( and indeed we have seen in paper i that at least half the sample of galaxies without detected grand design structure have @xmath6 values of @xmath9 or above ) . on the other hand , figure 17 presents convincing evidence that where there is strong interaction the _ amplitude of spiral structure is enhanced_. we have seen that the pitch angles of such ` induced ' spirals are indistinguishable from spiral patterns in apparently isolated galaxies . we have also seen that where we have tentatively identified trends in the amplitude of spiral features ( i.e. the positive correlations with @xmath10 and stellar mass shown in figures 4 and 5 and the mild anti - correlation with gas content shown in figure 10 ) then the strongly interacting galaxies ( i.e. those with @xmath11 ) roughly follow these trends but at enhanced amplitude . once again , we however caution against over - interpretation given our limited sample size . [ hubb_m3m2 ]    [ hub_pitch2 ]    [ rotgrad_amp2 ] the first half of this paper studied the effects of galaxy type , morphology and environment on the stellar spiral structure . however , as well as being influenced by the galaxy , a global spiral is expected to influence the galaxy in which it is located . the following section considers the effect of the spiral on the structure of the gas and the distribution of star formation in galaxies . here we briefly recapitulate the results of paper i with regard to the location of shocks in the gaseous component ( as traced by @xmath12 m emission ) with respect to the spiral in the stellar mass density ( as traced at @xmath13 m ) . of the @xmath1 galaxies in the detailed sample considered here , @xmath14 galaxies ( ngc 2403 , ngc 2841 , ngc 6946 , ngc 7793 ) exhibited a flocculent response at @xmath12 m so it is not possible to talk meaningfully about the azimuthal alignment between the spirals manifest in the gas and the stars . amongst the remaining @xmath15 galaxies there is a general tendency for the @xmath12 m spiral to be located on the trailing ( concave ) side of the stellar ( @xmath13 m ) spiral , i.e. ( at radii inwards of corotation ) being _ upstream _ of the stellar spiral . moreover the @xmath12 m spiral is generally somewhat more tightly wound . this is most clearly demonstrated in the case of ngc 3031 ( m81 ) where the azimuthal offset increased systematically with radius and whose behaviour was the subject of detailed analysis in kendall et al . four other galaxies ( ngc 3184 , ngc 3198 , ngc 3938 and ngc 4579 ) also demonstrated evidence of a systematic increase of offset angle with radius . in the remaining cases , the scatter in azimuthal offset is sufficiently large to prevent the identification of an obvious radial trend . the possible significance of an upstream azimuthal offset between the shock in the gas and the maximum stellar density is discussed in kendall et al . a scenario that predicts such an offset ( with magnitude increasing with radius ) is the case of a rigidly rotating spiral mode whose lifetime is sufficiently long to permit the gas to achieve a steady state flow in the imposed potential . in this case the solutions of roberts ( 1969 ) and shu et al . ( 1972 ) predict that the gas shocks upstream of the gas inwards of corotation , with the magnitude of this offset increasing with radius , a result that has been verified through hydrodynamic simulations ( gittins & clarke 2004 ) . kendall et al . ( 2008 ) used these results to infer the corotation radius in ngc 3031 and obtained a value that is in reasonable accordance with other estimates in the literature . on the other hand , clarke & gittins ( 2006 ) and dobbs & bonnell ( 2008 ) undertook hydrodynamic simulations of gas in the case that the stellar potential derived from a high resolution re - simulation of the n - body calculations reported in sellwood & carlberg ( 1984 ) . notably in these simulations with a ` live ' galactic potential , the gas shocks tended to trace the regions of instantaneous maximum stellar density and exhibited no systematic azimuthal offset with respect to the stellar spiral . a clear difference between this simulation and the idealised case involving a single long lived spiral mode is that , in the n - body calculations , spiral features come and go on a timescale of a few orbits . clarke & gittins related the gas morphology in this case to the fact that the gas did not have time to adjust to a steady state configuration on the timescale on which individual spiral features established and dissolved . in the case of spiral structure induced by a strong galactic encounter ( as in the modeling of m51 - ngc 5194 - by dobbs et al . 2010 ) no azimuthal offset was found between the gas and stellar spirals . how do these theoretical results relate to the range of results that we find in our sample ? we can divide our sample ( as detailed above ) into @xmath14 galaxies with flocculent @xmath12 m structure , @xmath16 with evidence for a systematic upstream offset with magnitude broadly increasing with radius and the remaining @xmath14 as having no simple offset pattern ( we denote these groups f(locculent ) , r(egular ) and c(omplex ) : see paper i for the phase diagrams on which this classification is based . we find no systematic differences between the galactic properties in these three groups apart from the fact that all the galaxies in the f(locculent ) category are isolated ( @xmath8 ; see table 2 ) ; the r(egular ) and c(omplex ) galaxies occupy a wide range of @xmath6 values , although notably the galaxy which is undergoing by far the strongest interaction ( ngc 5194 ) is in the c(omplex ) category , in line with the simulations of this encounter by dobbs et al . ( 2010 ) . we conclude from this that many of the isolated galaxies exhibit the kind of gas response seen in simulations based on live galactic potentials and that the most strongly interacting galaxy also agrees with simulations inasmuch as it exhibits no clear offset ; the r(egular ) group however exhibit a trend suggestive of a more long lived spiral pattern . we re - emphasise that this group is heterogeneous in terms of its galactic properties . to date there has been no detailed quantitative study of how tightly the offset data can be used to constrain the lifetime and nature of spiral structure in simulations . in addition to the position of the shock , it is instructive to look at the strength of the gas shocks as a function of the strength of the stellar spiral , since a larger spiral amplitude is expected to trigger a larger shock . we compared these properties for all galaxies in the sample , not just those for which the offset work was carried out . in contrast to the stellar mass maps , the gas response can not be normalised by an axisymmetric component , so instead we simply divide the 8@xmath5 m azimuthal profile into ` on ' and ` off ' arm regions ; ` on ' arm is defined as being angles @xmath17 of the m=2 peaks ( ` off ' arm is the other @xmath18 degrees , in two arcs ) . the response is calculated by taking the ratio of the highest peaks in the ` on ' arm region to the average ` off ' arm flux . note that we use azimuthal profile data for analysis of the @xmath12 m emission for all galaxies in the sample , including those whose amplitude at @xmath13 m is obtained through radial profile analysis . m gas shock strength indicator ( defined as the average of the four highest peaks in the 8@xmath5 m profile in the ` on ' arm region , normalised by the average value of the 8@xmath5 m profile in the ` off ' arm region ) plotted against average m=2 amplitude at @xmath13 m . see figure 2 for an explanation of the symbols.,width=313 ]    m excess ratio ( see text for definition ) . see figure 2 for an explanation of the symbols.,width=313 ]    figure 18 shows the trends in the gas response as a function of m=2 amplitude and show that , as expected , the gas response will tend to be stronger for a larger stellar wave amplitude . ( note that we expect such a correlation to be weakened somewhat by the fact that @xmath12 m pah features may be depleted by the strong uv emission associated with the enhanced star formation in shocked gas . )      in this section we apply the same method as in section 4.1.2 to define ` on ' and ` off ' arm regions , using the _ spitzer _ mips ( rieke et al . 2004 ) 24@xmath5 m wavelength as a tracer of star formation ( calzetti et al . @xmath19 m emission is dominated by reprocessed radiation from young massive stars as demonstrated by its excellent spatial correlation with other diagnostics associated with such stars ( see relano & kennicutt 2009 for such a demonstration in the case of h@xmath3 emission in m33 ) . helou et al . ( 2004 ) argued that the stellar continuum contributes at most a few per cent at @xmath19 m    in figure 19 we plot the ratio of the _ excess _ emission at 24@xmath5 m in the arm compared with the inter - arm region ( where the excess is defined as the mean emission levels in that region minus the minimum emission level in that radial ring : note that in the case of the @xmath12 m data this minimum level is essentially zero in many regions and so in contrast figure 18 presents the ratio of unsubtracted values . ) figure 19 is qualitatively similar to figure 18 and shows a clear dependence of the star formation response ( as measured at @xmath19 m ) and the amplitude of variations in the stellar potential . the levels of contrast can be used to interpret figure 11 , noting that the dispersion in specific star formation rates and in the @xmath19 m excess ratios are comparable . this suggests that strong spiral structure can contribute significantly to the dispersion in specific star formation rates in our sample . nevertheless , it can not be the only driver since there are some notable counter - examples ( for example , ngc 3031 with its strong @xmath19 m excess ratio and low specific star formation rate . by far the strongest correlation that we have identified ( figure 17 ) is between the _ amplitude _ of spiral structure and the strength of tidal interaction as measured by the @xmath6 parameter ( equation ( 1 ) ) . we conclude - in line with our conclusions in paper i - that a close tidal encounter appears to be a sufficient but not necessary condition for prominent spiral features . ngc 3198 is a good example of a galaxy that is unbarred and isolated and yet shows the ` grand design structure ' that permits its inclusion in the detailed sample studied in this paper . nevertheless , the amplitude of its spiral features is admittedly lower than the strongly interacting galaxies in our sample and moreover the @xmath20 mode is less dominant than in these galaxies . we have indicated the galaxies with the three largest values of @xmath6 ( the strongly interacting galaxies ) by the large crosses in all the other correlation plots . these galaxies are not atypical for the sample with regard to their pitch angles and the radial extent of spiral structure ( as measured by the azimuthal angle subtended by the arms in figures 8 and 9 ) and so do not strongly influence these plots . they are however all rather massive and gas deficient galaxies of type sbc or earlier and we must therefore be mindful that apparent correlations may be driven by this association . this may be particularly relevant to the interpretation of figures 4 , 5 and 10 . our sample covers a rather narrow range of pitch angles and it is thus unsurprising that we have not identified any correlations between pitch angles and other galactic parameters . figure 12 may throw some light on previous comparisons between the pitch angles of spiral structure in the optical and nir : for example in the most recent and sophisticated analysis of this issue by davis et al . ( 2012 ) there is generally good agreement between pitch angles in the b and i bands . discrepant galaxies have a tendency to have larger pitch angles at longer wavelengths . figure 12 ( which compares nir pitch angles with pitch angles measured in h@xmath3 , an optical star formation indicator ) suggests that in our sample a discrepancy of this sign is evident only for earlier hubble types . this would suggest - if borne out in larger samples - that an association between early hubble types and tightly wound arms may only be manifest in the optical ( clarke et al . 2010 ) . in the most discrepant galaxies ( ngc 3031 and ngc 4579 ) this difference in pitch angles can be explained by our previous demonstration ( kendall et al . 2008 , 2011 ) that in these galaxies the spiral pattern traced at @xmath12 m in shocked gas is significantly displaced upstream with respect to the nir ( stellar ) arms . in a differentially rotating galaxy it would then be expected that the spiral traced by star formation indicators should be more tightly wound than the stellar arms . the origin of azimuthal offsets is discussed in our previous papers ; it is currently unclear whether there is a more general association between such offsets and galaxies of early hubble type . we also draw attention to the other strong correlation that we have found , i.e. the demonstration in figures 18 and 19 of a clear dependence of the amplitudes of spiral structure in star formation ( @xmath19 m ) and shock ( @xmath12 m ) tracers on the amplitude of potential variations as traced in the nir . the amplitude of spiral structure at @xmath19 m demonstrates that the star formation within the arms can not on its own account for all the variation in specific star formation rate exhibited by the galaxies in our sample and there are indeed objects with strong arm contrast yet low specific star formation rates . however the range of arm amplitudes at @xmath19 m suggest that variation of strength of spiral features plays a significant contributory role in setting the specific star formation rates in disc galaxies . many thanks to bob carswell for much useful advice regarding the use of idl , in particular the curvefit program . we are indebted to hans - walter rix for the idea of using colour - correction to make mass maps from optical data . many thanks also go to jim pringle , clare dobbs , stephanie bush , giuseppe bertin and jerry sellwood for useful discussions over the course of this work and to phil james and marc seigar for valuable observational input . we also thank the referee for constructive comments . this work makes use of iraf . iraf is distributed by the national optical astronomy observatories , which are operated by the association of universities for research in astronomy , inc . , under cooperative agreement with the national science foundation . baba , j. , saitoh , t. , wada , k. , 2013 . apj 763,46 bendo , g. , calzetti , d. , englebracht , c. et al . , 2007 . mnras 380 , 1313 bertin , g. , lin , c. , 1996 . in spiral structure in galaxies a density wave theory , publisher : cambridge , ma mit press . block d. l. , bertin g. , stockton a. , grosbol p. , moorwood a. f. m. , peletier r. f. , 1994 , a&a , 288 , 365 block d. l. , bertin g. , stockton a. , grosbol p. , moorwood a. f. m. , peletier r. f. , 1994 , a&a , 288 , 365 . block d. l. , puerari i. , frogel j. a. , eskridge p. b. , stockton a. , fuchs b. , 1999 , ap&ss , 269 , 5 block d. l. , wainscoat r. j. , 1991 , nature , 353 , 48 byrd g. g. , howard s. , 1992 , aj , 103 , 1089 calzetti , d. , kennicutt , r. , engelbracht , c. et al . , 2007 . ap j 666 , 870 cepa j. , beckman j. e. , 1990 , apj , 349 , 497 clarke , c. , gittins , d. , 2006 . mnras 371,530 clarke , c. , kendall , s. , kennicutt , r. , 2010 in aip conf . proceedings 1242,129 davis , b. , berrier , j. , shields , d. et al . , 2012 . apj s 199,33 dobbs c. l. , bonnell i. a. , 2008 , mnras , 385 , 1893 dobbs , c.l . , theis , c. , pringle , j. , bate , m. , 2010 . mnras 403 , 625 elmegreen d. m. , chromey f. r. , bissell b. a. , corrado k. , 1999 , aj,118,2618 elmegreen d. m. , elmegreen b. g. , 1984 , apjs , 54 , 127 - 1987 , apj , 314 , 3 elmegreen , d.m . et al 2011 , apj , 737 , a32 fuchs , b. , 2000 . in dynamics of galaxies from the early universe to the present , 15th iap meeting , eds . combes , f. , mamon , g. , charmandaris , v. , p.53 gerola h. , seiden p. e. , 1978 , apj , 223 , 129 gittins d. m. , clarke c. j. , 2004 , mnras , 349 , 909 grand , m. , kawata , d. , cropper , m. , 2013 . a & a 553 , a77 grosbol p. , patsis p. a. , pompei e. , 2004 , a&a , 423 , 849 helou , g. , roussel , h. , appleton , p. et al . , 2004 . apjs 154,253 kendall s. , kennicutt r. c. , clarke c. , thornley m. d. , 2008 , mnras , 387 , 1007 kendall , s. , kennicutt , r. , clarke , c. , 2011 . mnras 414,538 ( paper i ) kennicutt , jr r. c. , hodge p. , 1982 , apj , 253 , 101 - 1981 , aj , 86 , 1847 kennicutt jr . r. c , 1998 , apj , 498 , 541 kennicutt jr . r. c. , armus l. , bendo g. , et al . pasp 115,928 kormendy j. , norman c. a. , 1979 , apj , 233 , 539 leroy , a. , walter , f. , brinks , e. , et al . 2008 . aj 136,2782 lin c. c. , shu f. h. , 1964 , apj , 140 , 646 - 1966 , proceedings of the national academy of science , 55,229 lindblad b. , 1964 , astrophysica norvegica , 9 , 103 ma , j. , 2001 . chjaa , 1 , 395 michikoshi , s. , kokubo , e. , 2014 . apj 787,174 . mueller m. w. , arnett w. d. , 1976 , apj , 210 , 670 oh s. h. , kim w .- t . , lee h. m. , kim j. , 2008 , apj , 683 , 94 perez - gonzalez p. g. , kennicutt jr . r. c. , gordon k. d. , et al . , 2006 . apj 648,987 relano , m. , kennicutt , r. jr . , 2009 . apj 699,1125 rieke g. h. , young e. t. , engelbracht c. w. , et al . apj s 154,25 roberts , w. w. , 1969 , apj , 158 , 123 roberts jr . w. w. , roberts m. s. , shu f. h. , 1975 , apj , 196 , 381 schweizer , f. , 1976 . apjs , 31 , 313 seigar m. s. , block d. l. , puerari i. , chorney n. e. , james p. a.,2005 , mnras , 359.1065 seigar , m. , bullock , j. , barth , a. , ho , l. , 2006 . apj 645,1012 seigar m. s. , chorney n. e. , james p. a. , 2003 , mnras , 342 , 1 seigar m. s. , james p. a. , 1998 , mnras , 299 , 685 - , 2002 , mnras , 337 , 1113 sellwood , j. , carlberg , r. , 1984 . apj 282,61 - , 2014 . apj 785,137 sheth , k. et al . 2010 , pasp , 122 , 1397 shu , f. , millone , v. , gebel , w. , yuan , c. , goldsmith , d. , roberts , w. , 1972 . apj 173,557 thornley m. d. , 1996 , apj , 469 , l45 toomre a. , 1981 , in structure and evolution of normal galaxies , fall s. m. , lynden - bell d. , eds . 111 - 136 tully r. b. , fisher j. r. , 1977 , a&a , 54 , 661 van den bergh , s. , 1960a , apj 131,215 van den bergh , s. , 1960b .
Abstract:

Answer: this paper presents a detailed analysis of two - armed spiral structure in a sample of galaxies from the spitzer infrared nearby galaxies survey ( sings ) , with particular focus on the relationships between the properties of the spiral pattern in the stellar disc and the global structure and environment of the parent galaxies . following paper i we have used a combination of spitzer space telescope mid - infrared imaging and visible multi - colour imaging to isolate the spiral pattern in the underlying stellar discs , and we examine the systematic behaviours of the observed amplitudes and shapes ( pitch angles ) of these spirals . in general , spiral morphology is found to correlate only weakly at best with morphological parameters such as stellar mass , gas fraction , disc / bulge ratio , and v@xmath0 . in contrast to weak correlations with galaxy structure a strong link is found between the strength of the spiral arms and tidal forcing from nearby companion galaxies . this appears to support the longstanding suggestion that either a tidal interaction or strong bar is a necessary condition for driving grand - design spiral structure . the pitch angles of the stellar arms are only loosely correlated with the pitch angles of the corresponding arms traced in gas and young stars . we find that the strength of the shock in the gas and the contrast in the star formation rate are strongly correlated with the stellar spiral amplitude . [ firstpage ]    galaxies : individual  galaxies : spiral  galaxies : structure  infrared : galaxies