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: single - ended optical fiber link characterization techniques have been a point of major focus since the late seventies when barnosky and others demonstrated that , by measuring the rayleigh backscattered signal from a pulse sent into the optical fiber , one could extract its physical properties such as attenuation coefficient and eventual power losses along its length @xcite . presently , different techniques for optical fiber characterization are available and standardized such as the optical time domain reflectometry ( otdr ) and frequency domain reflectometry ( ofdr ) @xcite . in any case , embedding any of these fiber - monitoring techniques in passive optical networks ( pon ) systems remains a challenge in terms of technology and cost - effectiveness not to mention the in - service monitoring paradigm @xcite . otdr makes use of high - power short optical pulses while measuring the rayleigh backscattered power as a function of time , its spatial resolution being defined by the duration of the optical pulse . ofdr may present two different flavors : coherent - ofdr ( c - ofdr ) , where the wavelength of a low - power continuous - wave coherent light source is swept and the heterodyne beat between the backscattered light signal and a reference reflection is measured ; and incoherent - ofdr ( i - ofdr ) , where a radio - frequency optical sub - carrier frequency is swept and , once again , a heterodyne beat between the detected backscattered signal and the modulating signal is measured . in both c- and i - ofdr , the spatial profile of the fiber is obtained from the fourier transform of the detected ofdr signal and the frequency range covered by the optical sweep limits its resolution . a very complete review of frequency domain optical network analysis , where the distinctions between c- and i - ofdr , their pro s , and con s are highlighted , may be found in @xcite .    in terms of the apparatus necessary to guarantee reasonable high dynamic range and spatial resolution , otdr and ofdr differ greatly : on one hand , otdr must employ high bandwidth detection in order to increase spatial resolution since it is directly proportional to the probe pulse width , which , in general , has a trade off associated to a reduction of dynamic range since the energy is reduced along with pulse duration @xcite . frequency domain techniques usually make use of a frequency sweep at the source level and a fourier transform which provides the conversion from frequency to distance based on the sweeping rate . c - ofdr does not require high frequency response detectors since the detection is performed in the base band , but highly coherent sources are needed in order to achieve long range in optical fibers @xcite . on the other hand , the resolution of the frequency domain techniques , such as in i - ofdr , are given by the frequency sweep range , which imposes high bandwidth detection bringing up again the trade - off between resolution and dynamic range . high resolution is bandwidth demanding whereas high sensitivity imposes low resolution @xcite . incoherent ofdr measurements have been extensively studied in @xcite and rayleigh scattering has been successfully observed . the mathematical model developed in the time - domain concludes that the amplitude of the rayleigh scattering signal within an optical sub - carrier is inversely proportional to the sub - carrier s frequency ( hereby referred to as the _ high - frequency issue _ ) , which complicates signal acquisition when there is a requirement for high spatial resolution . recently , however , urban _ et al . _ @xcite showed that fiber monitoring can be accomplished in a sub - carrier multiplexed pon ( scm - pon ) context by measuring the rayleigh backscattered signal from a high - frequency ( up to 220mhz was reported ) inside a 20mhz band , which translates into a 10 meter achievable spatial resolution . the method resembles the step - frequency method presented in the late eighties by nakayama _ @xcite however , instead of disregarding the rayleigh backscattering and focusing on reflective events as in the latter , the effect of the backscattered signal was exploited at its full potential in @xcite . an alternative to the _ high - frequency issue _ of i - ofdr rayleigh backscattering measurements would be to perform the fiber monitoring in the low - frequency range , up to a few hundred kilohertz , where the @xmath0 hindering effect is not too harsh . the downside of this approach would be the diminished spatial resolution associated to one s measurement when the fourier transform of the acquired signal is taken . in order to overcome this preclusion and , at the same time , improve the signal - to - noise ratio at the receiver s end , we tackle the fiber characterization procedure in a different manner : instead of taking the fourier transform we attempt to fit the resulting backscattered signal with its associated mathematical model , which , in order to simplify its description , we model as a spatial - dependent phasor and arrive at an equivalent expression to the one found in @xcite , but in the frequency domain . by employing a fitting procedure developed in matlab@xmath1 , we show that not only a fiber fault s position but also its magnitude can be identified with high precision . the remainder of the article is laid out as follows . in section [ expset ] , we present the experimental apparatus , its characteristics , and limitations . section [ matmodel ] attempts to develop an intuitive mathematical model associating the acquired signal to an analytic frequency - dependent function . the relation between a fault s magnitude and the signal s amplitude as well as the limitation associated to eventual reflection events along the fiber are also discussed . the experimental results are presented and discussed in section [ results ] , section [ impact ] presents the technique s limitations and the impact on in - service data transmission , and section [ conclusion ] concludes the article and details our future points of investigation . the experimental apparatus is rather simple and a block diagram identifying its main structures is presented in fig . [ fig : setup ] . in order to achieve even higher signal - to - noise ratio , we perform the detection in a low - frequency lock - in amplifier which is also responsible for generating the modulating optical sub - carrier tone . according to @xcite , forcing the measuring time to be higher than the round trip time of the sub - carrier signal inside the fiber is equivalent to measuring the steady state response of the system for a given frequency . this way , determining the relative values of magnitude and phase of the backscattered signal for each frequency value in the range is the same as tracing the system s transfer function . by setting the time constant of the lock - in amplifier to @xmath2 milliseconds is sufficient to account for the steady state measurement of any practical optical fiber link available ( shorter than @xmath3 kilometers ) and translates into reasonable measuring periods ( less than 5 minutes ) given that the frequency is swept from 100 hz to 100 khz in 100 hz steps . furthermore , the sweep time is long enough so that each measurement is acquired after one time constant of the lock - in amplifier , which was verified to be short enough to avoid dragging effects on the measured curve . the noise level at each step frequency measurement is @xmath4 for the sr530 analog lock - amplifier model employed in the measurements . hz to @xmath5 khz ) modulates the laser diode . the backscattered signal corresponding to the base tone is amplified by the low - frequency photodetector s pre - amplifier and sent to the lock - in amplifier . ] the lock - in amplifier generates the sweeping tone which falls inside the 100 hz to 100 khz band and is responsible for determining the relative amplitude and phase parameters of the backscattered signal which composes the system s transfer function . note , however , that if we take into account the transfer functions of each element apart from the fiber itself and then proceed to calibrate the resulting signal , we get rid of the contribution from the laser diode , the photodetector , the electronic amplifier , and the lock in , and are left with the fiber s transfer function . the system s calibration is achieved by inducing a reflective event at port @xmath6 of the optical circulator and measuring the system s response . the laser diode operates linearly within a @xmath7 range centered at @xmath8 bias current . the low - frequency tone with angular frequency @xmath9 ( @xmath10 ) at full modulation depth has a peak voltage amplitude of @xmath11 since the laser diode s input impedance is @xmath12 . the optical signal launched into the fiber under test ( assumed to have a length @xmath13 ) will carry its modulation properties as will the backscattered signal . upon detection , the modulated backscattered signal will produce an electrical signal with current amplitude proportional to the field amplitude squared and the detector s detectivity @xcite . if we disregard the dc term after the detection , which is a good assumption since we are performing a steady state matched detection at lock - in amplifier , the electric signal entering the lock - in amplifier as a function of time can be written as @xmath14 , where @xmath15 is the detected signal s current amplitude and @xmath16 is the photodetector s detectivity .    for simplicity purposes , we first neglect possible reflective events throughout our mathematical modeling so the total signal intensity arriving at the lock - in amplifier is given solely by the incoherent sum of light scattered from all positions of the fiber . such reflective events , associated to defective connectors and fresnel reflections due to an abrupt change in the refractive index , will be input to the model after these first considerations are laid out . if we assume the distribution of scattering centers is uniform along the fiber , we can write the signal entering the lock - in amplifier as a spatial - dependent phasor in integral form : @xmath17 where @xmath18 is the launched probe signal power at the fiber s input and @xmath19 is the rayleigh scattering coefficient , assumed to be constant along the whole fiber length . @xmath20 is the complex wave vector of the modulating signal which includes the fiber s attenuation coefficient @xmath21 in its imaginary part . @xmath22 is the group index of refraction of the fiber . we now suppose a faulty _ non - reflective _ event at position @xmath23 induces an optical loss @xmath24 , where @xmath25 . to take into account this fault , we re - write eq . [ eq : fullfiber_intform ] as : @xmath26 where the first integral corresponds to power being scattered back from the portion of the fiber before the loss event whereas the second integral corresponds to power scattered back from the last portion of the fiber , after the fault event at @xmath23 . the loss factor @xmath24 appears multiplying the second integral twice ( hence , squared ) due to the fact that light passes two times through the fault in its round trip and we assume that the transmission coefficient is symmetric ( it is the same whether light is propagating from the generator to the load or in the opposite direction ) . to simplify the expression , we condensed the constants into @xmath27 . note that the second term in eq . [ eq : fiberbreak_intform ] can be dismembered into two integrals with different limits @xmath28 which , when substituted back into eq . [ eq : fiberbreak_intform ] , yields @xmath29 the signal arriving at the lock - in amplifier can , thus , be interpreted as the combination of two parcels : the modulated backscattered signal from a fiber with length @xmath23 , which contributes with a factor @xmath30 to the overall optical power ; and the modulated backscattered signal from a fiber with length @xmath13 , which contributes with a factor @xmath31 to the overall optical power .    splitting the expression of @xmath32 into two parcels and dealing with the integrals separately , we find that the results are complex sinusoidal functions with the oscillatory parameter ( the product @xmath33 which governs the oscillation of the sinusoidal function ) being the distance of the fiber : @xmath34 where we have condensed the factors @xmath35 and @xmath36 into @xmath37 and @xmath38 , respectively , and @xmath39 . we would like to highlight that this modelling approach is rather intuitive as it has a somewhat direct correspondence to otdr - like traces , i.e. , the signal contributions from different portions of the fiber are independently weighted and piled up forming the resulting trace . the difference here is that , instead of step functions as in the standard otdr , the signal contributions are sinusoidal functions that oscillate more or less within the frequency range of the sweep depending on the fiber length to which they are associated .    in @xcite , a similar expression has been derived following a distinct mathematical framework . there , however , the phase factor of the monitoring signal has been disregarded and the authors focused on the magnitude of the phasor in order to obtain an otdr - like trace by employing a fourier transform . as we show in the next section , the phase factor is rather sensitive to changes in the fiber and , if used in conjunction with the magnitude information , may raise the precision of fault location . in order to take both the phase and magnitude informations into account , a fitting procedure in the complex frequency domain is necessary and , although the result is not presented as an otdr - like trace but , rather , as an event list , it is possible to extend the spatial resolution beyond the limit imposed by the fourier transform . although the fourier transform operator is unambiguous , there are some pitfalls that may arise from computational methods : the non - localizability of the fourier transform operator causes it to interpret all frequency components with the same weight ; for that reason , there are several methods in the literature that focus on enhancing the localizability of the fourier transform operator through signal processing in order to enhance the spatial resolution of frequency measurements ( @xcite and references therein ) . such methods , however , tend to be very expensive in computational complexity even though the result is equivalent to the fitting procedure demonstrated in this document . the fitting procedure on the frequency domain is much more attractive in our application since the phase information is highly dependent on the position of the break which translates into a computationally cheap and accurate method for determining the position of a fault thus enhancing the achievable spatial resolution with respect to a straightforward fourier transform . the result conveyed by eq . [ eq : int2sinc ] is that @xmath32 is the sum of two sinusoidal terms which take different arguments corresponding to the length of fiber traversed by the backscattered optical signal . from the expressions of eq . [ eq : int2sinc ] , we may attempt to write a generic term @xmath40 for the backscattered signal from a fiber of arbitrary length @xmath41 . this generic term is a function of both the frequency and the fiber length and has the form presented below , where we took the liberty of expanding and simplifying all terms : @xmath42 note that , for an arbitrary number of non - reflective faults along the fiber link , the procedure presented in eq . [ eq : dismember ] can be used indefinitely . moreover , the backscattered signal from an optical fiber link which contains a number @xmath43 of fault events along its length will be represented as the sum of @xmath43 generic terms presented in eq . [ eq : genericterm ] , each dependent on the fault position @xmath44 and with different amplitude factor @xmath45 , so we may write : @xmath46    the expression of @xmath45 gets progressively intricate as the number of fault event rises as can be perceived from the right - hand side of eq . [ eq : integral_form2 ] : with @xmath43 terms , for instance , the value of @xmath47 depends on the values of all @xmath48 from @xmath49 downto @xmath50 which also holds true for @xmath51 and so forth . extending the expression of @xmath45 for any number of fault events yields the following : @xmath52    once the non - reflective events have been considered and modeled , we turn our attention to the parcel associated to reflective events . observing reflections along an optical fiber link has been a matter of extensive research and , as previously commented , nakayama _ et al . _ were able to detect such reflections with a similar experimental apparatus @xcite . in a way , therefore , reflection events are manageable and the method hereby proposed offers a means of not only identifying these specific events , but also the non - reflective ones . as we know , reflective events are always accompanied by a loss of intensity due to energy conservation conditions ( namely , the @xmath53 condition at any partially reflective interface ) but the converse does not hold , i.e. , an intensity loss may or may not be accompanied by a reflection . generally , loss events that are not associated to a reflection peak in standard otdr measurements are due to absorption or irradiation effects , such as macro- and micro - bend losses @xcite . reflection events are more easily observed since the _ high - frequency issue _ does not arise in modeling its amplitude and a sinusoidal function which reflected back will arrive at the lock - in amplifier as a sinusoidal function with lesser magnitude and constant phase . hence , taking reflective events into account is as straightforward as including a constant magnitude and constant phase factor to the model . we , then , rewrite eq . [ eq : s_final ] as @xmath54 with @xmath55 , meaning that the set of reflective events is a subset of all the events but not necessarily equal to the latter .    in the same fashion as nakayama _ considered the rayleigh backscattering some sort of noisy contribution that could spoil the reflection measurements performed @xcite , we consider the reflection events as noisy contributions that can spoil our fiber characterization method . suppose , for instance , that a optical fiber link presents a non - reflective event at position @xmath23 and a reflective fiber - end at position @xmath13 . if the reflective fiber - end contribution rises above the backscattering power level , the non - reflective event wo nt be observed even though it is there . we first wish to validate the mathematical modelling of the signal arriving at the lock - in amplifier . for that , we simulated the results from eq . [ eq : s_final ] using @xmath56 , @xmath57 km , and a frequency range of @xmath5hz to @xmath5khz . we then compared the amplitude and phase data with the results acquired from measuring a @xmath58-km standard optical fiber within the same frequency range using our proposed method . the amplitude of the modulating signal throughout all experiments were kept at one eighth of the full modulation depth . the result is depicted in fig . [ fig : ampphasesimulationcomp ] . -km standard optical fiber . a ) amplitude data series ; b ) phase data series . ]    even though the results show good agreement between simulated and experimentally acquired data , there is still a slight discrepancy , specially in the phase data series . note , however , that the generic term presented in eq . [ eq : genericterm ] depends on @xmath59 and @xmath41 but also on @xmath21 and @xmath22 , the fiber s attenuation coefficient and the group index of refraction of the fiber , respectively . since , from an optical link monitoring perspective , there is usually no _ a priori _ access to such information , we utilized standard values within the international telecommunication union ( itu - t ) specification range , i.e. , @xmath60db / km and @xmath61 @xcite . nevertheless , these values may suffer variations from one fiber to the other and the discrepancy perceived between simulated and experimental data may be attributed to that . observe that , from a systematic point of view , if one was capable of identifying which frequency - dependent sinusoidal functions are present in a data series originated by the monitoring procedure described  which is equivalent to finding all the generic terms of eq . [ eq : genericterm ] that compose the signal  , the eventual positions of faults could be determined . before that , and in order to further investigate the correspondence between the mathematical model and the experimentally acquired data , as well as determining the availability of this fiber characterization method , we conducted laboratory tests referenced by a standard otdr device with previously unknown fibers ( unknown length and attenuation coefficient ) . first , the amplitude and phase profiles were extracted from the fiber in its original configuration using the lock - in amplifier . afterwards , a fiber bend loss was induced in an arbitrary position also with unknown magnitude . the goal is to identify how a fiber bend alters the transfer function s profile , so the amplitude and phase profiles of the modified fiber were extracted using the lock - in amplifier .        from the results of fig . [ fig : ampphasecomp]-a and fig . [ fig : ampphasecomp]-b , it is clear that tempering with the fiber has a direct effect on the signal measured with the lock - in amplifier , so the possibility of quantifying the effect by employing our knowledge of the mathematical model is a straightforward assumption .    to assess the capability of identifying a fault event along an optical fiber link , the otdr trace of the unknown fiber was acquired in its original and fault - induced conditions . at the same time , a model fit of the amplitude and phase profiles was ran based on an _ extensive search _ heuristic developed in matlab@xmath62 : the algorithm sorts an arbitrary fiber length and calculates the generic component ( eq . [ eq : genericterm ] ) of our mathematical model ; it then calculates the minimum square error between the original profile ( the fiber without the induced fault ) summed with this generic component and the profile after the induced fault ; the position ( @xmath44 ) , amplitude factor ( @xmath45 ) , and attenuation coefficient parameters are adjusted to find the candidate that best fits the new result . note that the nature of the mathematical model enables one to test an arbitrary number of generic components in order to best estimate the fiber fault position . in a sense , the spatial resolution is limited by the amount of generic components one is inclined to test since the higher this number , the longer the algorithm will take to come up with a result . in the examples that follow , the generic components were created each separated by 1 meter , i.e. , for a 4-km fiber , four thousand components will be tested and the one that best fits the acquired signal will be selected . the non - linear search algorithm employed in matlab@xmath1 , however , must deal with four parameters at a time ( position , magnitude , @xmath21 , and @xmath22 ) so the average processing time revolves around @xmath5 milliseconds . increasing this separation over 1 meter will be computationally costly and ever slower , for which reasons we set 1 meter as our default separation between generic components achieving a @xmath63 minutes computing time . the developed algorithm along with a testbench data series is available at @xcite . we tested two different optical links : link @xmath64 has its amplitude versus phase profile presented in fig . [ fig : ampphasecomp ] and corresponds to a @xmath65-km long composition of two fibers with the induced fault at @xmath66 kilometers . link @xmath6 corresponds to a @xmath67-km composition of two fibers with the induced fault at @xmath68 kilometers . the results are shown in table [ table : comp1 ] , for which we found a less than @xmath69 meters difference between the standard otdr and the low - frequency sweep measurements in both cases . this results are a clear indication that the fourier transform limitation of i - ofdr has indeed been overcome since a 100khz frequency sweep would allow for a maximum spatial resolution of 2 kilometers .    .standard otdr vs low - frequency sweep [ cols="^,^,^,^,^ " , ]     [ table : comp2 ]    the fiber end contribution varies inversely , as expected . according to these results , the dynamic range is limited to a maximum of @xmath70 db , i.e. , the maximum attenuation throughout the link which still permits one to identify the fiber end . in the access systems of rof , the optical link distribution does not , usually , exceed 15 km . therefore , even though the technique achieves low dynamic range , it is an interesting candidate for short range mobile fronthaul applications with low optical splitting ratios @xcite . the possibility of in - service network monitoring has been assessed by measuring the resulting error vector magnitude ( evm@xmath71 ) for different data signal amplitude under two conditions : with ; and without the monitoring tone . once again , the monitoring tone amplitude was kept at one eighth of the laser s full modulation depth in order to simulate an eight - channel scm network . the result indicates that the monitoring signal deteriorates the quality of the communication , even if by a very small factor . under these conditions , however , the required evm@xmath71 for lte communication is respected with a large margin @xcite . it can be observed that , above a certain value of the data signal s amplitude , the quality of communication is deteriorated . this behavior is associated to the laser s non - linear response region , achieved when the modulation signal violates the laser s full modulation depth at around 1.5 [ vpp ] , or 30 [ ma ] . the recent interest in the deployment of scm - pon and cost restrictions imposed by modern commercial needs bring up low cost embedded monitoring techniques as an essential tool for operational success of access networks . although a method capable of replacing one optical sub - carrier channel by a monitoring probe has been presented in @xcite , dynamic range is limited by the 1/f factor impairing backscattering measurements inside an optical sub - carrier . the results presented here shows that an embedded in - service short - reach link monitoring can be performed and the impact over data transmission is expected to be not higher than the cross contribution from any two sub - carrier channels due to the probe signal s very nature . using the baseband for monitoring purposes , the method benefits of the high - sensitivity low - frequency detection and further overcomes the _ high - frequency issue _ , fully preserving the rf sub - carrier band for data transmission . fulfilling the steady state condition over the round trip path along the fiber during the measurement is imperative for the performance of the technique and a feature that distinguishes it from the technique presented in @xcite . nevertheless , minute - long measurement times are achieved in fibers a few kilometers long . the mathematical model presented here is similar to the one previously developed in @xcite but with a different approach , namely the development in the frequency domain when one represents the incoming signal as an integral of phasors originated from the modulated light coming from distinct portions of the fiber . the intuitiveness of the model is useful since the signal is not transformed back to the spatial domain as in standard c- or i - ofdr and the results are given in the form of an event table as opposed to an otdr - like trace . the algorithm developed to treat the signal , even though simple and still limited enables highly reliable fault location , as backed by our experimental results contrasted to standard otdr measurements . this is possible due to the fact that several signal component candidates may be tested and the one that best approximates the acquired complex signal in the frequency domain is selected . using the phase combined with the magnitude information of the backscattered signal , we could overcome the fourier transform spatial resolution limitation of low - frequency measurements . we believe that , as the method developed by @xcite , the method presented here may also figure as a good candidate for embedded in - service short - reach scm - pon monitoring , not demanding a sub - carrier channel for that purpose . assuming that the data channels are delegated to highest frequencies ( in the order of hundreds of megahertz ) , the base band would be unoccupied and could be reserved for the low - frequency probing tone herewith described . we also believe that the standard deviation reported in figure [ fig : estimatevsreal ] is the limitation of localizability of the phase measurement even when the fitting procedure in the frequency domain is employed . even though we surpass the standard fourier transform limitation , there is still a fourier transform limitation to the technique , i.e. , we can not retrieve perfect information of the fiber unless an infinite frequency span is performed . unfortunately , the lock - in amplifier employed in the measurements is incapable of reaching higher than @xmath5 khz frequencies and , therefore , we are limited to this spatial resolution of @xmath72 meters . higher frequency spans should improve the resolution , so conducting the experiments with a device specified for higher frequencies is a goal for future developments . location of multiple faults is , unfortunately , still impracticable with the present signal processing technique . we believe , however , that it could be achieved and leave it as our major future point of investigation . the innovation centre , ericsson telecomunicaes s.a . , brazil and brazilian agencies cnpq and faperj supported this work . g.  c. amaral , l.  e.  y. herrera , d.  vitoreti , g.  p. temporo , p.  j. urban , and j.  p. von  der weid , `` wdm - pon monitoring with tunable photon counting otdr , '' _ photonics technology letters _ , vol . 26 , no .  13 , pp . 49214926 , 2014 .        g.  c. amaral , j.  d. garcia , l.  e. herrera , g.  p. temporao , p.  j. urban , and j.  p. von  der weid , `` automatic fault detection in wdm - pon with tunable photon counting otdr , '' _ lightwave technology , journal of _ , vol . 33 , no .  24 , pp . 50255031 , 2015 . j.  p. von  der weid , r.  passy , g.  mussi , and n.  gisin , `` on the characterization of optical fiber network components with optical frequency domain reflectometry , '' _ lightwave technology , journal of _ , vol . 15 , no .  5 , pp . 11311141 , 1997 . p.  j. urban , g.  c. amaral , and j.  p. von  der weid , `` fiber monitoring using a sub - carrier band in a sub - carrier multiplexed radio - over - fiber transmission system for applications in analog mobile fronthaul , '' _ lightwave technology , journal of _ , vol . 34 , no .  13 , pp . 31183125 , 2016 . j.  a. tropp , j.  n. laska , m.  f. duarte , j.  k. romberg , and r.  g. baraniuk , `` beyond nyquist : efficient sampling of sparse bandlimited signals , '' _ ieee transactions on information theory _ , vol . 56 , no .  1 , pp . 520544 , 2010 . l.  herrera , g.  amaral , and j.  von  der weid , `` investigation of bend loss in single mode fibers with ultra - high - resolution photon - counting optical time domain reflectometer , '' _ applied optics _ , vol . 55 , no .  5 , pp . 11771182 , 2016 . , `` a low - frequency tone sweep method for fault location in optical fiber links  matlab source files and testbench data series , '' https://github.com/tungstenio/tonesweepfibarlocation , accessed : 2016 - 05 - 07 .
Abstract:

