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: decoherence is the process which limits our ability to maintain pure quantum states , or their linear superpositions . it is the phenomenon by which the classical world appears from the quantum one @xcite . in more physical terms it is described as the rapid destruction of the phase relation between two , or more , quantum states of a system caused by the entanglement of these states with different states of the environment . the present widespread interest in decoherence is due to the fact that it is the main limiting factor for quantum information processing . we can store information , indeed , in two - level quantum systems , known as qubits , which can become entangled each other , but decoherence can destroy any quantum superposition , reducing the system to a mixture of states , and the stored information is lost . for this reason decoherence control is now becoming a rapidly expanding field of investigation .    in a series of previous papers @xcite we have faced the control of decoherence by actively modifying the system s dynamics through a feedback loop . this procedure turned out to be very effective , in principle @xcite , to slow down the decoherence of the only one experiment @xcite , up to present , in which the decoherence of a mesoscopic superposition was detected . the main limiting aspect of this procedure is connected with the need of a measurement . in order to do the feedback in the appropriate way , one has first to perform a measurement and then the result of this measurement can be used to operate the feedback . however , any physical measurement is subject to the limitation associated with a non - unit detection efficiency . we have shown @xcite that with detection efficiency approaching unity the quantum superposition of states stored in a cavity can be protected against decoherence for many decoherence times @xmath0 , where @xmath0 is defined as the cavity relaxation time divided by the average photon number @xcite . we wish now to face the problem of eliminating the measurement in controlling the decoherence . we show , here , how it is possible to inhibit decoherence through the application of suitable open - loop control techniques to the system of interest , that is , by using appropriately shaped time - varying control fields . to be more specific , decoherence can be inhibited by subjecting a system to a sequence of very frequent _ parity kicks _ , i.e. , pulses designed in such a way that their effect is equivalent to the application of the parity operator on the system .    the paper is organized as follows : in section ii the parity kicks method is presented in its generality and it is shown how decoherence and dissipation are completely frozen in the limit of infinitely frequent pulses . in section iii the method is applied to the case of a damped harmonic oscillator , as for example , a given normal mode of a system of trapped ions . in section iv the numerical results corresponding to this case are presented showing that a considerable decoherence suppression is obtained when the parity kicks repetition rate becomes comparable to the typical timescale of the environment . in section v the possibility of applying this scheme to harness the decoherence of the center - of - mass motion in ion traps are discussed . let us consider a generic open system , described by the hamiltonian @xmath1 where @xmath2 is the bare system hamiltonian , @xmath3 denotes the reservoir hamiltonian and @xmath4 the hamiltonian describing the interaction between the system of interest and the reservoir , which is responsible for dissipation and decoherence . we shall now show that if the hamiltonian ( [ htot ] ) possesses appropriate symmetery properties with respect to parity , it is possible to actively control dissipation ( and the ensuing decoherence ) by adding suitably tailored time - dependent external forcing acting on system variables only . the new hamiltonian becomes @xmath5 where @xmath6 ( @xmath7 is the usual step function ) is a time - dependent , periodic , _ system operator _ with period @xmath8 describing a train of pulses of duration @xmath9 separated by the time interval @xmath10 . the stroboscopic time evolution is therefore given by the following evolution operator @xmath11^{n } \label{evolu}\ ] ] where @xmath12 describes the evolution during the external pulse and @xmath13 gives the standard evolution between the pulses . we now assume that the external pulse is so strong that it possible to neglect the standard evolution during the pulse , @xmath14 , and then we assume that the pulse hamiltonian @xmath15 and the pulse width @xmath9 can be chosen so as to satisfy the following _ parity kick condition _ @xmath16 where @xmath17 is the system parity operator . it is now possible to see that such a time - dependent modification of the system dynamics is able to perfectly protect the system dynamics and completely inhibit decoherence whenever the following general conditions are satisfied : @xmath18 that is , the system hamiltonian is parity invariant and the interaction with the external environment anticommutes with the system parity operator . to be more specific , one has that in the ideal limit of continuous parity kicks , that is @xmath19 the pulsed perturbation is able to eliminate _ completely _ the interaction with the environment and therefore all the physical phenomena associated with it , i.e. , energy dissipation , diffusion and decoherence . to see this it is sufficient to consider the evolution operator during two successive parity kicks , which , using the parity kick condition of eq . ( [ parki ] ) , can be written as @xmath20 since @xmath21 one has that the time evolution after two successive kicks is driven by the unitary operator @xmath22 this expression clearly shows how the pulsed perturbation is able to `` freeze '' the dissipative interaction with the environment : the application of two successive parity kicks alternatively changes the sign of the interaction hamiltonian between system and reservoir . therefore one expects that , in the limit of continuous kicks , this sign inversion becomes infinitely fast and the interaction with the environment averages exactly to zero . this fact can be easily shown just using the definition of continuous kicks limit . in fact , in this limit @xmath23^{\frac{t}{2t+2\tau_{0 } } }   \ ; , \label{proof}\ ] ] which , using just the definition of the exponential operator , yields @xmath24 this means that in the ideal limit of continuous parity kicks , the interaction with the environment is completely eliminated and only the free uncoupled evolution is left . let us briefly discuss the physical interpretation of this result . the continuous kicks limit ( [ conti ] ) is formally analogous to the continuous measurement limit usually considered in the quantum zeno effect ( see for example @xcite ) in which a stimulated two - level transition is inhibited by a sufficiently frequent sequence of laser pulses . however this is only a mathematical analogy because at the physical level one has two opposite situations . in fact , during the pulses , in the quantum zeno effect the interaction with the environment ( i.e. , the measurement apparatus ) prevails over the internal dynamics , while in the present situation the interaction with the reservoir is practically turned off by the externally controlled internal dynamics ( see eq . ( [ parki ] ) ) . this parity kick idea has instead a strong relationship with spin echoes phenomena @xcite in which the application of appropriate rf - pulses is able to eliminate much of the dephasing in nuclear magnetic resonance spectroscopy experiments . these clever rf - pulses realize a sort of time - reversal and a similar situation takes place here in the case of parity kicks . in fact the relevant evolution operator is given by the unitary operator governing the time evolution after two pulses @xmath25 of eq . ( [ 4kiks ] ) , in which the standard evolution during a time interval @xmath10 is followed by an evolution in which the interaction with the environment is time - reversed for the time @xmath10 . if the time interval @xmath10 is comparable to the timescale at which dissipative phenomena take place , this time - reversal is too late to yield appreciable effects , while if @xmath10 is sufficiently small the very frequent reversal of the interaction with the reservoir may give an effective freezing of any dissipative phenomena . two very recent papers @xcite have considered a generalization of this parity kicks method and have showed that a complete decoupling between system and environment can be obtained for a generic system if one considers an appropriate sequence of infinitely frequent kicks , realizing a symmetrization of the evolution with respect to a given group . the present parity kicks method is in fact equivalent to symmetrize the time evolution of the open system with respect to the group @xmath26 , composed by the identity and the system parity operator @xmath17 . here we focus only on this case because the experimental realization of the parity kicks of eq . ( [ parki ] ) is quite easy in many cases and moreover the applicability conditions of the present method , eqs . ( [ condit1 ] ) and ( [ condit2 ] ) , are satisfied by many interesting physical systems . the continuous kicks limit is only a mathematical idealization of no practical interest . however it shows that there is in principle no limitation in the decoherence suppression one can achieve using parity kicks . moreover it clearly shows that one has to use the most possible frequent pulses in order to achieve a significative inhibition of decoherence . therefore the relevant question from a physical point of view is to determine at which values of the period @xmath8 one begins to have a significant suppression of dissipation and decoherence . we shall answer this question by considering in the next section , a specific example of experimental interest , that is , a damped harmonic oscillator , representing for instance a normal mode of a system of trapped ions , or an electromagnetic mode in a cavity . let us consider a harmonic oscillator with bare hamiltonian @xmath27 describing for example a given normal mode of a system of ions in a paul trap or an electromagnetic mode in a cavity . in these two cases the relevant environmental degrees of freedom can be described in terms of a collection of independent bosonic modes @xcite @xmath28 representing the elementary excitations of the environment ( in the cavity mode case they are simply the vacuum electromagnetic modes ) . moreover , the interaction with the environment is usually well described by the following bilinear term in which the `` counter - rotating '' terms are dropped @xcite @xmath29 in these cases time evolution is usually described in the frame rotating at the bare oscillation frequency @xmath30 in which the effective total hamiltonian of eq . ( [ htot ] ) becomes @xmath31 where @xmath32 it is immediate to see that in this example the conditions of eqs . ( [ condit1 ] ) and ( [ condit2 ] ) for the application of the parity kicks method are satisfied . more generally , conditions ( [ condit1 ] ) and ( [ condit2 ] ) are satisfied whenever @xmath2 is an even function of @xmath33 and @xmath34 and @xmath4 is an odd function of @xmath33 and @xmath34 . in the harmonic oscillator case the pulsed perturbation realizing the sequence of parity kicks can be simply obtained as a train of @xmath35-phase shifts , that is @xmath36 where @xmath37 and the pulse height @xmath38 and the pulse width @xmath9 satisfy the condition @xmath39    to see how large the pulses repetition rate has to be in order to achieve a relevant decoherence suppression , first of all one has to compare the pulsing period @xmath8 with the relevant timescales of the problem . the harmonic oscillator dynamics is characterized by the free oscillation frequency @xmath30 and by its energy decay rate @xmath40 ( which is a function of the couplings @xmath41 ) ; the bath is instead characterized by its ultraviolet frequency cutoff @xmath42 which essentially fixes the response time of the reservoir and generally depends on the system and bath considered . a reservoir is usually much faster than the system of interest and this means that one usually has @xmath43 . typically the reduced dynamics of the system of interest is described in terms of effective master equations which are derived using the markovian approximation ( see for example @xcite ) which means assuming the limit @xmath44 ; however , the existence of a finite cutoff @xmath42 is always demanded on physical grounds @xcite and this parameter corresponds for example to the debye frequency in the case of a phonon bath .    from the preceding section , it is clear that the pulsing period @xmath8 has to be much smaller than @xmath45 , otherwise the change of sign of the interaction hamiltonian is realized when a significant transfer of energy and quantum coherence from the system into the environment has already taken place . at the same time it is easy to understand that the condition @xmath46 is a sufficient condition for a significant suppression of dissipation , since in this case the sign of the interaction @xmath4 changes with a rate much faster than every frequency of the bath oscillators ; each bath degree of freedom becomes essentially decoupled from the system oscillator and there is no significant energy exchange , i.e. no dissipation . therefore the relevant question is : for which value of the pulsing period @xmath8 within the range @xmath47 $ ] one begins to have a relevant inhibition of dissipation ? we shall answer this question by studying in particular the time evolution of a schrdinger cat state , that is , a linear superposition of two coherent states of the oscillator of interest @xmath48 where @xmath49 we consider this particular example because these states are the paradigmatic quantum states in which the progressive effects of decoherence and dissipation caused by the environment are well distinct and clearly visible @xcite . from the general proof of the preceding section it is evident that whenever one finds a significant suppression of decoherence in the schrdinger cat case , this implies that the system - reservoir interaction is essentially averaged to zero and that one gets a significant system - environment decoupling in general . describing the evolution of such a superposition state in the presence of the dissipative interaction with a reservoir of oscillators which is initially at thermal equilibrium at @xmath50 is quite simple . in fact , it is possible to use the fact that a tensor product of coherent states retains its form at all times when the evolution is generated by the hamiltonian ( [ hrwa ] ) , that is @xmath51 where the time - dependent coherent amplitudes satisfy the following set of linear differential equations @xmath52 therefore , using eq . ( [ evocohe ] ) , one has the following time evolution for the schrdinger cat state ( [ cat ] ) @xmath53 the corresponding reduced density matrix of the oscillator of interest at time @xmath54 , @xmath55 , is given by @xmath56 \right\ } \ ; , \label{matocat}\end{aligned}\ ] ] i.e. is completely determined by the coherent amplitude @xmath57 which is a decaying function describing dissipation of the oscillator s energy and by the function @xmath58 describing the suppression of quantum interference terms , whose explicit expression is @xmath59    these results for the time evolution of the schrdinger cat are valid both in the presence and in the absence of the external pulsed driving . in fact the only difference lies in the fact that in the absence of kicks one has the standard evolution driven by @xmath60 , where @xmath61 is given by ( [ hrwa ] ) , while in the presence of parity kicks one has a stroboscopic - like evolution driven by a unitary operator analogous to that of eq . ( [ 4kiks ] ) @xmath62^{n}\;. \label{rwakiks}\ ] ] ( we have neglected only the system - reservoir interaction during the parity kick and therefore we have added two terms @xmath63 in eq . ( [ rwakiks ] ) with respect to eq . ( [ 4kiks ] ) ) .    due to the general result ( [ evocohe ] ) , it is convenient to express the state of the whole system in terms of a vector @xmath64 whose zero - th component is given by the amplitude @xmath65 and whose k - th component is given by the amplitude @xmath66 , @xmath67 . in this way the formal solution of the set of linear equations ( [ eqlin ] ) can be expressed as @xmath68 where @xmath69 is a matrix whose first matrix element @xmath70 is given by the following inverse laplace transform @xmath71 } \ ; , \label{azz}\ ] ] where @xmath72 all the other matrix elements can be expressed in terms of this matrix element @xmath70 in the following way : @xmath73 @xmath74    in the presence of kicks , the stroboscopic evolution is determined by the consecutive application of the evolution operator for the elementary cycle lasting the time interval @xmath75 ( see eq . ( [ rwakiks ] ) ) , during which one has a standard evolution for a time @xmath10 , followed by an uncoupled evolution for a time @xmath9 , then a `` reversed '' evolution for a time @xmath10 in which the interaction with the environment has the opposite sign and another uncoupled evolution for the time @xmath9 at the end . it is easy to understand that in terms of the vector of coherent amplitudes , the evolution in the presence of parity kicks can be described as @xmath76^{n }          \left(\matrix {          \alpha_{0 } \cr          \vdots \cr          0   \cr          \vdots \cr}\right )   \ ; , \label{matkik}\ ] ] where @xmath77 is the diagonal matrix associated to the free evolution driven by @xmath78 and is given by @xmath79 . we have seen that the quantities of interest are ( see eq . ( [ matocat ] ) ) @xmath65 and @xmath58 ; for the amplitude of the two coherent states one has @xmath80 in the absence of kicks and @xmath81^{n}\right\}_{00 } \alpha_{0 }      \label{aki}\ ] ] in the presence of kicks ( @xmath82 ) . as for the function @xmath58 , it is clear that each amplitude @xmath66 in ( [ decod ] ) is proportional to @xmath83 and therefore one can write @xmath84 where @xmath85 in the absence of kicks and @xmath86^{n}\right\}_{k0 } |^{2 }       \label{etaki}\ ] ] in the presence of kicks ( @xmath82 ) . in the description of dissipative phenomena one always considers a continuum distribution of oscillator frequencies in order to obtain an irreversible transfer of energy from the system of interest into the reservoir . moreover , most often , also the markovian assumption is made which means assuming an infinitely fast bath with an infinite frequency cutoff @xmath42 . this case of a standard vacuum bath in the markovian limit is characterized by an infinite , continuous and flat distribution of couplings @xcite @xmath87 which implies @xmath88 using eqs . ( [ anoki ] ) , ( [ azk ] ) and ( [ etanoki ] ) , we get @xmath89 the last sum in the expression for @xmath90 has to be evaluated in the continuum limit , that is , replacing the sum with an integral over the whole real @xmath91-axis and using @xmath92 . the result is the standard vacuum bath expression for the decoherence function @xmath93 @xcite @xmath94    from the general expressions of the above section it is clear that it is not possible to solve the dynamics in the presence of parity kicks in simple analytical form . we are therefore forced to solve numerically the problem , by simulating the continuous distribution of bath oscillators with a large but finite number of oscillators with closely spaced frequencies . to be more specific , we have considered a bath of @xmath95 oscillators , with equally spaced frequencies , symmetrically distributed around the resonance frequency @xmath30 , i.e. @xmath96 and we have considered a constant distribution of couplings similar to that associated with the markovian limit @xmath97    considering a finite number of bath oscillators has two main effects with respect to the standard case of a continuous markovian bath . first of all , the adoption of a discrete frequency distribution with a fixed spacing @xmath98 implies that all the dynamical quantities are periodic with period @xmath99 @xcite which can be therefore considered a sort of `` revival time '' . it is clear that our numerical solution will correctly describe dissipative phenomena provided that we focus only on not too large times ( say @xmath100 ) . moreover , the introduction of a finite cutoff ( @xmath101 in our case ) implies a modification of the coupling spectrum @xmath102 at very high frequency with respect to the infinitely flat distribution of the markovian treatment ( see eq . ( [ fla ] ) ) . this fact will manifest itself in a slight modification of the exponential behavior shown by eqs . ( [ mar1 ] ) and ( [ decomar ] ) at very short times ( @xmath103 ) @xcite . we have verified both facts in our simulations . however , to facilitate the comparison between the dynamics in the presence of parity kicks with the standard case of a markovian bath , we have chosen the parameters of our simulation so that , within the time interval of interest , @xmath104 say , we found no appreciable difference between the standard markovian bath expressions ( [ mar1 ] ) and ( [ decomar ] ) and the corresponding general expression for a discrete distribution of oscillators ( [ anoki ] ) and ( [ etanoki ] ) . let us now see what is the effect of parity kicks on decoherence , by studying first of all the behavior of the decoherence function @xmath93 for different values of the pulsing period @xmath8 . as discussed above , we expect an increasing decoherence suppression for decreasing values of @xmath8 and this is actually confirmed by fig .  1 , showing @xmath93 for different values of @xmath10 ( we have chosen @xmath105 for simplicity in the numerical simulation ) . in fact the increase of @xmath93 is monotonically slowed down as the time interval between two successive kicks @xmath10 becomes smaller and smaller . moreover fig .  1 seems to suggest that decoherence inhibition becomes very efficient ( @xmath106 ) when the kick frequency @xmath107 becomes comparable to the cutoff frequency @xmath42 . to better clarify this important point we have plotted in fig . 2 the value of the decoherence function at a fixed time as a function of the time between two successive kicks @xmath10 . in fig .  2a the decoherence function at half relaxation time @xmath108 is plotted as a function of @xmath109 , while in fig . 2b the decoherence function at one relaxation time @xmath110 is plotted again as a function of @xmath109 . in both cases one sees a quite sharp transition at @xmath111 : _ decoherence is almost completely inhibited as soon as the kick frequency @xmath112 becomes larger than the cutoff frequency @xmath113_. therefore we can conclude that by forcing the dynamics of an harmonic oscillator with appropriate parity kicks one is able to inhibit decoherence almost completely . this decoherence suppression becomes significant when the kick frequency @xmath107 becomes of the order of the typical reservoir timescale , i.e. the cutoff frequency @xmath113 . a similar conclusion has been reached by viola and lloyd @xcite , who have considered a model for decoherence suppression in a spin - boson model very similar in spirit to that presented here . in their paper they have considered a single 1/2-spin system coupled to an environment and they have shown how decoherence can be suppressed by a sequence of appropriately shaped pulses . they have shown that the 1/2-spin decoherence can be almost completely suppressed provided that the pulses repetition rate is at least comparable with the environment frequency cutoff . similar results have been obtained in a more recent paper @xcite in which a generalization to more qubits and more general sequences of pulses is considered . however , despite the similarity of our approach to decoherence control to that of these papers , there are important differences between the present paper and refs .  @xcite . first of all , we have considered a harmonic oscillator instead of a @xmath114-spin system , but above all , we have considered a _ dissipative _ bath , that is a zero temperature reservoir inducing not only decoherence ( i.e. quantum information decay ) , but also dissipation ( i.e. energy decay ) . in ref . @xcite on the contrary , there is only decoherence and the @xmath114-spin system energy is conserved . the more general nature of the present model helps in clarifying a main point of this impulsive method to combat decoherence : parity kicks do not simply suppress decoherence but tend to completely inhibit _ any interaction between system and environment_. this means that also energy dissipation is suppressed ; in more intuitive terms , the external pulsed driving is perfectly `` phase matched '' to the system dynamics , so that any transfer of energy and quantum coherence between system and bath is inhibited . the fact that decoherence suppression in our model is always associated with the suppression of dissipation can be easily shown using the fact that the evolution of the whole system is unitary . in particular this implies that the matrix @xmath69 , driving the evolution in absence of kicks ( see eq . ( [ matq ] ) ) , and the matrix @xmath115^{n } $ ] , driving the evolution in the presence of kicks ( see eq . ( [ matkik ] ) ) , are unitary matrices and therefore subject to the condition @xmath116 where @xmath117 denotes any of the two above mentioned matrices . now , using the reduced density matrix of eq . ( [ matocat ] ) , it is easy to derive the mean oscillator energy @xmath118 from eqs . ( [ anoki ] ) and ( [ aki ] ) one therefore derives that the normalized oscillator mean energy is given by @xmath119 in the absence of kicks and @xmath120^{n}\right\}_{00 } \right|^{2 }       \label{normeaki}\ ] ] in the presence of kicks ( @xmath82 ) . now , by considering eqs . ( [ etanoki ] ) and ( [ etaki ] ) and the unitarity condition ( [ unita ] ) , it is immediate to get the following simple relation between decoherence and dissipation @xmath121 which is valid both in the presence and in the absence of kicks . this equation simply shows that when decoherence is suppressed ( @xmath122 ) the oscillator energy is conserved .    a qualitative demonstration of the ability of the kick method to suppress dissipation and decoherence is provided by fig .  3 . ( a ) shows the wigner function of an initial odd cat state with @xmath123 , @xmath124 ; ( b ) shows the wigner function of the same cat state evolved for a time @xmath125 in the presence of parity kicks with @xmath126 and ( c ) the wigner function of the same state again after a time @xmath127 , but evolved in absence of kicks . this elapsed time is ten times the decoherence time of the schrdinger cat state , @xmath128 , the state is almost indistinguishable from the initial one and that the quantum wiggles of the wigner function are still well visible . the numerical results presented above provides a clear indication that perturbing the dynamics of an oscillator with an appropriate periodic pulsing can be a highly efficient method for controlling decoherence . however we have seen that a significant decoherence suppression is obtained only when the pulsing frequency @xmath129 becomes comparable to the cutoff frequency of the reservoir . this fact poses some limitations on the experimental applicability of the proposed method . for example the parity kick method is certainly unfeasible , at least with the present technologies , in the case of optical modes in cavities . in this case , realizing a parity kick is not difficult in principle since impulsive phase shifts can be obtained using electro - optical modulators or a dispersive interaction between the optical mode and a fast crossing atom @xcite . the problem here is that it is practically impossible to make these kicks sufficiently frequent . in fact , we can assume , optimistically , that the frequency cutoff @xmath42 , even though larger , is of the order of the cavity mode frequency @xmath30 ( this is reasonable , since the relevant bath modes are those nearly resonant with the frequency @xmath30 ) . this means @xmath130 hz and it is evident such high values for @xmath129 are unrealistic both for electro - optical modulators and for dispersive atom - cavity mode interactions . a similar situation holds for an electromagnetic mode in a microwave cavity as that studied in the schrdinger cat experiment of brune _ et al . _ @xcite . in this experiment the dispersive interaction between a rydberg atom and the microwave mode is used for the generation and the detection of the cat state and therefore the realization of @xmath35-phase shifts is easy . the problem again is to have sufficiently fast and frequent @xmath35-phase shifts : in this case the mode frequency is @xmath131 ghz and therefore one would have a good decoherence suppression for @xmath132 sec . this is practically impossible because it implies quasi - relativistic velocities for the crossing rydberg atom and an unrealistic very high dispersive interaction in order to have a @xmath35-phase shift . the situation is instead different in the case of ions trapped in harmonic traps . in this case in fact , the free oscillation frequency @xmath30 , and therefore the cutoff frequency @xmath42 too , are usually much smaller and it becomes feasible to realize fast and frequent @xmath35-phase shifts . let us consider for example the case in which the oscillator mode of the preceding sections is the center - of - mass motion of a collection of trapped ions in a paul trap , as that considered in the experiments at nist in boulder @xcite . in this case the free oscillation frequency is of the order of @xmath133 mhz and a frequent sequence of parity kicks could be obtained in principle by appropriately pulsing at about 10 mhz the static potential applied to the end segments of the rods confining the ions along the z - axis @xcite . the duration @xmath9 and the intensity @xmath38 of the pulses have to be tailored so to have the parity kick condition @xmath134 ( see eq . ( [ parity ] ) ) and this implies having a very good control of the pulse area .    controlling the decoherence of the center - of - mass mode is crucial for the realization of quantum imformation processing with trapped ions . in fact , even though only the higher frequency vibrational modes will be used for quantum gate transitions involving the ions motion ( as it has been already done in @xcite where the deterministic entanglement of the internal states of two ions has been achieved by manipulating the stretching mode ) , suppressing decoherence of the center - of - mass motion is still important because the heating of the center - of - mass motion partially couples also with the other vibrational stretching modes @xcite . the parity kicks method could be used even in the case in which the center - of - mass mode itself is used as quantum bus for the realization of quantum gates . however in this case , the use of parity kicks can not be applied to protect against decoherence all kinds of quantum gates . in fact the parity kicks tends to average to zero any term in the hamiltonian which is an odd function of the bosonic operators @xmath33 and @xmath34 , and this means that the parity kicks could average to zero just the system dynamics we want to protect from decoherence , as for example , some gate operations in an ion - trap quantum computer . for example , it is easy to see that all the gate operations involving _ first _ red or blue sideband transitions , which implies having _ linear _ terms in @xmath33 and @xmath34 in the system hamiltonian , as for example the cirac - zoller controlled - not ( c - not ) gate @xcite and the c - not gate experimentally realized at nist by monroe _ et al . _ @xcite , tend to be averaged to zero by frequent parity kicks . to overcome this problem , it is however sufficient to restrict to quantum gates based on carrier transitions , which involve functions of @xmath135 only , and are not therefore affected by parity kicks , as for example the one - pulse gate proposed by monroe _ et al _ @xcite . therefore , the parity kick method presented here could be very useful to achieve a significant decoherence control in ion traps designed for quantum information processing . however we have to notice that our considerations are based on the model of section iii which is the one commonly used to describe dissipation , even though the actual sources of decoherence in ion traps have not been completely identified yet @xcite . for example refs . @xcite have considered the effect of fluctuations in the various experimental parameters , yielding no appreciable energy dissipation in the center of mass motion but only off - diagonal dephasing . this work has been partially supported by infm ( through the 1997 advanced research project `` cat '' ) and by the european union in the framework of the tmr network `` microlasers and cavity qed '' . we acknowledge useful discussions with g.j . milburn , s. schneider and q.a . turchette . 99 w. h. zurek , phys . d * 24 * , 1516 ( 1981 ) ; _ ibidem _ * 26 * , 1862 ( 1982 ) ; phys . today * 44 * ( 10 ) , 36 ( 1991 ) . p. tombesi and d. vitali , phys . a * 51 * , 4913 ( 1995 ) . p. tombesi and d. vitali , appl . b * 60 * , s69 ( 1995 ) ; p. goetsch , p. tombesi and d. vitali , phys . a * 54 * , 4519 ( 1996 ) . d. vitali , p. tombesi and g.j . milburn , phys . rev . lett . * 79 * , 2442 ( 1997 ) . d. vitali , p. tombesi and g.j . milburn , phys . a * 57 * , 4930 ( 1998 ) . m. brune , e. hagley , j. dreyer , x. maitre , a. maali , c. wunderlich , j.m . raimond and s. haroche , lett . * 77 * , 4887 ( 1996 ) . a.o . caldeira and a.j . leggett , phys . rev . a * 31 * , 1059 ( 1985 ) ; d.f . walls and g.j . milburn , phys . a * 31 * , 2403 ( 1985 ) . r. cook , phys . scripta * t21 * , 49 ( 1988 ) . itano , d.j . heinzen , j.j . bollinger and d. wineland , phys . a * 41 * , 2295 ( 1990 ) . hahn , phys . rev . * 80 * , 580 ( 1950 ) . p. zanardi , e - print quant - ph/9809064 ; l. viola , e. knill and s. lloyd , e - print quant - ph/9809071 . a.o . caldeira and a.j . leggett , ann . phys . ( n.y . ) * 149 * , 374 ( 1983 ) . gardiner , _ quantum noise _ , ( springer , berlin , 1991 ) . wineland , c. monroe , w.m . itano , d. leibfried , b.e . king , d.m . meekhof , j. res . technol . * 103 * , 259 ( 1998 ) . h. grabert and u. weiss , z. phys . b * 55 * , 87 ( 1984 ) ; f. haake and r. reibold , phys . a * 32 * , 2462 ( 1985 ) ; b.l . paz and y. zhang , phys . d * 45 * , 2843 ( 1992 ) . cukier , k.e . shuler and j.d . weeks , j. stat . phys . * 5 * , 99 ( 1972 ) ; j.l . gruver , j. aliaga , h.a . cerdeira , a.n . proto , phys . e * 51 * , 6263 ( 1995 ) . l. viola and s. lloyd , phys . a * 58 * , 2733 ( 1998 ) . duan , g - c . guo , e - print quant - ph/9807072 . turchette , c.j . hood , w. lange , h. mabuchi and h.j . kimble , phys . lett . * 75 * , 4710 ( 1995 ) . turchette , c.s . wood , b.e . king , c.j . myatt , d. liebfried , w.m . itano , d.j . wineland , phys . lett . , * 81 * , 3631 ( 1998 ) . cirac , p. zoller , phys . rev * 74 * , 4091 ( 1995 ) . c. monroe , d.m . meekhof , b.e . king , w.m . itano , d.j . wineland , phys . rev . lett . * 75 * , 4714 ( 1995 ) . c. monroe , d. leibfried , b.e . king , d.m . meekhof , w.m . itano , d.j . wineland , phys . a. * 55 * , r2489 ( 1997 ) . s. schneider , g.j . milburn , phys . rev . a * 57 * , 3748 ( 1998 ) . m. murao , p.l . knight , phys . a * 58 * , 663 ( 1998 ) .
Abstract:

