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: synchronization of neural activity appears in different parts of the mammalian cerebral cortex @xcite , and underlies different neural processes in both normal and anomalous brain functions @xcite . it has been suggested that synchronization plays a vital role in information processing in the brain , e.g. , processing information from different sensory systems to form a coherent and unified perception of the external world @xcite . on the other hand , synchronization has been detected in pathological conditions such as parkinson s disease @xcite . and epileptic seizures have long been considered resulting from excessive synchronized brain activity @xcite , though some recent studies suggest that this picture may be an over - simplification @xcite . therefore understanding the mechanisms of synchronization may be a critical step in elucidating how neural systems work @xcite . it has stimulated a great deal of theoretical and numerical works , such as the studies on the effects of the topological properties of underlying networks @xcite and the dynamical properties of synaptic coupling @xcite . it was recently shown that the response time of synaptic couplings influences the stability of synchronized oscillation in the nonlocally coupled hodgkin - huxley ( hh ) equations @xcite . if the response time of synaptic coupling is slower , synchronized activity of the systems is instable for excitatory coupling . however , the underlying dynamical mechanism of the influence is not clear . in experimental studies @xcite , it has been suggested that the generation of prolonged epileptiform neuronal synchronization is favored by lower efficacy of synaptic transmission . the numerical studies @xcite in a detailed computational model revealed that seizure - like activity occurs when the excitatory synapses are weakened , and the results were confirmed experimentally in mouse neocortical slices . according to the common accepted assumption that synchronization of neuronal activity underlies seizures , the dynamical mechanism of synchronization may be useful for understanding the way the biological neural system works . in this work , we numerically investigated the dynamical mechanism underlying the influence of synaptic efficacy on firing synchronization in hh neuron networks . to do this , we first studied the dynamics of the response of hh neuron to excitatory synaptic current . when the efficacy of synapse is low , namely , strength is weak and duration is short , the limit cycle is stable to the perturbation of the synaptic current . when synaptic efficacy is high , synaptic current can induce the transition of the neurons from limit cycle to fixed point or transient state . the transition is determined by dynamics of neuron s ionic channel . the decrease of synaptic current depresses the feedback of sodium ionic current which is responsible for the initiation of the spike . for simplicity we will refer to the transitions as spike death .    in neuronal networks , synaptic input of a neuron is the accumulation of the currents received from all presynaptic neurons . when the coherence of firing time of neurons is enhanced by the excitatory interaction , the synaptic input of neurons transforms from the fluctuating waveform into the pulse shape like the signal produced by one synapse . if synaptic efficacy is high , the input signal can induce spike death of the neuron . then spike death disorders the adjustment of the rhythm of neurons and prevents neurons from firing spikes synchronously . in contrast , for synapses of lower efficacy , the duration of synaptic current is too short to induce spike death of neurons . additionally , the firing synchronization is different from synchronous activity of oscillators for the existence of the transitions of neuron s state . the paper is organized as follows . the hh neuron model and the synaptic coupling were introduced in sec . the response of a hh neuron to synaptic current was investigated in sec . the influence of the dynamics of neurons on firing synchrony was shown in sec . discussion and conclusion were given in sec . to investigate the dynamics of a neuron under the perturbation of synaptic stimulus , we adopted a system consisting of a hh neuron and a synapse . the hh neuron was originally proposed for the giant axon in a squid @xcite . it serves as a paradigm for the spiking neuron models that based on nonlinear conductances of ion channels . the model describes the evolution of the membrane potential @xmath0 and can be written as @xmath1 where @xmath2 is the ionic current , @xmath3 is the external current and @xmath4 is the synaptic current . the ionic current describes the ion channel on the membrane and is defined as @xmath5 where @xmath6 , @xmath7 and @xmath8 are the maximum conductances for the sodium , potassium and leak currents , @xmath9 , @xmath10 and @xmath11 are the corresponding reversal potentials . @xmath12 and @xmath13 are the activation and inactivation variables of the sodium current and @xmath14 is the activation variable of the potassium current . the gating variables @xmath15 @xmath12 , @xmath13 , @xmath14 satisfy the differential equation @xmath16-\beta_yy(t ) , \label{eq : dy}\ ] ] with nonlinear functions @xmath17 and @xmath18 given by @xmath19\ } , \\ \beta_m & = & 4\exp[-(v+65)/18 ] , \\ \alpha_h & = & 0.07\exp[-(v+65)/20 ] , \\ \beta_h & = & 1/\{1+\exp[-(v+35)/10]\ } , \\ \alpha_n & = & 0.01(v+55)/\{1-\exp[-(v+55)/10 ] \ } , \\ \beta_n & = & 0.125\exp[-(v+65)/80].\end{aligned}\ ] ] the parameter values are @xmath20 mv , @xmath21 mv , @xmath22 mv , @xmath23 ms/@xmath24 , @xmath25 ms/@xmath24 , @xmath26 ms/@xmath24 , and @xmath27 @xmath28f/@xmath29 @xcite . the external current @xmath3 determines the firing rate of the neuron . in the absence of synaptic coupling @xmath4 , the hh neuron has the following bifurcation diagram as a function of @xmath3 : in the parameter regions @xmath30 and @xmath31 fixed point is the global attractor . for @xmath32 the neuron possesses coexisting stable attractors , the fixed point and the limit cycle which are separated by an unstable limit cycle . for @xmath33 fixed point becomes unstable . the values of bifurcation points are @xmath34 @xmath28a/@xmath24 , @xmath35 @xmath28a/@xmath24 and @xmath36 @xmath28a/@xmath24 @xcite . we adopted the synaptic current @xmath4 described by an alpha function @xcite . the alpha function synapse is a phenomenological model based on an approximate correspondence of the time course of the waveform to physiological recordings of the postsynaptic response @xcite . the equation of the synapse is like : @xmath37 with @xmath38 where @xmath39 is the characteristic time of the interaction , @xmath40 is the heaviside step function , and @xmath41 is the beginning time of the synaptic interaction , i.e. , the firing time of the presynaptic neuron ( all delays are neglected ) . the synaptic effect is traditionally classified as excitatory or inhibitory depending on the value of @xmath42 . here , we took @xmath43 mv for excitatory synapses , and @xmath44 mv for inhibitory ones . ( [ alpha ] ) yields pulses with the maximum value of @xmath45 at @xmath46 and with the half - width of 2.45@xmath39 @xcite . so @xmath39 characterizes the duration of the synaptic interaction . for the alpha function synapse , we equated the synaptic efficacy to the maximum synaptic conductances @xmath47 and the characteristic time @xmath39 . the high efficacy means that synaptic current possesses strong strength and long duration . ( color online ) ( a ) the response of bistable hh neuron to excitatory synaptic current . the symbol @xmath48 denotes the sum of the external current @xmath3 and the synaptic current @xmath4 . ( b ) the corresponding phase portrait . the values of parameters are @xmath49 msec , @xmath50 ms/@xmath24 and @xmath51 @xmath28a/@xmath24 . ] we focused on the dynamics of the system in the parameter region near the bifurcation point @xmath52 . firstly , we studied the response of bistable neuron ( @xmath53 ) to excitatory synaptic current . in simulations , firing was identified as membrane potential @xmath54 is over 20 mv . when the neuron fired a spike , we triggered a pulse of synaptic current into it . we observed two types of dynamics of the response , which depended on the efficacy of synapse . for slow response time and strong synaptic strength , the neuron transited from the limit cycle to the fixed point . the transition of neuron s state is shown in figs . [ fig1 ] . in fig . [ fig1](a ) the response of the neuron to the synaptic current is represented by the membrane potential @xmath54 . in the figure , the synaptic current is illustrated by the pulse added on the external current . one can see that the periodic firing eliminated after the synaptic current was injected . the value of membrane potential tended to -61.15 mv through subthreshold oscillation . in fig . [ fig1](b ) the transition of neuron s state is shown in a three - dimensional space ( @xmath54 , @xmath13 , @xmath12 ) which is a projection of the phase space ( @xmath54 , @xmath13 , @xmath12 , @xmath14 ) . the trajectory left the limit cycle and was attracted into the basin of the fixed point . at last , through transient process , the trajectory stopped at the fixed point which is indicated by the dashed lines in the figure . the fixed point was @xmath55=(-60.15 , 0.423 , 0.092 , 0.394 ) as the parameter @xmath51 @xmath28a/@xmath24 . on the other hand , for systems with quick response time and weak synaptic strength , the transition did not occur and the trajectory was attracted back the limit cycle from weak perturbation . ( color online ) ( a ) the response of the neuron with @xmath56 @xmath28a/@xmath24 to excitatory synaptic current . ( b ) the corresponding phase portrait . the synaptic parameters are the same as fig . [ fig1 ] . ] when the external signal @xmath3 is larger than but near @xmath52 , the neuron possesses a stable limit cycle and an unstable fixed point . in this case there were also two types of dynamics of response . for high synaptic efficacy , the neuron exhibited a transient behavior when it received the synaptic current . [ fig2 ] ( a ) shows that the membrane potential responded to the synaptic current with a transient subthreshold oscillation . the subthreshold oscillation interrupted the periodic firing of the neuron . in phase space the transient behavior was a motion around the unstable fixed point . this is explicitly shown in the corresponding 3d phase space ( @xmath54 , @xmath13 , @xmath12 ) . in fig . [ fig2](b ) one can see that the trajectory of the neuron left the limit cycle and transiently moved around the unstable fixed point which is indicated by the dashed lines . the unstable fixed point was @xmath55=(-58.704 , 0.374 , 0.108 , 0.417 ) as the parameter @xmath56 @xmath28a/@xmath24 . on the other hand , for low synaptic efficacy , the coupling can not induce the transient motion around the unstable fixed point . like bistable neurons , the trajectory returned to limit cycle from weak perturbation . [ fig3 ] shows the boundary between the two types of dynamics on the parameter plane @xmath47 vs @xmath39 . for the bistable neuron of the external current @xmath51 @xmath28a/@xmath24 , the boundary between the two types of dynamics is represented by squares . above the curve , neurons responded to synaptic currents with transitions between attractors . inversely , the limit cycle of neurons was stable to the perturbation . for the neuron of the external current @xmath57 @xmath28a/@xmath24 , the boundary is shown by circles . it is notable that stronger strength and longer duration of synaptic current were needed by the transitions of neuron s state . based on numerical simulations , we obtained that the transitions of neuron s state can be induced by synaptic current in the parameter region @xmath58 @xmath28a/@xmath24 @xmath59 @xmath28a/@xmath24 . above the upper boundary of the region , the transient motion around the unstable fixed point can not be induced by the synaptic current . ( color online ) the dynamics diagram of minimal synaptic intensity @xmath47 vs the characteristic time @xmath39 for the neurons with the external current @xmath51 @xmath28a/@xmath24 ( squares ) or @xmath56 @xmath28a/@xmath24 ( circles ) . ]    in the following , we referred to the phenomena shown in both fig . [ fig1 ] and fig . [ fig2 ] as spike death . different from oscillator death which is the quenching of oscillation of coupled systems , spike death is the behavior of a single neuron . and spike death includes both the transition between stable attractors and the transient behavior . next we gave a qualitative interpretation of the spike death phenomena . although the synaptic current consists of a rise and a decay stage , we want to show that the decrease of the synaptic current induces the transition of neuronal activity . we illustrated the role of the decreased current in fig . [ fig4 ] . in the simulations , external currents linearly decreased from 20 @xmath28a/@xmath24 to 8.5 @xmath28a/@xmath24 , and synaptic current was absent . the time of beginning to decrease and the different slopes of the decrease current were chosen to represent the distinct duration and rate of decay of the current . in fig . [ fig4 ] ( a ) , the activity of the neuron transited from the periodic firing to silent state through transient subthreshold oscillation . the decrease of current induced the spike death of the neuron . in fig . [ fig4 ] ( b ) , the decay of the current occurred in the refractory period of the neuron . the decreased current did not depress the next spike . this is similar to the scenario of the synaptic current of short duration . so spike death requires that the decrease of current occurs at the end of refractory period and the stage of the initiation of next spike . in fig . [ fig4 ] ( c ) , the external current decreased with the slope -0.5 @xmath28a/(@xmath24 msec ) . comparing with fig . [ fig4 ] ( a ) , the rate of decrease was small . although the current decreased at the stage of initiating a spike , it did not depress the firing of the neuron . thus spike death requires that the rate of decrease of current is large at the stage of initiating spikes .    to interpret the effect of the decrease of current on the oscillation of neurons , we reviewed the generation of spikes . for a fixed value of membrane potential @xmath54 , the variable @xmath60 (= @xmath12 , @xmath13 , @xmath14 ) approaches the value @xmath61 $ ] with the time constant @xmath62^{-1}$ ] . the variable @xmath63 increases with @xmath54 , and the corresponding time scale @xmath64 is littler than @xmath65 and @xmath66 . if external current injects into the cell and raises the membrane potential @xmath54 , the conductance of sodium channels @xmath67 increases due to increasing @xmath12 . then sodium ions flow into the cell and raise membrane potential even further . if this positive feedback is large enough , a spike is initiated @xcite . when external current decreases at the onset of the generation of a spike , the raise of membrane potential slows down . then the increase of @xmath12 is slowed down . so the positive feedback is weakened . if the current decreases quickly , the spike of neuron can be depressed . therefore , for depressing the positive feedback of sodium ionic current , the decrease of current must occur at the onset of generation of spike , and the current must decrease in a large rate . this interpretation is consistent with above simulated results that the strong strength and the long duration of synaptic current are necessary for spike death . ( color online ) the change of the activity of hh neuron as external current decreases . the external current is reduced from 20 @xmath28a/@xmath24 to 8.5 @xmath28a/@xmath24 with the slop ( a ) -1 @xmath28a/(@xmath24 msec ) , ( b ) -2 @xmath28a/(@xmath24 msec ) , or ( c ) -0.5 @xmath28a/(@xmath24 msec ) . ] next we investigated the influence of spike death on the firing synchronization in neuronal networks . we considered a directed random network consisting of @xmath68 nonidentical hh neurons . the network was generated as following : with a probability @xmath69 , we connected each of the probable directed couplings ( such as the one from the @xmath70th neuron to the @xmath71th neuron ) . the network is described by adjacency matrix @xmath72 which entry @xmath73 is equal to 1 when the coupling from @xmath70 to @xmath71 exists , and zero otherwise . in the directed network , @xmath73 can be not equal to @xmath74 , and the signal travels in only the direction from @xmath70 to @xmath71 if @xmath73=1 and @xmath74=0 . the signal received by the neuron @xmath71 is the accumulation of all input synaptic currents , which is defined as @xmath75 where @xmath76 is the rescaled factor which equals the number of inputting synapses of neuron @xmath71 ( @xmath77 ) . @xmath78 is the latest firing time of presynaptic neuron @xmath70 . @xmath79 is the reverse potential of the synapse connecting @xmath70 to @xmath71 .    to study the global behavior of neuronal networks we computed the average activity @xmath80 of the network . the amplitude of average activity can intuitionally reveal the coherence of activity of neurons , which is defined as @xcite @xmath81 ^ 2 dt,\ ] ] where the angle brackets denote temporal average over the integration interval . in this work we studied the firing synchrony for the reason that neuronal states may transit away from the limit cycle and the synchronization of oscillators will be disrupted . here we focused on the coherence of firing time of neurons . we adopted the average cross correlation of firing time of neurons @xcite to quantify the degree of firing synchronization . average cross correlation is obtained by averaging the pair coherence @xmath82 between neuron @xmath71 and @xmath70 , i.e. , @xmath83 the pair coherence @xmath82 is defined as @xmath84^{\frac{1}{2}}},\ ] ] which is measured by the cross correlation of spike trains at zero time lag within a time bin @xmath85 . to transform the neuronal activity into spike train , the interval @xmath86 is divided into @xmath87 bins of @xmath88 msec . then spike trains of neurons @xmath71 and @xmath70 are given by @xmath89 or 1 and @xmath90 or 1 ( @xmath91 ) , where 1 represents a spike generates in the bin and 0 otherwise . ( color online ) ( a ) the amplitude of average activity @xmath92 versus the fraction of excitatory neuron @xmath93 . squares ( cycles ) represented the relations in networks with @xmath39=2 msec ( 1 msec ) . the error bars were the standard deviation across 20 realizations . ( b ) the relation between average cross correlation @xmath94 and the fraction of excitatory neuron @xmath93 . ]    we numerically investigated the collective activity of networks versus the fraction of excitatory neurons @xmath93 in the networks . in general , excitatory synapses tend to synchronize the activity of neurons in networks . however , it was shown that in neural systems different types of dynamics of response of neurons to couplings produce different synchrony property @xcite . here we studied the influence of the new type of response , spike death , on the formation of synchrony as the excitatory interactions increase . in simulations , the networks consisted of 1000 neurons and connecting probability was @xmath69=0.01 . the synaptic conductance took the value @xmath47=1 ms/@xmath24 . as an example , we used the synapses with the characteristic time @xmath39=2 msec or @xmath39=1 msec to generate neuronal networks with or without spike death , respectively . to ensure that the neurons had nonidentical properties , the external currents @xmath95 were in @xmath96 @xmath28a/@xmath24 and were generated at random .    in fig . [ fig5](a ) we plotted the amplitude of average activity @xmath92 versus the fraction @xmath93 of excitatory neurons . squares represents the results obtained in networks with @xmath39=2 msec , and circles represents the results with @xmath39=1 msec . in fig . [ fig5](b ) we did the same for average cross correlation @xmath94 . the relations showed that the coherence of activity increases with the fraction of excitatory neurons . it is notable that the networks with @xmath39=2 msec , i.e. , with the spike death property , had obviously lower values of @xmath92 and @xmath94 than networks with @xmath97 msec . so the degree of synchrony in networks with spike death property was obviously lower than networks without the property . average membrane potential of the networks with @xmath49 msec ( left ) and @xmath97 msec ( right ) . the fraction of excitatory neurons is @xmath98 , 0.8 , 1.0 from top to bottom . ] the difference of the degree of synchrony between two types of networks is shown more intuitionally by the average activity in fig . we made neurons in the networks beginning to oscillate with a high degree of coherence . in simulations the first firing time of neurons randomly distributed in ( 0 , 5 ) msec . when @xmath93=0.5 , the average activity of the network with @xmath49 msec decreased obviously and the network changed quickly to oscillating randomly , as shown in fig . [ fig6](a ) . in the network without spike death property , coherent oscillation @xcite appeared , as shown in fig . [ fig6](b ) . when the fraction of excitatory neurons increased , the excitatory interaction enhanced the average activity in the networks both with and without the spike death property , while the average activity in the latter was obviously larger . furthermore , the frequency of the average activity @xmath99 in networks with @xmath97msec was higher than the frequency in networks with @xmath39=2 msec . ( a ) the input synaptic current of a randomly chosen neuron from the network with @xmath49 msec , @xmath98 . ( b ) the input synapse current when @xmath93 increased to 0.8 . ( c ) the histogram of spike death events in the network with @xmath49msec and @xmath100 . ]    the synchrony properties can be qualitatively explained by spike death . in networks , if the fraction of excitatory neuron is increased , the coherence of firing time of neurons is enhanced . then the input synaptic current of a neuron may transform from fluctuating waveform into a smooth pulse . if the pulse of the accumulated synaptic current has long enough duration to depress the next spike , the synaptic current may lead to spike death . therefore the synaptic current disorders the adjustment of the rhythm of neurons firing and prevents synchronization . to demonstrate this , we showed the synaptic currents and spike death events in fig . 7 . in fig . 7(a ) and 7(b ) we plotted the input synaptic currents of a neuron , which was randomly chosen in networks with @xmath93=0.5 and @xmath93=0.8 respectively . the transform from fluctuate waveform to smooth pulse was observed . in simulations we took the peak of subthreshold oscillation of a neuron as a spike death event . 7(c ) shows the histogram of frequency of spike death events , i.e. , the number of spike death events in 1 msec . the spike death events periodically appeared with the same rhythm as the average activity of the network . in contrast , spike death events can not be obtained in the networks of low synaptic efficacy . thus the spike death can explain the above synchrony property . if network consists of identical neurons , the influence of spike death also exists and is more remarkable . 8(a ) shows the relation of @xmath94 to @xmath93 of the networks in which the external currents of all neurons was @xmath3=10.0@xmath28a/@xmath24 . one can see that the values of @xmath94 of the networks with @xmath39=2 msec ( squares ) was similar to fig . 5 . for networks with @xmath39=1 msec ( circles ) , however , the value of @xmath94 tended to 1 as @xmath93 increases . we calculated the degree of synchrony as a function of the heterogeneity of neurons in purely excitatory networks . in simulations the value of external currents @xmath3 distributed in the region @xmath101 @xmath28a/@xmath24 . 8(b ) shows the relation between @xmath94 and the width @xmath102 of the parameter region of @xmath3 . for networks with @xmath49msec ( squares ) , the degree of synchrony was unsensitive to the heterogeneity of neurons . in contrast , the degree of synchrony in networks with @xmath97msec ( circles ) remarkably increased and tended to 1 as the heterogeneity of neurons decreased . for networks with @xmath97msec , the relation between @xmath94 and the width @xmath102 of the parameter region was fitted to the first order exponential decay curve . the fitted curve ( dashed line ) is @xmath103 with @xmath104 = 0.362 @xmath105 0.005 , @xmath106 = 0.595 @xmath105 0.007 and @xmath107 = 1.017 @xmath105 0.030 . the remarkable difference between two kinds of networks shows that spike death can effectively prevent the firing synchronization in neuronal networks . ( color online ) ( a ) : the average cross correlation @xmath94 versus the fraction of excitatory neurons for the networks consisting of identical neurons . squares ( circles ) represent the relation in networks with @xmath49msec ( 1msec ) . the error bars were the standard deviation across 20 realizations . ( b ) : the average cross correlation @xmath94 changes with the width @xmath102 of parameter region of @xmath3 . dashed line is the fitted exponential curve . ] we have used hh neuron networks to investigate the dynamical origin of the influence of synaptic efficacy on the firing synchronization . a new dynamics of response of neurons to coupling , spike death , was suggested as a possible mechanism underlying the influence . when the firing time of neurons is so coherent that synaptic currents have a pulse waveform in excitatory networks , synaptic current induces the transition of neuron s state from limit cycle to fixed point or transient state . the transitions disrupt the adjustment of rhythm of neurons oscillations and prevent further increase of firing synchronization . we studied the dynamics of the hh neuron responding to the excitatory synaptic perturbation . we numerically demonstrated that the synapse of high efficacy , i.e. , large characteristic time and strong strength , induces spike death of the neuron . for bistable neurons , spike death means the transition from limit cycle to fixed point . for the neuron with unstable fixed point , spike death means the transition from limit cycle to transient state . the transient state is the motion around the unstable fixed point in phase space . spike death of neurons results from the decrease of synaptic current , which depresses the feedback of sodium ionic current at the stage of initiating a spike . we demonstrated the influence of spike death on the degree of firing synchronization . in simulations we considered the networks with or without the spike death property , which were generated using synapses of characteristic time @xmath39=2 msec and @xmath39=1 msec , respectively the degree of synchrony of the former was lower . this is consistent with results of @xcite that synchronous state is not stable for the excitation coupling of slow response time . however , we also showed that for slow response time , the degree of synchrony increased with the fraction of excitatory , and oscillation rate of whole network slowed down . our main results are that in the network with @xmath39=2 msec , the spike death events were observed . and spike death can explain the mechanism of preventing the raise of the degree of synchrony . related synchrony properties was found also in weakly coupled hh neurons @xcite . in the case of weak coupling , the phase of neuron was perturbed by couplings , but the oscillation of neuron was not destroyed . in contrast , the new dynamical mechanism we suggested is proper for strong coupling and underlies the synchrony of the interrupted oscillations . additionally , it is notable that , for the existence of spike death , the firing synchronization of neuronal networks is different from the usual oscillator synchronization in which each oscillator is stable to perturbations @xcite . our work relates to that of drover et al . @xcite . in the elegant work , using a simplified neuron model of two variables , they proposed a mechanism for slowing firing down . they found that slow decay synaptic variable induces that trajectory is attracted toward the unstable fixed point of the simplified model . this is similar with the transient behavior of hh neuron we proposed here . however , with their mechanism , synaptic excitation is strongly synchronizing in networks in contrast with that spike death prevents synchrony . in neural networks , collective behaviors sensitively depend on the intrinsic dynamics of neurons @xcite , and many types of response of neuron to synaptic coupling may exist @xcite . it is interesting to make further studies on the relations among different responses and their effects on the collective behaviors of networks . the variability of strength of synapses was not involved in the present investigation . the synaptic strength is often affected by the activity of neurons through synaptic plasticity @xcite and synaptic adaption @xcite . the effect of changes of synaptic strength will be studied elsewhere . as mentioned above , the phenomena that low efficacy of synapses favors the generation of neuronal synchronization underlying seizure were obtained in experiments and numerical simulations @xcite . here we equated the strength and especially the characteristic time of synapses to the synaptic efficacy , and studied the mechanism by which synapses influence firing synchronization . the mechanism of the influence may have potentially values for understanding the way the realistic neural system works . we thank l. wang for many helpful discussions and c .- f . feng for comment reading . this work was supported by the nsf of china , grant no . j.x acknowledges financial support from fct ( portugal ) , grant no . sfrh / bpd/30425/2006 . d. h. zanette and a. s. mikhailov , phys . e * 58 * , 872 ( 1998 ) ; q. li , y. chen , and y. h. wang , phys . e * 65 * , 041916 ( 2002 ) ; s. j. wang , x. j. xu , z. x. wu , and y. h. wang , phys . e * 74 * , 041915 ( 2006 ) .
Abstract:

