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: resolving in time the energy relaxation in molecules both at the vibrational and electronic levels is theoretically and experimentally challenging . there are a variety of experimental techniques for probing molecular dynamics , whose utility can vary from one polyatomic molecule to another one @xcite . among the most commonly used pump - probe signals , we can highlight photoelectron energy @xcite and angular distributions ( pad ) @xcite , electron - ion coincidence in the perturbative regime @xcite that can provide molecular frame - pad @xcite , or in the non - perturbative regime like coulomb explosion @xcite or above threshold ionization @xcite . all these detections are based on collecting charged species . an alternative experimental technique developped in the last decade is high harmonic spectroscopy , based on the high order harmonic generation process . this method has shown a relevant sensitivity to nuclear dynamics @xcite and an inherent connection to the symmetries of molecular orbitals @xcite .    the high harmonic generation ( hhg ) process , which occurs when a molecule is submitted to a strong probe laser field , can be described in a first approximation by a three step model @xcite . a valence electron first tunnels out of the molecule under the influence of the strong field . this electron is then driven away by the electric field and can finally recombine onto its ionic core to emit extreme ultraviolet radiation . this process occurs every half period of the probe laser , so that the resulting radiation is constituted of odd high - order harmonics . the amplitude , polarization and phase of the xuv - vuv emission encode information about the electronic configuration of the target molecules . since this pump - probe signal is purely based on optical detection , its contrast can be drastically reduced due to the harmonic emission from unexcited molecules . this can be improved by using for instance polarization resolved spectroscopy @xcite or transient grating spectroscopy @xcite . in the latter case , a spatial grating of molecular excitation is created by optical interferences between two identical non - colinear pump pulses . the sinusoidal spatial modulation of the pump intensity across the molecular beam results in a partial diffraction of the high harmonic emission . this approach has been used in the present paper to revisit the relaxation dynamics in no@xmath0 molecules electronically excited in the first optically - allowed electronic state @xmath2b@xmath0 close to the first dissociation limit ( @xmath1 3 ev , fig . [ figsurfaces ] ) . the dipole of this transition belongs to the b@xmath0 irreducible representation of the c@xmath3 point group , namely the pump pulses preferentially excite no@xmath0 molecules with the c@xmath3 axis perpendicular to the pump polarization . by reducing the point group symmetry to c@xmath4 , the 1@xmath2a@xmath5 and 1@xmath2b@xmath0 states become the 1@xmath2a and 2@xmath2a states respectively that are vibronically coupled . we call these diabatic states 1 and 2 , respectively , as labelled in fig this strong coupling results in a high density of states ( 0.754 levels/@xmath6 ) @xcite and leads to the loss of both the electronic and vibration / rotation signatures . this has made of no@xmath0 a kind of hellish grail of molecular physicists . time - resolved photoionization has revealed unassigned regular oscillations with a 500 fs period for a pump centered at 400 nm @xcite .    on the contrary , a recent high harmonic spectroscopy experiment using @xmath7 nm pump pulses has shown a few oscillations of @xmath1 100 fs characteristic time that have been assigned to the @xmath2b@xmath0-@xmath2a@xmath5 vibronic coupling @xcite , in qualitative agreement with theoretical predictions @xcite . in the present paper , we reproduce this experiment and find subtle but significant differences in the femtosecond transient with respect to previous results @xcite . we investigate the influence of the pump intensity and wavelength , and of the rotational temperature of the molecular beam . simulations based on trajectory surface hopping ( tsh ) reveal that the diffracted signal from the excitation grating directly reflects the vibrational dynamics of the first excited diabatic state around the conical intersection . the remainder of the paper is organized as follows . calculations of the molecular dynamics are presented in section ii . the experimental results of high harmonic transient grating spectroscopy are given in section iii . theory and experiments are combined to deliver an interpretation to the observed dynamics in section iv . a and first excited state 2@xmath2a of no@xmath0 as a function of the bending angle and no strecth distance r@xmath5 . these diabatic states are referred to as 1 and 2 in the text , respectively . the second no distance @xmath8 is fixed to 1.19 .,scaledwidth=40.0% ]    no@xmath9 has received much attention because it is a stable and small molecule for which the effect of a conical intersection can be observed and calculated . most of the theoretical studies were focused on the analysis of the experimental spectra and were performed in a time - independent approach . more recently , time - dependent simulations of the intramolecular nonadiabatic dynamics were also performed @xcite . these studies were based on a diabatic representation of the electronic hamiltonian obtained from _ ab initio _ calculations @xcite or from the experimental spectrum @xcite . by propagating a quantum wavepacket or a swarm of surface hopping trajectories , the time dependent populations of diabatic or adiabatic states were calculated and regular recurrences were observed in some of these populations . for a range of total energy going from @xmath1 13000 @xmath6 up to @xmath1 28000 @xmath6 , the period of these recurrences is in the range 50 - 100 fs with a propensity to increase when total energy increases . a good agreement was observed between quantum time - dependent ( qtd ) and tsh results @xcite , indicating that the semiclassical approach is relevant for this system . since experiments are carried out around 3 ev , the comparison with theory requires accurate pes including the o+no dissociation channel . among the three existing models of pes@xcite satisfying these requirements , the most reliable one was built by kurkal , fleurat - lessard and schinke@xcite ( kfs ) from _ ab initio _ data calculated with multireference configuration interaction and a quadruple zeta basis set . this last work provides both the adiabatic and diabatic representations of the ground and the first excited pes . the bound levels of the ground electronic state were calculated up to @xmath1 7050 @xmath6 and a good qualitative agreement with the experimental spectrum was obtained . however the kfs diabatic representation is incomplete because neither the diabatic electronic coupling nor the nonadiabatic derivative coupling have been calculated . nevertheless the square of the electronic coupling @xmath10 can be extracted from the available diabatic ( @xmath11 and @xmath12 ) and adiabatic ( @xmath13 and @xmath14 ) energies by @xmath15\ ] ] let @xmath16 = @xmath17 be the antisymmetric stretch coordinate , @xmath18 and @xmath19 being the n - o distances . when @xmath16 = 0 , the molecule has the @xmath20 symmetry , the electronic coupling should be equal to zero and the adiabatic energies should be equal to the diabatic energies . if @xmath21 0 , the symmetry is reduced to @xmath22 , the electronic coupling can be different from zero and the following inequality must be fullfiled : @xmath23 ) is not always satisfied due to the small errors introduced by the fitting procedures . since this occurs in regions where @xmath10 is close to zero , we safely reduce it to exactly zero . the sign of @xmath10 can not be known from equation ( [ eq : v12carre ] ) , but at the first order , @xmath10 is proportional to @xmath16 @xcite . for this reason which is a particular case of a more general symmetry reason @xcite , @xmath10 takes the sign of @xmath16 . finally , let us remind two important energy quantities of the pes . the o+no dissociation energy of the ground electronic state is @xmath24 = 3.00 ev while the experimental value is @xmath24 = 3.11 ev @xcite . the difference of potential energy from the ground to the excited electronic state , both at equilibrium configuration , is @xmath25 = 1.31 ev , while the experimental value is 1.21 ev @xcite . although the nonadiabatic dynamics of a triatomic molecule can be simulated nowadays with a qtd method without difficulties , we have preferred the semiclassical approach . the obvious benefit of using semiclassical methods is the possibility to study large polyatomic systems for which quantum dynamics calculations are out of reach . validating the semiclassical approach for the small system no@xmath9 will pave the way to the study of larger systems . while tsh simulations of the dynamics of no@xmath9 have been already successfully compared to qtd simulations @xcite , the present work provides the first comparison of the former with experimental data . the tsh method used here is based on the tully s fewest switches algorithm @xcite and some implementation details can be found in a previous work @xcite . the starting point in the phase space and the starting time of the trajectories are obtained by applying the classical treatment proposed by meier and engel@xcite . the initial phase space distribution is obtained by a two - step procedure . first , the initial distribution for the vibrational ground state of the electronic ground state , represented by the wigner function @xmath26 , is filtered by the resonant one - photon transition @xmath27 where @xmath28 is the potential energy difference of the transition , @xmath29 the central frequency of the pulse , and @xmath30 the parameter of the pulse envelope defined by @xmath31 . second , the random starting time of each trajectory is drawn from the distribution function @xmath32 . the quantities @xmath29 and @xmath30 have been defined by considering a pump pulse centred at 400 nm and with a 30 fs fwhm gaussian envelope . figure [ fig_pop2 ] shows two different evolutions of the population @xmath33 of the excited diabatic state : the first one is calculated by applying only the first step of the initial procedure , and the second one is calculated with the complete initial procedure . as expected , using random starting time tends to smooth out the details of the dynamics . these calculations have been done by propagating a batch of @xmath34 trajectories . smaller batches of @xmath35 and @xmath36 result in no observable difference in the @xmath33 population behavior , except that some numerical noise increases progressively as the number of trajectories decreases . the influence of the pump wavelength was also investigated by computing @xmath33 for a laser pulse centred at 395 nm and at 392 nm . no significant difference could be observed with respect to the result at 400 nm shown in fig [ fig_pop2 ] . all three curves were almost indistinguishable until @xmath37 fs . diabatic population calculated with a batch of trajectories starting at the same time ( green ) or starting at gaussian random time ( red ) , in the case of a 30 fs pump pulse centred at 400 nm.,scaledwidth=50.0% ]     excitation probability.,scaledwidth=40.0% ]    , with the second no distance @xmath19 fixed to 1.19  . the snapshots are obtained by representing only two coordinates ( bending angle and @xmath18 ) and ignoring the third one ( @xmath19).,scaledwidth=60.0% ]    let us now examine the intramolecular nonadiabatic dynamics revealed by these calculations . figure [ figpopangles ] shows the evolution of the populations @xmath38 and @xmath33 on the ground and excited diabatic pes respectively , along with the evolution of the average bond angle on each pes . additionally , figure [ fig_snapshot ] displays some snapshots of the trajectories swarm on both diabatic surfaces at selected time steps . because the details of the dynamics are more easily observed in the case where random initial delay is not applied , the snapshots where obtained from a batch of trajectories starting at the same time on the excited pes . only @xmath35 trajectories are included in order to avoid overloading the figure . at the beginning , @xmath39 = 0 and @xmath40 = 1 . then the swarm goes downhill and encounters the crossing seam at @xmath41 fs . because the swarm is still quite close to a @xmath20 geometry , the electronic coupling is inefficient , with few trajectories jumping on the ground pes . the swarm reaches the turning point at small bond angles at @xmath42 fs , spreads along q@xmath43 before returning to the crossing seam . due to this q@xmath43 dispersion , the electronic coupling is effective at this second crossing of the seam . therefore , a large number of trajectories jumps on the ground state . the main effect is the drastic decrease of @xmath33 observed between @xmath44 = 30 fs and @xmath44 = 60 fs on figure 3 ( a ) . the variation of @xmath33 is now controlled by the bending motion of the swarm on the ground pes . this swarm moves to the o - n - o linear configuration ( @xmath44 = 60 fs ) and then goes back towards small angles , crosses again the seam and thus repopulates the excited pes : this is the first recurrence of @xmath33 at @xmath45 fs . another period of bending occurs , indicated by a second recurrence which is hardly visible around @xmath46 fs . then both swarms are completely scattered in the whole phase space and no more change of the populations can be seen .    in the adiabatic representation , the dynamics is somewhat simpler . the trajectories start on the excited adiabatic state . around t=10 fs , most of the trajectories jump on the ground adiabatic state and then remain in that state , in which a large amplitude bending motion occurs subsequently . thus , the same dynamics is seen as a variation of the population in the diabatic representation or as a variation of the average bending angle in the adiabatic representation . the time dependencies of the diabatic populations ( fig . [ figpopangles](a ) ) as well as the snapshots of the trajectories ( fig [ fig_snapshot ] ) will be used in section iv to analyze the results of the experiment presented in the following section . transient grating spectroscopy is widely used in non - linear optical spectroscopy to perform background - free measurements . in this configuration , the molecular medium is excited by two synchronized non - colinear beams which optically interfere , forming a grating of molecular excitation . this technique can be directly transposed to the extremely non - linear regime of high - order harmonic spectroscopy @xcite . the probe beam generates high harmonics in this structured medium . the harmonics are diffracted by the grating of excitation , leading to the appearance of one or several diffraction peaks around the main harmonic beam . in the present experiment , the fringes of the excitation grating being horizontal the diffraction peaks appear above and below the main harmonic beam . this diffraction pattern can be used to perform high contrast pump - probe measurements , since there is no signal in the direction of the diffracted peaks in the absence of molecular excitation . in addition , since the diffraction pattern is formed by the interference of emission from excited and unexcited sources , it encodes the phase difference between harmonic emission from excited and ground state molecules @xcite .    high harmonic transient grating spectroscopy is particularly valuable in a context where the pump can induce dynamics both through one - photon and multi - photon absorption . the spatial modulation of the pump beam intensity is sinusoidal . for a pure one - photon process , the excitation grating is also sinusoidal and only first order diffraction can be observed in the far field ( fig . [ figtgs](a ) ) . if the excitation results from a two - photon process , then the grating is determined by the square of the sinusoidal modulation of the pump beam , which contains the second harmonic of the exciting beam spatial modulation frequency . consequently , second order diffraction peaks appear ( fig . [ figtgs](b ) ) . note that second ( and higher ) order diffraction peaks can also appear in the case of a saturated one - photon transition : the saturation of the molecular excitation induces an anharmonicity of the grating which is reflected in the far field diffraction pattern ( fig . [ figtgs](b ) ) . the experiment was performed using the aurore ti : sa source at celia , delivering 7 mj , 27 fs , 800 nm pulses at 1 khz . the laser beam is split into a pump ( 20@xmath47 in energy ) and probe ( 80@xmath47 in energy ) beams . the pump beam is frequency doubled using a 200 @xmath48 m sh - bbo crystal , which provides wavelength tunability around 400 nm without reducing the temporal resolution ( @xmath1 35 fs corresponding to a 7 nm bandwidth ) . the pump is further split in two parts by a 50 - 50 beamsplitter . the 800 nm light is filtered out from the pump beams using two dichroic dielectric mirrors ( in each arm ) with high reflectivity around 400 nm and high transmission around 800 nm . an adjustable aperture on the pump beams is used to reduce their energy and to ensure that the excitation of the molecular medium will be uniform over the area probed by high - harmonic generation . the pump energy used in the experiment typically varied between 10 and 35 @xmath48j / p in each arm . the two pump beams are parallel and vertically shifted with respect to each other by 16 mm . the 800 nm probe beam is centered between the two pumps . the three beams are focused onto the gas target using 750 mm radius spherical silver mirror placed under vacuum in order to minimize self phase modulation in the entrance window of the vacuum chamber . the pump beams produce an interference pattern at focus with a fringe spacing of 9.4 @xmath48 m . the probe pulse waist is @xmath135 @xmath48 m . the linear polarization of the probe pulse can be continuously rotated relative to the s - polarized pump pulses . the typical probe intensity determined from the harmonic cut - off ( see fig [ figtgs ] ) is @xmath11.2@xmath4910@xmath50 w / cm@xmath2 . the high - harmonic emission is produced a few 100 @xmath48 m downstream from a 100 @xmath48 m nozzle gas jet ( general valve ) operated at 100 hz with a backing pressure of pure no@xmath0 at 780 mbar . to avoid the dimer formation , the nozzle temperature is fixed at 80@xmath52 where only @xmath1 10 @xmath47 species are n@xmath0o@xmath53 . the latter has an extremely low excitation probability at 400 nm ( @xmath54 ) @xcite , so that it does not contribute to the detected signal . the harmonic signal is dispersed horizontally by a grazing incidence spherical grating with variable groove spacing and imaged by a detector consisting of microchannel plates , a phosphor screen and a charge - coupled device camera . j / p , at a pump - probe delay of 100 fs and with a probe pulse polarized orthogonally to the pump pulses . at higher pump intensity ( 35 @xmath48j / p , ( b ) ) , second order diffraction peaks appear in the far field . they can result from two - photon transitions or from the saturation of the one - photon absorption , the two effects leading to anharmonicities of the near - field harmonic spatial modulation.,scaledwidth=45.0% ]      as a first step we studied the picosecond dynamics of the diffracted and undiffracted signal . similarly to ref . @xcite , we observed picosecond transients reflecting the one - photon dynamics of no@xmath0 with a 1.8 @xmath55 0.3 ps dissociation time . this indicates that the dominant process detected in our measurement is dictated by one - photon excitation .    in order to study the electronic relaxation of photoexcited no@xmath0 molecules shown in fig . [ figpopangles ] and [ fig_snapshot ] , we have performed pump - probe measurements in the first few hundred femtoseconds following excitation . the pump pulses were 20 @xmath48j each and centered at 400 nm . first we determine angle between pump and probe polarizations that optimizes the contrast . figure [ figfemto_angle ] shows the first order diffraction efficiency @xmath56 for harmonic 15 as a function of delay and pump - probe polarization angle . the diffraction efficiency is defined as the fraction of the harmonic light contained within a single diffraction order : @xmath57 monitoring this quantity rather than the raw diffracted signal enables lowering the effect of laser fluctuations and provides smoother signals . the diffraction efficiency shows a maximum around 120 fs . the contrast of this feature is maximized when the pump and probe polarizations are perpendicular to each other . all of the following discussion will therefore be with reference to this configuration . j each , centered at 400 nm and polarized orthogonally to the probe . the cross correlation signal ( panel ( a ) : dots for the experimental points and green line for the gaussian fit ) from harmonic 16 allows an accurate determination of the zero delay.,scaledwidth=50.0% ]    figure [ figfemto ] shows the evolution of the total signal ( blue ) and first order diffraction efficiency ( red ) as a function of the pump - probe delay ( 10 fs steps ) , for harmonics 15 to 21 . this result was obtained by adding up 5 consecutive scans , resulting in an accumulation over 6250 laser shots per delay . the zero pump - probe delay is determined from the even harmonics signal . these harmonics result from wave mixing between the pump ( 400 nm ) and probe ( 800 nm ) pulses , and thus provide a cross - correlation signal . the influence of the dynamics initiated by the pump pulse on even harmonics is taken into account to determine the zero delay . we assume that the even and odd harmonics are similarly affected by the change of the molecular geometry . the cross correlation signal is thus obtained by normalizing the signal of even harmonics by that of adjacent odds . the resulting signal has a symmetric shape , which can be fitted by a gaussian function to extract the zero delay . note that this normalization shifts the zero delay by + 10 fs . all harmonics experience a fast intensity decay within the first few tens of femtoseconds of interaction with the pump . after @xmath58 fs , the signal oscillates , presenting a maximum at 65 fs and a minimum at 130 fs . after that , it shows a slow increase , characteristic of the ongoing dissociation of a fraction of the excited molecules , since at 400 nm , with a 7 nm laser pump bandwidth , 25@xmath47 of excited molecules lie above the dissociation limit . the diffraction efficiency shows an inverted behavior , with maxima at 15 fs and 120 fs , and a minimum at 58 fs . the contrast of diffraction efficiency is better than that of the total signal , in particular for highest harmonics . in parallel to the harmonic signal measurements , we measured the ionization yield from the generation jet ( not shown here ) . the ion signal increases by a factor 1.3 between @xmath59 fs and 60 fs and is stable afterwards . overall , these results are similar to the previously reported observation @xcite , nonetheless they present a few differences . first , the diffraction efficiency measured here is 6 times lower than in the previous experiment . second , dynamics appear as well in the undiffracted harmonic signal and not only in the diffracted light . third , only two oscillations of the diffracted light are observed : a first sharp peak , and a second one with a steep leading edge and slow falling edge . last but not least , the exact timing of the oscillations is different . in order to investigate the possible origins of these differences , we conducted a systematic study of the transient .    when the pump energy is increased to 30 @xmath48j , second order diffraction peaks appear around the harmonics . figure [ figfemto_m1m2 ] compares the evolution of the first and second order diffraction efficiencies as a function of pump - probe delay . the dynamics observed for the second order diffraction is clearly the same as in the first order . this indicates that the second order diffraction peak originates from the saturation of the excitation rather than from a two - photon absorption process . indeed , with a @xmath2b@xmath60a@xmath5 absorption cross - section @xcite @xmath1 6@xmath4910@xmath61 cm@xmath2 and a pump waist size of @xmath62 @xmath48 m , saturation is expected as soon as the pump energy reaches @xmath125 @xmath48j per arm . for the first order , the results are very similar to those obtained at lower pump intensity in fig . [ figfemto ] . while the relative heights of the peaks at 20 fs and 130 fs slightly differ between the two measurements , their positions are identical . this shows that increasing the pump intensity does not affect the detected dynamics . we performed measurements varying the pump energy between 10 @xmath48j and 35 @xmath48j and did not observe significant changes in the temporal behavior of the first order diffracted signal . the first peak in the second order diffraction is slightly shifted to longer delays ( fig . [ figfemto_m1m2](d ) ) , which could reflect the dynamics of the saturation process . j each , are centered at 400 nm and are polarized orthogonally to the probe.,scaledwidth=50.0% ]        as the rotational temperature of the molecules may play a role in the measurement , in particular through the dynamics of the anisotropy induced by the pump excitation @xcite , or the rovibrational and rovibronic couplings which are especially significant at long distances , we have compared the transients measured with two different gas sources . the first source is that used for all results presented so far : a 100 hz general valve pulsed source with a 100 @xmath48 m nozzle diameter , backed by @xmath63 mbar of pure no@xmath0 , and heated at 80@xmath64c . in this configuration the laser is focused at a few hundred microns away from the nozzle . the second source is a 1 khz even - lavie valve with a 250 @xmath48 m conical nozzle , backed by a mixture of 16 bars of he and @xmath65 mbar of no@xmath0 , and heated at 120@xmath64c . in that configuration the laser is focused at @xmath66 mm away from the nozzle . we expect the rotational temperature to be significantly lower in the latter case . figure [ figtlambda](a ) compares the evolution of the diffraction efficiency for harmonic 17 using these two sources . while the ratio between the first and second maximum is different , which probably reflects slightly different laser pump conditions , the results are remarkably similar . last , we studied the evolution of the transient as a function of the pump wavelength , by varying it between 392 and 400 nm ( fig . [ figtlambda](b ) ) . the results are robust against such wavelength variations .    in summary , none of the external parameters like pump intensity , wavelength and rotational temperature , can explain the differences observed between the data shown here and the ones published in ref . one possibility to explain the absence of oscillatory structure in the undiffracted hhg signal and more persistent oscillations seen in the earlier work could be the noise level , the statistics being improved in the present work due to the higher repetition rate of the laser system . the experimental results show that the evolution of the diffraction efficiency presents robust characteristics : a maximum around 20 fs , a minimum around 60 fs , and a broader maximum around 120 fs . this behavior can be compared to that of the molecular dynamics described in fig . [ figpopangles ] . the minimum diffraction efficiency at 60 fs corresponds to a minimum in the population of the excited diabatic state , as well as to a maximum average bond angle of the molecule . in order to distinguish which of the two plays a dominant role , we performed simple calculations of the harmonic signal and diffraction efficiency . since the tunnel ionization probabilities and recombination dipole moments have been calculated only as a function of the bending angle @xcite , our model is one - dimensional : we only consider the influence of the motion along the bending coordinate by integrating the trajectory swarms over the stretching coordinates . the resulting evolutions are shown in fig.[figcomparison](a - b ) .        in a three - step picture of hhg , the harmonic signal is determined by the strong field ionization rate of the molecule , the propagation of the electron in the continuum and the recombination dipole moment of this electron to the ion core . the ionization rates and complex recombination dipole moments were calculated in ref @xcite for the fundamental and first excited diabatic states , for three bond angles of the no@xmath0 molecules ( 85@xmath64 , 102@xmath64 and 134@xmath64 ) , taking into account the different possible ionization channels . we use these results in our model and interpolate them to intermediate bond angles . the propagation of the electron in the continuum induces an additional phase shift in the harmonic emission from excited molecules with respect to ground state molecules , which corresponds to the phase @xmath67 accumulated by the molecule between ionization and recombination : @xmath68 @xcite , where @xmath69 is the difference between the ionization potential of the considered state and the unexcited ground state , and @xmath70 is the electron travel time , typically 1.2 fs for harmonic 17 for the experimental conditions . as the bond angle varies , the ionization potential changes because of the dramatic difference between the neutral and cationic potential surfaces . this can induce important variations of the harmonic phase . we calculate this phase by using the ionic potential curves from ref @xcite . as shown in fig . [ figiptau ] , the phases vary linearly as a function of angle within the narrow range explored by the wavepacket , namely 120 - 142@xmath64 for the ground electronic state and 90 - 110@xmath64 for the excited state ( shaded areas in fig . [ figiptau ] ) .     in units of @xmath71 , as a function of bending angle , for emission from pes 1 ( @xmath72 ) ( blue ) and pes 2 ( @xmath73)(red ) . the shaded areas correspond to the main part explored by the trajectories . , scaledwidth=50.0% ]    for each pump - probe delay @xmath44 , the harmonic signal is calculated by coherently summing the contributions from the different bending angles @xmath74 , for each molecular state . the undiffracted and first order diffracted signals are given by @xcite : @xmath75 where @xmath76 is the fraction of excited molecules , @xmath77 the harmonic dipole moment for the non excited molecules , @xmath78 the complex dipole moment for emission from surface @xmath79 for the excited molecules , and @xmath80 is the time - dependent bending trajectory packet on the @xmath79 surface as depicted in fig . [ figcomparison ] ( a - b ) . note that since the high harmonic signal is a coherent macroscopic emission , its amplitude is quadratically dependent on the population . the total signal and diffraction efficiency are @xmath81 and @xmath82 . figure [ figcomparison](c ) shows the results obtained by assuming that @xmath83 of the molecules are excited , together with the experimental results . the results of the calculations have been convoluted by a gaussian function of 15 fs fwhm to take into account the duration of the harmonic generation process . the total signal is normalized to 1 at negative delays . the agreement is remarkably good , especially for the contrast , considering the simplicity of the model used for hhg calculations and the fact that we neglect the stretching motion . both the fast dynamics and the average level of signals at long delays are well reproduced . besides the slight temporal shift in the position of the first peak in the diffraction efficiency , the main difference between experiment and theory lies in the width of the second oscillation : the experimental results show a single broad secondary peak around 120 fs , associated to a dip in the total signal , while the calculated signal shows a peak at 100 fs followed by weaker replicas at 170 fs and 240 fs . this is probably an effect of the spreading of the packet along the stretching coordinates , which smoothes out the evolution at longer delays .        in order to ascertain the origin of the transients observed on diffracted signals , we cancel out the effect of the population transfers by normalizing the packets at each delay : @xmath84 and @xmath85 this enables us to isolate the influence of the bending motion . the results are remarkably similar to those obtained with the full calculation ( fig . [ figanalysis ] ) . the influence of the packet dynamics can be further disentangled by calculating the contribution of each packet alone ( fig . [ figanalyzebending ] ( a ) ) , setting @xmath86 and @xmath87 ( blue ) or @xmath88 and @xmath89 ( red ) . this reveals that the main origin of the observed experimental modulations is the bending motion of the trajectory packet on the excited diabatic surface 2 . this is surprizing , given the fact that the extension of the bending motion is more important on the diabatic surface 1 than on 2 . as the packets move , the harmonic emission is mainly modulated by the phase accumulated by the molecules during the electron travel in the continuum , @xmath90 ( fig . [ figiptau ] ) . this phase , which appears in the complex dipole moments , typically varies between @xmath91 and @xmath92 during the packet motion on surface 1 , and between @xmath93 and @xmath94 during the packet motion on surface 2 . figure [ figanalyzebending](b ) shows the evolution of the diffraction efficiency as a function of @xmath90 . it is clear from this graph that the packet motion on surface 2 will induce a much stronger diffraction signal as well as a deeper modulation . this shows that similar variations of @xmath90 as a function of the reaction coordinates can have very different impact on the signal depending on the absolute value of the phase . for instance an additional recombination phase in the harmonic emission could strongly modify the situation and make the dynamics on surface 1 dominate the signal . this issue is for instance illustrated in the different contrast observed as a function of the harmonic order for the first and second order diffraction efficiency . the sensitivity of high harmonic spectroscopy to specific dynamics will thus depend on the considered system and must be carefully studied case by case . since the signal is dominated by the vibrational motion on the excited diabatic surface , one may wonder what is the influence of the conical intersection on this signal . to estimate this , we calculated the molecular dynamics on surface 2 in the absence of electronic coupling to surface 1 . the results are shown in fig . [ figcoupling](b ) and compared to the results including the coupling but keeping the population constant @xmath89 ( fig.[figcoupling](a ) ) . the results are similar until 20 fs delay . this means that the first passage through the conical intersection , around 10 fs delay , has only little effect on the structure of the trajectory packet . by contrast at 60 fs there is a strong transfer to the ground diabatic state ( fig [ fig_snapshot ] ) which alters the packet , depleting the highest bending angles . this is followed by a transfer back to the excited state around 100 fs . figure [ figcoupling](c ) compares the high harmonic diffraction efficiencies obtained with ( red ) and without ( blue ) the coupling . the results are quite different . in particular they indicate that the second oscillation observed around 100 fs in the experimental signal is due to the transfer through the conical intersection from the ground to the first excited diabatic state . while the change in the overall population of the excited state has only little effect on the signal ( fig . [ figanalysis ] ) , the packet structure and dynamical changes induced by transfers through the conical intersection are clearly reflected in the signal . the present results have to be compared to recent calculations based on wavepacket propagation @xcite . in a first 1d calculation , the bending wavepackets were modelled as gaussian functions of constant width , with their centers following the angular position of the full 3d wavepackets center . the harmonic emission was calculated using the same procedure as described here , i.e. taking into account the phase modulations induced by the vertical ionization potential changes , and the changes in the harmonic dipole moments due to the change in molecular geometry . this calculation resulted in strongly contrasted modulations of the diffracted signal as a function of delay . in a second calculation , the wavepacket motion was neglected and only the influence of the populations was retained . the resulting modulation of the diffracted light was weaker than in the 1d case . last , a 3d calculation taking into account the evolution of the molecular wavepackets along the three dimensions was performed , assuming constant harmonic dipole moments and phases except for from the @xmath95 phase . the results hardly show any modulation in the harmonic signal @xcite . since the model neglecting wavepacket motion and taking into account only the population evolutions gave the best qualitative agreement with their experiments , kraus _ et al . _ concluded that the experimental signal reflected the variation of the diabatic state populations , the phase modulations @xmath95 associated to the vibrational motion being quickly washed out by wavepacket spreading in the three dimensions . this conclusion is quite different from that drawn from our 1d model , in which we find that the population transfer only has a weak effect on the detected signal . it is interesting to note that for the 1d models , the proper treatment of the extension of the wavepacket along the bending angle ( as done in equation 9 and 10 in the present paper ) without the assumption of a gaussian shape as done in @xcite leads to a better agreement with experiments . if only the populations or vertical ionization potential were probed by high harmonic spectroscopy , the measurement should reveal the same dynamics irrespective to the pump - probe angle . this is not the case : the dynamics is much more contrasted for orthogonal pump and probes , as shown in fig . [ figfemto_angle ] . this is a consequence of the angular dependence of the harmonic signal with respect to the z - molecular axis : the harmonic dipole moments can indeed play a major role in the outcome of the measurement . the fact that the 1d model used here reproduces better the experiment than the 3d model previously used in @xcite assuming constant harmonic dipole moments for all bending angles and bond lengths suggests that the harmonic dipole moments present significant variations along the stretching coordinates , such that the complete 3d averaging is closer to the complete 1d calculation than to the partial 3d calculation . a complete 3d calculation , including the evolution of the harmonic dipole moments along all molecular coordinates , is required to confirm this statement . , @xmath87 , blue ) or surface 2 ( @xmath88 , @xmath89,red ) . the black dashed line is the result of the full calculation . ( b ) normalized diffraction efficiency as a function of the phase modulation @xmath90 in units of @xmath71 . the shaded areas correspond to the main part explored by the trajectories.,scaledwidth=50.0% ] high harmonic transient grating spectroscopy experiments in photoexcited no@xmath0 molecules show modulations of both the total harmonic signal and first order diffraction efficiency within the first 200 fs following the pump pulse . these can be remarkably well reproduced by combining trajectory surface hopping calculations of the reaction dynamics around the conical intersection and a simple model of high harmonic generation . the analysis reveals that the dominant information encoded in the transient grating signal is the molecular vibration on the first excited diabatic state . this sensitivity results from the modulation of the vertical ionization potential as the molecule vibrates , which induces a phase modulation of the harmonic emission . the vibrational motion in the ground diabatic state also induces a significant phase modulation . however this modulation is centered on zero ( the zero phase shift corresponding to the equilibrium geometry of the molecule in the ground state ) , such that it has only little effect on the measured signal . these results show the importance of the harmonic phase in the nature of the information encoded in high harmonic transient grating spectroscopy experiments . indeed , an additional phase shift in the harmonic emission could invert the situation and make the signal most sensitive to the vibration in the ground diabatic state . since the phase of the recombination dipole moments in hhg can show significant structures as a function of harmonic order @xcite , it would be interesting to extend this study to longer driving wavelengths which produce much broader harmonic spectra . the transient grating signal could then encode different dynamics ( ground diabatic state vibrations , first excited diabatic state vibrations , or population transfer through the conical intersection ) depending on the harmonic order , giving a full picture of the ongoing dynamics . we acknowledge r. schinke for providing the pess used in the tsh calculations . we thank r. bouillaud and c. medina for their technical assistance . we acknowledge financial support by the anr ( anr-08-jcjc-0029 harmodyn ) , the conseil regional aquitaine ( 20091304003 attomol and 2.1.3 - 09010502 cola2 projects ) , and the european union through the molcotuv and itn - iconic contracts .
Abstract:

Answer: we study theoretically and experimentally the electronic relaxation of no@xmath0 molecules excited by absorption of one @xmath1 400 nm pump photon . semi - classical simulations based on trajectory surface hopping calculations are performed . they predict fast oscillations of the electronic character around the intersection of the ground and first excited diabatic states . an experiment based on high - order harmonic transient grating spectroscopy reveals dynamics occuring on the same timescale . a systematic study of the detected transient is conducted to investigate the possible influence of the pump intensity , pump wavelength , and rotational temperature of the molecules . the quantitative agreement between measured and predicted dynamics shows that , in no@xmath0 , high harmonic transient grating spectroscopy encodes vibrational dynamics underlying the electronic relaxation .


Instruction: Write an abstract for the following acadamic paper.

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

Paper: this paper is devoted to the survival analysis of a system subject to competing failure modes within an external stressing environment . the external environment is assumed to stress the system at random and isolated times according to a random shock model . such a model can represent external demands e.g. , which put some stress on the system at their arrivals . shock models have been the subject of an extensive literature . following @xcite , shock models may be classified into different categories , according to whether arrival times and shock magnitudes are correlated @xcite , or independent @xcite , and according to the assumption put on the shocks arrival process : homogeneous or non - homogeneous poisson process @xcite , renewal process @xcite or non - stationary pure birth process @xcite . according to the influence of shocks on the system , shock models may be further classified into three types : extreme shock models ( a shock can cause the system immediate failure , see @xcite ) , cumulative shock models ( a shock increases some intrinsic characteristic of the system such as its deterioration , failure rate , age , number of already endured shocks , ... , see @xcite ) and mixed shock models ( a shock can either cause the system immediate failure or increase some intrinsic characteristic , see @xcite ) . more references can be found in @xcite ; @xcite ; @xcite ; @xcite .    in this paper , a shock model with mixed effects is considered , where the occurrence of shocks is classically modelled through a non - homogeneous poisson process and where each shock may result in the system immediate failure through a bernoulli trial , independent of the system intrinsic behaviour . apart from shocks , the system suffers from competing soft and sudden failures , where soft failures refer to the reaching of some given threshold for the degradation level , and sudden failures to accidental failures , characterized by a failure rate . a possible representation corresponds to a two - component series system , where the first component is subject to soft failures and the second one to sudden failure . the system may hence fail through three different competing modes : traumatic failure due to a fatal shock ; soft failure ; sudden failure . each non fatal shock induces some increase of both deterioration level and failure , with possible dependence in - between .    in the oldest literature , most models considered one single possible type of failures for the system : e.g. soft failures in @xcite or traumatic failures due to shocks in @xcite , both in the multivariate setting . more recently , different models have been developed , which consider two different types of failures . for instance , competition between soft and sudden failures is studied in @xcite . several application cases are proposed in the paper ( see also references therein ) . an industrial example is also provided in @xcite , which studies the reliability of an aircraft engine . note that @xcite assumes soft and sudden failures to be independent whereas the stressing environment of the present paper makes them dependent . competition between soft failures and occurrence of a traumatic event are considered in degradation - threshold - shock models ( dts - models , denomination of @xcite ) which have been proposed by @xcite and further studied by @xcite . a case study is provided in @xcite for the analysis of fatigue crack growth , where the effects of shocks on the degradation are put into evidence . note that , contrary to the present paper , dts - models consider some possible influence of the deterioration level of the system on the shocks arrival rate . however , the first shock of a dts - model is always fatal to the system ( leading to a single shock possibly endured by the system ) , whereas successive non fatal shocks are here envisioned . we could not find any paper which takes into account three competing failure modes , as in the present paper . however , based on the case studies from the previous literature , one can think that our model can reflect lots of systems subject to competing soft and hard failures , within a stressing environment . for instance , one can think of an aerial cabin hooked to a cable : the supporting cable is deteriorating due to corrosion and fatigue ; the linking pulley is subject to sudden failures ( and maybe also to deterioration ) ; both cable and pulley endure shocks at each cabin travel , which increase jointly their respective deterioration level and failure rate . the present paper consider several kinds of dependence between the three competing failure modes : at each non fatal shock , the increase of deterioration and failure rate is simultaneous . this induces a first type of dependence between soft and sudden failures . each shock may induce a failure , either because the shock is fatal , or because the deterioration is suddenly increased beyond the threshold level . this induces a second type of dependence between soft and traumatic failures , which may be simultaneous . also , some possible dependence is envisioned between the increments of failure rate and of deterioration at each non fatal shock . this induces a third type of dependence between soft and sudden failures . finally , following @xcite , the probability for a shock to be fatal depends on the shock arrival time , which induces a last type of dependence . up to our knowledge , all these kinds of dependence have not been yet considered altogether and , as will be seen all along the text , this model enlarges several ones from the previous literature . the paper is organized as follows : the model is specified in section [ model ] . the system reliability is computed through different methods in section 3 . sufficient conditions are provided in section 4 * *  * * for the system lifetime to be new better than used . the influence of various parameters of the shock environment on the lifetime is studied in section 5 . numerical experiments are proposed in section [ section num ] and concluding remarks end the paper in section [ concl ] . to make the model clear , a two - component series system is considered , where the first component is subject to sudden failure and the second one to soft failure . this two - unit system is just a representation for the competing soft and sudden failure modes , with no restriction . in the ideal condition ( for example in a laboratory environment ) , the lifetime of the first component is characterized by its intrinsic hazard rate @xmath0 , @xmath1 while the second one is subject to some accumulative deterioration modeled by an increasing stochastic process @xmath2 ( e.g. a gamma process ) . the second component fails once its deterioration level exceeds a failure threshold @xmath3 . the lifetimes of the two components are made dependent by their common stressing environment . this environment is modelled by a random shock process , where the shocks arrive according to a non - homogeneous poisson process @xmath4 with intensity @xmath5 ( or cumulated intensity @xmath6 ) . to avoid useless technical details , we assume that @xmath7 for all @xmath8 . more generally , one might consider that @xmath7 only for @xmath9 greater than some @xmath10 , which would mean that the shocks would arrive only after time @xmath11 . the points of the poisson process are denoted by @xmath12 , ... , @xmath13 , ... with @xmath14 almost surely . a shock at time @xmath9 may cause the system immediate failure ( fatal shock ) with probability @xmath15 $ ] , which depends on the age @xmath9 of the system at the shock arrival . a shock at time @xmath9 is non fatal with probability @xmath16 . a non fatal shock at time @xmath17 increases the deterioration of both components in a different way :    * for the first component , its hazard rate is increased of a non negative random amount @xmath18 , * for the second one , its accumulated deterioration is increased of a non negative random amount @xmath19 . the random vectors @xmath20 , @xmath21 , @xmath22 , ... are assumed to be independent and identically distributed ( i.i.d . ) with common distribution @xmath23 , and independent of the shocks arrival times @xmath24 ( and hence independent of the poisson process @xmath4 ) . at each shock , the increments @xmath18 and @xmath19 are possibly dependent . when subscript @xmath25 is unnecessary , we drop it and set @xmath26 to be a generic copy of @xmath20 . for @xmath27 , the distribution of @xmath28 is denoted by @xmath29 . we set @xmath30 to be the bivariate compound poisson process defined by @xmath31 with@xmath32where @xmath33 .    the processes @xmath34 and @xmath2 are assumed to be independent . provided that the system is functioning up to time @xmath9 , the random variables @xmath35 and @xmath36 stand for the cumulated increments on @xmath37   $ ] of the failure rate of the first component and of the deterioration of the second component due to the external environment , respectively . setting @xmath38 to be the @xmath39field generated by @xmath40 and provided that the system is still up at time @xmath9 , the conditional hazard rate of the first component given @xmath41 is@xmath42 and the conditional deterioration of the second component given @xmath41 is@xmath43    to make it clearer , we introduce @xmath44 , @xmath45 to be the lifetime of the @xmath46 component under the external environment , without taking into account the possibility of fatal shocks for the system . to simplify the writing , we denote by @xmath47 the conditional expectation @xmath48 for any measurable set @xmath49 , where @xmath50 stands for the indicator function ( @xmath51 if @xmath52 , 0 elsewhere ) . we then have:@xmath53 where @xmath54 is the cumulated intrinsic failure rate of the first component and where @xmath55 stands for the cumulative distribution function ( c.d.f . ) of @xmath56 . ( recall that @xmath57 is independent of @xmath41 ) . we now let @xmath58 to be the time to the first fatal shock for the system with @xmath59 and we assume that the bernoulli trials ( fatal shocks or not ) which happen at each shock arrival are independent one with each other , and that they depend on @xmath41 only through the @xmath60 s , that is : @xmath61   } \left (   t_{i}\right ) = \prod_{i=1}^{n_{t}}q\left (   t_{i}\right ) \label{tau3}\ ] ] where @xmath62 . the system failure is induced either by a fatal shock or by a component failure ( soft or sudden failure ) , whatever arrives first . the lifetime of the system hence is @xmath63    we finally make the additional assumption that @xmath64 , @xmath65 and @xmath58 are conditionally independent given @xmath66@xmath67    to sum up , the whole model is specified by :    * @xmath68 the ( generic ) random increments in failure rate ( first component , @xmath69 ) and deterioration ( second component , @xmath70 ) , * @xmath71 the intensity of the non - homogeneous poisson process , * @xmath72 the intrinsic failure rate of the first component ( sudden failure ) , * @xmath73 the intrinsic deterioration of the second component ( soft failure ) , * @xmath74 the probability for a shock at time @xmath9 to be fatal at the system level ,    ( plus some independence assumptions ) .    by taking special cases for these five ingredients , we can see that our model extents some well - known models from the literature .    for instance , taking @xmath75 ( no fatal shocks ) , @xmath76 constant , @xmath77 and @xmath78 all @xmath1 ( no second component ) , one gets the `` stochastic failure model in random environment '' from @xcite . taking @xmath77 and @xmath78 all @xmath1 ( one single component ) , one gets the `` stochastic survival model for a system under randomly variable environment '' from @xcite . taking @xmath78 all @xmath1 , @xmath79 , the model resumes to a classical extreme shock model ( one single component ) , where system failures are only due to shocks arriving according to a non - homogeneous poisson process , with probability @xmath80 for a shock to be fatal ( and @xmath81 to be harmless ) . this model is interpreted as the brown - proschan model by @xcite ; see also @xcite where various properties of the model are explored . taking @xmath75 ( no fatal shocks ) , @xmath82 ( one single component ) , @xmath83 ( homogeneous poisson process ) , @xmath78 all @xmath1 ( no intrinsic deterioration for the second - and single - component ) , one gets the `` cumulative damage threshold model '' from ( * ? ? ? * section 4 , case of i.i.d . damage increments ) . taking @xmath84 ( constant ) , @xmath85 ( one single component ) , @xmath78 all @xmath1 , one gets the `` cumulative damage model with two kinds of shocks '' from ( * ? ? ? * case of i.i.d . damage increments ) . taking @xmath86 , @xmath70 exponentially distributed , @xmath87 , all @xmath1 , one gets the model from subsection * 3.b * in @xcite . all these models are summed up in table [ table1 ] . note that we do not pretend at any exhaustibility and our model will include lots of other previous models which are not provided here . [ c]ccbrown - proschan model from @xcite & @xmath88 , + & @xmath78 , @xmath89 + deterministic boundary in @xcite & @xmath90 , @xmath87 + & @xmath70 exponentially distributed + @xcite & @xmath91 , @xmath92 is a constant , + & @xmath77 , @xmath78 , @xmath89 + @xcite & @xmath77 , @xmath78 , @xmath93 + cumulative damage threshold models & @xmath91 , @xmath77 , @xmath94 is a constant + in section 4 of @xcite & @xmath95 , @xmath78 , @xmath96 + @xcite & @xmath97 is a constant , @xmath98 , + & @xmath78 , @xmath89 + the objective of this section is to calculate the reliability @xmath99 of the system at time @xmath9 , with@xmath100 where we recall that the system lifetime @xmath101 is defined by @xmath102 . a first way to compute @xmath99 is to use classical monte - carlo simulations and to simulate a large number of independent histories for the system up to time @xmath9 ( * method 1 * ) . this method will serve as a comparison tool in the numerical experimentations in section [ section num ] . this method requires the simulation of a random variable with conditional hazard rate @xmath103 ( see algorithm [ algo ] in section [ section num ] for details ) and may imply long computational times for the system reliability . we provide below a few alternate methods which may be quicker and also easier to implement . [ method 2][prop1]the reliability is given by @xmath104 where @xmath105 is provided by @xmath106 and where @xmath107    due to the conditional independence of @xmath64 , @xmath65 and @xmath58 given @xmath41 , we have : @xmath108 using @xmath109 , we get:@xmath110 which provides @xmath111 and next @xmath112 , due to @xmath113 and @xmath114    based on the previous result , the only point to get the reliability @xmath115 is to compute @xmath116 . the remaining of the section is hence devoted to the computation of @xmath116 . starting from @xmath111 ( or @xmath117 ) , a possibility is to compute @xmath116  through monte - carlo simulations of @xmath118 and @xmath30 , which is simpler and quicker than simulating trajectories of the system according to the initial model . this method is called * method 2 * in section [ section num ] .    following @xcite and @xcite , one may also use a series expansion of @xmath116 , as provided by the following proposition .    [ method 3 , general case][prop2 ] @xmath119 where@xmath120 and @xmath121 are i.i.d . random variables with probability density function ( p.d.f . ) @xmath122}(x)$ ] and independent of @xmath123 . conditioning on the poisson process @xmath124 , we have:@xmath125 with@xmath126 now , given that @xmath127 , the conditional joint distribution of @xmath128 is the same as the joint distribution of the order statistics @xmath129 of @xmath130 i.i.d . random variables @xmath131 with p.d.f . @xmath132}(x)$ ] ( see @xcite e.g. ) . using the fact that @xmath133 is independent of @xmath124 , we get:@xmath134 noting that the expression within the expectation is invariant through permutation of the @xmath135 s , we derive that:@xmath136 which provides the result . based on the previous result , one can see that our model is equivalent to a classical shock model , where the shocks arrive according to a non - homogeneous poisson process with intensity @xmath137 and with conditional probability of survival at time @xmath9 equal to @xmath138 , given that there has been @xmath130 shocks up to time @xmath9 . [ method 3 , independent case][cor1]in the special case where @xmath69 and @xmath139 are independent , we get:@xmath140 with @xmath141 and @xmath142 where @xmath143 stands for the convolution operator , @xmath144 , all @xmath145 , and @xmath146 stands for the laplace transform of the distribution @xmath147 of @xmath148 , with@xmath149    starting from @xmath150 and using the independence of all @xmath151 s , @xmath152 s and @xmath135 s , and the identical distributions of all @xmath151 s and of all @xmath135 s , we get:@xmath153 with@xmath154 substituting this expression into eq . @xmath155 and next into eq . @xmath150 provides the result . [ ex1]let @xmath139 be an exponentially distributed with mean @xmath156 and @xmath148 an independent random variable . in that case , @xmath157 is gamma distributed with parameter @xmath158 . this provides@xmath159 using @xmath160 and eq . @xmath161 , we get:@xmath162 taking @xmath87 for all @xmath1 , @xmath163 and @xmath164 as a special case , we get@xmath165 and , for @xmath166@xmath167 using successive integrations by parts for the last integral . this last expression is the result of theorem 2 in @xcite , which hence appears as a special case of the previous results .    from proposition [ prop2 ] and corollary [ cor1 ] , one may derive the following approximation for @xmath168 . [ method 3 , approximation][cor method3]for @xmath169 , let @xmath170 be defined by @xmath171 in the general case , and by @xmath172 in case @xmath148 and @xmath139 are independent , where @xmath138 and @xmath173 are provided by @xmath150 and @xmath174 , respectively . then , for all @xmath1 , the sequence @xmath175 increases to the limit @xmath116 when @xmath176 and for all @xmath169 , we have @xmath177 where@xmath178 @xmath179 is defined by @xmath180 and @xmath181 is poisson distributed with mean @xmath179 . we just look at the general case . from proposition , we have @xmath182 due to @xmath183 , we get:@xmath184 based on the proof of corollary [ cor1 ] . this provides:@xmath185 and the result . the previous proposition provides numerical bounds for @xmath116 , which may be adjusted as tight as necessary , taking @xmath186 large enough . also , the required number of terms in the truncated series is given , to get a specified precision . this method is quite adapted as soon as it is possible to compute the @xmath138 s ( or the @xmath187 s ) . this mostly requires the distribution of @xmath188 to be known in full form , which is the case e.g. when the @xmath189 s are constant or gamma distributed . an example is provided in example [ ex1 ] , in the special case of an exponential distribution . in the most general case , the computation of the @xmath138 s ( or of the @xmath173 s ) may be as difficult as the initial problem of computing @xmath168 , so that the previous method is not always adapted . we finally present another method based on laplace transform , which does not suffer from the same restriction .    [ method 4][prop3]we have:@xmath190 or equivalently @xmath191 where @xmath192 is provided by its laplace transform @xmath193 with @xmath194 the bivariate laplace transform of the distribution @xmath195 of @xmath68@xmath196 and @xmath197 , all @xmath198.in the special case where @xmath148 and @xmath139 are independent , @xmath199 may be simplified into:@xmath200 where @xmath179 is provided by @xmath180 and where @xmath201 is the univariate laplace transform of @xmath202 .    remembering that @xmath203 , we get from @xmath117 that:@xmath204   ~dl\nonumber\\ &   = \mathbb{e}\left [   \left (   \int_{0}^{\infty}e^{-sl}f_{g_{t}}\left ( l - a_{t}^{(2)}\right )   ~dl\right )   e^{-\sum_{i=1}^{n_{t}}v_{i}^{(1)}(t - t_{i})}\prod_{i=1}^{n_{t}}q(t_{i})\right ] \label{lap fi}\ ] ] with@xmath205 because @xmath57 is non negative . setting @xmath206 , we obtain @xmath207 substituting this expression into @xmath208 provides:@xmath209 with@xmath210 \\ &   = \mathbb{e}\left [   e^{-\sum_{i=1}^{n_{t}}\left (   v_{i}^{(1)}(t - t_{i})+sv_{i}^{(2)}-\ln q(t_{i})\right )   } \right ]   .\end{aligned}\ ] ] because @xmath211 is equivalent to @xmath212 , we get:@xmath213 with@xmath214 noting that the sequence @xmath215 are the points of a poisson random measure @xmath216 with intensity @xmath217 , the function @xmath218 may be interpreted as a laplace functional with respect of @xmath219 @xmath220    the formula for laplace functionals of poisson random measures ( * ? ? ? * theorem 2.9 ) next provides:@xmath221 with@xmath222 substituting @xmath223 by its expression @xmath224 , we get:@xmath225   ( t).\end{aligned}\ ] ] substituting this expression into @xmath226 and next into @xmath227 provides @xmath228 with @xmath199 given by @xmath229 . equation @xmath230 is a direct consequence .    finally , in case @xmath148 and @xmath231 are independent , we have:@xmath232 and@xmath233 which ends this proof .    based on the previous result , one can compute @xmath116 by inverting its laplace transform with respect of @xmath3 . looking at equation @xmath234 , the key point is the inversion of the laplace transform @xmath199 . we next provide an example where the inversion is possible in full form . in the most general case , this can be done numerically using some laplace inversion software . let @xmath95 , @xmath78 ( all @xmath1 ) , @xmath94 and @xmath97 constant , @xmath235 identically exponentially distributed with mean @xmath156 ( so that @xmath148 and @xmath70 are completely dependent ) . then:@xmath236 where @xmath237 stands for the dirac mass at @xmath238 . we easily get:@xmath239 and@xmath240 for @xmath241 , we have @xmath242 where @xmath243 inverting the laplace transform @xmath199 , we obtain : @xmath244 as @xmath245 , we get the following full form for the reliability:@xmath246 where @xmath247 is the cumulative distribution function of a gamma distributed random variable with parameter @xmath248 .    in the special case where @xmath148 and @xmath139 are independent , the laplace inversion of @xmath199 is reduced to inverting @xmath249 , or equivalently to inverting @xmath250 , where @xmath251 is a constant . this is hence easier than in the most general case of correlated @xmath69 and @xmath139 .    to sum up the section , we have at our disposal four different methods for computing the reliability :    * method 1 * ( direct mc simulations ) : :    the main drawbacks of this method are that it suffers from long    computation times and that its implementation is less direct than for    the other methods . * method 2 * ( computing @xmath168 through formula @xmath117 and mc simulations ) : :    this method is much quicker and much easier to implement than method    1 . besides , it is always possible to use it . however , method 3 ( when    possible ) and method 4 are quicker . * method 3 * ( truncated series expansion + control of the truncation error through corollary [ cor method3 ] ) : :    this method provides very good results as soon as the    @xmath138 s ( or the    @xmath173 s ) are available in full form . * method 4 * ( laplace transform inversion ) : :    this method is the best when it is possible to inverse the laplace    transform @xmath252 in full form . numerical    laplace inversion also provides quite good results . let us recall that a random variable @xmath253 ( or @xmath254 ) is new better than used ( nbu ) if @xmath255 all @xmath256 . we here provide sufficient conditions under which @xmath101 is nbu . [ propnbu]assume that the intrinsic lifetimes of both components are nbu , which means that : @xmath257 where the second condition is true as soon as @xmath2 is a univariate non negative lvy process . then , @xmath101 is nbu if one among the two following conditions is satisfied :    1 . @xmath97 is non increasing and @xmath94 is constant , 2 . @xmath97 is constant and @xmath258 is super - additive ( @xmath259 , all @xmath260 ) . let us first note that , in case @xmath2 is a univariate non negative lvy process , we have:@xmath261 due to the independent and homogenous increments of @xmath2 for the third line . assumption @xmath262 is hence true . starting again from @xmath263 and based on the nbu assumption @xmath264 , it is sufficient to show that @xmath116 is nbu . now : @xmath265    due to the second nbu assumption @xmath262 .    under each of the two provided conditions , @xmath97 is non increasing so that @xmath266 , all @xmath267 . using @xmath268 and splitting the exponential and the product into two parts , one gets:@xmath269   .\label{edfi}\ ] ] setting @xmath270 then @xmath271 are points of the poisson process @xmath272 with admits @xmath273 for intensity . equation @xmath274 now writes:@xmath275   , \end{aligned}\ ] ] or equivalently:@xmath276   .\end{aligned}\ ] ] as @xmath277 is independent on @xmath278 and as the @xmath279 s are i.i.d . and independent on @xmath278 , one gets:@xmath280 with@xmath281 \\ b_{n } &   = \mathbb{e}\left [   f_{g_{s}}\left (   l-\sum_{j=1}^{n_{s}^{\left ( t\right )   } } v_{j+n}^{(2)}\right ) e^{-\sum_{j=1}^{n_{s}^{\left (   t\right )   } } v_{j+n}^{(1)}(s - t_{j}^{\left (   t\right )   } ) } \prod_{j=1}^{n_{s}^{\left ( t\right )   } } q(t_{j}^{\left (   t\right )   } ) \right ]   \text{.}\ ] ] noting that @xmath282 is independent on @xmath130 ( @xmath283 , all @xmath284 ) and that @xmath285 , we finally have@xmath286 where @xmath287   .\ ] ] the point now is to prove that @xmath288 under the two different assumptions .    1 . if @xmath94 is a constant , then @xmath289 is identically distributed as @xmath290 . we hence have@xmath291 and the result is clear . 2 .   if @xmath97 is constant , then:@xmath292   , \ ] ] where @xmath293 for @xmath294 . moreover , the respective cumulated intensities of @xmath290 and @xmath295 are @xmath296 and @xmath297 , with@xmath298 due to the super - additivity of @xmath296 . we derive from ( * ? ? ? * theorem 6.b.40 , example 6.b.41 ) that @xmath299 for all @xmath300 , where @xmath301 stand for the standard stochastic order . as @xmath302 is non decreasing with respect to each @xmath303 , we get that:@xmath304   \geq\mathbb{e}\left [   f_{n}\left (   \left (   t_{i}^{(t)}\right ) _ { i=1}^{n}\right )   \right]\ ] ] for each @xmath305 . setting @xmath306 , we derive by lebesgue s dominated convergence theorem that@xmath307   = \phi_{s}\left (   l\right )   \geq \lim\limits_{n\rightarrow+\infty}\mathbb{e}\left [   f_{n}\left (   \left ( t_{i}^{(t)}\right )   _ { i=1}^{n}\right )   \right ]   = \phi_{s}^{\left (   t\right ) } \left (   l\right )   , \ ] ] which achieves the proof . the conditions of the previous theorem means that    1 .   the probability for a shock to be non fatal decreases with time ( and @xmath94 is constant ) , 2 . the cumulated rate of shocks arrivals is larger at time @xmath9 than at time @xmath308 ( and @xmath97 is constant ) . such conditions hence mean that the environment is more and more stressing , or that it is more stressing after a while than at the beginning . such conditions are quite natural . we here study the influence on the lifetime @xmath101 of different parameters of the stressing environment : we study the influence of probability @xmath309 , of the dependence between @xmath148 and @xmath139 and of the cumulated intensity function @xmath258 . we hence look at the influence on the lifetime @xmath101 of all characteristics of the stressing environment which make the components dependent . the influence of @xmath310 is straightforward . we mention it for sake of completeness . let us consider two different systems with identical parameters except from @xmath310 ( first system ) and @xmath311 ( second system ) and such that @xmath312 for all @xmath313 . then , adding a tilde ( @xmath314 ) to any quantity referring to the second system , we directly get from @xmath111 that @xmath315 or equivalently that @xmath101 is smaller than @xmath316 in the sense of the standard stochastic order @xmath317@xmath318 as expected , the lifetime @xmath101 is hence increasing with the probability @xmath310 for a shock to be non fatal . we here study the influence of the dependence between the two marginal increments @xmath148 and @xmath139 on the lifetime @xmath101 . to measure the dependence level between @xmath148 and @xmath139 , we use the lower ( or upper ) orthant order , where we recall that @xmath319 is said to be smaller than @xmath320 in the lower orthant order ( @xmath321 ) if @xmath322 or equivalently if @xmath323    [ prop_dependence]let us consider two different systems , with identical parameters except from @xmath324 ( first system ) and @xmath325 ( second system ) . as previously , a tilde ( @xmath314 ) is added to any quantity referring to the second system . assume that @xmath326 . then @xmath101 is smaller than @xmath316 in the sense of the standard stochastic order @xmath327 . the point is to show that @xmath328 . starting from @xmath117 and conditioning on @xmath329 , we have @xmath330 where @xmath331 for all @xmath332 , all @xmath333 , with @xmath334 based on the independence between @xmath324 and @xmath124 , it is sufficient to show that @xmath335 all @xmath332 , all @xmath333 to get @xmath328 . as @xmath336 is non decreasing in @xmath337 and @xmath338 , we first know from ( * ? ? ? * theorem 6.g.3 . ) that @xmath339 all @xmath340 . as @xmath341 and @xmath342 are two sequences of independent random vectors , we derive from the same theorem that @xmath343 as both functions @xmath344 and @xmath345 are non decreasing , we derive ( same theorem ) that @xmath335 which achieves this proof . the previous result shows that the more dependent @xmath148 and @xmath139 are , the larger the system lifetime is . we finally study how the lifetime of the system depends on the frequency of shocks . [ influence_lambda]let us consider two different systems , with identical parameters except from @xmath258 ( first system ) and @xmath346 ( second system ) . as previously , a tilde ( @xmath314 ) is added to any quantity referring to the second system . assume that @xmath347 , and that @xmath97 is non decreasing . then @xmath101 is smaller than @xmath316 in the sense of the standard stochastic order @xmath327 . using a similar method as for the proof of theorem [ propnbu ] , we can write @xmath116 as @xmath348 where @xmath349 is non decreasing with respect to each @xmath303 , all @xmath340 . as @xmath347 , we derive in the same way that @xmath350 for all @xmath300 , which allows to conclude . this result is very natural . the more frequent the shocks occur , the shorter the lifetime is . as already mentioned in subsection [ rel ] , a first possibility for computing the system reliability @xmath99 is to use classical monte - carlo ( mc ) simulations ( method 1 ) and simulate a large number of independent histories for the system up to time @xmath9 . we here provide the algorithm that we have used , considering the case where @xmath351 , all @xmath1 and where @xmath148 is almost surely positive ( @xmath352 ) ( not essential assumptions ) . [ algo]repeat @xmath216 times with @xmath216 large enough :    1 . simulate @xmath57 with given distribution . 2 .   simulate @xmath353 according to the poisson distribution with parameter @xmath6 . 3 .   simulate @xmath353 i.i.d . random variables @xmath354 with p.d.f . @xmath122}(x)$ ] . the shock arrival times are given by @xmath355 where @xmath356 is the order statistics of @xmath357 . 4 .   simulate @xmath135 as the result of a bernoulli trial between a fatal ( 0 ) and a non fatal shock ( 1 ) at time @xmath17 , with probability @xmath358 for a shock to be non fatal , all @xmath359 . simulate @xmath353 i.i.d . random vectors @xmath360 , @xmath359 according to distribution @xmath195 . 6 . simulate the lifetime @xmath361 of the first component with conditional hazard rate @xmath362 given @xmath363 . with that aim , setting @xmath364 ( see @xmath365 ) , @xmath366 is a one - to - one increasing function from @xmath367 into @xmath368 and , for @xmath369 with @xmath370 , we have:@xmath371 it is then known that , if @xmath372 is uniformly distributed on @xmath373 $ ] , then @xmath374 is identically distributed as @xmath361 , see ( * ? ? ? * proposition 1.20 ) . 7 . compute @xmath375 where @xmath376 refers to the @xmath377th mc history , with @xmath378 .    at the end of the algorithm , symbol @xmath379 stands for the realization of a bernoulli trial @xmath380 between an up ( 1 ) or down ( 0 ) system at time @xmath9 , with probability @xmath381 for the system to be up . the reliability @xmath381 is then classically approximated by the empirical mean @xmath382 of the @xmath379 s and a 95% asymptotic confidence interval is computed . method 2 is based on mc simulations of trajectories of @xmath383 and of @xmath124 ( see section [ rel ] ) . for both methods 1 & 2 , mc simulations are based on @xmath384 histories . methods 3 & 4 are described in section [ rel ] . the four methods are compared on a few specific examples . in all these examples , we compare the reliability at time @xmath385 ( @xmath386 ) and we suppose that the shock are due to a homogeneous poisson process with parameter @xmath387 , @xmath388 , @xmath95 and @xmath2 is a null process . all other parameters are provided in table [ table4 ] , where @xmath389  means that the random variable @xmath390 is exponentially distributed with mean 1 . [ c]cccccinput & @xmath391 & method & @xmath386 & 95 % ci + @xmath392 & & 1 & 0.5196 & [ 0.5181 0.5212 ] + @xmath393 , @xmath394 & independence & 2 & 0.5195 & [ 0.5184 0.5205 ] + & & 3 , 4 & 0.5198 & + @xmath395 & & 1 & 0.5049 & [ 0.5033 0.5064 ] + @xmath396 & complete dependence & 2 & 0.5049 & [ 0.5039 0.5059 ] + & & 4 & 0.5054 & + @xmath392,@xmath397 & & 1 & 0.4809 & [ 0.4793 0.4824 ] + @xmath92 , @xmath398 independent & dependence & 2 & 0.4813 & [ 0.4801 0.4825 ] + @xmath393 , @xmath399 & & & & +    methods 1 and 2 may be used in any case . method 2 is more effective than method 1 ( shorter c.p.u . time and tighter 95% confidence interval - ic - ) . method 3 is more practical when @xmath92 and @xmath139 are independent with some specific distribution . method 4 is adapted to the case where @xmath92 and @xmath139 are dependent . we here illustrate several properties from a numerical point of view on a few examples . examples parameters are provided in table [ table3 ] . [ c]ccccccccc & @xmath3 & @xmath400 & @xmath57 & @xmath401 & @xmath94 & @xmath92 & @xmath139 & @xmath402 + ex.[ex9 ] & - & 0 & 0 & - & 1 & @xmath403 & @xmath403 & independent + ex.[ex8 ] & 2 & 0 & 0 & 1 & 1 & @xmath403 & @xmath403 & - + ex.[ex6 ] & 2 & 0 & 0 & @xmath404 & - & 1 & @xmath403 & independent +    [ ex9 ] this example illustrates the nbu property of the lifetime when @xmath94 is constant and @xmath392 is non increasing , see theorem [ propnbu ] . taking @xmath388 , fig . [ fig7_1 ] indeed shows that the remaining lifetime of a system with age @xmath405 is stochastically smaller than the lifetime of a new system . on the contrary , when @xmath94 is still constant but @xmath406 is non decreasing , the remaining lifetime of a system with age @xmath405 is not comparable with that of a new system , see fig . [ fig7_2 ] with @xmath407 . so the nbu property does not hold anymore in that case .    , nbu case ]    , not nbu case ]    [ ex8]we here consider two extreme cases for the dependency between @xmath148 and @xmath139 : independent or completely dependent ( here @xmath408 ) . the reliability in the completely dependent case is always greater than in the independent case ( fig . [ fig6 ] ) . this result is coherent with proposition [ prop_dependence ] .     and @xmath139 , example [ ex8 ] ]    [ ex6]this example shows the monotony of the reliability with respect to the intensity of the poisson process @xmath409 when @xmath97 is increasing , see proposition [ influence_lambda ] . the more frequently the shocks occur , the lower the reliability is ( fig . [ fig3 ] ) .    , example [ ex6 ] ] we here proposed a random shock model with competing failure modes , which enlarges several models from the previous literature . the model takes into account different types of dependence between competing failures modes , where the dependence is induced by a common external shock environment . the reliability has been calculated by several different methods and conditions have been provided under which the system lifetime is new better than used . due to this ageing property , it might be of interest to propose and study some maintenance policy to enlarge the system lifetime . several versions might be proposed , according to the available information .    also , the influence of the characteristics of the stressing environment on the lifetime @xmath101 has been studied . as expected , we saw that the lifetime was stochastically increasing with the probability @xmath310 for a shock to be non fatal . besides , and that result was not necessarily so clear at first sight , we saw that the lifetime was also stochastically increasing with the dependence between the two marginal shock sizes . finally , in case of a non decreasing function @xmath310 , we saw that the lifetime was stochastically decreasing with the cumulated frequency of shocks . this means that the more frequent the shocks occur , the shorter the lifetime is . this result is natural but our proof is limited to the case of a non decreasing function @xmath310 . in the special case of @xcite , the survival function of @xmath101 is however given by @xmath410 and it is easy to check that if @xmath411 ( stronger assumption than @xmath347 ) then @xmath412 , without any special condition on @xmath310 . so , the stochastic monotonicity of the lifetime with respect of the ( cumulated ? ) frequency of shocks might be true under a more general setting than in the present paper , without assuming any monotonicity condition on @xmath309 . we however have not been able to conclude on this point , and whether it is true or not remains an open question . both authors thank the referees for their carefull reading of the paper and their constructive remarks , which lead to a better introduction and justification of the model , and to a clearer paper . this work has been initiated during hai ha pham s phd studies in pau ( france ) , and has been supported by the conseil rgional daquitaine ( france ) . this work has also been supported for both authors by the french national research agency ( ammsi project , ref . anr 2011 bs01 - 021 ) . hao , c.  su , and z .- z . reliability analysis for mechanical components subject to degradation process and random shock with wiener process . in _ 19th international conference on industrial engineering and engineering management _ , pages 531543 , 2013 . a.  lehmann . joint modeling of degradation and failure time data . _ journal of statistical planning and inference _ , 1390 ( 5):0 1693  1706 , 2009 . special issue on degradation , damage , fatigue and accelerated life models in reliability testing .                c.  qian , s.  nakamura , and t.  nakagawa . cumulative damage model with two kinds of shocks and its application to the backup policy . _ journal of the operations research society of japan - keiei kagaku _ , * * 42**0 ( 4):0 501511 , 1999 .
Abstract:

Answer: a system is considered , which is subject to external and possibly fatal shocks , with dependence between the fatality of a shock and the system age . apart from these shocks , the system suffers from competing soft and sudden failures , where soft failures refer to the reaching of a given threshold for the degradation level , and sudden failures to accidental failures , characterized by a failure rate . a non - fatal shock increases both degradation level and failure rate of a random amount , with possible dependence between the two increments . the system reliability is calculated by four different methods . conditions under which the system lifetime is new better than used are proposed . the influence of various parameters of the shocks environment on the system lifetime is studied . * keywords : * reliability ; bivariate non homogeneous compound poisson process ; hazard rate process ; poisson random measure ; stochastic order ; ageing properties ; two - component series system .    * ams msc : * 60k10 & 60g51