Answer: we show how it is possible to suppress decoherence using tailored external forcing acting as pulses . in the limit of infinitely frequent pulses decoherence and dissipation are completely frozen ; however , a significant decoherence suppression is already obtained when the frequency of the pulses is of the order of the reservoir typical frequency scale . this method could be useful in particular to suppress the decoherence of the center - of - mass motion in ion traps .


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: extreme events such as heat or cold waves have major impacts on human health and ecological sustainability . motivated by the pressing need for better risk assessment , inference for spatial extremes has developed rapidly during the last decade ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? max - stable processes are the only possible non - degenerate limits for rescaled pointwise maxima of random processes , thus motivating their use in modeling spatial extremes . flexible max - stable models have been developed for inference on spatial extremes . over a variety of applications , the brown  resnick process  @xcite has proven to be a valuable model for extremes of environmental variables , such as rainfall , temperature , wind speed and river discharge ( see , e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? it is built using gaussian processes and its dependence structure is determined by a variogram , so well - established ideas from classical geostatistics may be exported to the extremal context . the greater challenge is inference , since the full likelihood of the brown  resnick process is tractable only in relatively low - dimensional settings @xcite . composite likelihood techniques have been successfully used to fit max - stable models to block maxima or threshold exceedances , but with reduced estimation efficiency . more efficient procedures based on full , possibly censored , likelihoods @xcite have been developed for threshold exceedances , allowing the construction of more elaborate models for spatial extremes . bayesian hierarchical models are appealing in spatial statistics because of the flexibility with which they can track the spatial variation of marginal distributions . @xcite , @xcite and @xcite constructed hierarchical models for extremes assuming conditional independence of extreme occurrences . such models may have flexible margins but can not produce realistic simulations of extreme events @xcite . @xcite and @xcite used hierarchical models that model spatial dependence using gaussian processes , a somewhat restrictive class of models that are asymptotically independent and thus may underestimate potential dependence at the most extreme levels @xcite . @xcite and @xcite proposed the use of composite likelihoods in markov chain monte carlo ( mcmc ) algorithms and applied them for bayesian inference on max - stable models , but these methods present computational challenges and do not yield exact simulation from the posterior distribution . @xcite constructed a max - stable spatial model based on a hierarchical representation of logistic extreme - value distributions @xcite .    in this paper we propose a bayesian hierarchical model based on the brown  resnick process , which exploits a special case for which the full likelihood of such a process can be calculated , using the joint distribution of componentwise maxima and their partitions in terms of individual events @xcite . we discuss the difficulties of using the partitions in practice , and we propose a new bayesian approach to fitting max - stable processes without fixed partitions . our work is motivated by the changing climate of northern fennoscandia . in northern regions of norway , finland , russia and sweden , extreme low temperatures play a central role in ecological sustainability by limiting the extent of outbreaks of forest pests such as the autumnal moth _ epirrita autumnata _ , whose larvae have caused serious damage by defoliating mountain birch . forests in the region have mostly been protected from this pest by extreme low winter temperatures , since its eggs can not survive below @xmath0c @xcite . winter minimum temperatures are often below this threshold , but in recent years they have increased and projections suggest they will increase even more in the future ( * ? ? ? *  e.1 ) . although temperature is not the only factor limiting the spread of the moth , it is one of the most important , and thus plays a central role in the ecology of northern fennoscandia . in  [ sec : forecasts ] we discuss the possible changes in winter minimum temperatures in the region , focusing on the probability that they will be less than @xmath0c . motivated by recent disastrous heat wave events over the world , max - stable processes and related methods have been used to model temperature extremes and their potential changes in time . @xcite used the so - called @xcite process to model annual maximum temperatures in switzerland . they used latitude , longitude and elevation to estimate response surfaces for the marginal parameters of extremes . @xcite describe a bayesian spatial model for annual maximum temperatures in the us based on mixtures of gaussian distributions . @xcite used a bayesian max - stable model to model extreme temperature data over europe . section  [ sec : setting ] describes the data we study . section  [ sec : evt ] summarizes some extreme value ideas and introduces the inference methods for the brown  resnick model . the hierarchical model and the application are discussed in section  [ sec : bhmodel ] . we use temperature data from 20 meteorological stations located in northern regions of finland , norway and sweden ; see fig . [ fig : mapfinland ] . the smallest and largest distances between two stations are 26 and 471  km . daily minimum temperatures are available from 1960 to 2013 , but with data missing . in view of the motivation for our application , inference is needed only for the winter minima . we focus on modeling minimum temperatures in the months december  march and we use the daily data to investigate their co - occurrences . weather stations in the northern regions of finland , norway and sweden . the dashed lines shows the convex hull of the stations , which will be our region of interest for simulation.,scaledwidth=81.0% ]    missing daily data censor the winter minima . winter minimum values for the 20 stations are shown in fig . [ fig : annualmin ] . our dataset has at least 26 winter minima for each station , and each year has values of winter minimum for at least 15 of the 20 stations . c ) at the @xmath1 stations ( a different color for each station ) . the horizontal line shows the @xmath0c threshold . ] mean winter minimum temperatures depend on local geography and topography , and we shall model this non - stationarity in their marginal distributions using covariates and random effects , which will also allow us to extrapolate results to ungauged sites . we derived covariates that are known to influence minimum temperature in northern fennoscandia @xcite . we supplemented the coordinates of the stations with their absolute elevation above sea level and elevation relative to the surrounding areas ( calculated on a 2  km  @xmath2  2  km spatial moving average ) , a proximity index to the arctic ocean , calculated with a 160  km  @xmath2  160  km moving window assuming that the effect of the ocean reaches 80  km inland , and a lake cover index , which expresses the proportion of area covered by lakes in a 3  km  @xmath2  3  km window ( see * ? ? ? * ) .      in fig . [ fig : annualmin ] it is difficult to see the presence of a temporal trend with only 53 years of noisy data , but the figure suggests that minima may generally be slightly less extreme from 1990 onwards . the inclusion of a temporal effect in the model will help to assess this .      since we focus on modeling the winter minima , temporal dependence in the daily temperatures is not of interest in itself . but our approach requires us to partition minima in terms of independent generating events , which is complicated by the presence of temporal dependence in the daily data ; see   [ sec : declustering ] . cold waves stay over the region for several days , so minimum temperatures tend to be very dependent from day to day . preliminary analysis ( see the supplementary material , * ? ? ? * ) showed that the dependence is strong for the first few days and vanishes after approximately five to ten days at each location . we expect winter minimum temperatures to be dependent over large regions , and this is the case in our data . [ fig : annualmin ] suggests that the minima for the @xmath1 stations behave similarly for the @xmath3 years : the years 1966 and 1999 were very cold at all @xmath1 locations , whereas 1973 , 1993 and 2007 are warmer at all locations . it appears that most of the minima occurred on the same day or only a few days apart , so they may correspond to the same event , which suggests strong spatial dependence . results from a preliminary analysis ( see * ? ? ? * ) suggest that pairwise distributions are asymptotically dependent and max - stable . after accounting for differences in the marginal distributions , we will use a stationary isotropic max - stable model for modeling spatial dependence . let @xmath4 be a random process with continuous sample paths defined over a compact subset @xmath5 of @xmath6 , and let @xmath7 be independent replicates of @xmath8 . if there exist sequences of continuous functions @xmath9 and @xmath10 such that we have the distributional convergence @xmath11\right\}_{\bm s\in s } \rightarrow \{z(\bm s)\}_{\bm s\in s},\quad n\rightarrow\infty,\ ] ] where @xmath12 is a non - degenerate random process , we say that @xmath8 lies in the max - domain of attraction of the process @xmath13 . then the process @xmath13 must be max - stable @xcite and its marginals are of generalized extreme - value ( gev ) form , @xmath14_+^{-1/\xi(\bm s)}\right ) , \quad \bm s\in s,\ ] ] denoted @xmath15 , with @xmath16 , and real location @xmath17 , scale @xmath18 and shape @xmath19 parameters ; throughout the paper , the case @xmath20 is defined as the limit for @xmath21 . every max - stable process with marginal parameters @xmath17 , @xmath22 and @xmath19 can be written as @xmath23 where @xmath13 is a simple max - stable process with unit frchet margins , i.e. , @xmath24 . thus it suffices to consider simple max - stable processes , henceforth denoted by @xmath13 .    as a consequence of max - stability , for every set of locations @xmath25 , @xmath26 where the so - called exponent function @xmath27 is homogeneous of order @xmath28 , i.e. , @xmath29 , @xmath30 , and satisfies @xmath31 for any permutation of its arguments . in particular , @xmath32 where @xmath33 $ ] is known as the extremal coefficient of @xmath13 relative to @xmath34 . the extremal coefficient can be interpreted as the effective number of independent variables among @xmath35 and is useful for model checking , since it can be estimated non - parametrically .    as they appear as limits for rescaled pointwise maxima of random processes , max - stable processes are natural models for observed block maxima . moreover , the equality @xmath36 shows that such processes also provide natural models for block minima . several parametric max - stable models have been proposed . the @xcite process was the first to be considered but its realizations are too smooth for realistic modeling of complex extremes . the @xcite process is constructed from stationary gaussian processes and provides more realistic realizations controlled by different correlation functions . the extremal-@xmath37 process extends the schlather process by introducing an additional shape parameter @xcite . here we focus on the brown  resnick process @xcite which appears to be a good model for extremes of environmental processes . @xcite constructed a stationary max - stable process , called the brown  resnick process , as @xmath38 where the @xmath39 s are the points of a unit rate poisson process on @xmath40 and the @xmath41 are independent replicates of an intrinsically stationary centered gaussian process @xmath42 with variance function @xmath43 . the dependence structure of @xmath13 is inherited from that of @xmath42 , which is controlled by the variogram @xmath44 $ ] . without loss of generality , we suppose that @xmath45 with probability one , so @xmath46 . a large class of models may be obtained by choosing different variograms . later we use the stable variogram @xmath47 , with @xmath48 and @xmath49 $ ] . the exponent function of a brown  resnick process is @xcite @xmath50 where @xmath51 is the cumulative distribution function of a @xmath52-dimensional gaussian distribution with mean vector @xmath53 and covariance matrix @xmath54 , @xmath55 , @xmath56 , @xmath57 denotes the @xmath58 vector consisting of the @xmath59th column of @xmath60 without its @xmath59th component , @xmath61 denotes the matrix @xmath60 without its @xmath59th row and @xmath59th column , and so forth . let @xmath62 be a vector stemming from a simple max - stable process @xmath13 . differentiation of   shows that the full likelihood contribution for one realization @xmath63 of @xmath64 may be written as @xmath65 \exp\{-v(z_1,\ldots , z_d)\},\ ] ] where @xmath66 denotes the set of all partitions of @xmath67 and the subscript @xmath68 on @xmath69 denotes the partial derivative of the function @xmath27 with respect to the components in the subset @xmath68 . although the likelihood has a closed - form expression for the brown  resnick process @xcite , the size of @xmath66 equals the bell number and is larger than @xmath70 for @xmath71 , making full likelihood inference computationally unattainable in many situations . hence inference for max - stable processes has typically been based on composite likelihoods using lower - order densities @xcite . not only do composite likelihoods result in a loss of efficiency compared to full likelihood estimation , they also increase the computational burden when employed in a bayesian analysis , because the composite likelihood does not directly yield a posterior which appropriately accounts for parameter uncertainty @xcite . @xcite , @xcite and @xcite developed full likelihood inferences for threshold exceedances of processes in the max - domain of attraction of a brown  resnick process , which , while challenging , have lower computational burdens than does  .    for our likelihood , we employ the result of @xcite , who showed that , for the brown  resnick process , partial derivatives of @xmath27 with respect to the first @xmath72 indices , denoted @xmath73 , satisfy @xmath74,\end{gathered}\ ] ] with @xmath75 , @xmath76 , @xmath77 , @xmath78 the @xmath79 matrix consisting of the first @xmath72 rows and columns of @xmath80 , @xmath81 , @xmath82 , @xmath83 , @xmath84 , @xmath85 and @xmath86 , where @xmath87      @xcite derived a likelihood based on componentwise block maxima and their occurrence times that avoids direct differentiation of  . let @xmath88 denote independent replicates of a random vector @xmath89 . without loss of generality , suppose that the process @xmath8 has unit frchet margins and lies in the max - domain of attraction of a simple max - stable process @xmath13 . let @xmath90 denote the componentwise maxima , and let @xmath91 denote the unique partition of @xmath67 such that @xmath92 comes from the same individual observation , i.e. , for each @xmath93 , there exists @xmath94 such that @xmath95 . then @xcite showed that , as @xmath96 , the block maximum and its associated partition @xmath97 converge in distribution to @xmath98 , whose joint density is @xmath99 this suggests using   as a likelihood for inference on the joint distribution of block maxima and their occurrence times . this approach only requires the maximum and the partition to be observed , and thus represents a compromise between the block maximum and threshold exceedance approaches . the relative efficiency of stephenson  tawn and pairwise maximum likelihood estimators has been investigated by @xcite . because the partition contains information about the dependence , the stephenson  tawn estimator has smaller variance than that based only on the block maxima . however , for finite block length , the former may be more biased , because the model @xmath98 is misspecified for @xmath97 , i.e. , both for the maximum and the partition @xcite . use of the limiting distribution for @xmath97 represents a strong assumption that can lead to poor fit , whereas approaches based on the maximum @xmath100 alone are more robust . a further difficulty arises from the definition of the partition for serially dependent data . the partition in the stephenson  tawn model is defined in terms of independent generating events . while it is natural to define the partition in terms of the days of occurrence of the maxima when the daily observations are independent , its estimation is harder when data show temporal dependence . in  [ sec : declustering ] we propose a space - time declustering procedure to estimate the partition . although the use of the estimated fixed partition is theoretically appealing , it represents a practical challenge . first , the model is misspecified , which affects the estimation of the max - stable model by the introduction of potentially biased information . second , the partition must be estimated in a preliminary step , and inference based on the stephenson  tawn model is conditional on the fixed partition and ignores the uncertainty associated to its estimation . the partition is not always very well defined ( see  [ sec : declustering ] ) and its estimation is sensitive to measurement error . this motivates us to propose an approach to estimating the parameters of the brown  resnick process , without fixing the partition , which uses only the values of the maxima and can be applied when the occurrence dates are unavailable . let @xmath62 be a vector stemming from a brown  resnick process @xmath13 . to overcome the intractability of the likelihood , we propose to augment the data with the unknown partition @xmath101 that generates @xmath64 in terms of the functions @xmath102 in  , i.e. , @xmath103 is defined such that the points @xmath104 in the same set @xmath68 come from the same function @xmath105 . the joint distribution of @xmath98 is  . we impute the unknown partition @xmath101 given the observed value of @xmath64 in an mcmc algorithm , using ideas similar to @xcite , who considered the partition in the context of conditional simulation from max - stable processes . the partition is imputed from its conditional distribution given the data . from   we have @xmath106 however , the number of partitions of the set @xmath67 is the @xmath107-th bell number , so   is intractable for large @xmath107 . we use a gibbs sampler to update the partition . for a partition @xmath108 , we choose a location @xmath109 at random and let @xmath110 denote the restriction of @xmath111 to @xmath112 . the size of @xmath110 is @xmath113 , or @xmath114 if @xmath115 is itself a partitioning set of @xmath111 . given @xmath110 , the partition is updated by attributing the location @xmath109 either to an existing set of @xmath110 or to a new set , giving at most @xmath116 choices . the conditional probability for a new partition @xmath117 with @xmath118 is @xmath119 the discrete distribution   has support on a set of size at most @xmath116 , so a new partition @xmath117 can easily be simulated .    in practice , inference is based on independent replicates @xmath120 , and the mcmc algorithm independently updates the corresponding latent partitions @xmath121 . given the observed values of the @xmath122 and the simulated @xmath123 , the marginal and dependence parameters of the model are updated from the joint distribution of the @xmath124 , which is a product of terms of the form  . we first fitted gev distributions separately at each of the 20 stations ; see @xcite . the estimates for the scale and shape parameters were not significantly different from site to site , but variation in the location parameter estimates suggested use of a non - stationary gev model with spatially varying location parameter and constant scale and shape parameters . we then fitted non - stationary marginal models at the 20 locations by maximizing the independence likelihood , which assumes the independence of the temperatures from site to site . non - stationarity was modeled using trend surfaces defined from the covariates mentioned in   [ sec : dataset ] and temporal trends in the location and scale parameters . we used the r package spatialextremes @xcite to fit these models , which were compared using the takeuchi information criterion , an equivalent of the akaike information criterion for misspecified models . our best model had all six covariates for the location surface , constant scale and shape parameters , and a temporal trend in the location only . its fit was assessed using marginal qq - plots : departure from the @xmath125 confidence intervals was observed for five stations . this suggests that the trend surface model is not flexible enough to model the variability in the location parameters over the 20 stations . inclusion of additional covariates did not improve the fit , motivating the use of a more flexible approach using random effects in a bayesian hierarchical model . we consider two bayesian hierarchical models based on the brown  resnick process . the first uses the stephenson  tawn model and estimates the partitions in a preliminary step through space - time declustering ; see    [ sec : stmodel ] ,  [ sec : declustering ] . the second imputes the partitions in the mcmc algorithm ; see   [ sec : estpart ] . here we describe the part of the hierarchical model for the data conditional on the partitions , which is common to both approaches . let @xmath126 ( @xmath127 ) denote the observed minimum temperatures for the @xmath128 years of data , let @xmath34 denote the @xmath129 locations of the stations such that @xmath130 , let @xmath131 denote the covariate for the day of occurrence of the winter minimum @xmath132 , defined as the number of years from 1 january 2000 , and finally let @xmath133 denote the partitions of @xmath67 for each year . when data are missing , the dimension of @xmath134 is lower than @xmath107 and so is the size of the corresponding partition . the first level of the hierarchical model gives the probability density of the observed minima and the partitions @xmath135 , conditional on the marginal and dependence parameters @xmath136 , as @xmath137 \prod_{k=1}^{k_i } \left[-v_{\pi_{ik}}\{f_{i1}(-y_{i1}),\ldots , f_{id}(-y_{id})\ } \right ] \prod_{j=1}^d f_{ij}'(-y_{ij } )   \right ) , \end{gathered}\ ] ] where @xmath138 denotes the marginal transformation to the standard frchet scale for @xmath139 and the prime denotes derivative with respect to @xmath140 ; the negative signs for the @xmath132 s in   correspond to modeling negated winter minima . the second level of the hierarchical model specifies the models for the exponent measure @xmath27 and the gev parameters . we define a space - time model for @xmath141 with the covariates discussed in   [ sec : dataset ] , with a linear temporal trend , and with a gaussian prior distribution as @xmath142 where @xmath143 denotes the vector of the intercept and covariates at location @xmath144 , and @xmath145 denotes a gaussian process with mean @xmath146 , variance @xmath147 , and correlation function @xmath148 . the random effect allows the location parameters to vary smoothly over space , in addition to any effect of the covariates . we kept @xmath149 and @xmath150 constant over space and time because this appeared to suffice for a good fit . thus @xmath151 , @xmath152 , and @xmath153 for each pair of replicate and location indices @xmath154 , and the marginal gev distributions of the model at the observed data can be written as @xmath155 where @xmath156 is a gaussian vector with mean @xmath157 , where @xmath158 is a @xmath159 matrix with rows @xmath160 , and covariance matrix @xmath161 . we use a brown  resnick process with variogram @xmath162 as a model for the exponent function @xmath27 . prior distributions are assigned to the variables at the last level of a hierarchical model . we used vague priors for all parameters except those of the latent gaussian field ; see @xcite . the prior distributions for @xmath163 and @xmath147 were chosen to obtain conjugate full conditional distributions for these parameters . the use of the stephenson  tawn model involves choosing the partitions in  , whose sets should correspond to independent individual events . several approaches have been proposed to identify clusters in univariate time series ( e.g. ,  ( * ? ? ? * chapter  10 ) ; @xcite ) . in the multivariate case , declustering techniques have been proposed by @xcite and @xcite . here we follow @xcite , who identified clusters of storm surges by assuming that events separated by a fixed time lag are independent .    in our application , a difficulty arises from the use of the dates of occurrences . because of the measurement precision of the temperatures ( @xmath164c ) , the same annual minimum may occur on different days , affecting the estimation of the partition ; this is the case for 7% of the minima , with dates of occurrences separated by more than two days in 3% of cases . when a minimum has multiple days of occurrences , we choose one of these days at random to define our partition .    for our temperature data , preliminary analysis suggests choosing a time lag of five days : two observations at the same location and separated by more than five days are assumed to be independent . to identify clusters for the 20 time series , we use a single linkage ( or friend - to - friend ) method whereby a cluster is composed of data that are dependent in pairs . thus if the minima occurred at time 1 for the first series , time 4 for the second and time 7 for the third , these three minima are considered to lie in the same cluster . with this procedure , the size of the partition for each year is at most five , and winter temperature minima usually occurred from widespread events that lasted for between five and ten days .      the mcmc algorithm for our brown  resnick model is nonstandard because the likelihood   can not be calculated exactly but must be estimated . the likelihood involves the function @xmath27 and its partial derivatives , which themselves involve the gaussian cumulative distribution function in their expressions , and thus can only be estimated . we used the approach of @xcite to use a likelihood estimator in the mcmc algorithm . under the assumption that the likelihood estimator is unbiased , that algorithm generates a markov chain that has the targeted posterior distribution as stationary distribution . we estimated the likelihood   by replacing the gaussian cumulative distribution functions in its expression by monte carlo estimators . our estimator is not unbiased , but we constructed the estimator to minimize the potential bias in the estimation of the posterior and obtain good mixing properties for the markov chain in a reasonable time ; see @xcite . for our dataset and the fixed partitions estimated in   [ sec : declustering ] , a single estimation of the likelihood   takes about one second on a single cpu . we used a gibbs sampler with metropolis  hastings steps to update the marginal and dependence parameters of the model . to minimize the number of likelihood estimations in the mcmc algorithm , and because some parameters were found to be very dependent , we updated the parameters of the model in blocks ; see @xcite . with our choice of blocks , one iteration of the mcmc algorithm requires four likelihood estimations for the stephenson  tawn model . for the random partition model , the gibbs sampler for the partitions requires estimating the probabilities  , which increases the computational burden a bit more . since one iteration of the mcmc algorithm takes several seconds on a single cpu , we used parallel computing , fitting 50 independent markov chains each of length 15,000 and discarding the first 5000 iterations of each chain . the resulting chains were merged to obtain a final chain of length 500,000 . the mcmc diagnostics are discussed in @xcite . the mean of the posterior distributions were used as pointwise estimators . credible intervals were calculated from the quantiles of the markov chains . the results of a simulation study suggest that the mcmc algorithm provides sensible inferences on the model parameters : the posterior means were close to the true values and the empirical coverages of the credible intervals were close to the nominal values ; see @xcite . we compare three non - stationary brown  resnick models . the first , m1 , uses the same trend surfaces as does the best marginal model discussed in   [ sec : prelim ] and is fitted using pairwise likelihood as implemented in the r package spatialextremes ; this is our baseline model , to which we compare our new approaches . the second , m2 , is the bayesian hierarchical model that uses the fixed partitions from the declustering procedure . the third , m3 , is the bayesian hierarchical model with the random partitions . .estimated posterior mean parameters and @xmath125 confidence / credible intervals ( in parentheses ) from the trend surface model m1 , the bayesian model based on declustering m2 and the random partition bayesian model m3 . the second column is the estimated average location parameter @xmath165 over the 20 stations . [ cols="<,^,^,^,^,^,^",options="header " , ]     we discuss the estimates in terms of the negated location parameter @xmath141 , and negated trend @xmath166 , which correspond to the original scale of minimum temperatures . the estimates for the three models are given in table  [ tab : estimates ] . we assessed the fit of the marginal distributions using qq - plots ( see * ? ? ? we assessed the fit of the dependence model in terms of empirical pairwise extremal coefficients , computed using the empirical marginal distributions of the data to avoid potential bias due to the estimation of the marginals @xcite ; see fig .  [ fig : extremalcoef ] . as a global assessment of the fit in terms of the marginals and dependence , we computed the maximum of winter minimum temperatures over four different sets of stations , for the 53 years of the dataset , and compared the observed values with the predictions from the three brown  resnick models , computed by simulation using the algorithm of @xcite ; see fig . [ fig : modelchecking_max ] . further diagnostics for the prediction of the minimum and the averages of winter minimum temperatures are given in @xcite . confidence intervals in gray ) . curves are from the fitted brown  resnick models ( see table  [ tab : estimates ] ) : m1 ( dotted ) , m2 ( dashed ) , and m3 ( plain ) . ]    ): m1 ( top row ) , m2 ( middle row ) and m3 ( bottom row ) . the qq - plots compare the predictions and the observed maximum of winter minimum temperatures ( @xmath167c ) over some groups of stations , @xmath168 ( @xmath169 ) . the columns , from left to right , correspond to the group of stations @xmath170 ( western region ) , @xmath171 ( eastern region ) , @xmath172 ( northern region ) , and all 20 stations . grey dashed lines are @xmath125 overall confidence bands , and the solid grey line indicates a perfect fit.,scaledwidth=100.0% ]    for the trend surface model , m1 , the marginal qq - plots show poor fits for six stations . moreover , the model severely underestimates the dependence in terms of the pairwise extremal coefficient . this can be explained by the inadequacy of the fitted marginal distributions , which bias the dependence estimates ; @xcite also observed this when modeling wind gust data . in terms of predicting the maxima over groups of stations , the model performs poorly : many points lie outside the @xmath125 confidence bounds . in general , this model does not predict extreme temperatures well .    for the bayesian model based on the declustering , m2 , the marginal qq - plots show that the sample quantiles are outside the @xmath125 confidence bands for three stations , and underestimate them systematically at all stations ; the estimated scale parameter @xmath149 seems too large . the model overestimates the pairwise dependence somewhat , and tends to overestimate return levels for the maxima over different groups of stations .    for the random partition model , m3 , the marginal qq - plots show a good fit of the gev distributions : the empirical quantiles are in the @xmath125 confidence bands for all the stations . the empirical pairwise extremal coefficients agree adequately with the curve from the model . in terms of the multivariate predictions , m3 can predict the maximum , minimum and average of the minima better than does m2 , with all the observations within the @xmath125 confidence bounds ; see fig . [ fig : modelchecking_max ] and @xcite . this model appears to represent both the marginal and spatial behavior of extreme low temperatures adequately . we further discuss the results for the random partition bayesian model m3 . [ fig : mu_mean2016](a ) shows the posterior mean of the random effects for the location parameter over the convex hull of the meteorological stations ( see fig . [ fig : mapfinland ] ) , and the mean of the gev distribution for the year 2016 . we restrict our analysis to the convex hull of the stations because extrapolation outside this region is too uncertain . the posterior mean of the random effect in the location parameters shows the variability in this parameter that can not be explained by the covariates . the posterior mean for the fitted gev distributions shows lower temperatures in the inland region . the minima decrease as latitude , longitude , relative elevation and distance from the arctic ocean increase , and increase as lake cover and elevation increase . similar covariate effects were found by @xcite for extreme low temperatures in northern fennoscandia . the estimated shape parameter , corresponding to weibull  so - called light - tailed  distributions , agrees with previous studies on extreme temperatures . the dependence parameters correspond to strong dependence , with pairwise extremal coefficient @xmath173 at 400  km , while @xmath174 implies that the realizations from the brown  resnick process are continuous but nowhere differentiable . c).,scaledwidth=100.0% ]      we base our forecasts on the random partition bayesian model m3 . changes in extremes are entirely based on the linear trend @xmath166 in the location parameter . this translates to changes in the mean of the gev distribution at location @xmath144 and year @xmath37 , which is @xmath175 the estimate @xmath176 corresponds to an increase of @xmath177c per decade in the mean winter minimum temperature . fig . [ fig : mu_mean2016](b ) shows the predicted mean winter minimum temperature for 2016 . fig . [ fig : prob36 ] shows the probability that winter minimum temperatures are higher than @xmath0c , the threshold needed for protecting northern fennoscandia forests from outbreaks of _ epirrita autumnata _ ,  for 1980 , 2016 and 2030 . there is a strong change in these probabilities . over the region defined by the stations , the winter minimum temperatures had low probability of exceeding @xmath0c before 1980 , rising to about @xmath178 for 2016 , and then to @xmath179 for 2030 , according to the fitted model . the predictions must be interpreted with care , owing to the large uncertainty in the estimation of @xmath166 . c for 1980 , 2016 and 2030.,scaledwidth=100.0% ]    spatial simulation of winter minimum temperatures is useful for estimating joint return levels over certain regions , which could be useful for assessing the risks of pest outbreaks over a region for example . [ fig : simubr ] shows three simulations of winter minima fields from m3 , which show how the strong dependence in the brown  resnick model results in winter minima that are either low or high throughout the region . the spatial heterogeneity is mostly due to non - stationarity in the marginal distributions . conditional simulations are possible @xcite and could be useful as input to models to link past pest outbreaks with extreme temperatures , for example . c ) for 2016.,scaledwidth=100.0% ] although the stephenson  tawn approach was found to perform well on simulated data when the true partitions are known , the results of our application show that it can be inappropriate in practice . the hierarchical model for the marginals does not give a good fit to the marginal gev distributions , notwithstanding the additional flexibility provided by the random effects . we explain this by the use of the joint model for the minima and the partitions , and the declustering . the first of these requires good agreement of the data with the distribution of the block maxima and the partitions , and this is a strong assumption in practice . although the brown  resnick dependence structure is quite flexible for modeling stationary isotropic dependence structures , the distribution of the underlying partitions can be very different from the partitions obtained through the declustering . in the univariate case , previous studies have shown that declustering can bias estimators of the marginal parameters if the clusters are poorly - defined , and we can expect similar problems in the multivariate case . in our application , assuming that data separated by more than five days are independent in pairs results in partitions of small sizes each year ; see fig . [ fig : clust ] . a similar size for the true partitions for a brown  resnick process can only be obtained if the dependence in the process is very strong , and this is incompatible with the empirical extremal coefficients shown in fig . [ fig : extremalcoef ] . fig . [ fig : clust ] also shows the sizes of the random partitions in the mcmc algorithm for model m3 ; in several cases the size of the partitions from the declustering procedure have low probability for the random partition model . to further compare the partitions from the mcmc algorithm with that from the declustering procedure , we computed the rand index @xcite , which is , for each year , the proportion of pairs of locations that are simultaneously either together or separate in both partitions ; we then averaged the rand index over the years . it varies from @xmath180 to @xmath181 , indicating that pairwise classification agrees for less than 63% of the pairs between m2 and m3 . the size of the partitions is typically larger in the random partitions , and the relation with the occurrence dates is unclear . for example , in the winter 19731974 , the winter minima at the 16 stations with no missing data occurred on four consecutive days and the declustering procedure finds a unique partitioning set , whereas just 7% of the partitions visited in the mcmc run have a unique set and 37% have four partitioning sets or more . .,scaledwidth=100.0% ]    these results suggest that a different declustering procedure must be used . we varied the number of days used in our declustering scheme , and also considered another scheme that forces the stations that are at least 150  km apart to be in different clusters . the parameter estimates changed with the scheme but the diagnostic plots were never very satisfactory . we have shown how bayesian inference can be conducted for the brown  resnick process , either using the occurrence times of maxima and a declustering procedure , or by including the partitions in an mcmc algorithm . the partitions add powerful information into the model about the dependence structure . however our application shows that using fixed partitions can badly bias the estimation of the marginal and dependence parameters of a max - stable process . in a second approach we included the partitions in the mcmc algorithm without using information on the occurrence times , and this produced superior results . a more complex approach could use the occurrence times to construct a prior distribution on the space of the partitions , to add more information into the model and improve the efficiency of the estimation , while accounting for the uncertainty on the partitions , but this might be rather complex . @xcite proposed using second - order terms in the stephenson  tawn model to reduce potential bias due to the partitions , but this seems difficult in our case because the number of terms added to the likelihood would be huge and the mcmc algorithm would become even slower .    our bayesian model can be generalized in several ways . we focused on a stationary isotropic brown  resnick process for simplicity , and it appears to be sufficient for our application . a similar bayesian hierarchical model can be constructed for the extremal-@xmath37 process using ideas of @xcite , anisotropy would be easily incorporated using the idea of @xcite , and so perhaps could the work of @xcite on non - stationary dependence models for extremes . however , the generalization of the present method to more complex models and larger datasets may be limited by computational feasibility . monte carlo estimation of the likelihood is critical to keep the running times reasonable , and is relatively fast and accurate in low dimensions , but with more sites , more simulations may be needed to obtain reliable inference . the model predicts an increase of @xmath177c per decade in the mean winter minimum temperatures in northern fennoscandia . further work to validate the predictions from our model and to quantify the possible effects of global change on boreal forests would be valuable . the linear trend in time provides sensible short - term predictions , but for longer periods more complex approaches would be necessary , though limited by the availability of reliable data . one possibility would be to link the changes in extreme temperatures to mean temperatures and use predictions from global climate models under different scenarios as inputs to forecasting . this work was undertaken while e.t . was at epfl , switzerland , with support from the swiss national science foundation . e.t . and d.s.c.s work at csu was partially supported by us national science foundation grant dms-1243102 . the temperature data are from the finnish meteorological institute database for finland , and from the european climate assessment & dataset database ( http://www.ecad.eu ) for sweden and norway . this work was motivated by discussions with seppo neuvonen , who pointed out to us the relevance of winter minimum temperatures to insect outbreaks .
Abstract:

Answer: the brown  resnick max - stable process has proven to be well - suited for modeling extremes of complex environmental processes , but in many applications its likelihood function is intractable and inference must be based on a composite likelihood , thereby preventing the use of classical bayesian techniques . in this paper we exploit a case in which the full likelihood of a brown  resnick process can be calculated , using componentwise maxima and their partitions in terms of individual events , and we propose two new approaches to inference . the first estimates the partitions using declustering , while the second uses random partitions in a markov chain monte carlo algorithm . we use these approaches to construct a bayesian hierarchical model for extreme low temperatures in northern fennoscandia .    ,    ,    ,