Answer: we demonstrate an optical fiber fault location method based on the frequency response of the modulated fiber optical backscattered signal in a steady state low - frequency step regime . careful calibration and measurement allows for the reconstruction of the fiber transfer function , which , associated to its mathematical model , is capable of extracting the fiber characteristics . the technique is capable of identifying non - reflective fault events in an optical fiber link and is perfectly compatible with previous methods that focus on the reflective events . the fact that the recuperation of the complex signal is performed in the frequency domain and not via a fourier transform enables the measurements to overcome the spatial resolution limitation of fourier transform incoherent - ofdr measurements even with frequency sweep ranges down to 100 - 100000 hz . this result is backed up by a less than 10 meters difference in fault location when compared to standard otdr measurements .


Instruction: Write an abstract for the following acadamic paper.

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

Paper: in spite of quantum mechanics has demonstrated the amazing power to describe micro - phenomena , its behavior is so different from macro - intuitions that attempts to recover ( quasi- ) deterministic movement of particles never give up.@xcite however , due to non - local correlation of quantum states , these attempts all suffer quite difficulties to achieve self - consistent results and received a despaired strike by aspect s experiments@xcite on epr correlations.@xcite many evidences have revealed that non - locality ( or entanglement ) is a fundamental property of quantum process , and to restore the motion of particles back to deterministic manner is in fact impossible . nevertheless , one may ask a further question : is such kind of non - local correlation globally ( i.e. , simultaneously around the whole space ) correlated or partially localized ? from a point of view of standard non - relativistic quantum mechanics , since wave function is defined on the whole space , as well as all operators , the theory must be globally correlated . in terms of feynman s path integral formalism , that is to say that when calculating the transition function between two quantum states @xmath0\mathcal{d}[p(t ) ] \nonumber\\ & & \times \exp\left[\frac{i}{\hbar}\int_{t_a}^{t_b}\mathrm{d}t\left\{\sum_{r } p_{r}(t)\frac{\mathrm{d}}{\mathrm{d}t}q_{r}(t)-h\left(q(t),p(t)\right)\right\}\right ] , \label{eq : fpi}\end{aligned}\ ] ] all paths across the space are allowed,@xcite which results in globally and instantaneous corrections . here path integral @xmath1\mathcal{d}[p(t)]$ ] means to sum over the whole coordinate and momentum spaces along time - sliced paths with fixed end points at @xmath2 and @xmath3 ( see chapter 1 of ref .   for details ) . thus one may wonder how can we separate an atomic ( quantum ) system from the universe and treat it isolatedly as we did heaps of times before ? a consistent theory should explain why to ignore effects due to supernova explosions or other processes occurring at somewhere in the cosmos is reasonable . it is difficult within the framework of standard quantum mechanics within the framework of path integral formalism , one may would like to employ the stationary phase condition to argue that contributions from those paths with large deviation from classical one s canceled completely and results in required isolation conditions . but notice ( a ) this is only true for quasi - classical approximations ; ( b ) obviously different paths would have different variations of the contribution due to changes of the dynamic environment along each course , there is no guarantee that all of these variations with time can be completely and simultaneously canceled all the time . in a sense of that , we reach a point that some paths in eq.([eq : fpi ] ) should be forbidden in a self - consistent quantum theory in order to satisfy the isolation hypothesis and to remove instantaneous correlations . actually , it will be seen in this paper that to generalize the feynman s formalism slightly to covariant relativistic case will produce partial locality to quantum mechanics . this property is crucial to ensure that experiments on atomic system are really corresponding to the behavior of that system but not those relating to the external cosmos . it is necessary to point out that several attempts to set up a relativistic particle path integral have been proposed before , but they failed to produce the required isolation condition . the employed formalisms are rather unsatisfactory due to lack of explicit lorentz covariance and treated time and spatial coordinates on different footings.@xcite the physical implications associated with the derivations also are not clear enough . therefore these attempts are in fact a kind of mathematical techniques for conventional quantum mechanics , and can not be taken as the foundation for developing a new theory . in subsequent sections , a variational principle on world lines is developed , which is necessary in order to generalize feynman s formalism to lorentz covariance . the relationship with standard quantum mechanics is discussed in section [ sec - cpi ] . some understandings arised from this new point of view are given , followed by a comparison with previous derivations . thanks to the theory of relativity , time is deprived its privilege in motion equations and became an ordinary dimension of space - time . if one still would like to use lagrangian variational principle to recover the position of action principle in physics , finding out an other parameter to take the place of time in newtonian mechanics is necessary . considering the movement of a particle with non - zero rest mass in space - time , we always can choose a frame of reference in which the particle is at rest . usually it is non - inertial . however , it is reasonable to assume that instantaneously one can employ an inertial frame to approximate it ( in terms of the general theory of relativity , the space - time manifold can be approximated by a series of minkowski spacetimes locally@xcite ) . in this way the movement of a particle can be described by a series of inertial frames , which relating to the observer by respective lorentz transformations @xmath4 here @xmath5 is the ratio of instantaneous velocity to light speed , i.e. , @xmath6 . under these transformations , the length of proper time @xmath7 is invariant and satisfies @xmath8 , where @xmath9 is the metric of minkowski space . usually @xmath7 can be expressed as a functional as @xmath10 . thus the variation of frames characterized by parameter @xmath5 can also be characterized by proper time @xmath7 . namely , with eq.([eq : lorentz - t ] ) we have @xmath11 . in this way the proper time takes the privilege position of time as in newtonian mechanics , and the motion equations of @xmath12 along @xmath7 can be obtained by the generalized principle of least action that defined on world lines . ) fails and the equations of motion is independent from proper time . it corresponds to a static case in newtonian framework . ] suppose a relativistic dynamical system that can be determined by a characteristic functional @xmath13 , where @xmath14 , then one may have a formal action defined on world lines @xmath15=\int_{\tau_{i}}^{\tau_{f}}l(x_{\mu},\dot{x}_{\mu})\mathrm{d}\tau \label{eq : action-1}\ ] ] with fixed end points at @xmath16 and @xmath17 . the variation of this action with respect to world lines with fixed boundary conditions @xmath18 and @xmath19 gives @xmath20 when @xmath21=0 $ ] . this is just the generalized covariant formalism of least action principle , and one can employ eq.([eq : c - le ] ) to reproduce the covariant motion equation of particles .    using the formal conjugate momentum @xmath22 , we can define a formal hamiltonian by a legendre transform @xmath23 usually it is lorentz invariant . then the equation of motion can also be given by @xmath24 it is evident from above discussion that if @xmath25 is not explicitly @xmath7-dependent , then the formal hamiltonian is a conservative quantity . when @xmath12 and @xmath26 varied independently , the formal action of eq.([eq : action-1 ] ) can be expressed as @xmath27=\int_{\tau_{i}}^{\tau_{f}}\left[\sum_{\mu}p^{\mu}\dot{x}_{\mu}-m\left(x_{\mu},p^{\mu}\right)\right]\mathrm{d}\tau , \label{eq : action-2}\ ] ] which is lorentz invariant too .    for a free particle , one has @xmath28 for 4-momentum @xmath26.@xcite it is easy to prove that the corresponding formal hamiltonian should be @xmath29 then one gets @xmath30 where @xmath31 is the rest mass because @xmath32 . the conservation law of @xmath25 in fact becomes the conservation of rest mass or rest energy . in this formalism we see that mass no longer appears as an _ aprior _ parameter but a dynamical variable analogous to hamiltonian in classical physics . as discussed in introduction section , self - consistent quantum mechanics should be at least partial localized . if it is true , there is no any reason to save the conception that wave function ( or quantum state ) exclusively describes just the motion of that particles . this clew of thought eventually leads to a theory of quantum fields . however , here we would like to restrain our discussions within quantum mechanics .    ) where paths must satisfy @xmath33 ( domain @xmath34 ) ; ( b ) spatial summation lattice between time @xmath35 and @xmath36 ( domain @xmath37 ) , not only all paths belonging to end points @xmath38 and @xmath39 are restricted by the dot - dashed rectangle , every path also should be within the instantaneous light cones all the time to ensure a properly definition of proper time . [ fig : cpi],title="fig:",width=230 ] ) where paths must satisfy @xmath33 ( domain @xmath34 ) ; ( b ) spatial summation lattice between time @xmath35 and @xmath36 ( domain @xmath37 ) , not only all paths belonging to end points @xmath38 and @xmath39 are restricted by the dot - dashed rectangle , every path also should be within the instantaneous light cones all the time to ensure a properly definition of proper time . [ fig : cpi],title="fig:",width=172 ]    following feynman s consideration , it is quite natural to assume that the transition probability amplitude between two quantum states is completely determined by a phase factor given by the formal action of particle eq.([eq : action-2 ] ) , which results in the covariant version of eq.([eq : fpi ] ) @xmath40\mathcal{d}[p(\tau ) ] \nonumber\\ & \times \exp\left[\frac{i}{\hbar}\int_{\tau_a}^{\tau_b}\mathrm{d}\tau\left\{\sum_{\mu,\nu}p^{\nu}(\tau)\frac{\mathrm{d}}{\mathrm{d}\tau}x_{\mu}(\tau)\delta_{\mu\nu}-m\left(x_{\mu}(\tau),p^{\mu}(\tau)\right)\right\}\right ] . \label{eq : c - fpi}\end{aligned}\ ] ] evidently , this is just a natural special relativity generalization of feynman s path integral formalism . the proper time sliced expression is given by ( differs from other path integral formalisms where a time - slicing scheme is used ) @xmath41,\end{aligned}\ ] ] where @xmath42 and the integral domain @xmath43 of space - time coordinates is not only determined by boundary conditions @xmath44 and @xmath45 , it also has to satisfy a condition of that @xmath46 should be defined properly which constrained by physical conditions @xmath33 and @xmath47 . thus we have @xmath48 , corresponding to an intrinsic property of minkowski space - time and the non - spatial like motions of particle , respectively . it is clear that the most significant difference between this and the conventional path integral formalism lies on that in the latter case , all paths across the space are possible , while eq.([eq : c - fpi ] ) allows just some special paths , as showed in fig.[fig : cpi ] . for simplicity , only ( 1 + 1)-dimension case is shown here , where time - component of proper time sliced paths are restricted by @xmath33 and spatial component belonging to end points @xmath38 and @xmath39 are not only limited within the dot - dashed rectangle ( relativistic causality ) , but also each path should be time - like , i.e. , contained in the instantaneous light cones all the time , to guarantee @xmath49 at classical limit ( see eq.([eq : c - fpi ] ) ) . evidently , this property recovers the global causality of relativity and partial locality to the theory , the latter eventually leads to the isolation condition for massive particles . it is easy to see that eq.([eq : c - fpi ] ) contains conventional quantum mechanics as approximation as feynman case shown . the superposition and composition law of the transition probability amplitude is stated in a form of @xmath50\mathcal{d}[p]=\int_{\mathbb{c}}\mathrm{d}^{4}x'\int_{x'_{\mu}\left ( \tau'\right)}^{x_{\mu}\left ( \tau_{b}\right)}\mathcal{d}[x_{\uppercase \expandafter { \romannumeral 2 } } ] \mathcal{d}[p_{{\uppercase \expandafter { \romannumeral 2 } } } ] \int_{x_{\mu}\left ( \tau_{a}\right)}^{x'_{\mu}\left ( \tau'\right)}\mathcal{d}[x_{{\uppercase \expandafter { \romannumeral 1 } } } ] \mathcal{d}[p_{{\uppercase \expandafter { \romannumeral 1 } } } ] .\ ] ] it is worthwhile to note that j. schwinger once introduced proper time depended green functions to trick the gauge invariant problem,@xcite which is different from our considerations here completely .    consider the transition function of eq.([eq : c - fpi ] ) with an infinitesimal increase in proper time @xmath7 , the wave function can be simplified as @xmath51 , \label{eq : c - fpi - small - t}\end{aligned}\ ] ] with @xmath52 . as matrix element of operator @xmath53 is given by@xcite @xmath54\mathcal{d}[p(\tau)]o(x_{\mu}(\tau))\exp\left[\frac{i}{\hbar}s(\tau_{b},\tau_{a})\right],\ ] ] it is quite easy to prove that the corresponding form of operators acting on wave function eq.([eq : c - fpi - small - t ] ) are @xmath55 to the first order of @xmath56 when @xmath56 approaches zero . note here that the @xmath57-momentum operators have an oppositive sign comparing to traditional definition , but it does not matter . using eq.([eq : xt - operat ] ) we have the motion equation for wave function as @xmath58 for a formal hamiltonian eigenstate @xmath59 , one gets @xmath60 with @xmath61 . if studied system only includes gauge interactions , then the formal hamiltonian of particle part still keeps the form of eq.([eq : form - h ] ) except for a slight modification @xmath62 . so generally we can use the form of eq.([eq : form - h ] ) for discussion . applying one more @xmath63 on eq.([eq : eigeq - m ] ) and making use of eqs.([eq : form - h],[eq : xt - operat ] ) and @xmath64 , we reach the klein - gordon equation . a detailed prove is given in appendix a to show this formal inference is true . in order to obtain a linear equation , let s consider a mapping from @xmath57-vector to a lorentz scalar @xmath65 . conservation of the length requires @xmath66 with anticommutation defined by @xmath67 . obviously @xmath68 are exactly the dirac matrices . in this sense @xmath69 form a set of orthonormal basis with inner product given by anticommutation operation , and @xmath57-vector @xmath12 can be viewed as the components of lorentz scalar @xmath70 on this basis set @xmath71 . then the formal action can be reexpressed as @xmath72=\int_{\tau_{i}}^{\tau_{f}}\left[\frac{1}{2}\{p,\dot{x}\}-m\left(x , p\right)\right]\mathrm{d}\tau , \label{eq : action-3}\ ] ] with @xmath25 in eq.([eq : form - h ] ) given by @xmath73 , which is a lorentz scalar but usually is a matrix too . then the dirac equation is obtained directly from eq.([eq : eigeq - m ] ) . a derivation of non - relativistic approximation is given in appendix b , which reproduced the standard non - relativistic path internal formalism for free particles , demonstrating the validity of our proposal .    by comparing with previous attempts to the same issue , our proposed scheme is more systematic and elegant with high space - time symmetry . though all proposed methods have ability to reproduce relativistic equation of motion , the works of d. felea@xcite and p. gosselin@xcite are not lorentz covariant explicitly . without introducing proper time , they treated time and spatial coordinates non - equivalently , which is contrary to the spirit of relativity . t. padmanabhan@xcite employed a quite different approach and introduced proper time , unfortunately it is still non - covariant . in pazma s method even the action is non - invariant.@xcite the only covariant version is originated from feynman , but it suffers a difficulty that the physical meanings of the procedure and the conjugated variable of proper time are unclear . in a sense of that , these path integral schemes are actually just convenient mathematical techniques for quantum mechanics , and can not be treated as the basis for developing a new theory . by contrast , our proposed formalism is not only lorentz covariant , but also bears clear physical implications by developing the rest energy operator @xmath63 , which results in the evolution of states controlled by @xmath63 besides the conventional hamiltonian , therefore has the potential to treat variable rest mass problem . moreover , all previous proposed formalisms are just defined on coordinate space , it is well known that this kind of path integral is not always valid , and including phase - space integrals , as in our scheme , is more appropriate.@xcite it would be interesting to show how our formalism can produce the partial locality for quantum process , for example , measurements , which is beyond the scope of application of other proposals . for simplicity , let s assume the measurements are just perturbations to the state , and the corresponding action can be characterized by a dirac delta function as @xmath74 . then the wave function with two measurements performed respectively at space - like points of @xmath75 and @xmath76 is given by @xmath77\exp[\frac{i}{\hbar}\{s_{0}+\alpha s'\delta(x_{u}-x'_{\mu})+\beta s''\delta(x_{u}-x''_{\mu})\}].\ ] ] here we have suppressed unrelated expressions for brief , for example the path integral in phase space and the detailed information about actions . parameters @xmath78 and @xmath5 must be small quantities to ensure measurements are perturbation . making taylor expansion of the exponential to first order of @xmath78 and @xmath5 , we get @xmath79 the first term at the right hand side is the wave function without perturbations and the last two terms corresponding to influences of measurements whose subscripts @xmath80 and @xmath81 indicate the domains of definition , which given by @xmath82\psi(x_{\mu}^{0},\tau_{0})\alpha s'\delta(x_{u}-x'_{\mu})\exp[\frac{i}{\hbar}s_{0 } ] . \label{eq : local-1}\ ] ] similar expression holds for @xmath83 . it is interesting to notice that the paths with non - zero contributions to eq.([eq : local-1 ] ) must pass through point @xmath84 . on the other hand , the integral domain @xmath43 is determined by non - spatial condition and @xmath33 condition , the latter allows us to divide the integral as @xmath85=\int_{\mathbb{c}'}\mathcal{d}[x]+\int_{\mathbb{c}''}\mathcal{d}[x]\ ] ] where @xmath86 only contains paths with @xmath87 and other paths belong to @xmath88 . evidently , the paths of @xmath88 must contain contributions arised from @xmath89 , which describes a particle moving along reverse time , that is to say , it corresponds to an anti - particle .    at low energy approximation where no anti - particles are excited , we can ignore the paths belonging to @xmath88 safely . in this case , at time @xmath90 , the value of wave function are determined by integral over those paths belonging to @xmath86 and passing trough @xmath84 which constrained by non - spatial condition @xmath91 , see fig.[fig : cpi ] . this means that @xmath92 does not be defined around the whole space - time . to enlarge its definition domain to the whole space - time , we have to use a step function @xmath93 that has a value of 1 if @xmath94 otherwise 0 to characterize the definition domain @xmath80 explicitly . the result reads @xmath95 where @xmath96 . similar expression holds for @xmath83 . thus if define @xmath97 , then when @xmath98 one gets @xmath99 that is , the influence of two measurements occured at @xmath84 and @xmath100 are physically uncorrelated . however , when @xmath101 , it will become correlated since the intersection of @xmath80 and @xmath81 usually non - empty in this case . it is quite clear that in this case the propagation of quantum correlation of measurements also has the same speed as light ( if the space - time coordinates of one measurement lies within the influence region of another measurement , then their results will be correlated , otherwise are independent ) . however , if paths in @xmath88 are involved , i.e. , with anti - particle excitations , the segments of a path with @xmath102 allow the path sweeping a larger time interval than actually happened from the point of an observer . this implies that the permitted spatial region in fig.[fig : cpi](b ) are broadened , which would lead to a superluminal propagation of influence region of measurements . at high energy limits where excited anti - particle can exist permanently , quantum correlation then will approach globally and instantaneously . this novel mechanism for propagation of quantum correlation makes it possible to interpret the experiments on epr effect@xcite with a new point of view . a. aspect s experiments shows a contradiction between the reality and locality in a physical theory.@xcite however since there the locality is defined by the principle of maximum speed of light ( i.e. , the velocity of quantum correlation propagation is limited by @xmath103 ) , the propagation mechanism of perturbations arised from our formalism releases this definition and solve the contradiction automatically for massive matters . as long as energy does not approach infinite , the life - time of excited anti - particles must be finite , then we have a finite correlation speed and the required isolation condition for quantum mechanics . if we define this property as locality , which has the same spirit as the original one  to get rid of action at a distance from physical theory , then bell s inequality@xcite loses its availability to characterize the relationship of physical reality and locality because we do nt know exact information about anti - particle excitations , as well as the exact value of correction speed , despite it should larger than the velocity of light usually . we see that in this case we still have a possibility to preserve the reality and locality from epr s experiments . of course this also leaks a light of hope for hidden variable theories if they prefer to include anti - particle effects in their framework to reproduce superluminal mechanism , though i do not believe it could be a right choice for the truth . it is necessary to point out that in aspect s experiments , since correlation between photons are examined , which corresponding to a singular case in our formalism because photon is massless and its anti - particle is just itself , the predicted correlation velocity will approach infinite . this should be the only one case where action at a distance may be possible in our formalism and results in a maximum violation of bell s inequality . however , to eliminate this possibility is far beyond the scope of this paper . regards this , we suggest to perform analogous experiments on massive particles to examine the effects of anti - particle excitations . using developed variational principle defined on world lines , we obtained the formal action for a relativistic system . it permits us to generalize the feynman s path integral formalism of quantum mechanics to special relativity case elegantly , which reproduces the klein - gordon and dirac equation for infinitesimal intervals very well . some novel understandings are obtained from this new point of view . though quantum state can be defined on the whole space - time , quantum processes are always localized at low energy limits , while large - scale correlations are due to anti - particles excitations , which makes superluminal corrections become possible even within the framework of special relativity . with this mechanism for correlation propagation , quantum mechanics becomes more completeness and self - consistency . also it shows at the first time that we still have opportunity to save the concept of physical reality without lost the locality of a theory . sometimes this kind of correlation mechanism will bring a little inconvenience to calculation . for example , when considering a process described by a superposition state , one then need to check that if its sub - states are really correlated or not within studied time interval in order to give physical meaning to the final superimposed state . at last , it is necessary to note that most discussions in this paper are for single particle , but the same conclusions hold for many - particles case , which is rather straightforward to derive and do nt need to repeat here . the author thankfully acknowledges the help and kind hospitality from prof . ying chen at the university of tokyo . here we would like to derive the klein - gordon equation of free particle explicitly by starting from eq.([eq : c - fpi - small - t ] ) , which can be rewritten as @xmath104 , \label{eq : a - fpi}\end{aligned}\ ] ] here we have used a condition of @xmath105 to remove the square root in action , which means that we are considering a particle with fixed rest mass . note the 4-momentum in eq.([eq : a - fpi ] ) should be treated independently when making path integral , namely , it does not satisfy the mass - energy relation and corresponding to none real physical state . integral out the momentum , one gets @xmath106 . \label{eq : a - coord - fpi}\ ] ] here use has being made of @xmath107 . as @xmath108 is small and making the taylor expansion of @xmath109 to its second order , we have @xmath110 while considering that the integral domain @xmath43 is central symmetry with respect to the origin of coordinates , non - zero contributions to the integral are just come from terms of @xmath111 . thus we can write the right hand side of eq.([eq : a - coord - fpi ] ) as a sum of two terms , one depends on @xmath109 and the other depends on second order derivatives of @xmath109 . let s calculate the first term firstly , i.e. , @xmath112\nonumber\\ & = \frac{im_{0}^{2}}{4\pi^{2}\hbar^{2}\varepsilon^{2}}\psi(x_{\mu},\tau ) \int_{\mathbb{c}_{1}}\mathrm{d}(c\delta t)\int_{\mathbb{c}_{2}}\mathrm{d}^{3}(\delta \vec x)\,\exp\left[\frac{i}{\hbar}\frac{m_{0}}{2\varepsilon}\sum_{\mu}\delta x_{\mu}\delta x^{\mu}\right]\nonumber\\ \label{eq : a - ft}\end{aligned}\ ] ] where @xmath34 is given by @xmath113\cup[c\varepsilon,\infty]$ ] and @xmath37 is determined by @xmath114 ( see fig.[fig : cpi ] ) . define @xmath115 and integral out the spatial components , we have @xmath116 taking account of @xmath117 approaches zero when @xmath118 , we can make an approximation of the square root simply as @xmath119 to get accurate enough result on @xmath120 . thus integral over the time component and making use of @xmath121 , we finally get @xmath122 in the last equality we just keep to the first order of @xmath56 . the calculation of the second term are quite similar , @xmath123\exp\left[i\alpha\sum_{\mu}\delta x_{\mu}\delta x^{\mu}\right]\nonumber\\ & = \frac{im_{0}^{2}}{8\pi^{2}\hbar^{2}\varepsilon^{2}}\left[\sum_{\mu}\partial^{\mu}\partial_{\mu}\psi\right]\left(\frac{i\pi}{\alpha}\right)^{\frac{3}{2 } } \int_{\mathbb{c}_{1}}(c\delta t)^{2}e^{-i\alpha ( c\delta t)^{2}}\mathrm{d}(c\delta t ) . \label{eq : a - ft2}\end{aligned}\ ] ] integral term gives @xmath124 regardless the exponential contributions , the first term is in proportion to @xmath125 and the second term in proportion to @xmath126 . at the first approximation , we can neglect the first term . then we have @xmath127\left(\frac{i\pi}{\alpha}\right)^{\frac{3}{2 } } \frac{1}{i2\alpha }    \sqrt{\frac{\pi}{i\alpha}}e^{-i\alpha(c\varepsilon)^{2}/2}\nonumber\\    & = \frac{i}{4\alpha}\,e^{-i\alpha(c\varepsilon)^{2}/2}\sum_{\mu}\partial^{\mu}\partial_{\mu}\psi    \simeq\frac{i\hbar\varepsilon}{2m_{0}}\sum_{\mu}\partial^{\mu}\partial_{\mu}\psi.\end{aligned}\ ] ]    in this way we get @xmath128 expand the left hand term to first order of @xmath56 and comparing with right hand terms , one has one trivial identical equality and @xmath129 for first order of @xmath56 . define @xmath130 , one can rewrite above equation as @xmath131 or , @xmath132 i.e. , @xmath133 with @xmath134 . for space - time component @xmath135 , it satisfies @xmath136 which is exactly the klein - gordon equation . to obtain the non - relativistic approximation , it is natural to set light speed @xmath137 . then we have @xmath138 at classical limit for massive particles when time - arrow ( i.e. , the direction of proper time is identical with time ) is employed . this indicates that the domain @xmath34 in figure [ fig : cpi](a ) becomes a straight line which is given by @xmath139 , and the proper time is dependent and superfluous . the summation lattice in figure [ fig : cpi](b ) also becomes equivalent to the standard path integral of non - relativistic quantum mechanics because @xmath140 .    under this approximation , the action can be expressed as @xmath141 namely , the non - relativistic action between @xmath142 and @xmath143 plus a constant because @xmath144 and @xmath145 . take a free particle as example , the lagrangian is @xmath146 where subscript @xmath147 refers to standard non - relativistic lagrangian . then the action can be written as @xmath148 . therefore from eq.([eq : c - fpi ] ) we have the transition probability amplitude between two quantum states as @xmath149 taking out the constant @xmath150 where @xmath151 we get the transition function as @xmath152 the right side is exactly the standard non - relativistic path integral formalism for free particles.@xcite when potentials due to pure gauge fields are present , the conclusion still holds where instead eq.([eq : action-2 ] ) should be used . 0 d. bohm , _ phys . rev . _ * 85 * , 166 ( 1952 ) ; _ ibid _ 180 . e. nelson , _ phys . rev . _ * 150 * , 1079 ( 1966 ) . l. s. f. olavo , _ physica a _ * 262 * , 197 ( 1999 ) . a. aspect , p. grangier , and g. roger , _ phys . lett . _ * 47 * , 460 ( 1981 ) . a. aspect , p. grangier , and g. roger , _ phys . * 49 * , 91 ( 1982 ) . a. einstein , b. podolsky , and n. rosen , _ phys . _ * 47 * , 777 ( 1935 ) . r. p. feynman , _ rev . phys . _ * 20 * , 367 ( 1948 ) . m. masujima , _ path integral quantization and stochastic quantization _ ( springer , berlin , 2000 ) . d. felea , c. besliu , and i. s. zgura _ phys . j. c _ * 20 * , 713 - 716 ( 2001 ) . v. pazma , czech . b * 29 * , 481 ( 1979 ) . a. kull and r. a. treumann , _ int . phys . _ * 38 * , 1423 ( 1999 ) . p. gosselin and j. polonyi , _ ann . phys . _ * 268 * , 207 - 224 ( 1998 ) . t. padmanabhan , _ foundations of physics _ * 24 * , 1543 ( 1994 ) . j. r. fanchi , _ foundations of physics _ * 35 * , 1453 ( 2005 ) . s. weinberg , _ gravitation and cosmology : principles and applications of the general theory of relativity _ ( john wiley & sons , new york , 1972 ) . j. schwinger , _ phys . * 82 * , 664 ( 1951 ) . this is just based on a naive relationship between classical action and quantum motion equations , which is not always hold when considering variations of wave function with respect to internal symmetry operations and far beyond the scope of this paper . in fact , massive gauge bosons satisfy an equation different from klein - gordon one . however , from the point of view of point particles , it is still quite strange to have a klein - gordon field besides dirac field . j. s. bell , _ physics _ * 1 * , 195 ( 1965 ) .
Abstract:

Answer: a space - time symmetric and explicitly lorentz covariant path integral formalism of relativistic quantum mechanics is proposed , which produces partial locally correlations of quantum processes of massive particles with the velocity of light at low energy limits . a superluminal correlation is also possible if anti - particles that moving along reverse time direction are excited . this provides a new point of view to interpret epr experiments , also leaks a light of hope for hidden variable theories .


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 existence of two - level tunneling defects as a generic property in amorphous systems was postulated four decades ago@xcite in an attempt to explain the remarkable universality in the low energy characteristics of amorphous solids as were found earlier by zeller and pohl@xcite . the exact nature of the tunneling states is not known . yet , their characteristics at low energies has been thoroughly studied through measurements of properties such as specific heat , thermal conductivity , and internal friction ( see refs . @xcite and references therein ) , as well as their relaxation and dephasing times via echo experiments ( see e.g. @xcite ) . at very low temperatures interactions between the tunneling two - level systems ( ttlss ) lead to glassy characteristics and slow relaxation@xcite . recently , in remarkable experiments the coupling between superconducting qubits and ttlss was used to study single ttlss , their relaxation times , and the dependence of their bias energy on applied strain@xcite . many of the above mentioned properties can be explained using what is now referred to as the `` standard tunneling model''(stm)@xcite . in the stm each ttls is characterized by the energy bias between its two states @xmath3 , the tunneling amplitude between the two states @xmath4 , and its coupling to the phonon field denoted by @xmath5 . a central assumption of the stm for the ensemble of ttlss in a given system is that the energy biases and the magnitudes of the barriers are homogenously distributed , resulting in the distribution @xmath6 . the stm further assumes that @xmath5 has a narrow distribution . with regard to the energy biases , the assumption of homogeneity at low energies rests on the argument that the biases are dictated by the large energy related to the disorder , @xmath1ev , and are therefore homogenous at much lower energies . although very successful in explaining many of the universal properties at low temperature , the stm falls short in describing the nature of the tunneling states , the quantitative universality of phonon attenuation , and the energy scale of @xmath2k , only below which universality is observed . as the stm has no energy scale except that of the glass transition , it can not also explain the drastic change of behavior above @xmath2k , and phenomena such as the universal plateau in thermal conductivity at @xmath7k , and the boson peak at higher energies@xcite . thus , the mechanism , generic to the disordered state , that leads to the universality of phonon attenuation at low temperatures , and its relation to the microscopic structure of the amorphous state , remain long standing open questions . the low energy phenomena dictated by the ttlss are observed equally in disordered lattices@xcite . it was carefully shown@xcite that the phenomena are equivalent between the two systems , strongly suggesting that it is the same mechanism leading to the universal phenomena in all disordered systems . it is therefore expected that the central characteristics of the low energy excitations , such as the distribution of their coupling constants and the structure of their density of states ( dos ) , will be similar between the different systems showing universality . it was thus concluded@xcite that the defects in the crystal should be used to model the excitations in amorphous solids , rather than the amorphous structure itself .    in this paper we present numerical results for the density of tunneling states in kbr@xmath8(cn)@xmath9 , where the system is modeled by its bare inter - atomic potentials . kbr@xmath8(cn)@xmath9 is the archetypal disordered lattice showing universality at concentrations @xmath10 , and the magnitude of its energy disorder at the relevant concentration range is @xmath11ev , similar to that of amorphous solids . thus , we could expect a homogenous ttls dos for energies smaller than @xmath12k also in kbr : cn .    using a combined monte carlo ( mc ) and molecular statics technique , we relax fragments of kbr : cn in both 2 dimensions ( 2d ) and 3 dimensions ( 3d ) to a low energy state . we then calculate the energy of all states resulting from a tunneling of a single cn defect . remarkably , a sharp decrease of the dos of the tunneling states is found at energy @xmath2k in both 2d and 3d . this energy is much smaller than the energy scale of the disorder , as is reflected e.g. in the glass temperature . yet , it is the same energy scale that dictates the temperature below which universality is observed . we further find that it is cn flips , which constitute two states related to each other by local inversion symmetry , that overwhelmingly dominate the low energy excitations . at the same time , cn rotations ( asymmetric with respect to local inversion ) have a much larger energy scale ( @xmath12k in 2d , @xmath13k in 3d ) , and a diminishing ( non - homogenous ) dos at low energies .    our results here add up to the bimodality of the coupling strengths of ttlss to the phonon field , as was obtained numerically in @xcite . both these findings are in sharp contradiction to the assumptions of the stm , but in agreement with a recent  two - ttls " theory@xcite for the low temperature universality . our results thus suggest that it is the inhomogeneity of the dos of ttlss at low energies , and its dependence on the symmetry of the ttlss , that dictates quantitative universality of the acoustic properties of glasses , and the energy scale related to it . for a given realization of cn impurities , the low energy state of a kbr@xmath8(cn)@xmath9 sample corresponds to a specific orientation of each of the cn impurities in the sample , which is dictated by the impurity - impurity interactions . at the same time , the exact positions of all ions in the sample giving energy minimization have to be calculated . we therefore use a technique that hybridizes local energy minimization ( conjugate gradients ) and mc simulation . we start by creating a 3d grid of volume @xmath14 ( @xmath15 in 2d , @xmath16 is even ) of @xmath17^+^ and @xmath18 ^ -^ ions , having distance @xmath19 ( @xmath20 in 2d ) between the ions . these distance values were calculated by the energy minimization procedure of pure kbr grid . we replace randomly some of the @xmath18 ^ -^ ions by @xmath21 ^ -^ ions , according to our chosen concentration @xmath22 . the charges of @xmath17^+^ and @xmath18 ^ -^ ions are taken as + 1 and -1 respectively , and the charge of the @xmath21 ^ -^ ion is represented by fractional charges q@xmath23 , q@xmath24 and q@xmath25 placed respectively on the carbon atom , nitrogen atom and at the center of mass @xcite . the distance of the carbon and nitrogen atoms from the center of mass are taken as @xmath26 and @xmath27 respectively @xcite . interatomic potential is calculated by the formula :    @xmath28    the interatomic potential parameters a@xmath29 , @xmath30 and b@xmath29 , see table  [ tab : params ] , are taken from ref . @xcite , and @xmath31kj / mol . .interatomic potential parameters ( taken from ref . cross - interaction parameters were calculated by @xmath32 , @xmath33 , @xmath34 . [ cols="^,^,^,^",options="header " , ]     for a single cn impurity in an otherwise pure kbr lattice we find that its preferred orientation is along the eight - fold degenerate in - space diagonals , in agreement with ref .  [ ] . in two - dimensions the degeneracy is four - fold , along the plane diagonals . these degeneracies are lifted by tunneling . however , in strongly disordered systems , of interest to us here , the bias disorder is much larger than the tunneling amplitudes@xcite . calculation of the bias energy can therefore be carried out without consideration of the tunneling . the tunneling amplitude , which for almost all ttlss is much smaller than the bias energy , can then be considered after the distribution of biases , which is of interest for us here , has been established . we thus start by orienting each @xmath21 ^ -^ ion in one of the inspace diagonals ( plane diagonals in 2d ) randomly . we then relax the system for the given initial orientation of the @xmath21 ^ -^ impurities to the local energy minimum , optimizing all ion locations as well as @xmath21 ^ -^ orientations using the non - linear fletcher - reeves conjugate gradients method . periodic boundary conditions are used to simulate the infinite crystal . the conjugate gradient method , however , does not allow the crossings of the high barriers separating different single cn states . such crossings are required , though , to find a low energy state of the full system , as is dictated by the impurity - impurity interactions . we thus combine the conjugate gradients method with intervening mc steps . after the initial energy minimization , the orientation of one @xmath21 ^ -^ ion is changed randomly to another in space diagonal ( in plane diagonal in 2d ) , and the system s energy is minimized again using conjugate gradient method . the new configuration is accepted or rejected according to the standard metropolis algorithm and the assigned temperature . these steps are carried at 40 different temperatures from @xmath35k to @xmath36k , with 256 mc steps for each temperature , where a single mc step involves a single orientation change of each @xmath21 ^ -^ ion in the sample . we then do 256 mc steps at @xmath37 k , and reach a low energy metastable state . the final state of the system is not necessarily the ground state . however , such metastable states are expected to represent well the real low temperature glassy state of kbr : cn and its low energy dos . the latter is dictated by the efros shklovskii gap@xcite , which relies only on stability to single and double @xmath21 ^ -^ flips or rotations@xcite . once we reach the final state of the simulation we measure all single particle excitation energies . we also keep track of the symmetry of the excitations . the symmetric excitations , i.e. those where the final orientation of the @xmath21 ^ -^ ions is nearly inverse to its initial orientation ( within @xmath38 ) are denoted @xmath39-ttlss , and the asymmetric excitations are denoted @xmath40-ttlss@xcite . this procedure is then repeated for 3000 samples in 3d and 5000 samples in 2d . in fig.[fig : dos3d](a ) and fig.[fig : dos2d](a ) we plot the full dos of single impurity tunneling in 3d ( @xmath41 unit cells ) and 2d ( @xmath42 unit cells ) respectively . the dos has a maximum at low energies , but reduces considerably at an energy @xmath2k . we now plot separately the single particle excitations of cn flips ( @xmath39-ttlss ) and rotations ( @xmath40-ttlss ) in fig.[fig : dos3d](b),(c ) for the 3d samples and in fig.[fig : dos2d](b),(c ) for the 2d samples . the difference between the dos of the two types of excitations is striking . the inversion symmetric flip excitations are peaked at low energies , with typical energies of @xmath43k for both 3d and 2d . the inversion non symmetric rotation excitations , however , are broadly distributed , with typical energies of @xmath13k for 3d and @xmath12k for 2d . the peak value of the @xmath39-ttls dos is @xmath44 and @xmath45 times larger than the peak value of the @xmath40-ttls dos in 3d and 2d respectively . the difference between the peak values of the @xmath40-ttlss in 3d and in 2d is a consequence of each cn impurity having two @xmath40 excitations in 2d , but six @xmath40 excitations in 3d . focusing on the cn flips ( @xmath39-ttlss ) , their dos is quite close to a gaussian , except for a small dip in the dos at very low energies [ see inset of figs .  [ fig : dos3d](b ) ,  [ fig : dos2d](b ) ] , and some structure at @xmath46 k. the former is generic to the two - ttls model ( see below ) whereas the latter is a result of details of the lattice structure and form of the interaction . with regard to the @xmath40-ttlss , their dos is well fit to a log function at low temperatures , except below @xmath47 k , the energy at which the @xmath39-ttlss appear , see further discussion below .      our results for the dos are intimately connected with the fact that the coupling strength to the phonon field of inversion symmetric excitations @xmath48 is much smaller than that of the asymmetric excitations @xmath49 . for kbr : cn @xmath50 , whereas @xmath51@xcite . this bi - modality in the coupling strengths affects the phonon mediated acoustic interactions , which then result in an effective ttls - ttls interaction hamiltonian of the form@xcite    @xmath52    where in three dimensions ( two dimensions ) all interactions decay as @xmath53 ( @xmath54 ) at distances @xmath55 , with a short distance cutoff @xmath56 , and their typical values at near neighbor distance are related by@xcite    @xmath57    here @xmath58@xcite . this effective hamiltonian was shown to lead to quantitative universality and to account for the energy scale of @xmath59k related to it@xcite . the dos for the two types of ttlss within the hamiltonian  [ effectivestau ] was analyzed in detail in ref.@xcite . indeed , it was found that the weakly interacting ( @xmath39 ) ttlss dos are well fit by a gaussian , except for a small dip , at very low energies , of relative magnitude @xmath60 . the gaussian width is dictated by the strength of the @xmath61 interaction @xmath62 . the strongly interacting @xmath40-ttlss were found to have a typical energy dictated by @xmath63 , a logarithmic gap at intermediate energies@xcite , and a power - law gap below @xmath64 , the energy where the dos of the @xmath39-ttlss becomes appreciable . as a consequence of this gapping of the @xmath40-ttls dos @xmath65 , one finds below @xmath2 k that @xmath66 ( i.e. @xmath67 ) . this defines the condition for the @xmath39-ttls to dominate phonon attenuation , resulting in universality below this energy scale . our results here validate the applicability of the above two - ttls model to real systems , as we verify , in a calculation relying only on the bare interatomic interactions of a real glass , the central features predicted by the model : ( i ) the width of the distributions ( typical energies ) of the @xmath40-ttlss and @xmath39-ttlss - both their relative magnitudes and the absolute magnitude of each one . we further note that the value of @xmath68 inferred from the ratio between the typical energies is in good agreement with that obtained from the ratio of the interaction constants in ref.@xcite ; ( ii ) the gapping of the @xmath40-ttlss but not the @xmath39-ttlss at low energies , and the small dip , of relative magnitude @xmath60 , of the @xmath39-ttlss at the lowest energies : ( iii ) the change in the functional form of the gap of the @xmath40-ttlss at the energy scale @xmath46 k where the @xmath39-ttlss appear : ( iv ) the relation @xmath67 is fulfilled at energies smaller than @xmath69 k. we note that the energy below which the above condition is fulfilled in kbr : cn is slightly larger than what is inferred by our results , because as a consequence of finite size effect our results over estimate the @xmath40-ttls dos at the lowest energies . in contrast to the assumption of the stm , we find that in kbr : cn ttlss are divided to two distinct classes by their local inversion symmetry , with a dos distribution which is particular to each class of ttlss , and has a strong energy dependence at low energies , of order @xmath47 k. many disordered lattices , of which kbr : cn is a primary example , share the low temperature universal characteristics with amorphous solids . furthermore , experiments strongly support the notion that it is the same mechanism , pertaining to the disordered state of matter itself , that dictates universality in amorphous systems and disordered lattices alike . the structure we find for the dos of the symmetric and asymmetric ttlss is not only very different from the flat dos suggested by the stm , but it is this very structure of dos that provides an explanation of the quantitative universality of phonon attenuation and the energy scale of @xmath2k below which universality is observed@xcite . as it is expected that the same model would be relevant to both amorphous solids and disordered lattices showing universality , our results suggest that also in amorphous solids two types of ttlss exist , with a similar structure of dos to the one found here . verification of this idea could lead not only to a resolution of the long standing problem of the low temperature universality in disordered systems , but also to an enhanced understanding of the microscopic structure of amorphous solids , and its relation to the physical properties of amorphous solids in general . the standard tunneling model uses only one feature of the complicated structure of the density of states found here - the rather homogeneous dos of the inversion symmetric ttlss below @xmath70k . it completely neglects not only the sharp reduction of the dos at higher energies , but also the existence of the asymmetric ttlss with their much larger interaction with the phonon field . these characteristics of the dos , however , may lead to the explanation of further intriguing properties of disordered lattices and amorphous solids at low and intermediate temperatures . some examples at intermediate temperatures are the full temperature dependence of the thermal conductivity including the rather universal plateau between @xmath7k@xcite , and the boson peak . at low temperatures , the existence of the strongly interacting @xmath40-ttlss , although small in number , may be useful in explaining equilibrium and non - equilibrium properties of phonon attenuation , as well as the relaxation and decoherence of the prevalent @xmath39-ttlss . in addition to the standard ttlss responsible to the universal phonon attenuation at low temperatures , a second type of local excitations was introduced previously , phenomenologically , by black and halperin@xcite and by yu and freeman@xcite . black and halperin suggested the existence of ttlss with a much smaller interaction with the strain , thus affecting the specific heat but not phonon attenuation . yu and freeman suggested the existence of local phononic modes with a similar interaction with the phonon field to that of the ttlss , but with a dos which is gapped below @xmath71 k , to fit experimental data for thermal conductivity at intermediate and high temperatures . our results differ from the above two approaches in that ttlss of the second type have a much larger interaction with the strain , and are softly gapped below @xmath69 k.    experimentally , ttlss are observed via a plethora of techniques , including echo@xcite , coupling to superconducting qubits@xcite , specific heat@xcite and dielectric response@xcite . in general , our results suggest that care should be taken in experiments to characterize the ttlss according to their interaction strength . whereas most experiments are susceptible to the weakly interacting ttlss abundant at low energies , strongly interacting ttlss can dominate experimental observations susceptible to their high bias energy or strong interaction with the strain . specifically , our results disagree with previous estimates for typical energy biases of cn flips in kbr : cn@xcite , and suggest that experimental values of @xmath72 k obtained for ttls - ttls interaction at short distances , and for ttls bias energies , in kbr : cn@xcite concern cn rotations rather than cn flips .    the use of the empirical potential of the ( buckingham - coulomb ) form given in eq . ( [ potential ] ) rather than the more accurate dft or ab - initio formalisms is dictated by the required complexity of the problem : the need to consider a relatively large system , to find its low energy state using monte carlo , and to repeat the calculation for many samples for good statistics . this use of the empirical potential has been shown to be useful in obtaining qualitative , and to some degree also quantitative , predictions for ionic crystals ( for details see refs . @xcite ) . for the properties which are of interest to us here , i.e. the dos of the @xmath21 ^ -^ flips ( @xmath39-tlss ) and @xmath21 ^ -^ rotations ( @xmath40-tlss ) , all qualitative features are dictated by symmetry and are therefore robust . comparison to results for the model hamiltonian in eq.([effectivestau])@xcite suggest that all qualitative features of the dos for both @xmath21 ^ -^ flips and @xmath21 ^ -^ rotations are indeed retrieved by our calculation . quantitatively , the value of @xmath68 inferred from our results for the typical energies of @xmath21 ^ -^ flips and @xmath21 ^ -^ rotations is in very good agreement with that found by the ratio of their corresponding couplings to the phonon field using dft and ab - initio calculations@xcite . we also find good quantitative correspondence with the calculations for the model hamiltonian in eq.([effectivestau ] ) as obtained in ref . this suggests that our results provide also reasonable quantitative accuracy , except for the over estimate of the dos of the @xmath40-tlss at the very low energies , as mentioned above . we are concerned here with the calculation of the bias energies of the @xmath40-tlss and the @xmath39-tlss , but not their tunneling amplitudes . the latter are crucial for the discussion of dynamic and non - equilibrium properties . a plausible assumption for the distribution of the tunneling amplitudes @xmath73 is that suggested in refs.@xcite , @xmath74 . for the @xmath39-tlss which dominate the low energy universal phenomena , this assumption has been shown to be consistent with the various experimental findings . however , its applicability for the @xmath40-tlss needs to be checked . such a verification may become possible by using our results for the bias energies here , and for the coupling constant of the @xmath40-tlss in ref.@xcite to calculate measurable quantities dominated by @xmath40-tlss and comparing to experimental results . _ acknowledgments _  we thank a. burin , a. gaita - ario , niels grnbech - jensen , and joerg rottler for very useful discussions . this research was supported by the israel science foundation ( grant no .
Abstract:

Answer: amorphous solids , and many disordered lattices , exhibit a remarkable qualitative and quantitative universality in their acoustic properties at temperature @xmath0k . this phenomenon is attributed to the existence of tunneling two - level systems ( ttlss ) , characterized by a homogenous density of states ( dos ) at energies much lower than the disorder energy ( @xmath1ev ) . here we calculate numerically , from first principles , the dos of kbr : cn glass , the archetypal disordered lattice showing universality . in contrast to the standard tunneling model , we find that the dos diminishes abruptly at @xmath2k , and that tunneling states differ essentially by their symmetry under inversion . this structure of the ttlss dictates the low temperature below which universality is observed , and the quantitative universality of the acoustic properties in glasses . consequences to the properties of glasses at intermediate temperatures , as well as to the microscopic structure of amorphous solids , are discussed .