Answer: we investigated the influence of efficacy of synaptic interaction on firing synchronization in excitatory neuronal networks . we found spike death phenomena , namely , the state of neurons transits from limit cycle to fixed point or transient state . the phenomena occur under the perturbation of excitatory synaptic interaction that has a high efficacy . we showed that the decrease of synaptic current results in spike death through depressing the feedback of sodium ionic current . in the networks with spike death property the degree of synchronization is lower and unsensitive to the heterogeneity of neurons . the mechanism of the influence is that the transition of neuron state disrupts the adjustment of the rhythm of neuron oscillation and prevents further increase of firing synchronization .


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: wasp-12b is among the largest transiting planets discovered so far . first identified in an optical photometric transit survey ( * ? ? ? * hereafter , h09 ) , its mass is @xmath3 and radius @xmath4 , where @xmath5 and @xmath6 are the mass and radius of jupiter , respectively . after its discovery , additional observations have been acquired ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? interestingly , transit observations in the near - uv revealed a longer transit duration than in the optical @xcite . while the time of the egress occurs almost simultaneously for both set of observations , the ingress of the transit is first seen in the near - uv wavelength range . this asymmetric behavior has been explained by the presence of asymmetries in the planetary atmosphere . close - in giant gas planets are rather inflated and most have developed an exosphere ( hd  209458b : @xcite ; hd  189733b : @xcite ; corot - exo-1b : @xcite ; wasp-12b : @xcite ) that can fill or even overflow the planet s roche lobe @xcite . this may result in mass transfer through a lagrangian point to the star that could cause an asymmetry in the appearance of the transiting planet - star system as seen from the earth ( * ? ? ? * hereafter , lhv10 ) . asymmetries could also be produced by cometary tails . however , @xcite demonstrated for hd  209458b that a radiation - driven cometary tail would produce a late egress of the planetary transit light curve , instead of an early ingress . investigated the formation of a bow shock around the planet due to the interaction of the planet s magnetosphere with a stellar wind as the cause of the early ingress . assuming a typical solar wind mass - loss rate and adopting solar wind properties , derived the wind velocity at the distance of wasp-12b using a thermally driven wind model @xcite . they found that the wind is still subsonic at wasp-12b orbital distance , concluding that no bow shock is expected to form in the planet - wind interaction zone . here , we demonstrate that a bow shock can actually be formed around the planet if the relative azimuthal velocity between the planetary orbital motion and the ambient medium is taken into account . wasp-12b orbits its host star ( a late - f main - sequence star , with mass @xmath7 and radius @xmath8 ) at an orbital radius of @xmath9 , with an orbital period of @xmath10  d . due to its close proximity to the star , the flux of coronal particles impacting on the planet comes mainly from the azimuthal direction , as the planet moves at a keplerian orbital velocity of @xmath11  km  s@xmath12 around the star . therefore , stellar coronal material is compressed ahead of the planetary orbital motion , possibly forming a bow shock ahead of the planet . if such compressed material is optically thin , the planetary transit light curve is symmetrical with respect to the ingress and the egress of the planet . indeed , this is the case when the transit is observed at optical wavelengths . however , if the shocked material ahead of the planet can absorb enough stellar radiation , the observer will note an early ingress of the planet in the stellar disk , but no difference will be seen at the time of egress , as the shocked material is present only ahead of the planetary motion .    here , we investigate under which conditions the interaction of a planet with the stellar coronal plasma could lead to the formation of a bow shock ahead of the planetary orbital motion , and therefore , explain the early ingress observable in the near - uv . our shock model is described in section  [ sec.shockbasics ] . although we know the orbital radius of the planet , we do not know if at this radius the stellar magnetic field is still capable of confining the hot gas of the stellar corona , or if this gas is escaping in a wind . therefore , we investigate the validity of our model with respect to different stellar coronal conditions in sections  [ sec.corona ] and [ sec.wind ] . conclusions are presented in section  [ sec.conclusions ] . a bow shock around a planet is formed when the relative motion between the planet and the stellar corona / wind is supersonic . the shock configuration depends on the direction of the flux of particles that arrives at the planet . we illustrate two different limits of the shock configuration and an intermediate case in figure  [ fig.shockconditions ] , where @xmath13 is the deflection angle between the azimuthal direction of the planetary motion and @xmath14 , a vector that defines the outward direction of the shock . as seen from the planet , @xmath15 is the velocity of the impacting material .        the first shock limit , a `` dayside - shock '' , occurs when the dominant flux of particles impacting on the planet arises from the ( radial ) wind of its host star . for instance , the impact of the supersonic solar wind forms a bow shock at the dayside of earth s magnetosphere ( i.e. , at the side that faces the sun ) . this condition is illustrated in figure  [ fig.shockconditions](a ) and is met when @xmath16 , where @xmath17 and @xmath18 are the local radial stellar wind velocity and sound speed , respectively .    a second shock limit , an `` ahead - shock '' , occurs when the dominant flux of particles impacting on the planet arises from the relative azimuthal velocity between the planetary orbital motion and the ambient plasma . this condition is especially important when the planet orbits at a close distance to the star , and therefore , possesses a high @xmath19 . in this case , the velocity of the particles that the planet ` sees ' is supersonic if @xmath20 , where @xmath21 is the azimuthal velocity of the stellar corona . this condition is illustrated in figure  [ fig.shockconditions](b ) .    for intermediate cases , both the wind and the azimuthal relative velocities will contribute to the formation of a shock around the planet ( figure  [ fig.shockconditions](c ) ) and the deflection angle @xmath13 is given by @xmath22 in general , for a planet orbiting its host star at a close distance , the stellar wind is still accelerating and subsonic ( @xmath23 ) . in this case , conditions for an ahead - shock will more probably be met .    here , we adopt the estimates of to derive the minimum density of the stellar plasma at the orbital radius of wasp-12b . from the observations presented in @xcite , estimated the column density of the absorbing gas around the planet to be @xmath24 . for that , they assumed @xmath25 in the absorption lines of mg . because the time resolution of the transit measurement is sparse , it is difficult to place a firm constrain on the thickness of the absorbing material , which could either extend all the way to the planetary surface or be limited to a geometrically thin shocked region . through time differences of the optical / near - uv ingresses , estimated that the stand - off distance between the shock and the center of the planet is about @xmath26 , leading to a mg density @xmath27 . here , we convert this minimum density to hydrogen number density by using the observed metallicity of wasp-12 , @xmath28=0.3 $ ] @xmath29 where @xmath30 and @xmath31 are the solar element abundances for mg and h , respectively @xcite . therefore , the minimum density of the shocked stellar coronal plasma is @xmath32 the density behind such a shock would , in the adiabatic limit , be at most four times the density ahead of the shock . from equation  ( [ eq.shock-density ] ) , this would then require a pre - shock coronal density of @xmath33 in the external ambient medium at the orbital radius of the planet . the temperature of the absorbing material must be such as to allow for the presence of mg ii lines . for an adiabatic shock , the temperature immediately behind the shock depends on the squared mach number of the flow impacting on the planet . to determine the temperature and ionization structures of the atmosphere of the planet , and therefore quantitatively evaluate the optical depth of mg ii lines , detailed radiative transfer models are needed . in this section , we assume that the orbit of the planet lies within the stellar magnetosphere . in this case , the planet moves through a medium that is confined by a rigid stellar magnetic field and so rotates at the stellar rotation rate . the relative azimuthal velocity between the planet and this medium is therefore @xmath34 where @xmath35 is the velocity of the medium corotating with the star , and @xmath36 is the stellar period of rotation . if @xmath37 , a bow shock forms around the planet . we may therefore define a critical stellar rotation period such that @xmath38 and @xmath39 presents @xmath40 as a function of the coronal temperature @xmath41 . the shaded areas show regions on the @xmath42 parameter space where no shock is formed . in our notation , @xmath43 refers to cases where the star is counterrotating with respect to the planetary orbital rotation ( e.g. , wasp-8b ; * ? ? ? * ) . the vertical line in figure  [ min_prot_t ] represents a critical coronal temperature where @xmath44 and @xmath45 ( cf . equation  ( [ eq.pcrit ] ) ) . for wasp-12 , we obtain @xmath46 where @xmath47 is the mean particle mass , @xmath48 is the boltzmann constant , and @xmath49 is the proton mass . we adopt @xmath50 . above this temperature , no shock formation is possible at any positive @xmath36 . for all the remaining areas of the plot an ahead - shock will develop . , or @xmath51  d for wasp-12b . in such situation , the planet orbits beyond the keplerian corotation radius . the coronal plasma lags behind the planetary motion and the shock trails the planet . observationally , we would detect a late egress , instead of an early ingress . however , figure  [ min_prot_t ] shows that for @xmath52  d ( wasp-12b ) , the relative azimuthal velocity @xmath53 is subsonic , therefore , no shock will be formed . ] although @xmath36 is unknown , wasp-12 is likely to be a slow rotator @xcite . with this information , the presence of an ahead - shock for wasp-12b constrains the coronal temperature to be @xmath54  k.    ) ) required for the formation of a bow shock as a function of the environment temperature . negative periods refer to cases where the planet is in a retrograde orbit.,width=317 ]    for an isothermal corona , the density is @xmath55 , \end{aligned}\ ] ] where @xmath56 is the density at the coronal base . from equation  ( [ eq.dens-hyd ] ) , we note that the more slowly the star rotates , the lower is the density at any given radius . for the very slow stellar rotation rates expected for stars with detected transiting planets , @xmath36 has little influence on the density . the solid lines in figure  [ fig.wind-density ] show the density variation for a range of orbital radius for three different temperatures and a solar value @xmath57  d. the density was scaled to match the observationally derived value of @xmath58 ( equation  ( [ eq.coronal-density ] ) ) , which is the minimum density required for the detection of the early ingress . except for the very low temperature case ( black solid line ) , a typical solar coronal base density ( i.e. , @xmath59 ; * * ) provides the condition for the existence of a shock . , width=317 ]    we note that the observationally derived stand - off distance from the shock to the center of the planet can be taken as approximately the extent of the planetary magnetosphere @xmath60 . pressure balance between the coronal total pressure and the planet total pressure requires that , at @xmath61 , @xmath62 where @xmath63 , @xmath64 and @xmath65 are the local coronal mass density , thermal pressure , and magnetic field intensity , and @xmath66 and @xmath67 are the planet thermal pressure and magnetic field intensity at @xmath61 . neglecting the kinetic term and the thermal pressures  g ) , the magnetic pressure is @xmath68  dyn  @xmath69 . the kinetic term @xmath70 in equation  ( [ eq.equilibrium ] ) is dependent on the period of the star , which for wasp-12 is unknown . for example , for @xmath57  d , @xmath71  dyn  @xmath69 . for @xmath72 k , @xmath73  dyn  @xmath69 . to estimate the thermal pressure of the planet , we use values derived from @xcite for planetary density ( @xmath74  @xmath75 ) and temperature ( @xmath76  k ) , which result in @xmath77  dyn  @xmath69 . because @xmath78 is much larger than the aforementioned pressures , the magnetic terms in equation  ( [ eq.equilibrium ] ) will dominate . ] , we may therefore rewrite previous equation as @xmath79 from which we can estimate an upper limit for the planetary magnetic field intensity .    for a stellar dipolar magnetic field , @xmath80 . assuming that the planetary magnetic field is also dipolar , then @xmath81 . @xmath82 and @xmath1 are the magnetic field intensities at the stellar and planetary surfaces , respectively . from equation  ( [ eq.equilibrium1 ] ) , we have @xmath83 where we have used the observationally derived characteristics for the wasp-12 system : @xmath84 and @xmath85 . so far , tentative measurements of the stellar magnetic field have not provided significant detections ( for theoretical predictions , see * ? ? ? * ) , but have yielded an upper limit of @xmath86  g in the longitudinal component @xcite . adopting this upper limit as @xmath82 , our model predicts a maximum planetary magnetic field of about @xmath87  g. in this section , we investigate the scenario where the stellar magnetic field is not strong enough to confine the coronal plasma , which expands in the form of stellar wind . for a hot corona , the thermal pressure gradient is able to drive a wind @xcite . we consider the simplest case of a thermally - driven wind which is not magnetically channeled . the wind radial velocity @xmath17 is derived from the integration of the differential equation @xmath88 figure  [ fig.wind-velocity ] presents the wind velocity profile at different orbital radius up to a few stellar radii for three isothermal wind temperatures . in the distance range @xmath89 , the wind is still accelerating , and the radial velocities are considerably smaller than the terminal velocities @xmath90 achieved : @xmath91 , @xmath92 and @xmath93  km  s@xmath12 for @xmath94 , @xmath95 , and @xmath96k , respectively . the filled circles are the sound speed of the wind : @xmath97 and @xmath98  km  s@xmath12 for @xmath99 and @xmath96k . dashed lines represent the regions of supersonic velocities . for the lowest - temperature case ( @xmath100k ) , the wind becomes supersonic at a much larger radii , not shown in figure  [ fig.wind-velocity ] . the vertical dotted line shown in figure  [ fig.wind-velocity ] represents the orbital radius of wasp-12b .    ) and dashed lines represent the regions of supersonic velocities , which , for @xmath101k , is outside the range of orbital radii shown above . [ fig.wind-velocity ] , width=317 ]    the density structure of the wind is derived from the conservation of mass , where @xmath102 is a constant of the steady - state wind . the dashed lines in figure  [ fig.wind-density ] present the density profile for an isothermal wind , for three coronal temperatures . in contrast to the hydrostatic corona ( solid lines ) , @xmath103 does not depend on @xmath36 . except for the lowest temperature adopted , where the density is well described by a hydrostatic density profile , the coronal base density @xmath56 required for the formation of a shock in the wind case is larger than the one required by the hydrostatic corona case . nevertheless , as in the static scenario , a typical solar coronal density still allows for the formation of a shock .    in the stellar wind case , there are no azimuthal forces acting on the flow . therefore , through conservation of angular momentum of the particles leaving the star ( @xmath104 ) , the azimuthal velocity of the wind is @xmath105 . this implies that @xmath53 for the wind case is @xmath106 which is independent of the wind temperature or radial velocity . for slow rotators , @xmath53 becomes almost independent of @xmath36 and is given by @xmath107 . the condition for shock formation requires that @xmath108 , implying that @xmath109 ) . this means that the critical temperature derived in the previous section ( equation  ( [ eq.tcrit ] ) ) is enhanced by a temperature @xmath110  k , implying that for @xmath111 , a shock will be formed .    in the wind case , figure  [ fig.wind-velocity ] shows that @xmath112 at the orbital radius of wasp-12b , implying that the shock will not form directly ahead of the planet , but at an intermediate angle @xmath13 . because @xmath13 depends on @xmath17 ( equation  ( [ eq.theta ] ) ) , which depends on the temperature , the orientation of the shock is temperature - dependent . figure  [ fig.theta ] shows how the angle @xmath13 depends on the wind temperature for @xmath57  d. while in the hydrostatic case , we found that an ahead - shock ( @xmath113 ) is always formed , the wind case requires a very low temperature to form an ahead - shock . for higher wind temperatures , the shock normal gets closer to the line of centers between the planet and the star and @xmath13 approaches @xmath114 .     makes to the relative azimuthal velocity of the planet . [ fig.theta],width=317 ] motivated by the recent observation on the light curve asymmetry of the transit hot - jupiter wasp-12b @xcite , we proposed a model where the interaction of the stellar plasma with the planet results in the formation of a bow shock around the planet , which could explain the early transit ingress observed in the near - uv . although we know the orbital radius of the planet , we do not know if at this radius the stellar magnetic field is still capable of confining the hot gas of its corona , or if this gas is escaping in a wind . therefore , we investigated the physical conditions of the external ambient medium around the planet that could allow for the formation of such a shock . in either case , for plasma temperatures @xmath115k we expect that a shock capable of compressing the plasma to the observed densities will form around the planet . in the case where the coronal plasma is confined , we showed that the geometry of this shock is independent of the plasma temperature . the shock forms ahead of the planet in its orbital path . in the unconfined ( wind ) case , the angle between the shock normal and the direction of planetary motion depends on the ratio of the radial wind speed and the azimuthal speed of the planet relative to the stellar wind . since the wind speed depends on the temperature , the orientation of the shock is temperature - dependent . if the planet s orbit takes it through both regions of confined plasma and also wind plasma , the orientation and density of the shock may change , producing time - dependent absorption and duration of transit . given the large range of stellar rotation rates and coronal / wind temperatures at which a shock is capable of producing the densities estimated by , we conclude that this is likely to be a common feature of transiting systems and may prove to be a useful tool in setting limits on planetary magnetic field strengths . for the case of wasp-12b , we derived an upper limit for the planet s magnetic field of about @xmath87  g.    22 natexlab#1#1 , p. , baglin , a. , & auvergne , m. et  al . 2008 , , 482 , l17 , c.  j. , harrington , j. , & hardy , r.  a. et al . 2010 , arxiv:1003.2763 christensen u.  r. , holzwarth v. , reiners a. , 2009 , nature , 457 , 167 , d. , lecavelier des etangs , a. , hbrard , g. et  al . 2008 , , 483 , 933 , l. , haswell , c.  a. , & froning , c.  s. et al . 2010a , , 714 , l222 , l. , bagnulo , s. , & elmasli , a. et al .  2010b , , 720 , 872 , n. , asplund , m. , & sauval , a.  j. 2007 , space science reviews , 130 , 105 , p. , lin , d.  n.  c. , & bodenheimer , p.  h. 2003 , , 588 , 509 , l. , collier - cameron , a. , & loeillet , b. et  al . 2009 , , 693 , 1920 , h09 , n. , pont , f. , hebrard , g. , et  al . 2010 , arxiv:1004.1809 , l. , burrows , a. , & spiegel , d.  s. 2010 , , 713 , 751 lai d. , helling c. , van den heuvel e.  p.  j. , 2010 , apj , 721 , 923 , lhv10 , a. , ehrenreich , d. , vidal - madjar , a. , et  al . 2010 , , 514 , a72 , s. , miller , n. , lin , d.  n.  c. , & fortney , j.  j. 2010 , , 463 , 1054 , j.  l. , yang , h. , france , k. , et  al . 2010 , , 717 , 1291 , m. , coughlin , j.  l. , sing , d.  k. , et  al . 2010 , , 716 , l36 , r.  a. , chiang , e.  i. , & murray , n. 2009 , , 693 , 23 , e.  n. 1958 , , 128 , 664 , d. , anderson , d. , & collier cameron , a. et al . 2010 , , 517 , l1 vidal - madjar , a. , lecavelier des etangs , a. , & dsert , j .- m . et al . 2003 , , 422 , 143 vidal - madjar , a. , lecavelier des etangs , a. , & dsert , j .- m . et al . 2008 , , 676 , l57 , g.  l. 1988 , , 325 , 442
Abstract:

Answer: recently , @xcite observed that the uv transit of wasp-12b showed an early ingress compared to the optical transit . we suggest that the resulting early ingress is caused by a bow shock ahead of the planetary orbital motion . in this letter we investigate the conditions that might lead to the formation of such a bow shock . we consider two scenarios : ( 1 ) the stellar magnetic field is strong enough to confine the hot coronal plasma out to the planetary orbit and ( 2 ) the stellar magnetic field is unable to confine the plasma , which escapes in a wind . in both cases , a shock capable of compressing plasma to the observed densities will form around the planet for plasma temperatures @xmath0k . in the confined case , the shock always forms directly ahead of the planet , but in the wind case the shock orientation depends on the wind speed and hence on the plasma temperature . for higher wind temperatures , the shock forms closer to the line of centers between the planet and the star . we conclude that shock formation leading to an observable early uv ingress is likely to be a common feature of transiting systems and may prove to be a useful tool in setting limits on planetary magnetic field strengths @xmath1 . in the case of wasp-12b , we derive an upper limit of about @xmath2  g.


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 concept of self - organized criticality ( soc ) @xcite presents a simple way of modelling slowly - driven out of equilibrium systems . the interesting natural systems thought to exhibit soc , such as rain - fall @xcite , earthquakes @xcite , economic markets @xcite , biological evolution @xcite , or the brain @xcite , are difficult to study in a controlled experiment , and more simple systems have to be found with which the predictions of soc can be tested quantitatively . among the first toy systems studied , were sand - piles , however it has been shown that due to the appearance of e.g. kinetic effects , real sand does not behave as an soc system @xcite . another system , that has shown indications of soc in a controlled environment is the dynamics of vortices in type - ii superconductors . here , power - law scaling has been observed in the past @xcite , however not all authors have found avalanche scaling @xcite . this is probably due to the fact that internal avalanches are not measured in the usual setup @xcite . recently behnia _ et al . _ @xcite have measured internal avalanches using arrays of hall - probes , where they do indeed find power - law behavior of the avalanches . the ultimate hall - mark of soc , however , is the observation of finite size - scaling , indicating a true critical dynamics . so far , there are only two experimental systems , a one dimensional ( 1d ) pile of rice @xcite or a 1d pile of steel balls with a random arrangement at its bottom layer @xcite , which have been shown to exhibit finite size scaling in their avalanche distributions . the critical exponents obtained from a finite size scaling analysis , would yield more information on a soc system , as was discussed by paczuski _ @xcite for models in a stationary critical state . a system evolving through avalanches , which are distributed according to a power - law , will also show roughening dynamics of its surface . this can be seen as an indication to the origin of the abundance of fractal or self - affine structures in nature . roughening dynamics has been studied extensively in the past , including the study of interfaces in porous media @xcite , the growth of bacterial colonies @xcite , the slow combustion @xcite or the rupture @xcite of paper . these systems are well characterized experimentally and in many cases , the interface dynamics can be modelled analytically by the kardar - parisi - zhang ( kpz ) equation @xcite . in one spatial dimension , the kpz equation is solved exactly for the case of white noise , such that good comparison between experiment and theory is possible . however , e.g. in the burning of paper , there are avalanches observed in the propagation of the front , which are reminiscent of soc dynamics @xcite . this has recently been addressed by alava _ @xcite , who proposed a mapping between soc models and the kpz equation similar to that of paczuski and boettcher @xcite . with this mapping , they have been able to obtain soc from the kpz equation , with the proper kind of noise @xcite .    moreover , roughening processes have also been found in the experimental soc systems discussed above . the profile of a 1d rice pile is a rough interface @xcite , as is the front of magnetic flux penetrating a thin film superconducting sample @xcite . however , the universal scaling relations derived by paczuski _ _ @xcite have not yet been tested experimentally . but the fact that both roughening phenomena and avalanche dynamics can be observed in the same experimental system , such as a 3d rice - pile , means that the universal connections between soc and surface roughening can be tested experimentally .    here we present a two - dimensional ( 2d ) system , the surface of a 3d pile of rice , showing both roughening behavior and avalanche dynamics . the avalanche dynamics is studied in terms of the avalanche size distribution for different sizes of the field of view , @xmath7 . this allows the observation of finite size scaling and the determination of the critical exponents @xmath1 , describing the dependence of the cutoff scale on @xmath7 , and @xmath0 , the exponent of the avalanche size distribution . the roughening behavior is studied via the power spectrum of the surface @xcite , in both space and time , resulting in a determination of the roughness and growth exponents respectively @xcite . the connection between these two phenomena is shown from the derivation of scaling relations between the different exponents , where we obtain excellent agreement using the experimentally determined values .    the experimental setup for the growth of the rice - pile , as well as the reconstruction method used to determine the structure of the pile and the size of the avalanches is discussed in section ii . the avalanche size distributions , and their finite size scaling is presented in section iiia , together with the determination of the critical exponents . in section iiib , the surface roughness is analyzed and the necessary techniques are briefly introduced . still in section iiib , the scaling relations between the roughness exponents and the critical exponents are introduced . these results are also put into a wider context and compared with results from kpz roughening systems @xcite . the experiments were carried out on long grained rice with dimensions of typically @xmath82x2x7 mm@xmath9 , similar to rice a of ref . the pile was grown from a uniform line source 1 m wide . this uniform distribution was achieved via a custom built mechanical distributor based on a nail - board producing a binary distribution @xcite . the actual setup consists of a board with an arrangement of triangles as shown in fig . this means that the possible pathways of the rice are continuously split at each level , such that at the bottom we end up with a row of 64 uniformly distributed compartments . the uniformity of the distribution was measured to be of the order of 5@xmath10 .    at the bottom of the distributor , the grains are slowed down by a sheet of plastic before they hit the top of the pile . the rice is fed to the top of the distribution board , from a point source , which drops rice at an average rate of @xmath85g / s . this corresponds to 1500 grains per image , which means that over the length of the line there are 1 to 2 grains dropped per place . this rate is uniform over the time scale between two images . once a rice pile is grown , we measure the further evolution of the surface coordinates using a specially developed real - time technique based on the projection of a set of lines . in order to increase the spatial resolution in the direction of growth , we use a set of lines in the base colors ( red , green and blue ) , which can be easily separated digitally for the analysis . the lines are projected at a right angle to the average surface of the pile and observed at 45@xmath11 to both the direction of projection as well as the surface . this is illustrated in fig . [ lines ] , where a raw image is shown together with the lines extracted from it . from the distortions of the lines observed in this way , the coordinates in 3d space of the pile surface can be calculated by simple geometry @xcite . with this reconstruction technique , the 3d coordinates of the whole field of view can be determined with an accuracy of 2 - 3 mm , as we have checked on a number of test surfaces . using both structured and smooth surfaces , the resolution and accuracy were determined independently . since the resolution is roughly matched with the size of the grains , the method is well suited for the present purpose .    in a single experimental run , the growth is studied for a period of @xmath84 hours , with a picture taken every 30 s. thus an experiment consists of @xmath8 480 images . the pictures are taken with a high resolution digital camera , having a resolution of 2048x1536 pixels . for every time step , the surface structure is reconstructed , which gives information about the roughening process , as well as the avalanche dynamics . the size and shape of the avalanches can be determined from the height difference of the surface between two consecutive images . the overall growth of the pile is subtracted from this difference , however this correction is negligible and does not influence the results . this is shown in fig.[diff ] , where the height difference @xmath12 is shown for a medium size avalanche . this allows the study of internal avalanches instead of just the off - edge ones , which has proven to be important in previous studies ( for instance on superconducting vortex avalanches @xcite and 1d piles of rice @xcite ) . in order to obtain the size of such an avalanche , @xmath13 is integrated over the area , which yields the displaced volume , @xmath14 , corresponding to the size of the avalanche in mm@xmath9 . in order to use natural units , we will in the following measure the size of the avalanche in terms of the number of toppled grains , which is @xmath15 , where @xmath16 35 mm@xmath9 is the volume of a single grain of rice . the data discussed in this article was obtained in three separate experiments , with a total of 1330 images . in this section the experimental results are presented . we first show that there is finite size scaling in the avalanche statistics , which indicates the appearance of a critical state in the system . second , we characterize the surface roughness of the rice pile in 2d , both in space and time . the two characteristic exponents , @xmath5 for the spatial behavior , and @xmath6 for the dynamics @xcite can also be obtained from scaling relations of the critical exponents of the avalanche statistics and are in excellent agreement with those determined from a standard power - spectrum analysis . the time dependence of the displaced volume , @xmath14 , in each time step is shown in fig . [ time ] , for all three experiments discussed here . in the figure , the data are shown for avalanche sizes integrated over the total field of view of 600x600 mm@xmath17 . a histogram of these data , giving the avalanche size distribution is shown by the diamonds in fig .  [ hist]a ) . in that figure , size distributions of subsets of the data corresponding to smaller fields of view ( l = 50 , 150 and 300 mm ) are also shown . due to experimental resolution , the smallest avalanches measured for each subset depend on the size of that subset . the range of sizes from 50x50 mm@xmath17 to 600x600 mm@xmath17 spans more than a decade and all of the data taken together show a power - law scaling of the size distribution over three decades , with an exponent of @xmath0 = 1.20(5 ) . these data can however be scaled to fall onto a single curve , as shown in fig  [ hist]b ) . here the avalanche sizes are scaled with @xmath18 and the probabilities are scaled with @xmath19 . the good curve collapse visible in fig  [ hist]b ) indicates the presence of finite size scaling in the data , one of the hall - marks of critical behavior in a system and thus evidence of the appearance of soc in the 3d rice pile .    due to the finite size scaling found from the curve collapse , the avalanche size distribution as a function of system size can be written as a function of one parameter only : @xmath20 where @xmath21 is constant up to a cutoff scale @xmath22 the exponents @xmath0 and @xmath1 used to obtain the curve collapse in fig .  [ hist]b ) were @xmath0 = 1.21(2 ) and @xmath1 = 1.99(2 ) . we note here that in usual finite size scaling , separate experiments are used instead of subsets . however , as we have tested on simulations of a 2d version of the oslo - model @xcite , finite size scaling using subsets yields the same exponents , but is much more easily implemented experimentally . as an additional test , the avalanche dimension can be determined directly using a box - counting method in 3d . this yields @xmath1 = 2.05(10 ) consistent with the result from finite size scaling . the roughness of a surface can be characterized in different ways . most commonly , this is done via the width of the interface as given by its standard deviation @xcite . for a 2d surface , this is given by @xmath23 where @xmath24 is the average profile of the surface height . when the roughness of the surface has reached saturation , its value @xmath25 , will scale with the system size as a power - law , @xmath26 , where @xmath5 is the roughness exponent . a similar type of analysis can be achieved via the power spectrum of the surface . here , one has the advantage that better statistics can be obtained , since the whole surface can be used at once @xcite . we determined the power spectrum from a radial average of the 2d fourier transform of the surface @xmath27 where @xmath28 . the fourier transforms where performed via an fft algorithm , which is why a data subset corresponding to a power of two ( 512x512 ) was used . from the power spectrum we then determine the distribution function @xmath29 @xmath30 from which it can be shown @xcite that @xmath31 , where @xmath32 is the length scale over which the width is calculated . the roughness exponent @xmath5 can therefore be reliably obtained from a determination of @xmath29 . an alternative way of obtaining the roughness exponent is via the two point correlation function@xcite . @xmath33 the radial average of which obeys the same scaling behavior as the distribution function @xcite .    in fig . [ alpha ] , the correlation and distribution functions for the rice - pile surface is shown , averaged over all time steps of all experiments ( 1330 images ) . the part of the surface studied consists of 512x512 mm in the center of the surface laterally ( along the horizontal direction ) and starting from a height of 200 mm . as can be seen from the double logarithmic plot , there is power - law scaling with an exponent , @xmath5 = 0.42(3 ) . the scaling behavior of the interface width can , however , also be determined from the cut - off size of the avalanches in a system with size @xmath7 . the size of such an avalanche will be proportional to the maximal height difference , given by the saturation width @xmath25 , times the maximal area of the avalanche . as can be seen from fig . [ diff ] , the avalanche - shape is a fractal , which means that @xmath34 , where @xmath2 is the fractal dimension of the cluster of displaced grains . equating this expression with eq . [ cutoff ] , one obtains @xmath35 , which implies the scaling relation @xmath36 such a scaling relation between the roughness exponent and the exponents characterizing the avalanches was also derived by paczuski _ @xcite from different arguments .    in order to determine @xmath2 in an accurate manner , we studied avalanches more than two standard deviations bigger than the average size . this corresponds to a subset of @xmath8 100 images from all experiments . after thresholding the height difference fields @xmath12 at their mean value , we applied a simple box counting method @xcite to the resulting clusters . this is shown in fig . [ frac ] , where the number active pixels in a box of given size is shown as a function of the box size . the result is a power law with an exponent of -1.58(2 ) , indicating a fractal dimension of @xmath2 = 1.58(2 ) . inserting this value into scaling relation ( [ alpha2 ] ) , the roughness exponent is found to be @xmath5 = 0.41(3 ) , in excellent agreement with that determined from the roughness analysis . the dynamics of the roughening process is likewise analyzed via the distribution function . in contrast to the roughness analysis above , we now determine the distribution function for the time dependence @xmath37 for each pixel , where @xmath38 denotes the average over the duration of the experiment . again , the distribution function , given by @xmath39 is equal to the momentary width , @xmath40 . thus the growth exponent @xmath41 , describing the scaling of the width with time , can be determined from the distribution function . again , the correlation function shows the same scaling behavior in time , allowing a separate determination of @xmath41 . both results are shown in fig . [ beta ] , where it can be seen that there is good power - law scaling of the distribution function , as well as the correlation function , over two decades with an exponent of @xmath41 = 0.28(3 ) . another way of describing the dynamics of an interface is via the dynamic exponent @xmath6 , which describes the scaling of the saturation time with the system size @xmath42 it can be easily obtained @xcite from the roughness and growth exponents via @xmath43 from the values determined above we thus obtain @xmath6 = 1.5(1 ) .    the dynamic exponent , @xmath6 , can also be derived from a scaling relation using the critical exponents of the avalanche dynamics . the saturation time , @xmath44 , will be roughly given by the time it takes for an avalanche of the cut - off size , @xmath45 , to appear . since material is added to the system at a constant rate , the number of grains added until a cut - off avalanche occurs is proportional to the cross - over time @xmath44 . on the other hand , the material added not only result in an avalanche of size @xmath45 , but will also be lost by smaller avalanches , such that the total material necessary to obtain an avalanche of size @xmath45 can be estimated from : @xmath46 from eq . [ cutoff ] , we obtain for the cross - over time @xmath47 . with eq . [ tx ] , this immediately leads to the scaling relation @xmath48 which was also derived by paczuski _ @xcite from more general arguments . using the values for @xmath1 and @xmath0 determined above from the curve collapse of the avalanche size distributions , we obtain @xmath6 = 1.56(8 ) , which is again in very good agreement with the result of the power spectrum analysis , @xmath6 = 1.5(1 ) . note also that the exponents @xmath5 and @xmath6 fulfill the kpz scaling relation@xcite @xmath49 = 2 , which is valid independent of the dimension of the system and depends only on the fact that the growth is driven by height gradients . furthermore , the roughness exponent @xmath5 = 0.42(3 ) is in good agreement with a numerical determination of the behavior of the 2d kpz equation @xcite . the connection between roughening and soc will be discussed in more detail below . in conclusion , we have presented results on both the avalanche and the roughening behavior of an experimental soc system . in addition , we presented simple arguments for universal scaling relations , derived by paczuski _ @xcite on general grounds , connecting the avalanche and roughening behavior . we obtain _ quantitative _ agreement of the experimental exponents characterizing the roughness and those describing the avalanche statistics . this means that the roughening of the surface of the pile is governed by its underlying avalanche dynamics , which was already conjectured from models of interface depinning . earlier experimental studies @xcite have considered the avalanche dynamics and the roughening behavior of the 1d version of the present system , where finite size scaling was also found @xcite . the above scaling relations , showing quantitatively the connection between the two phenomena were not previously tested . the fact that kpz dynamics is observed in a soc system , as indicated by the fact that the roughness and dynamic exponents obtained here fulfill the kpz scaling relation , is surprising . this is mostly because in soc systems , the critical state is self organized , whereas in kpz systems , the roughening is put into the dynamics by obeying the proper symmetries . however , it has been argued by alava _ @xcite , that there is a mapping of soc models , which naturally lead to the necessary symmetries . in this context , it should also be noted that experiments on kpz systems such as burning paper have also observed avalanches in the advance of fronts @xcite . this work was supported by fom ( stichting voor fundamenteel onderzoek der materie ) , which is financially supported by nwo ( nederlandse organisatie voor wetenschappelijk onderzoek ) . 10 p.  bak , c.  tang , and k.  wiesenfeld , phys . lett . * 59 * , 381 ( 1987 ) and phys . a * 38 * , 364 ( 1988 ) . o.  peters , c.  hertlein , and k.  christensen , phys . lett . * 88 * , 018701 ( 2002 ) ; o. peters , and k. christensen , preprint . see e.g. p.  bak , k.  christensen , l.  danon , and t.  scanlon , phys . * 88 * , 178501 ( 2002 ) ; a.  sornette , and d.  sornette , europhys . lett . * 9 * , 197 ( 1989 ) ; j.  m.  carlson , j.  s.  langer , and b.  e.  shaw , rev . phys . * 66 * , 657 ( 1994 ) . see e.g. j. a. scheinmann , and m. woodford , am . * 84 * , 417 ( 1994 ) ; p. bak , k. chen , j. a. scheinmann , and m. woodford , ricerche economiche * 47 * , 3 ( 1993 ) . see e.g. j. bascompte , and r. v. sol , trends ecol . * 10 * , 361 ( 1995 ) ; h. flyvbjerg , p. bak , and k. sneppen , phys . lett . * 71 * , 4087 ( 1993 ) ; s. a. kaufmann , and s. j. johnsen , j. theor . biol . * 149 * , 467 ( 1991 ) . d. stassinopoulos , and p. bak , phys . e * 51 * , 5033 ( 1995 ) . h.  m.  jaeger , s.  r.  nagel , and r.  p.  behringer , rev . * 68 * , 1259 ( 1996 ) ; s.  r.  nagel , _ ibid . _ * 64 * , 321 ( 1992 ) and references therein . s.  field , j.  witt , and f.  nori , phys . lett . * 74 * , 1206 ( 1995 ) ; c.  m.  aegerter , phys . e * 58 * , 1438 ( 1998 ) . see e.g. r.  j.  zieve , t.  f.  rosenbaum , h.  m.  jaeger , g.  t.  seidler , g.  w.  crabtree , and u.  welp , phys . b * 53 * , 11849 ( 1996 ) ; e.  r.  nowak , o.  w.  taylor , l.  liu , h.  m.  jaeger , and t.  i.  selinder , phys . rev . b * 55 * , 11702 ( 1997 ) . j. feder , _ fractals _ ( plenum , new york , 1989 ) . k.  behnia , c.  capan , d.  mailly , and b.  etienne , phys . b * 61 * , r3815 ( 2000 ) ; k.  behnia , c.  capan , d.  mailly , and b.  etienne , j. mag . mag . mat . * 226 * , 370 ( 2001 ) . v.  frette , k.  christensen , a.  malthe - s@xmath50renssen , j.  feder , t.  j@xmath50ssang , and p.  meakin , nature ( london ) * 379 * , 49 ( 1996 ) . e. altshuler , o.  ramos , c.  martinez , l.  e.  flores , and c.  noda , phys . lett . * 86 * , 5490 ( 2001 ) . m. paczuski , s. maslov , and p. bak , phys . e * 53 * , 414 ( 1996 ) . m.  a.  rubio , c.  a.  edwards , a.  dougherty , and j.  p.  gollub , phys . lett . * 63 * , 1685 ( 1989 ) ; v.  k.  horvth , f.  family , and t.  vicsek , j. phys . a * 24 * , l25 ( 1991 ) ; s. he , g. l. m. k. s. kahanda , and p .- z . wong , phys . * 69 * , 3731 ( 1992 ) . e. ben - jacob , o.  shochet , a.  tenenbaum , i.  cohen , a.  czirk , and t.vicsek , nature(london ) * 368 * , 46 ( 1994 ) ; s.  matsuura , and s.  miyazima , fractals * 1 * , 336 ( 1993 ) ; t.  vicsek , m.  cserz , and v.  k.  horvth , physica a * 167 * , 315 ( 1990 ) . j.  maunuksela , m.  myllys , o.  p.  khknen , j.  timonen , n.  provatas , m.  j.  alava , and t.  ala - nissila , phys . lett . * 79 * , 1515 ( 1997 ) . j.  kertsz , v.  k.  horvth , and f.  weber , fractals * 1 * , 67 ( 1993 ) . m.  kardar , g.  parisi , and y .- c . zhang , phys . lett . * 56 * , 889 ( 1986 ) . a. myllys , j.  maunuksela , a.  alava , t.  ala - nissila , j.  merikoski , and j.  timonen , phys . e * 64 * , 036101 ( 2001 ) ; and j.  timonen , private communication . m. alava , and k. b. lauritsen , europhys . lett . * 53 * , 569 ( 2001 ) ; g. szabo , m. alava , and j. kertesz , europhys . lett . * 57 * , 665 ( 2002 ) ; m. alava , j. phys . c * 14(9 ) * , 2353 ( 2002 ) . m. paczuski , and s. boettcher , phys . lett . * 77 * , 111 ( 1996 ) . a.  malthe - s@xmath50renssen , j.  feder , k.  christensen , v.  frette , and t.  j@xmath50ssang , phys . lett . * 83 * , 764 ( 1999 ) . r.  surdeanu , r.  j.  wijngaarden , e.  visser , j.  m.  huijbregtse , j.  h.  rector , b.  dam , and r.  griessen , phys . lett . * 83 * , 2054 ( 1999 ) . j.  schmittbuhl , j .- vilotte , and s.  roux , phys . e * 51 * , 131 ( 1995 ) ; m.  siegert _ ibid . _ * 53 * , 3209 ( 1996 ) ; j.  m.  lopez , m.  a.  rodriguez , and r.  cuerno , _ ibid . _ * 56 * , 3993 ( 1997 ) . a.  l.  barabasi and h.  e.  stanley _ fractal concepts in surface growth _ ( cambridge university press , 1995 ) . j. g. amar , and f. family , phys . a * 41 * , 3399 ( 1990 ) . r. gnther , masters thesis , vrije universiteit ( 2002 ) . c.  m.  aegerter , physica a * 319 * , 1 ( 2003 ) . b. b. mandelbrot , _ the fractal geometry of nature _ ( freeman , new york , 1983 ) .
Abstract:

Answer: we present a two - dimensional system which exhibits features of self - organized criticality . the avalanches which occur on the surface of a pile of rice are found to exhibit finite size scaling in their probability distribution . the critical exponents are @xmath0 = 1.21(2 ) for the avalanche size distribution and @xmath1 = 1.99(2 ) for the cut - off size . furthermore the geometry of the avalanches is studied leading to a fractal dimension of the active sites of @xmath2 = 1.58(2 ) . using a set of scaling relations , we can calculate the roughness exponent @xmath3 = 0.41(3 ) and the dynamic exponent @xmath4 = 1.56(8 ) . this result is compared with that obtained from a power spectrum analysis of the surface roughness , which yields @xmath5 = 0.42(3 ) and @xmath6 = 1.5(1 ) in excellent agreement with those obtained from the scaling relations . -1.5 cm 24.5 cm    2