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 full - energy , @xmath4200 gev , au@xmath1au collisions at the relativistic heavy - ion collider ( rhic ) , pulses of tens to hundreds of photons from @xmath5 decays stream outwards from the primary collision vertex at the speed of light and very close to the beam pipe . a detector capable of measuring the arrival times of these photons can thus provide important information for event triggering , especially on the location of the primary vertex along the beam pipe , and on the time that the event occurred relative to free - running master clocks used by other fast timing detectors in the experiment . in the solenoidal tracker at rhic ( star ) @xcite , the first detector capable of such measurements that was implemented was the  pseudo vertex position detector " ( pvpd ) @xcite . this system consisted of two identical assemblies of three readout detectors , one on each side of star ( east and west ) , with each assembly mounted immediately outside the beam pipe . the pvpd worked well in full - energy au@xmath1au collisions @xcite . however , for lighter beams such as p@xmath1p and also au@xmath1au beams at lower beam energies , the efficiency of the pvpd for providing the vertex location and the event time degraded due to the relatively lower multiplicities of very forward prompt particles in such collisions . to address this , the pvpd was upgraded to increase the readout detector channel count on each side of star in the same acceptance . this detector is called the vertex position detector ( vpd ) .    unlike the pvpd , the vpd has been fully integrated into the star trigger system . the vpd provides the primary detector input to the star minimum bias trigger in a@xmath1a collisions . each vpd assembly measures up to nineteen times in each event . these times are thus available , both online and offline , to measure the location of the primary vertex along the beam pipe , @xmath6 , via the equation , @xmath7 where @xmath8 and @xmath9 are the times from each of the two vpd assemblies and @xmath10 is the speed of light . the event  start time , " which is needed by the star time - of - flight ( tof ) and muon telescope detector ( mtd ) @xcite systems to perform particle identification at mid - rapidity , is given by , @xmath11 where @xmath12 is the distance from either assembly to the center of star . one motivation for the increased channel count from the pvpd to the vpd in the same angular acceptance was to increase the efficiency of the detector for recording hits in single events . another important motivation involves the fact that the experimental resolution on either @xmath8 or @xmath9 used in equations ( [ eq : zvtx ] ) and ( [ eq : tstart ] ) improves like 1/@xmath13 , where @xmath14 is the number of channels on each side that were hit by prompt particles . that is , taking @xmath8 or @xmath9 as the average over all channels hit by prompt particles in an event results in a resolution on , _ e.g. _ , @xmath6 , given by , @xmath15 where @xmath16 is @xmath8 or @xmath9 , @xmath17 is the resolution of the time difference @xmath8@xmath18@xmath9 , @xmath19 is the resolution on @xmath8 or @xmath9 , and @xmath20 is the time resolution of a single readout detector . the larger channel count in the vpd thus allows an averaging over a larger number of lit channels and hence an improved performance for measuring both @xmath6 and @xmath21 . the larger number of readout channels also provides more time values in the same acceptance , allowing one to better recognize which times are prompt and which times should be rejected as ( delayed ) outliers . the rates for such non - prompt hits are nearly negligible in full - energy au@xmath1au collisions , but are significant for p@xmath1p and lower energy au@xmath1au collisions . this paper is organized as follows . section [ sec : design ] describes the design of the detectors and the electronics . section [ sec : calibandperform ] describes the configuration of the system for a new rhic beam , the offline calibration of the detectors , and the start - timing and @xmath6 performance . finally , section [ sec : summary ] presents the summary and conclusions . the vpd exists as two identical detector assemblies , one on the east and one on the west of star . the design of these assemblies is described in section [ sec : dets ] . each of the nineteen detectors used in each assembly is composed of a pb converter followed by a fast , plastic scintillator which is read out by a photomultiplier tube ( pmt ) . the signals from the nineteen detectors in each assembly are digitized independently by two different sets of electronics which are described in section [ sec : electronics ] . the detector was first installed in advance of the 2007 rhic run . the bases for the pmts were revised before the 2008 rhic run . the detector has been used without additional modifications in every rhic run since then . each vpd assembly consists of nineteen detectors , a side view of which is shown in figure [ fig : det ] . each detector housing is a 2 inch outer diameter and 0.049 inch thick aluminum cylinder with 3/8 inch thick aluminum front and back caps . inside this cylinder is a 0.25 inch non - conducting spacer , then the active elements consisting of a 0.25 inch ( 1.13 radiation lengths ) pb converter , and a 1 cm thick scintillator ( eljen ej-204 ) coupled to a 1.5 inch diameter hamamatsu r-5946 mesh dynode pmt via rtv-615 optically transparent silicone adhesive . the pmts used in the vpd were taken from the tofp detector @xcite after it was decommissioned in 2005 . the pmt dynode voltages are provided by a conventional linear resistive base . an initial version of these bases was used in the 2007 rhic run , and the detectors performed well , but there was some  ringing " in the trailing edge of the pmt line shapes . this slightly complicated the offline slewing corrections , but did not affect the timing performance . in 2008 , these bases were revised . no changes were made to the component values or connections , but the placement of these components was revised to minimize the lengths of the traces inside each pmt base circuit board . the excessive trace lengths in the previous version of the bases resulted in enough distributed inductance to cause the trailing edge ringing , which was completely removed in the final bases . a wire connects the pmt cathode pin to a 0.001 inch thick aluminum cylinder which extends past the active elements . another wire is used to connect the aluminum cylinder to the lead converter . in this way , the active elements are enclosed on all but one side by an electrostatic shield . this shield is electrically isolated from the active elements inside and the aluminum outer cylinder outside by several layers of kapton tape . the output coaxial connector shield is isolated from the detector housing and the high voltage ground but is indirectly connected via a 1 k@xmath22 resistor . this prevents the ( inductive ) shield of the coaxial signal cable from forming an undesirable resonant circuit with the ( capacitive ) electrostatic shield and detector housing while maintaining the high voltage ground return path . each vpd assembly consists of two rings of readout detectors and is mounted to the i - beam that supports the star beam pipe . a front view of one of the vpd assemblies is shown in figure [ fig : assy ] . the outer diameter of the beam pipe at this distance is five inches . an assembly exists as two semi - annular  clam - shells " that enclose the beam pipe . these are bolted together and are held in place by delrin support blocks which attach to a horizontal mount plate which is clamped to the beam pipe support i - beam . the beam pipe and i - beam are at a different ( dirty ) electrical ground than the experiment , so the delrin support blocks both hold the assembly in place and electrically isolate it . the two assemblies are mounted symmetrically with respect to the center of star at a distance of 5.7 m. the nineteen detectors in each assembly subtend approximately half of the solid angle in the pseudo - rapidity range of 4.24@xmath23@xmath24@xmath235.1 . when viewed from the rear and looking towards the center of star , the detectors are numbered 1 - 10(11 - 19 ) counter - clockwise starting from the lower right in figure [ fig : assy ] in the inner(outer ) ring . the signals from the vpd detectors are digitized by two different sets of electronics . a schematic view of the components that are involved is shown in figure [ fig : electronics ] . the signals from the vpd detectors are connected via low - attenuation @xcite belden 9310 rg-58  short coaxial " cables to  tpmd " boards which are mounted inside an aluminum  start box " that is positioned within @xmath26 feet of the vpd detector assemblies . also mounted inside this box are the  tdig " and  tcpu " boards . digital data are sent over category-6 cables to  thub " boards , which then transmit the data over optical fiber to receivers in the star data acquisition system @xcite . additional details on these electronics boards can be found in refs . @xcite and @xcite . the tpmd boards split the detector signals after they pass a diode - based clamping circuit which limits the pulse height to 4v . one output undergoes a leading - edge discrimination of these signals with respect to a tunable threshold at typically 25 mv . the maxim max9601 comparator outputs are sent to the tdig boards which use the cern hptdc chip @xcite to digitize the arrival times of the signals with respect to an externally - input free - running clock . each tdig board uses three hptdc chips in  very - high resolution " mode . the digitized times are 21 bit data words with a dynamic range of 52 @xmath25s and a least significant bit ( lsb ) time conversion of 25ns/1024 @xmath26 24.4 ps . the pulse sizes are measured by the tdig boards as the width of the pulse at the threshold , which is called the  time over threshold " ( tot ) . all 2@xmath019 channels of the vpd are digitized in these electronics . the hptdc chips exhibit a timing cross - talk @xcite of up to three lsb when two signals arrive at the same hptdc chip within a few nanoseconds of each other . to avoid this , there are five tpmd@xmath1tdig board pairs in each start box , and each board pair handles only 3 - 4 detector channels out of the 24 possible . when two detector channels are connected to the same hptdc chip , the  short coax " input cables are staggered in length in 12  ns steps so as to remove the possibility of hptdc timing cross - talk .    the other signal from the passive split performed in the tpmd boards is passed on with negligible attenuation and distortion since the stub feeding the high impedance comparator is very short . this signal is sent over long ( @xmath2100 ft ) coaxial cables to electronics in the star trigger system called  qt boards " @xcite . short sections of coaxial cable are used after the tpmd boards and before these  long coax " cables to equalize the signal arrival times at the qt electronics to a few nanoseconds given the presence of the short coax cable staggering before the tpmd boards used to avoid the hptdc cross - talk . the qt boards are 9u vme boards and are part of the star trigger system @xcite . they perform pulse area and time measurements using an analog - to - digital ( adc ) and time - to - amplitude - conversion ( tac ) circuitry for all of the star trigger detectors . the adc and tac measurements are each 12-bit numbers . the tac time measurement is common - stop with respect to the 9.4 mhz rhic clock with a digital - to - time conversion of @xmath218 ps / lsb . for system design simplicity , only 2@xmath016 channels of the vpd are digitized in qt boards . the information needed from the vpd for use in equations ( [ eq : zvtx ] ) and ( [ eq : tstart ] ) are time values with the best possible experimental resolution . this first requires that the high voltage values used for each of the vpd pmts result in equalized gains , _ i.e. _  the pulses from the different detectors should have the same size for the same amount of light generation in the scintillators . the time values digitized by the electronics described in section [ sec : electronics ] include offsets resulting from propagation times inside the pmts , signal cabling , and trace lengths in the electronics . there is also a dependence of the digitized pulse times on the size of the detector pulses which is called  slewing " @xcite . the contributions from these sources of smearing are removed via offline calibrations of the vpd data . these calibrations , and the gain - matching of the vpd pmts , is described in this section . the timing performance of the vpd detectors following these calibrations will also be presented in this section . rhic is an extremely flexible machine . in the six runs ( one per year ) since the vpd was installed , it has provided p@xmath1p collisions at beam energies , , of 200 and 510 gev , as well as cu@xmath1au , au@xmath1au , and u@xmath1u collisions at a number of beam energies , , ranging from 7.7 to 200 gev . the species and the number of particles per event producing signals in the vpd detectors depend strongly on the rhic entrance channel and the collision centrality . in au@xmath1au collisions at 200 gev , the signal in a vpd detector in one event is produced by 10 photons from @xmath5 decays and an approximately equal number of charged pions . in the lowest energy au@xmath1au collisions , _ i.e. _  7.7 gev , a vpd hit results from one or two spectator protons . the beam energy at which these trends cross depends on the centrality but is near 30 - 50 gev according to urqmd @xcite simulations . in similarity to the full - energy au@xmath1au collisions , vpd hits in p@xmath1p collisions result from a mixture of photons from @xmath5 s and charged pions , but unlike the full - energy au@xmath1au collisions , a vpd hit results from generally only one such particle per event . in order to provide consistently useful information from the vpd in these very different beam environments , the gains of the vpd pmts must be set for each rhic beam separately .    when a new beam is made available to the rhic experiments , star sets up a simple interaction trigger generally based on the zero degree calorimeter @xcite . using this trigger , data are collected from the vpd using three different sets of high voltage values , one at best guess values and the other two at @xmath27100 v with respect to the best guess values . the average adc values from the trg qt boards , @xmath28 , and the average tot values from the tdig electronics , @xmath29 , are measured for each of the three gain sets . the dependence of both quantities on the pmt high voltage value is fit with a power law for each channel . for the 2@xmath016 vpd channels digitized by the trigger system qt boards , these functions are interpolated to determine the high voltage resulting in @xmath28@xmath30300 . the relationship between the best high voltage values obtained from the @xmath28 and @xmath29 values is linear , which allows the gains for the 2@xmath03 channels not digitized in the qt boards to be consistently set on the basis of the @xmath29 values . the high voltage set - points so calculated are @xmath2400 v lower in 200 gev au@xmath1au collisions than they are in 200 gev p@xmath1p collisions . the slopes of the power - law gain curves vary by 6 - 8% for the different detectors and are @xmath25.8(7 ) in 200 gev p@xmath1p(au@xmath1au ) collisions .    once the gains are set , data are collected to check the final gains and to set the basic timing offsets in the qt boards . these offsets are determined by plotting the time in a vpd channel minus the average over the times in all other lit channels on the same side of star in the same event . the average values of these  1-@xmath31 " distributions are then set as the tac offset values in the qt boards .    at this point , the vpd is commissioned for the new rhic beam . minimum bias triggers in star commonly use the vpd as the primary detector input . these triggers generally include constraints on the vpd time difference , _ i.e. _  equation ( [ eq : zvtx ] ) , in order to enhance the rate for collisions near the center of star which have the best particle measurement efficiencies and the lowest backgrounds from particles produced in various detector support materials . these @xmath6 constraints are defined via upper and lower limits on the difference between the uncalibrated times from the earliest vpd hit on each side of star . such an approach results in a @xmath6 resolution of @xmath25.5 cm in full - energy au@xmath1au collisions . the integrated non - linearity in the hptdc chips is corrected via code - density tests @xcite . the smearing of the detector time values resulting from channel - dependent propagation delays , _ _ offsets , and from the slewing are removed via offline calibrations . the procedure is iterative , and results in a table of values binned by the pulse size , tot , for each channel . in subsequent analyses , the event - by - event value of the pulse size in a lit detector in an event is used to look up a value from these tables which is subtracted from the raw value to remove the offsets and slewing . these calibration tables are determined as follows . the events used are from the minimum bias trigger , often based on the vpd data itself , which requires that at least one vpd channel on each side of star is lit in an event . the location of the collision vertex along the beam pipe as reconstructed from the primary tracks in the star time projection chamber ( tpc ) @xcite , @xmath32 , was required to be within @xmath2750 cm of the center of star . the transverse radius of the primary collision vertex was required to be less than 2 cm ( inside the 2.5 cm radius beam pipe ) . at least two primary tracks were also required . the first step involves the determination of initial relative offsets . these are obtained as the average values of the times in specific lit channels with respect to a reference channel ( channel 1 on the west ) . once these offsets are obtained , they are applied to the raw time values . the resulting time values from a vpd channel minus the average of the times of all other lit channels on the same side of star in the same event are plotted as a function of that channel s tot value . the dependence of the average  1-@xmath33@xmath14@xmath34 " times on each channel s tot value is used in subsequent passes through the data to reduce the slewing and offsets . the procedure is iterative and determines three such calibration curves for each channel . the two vpd assemblies , one east and one west , are treated in parallel . the range of the dependence of the 1-@xmath33@xmath14@xmath34 time differences on the tot values is generally a few ns before the first pass , a few hundred ps before the second pass , and is negligible after the third pass . at the end of the procedure , the initial offsets and the three slewing correction tables for each channel are added together to produce the final correction tables versus the tot values .    in full - energy au@xmath1au collisions , every vpd channel is hit by 10 s of prompt particles in every event . as the beam energy of au@xmath1au collisions is decreased , the number of channels lit by prompt particles strongly decreases , and the probability that vpd channels are lit by non - prompt particles increases . these  late hit " channels must be rejected . this outlier rejection is performed as follows . the location of the primary vertex along the beam pipe calculated from the vpd times , @xmath35 , is calculated using equation ( [ eq : zvtx ] ) from every possible pair of lit channels with one channel on the east and one channel on the west . if the @xmath35 value from a pair of channels is consistent with the value of @xmath32 , these two channels are flagged as having prompt hits . this consistency is with respect to an analysis cut that becomes tighter in each calibration pass . it is also required that each of the @xmath14 times used in the 1-@xmath33@xmath14@xmath34 relative time calculation are consistent with the average values , @xmath33@xmath14@xmath34 , with a @xmath24@xmath36 cut that also becomes tighter in each calibration pass . " time differences versus @xmath14 for a typical vpd detector channel in the 7.7 gev au@xmath1au data . the functional fit ( see text ) allowing the extraction of the single detector resolution , @xmath20 , is shown as the solid line . in the left frame , the full outlier rejection is performed , and in the right frame , no outlier rejection is performed.,scaledwidth=70.0% ]    following the calibration procedure , the time resolution of a vpd channel after this procedure is called the single detector resolution , @xmath20 , and is extracted as follows . two dimensional plots of the  1-@xmath33@xmath14@xmath34 " time differences versus @xmath14 are filled , and each bin of the @xmath14 x - axis is fit with a gaussian . the dependence of the standard deviations and mean values from these gaussian fits as a function of @xmath14 for a typical vpd channel is shown in the left frame of figure [ fig : sigma0_meansigm ] . the solid line in this frame is the fit to the standard deviations using the functional form @xmath36(1-@xmath33@xmath14@xmath37@xmath30@xmath38 , where the single detector resolution , @xmath20 , is the fit parameter . important requirements of a successful vpd calibration are that the mean values of the  1-@xmath33@xmath14@xmath34 " time differences ( open points ) are near zero for all values of @xmath14 , and that the standard deviations ( solid points ) follow the expected @xmath38 trend shown as the solid line in figure [ fig : sigma0_meansigm ] . these requirements indicate the extent to which the vpd times measured by the different detector channels are self - consistent and result only from prompt particles . the right frame of figure [ fig : sigma0_meansigm ] shows the same standard deviations and mean values from the gaussian fits to the 1-@xmath33@xmath14@xmath34 time differences as a function of @xmath14 but when no outlier rejection is performed . in this frame , the standard deviations ( solid points ) are @xmath2100  ps larger and no longer follow the expected @xmath38 trend , and the mean values ( open points ) deviate significantly from zero . this underscores the importance of careful vpd time outlier rejection especially in p@xmath1p and low energy au@xmath1au collisions .          the single detector resolution , @xmath20 , values extracted as described in the previous subsection are plotted versus the vpd channel number in figure [ fig : sigma0 ] for minimum bias data from a number of different rhic entrance channels . the average values of the single detector resolution is 94 ps for full - energy au@xmath1au collisions , and increases to @xmath2150 ps for intermediate energy au@xmath1au collisions and 510 gev p@xmath1p collisions . the resolution for 7.7 gev au@xmath1au collisions is approximately 200 ps , and is not shown as the vpd is generally not used to provide the start time to the tof system due to its relatively low efficiency per event for those beams . the increasing single detector resolution from full - energy au@xmath1au to the lower energy au@xmath1au and p@xmath1p collisions reflects the fact that the vpd is doing multiple - particle timing ( per detector channel ) in the highest energy au@xmath1au collisions and there is a gradual evolution to single - particle timing in p@xmath1p and lower energy au@xmath1au collisions .    the resolution by which the vpd measures the start time needed by the tof and mtd detectors , _ i.e. _  equation ( [ eq : tstart ] ) , goes like @xmath2@xmath20/@xmath39 where @xmath40 is the total number of vpd channels lit by prompt particles in an event . in full - energy au@xmath1au collisions , the start time resolution is observed to be 20 - 30 ps , which is essentially negligible compared to the stop time resolution of the tof(mtd ) detectors of @xmath280(100 ) ps @xcite . in p@xmath1p collisions , the vpd start time resolution is approximately 80 ps , as for those beams @xmath20 is @xmath2150 ps and @xmath40@xmath263 .    , versus the same position as obtained using the primary tracks reconstructed in the tpc , @xmath32 in 510 gev p@xmath1p collisions ( left frame ) and 200 gev au@xmath1au collisions ( right frame ) . the insets depict the difference @xmath41 @xmath30 @xmath35-@xmath32 which allows the extraction of the vpd s @xmath6 resolution as the standard deviation of the difference distributions.,scaledwidth=90.0% ]    the resolution by which the vpd measures the primary vertex location is determined by comparing the vpd s measurement , @xmath35 ( _ i.e. _  equation ( [ eq : zvtx ] ) , to that obtained from the primary tracks reconstructed in the tpc , @xmath32 . these two quantities are plotted in figure [ fig : zz ] for 510 gev p@xmath1p collisions ( left frame ) and 200 gev au@xmath1au collisions ( right frame ) . the insets in each frame depict the difference @xmath41 @xmath30 @xmath35 @xmath18 @xmath32 and gaussian fits to determine the vertex resolution . the standard deviations of these distributions obtained from the fits are typically @xmath22.4 cm and @xmath21 cm in 510 gev p@xmath1p collisions and 200 gev au@xmath1au collisions , respectively . the 2@xmath03 channel  pvpd " @xcite vertex and start - timing detector in the star experiment at rhic has been replaced by a 2@xmath019 channel detector in the same acceptance . this vertex position detector ( vpd ) exists as two identical assemblies , one on each side of star , very close to the beam pipe and @xmath25.7 m from the center of star . the readout channels in each assembly include a pb converter followed by a fast plastic scintillator and a mesh dynode pmt . the pmt signals are digitized by two different sets of electronics for use in the star level-0 trigger to select minimum bias collisions , to constrain the location of the primary collision vertex along the beam pipe , and to provide the start time needed by other fast timing detectors in star . the system must be configured for each rhic beam separately to provide a consistent performance despite the wide range of beam particles ( protons to au ) and beam energies ( 7.7 to 510 gev ) provided by the rhic . the slewing and offset corrections are performed using an iterative procedure and require a careful rejection of outlier times from non - prompt particles . the single - detector resolution of the vpd , @xmath20 , is approximately 95  ps in full - energy au@xmath1au collisions , and degrades to @xmath2150  ps in p@xmath1p and lower energy au@xmath1au collisions . the start time resolution of the vpd ranges from 20 - 30 ps in full - energy au@xmath1au collisions to @xmath280 ps in p@xmath1p collisions . the resolution by which the vpd can measure the location of the primary vertex in full - energy p@xmath1p and au@xmath1au collisions is @xmath22.5 cm and @xmath21 cm , respectively . we thank the star collaboration for the use of the experimental data shown in this paper and the operation of this system during rhic running periods as part of star standard shift crew operations . we thank allan schroeder and the members of the ut - austin machine shop for the machining of the structural parts of the detector assemblies . we appreciate the expert assistance of the bnl collider - accelerator department technicians charlie bloxson , matt ceglia , and robbie karl . we gratefully acknowledge funding from the us department of energy under grant numbers de - fg02 - 10er41666 and de - fg02 - 94er40845 . k. h. ackermann _ et al . _ , . w.j . llope _ et al . _ , . w.j . llope _ a. ljubii jr . _ et al . _ , ; + j.m . et al . _ , . j.j . schambach _ et al . _ , e * 16 * , 2496 ( 2007 ) . et al . _ , . et al . _ , . h.j . et al . _ , http://hena.lbl.gov/fms/cap-equip-qt.pdf f.s . et al . _ , . b. bengtson and m. moszyski , ; + b. bengtson and m. moszyski , ; + t. sugitate _ et al . _ , . s.a . bass _ et al . _ , . m. anderson , _ et al .
Abstract:

Answer: the 2@xmath03 channel pseudo vertex position detector ( pvpd ) in the star experiment at rhic has been upgraded to a 2@xmath019 channel detector in the same acceptance , called the vertex position detector ( vpd ) . this detector is fully integrated into the star trigger system and provides the primary input to the minimum - bias trigger in au@xmath1au collisions . the information from the detector is used both in the star level-0 trigger and offline to measure the location of the primary collision vertex along the beam pipe and the event  start time " needed by other fast - timing detectors in star . the offline timing resolution of single detector channels in full - energy au@xmath1au collisions is @xmath2100 ps , resulting in a start time resolution of a few tens of picoseconds and a resolution on the primary vertex location of @xmath21 cm . # 1#2#3#4#1 * # 2 * , # 3 ( # 4 )    vertex position detector , time resolution 25.75.cj , 29.40.cs    @xmath3 corresponding author . _ e - mail address : _ llope@rice.edu .


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: recent years have witnessed several initiatives on enabling internet access to the next three billion people . it is widely recognized by all the stakeholders involved in the game of universal service provisioning , that providing access to the internet has no one size fits all solution for enabling wider universal internet access but requires exploring a variety of solutions . this is evident from organizations like facebook and google who are on a ( noble ) mission to connect the next three billion through novel ways utilizing high altitude platforms such as drones , balloons and satellites . access to the internet necessarily translates to access to its services . this means that the goal of providing internet access requires access to its critical service infrastructure , which are currently hosted in the cloud . however , recent works have pointed out that the current cloud centric nature of the internet is a fundamental barrier for internet access in rural / remote areas as well as in developing regions . taking africa as an example , we recently conducted a study mapping the african web ecosystem @xcite . our study highlighted that majority of the african web infrastructure is placed outside the continent . this has direct impact on aspects of dns resolution times as well as http page load times causing users to experience higher latencies for downloading the webpages and poor quality of experience . our paper concluded that there is a need for localised service infrastructures . access to services are also extremely important ( even potentially life saving ) during humanitarian crisis where access to services are challenged - due to intermittent connectivity , heavy interference and congestion leading to increased latencies to global servers or in most cases no access to them . hence its not rocket science to understand that access to localized service infrastructures is of paramount importance to solve the global access problem . recent work on mobile edge computing @xcite aims to push computation right at the edge of mobile networks , enabling computations at the edge improving latencies and performance . the recent work on infrastructure mobility @xcite illustrates the interesting concept of making the access infrastructure mobile thus providing better and much more efficient coverage based on the need / demand from the users . so an obvious question or follow up to this exciting idea is rather than making the access infrastructure mobile , why can not build on such work to make the service infrastructure mobile i.e. rather than expecting the users to access the cloud services , why can not we directly provision the cloud service infrastructure to the user on demand even in the absence of a terrestrial infrastructure e.g. rural / remote areas or disaster zones ?    in this paper , we present @xmath0- a preliminary idea of deploying a lightweight micro cloud infrastructure in the sky using indigenously built low cost drones , single board computers and lightweight operating system virtualization technologies such as unikernels / dockers @xcite . our paper lays out the preliminary ideas on such a system that can be instantaneously deployed on demand . to enable a lightweight micro cloud infrastructure in the sky using drones , we bring together a couple of fascinating recent innovations in computing : single board computers such as the raspberry pis ( pi ) and lightweightos virtualisation technologies such as dockersand integrate them with lightweight quadcopters . there were two possibilities for us to build the cloudrone : integrate the pi with an off the shelf ready made quadcopter or design and build an entire quadcopter from scratch with the pi powering both the drone as well as acting as a micro cloud server . we decide to go for the latter since we wanted to have the drone as lightweight and low cost as possible . there are two main challenges in building a low cost quadcopter . it has to be light weight yet robust and it has to be economical . the cost can not be controlled on components such as the motors , electronic speed controllers ( escs ) and the battery . so the focus was on the other two major areas to reduce the cost . the frame and the rc controller . for building the frame , some materials like carbon fiber are light and strong but they are expensive compared to other materials like aluminium shafts or high quality reinforced plastic . the frame of the quadcopter consists of carbon infused plastic , the landing gear is made of light - weight square aluminium sections fitted with high quality foam sponge balls to withstand crashes . four brushless dc motors ( 1400kv , 800gms thrust ) connected with carbon fibre propellers ( 8x4.5 " ) drive the drone . the speed of the motors is controlled by four electronic speed controllers ( peak current 30a ) with in - built battery eliminator circuits that provide 5v , 2a for powering up other circuits like the kk2 flight controller board and the pi . these escs are instructed by control signals that are generated by the output pins of the kk2 flight controller board . the kk2 board has a variety of modes to control different uavs like tricopters , quadcopters , hexacopters and octocopters easily . so any extension to the current quadcopter design can be easily accommodated and realised quickly . the kk2 has input pins which are generally connected to rc receivers . but rc transmitters and receivers tend to be quite expensive and make up a large portion of the entire cost .    a good quality 4 channel rc transmitter costs more than  200 . to build a low cost quadcopter that can act as a micro cloud server , we eliminated the rc transmitter and receiver entirely and replaced it with the pi . the pi has two important functions i.e. , controlling the drone s movements and to run micro services . since controlling the quadcopter is now being done by the pi , it can be programmed to fly autonomously along with manual override . we used the pi 2 model b equipped with quad - core armv7 cpu 900 mhz . the ram memory is 1024 mb . it requires power feed about 800 ma . the pi runs an arm version of ubuntu 14.04 . we used a edimax wifi dongle ( ew-7811un ) for providing the wifi interface . 0.33     0.33     0.33       the pi generates a portable wi - fi access point to which the laptop is connected . a gaming joystick is connected to the laptop to control the drone manually . by using a python module called pygame  , the joystick emulates an rc transmitter with different buttons and 2-axis joysticks , acting as controls for roll , pitch , yaw and throttle . a graphical user interface was also created to monitor the orientation , altitude , battery level , gps co - ordinates . the pi and the base laptop are connected by a software framework called ros ( robot operating system )   which simplifies sending and receiving data with the concepts of publishing and subscribing data . data which may be originating from a controller or from a sensor is published onto a  topic " , from which any other component of the entire system can access it by subscribing to that topic .    in our design , the joystick values are published to a topic which is subscribed by the pi . the pi then uses that data to send pulse width modulated control signals to the kk2 board . the program that runs on the pi has been designed to give signals exactly as an rc receiver would . so , for the kk2 board the signals appear to be coming from an rc receiver . based on the signal it receives , the kk2 board sends corresponding signals to the escs , which in turn send the brushless dc ( bldc ) motors the appropriate currents to change their speeds . all the electronics , motors and escs are powered using a lithium polymer battery ( 4200mah 3s 35c ) .    in order to maintain stability the kk2 implements p - i ( proportional - integral ) controller . the user has to manually set values to constants that govern the flight control . in the kk2 board , there are two important constants : p - gain , and i - gain for roll , pitch and yaw . these constant values will be the same for both roll and pitch owing to symmetry . firstly considering only one axis , say roll , setting i - gain to zero , p - gain is increased until the quadcopter produces oscillations periodically . if p - gain is too low , the reaction to roll left or right will be sluggish . if the p - gain is too high the quadcopter produces high frequency oscillations . based on observations , the p - gain is set to an appropriate value . i - gain is then increased to ensure that the quadcopter does not drift from the set - point . if i - gain is too high , there will be a lot of overshoot from the set point and if it is too low there will be a lot of drift . the quadcopter will not reach the set point forever . in order to tune the roll axis , the quadcopter is tied along the pitch axis to two parallel cots which act as a stable test rig . the accelerometers and gyroscopes in the kk2 board are then calibrated . we place the quadcopter on a tested perfectly flat surface and run this calibration on the kk2 board to make sure that the pitch and roll angles are zero while hovering . the escs are then calibrated to start at the same time . this is done by pressing the back and enter buttons simultaneously while switching on the kk2 board with all the escs connected . the joystick connected to the laptop is set to 100% throttle . the escs beep to indicate that they have been calibrated . propellers , even the high quality ones may have some weight discrepancies which affect the quadcopter s flight when the motors spin at 14,000 rpm . so with a simple balancing test , we find out which side of the propeller is lighter than the other . to make the weights of both sides of the propeller equal , we stick cellophane tape on the lighter side . similarly to reduce the vibrations generated by the bldc motors we first measure the vibrations using an android app and stick cellophane tapes accordingly to all four motors and minimize these vibrations as they can corrupt the values of the accelerometer and gyroscope in the kk2 board . the quadcopter s final weight is approximately 1.5 kgs . the maximum altitude to which it was flown was 50 feet . with the current setup it can fly up to 100 feet theoretically . the effectiveness of the landing gear was also tested . the quad copter was dropped from a 2-storey building ( 45 feet high free fall ) and no components were damaged . the quadcopter s flight time has been measured in two different case scenarios . first with a 2200mah battery with 60c discharge rating and then with a 4200mah 35c battery . the flight duration and the throttle levels vary noticeably in these two scenarios . when we use the 2200mah battery the payload of the quadcopter was bout 50 gms higher than when the 4200mah battery was used . these are the observations :    * the quadcopter in the first scenario lifts only at 70% throttle level and the flight time is very poor owing to the smaller capacity of the battery and slightly heavier load . the estimated and realised flight times are only around 2 - 3 minutes . * on redesigning certain aspects of the quadcopter to compensate for the heavier 4200mah battery , the overall weight of the quadcopter was reduced by 50 gms . the quadcopter now lifts at 55% throttle and was able to have sustained flight for over 8 minutes . if a smaller battery is used in the lighter quadcopter , the performance is still not as good as using a bigger and higher capacity battery . * by using a higher mah rated battery , the flight time can be further increased to an estimated 15 minutes . the quadcopter s cpu consumption was measured during the flight operation by using _ sysstat _ tools . on average our quadcopter consumes 30% for the flight control of the quadcopter leaving 70% for other processes .    to measure the wifi coverage of the quadcopter while it was in the air , we measured the signal strength from various distances ( with clear line of sight ) using wigle . we observed signal strengths varying from -65 dbm to -78 dbm within the distances 10 to 50 meters north , east and west of the direction of the wifi dongle . south being the direction away from the dongle had relatively poor coverage ( between -75 dbm to -89 dbm ) . wifi coverage can be enhanced by using additional portable lightweight wifi hotspots that can provide wider coverage ( such as the tp - link mr3040 ) . our vision to build a micro cloud infrastructure in the air is to use a swarm of pis on drones acting as micro cloud servers running dockers enabling us to run several lightweight containers . the lightweight nature of docker containers significantly reduce the size of image compared to the heavyweight vms @xcite . for instance , we can build a minimal static web server with only 2 mb size . this can be further reduced with new lightweight os virtualisation technologies such as includeos .    to interconnect the pis as a swarm / mesh , each of these devices also act as mobile routers operating in two wifi communication modes : one operates in the wifi ad - hoc mode to allow optimised link state routing ( olsr ) protocol , constructing a mobile ad hoc network ( manet )  @xcite ; the other operates in the access point ( ap ) mode to allow user devices to connect with dhcp . the benefit of using olsr protocol compared to other manet routing protocols is it uses a special mechanism called multi - point relay ( mpr ) to reduce the number of flooded messages . only a few devices that are located in strategically better spatial positions are chosen to relay ( i.e. , re - transmit ) messages in the path from a source device to a destination device . the mpr mechanism helps reduce overall energy consumption .    to create a swarm of micro cloud servers , we use the docker swarm technology that allows us to create a cluster of multiple docker hosts ( pi running docker ) and migrate the service containers across the cluster . specifically , there are two types of nodes classified in the swarm . 1 . swarm manager - who manages overall resources ( e.g. , swarm members , number of running containers ) and decides where to place the service containers . 2 . swarm agent - nodes registered with the swarm manager . when the swarm manager and agents are created , they have to register with the discovery backend as members of the swarm . the discovery backend maintains an up - to - date list of swarm members and updates that list with the swarm manager . the swarm manager uses this list to assign tasks and schedule the service containers to the agents . to determine where to place the new container in the swarm , the swarm manager uses either _ spread _ or _ bin packing _ strategy to compute the rank regarding node s available cpu , ram and the number of running containers . with the spread strategy , the swarm manager gives a priority to the node who has the largest available memory or has the minimum number of running containers . on the other hand , the bin packing strategy tries to pack as many containers to a node until reaching its maximum capacity ( e.g. , ram , cpu ) . in this sense , a swarm can optimise the number of nodes running the containers and leave room for future assignment which may require large space of resources . @xmath0targets to provide localised service infrastructure in challenged network environments . such scenarios refer to the emergency and post - disaster situations wherein traditional communication services are completely inoperable . an example of @xmath0 s deployment is illustrated in figure  [ img : cloudrone ] . a base camp command center with backhaul internet connectivity ( e.g. , satellite link ) is set up close to the affected area . a cluster of drones can fly to cover the target area ( and could also land ) , forming a mesh network and provide localised services to the users ( immigrants or post - disaster victims ) on the ground via wifi . a variety of ( crucial ) services ( lightweight docker containers ) can be either pre - loaded onto the pi or on demand from the ground . the manet of drones , facilitates the swarm manager to communicate and control the cluster remotely from the base camp through the long haul link . the swarm manager can update the necessary services from the internet and disseminate throughout the cluster .        in some operations , the cluster can be out of contact with the swarm manager ( i.e. , the target area is far away from the command center ) . to deal with any interruption of service , we can create a primary swarm manager operating as the main point of contact and multiple replicas to be the backup swarm managers . using this feature , the replicas can seamlessly take over the functionalities from the primary swarm manager when it fails . if the cluster fails to contact the primary swarm manager , the most powerful replica automatically takes over the control . battery life time is one of the major challenges in our deployment , since the flight time is limited to 15 minutes . to mitigate this problem , we envision that the drone can either fly or land to provide service access . in case , the battery is low , our drone can decide to land on the ground to save the battery while the docker containers and wifi access point are still operating to provide service access . the preliminary benchmarking presented in this paper focusses on understanding the scalability issues of lightweight os virtualisation technology such as docker on a pi . we focus on two main questions 1 ) how many docker containers can a single pi support ? 2 ) how many user requests can a single container running on a pi support ? 0.24     0.24     0.24     0.24     0.33     0.33     0.33       the first evaluation is to explore the maximum number of containers that could be operated concurrently over a single pi . to scale up the containers , the docker image could be prepared as small as possible to minimize memory footprint . consequently , the memory allocation for kernel to handle a web server process can be optimized . for this , we use a nano web server image developed in assembly code ( size is less than 90 kb ( included index.html and a small jpg file ) ) . to benchmark the capability of a pi ( pi 2 ) , we base our evaluation on memory consumption , cpu utilisation and the creation time ( time taken to create a container ) by using _ sysstat _ , a collection of performance monitoring tools for linux . in our first attempt , we were able to spin up only 37 containers , even though the memory usage and cpu utilisation were only 40% and 18% respectively . we then hacked the docker daemon , to scale up the deployed containers to 2408 . the procedures that we used is summarised below :    * _ tweak docker command : _ we disabled some docker options that are unnecessary for running a web server . especially , running with a dedicated ip stack per container involve a huge resource usage , so we run each container with @xmath1 . we also disabled log driver ( @xmath2 ) to let docker use less resources . * _ unlock docker limit : _ linux distributions use systemd to start the docker daemon . we customised some parameters to unlock the limitation of maximum number of running containers . the number of processes ( limitnproc ) and number of queued signals ( limitsigpending ) were set to infinity . * _ _ tweak docker configuration:__with the default configuration , docker daemon is run with many options which is not required for a simple web server such as ipv6 , proxies , ip forwarding , log - driver . we disabled all these options to reduce the memory usage of docker daemon . the results for our optimised docker are depicted in figure  [ result : benchmark ] where we were able to scale up the number deployed containers to 2408 . specifically , memory usage is a key factor that limits the capability of running the containers on a pi . as shown in figure  [ result : mem ] , the memory usage increases gradually when a container is added . the initial memory usage before creating the first container was about 98 kb ( 9.89% ) . we hit the limit of 2408 containers where there is no space for available memory . figure  [ result : cpu ] shows the average utilisation of the quad - core cpu of the pi . over the first thousand of deployed containers , the average cpu usage is very low ( about 0.2% ) with a few spikes . however as expected , the cpu usage increases significantly in the high load state ( when the number of containers is larger than 2000 ) . as mentioned in section  [ benchmarking the drone ] , the quadcopter consumes about 30% of cpu resources . the available space can be around 70% which is sufficient to provide multiple services with docker containers . hence when the quadcopter is flying , we will not be able to have 2408 containers running in parallel . however , a @xmath0can still support a large number of concurrent containers . figure  [ result : time ] depicts the results of creation time where each point denotes the time it took for the nth container to start up . the creation time for each container varies from 0.62 s to 38.37 s depending on the current cpu load and memory usage . in addition , we also plot the cumulative creation time of the 2408 containers ( figure  [ result : cumulativetime ] ) . the pi spent approximately 1 hr and 50 minutes to spin out 2408 containers . on average each container requires about 2.79 s to start up the web server .    _ key takeaway message _ : a single pi ( pi 2 model b ) can support significant amount of concurrent lightweight services .      in order to investigate the feasibility of using lightweight containers as a platform for the @xmath0 , we aim to evaluate the scalability of each of these containers running on a single pi while serving a large number of requests . we deploy a minimal static web server using httpd docker base image using a similar configuration as the experiments in the previous section . the benchmarking scenario represents the _ cloudrone s _ operation where users on the ground can access the services provided by the _ cloudrone _ through a wireless interface . using the _ ab _ - apache http server benchmarking tool , we conduct stress tests on the _ cloudrone _ while scaling the number of concurrent users from 10 to 250 . the total number of request was set as 10000 transactions per experiment . for instance , in case of 10 users , 1000 of requests were sent by each user . the measured rtt via ping tests with a clear line of sight at 50 feet was between 8ms-10ms . figure  [ result : cdf_delay ] illustrates the cdf of response time from the web container running on pi while varying the number of simultaneous users accessing the service . as shown in the figure , a single deployed docker container is capable of serving a large number of concurrent users . as expected , the average response time increases when the number of concurrent users is scaled up . figure  [ result : cpu_utilization ] and  [ result : cpu_load ] plot the cpu utilisation and cpu load of the pi using _ sysstat _ tools . the cpu utilisation increases almost 20% when the number of users increase from 10 to 250 . this increases the response time for the processes . even though the utilization has not reached the capacity of cpu , processes still run slower as the cpu load increases . figure  [ result : cpu_load ] shows the average cpu load of the pi ( sampled in one min intervals ) . as the arrival rate of user requests increases , the amount of computational work need to process also increases this has impact on the response of time(figure  [ result : cdf_delay ] ) . _ key takeaway message _ : a docker container running on a single pi ( pi 2 model b ) can support significant amount of concurrent users . our preliminary benchmarking demonstrates that the pi is capable of functioning as a great micro cloud platform . our benchmarking was carried out using a lightweight web server serving a lightweight webpage . it is important that the scalability of the @xmath0should also be tested with heavier web servers serving applications such as openstreetmaps . each application will have different memory and cpu requirements and hence the number of containers that can be instantiated will vary depending on the type of applications which directly influences the size of the container ( e.g. , packages , data , library ) .    another scalability challenge is to support a larger number of users within an area . increasing the number of users causes an adverse influence on the response time which will cause service degradation . we envision , the different services will be provided by a swarm of drones and hence appropriate load balancing using techniques such as application layer anycast  @xcite could be used . _ cloudrones _ have two main challenges in terms of providing a reliable service . first , mobility poses a critical challenge . during mobility , ongoing sessions may break and sessions need to be reestablished . second , considering the distributed nature of the services across a mesh of drones , identifying the location of the service across a mobile adhoc network is challenging . to solve the latter , techniques such as multicast dns ( mdns )  @xcite or application layer anycast could be used  @xcite . another potential way to mitigate the problem of mobility and service discovery is to explore new architectures that utilise information centric networking )  @xcite . icn architectures such as named data networking ( ndn )  @xcite or scandex  @xcite decouple the content from the location thus removing the need for the current end - to - end client server model such that the service and/or content can be served directly by any host that currently has the service / content . icn thus integrates the provisioning of content with the locationless notion of information delivery in icn allowing different flavours of caching , from on - path caching to edge caching through a farm of surrogate micro servers running on the _ cloudrones _ that can be quickly integrated into the overall ( icn ) routing fabric without the need for dns redirection or other solutions of the current internet . this inherently addresses the issues of mobility and reliability . as ongoing work , we are in the process of developing an icn architecture for _ cloudrones _ via the eu umobile project  . although _ cloudrones _ demonstrate excellent potential for deploying localised service infrastructures in areas where access to services are crucial but are beyond reach - there are still major challenges that need to be surmounted . drones such as quadcopters have reduced flight times due to battery life . we envision this situation will change in the near future with better innovations in battery design and production or innovations in alternate sources of power e.g. hydrogen powered drones have flight times upto two hours  . _ cloudrones _ also need not be in the air for their entire flight duration to provide access to its services . we envision that _ cloudrones _ can be flown to an area and then can provide it s localised services from the ground ( ideally powered by an energy source on the ground ) . there are tight regulations in flying drones such as quadcopters . these rules have been laid out by civil aviation authority ( in the case of uk ) . hence these rules should be adhered to and in some cases may be restrictive e.g. land or fly in a congested area . however , we believe , @xmath0deployments will fall under the commercial aerial work , and hence special permission from the aviation authority will be required to fly . in this paper , we present @xmath0- a preliminary idea of deploying a lightweight micro cloud infrastructure in the sky using indigenously built low cost drones , single board computers and lightweight operating system virtualization technologies . we describe an initial design of the @xmath0and provide a preliminary evaluation of the proposed system mainly focussed on the scalability issues of supporting multiple services and users . as part of future work , we plan to conduct large scale evaluation trials benchmarking the @xmath0performance while in the air ( in terms of throughput , latencies and energy ) across a wide set of scenarios . we are also in the process of integrating docker with ndn and performance benchmarks will be carried out . finally , the current @xmath0design does not fly autonomously and hence is strictly limited in terms of distance it can cover without manual intervention . as part of future work , we intend to build autonomous flying capabilities .
Abstract:

Answer: recent years have witnessed several initiatives on enabling internet access to the next three billion people . access to the internet necessarily translates to access to its services . this means that the goal of providing internet access requires access to its critical service infrastructure , which are currently hosted in the cloud . however , recent works have pointed out that the current cloud centric nature of the internet is a fundamental barrier for internet access in rural / remote areas as well as in developing regions . it is important to explore ( low cost ) solutions such as micro cloud infrastructures that can provide services at the edge of the network ( potentially on demand ) , right near the users .    in this paper , we present @xmath0- a preliminary idea of deploying a lightweight micro cloud infrastructure in the sky using indigenously built low cost drones , single board computers and lightweight operating system virtualization technologies . our paper lays out the preliminary ideas on such a system that can be instantaneously deployed on demand . we describe an initial design of the @xmath0and provide a preliminary evaluation of the proposed system mainly focussed on the scalability issues of supporting multiple services and users .


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 geometric phases in quantum systems , ever since its discovery by berry @xcite , has attracted considerable interest both from the theoretical and experimental viewpoints . several applications of this phenomenon in different areas of physics have also been studied@xcite . comparatively much less work has been devoted to the question of the dynamical evolution of such systems in the presence of a weakly dissipative medium . basically the available results can be summarized as follows : nonhermitian operators lead to a modification of berry s phase@xcite , stochastically evolving magnetic fields produce both energy shift and broadening@xcite , phenomenological weakly dissipative liouvillians alter berry s phase by introducing an imaginary correction@xcite or causing damping and mixing of the density matrix elements@xcite . @xcite obtain their results studying the eigenmatrices of the complete liouvillian superoperator in a time - independent basis , while gamliel and freed@xcite , find closed formal expressions for the density matrix representation in the instantaneous frame of the hamiltonian in the adiabatic limit , and in the weak dissipation approximation . however , even in this regime the results in the cases studied can only be extracted through approximations ( in addition to the adiabatic approximation ) or numerical computation .    in the present work , we consider the celebrated example of a spin 1/2 particle coupled to a slowly evolving magnetic field in the presence of a weakly dissipative medium as represented by a lindblad form superoperator@xcite , incorporated in a convenient and physically motivated frame . in the absence of dissipation , a simple geometrical interpretation of the results emerges in terms of both the bloch s vector@xcite and the phase vector@xcite . the geometrical phase appears as a `` delay '' or an `` advance '' in the precession period of the bloch s vector with respect to the period dictated by the magnetic field s frequency . the precession of this vector occurs at a fixed angle with respect to a fixed axis about which the external magnetic field precesses , as usual . we will introduce the dissipation via a semi - group type dynamics , such that dissipation does not alter the precession frequency . however , the @xmath0 , @xmath1 and @xmath2 components of bloch s vector are altered in different ways by nonunitary effects . in particular , the modulus of bloch s vector will shrink . the relation between the reduction of bloch s vector modulus and the loss of coherence has been explored by stodolsky and collaborators@xcite . in order to make it explicit , we introduce the linear entropy ( or idempotency defect ) as a measure of purity loss@xcite    @xmath3    and write    @xmath4    where @xmath5 is the bloch s vector . we get    @xmath6    which confirms that the shrinking of the modulus of @xmath7 is a measure of coherence loss .    in order to discuss dissipative effects on interference patterns due to berry s geometrical phase the natural framework is to consider the evolution of the density matrix , which completely characterizes the interference effect due to the different geometric phases acquired by the eigenvectors of the hamiltonian of the system , even in the absence of dissipation@xcite . of course , this time we add a nonunitary liouville operator contribution to the dynamics .    in section 2 we present the two level quantum system of a spin 1/2 in the presence of an external magnetic field precessing with constant angular velocity around a fixed axis . the master equation of the model is written in what we call the diagonal frame where lindblad superoperators are introduced to describe the nonunitary part of the liouville operator . in subsection 2.1 we consider the adiabatic limit of the quantum system in thermal equilibrium with a reservoir of electromagnetic fields and in subsection 2.2 we consider the case of a dephasing process . in both cases , we treat the coupling of the quantum system with its environment in the weak regime . in section 3 we show how the geometric phases and the dissipation effects acts on the @xmath2 component of the magnetization of the spin 1/2 coupled to a reservoir of electromagnetic fields at thermal equilibrium . in section 4 we summarize our conclusions ; finally , in appendix a , we present the non - unitary part of the liouvillians of the models under consideration in the instantaneous basis of the hamiltonian . we consider a spin 1/2 variable ( two level model ) coupled to a time dependent magnetic field precessing around the z - axis . the unitary contribution for this evolution is given by the hamiltonian    @xmath8    written in the basis of the eigenstates of the @xmath2-component of the spin , where @xmath9 is the norm of the external magnetic field , @xmath10 its azimuthal angle , @xmath11 the precession frequency and the constant @xmath12 , being @xmath13 land s factor and @xmath14 the bohr magneton . we are using natural units ( @xmath15 ) . for the sake of later calculations , it is convenient to define two unitary transformations : the first one , @xmath16 , takes us to the rotating frame where the hamiltonian is no longer time dependent ; the second one , @xmath17 , diagonalizes the effective hamiltonian ( time independent ) that drives the dynamics of the final matrix representation of the density operator . after the first transformation , @xmath18 , the density matrix and the hamiltonian read    @xmath19    and    @xmath20    in analogous manner , after the second transformation we get , in the diagonal frame , the density matrix    @xmath21    and the effective hamiltonian    @xmath22    where @xmath23 . the rotation matrix @xmath24 is equal to    @xmath25    one possible way to add dissipative contributions to the above dynamics is to include a lindblad type superoperator in the evolution equation in the diagonal frame    @xmath26 + k { \cal l}_d \rho_d ( t ) , \label{6.2}\ ] ]    where @xmath27 is the dissipation constant . the weak coupling regime is characterized by the condition @xmath28 . our aim is isolating , in the density matrix , the effects of dissipation on interference due to geometric phases . the representation of the density matrix in the diagonal frame is not very enlightening for this purpose . this is better realized in some basis of the instantaneous eigenvectors of hamiltonian ( [ 4 ] ) , as was done by gamliel and freed @xcite . we define @xmath29 , the matrix density in a basis of the instantaneous eigenvectors of hamiltonian ( [ 4])@xcite . the relation between @xmath29 and @xmath30 is    @xmath31    where the matrix @xmath32 is equal to    @xmath33    the time evolution of @xmath29 is given by    [ 9 ] @xmath34                             + k { \cal l}_i \rho_i ( t)\label{9.1}\ ] ]    where    @xmath35    and @xmath36 is obtained from @xmath37 through a similarity transformation equivalent to ( [ 7 ] ) .    before we specialize our discussion to any particular liouvillian , we study the adiabatic limit of eqs.([9 ] ) for the coupling constant @xmath27 in the weak regime . in this regime , the matrix @xmath36 is written as a linear superposition of the elements @xmath38 . we remind that the density matrix of a two level model must satisfy two conditions : _ i _ ) @xmath39 @xmath40 and _ ii _ ) @xmath41 . as a consequence of those conditions , the density matrix has only two independent elements . we take the elements @xmath42 and @xmath43 as our two independent entries . the general form for the time equations of those two elements in the weak coupling regime is    [ 11 ] @xmath44    and    @xmath45    we make the change of variables :    @xmath46    with @xmath47 and @xmath48 . we may disclose the time scale @xmath49 in the differential equations by introducing the following transformation upon the time parameter@xcite :    @xmath50    where @xmath51 .    with the new variables , eqs.([11 ] ) become    [ 13.0 ] @xmath52    we take the differential equation for @xmath53 to exemplify the discussion of the adiabatic limit of eqs.([13.0 ] ) . at this point we will follow closely the references @xcite . we point out that the adiabatic approximation is not recovered by an @xmath11 expansion of the terms on the r.h.s . of eqs.([13.0 ] ) . the volterra equation@xcite obtained from eq.([13.2 ] ) is    @xmath54    in the limit @xmath55 , the riemann - lebesgue theorem @xcite gives that    @xmath56    if @xmath57 is a piece - wise continuous function in the interval [ @xmath58 and @xmath59 . as a consequence of this theorem the last three terms on the r.h.s . of eq.([13.3 ] ) vanish in the adiabatic limit ( @xmath60 ) .    integrating by parts eq.([14 ] ) we get    @xmath61   - \frac{1}{it } \frac{1}{\alpha \mu b }               \int_0^s \frac{d}{ds^\prime}\left ( f(s^\prime ) \right )                                                    e^{i\alpha \mu tb s^\prime } ds^\prime . \label{15}\ ] ]    from eqs.([9 ] ) and ( [ 11.2 ] ) , the coefficients @xmath62 and @xmath63 have the general dependence on @xmath11 and @xmath27 :    [ 16 ]    @xmath64    and    @xmath65    the integrands of all integrals on the r.h.s . of eq.([13.3 ] ) which contain the oscillatory function @xmath66 , with @xmath67 or @xmath68 , after integration by parts , each of those integrals has its order in t decreased by one unit , and acquires a multiplying constant of value @xmath69 or @xmath70 . from the conditions satisfied by the two simultaneous regimes : _ i _ ) adiabatic limit ( @xmath71 ) and _ ii _ ) weak coupling limit ( @xmath72 ) , we can neglect those terms in comparison to the first two terms on the r.h.s . of eq.([13.3 ] ) . the two previous inequalities do not impose any constraint to the ratio @xmath73 , though . the differential equation satisfied by @xmath74 in the adiabatic limit and weak coupling regime is    @xmath75   \tilde{\rho}_{12 } ( s ) . \label{17}\ ] ]    the term proportional to @xmath76 was dropped , since it is of higher order in @xmath77 . for @xmath78 , the off - diagonal elements of the density matrix vanish before the external magnetic field @xmath79 returns to its configuration at @xmath80 . the other uninteresting situation from the point of view of the appearance of an imaginary correction to the geometric phase is the condition @xmath81 when the dissipation effects can still be neglected after one period . in this work we discuss the case @xmath82 when both terms on the r.h.s . of eq.([17 ] ) contribute to the time evolution of @xmath83 .    by a similar discussion we obtain the time equation of @xmath84 in the adiabatic and weak coupling regimes ,    @xmath85 . \label{18.1}\ ] ]    from eqs.([9 ] ) we can affirm that    @xmath86    in order to verify if we can get imaginary phases from eqs . ( [ 17 ] ) and ( [ 18.1 ] ) due to the coupling of the quantum system to a dissipative medium , we write the density operator at @xmath80 in the instantaneous basis of hamiltonian ( [ 4 ] )    @xmath87    we have @xmath88 , @xmath89 and @xmath90 . being @xmath91 the time evolution of the eigenvector @xmath92 , we have    @xmath93    where the time evolution of @xmath91 is driven by @xmath94 . differently from gamliel and freed , we include the phases coming from the unitary evolution and the geometric phase in the dyadic product @xmath95 . in the adiabatic approximation , we get    @xmath96    where @xmath97 is the geometric phase acquired by the eigenvector @xmath98 , @xmath99 is the instantaneous eigenvector of @xmath100 ( @xmath101 ) and @xmath102 . the density operator at any time , in the adiabatic limit is    @xmath103    from eq.([22 ] ) we recognize that the phase @xmath104 in the element @xmath105 is just the difference of the geometric phases of the instantaneous eigenstates @xmath106 and @xmath107 in the absence of dissipation . the dynamics of the coefficients @xmath108 is ruled by the nonunitary evolution of the quantum system and it is independent of the particular choice for the instantaneous eigenstates of the hamiltonian , up to a multiplicative constant@xcite . in the model under consideration ( see hamiltonian ( [ 4 ] ) ) , the eigenvalues @xmath109 , @xmath110 and @xmath90 , are time - independent . it is simple to get the time equations of @xmath111 and @xmath112 from eqs.([18.1 ] ) and ( [ 17 ] ) , respectively    [ 23 ]    @xmath113    the constants @xmath114 , @xmath115 and @xmath116 depend on the particular master equation that describes the behaviour of the quantum system interacting with the dissipative medium . in the next sub - sections we consider two particular interactions of the two level model with a reservoir : _ i _ ) two level model in thermal equilibrium with a reservoir of electromagnetic fields ; _ ii _ ) dephasing process in a two level model . as discussed before , to incorporate the dissipative effects in the two level model we introduce the lindblad superoperator in the diagonal frame . the master equation of the spin 1/2 model coupled to a reservoir of electromagnetic fields in thermal equilibrium in the diagonal frame is@xcite    @xmath117 + k\left ( \overline{n}+1\right ) \left ( 2{\sigma } _ { -}\rho_d \left ( t\right ) { \sigma } _ { + } -\rho_d\left ( t\right ) { \sigma } _ { + } { \sigma } _ { -}-{\sigma } _ { + } { \sigma } _ { -}\rho_d\left ( t\right ) \right ) + \nonumber \\ & & + k\overline{n}\left ( 2{\sigma } _ { + } \rho_d\left ( t\right ) { \sigma } _ { -}-\rho_d\left ( t\right ) { \sigma } _ { - } { \sigma } _ { + } -{\sigma } _ { -}{\sigma } _ { + } \rho_d \left ( t\right ) \right ) ,             \label{24}\end{aligned}\ ] ]    where @xmath27 is the dissipation constant at zero temperature and @xmath118 is the average number of excitations of the weakly coupled thermal oscillators at inverse temperature @xmath119 . an important requirement for the introduction of this lindblad type superoperator is that it leads asymptotically to a thermal equilibrium . in appendix a we give the master equation of this physical process in the instantaneous basis of hamiltonian ( [ 4 ] ) for arbitrary value of @xmath11 .    from the master equation in the instantaneous basis of the hamiltonian we obtain the equations for @xmath120 and @xmath121 . these equations in the adiabatic and in the weak coupling limits become    [ 25 ]    @xmath122    the constant @xmath27 does not come up on the r.h.s . of eqs.([25 ] ) due to the time variation of any classical parameter that characterizes the reservoir . from its explicit definition@xcite it can not be written as : @xmath123 , where @xmath124 and @xmath125 are two regular time dependent functions . the same is true for @xmath118 . therefore the imaginary phase @xmath126 defined as : @xmath127 , with    @xmath128    is not geometric . the solution of eq.([25.1 ] ) is    @xmath129   e^{-2k\int_{0}^{t } ( 1 + 2\bar{n } ) dt^\prime } . \label{27}\ ] ]    the exponential decay on the r.h.s . of eq.([27 ] ) means that the population of the instantaneous eigenstates of hamiltonian ( [ 4 ] ) varies in time and consequently the adiabatic theorem is not valid for this dissipation mechanism . exactly soluble models are always important checks to approximation schemes . eq.([24 ] ) is exactly solved and the solutions are    [ sol - ex ]    @xmath130 + \rho_{11}^d(0)e^{-2k(2\overline{n}+1)t } , \label{sol - ex1 } \\ \nonumber\\ \nonumber\\ \rho_{12}^d(t)&= & \rho_{12}^d(0)e^{-(2i\lambda_1 + k(2\overline{n}+1))t}.    \label{sol - ex2}\end{aligned}\ ] ]    it is straightforward to obtain the adiabatic and weak coupling limit eqs.([25 ] ) from eqs.([sol - ex ] ) . another interesting process well studied in the standard textbooks@xcite is the phase destroying process which might appear due to elastic collisions . in general , those effects are incorporated in the master equation of the two level model , besides the energy dissipation process studied in subsection 2.1 . since we are studying the coupling of the spin 1/2 to a dissipative medium in the weak coupling limit , the inclusion of the phase destroying process in eq.([24 ] ) gives corrections to the coefficients @xmath131 in eqs.([25 ] ) . due to the linearity of the equations , the new imaginary phases coming from the dephasing process are added to the ones obtained previously . for the sake of simplicity we study the imaginary phases acquired by the variables @xmath111 and @xmath112 only due to the dephasing process . the master equation written in the diagonal frame is    @xmath132   + \frac{k}{2 }   \left ( \sigma_z \rho_d ( t ) \sigma_z   - \rho_d ( t ) \right ) . \label{28}\ ] ]    in appendix a we present the master equation of the spin 1/2 with the dephasing effect included for arbitrary value of the angular velocity @xmath11 of the external magnetic field . in this subsection we study the time equations of coefficients @xmath120 and @xmath121 in the adiabatic and weak coupling regimes . taking into account our discussion in section 2 , eqs.([a2 ] ) and ( [ 23 ] ) we obtain    [ 29 ]    @xmath133    that have the solutions :    [ 30 ]    @xmath134    from eq.([30.1 ] ) we conclude that the adiabatic theorem is valid in this process , since the population at each quantum state does not vary along the adiabatic process . by analogous reasons to the ones discussed in subsection 2.1 , the phase in eq.([30.2 ] ) is not geometric but a time dependent imaginary phase that destroys the off - diagonal elements of the density matrix . eq.([28 ] ) is exactly solved . the solutions ( [ 30 ] ) are easily recovered from the exact solutions when we calculate them in the adiabatic and weak coupling limits .    in order to understand why references @xcite-@xcite give an imaginary correction to the geometric phase and the models studied here do not , we compare eq.(2.3 ) of reference @xcite with our eqs.([26 ] ) and ( [ 30.2 ] ) . the coefficient that multiplies of variable @xmath135 in eq.(2.3 ) of reference @xcite has the form @xmath136 . since the vector states depend on a periodic external parameter @xmath137 , this coefficient , in the adiabatic approximation , corresponds to a correction to berry s phase written as a closed curve in parameter space . the coefficient @xmath138 that multiplies @xmath112 in eq.([25.2 ] ) does not arise from any variation of an external parameter . the same is true for the coefficient @xmath27 in eq.([29.2 ] ) . that is the reason that allows us to claim that the exponentials acquired from eqs.([26 ] ) and ( [ 30.2 ] ) have no geometric origin , which means that the suppression terms are a function of time instead of some path parameter in order to verify how the geometric phases and the dissipation affect the physical quantities , we return to the spin @xmath139 model in the presence of a reservoir of electromagnetic fields in thermal equilibrium . under the initial condition @xmath140 , with @xmath141 , in the adiabatic limit and weak coupling limit we end up with    [ 12.todas ]    @xmath142    we have introduced the berry s phase tracer @xmath143 that helps us identify the contribution of berry s phase to physical quantities . at the end , we take the tracer equal to one .    thus we see that , in this model , dissipation does not affect berry s geometrical phase , but only makes it harder to observe their interference effect : such information is contained in the off - diagonal terms of the density matrix in the instantaneous frame , which vanish . by the way , we can observe that the system considered in this contribution is analogous to the case of the classical foucault pendulum where dissipation diminishes the amplitude , but do not affect the rotation of the oscillation plane . the definition of the geometrical phase can also be given in terms of a phase vector as in reference@xcite , where the discussion is confined to pure states . it is , however , a relatively simple matter to extend the definition of the phase vector for mixed states . in this case we find that the loss of coherence will shorten the phase vector in a manner which is completely analogous to what happens to bloch s vector . this norm reduction , as in the present case , reflects the asymptotic vanishing of off - diagonal density matrix elements@xcite . since we are interested on interference effects due to geometric phases , bloch s vector , defined in eq.([2 ] ) , is suitable to provide a graphic visualization of the density matrix @xcite . from eq.([3 ] ) , we have that these effects are clearly described by bloch s vector in the instantaneous frame . in this frame the projection of bloch s vector sweeps the @xmath144 plane and makes an angle smaller than that of the magnetic field by an amount proportional to the solid angle @xmath145 , while its length decreases exponentially with a time rate of @xmath146 . we illustrate in figure [ figura01 ] the shrinking of the projection of bloch s vector in the @xmath144 plane , due to the presence of dissipation . even though we have discussed the case @xmath147 , we choose @xmath148 in order to show more clearly the decreasing of this projection in an interval @xmath149 . shortening in @xmath2 is faster than in the @xmath144 plane , causing a time dependent azimuthal angle . figure [ figura02 ] shows the plot of the time evolution of bloch s vector for @xmath150 . for @xmath151 bloch s vector has only non - zero @xmath2-component and @xmath152 .    for time @xmath153 , such that @xmath154 ( @xmath155 and @xmath156 is an integer ) the off - diagonal elements of the density matrix @xmath29 depends on the chosen condition satisfied by @xmath157 @xcite . for the sake of comparison with experiments it is necessary to calculate dissipation and geometrical phases effects on measurable quantities . in this model the natural candidates are the components of the magnetization vector @xmath158 . let us consider the @xmath2-component of magnetization whose average value has the expression @xmath159 . since the trace is independent of the particular basis applied to calculate it the result is independent of our particular choice of @xmath160 . in the adiabatic approximation and the weak coupling limit , using eqs.([12.todas ] ) , we get the fourier transform of @xmath161 ,    @xmath162 ,    \nonumber \\   &   &      \label{13}\end{aligned}\ ] ]    where @xmath163 is given by the initial condition , and the resonant frequency @xmath164 is equal to    @xmath165    the first term of eq.([13 ] ) corresponds to the constant component of the magnetic field , and the second one shows the dissipation effects on this field component . the last term on the r.h.s . of the above equation displays a real frequency shift which contains the contribution of the geometrical phase , as can be seen from eq.([gamma ] ) , due the presence of the tracer @xmath143 , and a line broadening caused only by the dissipative evolution . these effects on the magnetization agree with the ones derived in reference @xcite where the path integral formalism was applied . the expressions for the other components of magnetization are analogous , but somewhat lengthy . in summary , we have presented an analytical solution of the adiabatic limit of a spin @xmath139 in a precessing magnetic field embedded in a weakly dissipative medium , introduced phenomenologically . we consider two distinct nonunitary contributions that were accounted for by a lindblad type superoperator in the diagonal frame . we are able to derive analytical expressions for the geometric and imaginary phases in both cases in the presence of a weak dissipation in the adiabatic limit without further approximations .    from eqs.([26 ] ) and ( [ 30.2 ] ) we get that the nature ( path dependent or time dependent ) of the imaginary phase acquired by @xmath112 depends on the mechanism that introduces the dissipation in the quantum system . in both cases that we have studied , the dissipation is present due to the two level system being in contact with a reservoir . the constant @xmath27 on the r.h.s . of eq.([26 ] ) is associated to the time rate of population and not due to the variation of any external parameter . consequently the quantum geometric phase for @xmath166 is not modified by a complex value . an analogous argument is valid to explain why the imaginary phase acquired by the entry @xmath167 in the dephasing process is not geometric , either . in this last model the adiabatic theorem continues to be true while it is not true anymore for the spin 1/2 coupled to the electromagnetic fields at thermal equilibrium . differently from ellinas _ et al . _ in reference@xcite we do not call this complex phase as berry s phase . we reserve the name of `` berry s phase '' only to phases ( real or imaginary ) that are path dependent . decoherence effects are present and their manifestation is the shortening of the three components of the bloch s vector . the fact that the dissipation effect causes the suppression of interference patterns due to the geometric phase is not a particular result for the chosen liouvillian . it is rather general , stemming from the fact that the dissipation mechanism is not related to the variation of any set of external periodic parameters . the authors are in debt with m.c . nemes and j.g . peixoto de faria for useful discussions . the authors thank the referee for bringing to their attention the interesting model discussed in section 2.2 . a.c.a.p . and m.t.t . thank e.v . corra silva for the careful reading of the manuscript . aguiar pinto thanks cnpq for financial support . m.t . thomaz thanks cnpq for partial financial support . the master equation ( [ 24 ] ) in the instantaneous basis of the hamiltonian for arbitrary value of @xmath11 is        @xmath169               + ( 1+\lambda^2)\big[\sigma_+ ( t ) \rho_i ( t ) \sigma_-(t ) +                    \sigma_- ( t ) \rho_i ( t ) \sigma_+(t)\big ] -                              \nonumber \\ \nonumber \\               & - & \lambda\sqrt{1-\lambda^2 }                \big[e^{i\omega t } \big ( \sigma_n(t ) \rho_i ( t ) \sigma_-(t )                 + \sigma_- ( t ) \rho_i ( t ) \sigma_n(t ) \big ) +                    \nonumber\\                \nonumber\\               & + & e^{-i\omega t } \big ( \sigma_+ ( t ) \rho_i ( t ) \sigma_n(t ) +                   \sigma_n ( t ) \rho_i ( t ) \sigma_+(t ) \big)\big]\big\ } \ ; -                      \nonumber\\                   \nonumber \\               & - & \frac{1}{2 } \big\ { \{\rho_i ( t ) , \lambda \sigma_n(t ) +                      \sqrt{1-\lambda^2}(e^{-i\omega t}\sigma_+(t ) +                      e^{i\omega t}\sigma_-(t)\ } + 2\lambda\big [                      \sigma_+ ( t ) \rho_i ( t ) \sigma_-(t ) -                         \nonumber \\                      \nonumber\\               & - & \sigma_- ( t ) \rho_i ( t ) \sigma_+(t)\big ] + \sqrt{1-\lambda^2 }               \big [ e^{-i\omega t } \big ( \sigma_n(t ) \rho_i ( t ) \sigma_+(t ) -               \sigma_+(t ) \rho_i ( t ) \sigma_n(t ) \big ) +                   \nonumber \\                 \nonumber\\              & + & e^{i\omega t } \big ( \sigma_-(t ) \rho_i ( t ) \sigma_n(t ) -              \sigma_n(t ) \rho_i ( t ) \sigma_-(t ) \big)\big]\big\ } .              \label{a1.2}\end{aligned}\ ] ]                @xmath173 +                          ( 1-\lambda^2)\big [                          e^{-2i\omega t } \sigma_+(t ) \rho_i ( t ) \sigma_+(t ) +                             \nonumber \\ \nonumber\\                      & + & e^{2i\omega t } \sigma_-(t ) \rho_i ( t ) \sigma_-(t ) +                          \sigma_+(t ) \rho_i ( t ) \sigma_-(t ) +                          \sigma_-(t ) \rho_i ( t ) \sigma_+(t ) \big ] - \rho_i(t ) . \label{a2.2}\end{aligned}\ ] ]                                          like the case of no dissipation , the off - diagonal elements of the density matrix are not uniquely defined unless we impose a condition on @xmath174 , @xmath175 and @xmath90 . in order to get eq.([8 ] ) , we impose the conditions : @xmath176 and @xmath177 . if we make another choice for the instantaneous eigenstates of hamiltonian , @xmath178 , @xmath175 and @xmath90 , the dynamics of the coefficient @xmath179 in the new basis is related to to the dynamics of @xmath180 ( from eqs . ( [ 23 ] ) ) by @xmath181 .
Abstract:

Answer: we consider the spin 1/2 model coupled to a slowly varying magnetic field in the presence of a weak damping represented by a lindblad - form operators . we show that berry s geometrical phase remains unaltered by the two dissipation mechanism considered . dissipation effects are twofold : a shrinking in the modulus of the bloch s vector , which characterizes coherence loss and a time dependent ( dissipation related ) precession angle . we show that the line broadening of the fourier transformation of the components of magnetization is only due to the presence of dissipation . pacs numbers : 03.65.bz , 03.65.-w , 05.30 , 42.50.lc    keywords : berry s phase , lindblad - form operator , decoherence